An Integrated Energy-efficient Operation Methodology for Metro Systems Based on a Real Case of Shanghai Metro Line One

Metro systems are one of the most important transportation systems in people's lives. Due to the huge amount of energy it consumes every day, highly-efficient operation of a metro system will lead to significant energy savings. In this paper, a new integrated Energy-efficient Operation Methodology (EOM) for metro systems is proposed and validated. Compared with other energy saving methods, EOM does not incur additional cost. In addition, it provides solutions to the frequent disturbance problems in the metro systems. EOM can be divided into two parts: Timetable Optimization (TO) and Compensational Driving Strategy Algorithm (CDSA). First, to get a basic energy-saving effect, a genetic algorithm is used to modify the dwell time of each stop to obtain the most optimal energy-efficient timetable. Then, in order to save additional energy when disturbances happen, a novel CDSA algorithm is formulated and proposed based on the foregoing method. To validate the correctness and effectiveness of the energy-savings possible with EOM, a real case of Shanghai Metro Line One (SMLO) is studied, where EOM was applied. The result shows that a significant amount of energy can be saved by using EOM.


Introduction
Nowadays metros systems are playing important roles in our daily lives, especially in large cities.In developing countries like China, dozens of metro lines have been built in large cities and more are going to be built in the next few decades.For example, in Shanghai, the metro system has a total route length of over 500 km and carries over 8 million people each day.By the end of 2020, the total route length will be more than 1000 km [1].Since metro systems consume huge amounts of energy each day, saving even a few percent of the total energy can lead to great economic growth and environmental protection benefits.In the last few years, researchers have begun to pay more attention to finding the economic operation strategies for metro systems and have made some achievements.
Gonzá lez-Gil [2] provided a comprehensive appraisal of the main practices, strategies and technologies currently available to minimize the energy use of urban rail systems.In all these technologies, the two main methods are energy-efficient driving (eco-driving) and utilization of regenerative braking energy.The operation strategy of eco-driving is to optimize the speed profile of each train, while the utilization of regenerative braking energy is for one train to re-use the energy generated by the braking of other trains.The first topic has been well studied by many researchers [3][4][5][6][7][8][9][10][11][12][13][14][15], while the latter is a rather new approach.It is based on the practice that the acceleration process consumes energy from the power supply network whereas the braking process returns energy to the supply network.One way to use the returned energy from braking is to store the energy it generates.Up to now, several methods have been proposed to store the energy generated by braking, such as using supercapacitors [16][17][18][19][20][21], flywheels [22][23][24][25][26] and batteries [27][28][29][30][31].An exclusive review of these technologies can be found in [32].However, it is not easy to apply these methods in practice.On the one hand, high levels of capital investment are needed to build the new infrastructure.On the other hand, the implementation of the new infrastructure in and existing metro system is hard as well.
An alternative practical and economical approach to use the energy from regenerative braking in metro systems is to modify the driving strategy of each train in order to obtain an optimal energy-efficient timetable.Different from the forgoing energy-storage methods, a timetable optimization method does not incur any additional capital costs.It rearranges the departure and arrival times of the rolling stock so that the energy returned by the braking trains can be simultaneously consumed by those accelerating trains.This technology has been studied by many researchers.For example, Chen [33] optimized the timetable to reduce the maximum traction power.Nasri [34] got an optimal timetable for the test system and analyzed the effect of headway and reserve time.Albrecht [35] used running time modification instead of dwell time modification to reduce the power peak and energy consumption.Fournier [36] introduced a hybrid genetic/linear programming algorithm to reduce the energy consumption.Bocharnikov [37] set a weighted combination of energy consumption and journey time requirements as an objective function and changed the rates of train acceleration, braking and coasting to get an optimal driving strategy.Kim [38] developed a mixed integer programming model to minimize the maximum traction energy that occurs when the trains are running simultaneously.Peña-Alcaraz [39] maximized the use of regenerative-braking energy based on a power flow model of the electrical network and applied the scheduling strategy on Line Three of the Madrid underground system.A research group from Beijing Jiao Tong University has also contributed a lot to this issue in recent years.Yang [40] proposed a timetable optimization model to increase the utilization of regenerative energy and to shorten the waiting time of the passengers at the same time with a two-objective integer programming model.Li [41] combined the timetable optimization method and speed profile optimization method to achieve a better net energy consumption performance.Su [42] proposed a mathematical model to maximize the utilization of recovery energy.Then he also integrated the fleet size and cycle time decision, distribution of cycle time and driving strategy optimization to find a globally optimal schedule [43].
However, all the above researchers only considered how to generate an optimal timetable.None of them considered the solution on how to adapt the established timetable under disturbances, which is a very common situation, to a certain metro train.For example the doors of the train have to be reclosed if the train is too crowded or something got stuck in the doors.Although this kind of disturbance may only contribute to some short delay of the schedule, it may indeed change the subtle balance of the optimal timetable significantly, which will lead to extra energy costs.In this paper, a novel algorithm named Compensational Driving Strategy Algorithm (CDSA), which is motivated by the real demand of Shanghai Metro Line One, is developed especially for solving the disturbance problem.
This paper is arranged as follows: in Section 2, a new integrated Energy-efficient Operation Methodology (EOM) of metro systems is proposed.First, the timetable of the metro system is optimized to make full use of the regenerative braking energy.Then, the driving strategy of the metro train is changed automatically according to the CDSA under disturbance.In Section 3, a real case of Shanghai Metro Line One (SMLO) is studied and experiments were done to validate the correctness and effectiveness of the EOM energy-savings.A six-station, four-train and bidirectional pilot system, from Xujiahui Station to People's Square Station, is established in Simulink.Finally, the EOM is applied to the pilot system and the energy consumption is compared among different situations.The results show that 5.12% of the total energy is saved when there is no disturbance.Under disturbance conditions, 81.98% more energy on average can be saved with CDSA.

