1. Introduction
Recent years have seen an accelerated development of high-speed railway (HSR) in China and many other countries to meet the increasing demand for travel and the socio-economic developments in these countries. This rapid development poses huge challenges to railway infrastructure managers and operators to best manage large HSR networks and to produce the most effective operating plans. Railway planning generally follows a hierarchical top-down approach, starting from network planning, to line planning, timetable generation, train routing, rolling stock scheduling, and crew scheduling [
1].
Within a given network, line planning becomes the first stage of the railway planning process, which is the process that determines the train routes, frequencies, and stop patterns in the network. Generally, train routes are determined first. Then, train service frequencies can either be optimized with the pre-determined stops or be simultaneously optimized with stops [
2]. Alongside train route and frequency, train stop planning directly represents the service level to stations and for passengers.
An effective stop planning is to balance the trade-off between railway supply and passenger demand; the latter is usually characterized as passenger flow levels through a station, or between an origin-destination station pair. A stop plan with a large number of stops will clearly be able to serve more passengers more conveniently but at the expense of increased train travel times and energy consumption caused by acceleration and deceleration. On the other hand, a stop plan with fewer stops will lead to longer waiting times for passengers but shorter train running times. Such practical supply-demand trade-offs increase the complexity of line planning or stop planning.
Most of the studies on line/stop planning have attributed travel demand to stations and categorize the station importance (or level) by the volume of passengers through the station—stations with higher degrees of importance correspond to higher passenger volumes. Then, the stop planning problem is to address the different service requirements for different stations. A commonly adopted approach in line/stop planning for station categories is to select from a set of lines with pre-determined routes and stop patterns that satisfies certain objectives [
2,
3,
4,
5]. Fu et al. [
2] classified stations and trains into two levels and applied the concept of “possible line set” to define train routes and stop stations for line planning, where the intermediate stops were pre-set. Kaspi and Raviv [
5] presented a model that combined line planning, timetabling, and routing decisions of the passengers. Again, they only selected lines from a predefined set of lines, each of which had its route and stopping stations pre-determined. A key drawback of the line set approach is that the candidate lines have fixed stopping stations. This makes it inflexible in line planning, and it may not always find the desired service frequencies for origin-destination (O-D) pairs.
To ensure the required level of service to a station, a minimum service frequency can be set based on the degree of importance of the station [
6,
7,
8,
9,
10,
11,
12]. Yang et al. [
6] proposed a mix-integer linear programming model to integrate train stop planning and train scheduling problems. The constraint that each station must be served by a given number of trains was considered. With given numbers and types of trains, Chen et al. [
7] proposed a stop-schedule optimization model that jointly minimized the train stop cost and maximized passenger travel convenience. They considered the reasonable range of train stop number and the minimum service frequency for each station. In addition, there has also been extensive literature on the design of stop patterns to cater to different demand at stations [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22]. Patterns include all-stop pattern, skip-stop pattern, zonal pattern, and express/local pattern, which cater to different demand characteristics (see a comparison in reference [
13]). Each stop pattern has the specific applied area, so they are not suitable to be adopted for a general stop planning problem. Thus, some researchers also study a combination of several patterns [
23,
24,
25,
26,
27,
28]. However, the complexity and computation time required to solve a combined stop-pattern problem become significant as the number of variables to optimize increases, so this combined approach is not suitable to be applied to large networks unless an efficient algorithm can be found.
Overall, the above literatures in stop planning are concerned with the service levels at stations for simplification, i.e., the passenger volumes through stations. They do not specifically deal with passengers’ O-D travel.
A good stop plan should consider travel characters and satisfy the passengers’ direct travel requirements between their O-D stations. Different O-D demands need different levels of service, and a certain number of train services need to simultaneously stop at the origin and destination stations of each O-D pair, not single stations. Therefore, the ultimate purpose of stop planning is to provide a sufficient level of direct services for different O-D demand by determining the stop stations for each train on railway lines to reduce the transfer costs and travel times of passengers.
There have been studies, though limited, in the literature considering the direct choices of passengers in stop planning. Lin and Ku [
29] constructed an integer program model to determine stop patterns and stop intervals for trains without passenger transfer. The objective was to maximize the profit earned by the train company, and the model was solved by two genetic algorithms. Lai et al. [
30] designed an innovative stop planning network to formulate the stop sequences of trains and developed two binary integer programming models to tackle scenarios without and with passenger transfer. However, these two studies only considered a single O-D demand level and did not present a specific level of service (service frequency) for different O-D pairs. To ensure the defined level of service for O-D demand, Yue et al. [
31] incorporated the service for major O-D demand in train stop schedule optimization, in which the major O-D demand were defined by the importance of their O-D stations. They set a constraint to ensure minimum percentages of direct services for these major O-D pairs. The minimum service requirements for major O-D pair were given as input to the model, and the authors suggested that the values could be generated from empirical studies.
The decisions made in the above studies are in terms of whether a train is to stop at a station or not, of which the decision variables are discrete 0–1 variables from the point of view of service for single stations, and the required service level (i.e., the service frequency) is determined by empirical methods. Thus, there has not been a theoretical approach to analytically derive the service frequencies for different O-D demand levels. In practice, trains would not stop at each passing station. Stations on a higher level of importance would generally have more stopping trains; in other words, trains would have a higher probability to stop at a station with higher level of importance. Thus, we can formulate stop planning from another perspective. Different from the binary integer decision variables in the mentioned studies, we propose the concept of stop probability, i.e., the probability of a train stopping at a station that is approximately equal to the ratio of the number of stop trains to that of passing trains, which is a continuous variable to represent the relative importance of a station in stop planning. This service probability for stations or O-Ds is usually adopted in practical planning stage by empirical methods, when operators set required train percentages for major stations or O-D pairs [
31].
This paper formulates stop planning from the point of view of direct travel between O-D stations. Considering different O-D demand levels and train service types, we introduce the concept of stop probability to present the general principle and mathematical formulation for stop planning with the objective of minimizing the per capita travel time, and determine the service frequencies for each O-D pair according to demand characteristics. The comparison to existing literature is shown in
Table 1, and the key contribution of this paper is summarized below.
First, we characterize railway stations in the HSR network by their geo-political importance, which also represents the attraction effects to travel demand. We consider spatial O-D movements that cross different, and movements within the same geo-political areas to present the travel characters of different O-D demand, and then define the O-D pairs jointly by the respective O-D station characteristics and by their spatial movements.
Section 2 presents the detailed description of the geo-political classification of stations and the dual-level characterization of demand.
Second, we introduce the concept of stop probability of trains by station level and formulate the stop planning problem as to optimize trains’ stop probabilities at stations on different levels. We show that the service frequency for a specific O-D pair can be derived directly from the stop probabilities at the respective O-D stations. We also consider two types of trains with different service availabilities—type 1 trains can stop at any station, while type 2 trains would only stop at higher level stations.
Third, we consider two additional constraints in the optimization to obtain feasible solutions—equilibrium of train load and equilibrium of stop density. Equilibrium of train load means that the proportions of train kilometers of the two train types should be similar to those of their transported passenger kilometers to keep a balance between the loads of the two types. To ensure the well-distributed stop stations within a geo-political area for the two train types, we adopt the concept of stop density, i.e., the average stop number in a province/district, to describe the distribution uniform degree of the stop stations of the train types, and keep it in a reasonable range.
We present the analytical formulations for the derivation of the optimal service frequencies for different O-D demand categories on an HSR line, which are important inputs for line planning and stop planning. From the point of view of direct travel between O-D stations, we formulate the stop-probability based stop planning problem to minimize the per capita travel time. The resulting optimal stop probabilities are then used to calculate the number of trains required for each train type serving the different demand categories. A case study on three HSR lines in China validates the effectiveness and correctness of the approach, and the solved results show a more reasonable load division between the two train types.
The remainder of this paper is organized as follows. The stop probability optimization problem is formulated in
Section 2, and the classification of stations and demand is also introduced in this section.
Section 3 constructs a stop-probability based stop planning optimization model, and the solving method is also designed. The case studies of three HSR Lines are presented in
Section 4, and conclusions follow in
Section 5.
2. Problem Formulation
In this section, we first introduce the classification of stations, demand and trains, and then formulate the train stop probability optimization problem.
Table 2,
Table 3 and
Table 4 show the input parameters and variables in this paper.
2.1. Station Categorization
In this paper, the concept of stop probability represents the importance of stations on the network, which can also be presented by station levels. Thus, in order to present the station characters and reduce the scale of decision variables, we next characterize railway stations and categorize them into different levels, which also influence the attraction effects to travel demand.
Stations on an HSR network can be classified into different levels according to different criteria [
17,
27,
28,
29], including administrative level, location on the network, daily dispatched passenger number, the number of passing trains, etc. In China, an important classification criterion of stations is the geo-political (or geo-administrative) level of the locations of stations, which also represents the attraction to passenger demand. The government’s administrative division follows a four-level structure—national, provincial (including municipalities), district, and county. Generally, a province includes the provincial city and several districts, which are administrated by the provincial capital. Each district has its own district capital city and its administrated counties.
Figure 1 illustrates the standard structure of the station location within a province.
In this paper, we classify stations into three levels according to their administrative division: provincial, district, and county level, and denote them as level and , respectively. The proportions of stations on the three levels are respectively denoted by satisfying . The provincial station set includes the stations in the provincial capitals (different HSR stations in one city can be viewed as a combined single station), which is denoted by , and the district and county station sets are denoted by respectively. As for a provincial station , the set of the district stations located in the corresponding province is denoted by ; similarly, the set of the county stations located in the corresponding district to station is denoted by .
2.2. O-D Demand Categorization
The above administrative structure influences the travel of demand and results in the attraction characters of provincial/district cities to travel demand. The attraction effects can be demonstrated in two aspects:
(1) High demand for travel across provinces or districts. Due to the high level of service for provincial stations, more trains would stop at these stations and a large number of passengers across provinces are generally attracted to get on/off at provincial stations to board direct trains to the destinations or find transfer travel schemes. Similarly, if there are a large number of county stations on the railway line, district stations would also attract a number of passengers across different districts to gather. Railway companies are required to provide sufficient train services to satisfy these direct travel requirements.
(2) Travel demand within a province or within a district. When a HSR line has been operated in a certain region such as a province or district, passengers traveling within that region often desire more direct train services to travel to/from the center city, i.e., the provincial or district city, for some political or economic reasons, even if the proportions may not be large enough. In order to better serve demand within a certain region, railway companies usually need to design the stop plan by taking provincial capitals and district cities as centers and improving the level of service for these passengers, i.e., increasing the service frequencies within a province or district to provide more convenience.
Therefore, stop planning should consider the attraction effects of stations on different levels and the travel characters of different O-D demand. The level of service for demand between the following O-D pairs should be mainly focused on: O-Ds between provincial stations and O-Ds traveling within a province or district. In this paper, we defined three movement types according to the O-D station location of passengers—(i) cross-province travel: passengers’ origin and destination stations belong to different provinces; (ii) same-district travel: the O-D stations belong to the same district, including travel to/from the provincial station; and (iii) cross-district travel: the O-D stations belong to different districts of the same province.
Based on the above O-D movement types and station levels, we classify the travel demand into 13 categories of bi-directional travel demand. The 13 demand categories, denoted as
, can be expressed by a spatial-level matrix shown in
Table 5, where
represent the levels of the O-D stations. The proportions of passenger volumes in the 13 demand categories are respectively denoted by
, which satisfies
and
. These proportions can be collected in the set
.
2.3. Train Types
As the demand to/from a county stations are relatively low, in practice, not all trains need to stop at county-level stations. In this paper, we consider two types of trains: trains of type 1 can stop at any station, while those of type 2 would not stop at county stations.
Then the 13 demand categories can be further grouped in terms of the train service availability. Eight of the 13 demand categories are relevant to county stations and can therefore only be served by train type 1. We term this group of demand as the “single-service” demand category, and denote the set of these eight categories by . The remainder five demand categories are grouped into set , which is termed as “double-service” demand category. The proportions of the two kinds of demand are denoted by .
Double-service demand have a choice of which train type to take. The proportions of passengers who only choose train type 2, denoted as and satisfying , can be obtained by observations and are considered as inputs to the model.
The proportions of the train kilometers of the two train types are denoted by satisfying . Since the single-service demand can only choose train type 1,
should be a positive value, and the two proportions need to satisfy . From the view of average, the proportions of the numbers of passing trains of type at any station can also be viewed as .
2.4. Train Stop Probability Optimization Problem
Generally, the service-oriented HSR line planning can be defined as follows.
Definition 1. With the given HSR network allocation limit, passenger travel demand, and a predefined level of service, an optimal line plan should minimize the total travel time of passengers.
If routes, frequencies and train types are provided, the only part to optimize in line planning is the stop plan. To solve this sub-problem of line planning, we propose the following assumptions for the train stop probability optimization:
(A1) In this paper, the train stop probability optimization only refers to a single HSR line.
(A2) As for any county or district station on the HSR line, its corresponding administrative district or provincial station is also on this line, which is illustrated in
Figure 2. This assumption may not be absolute in reality, but it can be satisfied by merging stations, or upgrading/degrading the station level by the dispatched passenger volumes at each station.
(A3) All trains running on the line have the same capacity which can satisfy the predefined level of service.
(A4) The passenger demand of all O-D pairs follows the uniform distribution pattern during the whole operation period .
(A5) The headway between any two adjacent trains is uniform.
(A6) In the planning stage, the ultimate purpose of stop planning is to provide sufficient direct services for different O-D pairs. Thus, the direct travel scheme is the only considered scheme in this paper.
(A7) Trains have independent stop probabilities in any two provinces, and the stop probabilities at provincial stations are independent and identically distributed at any two provincial stations. When the stop at a provincial station is determined, the stop probabilities at any two administrated district stations in the corresponding province are independent and identically distributed. Similarly, after the stop at a district station is determined, the stop probabilities at any two county stations in the corresponding district are independent and identically distributed.
Thus, the train stop probability optimization problem can be defined as follows:
Definition 2. Given the HSR line, routes and frequencies of trains, passenger demand, and assumptions (A1)–(A7), the stop planning can be described by the stop probabilities of the two train types at stations on different levels. The train stop probability optimization is to determine the train stop probabilities at stations on different levels and the reasonable proportion of the train kilometers of each type with the objective of minimizing per capita travel time. The outcome of the solved stop probabilities can be used to obtain the required average number of train services for O-D pairs in each demand category.
The stop-probability approach follows a “stop probability
service probability
service frequency” process. As illustrated in
Figure 3, given the input data, we adopt the stop probabilities to formulate stop planning with the objective of minimizing the per capita travel time, subject to constraints such as equilibrium of train load and stop density, and obtain the optimized stop probabilities of train types at different station levels. For an O-D pair, a train service should simultaneously stop at the origin station and destination station. The service probability for a demand category is deducted by the stop probabilities at the O-D stations that are determined by their station levels. Given the average number of passing trains of the demand category, we can further obtain the required average frequency for the O-D pairs in this category.
3. A Stop Probability Based Stop Planning Optimization
3.1. Probability of Train Services by O-D Pair Categories
Let denote the stop probabilities of train type at stations on level . We derive the probability of train services for each O-D category.
Due to the attraction effects and high levels of center stations within a region, the stop planning of stations in a district (or county) in a province (or district) is dependent on the train stopping at the corresponding provincial (district) station, and we adopt the conditional probabilities to describe the stop planning process between each two neighboring station levels. As for a station
on a lower level, which satisfies
, the stop probability at
can be deducted in two situations—if a train of type
stops at the corresponding administrative station
, the conditional stop probability at
is denoted by
; otherwise, the conditional stop probability at
is denoted by
. Since trains of type 2 do not stop at county stations,
Thus, the stop probability of trains of type
at a station on a lower level
can be
The average stop probability of trains of type
is
In the range of a province, after the stop at the provincial station is determined, the stops at the district/county stations in that province have specific conditional probabilities, which can be deducted by the above conditional probabilities. If a train of type
has stopped at the provincial station, the conditional probability of stopping at a station on level
in that province is denoted by
; otherwise, it is
.
Similarly, if a train of type has stopped at a district station, the conditional probability of stopping at a county station in that district is just ; otherwise, it is .
We denote as the probability of train type serving the demand category . Next, we present the detailed deduction by the travel movement types of each demand category.
(1) Cross-province O-D movements
When
, the demand categories correspond to cross-province O-D movements, i.e.,
. Since the stop probabilities at the O-D stations are independent of each other, the probability of serving O-D pairs in each category is the product of the stop probabilities at the O-D stations.
(2) Same-district O-D movements
When
, the demand categories correspond to same-district O-D movements, i.e.,
. The stop probabilities at the O-D stations are not independent, and the stop probability at a station on a lower level in the O-D pair is related to the stop at the other station, we should consider the corresponding conditional probabilities of stopping at a station on level
in the range of a province or district.
(3) Cross-district O-D movements
When
, the demand categories correspond to cross-district O-D movements, i.e.,
. The stop probabilities at the O-D stations are independent, but their stop probabilities both depend on the stop at the corresponding provincial station, so we should consider the corresponding conditional probabilities of stopping at a station on level
in the range of a province.
With given train routes and frequencies, we can obtain the average numbers of passing train services for the O-D pairs in each demand category, and denote them by
, which can be collected in the set
. For each demand category, the average numbers of train services of the two types serving each O-D pair is
times of these numbers. Using Equations (5)–(7), the average number of train services for the O-D pairs in
, can be obtained by the following formula, and these service frequencies can be collected in the set
3.2. Per Capita Dwell Time
Generally, more stop stations would result in longer passenger travel time, so the concept of per capita dwell time, i.e., the average passenger dwell time at intermediate stop stations, is a significant index of stop planning. We denote
as the per capita travel distance of
,
.
Figure 4 illustrates the route of an O-D pair in
including the O-D stations and intermediate passing stations. If the demand in this category choose trains of type
, the corresponding per capita dwell time can be denoted by
, which is the product of the number of intermediate stop stations and the average stop time at a station
. In the view of average, using the average station distance
, we can obtain the number of intermediate passing stations is
. Then, the average number of intermediate stop stations of this O-D pair is the product of the number of intermediate passing stations and the average stop probability at any station.
Using Equation (2), we can obtain the average stop probability of trains of type
at any station is
. Thus, the per capita dwell time of passengers in
choosing trains of type
is
Since the single-service demand can only choose train type 1, the average dwell time of the single-service demand
can be obtained if we replace
in Equation (9) by the per capita travel distance of the single-service demand
, i.e.,
In each double-service demand category , except the passengers who only choose train type 2 (accounting for ), the remaining demand can choose either train type, so we need to calculate the division proportions of the two train types for these passengers. We denote the division proportions of passengers choosing the two train types in each double-service demand category by , which satisfies .
From assumption (A4), passenger demand of all O-D pairs in each category follow the uniform distribution pattern during the whole operation period , so the division proportions of double-service demand choosing the two train types are equal to those of attraction time lengths of the two train types in the whole operation period .
Generally, passengers prefer the train type with less per capita dwell time, and a larger number of trains will also attract more passengers. Initial attraction time is the average attraction times of the two train types to each O-D pair in , which are determined by the per capita dwell time, and denoted by respectively.
If , then the attraction times of the two train types are both 0, i.e., .
If , then the attraction time radius of each train of type 2 to each O-D pair in is . In view of situations of departing earlier or later, the attraction time of a train of type 2 is . The average number of train services of type 2 for O-D pairs in is . Assumption (A5) indicates that trains run with a uniform headway, so we can obtain the initial attraction times , and .
If , then we can obtain the attraction times similarly: .
Except from the total attraction time, the surplus time length is . We can roughly divide the surplus time by the proportions of the train services of the two types for each O-D pair in . The surplus attraction times of the two train types can be shown as .
According to the initial attraction times and surplus attraction times of the two types of trains, we can obtain the division proportions of train type
satisfying
.
Combined with the division proportions, we can obtain the per capita dwell time of the double-service demand.
3.3. Per Capita Boarding Deviation Time
3.3.1. Per Capita Boarding Deviation Time of Single-Service Demand
Since single-service demand can only choose type 1, and the whole O-D demand follows the uniform distribution, we assume the average headway between the train services for the corresponding O-D demand is
. As shown in
Figure 5, the attraction time radius of a train of type 1 is
, and the boarding deviation time follows a linear distribution in the attraction time radius, so the per capita boarding deviation time of single-service demand is equal to the integral of the boarding deviation time of the demand departing during the attraction time, i.e.,
Thus, the per capita boarding deviation time of single-service demand is
3.3.2. Per Capita Boarding Deviation Time of Double-Service Demand
As for the passengers who only choose trains of type 2, we obtain the per capita boarding deviation times
As for the passengers who can choose either train type, we can determine the service time of the two train types
. Then we can obtain the per capita boarding deviation times of passengers who choose the two train types
Thus, the per capita boarding deviation time of double-service demand is
If the denominators in formulas (15) and (16) are 0, then the per capita boarding deviation times can be replaced by a sufficiently large constant .
3.4. Train Load and Stop Density Equilibrium
3.4.1. Equilibrium of Train Load
The equilibrium of train load is that the proportions of the train kilometers of the two train types are close to those of their transported passenger kilometers, so as to keep a balance between the loads of the two train types. The average transported passenger kilometers of the total demand by the two types are
and
, respectively. Thus, the proportion of the transported passenger kilometers by type 1 is
. The equilibrium of train load can be expressed by the relative error between
and
, i.e.,
in which
, and a larger line capacity corresponds to a larger
.
3.4.2. Equilibrium of Stop Density
If a train stops at a provincial/district station, then the stop probability in the corresponding province/district will reach maximum; otherwise, it will be minimal. We define the concept of stop density to describe the stop distribution in a province/district, which is the average stop number in that region. The imbalanced stop probabilities would result in the deviation among stop densities of trains at one province or district, in other words, the deviation among travel speeds (the travel time of a train includes the stop time at intermediate stations). With given headways at the starting station, the utilization ratio of the line capacity will be improved if trains have approximate travel speeds because in this way the headways between trains in each section are close to the defined headways at the starting station.
Next, we discuss how to ensure the equilibrium of stop density by stop probabilities.
(1) In a province
If a train stops at a provincial station, then the average stop probability in that province is ; otherwise, it is . The difference between them, i.e., , is equal to . If the stop probability of trains of type 1 at provincial stations is 1, then .
Generally, there is only one provincial station in a province, so the average number of stations in a province is
, and the difference between average numbers of stop stations is
, of which the upper limit is
.
is determined by the utilization ratio of the line capacity: as for a busy line,
is usually no larger than
. Thus, the difference between the average stop probabilities of trains of type 1 must satisfy
(2) In a district
If a train stops at a district station, then the average stop probability in that district is ; otherwise, it is . The difference is . If the stop probability of train of type 1 at district stations is 1, then .
If there is only one district station in a district, then the average number of stations in a district is
, and the difference between average numbers of stop stations is
, of which the upper limit is
, which is usually no larger than
on a busy line. Thus, the difference between the average stop probabilities of trains of type 1 must satisfy
(3) The constraint of the stop probabilities of trains of type 2
Since trains of type 2 do not stop at county stations, the stop probabilities of trains of type 2 are far less than those of type 1. In order to ensure reasonable stop probabilities of trains of type 2, we propose a constraint that the stop probabilities of trains of type 2 at provincial or district stations should be no less than those of type 1, i.e.,
3.5. Optimization Formulation and Solution Algorithm
3.5.1. Model Formulation
As introduced earlier, the stop planning problem can be described as to optimize stop probabilities to minimize passengers’ direct O-D travel time. The travel time of a passenger consists of the boarding deviation time, running time and stop time at intermediate stations. Since the train routes on a line are predefined, we do not consider the running times between stations, and the per capita travel time is equal to the weighted average of the travel times of single/double-service demand, of which the weights are
and
respectively. Thus, the optimization model of stop probability is as follows:
The objective function (22) is to minimize the per capita travel time; constraint (18) aims to ensure a balance between the train loads of the two types; constraints (19) and (20) are aimed to improve the utilization ratio of the line capacity; constraint (21) is to ensure the reasonable demand allocation on the two train types.
3.5.2. Solution Algorithm
Since the optimization model is a non-convexity problem, we adopt the penalty function method and introduce large penalty factors to add the constraints (19)–(21) to the objective function as penalty functions. The proportion of transported passenger kilometer by type 1, i.e., , is not a continuous variable, so we keep the formula (18) in the constraint to avoid shocks of the objective function in the solving process.
After eliminating the intermediate variables, the remaining variables include and ( is equal to ). We set to ensure the feasibility of and Equation (18). Due to the categorization of stations and O-D demand, the number of decision variables is not influenced by the structure and scale of the studied HSR line, and the solving scale of decision variables is fixed.
Given that all of the variables in
have continuous ranges, we randomly generate several initial solutions, search for the corresponding local optimal solutions from each initial solution, and select an optimal solution from these local optimal solutions. The detailed procedure is described below.
Algorithm 1 Procedure for solving stop probability optimization problem |
Step 0: Initialization |
Deduct the following parameters from the given data, i.e., train routes and frequencies, O-D demand, and the line operation data: |
(i) HSR line parameters: |
(ii) demand parameters: |
(iii) operation parameters: |
Step1: Random generations of solutions |
Discretize the range of by the step length 0.001 to generate different values of , and randomly generate an initial feasible solution from the uniform distribution for each . Calculate the objective with determined and according to the following procedure. |
Solve by Equation (9); |
Obtain by solving Equations (10)–(12); |
Obtain by solving Equations (13), (14), and by solving Equations. (15)–(17); |
Solve the objective by Equation (22). |
Step 2: Local search |
As for each group of , set the local optimal solution . |
Local search process |
Set the initial search step length . |
Whiledo |
Search for the neighbor solution using the Gradient Descent Method and calculate the objective . If and the solution satisfies the constraint (18), update the solution . |
End |
Step 3: Local search repetition |
As for the group of local optimal solution , update the values of in the solutions of by , and repeat Step 2 to search a better local optimal solution for . |
Step 4: Select the optimal solution |
Compare the objectives of the local optimal solutions and select the solution with minimal objective as the optimal solution, and obtain by Equations (1)–(8). |
5. Conclusions
In this paper, we propose an analytical method that provides the theoretical deduction of required service frequencies for passenger O-D demand as specified by geo-political station levels and passenger demand levels. Considering the attraction effects of provincial and district stations, we classify demand into several categories that can reduce the scale of the number of O-D pairs and also present the travel characters of each demand category. Given the level of service, trains are classified into two types by service availabilities. We propose a concept of stop probability to describe the stop planning problem and formulate the optimization of stop probabilities of HSR trains to obtain the average numbers of train services for O-D pairs in each demand category. The proposed model can determine the train stop probabilities at stations on different levels and the reasonable proportion of the train kilometers of each type with the objective of minimizing per capita travel time.
In the model formulation, we consider two equilibrium conditions—(i) equilibrium of train load to keep a balance between the train kilometers and transported passenger kilometers of the two train types and (ii) equilibrium of stop density to improve the utilization ratio of the line capacity.
The proposed method is validated on three real-life HSR lines. The solved stop probabilities are shown to match closely to those in practice. The solved stop probabilities of the two train types at provincial stations are close to 1, which complies with the attraction effect of provincial stations. The stop probability of train type 1 at county stations should be increased to serve more single-service demand, while that at district stations is decreased to lead more double-service demand to choose the other train type, which can reach a more reasonable load division between the two train types. The station set with a smaller proportion (the provincial station set) generally needs a higher stop probability, which can ensure more serving trains and reduce the added travel time.
The proportion of passenger kilometers of single-service demand is the lower bound of the proportion of transported passenger kilometers of train type 1, so it has an influence on the minimal feasible proportion of train kilometers of train type 1. We obtain a regression relationship between the proportions of train kilometers and passenger kilometers of single-service demand; the result has good practical significance.
Our future research can be extended in three aspects: the first is to optimize the service frequency considering more O-D demand characteristics, such as passengers’ time value, different distribution patterns during a day, the influence of O-D distances, and so on; the second is to adopt the solved stop probabilities to design a more reasonable stop plan for different O-D levels; and the third is to extend the studied structure to HSR networks considering characteristics of different train services, such as in-line and cross-line trains.