An Optimization Strategy for Truncating Ultralong Bus Lines Integrated with Metro Networks

,


Introduction
Te current public transportation system comprises bus lines of varying lengths, designed to provide travel service for passengers in residential areas or the central business district (CBD).Bus lines that surpass a specifc length threshold are referred to as ultralong bus lines, with the threshold determined based on the city's development scale and the average speed of buses [1].Trough an analysis of route and passenger fow characteristics on ultralong bus lines, it was found that although these lines enhance the trafc radiation of the metropolis to the surrounding towns and reduce transfers, long-distance travel often leads to a low punctuality rate, increasing the probability of unknown security risks.In addition, it will also lead to an uneven distribution of passenger fow both in time and space, resulting in the frequent occurrence of bus bunching [2,3].Terefore, in order to avoid these negative efects, such ultralong bus lines in cities need to be optimized or adjusted.
Currently, the research on the operational strategy of ultralong bus lines mainly focuses on three aspects: line truncation, line shortening, and skip-stop service.For the frst strategy, a bus line is split into two or more lines with unchanged stops [4,5], while the second strategy involves creating a shorter line by selecting stops from a specifc section of the original route; the last strategy is commonly employed in a public transit system that a vehicle passes through some stops that may not be continuous in a line [6,7].Early research on the abovementioned operation strategies mainly focuses on passenger fow distribution (PFD).Zhang et al. [8] proposed splitting strategies for long bus lines based on PDF and the average distance of the bus line by introducing the alternating coefcient.Te skip-stop operation strategy is adopted when PFD is uniform, whereas for lines with nonuniform characteristics such as unimodal, symmetrical bimodal, and asymmetric bimodal PFD, the line-truncated operation strategy is commonly applied, and the truncated points are typically located near peak or valley points [5].In addition to relying primarily on PFD, researchers gradually focus on optimization models aimed at saving time or cost [9,10].Tang et al. [11] developed a bus skip-stop operation optimization model with the aim of minimizing system costs for passengers and operators while simultaneously enhancing the energy utilization rate of electric buses.
In terms of adjustment or optimization of ultralong bus lines, most scholars mainly focus on factors such as PFD and unbalance coefcient, as well as passenger or operator cost.However, the combination with the metro network is often overlooked in operating strategies, aimed at enhancing these bus lines just like the rail transit system in many megacities that have been well developed.Furthermore, since metro and bus are two integral components of urban public transport systems, it is necessary to strengthen the relationship between bus lines and metro networks to better solve public transit operation problems [12,13].Promoting the development and utilization of urban rail transit can also contribute in reducing greenhouse gas emissions and alleviating road congestion [14].However, to the best of our knowledge, few studies pay special attention to ultralong bus lines with extremely long travel time and low occupancy rate, further reducing passenger satisfaction and operating benefts.
With the aim of fully improving the efciency of urban public transport system, this paper provides a novel scheduling methodology for optimizing ultralong bus lines.We made the following threefold contributions: frst, we defne the ultralong line evaluation index to flter out target lines that need to be optimized.Second, we introduce a truncation optimization procedure to improve ultralong bus lines.A multiobjective optimization model is developed with the consideration of mutivehicle type and cost of both passenger and operator, combining with the existing coline metro networks, and NSGA-II is employed to solve the model in order to generate the optimized truncated bus stop.Tird, we perform extensive numerical experiments on a real-world instance based on IC card data, which yield valuable insights into the capability of the proposed optimization strategy.
Te rest of the paper is listed as follows.Section 2 presents a review of the literature.Section 3 describes the analysis of ultralong bus lines.Section 4 introduces the multiobjective optimization model, and Section 5 represents the specifc algorithm according to the model.Section 6 is case analysis.Section 7 is the conclusion, which summarizes the results of case analysis.

