GET THE APP

Optimized Train Dispatching and Rescheduling During a Disruption in a Bottleneck Section
..

Journal of Civil and Environmental Engineering

ISSN: 2165-784X

Open Access

Review Article - (2023) Volume 13, Issue 5

Optimized Train Dispatching and Rescheduling During a Disruption in a Bottleneck Section

Danson Byegon*, Sosina M Gashaw and Birhanu Reesom Bisrat
*Correspondence: Danson Byegon, Department of Environmental Engineering, Addis Ababa Institute of Technology, Addis Ababa, Ethiopia, Tel: 251942279643, Email:
Department of Environmental Engineering, Addis Ababa Institute of Technology, Addis Ababa, Ethiopia

Received: 28-Sep-2022, Manuscript No. JCDE-22-76181; Editor assigned: 03-Oct-2022, Pre QC No. JCDE-22-76181(PQ); Reviewed: 17-Oct-2022, QC No. JCDE-22-76181; Revised: 16-Dec-2022, Manuscript No. JCDE-22-76181(R); Published: 05-Jan-2023 , DOI: 10.37421/2165-784X.2023.13.520
Citation: Byegon, Danson, Sosina M Gashaw and Birhanu Reesom Bisrat. "Optimized Train Dispatching and Rescheduling During a Disruption in a Bottleneck Section." J Civil Environ Eng 13 (2023): 489.
Copyright: © 2023 Byegon D, et al. This is an open-access article distributed under the terms of the creative commons attribution license which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited.

Abstract

Railway transportation is nowadays becoming one of the most preferred mode of transport due to its safety, capacity and reliability; the capital cost for the construction of the railway infrastructure is however very high and is characterized by high rigidity as the track layout is fixed; therefore there is need to optimally use the available infrastructure. Minor delays arising from a simple disruptions or even a single train failure can have massive impacts in terms of overall delays for subsequent trains using the track facility if not solved amicably. Disruptions can be attributed to power outages, mechanical failures, derailments, accidents or even environmental factors.

In a case of multiple uncertain perturbations happening in a busy complex railway network, where there are many trains requesting to use the available track resources concurrently, there will be massive delays which has a lot of negative operational and economic implications as well as passengers’ dissatisfaction. A mathematical model that is; a mixed integer linear programming formulation is modelled to minimize total time delays in case of a set of multiple disruptions occurring on a busy track section i.e. bottleneck section, the model is formulated with consideration of sets of constraints factoring in feasible routes and safety margins and other operational dynamics such as dwell times to achieve optimal use of the available infrastructure.

A number of numerical experiments based on arbitrary data and real network data are carried out to verify the effectiveness and efficiency of the proposed model. Performance of the designed model is evaluated and results are validated, the results obtained shows that the model offers an efficient rescheduled trains operation plan during disruptions, furthermore the performance of Fmincon solver and Genetic Algorithms (GA) are compared and their robustness evaluated, GA shows better performance during multiple disruption scenarios.

Keywords

Delays • Disruptions • Infrastructure • Optimization • Scheduling • Train dispatching

Introduction

  From railway planning and operational point of view, research towards the development of optimal dispatching model for the transportation industry is deemed indispensable. Railway transportation is deemed to be one of the most reliable, efficient, sustainable and competitive mode of transport as compared to other means.

Despite these merits, at times the trains might not always arrive and depart as scheduled or even sometimes train operations are even cancelled due to uncertain and unforeseen events which can be attributed to external or internal failures in the railway system such as railway line derailment, severe weather conditions, locomotive breakdowns etc., this events may lead to underutilization of the available railway infrastructure. In such situations, Train dispatchers need to put in place some mitigation measures such as reordering the trains movements and routes utilization, rescheduling of train plans and retiming within a very short time to avoid impacts of delays caused from spreading further and avert conflicting route scenario from occurring.

During severe disruptions, the stopping pattern of some trains may also be changed to offer more options to the disrupted passengers to travel to their destinations [1,2]. The rescheduling activities should be done in a manner which considers passenger weights and distributions along the platforms and the allocation priorities given to the highly demanded trains basing on such factors like passenger demand and general passenger flow, other remedies to optimize operations in case of a disruption includes insertion of train services and even cancelations.

Delays encountered in railway systems can be categorized as either primary or secondary delays. Warg et al. defines Primary delays as the types of delays affecting a particular train schedule and those delays are not as a result of malfunction of other trains using the line. Whereas secondary delays are defined as the delays caused by earlier delays of other trains within the network sharing the same infrastructures, passenger transfers, transfers and reallocation of crew and rolling stock connections among others [3].

