# An Operation Plan Optimization Model for Integrated Customized and Conventional Bus Services Based on Cost Analysis

Authors: Xi Lu; Junru Han; Peng Xu; Haode Liu; Xumei Chen
Journal: Journal of Advanced Transportation
 (2022)

---

## Abstract

A proper public transportation system should be economically efficient to improve its sustainability and competitiveness. The superfluous provision of conventional bus services in areas with a weak or scattered travel demand distribution might result in a waste of transport capacity. This study proposes an operation plan optimization model based on the characteristics and cost composition of conventional and customized bus services, which systematically analyzes the feasibility of fully or partially replacing conventional bus services with customized bus services. A mixed-integer planning model and a hybrid genetic algorithm are used to optimize the operation plan for specific travel demands and supply capacities. Then, a sensitivity analysis of key factors is carried out to identify the optimal operation mode under different circumstances. The results show that with an optimized fleet size, service frequency, and service area, integrated operation modes can reduce operating costs without sacrificing the passenger time.


---

## Body

## 1. Introduction

In recent years, the urbanization process of major economies worldwide has continued to develop rapidly, as has the scale of cities. Furthermore, the urban traffic structure has also undergone significant changes; motorized travel has escalated rapidly, and urban traffic congestion has become increasingly serious. These factors have put tremendous pressure on the urban transportation system. Public transportation (PT) has been proven as an effective way to alleviate urban traffic pressure and respond to the development of urbanization. Therefore, many countries have put forward the clear strategy of giving priority to the development of public transport as an important path to the alleviation of traffic congestion. For example, China has proposed a binding national development goal of 100% coverage of public transport stations in built-up areas in cities with a population of more than 1 million [1].A proper PT system should not only increase the attractiveness of PT and promote the coordinated development of urban and rural regional transportation but should also be economically efficient to improve its competitiveness. Due to the basic public service attributes of PT, the PT management department usually opens new bus lines or extends the original bus lines in remote areas of the city, even in sparsely inhabited areas, to increase PT coverage. However, the provision of fixed PT services in areas with low travel demand and scattered distribution would inevitably lead to a waste of supply capacity to ensure that the fixed departure interval is within a proper range. Thus, customized public transport is one of the potential ways to solve this issue.Customized bus services are one form of demand-responsive transport (DRT), which also exists as flexible transport services (FTS), dial-a-ride services, or paratransit services; these services represent flexible bus services for which the travel needs of passengers are investigated in advance, and the bus route is then optimized according to a set of time windows and capacity restrictions [2]. Some studies in the field of customized PT have explored the cost difference between conventional and variable PT methods. These studies agree that, in some cases, the replacement of conventional PT with variable PT can achieve the goals of not affecting the use and reducing costs [3, 4]. However, few studies have systematically analyzed the feasibility of the partial replacement of conventional bus routes with customized bus services in different operating environments.This study aims to fill this gap via the use of a mixed-integer planning model with an improved genetic algorithm, and the following is achieved. (1) A bus system operating cost model based on the operating characteristics and cost element composition of conventional and customized bus services is established. (2) An operation plan selection and optimization model, as well as a solution algorithm, is proposed for the identification of an appropriate operation plan to use customized buses to fully or partially replace conventional bus lines for specific travel demands and supply capacities. (3) A sensitivity analysis is carried out, the best operation plan for corresponding travel demand is determined, and the potentially applicable environment in which a particular operation mode has advantages is identified.The remainder of this paper is organized as follows. A literature review is presented in Section2, and the proposed model is presented in Section 3. Then, the solution algorithm of the problem is detailed in Section 4. The details of numerical experiments and the corresponding results are presented and discussed in Section 5. Finally, Section 6 concludes this article.


## 2. Literature Review

With the continuous deepening of PT research, innovative flexible bus transit services, such as customized buses, community buses, circular routes, and feeder services, have emerged. Kirby and Bhatt [5] analyzed and identified the characteristics of subscription bus services, an original type of customized bus service, such as the increased focus on long-distance trips, high service reliability, and satisfactory personalized features. Daganzo [6] conducted a preliminary study on dial-a-ride services, and the results showed that as the demand level decreases, the demand-responsive system becomes relatively more attractive than the fixed-route system, and the checkpoint system may become more cost-effective. Khattak and Yim [7] explored the demand for a personalized DRT service via a survey conducted in the San Francisco Bay Area, which revealed the potentially high demand of people willing to pay for the service. Davison et al. [8] explored factors that affect the competitiveness of the demand-responsive transit systems in Great Britain as well as their underlying mechanisms. Wang et al. [9] conducted a comprehensive analysis of a successful customized bus system based on practical subscription data. Later, Wang et al. [10, 11] analyzed the mechanism of customized bus user loyalty modeling as well as the spatial dependence and spillover effects in customized bus demand; the results demonstrated that areas with poor accessibility have a larger market for customized bus service.Some studies have compared flexible bus services with regular bus services with fixed routes and have explored the feasibility of replacing the fixed-route policy with the flexible-route policy of a branch line system without disturbing the existing coordination between the main transit and branch line services [2, 12, 13]. With the goal of minimizing the average cost (operating and passenger costs), Chien et al. [14] concluded that operators can reduce daily operating costs by providing conventional services during periods of high demand and operating customized services during off-peak periods. Li and Quadrifoglio [15] explored the operation model of switching between a demand-responsive and a fixed-route operating policy and concluded that the demand-responsive policy is more preferred when the demand falls between 10 and 50 passengers/mile2/h. Velaga et al. [16] explored public transport provisions in rural and remote areas illustrated with experience from Scotland and demonstrated that flexible transportation services provide considerable potential for PT systems in rural areas. Nourbakhsh and Ouyang [17] proposed an alternative flexible-route bus system within a predetermined area by seeking the best network layout, the service area of each bus, and the headway of the bus to minimize the total cost of the system; they found that the proposed flexible-path system tends to have the lowest system cost under low- to medium-demand levels. Qiu et al. [18] proposed a model with which to choose between fixed-route and flexible-route policies for transit services based on the assessment of the service quality function under varying demand. Ronald et al. [19] studied the performance changes of a fixed-time flexible-route scheme and a fully self-organizing scheme using MATSIM and demonstrated that the two programs have different effects on operators and passengers.Although high-quality research has been conducted in recent decades to optimize the timetabling and vehicle scheduling of conventional and customized bus services, it is difficult to formulate a cohesive solution for the planning problem of the entire bus network via sequential approaches. Kim and Schonfeld [20, 21] employed optimization models to analyze conventional-only, flexible-only, and variable-type (can switch between conventional and flexible services as demand changes over time) services by jointly considering the bus size and service type; the results confirmed that variable-type operation can reduce the system cost by changing the operation type as the demand density changes. Later, Kim and Schonfeld proposed probabilistic optimization models to integrate and coordinate original and flexible bus services to connect one terminal and multiple local regions [22]. Liu et al. [23] compared the constant-frequency conventional bus dispatching mode with the dynamic-frequency demand-responsive bus dispatching mode; they then proposed a hybrid operating mode that combines fixed- and dynamic-frequency services to minimize the average waiting time of passengers and maximize the operator’s profit. In summary, most previous studies have separately focused on the optimization of the schedule or fleets of only conventional bus services, customized bus services, or variable-type bus services; research on the joint optimization of conventional buses, customized feeder buses, and customized point-to-point buses is comparatively lacking. Therefore, the methods for the optimization of integrated customized and conventional bus services require further exploration.


## 3. Method

Assume that a bus line provides travel services between the city center and the suburbs. This line that has been built or is about to be constructed hasn bus stations with known departure intervals and operating service times. Moreover, the travel demands are known and stable.Table1 lists a portion of the sets, indices, and parameters used in this study.Table 1
The sets, indices, and parameters.
SymbolDefinitionPSet of passengersUkSet of customized busk route sectionsKconSet of conventional busesKcusSet of customized busesQSet of passenger origins in customized busESet of intervalse between adjacent stations in fixed routePIndex of passengerskIndex of vehiclei,jIndices of stop locationsi,jkIndices of route section of customized busk from point i to point jeIndex of interval in fixed routeaIndex of operation planωIndex of operation modeNNumber of original conventional bus stationsMNumber of stations replaced


### 3.1. Operation Mode

#### 3.1.1. Conventional Operation Mode

The conventional operation mode (Figure1) is the most widely used operation mode of conventional bus systems, in which buses always run on urban roads according to fixed routes and fixed schedules. Passengers walk to the appropriate station from various origins, take the bus, get off at the appropriate station, and then walk to their destinations.Figure 1
The conventional operation mode.


#### 3.1.2. Integrated Operation Mode

In the integrated operation mode, PT services are jointly completed by two types of buses, namely, conventional and customized buses. Specifically, in low-demand districts, customized bus services are used to replace portions of the original conventional bus services, and the remaining stations are still served by conventional buses. In this study, the integrated operation mode is divided into two types, namely, (1) integrated operation mode I: customized buses (feeder transit) and conventional buses and (2) integrated operation mode II: point-to-point (p2p) customized buses and conventional buses.(1) Integrated Operation Mode I: Customized Buses (Feeder Transit) and Conventional Buses. In this mode (Figure 2), a customized bus provides a connection service in its service area, picks up passengers according to the needs in the area, and transfers them to the transfer station; a conventional bus completes the rest of the transportation service. Specifically, a customized bus starts from the outer terminal station of the existing conventional bus route, picks up passengers according to their specific location within the m stations, and transfers them back to the transfer station to complete connection services. Then, the entire service is completed by the conventional bus covering the remaining n−m stations.Figure 2
Integrated operation mode I customized buses (feeder transit) and conventional buses.(2) Integrated Operation Mode II: Customized Buses (p2p) and Conventional Buses. In this mode (Figure 3), customized buses transfer some passengers from their origin to their destination, which is a type of “direct” public transport service. Passengers are no longer required to go to the departure bus station or get off at the arrival bus station. In other words, the p2p customized bus departs from the depot, picks up passengers according to their specific location within the m stations, transfers them to their destination, and finally returns to the depot. Other passengers continue to be served by a conventional bus.Figure 3
Integrated operation mode II: customized buses (p2p) and conventional buses.


### 3.2. Cost Element

The total operating cost of a bus service includes the operating costs of both the conventional and customized buses, i.e.,(1)Ctotal=Ccon+Ccus,where Ccon is the conventional bus operating cost, Ccus is the customized bus operating cost, and Ctotal is the total operating cost. The conventional bus operating cost includes the fixed input cost Ccon,f and variable operating cost Ccon,m, and the same is true for the customized bus operating cost.


#### 3.2.1. Fixed Input Cost

The fixed input cost is the sum of the fixed cost of each vehicle in the fleet. The fixed cost is generally related to the size of the fleet, the operating time, the drivers’ salaries, vehicle purchase fees, vehicle depreciation, insurance fees, etc. Formally, the fixed input costs of conventional busesCcon,f and customized buses Ccus,f can, respectively, be expressed as(2)Ccon⋅f=αcon⋅f×Gcon×T,Ccus⋅f=αcus⋅f×Gcus×T,where αcon,f is the fixed input cost per conventional bus hour, Gcon is the size of the conventional bus fleet, αcus,f is the fixed input cost per customized bus hour, Gcus is the size of the customized bus fleet, and T is the length of the bus service time. Gcon can be calculated based on the departure interval hcon by equation (4), where vbus is the average driving speed of the bus during the bus service time, and L is the total mileage of one cycle. Moreover, Gcus is related to the number of customized bus departures and each departure time.(3)Gcon=Lhcon×vbus.


#### 3.2.2. Variable Operating Cost

The variable operating cost is generally related to the fuel cost, which is associated with the total mileage of the fleet. The variable operating costs of conventional busesCcon,m and customized buses Ccus,m can, respectively, be expressed as(4)Ccon,m=αcon,m×∑k∈KconLk,Ccus,m=αcus,m×∑k∈Kcus∑i,j∈UkIi,jk,where αcon⋅m is the fuel cost of conventional buses per mile, Lk is the total mileage of conventional bus k, Kcon is the set of conventional buses, αcus⋅m is the fuel cost of customized buses per mile, Kcus is the set of customized buses, Ii,jk is the mileage of customized bus k driving from point i to point j, and Uk is the set of route sections of customized bus k. Thus, ∑i,j∈UkIi,jk is the total mileage of customized bus k.


### 3.3. Assumption

As it would be difficult to determine the specific departure and destination locations of passengers during the operation of the conventional mode, certain assumptions are made regarding the characteristics of passenger travel demand based on actual conditions:(1)
The operation mode has no effect on bus travel demand, which means that conventional bus travel demand would be fully transferred by customized bus services in areas where the conventional operation mode has been changed.(2)
The location and travel time of passengers are randomly generated based on an assumed distribution.(3)
Customized bus provides pickup and drop-off service; therefore, the passenger’s walking time is ignored.


### 3.4. Objective Formulation

The research goal is to determine the optimal bus operation plan and the related operating parameters after clarifying the corresponding travel needs. Therefore, 0-1 integer programming was used to construct the objective function of the operation mode selection and optimization model. The objective function is divided into three parts, and its expression is(5)minC=∑ω=02Zω×Ctotal,ω,where ω represents the bus operation mode. When ω=0, Ctotal,0 is the total operating cost of the conventional operation mode; when ω=1, Ctotal,1 is the total operating cost of integrated operation mode I; and when ω=2, Ctotal,2 is the total operating cost of integrated operation mode II. Moreover, C is the operating cost of the optimal bus operation plan, and zω is a binary variable.As the conventional operation mode does not include customized buses, the operating cost of customized buses is considered as 0, and the expression ofCtotal,0 is(6)Ctotal,0=αcon,f×Gcon,0×T+αcon,m×∑k∈Kcon,0Lk,where Gcon,0 is the size of the conventional bus fleet in the conventional operation mode and Kcon,0 is the set of conventional buses in the conventional operation mode.The operating costs of the integrated operation modes, namely,Ctotal,1 and Ctotal,2, include two parts, namely, the operating cost of conventional bus operation and the operating cost of customized bus operation. The corresponding expressions are as follows:(7)Ctotal,1=αcon,f×Gcon,1×T+αcus,f×Gcus,1×T+αcon,m×∑k∈Kcon,1Lk∑k∈Kcus,1∑i,j∈UkIi,jk,Ctotal,2=αcon,f×Gcon,2×T+αcus,f×Gcus,2×T+αcon,m×∑k∈Kcon,2Lk∑k∈Kcus,2∑i,j∈UkIi,jk.


### 3.5. Constraints

Because the goal of the model is to reduce the operating costs while ensuring that the needs of passengers are met, the main constraints of this model are given as follows:(8)z0+z1+z2=1,(9)loadi,jk≤Capk∀i,j∈k,i≠j,∀k∈Kcus,(10)T¯user,a≤T¯user,0,(11)1≤m<n,m∈N,(12)∑k∈Kcus∑j∈Qxi,jk=1,∀i∈Q,i≠j,(13)∑k∈Kcus∑j∈Qxj,ik=1,∀i∈Q,i≠j,(14)∑j∈Qxi,jk−∑j∈Qxj,ik=0,∀k∈Kcus,i≠j.(8) is the mode selection constraint; there is only one variable, which corresponds to the operation mode of the optimal bus operation plan. (9) is the vehicle load constraint, which stipulates that the number of served passengers in any route section of customized bus k cannot exceed its capacity. (10) is the average passenger time constraint, which guarantees that the average passenger time T¯user,a of plan a in the integrated operation mode does not exceed the average passenger time T¯user,0 in the conventional operation mode. (11) restricts the number of replacement stations m to not exceed the total number of original conventional bus stations n. (12) and (13) indicate that for each passenger who uses a customized bus, there is only one customized bus to provide service, where Q is the set of origins of passengers using a customized bus. Moreover, xi,jk=1 represents that customized bus k travels from point i to point j; otherwise, xi,jk=0. Finally, (14) indicates that the customized bus that arrives at point i and departs from point i must be the same bus.


#### 3.5.1. Total Passenger TimeTuser

(1) Total Passenger Time. Tuser,0in the Conventional Operation Mode. The total passenger time on a fixed route is Tfixed, which includes the passenger walking time Tcon,wk, the passenger waiting time Tcon,wt, and the passenger travel time Tcon,tl.(15)Tfixed=Tcon,wk+Tcon,wt+Tcon,tl,where Tcon,wk is the passenger walking time for the conventional bus mode, which is equal to the total passenger walking distance divided by the walking speed, i.e., Tcon,wk=Dcon/vpassenger, and Dcon is the summation of the walking distance of all passengers P from the origin to the departure station and from the arrival station to the destination, namely, Dcon=∑p∈PDu,p+Dd,p.Tcon,wt is the passenger waiting time for conventional buses. Assuming that the time for passengers to arrive at the departure stations is uniformly distributed, Tcon,wt can be calculated from the departure interval hcon, namely, Tcon,wt=Ncon,wt×1/2hcon, where Ncon,wt is the number of passengers waiting for conventional buses.Tcon,tl is the passenger travel time on conventional buses. The total travel time of passengers is the sum of the in-vehicle time of all passengers, which is equal to the multiplication of the number of passengers in each section by the section travel time. Hence, the calculation of Tcon,tl is converted into the calculation of the sum of the travel time of passengers in each interval between adjacent stations, where E represents the set of intervals e between adjacent stations on the fixed route, rk,e is the number of people within conventional bus k in interval e, and tk,e represents the travel time of conventional bus k through interval e.(16)Tcon,tl=∑k∈Kcon∑e∈Erk,e×tk,e.In summary, the expression of the total passenger timeTuser,0 in the conventional operation mode is as follows:(17)Tuser,0=Dconvpassenger+Ncon,wt×12hcon,0+∑k∈Kcon,0∑ei∈Erk,ei×tk,ei.(2) Total Passenger Time. Tuser,ain the Integrated Operation Mode. The total passenger time in the integrated operation mode includes the total passenger time on a fixed route Tfixed and the total passenger time on a customized route Tcus.Tcus is related to the type of customized bus in the integrated operation mode. The total passenger time for integrated operation mode II includes the passenger waiting time Tcus,wt and passenger travel time Tcus,tl. In addition, the passenger transfer time Tcus,tr must also be calculated for integrated operation mode I because passengers may need to wait for a conventional or customized bus at a transfer station at this time.Tcus,wt is the summation of the customized bus arrival times exceeding the origin time window, namely, Tcus,wt=∑j∈Qovertj. When a customized bus arrives at a passenger’s origin within a reasonable time window, the passenger does not generate a waiting time. Moreover, a passenger waiting time will only occur when a customized bus arrives at an origin outside its time window.(18)overtj=sj,min−tjk,iftjk<sj,min,0,ifsj,min<tjk<sj,max,tjk−sj,max,iftjk>sj,max,where overtj is the waiting time at point j, tjk is the arrival time of customized bus k, the reasonable time window is sj,min,sj,max, and the length of the time window is tw, namely, tw=sj,max−sj,min.Moreover,Tcus,tl is related to the routes of customized buses. Assuming that customized buses will depart Kcus times, Ni,jk represents the number of passengers on the bus when customized bus k passes through i,j, and Ti,jk represents the travel time of customized bus k passes through i,j. Then, the expression of Tcus,tl is as follows:(19)Tcus,tl=∑k∈Kcus∑i,j∈kNi,jk×Ti,jk.Tcus,wt is the sum of the waiting times when passengers wait for a customized/conventional bus at a transfer station. Therefore, it is composed of the total waiting time of passengers transferring from a customized bus to a conventional bus Tcus,wt,1 and the total waiting time of passengers transferring from a conventional bus to a customized bus Tcus,wt,2, which is denoted as Tcus,tw=Tcus,tw,1+Tcus,tw,2.As discussed previously, if the operation mode of integrated operation modea is ω and a customized bus replaces the service for m stations, the total passenger time Tuser,a of mode a in integrated operation mode ω can be expressed as follows:(20)Tuser,a=Tfixed+Tcus,Tfixed=Dconvpassenger+Ncon,wt×12hcon,ω+∑k∈Kcon,a∑e∈Erk,e×tk,e,Tcus=∑j∈Qovertj+∑k∈Kcus,a∑i,j∈kNi,jk×Ti,jk+z1×Tcus,wt.


## 3.1. Operation Mode

### 3.1.1. Conventional Operation Mode

The conventional operation mode (Figure1) is the most widely used operation mode of conventional bus systems, in which buses always run on urban roads according to fixed routes and fixed schedules. Passengers walk to the appropriate station from various origins, take the bus, get off at the appropriate station, and then walk to their destinations.Figure 1
The conventional operation mode.


### 3.1.2. Integrated Operation Mode

In the integrated operation mode, PT services are jointly completed by two types of buses, namely, conventional and customized buses. Specifically, in low-demand districts, customized bus services are used to replace portions of the original conventional bus services, and the remaining stations are still served by conventional buses. In this study, the integrated operation mode is divided into two types, namely, (1) integrated operation mode I: customized buses (feeder transit) and conventional buses and (2) integrated operation mode II: point-to-point (p2p) customized buses and conventional buses.(1) Integrated Operation Mode I: Customized Buses (Feeder Transit) and Conventional Buses. In this mode (Figure 2), a customized bus provides a connection service in its service area, picks up passengers according to the needs in the area, and transfers them to the transfer station; a conventional bus completes the rest of the transportation service. Specifically, a customized bus starts from the outer terminal station of the existing conventional bus route, picks up passengers according to their specific location within the m stations, and transfers them back to the transfer station to complete connection services. Then, the entire service is completed by the conventional bus covering the remaining n−m stations.Figure 2
Integrated operation mode I customized buses (feeder transit) and conventional buses.(2) Integrated Operation Mode II: Customized Buses (p2p) and Conventional Buses. In this mode (Figure 3), customized buses transfer some passengers from their origin to their destination, which is a type of “direct” public transport service. Passengers are no longer required to go to the departure bus station or get off at the arrival bus station. In other words, the p2p customized bus departs from the depot, picks up passengers according to their specific location within the m stations, transfers them to their destination, and finally returns to the depot. Other passengers continue to be served by a conventional bus.Figure 3
Integrated operation mode II: customized buses (p2p) and conventional buses.


## 3.1.1. Conventional Operation Mode

The conventional operation mode (Figure1) is the most widely used operation mode of conventional bus systems, in which buses always run on urban roads according to fixed routes and fixed schedules. Passengers walk to the appropriate station from various origins, take the bus, get off at the appropriate station, and then walk to their destinations.Figure 1
The conventional operation mode.


## 3.1.2. Integrated Operation Mode

In the integrated operation mode, PT services are jointly completed by two types of buses, namely, conventional and customized buses. Specifically, in low-demand districts, customized bus services are used to replace portions of the original conventional bus services, and the remaining stations are still served by conventional buses. In this study, the integrated operation mode is divided into two types, namely, (1) integrated operation mode I: customized buses (feeder transit) and conventional buses and (2) integrated operation mode II: point-to-point (p2p) customized buses and conventional buses.(1) Integrated Operation Mode I: Customized Buses (Feeder Transit) and Conventional Buses. In this mode (Figure 2), a customized bus provides a connection service in its service area, picks up passengers according to the needs in the area, and transfers them to the transfer station; a conventional bus completes the rest of the transportation service. Specifically, a customized bus starts from the outer terminal station of the existing conventional bus route, picks up passengers according to their specific location within the m stations, and transfers them back to the transfer station to complete connection services. Then, the entire service is completed by the conventional bus covering the remaining n−m stations.Figure 2
Integrated operation mode I customized buses (feeder transit) and conventional buses.(2) Integrated Operation Mode II: Customized Buses (p2p) and Conventional Buses. In this mode (Figure 3), customized buses transfer some passengers from their origin to their destination, which is a type of “direct” public transport service. Passengers are no longer required to go to the departure bus station or get off at the arrival bus station. In other words, the p2p customized bus departs from the depot, picks up passengers according to their specific location within the m stations, transfers them to their destination, and finally returns to the depot. Other passengers continue to be served by a conventional bus.Figure 3
Integrated operation mode II: customized buses (p2p) and conventional buses.


## 3.2. Cost Element

The total operating cost of a bus service includes the operating costs of both the conventional and customized buses, i.e.,(1)Ctotal=Ccon+Ccus,where Ccon is the conventional bus operating cost, Ccus is the customized bus operating cost, and Ctotal is the total operating cost. The conventional bus operating cost includes the fixed input cost Ccon,f and variable operating cost Ccon,m, and the same is true for the customized bus operating cost.


### 3.2.1. Fixed Input Cost

The fixed input cost is the sum of the fixed cost of each vehicle in the fleet. The fixed cost is generally related to the size of the fleet, the operating time, the drivers’ salaries, vehicle purchase fees, vehicle depreciation, insurance fees, etc. Formally, the fixed input costs of conventional busesCcon,f and customized buses Ccus,f can, respectively, be expressed as(2)Ccon⋅f=αcon⋅f×Gcon×T,Ccus⋅f=αcus⋅f×Gcus×T,where αcon,f is the fixed input cost per conventional bus hour, Gcon is the size of the conventional bus fleet, αcus,f is the fixed input cost per customized bus hour, Gcus is the size of the customized bus fleet, and T is the length of the bus service time. Gcon can be calculated based on the departure interval hcon by equation (4), where vbus is the average driving speed of the bus during the bus service time, and L is the total mileage of one cycle. Moreover, Gcus is related to the number of customized bus departures and each departure time.(3)Gcon=Lhcon×vbus.


### 3.2.2. Variable Operating Cost

The variable operating cost is generally related to the fuel cost, which is associated with the total mileage of the fleet. The variable operating costs of conventional busesCcon,m and customized buses Ccus,m can, respectively, be expressed as(4)Ccon,m=αcon,m×∑k∈KconLk,Ccus,m=αcus,m×∑k∈Kcus∑i,j∈UkIi,jk,where αcon⋅m is the fuel cost of conventional buses per mile, Lk is the total mileage of conventional bus k, Kcon is the set of conventional buses, αcus⋅m is the fuel cost of customized buses per mile, Kcus is the set of customized buses, Ii,jk is the mileage of customized bus k driving from point i to point j, and Uk is the set of route sections of customized bus k. Thus, ∑i,j∈UkIi,jk is the total mileage of customized bus k.


## 3.2.1. Fixed Input Cost

The fixed input cost is the sum of the fixed cost of each vehicle in the fleet. The fixed cost is generally related to the size of the fleet, the operating time, the drivers’ salaries, vehicle purchase fees, vehicle depreciation, insurance fees, etc. Formally, the fixed input costs of conventional busesCcon,f and customized buses Ccus,f can, respectively, be expressed as(2)Ccon⋅f=αcon⋅f×Gcon×T,Ccus⋅f=αcus⋅f×Gcus×T,where αcon,f is the fixed input cost per conventional bus hour, Gcon is the size of the conventional bus fleet, αcus,f is the fixed input cost per customized bus hour, Gcus is the size of the customized bus fleet, and T is the length of the bus service time. Gcon can be calculated based on the departure interval hcon by equation (4), where vbus is the average driving speed of the bus during the bus service time, and L is the total mileage of one cycle. Moreover, Gcus is related to the number of customized bus departures and each departure time.(3)Gcon=Lhcon×vbus.


## 3.2.2. Variable Operating Cost

The variable operating cost is generally related to the fuel cost, which is associated with the total mileage of the fleet. The variable operating costs of conventional busesCcon,m and customized buses Ccus,m can, respectively, be expressed as(4)Ccon,m=αcon,m×∑k∈KconLk,Ccus,m=αcus,m×∑k∈Kcus∑i,j∈UkIi,jk,where αcon⋅m is the fuel cost of conventional buses per mile, Lk is the total mileage of conventional bus k, Kcon is the set of conventional buses, αcus⋅m is the fuel cost of customized buses per mile, Kcus is the set of customized buses, Ii,jk is the mileage of customized bus k driving from point i to point j, and Uk is the set of route sections of customized bus k. Thus, ∑i,j∈UkIi,jk is the total mileage of customized bus k.


## 3.3. Assumption

As it would be difficult to determine the specific departure and destination locations of passengers during the operation of the conventional mode, certain assumptions are made regarding the characteristics of passenger travel demand based on actual conditions:(1)
The operation mode has no effect on bus travel demand, which means that conventional bus travel demand would be fully transferred by customized bus services in areas where the conventional operation mode has been changed.(2)
The location and travel time of passengers are randomly generated based on an assumed distribution.(3)
Customized bus provides pickup and drop-off service; therefore, the passenger’s walking time is ignored.


## 3.4. Objective Formulation

The research goal is to determine the optimal bus operation plan and the related operating parameters after clarifying the corresponding travel needs. Therefore, 0-1 integer programming was used to construct the objective function of the operation mode selection and optimization model. The objective function is divided into three parts, and its expression is(5)minC=∑ω=02Zω×Ctotal,ω,where ω represents the bus operation mode. When ω=0, Ctotal,0 is the total operating cost of the conventional operation mode; when ω=1, Ctotal,1 is the total operating cost of integrated operation mode I; and when ω=2, Ctotal,2 is the total operating cost of integrated operation mode II. Moreover, C is the operating cost of the optimal bus operation plan, and zω is a binary variable.As the conventional operation mode does not include customized buses, the operating cost of customized buses is considered as 0, and the expression ofCtotal,0 is(6)Ctotal,0=αcon,f×Gcon,0×T+αcon,m×∑k∈Kcon,0Lk,where Gcon,0 is the size of the conventional bus fleet in the conventional operation mode and Kcon,0 is the set of conventional buses in the conventional operation mode.The operating costs of the integrated operation modes, namely,Ctotal,1 and Ctotal,2, include two parts, namely, the operating cost of conventional bus operation and the operating cost of customized bus operation. The corresponding expressions are as follows:(7)Ctotal,1=αcon,f×Gcon,1×T+αcus,f×Gcus,1×T+αcon,m×∑k∈Kcon,1Lk∑k∈Kcus,1∑i,j∈UkIi,jk,Ctotal,2=αcon,f×Gcon,2×T+αcus,f×Gcus,2×T+αcon,m×∑k∈Kcon,2Lk∑k∈Kcus,2∑i,j∈UkIi,jk.


## 3.5. Constraints

Because the goal of the model is to reduce the operating costs while ensuring that the needs of passengers are met, the main constraints of this model are given as follows:(8)z0+z1+z2=1,(9)loadi,jk≤Capk∀i,j∈k,i≠j,∀k∈Kcus,(10)T¯user,a≤T¯user,0,(11)1≤m<n,m∈N,(12)∑k∈Kcus∑j∈Qxi,jk=1,∀i∈Q,i≠j,(13)∑k∈Kcus∑j∈Qxj,ik=1,∀i∈Q,i≠j,(14)∑j∈Qxi,jk−∑j∈Qxj,ik=0,∀k∈Kcus,i≠j.(8) is the mode selection constraint; there is only one variable, which corresponds to the operation mode of the optimal bus operation plan. (9) is the vehicle load constraint, which stipulates that the number of served passengers in any route section of customized bus k cannot exceed its capacity. (10) is the average passenger time constraint, which guarantees that the average passenger time T¯user,a of plan a in the integrated operation mode does not exceed the average passenger time T¯user,0 in the conventional operation mode. (11) restricts the number of replacement stations m to not exceed the total number of original conventional bus stations n. (12) and (13) indicate that for each passenger who uses a customized bus, there is only one customized bus to provide service, where Q is the set of origins of passengers using a customized bus. Moreover, xi,jk=1 represents that customized bus k travels from point i to point j; otherwise, xi,jk=0. Finally, (14) indicates that the customized bus that arrives at point i and departs from point i must be the same bus.


### 3.5.1. Total Passenger TimeTuser

(1) Total Passenger Time. Tuser,0in the Conventional Operation Mode. The total passenger time on a fixed route is Tfixed, which includes the passenger walking time Tcon,wk, the passenger waiting time Tcon,wt, and the passenger travel time Tcon,tl.(15)Tfixed=Tcon,wk+Tcon,wt+Tcon,tl,where Tcon,wk is the passenger walking time for the conventional bus mode, which is equal to the total passenger walking distance divided by the walking speed, i.e., Tcon,wk=Dcon/vpassenger, and Dcon is the summation of the walking distance of all passengers P from the origin to the departure station and from the arrival station to the destination, namely, Dcon=∑p∈PDu,p+Dd,p.Tcon,wt is the passenger waiting time for conventional buses. Assuming that the time for passengers to arrive at the departure stations is uniformly distributed, Tcon,wt can be calculated from the departure interval hcon, namely, Tcon,wt=Ncon,wt×1/2hcon, where Ncon,wt is the number of passengers waiting for conventional buses.Tcon,tl is the passenger travel time on conventional buses. The total travel time of passengers is the sum of the in-vehicle time of all passengers, which is equal to the multiplication of the number of passengers in each section by the section travel time. Hence, the calculation of Tcon,tl is converted into the calculation of the sum of the travel time of passengers in each interval between adjacent stations, where E represents the set of intervals e between adjacent stations on the fixed route, rk,e is the number of people within conventional bus k in interval e, and tk,e represents the travel time of conventional bus k through interval e.(16)Tcon,tl=∑k∈Kcon∑e∈Erk,e×tk,e.In summary, the expression of the total passenger timeTuser,0 in the conventional operation mode is as follows:(17)Tuser,0=Dconvpassenger+Ncon,wt×12hcon,0+∑k∈Kcon,0∑ei∈Erk,ei×tk,ei.(2) Total Passenger Time. Tuser,ain the Integrated Operation Mode. The total passenger time in the integrated operation mode includes the total passenger time on a fixed route Tfixed and the total passenger time on a customized route Tcus.Tcus is related to the type of customized bus in the integrated operation mode. The total passenger time for integrated operation mode II includes the passenger waiting time Tcus,wt and passenger travel time Tcus,tl. In addition, the passenger transfer time Tcus,tr must also be calculated for integrated operation mode I because passengers may need to wait for a conventional or customized bus at a transfer station at this time.Tcus,wt is the summation of the customized bus arrival times exceeding the origin time window, namely, Tcus,wt=∑j∈Qovertj. When a customized bus arrives at a passenger’s origin within a reasonable time window, the passenger does not generate a waiting time. Moreover, a passenger waiting time will only occur when a customized bus arrives at an origin outside its time window.(18)overtj=sj,min−tjk,iftjk<sj,min,0,ifsj,min<tjk<sj,max,tjk−sj,max,iftjk>sj,max,where overtj is the waiting time at point j, tjk is the arrival time of customized bus k, the reasonable time window is sj,min,sj,max, and the length of the time window is tw, namely, tw=sj,max−sj,min.Moreover,Tcus,tl is related to the routes of customized buses. Assuming that customized buses will depart Kcus times, Ni,jk represents the number of passengers on the bus when customized bus k passes through i,j, and Ti,jk represents the travel time of customized bus k passes through i,j. Then, the expression of Tcus,tl is as follows:(19)Tcus,tl=∑k∈Kcus∑i,j∈kNi,jk×Ti,jk.Tcus,wt is the sum of the waiting times when passengers wait for a customized/conventional bus at a transfer station. Therefore, it is composed of the total waiting time of passengers transferring from a customized bus to a conventional bus Tcus,wt,1 and the total waiting time of passengers transferring from a conventional bus to a customized bus Tcus,wt,2, which is denoted as Tcus,tw=Tcus,tw,1+Tcus,tw,2.As discussed previously, if the operation mode of integrated operation modea is ω and a customized bus replaces the service for m stations, the total passenger time Tuser,a of mode a in integrated operation mode ω can be expressed as follows:(20)Tuser,a=Tfixed+Tcus,Tfixed=Dconvpassenger+Ncon,wt×12hcon,ω+∑k∈Kcon,a∑e∈Erk,e×tk,e,Tcus=∑j∈Qovertj+∑k∈Kcus,a∑i,j∈kNi,jk×Ti,jk+z1×Tcus,wt.


## 3.5.1. Total Passenger TimeTuser

(1) Total Passenger Time. Tuser,0in the Conventional Operation Mode. The total passenger time on a fixed route is Tfixed, which includes the passenger walking time Tcon,wk, the passenger waiting time Tcon,wt, and the passenger travel time Tcon,tl.(15)Tfixed=Tcon,wk+Tcon,wt+Tcon,tl,where Tcon,wk is the passenger walking time for the conventional bus mode, which is equal to the total passenger walking distance divided by the walking speed, i.e., Tcon,wk=Dcon/vpassenger, and Dcon is the summation of the walking distance of all passengers P from the origin to the departure station and from the arrival station to the destination, namely, Dcon=∑p∈PDu,p+Dd,p.Tcon,wt is the passenger waiting time for conventional buses. Assuming that the time for passengers to arrive at the departure stations is uniformly distributed, Tcon,wt can be calculated from the departure interval hcon, namely, Tcon,wt=Ncon,wt×1/2hcon, where Ncon,wt is the number of passengers waiting for conventional buses.Tcon,tl is the passenger travel time on conventional buses. The total travel time of passengers is the sum of the in-vehicle time of all passengers, which is equal to the multiplication of the number of passengers in each section by the section travel time. Hence, the calculation of Tcon,tl is converted into the calculation of the sum of the travel time of passengers in each interval between adjacent stations, where E represents the set of intervals e between adjacent stations on the fixed route, rk,e is the number of people within conventional bus k in interval e, and tk,e represents the travel time of conventional bus k through interval e.(16)Tcon,tl=∑k∈Kcon∑e∈Erk,e×tk,e.In summary, the expression of the total passenger timeTuser,0 in the conventional operation mode is as follows:(17)Tuser,0=Dconvpassenger+Ncon,wt×12hcon,0+∑k∈Kcon,0∑ei∈Erk,ei×tk,ei.(2) Total Passenger Time. Tuser,ain the Integrated Operation Mode. The total passenger time in the integrated operation mode includes the total passenger time on a fixed route Tfixed and the total passenger time on a customized route Tcus.Tcus is related to the type of customized bus in the integrated operation mode. The total passenger time for integrated operation mode II includes the passenger waiting time Tcus,wt and passenger travel time Tcus,tl. In addition, the passenger transfer time Tcus,tr must also be calculated for integrated operation mode I because passengers may need to wait for a conventional or customized bus at a transfer station at this time.Tcus,wt is the summation of the customized bus arrival times exceeding the origin time window, namely, Tcus,wt=∑j∈Qovertj. When a customized bus arrives at a passenger’s origin within a reasonable time window, the passenger does not generate a waiting time. Moreover, a passenger waiting time will only occur when a customized bus arrives at an origin outside its time window.(18)overtj=sj,min−tjk,iftjk<sj,min,0,ifsj,min<tjk<sj,max,tjk−sj,max,iftjk>sj,max,where overtj is the waiting time at point j, tjk is the arrival time of customized bus k, the reasonable time window is sj,min,sj,max, and the length of the time window is tw, namely, tw=sj,max−sj,min.Moreover,Tcus,tl is related to the routes of customized buses. Assuming that customized buses will depart Kcus times, Ni,jk represents the number of passengers on the bus when customized bus k passes through i,j, and Ti,jk represents the travel time of customized bus k passes through i,j. Then, the expression of Tcus,tl is as follows:(19)Tcus,tl=∑k∈Kcus∑i,j∈kNi,jk×Ti,jk.Tcus,wt is the sum of the waiting times when passengers wait for a customized/conventional bus at a transfer station. Therefore, it is composed of the total waiting time of passengers transferring from a customized bus to a conventional bus Tcus,wt,1 and the total waiting time of passengers transferring from a conventional bus to a customized bus Tcus,wt,2, which is denoted as Tcus,tw=Tcus,tw,1+Tcus,tw,2.As discussed previously, if the operation mode of integrated operation modea is ω and a customized bus replaces the service for m stations, the total passenger time Tuser,a of mode a in integrated operation mode ω can be expressed as follows:(20)Tuser,a=Tfixed+Tcus,Tfixed=Dconvpassenger+Ncon,wt×12hcon,ω+∑k∈Kcon,a∑e∈Erk,e×tk,e,Tcus=∑j∈Qovertj+∑k∈Kcus,a∑i,j∈kNi,jk×Ti,jk+z1×Tcus,wt.


## 4. Solution Algorithm

The optimization of the proposed model is fundamentally an NP-hard problem, for which the use of intelligent algorithms is currently an efficient method.In this research, a hybrid genetic algorithm with a local search is used to solve the proposed model. In addition to the known constants, the fleet size and mileage of conventional and customized buses are variables that directly affect the results of the model and are the key variables to be optimized under the premise of satisfying the model constraints. The algorithm is designed as follows.Algorithm 1:Solution Algorithm.
(i)
Input: (1) travel demandsDEMAND, (2) the maximum number of iterations MAXGEN, (3) population size NIND, (4) number of original bus stations n(ii)
Output: (1) minimal operating costminC, (2) the optimal operating plan am,ω with minC(1)
Calculate the operating cost of the conventional operation modeCtotal,0(2)
Whilem<n(3)
Forω=1to 2(4)
Ifω=1(5)
ConstructQ according to operation principles of Integrated operation mode I(6)
Else(7)
ConstructQ according to operation principles of Integrated operation mode II(8)
End(9)
Fori=1toNIND(10)
EncodeQ into the genes of chromosome i(11)
i=i+1(12)
End(13)
Establish a primary populationChrom(14)
Whilegen<MAXGEN(15)
Determine the parent populationChrom(16)
Calculate the fitness valueobj of chromosome in Chrom(17)
Select some chromosomes through selection, crossover and mutation(18)
Optimize selected chromosomes by large-scale neighborhood search, form a newChrom(19)
IfNINDChrom<NIND(20)
Insert new chromosomes to addNINDChrom to NIND(21)
End(22)
gen=gen+1(23)
Find the chromosomeMaxCh with the highest fitness value in Chrom, then decode it to get optimal plan am,ω(24)
ω=ω+1(25)
End(26)
m=m+1(27)
Find the minimum operating costCtotal,1 among all am,ω whose ω=1(28)
Find the minimum operating costCtotal,2 among all am,ω whose ω=2(29)
SubstituteCtotal,0, Ctotal,1, Ctotal,2 into objective formulation and solve it(30)
Return minimal operating cost minC, and the optimal operating plan am,ω


## 5. Case Study

### 5.1. Numerical Information

As shown in Figure4, in this experiment, the total length of the original bus line was 8.4 km, the line passed 15 stations, the distance between adjacent stations was 600 m, and the departure interval of conventional buses was 15 min. Furthermore, to be consistent with actual circumstances, the passengers’ origin and destination, departure and arrival stations, and ideal travel time met the following experimental assumptions:(1)
The number of boarding passengers decreases linearly from the urban district to the low-demand district, and the demand of the farthest station is 20% of the first station in the urban district(2)
The average travel distance of passengers is 3 km(3)
The passengers’ origin and destination locations are randomly located within 500 m of their departure and arrival stations(4)
The ideal travel time of a passenger is randomly generated within the operating service timeFigure 4
The bus route and passenger distribution area.


### 5.2. Results Analysis

The experiment was carried out in the MATLAB 2018b environment. The total travel demands were set to 200 passengers per hour, and the values of the known constants are reported in Table2. An experiment dealt with 200 demands, as the service time T was 1 hour.Table 2
The values of known constants.
Constantαcon,fαcus,fTαcon,mαcus,mCapkvbusvpassengertwValue30 yuan/h25 yuan/h1 h15 yuan/km8 yuan/km1520 km/h10 km/h40 minPassengers/vehicleBecause the passenger demand was randomly generated, the experiment was repeated five times, and the average value was taken as the final result. The operation plans under each operation mode are exhibited in Table3. The optimal operation plan was integrated operation mode I (feeder transit and conventional buses) that 4 of 15 stations were replaced and the passengers were transfered with 4 customized buses. Compared with the conventional operation mode, the operating cost of this plan was found to be reduced by 124.08 yuan per hour.Table 3
The optimal bus operation plans in different modes (travel demand = 200 passengers per hour).
ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Headway/minFleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)Conventional operation modeω=0—15.00667.20———1188.00Combined operation mode Iω=1415.00448.0072415.491063.92Combined operation mode IIω=2322.50435.6454442.051090.97Figure5 exhibits the change in the operating cost per capita with the increase in the number of replaced stations. The following can be observed: (1) the operating costs per capita of the two integrated operation modes were both found to decrease and then increase with the increase of the number of replaced stations and (2) the numerical range of the reduction of the operating cost per capita of integrated operation mode I was found to be more affected than that of integrated operation mode II.Figure 5
The reductions in the operating cost per capita in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)Figure6 presents the comparison of the operating costs between the conventional and customized buses in the two integrated operation modes. It was found that the steady increase in operating costs with the increase in the number of replaced stations influenced by the cost of operating customized buses. The cost proportion of conventional bus operating costs in integrated operation mode I was found to increase at a slower rate than that in integrated operation mode II, i.e., at the same number of replaced stations, the customized bus operating cost in integrated operation mode I was lower than that in integrated operation mode II.Figure 6
The compositions of operating costs in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)


### 5.3. Sensitivity Analysis

Governments and bus companies are facing diverse operating environments. Therefore, it is necessary to explore the impact of changes in some key parameters on the results.


#### 5.3.1. Travel Demand

Travel demand is a key parameter that affects the change of the optimal bus operation plan. To study the impact of travel demand, the optimal solutions for different values of passenger demand were compared.Figure7 exhibits the changes in the operating cost per capita of the two integrated operation modes under changes in both the travel demand and the number of replaced stations. It should be noted that although the degrees of change in the two integrated operation modes were different, the trends were the same. Regarding the travel demand, under low-demand conditions (less than 150 passengers per hour), the reduction in the operating cost was found to increase with the increase of the number of replaced stations, especially for integrated operation mode I; this indicates that integrated operation mode I exhibited an obvious advantage in reducing operating costs in low-demand districts. Regarding the medium-demand conditions (greater than or equal to 150 and less than 300 passengers per hour), integrated operation mode I also exhibited stable economic benefits when few stations were replaced.Figure 7
The reductions in the operating cost per capita in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)Table4 summarizes the optimal bus operation plans for different travel demands. With the increase in the number of passengers, the change of the optimal bus operation plan exhibited three obvious characteristics: (1) the mode changed from integrated operation mode I to integrated operation mode II and then to the conventional operation mode, (2) the number of replaced stations was gradually reduced, and (3) the reduction in the cost per capita was gradually reduced. The analysis results support the use of the integrated operation mode to replace the conventional operation mode in some cases and reveal the economic benefits of both integrated operation modes.Table 4
The optimal bus operation plans with different travel demands.
Travel demands (passengers/per hour)ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)5011127.5018.4948328.73462.1114.5210011020.00214.4093328.61579.886.081301618.13432.8877525.75944.221.881501617.50433.6086631.971029.741.061801617.50433.60101631.391025.100.912001415.00448.0072415.491063.920.622302324.00433.2657442.191056.470.572502323.00434.8559342.981061.570.512802324.00433.2666449.031111.210.273001315.00552.8078412.821144.550.143501215.00557.605136.501140.970.134001215.00557.606247.901177.180.034501215.00557.606548.571182.530.015000—15.00667.20———1188.000.005500—15.00667.20———1188.000.006000—15.00667.20———1188.000.006500—15.00667.20———1188.000.007000—15.00667.20———1188.000.008000—15.00667.20———1188.000.00Furthermore, Tables5 and 6 report the optimal plan indicators of the two integrated operation modes under different travel demands. It is evident that the two integrated operation modes exhibited similar trends in terms of the reduction in the cost per capita, but the operating cost of integrated operation mode I tended to remain unchanged after 500 passengers per hour whereas that of integrated operation mode II continued to increase.Table 5
The optimal bus operation plans with different travel demands in integrated operation mode I.
Travel demands (passengers/per hour)Total replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)501127.5018.4948328.73462.1114.521001020.00214.4093328.61579.886.08130618.13432.8877525.75944.221.88150617.50433.6086631.971029.741.06180617.50433.60101631.391025.100.91200415.00448.0072415.491063.920.62230315.00552.805037.541077.310.48250315.00552.806049.011114.050.30280315.00552.806549.601118.770.25300315.00552.8078412.821144.550.14350215.00557.605136.501140.970.13400215.00557.606247.901177.180.03450215.00557.606546.571166.530.05500115.00662.403532.681212.45−0.05550115.00662.403732.801213.42−0.05600115.00662.404032.991214.95−0.04650115.00662.404543.861246.85−0.09700115.00662.404944.031248.27−0.09800115.00662.405344.221249.77−0.08Table 6
The optimal bus operation plans with different travel demands in integrated operation mode II.
Travel demands (passengers/per hour)Total replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)501330.0012.4050449.90565.2312.461001035.0018.2393463.97765.214.23130535.00218.8761659.26967.061.70150529.38322.5669658.531046.740.94180321.67436.9642436.561066.900.67200322.50435.6454442.051090.970.49230324.00433.2657442.191056.470.57250323.00434.8559342.981061.570.51280324.00433.2666449.031111.210.27300220.00443.2048436.751162.030.09350219.00546.0857437.011237.25−0.14400218.75546.8056439.331266.67−0.20450115.00662.4031324.361385.87−0.44500115.00662.4035324.461386.71−0.40550115.00662.4037327.631412.06−0.41600115.00662.4040333.361457.86−0.45650115.00662.4045433.521484.19−0.46700115.00662.4049435.191497.55−0.44800115.00662.4053643.121610.98−0.53


#### 5.3.2. Service Coverage Radius of a Station

The maximum walking distance of passengers for a new conventional bus line is limited in most large cities with developed PT [24]. However, this criterion is clearly excessively stringent and does not reflect reality in low-demand areas. Therefore, the influence of the service coverage radius is discussed in this section. The service range radius r was respectively set to 250 m, 500 m, 750 m, and 1 km. Each experiment was carried out three times, and the average value was used as the outcome.Table7 presents the optimal operation plans for different radii. It was found that more stations can be replaced as the service coverage radius limit of a station increased. This means that the integrated operation mode has a distinct advantage in areas with a limited number of bus stops. Furthermore, in all scenarios, integrated operation mode I was found to be the best option, thereby demonstrating its versatility in a variety of situations.Table 7
The optimal bus operation plans with different service range radii (travel demand = 200 passengers per hour).
Service coverage radius (m)ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)2501415.00448.0068411.281030.280.795001415.00448.0072415.491063.920.627501517.50437.8085524.521008.160.9010001517.50437.8096728.151087.200.50


#### 5.3.3. Departure Interval of Conventional Bus

The departure interval of a conventional bus is the key factor affecting the operating cost. Figure8 reveals that the operating cost of the conventional operation mode was found to gradually decrease with the increase of the departure interval, whereas the operating costs of the other two integrated operation modes began to increase after failing to a certain degree due to the overload demand of customized buses.Figure 8
The total operating costs of different modes (the number of replaced stations is 5).


## 5.1. Numerical Information

As shown in Figure4, in this experiment, the total length of the original bus line was 8.4 km, the line passed 15 stations, the distance between adjacent stations was 600 m, and the departure interval of conventional buses was 15 min. Furthermore, to be consistent with actual circumstances, the passengers’ origin and destination, departure and arrival stations, and ideal travel time met the following experimental assumptions:(1)
The number of boarding passengers decreases linearly from the urban district to the low-demand district, and the demand of the farthest station is 20% of the first station in the urban district(2)
The average travel distance of passengers is 3 km(3)
The passengers’ origin and destination locations are randomly located within 500 m of their departure and arrival stations(4)
The ideal travel time of a passenger is randomly generated within the operating service timeFigure 4
The bus route and passenger distribution area.


## 5.2. Results Analysis

The experiment was carried out in the MATLAB 2018b environment. The total travel demands were set to 200 passengers per hour, and the values of the known constants are reported in Table2. An experiment dealt with 200 demands, as the service time T was 1 hour.Table 2
The values of known constants.
Constantαcon,fαcus,fTαcon,mαcus,mCapkvbusvpassengertwValue30 yuan/h25 yuan/h1 h15 yuan/km8 yuan/km1520 km/h10 km/h40 minPassengers/vehicleBecause the passenger demand was randomly generated, the experiment was repeated five times, and the average value was taken as the final result. The operation plans under each operation mode are exhibited in Table3. The optimal operation plan was integrated operation mode I (feeder transit and conventional buses) that 4 of 15 stations were replaced and the passengers were transfered with 4 customized buses. Compared with the conventional operation mode, the operating cost of this plan was found to be reduced by 124.08 yuan per hour.Table 3
The optimal bus operation plans in different modes (travel demand = 200 passengers per hour).
ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Headway/minFleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)Conventional operation modeω=0—15.00667.20———1188.00Combined operation mode Iω=1415.00448.0072415.491063.92Combined operation mode IIω=2322.50435.6454442.051090.97Figure5 exhibits the change in the operating cost per capita with the increase in the number of replaced stations. The following can be observed: (1) the operating costs per capita of the two integrated operation modes were both found to decrease and then increase with the increase of the number of replaced stations and (2) the numerical range of the reduction of the operating cost per capita of integrated operation mode I was found to be more affected than that of integrated operation mode II.Figure 5
The reductions in the operating cost per capita in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)Figure6 presents the comparison of the operating costs between the conventional and customized buses in the two integrated operation modes. It was found that the steady increase in operating costs with the increase in the number of replaced stations influenced by the cost of operating customized buses. The cost proportion of conventional bus operating costs in integrated operation mode I was found to increase at a slower rate than that in integrated operation mode II, i.e., at the same number of replaced stations, the customized bus operating cost in integrated operation mode I was lower than that in integrated operation mode II.Figure 6
The compositions of operating costs in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)


## 5.3. Sensitivity Analysis

Governments and bus companies are facing diverse operating environments. Therefore, it is necessary to explore the impact of changes in some key parameters on the results.


### 5.3.1. Travel Demand

Travel demand is a key parameter that affects the change of the optimal bus operation plan. To study the impact of travel demand, the optimal solutions for different values of passenger demand were compared.Figure7 exhibits the changes in the operating cost per capita of the two integrated operation modes under changes in both the travel demand and the number of replaced stations. It should be noted that although the degrees of change in the two integrated operation modes were different, the trends were the same. Regarding the travel demand, under low-demand conditions (less than 150 passengers per hour), the reduction in the operating cost was found to increase with the increase of the number of replaced stations, especially for integrated operation mode I; this indicates that integrated operation mode I exhibited an obvious advantage in reducing operating costs in low-demand districts. Regarding the medium-demand conditions (greater than or equal to 150 and less than 300 passengers per hour), integrated operation mode I also exhibited stable economic benefits when few stations were replaced.Figure 7
The reductions in the operating cost per capita in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)Table4 summarizes the optimal bus operation plans for different travel demands. With the increase in the number of passengers, the change of the optimal bus operation plan exhibited three obvious characteristics: (1) the mode changed from integrated operation mode I to integrated operation mode II and then to the conventional operation mode, (2) the number of replaced stations was gradually reduced, and (3) the reduction in the cost per capita was gradually reduced. The analysis results support the use of the integrated operation mode to replace the conventional operation mode in some cases and reveal the economic benefits of both integrated operation modes.Table 4
The optimal bus operation plans with different travel demands.
Travel demands (passengers/per hour)ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)5011127.5018.4948328.73462.1114.5210011020.00214.4093328.61579.886.081301618.13432.8877525.75944.221.881501617.50433.6086631.971029.741.061801617.50433.60101631.391025.100.912001415.00448.0072415.491063.920.622302324.00433.2657442.191056.470.572502323.00434.8559342.981061.570.512802324.00433.2666449.031111.210.273001315.00552.8078412.821144.550.143501215.00557.605136.501140.970.134001215.00557.606247.901177.180.034501215.00557.606548.571182.530.015000—15.00667.20———1188.000.005500—15.00667.20———1188.000.006000—15.00667.20———1188.000.006500—15.00667.20———1188.000.007000—15.00667.20———1188.000.008000—15.00667.20———1188.000.00Furthermore, Tables5 and 6 report the optimal plan indicators of the two integrated operation modes under different travel demands. It is evident that the two integrated operation modes exhibited similar trends in terms of the reduction in the cost per capita, but the operating cost of integrated operation mode I tended to remain unchanged after 500 passengers per hour whereas that of integrated operation mode II continued to increase.Table 5
The optimal bus operation plans with different travel demands in integrated operation mode I.
Travel demands (passengers/per hour)Total replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)501127.5018.4948328.73462.1114.521001020.00214.4093328.61579.886.08130618.13432.8877525.75944.221.88150617.50433.6086631.971029.741.06180617.50433.60101631.391025.100.91200415.00448.0072415.491063.920.62230315.00552.805037.541077.310.48250315.00552.806049.011114.050.30280315.00552.806549.601118.770.25300315.00552.8078412.821144.550.14350215.00557.605136.501140.970.13400215.00557.606247.901177.180.03450215.00557.606546.571166.530.05500115.00662.403532.681212.45−0.05550115.00662.403732.801213.42−0.05600115.00662.404032.991214.95−0.04650115.00662.404543.861246.85−0.09700115.00662.404944.031248.27−0.09800115.00662.405344.221249.77−0.08Table 6
The optimal bus operation plans with different travel demands in integrated operation mode II.
Travel demands (passengers/per hour)Total replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)501330.0012.4050449.90565.2312.461001035.0018.2393463.97765.214.23130535.00218.8761659.26967.061.70150529.38322.5669658.531046.740.94180321.67436.9642436.561066.900.67200322.50435.6454442.051090.970.49230324.00433.2657442.191056.470.57250323.00434.8559342.981061.570.51280324.00433.2666449.031111.210.27300220.00443.2048436.751162.030.09350219.00546.0857437.011237.25−0.14400218.75546.8056439.331266.67−0.20450115.00662.4031324.361385.87−0.44500115.00662.4035324.461386.71−0.40550115.00662.4037327.631412.06−0.41600115.00662.4040333.361457.86−0.45650115.00662.4045433.521484.19−0.46700115.00662.4049435.191497.55−0.44800115.00662.4053643.121610.98−0.53


### 5.3.2. Service Coverage Radius of a Station

The maximum walking distance of passengers for a new conventional bus line is limited in most large cities with developed PT [24]. However, this criterion is clearly excessively stringent and does not reflect reality in low-demand areas. Therefore, the influence of the service coverage radius is discussed in this section. The service range radius r was respectively set to 250 m, 500 m, 750 m, and 1 km. Each experiment was carried out three times, and the average value was used as the outcome.Table7 presents the optimal operation plans for different radii. It was found that more stations can be replaced as the service coverage radius limit of a station increased. This means that the integrated operation mode has a distinct advantage in areas with a limited number of bus stops. Furthermore, in all scenarios, integrated operation mode I was found to be the best option, thereby demonstrating its versatility in a variety of situations.Table 7
The optimal bus operation plans with different service range radii (travel demand = 200 passengers per hour).
Service coverage radius (m)ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)2501415.00448.0068411.281030.280.795001415.00448.0072415.491063.920.627501517.50437.8085524.521008.160.9010001517.50437.8096728.151087.200.50


### 5.3.3. Departure Interval of Conventional Bus

The departure interval of a conventional bus is the key factor affecting the operating cost. Figure8 reveals that the operating cost of the conventional operation mode was found to gradually decrease with the increase of the departure interval, whereas the operating costs of the other two integrated operation modes began to increase after failing to a certain degree due to the overload demand of customized buses.Figure 8
The total operating costs of different modes (the number of replaced stations is 5).


## 5.3.1. Travel Demand

Travel demand is a key parameter that affects the change of the optimal bus operation plan. To study the impact of travel demand, the optimal solutions for different values of passenger demand were compared.Figure7 exhibits the changes in the operating cost per capita of the two integrated operation modes under changes in both the travel demand and the number of replaced stations. It should be noted that although the degrees of change in the two integrated operation modes were different, the trends were the same. Regarding the travel demand, under low-demand conditions (less than 150 passengers per hour), the reduction in the operating cost was found to increase with the increase of the number of replaced stations, especially for integrated operation mode I; this indicates that integrated operation mode I exhibited an obvious advantage in reducing operating costs in low-demand districts. Regarding the medium-demand conditions (greater than or equal to 150 and less than 300 passengers per hour), integrated operation mode I also exhibited stable economic benefits when few stations were replaced.Figure 7
The reductions in the operating cost per capita in (a) integrated operation mode I and (b) integrated operation mode II.
(a)(b)Table4 summarizes the optimal bus operation plans for different travel demands. With the increase in the number of passengers, the change of the optimal bus operation plan exhibited three obvious characteristics: (1) the mode changed from integrated operation mode I to integrated operation mode II and then to the conventional operation mode, (2) the number of replaced stations was gradually reduced, and (3) the reduction in the cost per capita was gradually reduced. The analysis results support the use of the integrated operation mode to replace the conventional operation mode in some cases and reveal the economic benefits of both integrated operation modes.Table 4
The optimal bus operation plans with different travel demands.
Travel demands (passengers/per hour)ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)5011127.5018.4948328.73462.1114.5210011020.00214.4093328.61579.886.081301618.13432.8877525.75944.221.881501617.50433.6086631.971029.741.061801617.50433.60101631.391025.100.912001415.00448.0072415.491063.920.622302324.00433.2657442.191056.470.572502323.00434.8559342.981061.570.512802324.00433.2666449.031111.210.273001315.00552.8078412.821144.550.143501215.00557.605136.501140.970.134001215.00557.606247.901177.180.034501215.00557.606548.571182.530.015000—15.00667.20———1188.000.005500—15.00667.20———1188.000.006000—15.00667.20———1188.000.006500—15.00667.20———1188.000.007000—15.00667.20———1188.000.008000—15.00667.20———1188.000.00Furthermore, Tables5 and 6 report the optimal plan indicators of the two integrated operation modes under different travel demands. It is evident that the two integrated operation modes exhibited similar trends in terms of the reduction in the cost per capita, but the operating cost of integrated operation mode I tended to remain unchanged after 500 passengers per hour whereas that of integrated operation mode II continued to increase.Table 5
The optimal bus operation plans with different travel demands in integrated operation mode I.
Travel demands (passengers/per hour)Total replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)501127.5018.4948328.73462.1114.521001020.00214.4093328.61579.886.08130618.13432.8877525.75944.221.88150617.50433.6086631.971029.741.06180617.50433.60101631.391025.100.91200415.00448.0072415.491063.920.62230315.00552.805037.541077.310.48250315.00552.806049.011114.050.30280315.00552.806549.601118.770.25300315.00552.8078412.821144.550.14350215.00557.605136.501140.970.13400215.00557.606247.901177.180.03450215.00557.606546.571166.530.05500115.00662.403532.681212.45−0.05550115.00662.403732.801213.42−0.05600115.00662.404032.991214.95−0.04650115.00662.404543.861246.85−0.09700115.00662.404944.031248.27−0.09800115.00662.405344.221249.77−0.08Table 6
The optimal bus operation plans with different travel demands in integrated operation mode II.
Travel demands (passengers/per hour)Total replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)501330.0012.4050449.90565.2312.461001035.0018.2393463.97765.214.23130535.00218.8761659.26967.061.70150529.38322.5669658.531046.740.94180321.67436.9642436.561066.900.67200322.50435.6454442.051090.970.49230324.00433.2657442.191056.470.57250323.00434.8559342.981061.570.51280324.00433.2666449.031111.210.27300220.00443.2048436.751162.030.09350219.00546.0857437.011237.25−0.14400218.75546.8056439.331266.67−0.20450115.00662.4031324.361385.87−0.44500115.00662.4035324.461386.71−0.40550115.00662.4037327.631412.06−0.41600115.00662.4040333.361457.86−0.45650115.00662.4045433.521484.19−0.46700115.00662.4049435.191497.55−0.44800115.00662.4053643.121610.98−0.53


## 5.3.2. Service Coverage Radius of a Station

The maximum walking distance of passengers for a new conventional bus line is limited in most large cities with developed PT [24]. However, this criterion is clearly excessively stringent and does not reflect reality in low-demand areas. Therefore, the influence of the service coverage radius is discussed in this section. The service range radius r was respectively set to 250 m, 500 m, 750 m, and 1 km. Each experiment was carried out three times, and the average value was used as the outcome.Table7 presents the optimal operation plans for different radii. It was found that more stations can be replaced as the service coverage radius limit of a station increased. This means that the integrated operation mode has a distinct advantage in areas with a limited number of bus stops. Furthermore, in all scenarios, integrated operation mode I was found to be the best option, thereby demonstrating its versatility in a variety of situations.Table 7
The optimal bus operation plans with different service range radii (travel demand = 200 passengers per hour).
Service coverage radius (m)ModeωTotal replaced stationsConventional busCustomized busTotal operating cost (yuan)Saving operating cost per capita (yuan)Headway (min)Fleet size (vehicles)Mileage (km)The number of passengers servicedFleet size (vehicles)Mileage (km)2501415.00448.0068411.281030.280.795001415.00448.0072415.491063.920.627501517.50437.8085524.521008.160.9010001517.50437.8096728.151087.200.50


## 5.3.3. Departure Interval of Conventional Bus

The departure interval of a conventional bus is the key factor affecting the operating cost. Figure8 reveals that the operating cost of the conventional operation mode was found to gradually decrease with the increase of the departure interval, whereas the operating costs of the other two integrated operation modes began to increase after failing to a certain degree due to the overload demand of customized buses.Figure 8
The total operating costs of different modes (the number of replaced stations is 5).


## 6. Conclusion

This paper proposed an optimization model for bus operation planning based on cost analysis for low-demand districts. In addition to the conventional operation mode, two integrated operation modes were proposed to simultaneously optimize key operation variables such as the headway, fleet size, and service area under certain demand. The results of experiments revealed that the proposed model and solution algorithm can provide reliable results under different circumstances. Nevertheless, this study faced some limitations as follows: (1) the impacts of passenger loyalty and continuous usage intention were not considered, which could have resulted in differences in the experimental environment as a result of the reduced customized bus demand, and (2) in this study, the Euclidean distance between passengers and stations was used as the passenger walking distance to simplify the calculation. More accurate computation results could be obtained by applying the real walking distance of passengers according to the road network geography. Future studies should also consider the use of real data to validate the proposed model.



---

*Source: 4142517-2022-06-03.xml*