Bus timetabling considering passenger satisfaction: An empirical study in Beijing

Bus timetabling considering passenger satisfaction: An empirical study in Beijing

Accepted Manuscript Bus Timetabling Considering Passenger Satisfaction: an Empirical Study in Beijing Hua-Yan Shang, Hai-Jun Huang, Wen-Xiang Wu PII: ...

2MB Sizes 0 Downloads 67 Views

Accepted Manuscript Bus Timetabling Considering Passenger Satisfaction: an Empirical Study in Beijing Hua-Yan Shang, Hai-Jun Huang, Wen-Xiang Wu PII: DOI: Reference:

S0360-8352(19)30069-5 https://doi.org/10.1016/j.cie.2019.01.057 CAIE 5683

To appear in:

Computers & Industrial Engineering

Received Date: Revised Date: Accepted Date:

16 December 2017 4 December 2018 27 January 2019

Please cite this article as: Shang, H-Y., Huang, H-J., Wu, W-X., Bus Timetabling Considering Passenger Satisfaction: an Empirical Study in Beijing, Computers & Industrial Engineering (2019), doi: https://doi.org/ 10.1016/j.cie.2019.01.057

This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Bus Timetabling Considering Passenger Satisfaction: an Empirical Study in Beijing

Hua-Yan Shanga,, Hai-Jun Huangb, Wen-Xiang Wuc a

Information College, Capital University of Economics and Business, Beijing 100070, China

b

School of Economics and Management, Beihang University, Beijing 100191, China and Key Lab of

Complex System Analysis Management and Decision, Ministry of Education, China c

Beijing Key Lab of Urban Intelligent Traffic Control Technology, North China University of Technology, Beijing 100144, China



Corresponding author. E-mail address: shanghuayan@126. com (H. Y. Shang) 1

Bus Timetabling Considering Passenger Satisfaction: an Empirical Study in Beijing

Abstract: Public transport is the key component in sustainable transportation. In China, large subsidies are provided by the national and local governments to encourage the use of public transport systems. However, high commuting demand leads to such incredible phenomenon that passengers have to climb windows for riding in buses during peak hours. Thus, a case study was conducted in a Beijing’s transport hub which covers 116 bus lines and ten stations. The existing bus timetabling is surprisingly manually made and 26.53% of the bus supply undertakes 59.19% of the bus demand, which results in unbearable long waiting time and on-board crowding. To find a bus timetabling method which is practical and applicable in China, three procedures are carried out. Firstly, the passenger satisfaction and the bus transit efficiency are simplified and formulated. Secondly, a preliminary bus timetabling method with consideration of passenger satisfaction is proposed for optimizing the bus frequency and headway. Thirdly, the bus timetabling is optimized by embodying a balance between the passenger satisfaction and the bus transit efficiency, subject to the constraints of load factor. This method is verified by field data and the results demonstrate that it has the advantage of reducing waiting time and lessening on-board discomfort. Keywords: Public transport; Bus timetabling; Waiting time; Congestion effect

1. Introduction Population growth exerts considerable pressure on infrastructure and natural resources in urban regions. For example, China has 102 extra large cities which have populations of more than ten million. By 2016, its urbanization level has risen to 56.1%. One of the most obvious problems in daily life is the transportation system, including environmental pollution, road accidents, and road congestions. How additional transportation demand will be served in growing urban regions is an important issue for achieving sustainability. Indeed, the relationship between the transportation system, urban form, trip demand, and energy use is paramount in addressing the challenges presented by urban growth. This may be attributed to considerable economic inefficiency and environmental degradation associated with excessive private vehicle travel. Thus, public transportation is recognized as a key component in the management and planning of urban regions (Otto and Boysen, 2014). With rapid economic growth of China, traffic congestion can be found everywhere in 2

China’s large cities. To alleviate traffic congestion, the Chinese government is insistently and constantly constructing a comprehensive transportation system with the use of a policy priority to public transport. The policy of vigorously developing the public transportation has been listed in The 13th Five-year Plan for National Economy and Social Development, and issued by the State Council (13th Five-Year Plan, 2015–2020). The general objective of this policy is to build a modern public transport system for each city so as to provide people fast, efficient, safe, green, and economical services. Accordingly, large subsidies are provided by the national and local governments to encourage the use of public transport systems. In Beijing, China, the government subsidies have grown by seven times from 2006 to 2013. The welfare was 3,672 million yuan (US$ 530 million) in 2006, and increased to 250,000 million yuan (US$ 3,600 million) in 2013. In 2014, Beijing government tried to adjust the bus fare, but the fare was still defined well below the system cost. The government bears over 60% of the total cost of public transport now (Statistical Yearbook, 2017). Unfortunately, the current public transport systems in most Chinese cities are not satisfactory. Low speed, low punctuality, poor comfort, inconvenient transferring, and low coverage, all these are problems being faced by the Chinese public transportation systems. Thus, we have conducted a series of field surveys to gain an in-depth understanding of the bus operations in real traffic, including a survey on public transport in Guomao Bridge transport hub (Shang et al., 2014). Nearly 1,500 buses are dispatched from this hub or pass by it with stopping during evening peak hours (17:00-19:30) and the bus fleet is called the world’s largest fleet by Chinese netizens. However, it has a large number of passengers queue at bus stops (see Fig. 1a). Some passengers even have to climb up the bus windows for riding in a bus (see Fig. 1b). Guomao Bridge hub is a classic transit-oriented transport hub. It is located at a central business district (CBD) of Beijing. A variety of public transportation services (e.g., metro trains, urban and suburban bus services) are provided. Tidal-traffic patterns and suburban bus services are prominent. Yanjiao Town is a large residential district in the suburbs. Only one corridor, named Beijing-Tongzhou Expressway, connects the work zone and the home zone (Fig. 2). On each weekday morning, over 300,000 commuters leave homes from Yanjiao Town and travel to Guomao Bridge along the corridor, and return in the afternoon. The suburban bus routes are operated in the region between Guomao Bridge hub and Yanjiao Town district, with fixed routes and a unique route number/name (e.g. Line 813). These routes are very important for the suburban commuters with stable passenger flows. Although they operate similarly with urban bus routes, this paper only focuses on the headway design of these trunk bus routes. 3

The objective of this paper is to find a bus timetabling method which is practical and applicable in China. This method relies heavily on the research background of Guomao Bridge transport hub (detailed in Section 3). Three procedures are carried out: Step 1: Formulate the passenger satisfaction and the bus transit efficiency. There are various important attributes that cannot easily be quantified and measured (Ceder, 2016). For simplicity, waiting time and bus comfort are considered to be the main attributes of the passenger satisfaction according to our survey data. Step 2: Develop a preliminary bus timetabling. The max load method is used to determine the frequencies and headways with smooth transitions are proposed. Step 3: Optimize the bus timetabling. An optimization model is developed with regard to the trade-off between the passenger satisfaction and the bus transit efficiency, subject to the constraints of load factor. In this way, the scheduler will be in a position not only to expedite manual/semi-manual tasks, but also to compare procedures and methods in regard to the trade-off between the passengers’ comfort and the operating efficiency. The rest of the paper is organized as follows: In the next section, a literature review of bus scheduling and passenger satisfaction is provided. The survey and some analyses are presented in Section 3. In Section 4, we formulate the passenger satisfaction and develop the new bus timetabling method. The practical application of the method is illustrated in Section 5. Section 6 concludes the paper.

2. Literature review In the last decades, a fruitful development of models and solution techniques were addressed in bus transport systems. Earlier researches included the design and evaluation of vehicle holding strategies (Osuna and Newell, 1972; Barnett and Kleitman, 1973; Barnett, 1974). Models proposed in these studies aimed to minimize the number of buses required under rather restrictive conditions. At that time, there were not good computational techniques to solve these models and survey data to verify them. Afterwards, the research interest was extended to investigate other strategies such as stop skipping, signal priority, reduction of bus dispatch uncertainty, short-turning, and headway control (Turnquist and Bowman, 1980; Furth and Wilson, 1981; Abkowitz and Engelstein,1984; Wilson et al., 1992; Khasnabis et al., 1999). These studies were empirically validated with data from actual transit operations. For example, Turnquist and Bowman (1980) utilized the data from a bus route in Evanston of Illinois, and Abkowitz and Engelstein (1984) employed the data from Los Angeles of California. It was generally supposed that the strategies proposed were expected to be 4

effective in reality as they were examined by real data (Strathman et al., 2002). Later, passenger waiting time at a bus station and load factor were paid much attention and deeply studied (Rajbhandari et al., 2003; Mishalani, 2006; Badami and Haider, 2007; Toledo et al., 2010; Ceder, 2013). Timetable synchronization is regarded as a useful strategy utilized to reduce passenger waiting time and improve service connectivity. Fonseca et al. (2018) studied the integration of timetabling and vehicle scheduling to reduce passenger-associated costs and operating costs. However, most of the studies on timetable synchronization design focused only on minimizing the transfer waiting time (Liu and Ceder, 2017). Load factor was widely used for frequency and vehicle-size determination. Max-load method is still popular now. It is a general agreement that the in-vehicle and out-of-vehicle times cannot be equally treated since the out-of-vehicle time is much more onerous for determination (Ei-Geneidy et al., 2006). Thus, the load factor can be regarded as an un-direct variable affecting passenger waiting time. These studies consciously inferred that high load factor could reduce waiting time. However, they ignored that higher load factor means higher level of on-board crowding and lower passenger satisfaction. The global public transit problem is computationally intractable and can hardly be tackled at once. Desaulniers and Hickman (2007) divided it into a set of sub-problems that were usually solved sequentially at various stages of the planning process (strategic, tactical, and operational) and during operations (real-time control). They provided a systematic way of reviewing state-of-the-art models and approaches for solving the public transit problems. Ibarra-Rojas et al. (2015) followed this classification. Guihaire and Hao (2008) presented a global review of the crucial strategic and tactical steps of transit planning: the design and scheduling of the network. They followed a five-step planning process including network design, frequencies setting, timetable development, bus scheduling, and driver scheduling (Ceder, 2016). Some studies aim to get a bus schedule which can implement “fewer buses, no waiting time” (Ceder, 2016), but congestion cost or crowding effects have not been fully elucidated (Prud’homme et al, 2012; Tirachini et al. 2013). From the perspective of a transit provider, optimal transit service can be characterized by predictable passenger activity and few service reliability problems. Passengers expect that buses arrive promptly at stops that access, egress, and waiting times are minimized. Also, they seek to minimize their combined in-vehicle and out-of-vehicle travel times. Thus, many studies tried to reveal the quality paradox between service supply and passenger satisfaction (Michaelis and Schöbel, 2009; Friman and Fellesson, 2009; Lai and Chen, 2011). These works mainly focused on such aspects as qualitatively evaluating the reliability, frequency, travel time, fare, comfort, cleanliness, network coverage, safety, and the information released to passengers (Friman and Fellesson, 2009). However, for 5

passenger satisfaction, there are too many quality-based attributes and the importance of attributes being negatively or positively perceived by individual passengers depends on their preferences. It is hard to quantify passenger satisfaction thus qualitative analyses are preferred. In all the attributes of passenger satisfaction, congestion/crowding effects in transit systems have attracted much attention. Congestion leads to on-board discomfort, denied boarding, and low service reliability. Kurauchi et al. (2003) considered the failure-to-board probability and Schmöcker et al. (2008) developed a quasi-dynamic frequency-based model. In addition to capacity constraints, the on-board discomfort effect is presented in Trozzi et al. (2013). Sumalee et al. (2009) and Schmöcker et al. (2011) introduced fail-to-sit probability to satisfy the set of priority rules and the seat capacity constraint. In-vehicle seat priorities were also incorporated by Leurent et al. (2014). Tirachini (2013, 2014) identified and analyzed the interplay between congestion and crowding externalities in the design of urban bus systems. Nuzzolo et al. (2012) and Pel et al. (2014) proposed models that estimate on-board discomfort based on average crowding levels (e.g. volume/capacity or load/seats ratios). These models would underestimate the congestion effect since compared to less crowded vehicles, passengers experienced more discomfort in overcrowded vehicles. Cats et al. (2016) took three distinct congestion effects into account: on-board discomfort, denied boarding, and irregular vehicle arrivals. However, they focused on modeling these effects and the bus timetabling was not considered in their paper. On the other hand, the bunching phenomenon became the most noticeable phenomenon (Cats et al., 2016) and He (2015) suggested an anti-bunching strategy to improve bus schedule and headway reliability by using available accurate information. And yet, up to date, there has not been an explicit function which is calibrated by survey data and can be used to measure the passengers’ discomfort cost caused by crowding in vehicle. Huang et al. (2007) and Tian et al. (2007) quantitatively considered the congestion cost which was formulated as a function of congestion degree and travel time. They conducted a survey at a subway station in Beijing. However, they assumed that the function was linear for deriving an analytical solution. The impact of congestion/crowding effects has seldom been considered in scheduling bus timetable. To our knowledge, discomfort due to crowding in vehicle is usually neglected in the objective functions. Especially, long travel time due to long distance in a crowding vehicle always causes additional psychological discomfort. Jiang et al. (2014) proposed an optimization model to reach the minimization of the waiting time cost, the crowding cost, and the operating cost. But the passenger demand was assumed to be relatively low. With rapid urbanization in China, many transport hubs have sprung up to serve the transportation between CBD and outskirts, similar to Guomao Bridge hub. A variety of public 6

transportation services are provided, among which the suburban bus services account for a large proportion. However, different from those traditional suburban ones, the demands for these outskirt bus services are high. They operate similarly with the urban bus routes, but the trip times are usually long. Thus, passengers expect more efficient service (e.g., high bus frequency) to reduce the waiting time and improve the in-vehicle comfort. In contrast, bus companies are unwilling to operate the routes with low ridership. They always expect a longer headway to reduce operating costs. Moreover, since the passenger flow decreases from downtown to outskirts, bus companies expect higher ridership at the departure stops. Hence, it is important to balance the interests of both passengers and bus companies. In this paper, our model tries to reach the balance of the passenger satisfaction and the bus transit efficiency. The waiting time cost, the discomfort due to crowding in vehicle, and the load factor are considered. The objective of this paper is to find a bus timetabling method which is practical and applicable in China. Thus, some simplified formulations are selected and a quantified modeling approach is adopted to account for the on-board crowding effects.

3. Survey In this section we first introduce the survey on public transport in Guomao Bridge hub and then address the initial findings related to our studies. This section is the model basis of Section 4 and the data source of Section 5. 3.1. Data collection Guomao Bridge transport hub of Beijing was selected for data collection during the evening peak hours (17:00-19:30). Fig. 3 shows the map of this hub with 1.1 square kilometers. All the parking lots, subways, bus stops, and surrounding buildings within this hub are our objects for investigation. We recorded the passenger flows generated from these sites. The operations of 116 bus lines at ten stations, including passenger arrival rates, bus arrival time intervals, passenger waiting times, and queuing lengths, were observed and statically analyzed. The survey lasted for three months from April 2012 to June 2012, and 60 undergraduate students were recruited. Direct observation and video tape recording were simultaneously employed. All field works were finished on weekdays with good weather conditions. 3.2. Initial findings (1) About 80,000 commuters go home during the peak hours of each workday in this hub,, 7

among which 5,857 are transported by private car or taxi, 27,660 by subway, and 42,238 by bus. 1,500 buses of 116 lines serve the commuters from ten stations. (2) 19 bus lines (16.38% of the total lines) serve on suburban bus routes. They dispatch 398 buses (26.53% of the total bus supply) during evening peak hours while transport nearly 25,000 passengers (59.19% of the total bus demand), which means 26.53% of the bus supply undertakes 59.19% of the bus demand. Thus, overload is inevitable. (3) 97 bus lines (83.62%) served on the urban bus routes. 1102 buses of these lines pass by Guomao Bridge transport hub with low bus loads. (4) The special geographical condition leaves the suburban bus commuters with no choice. The sketch map of the suburban bus routes is shown in Fig. 4. These routes depart from Guomao Bridge hub. Each has an independent bus line and a unique route number/name. The passenger flow decreases from the hub to the home zone. So, optimal timetable at the departure stops in the hub is highly important. (5) All vehicles of the suburban bus lines seem to be scheduled manually. Total number of

buses is basically constant but headways fluctuate largely for each line. Bunched and delayed bus arrivals are common. Fig. 5 shows the bus headways for Line 807 in the direction of vertical axis (in minutes), whereas the horizon axis represents the bus order from the first vehicle to the last one. The average headway was three minutes during the surveyed periods. However, the vehicles did not arrive evenly. A vehicle might be late for over ten minutes while other headways might be less than one minute (delayed arrivals vs. bunched arrivals in Fig. 5).

Consequently, vehicles’ uneven arrivals lead to long passengers’ waiting times.

(6) There are significant differences among actual, perceived, and tolerable waiting times for riding in the suburban buses1, as shown in Table 1. The average tolerable waiting time is nearly half an hour because the suburban passengers have no other choice but to wait. The average perceived waiting time is much larger than the actual one, i.e., the passengers over-evaluate their waiting times. In other words, long actual waiting times make the passengers anxious and cause higher perceived waiting times. 1

The actual waiting time is the time that a passenger actually experiences, i.e., the real gap between times of reaching and leaving a bus station. The perceived waiting time which was first defined by Mishalani et al. (2006), is the time that a passenger subjectively feels and evaluates. This time is usually larger than the actual waiting time since passengers generally over-evaluate their waiting times due to impatient mood. The tolerable waiting time is the time that a passenger is maximally willing to keep waiting for a bus. Over this time, he/she will change to other bus lines or transport modes. Therefore, the actual waiting time is objective and can be obtained by recording. The perceived and tolerable waiting times are subjective and can be obtained by questionnaire. 8

In short, the traffic at Guomao Bridge transport hub is abominable because of congestion, disorder and being out of control. One of urgent works for the hub is to optimize the timetable of the suburban bus routes.

4. A model to optimize the bus scheduling This section presents the procedures of the bus timetabling method. Firstly, after model assumptions in Section 4.1, the passenger satisfaction and the bus transit efficiency are formulated in Section 4.2. Secondly, a preliminary bus timetabling is developed in Section 4.3 and thirdly, the bus timetabling is optimized in Section 4.4. Finally, a solution algorithm is given in Section 4.5. 4.1. Model assumptions The model proposed in this study is based on the real traffic of Guomao Bridge transport hub with special geographical structure as shown in Fig. 2. It is about a work zone in one CBD of Beijing and a home zone in the suburban. Only one corridor connects the two zones. The map of the transportation system is sketched in Fig. 4. We only consider the suburban bus routes and other bus services are not be involved in the model. The proposed model is developed based on some dedicated suburban bus routes (Fig. 4). Each route can operate independently. The trip demand for each route is assumed to be constant and known. Arrivals of passengers at the bus stops are considered to be random. The waiting time that passengers are willing to bear is not long. For simplicity, the evening peak period is discretized into some equal intervals. In the case study, we discretize the period 17:00-19:30 into five equal intervals. The following notations are used in formulating the model: (a) Indices

i

bus line

j time interval, j  J k point of time in minute, k  1,

(b)

, K

Parameters

gi

vehicle size of line i

ij

passenger arrivals of bus line i within interval j

ijk

passenger arrivals of line i within k minute of interval j

9

t1 ( t2 )

critical values of waiting time

 min

minimal load factor (  min = 0.5 in this case study)

 max

maximal load factor (  max = 2 in this case study)

K length of an interval ( K  30 minutes in this case study) (c) Variables to be determined

i

load factor of bus line i

 ij

load factor of bus line i within interval j

zij

average load of line i within interval j

 ij

vehicle headway of line i within interval j

 ij

average headway of line i within interval j

nijF

number of buses dispatched with headway  ij   

nijC

number of buses dispatched with headway  ij   

nij

number of buses required for line i within interval j

Ni

total number of buses required for line i within the evening peak hours

wij

passenger waiting time of bus line i within interval j

sijpw

passenger satisfaction attributed to waiting time

sijpc

passenger satisfaction attributed to in-vehicle discomfort

sijp

satisfaction of passengers who choose bus line i within interval j

sip

average satisfaction of passengers who choose bus line i

eijb

efficiency of bus line i within interval j

eib

average efficiency of bus line i

4.2. Step 1: Passenger satisfaction and bus transit efficiency This subsection formulates the passenger satisfaction and the bus transit efficiency mathematically. 10

4.2.1. Passenger satisfaction Passenger satisfaction is hard to be quantified since there are too many quality-based attributes and the importance of attributes being negatively or positively perceived by individual passengers depends on their preferences. Perhaps a passenger satisfaction function can be mathematically expressed as a function of some explanatory variables such as waiting time, boarding time, seat availability, in-vehicle time, alighting time, total travel time, transfer time, and so on (Cats et al., 2016). However, we do not know which attribute is suitable and its importance. Our survey reveals that waiting time and in-vehicle bus comfort are the most important attributes. Then we present a model involving these two key aspects. These two attributes are both related to the load factor  ij which is measured by

ij 

zij

.

(1)

gi

Where zij denotes an average load and g i is the vehicle size. We can see that, if ij  1 , all passengers have seats and their requirement for on-board comfort can be ignored. If ij  1 , g i passengers have seats and expect to shorten their waiting times while (ij  1) gi

passengers have to experience in-vehicle crowding with no seats. pw pc Let sij and sij be the passenger satisfaction caused by waiting time and in-vehicle p bus comfort, respectively. Then, the passenger satisfaction sij is

sijp  (sijpw , sijpc ) A ,

(2)

 aiw  pw pc where A is the weight set of sij and sij . Let A   c  , here aiw and aic respectively  ai  pw pc denote the impacts of sij and sij . They are non-negative and satisfy

aiw  aic  1, aiw  0, aic  0 .

(3)

Paired comparison method is suitable here to handle the weights. If ij  1 , aiw  1 and aw

1

i a  0 holds. If ij  1, a c    1 holds, which represents that g i seating passengers i ij c i

pay attention to the waiting times while (ij  1) gi standees attach importance to the comfort. p Thus, sij can be computed by

11

 sijpw , if ij  1  p sij   1 pw ij  1 pc  sij   sij , if ij  1 ij  ij

(4)

pw pc Semi-trapezoidal distributions can be introduced to express sij and sij (Cats et al., 2016;

Ceder, 2016), as follows.

sijpw

1, wij  t1   t  wij  2 , t1  wij  t2  t1 0, wij  t2 

(5)

zij  gi 1,  s   max gi  zij  g  g , gi  zij   max gi i  max i

(6)

pc ij

where t1 and t2 are two critical values of waiting time. Eq. (5) states that it is satisfying if

wij (waiting time) is less than t1 , acceptable if wij is between t1 and t2 , and intolerable if wij is longer than t2 . Based on our survey data, t1 and t2 are assumed to be five and ten minutes, respectively. Eq. (6) shows that if there is no seats, the more the remaining space is, the more comfortable passengers feel. (max gi  gi ) denotes the maximum allowable standees and (max gi  zij ) represents the remaining space inside the vehicle. Then the weighted average satisfaction sip is

   J

p i

s

j 1 J

ij sijp

j 1

ij

.

(7)

Eqs. (2)-(7) comply with a certain level of service for that trip: waiting time, degree of comfort, seat availability, and other operational features. In this study, we propose a data-driven approach to obtain the waiting time wij in Eq. (5) (Mishalani et al., 2006; Turnquist and Bowman, 1980; Özekici, 1987; Ceder, 2016). This method is presented under the assumption that: (all actual waiting times) = (all assumed waiting times)-(all actual no-waiting times). At each time interval, the first bus departs at time

 ij1 . We assume  ij1   for ease of description. If all ij passengers reach the i th bus line at the start time of the j th interval, which is called Case 1, the waiting time for all ij 12

passengers can be obtained by: wij1  zij  ij1  zij  ( ij1   ij 2 )  +zij  ( ij1   ij 2 

 zij  ( ij1   ij 2 

  ijn 

  ijn ) 



  ijN )

(8)

 1 Where n  Z ,n [1, N ], N  nij , wij denotes the waiting time for all ij passengers under

Case 1. Eq. (8) indicates “all assumed waiting times” and can be illustrated in Fig. 6a. We can see that each term on the right-hand side of Eq. (8) corresponds to a rectangle’s area in Fig. 6a. 1 In other words, wij of Eq. (8) equals to the sum of N rectangles’ areas in Fig. 6a.

In real life, passengers arrive at the i th bus line one by one (Case 2). They may not wait at the start time of the j th interval. The waiting times before they arrive are assumed to be wij2 and can be obtained by wij2  ij ( 1) 1  ij ( 2)  2 

 ijK  ( K  )

(9)

So Eq. (9) indicates “all actual no-waiting times” and should be subtracted. Similarly, each 2 term on the right-hand side of Eq. (9) corresponds to a rectangle’s area in Fig. 6b and wij of

Eq. (9) equals to the sum of ( K  ) rectangles’ areas in Fig. 6b. It should be noted that

   ij1 and  is not a constant because of various j . Therefore, more accurate waiting time for ij passengers is ( wij1  wij2 ) and the average waiting time per passenger is wij 

wij1  wij2

(10)

ij

4.2.2. Bus transit efficiency All transit agencies are willing to improve the efficiency, productivity, and effectiveness of their systems. However, it is difficult to define the bus transit efficiency, which involves too many attributes including fleet size, empty-seat hours, length of routes, and number of transfers. Moreover, the optimization criteria and its calculation are cumbersome and complex. In this study, the operating benefit (ticket revenue) can only be found at the departure stop for a suburban bus route. Thus, greater load means more income. For simplicity, the bus transit efficiency eijb is defined as the ratio of traffic load over the maximum load, i.e.

eijb 

zij

 max gi

.

(11)

Since zij  giij from Eq. (1), we have

13

eijb 

ij gi   ij .  max gi  max

(12)

The average efficiency eib is given by

   J

b i

e

b j 1 ij ij J

ne

.

(13)

n

j 1 ij

4.3. Step 2: Preliminary bus timetabling This subsection introduces the way to determine the preliminary bus frequency and time headway. 4.3.1. Bus frequency determination The max load method has been widely used to determine the frequency of dispatching buses (Ceder, 2016). This method can ensure adequate space to accommodate the maximum number of on-board passengers. It is also used with a little revision in this study. Let  i be the bus load factor of line i . Suppose a bus company dispatches nij buses to transport ij passengers within interval j . Then, nij is computed by   ij  ij  Z   , if  g  g  i i nij    i i   ij , otherwise  g  i i

(14)

It should be noted that in many of the previous studies, two methods were adopted: both frequency and headway were rounded to the nearest integer; or, frequency were retained to be non-rounded but timetables were determined by smooth transition or some computation algorithms. Eq. (14) is different from these methods since nij is rounded off to the next highest integer if

ij is not an integer. In China, many transit agencies routinely round off i gi

the frequency nij to the next highest integer. The level of passenger comfort is increased at the expense of unnecessary operating costs. However, the “rounding-off” procedure is probably reasonable due to high load and overcrowding may be alleviated. When i  max , nij approaches the minimum. The average load zij can be obtained by

14

zij 

ij nij

,

(15)

and the total number of buses required for line i is J

Ni   nij .

(16)

j 1

4.3.2. Time headway Similarly, the “rounding-off” procedure is used here. Each time interval has K minutes. The even time headway  ij can be computed by

 ij =

K . nij

(17)

In reality, the smallest unit of a feasible bus timetable is integer in minute. If  ij is an integer, the bus company can schedule nij buses. However, the fractional  ij is inevitable and integer processing is necessary. A common headway smoothing rule is to use an average headway. However, it can result in either undesirable overcrowding or underutilization. Some headway-smoothing techniques are developed in the transition segments between adjacent time periods. For example, Ceder (2013) proposed an even-load timetable with smoothing transitions which move per load factor horizontally until intersecting the cumulative-load curve and then vertically. Later, they presented an even-load and even-headway smoothing techniques. These techniques are advanced. However, passenger demand in their studies is far less than that in China and the time headway is much longer. Ten-minute headway at peak hours may be usual in their studies but can cause riots of waiting crowds in China. Thus, a more practical but much simpler method is needed in this paper. C F Here,  ij    ij   ij  holds. Let nij and nij be the numbers of buses dispatched

with headways  ij  and  ij  , respectively. We have

nijF  ij   nijC  ij   K      F C nij  nij  nij

(18)

C F Solving Eq. (18), we can obtain nij and nij for all intervals of each bus line.

Thus, the bus headway can be obtained by

 ijn

   ij , if  ij  Z          ij  or  ij  , if  ij  Z

(19) 15

4.4. Step 3: Optimization of bus timetabling The purpose of this section is to propose a formulation for a trade-off between the passenger satisfaction and the bus efficiency for the determination of load factor. The transit agencies wish to utilize their resources more efficiently by minimizing the number of required vehicles and crew costs. Nevertheless, they also need to accommodate the observed passenger demand as well as possible. A cost-effective and efficient transit timetable embodies a balance between passenger satisfaction and service cost. Since it is hard to explicitly formulate the practical operational cost, an optimization model is developed with regard to the trade-off between the passenger satisfaction and the bus transit efficiency, subject to the constraints of load factor. In this way, the scheduler will be in a position not only to expedite manual/semi-manual tasks, but also to compare procedures and methods in regard to the trade-off between the passengers’ comfort and the operating efficiency. A balance is made here to reach “fewer buses, shorter waiting times, and less on-board discomfort” (Cats et al., 2016). This can be represented by maximizing the product of two dimensionless indexes, namely efficiency and satisfaction. We solve the following optimization problem for getting an optimal load factor  i* for each bus line: max eib sip

(20)

s.t .  min  i   max

Where, the objective function demonstrates the trade-off between the passenger satisfaction and the bus efficiency, and the constraint is to avoid undesirably low or high load. Note that the connection between the objective function and the constraint is not clear at p the first look in model (20), just as the models in Ceder (2016). Passenger satisfaction si is b computed by Eqs. (1)-(10) and bus transit efficiency ei is computed by Eqs. (11)-(13). Load

factor  i is an input variable. Via the model (20), the optimal load factor  i* can be obtained by a numerical search approach in the next subsection. Once  i* is obtained, we can * * determine the optimal nij from Eq. (14), then the optimal  ij from Eq. (17), the optimal

nijF * and nijC * from Eq. (18), and finally, the optimal bus timetable.

Specially, we should point out that the timetabling only considers the arrival at the original stop in CBD, as shown in Fig. 4. The arrivals at downstream stops are not considered. Thus, the load factor and waiting times are only for the passengers at the very first stop of the bus line. Thus, this is not the definition of “average load factor” used in most papers. 16

4.5. Solution algorithm The proposed maximization model is a nonlinear programming model. A numerical search approach was used, as shown in Fig. 7. It is calculated by a few simple mathematical procedures of Matlab R2012a and implemented on a personal computer with Intel(R) Core(TM)2 Duo CPU P8600 @2.40GHz 2.40GHz and 4.00G RAM. Fig. 7 represents all the steps of our algorithm.  ij can be obtained from Eq. (1), and wij p can be obtained from Eqs. (8)-(10). Then passenger satisfaction si can be obtained from Eqs. b (2)-(7). Meanwhile, bus transit efficiency ei is computed by Eqs. (11)-(13). If we let

i  i  0.01% and substitute sip , eib into the model (20), the optimal eib sip and the * * optimal  i* can be calculated. With  i* , the optimal nij from Eq. (14), the optimal  ij F* C* from Eq. (17), the optimal nij , nij from Eq. (18), and the optimal bus timetable are all

determined. It should be noted that we cannot guarantee the global or local optimum of the solutions, since we let i  i  0.01% in the iteration for simplicity and the continuity of the target function is unknown. The model in this paper is simple and the method can be applicable. However, more research should be conducted in future.

5. Model implementation and validation All data are taken from the survey described in Section 2. We select four suburban bus lines 813, 807, 817, and 818 whose average waiting times are over ten minutes. All vehicles of these four lines have the same size, i.e., g813  g807  g 818  g 817  59 . Fig. 8 shows the five-minute passenger arrivals. Obviously, the passenger-arrival rates are not uniform. The buses should not be dispatched with fixed headway for matching with the time-varying passenger arrivals. The algorithm in Section 4.5 is used to solve the model (see Fig. 7). The optimal load factors for these four lines are obtained.

17

* 813  * 807  * 817  *  818

 145.25%  145.48%

(21)

 67.05%  75.83%

As shown in Eq. (21), we can see that standees are inevitable with high demands (Line 813, Line 807) and shorter waiting times are expected with low demands (Line 817, Line 818). And we obtain the following relationships between bus transit efficiency eib and load factor  i , i  813, 807, 817, and 818. 2 b e813  0.051  0.580813  0.051813   2 b e807  0.192  0.772807  0.111807   b 2 e817  0.012  0.489817  0.042 817   b 2 e818  0.072  0.366818  0.005 818 

(22)

The relationships between the passenger satisfaction sip and load factor  i are as follows: 2 p  s813  1.000  0.316813  0.292 813   2 p  s807  0.851  0.535807  0.371807   p 2  s817  1.427  1.402817  0.383 817   p 2  s818  0.740  0.627818  0.163 818 

(23)

Fig. 9 depicts the functions (22) and (23) against load factor. It can be seen that the bus efficiency eib basically increases with  i but the passenger satisfaction sip decreases with

 i . These two attributes are indeed incompatible. Transit agencies wish to utilize their resources more efficiently by minimizing the number of required vehicles and crew costs. The higher the load factor is, the higher the bus efficiency is. Nonetheless, higher load factor means poorer on-board bus comfort and lower passenger satisfaction. Moreover, from Fig. 9, we can find that not all the meeting points appear at i  1 . There is an intersection point between the two lines of the bus efficiency eib and the passenger satisfaction sip for Line 813 and Line 807. However, there is no intersection point for Line 18

817 and Line 818. The results relate to the passenger demands in Fig. 8. Substituting Eq. (21) into Eqs. (14)-(19), we obtain the optimal bus timetables of four lines which are depicted in Fig. 10. Passenger arrival rates are not uniform and the headways are not fully even. The beginning and end of the schedule horizon are set to be at 17:00 and 19:30. The left-side of the long vertical line is time headways and the right-side of the line is departure times. For instance, to comply with the surveyed demands in Fig. 8, Line 813 will dispatch 41 vehicles. The first departure time after 17:00 will be 17:04. Therefore, the complete timetable is: 17:04, 17:08, …, 19:25, and 19:30. Similarly, Lines 807, 817, and 818 will dispatch 41, 38, and 30 vehicles. Table 2 compares the existing (surveyed) and optimal bus timetables of the four lines. The average waiting times are from the surveyed data and the optimum ones are obtained according to the algorithm in Fig.7. From Table 2, we can find: (1) An interesting observation is that all the bus transit efficiency scores are reduced for the optimal load factor (the 4th row), which shows that the existing bus timetables pay more attention to the bus efficiency (the 4th row) than the passenger satisfaction (the 5th row). After optimization, all the efficiencies decrease but the satisfactions greatly increase, consequently the products improve (the 6th row). (2) Existing supply can meet demand for Line 813 and Line 807 (the 1st row). However, uneven headways (see Fig. 5) lead to low satisfaction with long waiting time (the 2nd row) and poor bus comfort (the 3th row). Thus, only one more bus is required to enhance the satisfaction by 472.0% for Line 813 and one less bus is required to promote the satisfaction by 435.4% for Line 807. Due to the optimal headways, the waiting times are significantly shortened and the in-vehicle comforts are greatly improved. (3) The buses for Line 817 and Line 818 increase since the existing supply cannot meet demand (the 1st row). Insufficient vehicles lead to long waiting times (the 2nd row). More vehicles are suggested to be provided. This will decrease the efficiency but increase the satisfaction much more.


6. Conclusions Public transport priority has become an important policy in China. However, long waiting time and poor in-vehicle comfort make passengers very dissatisfied with the existing bus systems. Previous studies mostly aimed to get fewer buses and shorter waiting time, seldom addressing the in-vehicle comfort. The objective of this paper is to find a bus timetabling 19

method which is practical and applicable in China. We started from a survey on Guomao Bridge transport hub in Beijing and highlighted passenger satisfaction in developing a new bus timetabling model. The passenger satisfaction and the bus transit efficiency are firstly formulated. For simplicity, waiting time and bus comfort are considered to be the main attributes of the passenger satisfaction according to our survey data. Waiting time and on-board comfort are two basic elements considered and their semi-trapezoidal distributions are introduced to formulate the passenger satisfaction. The paired comparison method is adopted for determining the weights. Secondly, a preliminary bus timetabling is developed. The max load method is used to determine the frequencies and even headways with smooth transitions are proposed. Thirdly, an optimization model is developed with regard to the trade-off between the passenger satisfaction and the bus transit efficiency, subject to the constraints of load factor. A numerical search approach is used. Finally, this method is verified by field data and the results demonstrate that our model can greatly reduce passengers’ waiting times and increase their in-vehicle comforts. Our model shows that even fewer vehicles can greatly promote passenger satisfaction for those lines whose supply can meet demand. It should be noted that the model in this paper is simple and only considers one directional (outbound) demand during the afternoon peak periods. The number of buses needed for a line should consider the round travel time and dwelling time at both ends. Moreover, if all the bus lines are operated by a single company, the total number of buses available could be an important constraint. Another important factor is the number of seats in a vehicle. The vehicle seat capacity could be another decision variable for the timetabling design. Our next work is to extend the model to cover more factors associated with passengers and transit agency. In addition, optimizing the bus crew scheduling and bus stop placement are also urgent to Beijing and worthy of studying deeply and extensively.

Acknowledgements This research was jointly supported by grants from the National Natural Science Foundation of China (71890971, 71371128, 71271004) and the Key Program for Social Science Fund of Beijing Municipal Education Commission.

References [1] 13th Five-Year Plan for National Economy and Social Development, 2016–2020. Chinese version. Retrieved from: 20

http://www.stdaily.com/zhuanti01/guojia/2018-01/04/content_617866.shtml [2] Abkowitz, M. D., & Engelstein, I. (1983). Factors affecting running time on transit routes. Transportation Research Part A, 17(2), 107-113. [3] Abkowitz, M. D., & Engelstein, I. (1984). Methods for maintaining transit service regularity. Transportation Research Record, 961, 1-8. [4] Badami, M. G., & Haider, M. (2007). An analysis of public bus transit performance in Indian cities. Transportation Research Part A, 41(10), 961-981. [5] Barnett, A. (1974). On controlling randomness in transit operations. Transportation Science, 8(2), 102-116. [6] Barnett, A., & Kleitman, D. J. (1973). Optimal scheduling policies for some simple transportation systems. Transportation Science, 7(1), 85-99. [7] Cats, O., West, J., & Eliasson, J. (2016). A dynamic stochastic model for evaluating congestion and crowding effects in transit systems. Transportation Research Part B, 89, 43-57. [8] Ceder, A. (2016). Public transit planning and operation: modeling, practice and behavior. CRC press. [9] Ceder, A., Hassold, S., & Dano B. (2013).Approaching even-load and even-headway transit timetables using different bus sizes. Public Transport, 5(3), 193-217. [10] Ei-Geneidy, A. M., Strathman, J. G., Kimpel T. J, & Crout, D. T. (2006). Effects of bus stop consolidation on passenger activity and transit operations. Transportation Research Record, 1971(1): 32-41. [11] Friman, M., & Fellesson, M. (2009). Service supply and customer satisfaction in public transportation: The quality paradox. Journal of Public Transportation, 12(4): 57-69. [12] Fonseca, J. P., Van Der Hurk E., Roberti R., & Larsena A. (2018). A matheuristic for transfer synchronization through integrated timetabling and vehicle scheduling. Transportation Research Part B, 109: 128-149. [13] Furth, P. G., & Wilson, N. H. M. (1981). Setting frequencies on bus routes: Theory and practice. Transportation Research Record, 818, 1-7 [14] He, S. X.(2015). An anti-bunching strategy to improve bus schedule and headway reliability by making use of the available accurate information. Computers & Industrial Engineering, 85:17-32 [15] Huang, H. J., Tian, Q., Yang, H., & Gao, Z. Y. (2007). Modal split and commuting pattern on a bottleneck-constrained highway. Transportation Research Part E, 43(5): 578-590. [16] Jiang, X., Guo, X., & Ran, B. (2014). Optimization model for headway of a suburban bus route. Mathematical Problems in Engineering, 979062, 1-6. [17] Khasnabis, S., Rudraraju, R. K., & Baig, M. F. (1999). Economic evaluation of signal preemption projects. Journal of transportation engineering, 125(2), 160-167. [18] Kurauchi, F., Bell, M. G. H., & Schmöcker, J-D. (2003). Capacity constrained transit 21

assignment with common lines. Journal of Mathematical Modelling and Algorithms, 2, 309–327. [19] Lai, W. T., & Chen, C. F. (2011). Behavioral intentions of public transit passengers—The roles of service quality, perceived value, satisfaction and involvement. Transport Policy, 18(2): 318-325. [20] Leurent, F., Chandakas, E., & Poulhes, A. (2014). A traffic assignment model for transit on a capacitated network: bi-layer framework, line sub-models and large-scale application. Transportation Research Part C, 47, 3–27. [21] Liu T, Ceder A. (2017). Integrated public transport timetable synchronization and vehicle scheduling with demand assignment: A bi-objective bi-level model using deficit function approach. Transportation Research Part B, 23, 341-361. [22] Michaelis, M., & Schöbel, A. (2009). Integrating line planning, timetabling, and vehicle scheduling: a customer-oriented heuristic. Public Transport, 1(3), 211-232. [23] Mishalani, R. G., McCord, M. M., & Wirtz, J. ( 2006). Passenger wait time perceptions at bus stops: empirical results and impact on evaluating Real-Time Bus Arrival Information. Journal of Public Transportation, 9(2), 89-106. [24] Nuzzolo, A., Crisalli, U., & Rosati, L. (2012). Schedule-based assignment models with explicit capacity constraints for congested public transport networks. Transportation Research Part C, 20, 16–33. [25] Osuna, E. E., & Newell, G. F. (1972). Control strategies for an idealized public transportation system. Transportation Science, 6(1), 52-72. [26] Otto, B., & Boysen, N. ( 2014). A dynamic programming based heuristic for locating stops in public transportation networks. Computers & Industrial Engineering, 78, 163-174. [27] Özekici, S. (1987). Average waiting time in queues with scheduled batch services. Transportation science, 21(1): 55-61. [28] Pel, A., Bel, N. H., & Pieters, M. (2014). Including passengers’ response to crowding in the Dutch national train passenger assignment model. Transportation Research Part A, 66, 111–126. [29] Prud’homme, R., Koning, M., Lenormand, L., & Fehr, A. (2012). Public transport congestion costs: the case of Paris subway. Transport Policy, 21, 101–109. [30] Rajbhandari, R., Chien S. I., & Daniel, J. R. (2003). Estimation of bus dwell times with automatic passenger counter information. Transportation Research Record, 1841(1), 120-127. [31] Schmöcker, J-D. , Bell, M. G. H., & Kurauchi, F. (2008). A quasi-dynamic capacity constrained frequency-based transit assignment. Transportation Research Part B, 42 (10), 925–945. 22

[32] Schmöcker, J-D., Fonzone, A., Shimamoto, H., Kurauchi, F., & Bell, M.G.H. (2011).

Frequency-based transit assignment considering seat capacities. Transportation Research Part B, 45, 392–408. [33] Shang, H. Y., Huang, H. J., & Xiao, L. L. (2014). Investigating a transit-oriented transport

hub in Beijing. Journal of Transportation Systems Engineering and Information Technology, 14(3), 201-206+220. [34] Statistical Yearbook, 2017. Beijing Statistical Yearbook, 2017, Beijing, China. Chinese version. Retrieved from: http://www.bjstats.gov.cn/nj/main/2017-tjnj/zk/indexch.htm [35] Strathman,J. G., Kimpel, T. J., Dueker, K. J., Gerhart, R. L., & Callas, S. (2002). Evaluation of transit operations: data applications of Tri-Met's automated Bus Dispatching System. Transportation, 29(3), 321-345. [36] Sumalee, A., Tan, Z., & Lam, W.H.K. (2009). Dynamic stochastic transit assignment with explicit seat allocation. Transportation Research Part B, 43, 895–912. [37] Tian, Q., Huang, H. J., & Yang, H. (2007). Equilibrium properties of the morning peak-period commuting in a many-to-one mass transit system. Transportation Research Part B, 41(6): 616-631. [38] Tirachini, A., Hensher, D. A., & Rose, J. M. (2013). Crowding in public transport systems: effects on users, operation and implications for the estimation of demand. Transportation Research Part A, 53, 36–52 . [39] Tirachini, A., Hensher, D. A., & Rose, J. M. (2014). Multimodal pricing and optimal design of urban public transport: The interplay between traffic congestion and bus crowding. Transportation Research Part B, 61, 33-54 [40] Toledo, T., Cats, O., Burghout, W., & Koutsopoulos Haris, N. (2010). Mesoscopic simulation for transit operations. Transportation Research Part C, 18(6): 896-908. [41] Trozzi, V., Gentile, G., Bell, M., & Kaparias, I. (2013). Dynamic user equilibrium in public transport networks with passenger congestion and hyperpaths. Transportation Research Part B, 57, 266–285. [42] Turnquist, M. A., & Bowman, L. A. (1980). The effects of network structure on reliability of transit service. Transportation Research Part B, 14(1), 79-86. [43] Wilson, N. H. M., Macchi, R. A., Fellows, R. E., & Deckoff, A. A. (1992). Improving service on the MBTA Green Line through better operations control. Transportation Research Record, 1361, 10-15.

23

Table 1 Waiting times (minutes) for riding in buses. Line 813

Line 807

Line 817

Line 818

Actual

Perceived

Tolerable

Actual

Perceived

Tolerable

Actual

Perceived

Tolerable

Actual

Perceived

Tolerable

Sample size

878

874

561

372

374

374

342

338

333

476

476

475

Mean value

12.33

22.35

22.95

9.62

23.09

29.26

12.83

15.43

30.26

20.72

23.82

30.64

Std. deviation

7.67

10.46

10.71

6.76

11.73

15.42

6.77

7.50

15.79

11.58

12.17

16.11

Minimum

0.00

3.00

1.00

0.00

2.00

1.00

0.00

1.00

10.00

0.00

1.00

5.00

Maximum

42.00

60.00

60.00

38.00

60.00

90.00

35.00

30.00

60.00

47.00

60.00

60.00

24

Table 2 Optimal and existing (surveyed) bus schedules of four lines. Line 813

Line 807

Line 817

Line 818

Surveyed

Optimal

40

41



42

41



24

38



20

30



Average waiting time

12.3292

4.5415



9.6210

4.4307



12.8275

6.4369



20.7185

6.9336



Average in-vehicle comfort

0.4398

0.5962



0.5541

0.5930



0.7377

1.0000



0.6976

1.0000



0.7801

0.7003



0.7205

0.7040



0.6090

0.3196



0.6088

0.3452



Passenger satisfaction ( si )

0.1476

0.8443



0.1609

0.8614



0.1107

0.7369



0.1055

0.6373



eib sip

0.1151

0.5913



0.1159

0.6064



0.0674

0.2355



0.0642

0.2200



No. of vehicles required

b

Bus transit efficiency ( ei ) p

Rise/Fall Surveyed Optimal Rise/Fall Surveyed Optimal Rise/Fall Surveyed Optimal Rise/Fall

25

Fig. 1. Crowding at the bus stop. (a) People queuing for Line 813, the upper is queues inside the barrier and the lower is queues outside. (b) People climbing windows. The photos are taken from our survey.

26

Fig. 2. Spatial distribution of the corridor. The picture is taken from the website: http://map.baidu.com/. Two yellow circles are used to highlight Guomao Bridge and Yanjiao Town. Only one expressway connects them.

27

Fig. 3. Surveyed zone with an area size of 1.1 square kilometers around Guomao Bridge. The picture is taken from the website: http://map.baidu.com/.

28

Fig. 4. Sketch map of the suburban bus routes between work zone and home zone.

29

Fig. 5. Bus headways of Line 807.

30

zij zij

Bus

zij

 ij1  ij 2  ij 3

Headway

… …

zij

zij

 ijN

ij ( 1) ij (  2) ij ( 3) Passenger Time



ijK … … K-1 K

 +1+2

ijK ij ( K 1)

zij zij zij





zij zij

 ij1  ij 2  ij 3 0

1

2

ij ( 3) ij (  2) ij ( 1)

 ijN

3 … N-1 (a)

N Schedule

1 1 1  +1+2

1

… K-1 K Time (b)

Fig.6. Waiting times: (a) is for Eq.(8) and (b) is for Eq.(9).

31

Fig. 7. Flow chart of the solution algorithm

32

Fig. 8. Passenger arrival rates (per 5 minutes) of four lines.

33

Fig. 9. Passenger satisfaction/bus transit efficiency against different load factors.

34

Fig. 10. Optimal bus timetables of four lines originated from Guomao Bridge transport hub.

35

Highlights:  Conduct an empirical study in Beijing.  Formulate passenger satisfaction addressing waiting time and on-board discomfort. 

Propose a practical timetabling with a compromise of efficiency and satisfaction.

36