ISRN Applied Mathematics

Volume 2013 (2013), Article ID 981692, 8 pages

http://dx.doi.org/10.1155/2013/981692

## Nonlinear Second-Order Model and Equilibrium Point Characteristic Analysis of DC Traction Power Supply System

^{1}School of Automation Science and Electrical Engineering, Beihang University, Beijing 100191, China^{2}College of Information Engineering, Beijing Institute of Petrochemical Technology, Beijing 102617, China

Received 21 April 2013; Accepted 9 June 2013

Academic Editors: E. Kita and H.-T. Yau

Copyright © 2013 Y. H. Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

According to the different features of the DC traction network transient current in the Metro, the electromagnetic transient process of traction power supply system that should be divided into the two forms of short-circuit fault and the low-frequency oscillation was proposed. While traction network short-circuit fault happens, the system model is a weakly nonlinear first-order differential equation, which feeder current is increasing exponentially and eventually stabilized. While the load current becomes low-frequency oscillation current, the system model is a nonlinear second-order differential equation. After linear processing in unique equilibrium point neighborhood, the state space equation of linear system is built, and the structural reason of low-frequency oscillation is revealed by eigenvalue analysis method. The simulation result shows that the simulation current waveforms are consistent with the recorded current waveforms and the division of transient process is reasonable.

#### 1. Introduction

Metro becomes the favor of big cities in the world because of its transport features of rapid, larger capacity, three dimensions, and environment protection effect of green, pollution-free [1]. Because Metro vehicle adopts electric traction, their safe operation state depends on stability of the DC traction power supply system (TPSS) and its fault protection measures. Traditionally, the main factor which influences Metro operation safety is the traction network short-circuit fault of high incidence and heavy harmfulness. However, with the continuous increase of transport capacity and the extensive application of regenerative braking technology, the low-frequency oscillation current of high incidence appears in traction network, which has the characteristic of faster start up and short duration, and is not characterized as fault. However, as the similarity of divergent oscillation current and remote short-circuit fault current in the amplitude and slope feature, the protection device often regards oscillation current as fault current and causes the protection system malfunction. Both determination of traction network fault type and reveal of oscillation generation mechanism rely on the analysis of electromagnetic transient process of TPSS; thus it is very significance to establish accurately system model and study its stability for Metro operation safety.

While the traction network short-circuit fault happens, the electromagnetic transient model of TPSS is presently the zero-state response of first-order *RL* series circuit [2], whose corresponding main protection algorithm is the combination algorithm of (current rate-current increment) [3, 4], thereby the feeder current increases by index law and reaches eventually a stable status during the traction network short-circuit fault. Against the fact that there exist oscillation phenomenon in the transient current of traction network, the transient model considering AC side of power supply effect to short-circuit current waveform is presented in [5], and the simulation waveform with “overshoot” phenomenon and the important conclusion of positive correlation between overshoot degree and fault distance are got. The transient physical model of TPSS that is a second-order *RL* hybrid circuit is presented in [6], the short-circuit current is provided from many nearby substations, and the same conclusion is concluded. However, it can be learned from actual recorded current waveform of Beijing Metro that the oscillation phenomenon does not appear in feeder current during short-circuit fault of traction network, while the obvious low-frequency oscillation process appears in normal operation stages, so the electromagnetic transient process of TPSS should be divide into the two forms of short circuit fault and the low-frequency oscillation.

It can be learned from the waveform analysis of short-circuit and oscillation current that the electrical drive structure of vehicle and operation status have great influence to oscillation process, so the transient model of TPSS should take vehicle factor into consideration. Firstly, the transient model of TPSS is established during the traction network short-circuit, and the fault type is determined by the feeder current analysis of transition process and change feature. Secondly, the nonlinear second-order dynamic system model considering vehicle factor and the corresponding linear system model in balance-point neighborhood are established, and the structural reasons of low-frequency oscillation is concluded by analysis of the state space matrix of linear system.

#### 2. Transient Current and Waveform Feature of Traction Network

The function of DC TPSS is to change three-phase AC voltage into DC 750 V (or DC 1500 V) voltage by step-down and rectifier and then provide power to vehicle by the overhead contact line (or the contact rail). The TPSS is mainly composed of DC traction substation, feeder, DC traction network (contact line and travel rail), and return line, and its structure features are bilateral power supply and shorter transmission line.

