#### Abstract

Different configurations of Bus Rapid Transit (BRT) system may cause different additional riderships. In this paper, in terms of network traffic equilibrium assignment principle, the additional riderships estimation model based on Variational Inequality (VI) model is presented. The bus frequency is related to variables including the travel time, the residence time in terminals, and the dwelling time at the stops. The additional riderships are translated into network additional traffic flow firstly. Given the bus frequency, VI model can be turned into Stochastic User Equilibrium (SUE) model to calculate the other variables. The similarity diagonalization method is used to calculate the elastic bus frequency and finally the network additional traffic flow can be computed. The additional riderships under different configurations of BRT system are compared in the numerical test. The results show that the additional riderships under different configurations have large differences and occupy a high percentage of the total ridership.

#### 1. Introduction

Although the foundation ridership can be generated by the increase of the population and the growth of the economy, additional ridership can be generated by Bus Rapid Transit (BRT) system. Different configurations of BRT system induce the diverse scale of additional riderships. Good configurations of BRT system cause the increase of traffic demand and make additional ridership account for a great percentage of the total ridership. The additional riderships induced by different configurations of BRT system are proposed in this paper, which should be considered in the ridership estimation.

The different configurations of BRT system have been attempted by lots of researchers. Laporte et al. [1] reviewed the main optimization methods for the BRT planning. Abdelghany et al. [2] developed a modeling framework for the planning of BRT services in urban transportation network. Aiming at minimizing the total travel time of passengers, Li et al. [3] presented an optimization model for the BRT planning. Schmid [4] proposed a hybrid metaheuristic approach based on large neighborhood search and liner programming to solve the bus rapid transit design problem. In order to capture the effects of BRT services on urban transportation, various methods have been developed, and a majority of those methods are based on simulation. Salem et al. [5] used a CORSIM model to conduct the benefit/cost analysis of the BRT service. Yagi and Mohammadian [6] simulated the BRT development to study the variation of modal split between automobiles and bus rapid transit ridership. Gunawan et al. [7] presented a numerical simulation based on discrete-event approach to identify variables which affects the performance of BRT system. The microscopic simulation techniques are also widely used in the research field of BRT. Yu et al. [8] applied GPS data in the VISSIM-based simulation for BRT system in Beijing. Godavarthi et al. [9] used a microsimulation to find the optimum volume/capacity ratio on BRT routes. Cervero and Kang [10] measured the impacts of implementing BRT on adjacent land usage and land value, by using multilevel models. Based on six characters of BRT system including infrastructure, transport capacity, service level, economic results, safety and emergency management, and energy saving and emission reduction [11], established an evaluation methodology for BRT operation. To study the modal shifts to BRT from other modes like automobiles, normal buses, Wang et al. [12] developed a binary logistic analysis approach. Falbel et al. [13] studied the impacts of implementation of BRT on the traffic improvement over the urban transport network, including the additional ridership to the transit system. Currie and Delbosc [14] used several regression models to explore the relationship between BRT design features and increment of ridership. Sun et al. [15] established a numerical model for the headway optimization as well as scheduling combination of BRT vehicles to increase ridership and improve operation performance of BRT. There are few previous literatures that studied the additional ridership under different configurations of BRT system. This paper developed a numerical test based on Variational Inequality (VI) model for additional ridership estimation under different configurations of BRT system. The paper is organized as follows. Section 2 is the problem description. Section 3 is the model development and algorithm for the proposed model. A numerical test is given in Section 4. Section 5 presents the conclusions.

#### 2. Different Configurations of BRT System Description

In transportation system, there are many literatures about traffic assignment. Some equilibrium analysis for urban transportation networks were proposed in literature [16, 17]. A traffic equilibrium assignment principle was presented, in the principle, passenger’s selection behavior of a transit note was considered. The obtained results showed that the method was effective for network traffic equilibrium assignment [18]. In this paper, the additional riderships estimation under different configurations of BRT system is based on the network traffic equilibrium assignment principle. In terms of the network, the addition ridership account for less than 25% of the total ridership which is the main restriction in the model [19].

The bus frequency which can reflect the travel time has a function relation with the traffic flow on the road link [20]. It is the main variable in the model, which can represent the requirements of the different hardware configurations and the standards of operational service.

The traffic flow on the road link of BRT system is a function of the bus frequency, which is determined by the bus operation cycle and the bus number of each line. denotes the travel time on the road link of the bus route , and the bus operation cycle is determined by the traffic conditions of network and the dwelling time at the stop :where and represent that the road link and transfer stop belong to the bus route , respectively. represents the residence time in the terminal. denotes the number of terminals and indicates there is only a terminal, while indicates there are two terminals.

denotes the bus number of the route ; then bus frequency of the route can be calculated with the following equation:The travel time has a relationship with adopting signal priority control and the configurations of route. If the signal priority control is adopted, the intersection delay and the travel time will be shortened. When the route is a part of BRT system, the security and the speed of driving will increase. According to the travel time survey in Changzhou BRT, for 1 kilometer section, the variation range of the travel time of different configurations is as shown in Table 1.

If the bus frequency is held constant, the residence time in the terminal is related to the number of buses on the route. And the existence of more buses means longer residence time in the terminal. In terms of the BRT system under different configurations, the bus capacity varies with the bus model. For the convenience of calculation and comparison, the different BRT bus models are replaced by a standard BRT bus model, the size of which is 12 m, and the conversion coefficient is the ratio of bus capacity.

The dwelling time at a stop is consisted of the time of deceleration and acceleration, the time of switching door, and the time of boarding and alighting. When the bus model, the platform, and the lane are constant, the dwelling time at a stop is determined by the number of the passengers boarding and alighting. The boarding and alighting time are determined by passenger number and the average time of per boarding and alighting passenger , , respectively. Wirasinghe and Szplett [21] presented that if the given road link traffic flow is , the number of passengers alighting and boarding at the stop on the bus route can be yielded as follows:where and denote the number of passengers alighting and boarding at the stop on the bus route , respectively. and denote the average time of per boarding and alighting passenger, respectively. Thus, denotes the boarding time and denotes the alighting time at the stop . denotes the percentage of riderships who choose route on the road link , and denotes the ridership on the road link . , , and are defined as exponential functions, where the superscripts with plus sign and minus sign indicate the departure terminal and the destination terminal on the road link . The dwelling time at the stop can be yielded as follows:where denotes the time of deceleration, which means the time from vehicle slowing down to a full stop. denotes the time of acceleration, which means the time from starting acceleration to moving in a steady speed. denotes the time of opening and closing door. Passengers alighting are assumed to be completed before passengers boarding at the stop. The main factors influencing the dwelling time at a stop are the average time of per boarding and alighting passenger , , respectively. The average time of per boarding and alighting passenger is influenced by fare collection method, the density of passengers inside and outside the vehicle, width of the bus door, whether to keep the platform and the vehicle floor in the same horizontal plane, and other uncertain conditions like weather.

Thus, the bus frequency varies with the different configurations of BRT system, which reflects the hardware configuration and the operation state significantly. When the vehicle model, the running ways form, and the platform style are fixed, the dwelling time of a vehicle can be presented as a function of the traffic flow as shown in formula (5). denotes the traffic flow on the road link in the public transportation network, which can be transformed to ridership based on the average capacity and the vehicle conversion coefficientWhen the configuration is determined, bus frequency is also a function of the traffic flow ; the elastic bus frequency can be yielded as follows:where denotes the bus number of the route and is the operation cycle related to the traffic flow on bus line . is the traffic flow in the transit network. The traffic flow can be converted to ridership by bus capacity and vehicle conversion coefficient.

#### 3. The Additional Riderships Estimation Based on VI Model

In the model, the additional ridership is translated into additional traffic flow in the network . The relationship between traffic generation and traffic distribution can be represented as follows:where denotes the additional traffic distribution. The relationship between the traffic distribution and the path flow can be formularized as follows:where denotes the route flow between and on route . The relationship between the traffic flow in the network and the path flow can be defined as follows:where denotes an incidence matrix regarding road link-path. If the road link is on path , it is equal to 1; otherwise it equals 0. denotes the path flow driven by the current traffic demand, while denotes the additional path flow driven by the additional traffic demand. Thus,where represents the number of buses on the route , which is often held constant:The passenger flow at the stop is a function of the bus frequency. It makes it hard to obtain the Karush-Kuhn-Tucker (KKT) condition of the Stochastic User Equilibrium (SUE) model [22] and leads to the difficulty of solving the model. When a variable (such as bus frequency) is held constant, the VI model is used to avoid this problem. The VI model can be transformed to the UE optimization model, which has a unique solution.

##### 3.1. VI Model Formulation

It has been proved that the network equilibrium model and the VI model are equivalent. Assuming that the traffic cost and the additional traffic demand , let , where is a vector of parameters; denotes the path flow driven by the current traffic demand; denotes the additional path flow driven by the additional traffic demand; denotes the additional traffic distribution; is the bus frequency vector of route. According to the SUE model, a continuous function of each vector variable can be referred to as formula (12) by the letters , , , and in vector:where , , .

The VI model can be expressed as follows:where there is a constraint . In formula (12), represents the travel cost on the road link ; denotes the user waiting time on the road link and denotes the destination attractive function. The passenger waiting time on the road link is determined as . denotes the bus frequency on the road link and is equal to the bus frequency of route , . The parameter is determined by the distribution value of bus headway, where represents an exponential distribution and represents a uniform distribution [23]. is commonly used in the practical analysis.

Assuming that is a solution of the VI model (13) and the linear independence condition is satisfied at , then the GKKT (General Karush-Kuhn-Tucker) conditions are satisfied.

The constraint can be yielded as follows:and the constraint can be yielded as follows:Formula (17) is based on the condition that the addition ridership accounts for less than a percent of the total ridership. According to the TCRP118, the percentage was suggested as 25%.

In formula (21), the part “” is based on the practices. Usually the time interval of two consecutive buses is in the range of 1 min to 6 min; thus veh/h. The GKKT conditions are defined as follows:where, , , .

By invoking (14), (16) is solved as follows:Let present the travel cost on route and present the user waiting time on route . Equations (15) and (16) can be derived as follows:where is the logit model based on the path selection probability.

##### 3.2. Solution Algorithm

When is given, let ; the objective function can be written as follows:which subjects to constraints (14)–(21). This model is a SUE assignment problem and satisfies the KKT conditions.

In public transport problems, the passenger travel time is independent of the traffic flow [24–27]. Once the bus frequency is certain, the expected passenger travel time is independent of the path flow and the bus frequency: denotes the travel time at the road link . indicates the waiting time at the stop, and is the random error. As long as is fixed, the optimization problem (25) can be solved without iteration. Meanwhile, the VI model can be transformed into the optimization SUE model to solve the other two variables and . This characteristic can be used in the similarity diagonalization method to iterate and solve the VI model based on the elastic frequency.

#### 4. Numerical Test

##### 4.1. Network Basic Information

In order to demonstrate the performance of the proposed model and solution algorithms, numerical tests were carried out on a real transit network in Changzhou city (see Figure 1). In this study, Wujin district (node 1) and New North district (node 2) are considered as the origins, while Jin Ling road (node 3) and River road (node 4) are considered as the destinations. (from node 1 to node 5), (from node 5 to node 6), and (from node 6 to node 3) are road links of route , while and are the road links of route .

The road link information of the transit network shown in Figure 1 is given in Table 2. , , and are BRT routes. The traffic demand is shown in Table 3.

The dwelling time at a stop under different configurations can be yielded as follows: primary configuration: intermediate configuration: advanced configuration: the dwelling time for a regular bus:

##### 4.2. Results

Through the numerical test, it is concluded that the bus frequency under different configurations of BRT line 1 ( in the network) increases with the traffic flow. The increase of the bus frequency under primary configuration is bigger than that under advanced configuration (Figure 2).

Figure 3 shows the change of the bus frequency on each route with the increase of traffic flow, where the BRT system is under intermediate configuration. The bus frequency variation of BRT lines (, , and ) is more uniform than that of common ones (, , , and ). The bus frequency of is the most sensitive to the additional traffic flow, while the frequency of is lower than others because of the shorter travel distance. In addition, Figure 4 shows the additional traffic flow occupy. From Figure 4, it can be found that the total traffic flow increases obviously and the percentages of total traffic flow under different configurations are 11.9%, 18.8% and 22.8%, respectively compared with the original traffic flow (5000/h). The additional traffic flow occupies a percentage of total traffic flow that is 11.9%, 18.8%, and 22.8%, respectively.

#### 5. Conclusions

In this study, the additional riderships prediction method based on VI model under different BRT system configuration are presented, as well as the solution algorithm. The VI model is turned into the SUE model, in order to avoid the solving difficulty resulted from the elastic frequency. A transform between the VI model and the UE model is proposed in the solution, where the similarity diagonalization method can be used. The additional traffic flow under different configurations of BRT system is computed in the numerical test. The results show that the potential riderships of different configurations are obviously different, and the bus frequency under primary configuration has the largest increase with the rise of network traffic flow. It also can be seen that the potential riderships have a big proportion in the total ridership.

Considering those differences in the construction of BRT, different elements of the system should be chosen in different development stages to achieve the objective of TOD, in addition to saving the cost of investments.

#### Conflict of Interests

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

#### Acknowledgment

This work is supported by National Natural Science Foundation of China 51408314.