Next Article in Journal
A Quality Control Methodology for Heterogeneous Vehicular Data Streams
Previous Article in Journal
Learning to Reason on Tree Structures for Knowledge-Based Visual Question Answering
Previous Article in Special Issue
FMCW Radar Estimation Algorithm with High Resolution and Low Complexity Based on Reduced Search Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Objective Real-Time Tuning of SVC Used in Electrified Traction Systems

by
Mohammad Hossein Bigharaz
1,2,
Mehdi Amiri Dehcheshmeh
1,
Hadi Givi
3 and
Štěpán Hubálovský
4,*
1
Department of Electrical Engineering, Amirkabir University of Technology (Tehran Polytechnic), Tehran 159163-4311, Iran
2
Monenco Iran Consulting Engineers Company, Tehran 1994643315, Iran
3
Department of Electrical Engineering, Shahreza Campus, University of Isfahan, Iran
4
Department of Applied Cybernetics, Faculty of Science, University of Hradec Králové, 500 03 Hradec Králové, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 10 January 2022 / Revised: 11 February 2022 / Accepted: 11 February 2022 / Published: 17 February 2022

Abstract

:
Electric train system is a very large load for the power network. This load consumes a large amount of reactive power. In addition, it causes a massive unbalance to the network, which results in many problems such as voltage drops, high transmission losses, reduction in the transformer output ability, negative sequence current, mal-operation of protective relays, etc. In this paper, a novel real-time optimization approach is presented to adjust the static VAR compensator (SVC) for the traction system to realize two objectives; current unbalance reduction and reactive power compensation. A multi-objective optimization technique entitled non-dominated sorting genetic algorithm (NSGA-II) is used to fulfill the regarded objectives simultaneously. A comprehensive simulator has been designed for electric train network modeling that is able to adjust the parameters of SVC in an optimum manner at any time and under any circumstances. The results illustrate that the provided method can efficiently reduce the unbalancing in current as well as supply the demanded reactive power with acceptable precision.

1. Introduction

Nowadays, using railway networks, especially 2 × 25 kV autotransformer type electric traction systems, is increasing in some countries such as Spain, France, Belgium, Italy, Netherland, Russia, etc. [1]. So the power demand for the traction systems is also increasing, which in turn leads to more current unbalance as well as more reactive power consumption [2]. In [3], while the traction vehicle (with three-phase induction motors) is operated, the effects of the operation on power quality in a 110 kV transmission network were studied. Traction systems are generally fed by only two phases of the electrical network. This arrangement causes current unbalance and some problems to the network, such as increasing power losses reducing the efficiency and productivity of the generation, transmission, and distribution of electric energy [4]. Moreover, the current unbalance may lead to system instability [5,6]. To solve the mentioned problems, some solutions have been proposed, such as using Scott [7] and V-V configuration of the traction system transformer for current balancing [8,9], shunt reactors [10], and flexible AC transmission system (FACTS) devices for compensation of the reactive power. Static VAR compensator (SVC) is a widely-used FACTS device in this regard [11]. Reduction in the co-phase traction power supply (CTPSS) capacity and improving the power quality of freight railways can be realized through suppression of negative-sequence current and compensation of reactive power. To achieve these goals, voltage imbalance compensation was carried out in [12] for the stator of doubly-fed induction generators (DFIGs) utilized in power networks with high penetration of wind generators.
Many investigations in the literature have focused on SVC implementation for the traction system applications. In [5], a combination of thyristor-switched shunt reactor and current source converter was used for reactive power compensation in the traction system. In [13], SVC was used for current balancing through the elimination of negative sequence current. The effects of SVC implementation on the power quality improvement in the Shanghai electric railway network were analyzed in [14] using PSCAD software. In [15], a new method was proposed to predict the load current unbalance with negative sequence current in the AC traction network. Power quality problem in the DC traction system was analyzed in [16], and by using a modified SVC, the network distortion was reduced, which led to voltage profile improvement. Performance of the SVC in the traction substation of the Beijing–Shanghai electrified railway was studied in [17]. Characteristics of harmonic and negative sequence currents were analyzed in [18] to provide the optimum arrangement of SVC installation for low-voltage side of traction substation. Combination of SVC and hybrid power quality conditioner (HPQC) was proposed in [19] to improve the power quality compensation capability in AC electrified railway systems. In this technique, the reactive power demand of the load is met by SVC, while HPQC helps to supply the remaining required reactive power to mitigate the current unbalance problem. In [20], a new electromagnetic hybrid system and a magnetic static VAR compensator (MSVC) were provided to compensate for the power quality parameters of the network. In [21], the authors adopted an integrated compensation technique to mitigate power quality problems encountered in railway power distribution networks. In this approach, DQ decoupling technique was utilized to achieve a double closed-loop control strategy, by which harmonic current detection, as well as reactive power compensation, can be accomplished. In [22], the operation of a half-bridge railway power conditioner (HB-RPC) was analyzed for decreasing power quality problems including reactive power, harmonics, and negative sequence current in railway system. Moreover, a Genetic Algorithm (GA)-based control approach was adopted to achieve an optimized performance. In [23], a novel type of power quality compensation system with two sets of thyristor-controlled reactors and two sets of a thyristor-controlled third filter was presented. The authors claimed that the proposed power quality compensator has some features such as small capacity, low cost, reduced capacity of the active part, and favorable performance for compensating the negative sequence current. To overcome voltage dip effectively, a static VAR compensator was proposed in [24], which was paralleled with a 50 km section of a traction transmission line supplied by local substation. A nonuniform power factor partial compensation approach was presented in [25] to reduce the current rating of railway power conditioner (RPC). This purpose was realized by using particle swarm optimization (PSO) algorithm to specify lower boundary of the grid side power factor. In most of the mentioned works, the authors focused on only one objective for compensation and neglected other advantages.
The monitoring power supply of the electric traction system is an important concern for transportation infrastructures. In this regard, a monitoring system has been proposed in [26], which is capable of detecting short circuit fault as well as electro-erosion of the circuit breakers. This system acquires electrical parameters of the railway through a hardware-software structure combined with a global system for mobile communications (GSM).
In this paper, an SVC is used in the traction substation in order to improve two independent objectives simultaneously. The first objective is reducing the network current unbalance and the second objective is considered reactive power compensation. A multi-objective evolutionary optimization approach entitled NSGA-II is utilized for this purpose. After defining the desirable Pareto front, a Lagrange equation is utilized for selecting the final individual. As a case study, a simulator of the 2 × 25 kV autotransformer type electric traction system, including some moving trains with various operational modes (motoring, coasting, and braking), is used.
The rest of this paper is organized as follows. Section 2 presents a concise overview of the simple 2 × 25 kV AC traction system. In Section 3, the general configuration of the SVC is explained. Section 4 is assigned to the proposed multi-objective optimization. In Section 5, the simulation results of this investigation are presented. At the end, the main conclusions are given in Section 6.

