# Zone Fare System Design in a Rail Transit Line

Authors: Yi Yang; Lianbo Deng; Qing Wang; Wenliang Zhou
Journal: Journal of Advanced Transportation
 (2020)

---

## Abstract

As a widely existing form of public transit fare structures, zone fare system is traditionally designed in a separate way, which may lead to suboptimal results. This paper aims to concurrently address the zone division and fare calculation issues of the zone fare system design in a rail transit line. It is necessary to consider passengers as well as operators to find an impartial zone fare system. A fair zone fare system is one where the zone fares are as close as possible to the distance-based fares, for the fares of the distance-based fare system are highly correlated with the actual distance of trips. Thus, the fare deviations for the trips between the zone fare system and the distance-based fare system are utilized as the evaluation metric of fairness. To achieve the goal of minimizing fare deviations for all trips, we introduce three indexes: average absolute deviation, average squared deviation, and maximum mileage fare deviation. With the three indexes as the objective functions, we develop a joint optimization algorithm where a novel zone boundary adjustment scheme is proposed as the key technique. Numerical results show that the proposed algorithm can effectively provide a joint optimal scheme and the optimal number of planned zones is 5 for Changsha Metro Line 2. The proposed algorithm can provide guidance for the practical design and adjustment of the zone fare system.


---

## Body

## 1. Introduction

Fare is crucial to the financial situation of the rail transit. The price unit of tickets is directly determined by fare structure, which is the base and premise of the price formation. According to the relationship between the fare and the distance, fare structure generally includes the flat fare and the graduated fare.The flat fare assigns the same fare to all trips, which is the easiest and cheapest policy to implement. However, the flat fare is extremely unfair to the short-distance passengers, who pay the same price as the long-distance passengers. Generally, the public does not accept the flat fare. The graduated fare, by contrast, can be calculated respectively based on the distance and the number of crossed zones. Thus, the graduated fare is more equitable than the flat fare. Two types of the graduated fare are therefore commonly used: distance-based fare and zone fare.The most frequently used approach is the distance-based fare system, where the price only depends on the travelled distance. The longer the trip, the higher the fare will be. Therefore, the distance-based fare system is often considered as fair [1–3]. However, in order to calculate a specific fare, the distance for each pair of stations needs to be available for passengers and operators. This makes the system inconvenient, especially for the big city with a large number of stations [1, 2].There are two types of zone fare systems. One is the zone fare system with arbitrary prices, such as the system in the north of San Francisco [1], where the fare is given arbitrarily for each pair of zones. The other is the counting zone fare system. As the most popular type of zone fare system, it has been widely accepted. Examples of such a zone fare system can be found in Oslo [4] and in the south of San Francisco [1], where the fare of a trip depends only on the crossed zones on the trip. In other words, crossing the same number of zones, passengers must pay the same price, even if they have different starting zones and different ending zones. Because of these unique properties, zone fare system is easy to handle and can significantly facilitate passengers’ travel and operators’ management.In brief, the graduated fare is more acceptable to the public than the flat fare. In terms of the graduated fare, as a very popular system in Germany and Switzerland [5], the zone fare system is simpler than the distance-based fare system [1–3]. Over the past several decades, the distance-based fare system has been extensively studied [6–15], but the research on the design of zone fare system has been neglected. To the best of our knowledge, the literature on zone fare system is scarce.Several studies evaluated the application effect of zone fare system from the view of economics. Jansson and Angell [4] indicated that the ticket prices of zone fare system could approximately reflect the social marginal costs of the trips, because zones vary with size and location. Andrle et al. [16] showed that zone fare system could be introduced effectively in a small transit system. After the historically complex fare system of Haifa city was simplified to a five-zone fare system, Sharaby and Shiftan [17] found that the system could reduce the travel cost of passengers and attract more passengers. By using the analytic hierarchy process, Nassi and Costa [18] found that the zone fare system was the optimal mode compared to other fare systems. Moreover, with regard to exploiting the passengers’ willingness to pay, Otto and Boysen [2] showed that zone fare system is quite robust compared to the flat and distance-based fare systems.According to these studies, we find that zone fare system is an economical and effective system in public transit. Thus, it is necessary to explore how to design the zone fare system. However, the literature on zone fare system design is rare. The problem of designing zone fare system is NP-hard [2, 19], which includes two subproblems: zone division and fare calculation. In terms of the implemented zone fare system in Potsdam (Germany), Borndörfer et al. [20] proposed a nonlinear fare optimization model and obtained the optimal fares. Schöbel [21] and Hamacher and Schöbel [1] provided three heuristics to obtain a fair zone design scheme in linear station networks, and they successfully redesigned the fare structure in one of Germany’s transport systems by processing zone division and fare calculation individually. Koháni [22] introduced a mathematical model based on counting zones and developed two approaches to determine the quality of solution. Then, Koháni [3] made the case study for Žilina Municipality with that mathematical model. But the computation time varies from few seconds up to 2 hours for one particular parameter setting, which is not efficient enough. In order to minimize the maximum fairness deviation, Babel and Kellerer [19] proposed several heuristic algorithms and presented theoretical results for special networks. Conejero et al. [23] put forward a method of the zone fare system based on finding the longest path in a weighted root tree, and they considered all possible distributions in turn to obtain the optimal solution. Pratelli [5] developed a two-level model to deal with the problems of zone division and fare calculation separately.In summary, the previous studies mainly focused on evaluating the application effect of the zone fare system based on the cities that have implemented this fare system. In terms of the system design, the joint optimization of zone division and fare calculation has been largely ignored in the previous studies, which may undermine the benefit of the zone fare system. In light of this, we develop a novel optimization algorithm to simultaneously solve the problems of zone division and fare calculation in designing the zone fare system. We should ensure the fairness when designing a zone fare system. Thus, the zone fares of all trips should be as close as possible to the distance-based fares, for the distance-based fare system highly relates the price to the actual distance of a trip and is mostly considered as fair [1–3]. In view of this, the fare deviations between the zone fare system and the distance-based fare system can be adopted as the objective function of the algorithm [1, 19–21]. Therefore, three objective functions are introduced to evaluate the fairness: average absolute deviation, average squared deviation, and maximum mileage fare deviation. Then, with the zone boundary adjustment as the key technique, a joint optimization algorithm is proposed based on the zone division and the fare calculation. The paper not only provides a joint optimization algorithm to design the zone fare system but also a decision-making basis for the future of this system.The main contributions of the paper are threefold: (1) we propose an innovative algorithm to concurrently process the zone division and fare calculation for the zone fare system design; (2) we introduce three indexes to evaluate the fairness of the zone fare system scheme; (3) a novel zone boundary adjustment method is put forward as the key technique in the algorithm.The remainder of this paper is organized as follows. Section2 analyzes the two subproblems: zone division and fare calculation. To select the objective function, the evaluation indexes of the scheme are analyzed in Section 3. In Section 4, with zone boundary adjustment as the key technique, we present a novel heuristic algorithm; In Section 5, a numerical example is conducted on Changsha Metro Line 2. Finally, conclusions are drawn in Section 6.


## 2. Analysis

The problem can be divided into two subproblems: zone division and fare calculation. For the rail transit lines, their operating policies are independent of each other, and so are the fare policies. Therefore, the study scope can be limited to a rail transit line. Based on a line withn stations, the set of stations is denoted as S=s1,s2,…,sn.


### 2.1. Zone Division

Zones are divided along the line sequentially and fares are the same for both directions. In the zone fare system,K represents the number of planned zones, and the zone is expressed as zk, 1≤k≤K. zk denotes the number of stations within zk. Besides, each zone contains at least one station, namely, zk≥1. For any station si, it must belong to only one zone that can be denoted as zk. Then, zsi=zk, 1≤k≤K, and 1≤i≤n. zsi,sj is the number of zones for the trip from si to sj. Thus, the first and last stations on the line respectively belong to the zones z1 and zK, namely, zs1∈z1, zsn∈zK. To ensure there is at least one station in each zone, we set zk≥1, for 1≤k≤K.We definepij as the number of passengers from si to sj, 1≤i,j≤n,i≠j. As discussed above, the fare deviations between the zone fare system and the distance-based fare system are used to evaluate the fairness of zone fare system. For this reason, in the study period, the impact of the fare level on the number of passengers can be ignored. That is, pij is a given constant based on the real data without considering the demand elasticity. The fare of pij can be expressed as hij=hza,zb, where zsi=za, zsj=zb.


### 2.2. Fare Calculation

In zone fare system, the price of a trip is determined by the starting and ending zones of the trip. In general, for passengers with the same number of crossed zones, according to whether their fares are different, two methods of the ticket prices can be calculated as follows: the zones-differential pricing and the zones-uniform pricing. For passengers with the same zone fare, their fare levels can be determined in accordance with the mean of the maximum and minimum fares, the median price, the average fare, and other ways under the distance-based fare system. The distance-based fare is formulated as follows:(1)fij=fw0+fwwij,where fw0 and fw are the fixed fare and variable fare coefficients, respectively, and wij is the distance from station i to station j, 1≤i,j≤n,i≠j.


#### 2.2.1. The Zones-Differential Pricing

The differential pricing is determined by the starting and ending zones. The fares may be different for passengers crossing the same number of zones in the trips. For the zone fare system with arbitrary prices, utilizing the average distance-based fares, the zone fare from zoneza to zone zb can be computed by(2)hza,zb=∑si∈za∑sj∈zbfijpij∑si∈za∑sj∈zbpij.


#### 2.2.2. The Zones-Uniform Pricing

The uniform pricing is determined only by the number of crossed zones. For the counting zone fare system, the zone farehk is expressed as(3)hk=12maxzsi−zsj=kfw0+fwwij+minzsi−zsj=kfw0+fwwij,where zsi−zsj=k indicates the trips from station si to station sj with k crossed zones.In general, the zones-differential pricing method is used to calculate the zone fares for a given zone division scheme, and the zone fare level is determined according to the average distance-based fare in the following sections. Moreover, the following optimization methods can be easily extended to the zone-uniform pricing methods and other pricing means.


## 2.1. Zone Division

Zones are divided along the line sequentially and fares are the same for both directions. In the zone fare system,K represents the number of planned zones, and the zone is expressed as zk, 1≤k≤K. zk denotes the number of stations within zk. Besides, each zone contains at least one station, namely, zk≥1. For any station si, it must belong to only one zone that can be denoted as zk. Then, zsi=zk, 1≤k≤K, and 1≤i≤n. zsi,sj is the number of zones for the trip from si to sj. Thus, the first and last stations on the line respectively belong to the zones z1 and zK, namely, zs1∈z1, zsn∈zK. To ensure there is at least one station in each zone, we set zk≥1, for 1≤k≤K.We definepij as the number of passengers from si to sj, 1≤i,j≤n,i≠j. As discussed above, the fare deviations between the zone fare system and the distance-based fare system are used to evaluate the fairness of zone fare system. For this reason, in the study period, the impact of the fare level on the number of passengers can be ignored. That is, pij is a given constant based on the real data without considering the demand elasticity. The fare of pij can be expressed as hij=hza,zb, where zsi=za, zsj=zb.


## 2.2. Fare Calculation

In zone fare system, the price of a trip is determined by the starting and ending zones of the trip. In general, for passengers with the same number of crossed zones, according to whether their fares are different, two methods of the ticket prices can be calculated as follows: the zones-differential pricing and the zones-uniform pricing. For passengers with the same zone fare, their fare levels can be determined in accordance with the mean of the maximum and minimum fares, the median price, the average fare, and other ways under the distance-based fare system. The distance-based fare is formulated as follows:(1)fij=fw0+fwwij,where fw0 and fw are the fixed fare and variable fare coefficients, respectively, and wij is the distance from station i to station j, 1≤i,j≤n,i≠j.


### 2.2.1. The Zones-Differential Pricing

The differential pricing is determined by the starting and ending zones. The fares may be different for passengers crossing the same number of zones in the trips. For the zone fare system with arbitrary prices, utilizing the average distance-based fares, the zone fare from zoneza to zone zb can be computed by(2)hza,zb=∑si∈za∑sj∈zbfijpij∑si∈za∑sj∈zbpij.


### 2.2.2. The Zones-Uniform Pricing

The uniform pricing is determined only by the number of crossed zones. For the counting zone fare system, the zone farehk is expressed as(3)hk=12maxzsi−zsj=kfw0+fwwij+minzsi−zsj=kfw0+fwwij,where zsi−zsj=k indicates the trips from station si to station sj with k crossed zones.In general, the zones-differential pricing method is used to calculate the zone fares for a given zone division scheme, and the zone fare level is determined according to the average distance-based fare in the following sections. Moreover, the following optimization methods can be easily extended to the zone-uniform pricing methods and other pricing means.


## 2.2.1. The Zones-Differential Pricing

The differential pricing is determined by the starting and ending zones. The fares may be different for passengers crossing the same number of zones in the trips. For the zone fare system with arbitrary prices, utilizing the average distance-based fares, the zone fare from zoneza to zone zb can be computed by(2)hza,zb=∑si∈za∑sj∈zbfijpij∑si∈za∑sj∈zbpij.


## 2.2.2. The Zones-Uniform Pricing

The uniform pricing is determined only by the number of crossed zones. For the counting zone fare system, the zone farehk is expressed as(3)hk=12maxzsi−zsj=kfw0+fwwij+minzsi−zsj=kfw0+fwwij,where zsi−zsj=k indicates the trips from station si to station sj with k crossed zones.In general, the zones-differential pricing method is used to calculate the zone fares for a given zone division scheme, and the zone fare level is determined according to the average distance-based fare in the following sections. Moreover, the following optimization methods can be easily extended to the zone-uniform pricing methods and other pricing means.


## 3. Objective Function

As a simplified form of the distance-based fare system [1–3], the rationality of a zone fare system scheme can be reflected in the difference between the two fares. Inspired by the previous studies [1, 19–21], the evaluation indexes are shown below. P=∑si,sj∈Spij represents the total passenger flow. The scheme of zone fare system becomes better as the index becomes smaller and more stable.


### 3.1. Average Absolute Deviation (AAD)

The average absolute deviation can be calculated as follows:(4)b1=1P∑si,sj∈Spijhij−fw0+fwwij.


### 3.2. Average Squared Deviation (ASD)

In addition to the average absolute deviation, the average squared deviation can be used as another evaluation index. It is calculated as follows:(5)b2=1P∑si,sj∈Spijhij−fw0+fwwij2.


### 3.3. Maximum Mileage Fare Deviation

Passengers who travel fromza and zb pay the same price in zone fare system, while the fares may be different for the passengers in the distance-based fare system due to their different travel distances. The maximum value of these fare deviations is denoted as the maximum mileage fare deviation bmax−min. For the passengers travelling from za and zb, we calculate the maximum difference between the highest distance-based fares and the lowest distance-based fares. The maximum difference can be used to compute bmax−min:(6)bmax−min=12maxza,zb∈Zmaxi∈za,j∈zbfw0+fwwij−mini∈za,j∈zbfw0+fwwij.We takeb1, b2, and bmax−min as the objective function, respectively. When one of them reaches the minimum, the optimal solution is obtained.


## 3.1. Average Absolute Deviation (AAD)

The average absolute deviation can be calculated as follows:(4)b1=1P∑si,sj∈Spijhij−fw0+fwwij.


## 3.2. Average Squared Deviation (ASD)

In addition to the average absolute deviation, the average squared deviation can be used as another evaluation index. It is calculated as follows:(5)b2=1P∑si,sj∈Spijhij−fw0+fwwij2.


## 3.3. Maximum Mileage Fare Deviation

Passengers who travel fromza and zb pay the same price in zone fare system, while the fares may be different for the passengers in the distance-based fare system due to their different travel distances. The maximum value of these fare deviations is denoted as the maximum mileage fare deviation bmax−min. For the passengers travelling from za and zb, we calculate the maximum difference between the highest distance-based fares and the lowest distance-based fares. The maximum difference can be used to compute bmax−min:(6)bmax−min=12maxza,zb∈Zmaxi∈za,j∈zbfw0+fwwij−mini∈za,j∈zbfw0+fwwij.We takeb1, b2, and bmax−min as the objective function, respectively. When one of them reaches the minimum, the optimal solution is obtained.


## 4. Algorithm

According to the characteristics of the zone fare system, we analyzed at the beginning, a joint optimization algorithm is designed to deal with the problems of zone division and fare calculation.For a zone division scheme, zone farehza,zb (za,zb∈Z) is calculated with equation (2).When moving the zone boundary, the number of crossed zones changes. Because the affiliation’s change of a boundary station between two adjacent zones will change the zone division scheme, the improved solutions can be achieved by adjusting the zone boundary. The most suitable objective functionb for a specific scenario will be selected as the goal when adjusting zone boundary.In the left-moving operation of the boundary, passengers who have a change in the number of their crossed zones are shown in the shadow part of Figure1(a). The boundary of zk and zk+1 is si,si+1. When the boundary moves to si−1,si along the left side of the line, station si changes from zsi=k to zsi=k+1. Thus, the number of crossed zones will be decreased by 1 for the passengers travelling from si to si+1,sn, and the number of crossed zones will be increased by 1 for the passengers who travel from s1,si−1 to si. Similarly, the right-moving operation that the boundary si−1,si moves to si,si+1 along the right side is shown in Figure 1(b).Figure 1
The range of the affected passengers in the moving operation. (a) The range of the affected passengers in the left-moving operation. (b) The range of the affected passengers in the right-moving operation.
(a)(b)In the moving operation, we only consider the passengers whose number of crossed zones changes. Then, we calculate the change of the objective function value under the current level of zone fares. A boundary adjustment scheme which mostly improves the objective function can be selected as a new scheme of zone division. The fares are updated according to the new zone division scheme. An iteration is finished. In this way, the optimal scheme is finally formed. In the rail transit line, we can get the initial fare scheme when the initial zone division scheme is obtained based on the passenger-kilometers (PKM) and passenger flow (PF). Thus, the initial scheme of the zone fare system is obtained.To illustrate the algorithm more explicitly, the algorithm where the initial solution is generated based on PKM is designed as follows:Input: station (S), distance (w), passenger demand (P), and the number of planned zones (K).Output: zonezsi and zone farehza,zb.(1)
Calculate the PKM from the starting station to the second-to-last station on the line. That is,Tsi,si+1=∑i′=1i∑j′=i+1nwi,i+1pi′j′+pj′i′.The zones are divided by(7)zsi=∑i′=1i−1Tsi′,si′+1∑i′=1n−1Tsi′,si′+1/K+1,1≤i≤n−1,K,i=n.Then, form the initial zone division schemeZ=zsi|si∈S.(2)
Calculate and update the zone fareshzk,zk′=∑si∈zk∑sj∈zk'fijpij/∑si∈zk∑sj∈zk'pij, 1≤k≤K.(3)
Setk=1, get the current scheme Z, and set the best scheme so far Z¯=Z.(4)
Explore the neighborhood of current scheme by moving left and right. Ifzsi=k, zsi+1=k+1, let si,si+1 be the boundary point between zone zsi and zone zsi+1. All boundary points move according to the following rules:(4.1)
Left-moving operation of the boundary point is(8)z−si,si+1=z¯si=zsi+1,z¯si+1=zsi+1.In this operation,z¯si and z¯si+1 represent the zones that si and si+1 belong to, respectively. For the adjusted zone-division scheme after the left-moving operation, the number of crossed zones from si' and sj' is(9)z′sj′,si′=zsj′,si′+1,i′=i;j′=1,…,i−1,zsj′,si′−1,i′=i;j′=i+1,…,n,zsj′,si′,i′≠i.After the left-moving operation, compared withZ, the change of the objective function is calculated by the formula Δbz′sj,si=bz′sj,si−bZ.(4.2)
Right-moving operation of the boundary point is(10)z+si,si+1=z¯si=z¯si,z¯si+1=z¯si+1−1.After the right-moving operation, the number of crossed zones fromsi' and sj' is(11)z″sj′,si'=zsj′,si′−1,i′=i+1;j′=1,…,i,zsj′,si′+1,i′=i+1;j′=i+2,…,n,zsj′,si′,i′≠i+1.After the right-moving operation, compared withZ, the change of the objective function is calculated by the formula Δbz″sj,si=bz″sj,si−bZ.(5)
Update the solution.(5.1)
WhenΔbz″sj,si≤Δbz′sj,si, if Δbz″sj,si<0 and z¯si+1>1, update z+si,si+1 inZ¯; otherwise, if Δbz′sj,si<0 and z¯si+1=1, z¯si>1, update z−si,si+1 in Z¯.(5.2)
WhenΔbz′sj,si≤Δbz″sj,si, if Δbz′sj,si<0 and z¯si>1, update z−si,si+1 in Z¯; otherwise, if Δbz″sj,si<0 and z¯si=1, z¯si+1>1, update z+si,si+1 in Z¯.(6)
Ifk<K, then k=k+1, and go back to step (4). Otherwise, if the scheme Z=Z¯, then the algorithm terminates and output the optimal results; if Z≠Z¯, then set Z=Z¯, and go back to step (2).The advantage of the joint optimization algorithm is that it can fill the gap of the previous studies where the zone division and the fare calculation are usually solved, respectively. Using the proposed algorithm, we can jointly address the two problems and achieve an optimal solution.The complexity of the algorithm depends on the time required in fare calculation and the zone boundary adjustment time. The complexities of constructing the initial schemes of zone division and fare calculation areon and on2, respectively. There are K−1 boundary points in this algorithm, and the total number of the left-moving and right-moving operations is 2K−1. The complexity of comparing different schemes is nK−1. The complexity of updating the new ticket prices is on2. Therefore, the sum of the algorithm complexity is ολKn3, where λ represents the number of iterations when the algorithm finally converges. The algorithm efficiency can be calculated according to λ.


## 5. Case Study

### 5.1. Background

As shown in Figure2, there are 19 stations in the Changsha Metro Line 2, and the length of the line is ws1,s19 = 21.36 km. The number in the gap between each pair of stations indicates the station distance in meters. In the current distance-based fare system, the fare can be computed by fi,j=2+0.24wi,j. The fare unit is yuan (¥) in China.Figure 2
Changsha Metro Line 2.


### 5.2. Scenario of 5 Planned Zones under the Objective Functionb2

To clearly illustrate the algorithm, we first study the scenario where the goal is to design 5 zones in a zone fare system with arbitrary prices. The initial solution is constructed based on the passenger-kilometers (PKM), and we takeb2 as the objective function, i.e., b=b2. In the optimal scheme, the value of the objective function is b2=0.15. Figure 3 shows the optimal zone division, where stations 1–4, 5–9, 10–12, 13–15, and 16–19 belong to zones 1, 2, 3, 4, and 5, respectively. The corresponding optimal fare scheme is shown in Table 1.Figure 3
Optimal zone division.Table 1
Optimal zone fares.
Zone1234512.573.103.814.635.9523.102.312.673.474.8033.802.672.292.804.1144.633.482.802.533.2955.954.804.113.292.60The indexes’ values of the optimal scheme are shown in Table2. The fare range is 2.15–7.13 yuan in the distance-based fare system, while it is 2.29–5.95 yuan in the zone fare system. Besides, the values of AAD and ASD are very small, which are respectively 0.31 yuan and 0.15 yuan. The indexes indicate that the zone fares are fair close to the corresponding distance-based fares for most trips. Thus, we can conclude that this scheme of the zone fare system is feasible. Moreover, figure 4 shows that the objective function can converge rapidly as the number of iterations increases. Therefore, the algorithm is proved to be convergent and efficient.Table 2
Evaluation of the optimal zone fare system scheme.
Evaluation indexAverage fareAverage absolute deviation (AAD)Average squared deviation (ASD)Highest fareLowest fareValue3.520.310.155.952.29Figure 4
Convergences of the initial schemes under PKM and PF. (a) The value ofb2 under PKM. (b) The value of b2 under PF.
(a)(b)To explore the application of the proposed algorithm for different initial schemes, we construct an initial zone division scheme based on the passenger flow (PF). By takingb2 as the objective function, it can be found that the optimal zone division scheme obtained from the initial zone division scheme under PF is the same as that under PKM. Figure 4 indicates that the algorithm efficiency under PF is similar to that under PKM, and the optimal objective function values are almost equal. Therefore, it can be concluded that the algorithm is applicable for different initial schemes.


### 5.3. Analysis underK Planned Zones

In order to further investigate the practical application, the number of planned zones is expanded toK∈3,13. Figure 5 shows the optimal zone division schemes obtained by using the three objective functions. Figure 5(a) depicts that the optimal schemes obtained by taking b1 and b2 as the objective functions are exactly the same. Besides, as the number of planned zones K changes from 3 to 13, the corresponding optimal zone division schemes under bmax−min are presented in Figure 5(b). Generally, Figure 5 indicates that the method is available for various numbers of planned zones.Figure 5
Optimal zone division schemes underb1, b2, and bmax−min for K∈3,13. (a) Optimal zone division schemes under b1 and b2. (b) Optimal zone division schemes under bmax−min.
(a)(b)The main results under different numbers of planned zones are analyzed as follows.


#### 5.3.1. Range of the Fares in the Optimal Solutions

Figure6(a) presents the highest and lowest zone fares obtained by respectively taking b1 and b2 as the objective functions. The highest and lowest distance-based fares are 7.13 yuan and 2.15 yuan. Figure 6(a) indicates that the curves corresponding to b1 and b2 are exactly same in the zone fare system for K∈3,13, where K=5 is an obvious inflection point. Figure 6(b) shows the curves of the highest and lowest fares obtained by taking bmax−min as the objective function, where K=5 is also a significant inflection point. In general, the zone fares are both gradually approaching to the distance-based fares as K increases, which conforms to an expected result of the algorithm that fare deviations between the zone fares and the distance-based fares should be small.Figure 6
Analysis under objective functionsb1, b2, and bmax−min for K∈3,13. (a) Fare ranges under b1 and b2. (b) Fare range under bmax−min. (c) Objective functions. (d) AAD and ASD under b1 and b2. (e) AAD and ASD under bmax−min.
(a)(b)(c)(d)(e)


#### 5.3.2. Values of the Objective Functions

As shown in Figure6(c), the three curves denote the values of b1, b2, and bmax−min for K∈3,13. The trend of curve b1 is similar to that of curve b2, and both of them are smoother than that of bmax−min. All the three curves have a downward trend from fast to slow, and the inflection points of them are K=5.


#### 5.3.3. Average Absolute Deviation (AAD) and Average Squared Deviation (ASD)

In Figure6(d), the blue and red curves show the values of AAD and ASD, respectively; these values are obtained by taking b1 and b2 as the objective functions for K∈3,13. Similarly, Figure 6(e) presents the curves of AAD and ASD when the objective function is bmax−min. Figures 6(d) and 6(e) show that, as the number of planned zones increases, AAD and ASD corresponding to the three objective functions gradually decrease. Besides, the curves in Figures 6(d) and 6(e) have the downward trend from fast to slow, and the spacing between the two curves becomes obvious when K=5. Thus, the optimal number of planned zones is 5. Figures 6(d) and 6(e) show that both b1 and b2 share the same value of AAD when K=5, and the value is slightly smaller than that under bmax−min. Therefore, compared with bmax−min, b1 and b2 are speculated to be more suitable to be the objective function in practice.In summary, Figures6(a)–6(e) show that there is a significant inflection point when in this scenario K=5. In terms of fares concerned by operators and passengers, compared with the distance-based fares, the highest and lowest zone fares are decreased by 1.18 yuan and increased by 0.14 yuan, respectively, and the average absolute deviation is 0.31 yuan. Therefore, there will be no obvious change in the interests of passengers and operators by the transition from the current distance-based fare system to a zone fare one with 5 planned zones, and it can be speculated that this zone fare system scheme is acceptable to operators and passengers. Meanwhile, as mentioned above, zone fare system has the advantage of greatly facilitating the travel of passengers and the management of operators. Therefore, we recommend that the optimal number of planned zones is 5 for Changsha Metro Line 2 if the operator plans to implement a zone fare system, and the optimal results of zone division and fare calculation for 5 planned zones are presented in Section 5.2.


## 5.1. Background

As shown in Figure2, there are 19 stations in the Changsha Metro Line 2, and the length of the line is ws1,s19 = 21.36 km. The number in the gap between each pair of stations indicates the station distance in meters. In the current distance-based fare system, the fare can be computed by fi,j=2+0.24wi,j. The fare unit is yuan (¥) in China.Figure 2
Changsha Metro Line 2.


## 5.2. Scenario of 5 Planned Zones under the Objective Functionb2

To clearly illustrate the algorithm, we first study the scenario where the goal is to design 5 zones in a zone fare system with arbitrary prices. The initial solution is constructed based on the passenger-kilometers (PKM), and we takeb2 as the objective function, i.e., b=b2. In the optimal scheme, the value of the objective function is b2=0.15. Figure 3 shows the optimal zone division, where stations 1–4, 5–9, 10–12, 13–15, and 16–19 belong to zones 1, 2, 3, 4, and 5, respectively. The corresponding optimal fare scheme is shown in Table 1.Figure 3
Optimal zone division.Table 1
Optimal zone fares.
Zone1234512.573.103.814.635.9523.102.312.673.474.8033.802.672.292.804.1144.633.482.802.533.2955.954.804.113.292.60The indexes’ values of the optimal scheme are shown in Table2. The fare range is 2.15–7.13 yuan in the distance-based fare system, while it is 2.29–5.95 yuan in the zone fare system. Besides, the values of AAD and ASD are very small, which are respectively 0.31 yuan and 0.15 yuan. The indexes indicate that the zone fares are fair close to the corresponding distance-based fares for most trips. Thus, we can conclude that this scheme of the zone fare system is feasible. Moreover, figure 4 shows that the objective function can converge rapidly as the number of iterations increases. Therefore, the algorithm is proved to be convergent and efficient.Table 2
Evaluation of the optimal zone fare system scheme.
Evaluation indexAverage fareAverage absolute deviation (AAD)Average squared deviation (ASD)Highest fareLowest fareValue3.520.310.155.952.29Figure 4
Convergences of the initial schemes under PKM and PF. (a) The value ofb2 under PKM. (b) The value of b2 under PF.
(a)(b)To explore the application of the proposed algorithm for different initial schemes, we construct an initial zone division scheme based on the passenger flow (PF). By takingb2 as the objective function, it can be found that the optimal zone division scheme obtained from the initial zone division scheme under PF is the same as that under PKM. Figure 4 indicates that the algorithm efficiency under PF is similar to that under PKM, and the optimal objective function values are almost equal. Therefore, it can be concluded that the algorithm is applicable for different initial schemes.


## 5.3. Analysis underK Planned Zones

In order to further investigate the practical application, the number of planned zones is expanded toK∈3,13. Figure 5 shows the optimal zone division schemes obtained by using the three objective functions. Figure 5(a) depicts that the optimal schemes obtained by taking b1 and b2 as the objective functions are exactly the same. Besides, as the number of planned zones K changes from 3 to 13, the corresponding optimal zone division schemes under bmax−min are presented in Figure 5(b). Generally, Figure 5 indicates that the method is available for various numbers of planned zones.Figure 5
Optimal zone division schemes underb1, b2, and bmax−min for K∈3,13. (a) Optimal zone division schemes under b1 and b2. (b) Optimal zone division schemes under bmax−min.
(a)(b)The main results under different numbers of planned zones are analyzed as follows.


### 5.3.1. Range of the Fares in the Optimal Solutions

Figure6(a) presents the highest and lowest zone fares obtained by respectively taking b1 and b2 as the objective functions. The highest and lowest distance-based fares are 7.13 yuan and 2.15 yuan. Figure 6(a) indicates that the curves corresponding to b1 and b2 are exactly same in the zone fare system for K∈3,13, where K=5 is an obvious inflection point. Figure 6(b) shows the curves of the highest and lowest fares obtained by taking bmax−min as the objective function, where K=5 is also a significant inflection point. In general, the zone fares are both gradually approaching to the distance-based fares as K increases, which conforms to an expected result of the algorithm that fare deviations between the zone fares and the distance-based fares should be small.Figure 6
Analysis under objective functionsb1, b2, and bmax−min for K∈3,13. (a) Fare ranges under b1 and b2. (b) Fare range under bmax−min. (c) Objective functions. (d) AAD and ASD under b1 and b2. (e) AAD and ASD under bmax−min.
(a)(b)(c)(d)(e)


### 5.3.2. Values of the Objective Functions

As shown in Figure6(c), the three curves denote the values of b1, b2, and bmax−min for K∈3,13. The trend of curve b1 is similar to that of curve b2, and both of them are smoother than that of bmax−min. All the three curves have a downward trend from fast to slow, and the inflection points of them are K=5.


### 5.3.3. Average Absolute Deviation (AAD) and Average Squared Deviation (ASD)

In Figure6(d), the blue and red curves show the values of AAD and ASD, respectively; these values are obtained by taking b1 and b2 as the objective functions for K∈3,13. Similarly, Figure 6(e) presents the curves of AAD and ASD when the objective function is bmax−min. Figures 6(d) and 6(e) show that, as the number of planned zones increases, AAD and ASD corresponding to the three objective functions gradually decrease. Besides, the curves in Figures 6(d) and 6(e) have the downward trend from fast to slow, and the spacing between the two curves becomes obvious when K=5. Thus, the optimal number of planned zones is 5. Figures 6(d) and 6(e) show that both b1 and b2 share the same value of AAD when K=5, and the value is slightly smaller than that under bmax−min. Therefore, compared with bmax−min, b1 and b2 are speculated to be more suitable to be the objective function in practice.In summary, Figures6(a)–6(e) show that there is a significant inflection point when in this scenario K=5. In terms of fares concerned by operators and passengers, compared with the distance-based fares, the highest and lowest zone fares are decreased by 1.18 yuan and increased by 0.14 yuan, respectively, and the average absolute deviation is 0.31 yuan. Therefore, there will be no obvious change in the interests of passengers and operators by the transition from the current distance-based fare system to a zone fare one with 5 planned zones, and it can be speculated that this zone fare system scheme is acceptable to operators and passengers. Meanwhile, as mentioned above, zone fare system has the advantage of greatly facilitating the travel of passengers and the management of operators. Therefore, we recommend that the optimal number of planned zones is 5 for Changsha Metro Line 2 if the operator plans to implement a zone fare system, and the optimal results of zone division and fare calculation for 5 planned zones are presented in Section 5.2.


## 5.3.1. Range of the Fares in the Optimal Solutions

Figure6(a) presents the highest and lowest zone fares obtained by respectively taking b1 and b2 as the objective functions. The highest and lowest distance-based fares are 7.13 yuan and 2.15 yuan. Figure 6(a) indicates that the curves corresponding to b1 and b2 are exactly same in the zone fare system for K∈3,13, where K=5 is an obvious inflection point. Figure 6(b) shows the curves of the highest and lowest fares obtained by taking bmax−min as the objective function, where K=5 is also a significant inflection point. In general, the zone fares are both gradually approaching to the distance-based fares as K increases, which conforms to an expected result of the algorithm that fare deviations between the zone fares and the distance-based fares should be small.Figure 6
Analysis under objective functionsb1, b2, and bmax−min for K∈3,13. (a) Fare ranges under b1 and b2. (b) Fare range under bmax−min. (c) Objective functions. (d) AAD and ASD under b1 and b2. (e) AAD and ASD under bmax−min.
(a)(b)(c)(d)(e)


## 5.3.2. Values of the Objective Functions

As shown in Figure6(c), the three curves denote the values of b1, b2, and bmax−min for K∈3,13. The trend of curve b1 is similar to that of curve b2, and both of them are smoother than that of bmax−min. All the three curves have a downward trend from fast to slow, and the inflection points of them are K=5.


## 5.3.3. Average Absolute Deviation (AAD) and Average Squared Deviation (ASD)

In Figure6(d), the blue and red curves show the values of AAD and ASD, respectively; these values are obtained by taking b1 and b2 as the objective functions for K∈3,13. Similarly, Figure 6(e) presents the curves of AAD and ASD when the objective function is bmax−min. Figures 6(d) and 6(e) show that, as the number of planned zones increases, AAD and ASD corresponding to the three objective functions gradually decrease. Besides, the curves in Figures 6(d) and 6(e) have the downward trend from fast to slow, and the spacing between the two curves becomes obvious when K=5. Thus, the optimal number of planned zones is 5. Figures 6(d) and 6(e) show that both b1 and b2 share the same value of AAD when K=5, and the value is slightly smaller than that under bmax−min. Therefore, compared with bmax−min, b1 and b2 are speculated to be more suitable to be the objective function in practice.In summary, Figures6(a)–6(e) show that there is a significant inflection point when in this scenario K=5. In terms of fares concerned by operators and passengers, compared with the distance-based fares, the highest and lowest zone fares are decreased by 1.18 yuan and increased by 0.14 yuan, respectively, and the average absolute deviation is 0.31 yuan. Therefore, there will be no obvious change in the interests of passengers and operators by the transition from the current distance-based fare system to a zone fare one with 5 planned zones, and it can be speculated that this zone fare system scheme is acceptable to operators and passengers. Meanwhile, as mentioned above, zone fare system has the advantage of greatly facilitating the travel of passengers and the management of operators. Therefore, we recommend that the optimal number of planned zones is 5 for Changsha Metro Line 2 if the operator plans to implement a zone fare system, and the optimal results of zone division and fare calculation for 5 planned zones are presented in Section 5.2.


## 6. Conclusions

In the paper, we study the problems of zone fare system design in a rail transit line. Considering the interaction between zone division and fare calculation, a joint optimization algorithm is proposed. In order to test the algorithm, the paper conducts a case study on Changsha Metro Line 2 based on real data. The conclusions are drawn as follows:(1)
The optimal number of planned zones: through the theoretical and empirical research, the analysis of the solutions under the three objective functions shows that the optimal number of planned zones is 5 for Changsha Metro Line 2.(2)
Analysis of the three objective functions: the experimental results indicate that the solutions obtained by takingb1 and b2 as the objective functions are quite similar, and the values of AAD under b1 and b2 are both slightly smaller than that under bmax−min. In general, b1 and b2 are more suitable to be adopted as the objective function than bmax−min in practice.(3)
The zone fare system we designed is acceptable and the algorithm is convergent and applicable. Comparing with the fares in the current fare system of the Changsha Metro Line 2, AAD and ASD indicate that there is no significant change in fares for most trips in the zone fare system. Thus, the algorithm generates an impartial zone fare system. Then, with regard to the property of simplicity as mentioned above, both operators and passengers will tend to accept the reform when the current distance-based fare system is transformed to the zone fare system. Therefore, the zone fare system we designed is acceptable. The algorithm is proved to be convergent and the algorithm complexity isολKn3. Moreover, as shown in Sections 5.2 and 5.3, respectively, the algorithm proffered can not only be applied to different initial schemes under PKM and PF but also is available for various numbers of planned zones. Thus, the algorithm has a fairly good applicability for different contexts. Besides, it is easy to implement with a short computing time.To the best of our knowledge, it is the first time that the zone division and fare calculation issues are concurrently addressed with a zone boundary adjustment as the key technique. It is worth noting that the choice of the fare system is not constant, which depends on the specific scenario of the rail transit systems and is not the focus of this study. In the future, we will explore the factors that affect the fare system and the corresponding application environment. In terms of the practical application, the study can provide a useful reference for future studies and cities that plan to implement the zone fare system.



---

*Source: 2470579-2020-12-04.xml*