Timetable Optimization
Timetable optimization is a method that does not incur in additional costs when regenerative braking energy is used, compared with other energy-storage methods such as using supercapacitors or flywheels.The main idea is to synchronize the dwell time of each train at every station and make the braking and acceleration happen simultaneously so that the regenerative braking energy can be used to the greatest extent, which is illustrated in Figure 1.
Figure 1 shows the power consumption of two metro trains in the same metro system.The positive value of the power means the train is consuming energy (accelerating) while the negative implies the train is feeding regenerative energy back to the system (braking).If one train is braking while some other trains are accelerating at the same time, part of the regenerative energy can be used by the accelerating trains.The saved energy is highlighted in red in the figure.Otherwise, this part of the energy will be wasted on some resistors and transformed into heat.In the upper diagram it can be seen that only a little bit of the regenerative energy is used because the dwell time of the 2nd train is not optimized.However, if more dwell time is added to the 2nd train, illustrated on the lower diagram, the red part increases dramatically, which means that more energy is saved in the second situation.To get an optimal timetable, we set the dwell time of the train as decision variables.Suppose there are N periods of dwell time that can be changed for each train and there are M trains in the whole system, so there are K variables in total: In real situations, the dwell time can be modeled as two types: the fixed dwell time and the flexible dwell time.The former one, which is the necessary time for passengers to board the train, is fixed as a constant.The latter one, which has an upper bound of xmax, is set as the variables xi (i = 1,2,…,K).So the optimization problem can be formulated as follows: . .0 1, 2,..., K is integer 1, 2,..., K where E represents the total energy consumption of the system, Q is the number of substations of power supply for the metro system, tf is the total travelling time of the system, and Uj and Ij are the Direct Current (DC) voltage and current of the jth substation, respectively.
The optimization problem of Equation ( 2) is an Integer Programming Problem (IPP).In this paper, it is solved by the Mixed Integer-Laplace Crossover Power Mutation (MI-LXPM) genetic algorithm that was proposed in [44].It is a real coded genetic algorithm, which means the decision variables are encoded with real numbers instead of binary code.The whole algorithm is made up of Laplace crossover, power mutation, tournament selection, constraint handling techniques, and truncation procedure for integer restrictions.

