A Fuzzy Optimization Model for High-Speed Railway Timetable Rescheduling

. A fuzzy optimization model based on improved symmetric tolerance approach is introduced, which allows for rescheduling high-speed railway timetable under unexpected interferences. The model nests di ﬀ erent parameters of the soft constraints with uncertainty margin to describe their importance to the optimization purpose and treats the objective in the same manner. Thus a new optimal instrument is expected to achieve a new timetable subject to little slack of constraints. The section between Nanjing and Shanghai, which is the busiest, of Beijing-Shanghai high-speed rail line in China is used as the simulated measurement. The fuzzy optimization model provides an accurate approximation on train running time and headway time, and hence the results suggest that the number of seriously impacted trains and total delay time can be reduced signiﬁcantly subject to little cost and risk.


Introduction
The infrastructures of high-speed railway have been extensively developed in China for the past several years.The network topology structure and operation mode of the railway are changing profoundly.The target is to cover its major economic areas with a highspeed railway network, which consists of four horizontal and four vertical lines 1 , in the following several years.The network scale is much larger than any existing ones in the world.However, passengers have to face huge challenging variables generated by unexpectedness, for example, weather, equipment failure.On June 20, 2010, for example, the southern part of the existing Beijing-Shanghai line experienced speed restriction for a heavy rain storm.As announced, 18 trains departing from Shanghai were canceled, thereby more than 20,000 passengers being affected.Also, the rest in-service trains were more or less delayed.Considering the high speed and frequency of bullet train, the impact of train delay and with fuzzy passenger demand on a single line railway.Since the number of passengers getting on/off the train at each station is assumed to be a fuzzy variable, the total passengers time is also a fuzzy variable.An expected value goal-programming model is constructed to minimize the total passenger time and the total delay.A branch-and-bound algorithm based on the fuzzy simulation is designed in order to obtain an optimal solution.Then they extend the uncertainty in 31 by considering stochastic and fuzzy parameters synchronously to solve the railway freight transportation planning problem.Based on the chance measure and critical values of the random fuzzy variable, three chance-constrained programming models are constructed for the problem with respect to different criteria.A fuzzy periodic job shop scheduling model is introduced to address the framework of the periodic robust train scheduling problem in 32 .Fuzzy approach is used to reach a tradeoff among the total train delays, the robustness of schedules, and the time interval between departures of trains from the same origins in this paper.Cucala et al. 33 proposes a fuzzy linear programming model to minimize energy consumption with uncertain delays and drivers behavioral response.The method is applied to a real Spanish high speed line to optimize the operation, and comparing to the current commercial service evaluates the potential energy savings.Wang et al. 34,35 and Guo et al. 36 analyzed the driver's safety approaching behaviour and pedestrian safety crossing behaviour in the urban traffic environment.The traffic participants revealed different behavioral decisions with various personal characteristics and can be described as fuzzy parameters.Overall, none of the previous papers analyzed the uncertainty of all factors in train rescheduling, although they considered some constraints as fuzzy member, such as passenger time and delay time.Furthermore, the fuzzy constraints are defined as triangular or trapezoidal fuzzy numbers, which are very difficult to determine the distribution.
In this paper, we attempt to achieve optimal timetable rescheduling under the uncertainties, for example, constraints and/or unexpected parameters, by means of proposing the fuzzy optimization model as discussed in the above.In the rest of this paper, typical rescheduling model will be discussed in the section of Timetable Rescheduling Problem.In the following section, we will describe in details of the fuzzy optimization model based on improved tolerance approach to timetable rescheduling, including the fuzzy membership functions of the original objective and soft constraints.A case study on the busiest section of Beijing-Shanghai high speed line will be illustrated in the section of Case Study.The final section concludes the results of the paper and suggests for further research.

Timetable Rescheduling Problem
The aim of train rescheduling is to get a new timetable that adjusts the train movements to be consistent with the planned schedule as much as possible under some interference 37 .The following model focuses on minimizing the total delay as well as the number of seriously impacted trains.