Literature Review
2.1.Bus Line and Vehicle Scheduling Optimization.Te bus service planning process roughly includes the following fve stages: route planning, frequency setting, schedule making, vehicle scheduling, and personnel arrangement.A highquality bus service requires all stages to be optimized.Among them, the optimization of route planning is the basis of other stages, and it has been extensively studied by many scholars.Lin et al. [15] considered improving bus lines to enhance the quality of passenger travel services by shortening the total length of bus routes under the premise of fully meeting the travel needs of passengers.Since an ultralong bus line may afect the punctuality rate of vehicles and cause vehicle bunching, a route optimization method referring to long-distance commuter buses is proposed to determine bus route length and each stop location and to minimize the total operating distance of the bus line under the limitation of the maximum capacity [16].Zhang et al. [17] proposed dividing the long through-type bus lines into several independent bus lines by adjusting the stop setting and headway of bus lines, which not only efectively reduced the operational costs but also improved the unbalanced distribution of passenger fow on bus lines.Te abovementioned studies mainly consider the benefts or costs of passengers and operators separately.Wu et al. [18] constructed a model with the objective of minimizing the costs of both passengers and operators using continuous functions, which saved the operating costs of enterprises and improved the travel satisfaction of passengers concurrently.Furthermore, the vehicle type is also considered in the route optimization problem [19].For regular bus lines, previous studies mainly considered the optimization of stop location and route length, combined with the distribution of passenger fow, with the aim of reducing route length and minimizing the cost to the operators and the average travel time of passengers.
Regarding the scheduling optimization of public transport vehicles, the existing studies take the uncertainty of vehicle travel time and passenger arrival into account [20][21][22][23].Zhang et al. [20] and Li et al. [21] considered the uncertainty of bus running time; the former optimized the existing bus schedule by fne-tuning vehicle departure time and parking vehicle at other stops, and the latter proposed a schedule optimization model containing fuzzy travel time by considering the bus running time as a fuzzy variable subject to various disturbances such as weather and trafc conditions.On the other hand, Liu et al. [22] considered the uncertainty of passenger arrival rate and built a robust optimization model of dynamic bus departure considering diferent passenger arrival rates under diferent scenarios by setting three scenarios of baseline passenger fow, high passenger fow, and low passenger fow so as to achieve higher operational efciency.Xue et al. [23] simultaneously considered the uncertainty of bus running time and passenger arrival rate and established a bilevel programming model in order to tackle the bus scheduling problem.Te abovementioned studies consider the efect of uncertain factors individually or simultaneously in determining an optimal vehicle schedule without combining associated data.Ma et al. [24] developed a datadriven bus schedule utilizing bus GPS data and IC card data to replace the current schedules based on human experience.In the abovementioned research, when combining uncertain factors or historical data to optimize vehicle scheduling, the vehicle capacity is set as fxed, while the vehicle scheduling problem referring to multiple vehicle types is not considered.Sun et al. [25] proposed a fexible scheduling optimization method considering 2 Journal of Advanced Transportation mixed vehicle types.Te abovementioned studies mostly solve the bus route design or vehicle scheduling problem separately, but less consider jointly optimizing the above two problems.
In the existing literature, the joint optimization for both bus routes and vehicle scheduling is mainly applied in various scenarios and also considers relative factors.In terms of the objective function, Yu et al. [26] considered the passenger cost, and with the objective of minimizing the total travel time of passengers, a vehicle scheduling model was proposed to determine vehicle frequency and bus travel route.Besides, Kenyon et al. [27] constructed an optimization model with the objective of minimizing the total vehicle travel time in view of path issues of vehicles.Te abovementioned studies establish an optimization model based on a single-objective function such as passenger benefts or operator cost.However, Chen et al. [28] and Israeli et al. [29] considered both operator and passenger costs to build a model and, respectively, solved the model by using a forbidden search heuristic algorithm and an exhaustive search algorithm so as to improve the efciency of the bus service while reducing the travel time of passengers.Li et al. [30] and Ahern et al. [31] took both operator and passenger benefts into account and established a biobjective optimization model.Te former established a stochastic expected value model for bus routes in response to the random behaviors of passengers and bus vehicles and used a hybrid intelligent algorithm combining stochastic simulation and genetic algorithms to solve it.Te latter proposed a model combining vehicle scheduling and the assignment of passengers to reduce the operating costs while ensuring that passengers reached their destinations without afecting vehicle travel time.Yu et al. [32] introduced headway in modeling while considering operator and passenger costs and proposed a headway optimization model for bus routes, which was solved by a parallel genetic algorithm.In the abovementioned studies, the objectives of joint optimization are divided into single objective and multiobjectives.Te establishment of a multiobjective joint optimization model can take into account the benefts of both passengers and operators and the compromised proposal obtained through the solution can improve the satisfaction of both parties, which is more applicable than the single-objective optimization model.Some studies take the headway into account, and the established model referring to diferent types of vehicle capacity is less analyzed.In terms of model solving, it mainly consists of designing multiobjective algorithms [28][29][30], which is generally done by converting a multiobjective problem into a single-objective problem to be solved by setting weighting values for subobjectives [31,32].However, the weights of diferent objective values are often difcult to be determined, and if the weights are not set properly, it will have a greater impact on the optimization results.In contrast, the multiobjective optimization algorithm can efectively avoid the problem of improperly setting the weights of subobjectives; thus, it is easier to obtain a better solution.In this paper, with the objective of minimizing the cost for both operators and passengers, a model with consideration of diferent vehicle types is developed to simultaneously tackle bus routes designing and vehicle scheduling problem, and it is solved by using the fast nondominant sorting genetic algorithm II (NSGA-II).