2. Arrangement of AC Traction System (2 × 25 kV)

The general arrangement of 2 × 25 kV AC traction system is shown in Figure 1. This system includes a multi-transmission line (MTL), autotransformer (AT), traction substation (TSS), trains, and sectioning cabin (SC). The traction system is powered by two phases of the electrical network (e.g., 230 kV) and with a center-tapped transformer, which energizes the electric railway conductors means Catenary (+25 kV), Rail (0 V), and Feeder (−25 kV). The AT is installed at every 8–15 km of the MTL. This method of supply imposes a big unbalance on the network, and due to the moving trains, the system demands variable reactive power in each time period [27,28,29]. Every part of the electric railway system is described as follows.

2.1. Train Modeling

In this paper, the train is considered as a current source, and its value would be updated at every moment based on some parameters according to the following equation:
I T r = m o d e f v η e 1 + j tan ( ϕ ]
η = P m e c P e l e c
Newton’s second law can be applied to the train’s movement:
M e f f d 2 x d t 2 = f R G x , v
The parameters in the above equations are described as:
m o d e Control mode of train;
f Tractive effort (N·m);
v Speed of train (m/s);
η Factor of efficiency;
e Voltage of catenary (V);
ϕ Current and voltage phase difference;
P m e c Mechanical power (Watt);
P e l e c Electrical power (Watt);
M e f f Train’s mass (effective) (kg);
R Resistance of train (N);
g Gravity force (N);
θ track slope;
x Train’s position (m);
G Track resistance (N).
If e = 0 , the train is in a coasting state, and if 0 < m o d e 1 , the train is in powering and cruising states, and if 1 m o d e < 0 , the train is in the regenerative braking state. The maximum acceleration of a train can be defined by the following formula:
a m a x = g P m . μ
where Pm is the ratio of motoring axles, and μ is the frictional coefficient. In railway’s rolling stocks, μ is based on the rail frictional condition and is generally limited to (0.05, 0.5). It may be bigger than 0.5 while the rail is dry or sandy.
For train resistance, some investigations have proposed the empirical relationship (5):
R = A + B v + C v 2
where A (kN), B (kNs/m), and C (kNs2/m2) are known as Davis coefficients. These empirical coefficients are determined by the train’s manufacturer.
In this paper, the following state equations are considered for train’s movement under all possible conditions throughout the track:
x j v j = v j t 1 M e f f f m o d e t , v t R G
m o d e t 1   ,   ,   1

2.2. Multi-Transmission Line (MTL) Modeling

According to Figure 2, because of the existing many conductors, electrical modeling of the MTL will be difficult. Although various methods have been presented for MTL modeling, a reduction method provided by A. Mariscotti [30] is used in this paper. In this method, three main categories are determined based on the voltage and all conductors are divided into these categories.

2.3. Modeling of Autotransformer (AT) Substations

This is a special transformer with equal winding turns ratios (a = 1) and coupled windings with three terminals, which are connected to the conductors of catenary, rail, and feeder throughout the line. Using ATs brings many advantages to the electric railway system, such as: reducing the number of required traction substations along the line, minimizing the electromagnetic interfaces (EMIs) between power electronic equipment and nearby telecommunication circuits, reducing voltage drop of the line, and decreasing the current that may be strayed to the underground’s metal constructions of the rail lines, and eventually reducing the power loss. In this paper, an impedance model of AT is used. As seen in Figure 3a, by using Kirchhoff laws, the following impedance equation for AT can be determined.
V C V R V F = Z a Z a Z a Z a Z m 0 Z a 0 Z a I C I R I F
Z a = Z 1 2 Z m
In these equations, Z1 and Z2 are the leakage impedances of the windings ( a = 1 Z 1 = Z 2 ) and Zm is the magnetizing impedance.

2.4. Modeling of Traction Substations (TSS)