At present, the two types of transient current waveform of traction network short-circuit fault and system low-frequency oscillation are recorded in line 1 of Beijing Metro, as shown in Figure 1. In Figure 1(a), the protection started at 100 ms, and the fault is cut off before it reaches steady state, while in Figure 1(b), the protection started at 300 ms, and power supply is cut off before it reaches peak.

It can be seen from Figure 1(a) that the short-circuit fault current waveform has the following features: (1) the current pulsation amplitude is very small before short-circuit fault, while the current rises exponentially and reaches peak (or saturation value) in shorter time after fault; (2) the current has a certain initial value, which indicates that traction network is firstly in transmission condition and then happens to short-circuit fault, and its change rule does follow the zero-state response of first-order *RL* serials circuit. These two features have been confirmed by the recorded short-circuit fault cases of traction network of Beijing Metro in the last two years.

Figure 1(b) shows that the oscillation current has the following features: (1) the startup of oscillation is quick, the duration is generally less than 1 s, and the oscillation period is about 100 ms; (2) the start-up point at the negative current indicates that vehicle is in a regenerative braking state; (3) the amplitude of oscillation current is small, even smaller than the traction load current.

It can be concluded by the previous analysis that these two types of transient current have essential difference in waveform morphology and generation mechanism. Therefore, the modeling and analysis of TPSS will be carried out in accordance with two forms of traction network short-circuit fault and system low-frequency oscillation.

#### 3. Modeling and Analysis of TPSS during the Traction Network Fault

##### 3.1. Physical Model of Traction Substation

Generally, it is suitable that the equivalent modeling of traction substation adopts the theorem of L.C. Thevenin in study and analysis of TPSS electromagnetic transient process. The main electrical equipment of traction substation is 24-pluse rectifier unit (combination of traction transformer and rectifier), and the work feature of rectifier unit is that its external characteristic (function relation between DC output voltage and load current) is a weakly nonlinear relationship. What is more, the AC side equivalent inductance makes little effect on operation of TPSS because of large leakage reactance of rectifier transformer. Therefore, the equivalent physical model of traction substation is the series circuit of ideal voltage source and internal resistance , as shown in Figure 2.

Obviously, the absolute value of characteristic curve slope is the internal resistance , and the intercept of *y* axis is the voltage value of ideal voltage source , thereby the output equation of substation is

##### 3.2. System Modeling and Qualitative Analysis

The transmission lines of TPSS are uniform structure, their total length is generally less than 3 km, and voltage is low; thereby the equivalent resistance model of short-circuit line can be used in modeling of transmission line. Thus, the electromagnetic transient physical model of TPSS during traction network short-circuit fault can be got according to the model of substation and transmission line, as shown in Figure 3.

It can be learned from the circuit structure that the electromagnetic transient physical model of TPSS is a weakly nonlinear first-order *RL* series circuit, and the two series circuits are independent. In Figure 3, is the distance between substation A and substation B, is the distance between short-circuit point and substation A, is the arc voltage (only relevant to the arc length), voltage value of ideal voltage source is , and and are, respectively, feeder currents provided by substations A and B. In addition, if we assumed that the resistance and inductance of unit length for contact line and the resistance and inductance of unit length for travel rail are constant value, while the current changes, then the contact line impedance (, ) and travel rail impedance (, ) are proportional to distance of power supply.

According to the mesh current method and round direction shown in Figure 3, the first-order differential equations of feeder current based on small signal analysis method are where and .

The time-domain solution of (2) is

It is shown in formula (3) that the small signal of feeder currents and contains two components of steady-state and transient if the traction network short-circuit fault occurs in . That is, the feeder current contains steady-state component and transient components , and the feeder current contains steady-state component and transient component . As time goes on, the transient component and will attenuate in accordance with the law of exponential function, while the feeder current increases exponentially and eventually stabilizes to a steady-state current.

Through the real-time simulation of transient model of TPSS carried in the digital simulation platform of RT-Lab, the time-domain transition process of feeder current is obtained, and the simulation result is shown in Figure 4.