In solving resultant conflicts by train dispatcher; simple, reliable and efficient cost effective model needs to be defined and employed. This dispatching model will be helpful for dispatchers in ensuring that they make reliable decisions within a very short time frame and ensuring that there are no conflicting route scenarios. A lot of efforts have been put in place in regards to solving train dispatching or scheduling and rescheduling problems, however, the preference of use of railway system is increasing in demand and thus requires better solutions to achieve more reliable railway systems.

This article focuses on the integration problems associated with train rescheduling, train reassignment under multiple disruptions in a railway system particularly in a constricted or a “bottleneck” section, timetabling and delays experienced to come up with an optimum dispatching model to achieve the best use of the available infrastructure with minimal delays.

Literature Review

Problem description

Figure 1 below shows a disrupted route network, there are two disruptions which occur concurrently along the route, the causes of the disruptions could be different and uncertain or same causes of disruption but generally it affects the smooth flow and routine operations of the trains. To have a clear understanding of the context, we consider n number trains using the route as shown above, the trains are given the notations T1,T2,T3,T4,Txn,Tyn respectively. The track network has 4 stations in total with notation subscript X for the stations along eastwards flow route and script Y for the stations along westwards direction flow.

civil-environmental-illustration

Figure 1. Problem illustration.

The disruptions in the problem illustration, occurring in the line affects the scheduling operations and propagates the delays, for example, the train T3 cannot be dispatched from station By until the train T1 has passed through the disrupted section and a certain clearance time has been covered. To solve this problem, possible solutions would be rescheduling, reordering or even cancelling some train operations, any of the approach chosen should have sensible economical cost benefit resolutions and proper priority considerations.

The Bottleneck sections of mainline railways are usually characterized by high density traffic with many resource requests, where a disturbance from one single train could result in heavy impactful delays across a number of trains, multiple disruptions occurring concurrently in the busy section will in fact have detrimental impacts if no mitigation measures are taken may lead to underutilization of its infrastructural capacity and unprecedented delays. The model formulation is constrained to ensure that such conflicting situations does not occur after rescheduling the trains track allocations. Figure 2 presented below access the impact of disruption on a normal traffic flow.

civil-environmental-occupancy

Figure 2. Route occupancy during normal and disrupted scenarios.

A train i moving from station A to station E is shown in the Figure 2 above, the first blue line shows the movement of the train i as guided by the scheduled time table, whereas the orange line represents the rescheduled movement of the same train i in a path having multiple disruptions occurring between stations A,B,C and D, the duration of the disruptions are not uniform as shown by the time delays each train has to wait in the disrupted section, this shows the uncertainty index and weight of the disruption.

Related work

Train dispatching optimization is one of the fields of interest in terms of railway planning and operations and it has attracted extensive research in past few decades; several optimization models have been proposed, some of the cutting edge solutions proposed includes: Demand forecasting modules, discrete event simulation models, integrated models, multiple agent system, linear programming models and computer aided dispatching systems among others.

A number of previous works done on train dispatching problems has been able to assist in trying to solve problems faced by train dispatchers optimally. Some of the earliest works uses linear programming model and mixed integer linear programming as discussed by Higgins et al. in determining the best crossing and overtaking positions for trains on a single line track using fixed velocities and departure times of the trains [4]. Later on, a computer aided dispatching system for generating effective train schedules on single-line track using branch and bound algorithm was formulated. In addition, a high level petri nets simulation model that detects and manages train path conflicts on a single track railway line is presented by Jeremic D, et al. [5].

Over time several dispatching models have been developed employing greedy approach (heuristic method) in solving dispatching problems for timetabling and dispatching trains on single line track, heuristic algorithm for solving conflicts encountered by trains on single-line track that appear in time; Afeez et al. and Shakibayifar et al. proposed detailed information on the implementation of different heuristic algorithms for train dispatchers in solving train dispatching problems [6,7].

Bettineli et al. furthermore proposed a fast and effective parallel algorithm that is based on an iterated greedy scheduling of trains on a time space network [8]. Given a disturbed railway operations or timetable, a real time train dispatching and rescheduling solutions are suggested by Ariano et al. and who designed a dispatching system of decision support system for traffic controllers. Khwanpruk et al. also proposed a Timetable Optimizer (TO) consisting of demand forecasting module, train optimization module and timetable generator module [9,10]. Whereas developed a discrete event simulation model that employs a variable neighbourhood search algorithm to maintain the service level under infrastructure elements’ unavailability [11,12].