Input Data
Take a rail line with n trains and m stations for example.The numerical inputs are described as follows: S j : the station j according to original timetable.
x * i,j : the departure time of train i at station j according to original timetable.
y * i,j : the arrival time of train i at station j according to original timetable.a i,j : the minimum running time of train i on section j, j 1 .τ f j : the time interval between two adjacent trains at station j. τ d j : the headway of section j, j 1 .
τ tr j : the track number of station j.T i,j : the minimum dwell time of the train i at station j.
h i,j : if train i stop at station j according to original timetable, then h i,j 1; otherwise, h i,j 0.
: the cost per time unit delay for train i.
T D i : delay tolerance for train i. θ: weight for objective.
M: a very big integer, for example, 100000.

Decision Variables
The decision variables are described as follows: x i,j : the new departure time of train i at station j after adjustment.y i,j : the new arrival time of train i at station j after adjustment.d i,j : delay of train i at station j, which is defined as the difference between the arrival time after adjustment and the planned arrival time in the original timetable.b i : if train reaches its final considered stop with a delay larger than w i , then b i 1; otherwise, b i 0. η i,j,l : if train i use track l at station j, then η i,j,l 1; otherwise, η i,j,l 0. α i,j,i ,j : if the departure time of train i at station j is earlier than the departure time of train i at station j to original timetable, then α i,j,i ,j 1; otherwise, α i,j,i ,j 0. β i,j,i ,j : if the departure time of train i at station j is changed to occur after the departure time of train i at station j , then β i,j,i ,j 1; otherwise, β i,j,i ,j 0.

Objective Functions
i To minimize the delay cost: 2.1 ii To minimize the number of seriously impacted trains: We set the final objective function as:

Constraints
i Section running time restrictions:

2.4
The real departure time cannot be earlier than the original departure time:

2.5
ii Station dwell time restrictions: x i,j − y i,j T i,j .

2.6
iii Track restrictions: iv Station headway restrictions.
For each station, if two trains use the same track, at least one of α and β is forced to be 1:

2.10
v Section headway restrictions.
For each section, at least one of α and β is forced to be 1 because there is only one track: α i,j,i ,j β i,j,i ,j 1, 2.11 x i ,j − x i,j τ d j α i,j,i ,j − M 1 − α i,j,i ,j , 2.12 x i,j − x i ,j τ d j β i,j,i ,j − M 1 − β i,j,i ,j .

2.13
vi Auxiliary restrictions: In practice, some constraints of the model are not strictly satisfied due to the inexact operation time.Thus four of the above constraints need to be changed as below.
(i) Section Running Time Restrictions a i,j is the minimum running time of train i in section j, j 1 , including the pure running time in section and the additional time for train stop or departure at the station.The minimum running time is decided by the length of the section, the infrastructure characters of the section, and the train type.In the actual operation environment, the train running speed is not a constant value because of the various infrastructure characters of the railway, that is, bridges, tunnels, and culverts.Furthermore, some factors i.e., railway equipment statuses, technological level of crew, and weather condition also increase the uncertainty of running time.Thus it is very important to find a safe and reasonable average speed that is far below the limited speed and full utilization of railway capacity, which will greatly improve the optimization result especially when some trains are delayed due to some interference.Since the minimum running time of the train in the section usually changes within a certain range, there should be a tolerance for a i,j .Thus 2.4 does not need to be strictly satisfied and can be changed as below:

ii) Station Dwell Time Restrictions
T i,j is the minimum dwell time of the train i at station j, including the pure operation time at station, passengers on and off time, crew setup time, and some additional time, like waiting for other trains.The minimum dwell time is decided by station operation type, station level, and the train type.In the actual operation environment, railway equipment statuses, driving behavior of drivers, and environment condition may increase the uncertainty of running time.
Thus it is also important to find a safe and reasonable average dwell time that can be full utilization of station capacity.Similar to section running time, there also should be a tolerance for T i,j , and the dwell time constraint can be changed as below: 2.16

(iii) Station Headway Restrictions
The headway time for each station is decided by the number of receiving-departure track, the operation time of turnout, and holding time of the track, which also faces uncertainty problems due to the factors of equipment status and human technological level.So there should be a tolerance for τ f j .Equations 2.9 and 2.10 do not need to be strictly satisfied and can be changed as below: 2.17