Compensational Driving Strategy Algorithm
With the optimal timetable, a certain percentage of the total energy consumption can be reduced if all the metro trains strictly stick to it.However, this is not necessarily the case in real situations.When the trains are in service, there are a variety of disturbances that will keep the train off the established timetable, such as doors that cannot open or reopen/reclose for safety purposes.Normally this kind of disturbance does not last very long and can be fixed automatically.Nevertheless, since the energy-efficient timetable is arranged for the optimized simultaneous use of accelerating and braking among all the metro trains, such a time delay may cause a suboptimal use of the regenerative energy and influence the later operation of the disturbed train.Consequently, a little amount of time delay would significantly break the balance of the use of regenerative energy and reduce the energy-saving effect.In this section, a novel algorithm named Compensational Driving Strategy Algorithm (CDSA) is introduced to save the additional energy when disturbances happen.
First, the driving strategy of a single metro train has to been determined.By the application of Pontryagin's Maximum Principle, Horn [45] has proven that the optimal driving scheme consists of maximum five ordered phases: 1 Driving with maximal acceleration; 2 Traveling with constant speed (cruising); 3 Coasting; 4 Braking to target with maximal deceleration; 5 Waiting for passengers to board the train (dwelling).
Moreover, as for the specific metro system such as SMLO, only four phases of acceleration, coasting, braking and dwelling are needed in normal operation conditions because the metro system always has a relatively short station spacing (normally 1 to 2 km) [3,42,46].Note that this driving strategy modeling is a simplified model that is based on the real case of SMLO.Whether there is a cruising process or not will depend on the maximum speed and the station spacing of the line.If this maximum speed is not high enough or the station spacing is too long, it is necessary to travel at a constant speed to reach the next station or even re-motor the train in the middle of the journey.The schematic diagram of the driving strategy is shown in Figure 2.
In Figure 2, t1,n, t2,n, t3,n and t4,n represent the durations of acceleration, coasting, braking, and dwelling between the nth and (n + 1)th stations, respectively.a1,n, a2,n, and a3,n are the accelerations of accelerating, coasting, and braking between the nth and (n + 1)th stations, respectively.v1,n and v2,n are the instant speeds at the end of the acceleration and the coasting between the nth and (n + 1)th stations, respectively.In order to maintain the safety, the speed of the train at each stage cannot exceed the limit set by the red line.The blue line represents a typical Automatic Train Operation (ATO) system, whose function is to drive the train automatically according to a set of pre-programmed speed profile.As shown in Figure 2, the proposed driving strategy for a single train seems to be diverse.However, not all the variables can be changed for a certain metro train.In fact, a1,n and a3,n are always set to be a constant value in real situations, which is shown by the experiment in Section 3. The variable a2,n, however, has a tricky value.In fact, a2,n is determined by the running resistance Fr, which has the following equation: Equation ( 3) is called the Davis equation [47], where A, B, and C are coefficients specific to certain rolling stock.Usually B and C are much smaller than A, so we can also treat a2,n as a constant value of a2 for simplicity if the gradient is assumed to be constant (this also can be proved in the real case in Section 3).So, the driving process can be treated as three segments of uniformly accelerated motion.Thus, as long as v1,n and t2,n are determined, the travelling distance can be calculated, which means the distance Ln is the function of v1,n and t2,n.According to the basic formulation of dynamics, the relationship can be derived as follows: start, where tstart,n is the start time of the train leaving the nth station and can be calculated as: In addition, we have: Then substituting v1,n and v2,n in Equation ( 4) from Equations ( 6) and ( 7), and treating v1,n and t2,n as variables, we get Equation (8): where: Let's say that a disturbance causes a time delay at the nth station.For convenience, all the variables used in CDSA are symbolized the same as before except that there is a prime on each variable, e.g., the instant speed of v1,n is symbolized as v1,n' in CDSA.The very instinct idea of CDSA is to drive the train faster to make up the delay, which is shown in Figure 3.  Figure 3 shows two driving strategies between the nth and (n + 1)th stations.The first strategy, the original one, is depicted as the green line.The second strategy, which is shown by the red line, has a greater maximal speed and consequently takes less travelling time.Thus, the time difference between the two strategies can be used to make up the delay.However, according to Equation (4), the yellow area covered by the two driving strategies should be the same because the travelling distances of the two strategies are the same.
The next question should be foreseen and determined is how fast the train should be driven and how much time delay can be compensated at most.In fact, there has to be a tcmax,n that can compensate at most for the previous delay.Obviously, it is the case that the metro train can be driven at its maximal limited speed, so given the time delay of the previous station as td,n-1, and this time delay is to be compensated at the next station, we have: (11) Solving Equation ( 8) for t2,n, we have: where: In addition, according to Figure 2, we have: where tn is the total running time between the nth and (n + 1)th station.From Equations ( 11), ( 12) and ( 14), we can get the minimum travelling time allowed as follows: Then the maximal compensable time delay is: Note that if td,n-1 > tcmax,n, it is impossible to compensate the time delay within only one following station.The train should be driven at the maximal speed, leaving the uncompensated part of the delay as the new time delay between the nth and (n + 1)th stations and then try to compensate it in addition to the real possible time delay of this segment in the next period of travelling.Thus, the system can return to the original energy-saving status within a couple of stations.The mechanism of CDSA is shown in Figure 4. (a) State of optimal timetable (Ideal energy-saving status) The system returns to energy-saving status within two stations (c) Process of applying CDSA (CDSA returns the system to energy-saving status).