Due to the unpredictable nature of the incidents, the uncertainty associated with obstruction duration was experimentally analyzed through probabilistic scenarios. Barbour et al. proposed a methodological tool set, called the dispatch analysis problem that can analyse recent, empirical train dispatching data against an optimal dispatching plan [13]. In his work, Wang et al. presented an integrated model to optimize the train schedule and circulation plan simultaneously based on a given service pattern generated by the demand analysis and line planning [14].

Furthermore, Josyula et al. presented a parallel algorithm to efficiently solve the real time railway rescheduling problem on a multicore parallel architecture [15]. With the aim of improving simulation efficiency by directly confirming the next discrete time, Zhang et al. focused on discrete dynamic optimization on train rescheduling on single-track railway with the consideration of train punctuality and station satisfaction degree [16]. Takagi et al. also proposed the use of a Genetic Algorithm (GA) with Object Oriented Multi Train Simulator (OOMTS) developed by Birmingham University, as an embedded simulator to optimize the order of route setting [17].

Wang et al. and Gong et al. studied an integrated train rescheduling and rolling stock circulation planning problem for the complete blockage situations in a metro line, by modelling the problem as a complex multi-objective MILP formulation [18,19]. Withal Leonardo et al. also suggested the formulation of a mixed integer linear programming model using Benders'-like decomposition approach to minimize the number of (heavily) delayed trains in cases of disruptions [20]. Guo et al. proposed a Multi-Agent System (MAS) to reduce rail transit energy consumption when disturbances occur. Furthermore, Huang et al also discussed three models in a progressive fashion to optimize the last train timetable incorporating multimodal coordination.

Railway traffic management

There is a need to offer quick railway traffic rescheduling to ensure least disturbance along the line operations when a disruption occurs. One of the mitigation measures to this problem was proposed by Shakibayifar et al. who designed a simulation based optimization model that incorporated dynamic dispatch priority rules with the objective of minimizing the total delay time of trains and a variable neighbourhood search meta-heuristic.

A railway transport system is a very complex system with integration of human controls and behavior, technical devices, and external environment integrated together and thus contains many risks of disturbances. Railway traffic performance are usually managed by a carefully designed plan of operations also referred to as train timetable, the timetable which incorporates orders and the timing of the trains several months in advance is referred to as offline timetable. There is however usually a secondary remedy to offer real time measures to manage traffic in cases of disturbances.

Railway bottleneck

A bottleneck section is a section in a railway network consisting of two or more junctions points at each end and it is usually a very busy section as many trains going in different direction uses the section. This definition of the bottleneck shall be used in our problem modelling in the subsequent sections of this paper. The rescheduling of approaching trains for one portal junction needs, therefore, to consider the movements of the trains approaching this portal junction from another portal junction, looked at the junction designs and efficiency.

Railway timetabling

A timetable is the foundation of any public transportation system. It is critical in determining the paths or routes to be taken and the travel times from the points of origins to the desired destinations. Designing of train timetable is a complex process which must take into consideration the feasible infrastructural capacity, safety assurance, and resource constraints. The extent to which disturbances such as delays propagate through the system is highly depends on the robustness of the timetable.

The urge to come up with a good and reliable timetable has triggered many researchers to study route planning and the field of public transport optimization Lieshout et al. Rail traffic planning process involves strategic and tactical levels. Strategic planning entails railway network planning and line planning. Whereas, tactical level is aimed at ensuring efficiency and reliability and it entails timetabling, rolling stock planning and crew scheduling.

Railway timetables can be grouped into two categories as described by Gao et al. namely:

• A regular timetable which allows passengers to easily remember the exact departure times at stations.

• An irregular timetable, which can better match passenger demand and reduce operating costs and minimize passengers’ cost for their trips, and enhance passengers train-travel experience Xie et al. Skip-stop plan; is a strategy commonly used to increase vehicle speeds and reduce passengers travel times as explained by Mei et al. and Tang et al.

Route disturbance

Railway traffic is operated according to a detailed off-line plan of operations, specifying for each train its path through the network plus arrival and departure times at its scheduled stops and studied how to support dispatchers in the management of infrastructure disruptions by considering a seriously disturbed traffic conditions on double track railway lines where some block sections of one track are unavailable for traffic, e.g., due to a temporary track blockage.

