#### Abstract

In the moving block signalling (MBS) system where the tracking target point of the following train is moving forward with its leading train, overload of the substations occurs when a dense queue of trains starts (or restarts) in very close distance interval. This is the peak power demand problem. Several methods have been attempted in the literature to deal with this problem through changing train’s operation strategies. However, most existing approaches reduce the service quality. In this paper, two novel approaches—“Service Headway Braking” (SHB) and “Extending Stopping Distance Interval” (ESDI)—are proposed according to available and unavailable extra station dwell times, respectively. In these two methods, the restarting times of the trains are staggered and traction periods are reduced, which lead to the reduction of peak power demand and energy consumption. Energy efficient control switching points are seen as the decision parameters. Nonlinear programming method is used to model the process. Simulation results indicate that, compared with ARL, peak power demands are reduced by 40% and 20% by applying SHB and ESDI without any arrival time delay, respectively. At the same time, energy consumptions are also reduced by 77% and 50% by applying SHB and ESDI, respectively.

#### 1. Introduction

Moving block signalling (MBS) [1] was proposed a few decades ago to reduce headway among successive trains in a track line. Theoretically, two successive trains are separated by a short distance, which is equivalent to the braking distance of the following train from its current speed, as well as a safety margin. This distance can be changed with the limit for the given operating speed and train characteristics, such as train length and braking rate. In MBS, when a leading train stops for a long time, the following trains will stop at the tail of the leading train. When this leading train restarts, the following trains will start almost simultaneously. It could cause synchronization of the peak demand of trains and increases the total peak power demand significantly, which is highly energy consumed and may lead to overload of the nearby substations. How to reduce peak power demand is called “Peak Demand Reduction” (PDR) problem.

This problem could be solved by improving the infrastructure or the train operation strategies. Energy storage system (ESS) is a very important component in modern railway power supply system. Advanced control and manufacturing technology improve the stability, capacity, and weight of ESS [2–5], which makes it suitable to be equipped in the substations or in trains. It could store the regenerative braking energy and assist the main power source during train’s acceleration period. However, it is costly to equip ESS since the whole system needs a large number of them, for example, in urban railway system, nearly each of the stations or trains needs an ESS. At the same time, the cost of maintenance for ESS is also huge, especially when the frequency of acceleration*∖*braking is high (which means the high frequency of charge and discharge and leads to life time reduction). Furthermore, it is inconvenient to be adopted by existing lines since the reconstruction project is costly and complex. Therefore, changing train operation strategies is a more convenient and economical way to reduce peak power demand.

There are two kinds of traditional PDR techniques based on changing operation strategies; one is called starting time delay (STD), which introduces a starting time delay to each of the following trains. Under this category, there are two specific techniques called single STD and grade STD. The difference between them is the introduced starting time delay to each of the following trains, which are the same in single STD technique but different in grade STD (showing a deceasing trend). The other one is called acceleration rate limit (ARL), which means the acceleration of the following trains is limited to a certain extent (or different extents). Under this category, there are two specific techniques called single ARL and grade ARL. The difference between them is the limited acceleration rate to each of the following trains, which are the same in ARL technique but different in grade ARL (showing a deceasing trend). In addition to the above techniques, there is also a PDR technique called coordinated PDR. It is a combination of the STD and ARL techniques with feeding the regenerated power of decelerating trains to accelerating trains in the same queue, by coordinating the movement of queued trains.

Takeuchi and his colleagues discussed these techniques in [6–8]. Simulation results show that the graded ARL technique has the best performance in reducing peak power demand among single STD, grade STD, single ARL, grade ARL, and coordinated PDR techniques. In the traditional techniques, time delay is introduced and quality of service is degraded. Ho and Wong [9] use an expert system to help the operators for decision making, and it focuses on the balance between time delay and peak power demand. In [10], by braking and powering the trains simultaneously, Albrecht proposed a way to reduce the peak energy consumption and maximize the regenerative energy. Chen and his colleagues [11] proposed a method to minimize the peak energy consumption by adjusting the train dwell times at each station in MRT systems. Kim et al. [12] developed a mixed integer programming model to minimize the peak power energy demand that occurs when trains are running simultaneously, and the basic principle of the method is the same with STD. All of the methods proposed in [6–12] are implemented after the trains’ restarting, by adjusting the timetable or driving strategy.