v1,n
Delay tcmax Delay Figure 4 illustrates the power consumption of two trains with compensation of time delay by CDSA. Figure 4a indicates the status where the metro system sticks to the optimal timetable, which is the ideal energy-saving status.It can be seen that energy is saved due to timetable optimization.Figure 4b shows the time delay caused by an unexpected disturbance.In this case a small amount of time delay eliminates the use of regenerative energy.Figure 4c represents the process of applying CDSA.When the delay happens, CDSA immediately takes effect and changes the speed of the train after leaving the nth station, which is indicated by the green line.In the first segment, v1,n of the second train is set to be vmax.In this way the travelling time between the nth station and the (n + 1)th station is saved at the most extent, which is used to compensate for the time delay.In the next segment, same thing happens except that v1,n+1 is set to be a certain value of v1,n+1' instead of vmax according to CDSA.By the end of the (n + 1)th station, the system is returned to energy-saving status.
When the time delay at the previous station td,n-1 is less than tcmax,n, a certain value of v1,n' should be determined so that the time delay can be exactly compensated at the next stop.In this case, we have: According to Equations ( 14) and ( 17), t2,n' can be determined as follows: 2, where p and q are the coefficients defined as follows: Substituting Equation (18) into Equation ( 8), we have: where: So we can determine v1,n' as: In conclusion, the whole algorithm is implemented according to the flowchart shown in Figure 5.

Case Study: Shanghai Metro Line One
To validate the correctness and evaluate the energy-saving effects of EOM, experiments were done on a traction box of a metro car in SMLO.Models were established in Simulink.All the data was acquired from Shanghai Shentong Metro Group Co., Ltd.

Experiment Data and Parameters
Data measurement circuits were set up to acquire the online data from the traction part of a moving train.The diagram of the circuit is shown in Figure 6.The HIOKI (Nagano, Japan) 8860-50 data logger, which has 16 channels of analog inputs, up to 1 MS/s sampling rate, was used to record the pantograph's DC supply voltage & current, breaking chopper current, and the train's speed.Some issues should be noted for the experiment: 1 The experiment was done on a separate power supply system designed especially for the test and maintenance, so the supply capacity is limited. 2 There were no other trains in the same supply network.3 During the experiment, the test metro train was unloaded.That means the mass of the whole test metro train is just its own mass, which is 296 tons.4 Most of the auxiliary electric devices, like air-conditioners, were switched off during the experiment.5 The metro train is composed of six cars.All the measurements were taken from one of the cars.
Figure 7 shows the experiment data of the test metro car. Figure 7a records the speed profile of the test train.It shows three periods of movement.From 3 to 23 s, the controller gave the full output to accelerate the train from 0 to 60 km/h.Then the traction motors were switched off for 7 s.The train started braking at 30 s. Figure 7b,c show the DC bus voltage and current of the test car, respectively.From Figure 7a it can be seen that the modeling of the driving strategy proposed in this paper is true and has already been applied on the real metro system.It could also be seen that when the train is accelerating, the voltage of the DC supply network was decreasing and the current demand of the train was growing fast because the train was in the process of changing from electric energy to kinetic energy.In addition, we can approximate the voltage and current as linear functions of time according to the diagrams.In the coasting period, the current demand decreased to zero because the power was switched off.The train speed dropped a little for the sake of extremely low resistance.In the braking process, the train started to feed the regenerative energy back.However, the regenerative braking was an effective method to save energy only when it is simultaneously used by other metro trains.If not, it would otherwise raise the voltage of the supply network, which is not permitted.Thus, a braking chopper was then used to limit the DC bus voltage by switching the braking energy to a resistor where the braking energy was converted to heat.The breaking chopper would take effect when the DC line voltage exceeds 1800 V in SMLO. Figure 7d shows the current of the braking chopper.From these diagrams, we can easily derive the power of different part of energy consumption by multiplying the corresponding voltages and currents.Note that the braking chopper power was not continuous because the braking resistor was switched on discretely.It is denser in the forepart and sparser in the rear.This is because that at the beginning of braking process, the velocity of the train is greater and thus more kinetic energy is turned into heat.If we integrate the braking chopper energy separately and add the feedback energy, we can get the absorbed power and the regenerative power of one metro car as shown in Figure 8. Recall that in Section 2, we analyzed and derived EOM based on the modeling of the power diagram with two triangles, which means linearity.Figure 8 validates the use of this power network model of metro systems.