(iv) Section Headway Restrictions
The headway time for each section is decided by the number and length of the block between two tracking trains and the speeds of the trains, which ranges from 2.4 minutes to 3 minutes according to Shi 38 .Similar to station headway constraint, 2.12 and 2.13 can be changed as below: g 3 x i ,j − x i,j τ d j α i,j,i ,j − M 1 − α i,j,i ,j , g 4 x i,j − x i ,j τ d j β i,j,i ,j − M 1 − β i,j,i ,j .

2.18
Since the four constraints are changed as above, the model turns to be not representative in the sense of mathematical viewpoints.To construct a reasonable mathematical model under the uncertain environment, the tolerance approach based timetable rescheduling model will be introduced in the next section.

Fuzzy Optimization Model for Timetable Rescheduling
In the paper, we use the fuzzy optimization based improved tolerance approach to solve the uncertainty program, and some necessary backgrounds and notions of the approach are reviewed.

Improved Tolerance Approach
Tolerances are indicated in any technical process, that is, the admissible limit of variation around the object value and the deviations allowed from the specified parameters 39 .
A general model of a fuzzy linear programming problem is presented by the following system 40 .A i,j , B i , and The symbol ⊕ represents the extended addition.Each real number can be modeled as a fuzzy number. 3.1 Here we only discuss the special case 3.2 where the objective function is crisp, some constraints have the soft form, and the rest constraints are crisp.

3.2
The m 1 soft constraints may be described more precisely by the fuzzy set with the support b i −d i , b i , and d i is the tolerance according to b i .Moreover the membership function of g i x can be specified as 3.3 , and its graphics is Figure 1 a .Then we can directly assigns a measure of the satisfaction of the ith constraint to the solution X x 1 , x 2 , . . ., x n .

3.3
Fuzzy constraints will inevitably lead to the fuzzy objective based on the ideology of symmetric model.Thus, for the objective function z x , there is a fuzzy set and z Min x∈X L z x .So the membership function μ z z of Z is given by 3.4 , and its graphics is Figure 1 b .

3.4
In order to determine a compromise solution, it is usually assumed that the total satisfaction of a decision maker may be described by λ x min μ z x , μ 1 x , μ m 1 x .Since not all the constraints are equally important to the solution in practice, we take λ i to describe the priority of the ith soft constraint to the solution, and λ x can be described by λ x m 1 1 w i λ i x .w i is the weight of λ i , which shows the importance of different soft  constraints to the optimization purpose.Generally, there are many complex constraints for different optimization systems in practice.It is very hard to handle all constraints during the system optimization, and therefore the systems always treat the constraints with different priorities.Moreover the decision makers also hold different views on the importance of various factors about the same issue.The timetable rescheduling problem, for example, will take different optimization objectives and constraint priorities in different emergencies.The dispatcher will focus on the reducing of train delay time and train headway time if few trains are disturbed by equipment failure.However, minimizing the number of seriously impacted trains and the train interval running time as much as possible are more effective measures to return to the normal states, when some accidents or natural disasters lead to speed restriction on lots of sections for a long time.
Since we treat the objective in the same manner as the soft constraints, this is also a symmetric model 41 .More information about symmetric fuzzy model can refer to 42 .The optimization program is clearly equivalent to 3.5 .Using linear membership functions or piecewise linear, concave membership function, the system can easily be solved by wellknown algorithms.

3.5
Figure 2 shows the principle of the tolerance based symmetric fuzzy model.When the constraint relaxes to b 1 , the objective value minimized to Z 1 a significant improvement .However, when the constraint relaxes to b 2 , the objective value improves little Z 1 − Z 2 compare to the slack of constraint b 2 − b 1 .The tolerance-based symmetric fuzzy model is to get the maximum value of the fuzzy member λ, such that the system can obtains a great improvement of the objective on the conditions of less relaxation of the constraints.This means, for the train adjustment problem, we can get a new timetable to eliminate interference as much as possible with little slack of train rescheduling constraints.
The principle of the fuzzy model.