According to Figure 3b, traction substations energize the electric railway system. The transformer of this substation is a center-tab two-phase transformer, for which the primary side is connected to two phases of the electric network, and the secondary side with a 180-degree phase difference is directly connected to the electric railway system. The impedance equation of the TSS can be determined as follows:
I C I R I F = I S S 1 0 1 + ( Z S S ) 1 . V C V R V F
Z S S = y a + y b 2 y b y a y b 2 y b 4 y b 2 y b y a y b 2 y b y a + y b 1
y a = 1 Z 1 T + Z 2 , y b = 1 Z 2 + 4 Z e
I s s = V P T 2 a
Z 1 T = Z 1 2 a 2
VPT represents the voltage of the primary side, a is the transformer ratio, Z1T is the leakage impedance of the primary side, Z2 is the leakage impedance of the secondary side, and Ze denotes the grounding impedance.
After defining the impedance models of all parts of the systems, means trains, autotransformer substations, traction substations, and multi-transmission lines, by composing these models as seen in Figure 4, a dynamic impedance model of the electric railway system can be determined. This dynamic impedance model is a central core for electric load flow calculation of the simulator developed for electric railway system analysis. The developed simulator is briefly described in Appendix A and in Figure A1.

3. SVC Implementation

SVC could be considered as the most popular FACTS device with two specific properties: (a) despite all its constructed elements being passive, it can distribute active power between three phases correctly, (b) although SVC is made of the elements with constant impedance, it is able to establish variable impedance in the network [31,32]. Specifically, in the traction system applications, SVC is required in order to fulfill these purposes (a) load current balancing and power factor correction, (b) adjusting the system positive sequence current in an acceptable range, (c) controlling the negative sequence current magnitude in a selectable range [33,34,35].
From an economic point of view, traditionally, in order to improve the unbalancing in the current and also to compensate reactive power in a traction system, two separate pieces of equipment are needed. For example, for improving the current unbalancing, a specific device like Scott, V-V, or open delta transformer is utilized and for reactive power compensation, using a FACTS device is inevitable. Equipping a traction system with both devices is costly. In the current research, both aims, including reducing unbalancing in the current and reactive power compensation are realized with only an advanced SVC that has remarkably lower price than two separate equipment. It should be noted that, among FACTS devices, SVC has lower price compared with other complex FACTS devices in combined series-shunt controllers category.

3.1. Network Current Balancing

As mentioned already, the traction system has some trains with various operational modes. For the sake of simplicity, equations of the currents are determined in this part only for one train, but in the simulations, some trains with various conditions will be considered in the comprehensive simulator. All electrical parameters for SVC adjustment are determined by using an advanced load flow analyzer. According to Figure 5, the train has been located between phases A and B, then the SVC and the network currents are calculated as follows [36,37]:
I A , S V C I B , S V C I C , S V C = j 3 V 1 A ξ 0 ξ 1 ξ 5 ξ 3 0 0 ξ 3 ξ 5 B A B B B C B C A
I A I B I C = j 3 V 1 A ξ 0 ξ 1 ξ 5 ξ 3 0 0 ξ 3 ξ 5 B A B B B C B C A + I T r I T r 0
where IA,SVC, IB,SVC, and IC,SVC denote the 3-phase line currents of the SVC. BAB, BBC, and BCA are the variable susceptances of the delta-connected SVC, V1A is the positive sequence voltage of phase A, which is considered as the reference voltage, and ξ is equal to 1 < 30 . IA, IB, and IC represent the 3-phase network currents, and ITr indicates the train current. By decomposing the three-phase currents to positive, negative, and zero sequences, it can be shown that the load current includes only positive and negative sequence currents. The negative sequence current is a significant factor in the network current unbalance. Though, in order to balance the load current, it is necessary to limit the negative sequence current in a defined range. This procedure is carried out through the regulation of SVC susceptances for any condition of the system.
Assuming the train current to be I T r = I T r < φ and a = 1 < 120 , the following equations are used to obtain the triple sequences currents:
I 1 I 2 I 0 = 1 3 1 a a 2 1 a 2 a 1 1 1 I A I B I C
where I1, I2, and I0 are the positive, negative, and zero sequence currents, respectively. Subsequently:
I 1 I 2 I 0 = V 1 A ξ 3 ξ 3 ξ 3 ξ 5 ξ 3 ξ 0 0 0 B A B B B C B C A + 1 3 I T r 1 ϕ 30 1 ϕ + 30 0
By separating I1 and I2 into their real and imaginary components, the following equation is derived:
Re ( I 1 ) I m ( I 1 ) Re ( I 2 ) I m ( I 2 ) = V 1 A 0 0 0 0 0 1 1 1 0 3 2 0 3 2 0 1 2 1 1 2 0 B A B B B C B C A + 1 3 I T r cos ( ϕ 30 ) sin ( ϕ 30 ) cos ( ϕ + 30 ) sin ( ϕ + 30 )
Consequently, the negative sequence current (I2) can be controlled in a specified range. The above equation is used to obtain the SVC succeptances. The equivalent succeptances are a function of the trigger angles and can be expressed as [38]:
δ i 3 3 ! δ i 5 5 ! + δ i 7 7 ! = B i π ω t
where i   ϵ   A B ,   B C ,   C A and δ denotes the trigger angle of the thyristor.

3.2. Reactive Power Compensation

The reactive power needed for supplying the traction system changes alternatively. This is due to the variable nature of the load current based on the trains’ operating condition. For example, when many trains operate in the accelerating mode simultaneously, more reactive power is demanded from the network. This reactive power causes some difficulties such as power losses, limiting the network capacity, etc. Therefore, it is desirable to supply the required reactive power locally through implementing an SVC, as shown in Figure 6.