Modeling of the Pilot Metro System
SMLO is one of the most important and busiest transportation systems of the city.It has 28 stations with a total distance of 37 km.For time-saving purpose, we only simulated the system with the most important segment, which contains six stations.The information of the pilot system is listed in Table 1.In the simulation, we number the six stations in sequence from 1 to 6.Each metro train departs at the 1st station and travels in sequence to the 6th station, which is called forward direction.Then, a turning of 60 s in duration occurs at the 6th station.After that, each train drives from the 6th station back to the 1st station, which is called backward direction.In the backward direction, the geographically nth station is defined as the (12 − n)th station in consistency with the forward direction.The segment between the nth and (n + 1)th station is defined as the nth segment.For simplicity, we modeled four trains in the system.Two of each depart from both ends of the track and turn around at the other end.The two trains have a headway time of 120 s.The schematic diagram is shown in Figure 9.The train line of the six stations is powered by two substations at both ends of the system.According to reference [48], we model each train as a current source.In addition, from Figure 8, the DC current should have the same shape (linear) as the power because the DC voltage is also linear.This can easily be shown according to the acceleration process of Figure 7b-c and Figure 8.The equivalent circuit of the system is shown in Figure 10a.(a) The equivalent circuit (b) The correspondence of the speed and current The equivalent current of each train can be determined by the train's speed and the acceleration rate according to Figure 7a,c.When the train is accelerating, the current is a positive value.When the train is coasting, it does not consume any energy and the current is zero.When the train is braking, which implies feeding the energy back, the current is negative.Note that the value of the negative current is determined according to Figure 8.In this case we assume that all the regenerative energy can be utilized by other trains.
The corresponding relationship between the current and the speed of the train is shown in Figure 10b.It should be noted that the data in Figure 7c shows the current value of only one metro car.In simulation, we should multiply it by six since a real metro train consists of six such cars.The voltage of substation one and the current of all the four trains in simulation are shown in Figure 11.

Specifying the Driving Strategy
In simulation, we got the accelerations of the accelerating and braking from Figure 7a.According to the figure, the absolute value of acceleration of coasting is very small compared with the other two.In fact, the specific Davis equation of SMLO is as follows: 2 r 16462.044.07 11.18 From Equation ( 26) it can be seen that different velocity values, from 0 to vmax, will change the running resistance by less than 30%.Consequently, it will have even less influence on the train's speed if we assume a suitable value of v in Equation ( 26), so taking a2,n as a constant of a2 for simplicity in both modeling and simulation is a reasonable assumption.Note that this assumption is acceptable only in particular situation that the gradient does not vary on the track segment considered.From all above, we get the acceleration values of the pilot system as shown in Table 2. Also, Equation (8) will have a specific form for the pilot system, which is: 1.0265 0.0176 0.9690 From Equation ( 27) and Table 1, the specified driving strategy for the pilot system is shown in Table 3.  2), we set the fixed dwell time to be 20 s and the flexible dwell time maximum 10 s in the pilot system by considering the real situation.For each train, there are eight dwell times that can be changed, that is at station 2, 3, 4, 5, 7, 8, 9, and 10.So there are 32 variables in total.The travelling time of the system, i.e., the simulation time, is set to be 1600 s, so the specified optimization problem of Equation ( 2) is as follows: . .0 10 1, 2,...,32 is integer 1, 2,..., 32 When MI-LXPM is used to solve Equation ( 28), we set the population size as 40, crossover fraction as 0.8, and elite count as 5.The algorithm finds the best solution in 50 generations.The results are shown in Figure 12.With the best solution, we can get the optimal energy-efficient timetable of departure time of the pilot system as shown in Table 4.