Although the existing techniques can reduce peak power demand at different degrees, they increase the travel time between the successive stations and decrease the service quality. And the energy consumption is increased since there are more traction periods.

In this paper, in order to reduce peak power demand and energy consumption, we first analyze the reasons of the formation of the peak power demand, and, then, two novel approaches are proposed based on the main reasons with considering the energy efficient driving strategies. One is for available extra station dwell time, named Service Headway Braking (SHB), and the other one is for unavailable extra station dwell time, named Extending Stopping Distance Interval (ESDI). Both of them are real-time adjustment methods and can be implemented before the leading train’s restarting. Therefore, there is no need to change timetable. Considering energy saving driving strategy could be seen as a kind of driving mode switching process [13–16]; a nonlinear programming model is contributed and the simulation results show that, compared with the best traditional PDR techniques, both of these two methods could reduce the peak power demand and energy consumption without any arrival time delay increasing. And when the available extra station dwell time is the same as the unavailable extra station dwell time, SHB performs better than ESDI since more peak power and energy consumption could be reduced.

#### 2. Tracking Dynamics and Peak Power Demand in MBS

##### 2.1. Tracking Model in MBS

Under MBS, the tracking target point of the following train moves forward continuously as the leading train travels. The instantaneous distance of two successive trains is calculated as where is the position of the leading train’s head and is the position of the following train’s head.

The distance between two successive trains must be larger than the safety margin at any moment even if the leading train comes to a sudden halt, so we have where is the length of the train, is the length of safety margin, is the instantaneous speed of the following train, and is deceleration rate.

Based on (1) and (2), the relation between the leading train and following train should satisfy which implies that the instantaneous speed and position of the following train should satisfy

##### 2.2. Reasons of the Formation of Peak Power Demand

The reason of the formation of the peak power demand is the restarting of the dense queue, and the reasons for the formation of the dense queue are listed as follows.(1)Features of moving block signaling system. (Two trains will start simultaneously if the distance interval between them is (2)Extra dwell time in station.

There are two kinds of extra station dwell time: one is available and the other is unavailable. In daily railway operation, there may be some exceptions, such as a passenger may be caught in the door of the train or a short-term surge in passenger flow (i.e., passenger flow increases sharply after a football match). In these circumstances, adjusting the whole timetable is not convenient, because the circumstances only exit in a short period. In this case, the operator will arrange the train to stop a little longer and this extra station dwell time is available. However, if a train is broken in a station and we are not sure how long we need to fix it, in this situation, the extra station dwell time is unavailable.

Based on the analysis above, it is known that peak power demand could be reduced by avoiding the dense queue. In order to achieve this goal, we first analyze the relation between extra station dwell time and the number of delayed trains.

##### 2.3. Station Delay Propagation

Generally speaking, each train has a required dwell time at a station. If a train stops longer than the required dwell time, we call the extra time as *delay time*. In this section, we focus on the relation between the *delay time* and the number of delayed trains. In MBS, the *delay time* may impact the following trains and cause a dense queue. Figure 1 shows the formation of the dense queue.

As it is shown in Figure 1, there are trains in the track. Train 1 stops at station A. The position of station A is . For each train, the dwell time is . The target speeds of the trains are the same and constant. According to the normal condition, each train arrives at station A and stops for and then starts to run.

When train 1 starts, the positions of following trains are , . The tracking time (the running time of the th train running from to ) between two successive trains is ; it is also the scheduled departure time interval. However, the following trains may become a dense queue if train A does not run immediately after . Defining is the *delay time* of the leading train after and is the number of the delayed trains caused by . Based on (5), train which is delayed should stop at point ; in other words, is the stop position of train in the dense queue:

Let be the running time in which train arrived at . Based on (2), we have , , , and .

Letting , we have

then the number of delayed trains should satisfy

##### 2.4. Peak Power Demand Calculation

The power demand of the th train is calculated by where is the traction force of the th train and is the speed of the th train after restarting.

The peak power demand of the th train is calculated by where is the traction force of the th train and is the highest speed of the th train after restarting.

The total peak power demand of all delayed following trains is calculated by where is the number of delayed following trains.

##### 2.5. Energy Efficient Control Switching Points

Because there are more traction periods, traditional PRD techniques cause energy consumption increasing. Energy conservation is the research interest in many fields [17, 18], especially in rail transport [13–16]. Therefore, it would be better to reduce peak power demand and energy consumption at the same time. Train operation is a switching process; the operation modes are switching among traction, braking, and coasting (with no traction and braking force). There are lots of researches on switching system [19–21]. For optimal control of switching system, we could apply the methodologies proposed in [22, 23]. However, the computation complexity will be a big problem when we adopt the technologies above. In this paper, we first apply Pontryagin maximum principle to find the discrete optimal control modes and then treat the control switching speed as the decision parameters to obtain a reference trajectory.

For electric traction systems, the motion equations of train have the following forms: where is the position of the train, is the speed of the train at position , is the time at which the train is located at the position , is the relative traction force, is the specific maximum traction force per mass unit, is the relative braking force, is the specific maximum braking force per mass unit, and is the specific basic resistance.

Note that the basic resistance is usually given by Davis equation: where , , and are regression coefficients obtained by fitting test data to the Davis equation. accounts for air resistance. account for mass and mechanical resistance. And has following characteristics:

The energy efficient operation problem is modeled as follows: where is the specific mechanical work of the traction force, is the speed limit. The boundary conditions are

By applying Pontryagin maximum principle to solve the problem as specified in (15) with constraints (16), the optimal solution should maximize the Hamiltonian function: where should satisfy the differential equation:

It is easy to prove that the Hamiltonian reaches the maximum with respect toand. And there are five energy efficient control modes as follows.(i) Full power (, ) if .(ii) Partial power (, ) if .(iii) No power and no braking (, ) if .(iv) Partial braking (, ) if .(v) Full braking (, ) if .

They could be seen as four possible driving phases, which are (i) acceleration with full power;(ii) speed holding with partial power or braking;(iii) coasting with no power and braking;(iv) braking with full braking.

Based on the above analysis, it is seen that the energy saving strategy relies on these optimal controls. Train energy saving driving process is the process of switching from an optimal control to another, and the switching sequence is acceleration, speed holding, coasting, and braking.

Then, the optimal control switching speeds and running phases duration are treated as decision parameters. The energy and power can be expressed by these decision parameters. The power demand is used as the objective function of a nonlinear programming problem. The constraints include running time and distance. In actual train operation, maximum power and braking cannot be applied by considering the ride comfort. Instead, a service acceleration/braking rate is applied. Therefore, by solving the nonlinear programming model, a reference trajectory leading to less power and energy consumption can be obtained.

#### 3. Peak Demand Reduction Techniques

##### 3.1. Service Headway Braking Strategy

Based on the analysis above, we know the restarting of the dense queue in a small area leads to peak power demand and both of the two traditional PDR techniques are carried out after the formation of the dense queue. In this section, we propose a novel operation strategy to reduce the peak power by avoiding the formation of a dense queue based on the available extra station dwell time. In the following parts, we use *reaccelerate* to indicate train accelerates from a nonzero initial speed and restart to indicate train starts from 0 m/s.

###### 3.1.1. SHB Strategy Analysis

Figure 3 shows the new operation strategy; the leading train stops at station A and the position is . The black thick solid line is calculated by (2), and we call it braking curve in this paper. Under normal circumstance, the following train should brake when it touches the brake curve. However, if the leading train starts to run when the following train touches braking curve, the following train will not brake but reaccelerate. Therefore, in the new strategy, if the following train touches the braking curve when the leading train starts to run, and the speed at this moment is not 0 m/s, the moment in which each train reaches the peak power demand can be staggered. Meanwhile, since the following trains brake at first, they will not stop too to cause a dense queue. By these two reasons, the total peak power demand can be reduced. Based on the analysis above and considered the energy saving operation process, a new operation strategy is proposed as follows.

As shown in Figure 2(a), the target speed of the trains in the track is ; the leading train has an extra dwell time, which is . If is long enough to cause a dense queue with trains, then, after , let train brake with service braking deceleration to (the position at this time is ); if , then it stops at for , then accelerates to with service acceleration (the position at this time is ), and then runs with , () for to ; after that it reduces traction force and keeps the train moving with a constant deceleration to (the position at this time is and are a pair of points in the braking curve). At this time, train -*1* restarts, therefore, train starts to track train *-1 *according to the moving block tracking distance interval until arrives at station A.

**(a)**

**(b)**

Figure 2(b) shows the position of train 3 when train 2 starts to move from station A. It is observed that train 3 will reaccelerate from the braking curve () to arrive at station A (); we use to represent this time period and all of the delayed following trains need this time period to arrive at station A after their leading train’s restarting. In this circumstance, for train 2, the running time between and is ; for train 3, the running time between and is ; for train , the running time between and is .

In order to obtain , a simulation is done as follows: (m), (m). The leading train stops at station A where the position is 400 (m): (m). Let the following train start from the point in braking curve; the speed and position of the following train is calculated by (4) or (5); the results are shown in Table 1.

From Table 1, we can see that the arrival times of following trains seem constant when they reaccelerate from the braking curve. In order to show the trend of this tracing process, Figure 3 gives the *V*-*T* profile when the reaccelerating speed and position are 8 m/s and 178 (m).

Based on the analysis above, the whole strategy is consisting of 6 steps, braking, waiting, traction, speed holding, coasting, and finally tracking the front train. The indexes of them are .

It is worth noting that this operation strategy is energy efficient, because the duration of some steps may be 0 based on an appropriate target function. For example, if we get based on a certain target function, there will be no traction process before reaccelerating and the train will brake and coast until reaccelerating.

In this new strategy, all the following trains brake at first; therefore, they will not stop too to cause a dense queue; at the same time, the other trains are coasting when one train reaccelerates, so the reaccelerating times of the following trains are staggered. Therefore the peak power demand is avoided.

###### 3.1.2. Model of SHB

In this section, we formulate the mathematical model of the operation process. For train , define as the starting time of each step; and are the running distance and time of each step, respectively. In order to stagger the reaccelerating times of the following trains, we should take an appropriate value to . Because if , the successive two trains will reaccelerate simultaneously again and peak power demand cannot be reduced. On the contrary, if is close to , the effect of energy saving will be reduced. Therefore, we recommend . At the same time, for train , the total running time and distance should satisfy the following equations:

Since the traction phase is highly energy consumed and power supported, we minimize the peak power. Therefore, based on the analysis above, the problem could be seen as a constrained nonlinear programming problem as follows:
where
where are penalty factors, ; are service acceleration and deceleration of the train; is a very small deceleration; is the index of the delayed following trains ; is the index of operation steps, ; is the starting speed of each step for train ; is the position of the leading train (the leading train); is the position of the th following train; is the running distance of the th step for train ; is the duration of of the th step for train ; is the braking distance if train brakes from ; is the scheduled station dwell time; is the tracking time in which train tracks train *-1* from to .

During the running process of successive following trains, the time interval among them is . When train starts to move from , the distance between train and train *-1* is shortest; at this time, the speed of train is larger than and the distance interval between them is at least ; according to the general condition in mass transit system, , , , , , it is easy to know that ; this is satisfied by (3). Therefore, it is no need to consider the distance constraint shown in (3) between successive following trains in (22) again.

##### 3.2. Extending Stopping Distance Interval Strategy

The previous section shows SHB strategy; it could reduce the peak power demand if the *delay time* is available. In this section, we will show a strategy which could reduce the peak power demand by decreasing the density of the dense queue. By adopting this strategy, the peak power demand could be reduced if the extra station dwell time is unavailable. Because, in the same , the longer the stopping distance interval of the following trains, the lower densities of the stopped train queue. Therefore, when they restart, the peak power demands are reduced.

###### 3.2.1. ESDI Strategy Analysis

According to the analysis in Section 3.1.1, it can be seen that no matter the train reaccelerates from which point in the braking curve, is nearly the same. Based on , we can extend the stopping distance interval of the following trains to reduce the density of the dense queue. The new operation strategy is proposed as follows.

Figure 4 shows the new operation strategy. The operation process of the following trains when the leading train dose does not move after could be divided into 6 steps, speed holding, braking, traction, speed holding, coasting, and braking to stop. Define as the starting time of each step,, and , ; define the running time and distance of each step as and , respectively. When the leading train has an extra dwell time, each of following trains will run according to these 6 steps. Define as the traveling distance during ( is also stopping distance interval) for the following trains. During this process, the running time of steps 3 to 6 should be equal to .

In this new strategy, the stopping distance interval of the following trains is increased, so the density of the queue is decreased, and the peak power demand is reduced.

###### 3.2.2. Model of ESDI

In this section, we formulate the mathematic model of the operation process. According to the last section, it is known that we know that the stopping intervals of the trains are extended in the new operation strategy; therefore, peak power demand can be reduced when the trains restart. Define *MI* as the stopping interval of two successive trains; we have

In order to grantee the feasibility of the new operation strategy, if leading train stops at the station, after , the distance between the first following and the leading train should satisfy the following equation: where is the second train’s position when the leading train stops at station, and it is calculated by

Since the traction phase is highly energy consumed and power supported, we minimize the peak power. Based on the analysis above, we use nonlinear programming to model the problem as follows:

where , are penalty factors, ; , are service acceleration and deceleration of the train; is the starting speed of step ; is the position of the leading train (the leading train); is the position of the following train; is the second train’s position when the leading train stops at station; is the scheduled station dwell time.

#### 4. Simulation and Discussion

In this section, a simulation is used to test and verify the new strategies. The length of train () is 140 m, safety margin () is 50 m, service tracking headway is 120 seconds, dwell time () is 10 seconds, target speed () is 16 m/s, service acceleration rate () is 1 m/s^{2}, service braking deceleration rate () is 1 m/s^{2}, and position of station A is 3710 m ().

In order to choose , we analyze the practical data of coasting phase from Dalian Fast Track [24]. Because the speed in coasting phase declines very slowly, so the profile of speed time could be seen as a straight line and the slope of the line is the deceleration rate of coasting phase. We use least-square procedure to fit the speed-time date sectional and the results are shown in Table 2.

From Table 2, it is clear that the higher the coasting starting speed, the greater the deceleration. In order to supply a small traction force to keep the train moving in a constant deceleration , we choose m/s^{2} as appropriate.

##### 4.1. Service Headway Braking Strategy

If is 250 seconds, according to (8), 3 trains will be delayed (including the leading train). Because we choose , therefore, is selected from Table 1. Based on (21)-(22), we have

Figures 9 and 10 and Table 3 show the simulation results.

In order to show the peak power demand without any PDR technique and compare it with the performance of traditional PDR technique, Figures 5–8 are given. According to Figures 5 and 6, we can see that the two following trains are starting simultaneously from 250 s; the peak power demand is 25.1 kw/t. The arrival times of train 2 and train 3 are 289.38 s and 338.76 s, respectively.

Figures 7 and 8 show the performance of graded ARL technique. The acceleration of train 2 and train 3 is 0.5 m/s^{2} and 0.3 m/s^{2}. By applying different acceleration, the peak power demand is reduced to 22.07 kw/t. However, the time delay is increased. The arrival times of train 2 and train 3 are 297.78 s and 400.36 s, respectively. That means the arrival times of the two trains are 8.4 s and 61.6 s later than their times without PDR techniques, respectively.

The performance of applying SHB technique is shown from Figures 9 and 10. As can be seen, train 2 has a waiting time 65.66 s at the position of 2048 m and then restarting at 81.66 s. At this time, train 3 is coasting; therefore, there is no peak power in this time. The two trains restart with different speeds in different time points, which can be seen form the red and blue circles in Figure 9. That means the reaccelerating times of the two trains are staggered. Furthermore, the two trains reaccelerate with nonzero initial speeds, which means the reacceleration periods are reduced. Therefore, the peak power demand is reduced to 13.42 kw/t. According to Figure 9, the two trains have no speed holding phases, and train 3 has no traction phases before step 6. Therefore, the energy consumption is reduced. The arrival times of train 2 and train 3 are 289.5 s and 339.16 s, respectively.

Table 3 shows the result data of the simulation. Graded ARL technique can reduce the peak power demand; however, the arrival times of the two trains are delayed significantly. SHB technique has great advantages. According to Table 3, the delay times of the two trains are very short; at the same time, the energy consumption is also reduced to a low level, even less than value without any PDR technique. In addition, the stopping times before arrival of the two trains are reduced to 65.66 and 0, respectively.

##### 4.2. Extending Stopping Distance Interval Strategy

According to Section 3.1.2, we choose s; when the leading train stops at the station and dose not start to run after , based on (27), we have

If the extra dwell time is 250 seconds, in graded ARL, 3 trains will be delayed (including the leading train) according to (8). And the length of the dense queue is 380 m. The performance of graded ARL technique can be seen in Figure 5.

In ESDI, based on (24) the stopping distance interval for the following trains () is 383.68 m. Based on (8), 2 trains (besides the leading train) will be delayed. However, the length of the dense queue is 767.36 m. In the same time period and distance area with graded ARL technique, we calculate the peak power demand, and the performance is shown in Figures 11 and 12.

As we see, using ESDI technique, the arrival times of the following two trains are almost the same as their times with no PDR technique. Without any PDR technique, the stopping distance interval is 190 m, so the dense queue caused by the extra station time (250 s) is 380 m. The peak power demand after 250 s within 380 m is 25.1 kw/t. In ESDI strategy, stopping distance interval is extended to 383.86 m, which is two times longer than the distance in ARL. That is to say, the restarting position of train 2 belongs to the nearby substation. Thus the power demands by train 2 and train 3 are not afforded by the same substation. Therefore, the peak power demand after 250 s within 380 m is reduced to 18.18 kw/t. Compared with non-PDR technique, the peak power demand is reduced by 18%. In ESDI, train 2 has a stopping time 153.98 s at the position of 2944.6 m. Train 3 has a stopping time 57.96 s at the position of 3328.32 m. Compared with non-PDR technique, the stopping time before arrival is increased by 18%.

According to above analysis, we can see that the new strategy reduces peak power demand by sacrificing the stopping time before arrival. However, the advantages of ESDI are still obvious, because it is more efficient on energy saving; it can reduce energy consumption by 50% compared with graded ARL technique.

From Table 4, it is seen that, although the arrival times of trains when applying SHB and the arrival times of trains when applying ESDI are the same, SHB reduces more peak power demand and energy consumption, which means SHB performs than ESDI when the available and unavailable extra station dwell times are the same.

#### 5. Conclusion

Peak power demand reduction strategies are discussed in this paper. The reasons of peak power demand problem are analyzed deeply and two main reasons are given.

Based on the reasons and according to different situations, two new peak demand reduction techniques are proposed. One is Service Headway Braking (SHB) strategy, which is used to reduce peak power demand when the extra station dwell time is available. The other is Extending Stopping Distance Interval (ESDI) strategy, which is used to reduce peak power demand when the extra station dwell time is unavailable. Nonlinear programming approach is introduced to model the operation strategy. The simulation results show that, compared with the best traditional PDR techniques, SHB can reduce 40% of peak power demand and 77% of energy consumption without increasing the arrival time delay. ESDI can reduce 20% of peak power demand and 50% of energy consumption without increasing the arrival time delay. Therefore, SHB has a better performance than ESDI when the available and unavailable extra station dwell times are the same. The basic principle of the two strategies is to avoid the formation of the dense queue. Therefore, both of the two strategies are implemented before the formulation of the dense queue.

#### Conflict of Interests

The authors declare that there is no conflict of interests.

#### Acknowledgments

This work was supported by the National High Technology Research and Development Program of China (no. 2011AA110502) and the National Natural Science Foundation of China (no. 61134001).