4. Multi-Objective Optimization

Generally, a multi-objective optimization problem (MOP) is considered as a set of ‘P’ objective functions f= (f1,…,fP) and a set of ‘U’ possible solutions in the decision space. Assume that Z = f (U) is an objective space in which every objective must be minimized. For every solution that means u U, an objective string z Z has been allocated. Z = {z1,…,zP } is calculated based on the function f: U→Z such that Z = f (U) = (f1 (u),…,fP (u)). The objective string z ϵ Z dominates z’ Z if zi ≤ zi, ∀ i {1,…, P} and also ∃ j {1,…,P}, zj < zj. The decision string u U dominates u’ U if f (u) dominates f (u’). On the other hand, a solution u U is named an efficient or Pareto optimal member if it is mapping on the objective space is a non-dominated value. The main goal is to identify a proper set ‘U’ of efficient solutions that satisfies a minimum deviation with the ideal Pareto front. Recently, evolutionary algorithms have been widely utilized for solving such problems [27].

4.1. Problem Definition

As mentioned in the previous parts, two objectives are tracked in this paper, network current balancing (f1) and reactive power compensation (f2), as expressed in (21) and (22). In this regard, some constraints are taken into account according to (23) and (24).
f 1 = ( I A + a 2 I B + a I C ) 2 + ( 1 I A ) 2 + ( 1 I B ) 2 + ( 1 I C ) 2
f 2 = Q d e m a n d Q S V C
λ < I 2 < λ
B m i n < B < B m a x
where Qdemand denotes the reactive power demanded by the traction system, QSVC is the reactive power generated by the SVC, I2 represents the negative sequence current, λ is the marginal amount of I2, and B is the SVC’s susceptance, which is restricted between Bmin and Bmax. These limits are defined according to the SVC’s construction.

4.2. Non-Dominated Sorting Genetic Algorithm

Non-dominated sorting genetic algorithm (NSGA-II) refers to a multi-objective optimization strategy that was introduced first by K. Deb in 2002 [39]. This algorithm is similar to the common genetic algorithm (GA) but with two additional functions entitled non-dominated sorting and crowding distance. First, the non-dominated individuals by others are located in the first group, then this group is deleted, and the above process is repeated for other individuals until all of them are rated. Crowding distance is a control parameter, which is calculated for each member of each group. This parameter indicates a degree of closeness to the other members in a group, as expressed in (25).
d j ( k ) = i = 1 p f i ( k 1 ) f i ( k + 1 ) f i max f i min
In (25), k is the member number existing in the jth group, and dj(k) is the crowding distance of this member. fi represents the ith objective function for a problem with P objective functions. Cross-over and Mutation are done as classical genetic algorithms.
Two significant properties of NSGA-II that lead to the selection of this algorithm in this investigation are running speed and also a great diversity of solutions along the Pareto front. The running speed of this algorithm is too high in comparison to the other multi-objective algorithms, and this property is obtained by the elitism performance of this algorithm. In this study, due to using the optimization strategy in a real-time manner, the running speed of the optimization algorithm is of great importance. Moreover, the selection strategy of the algorithm, which is non-dominated sorting of population, leads to a fine diversity of solutions on the Pareto front.
In NSGA-II, the individuals selected for the next iteration are based on the non-dominated sorting and the crowding distance parameters. More information in this regard can be found in [40,41]. The Block diagram model of the proposed system for regulation of SVC is shown schematically in Figure 7.
In real use of the proposed system, the initial data required for the system will be provided by some sensors, including voltage, current, and phase sequence sensors. These sensors will be located on the main busbars of traction substations and sense and provide the required data at every predefined moment for the optimization block. Finally, by determining the susceptance and regulation, the SVC will be tuned.

5. Simulation Results