Applying CDSA
Under disturbance, CDSA is applied in order to recover the system to energy-saving status. Figure 13 shows the current and the corresponding energy consumption of different situations.Suppose that a disturbance of 20 s occurs on the 1st train at the dwell time of the 2nd station.The green line shows the optimal timetable of the 1st train.The red line and the blue line show the situations with and without CDSA, respectively.At the first two stations, the train is operating at its optimal status so that everything was fine and the three lines stuck together.From the 3rd station, the three lines separated.Both the red line and the blue line were delayed 20 seconds.However, with CDSA, the red line caught up with the green line and the system returned to energy-saving status from the 4th station.The ultimate energy consumption shows clearly that CDSA can compensate for a certain amount of energy.
To get an unbiased performance of CDSA, the Monte Carlo method [49] is used.In each experiment, a disturbance that has a normal distribution with expectation of 20 s is applied randomly to the four trains at a random station.Then CDSA is activated automatically on the delayed train to keep it on schedule.Finally, we compare the energy-saving effects of the two situations with and without CDSA.The result is shown in the next section.Note that CDSA is designed to be applied independently on each train on the line.It is focused on whether the specific train is delayed or not and then tries to compensate it, so if several trains are delayed, CDSA will take effect automatically and independently on each delayed train, which means there will be no significant difference if several disturbances happen.On the other hand, according to the current operating mechanism of SMLO, short delays are usually not being recovered or unable to be recovered until the end of a single journey as long as they are not significant enough to cause some safety problems.This existing mechanism is simulated and compared with the CDSA introduced in this paper in the next section.

The Energy-Saving Analysis of EOM
The energy-saving effect is remarkable with the optimal timetable.A total number of 100 experiments were done to get an unbiased statistical result.The result is shown in Table 5.From Table 5 it can be seen that with timetable optimization, 51.59 MJ (5.12%) of energy is saved in the pilot system.When a disturbance occurs, only 22.34 MJ (2.22%) of the total energy is saved on average.However, 40.66 MJ (4.04%) of the total energy can be saved if CDSA is applied.This shows that 81.98% more energy can be saved by using CDSA under disturbance conditions.Considering the small scale of the pilot system, the energy-saving effect is noticeable.If applied to a real metro system with much larger scale, i.e., hundreds of metro trains and stations, EOM will definitely have a much more significant result.
To get a more detailed analysis, the average energy-saving effect of each station where the disturbance happens is shown in Figure 14.It can be seen from the blue bar that the earlier the disturbance happens, the less energy would be saved if CDSA was not applied.This is because that TO synchronizes the accelerating and braking on a system level.The earlier the disturbance occurs, the earlier it will break the balance of the synchronization and the less energy will be saved consequently.On the other hand, from Figure 14, it can also be seen that more energy is saved with CDSA than without for earlier disturbances, which is shown by the difference of the red bar and the blue bar.For example, if the disturbance occurs at the 1st station, nearly 5% of the total energy can be saved with CDSA (red bar) while only 1% is saved without CDSA (blue bar).The difference is about 4%.If the disturbance occurs at the 8th station, nearly 5% of the total energy can be saved with CDSA.However, without CDSA, nearly 4% energy can still be saved, of which the difference is only about 1%.In other words, this shows that the CDSA has a better performance for earlier disturbances.
Note that the energy-saving percentage using CDSA is less when the disturbance occurs at the 9th station.This is because that using CDSA means raising the train's travelling speed temporarily, trading for sticking to the later system-level energy-efficient timetable.Faster speed means more energy consumption.This temporarily more energy consumption can usually be overwhelmed easily by the later system-level-saved energy, which is the case of stations 1 to 8.However, there is no more station after the 10th station.Thus, the using of CDSA caused an opposite effect at the 9th station.This phenomenon can also be seen from Figure 13.The application of CDSA caused more energy consumption temporarily at the beginning of the 3rd station.However, the eventual overall energy consumption is much less compared with the case without CDSA, so for a more realistic application of CDSA, whether it should be used or not at the end of a journey shall be considered seriously.