Bus Line Optimization Combined with Metro Networks.
Buses and metros are important components of the urban public transportation system, so it is of great signifcance to optimize bus lines in conjunction with metro networks.Te key to the design of the integrated system of bus and metro is how to optimize the bus routes, make them adapt to metro lines, and form an organic connection in time and space, so as to shorten the transfer time of passengers, thus improving the transfer efciency between lines, enhancing the transportation capacity of the entire integrated transportation system, and providing urban residents with more convenient and efcient travel experience.Meanwhile, it also helps in reducing trafc congestion and environmental pollution [33][34][35].
In recent years, research on the joint optimization of bus and metro has attracted much attention, and many researchers have conducted relevant research on the design and optimization of bus routes connecting with metro stations.Zhu et al. [36] established a route model for feeder bus service with the objective of maximizing potential demand and used genetic algorithms to generate optimal routes that start and end at rail transit stations.Wei et al. [37] proposed a feeder bus route positioning process in order to attract the potential travel and transfer demand by analyzing the attraction range of the metro station.In the process of bus route optimization, in addition to considering passenger demand, some researchers took the interests of both passengers and operators into account.Lai et al. [38] proposed a multipleto-multiple optimization model between metro stations and feeder bus stops in the feeder service area with the objective of minimizing feeder bus operating cost and total passenger travel time; however, they overlooked transit transfers.Cao et al. [39] took into consideration the impact of passenger transfer on their travel time cost and established the optimization model in view of transfer time and the number of transfers.In recent studies, the integrated model connecting bus and metro lines is established to tackle both route optimization and vehicle scheduling problems.Xiong et al. [40,41] studied the optimization problem of community bus routes connected to metro stations with the objective of minimizing the total cost of both operator and passenger costs.Zhang et al. [17] studied long through-type bus lines near metro stations and established a bilevel programming model considering passenger and operator cost, as well as passenger fow allocation.Te experimental results show that the operating cost could be evidently reduced and the passenger fow balance coefcient of long-distance bus lines near metro stations could be improved.
Te multimodal public transport network realizes the organic combination of metro, urban rail, bus, and other transportation modes, which can expand the accessibility of rapid transit lines and improve the operating efciency of public transportation [42].Some researchers have studied the optimization of bus routes in the integrated system of bus and metro.Synchronization of the feeder bus network system and rail transit system was considered by Almasi et al. [43,44], and a model with the goal of minimizing the total operational cost was established.Costescu et al. [42] considered not only the minimization of the operational cost but also the minimization of passenger travel time in the bus route planning model.Huang et al. [45] proposed a threestage design method for the integrated bus and metro system.Hub stations and radiating nodes were selected from a clustering method, and a bilevel programming model was established with the minimum system cost to obtain the optimal passenger route.In the joint optimization of rail transit and bus networks within the transportation corridor, Chien et al. [46] gave a classifcation method of bus routes under the condition of rail transit and established a model to optimize the bus headway, bus stop location, and bus route distance, whereas Sun et al. [35] built a multiobjective model with the objective of maximizing rail trafc fow and minimizing total passenger travel time.
In studies on joint optimization of bus and metro networks, not all bus routes need to be optimized and only a few scholars have mentioned how to select bus lines that need to be optimized.Wei et al. [d7] used the smart card data in the urban public transportation system to generate a passenger origin-destination (OD) matrix and determined the routes and stops to be adjusted by analyzing boarding and alighting numbers at each stop.Cui et al. [48] determined potential bus lines to be optimized by evaluating the utility of colocation segments of bus lines and metro networks.In addition, some scholars selected bus lines along metro lines that need to be adjusted based on passenger travel preferences.Cao et al. [49] determined the bus line to be optimized by establishing a travel choice logit model based on residents' travel preferences.
In the abovementioned studies on combined bus lines and metro networks, few studies focused on the optimization and adjustment of ultralong bus lines.To improve the transportation efciency of ultralong bus lines and reduce the negative impact of factors such as bus bunching and low punctuality rates, a multiobjective optimization model with the aim of minimizing passenger travel cost and operator cost, respectively, is established to optimize bus lines combined with metro networks and selected based on IC card data.NSGA-II is used to solve the model to obtain the optimal truncated stop, feet size, vehicle type, and schedule.