Specifications of the traction substation (TSS), autotransformer (AT), and the trains in the simulation are presented in Table 1. In order to illustrate the current unbalance properly and to consider traction systems with considerable impact on the network, it is assumed that the traction system consumes 10 percent of the network power. This helps to validate the capability of the proposed SVC regulation strategy. In this paper, the priority of the SVC operation is assigned to the current balancing before power compensation. Therefore, in each time period (1 sec), the best solution is selected from the Pareto front using (26).
min ( λ 1 . O b j 1 i + λ 2 . O b j 2 i ) , i = 1 , , n p o p
where Obj1 indicates the first objective (f1 according to (21)) and Obj2 represents the second objective (f2 according to (22)). The coefficients λ1 and λ2 are weighting factors, which are applied in order to allocate preference between the objectives. Here λ1 is the weighting factor of current balancing and λ2 is the weighting factor of reactive power compensation, which are considered to be 0.9 and 0.1, respectively. It should be noted that λ1 and λ2 are generally determined by the managing parties of railway and electric network organizations. These factors may be changed during the day based on the parties’ decisions. In (26), npop represents the number of population members.
Due to real-time operation and for better convergence of the optimization algorithm, Equations (15)–(19) were used for generating the initial population in the decision space. As mentioned, two separated objectives were considered that should be minimized. For this purpose, NSGA-II was utilized to determine the Pareto front. Figure 8 shows the optimal solutions in the objective space for a specific moment (for example, the tenth sec). In the figure, the optimal solutions of the objective functions are identified. Users can change λ1 and λ2 coefficients to change the preference between the two objectives based on the circumstances. Table 2 presents some feasible solutions derived from Pareto front individuals in both objective and decision spaces.
It is necessary to note that sampling from the system was carried out every one second, and the optimization algorithm was always performed during a time shorter than one second. It can be seen from Table 2 that the total running time of the algorithm was 0.839 s, and this time satisfies the real-time requirement since the sampling time of the system is one second. This simulation was executed on a simulator of 2 × 25 kV autotransformer type electric traction system based on a timetable shown in Figure 9a,b. The simulation was performed with 250 sec time duration for various operational modes of some moving trains (Figure 9c,d) along a Multi-Transmission Line (MTL).
Some results such as network currents (before and after SVC implementation), susceptance, and reactive power profiles are illustrated in Figure 10, Figure 11, Figure 12 and Figure 13, respectively. Figure 10 illustrates the currents of the network before applying the SVC. It is obvious that the electric traction system imposes a massive impact on the upstream network, which leads to a large unbalance. Figure 11 shows the SVC capability for reducing the mentioned problem. Table 3 shows the maximum deviations in amplitude and phase from desirable values before and after SVC installation. In this table, AmA, AmB, and AmC denote the amplitude of the three-phase currents, and φAB, φBC, and φCA are the angle difference between the three-phase currents. Considerable mitigation of the unbalance in three-phase currents is observed due to the installation of SVC with proper parameters tuning. This leads to improvement in the power factor of the load as well as a reduction in the losses of the traction system. For instance, considering phase A, the percent of amplitude deviation is decreased from 23% to 1.8%. Moreover, a reduction in angle deviation from 8.5% to 1% is observed between phases A and B as a result of SVC installation. It should be noted that such deviations in amplitude and angle are very undesirable in high-load conditions.
As mentioned above, two independent objectives were considered. The first objective was mitigating the current unbalance that has been explained. The next objective was to provide reactive power to the electric traction system appropriately. Figure 12 shows that the SVC reactive power tracks the reactive power demanded by the traction system, as well. As shown in Figure 13, the values of the SVC succeptances were updated at every moment, which implies the system’s capability in adaptation to different conditions. This also indicates the flexibility of the proposed technique for tunning the SVC.

6. Conclusions

Due to imposing large current unbalance and consumption of variable reactive power, the traction system has a considerable impact on the whole network. In this paper, an SVC was considered in the traction substation to mitigate these problems. To achieve a favorable performance, it was proposed to determine the SVC optimum parameters using an intelligent optimization algorithm entitled NSGA-II. In this regard, a multi-objective optimization problem was defined mathematically to minimize network current unbalance and the difference between the demanded reactive power and the SVC reactive power. Two weighting factors were considered so that one can set the priority of the two objectives. Using the output solutions of NSGA-II, appropriate tunning of SVC parameters was realized. To evaluate the performance of the proposed approach, a comprehensive simulator for 2 × 25 kV AC traction systems was developed based on a mathematical model of different parts of the AC traction system. As the simulation scenario, a complete travel which means starting from a point and stopping at a destination by considering all movement modes (motoring, coasting, and braking), was allocated to all trains. The simulation results confirmed that even under the worst condition where all trains are in the accelerating mode, reduction in current unbalances and reactive power compensation can be realized effectively using the SVC tunned by the proposed approach. In other words, succeptances of the SVC are updated at every moment to adapt to the corresponding operating condition. Specifically, adjustment of SVC parameters can be performed such that the traction system is considered as a relatively constant load to the system. Moreover, it was shown that the performance of the proposed algorithm is real-time since the total running time of the algorithm is less than the sampling time. Due to the high running speed of NSGA-II in comparison with other multi-objective optimization techniques and the diversity of the solutions along the Pareto Front obtained by NSGA-II, optimum tunning of SVC can be achieved efficiently in a real-time manner to satisfy current unbalance reduction and reactive power compensation in a realistic traction system.

Appendix A

The following figure presents a view of the developed simulator. This simulator is sectionalized into several parts. The upper part is allocated for training settings. Electrical parameters curves can be shown dynamically in the monitoring section. As some load-flow methods are developed, the method for performing the load-flow analysis can be selected in the load flow section. The measurement part shows every electric parameter of the system. In the optimization section, objective functions and other settings are selected, and ultimately, the optimization process and results are displayed dynamically, which provides a fine scheme for analyzing the results of optimization while the simulation of the overall traction system is carried on.
Figure A1. Simulator window specifications developed for electric traction system [28].
Figure A1. Simulator window specifications developed for electric traction system [28].
Sensors 22 01584 g0a1

Author Contributions