A New Timetable Rescheduling Model
In the rescheduling model mentioned above, four additional inputs are involved to describe the tolerances of minimum running time, dwell time, and headway time for section and station, respectively.An additional decision variable λ is used to describe the fuzzy number for the objective and constraints.All the additional parameters are listed as below: c a i,j : tolerance of a i,j .
c T i,j : tolerance of T i,j .
λ i : decision variable for the fuzzy number.Fuzzy membership functions of objective and soft constraints are given by μ 0 to μ 6 3.6 ; their graphics are similar to Figure 1.The timetable rescheduling model is changed as 3.7 .Some principles for the relaxation extent of the constraints are summarized as follows. 1 Collect statistics of train running states in all kinds of emergencies in reality. 2 Consult with the train drivers and dispatchers, who have a wealth of experience about train driving and operation.3 Refer to the train safety operation specifications under different emergency conditions.

3.7
Since M is a very integer, the headway time constraints can be changed as below: 3.8

Case Study
This model is simulated on the busiest part of Beijing to Shanghai high speed line, between Nanjing Ning for short and Shanghai Hu for short .In the rest of this paper, "Hu-Ning Section" is used to represent this part of the Beijing to Shanghai high speed line.There are seven stations on the Hu-Ning Section, thereby there are 6 sections whose lengths are 65110 m, 61050 m, 56400 m, 26810 m, 31350 m, and 43570 m.Its daily service starts at 6:30 am and ends at 11:30 pm.As currently planned, there are 52 trains 14 high speed trains and 38 medium speed trains from Beijing to Shanghai line and 8 extra medium speed trains from Riverside line that go from Ning to Hu by the Hu-Ning Section.The trains from Beijing to Shanghai line are called self-line trains, while those from Riverside line are called cross-line trains.In reality, self-line trains have higher priority than cross-line trains.As for self-line trains, the high speed trains have higher priority than the quasi-high speed trains.Since Beijing to Shanghai high speed line is double track, we only consider the direction from Ning to Hu without loss of generality.The experimental procedure is divided into two stages for the full proof of the validity of the fuzzy model.Firstly, we do research in six aspects with different fuzzy constraints of the same weight, which proves the effectiveness of the tolerance-based fuzzy model in different trains operation conditions.Then a sensitivity analysis of the weighing factors is realized based on the previous operation circumstances.
All the models are solved by Ilog Cplex 12.2.

Delay of One Train
In the simulation, we assume that the train G103 is late for 20 minutes in the section from Nanjing to Zhenjiang, and the trains can run at the normal speed in all the sections;  the headway is set from 3 minutes to 2.5 minutes, and the minimum separation time on track possession in each station is set to 1 min.
First we use the original model to solve the problem and get the follow results.When the speed of high-speed train is set to 360 km/h, the speed of medium high-speed train is set to 320 km/h, the headway time is defined as 3 minutes, and the solution value is 368784; there are 3 trains late with the delay time between 10 minutes to 20 minutes and 1 train late with the delay time between 20 minutes to 30 minutes; and the total delay time is 47.5 minutes.When the speed of high-speed train is set to 380 km/h, the speed of medium high-speed train is set to 350 km/h, the headway time is defined as 2.5 minutes, and the solution value is 350449; there are one train late with the delay time between 0 minutes to 10 minutes and 2 trains late with the delay time between 10 minutes to 20 minutes; and the total delay time is 31.85minutes.Then we use the above results as the inputs of the fuzzy model the paper described before and get the follow results.The fuzzy member is 0.86392, which means the average high speed is 365 km/h, the average medium high speed is 335 km/h, and headway time is 2.9319 minutes; there are also only one train late with the delay time between 0 minutes to 10 minutes and 2 trains late with the delay time between 10 minutes to 20 minutes; and the total delay time is 33.5 minutes.The adjustment strategy is to extend the dwell time for the train K101 and then make the train G105 overtake the train K101 at Changzhou North Station, which reflects the adjustment priority for high-grade train.Figure 3 shows detail rescheduling result by fuzzy model.The gray lines mean the initial timetable, and blue lines mean the rescheduling timetable.