In emergency perturbations occurrence, there is strict time limit feasibly available to compute a new timetable to satisfactorily gather for all the rail network needs particularly if the perturbation occurs on a very busy section and in most of the cases operators usually restrict themselves mostly to a few manual modifications of the timetable. Ariano et al. looks into how network dispatchers regulate the railway traffic. Existing dispatching systems are able to provide viable solutions only for small instances or for simple perturbations. However, the real-time traffic management is still under the control of human dispatchers, in most cases, who usually do not have precise information about the future evolution of train traffic and the chosen traffic control actions may be often sub-optimal.

Train rescheduling

A scheduler is used to generate feasible timetable for the railroad system, it allocates departure and arrival times for all trains within the network and also tracks needed for their journeys. Some of the major scheduling constraints are described in and they include; speed constraints ensures that a trains within the network should not exceed the maximum allowed or maximum safe speed specified for it to use on any track during their journey. Furthermore, station entrance constraints ensure that there must exist fixed time interval between two trains arriving a stations using same track line. Whereas, load or weight constraints checks on the weight of the locomotives and ensures that the weight of trains specified to run on any track within the network must not exceed the maximum allowed weight resistance of that track.

Timetable conflicts may occur when the trains within the network does not abide by the scheduling rules, possible causes of these problems can be breakdown of trains or stochastic events. Possible timetable conflicts are discussed by and are describes as; trains meeting conflict, station capacity conflict and trains headway conflict. These conflicts can cause problems within the network and needs to be solved to avoid accidents (Table 1).

Table 1. Review summary

No. Identified gap Objective function Methodology
1 Need to design a predictive model to determine the priority trains for a network and route optimization. Minimize total delay time of the trains Meta- heuristic scheme approach
2 Integrating optimization model with other meta- heuristics to enhance search procedures. Real time traffic scheduling and generate Pareto frontiers. Multi-objective simulation.
3 Modelling a demand- driven timetable which is dynamic and autonomous in nature. Joint optimization of timetable and vehicle circulation schedule MIP
4 There is need to develop a re-scheduling plan taking into account passengers response during random disruptions. Generating resilient operation plans during random disturbance. MINLP and GA approach.
6 Need to handle complex rail nodes. Minimize the number of (heavily) delayed trains MILP using Benders'-like decomposition
7 Looking into multiple disruptions occurring in a bottleneck section. Minimizing total train delay times. MILP using Fmincon and GA solvers

Rescheduling model development

Train rescheduling during a disruption is considered to be one of efficient mitigation remedy to ensure that train services continue operating using the available feasible routes, even though it can be a very complex and challenging task. The rescheduling model formulated for this purpose is a mathematical model as discussed below.

Mathematical model formulation

A mathematical model can be so helpful in coming up with an efficient and reliable rescheduling solution. In formulating a mathematical model; the first step creating a conceptual mathematical model that captures the important resources and constraints in the train dispatching and routing problem. The conceptual mathematical model can then be translated into a computer program that builds that model and solve the model using a math programming solvers (Tables 2-4).

Notations:

Table 2. Notations.

Symbol Definition
Qi Set of train i i.e. Q= {1,2,3,…,n}, i is the index for trains
Ek Representation of events, k is the index for events
Rj Routes representation, j is the index for routes sections
Gx,y Set of direction of the train, x is the index for trains in eastwards direction and y is the index for the trains in westwards direction.
x Directional notation for train i moving in eastwards direction.
y Directional notation for train i moving in westwards direction.
Sn,x Stations along the east-west direction.
Sn,y Stations along the west-east direction.
Rti,x Real departure time of train i in X direction during disruption.
Rdi Real distance to be travelled by train i in the disrupted section.
Sti,x Scheduled departure time of the train i in X direction during normal operations.
t0i,z Arrival time of train i at destination point (z) in the nominal timetable i.e. scheduled arrival time when there is no disruption.
ti, z Real arrival time of train i at the destination point (z) when there is a disruption.
tw Represents the waiting time the train i going in eastwards direction will have to wait for another train going in opposite direction(westwards) to clear off the shared section i.e. the feasible use the section route before train i can proceed to occupy and use the section
tm H Minimum headway for the trains moving in the same direction.
t D Trains dwell time at the respective stations.
Vds Speed of the train i when maneuvering through the disrupted section using the alternative feasible route. This speed is usually lower than the actual operating speed (Vr) during normal operations due to the curvatures of the switch points.
Vr Actual train speed during normal operations without disruption.
B Set of block sections i.e. bi={1,2,3,…,n}; i ∈ Q, b ∈ B
Mj Set of events on a section j. Mj= {(i1, k1), (i2, k2), (i3, k3)…}; i ∈ Q, k ∈ E, j ∈ R