Conceptualization, M.H.B.; methodology, M.A.D. and H.G.; software, M.H.B., M.A.D. and H.G.; validation, Š.H.; formal analysis, H.G.; investigation, M.H.B., M.A.D. and H.G.; resources, Š.H.; data curation, Š.H.; writing—original draft preparation, M.H.B. and M.A.D.; writing—review and editing, H.G. and Š.H.; visualization, H.G.; supervision, Š.H.; project administration, M.H.B.; funding acquisition, Š.H. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Excellence Project PřF UHK No. 2214/2022–2023.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data present in this study are available on request from the author M.H.B.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Givoni, M. Development and impact of the modern high-speed train: A review. Transp. Rev. 2006, 26, 593–611. [Google Scholar] [CrossRef]
  2. He, Z.; Zheng, Z.; Hu, H. Power quality in high-speed railway systems. Int. J. Rail Transp. 2016, 4, 71–97. [Google Scholar] [CrossRef] [Green Version]
  3. Milešević, B.; Uglešić, I.; Filipović-Grčić, B. Power quality analysis in electric traction system with three-phase induction motors. Electr. Power Syst. Res. 2016, 138, 172–179. [Google Scholar] [CrossRef]
  4. Ciccarelli, F.; Fantauzzi, M.; Lauria, D.; Rizzo, R. Special transformers arrangement for AC railway systems. In Proceedings of the 2012 Electrical Systems for Aircraft, Railway and Ship Propulsion, Bologna, Italy, 16–18 October 2012. [Google Scholar]
  5. Terciyanli, A.; Acik, A.; Cetin, A.; Ermis, M.; Cadirci, I.; Ermis, C.; Demirci, T.; Bilgin, H.F. Power quality solutions for light rail public transportation systems fed by medium-voltage underground cables. IEEE Trans. Ind. Appl. 2012, 48, 1017–1029. [Google Scholar] [CrossRef]
  6. Zhao, T.; Peng, X. The AC traction power supply system for urban rail transit based on negative sequence current compensator. In Proceedings of the 2018 International Power Electronics Conference (IPEC-Niigata 2018-ECCE Asia), Niigata, Japan, 20–24 May 2018. [Google Scholar]
  7. Barros, L.A.M.; Tanta, M.; Martins, A.P.; Afonso, J.L.; Pinto, J.G. STATCOM evaluation in electrified railway using V/V and Scott power transformers. In Proceedings of the International Conference on Sustainable Energy for Smart Cities, Braga, Portugal, 4–6 December 2019. [Google Scholar]
  8. Roudsari, H.M.; Jalilian, A.; Jamali, S. Flexible fractional compensating mode for railway static power conditioner in a V/v traction power supply system. IEEE Trans. Ind. Electron. 2018, 65, 7963–7974. [Google Scholar] [CrossRef]
  9. Tanta, M.; Pinto, J.G.; Monteiro, V.; Martins, A.P.; Carvalho, A.S.; Afonso, J.L. A Comprehensive Comparison of Rail Power Conditioners Based on Two-Level Converters and a V/V Power Transformer in Railway Traction Power Systems. In Proceedings of the 7th Transport Research Arena Conference, Vienna, Austria, 16–19 April 2018; pp. 1–10. [Google Scholar]
  10. Tümay, M.; Demirdelen, T.; Bal, S.; Kayaalp, R.I.; Doğru, B.; Aksoy, M. A review of magnetically controlled shunt reactor for power quality improvement with renewable energy applications. Renew. Sustain. Energy Rev. 2017, 77, 215–228. [Google Scholar] [CrossRef]
  11. Pavleka, J.; Nikolovski, S.; Marušić, A. Using a FACTS Device as a Power Conditioner Suitable for Dynamic Reactive Power Compensation in Railway Application. Teh. Vjesn. 2019, 26, 1686–1693. [Google Scholar]
  12. Xiao, D.; Chen, M.; Chen, Y. Negative sequence current and reactive power comprehensive compensation for freight railway considering the impact of DFIGs. CPSS Trans. Power Electron. Appl. 2021, 6, 235–241. [Google Scholar] [CrossRef]
  13. Guiping, Z.; Jianye, C.; Xiaoyu, L. Compensation for the negative-sequence currents of electric railway based on SVC. In Proceedings of the 2008 3rd IEEE Conference on Industrial Electronics and Applications, Singapore, Singapore, 3–5 June 2008. [Google Scholar]
  14. Meidan, C.; Mingli, W. The switching transient simulation of SVC used in electric railway traction substation. In Proceedings of the 2009 International Conference on Sustainable Power Generation and Supply, Nanjing, China, 6–7 April 2009. [Google Scholar]
  15. Chen, S.-L.; Li, R.-J.; Hsi, P.-H. Traction system unbalance problem-analysis methodologies. IEEE Trans. Power Deliv. 2004, 19, 1877–1883. [Google Scholar] [CrossRef]
  16. Martinez, J.; Ramos, G. Reactive power and harmonic distortion control in electric traction systems. In Proceedings of the 2010 IEEE/PES Transmission and Distribution Conference and Exposition: Latin America (T&D-LA), São Paulo, Brazil, 8–10 November 2010. [Google Scholar]
  17. Jianzong, M.; Mingli, W.; Shaobing, Y. The application of SVC for the power quality control of electric railways. In Proceedings of the 2009 International Conference on Sustainable Power Generation and Supply, Nanjing, China, 6–7 April 2009. [Google Scholar]
  18. Wang, H.; Liu, Y.; Yan, K.; Fu, Y.; Zhang, C. Analysis of static VAr compensators installed in different positions in electric railways. IET Electr. Syst. Transp. 2015, 5, 129–134. [Google Scholar] [CrossRef]
  19. Habibolahzadeh, M.; Roudsari, H.M.; Jalilian, A.; Jamali, S. Hybrid SVC-HPQC scheme with partial compensation technique in co-phase electric railway system. In Proceedings of the 2019 27th Iranian Conference on Electrical Engineering (ICEE), Yazd, Iran, 30 April–2 May 2019. [Google Scholar]
  20. Chen, B.; Zhang, C.; Zeng, W.; Xue, G.; Tian, C.; Yuan, J. Electrical magnetic hybrid power quality compensation system for V/V traction power supply system. IET Power Electron. 2016, 9, 62–70. [Google Scholar] [CrossRef]
  21. Bai, X.; Peng, J.; He, X. Design of Comprehensive Compensation System for Railway Power Distribution System. In Proceedings of the 2020 IEEE 9th International Power Electronics and Motion Control Conference (IPEMC2020-ECCE Asia), Nanjing, China, 31 May–3 June 2020. [Google Scholar]
  22. Amira, C.R.; Ramdane, B.; Hamza, B.; Meriem, A.; Imen, M. Genetic Algorithm Approach of HB-RPC for Adaptive Power Quality Improvement in Railway Traction Chain. In Proceedings of the 2021 International Conference on Electrical, Communication, and Computer Engineering (ICECCE), Kuala Lumpur, Malaysia, 12–13 June 2021. [Google Scholar]
  23. Fang, L.; Xu, X.; Fang, H.; Xiao, Y. Negative-sequence current compensation of power quality compensator for high-speed electric railway. In Proceedings of the 2014 IEEE Conference and Expo Transportation Electrification Asia-Pacific (ITEC Asia-Pacific), Beijing, China, 31 August–3 September 2014. [Google Scholar]
  24. Sekhar, G.C.; Kale, V.; Krishna, G.V. Application of DVR to improve voltage profile of indian railway traction system. In Proceedings of the 2014 6th IEEE Power India International Conference (PIICON), Delhi, India, 5–7 December 2014. [Google Scholar]
  25. Liu, L.; Dai, N.; Lao, K.W.; Song, Y. Nonuniform Power Factor Partial Compensation for Compensating Current Reduction Using Particle Swarm Optimization in Traction Power Supply System. IEEE Trans. Ind. Electron. 2021, 69, 6140–6151. [Google Scholar] [CrossRef]
  26. Andrusca, M.; Adam, M.; Dragomir, A.; Lunca, E. Innovative Integrated Solution for Monitoring and Protection of Power Supply System from Railway Infrastructure. Sensors 2021, 21, 7858. [Google Scholar] [CrossRef]
  27. Raygani, S.; Tahavorgar, A.; Fazel, S.; Moaveni, B. Load flow analysis and future development study for an AC electric railway. IET Electr. Syst. Transp. 2012, 2, 139–147. [Google Scholar] [CrossRef]
  28. Bigharaz, M.H.; Hosseinian, S.H.; Afshar, A.; Suratgar, A.A.; Dehcheshmeh, M.A. A comprehensive simulator of AC autotransformer electrified traction system. Int. J. Power Energy Convers. 2019, 10, 129–147. [Google Scholar] [CrossRef]
  29. Dehcheshmeh, M.A.; Hosseinian, S.H.; Bigharaz, M.H.; Mohseni, K. Analysis of lightning transient in 2× 25 kV AC autotransformer traction system. Int. J. Power Energy Convers. 2018, 9, 89–103. [Google Scholar] [CrossRef]
  30. Mariscotti, A.; Pozzobon, P.; Vanti, M. Simplified Modeling of 2x25-kV AT Railway System for the Solution of Low Frequency and Large-Scale Problems. IEEE Trans. Power Deliv. 2007, 22, 296–301. [Google Scholar] [CrossRef]
  31. Thomas, R. Reactive power compensation in electrical traction Using Active impedance concepts. In Proceedings of the 2013 International Conference on Circuits, Power and Computing Technologies (ICCPCT), Nagercoil, India, 20–21 March 2013. [Google Scholar]
  32. Saini, D.K.; Venkatasubramanian, B.V.; Yadav, M. Effect of Unaccounted Parameters on Reactive Power Compensation in Indian Electric Traction Line. IEEE Access 2020, 8, 182679–182692. [Google Scholar] [CrossRef]
  33. Chu, W.-S.; Gu, J.-C. A new hybrid SVC scheme with Scott transformer for balance improvement. In Proceedings of the American Society of Mechanical Engineers (ASME)/IEEE Joint Rail Conference, Atlanta, GA, USA, 4–6 April 2006. [Google Scholar]
  34. Lee, G.-S.; Hwang, P.-I.; Lee, B.-G.; Moon, S.-I. Loss Minimization of Electrified Railway Traction Systems Using SVC Based on Particle Swarm Optimization. IEEE Access 2020, 8, 219680–219689. [Google Scholar] [CrossRef]
  35. Lao, K.W.; Wong, M.C.; Santoso, S. Recent advances of FACTS devices for power quality compensation in railway traction power supply. In Proceedings of the 2018 IEEE/PES Transmission and Distribution Conference and Exposition (T&D), Denver, CO, USA, 16–19 April 2018. [Google Scholar]
  36. Vedam, R.S.; Sarma, M.S. Control of static Var compensators. In Power quality: Var Compensation in Power Systems, 1th ed.; CRC press: Boca Raton, FL, USA, 2008. [Google Scholar]
  37. Morais, V.A.; Afonso, J.L.; Carvalho, A.S.; Martins, A.P. New Reactive Power Compensation Strategies for Railway Infrastructure Capacity Increasing. Energies 2020, 13, 4379. [Google Scholar] [CrossRef]
  38. Wang, D.; Wang, J.; Yang, C.; Zhang, X.; Li, G. Research on Application of TCR+ FC Typed SVC in Power Quality Integrated Management for Power Traction System. In Proceedings of the International Conference on Sustainable Power Generation and Supply (SUPERGEN), Hangzhou, China, 8–9 September 2012. [Google Scholar]
  39. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  40. Li, Y.; Wang, S.; Duan, X.; Liu, S.; Liu, J.; Hu, S. Multi-objective energy management for Atkinson cycle engine and series hybrid electric vehicle based on evolutionary NSGA-II algorithm using digital twins. Energy Convers. Manag. 2021, 230, 113788. [Google Scholar] [CrossRef]
  41. Xing, C.; Li, K.; Zhang, L.; Li, W. Optimal compensation control of railway co-phase traction power supply integrated with renewable energy based on NSGA-II. IET Renew. Power Gener. 2020, 14, 3668–3678. [Google Scholar] [CrossRef]