Two Trains Delay
In the simulation, we assume that the train G305 is late for 30 minutes in the section from Nanjing to Zhenjiang, then the train G103 gets further delay for 20 minutes in the section from Zhenjiang to Changzhou under the condition of the existing delay, and other conditions are identical to Section 4.1.1.First we use the original model to solve the problem and get the follow results.When we take the strict constraints, the solution value is 860720; there are 2 trains late with the delay time between 0 minutes to 10 minutes, 3 trains late with the delay time between 10 minutes to 20 minutes, and 3 trains late with the delay time between 20 minutes to 30 minutes; the total delay time is 122.43 minutes.When we relax the constraints, the solution value is 848871; there are 2 trains late with the delay time between 0 minutes to 10 minutes, 2 trains late with the delay time between 10 minutes to 20 minutes, and also 3 trains late with the delay time between 20 minutes to 30 minutes; and the total delay time is 109.18minutes.Then we solve the fuzzy model and get the follow results.The fuzzy member is 0.89425, which means the average high speed is 362 km/h, the average medium high speed is 332 km/h, and headway time is 2.9512 minutes; there are also 7 delay trains, 2 trains of which late with the delay time between 0 minutes to 10 minutes, 2 trains late with the delay time between 10 minutes to 20 minutes, and 3 trains late with the delay time between 20 minutes to 30 minutes; and the total delay time is 114.29 minutes.The adjustment strategy is to make the high level train G105 subsequently overtakes the train K115 L15 and K101 at Zhenjiang West Station, and train L15 overtakes the train G103 at the Wuxi East Station, which effectively avoid the high level train G105 being late and reduce the delay time of the train L15. Figure 4 shows detail rescheduling result by fuzzy model.

Speed Restriction in All Sections
When some natural hazards happened, like heavy storm and strong wind, railway will be greatly affected in a large area.In the simulation, we assume there is a speed restriction in all sections and the average limited speed ranges from 170 km/h to 150 km/h; other conditions are identical to Section 4.1.1.
First we use the original model to solve the problem and get the follow results.When we set the speed as 150 km/h and headway time as 3 minutes, the solution value is 7969213; there are 2 trains late for the time between 30 and 40 minutes, 18 trains late for the time between 40 and 50 minutes, 34 trains late for the time between 50 and 60 minutes, and 6 trains late over one hour; and the total delay time is 3052.95minutes.When the speed is set to 170 km/h, and headway time is set to 2.5 minutes, the solution value is 7473849; there are 5 trains late for the time between 20 and 30 minutes, 38 trains late for the time between 30 and 40 minutes, 17 trains late for the time between 40 and 50 minutes, and no train late over one hour; and the total delay time is 2248.26minutes.Although the latter result is very exciting, it is a significant risk to take the speed of 170 km/h, as the highest speed must bring the high operation cost and may arouse some new delays.Then we use the above results as the inputs of the fuzzy model the paper described before and get the follow results.The fuzzy member is 0.689117, which means the average speed is 156 km/h, and headway time is 2.8445 minutes; there are 19 trains late for the time between 30 and 40 minutes, 25 trains late for the time between 40 and 50 minutes, 16 trains late for the time between 50 and 60 minutes, and no train late over one hour; and the total delay time is 2655.16minutes.We can see that the fuzzy optimization result improved greatly with little slack of constraints that means little risk and operation cost.Figure 5 is the rescheduling timetable by fuzzy model.In addition, the result is a parallel diagram since all trains must comply with the same speed.