Table 3. Definition of terms.

Indices Definition
d ϵ D Set of distances travelled by train i i.e. D={d1,d2,d3,…,dn}
r ϵ R Set of routes.
i ϵ Q Set of trains, where Qx is a set of trains going towards eastwards direction and Qy is set of trains going towards west direction.
g ϵ G Set of directional priorities, this is guided by the notations x and y as defined above.
h ϵ H Minimum safety headway t m, d; distance between two subsequent trains. H
v ϵ V Set of acceptable speed profiles.
b ϵ B Set of block sections.
j ϵ R Route indices. Feasible Route is composed of block sections.
t ϵ T Set of running time intervals for train i.
s ϵ S Set of stations along a particular given route i.e. S= {s1,s2,s3,…,sn}

Table 4.Set of indices.

Parameter Definition
Rti,k begin Actual start time of an event k for train i during a disruption.
Rti,k end Actual end time of an event k for train i during a disruption.
Sti,k begin Planned start time of an event k for train i during normal train operations as scheduled in the timetable.
Sti,k end Planned end time of an event k for train i during normal train operations as scheduled in the timetable.
Δj Minimum clearance time for the train already using the requested route j.
Gi,k Direction of the train i during an event k
Rd,i Actual distance travelled by train i
Vrij Actual speed of train i using the feasible route indexed by j

Conceptual mathematical model entails:

• A set of decision variables, which capture the range of possible decisions to be considered are defined.

• A set of constraints, which capture global limits on the values of the decision variables are defined.

• An objective functions on the decision variables that captures the value to be minimized or maximized.

In creating a mathematical model to offer optimum rescheduling operations, the problem illustration in Figure 1 above is considered, the objective function is to minimize total train delay time, the constraints and other parameters such as decision variables shall be discussed in subsequent sections:

The diagrammatic representation of steps in optimization model development is shown in Figure 3 below.

civil-environmental-solving

Figure 3. Steps in solving the model.

Parameters terms definitions:

Routes/Available blocks: Feasible routes in the disrupted section represented by blocks represented by b1, b2, b3, b4.

Trains: Trains waiting to use the disrupted section represented by Qix, Qiy. i.e. train i.

Trains speed: The permissible speed profile allowable in the disrupted section i.e. Vd and actual train speed during normal operations, i.e. Vr.

Direction: The directional flow provisions of the trains. Eastwards directional flow is named X direction and westwards directional flow is named Y direction.

Passenger and priority train allocation weights: Represented by Wp.

Number of disruptions: represented as μk.

Weight μk: Impact of the disruptions on the routine trains operations and timetable plan.

Clearance time: It is the time the trains will have to wait for the feasible route to be available for their use, i.e. the time it will take a train currently occupying a given route section to clear off.

Weight q: Maximum number of trains waiting to use the available route in the disrupted bottleneck junction section.

Station dwell time: The time a train will stop at a given station for the passengers to board or alight from the train; this is determined by the passengers’ distribution weights.

Scheduled Time (St): specifies the minimum running time (computed prior to the problem solving) i.e. scheduled running time for event k if event k is an event for a segment between stations.

Decision variables

• Assigning route r to event E if train i request to use the section pi,E,1= ε {0, 1}; 1 if event E of train i uses the route section r; 0 otherwise.

• Select route r for set qi trains ϵ {0, 1}; 1 if the minimum safety headway tmH constraints, speed constraints (v), time constraints (t) and directional constraints G; are met; 0 otherwise.

• Assign qi, r ϵ {0, 1}; 1 if there is feasible route r available for trains s et q, clearance time provisions are met and only if there is no conflicting route scenarios; 0 otherwise.

• Vrij=Actual speed of train i using the feasible route indexed by j.

• Δj=Minimum clearance time for the trains moving in different directions to complete using the section so that the section can be available for use by other trains moving in different directions.

• For each event (i, k,) ϵ Mj, we have the decision variables pi, k,1, pi, k,2,… pi, k,n, where pi, k,1 is a binary variable that specifies whether route r (associated with j) is used by event E of train i.

Objective function

The objective function is to minimize the total time delays encountered by the trains using the disrupted section of the bottleneck junction (i.e. the sum of the final delays when trains arrive at their final destination), which is represented in equation (1) below.

Equation

No wait constraint: This constraint ensures that the resource allocations start immediately the previous one ends.

Equation

Direction constraints: This constraint serves to balance the traffic flow in the disrupted section to ensure that the flow is balanced in both directions while taking into account the priority considerations.

The directional constraints can be equated as:

Equation

Run time and dwell time constraints: This constraint ensures that the new event starts after the running time and the station dwell time of the previous event has been completed.

Equation

The running time of the train within a particular section of the railway is guided by the availability of feasible route with no other train using the section at the time of allocation.

Equation

Running time check: The main task of the running time check is to find the invalid solution individuals that do not satisfy the running time restriction shown in equation (8)

Equation

The running time check rule is shown in equation (9) below; If

Equation

At the moment of disruption, if the train is still not dispatched from the origin, it must satisfy the departure time interval from the origin. This is represented in the constraint below:

Equation

The travel time (or stop time) of the trains in the block (or station) are variable but should be defined by the minimum and maximum ranges

Equation

Track occupancy constraint: This constraint ensures that one train occupies a particular route section at any given time. This can be expressed as;

Equation

Route constraints: The trains can only proceed if there is no conflicting routes scenario.

If a single-track section is blocked, then the arrival time to the block for all incoming trains should be after the completion time of the blockage (t1+τ1) i.e.

Equation

If the blockage occurred in one of the lines of a double-track block section, then during the disturbance, and the partial service is accepted. This is expressed by the constraints below:

Equation

Where Sn is the next station after disruption point and Si is the station before disruption point.

The route setting control sequence can be defined as follows, using 0 and 1 to denote the two possible routings i.e. 1 show that there is a feasible route and 0 shows that there is no feasible route.

Routes in eastwards direction is represented by X and the routes in the westwards direction is represented by Y.

For the eastwards direction:

Equation

For the westwards direction:

Equation

To ensure that not more than two routes is allocated to one train i, the constraints equation (22) and (23) implements that condition

Equation

Optimal route setting for train i can be summarized as shown in equation (24) below

Equation

Headway constraints: The constraint ensures that the trains moving in the same direction must stick to the safe minimum headway requirements for safety purposes.

The minimum safety headway must not be exceeded.

Equation.

Rtix(sn-s1) is the running time of train i two subsequent blocks or stop overs i.e. stations Sn and S1.

Sequence and headway admissibility check: This constraint is used to ensure that solution individuals generated by the stochastic algorithm are valid in terms of the train sequence restriction at control boundaries. For example, if there is only one route for trains from point A to point B, then the train arrival sequence at point B must be the same as that at point A.

Equation

t s2-t s1 is the time taken between station s2 and station s1,(ts3−ts4) is the time taken by train between station s3 and s4.

Speed constraints: This constraint ensures that the maximum allowable speed profiles must not be exceeded. There is maximum allowable speed of the trains when running in the disrupted section must not be exceeded due to curvatures i.e.

Vr ≤ max permisible speed (Vp)

The speed of the train must also no be lower than a certain allowable threshold speed in order to utilize the infrastructure capacity in a relatively optimal manner.

Therefore, considering the two speed conditions, the speed constraint is presented as;

Equation

vpmin is minimum allowable speed for the train to cruise on the track.

For this particular formulation, there are some time delays which we assumed to be negligible and are not exclusively calculated, i.e. the time taken in setting and resetting the switch points, as well as stochastically varying delays at the stations which might be attributed to dwell time inconsistencies, we consider this time delays as very minimal and thus inconsequential in the overall delay encountered by the train running through the disrupted section.