It can be seen from the simulation results that the simulation waveform is consistent with the recorded fault current waveform in change process; fault steady-state current value and first-order *RL* circuit response characteristics are decided by the value of substation internal resistance and traction network resistance at fault moment; and feeder current value is positive correlation to fault distance; the change rule of feeder current follows the zero-state response of first-order *RL* series circuit.

#### 4. Modeling and Analysis of TPSS during Low-Frequency Oscillation

Currently, the AC speed regulation train with VVVF function and EMU (electrical multiple units) is widely used in Beijing Metro, each train has usually six marshalling (4 EMUs + 2 trailers), and each EMU is equipped with two sets of AC speed regulation device with a filter capacitor of 12 mF and a filter inductance of 6 mH. In addition, all inverters and AC motors of vehicle can be equivalent to the series circuit of nonlinear negative resistance and DC voltage source.

##### 4.1. Physical Model of Low-Frequency Oscillation

The aim of regenerative braking of Metro is to achieve the effect of energy saving by converting the part of vehicle kinetic energy into electric power and feedback it to DC traction network. Therefore, the operating status of regenerative braking is that the traction motors are in running condition and the voltage of traction network is higher than the voltage of rectifier unit. Considering that the low-frequency oscillation of TPSS occurs in regenerative braking condition, obviously the vehicle is an important part of TPSS at this time. Thus, the electromagnetic transient physical model of TPSS during the low-frequency oscillation is established by the model of electric drive, substation, and transmission, as shown in Figure 5(a).

It can be seen from Figure 5(a) that the transient network physical model of TPSS during low-frequency oscillation is a nonlinear second-order *RLC* hybrid series circuit network and the DC voltage source of vehicle which contains a certain low-frequency component can be regarded as the system disturbance source. In the vehicle model part, is the filter capacitor, is the filter inductance, and is the equivalent nonlinear resistor of electric drive (traction converters and motors).

To simplify the calculation and analysis, the equivalent transformation of the transient model is carried, and the transformation results are shown in Figure 5(b). Thus, according to the contrast analysis of circuit structure, the quantitative relation of main circuit parameters is

##### 4.2. Determination of the Equilibrium Point

It can be concluded from Figure 5(b) that the system physical model is a nonlinear second-order *RLC *hybrid circuit after equivalent transformation. Therefore, the mathematical model of feeder current is a second-order nonlinear differential equation, which can be described by two first-order differential equations as

Because it is difficult to discuss the phase diagram on the entire phase plane for formula (5), the qualitative information is usually obtained by the discussion of local properties in the equilibrium point neighborhood, and the solution properties in the equilibrium point neighborhood can be studied by linearization method.

While vehicle is in the regenerative braking status, the operation condition of its electric drive system is very complex. If we assumed that the output power is constant and ignore the energy loss of electrical equipment, the voltage-current characteristic of nonlinear resistance is a function relationship of inverse proportion. Therefore, the external characteristics of converter can be described as where and are, respectively, the output voltage and the output current of converter, is a positive real constant, and is the output power of inverter.

In formula (5), the equilibrium points meet conditions of and , so the second-order nonlinear differential equations of the feeder current can be simplified as

Formula (7) is a monotonic function, and the voltage-current characteristic of nonlinear resistance is still a monotonic function after linear transformation; thereby these two monotonic functions have only one intersection point which indicates that the circuit working point is the only equilibrium point. Therefore, after the inductances are replaced with short circuit and the capacitors are replaced with open circuit in Figure 5, the equilibrium working point can be obtained using KVL; that is,

##### 4.3. Equilibrium Point Properties of the General Second-Order Dynamic Circuit

Generally, the equation of nonlinear second-order dynamic circuit can be described by formula (9) as follows

So, formula (9) in the equilibrium point neighborhood can always be transformed into

In formula (10), and are the functions that is composed of the quadratic term and the more order item, and

Without considering the higher-order terms, then

Obviously, the formula (12) is the linear equation of formula (9) in the neighborhood of equilibrium point, and are the equilibrium points of formula (12). If the matrix does not exist eigenvalues with zero-real part, the equilibrium point of non linear equation is called the hyperbolic equilibrium point. The trajectory shape in the hyperbolic equilibrium point has the same characteristics with the trajectory shape of linear equations (12) in the equilibrium point. Therefore, when the equilibrium point of the linear system is, respectively, a saddle point, node, and focus, the equilibrium point of corresponding nonlinear system is also, respectively, saddle point, node and focus, and the trajectory of nonlinear system and linear system are tangent at the equilibrium point.

##### 4.4. Properties of Equilibrium Point during System Low-Frequency Oscillation

According to the previous linearization method of second-order dynamic circuit and designating the traction current and the filter capacitor voltage as state variables, the state space equation of the formula (5) in the neighborhood of equilibrium point is

The compact form of (13) is

In formula (14), is the state matrix, and is the input matrix.

If we assume to be the initial value of formula (14) of , its solution can be expressed as

The linear system has only one equilibrium point , and the property of the equilibrium point is decided by the eigenvalue of matrix .

The characteristic equation is

The new variables and are defined as

Thus, the eigenvalues of formula (16) are

When the parameters of matrix change, the eigenvalues will have lots of combination forms, and the results are shown in Table 1 after being discussed.

##### 4.5. Analysis and Simulation of Divergent Low-Frequency Oscillation

Considering the malfunction of relay protect caused by divergent low-frequency oscillation, it is necessary to carry detailed analysis of divergent low-frequency oscillation.

When the divergent low-frequency oscillation occurs, the system equilibrium point is an unstable focus, and the eigenvalues are a pair of conjugate roots of positive real part, that is, , so the conditions of divergent oscillation are

Many electrical parameters are involved in formula (19), and most of them are related to distance of power supply, so the results are difficult to obtain by direct calculation. Now enumeration method is used to verify the structural conditions of divergent oscillation.

The output power range of Metro vehicles converter is 0~2.4 MW, the output voltage is 825 V, and then the change range of resistance is

In Figure 5, the filter inductor is 1.5 mH, and the filter capacitor is 48 mF (the four sets of electrical drive units are parallel). The loop electrical parameters of unit length , , , and are, respectively, 0.019 Ω/km, 1.139 mH/km, 0.016 Ω/km, and 0.328 mH/km if the actual average of resistance and inductance is adopted. Suppose that the power supply distance between two substations () is 2 km, the internal resistance of substation is 0.04 Ω, and fault distance is 0.5 km in representative instance, and substitute these electrical parameters into formula (4), then

Then, the parameters , , , , and are substituted into formula (19), and the result is

Therefore, while , the system configuration is unstable.

Frequency and the damping ratio (this case is negative and which show diverging effect) are the two important indexes in describing oscillatory waveform morphology, so the relationship between indexes and electrical parameters can be expressed as

Generally, there exists condition of in formula (23), so

Now, in order to clarify the feature of system divergent oscillation, the frequency and the damping ratio are calculated by an actual example. Assuming that the actual output power of the converter is 80% of the maximum power and the output voltage is 825 V, then

Substitute and other electrical parameters (, , , and ) into formula (23), and the result is

From the calculation results, we can conclude the following: (1) the change of system oscillation frequency is small, and its calculation value is basically consistent with the actual oscillation frequency; (2) the damping ratio is relatively large in the neighborhood of equilibrium point, which indicates that the start up of system divergent oscillation is fast.

In addition, as can be learned from qualitative analysis of formula (16) and (19), the equivalent nonlinear negative resistance is the fundamental reason of system structure instability. Considering that in system model, the electrical parameters that can be easily adjusted are filter capacitor and filter inductor , and the system structure stability can be improved by reducing and increasing . However, it is required to reduce the cut off frequency of low-pass filter by increasing in order to prevent vehicle disturbance into the TPSS, so the adjustable margin of structure parameters is very smaller.

Through the real-time simulation of low-frequency oscillation transient model of TPSS carried in the digital simulation platform of RT-Lab, the time-domain transition process of feeder current is obtained, as shown in Figure 6.

It can be seen from the simulation results that the simulation waveform is consistent with the recorded divergent oscillation current waveform in change process, the oscillation frequency is mainly determined by filter capacitor and filter induction , and the positive damping ratio is mainly affected by the equivalent resistance of electric drive.

#### 5. Conclusion

Currently, there exist two types of transient current in the DC traction network, namely, short-circuit fault current with increasing exponentially and low-frequency oscillating current with faster startup and short duration. Because the short circuit fault of traction network can make great harm to system, the action of short-circuit is characterized as a serious fault. The divergent low-frequency oscillation does not belong to short-circuit fault, but it will cause malfunction of relay protection because of the imperfect of electromagnetic transient model of TPSS. Therefore, the two types of transient current have essential difference in waveform feature and generation mechanism.

When the short-circuit fault of DC traction network occurs, the electromagnetic transient physical model of TPSS is a weakly nonlinear first-order *RL *series circuit, and the mathematical model of feeder current is a weakly nonlinear first-order differential equations based on the small-signal analysis methods. It can be learned from analysis of the feeder current time-domain solution that the feeder current is composed of steady-state component and transient components, and the steady-state component is independent of time, while the transient component will decrease exponentially with time. Therefore, the feeder current will increase exponentially and eventually stabilize to a steady-state current, which indicates that the change of feeder current is a monotonically increasing process. In addition, it can be learned from simulation experiment that the feeder current value is positivly correlated with fault distance, so the short-circuit fault type can be judged by the value of feeder current.

It can be learned from waveform feature of transient current that the system low-frequency oscillation appears in vehicle regenerative braking state, and the structure of the system itself is a first-order *RL* series circuit, hence the vehicle must be incorporated in the TPSS in transient modeling. The equivalent model of vehicle electric drive is a series circuit; of DC voltage source and non-linear negative resistance, so the electromagnetic transient physical model of TPSS with vehicles is a nonlinear second-order *RLC *hybrid circuit, and the mathematical model of feeder current is a nonlinear second-order differential equations.

The transient model of TPSS with vehicle is a complex nonlinear second-order dynamic system, so it is a relatively simple method in discussing the stability of system structure that the linear process of model is carried in the neighborhood of equilibrium point. The operating point is the equilibrium point of the system; thus the transient model can be transformed into a linear system in the neighborhood of operating point. As can be learned by analysis of the state space matrix analysis of linear system, the equivalent resistance of electric drive is the main electrical parameters that affect system instability and positive damping process, while the system oscillation frequency depends mainly on filter capacitor and filter inductor. In addition the frequency is the essential attribute of system structure, and its change is small, so the relay protection algorithm of DC traction network based on the frequency characteristic can effectively distinguish oscillation current and fault current.

#### Acknowledgment

This paper was supported by Scientific Research and Development Planning Projects “Study of Automation Systems Solutions and Technical Specifications in Traction Power Supply” in Ministry of Railways, China, under Grant 2007J022.

#### References

- X. Li and L. Zuo, “Energy spectrum and eigenvector of DC traction power supply short-circuit current,”
*Transactions of China Electrotechnical Society*, vol. 25, no. 11, pp. 164–168, 2010. View at Google Scholar · View at Scopus - J. C. Brown, J. Allan, and B. Mellitt, “Calculation of remote short circuit fault currents for DC railways,”
*IEE Proceedings B*, vol. 139, no. 4, pp. 289–294, 1992. View at Google Scholar · View at Scopus - B. Lang, “Research on feeder protection for DC traction power supply system,”
*Journal of Beijing Jiaotong University*, vol. 33, no. 5, pp. 65–68, 2009. View at Google Scholar · View at Scopus - G. F. Wang, Y. K. Sun, and K. H. Chen, “DDL protection in subway DC traction power supply system,”
*Proceedings of the CSU-EPSA*, vol. 19, no. 1, pp. 59–62, 2007. View at Google Scholar - K. Kongwei, Q. Qinlijun, Y. Yangqixun, and D. Dingfuhua, “DC side short circuit transient simulation of DC traction power supply system,” in
*Proceedings of the International Conference on Power System Technology (POWERCON '04)*, pp. 182–186, November 2004. View at Scopus - C. L. Pires, S. I. Nabeta, and J. R. Cardoso, “Second-order model for remote and close-up short-circuit faults currents on DC traction supply,”
*IET Power Electronics*, vol. 1, no. 3, pp. 348–355, 2008. View at Google Scholar