Te Determination of Ultralong Bus Lines.
Te frst problem that needs to be addressed is how to specifcally identify ultralong bus lines from numerous regular bus lines before optimizing procedures.In this study, ultralong bus lines are determined using (1) from the standard [1].If a bus line's length exceeds the upper limit predefned in (1), it is defned as an ultralong line.
where L max is the theoretical maximum length of urban bus lines (km).When a city's regular bus line crosses this threshold, it is considered to be an ultralong bus line.S n is the the average running speed of buses in the line n (km/h).
T max is an average maximum one-way travel time for 95% local residents using public transit.For megacities, such as Beijing and Shanghai, the value can be set as 60 minutes [1].

Te Determination of Ultralong Bus Line to Be Adjusted.
Te ultralong bus lines that are attempted to be adjusted in this study are lines not only constrained by length limit in (1) but also with low punctuality rate and frequent bus bunching.Te bus bunching phenomenon refers to the meeting of two or more vehicles on the same bus route in operation, including the meeting at one stop or between bus stops.Meanwhile, this study discusses bus bunching at stop in the following two situations: (a) the bus j has not left stop k when j + 1 bus arrives at the stop k in the bus line and (b) j + 1 bus leaves stop k earlier than bus j.Te bus bunching rate is defned as the ratio of the number of vehicles occurring bus bunching to the total number of buses within a given time interval.A line is regarded as a high-frequent bunching line when the bunching rate is more than 5% [50].
In addition, the punctuality rate refers to the probability that a vehicle arrives at its scheduled stop within a given time window.Kho et al. [3] defned bus punctuality as the adherence to on-time performance and headway and proposed three related indicators to determine the punctuality rate.In this study, we say a vehicle is punctual when the vehicle arrives at the target stop less than 5 minutes early or late compared to the prescribed time.In this study, we defne a bus route is not punctual when its punctuality rate is less than 60%.

Model Formulation
4.1.Assumptions.For the convenience of model development, the following assumptions are proposed: A1.In the optimization process, the total passenger capacity of the network is set to a fxed value A2.Te truncation stop of the bus line is located near the metro stations A3.Passengers can only transfer to bus or metro A4.In the original bus line, 70% of the passengers passing through the truncated stop transferred to metro and 30% to the bus A5.Ignoring the deceleration and acceleration of the bus as it arrives and leaves the station, the average speed is taken as the travel speed A6.It is assumed that the time taken by the passengers in the boarding and alighting process is neglected; the time taken by the vehicle to reach the station is the boarding and alighting time of the passengers 4.2.Notations.Te parameters and variables used in the model are defned as shown in Table 1.

Problem Description.
In the study, optimization of an ultralong bus line relies on metro stops nearby.As shown in Figure 1, bus line #52 is an ultralong line, which intersects with several metro lines, and the intersection points are several stops, which can be used as potential truncation sites.Te ultralong bus line is divided into several segments by these potential truncation sites.Te main research in this paper is how to select the optimal truncation stop from the set of potential truncation stops, and how to obtain the optimal vehicle scheduling after the line is truncated at the site.

Objective Function.
Te objective function is a dual objective, where (2) and (3) represent the minimization of total costs for passengers and bus companies, respectively.
min P for each passenger's metro and bus trips can be expressed as follows: In ( 4), the frst term is the ticket purchase cost for all passengers on the bus, including the frst boarding and bus transfer; the second term is the ticketing cost of transferring to the metro at the truncation stop.
Te set of passengers at the truncation stop k transferring to a bus or metro can be expressed as follows: Equation ( 5) represents the set of passengers who pass through a truncated stop and transfer to the bus at that stop with a probability greater than 50%.Equation ( 6) expresses the set of passengers who pass through a truncated stop and transfer to the metro at that stop with a probability greater than 50%.
Te cost of passenger time including in-vehicle time and waiting time cost can be expressed as follows: In ( 7), the frst term represents the bus ride time cost, which includes the waiting time cost and in-vehicle time cost before the splitting stop.Te second term expresses the waiting time cost and in-vehicle time cost after the passenger transfers to the bus or metro at the splitting stop.Te number of passengers in bus j at stop k (person) 6 Journal of Advanced Transportation

Journal of Advanced Transportation
Te frst and second items of (8) mean that passenger i traveling by bus j boarding and alighting at stops k 1 and k 2 , respectively; the third and fourth items calculate the travel time before and after the splitting stop.Equation ( 9) calculates the passenger's waiting time based on the headway of the lines.Equation (10) calculates the travel time of passengers on the metro mainly based on the average travel speed of the local metro line.Te waiting time for passengers when they transfer to the subway is expressed by the authors in (11).
Te total penalty cost is comprised of the following three costs: the penalty cost for an additional transfer brought on by line truncation, the penalty cost for a long wait, and the penalty cost for in-vehicle congestion.
In ( 12), the frst term is the additional transfer penalty cost for passengers, which depends on the passengers passing through the truncated stop.Te second term is the congestion penalty cost, when in-vehicle passengers exceed the expected load value τV y max .Generally, τ relies on how passengers actually feel.Te last term is the waiting time penalty cost, when passengers' waiting time exceeds their expected value.

Bus Company Cost.
Te cost C 2 of the bus company consists of a fxed cost C F and an operating cost C V , as shown in (3).
Te frst term in (13) represents the cost of vehicle purchase, which assumes that the wage of each driver depends on the type of bus.Te second term expresses the labor cost, which relies on the number of drivers deployed on the route.
Equation ( 14) calculates the number of buses of type y needed in line n, and the number of drivers equipped is equal to the number of buses.Te operating cost of the bus company can be expressed as follows: In ( 15), the operating cost of the bus company depends on the length of the operating line L n and the vehicle type y.
Te length of the bus line is calculated by (16), which is mostly determined by the position of the truncation stop, the number of stops, and the average stop distance.(16) defnes that when only one truncated stop k exists in the line, namely, x k � 1, the upper n � 1 represents the part from the origin depot to the truncated stop and the lower n � 2 is the other part from the truncated stop to the terminal stop.

Constraints
Te waiting time limit for each passenger is represented by formula (17).Te limit number of truncated stops is expressed by formula (18).Formula ( 19) represents the bus capacity constraints, where V j k is determined by the number of passengers alighting and boarding at each stop.Equation (20) calculates the number of passengers in the bus j when it arrives at stop k in the ascending or descending direction, which is mainly determined by the number of passengers boarding and alighting at each stop.Formula (21) refects the headway limit, which is primarily determined by the maximum passenger fow in an hour and bus capacity, as well as applicable restrictions such as minimum and maximum headway limits.Formula ( 22) is the constraint for the number of truncated stops derived from the metro stop set.Te upper and lower limits of truncated stops can be represented by equations ( 23) and (24).Formula (25) expresses the quantity constraint of each vehicle type in each line.Equation ( 26) calculates the maximum number of buses allocated to the bus line determined by the length of the line, the average speed, and the headway.Formula ( 27) represents a constraint on the 0-1 variable.

Algorithm Solution
Te solution of the model is primarily divided into two steps.First, alternative truncation stops that satisfy relevant constraints are identifed as potential truncation stops.Second, an algorithm is employed to determine the optimal headway, feet size, and schedule for each alternative truncation stop.Te target line optimization process is shown in Figure 2.
Te data utilized in this research encompass the arrival time information of each passenger.Te algorithm's variables, including passengers' waiting time, maximum passenger count, and vehicle type, are determined by the vehicle's headway as depicted in Figure 3. Consequently, headway serves as the fundamental variable to be resolved, infuencing both other variables' outcomes and optimization performance.
Bus network design and frequency setting are nondeterministic polynomial (NP-hard) problems, and heuristic algorithms have obvious advantages in solving these problems [51].Te improved NSGA-II based on the genetic algorithm can efectively solve the difcult problem of weight setting in multiobjective problems [52].Among them, fast nondominated sorting, crowding distance assignment, and building partial ordering sets, which are all infuenced by the quantity of targets and individuals, are the key factors that afect the computing efciency of NSGA-II.In the algorithm, it is imperative to continuously calculate how many people board and alight at each stop with the constraints of a problem, the program's space complexity is large, as well as constraints on its computing efciency [53].In summary, NSGA-II is used to solve a multiobjective mathematical programming model for the optimization of ultralong bus line in this study, which yields a Pareto optimal solution through the process of fast nondominated sorting.

Chromosome Representation.
In this study, each chromosome represents the headway of two subroutes after being truncated.An 8-bit binary encoding is used to keep the headway of each line within the range of 8 to 15 minutes.Since the frst and ffth positions are set to fxed values, these two genes in the chromosome are ignored, as shown in Figure 4.

Journal of Advanced Transportation
Te crossover operation is shown in Figure 5.A partial segment crossover operation is utilized involving exchanging the two segments from two parent chromosomes to create new individuals according to a predetermined probability.
Te process of mutation is shown in Figure 6.
Te selection method is tournament selection, as shown in Figure 7.

Procedures. Te algorithm fow is as follows:
Step 1: set the algorithm parameter values, including the number of individuals N in the population, the maximum number of iterations Q, the probability of crossover and mutation, and the iteration number t � 0 Step 2: initialize, encode the individuals, and generate populations Step 3: calculate the ftness values of the individuals in the initial population Step 4: check each individual whether it satisfes the constraints; if not, its ftness value is set to infnity Step 5: the individuals in the population are sorted by fast nondominated sorting according to their ftness values Step 6: the crowding distance between individuals is calculated to measure the quality of each solution in the same Pareto front Step 7: using elitism, according to the results of fast nondominated sorting and crowding distance calculation, the individuals in the population are sorted and the ones are selected Step 8: generate a new population through crossover and mutation Step 9: calculate the ftness values of all individuals covering the parent population and the ofspring population, and then set, t � t + 1 Step 10: if t < Q, then return to Step 4; otherwise, end.Te optimal cost of each stop is compared to select the optimal truncation stop Te algorithm solves the optimal departure interval and cost for each alternative truncation stop Te best individual is selected from the Pareto front solutions as a representative in each iteration and it plots the iteration curve as shown in Figure 8 (only S15 is shown as an example).Te horizontal axis represents the number of iterations, and the vertical axis represents the minimum ftness value obtained by the weighted summation of the two objective function values of each individual on the current Pareto frontier.Te current optimal solution is determined by normalizing the objective values of all options in the Pareto frontier and weighting them.Terefore, the solution with the lowest value is then considered to be optimal.
Figure 8 demonstrates that the method reaches the convergence stage after 17 iterations when the weighted sum value of the ftness of the best individual no longer varies.

Case Study
In this case study, an ultralong bus line is determined based on public transport IC card data from Beijing in 2018, and then the selected bus line is utilized to show the optimization results and applicability of the proposed strategy on the ultralong bus line.

Determination of the Ultralong Bus Line.
Te IC card data mainly include bus number, line number, passenger boarding, and alighting time and stop.Te specifc felds are shown in Table 2.
In this study, Beijing Bus #52 is selected as the target line.Te total length of this line is 22 km, and the average bus speed is 16.9 km/h.Te population of Beijing exceeds 20 million, so T � 60 minutes, and the maximum length limit of the line can be calculated as 16.9 kilometers according to formula (1).After processing and analyzing the passenger fow data, it is found that numerous bunching phenomena exist on this line, and its punctuality rate is relatively low.Te partial bus bunching phenomenon will afect the overall service efciency of the whole public transport network, resulting in a serious imbalance of in-vehicle passengers and waiting time on several bus stops, as well as a low punctuality rate.
In Figures 9(a) and 9(b), 8 vehicle trips are listed both in the up and down directions between 7: 00 a.m. and 8: 00 a.m. at each stop in #52.Te bunching rate in the up direction is up to 44.1% and the punctuality rate is 50.0%, while in the contrary direction, the bunching rate is up to 58.5% and the punctuality rate is 43.9% during the peak travel hours from 7: 00 a.m. to 9: 00 a.m. and 5: 00 p.m. to 7: 00 p.m. Consequently, urgent adjustments to the bus line are necessary to address the frequent bunching and low punctuality rate during this time period.

Results' Analysis.
In this paper, #52 is taken as a specifc example, and the main parameters are presented in Table 3. Te distribution of metro lines and stops along #52 is shown in Figure 10.
In Table 3, the vehicle costs are calculated on a daily basis and determined by the purchase price and service life of each vehicle.Te fare of #52 is 2 CNY for 10 kilometers and adds 1 CNY for an additional 5 kilometers.In 2018, the per capita GDP in Beijing was approximately 140,000 CNY, so the travel time cost for residents during that period is about 34.5 CNY/h [54].Te selection of truncation stops in this study is primarily based on the analysis of passenger fow patterns and the presence of subway stops under specifc constraints.
each gene into a relative value according to the given probability, such as 0-1 and 1-0    In Figure 10, S13, S15, S16, S17, and S18 have been identifed as potential candidates for truncation stops.A detailed analysis of these candidates will be presented in the subsequent section.

Afer mutation
6.2.1.Total Cost Analysis.Te paper calculates the total cost for passengers and bus companies before and after splitting at bus stops S13, S14, S15, S16, S17, and S18.Tis information is presented in Table 4.Among these stops, S13, S14, S15, S17, and S18 are close to metro stops, while S16 is a nonmetro stop, meaning it is a bus stop without a metro station nearby.
In Table 4, compared with the average cost value at splitting metro stops, the total travel cost for passengers after the truncation operation is not signifcantly changed, but the total cost of bus companies is reduced by 65.1%.Tis indicates that the splitting strategy operation is efcient for the ultralong bus line.Compared with nonmetro stops, the average travel cost is decreased by 12.0% at metro stops, but no signifcant change in the bus company cost occurs.Moreover, the costs of both the bus company and passengers at S14 and S15 are relatively lower than those at S13, S17, and S18.Terefore, the optimized truncation stop is chosen from S14 and S15 in order to obtain better performance.After truncation, compared with the average value of splitting metro stops, the average passenger ticket purchase cost increases by 39.2% due to more expensive metro purchases, and the time cost decreases by 5.3%.Also, the reduction in travel time costs means less travel time is taken for passengers with slightly higher ticket costs.Since the average ticket purchase fee per passenger increased by 0.88 CNY, the maximum savings for passengers is about 20 minutes in a trip.Te results will be much better if the travel fee has a discount or the transfer fee is free, as seen in the analysis in the next section.However, for the nonmetro truncated stops, the ticketing cost is reduced by 8.6%, but the time cost and penalty cost are increased by 14.4% and 74.5%, respectively, compared with the average value.It is indicated that if the truncation stop is a nonmetro stop, passengers still have to transfer to bus travel, which may increase extra transfer time without much travel time savings.

Bus Company Cost Analysis.
Te same conclusion can be drawn from Figures 12 and 13 that after splitting both the operation cost and fxed cost reduce to relatively low values, 41.0% and 70.0%, compared to average values, 40.0% and 67.6%, compared to nonmetro stops.It may be because the length of the bus line becomes shorter after the truncation operation, and the required feet size on the line is reduced.At the same time, with the adoption of diferent bus types, both the fxed cost and the operating cost are reduced.
Compared to the nonmetro stops, the average fxed cost and operating cost of splitting stops near metro stops decrease by 1.7% and 7.3%, respectively.Figure 14 shows that the passenger fow at metro stops is signifcantly lower than the fow at nonmetro stops.It may indicate that part of the passengers will transfer to the metro, thus reducing the number of passengers on the new split line.However, when the truncated stop is a nonmetro stop, passengers can only transfer to the split bus line, resulting in a relatively large amount of passenger fow on the latter line.Terefore, it is necessary to reduce the headway and increase the number of vehicles in order to deliver more transfer passengers.
From the perspective of total cost, the optimized truncation stop is attempted to be generated from S14 and S15.Te costs for both passengers and bus company at S15 are relatively lower than those at other stops.Terefore, S15 is determined to be the optimal truncation stop, and the relative line information of the two parts is shown in Table 5.Part of the timetable after truncation is shown in Table 6, and the complete information is represented in the attachment.
6.3.Sensitivity Analysis 6.3.1.Impact of the Vehicle Type.In this paper, two types of buses with diferent capacities are utilized.Among them, the results of the bus company cost splitting at S15 stop with diferent bus type combinations are detailed in Table 7.In Table 7, the bus company cost for the combination of largeand small-capacity vehicle types is relatively lower than the other two single types.Although the vehicle purchase cost of the small capacity type is lower than the cost of the largecapacity type, more vehicles and drivers with a much higher purchase cost and driver salary are needed for this type.However, the load rate of vehicles with larger capacity is generally lower, especially during of-peak hours.In Figure 15, during the period from 9: 00 to 10: 00, the minimum load rate of the large-capacity bus type is 21.7%, and the average load rate is about 36.4%, which will lead to a waste of resources and a decrease in the transportation efciency of the whole bus system, violating the concepts of energy saving, emission reduction, and low-carbon travel.Terefore, the combination of large and small bus types can not only reduce the operating cost of the bus company associated with meeting the travel needs of passengers but also contribute to alleviating road trafc congestion and improving the transportation efciency of public transportation to a certain extent.

Results of Diferent Travel Mode
Expectations.Tis paper starts by assuming that 70% of the passengers choose to transfer to the metro and 30% choose to transfer to the bus at the truncated stop.Due to the diference in passengers' preferences for public transport, the proportion of transfers to the subway or bus may change.Terefore, three values of proportion are attempted to be utilized to verify the impact on passenger travel costs, as shown in Table 8.Te total passenger fow of #52 passing through S15 is 8021.Based on the results, as the number of bus-transferring passengers continues to increase, the total ticket purchase cost is gradually decreasing, caused by the fact that metro tickets are generally much more expensive than bus tickets.For 14 Journal of Advanced Transportation instance, in Beijing, the average cost per ticket for a metro ride is 3 CNY/person, whereas the average cost for a bus ticket is 2 CNY/person.However, the ticket cost will fall sharply if the transfer fee has a discount or is free for a period so as to encourage passengers to choose transit trips, which is now implemented in several cities in China such as Chengdu and Wuhan.In addition, the travel cost of passengers is decreasing, especially the travel time cost.It is obvious that transferring to the metro can signifcantly save passengers' travel time and improve the efciency of the public   In order to reduce the loss of time and money for the passengers when facing additional transfers, several preferential measures, such as transfer discounts or free in a period, are attempted to be adopted.With these measures in place, the cost of passenger travel is shown in Table 9.As can be seen from Table 9, compared to not taking any preferential measures, when transfer fares are discounted by half, the cost of passenger tickets is reduced by 17% and the total cost of travel is reduced by 6.2%; when transfer fares are free, the cost of passenger tickets is decreased by 34% and the total cost of travel is decreased by 9.3%.Terefore, in the truncation and optimization of ultralong bus lines, the adoption of transfer concessions at the truncated stops can efectively reduce passenger cost and the resistance generated by transfer.

Conclusion
Tis study develops a multiobjective model for optimizing the operation of long bus routes by truncating them.Te objectives are to minimize passenger travel costs and minimize bus company costs.Te model incorporates the infuence of subway networks and adopts a coordinated optimization approach between bus routes and subway networks.Furthermore, it considers the operational mode of coordinating multiple types of buses in order to enhance the   Journal of Advanced Transportation operational efciency of the bus company.Based on an analysis and verifcation using Beijing's #52 bus route, several key conclusions can be drawn.
(1) By implementing a truncated operation method for ultralong bus lines, signifcant travel time savings can be achieved at the expense of a slight increase in purchase ticket cost.In this scenario, the average ticket cost after truncation is raised by 0.88 CNY/ person, while passengers experience a decrease of 5.8% in overall travel time costs and the bus company achieves up to a 64.9% reduction in total costs.However, it should be noted that the maximum possible reduction in travel time after optimization is 20 minutes.Terefore, it is recommended to focus on improving ultralong lines with minimal operational benefts.
(2) Te selection of a splitting stop depends on its ability to signifcantly reduce costs for both bus companies and passengers.Based on this analysis, locating splitting stops near metro stations can efectively decrease the aforementioned expenditures.Terefore, in order to maximize the synergistic benefts between these two modes of transportation, it is essential to consider the placement of splitting stops along ultralong bus routes in relation to nearby metro stations.(3) Combining large and small capacity bus types in an optimization model's operation can reduce the bus company's operating costs while simultaneously increasing the full load rate and lowering certain carbon emissions and trafc congestion.
(4) Considering the additional transfers following the route truncation, it is advisable to implement preferential measures such as fare discounts or complimentary transfers in order to mitigate ticket purchasing costs.Tese measures can efectively address the operational defcit resulting from truncation and enhance the attractiveness of public transportation.
Although this work makes a fresh attempt to enhance the optimization of ultralong bus lines in both theory and reality, there are still a number of issues that need to be addressed.(1) Regardless of a person's preferred mode of transportation, a predetermined and fxed percentage of passengers will take a bus or metro.In the future study, a model of individual travel mode selection will be taken into account; (2) after truncation, the efect on passenger choice or the overall passenger fow of the bus trip is not taken into account.In order to gather pertinent data for the ensuing study, an SP survey based on the operation tactics on the extremely long bus route will be conducted; (3) on the ultralong bus line, only one line-truncated optimization technique is investigated.More tactics, including linetruncated, skip-stop, and line-shortened ones as well as combination strategies, will be used in the future study to optimize ultralong lines; (4) the implementation scope is intended to be broadened to include all relevant transit networks, not only the one bus or metro line at the splitting stop.

Figure 1 :
Figure 1: Part of #52 bus line and metro network.

Figure 8 :
Figure 8: Iteration curve: the horizontal axis represents the number of iterations, and the vertical axis represents the minimum ftness value obtained by the weighted summation of the two objective function values of each individual on the current Pareto frontier.

6. 2 . 2 .
Passenger Cost Analysis.Various types of passenger travel costs before and after truncation are shown in Figure 11.

Figure 11 :Figure 12 :Figure 13 :Figure 14 :
Figure 11: Various travel costs for passengers include passenger ticketing costs, travel time costs, and penalty costs at each truncated stop, respectively.

Figure 15 :
Figure 15: Maximum load ratio at each time points for two diferent feet types.

Table 1 :
Continued.Number of total passengers in an hour of line n (person/h) I kTe total number of passengers transferring at stop k (person) IB kTe total number of passengers transferring to bus at the truncation stop k (person)Binary variable, 1 if the vehicle j of type y is in the line n, and 0 otherwise

Table 2 :
Partial felds of passenger travel IC card data.

Table 4 :
Cost results at diferent truncation stops.

Table 5 :
New line information.

Table 7 :
Cost of bus company when using a single model.

Table 8 :
Passenger travel costs under diferent transfer proportion.