To solve the complete MILP problem presented above (composed of objective function (1) and constraints presented in the subsequent equations from equation (2) to equation (27), may require a lot of computational effort to simply this, the absolute term in the objective function can be linearized.

Validation test scenarios

To validate the formulated mathematical model, multiple disruptions scenarios occurring in the rail network are considered. For this purpose, a numerical example is used. Assume a rail network as illustrated in Figure 4 below.

civil-environmental-validation

Figure 4. Context for validation

The first and the second disruptions are induced in the segments between stations STx1 and STx2 and between stations STy3 and STy2 respectively in different directions. The third disruption happens between stations STx3 and STx4. The disturbance periods vary i.e. the first disruption starts at time T1 and the second disruption starts at time T2 and the duration of each disruption is

τ1 and τ2 respectively. The data required for numerical experiments is presented in subsequent sections. The departure time tolerance has been defined in order to achieve the best initial plan. However, in this example, the choice of departure time interval is arbitrary.

Without loss of generality, it is assumed that the segment distances are equal i.e. stations are 1 KM apart from each other. The running time between adjacent stations is calculated using velocity time formula taking into consideration speed limits restrictions and the network is a double rail and each section has switches and two rail tracks that permit trains to cross as illustrated in Figure 4 above.

The problems presented above are typical hybrid optimization problem that includes continuous variables (train arrival time) and discrete variables (route setting decisions). The discrete route setting decisions are identified once availability of a feasible route has been confirmed. If there are constraints for the train sequence on the converging routes, the route setting at the junction point must be decided on the basis of the train sequence demands. Also, the calculation of the train arrival time must take into account the route setting constraints.

Results and Discussion

Scheduled train timetable

The scheduled plan of the train operations during normal network operations without a disruption is presented in the state space chart below:

Figure 5 above show the movements of trains under normal scenario when there is no disruption, there are eight trains represented by the color codes and a total of six stations denoted by ST. The upstream graphical representation shows the trains moving in eastwards direction and the downstream graphical representation show the trains moving in westwards direction. During normal trains operations, the traffic flow is smooth and the traffic movements are guided by the scheduled timetable.

civil-environmental-timetable

Figure 5. Scheduled train timetable during normal operation.

Disrupted route scenario

Due to the uncertain nature of the disruptions, different possibilities for the start time, duration of delay and the location of the disturbance are probable. Different disruption scenarios are induced in the network for this purpose, the induced disruptions are shown in Table 5 below. A case of a single and double disruptions scenarios is considered, with a the pair of double disruptions occurring concurrently, and the disruption starts just before the most congested periods of the day, in the bottlenecks with the highest traffic as computed by a traffic congestion factor defined as the number of train conflicts in a schedule during a specified interval.

Table 5. Disruption scenarios.

Induced disruption Duration of blockage ( hours) Location of blockage Number of disruptions Direction of traffic disrupted
D1 0.5 Between points A-B 1 Eastwards direction
D2 0.75 At point B 1 Westwards direction
D3 1 Between points B-C 1 Eastwards direction
D1: D2 1.25 Between points A-B 2 Both directions
D1: D3 1.5 Between points A-B;B-C 2 Eastwards direction
D2: D3 1.75 Between points B-C 2 Both directions
D1: D2: D3 2.25 Between points A-B; B-C 3 Both directions

Choice of the algorithm

Matlab optimization algorithms i.e. Fmincon solver and genetic algorithm are used to optimize the continuous real time train arrivals during different disruption scenarios compared to scheduled arrivals times and their performance compared, taking into account the objective function, defined constraints and discrete route setting decisions variables.

Comparison between Fmincon and GA Solvers optimization during different disruption scenarios.

The performance and robustness of the two algorithms used are compared in Figure 6 and it is noted that Fmincon solver algorithm and GA solver offers almost similar optimal results in terms of simple or single disruption scenarios although there is some small variance. Whereas, for multiple and more complex perturbation cases i.e. D1D2, D1D3, D2D3 and D1D2D3, GA solver offers better results as compared to Fmincon solver, this is attributed to the nature in which the two solvers approach and defines the problem; Fmincon solver applies the use of a Successive Quadratic Programming (SQP) method. In this method, a Quadratic Programming (QP) sub problem is solved at each iteration.

civil-environmental-genetic

Figure 6. Comparison between Fmincon and genetic algorithm solvers

Genetic algorithm offers a global search solution beginning with a randomly chosen assortment of chromosomes, which serves as the first generation (initial population). Then each chromosome in the population is evaluated by the fitness function to test how well it solves the problem at hand. Now the selection operator chooses some of the chromosomes for reproduction based on a probability distribution defined by the user. The fitter a chromosome is, the more likely it is to be selected. This explains better performance and robustness of GA algorithm in a complex and uncertain disruptions.

Optimized reschedule during a single disruption

Figure 7 presented below show optimized trains movements in case of a single disruption occurring in a rail network, the Disruption case occurred between station STx1 and STx2 (Tables 6 and 7).

civil-environmental-disruption

Figure 7. Rescheduled plan during a single disruption.

Optimized reschedule in a case of double disruption: Figure 8 presented below show an optimized trains movements in case of double disruptions occurring in a rail network, the disruption cases occurs between station ST1 and ST2 and between ST3 and ST4 respectively with different start times and different disruption durations as shown graphically below. It can be seen that the train takes a longer time to move through the disrupted section as compared to unperturbed sections.

civil-environmental-rescheduled

Figure 8. Double disruption rescheduled traffic flow.

For this particular study, seven different random disruption scenarios presented in Table 3 are used to evaluate the mathematical model designed. Tables 6-7 above shows the optimal schedules for each train and optimized variable parameters respectively in each case scenario, the objective value for each disruption scenario is obtained.

Table 6. Individual train optimized schedule.

Train codes Objective value Travel time during normal scenario % Increase in travel time during disruption
T1 146 739.897 19.73
T2 152.121 857.831 17.73
T3 178.965 977.046 22.89
T4 346.954 1096.261 31.65
T5 238.09 738.91 36.84
T6 57.55 857.149 7.01
T7 40.654 976.346 4.16
T8 137.457 1095.543 12.6

Table 7. Optimized parameters.

Parameters Disruption case scenarios
D1 D2 D3
RTX A-B 91.145 51.43 102.164
RTX B-C 104.285 51.43 102.255
Headway X direction - - 180
Clearance time X direction - - 46.973
Stations Dwell time in X direction - - 40.467
RTY C-B 147.842 242.046 -
RTY B-A 147.842 234.286 -
Headway Y direction 180 180 180
Clearance time Y direction 141.706 93.437 -
Stations Dwell time in Y direction 39.47 3.931 -

Figure 5 shows a train traffic flow during normal operating conditions when there is no disruption, whereas Figure 7 and 8 shows a traffic flow state space graph during a single and double disruption respectively. From the observation in both cases, it is observed that more total travel time delays are encountered for the highly disrupted scenario with multiple disruptions as there are few feasible routes with main trains requesting to use the available resource. It is also observed that the nature and the location of the disruption along the railway line have an impact on the total time delays. The performance and robustness of the two algorithms used are compared in Figure 6.

Conclusion

In conclusion, a mathematical model that is able to minimize the total delay time is designed and tested, the results shows that the model is able to achieve its objective by adjusting and offering optimal time delays for each individual parameter contributing to the total train time delays for the train using the disrupted section, the model is also able to offer a rescheduled plan in a case of a single or multiple disruptions occurring concurrently at the railway network.

A small delay resulting from a single train or a simple disruption can result to massive impacts when it comes to utilization of the available railway capacity, the disruptions could be as a result of mechanical failures, power outages, derailments, accidents and even other uncertain and unprecedented causes, this disruptions generally alters the normal operational plans in terms of time table scheduling, passengers and freight transportations, therefore to mitigate the impact of this disruptions, there is need to minimize time delays resulting from the mentioned disruption. This article developed and validated a mathematical model to minimize time delays from the resultant disruption, the results from the validation confirms that the model offers a fast rescheduling and route allocation approach with minimized total weighted delays.

This model furthermore compares solutions using two different algorithms and from the results, there is no much variance in terms of the best solutions however each algorithm approaches the problem differently i.e. Fmincon solver looks at the scalar variables with respect to the constraints whereas the genetic algorithm solver modify the problem into a population approach and offers its search solutions using mutation approach, genetic algorithm is found to offer more optimal solutions in a context of a global problem search, whereas Fmincon solver is suitable when solving linear scalar problems.

It is recommended that the railway engineers should design and incorporate switch points at appropriate locations between the bottleneck sections to enable trains change of tracks in case of unprecedented disruptions so as to offer alternative feasible routes for the trains.

The future study aims at using artificial intelligence to offer faster and trained solutions to train rescheduling problems. Furthermore, looking at larger disruption contexts will be more ideal for a more realistic contexts as railway transport is becoming a more preferred mode of transport. Furthermore, the designed mathematical model can be used to analyze other performance parameters for instance stability and punctuality of rail operations. In addition, developing a simulation model that considers the main uncertain factors in practical applications and looking at the relationship between train volume and dispatch performance relative to an optimal schedule and its impact on capacity and infrastructure planning can be good research area in the future.

References

Google Scholar citation report
Citations: 1798

Journal of Civil and Environmental Engineering received 1798 citations as per Google Scholar report

Journal of Civil and Environmental Engineering peer review process verified at publons

Indexed In

 
arrow_upward arrow_upward