Speed Restriction in One Section
Sometimes equipment failure, like train signal failure, happens in some but not all the sections.In the simulation, we assume there is a speed restriction in first section, and the limited speed ranges from 60 km/h to 50 km/h; other conditions are identical to Section 4.1.1.
First we use the original model to solve the problem and get the follow result.When we take the strict constraints, the solution value is 9294339; there are 19 trains late for the time between 40 and 50 minutes, 30 trains late for the time between 50 and 60 minutes, and 11 trains late over one hour; and the total delay time is 3269.32 minutes.When the speed is 60 km/h, and headway time is 2.5 minutes, the solution value is 8577661; there are 20 trains late for the time between 30 and 40 minutes, 37 trains late for the time between 40 and 50 minutes, 3 trains late for the time between 50 and 60 minutes, and no train late over one hour; and the total delay time is 2486.32 minutes.Then we get the follow result using the fuzzy model the paper described before.The fuzzy member is 0.70606, which means the average speed is 53 km/h, and headway time is 2.85303 minutes; there are 5 trains late for the time between 30 and 40 minutes, 39 trains late for the time between 40 and 50 minutes, 16 trains late for the time between 50 and 60, minutes and no train late over one hour; the total delay time is 2753.17minutes.The adjustment strategy is to make the train G117 overtakes the train K105 at Wuxi East Station, and then let the train L7 overtake the train G143 at the Wuxi East Station.Figure 6 is the rescheduling timetable with speed restriction in the first section.

Speed Restriction in One Section with One Train Delay
In the simulation, we assume that there is a speed restriction in first section, and the limited speed ranges from 60 km/h to 50 km/h; the train G103 gets delay in the section from Zhenjiang to Changzhou.Other conditions are identical to Section 4.1.1.

Speed Restriction in All Sections with Two Trains Delay
In the simulation, we assume that the train G305 is late for 40 minutes in the section between Nanjing to Zhenjiang and then get further delay for 15 minutes in the section between Zhenjiang to Changzhou; the train G103 gets a delay for 30 minutes in the section from Zhenjiang to Changzhou; also there is a speed restriction in all sections, and the average limited speed ranges from 170 km/h to 150 km/h for all the trains; other conditions are identical to Section 4.1.1.
First we use the original model to solve the problem and get the follow results.When we set the speed as 150 km/h and headway time as 3 minutes, the solution value is 8042152; there are 1 train late with the delay time between 30 minutes to 40 minutes, 17 trains late with the delay time between 40 minutes to 50 minutes, 27 trains late with the delay time between 50 minutes to 60 minutes, and 15 trains late with the delay time more than one hour; and the total delay time is 3197.86 minutes.When the speed is 170 km/h, and headway time is 2.5 minutes, the solution value is 7588658; there are 5 trains late with the delay time between 20 minutes to 30 minutes, 35 trains late with the delay time between 30 minutes to 40 minutes, 14 trains late with the delay time between 40 minutes to 50 minutes, 3 trains late with the delay time between 50 minutes to 60 minutes, and 3 trains late with the delay time more than one hour; and the total delay time is 2397.15minutes.Then we use the above results as the inputs of the fuzzy model the paper described before and get the follow results.The fuzzy member is 0.792877, which means the average speed is 153 km/h, and headway time is 2.89 minutes; there are 17 trains late with the delay time between 30 minutes to 40 minutes, 23 trains late with the delay time between 40 minutes to 50 minutes, 15 trains late with the delay time between 50 minutes to 60 minutes, and 5 trains late with the delay time more than one hour; and the total delay time is 2704.43minutes.The adjustment strategy is to make the train  L1 overtakes the train G101 at the Wuxi East Station, the train L15 and the train K101 overtake K115, and then extend the dwell time of K115 at the Wuxi East Station for the stopover by the train G101, so the importance of the high level trains is also shown from the steps above.Figure 8 is the rescheduling timetable by fuzzy model.
These simulations are realized on the computer with Intel Core 2 Duo CPU E7500 and 2 sG Memory.All the optimization results for the 6 cases are given in Table 1, where "R1" denotes original optimization subject to high speed, "R2" denotes original optimization subject to low speed, "R3" denotes fuzzy optimization, "D1" denotes the number of trains with a delay less than 10 mins, "D2" denotes the number of trains with a delay between 10 and 20 mins, "D3" denotes the number of trains with a delay between 20 and 30 mins, "D4" denotes the number of trains with a delay between 30 and 40 mins, "D5" denotes the number of trains with a delay between 40 and 50 mins, "D6" denotes the number of trains with a delay between 50 and 60 mins, "D7" denotes the number of trains with a delay not less than 60 mins, "D8" denotes the objective value, "D9" denotes the total delay time, "D10" denotes the number of stopover trains, and "D11" denotes the calculation time.