Figure 1. General configuration of the traction system.
Figure 1. General configuration of the traction system.
Sensors 22 01584 g001
Figure 2. A typical view of MTL system (double lines): (a) conductors’ arrangement; (b) simplified impedance model.
Figure 2. A typical view of MTL system (double lines): (a) conductors’ arrangement; (b) simplified impedance model.
Sensors 22 01584 g002
Figure 3. (a) The impedance model of AT and (b) the impedance model of TSS.
Figure 3. (a) The impedance model of AT and (b) the impedance model of TSS.
Sensors 22 01584 g003
Figure 4. Block diagram model of the electric railway system.
Figure 4. Block diagram model of the electric railway system.
Sensors 22 01584 g004
Figure 5. SVC and single train view.
Figure 5. SVC and single train view.
Sensors 22 01584 g005
Figure 6. SVC location in the network.
Figure 6. SVC location in the network.
Sensors 22 01584 g006
Figure 7. Block diagram of the proposed system.
Figure 7. Block diagram of the proposed system.
Sensors 22 01584 g007
Figure 8. Determined Pareto Front.
Figure 8. Determined Pareto Front.
Sensors 22 01584 g008
Figure 9. (a,b): timetable of trains moving on two opposite rail lines considered in the simulation; (c,d): speed profiles of the trains moving on two opposite rail lines.
Figure 9. (a,b): timetable of trains moving on two opposite rail lines considered in the simulation; (c,d): speed profiles of the trains moving on two opposite rail lines.
Sensors 22 01584 g009aSensors 22 01584 g009b
Figure 10. Network currents without SVC.
Figure 10. Network currents without SVC.
Sensors 22 01584 g010
Figure 11. Network currents with SVC tunned by the proposed approach.
Figure 11. Network currents with SVC tunned by the proposed approach.
Sensors 22 01584 g011
Figure 12. Reactive power profile.
Figure 12. Reactive power profile.
Sensors 22 01584 g012
Figure 13. SVC’s succeptances profile based on the considered timetable.
Figure 13. SVC’s succeptances profile based on the considered timetable.
Sensors 22 01584 g013
Table 1. Specifications of the traction substation, autotransformer, and the trains in the simulation.
Table 1. Specifications of the traction substation, autotransformer, and the trains in the simulation.
TSSATTrain
Primary voltage230 kVRated Voltage25 kVRated Power4 MVA
Secondary voltage25 kVLeakage Impedances0.1564 +
j0.0997 Ω
Power Factor0.7
Rated Power16.5 MVARated Power5 MVAMass200 tons
Table 2. Some solutions derived from the determined Pareto front.
Table 2. Some solutions derived from the determined Pareto front.
Sol.Objective SpaceDecision SpaceRunning Time
Obj. IObj. II B A B B B C B C A 0.893 sec (in time interval (910) sec of simulator running)
10.0501350.3093360.1470190.1944590.228315
20.2453690.1823190.1820850.1977680.214491
30.1273700.2583600.1609460.1966070.222621
40.2523640.1777150.1855360.1976190.216171
50.0542330.3061070.1498500.1961430.229902
60.3120230.1398730.2011270.1974040.217196
70.3204610.1343440.1996560.1987890.213597
80.4657300.0390190.2232660.1997450.200516
90.3392990.1218270.1940020.1994690.203125
100.2717790.1659500.1936250.1972610.219731
Table 3. Maximum deviations in amplitude and phase of currents.
Table 3. Maximum deviations in amplitude and phase of currents.
Without SVCAmplitude
Deviation
Δ A m A −23%
Δ A m B −9.9%
Δ A m C −4.8%
Phase
Deviation
Δ φ A B −8.5%
Δ φ B C +11.6%
Δ φ C A + 4.1%
With SVCAmplitude
Deviation
Δ A m A −1.8%
Δ A m B −0.2%
Δ A m C +0.9%
Phase
Deviation
Δ φ A B +1%
Δ φ B C −0.5%
Δ φ C A +1.5%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bigharaz, M.H.; Dehcheshmeh, M.A.; Givi, H.; Hubálovský, Š. Multi-Objective Real-Time Tuning of SVC Used in Electrified Traction Systems. Sensors 2022, 22, 1584. https://0-doi-org.brum.beds.ac.uk/10.3390/s22041584

AMA Style

Bigharaz MH, Dehcheshmeh MA, Givi H, Hubálovský Š. Multi-Objective Real-Time Tuning of SVC Used in Electrified Traction Systems. Sensors. 2022; 22(4):1584. https://0-doi-org.brum.beds.ac.uk/10.3390/s22041584

Chicago/Turabian Style

Bigharaz, Mohammad Hossein, Mehdi Amiri Dehcheshmeh, Hadi Givi, and Štěpán Hubálovský. 2022. "Multi-Objective Real-Time Tuning of SVC Used in Electrified Traction Systems" Sensors 22, no. 4: 1584. https://0-doi-org.brum.beds.ac.uk/10.3390/s22041584

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop