CN110598942A - Community public transport network and departure frequency synchronous optimization method considering area full coverage for connecting subways - Google Patents

Community public transport network and departure frequency synchronous optimization method considering area full coverage for connecting subways Download PDF

Info

Publication number
CN110598942A
CN110598942A CN201910880061.4A CN201910880061A CN110598942A CN 110598942 A CN110598942 A CN 110598942A CN 201910880061 A CN201910880061 A CN 201910880061A CN 110598942 A CN110598942 A CN 110598942A
Authority
CN
China
Prior art keywords
line
path
network
scheme
operator
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910880061.4A
Other languages
Chinese (zh)
Other versions
CN110598942B (en
Inventor
熊杰
陈彪
陈艳艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201910880061.4A priority Critical patent/CN110598942B/en
Publication of CN110598942A publication Critical patent/CN110598942A/en
Application granted granted Critical
Publication of CN110598942B publication Critical patent/CN110598942B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047Optimisation of routes or paths, e.g. travelling salesman problem
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40Business processes related to the transportation industry
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Train Traffic Observation, Control, And Security (AREA)

Abstract

The invention discloses a community bus network and departure frequency synchronous optimization method for a connection subway considering area full coverage, and belongs to the field of public transport network planning. The invention divides the research area into a plurality of latticed cells, takes the minimization of the total cost as the optimization target, and also considers the problems of the layout of the micro-circulation public transportation network of the community and the setting of the departure frequency of the line, thereby constructing a mixed integer optimization model. In the solving process, the invention provides a series of constraint conditions based on reality, such as: line length constraints, net total service area constraints, vehicle capacity, fleet size constraints, and the like. The invention designs a plurality of crossover operators and mutation operators to optimize in a network scheme population; the invention also provides an adaptive weight adjustment algorithm, which can adjust the probability of executing the next generation of operators according to the execution condition of the genetic operators so as to improve the execution efficiency of the genetic algorithms.

Description

Community public transport network and departure frequency synchronous optimization method considering area full coverage for connecting subways
Technical Field
The invention belongs to the field of public transport network planning, and relates to a community public transport network and departure frequency synchronous optimization method for a connection subway considering area full coverage. The invention is suitable for optimizing the public transport network in medium and small-scale road network areas.
Background
With the development of urban traffic, conventional buses play an important role in urban passenger transport, but due to the limitation of vehicle types, the bus type is difficult to well play a connection function in an area with sparsely distributed subway stations (particularly in suburbs), and the micro-circulation buses in the community can more flexibly run in the area, so that more convenient and faster service can be provided for passengers. In the traveling process, passengers may need to transfer for many times to travel in a public transport-subway mixed transportation mode, and for better coordinating with the subway to complete passenger transport tasks, a bus operator needs to continuously adjust a microcirculation public transport network.
Due to the lack of sufficient basic theory support, the current microcirculation public traffic line network arrangement method also mainly depends on the experience of planners. Therefore, the invention aims to design an optimal community microcirculation public transportation network from a theoretical level, so that the optimal community microcirculation public transportation network plays an important role in subway connection, and simultaneously realizes the advantage complementation of two public transportation modes, reduces the total cost of a public transportation system and improves the total operation benefit.
Disclosure of Invention
A community bus net and departure frequency synchronous optimization method considering area full coverage for connecting subways comprises the following specific steps:
the method comprises the following steps: optimization model construction
Based on realistic constraint conditions, corresponding penalty functions are provided for the line length, the vehicle load and the service level, and a mixed integer optimization model is constructed by taking the minimization of the total cost of the system as an optimization target.
To more fully describe this problem, the present invention proposes the following constraints:
constraint 1: and (4) restricting the full coverage of the cell. In order to improve the accessibility of passengers in the community for traveling, the service range of the public traffic network needs to cover all the square cells. That is, each grid cell is within the maximum walking distance tolerable to passengers, i.e., each cell is within the service range of at least one bus stop or station in the network, and if the station or station closest to the cell exceeds the tolerable maximum walking distance, the tolerable maximum walking distance is regarded as the distance between the cell and the station.
Constraint 2: and (5) full coverage constraint of the subway station. Subway stations in a research area may not be all distributed on one subway line, and each node in the ST appears in the line network at least once to ensure the connection of two public transportation modes. Meanwhile, all nodes in S do not necessarily all appear in the public transportation network.
Constraint 3: total fleet size constraints. For a community, the public transportation line is often operated by the same operator, and the total fleet size of the public transportation network operation cannot exceed a given upper limit due to limited capital investment of the operator.
Constraint 4: vehicle capacity limitations. On each line, the section traffic volume per trip must not exceed the maximum traffic volume of the operating vehicles, which requires that each line must be set to a minimum departure frequency to ensure the corresponding service level.
Constraint 5: a line length constraint. According to the scale of the road network, the length of each line should be within a reasonable range, which is one of the most common constraints in the research.
Constraint 6: and (5) line trend constraint. There cannot be duplicate stations in each line to avoid a turn-around or local loop in the line.
For ease of study, the following assumptions were made: (1) the study area is divided into several equal sized square cells as in fig. 1. The public transportation requirements of each district are different and are uniformly distributed in the district; (2) the starting point/end point of each OD trip is at the center of the cell, and the end point/starting point is at a subway station (in reality, the trip demand among bus stations in a suburb community is negligibly small, especially in peak hours); (3) regular bus speed and station stop time are constant, since planning is the focus of the study, not operation; (4) for convenience of calculation, the walking distance for passenger transfer is 0; (5) all bus stations and subway stations are provided with conditions for regular bus turning; (6) node 0 represents a bus stop where passengers cannot take a ride or transfer, while all other stops can provide ride and transfer services.
The objective function of the present invention mainly includes user cost and operation cost, and can be expressed as:
CT=CU+CS (1)
in the formula, CTIs the total cost; cUThe cost for the user; cSIs the operating cost.
The user cost can be divided into three parts: cost of passenger walking CAPassenger cost in-vehicle CIPassenger waiting cost CW:
CU=CA+CI+CW (2)
In the formula, CA、CIBoth with respect to the net layout and the path selected by the passenger, they can be expressed as:
in the formula, Z is a square cell set; pathodIs a selectable path set between o and d; dod kThe walking distance from the cell o to the corresponding getting-on/off station on the kth path between o and d; l isod kIs the on-vehicle distance in the kth path between o and d; t is tsDelaying the time for the station; nsod kThe number of sites of the kth path between o and d; q. q.sodThe number of passengers between o and d (or d and o); p is a radical ofod kProbability of selecting the kth path for the passenger; vA、VBThe walking speed and the vehicle speed of the passenger; mu.sA、μIThe walk time value and the vehicle time value.
In formula (3), Dod kIndicating the walking distance from the area o to a certain bus stop on the k-th route between o and dAnd the station is within the maximum walking distance range that the residents of the cell can tolerate. If there are a plurality of stations on the route within the range, the residents are matched according to the minimum sum of the walking time and the in-vehicle time. In particular, for a cell close to a subway station, D may be within the range of cell o, which means that the resident will walk directly to the subway station instead of riding a regular bus, for which case the cell is matched directly to the subway station, Dod kRepresenting the walking distance from o to d.
The average waiting cost of passengers on a route is mainly related to the departure interval of the related line. That is, if the route passes through two or more lines, the waiting time for transfer is also included in CWIn (1). Thus, CWCan be expressed as:
in the formula, R is a line set of a line network scheme; nsrThe number of bus stations/subway stations (except bus stations) at which regular buses need to stop in one-way driving of the vehicles; HBrThe departure interval of the line r; nbs rThe number of passengers boarding route r at station s; z is a radical ofwA ratio of arrival time/departure interval desired for the passenger; mu.sWThe value of waiting time is obtained.
The operating cost is related to the net layout and departure interval for each line, which is equal to the total fleet size multiplied by the bus operating cost, and the fleet size for each line can be further expressed as equation (7).
In the formula, B is the operation cost of the vehicle per hour; frIs the fleet size of line r. L isrIs the unidirectional length of the line r.
A decision variable 0-1 is also determined, and if node j is a node subsequent to node i in line r, j belongs to SST, R belongs to R; if not, then,thus, Lr、nsrCan be further expressed as:
in the formula Iij rIs the distance on line r from i to j.
In the above equation, xij r、HBrAre the decision variables that need to be solved. For a given net plan and departure frequency set, Fr、Lr、nsrCan be directly obtained by calculation; and Dod k、Lod k、pod k、Nbs rNeeds to be determined by the passenger flow allocation method.
Because the line length repairing method can not ensure that each unqualified line can meet the length limit through repairing, a line length punishment is made, CP1As follows:
in the formula, muP1Are coefficients.
Considering the vehicle capacity constraint as a flexible constraint, an overload penalty is set. Before calculating the penalty, a penalty function is defined with respect to the vehicle capacity, g (x):
g(x)=ω(αx-1) (11)
in the formula, x is the cross section passenger flow volume; and omega and alpha are coefficients. By adjusting ω and α, g (x) approaches 0 when x is 0. ltoreq. P, and x approaches 0>P, rises rapidly. Overload penalty based on equation (21), net scheme, CP2Expressed as follows:
in the formula, Ldij rDetermining the section passenger flow quantity between the nodes i and j in the line r per hour according to the passenger flow distribution result; mu.sP2Are coefficients.
Under each net scheme, the bus service may not meet the travel requirement of each part, that is, there is no route from the starting point to the end point in the current net, or too many transfers are needed in the route (when the route needs more than two transfers, the current travel is regarded as being out of service). Thus, another penalty is built for demand that cannot be serviced, CP3As follows:
in the equation, if the current net can serve the travel demand, δod1, otherwise, δod=0;μP3Are coefficients.
C is to beP1、CP2、CP3Added to the initial objective function, the objective function can be expressed as:
min CT'=CA+CI+CW+CS+CP1+CP2+CP3 (14)
step two: net initialization
And generating a net scheme set meeting all the constraint conditions in the model by using a net initialization algorithm, wherein the net scheme set also comprises a sub-algorithm, namely a line length repairing algorithm.
The randomly generated net schemes can all meet the constraint condition. For the constraint of the line length, a length repair algorithm is provided, so that the length of each line is in a reasonable range. When the area full coverage constraint is realized, a corresponding algorithm is designed, namely, a new line is repeatedly generated to serve the current farthest uncovered area until all areas are covered. Finally, the number of subway stations is checked to ensure that all stations are connected into the network.
For the cell Z belonging to Z, the bus station S belonging to S and wdzsDenote the walking distance from z to s, let wdmin z=min{wdzs| S ∈ S }, Z0 denotes an uncovered cell, let wd bemax=max{wdmin z| Z ∈ Z0}, ST0 denotes the set of subway stations not connected to the net, NS denotes the currently generated net solution, WD denotes the maximum walking distance tolerable for the passenger, and the matrix RN denotes the connectivity of every two nodes of the road topology based network.
The steps of randomly generating a feasible net scheme are as follows:
in the above algorithm, the shortest path between two nodes always needs to be found (steps 4, 7, 12), and to solve this problem, the Dijkstra algorithm is repeatedly called, and the reason why the shortest path is taken as a path between two nodes is as follows: 1) the generated line is connected with some key nodes, namely target nodes which want to join the network; 2) there is no additional detour on the line and the line length generally does not exceed an upper limit, so line length repair-1 is primarily used to increase the length of shorter lines.
For line r, the length repair algorithm uses two ways to increase the length of the line: 1) for Lr<LminThe length of the line meets the requirement by extending the line; 2) for Lr≥LminAccording to a probability p associated with the current line lengthinsInserting some nodes into the line r:
the detailed steps are as follows:
the invention also provides another line length repairing method, which aims to shorten the longer line, namely line length repairing-2.
For line r, line length repair-2 shortens the line in two ways: 1) for Lr>LmaxThe shortest path between the nodes is used for replacing the sub-line, so that the line length is shortened; 2) for Lr≤LmaxAccording to a probability p associated with the current line lengthremDeleting nodes from r:
and pinsDifferent from, premThis operation results in longer line lengths, calculated by exponential functions, rather than linear functions, because a net consisting of a large number of short lines cannot be considered a good net solution.
The first line shortening method is that starting from a first node, detecting whether a sub-line between the node and other nodes is the shortest path, if not, replacing the sub-line with the shortest path; the second method for shortening the line is to detect whether the front and rear nodes of each node except the nodes at the two ends of the line can be directly connected (assuming that the triangle inequality rule is effective), and if so, deleting the middle node. Comparing the two methods, it was found that if a line cannot be shortened by the first method, it cannot be shortened by the latter, and that a line which cannot be shortened by the latter method can be tried to be shortened using the first method. Therefore, when L isr>LmaxThen, the first method is used to shorten the line, and the detailed steps are described as follows:
line length repair-2 is not suitable for the net initialization process because the initial line is generated mainly using the shortest path algorithm and each method of shortening the line may cause the loss of critical nodes on the line. However, it can be applied to the crossover operator, as will be described in detail in section 6.
Two line length repair heuristic algorithms ensure that each line is within a proper length range to the maximum extent, but the connectivity of a receiving network and the Lmin,LmaxAnd (4) calculating the penalty of the line length if the line still can not meet the length constraint, and adding the penalty into the objective function.
Step three: fitness evaluation
Determining alternative routes for passenger travel between each pair of ODs based on a public transit network, constructing a Logistic model taking walking time and vehicle time as effect functions, calculating the proportion of passenger flow on each route between the ODs, determining route departure frequency, providing a fleet scale adjustment algorithm to meet the limitation of the maximum fleet, and finally calculating an objective function value of a scheme.
PathodAnd (4) representing the path set with the least transfer between the o and the d, and determining the boarding station and the transfer station of the passenger through the matching relation between the path and the public traffic network, so that the cross-section passenger flow of each cross section on each line is obtained through calculation. Finding out all paths between each pair of ODs by using the method in graph theory, and judging the transfer times of each Path so as to determine Pathod
The current problem is how to quickly find all paths with the least number of transfers between each pair of ODs to determine the Pathod. In determining PathodIn the process of (1), firstly, whether the cell o is in the service range of the subway station d is detected, if yes, the Path between o and d is the shortest walking Path, and the shortest walking Path is recorded in the PathodSuch a demand does not require a regular bus ride (scenario 0); otherwise, whether a riding path without transfer exists between each pair of ODs (scene 1), a riding path with one transfer (scene 2), a riding path with two transfers (scene 3) is detected respectively,if there is no path not exceeding two transfers (scenario 4), this OD inter-pair demand is considered as an unserviceable demand, and the corresponding C will be calculatedP3And added to the objective function as shown in fig. 1. In fig. 1, there is no case where two scenarios occur simultaneously, and the occurrence of each scenario will terminate the algorithm.
In scenarios 1-3, PathodThe generation mechanism of (2) is similar, and the most significant difference between the three scenarios is that the number of cycles in the algorithm grows exponentially. In three scenarios, in order to remove too many detours and save the calculation time, the invention makes a screening strategy: will cost the effect more than PathodThe minimum value of all paths in the Path is 1.5 times from PathodIs deleted. Scenario 4 is used as a supplement to other scenarios, the routes of the scenarios exceed two transfers, the travel demand of the OD pair is identified as a non-service demand, and corresponding punishment is calculated.
In scenarios 2 and 3, the initial travel demand needs to be allocated to the bidirectional routes, since the routes are directional, which means that the directions of the sub-routes distributed on different routes are also considered while recording each route transfer station, so as to calculate the number of people getting on and off at each station during the round trip of the vehicle.
Due to PathodThe transfer time of all the paths is equal, the waiting time is not usually the main factor influencing the path selection of the passengers, and therefore, the proportion of the passengers assigned to each path is assumed to depend on the walking time and the time of the vehicles. Based on the above, the invention constructs a Logistic model taking walking time and vehicle time as effect functions, thereby calculating and obtaining the passenger flow proportion on each path between od:
in the formula, Vod kThe effect function for the kth path can be expressed as:
in the formula, a1,a2Are coefficients.
The cell o belongs to Z, the subway station d belongs to ST, NS represents the current network scheme, SSToRepresenting bus and subway stations, TRS, within the maximum walking distance of a cell oodRepresenting transfer stations for all paths.
Suppose that ns exists in a single pass of line rrThe round trip of a bus station/subway station (except for a station) can be divided into (2 ns)r-1) section forNbs rIndicates the number of passengers boarding the vehicle at station s, Nat rThe number of passengers getting off the vehicle at the station t can be expressed as:
in the formula (I), the compound is shown in the specification,all are variables from 0 to 1. If a passenger on the k-th route gets on board at station s on route r, thenIf the k-th path passenger transfers at station s on route r, thenIf the passenger on the k-th route stops at a stop t on route r, thenIf the k-th path passenger transfers at a stop t on line r, then
Passenger flow volume Ld of upper section of line rij rThe calculation method comprises the following steps:
setting a reasonable departure frequency for a route based on a particular fleet size is an important sub-problem in this study. In order to solve the problem, the invention provides a departure frequency setting method, which aims to determine the optimal departure frequency for each line and mainly comprises two stages: a preliminary calculation phase and a fleet adjustment phase.
In the preliminary calculation stage, for the route r, the waiting cost C of the passengerW rAnd an operating cost CS rCan be expressed as:
the departure interval of each line is in positive correlation with the waiting cost of passengers and in negative correlation with the operation cost, and the total cost C of the line r is solvedT rRelative to HBrThe first derivative of (a) yields:
the minimum C can be found by making the first derivative 0T rAnd corresponding departure interval HbrTo give formula (25):
upper limit HM of departure intervalrThe calculation method comprises the following steps:
optimized departure intervals, Hb, of line rr *The calculation method comprises the following steps:
round trip time of the line divided by Hbr *The value of (b) is rounded up to be the number of vehicles needed by the line r, namely the scale of the line fleet, and the calculation method comprises the following steps:
corresponding departure frequency fr *The calculation method comprises the following steps:
for all lines s ∈ R, Fr *The sum of (a) may exceed the fleet size upper limit TF, and when this occurs, the fleet size of the lines needs to be adjusted, so that the total fleet size meets the constraint, i.e., the fleet adjustment phase, by reducing the number of vehicles on certain lines. This may lead to some cross-sectional overloads, especially for Hbr *=HMrThe line of (2), at which time an overload penalty C needs to be calculatedp2And adds it to the objective function.
One net scheme includes a line set R and a fleet size set F ═ F1,…Fr,...,F|R|},ΔCT rRepresents the total cost change of the route r after removal of the vehicle, which can be expressed as Δ CT r=ΔCW r+ΔCS r+ΔCP2 r(ΔCW r,ΔCP2 r>0,ΔCS r<0) The algorithm is described in detail as follows:
step four: the design of an optimization algorithm is carried out,
the invention optimizes the line network by using the genetic algorithm, provides a plurality of crossover and mutation operators in the genetic algorithm, and embeds an adaptive operator selection strategy into the genetic algorithm, so that the whole iterative process is more intelligent. In addition, various repair operators are proposed to repair the infeasible schemes obtained after crossover and mutation.
The genetic algorithm optimization framework is shown in fig. 2. At the beginning of each iteration, the selection operator randomly selects a group of individuals to perform the following operation, and the individuals with higher fitness are considered as the better individuals and have higher probability of being selected. Therefore, the inverse of the objective function is taken as the fitness function: f 1/CT'。
The probability that each individual is selected can be expressed as:
in the formula, piProbability of being selected for individual i; f. ofiThe fitness function value of the individual i. Then, using roulette, NIND' individuals (NIND × GGAP) were selected as initial populations, and prepared for crossover and mutation.
The invention provides three kinds of crossover operators, namely site crossover, line crossover and scheme crossover.
In site crossing, Parent1, Parent2 represent two parents for crossing, as shown in fig. 3. The operator first randomly selects a line from the two individuals, respectively, then cuts each line at a single site and swaps subsequences of the two lines to generate two new lines. The purpose of this operator is to explore better site order, changing routes from a microscopic perspective.
This operator may be ineffective when two lines are split and switched at randomly selected sites, and this crossover operation is considered an effective crossover only when at least one feasible line is generated that is different from the other lines, and therefore the emphasis of operator design is on how quickly to find two reasonable split points to accomplish an effective crossover of lines. In order to solve the problem, the invention provides two single-point crossing methods as the sub-heuristic algorithms of the operator, namely single-point crossing-1 and single-point crossing-2.
The single point cross-1 is described in detail as follows:
the single point cross-2 is described in detail as follows:
the loop cancellation algorithm is used to cancel a local loop present in a line, and is described in detail as follows:
to compare the performance of the two intersection methods, the overlap ratio of the two lines is quantified by the ratio of the number of nodes of intersection to the number of nodes of union, and is expressed by equation (31):
in the formula, cr1,r2 1∈[1,1]Indicating the degree of coincidence of r1 and r 2. c. Cr1,r2 1A larger value of (a) indicates a higher coincidence degree of the two lines.
Curve of a single point crossing-1 with cr1,r2 1Gradually decreases with increasing; while the curve of the single-point cross-2 has a fitting interval, which can be expressed as fs,fe]. In [ fs,fe]Inner, single point cross-2 stable performance, but the curve drops sharply outside the interval, fs,feThe value of (c) may vary with the scale of the problem.
In terms of efficiency, single-point cross-2 operates at a slightly faster speed than single-point cross-1 because single-point cross-1 repeatedly runs the shortest path search requiring a large amount of computation time. The detailed steps of the site intersection operator are as follows:
unlike site crossing, line crossing aims to explore different line combinations from the macro level to create a new net scenario. In operator design, a line intersection separates an individual from one or more pairs of split points, which are all the end points of the line, as shown in fig. 4:
in fig. 4, the line crossing operator generates a new individual by swapping the same number of lines from two individuals, the present invention limits the number of swapped lines to:
in the formula, nsIs the number of lines switched; and | P1|, and | P2| are the number of lines of two individuals.
Notably, after selecting the route r1 from parent generation 1, the route r2 selected from parent generation 2 is not randomly selected. In order to ensure the coverage cell and the number of subway stations connected by the generated new individual to the maximum extent, a route with the maximum correlation coefficient r1 is selected from the parent generation 2 as r2, and the correlation coefficient calculation method comprises the following steps:
in the formula, c2 r1,r2Is the correlation coefficient of r2 and r1, Zr1 *Set of cells covered by the service area of r1, Zr2Set of cells covered by r2 service area, STr1 *Set of subway stations connected in r1, STr2For the set of subway stations connected in r2, | r2| is the number of nodes of r 2.
The line crossing operator is described in detail as follows:
a solution crossover is to split an individual from a single point and swap subsequences between two parents to form a pair of new individuals. The split points may be end points or intermediate points of the line. This operator is easy to implement and it can explore different line combinations and site sequences simultaneously, as shown in fig. 5.
Similar to a site crossing, if the splitting point in the solution crossing is located in the middle of a line, the operator involves the exchange of the site sequence between two lines, which can be done by the mixed use of single-point crossing-1 and single-point crossing-2, and after the new line is generated, the loop elimination algorithm and line length repair 1/2 need to be performed.
It is noted that the two lines used for crossing are not randomly selected, as shown in fig. 5, the selected lines divide the whole net scheme into two line groups, one before the other, and if the selected line position order is biased to one side, the number of lines included in the two line groups is unbalanced, which will make the descendants include too many or too few lines, which is not favorable for generating a competitive scheme. Therefore, the invention makes a weight coefficient for each line in the net scheme based on the gaussian function:
in the formula, wckIs the weight coefficient of line k, where 'k' is the sequential index of the line; the | R | is the number of lines in the line network; ps iskIs the probability that line k is selected.
The invention provides three mutation operators, namely deletion mutation, insertion mutation and shortening mutation. Deleting a variation is to change the scheme by deleting a line, and inserting a variation and shortening a variation can change the scheme by changing the node sequence of a line, as shown in fig. 6.
In fig. 6(a), the deletion variant is to remove one line from the parent to generate the child, and the deleted line is the line with the smallest importance coefficient, which can be quantified according to the number of coverage cells and connection sites:
in the formula, cr 3Is the importance coefficient of the line r; i Zr *L is the number of cells covered in the service area of the line r, | STr *And | is the number of subway stations connected by the line r. The line selection strategy reduces the damage to individual feasibility caused by line deletion, and the purpose of the mutation operator is to stabilize the number of lines in the scheme, because new lines can be added to the scheme during the scheme repair process.
The insertion mutation operator selects a line R (L) from the individualR≤Lmax) And attempts to insert a node (which is not present in the line R and can simultaneously connect adjacent nodes (bus stations or subway stations) in the line R) into the line R, as shown in fig. 6 (b). If all nodes can not be inserted into the line R, replacing one line and repeating the insertion operation until the nodes are successfully inserted or the line is selectedThe number reaches the maximum number of attempts max rs1, wherein | NSfL to satisfy the line length constraint (L)R≤Lmax) The detailed steps are as follows:
the probability of success of inserting a node using the shortest path search method is higher than the probability of success by detecting the connectivity of the node with the node in the line. When searching for the shortest path, some unexpected nodes may be introduced into the line, as shown in fig. 6(b), and node 3 is the unexpected node introduced when node 4 is inserted, and since the line has a local loop, node 4 is deleted later, and node 3 is still retained. This allows this to happen since the purpose of this operation is to accomplish as efficient an insert operation as possible without the need to insert the specified node.
As shown in FIG. 6(c), the shortening mutation is to select a line R (L) from the individualsR≥Lmin) And attempts to replace a particular sequence of nodes with the shortest path. If the entire route is the shortest route from the start point to the end point, the process is repeated for a route replacement until the route replacement is successfully performed or the maximum number of attempts for route selection max _ rs2 is reached. Similarly, letWherein | NSfI satisfies the length constraint (L) in individualsR≥Lmin) The number of lines of (2).
In the scheme after the crossover and mutation operators, infeasible lines still exist, and the infeasibility is mainly expressed in the following aspects:
(1) two or more identical lines exist in the line network scheme;
(2) there are areas not covered by the current line service area;
(3) not all subway stations are connected to the network.
In order to ensure the feasibility of the scheme, the invention provides a repair operator, which comprises the following detailed steps:
the repair algorithm detects the condition (1) first, and eliminates the repeated line (step 2); then detecting the condition (2), inserting some sites which can expand the service range of the net into the net scheme (step 4-7); similarly, the algorithm inserts unconnected subway stations into the current net (steps 8-10). It should be noted that finally, if there are still uncovered cells or unconnected subway stations, new lines will be generated to ensure the feasibility of the solution.
The basic idea of the adaptive operator selection strategy is: the score is given to each previous stage operator, and measures the performance of the operators so that the probability of selecting crossover and mutation operators is adaptively adjusted in the next stage according to their previous performance.
In the design of this strategy, the whole iterative process is divided into several stages, each stage being defined as a certain algebra. Each stage is the basic unit for updating the score of an operator according to the performance of the operator, and each stage can be further divided into multiple iterations of intersection and variation (and intersection probability P)cAnd the mutation probability PmCorrelated), that is to say, at the end of each phase, each newly generated individual is detected according to the different cases shown in table 1 and passes through σ16And scoring each operator.
TABLE 1 fraction adjustment parameters and description
When the operator is given the adjustment parameter, two aspects are mainly considered: one is the objective function value of the generation scheme; the other is whether it will generate a new wire net scheme. This strategy is designed to use as many operators as possible that can improve the solution, promoting solution diversity. If a new solution is generated after the crossover operator and mutation operator are executed, the two operators will get the same score because there is no way to distinguish how effective the two operators are for generating the new solution.
The weight calculation method of the operator i at the next stage comprises the following steps:
in the formula, wijFor the weights of algorithm i in run phase j,wi1=0;λ∈[0,1]a coefficient for adjusting the recent expression change reaction speed of an operator; beta is aijIs the total fraction, gamma, obtained after the operating phase j for the operator iijThe number of times operator i is used at the end of run phase j. At the beginning of the run phase j +1, the beta values of all operators are seti,j+1,γi,j+1Reset to 0. In the operating phase j +1, the probability that the operator is selected is:
in the formula, pi,j+1For the probability that operator i is selected in operating phase j +1,pi,1=1/3。
to avoid premature trapping of genetic algorithms into suboptimal, the present invention proposes an evolutionary stagnation detection strategy, i.e. if the optimal solution does not improve significantly in multiple successive operational phases, the population initialization algorithm is used to update the selected population at the beginning of the next operational phase to ensure that genetic algorithms can search for optimized net solutions in a wider solution space.
Step five: wire mesh optimization scheme
The method comprises the steps of optimizing a practical public traffic network by using a genetic algorithm of a plurality of intersection and mutation operators and a self-adaptive operator selection strategy combination, obtaining an optimized public traffic network scheme, and providing information such as the trend, departure interval and fleet scale of each line in the network scheme.
Drawings
FIG. 1 is a graph of path conversion discrimination.
FIG. 2 is a flow chart of a genetic algorithm.
Fig. 3 is a schematic cross-site diagram.
Fig. 4 is a schematic cross-line diagram.
FIG. 5 is a schematic cross-sectional view of the protocol.
FIG. 6 is a schematic diagram of a mutation operator.
Fig. 7 studies regional road networks and cell maps.
Fig. 8 studies a regional topology network map.
Fig. 9 studies traffic zone division and travel demand of a zone.
Figure 10 current bus net.
Figure 11 optimizes a bus network scheme.
FIG. 12 is a flow chart of an embodiment of the method of the present invention.
Detailed Description
The real case study area is the Tokyo Community in Beijing, which has an area of about 4.8km × 2.2km and is located in the suburban area of Beijing, as shown in FIG. 7. The area is divided into 264 grid cells, the side length of each cell is 200 meters, 106 nodes are totally arranged in a road network, the nodes comprise 1 station (node 0), 55 bus stations (nodes 1-55) and 6 subway stations (nodes 56-61), and the topological network is shown in fig. 8.
In order to measure and calculate the travel demand between each grid cell and each subway station, the public transportation IC card data is used for identifying a 'public transportation-subway' combined travel chain of each bus station in an early peak period in a research area, then the demand at the station is distributed to the cells covered by a service range according to a certain proportion, the proportion can be determined by the total travel demand in the current traffic cell, and the calculation formulas (39) - (40) are given out for the demand of each grid cell.
In the formula, Do,dThe public transportation trip demand (bidirectional demand) from the grid cell o to the subway station d; SDs,dThe current bus demand from the bus stop s to the subway stop d is obtained through the IC card data; soA station set within the maximum walking distance range of the grid cell o; TDoThe total travel demands of various traffic modes in the grid cell o are obtained; epsilons,zAs a coefficient, if the grid cell z is within the walking distance of the station s, then es,z1, otherwise, ∈s,z=0;TTDkThe total travel demand of various traffic modes in the traffic cell k is obtained through the mobile phone signaling data; etao,kThe ratio of the overlapping area of the grid cell o and the traffic cell k to the area of the traffic cell k.
Traffic cell division in the research area is shown in fig. 9, the number in each cell is the number of the traffic cell, and the demand of all traffic modes in each traffic cell is obtained through mobile phone signaling data.
Currently, there are seven lines serving the community during peak hours, each of which is connected to at least one subway station, as shown in fig. 10. Lines 101 and 102 are private lines serving community passengers, while the other lines are public lines that only partially pass through the area under study.
In calculating the fleet size and current operating costs, a factor defined by L '/L is considered, L representing the total link length and L' representing the length of the portion of the link within the study area, and the information for each link is shown in table 2.
Table 2 present information per line
aConverted fleet size, F ═ F × (L'/L);
bthe operating cost of the conversion, CS ═ CS × (L'/L).
The method comprises the steps that various costs and punishments are measured and calculated according to a given passenger flow distribution method, the converted fleet scale and the operation cost can be used for evaluating the current network scheme, and in the process of solving the optimization scheme, the converted fleet scale is used as the upper limit of the fleet scale, namely TF-57 veh. Other parameter settings in the model: vA=5.8km/h;VB=25km/h;ts1/90 h; p is 60 persons; l ismin=5km;Lmax14 km; $ 100/hour; mu.sA$ 10/hour; mu.sW$ 10/hour; mu.sI$ 5/hour; mu.sP1500 yuan/km; mu.sP2=1;μP3$ 10/person; z is a radical ofw=0.5,g(x)=0.01×(<0.001). The parameters of the genetic algorithm are set as: NIND-20; MAXGEN 1000; pc=0.9;Pm0.1; the length of the operation stage is 20 generations; sigma1=8,σ2=6,σ3=3,σ4=2,σ5=1,σ61 is ═ 1; λ is 0.5. For simplicity, the present disclosure treats a demand that requires more than two transfers as an unserviceable demand.
The data in table 3 are the costs and penalties of the current net and the costs and penalties of the optimized net.
TABLE 3 optimization scheme information
*The experiment was repeated 5 times to obtain results;
arun the number of detected protocols 5 times;
as shown in table 3, the optimized net scheme is significantly better than the existing net scheme, and the total cost of the searched net scheme is reduced by 57.09%. In practical application, by redesigning the community public transportation network under the scale of the existing fleet, the travel cost of passengers can be reduced by 62.8%, and the operation cost can be reduced by 18.3%. The optimization scheme is shown in table 4, and the route trend is shown in fig. 11.
TABLE 4 optimization scheme

Claims (6)

1. A community public transport network and departure frequency synchronous optimization method for a connection subway considering area full coverage is characterized by comprising the following steps:
step one, constructing a bus network optimization model
Providing a constraint condition based on reality, a line length constraint, a net total service range constraint, a vehicle capacity and a fleet scale constraint, taking the total cost, namely the minimum value of user cost and operation cost, as an optimization target, giving consideration to community microcirculation public traffic net layout and line departure frequency setting, and constructing a mixed integer optimization model;
step two, initializing the public traffic network
For the constraint of the line length, a length repair algorithm is provided, so that the length of each line is in a reasonable range; when the area full coverage constraint is realized, designing a corresponding algorithm, namely, repeatedly generating a new line to serve the current farthest uncovered area until all areas are covered; finally, detecting the number of the subway stations to ensure that all the subway stations are connected to the network;
step three, evaluating the fitness of the wire mesh
Fitness evaluation is used for evaluating the overall performance of a wire mesh scheme and comprises two parts: passenger flow distribution and fleet scale setting; the passenger flow distribution aims at distributing travel demands in the grid cells to bus stops or route sets of a bus network, so that the getting-on and getting-off passenger flow of each stop and the section passenger flow of each line are obtained; the motorcade size determination method is that a given number of vehicles are accurately distributed to each line under the condition of considering the total motorcade size and the vehicle capacity; finally, the fitness of the net scheme, namely the quality degree, is reflected by calculating a model objective function value;
step four, optimizing algorithm design
Optimizing the network scheme by using a genetic algorithm, and designing a plurality of crossover operators, mutation operators and repair operators to ensure that a new feasible scheme is obtained based on the original scheme;
step five, network optimization scheme
Based on the characteristics of different algorithm combinations, the final bus network optimization scheme is obtained by using the algorithm combination with the best performance.
2. The method for optimizing the community public transportation network and departure frequency of the connected subway in consideration of the full coverage of the area as claimed in claim 1, wherein said step one,
considering various constraint conditions based on reality, providing corresponding penalty functions for the line length, the vehicle load and the service level, and constructing a mixed integer optimization model by taking the minimization of the total system cost as an optimization target;
the proposed constraints are as follows:
constraint 1: cell full coverage constraint; in order to improve the accessibility of passengers in the community during traveling, the service range of a public traffic network needs to cover all grid cells;
constraint 2: the full coverage of the subway station is restricted; in the network, each subway station appears at least once to ensure the connection of two public transportation modes;
constraint 3: total fleet scale constraints; the fleet size of the line network cannot exceed the upper limit of the total fleet size;
constraint 4: a vehicle capacity limit; on each line, the section passenger flow of each journey cannot exceed the maximum passenger capacity of the operating vehicle;
constraint 5: line length constraints; the length of each line should be within a reasonable range;
constraint 6: line trend constraint; repeated stations cannot exist in each line, so that a turning back or a local loop is avoided;
the objective function contains user cost and operation cost, and is expressed as:
CT=CU+CS (1)
in the formula, CTIs the total cost; cUThe cost for the user; cSIs the operating cost;
the user cost is divided into three parts: cost of passenger walking CAPassenger cost in-vehicle CIPassenger waiting cost CW:
CU=CA+CI+CW (2)
In the formula, CA、CIBoth related to the net layout and the passenger selected route are represented as:
in the formula, Z is a square cell set; pathodIs a selectable path set between o and d; dod kThe walking distance from the cell o to the corresponding getting-on/off station on the kth path between o and d; l isod kIs the on-vehicle distance in the kth path between o and d; t is tsDelaying the time for the station; nsod kThe number of sites of the kth path between o and d; q. q.sodThe number of passengers between o and d (or d and o); p is a radical ofod kProbability of selecting the kth path for the passenger; vA、VBThe walking speed and the vehicle speed of the passenger; mu.sA、μITo the walk time value and the time-on-vehicle value;
in formula (3), Dod kOn the k-th path between o and dThe walking distance from the cell o to a bus stop within the maximum walking distance range tolerable by the residents of the cell; if a plurality of stations exist on the path within the range, matching the residents according to the minimum sum of the walking time and the time in the vehicle; for a cell close to a subway station, D may be within the range of cell o, meaning that the resident will walk directly to the subway station instead of riding a regular bus, matching the cell directly to the subway station, Dod kRepresents the walking distance from o to d;
the average waiting cost of passengers on one path is related to the departure interval of the related line; cWExpressed as:
in the formula, R is a line set of a line network scheme; nsrIn the one-way running of the vehicles, the number of bus stations/subway stations needing stopping of regular buses is except for bus stations; HBrThe departure interval of the line r; nbs rThe number of passengers boarding route r at station s; z is a radical ofwA ratio of arrival time/departure interval desired for the passenger; mu.sWThe value of waiting time is obtained;
the operation cost is related to the net layout and the departure interval of each line, the operation cost is equal to the total fleet size multiplied by the bus operation cost, and the fleet size of each line is further expressed as formula (7);
in the formula, B is the operation cost of the vehicle per hour; frIs the fleet size of line r; l isrIs the unidirectional length of the line r;
0-1 decision variable if node j is subsequent to node i in line rWhen node, xij r=1,j belongs to SST, R belongs to R; otherwise, xij r=0;Lr、nsrExpressed as:
in the formula Iij rIs the distance on line r from i to j;
line length penalty CP1The calculation method is as follows:
in the formula, muP1Is a coefficient;
penalty function for vehicle capacity, g (x):
g(x)=ω(αx-1) (11)
in the formula, x is the cross section passenger flow volume; omega and alpha are coefficients;
overload penalty C for net schemeP2The calculation method is as follows:
in the formula, Ldij rDetermining the section passenger flow quantity between the nodes i and j in the line r per hour according to the passenger flow distribution result; mu.sP2Is a coefficient;
unserviceable demand penalty CP3The calculation method is as follows:
where if the current net serves the travel demand, δod1, otherwise, δod=0;μP3Is a coefficient;
c is to beP1、CP2、CP3Added to the initial objective function, which is expressed as:
min CT'=CA+CI+CW+CS+CP1+CP2+CP3 (14)。
3. the method for optimizing the community public transportation network and departure frequency of the connected subway in consideration of the full coverage of the area as claimed in claim 1, wherein in the second step,
generating a net scheme set meeting all constraint conditions in the model by utilizing a net initialization algorithm, wherein the net scheme set also comprises a sub-algorithm, namely a line length repairing algorithm;
all the randomly generated wire net schemes can meet constraint conditions;
for the cell Z belonging to Z, the bus station S belonging to S and wdzsDenote the walking distance from z to s, let wdmin z=min{wdzs| S ∈ S }, Z0 denotes an uncovered cell, let wd bemax=max{wdmin z-Z ∈ Z0}, ST0 denotes the set of subway stations not connected to the wire network, NS denotes the currently generated wire network solution, WD denotes the maximum walking distance tolerable for the passenger, and the matrix RN denotes the connectivity of every two nodes of the road topology based network; the steps of randomly generating a feasible net scheme are as follows:
for line r, the length repair algorithm uses two ways to increase the length of the line: 1) to pairIn Lr<LminThe length of the line meets the requirement by extending the line; 2) for Lr≥LminAccording to a probability p associated with the current line lengthinsInserting some nodes into the line r:
the detailed steps are as follows:
for line r, line length repair-2 shortens the line in two ways: 1) for Lr>LmaxThe shortest path between the nodes is used for replacing the sub-line, so that the line length is shortened; 2) for Lr≤LmaxAccording to a probability p associated with the current line lengthremDeleting nodes from r:
the detailed steps are described as follows:
two line length repair heuristic algorithms ensure that each line is within a proper length range to the maximum extent, but the connectivity of a receiving network and the Lmin,LmaxAnd (4) influence of the value, if the line cannot meet the length constraint, calculating the line length penalty, and adding the line length penalty into the objective function.
4. The method for optimizing the community public transportation network and departure frequency of the connected subway in consideration of the full coverage of the area as claimed in claim 1, wherein said step three,
determining alternative paths for passenger travel between each pair of ODs based on a public traffic network, constructing a Logistic model taking walking time and vehicle time as effect functions, calculating the proportion of passenger flow on each path between the ODs, determining route departure frequency, providing a fleet scale adjustment algorithm to meet the limitation of the maximum fleet, and finally, calculating a target function value of a scheme;
Pathodrepresenting the path set with the least transfer between the o and the d, and determining the boarding station and the transfer station of the passenger through the matching relation between the path and the public traffic network, thereby calculating the cross-section passenger flow of each cross section on each line; finding out all paths between each pair of ODs by using the method in graph theory, and judging the transfer times of each Path so as to determine Pathod
In determining PathodIn the process of (1), firstly, whether the cell o is in the service range of the subway station d is detected, if yes, the Path between o and d is the shortest walking Path, and the shortest walking Path is recorded in the PathodSuch a demand does not require a regular bus ride — scenario 0; otherwise, whether a riding path without transfer, namely scene 1, a riding path with one transfer, namely scene 2, and a riding path with two transfers, namely scene 3 exist between each pair of ODs is respectively detected, if the path without more than two transfers, namely scene 4, does not exist, the demand between the OD pairs is regarded as the non-service demand, and the corresponding C is calculatedP3And adding it to the objective function; moreover, there is no case where two scenarios occur simultaneously, and the occurrence of each scenario will terminate the algorithm;
in scenarios 1-3, PathodThe generation mechanism of (2) is similar, and the most obvious difference among the three scenes is that the cycle number in the algorithm increases exponentially; in three scenarios, in order to remove too many detours, save the calculation time, a screening strategy is formulated: will cost the effect more than PathodThe minimum value of all paths in the Path is 1.5 times from PathodDeleting; scene 4 is used as a supplement of other scenes, the paths of the scenes exceed twice transfer, the travel demand of the OD pair is identified as the unserviceable demand, and corresponding punishment is calculated;
In scenarios 2 and 3, the initial travel demand needs to be allocated to the bidirectional routes, and since the routes are directional, while recording the transfer station of each route, the directions of the sub-routes distributed on different routes need to be considered, so as to calculate the number of people getting on or off at each station during the round trip of the vehicle;
constructing a Logistic model taking walking time and vehicle time as effect functions, and calculating to obtain the passenger flow proportion on each path between the od:
in the formula, Vod kThe effect function for the kth path is expressed as:
in the formula, a1,a2Is a coefficient;
suppose that ns exists in a single pass of line rrThe round trip of a bus station/subway station (except for a station) can be divided into (2 ns)r-1) section forNbs rIndicates the number of passengers boarding the station s (see equation (5)),the number of passengers getting off the vehicle at the station t can be expressed as:
in the formula (I), the compound is shown in the specification,all are variables from 0 to 1; if a passenger on the k-th route gets on board at station s on route r, thenIf the k-th path passenger transfers at station s on route r, thenIf the passenger on the k-th route stops at a stop t on route r, thenIf the k-th path passenger transfers at a stop t on line r, then
Passenger flow volume Ld of upper section of line rij rThe calculation method comprises the following steps:
a departure frequency setting method, which aims to determine the optimal departure frequency for each line, comprises two stages: a preliminary calculation stage and a motorcade adjustment stage;
in the preliminary calculation stage, the departure interval Hb of the line rr *The calculation method comprises the following steps:
Hbrthe calculation method comprises the following steps:
HMrthe calculation method comprises the following steps:
divide the round trip time of the line by Hbr *The value of (b) is rounded up to be the number of vehicles needed by the line r, namely the fleet scale of the line, and the calculation method comprises the following steps:
frequency f of departurer *The calculation method comprises the following steps:
for all lines s ∈ R, Fr *The sum of (2) may exceed the fleet scale upper limit TF, when this occurs, the fleet scale of the line needs to be adjusted, and the total fleet scale meets the constraint condition by reducing the number of vehicles of some lines, i.e. the fleet adjustment stage;
one net scheme includes a line set R and a fleet size set F ═ F1,…Fr,...,F|R|},ΔCT rRepresents the total cost change of the route r after removal of the vehicle, which can be expressed as Δ CT r=ΔCW r+ΔCS r+ΔCP2 r(ΔCW r,ΔCP2 r>0,ΔCS r<0) The algorithm is described in detail as follows:
5. the method for optimizing the community public transportation network and departure frequency of the connected subway in consideration of the full coverage of the area as claimed in claim 1, wherein said step four,
optimizing the line network by using a genetic algorithm, providing a plurality of crossover and mutation operators in the genetic algorithm, and embedding a self-adaptive operator selection strategy into the genetic algorithm, so that the whole iterative process is more intelligent; in addition, a plurality of repair operators are provided to repair the infeasible schemes obtained after intersection and mutation;
when each iteration starts, the genetic algorithm randomly selects a group of individuals to perform the following operations, the individuals with higher fitness are considered as the better individuals, and the selected probability is higher; taking the inverse of the objective function as the fitness function: f 1/CT';
The probability that each individual is selected can be expressed as:
in the formula, piProbability of being selected for individual i; f. ofiA fitness function value of the individual i; then, selecting NIND' individuals (NIND multiplied by GGAP) individuals to form an initial population by using a roulette method, and preparing for crossing and mutation;
three kinds of crossover operators are provided, namely site crossover, line crossover and scheme crossover;
in site crossing, a line is randomly selected from two individuals respectively, then each line is cut at a single site, and subsequences of the two lines are exchanged to generate two new lines; the operator aims at exploring a better site sequence and changing a line from a microscopic angle;
when two lines are split and switched at randomly selected sites, the operator may be invalid, and this crossover operation is considered to be an effective crossover only when at least one feasible line different from the other lines is generated, so the operator design is focused on how to quickly find two reasonable split points to complete an effective crossover of lines; in order to solve the problem, two single-point crossing methods are provided as the sub-heuristic algorithms of the operator, namely single-point crossing-1 and single-point crossing-2;
single point cross-1, as follows:
single point cross-2, as follows:
the loop cancellation algorithm is used to cancel the local loop present in the line, as follows:
the site intersection operator comprises two single-point intersection operators, as follows:
unlike site crossing, line crossing aims at exploring different line combinations from the macro level to generate a new line network scheme; in operator design, line crossing separates an individual from one or more pairs of split points, which are all end points of a line, and a new individual is generated by exchanging the same number of lines from two individuals; the number of lines switched is limited to:
in the formula, nsIs the number of lines switched; p1 and P2 are two individualsThe number of lines;
it should be noted that, when selecting a route, after selecting the route r1 from the parent1, the route r2 selected from the parent2 is not randomly selected, but a route with the largest correlation coefficient with r1 is selected from the parent2 as r2, and the correlation coefficient can be expressed as:
in the formula, c2 r1,r2Is the correlation coefficient of r2 and r1, Zr1 *Set of cells covered by the service area of r1, Zr2Set of cells covered by r2 service area, STr1 *Set of subway stations connected in r1, STr2The number of nodes is r2, | r2| is a subway station set connected in r 2;
the operation of the line crossing operator is described as follows:
a solution crossover is to split an individual from a single point and exchange subsequences between two parents to form a pair of new individuals, the split points can be end points or intermediate points of a line, which can simultaneously explore different line combinations and site sequences; similar to the site crossing, if the splitting point in the solution crossing is located in the middle of the line, the operator includes the exchange of site sequences between the two lines, the operation can be completed by the mixed use of the single-point crossing-1 and the single-point crossing-2, and after a new line is generated, a loop elimination algorithm and line length repair 1/2 need to be executed;
the two lines in the solution intersection are not randomly selected, because the selected lines divide the whole net solution into two line groups, one before the other, and if the position sequence of the selected lines is biased to one side, the number of lines included in the two line groups is not balanced, which will make the filial generation include too many or too few lines, which is not favorable for generating a competitive solution; therefore, a weighting factor is formulated for each line in the net scheme based on the gaussian function:
in the formula, wckIs the weight coefficient of line k, where 'k' is the sequential index of the line; the | R | is the number of lines in the line network; ps iskProbability of being selected for line k;
three mutation operators are provided, namely deletion mutation, insertion mutation and shortening mutation; deletion of a variation is to change the scheme by deleting a line, and insertion of a variation and shortening of a variation can change the scheme by changing the node sequence of a line;
deletion variant is to remove a line from a parent to generate a child, the deleted line being the line with the smallest importance coefficient, which can be quantified according to the number of coverage cells and connection sites:
in the formula, cr 3Is the importance coefficient of the line r; i Zr *L is the number of cells covered in the service area of the line r, | STr *L is the number of subway stations connected by the line r;
the insertion mutation operator selects a line R (L) from the individualR≤Lmax) And try to insert a node (not present in the line R and capable of connecting adjacent nodes (bus stations or subway stations) in the line R at the same time) into the line R; if all nodes cannot be inserted into the line R, the insertion operation is repeated by replacing one line until the node is successfully inserted or the number of line selections reaches the maximum number of attempts max _ rs1,wherein | NSfL to satisfy the line length constraint (L)R≤Lmax) The detailed steps are as follows:
the shortened variants are selected from the individuals for the line R (L)R≥Lmin) And attempting to replace the particular sequence of nodes with the shortest path; if the whole line is the shortest path from the starting point to the end point, replacing one line and repeating the process until the path replacement is successfully realized or the maximum attempt time max _ rs2 of line selection is reached; similarly, letWherein | NSfI satisfies the length constraint (L) in individualsR≥Lmin) The number of lines of (a);
the crossover and mutation operations result in unfeasible lines, which are expressed in the following aspects:
(1) two or more identical lines exist in the line network scheme;
(2) there are areas not covered by the current line service area;
(3) not all subway stations are connected into a wire network;
in order to ensure the feasibility of the scheme, a repair operator is provided, and the detailed steps are as follows:
the repair algorithm detects the condition (1) first, and eliminates the repeated line (step 2); then detecting the condition (2), inserting some sites which can expand the service range of the net into the net scheme (step 4-7); similarly, the algorithm inserts unconnected subway stations into the current net (steps 8-10); if there are still uncovered cells or unconnected subway stations, a new line is generated and added to the scheme (step 11);
the basic idea of the adaptive operator selection strategy is to score the scores of operators in each previous stage, the scores measuring the performance of the operators, so that the probability of selecting each pair of crossover and mutation operators is adaptively adjusted in the next stage iteration according to their previous performance;
in the design of this strategy, the whole iterative process is divided into several stages, each stage being defined as a certain algebra; each stage is a basic unit for updating the score of an operator according to the performance of the operator, and each stage can be further divided into multiple iterations of intersection and variation (and intersection probability P)cAnd the mutation probability PmCorrelated), that is to say, at the end of each phase, each newly generated individual is detected according to the different cases shown in table 1 and passes through σ16Scoring each operator;
TABLE 1 fraction adjustment parameters and description
The operator is endowed with an adjusting parameter, and two aspects are considered: one is the objective function value of the generation scheme; the other is whether it will generate a new wire net scheme;
the operator i weight of the next stage is updated as:
in the formula, wijFor the weights of algorithm i in run phase j,wi1=0;λ∈[0,1]a coefficient for adjusting the recent expression change reaction speed of an operator; beta is aijIs the total fraction, gamma, obtained after the operating phase j for the operator iijThe number of times operator i is used at the end of run phase j; at the beginning of the run phase j +1, the beta values of all operators are seti,j+1,γi,j+1Reset to 0; in the operating phase j +1, the operatorThe probability of being selected is:
in the formula, pi,j+1For the probability that operator i is selected in operating phase j +1,pi,1=1/3。
6. the method for optimizing the community public transportation network and departure frequency of the connected subway in consideration of the full coverage of the area as claimed in claim 1, wherein said step five,
a practical public transport network is optimized by using a genetic algorithm of a plurality of intersection and mutation operators and a self-adaptive operator selection strategy combination, an optimized public transport network scheme is obtained, and the trend, departure interval and fleet scale information of each line in the network scheme are given.
CN201910880061.4A 2019-09-18 2019-09-18 Method for synchronously optimizing community public transportation network and departure frequency of connection subways in consideration of full coverage of area Active CN110598942B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910880061.4A CN110598942B (en) 2019-09-18 2019-09-18 Method for synchronously optimizing community public transportation network and departure frequency of connection subways in consideration of full coverage of area

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910880061.4A CN110598942B (en) 2019-09-18 2019-09-18 Method for synchronously optimizing community public transportation network and departure frequency of connection subways in consideration of full coverage of area

Publications (2)

Publication Number Publication Date
CN110598942A true CN110598942A (en) 2019-12-20
CN110598942B CN110598942B (en) 2023-10-20

Family

ID=68860622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910880061.4A Active CN110598942B (en) 2019-09-18 2019-09-18 Method for synchronously optimizing community public transportation network and departure frequency of connection subways in consideration of full coverage of area

Country Status (1)

Country Link
CN (1) CN110598942B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111738895A (en) * 2020-07-20 2020-10-02 平安国际智慧城市科技股份有限公司 Multi-station passenger transport line generation method, device, medium and electronic equipment
CN112149878A (en) * 2020-08-31 2020-12-29 北京工业大学 Pure electric bus running plan and charging plan synchronous optimization method considering incomplete charging
CN112434844A (en) * 2020-11-10 2021-03-02 郑州天迈科技股份有限公司 New development and extension method for sorting net based on convex hull calculation and genetic algorithm
CN113139030A (en) * 2021-04-27 2021-07-20 北京交通发展研究院 Bus microcirculation line determination method and device, computer equipment and storage medium
CN113514071A (en) * 2021-06-29 2021-10-19 阿波罗智联(北京)科技有限公司 Method and device for determining public transport route
CN113570185A (en) * 2021-06-08 2021-10-29 深圳市综合交通设计研究院有限公司 Bus station supply and demand index system based on GIS and publishing method thereof
CN113962615A (en) * 2021-12-22 2022-01-21 华睿交通科技有限公司 Traffic resource management and control system, method and storage medium
CN115186049A (en) * 2022-09-06 2022-10-14 深圳市城市交通规划设计研究中心股份有限公司 Intelligent bus alternative station site selection method, electronic equipment and storage medium
CN115497326A (en) * 2022-11-21 2022-12-20 南京邮电大学 Multi-mode bus travel path planning method
CN115577818A (en) * 2022-12-01 2023-01-06 武汉好人科技股份有限公司 Passenger demand response type carpooling scheduling method and system for intelligent bus
WO2024045727A1 (en) * 2022-08-30 2024-03-07 比亚迪股份有限公司 Line network dispatching method and apparatus, medium, and electronic device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140149320A1 (en) * 2012-11-29 2014-05-29 International Business Machines Corporation Consistent price optimization in transportation networks
CN105205557A (en) * 2015-09-17 2015-12-30 华南理工大学 Design method for conventional urban public transit network
CN106097226A (en) * 2016-06-20 2016-11-09 华南理工大学 City Routine Transit Network Design method based on Hierarchical Programming
CN106779163A (en) * 2016-11-18 2017-05-31 华南理工大学 A kind of customization transit network planning method based on intelligent search
CN110119835A (en) * 2019-03-26 2019-08-13 杭州电子科技大学 A kind of public transport dynamic based on interval computation is dispatched a car method for optimizing scheduling
CN110245779A (en) * 2019-05-10 2019-09-17 杭州电子科技大学 A kind of public transport dynamic based on genetic algorithm is dispatched a car method for optimizing scheduling

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140149320A1 (en) * 2012-11-29 2014-05-29 International Business Machines Corporation Consistent price optimization in transportation networks
CN105205557A (en) * 2015-09-17 2015-12-30 华南理工大学 Design method for conventional urban public transit network
CN106097226A (en) * 2016-06-20 2016-11-09 华南理工大学 City Routine Transit Network Design method based on Hierarchical Programming
CN106779163A (en) * 2016-11-18 2017-05-31 华南理工大学 A kind of customization transit network planning method based on intelligent search
CN110119835A (en) * 2019-03-26 2019-08-13 杭州电子科技大学 A kind of public transport dynamic based on interval computation is dispatched a car method for optimizing scheduling
CN110245779A (en) * 2019-05-10 2019-09-17 杭州电子科技大学 A kind of public transport dynamic based on genetic algorithm is dispatched a car method for optimizing scheduling

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
罗孝羚等: "中小城市公交线网及发车频率同步优化", 《工业工程》 *
邢雪: "基于粒子群算法的城市接驳公交网络优化调度方法", 《北京工业大学学报》 *
陈巍等: "公交线网及发车频率同步优化研究", 《交通运输工程与信息学报》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111738895A (en) * 2020-07-20 2020-10-02 平安国际智慧城市科技股份有限公司 Multi-station passenger transport line generation method, device, medium and electronic equipment
CN112149878A (en) * 2020-08-31 2020-12-29 北京工业大学 Pure electric bus running plan and charging plan synchronous optimization method considering incomplete charging
CN112149878B (en) * 2020-08-31 2024-06-04 北京工业大学 Synchronous optimization method for pure electric bus driving plan and charging plan considering incomplete charging
CN112434844A (en) * 2020-11-10 2021-03-02 郑州天迈科技股份有限公司 New development and extension method for sorting net based on convex hull calculation and genetic algorithm
CN112434844B (en) * 2020-11-10 2024-01-26 郑州天迈科技股份有限公司 New opening and extension method of sorting wire net based on convex hull calculation and genetic algorithm
CN113139030A (en) * 2021-04-27 2021-07-20 北京交通发展研究院 Bus microcirculation line determination method and device, computer equipment and storage medium
CN113570185A (en) * 2021-06-08 2021-10-29 深圳市综合交通设计研究院有限公司 Bus station supply and demand index system based on GIS and publishing method thereof
CN113570185B (en) * 2021-06-08 2024-05-24 深圳市综合交通设计研究院有限公司 GIS-based bus station supply and demand index system and release method thereof
CN113514071A (en) * 2021-06-29 2021-10-19 阿波罗智联(北京)科技有限公司 Method and device for determining public transport route
CN113514071B (en) * 2021-06-29 2024-04-16 阿波罗智联(北京)科技有限公司 Public traffic route determining method and device
US11867518B2 (en) 2021-06-29 2024-01-09 Apollo Intelligent Connectivity (Beijing) Technology Co., Ltd. Method and apparatus for determining public transport route
CN113962615A (en) * 2021-12-22 2022-01-21 华睿交通科技有限公司 Traffic resource management and control system, method and storage medium
WO2024045727A1 (en) * 2022-08-30 2024-03-07 比亚迪股份有限公司 Line network dispatching method and apparatus, medium, and electronic device
CN115186049A (en) * 2022-09-06 2022-10-14 深圳市城市交通规划设计研究中心股份有限公司 Intelligent bus alternative station site selection method, electronic equipment and storage medium
CN115497326A (en) * 2022-11-21 2022-12-20 南京邮电大学 Multi-mode bus travel path planning method
CN115577818A (en) * 2022-12-01 2023-01-06 武汉好人科技股份有限公司 Passenger demand response type carpooling scheduling method and system for intelligent bus

Also Published As

Publication number Publication date
CN110598942B (en) 2023-10-20

Similar Documents

Publication Publication Date Title
CN110598942B (en) Method for synchronously optimizing community public transportation network and departure frequency of connection subways in consideration of full coverage of area
CN110648022B (en) Community public transportation network and departure frequency synchronous optimization method of connection subways considering full coverage of stations
CN110580404B (en) Network operation energy determining method based on urban multi-mode traffic network
CN108053062B (en) Customized bus route generation method based on multi-source data
CN108470444B (en) A kind of city area-traffic big data analysis System and method for based on genetic algorithm optimization
Ciaffi et al. Feeder bus network design problem: A new metaheuristic procedure and real size applications
CN112146673B (en) Expressway multipoint collaborative rescue path planning method based on improved ant colony algorithm
CN109034465A (en) Consider the charging station bi-level optimization method that charging station addressing is coupled with trip route
CN105513400A (en) Method for dynamically planning travel route
CN106779190A (en) A kind of urban track traffic passenger trip route suggesting method and system
JP2001504965A (en) Using multiple cost levels for route discovery computation
CN111754039A (en) Method for comprehensive integrated optimization design of pure electric bus network
CN111879329B (en) Customized public transport passable shortest path calculation method based on A-x algorithm
CN107180274A (en) A kind of charging electric vehicle facilities planning typical scene is chosen and optimization method
CN110599760A (en) Travel behavior simulation method under multi-mode traffic network
CN114037175B (en) Large-scale public transportation network hierarchical optimization method based on multi-scale clustering
Badia Comparison of bus network structures in face of urban dispersion for a ring-radial city
CN102567389A (en) Combined traffic network forming method and equipment as well as path searching method and equipment
Xiong et al. Route network design of community shuttle for metro stations through genetic algorithm optimization
CN115713207A (en) Hybrid bus service combination optimization method
CN115048576A (en) Flexible recommendation method for airport passenger group travel mode
CN115423168A (en) Customized public transport network optimization method and system considering service stability
CN111311002A (en) Bus trip planning method considering active transfer of passengers in transit
CN110162716A (en) A kind of influence power community search method and system based on community&#39;s retrieval
Zhang et al. Timetable synchronization of the last several trains at night in an urban rail transit network

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant