Abstract

Unsuitable locations of bus stops which provide feeder transportation connecting subways near urban intersections usually lead to the low efficiency of public transport and level of passenger service. A multiobjective optimization model to distribute such stop locations is proposed to attain the shortest total walk distance of passengers and minimum delay time of cars through intersections and travel time of buses. The Pareto frontier and optimal solutions for the proposed model are given by the distance-based and enumerative methods. The Xizhimen bus stop is selected to implement case studies for verifying the validity of the proposed model. The analysis of sensitivity on possible solutions is also carried out in the case studies. The results show that the proposed model is capable of optimizing the locations of bus stops connecting subways near intersections and helpful to improve the level of passengers service and operational efficiency of public transportation.

1. Introduction

The demand of urban traffic in China has been dramatically increasing with the rapid development of economy, urbanization, and mobilization for last decades. Moreover, severe traffic congestion is challenging to the operation and management of urban traffic. More people have realized that urban public transit is an effective way to mitigate traffic congestion and improve the sustainability of transport development.

Urban public transit is a highly complex system where urban railways and buses are crucial parts of the sustainably social development. Due to reliable running time, efficient transport service, and energy consumption, urban rail transit which can carry a great number of passengers has been playing a very important role as the skeleton of urban traffic system. The number of passengers served by the Beijing subway system has exceeded 10 million persons per day. The efficiency of urban rail transit depends not only on its operational performance but also on good cooperation with other transportation modes, especially with those buses which are able to provide feeder service for passengers arriving at or leaving urban rail stations.

Bus stops, as the connecting points between buses and urban rail transit services, are important spaces for passengers to alight or board [1]. According to the real situation of urban roads, the areas around intersections are suitable and convenient to arrival and leaving of passengers. Thus, the bus stops near intersections are often considerable for passengers transferring between urban rail transit and buses, which have a great impact on level of passenger service and attraction of public transit.

The bus stops connecting subways have some different characteristics, such as more bus routes and berths, longer bus stopping time for passengers to alight or board, larger ridership amount, and more complex passenger flow than common bus stops. If bus stops are not located correctly, they are likely to cause serious impacts on road traffic flow [2] and intersection capacity [3] and even become the bottlenecks of urban traffic systems [4].

Three major criteria, the efficiency, effectiveness, and equity of service, are employed to evaluate public systems [5]. The efficiency of service is defined as the ratio of LOS (level of service) to the cost of resource consumed and thus usually measured by its cost. The effectiveness of service shows what comfortable service public transport can provide. The equity of service requires an indiscrimination of transport supplies that different users can obtain. The above three criteria reflect the interests of different stakeholders involved in operators (bus companies), bus users (passengers), and other users (private cars, taxis), respectively, which are imported into this study.

2. Literature Review

The location of feeder bus stops near intersections is a hotspot in research of public transit which has been drawing the attention of researchers for decades. Those achievements have an important impact on the quality of service provided for passengers and the operation of bus vehicles.

There are two major types of studies relating tothe location optimization of bus stops. The first type investigates the layout of bus stops to cut down the cost of buses or overall transportation system. For example, by taking possible changes in demand into account, a bilevel optimization model for locating bus stops was developed to minimize the social cost of the overall transport system [6]. Hafezi and Ismail [7] analyzed the efficiency of bus operation where the studies on bus stop location at three points including near-side, far-side, and mid-block were considered. Li [8] established a model that lists the stops in the order that they appear after organizing bus stops by routes and then transformed the problem of matching bus stops into a shortest path problem. Similar research was also found in the work of Hu et al. [9] which attempted to analyze the characteristics of urban rail transit and conventional buses.

The second type of studies are mainly focused on the impact on the capacity of roads or intersections as well as delay incurred by bus stops. For instance, a weighted-least-squares regression model was used to estimate the dwell time of buses at stops that associated prediction interval to reduce the negative impacts of nearside bus stops [10]. Furth and SanClemente [11] worked on the impact of bus stop locations on bus delay where buses are running on the near-side, far-side, uphill, downhill of road. By dividing bus delay into service delay and nonservice delay, Xu et al. [12] proposed a delay estimation model for buses at a bus bay stop and a curbside bus stop. Lu et al. [13] concentrated on the delay of buses near a stop when mixed traffic flow was considered by a special cellular automation model. The relationship between delay time and distance from a bus stop to a stop line of intersection, arriving rates, and dwell time of buses, signal circle, was studied [14]. A kinematic wave theory-based model was used to determine where to place a near-side stop to achieve a target level of residual car queuing [15]. Chen et al. [16] developed a computation method on the bus delay at stops in Beijing through statistical analysis.

The existing efforts that considered the benefits from any one perspective of passengers [17], buses [18], and cars [19], and sometimes even the combination of any two, such as the cost of buses and delay time of cars [20], have provided important references for optimizing locations of bus stops. However, there is insufficient research that fully accounts for the influences from some vital factors at the same time, such as the walking distance of passengers, delay time of cars, and travelling time of buses. It is necessary and helpful to make up the gap by implementing a multiobjective analysis.

The rest of this study is organized as follows. A multiobjective optimization model is developed in Section 3 to aim at the shortest total walking distance of passengers, minimum delay time of cars through intersections, and least travelling time of buses. The solution method is given in Section 4. The empirical studies for the Xizhimen bus stops in Beijing and their sensitivity analyses are finally conducted to ensure that the proposed model is effective in Section 5.

3. Model Formulation

3.1. Objective Functions

The locations selection of bus stops mainly influences the walking distance of passengers and travelling speed of cars and buses through impact regions. The optimization of total walking distance of passengers, delay time of cars through intersections, and travelling time of buses plays a crucial role in determining an appropriate position for bus stops.

3.1.1. Total Walking Distance of Passengers

Different from common passengers, the walking distance of passengers who want to transfer at the feeder bus stops connecting subways is usually classified into two parts: transfer distance from subways to buses and direct distance to take buses. Thus, the total walking distance of transferring passengers can be expressed by where is the total walking distance of passengers transferring from adjacent subways to buses, is the total walking distance which passengers take buses through upstream or downstream intersections, is the number of passengers who transfer between subways and buses, is the number of passengers who come from the upstream intersection adjacent to bus stops, is the number of passengers who come from the downstream intersection adjacent to bus stops, is the distance from bus stops to subway entrances where the abscissa of bus stop center is set to , is the distance from bus stops to their adjacent upstream intersection, and is the distance from bus stops to their adjacent downstream intersection. and are the abscissas of subway entrance and its adjacent downstream intersection, respectively. The upstream intersection is selected as the origin of horizontal ordinate. The detailed layout including the relative positions of bus stops, adjacent intersections, and subway stations is represented in Figure 1.

3.1.2. Delay Time of Cars through Intersections

It will be convenient for passengers to transfer between subways and other transportation modes and gather or disperse in all directions if the locations of bus stops are close to intersections, which is however likely to generate the delay time of private cars and taxis.

The following factors are key parameters in modeling delay at intersections with a bus stop: is distance between bus stop and stop line (m), is traffic flow including vehicles and buses at the key lane (vehicles/s), is the average arrival rate of every bus line at the bus stop (vehicle/s), is the average boarding and alighting time for passenger at the multiple feeder bus stop (s), is proportion of time the bus is blocked at the bus stop, which is equal to the value that multiplies by , is effective green times (s), and is cycle times (s).

The delay of cars at signalized intersections is expressed by the following equation [21]:

The delay of cars at unsignalized downstream intersections is expressed by the following equation [21]: where is the proportion of effective green time with a cycle time and is calculated by (4). is equal to 1 if intersections are unsignalized:

is the degree of key entrance lanes at intersections and can be expressed by the following equation; is saturation flow: where is the traffic flow volume including vehicles and buses per lane (pcu/s). is the degree of saturation at bus stops which can be described as where is the proportion of time when buses are dwelling at bus stops, whose value is equal to the product of the average arrival rate of buses to bus stops, (vehicle/s), and the average boarding and alighting time of passengers at bus stops, (s). and are obtained by where is the arrival rate of bus berth ; is the average alighting and boarding time of passengers for bus berth ; and is the number of buses dwelling at bus berth .

Moreover, is the distance between bus stops and the vehicle stop line of downstream intersections (m). are coefficients and values are 106.5, −0.09, 0.07, 1.27, −0.53, and 0.57, respectively.

3.1.3. Travelling Time of Buses

The influenced area for bus stops is defined as the area that is located between upstream intersection and downstream intersection in this paper. The process of buses passing the area is divided into seven stages, including accelerating to pass upstream intersections, running at a constant speed, decelerating and stopping at bus stops, boarding or alighting of passengers at bus stops, accelerating to leave bus stops, running at a constant speed, and decelerating and stopping before downstream intersections. The speed variation of buses with their locations which are differently positioned between two adjacent intersections is illustrated in Figures 2, 3, and 4.

The total travelling time of buses within the influenced area is computed by the following equation in terms of the above analysis: where is the consumed time of buses within each stage as mentioned above, respectively.

The corresponding travelling distance for each stage as introduced before can be calculated by different methods. For the stages where buses are accelerating or decelerating, the travelling distance is obtained by where is the order number of the stages as shown in Figures 2 to 4 that buses run through (). is the travel distance of buses with the stage . is the initial travel velocity of buses at the beginning of process ; is the maximum travelling velocity of buses in the stage . is the accelerated velocity and is set to 1.4 m·s−2 [22].

For the stage when buses dwell at stops, the boarding and alighting time of passengers, , that is, , is obtained by (10) given in terms of [16]: where and are the dwelling time of buses during boarding and alighting passengers at curbside and bay-style stops, respectively. , , and are the number of passengers boarding and alighting at different bus doors. is the load factor of buses, which is usually calculated through the ratio of passenger number to bus capacity.

For the stages when buses are running at a constant speed, the corresponding travelling times and distances are calculated by (12). The constant operational velocity is set to 35 km·h−1, which is equal to 9.7 m·s−1 [22]:

The abscissa of bus stop center is required to satisfy

3.2. Constraints

The performance of vehicles at upstream intersections will not be affected by queuing buses. Thus, the location of the bus stops is required to be confined within a certain range and the corresponding constraint is expressed by where is the set of alternative locations of bus stops. is the vehicle length of buses. is the number of buses permitted to queue before bus stops.

The constraint indicates that the minimum distance from bus stops to the upstream intersection is more than the permitted maximum length of queuing bus while the maximum distance is less than the road length from upstream to downstream intersections.

The number of buses permitted to queue can be calculated through the following equation as the service system of buses is really an M/M/C queuing system: where is the number of bus berths. The probability without any bus at stops, , is expressed by

The total arrival rate of buses at bus stops, , is attained through where is the arrival rate of buses at bus stops along the bus route . The average service rate of buses at stops, , is obtained by

The service intensity of bus stops is expressed by

4. Solution of the Multiobjective Optimization Model

The mathematical model for the optimal locations of bus stops connecting subways near urban intersections is proposed as below where is a decision vector containing all continuous variables, denotes multiple objective functions, gives the nonlinear objective function, and denotes the set of feasible solutions.

For multiobjective problems, a set of Pareto optimal solutions form the Pareto frontier of them. Decision makers (DM) usually select a particular Pareto solution based on additional preference information about the objectives.

4.1. Multiobjective Analysis

The total walking distance of passengers, delay time of cars through intersection, and travel time of buses are given in different units so that they are unable to be measured and compared directly. A normalization formula on the basis of the original Pareto frontier is established as follows: where , , and are the normalized, maximum, and minimum values of the objective, respectively. The original Pareto frontier can be converted to the normative Pareto frontier after the normalization procedure.

The best alternatives under different weights on objectives by the distance-based method can be figured out. The distance-based method is expressed by where , , and are the weight coefficients of three objectives, which represent the preferences of managers or decision makers on the objectives, respectively. The weight coefficients also denote the relative importance of the objectives and are more important when they are more close to 1.

4.2. Solution Method

There are a large variety of methods to solve multiobjective optimization problems, including exact methods, such as the -constraint approach and Tchebycheff algorithm, and heuristic algorithms, such as evolutionary, Tabu search, ant colony, and particle swarm algorithms.

For the multiobjective model in (20), all feasible positions of bus stops, , should lie within the range between and . The model is thus converted into a nonlinear integer programming model as bus stops are permitted actually to set at those positions in meter in practice, which means that is an integer. The solution tool, Lingo, is employed here to obtain the Pareto frontier of the proposed model and optimal solutions.

5. Case Studies

One bus stop in the Xizhimen terminal of Beijing is selected to implement to verify the validity and practicability of the proposed method in case studies. Moreover, a sensitivity analysis of all alternatives on the weights of three objectives is conducted to provide decision-making supports for bus operators.

5.1. Set-Up

The Xizhimen terminal is a large transportation terminal in the northwestern region of Beijing. Passengers can transfer here among five transportation modes, including railway, subway, private cars, taxis, and bicycles. The daily amount of passengers served by the bus stops in the Xizhimen terminal is over 30000 according to the field survey. The layout of the Xizhimen terminal and its surrounding facilities is described in Figure 5.

The bus stop 1 with 3 bus berths is one of busiest bus stops in the Xizhimen terminal which is located on the eastern side of Gaoliang Bridge Street and between two intersections as shown in Figure 5. The stop is 70 and 170 meters away from the upstream and downstream intersection along the Gaoliang Bridge Street. There are two major passenger flows relating to bus stop 1: interchange passengers from/to Line 2, Line 4, and Line 13 of subway and those from/to other bus stops. A few of passengers from/to the downstream intersections are not considered in the case studies.

The filed data on traffic current situation was collected during the rush hour from 17:30 to 18:30 on December 27, 2013. The effective signal cycle and green time are 140 and 100 seconds at the intersection of Gaoliang Bridge Street and Beizhan North Street. The average length of buses serving passengers here is about 12 meters. The traffic volume and saturated flow rate including buses, cars, and taxis are 526 and 1260 pcu/lane per hour. Other data on buses and passenger flow are listed in Table 1.

The other basic parameters for the proposed model are calculated in light of the filed survey data and listed in Table 2.

The three objectives of the proposed model under current situation can be computed based on the above field data and calculated parameters. The total walking distance of passengers, delay time of cars through intersection, and travelling time of buses are 93.6 kilometers, 29.7 seconds, and 38.6 seconds, respectively.

5.2. Optimization Alternative Analyses

The suitable position of bus stops near intersections is given within a range of 80 to 150 meters away from their upstream intersection in terms of the Code for Design of Urban Road Engineering [23]. To ensure the completeness of this research, the range of possible positions of bus stops is set from 50 to 240 meters. The original Pareto frontier through the proposed model is illustrated in Figure 6.

The variation trend of three objectives with the positions of bus stops is listed in Table 3. “” and “” denote the values of three objectives increase and decrease with the growth of , respectively. The maximum and minimum values of three objectives from Table 3 are 435.24 and 84.24 km, 36.03 and 29.62 seconds, and 38.6 and 31.67 seconds, respectively.

The normalization of Pareto frontier is shown in Figure 7.

A group of optimal solutions are obtained and listed in Table 4 under the different combination of weight coefficients according to the preferences of decision makers or managers.

All Pareto solutions for 66 weight combinations are summarized into 5 groups due to values from Table 4. The details on the grouping and improvement rates of solutions comparing to the current position of the bus stop, that is, is 70, are listed in Table 5. It indicates that weight combinations have remarkable influences on optimal solutions. For example, if decision makers pay more attention to the reduction of the delay of cars through intersection and the travelling time of buses, the optimal position of bus stops will tend to be close to the upstream and downstream intersection, respectively. If the total walking distance of passengers is considered most, the optimal position will be located near the subway stations and keep a certain distance from adjacent intersections. Thus the balanced optimal position of bus stops is here recommended to be 90 meters away from the upstream intersection due to the average improvement rates of three objectives in Table 5.

6. Conclusion

Locations of bus stops connecting subways are much crucial to improve the efficiency and level of service of bus operations. A multiobjective optimization model to determine suitable locations of bus stops is proposed in this study considering the total walking distance of passengers, delay time of cars through intersections, and travelling time of buses between adjacent intersections. The proposed model is available and effective to determine locations of bus stops through case studies where an empirical study on the bus stop at Xizhimen in Beijing is carried out. The balanced optimal location of bus stops connecting to subways is recommended to be near its upstream intersection and subway station through this study when it meets the requirement from the Code for Design of Urban Road Engineering [23].

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are grateful to the National Basic Research Program of China (2012CB725406) and the National Natural Science Foundation of China (71001006 & 71131001) for their financial supports on this work.