Conclusions
Metro systems consume huge amounts of energy every day.Saving even a small percentage of this energy consumption can have significant economic growth and environmental protection benefits.For this purpose, a new integrated Energy-efficient Operation Methodology (EOM) of the metro system is proposed in this paper.
EOM can be divided into two parts.First, to get a system-level energy-saving effect, the Timetable Optimization (TO) method is used.Compared with other methods, TO is a method without additional cost.Instead of storing regenerative braking energy, it rearranges the departure and arrival times of rolling stock so that the energy returned by the braking trains can be consumed simultaneously by those accelerating trains.A genetic algorithm is used to modify the dwell time of every metro train at each stop so as to make full use of the regenerative energy.Then, to solve the problem of disturbance, a novel algorithm named Compensational Driving Strategy Algorithm (CDSA) is proposed in order to save additional energy when disturbances happen.The main idea is to adjust the driving speed and coasting time of the disturbed train in the next several stations after a disturbance so as to fix the incoherence with the optimal timetable and recover the system to energy-saving status as soon as possible.
To validate the correctness and evaluate the energy-saving effect of EOM, a real case of Shanghai Metro Line One (SMLO) is studied.Experiments were done on a test metro train of SMLO.From the experimental data and corresponding data processing, the modeling of the driving strategy is validated.Moreover, the modeling of the traction and braking power of a single metro train is achieved.Based on this, the equivalent circuit model of the metro system is established.A six-station, four-train and bidirectional pilot system is built in Simulink according to the circuit model.Then EOM is specified to the pilot system.The results show that 5.12% of the total energy is saved in the ideal situation.When a disturbance occurs, only 2.22% of the total energy can be saved.However, 81.98% more energy on average can be saved with CDSA.The overall energy-saving effect in real situation is expected to be much more significant since there are many more trains and stations in real metro systems.
The proposed methodology was based on a simplified model of the real metro system.Future work will be focused on how to modify the whole methodology on a more detailed and realistic model.Adding energy-efficient driving into the model may also be concerned.Moreover, how to get the optimal timetable in a real-time level would probably be studied in the future as well.

Figure 1 .
Figure 1.Change dwell time to use regenerative energy.

Figure 2 .
Figure 2. The schematic diagram of the driving strategy. v

Figure 3 .
Figure 3. Drive faster to make up the delay.
Situation of disturbance (Delay eliminates the use of braking energy).

Figure 6 .
Figure 6.Schematic diagram of the electrical data measurement.

Figure 7 .
Figure 7. Experiment data of the test metro car.

Figure 8 .
Figure 8.The experiment power diagram of one metro car.

Figure 9 .
Figure 9.The schematic diagram of the pilot system.

Figure 10 .
Figure 10. Circuit model of the metro system.

Figure 11 .
Figure 11.The voltage and current in simulation.

Figure 12 .
Figure 12.The results of the genetic algorithm.

Figure 13 .
Figure 13.Current curves and the corresponding energy consumption curve.

Figure 14 .
Figure 14.Average energy-saving effect with and without CDSA among different stations. ]

Table 1 .
The information of the pilot system.

Table 2 .
The acceleration values of the pilot system.

Table 3 .
The specified driving strategy of the pilot system.

Table 4 .
The optimal departure time at each station.

Table 5 .
Average energy-saving effect of different situations.

situations Average energy consumption (MJ) Average energy saving (MJ) Average energy-saving percentage (%)
Etotal energy-consumption of the system Q number of substations of power supply for the metro system Uj DC voltage of the jth substation Ij DC current of the jth substation Fr running resistance Ln distance between the nth and (n + 1)th station v1,n instant speed at the end of accelerating between the nth and (n + 1)th station v2,n instant speed at the end of coasting between the nth and (n + 1)th station vmax maximal speed of the train vlimit speed limitation a1,n acceleration of accelerating between the nth and (n + 1)th station a2,n acceleration of coasting between the nth and (n + 1)th station a3,n acceleration of braking between the nth and (n + 1)th station tstart,n starting time of leaving the nth station t1,n duration of accelerating between the nth and (n + 1)th station t2,n duration of coasting between the nth and (n + 1)th station t3,n duration of braking between the nth and (n + 1)th station tn total running time between the nth and (n + 1)th station t4,n duration of dwelling between the nth and (n + 1)th station td,n time delay caused by disturbance between the nth and(n + 1)th station tcmax,n maximal compensable time delay between the nth and (n + 1)th station tf total travelling time of the system