Different Weights
There are four kinds of fuzzy constraints as described in Sections 2.4 and 3.2.Since the minimum separation time on track possession in each station is set to 1 min in the previous experiments, we can ignore the uncertainty of station dwell time constraints.Furthermore, all the station and section headway times are set to the same value that from 3 minutes to 2.5 minutes; thereby the two constraints are considered equally important for the optimal objective and can use the same fuzzy member.Therefore two kinds of fuzzy constraints are analyzed in this section, section running time and headway time, and the weights of their fuzzy member are w 1 and w 2 .
Table 2 and Figure 9 are the comparison of the optimal objectives for the six cases.It can be seen from Figure 9 that the objective increases obviously as w 1 is set from 0.1 to 0.9 in the last four cases; however, the objective curves tend to flat in the former two cases, Case w 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 Case 1 0.860180 0.860920 0.861390 0.862250 0.863920 0.864470 0.865731 0.865020 0.864170 Case 2 0.889860 0.890440 0.891050 0.892280 0.894240 0.896030 0.895490 0.892661 0.890040 Case 3 0.671015 0.674356 0.678921 0.684357 0.689117 0.693258 0.698091 0.702655 0.705186 Case 4 0.697540 0.698920 0.700130 0.702480 0.706060 0.709138 0.712240 0.714521 0.715390 Case 5 0.730050 0.732160 0.734875 0.737280 0.740856 0.744030 0.747980 0.749721 0.751060 Case 6 0.779810 0.781530 0.784010 0.787903 0.792877 0.797460 0.801130 0.804655 0.806515 especially the curves decline slightly when w 1 is greater than 0.6.This may explain that the same constraints contribute differently for the optimal objective in different emergencies.Section running time constraint has a higher priority than headway time constraint when some accidents or natural disasters lead to speed restriction for a long time, and adjusting the train headway time is more effective if few trains delay by some interference.

Conclusion and Future Work
This paper presents a fuzzy optimization model based on improved tolerance approach for train rescheduling in case of train delay and speed restriction, which deals with the train running time at sections, the headway time at sections, and station as the fuzzy parameters.The simulations on Beijing to Shanghai high speed line reveal that the fuzzy optimization result improved greatly with little slack of constraints.This means we can get a new timetable with less total delay time as well as the number of seriously impacted trains in safe and lower average speed, little dwell time, and enough headway time.In addition, the sensitivity analysis of the weighing factors shows that the same constraints contribute differently to the optimal objective in different emergencies, thus the dispatchers should take different trains adjustment strategies to eliminate interference as much as possible.
There remain many interesting areas to explore around the uncertainty in timetable rescheduling problem.Firstly, the membership functions of fuzzy parameters used in the paper may be more complex form in practice than the linear function in the paper, so we can take some genetic functions, like gauss membership function, to model the fuzzy programming.The more accurate the membership function is, the better result the fuzzy optimization model gets.Secondly, the tolerance can also be described by fuzzy set-based schemes.Finally, fuzzy operators in the model can be improved to adapt well to the information processing mechanism of despatchers in dealing with rescheduling problems.Our ultimate goal is to develop a real-time rescheduling system to significantly improve operation management and scheduling efficiency in the future.

s 0 :
objective value of the original model based on original constraints.s * : objective value of the original model based on relaxed constraints.

Figure 3 :
Figure 3: Timetable of one train delay.

Figure 4 :
Figure 4: Timetable of two trains delay.

Figure 5 :
Figure 5: Timetable of speed restriction in all sections.

Figure 6 :Figure 7 :
Figure 6: Timetable of speed restriction in first section.

Figure 8 :
Figure 8: Timetable of speed restriction in all sections and two trains delay.

Table 1 :
Rescheduling results of different cases.

Table 2 :
Comparison of the optimal objectives.