CN113177320B - Simulation optimization method, system, equipment and medium for road dynamic congestion charging - Google Patents
Simulation optimization method, system, equipment and medium for road dynamic congestion charging Download PDFInfo
- Publication number
- CN113177320B CN113177320B CN202110494686.4A CN202110494686A CN113177320B CN 113177320 B CN113177320 B CN 113177320B CN 202110494686 A CN202110494686 A CN 202110494686A CN 113177320 B CN113177320 B CN 113177320B
- Authority
- CN
- China
- Prior art keywords
- simulation
- traffic demand
- traffic
- charging
- charging scheme
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 144
- 238000005457 optimization Methods 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000011156 evaluation Methods 0.000 claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 230000006870 function Effects 0.000 claims description 40
- 238000009826 distribution Methods 0.000 claims description 28
- 238000007726 management method Methods 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000013468 resource allocation Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 230000006872 improvement Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000005315 distribution function Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- HOWHQWFXSLOJEF-MGZLOUMQSA-N systemin Chemical compound NCCCC[C@H](N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(O)=O)C(=O)OC(=O)[C@@H]1CCCN1C(=O)[C@H]1N(C(=O)[C@H](CC(O)=O)NC(=O)[C@H](CCCN=C(N)N)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H]2N(CCC2)C(=O)[C@H]2N(CCC2)C(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)N)C(C)C)CCC1 HOWHQWFXSLOJEF-MGZLOUMQSA-N 0.000 description 1
- 108010050014 systemin Proteins 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
- Evolutionary Computation (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a distributed robust simulation optimization method, a system, equipment and a medium for road dynamic congestion charging, wherein the method comprises the following steps: step 1, for a road network to be subjected to dynamic congestion charging scheme optimization, creating a simulation road network in macroscopic traffic simulation software; step 2, constructing a dynamic congestion charging optimization model of the road by taking the total travel time as a target evaluation index; and 3, initializing the charging scheme set and the traffic demand parameter set, performing simulation evaluation on the created simulation road network by using macroscopic traffic simulation software, and solving the optimization model established in the step 2 by adopting a distributed robust simulation optimization algorithm to obtain a congestion charging scheme with optimal robustness. The method can solve the problems of simulation evaluation and optimization of the road dynamic congestion charging under the condition of traffic demand uncertainty.
Description
Technical Field
The invention relates to the field of road dynamic congestion charging, in particular to a distributed robust simulation optimization method, a system, equipment and a medium for road dynamic congestion charging.
Background
The traffic congestion charging, namely, the management and control measures for applying extra travel fees to the traffic travel in a specific time period, induces the travelers to make travel time and path selection again by increasing the cost of using traffic facilities in a charging area, realizes the transfer of part of the travelers from a congested area to a non-congested area, or selects to travel in a non-charging time period, thereby reducing the traffic demand in the charging area to a certain extent, relieving the contradiction between traffic supply and traffic demand, relieving traffic congestion, enabling road network resources to be more fully utilized, and having positive significance to the healthy development of cities.
The traffic simulation technology is to construct a model to simulate the existing or planning traffic system, so that the user can optimize the researched traffic network system or predict the future traffic state conveniently. However, the traffic simulation technology also has the problems of black box property, long time consumption and the like. Aiming at the problems, a simulation optimization algorithm is generated at the same time, and the calculation speed is greatly improved. The mainstream simulation optimization methods mainly include Gradient-based search methods (Gradient based search methods), Stochastic approximation methods (Stochastic approximation methods), Heuristic algorithms (Heuristic search methods), and Response surface methods (Response surface methods). The response surface method, namely the proxy model method, has the advantages of saving simulation cost and fitting complex and unknown system input and output relations, and the advantages are more prominent under the condition that simulation experiments are time-consuming.
In complex real-world engineering problems, even if the uncertainty effects of the data are ignored, the optimal solution, which is found according to the assumed conditions, is likely to be unstable or even infeasible. Therefore, compared with the global optimal solution, the solution with less sensitivity to disturbance, namely the distributed robust solution, in the practical problem has more realistic significance.
Disclosure of Invention
The invention provides a distributed robust simulation optimization method, a system, equipment and a medium for road dynamic congestion charging, which can solve the problems of road dynamic congestion charging simulation and evaluation under the condition of traffic demand uncertainty.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
a distributed robust simulation optimization method for road dynamic congestion charging comprises the following steps:
step 1, for a road network to be subjected to dynamic congestion charging scheme optimization, creating a simulation road network in macroscopic traffic simulation software;
step 2, taking the total travel time as a target evaluation index, and constructing the following road dynamic congestion charging optimization model:
wherein, S is a congestion charging scheme set, and x is any congestion charging scheme in S; j is the simulation times index, NxTotal number of simulations for x; t is a time segment index, and T is the total number of time segments; o is an initial node set, D is an end node set, OD is an origin and belongs to an OD pair set of which the O end belongs to D, and OD is any OD pair in the OD pair set; a is a road section set in the simulation road network, and a is any road section in the road section set A;obey normal distribution for the traffic demand of od in the time segment tAndmean and variance of the traffic demand of od at time segment t, respectively;a traffic demand matrix for a time segment t, consisting ofForming;the random traffic demand matrix obeys normal distribution in the jth simulation; theta is traffic demandThe set of fuzzy normal distributions that is satisfied,e is a parameter (traffic demand parameter for short) for controlling the traffic demand distribution variance, which represents the uncertainty of the traffic demand, e belongs to [ e [l,eu]Wherein e isl、euUpper and lower limits of e, respectively;representing macroscopic traffic simulation software xi according to simulation inputs x andcarrying out simulation evaluation on the simulation road network, and outputting the flow of the obtained road section a in the time slice tAnd travel time
Step 3, initializing to obtain N charging schemes xiI is 1,2, …, N forms a set S, and M traffic demand parameters e are obtained through initializationjJ ═ 1,2, …, M constitutes set E; each charging plan and traffic demand parameter combination (x)i,ej) Performing simulation evaluation for a plurality of times by using macroscopic traffic simulation software, and solving the optimization model established in the step 2 by using a distributed robust simulation optimization algorithm to obtain a congestion charging scheme with optimal robustness;
the distributed robust simulation optimization algorithm adopted in the step 3 comprises the following specific processes:
step 3.1, for each charging plan in combination with the traffic demand parameter (x)i,ej) Firstly, using macroscopic traffic simulation software to carry out N0Simulation evaluation of/M times, N0The number of times of the preset initial simulation evaluation is obtained; then, the objective function value corresponding to each simulation is calculated according to the obtained flow and travel time, and the mean value and the variance of the objective function values of all simulation evaluations are calculated
Step 3.2, setting the initial value of the iteration number Num to be 0, and setting the maximum iteration number Num to be 0maxIs (B-N.N)0)/N0B is a preset upper limit of simulation evaluation times, and an initial value N of the simulation times of the distribution stagea(0) Setting the initial value of the algorithm convergence index as 0; when the convergence index of the algorithm does not reach the convergence threshold value and the simulation residual resources exist, the following steps are executed:
step 3.2.1, for each charging scenario xiAccording to the mean value and the variance of the objective function values corresponding to the combination of the objective function values and the traffic demand parameters, a random proxy model approximating the response relation between the objective function values and the traffic demand parameters is constructed to be used as a charging scheme xiThe first random agent model of (1);
step 3.2.2, searching the corresponding traffic demand parameter when the objective function value is maximum based on the first random agent model, and recording the traffic demand parameter as word-caseopt,i;
Step 3.2.3, set S and e according to the charging schemeopt,iConstructing a random agent model of a response relation between an approximate charging scheme and a worst-case objective function value as a second random agent model according to the corresponding objective function value;
step 3.2.4, based on the second random agent model, searching for a new charging plan sample point x*Updating congestion charging scheme set S ← S ×*,N←N+1;
Step 3.2.5, applying the same method to x as step 3.1*Combining with traffic demand parameters by Ns(Num +1)/M times of simulation evaluation, Ns(Num+1)=N0-Na(Num +1) corresponding to x*To communicate withA mean and variance of objective function values of the combination of demand parameters;
step 3.2.6, adopting a simulation resource allocation method to allocate the total simulation times N of the stages in the Num +1 iterationa(Num +1) assigning to all charging schemes and combinations of traffic demand parameters;
step 3.2.7, for each combination (x) according to the additional simulation times assigned in step 3.2.6i,ej) Executing additional simulation evaluation, and counting the mean value and the variance of the corresponding objective function value;
step 3.2.8, if the convergence condition is reached and the simulation residual resources are still available, updating iteration number Num ← Num +1, and returning to execute step 3.2.1; otherwise, outputting a scheme corresponding to the minimum worst-case objective function value in the existing sample set evaluated by the second random agent model as the optimal robust charging scheme. In a more preferred embodiment, an improved expectation-improvement criterion (abbreviated MEI) is employed and x is searched by maximizing the MEI*The expression of the MEI is:
in the formula (I), the compound is shown in the specification,andpredicted by the second random agent model at xiThe mean and variance of phi is a standard normal distribution function, and psi is a standard normal probability density function;is the minimum target value evaluated by the second random agent model in the existing charging scheme.
In a more preferred embodiment, step 3.2.7 provides a total number of simulations N to be assigned in the Num +1 iterationa(Num +1) is:
wherein N isminIs the preset minimum simulation times.
In a more preferred technical solution, step 3.2.1 adopts a stochastic kriging model to approximate the response relationship between the traffic demand parameter and the objective function value for each charging scheme, and/or step 3.2.3 adopts the stochastic kriging model to approximate the response relationship between the charging scheme and the worst-case objective function value.
In a more preferred technical scheme, the N charging schemes and the M traffic demand parameters obtained by initialization are obtained by a space uniform sampling method such as Latin hypercube sampling.
In a more preferred technical solution, the macroscopic traffic simulation software employs VISUM.
A distributed robust simulation optimization platform for road dynamic congestion charging comprises macroscopic traffic simulation software, a control platform and a database management system;
the database management system is to: storing a road network information file, a traffic demand file and a congestion charging scheme set;
the macroscopic traffic simulation software is used for: acquiring a road network information file from a database management system, and creating a corresponding simulation road network; acquiring a traffic demand file and a congestion charging scheme from a database management system, performing simulation evaluation, and outputting the flow and the travel time of each road section;
the control platform is used for: constructing a dynamic congestion charging optimization model of the road by taking the total travel time as a target evaluation index; initializing N charging schemes xiI is 1,2, …, N forms a set S, and M traffic demand parameters e are obtained through initializationjJ ═ 1,2, …, M constitutes set E; controlling the combination of the macroscopic traffic simulation software for each billing plan and traffic demand parameter (x)i,ej) And carrying out simulation evaluation for a plurality of times, and solving the optimization model by adopting a distributed robust simulation optimization algorithm. Updating the congestion charging scheme obtained from each iteration to the database management systemIn the system, the algorithm is ended to finally obtain a robust optimal congestion charging scheme;
the road dynamic congestion charging optimization model comprises the following steps:
wherein, S is a congestion charging scheme set, and x is any congestion charging scheme in S; j is the simulation times index, NxTotal number of simulations for x; t is a time segment index, and T is the total number of time segments; o is an initial node set, D is an end node set, OD is an origin and belongs to an OD pair set of which the O end belongs to D, and OD is any OD pair in the OD pair set; a is a road section set in the simulation road network, and a is any road section in the road section set A;obey normal distribution for the traffic demand of od in the time segment tAndmean and variance of the traffic demand of od at time segment t, respectively;traffic demand matrix for time segment t, consisting ofForming;the random traffic demand matrix obeys normal distribution in the jth simulation; theta is traffic demandThe set of fuzzy normal distributions that is satisfied,e is a parameter (traffic demand parameter for short) for controlling the traffic demand distribution variance, which represents the uncertainty of the traffic demand, e belongs to [ e [l,eu]Wherein e isl、euRespectively the upper and lower limits of e;representing macroscopic traffic simulation software xi according to simulation inputs x andcarrying out simulation evaluation on the simulation road network, and outputting the flow of the obtained road section a in the time slice tAnd travel time
An electronic device comprising a memory and a processor, wherein the memory stores a computer program, and when the computer program is executed by the processor, the processor is enabled to implement any one of the above-mentioned technical solutions of the distributed robust simulation optimization method.
A computer readable storage medium, on which a computer program is stored, wherein the computer program, when executed by a processor, implements any one of the above-mentioned technical solutions of the distributed robust simulation optimization method.
Advantageous effects
The method can solve the problem of simulation evaluation and optimization of dynamic crowded charging of roads under the condition of uncertain traffic demands, aims to improve the efficiency of a road network, can quickly and accurately optimize and obtain a robust crowded charging scheme by intelligently distributing limited simulation resources, and can evaluate and analyze the charging scheme by various traffic state indexes output by simulation.
Drawings
FIG. 1 is a schematic flow diagram of a distributed robust simulation optimization method according to the present invention;
FIG. 2 is a schematic diagram of a distributed robust simulation optimization system according to the present invention.
Detailed Description
The following describes embodiments of the present invention in detail, which are developed based on the technical solutions of the present invention, and give detailed implementation manners and specific operation procedures to further explain the technical solutions of the present invention.
Example 1
The embodiment provides a distributed robust simulation optimization method for road dynamic congestion charging, as shown in fig. 1, including the following steps:
step 1, for a road network to be subjected to dynamic congestion charging scheme optimization, a simulated road network is created in macroscopic traffic simulation software, and in this embodiment, VISUM is preferred by the macroscopic traffic simulation software.
Step 2, taking the total travel time as a target evaluation index, and constructing the following road dynamic congestion charging optimization model:
wherein, S is a congestion charging scheme set, and x is any congestion charging scheme in S; j is the simulation times index, NxTotal number of simulations for x; t is a time segment index, and T is the total number of time segments; o is an initial node set, D is an end node set, OD is an origin and belongs to an OD pair set of which the O end belongs to D, and OD is any OD pair in the OD pair set; a is a road section set in the simulation road network, and a is any road section in the road section set A;obey normal distribution for the traffic demand of od in the time segment tAndmean and variance of the traffic demand of od at time segment t, respectively;traffic demand matrix for time segment t, consisting ofForming;the random traffic demand matrix obeys normal distribution in the jth simulation; theta is traffic demandThe set of fuzzy normal distributions that is satisfied,e is a parameter (traffic demand parameter for short) for controlling the traffic demand distribution variance, which represents the uncertainty of the traffic demand, e belongs to [ e [l,eu]Wherein e isl、euRespectively the upper and lower limits of e;representing macroscopic traffic simulation software xi according to simulation inputs x andcarrying out simulation evaluation on the simulation road network, and outputting the flow of the obtained road section a in the time slice tAnd travel time
In this embodiment, a distributed robust congestion charging optimization model aiming at improving the road network efficiency is constructed, and the optimization of the charging level and the charging position is specifically as follows: the distributed robust congestion charging optimization model aiming at improving the road network efficiency is established by taking a charging level and a charging position as decision variables, traffic demand uncertain parameters as environment variables and a total travel time as an objective function, and relates to traffic demand uncertainty. Assuming that traffic demand follows a normal distributionAnd defining a fuzzy set based thereon WhereinThe value of (b) is determined by historical traffic demand data, so that a specific normal distribution can be determined by determining e.
The charging scheme comprises a charging level and a charging position, the current stage comprises a road section-based charging scheme and a region-based charging scheme, and other charging scheme forms appearing later are also within the protection scope of the patent.
Step 3, obtaining N charging schemes x by Latin hypercube sampling initializationiI is 1,2, …, N forms a set S, and M traffic demand parameters e are obtained through initializationjJ ═ 1,2, …, M constitutes set E; each charging plan and traffic demand parameter combination (x)i,ej) Performing simulation evaluation for a plurality of times by using macroscopic traffic simulation software, and solving the optimization model established in the step 2 by using a distributed robust simulation optimization algorithmAnd obtaining a robust optimal congestion charging scheme.
The specific process of solving and obtaining the optimal congestion scheme by adopting the distributed robust simulation optimization algorithm in the step 3 is as follows:
step 3.1, for each charging plan in combination with the traffic demand parameter (x)i,ej) Firstly, using macroscopic traffic simulation software to carry out N0Simulation evaluation of/M times, N0The number of times of the preset initial simulation evaluation is obtained. Then, the objective function value corresponding to each simulation is calculated according to the obtained flow and travel time, and the mean value and the variance of the objective function values of all the simulation evaluations are calculated
Step 3.2, setting the initial value of the iteration number Num to be 0, and setting the maximum iteration number Num to be 0maxIs (B-N.N)0)/N0Wherein B is the upper limit of simulation evaluation times and the initial value N of the simulation times of the distribution stagea(Num) is 0, the convergence index of the algorithm is set to be the MEI value, and the initial value is 0.5; when MEI is not less than the set algorithm convergence value (in the embodiment, the convergence threshold is set to be 0.01), and simulation residual resources (B-N.N) are also available0-Na(Num +1) > 0), the following steps are performed:
step 3.2.1, for each charging scenario xjAccording to the mean value and the variance of the objective function values corresponding to the combination of the objective function values and the traffic demand parameters, a random kriging model is adopted to construct a random proxy model approximating the response relation between the objective function values and the traffic demand parameters, and leave-one-out cross validation is performed to serve as xiThe first random proxy model of (1), denoted SSM-I.
Step 3.2.2, based on SSM-I, searching the corresponding traffic demand parameter when the objective function value is maximum, and recording the traffic demand parameter as word-caseopt,i。
Step 3.2.3, set S and e according to the charging schemeopt,iCorresponding objective function value, adopting random Kriging model to construct random proxy model of response relation between approximate charging scheme and worst-case objective function value, and executing leave-one cross validationIs a second random proxy model, denoted SSM-II.
Step 3.2.4, based on SSM-II, search for new charging plan sample points x with improved expected improvement criterion (MEI) and by maximizing MEI*Updating congestion charging scheme set S ← S ×*N ← N + 1. Other methods suitable for searching for new sample points by the random proxy model are also within the scope of the present patent.
Step 3.2.5, applying the same method to x as step 3.1*Combining with traffic demand parameters by Ns(Num +1)/M times of simulation evaluation processing are carried out, and a charging scheme x is correspondingly obtained*Mean and variance of objective function values combined with traffic demand parameters. Wherein N iss(Num+1)=N0-Na(Num+1)。
Step 3.2.6, adopting a simulation resource allocation method to allocate the total simulation times N of the stages in the Num +1 iterationa(Num +1) is assigned to all charging schemes in combination with the traffic demand parameter, so that simulation resources can be assigned to more promising charging schemes, i.e. charging schemes with small mean or large variance. The simulation resource allocation method includes, but not limited to, toa (theoretical optimal allocation), ea (equivalent allocation), ptv (procedural to variance), ocba (optimal computing budget), and ROCBA (Ramdomized optimal computing budget).
Wherein the total simulation total times N distributed in the Num +1 th iterationa(Num +1) is:
wherein N isminA preset minimum number of simulations of the congestion charging scheme.
Step 3.2.7, for each combination (x) according to the additional simulation times assigned in step 3.2.6i,ej) Additional simulation evaluations were performed and the mean and variance of the corresponding objective function values were counted.
3.2.8, if the convergence condition is not reached, and furtherWhen simulation residual resources exist, updating iteration times Num ← Num +1, and returning to the step 3.2.1; otherwise, use each combination (x)i,ej) And outputting a scheme corresponding to the minimum worst-case objective function value in the existing charging scheme set evaluated by the second random agent model as the optimal robust charging scheme.
In this embodiment, a new charging plan is searched for using the improved desired improvement criterion, and therefore, the convergence condition may be set as whether the MEI value is smaller than a set value, for example, 0.01, and when the improved charging plan cannot be found, the algorithm is terminated.
Example 2
The embodiment provides a distributed robust simulation optimization platform for road dynamic congestion charging, which, as shown in fig. 2, includes a macroscopic traffic simulation software, a control platform and a database management system;
the database management system is to: storing a road network information file, a traffic demand file and a congestion charging scheme set;
the macroscopic traffic simulation software is used for: acquiring a road network information file from a database management system, and creating a corresponding simulation road network; acquiring a traffic demand file and a congestion charging scheme from a database management system, performing simulation evaluation, and outputting the flow and the travel time of each road section;
the control platform is used for: constructing a dynamic congestion charging optimization model of the road by taking the total travel time as a target evaluation index; initializing N charging schemes xiI is 1,2, …, N forms a set S, and M traffic demand parameters e are obtained through initializationjJ ═ 1,2, …, M constitutes set E; controlling the combination of the macroscopic traffic simulation software for each billing plan and traffic demand parameter (x)i,ej) And carrying out simulation evaluation for a plurality of times, and solving the optimization model by adopting a distributed robust simulation optimization algorithm. Updating the congestion charging scheme obtained by each iteration into a database management system, and finally obtaining the congestion charging scheme with optimal robustness after the algorithm is ended;
the road dynamic congestion charging optimization model comprises the following steps:
wherein, S is a congestion charging scheme set, and x is any congestion charging scheme in S; j is the simulation times index, NxTotal number of simulations for x; t is a time segment index, and T is the total number of time segments; o is an initial node set, D is an end node set, OD is an origin and belongs to an OD pair set of which the O end belongs to D, and OD is any OD pair in the OD pair set; a is a road section set in the simulation road network, and a is any road section in the road section set A;obey normal distribution for the traffic demand of od in the time segment tAndmean and variance of the traffic demand of od at time segment t, respectively;a traffic demand matrix for a time segment t, consisting ofForming;the random traffic demand matrix obeys normal distribution in the jth simulation; theta is traffic demandThe set of fuzzy normal distributions that is satisfied,e is a parameter (traffic demand parameter for short) for controlling the traffic demand distribution variance, which represents the uncertainty of the traffic demand, e belongs to [ e [l,eu]Wherein e isl、euUpper and lower limits of e, respectively;representing macroscopic traffic simulation software xi according to simulation inputs x andcarrying out simulation evaluation on the simulation road network, and outputting the flow of the obtained road section a in the time slice tAnd travel time
The control platform described in this embodiment obtains the optimal congestion charging scheme by solving the dynamic congestion charging optimization model of the road by using a distributed robust simulation optimization algorithm, and the specific solving method is the same as that described in embodiment 1, and is not described again in this embodiment.
In a more preferred embodiment, the macro traffic simulation software may adopt VISUM, the control platform may adopt MATLAB, Python, or the like, and the database management system may adopt ACCESS, and the present invention does not limit the specific types.
Example 3
The present embodiment provides an electronic device, which includes a memory and a processor, where the memory stores a computer program, and when the computer program is executed by the processor, the processor implements the distributed robust simulation optimization method described in embodiment 1.
Example 4
The present embodiment provides a computer-readable storage medium, on which a computer program is stored, which, when executed by a processor, implements the distributed robust simulation optimization method of embodiment 1.
The above embodiments are preferred embodiments of the present application, and those skilled in the art can make various changes or modifications without departing from the general concept of the present application, and such changes or modifications should fall within the scope of the claims of the present application.
Claims (9)
1. A distributed robust simulation optimization method for road dynamic congestion charging is characterized by comprising the following steps:
step 1, for a road network to be subjected to dynamic congestion charging scheme optimization, creating a simulation road network in macroscopic traffic simulation software;
step 2, taking the total travel time as a target evaluation index, and constructing the following road dynamic congestion charging optimization model:
wherein, S is a congestion charging scheme set, and x is any congestion charging scheme in S; j is simulation frequency index, NxTotal number of simulations for x; t is a time segment index, and T is the total number of time segments; o is an initial node set, D is an end node set, OD is an origin and belongs to an OD pair set of which the O end belongs to D, and OD is any OD pair in the OD pair set; a is a road section set in the simulation road network, and a is any road section in the road section set A;obey normal distribution for the traffic demand of od in the time segment t Andmean and variance of the traffic demand of od at time segment t, respectively;a traffic demand matrix for a time segment t, consisting ofForming;the random traffic demand matrix obeys normal distribution in the jth simulation; theta is traffic demandThe set of fuzzy normal distributions that is satisfied,e is a traffic demand parameter which represents the uncertainty of the traffic demand, and e belongs to [ e ∈ ]l,eu]Wherein e isl、euUpper and lower limits of e, respectively;representing the simulation input x and x of macroscopic traffic simulation software xi in a simulation road networkCarrying out simulation evaluation, and outputting the flow of the obtained road section a in the time segment tAnd travel time
Step 3, initializing to obtain N charging schemes xiI is 1,2, …, N forms a set S, and M traffic demand parameters e are obtained through initializationjJ ═ 1,2, …, M constitutes set E; each charging plan and traffic demand parameter combination (x)i,ej) Performing simulation evaluation for a plurality of times by using macroscopic traffic simulation software, and solving the optimization model established in the step 2 by using a distributed robust simulation optimization algorithm to obtain a congestion charging scheme with optimal robustness;
the distributed robust simulation optimization algorithm adopted in the step 3 comprises the following specific processes:
step 3.1, for each charging plan in combination with the traffic demand parameter (x)i,ej) Firstly, using macroscopic traffic simulation software to carry out N0Simulation evaluation of/M times, N0The number of times of the preset initial simulation evaluation is obtained; calculating objective function values corresponding to each simulation according to the flow and travel time obtained by simulation evaluation, and then calculating the mean value and variance of all objective function values
Step 3.2, setting the initial value of the iteration number Num to be 0, and setting the maximum iteration number Num to be 0maxIs (B-N.N)0)/N0B is a preset simulation evaluation frequency upper limit, and simulation frequency N of the distribution stagea(Num), setting an initial value of an algorithm convergence index; when the convergence index of the algorithm does not reach the convergence threshold value and the simulation residual resources exist, the following steps are executed:
step 3.2.1, for each charging scenario xiAccording to the mean value and the variance of the objective function values corresponding to the combination of the objective function values and the traffic demand parameters, a random proxy model approximating the response relation between the objective function values and the traffic demand parameters is constructed to be used as a charging scheme xiThe first random agent model of (1);
step 3.2.2, searching the corresponding traffic demand parameter when the objective function value is maximum based on the first random agent model, and recording the traffic demand parameter as worst-case trafficRequirement parameter eopt,i;
Step 3.2.3, set S and e according to the charging schemeopt,iConstructing a random agent model of a response relation between an approximate charging scheme and a worst-case objective function value as a second random agent model according to the corresponding objective function value;
step 3.2.4, search for a new charging scenario x based on the second stochastic proxy model*Updating congestion charging scheme set S ← S ×*,N←N+1;
Step 3.2.5, applying the same method to x as step 3.1*Combining with traffic demand parameters by Ns(Num +1)/M times of simulation evaluation, Ns(Num+1)=N0-Na(Num +1) corresponding to x*A mean and variance of objective function values combined with the traffic demand parameters;
step 3.2.6, adopting a simulation resource allocation method to allocate the total simulation times N of the stages in the Num +1 iterationa(Num +1) assigning to all charging schemes and combinations of traffic demand parameters;
step 3.2.7, for each combination (x) according to the additional simulation times assigned in step 3.2.6i,ej) Executing additional simulation evaluation, and counting the mean value and the variance of the corresponding objective function value;
step 3.2.8, if the convergence condition is not reached and there are simulation remaining resources, updating iteration number Num ← Num +1, and returning to execute step 3.2.1; otherwise, outputting a scheme corresponding to the minimum worst-case objective function value in the existing sample set evaluated by the second random agent model as the optimal robust charging scheme.
2. The method of claim 1, wherein an improved expected improvement criterion, abbreviated MEI, is employed and a new charging scheme is searched by maximizing the MEI, wherein an algorithm convergence index is a MEI value, and wherein the MEI is expressed as:
in the formula (I), the compound is shown in the specification,andpredicted by the second random agent model at xiThe mean and variance of phi is a standard normal distribution function, and psi is a standard normal probability density function;is the minimum target value evaluated by the second random agent model in the existing charging scheme.
4. A method according to claim 1, characterised in that step 3.2.1 approximates the response between the traffic demand parameter and the objective function value for each charging scheme using a stochastic kriging model, and/or step 3.2.3 approximates the response between the charging scheme and the worst-case objective function value using a stochastic kriging model.
5. Method according to claim 1, characterized in that the obtained N charging schemes x are initializediAnd M traffic demand parameters ejAll obtained by a spatially uniform sampling method.
6. The method of claim 1, wherein the macro traffic simulation software employs VISUM.
7. A distributed robust simulation optimization platform for road dynamic congestion charging is characterized by comprising macroscopic traffic simulation software, a control platform and a database management system;
the database management system is to: storing a road network information file, a traffic demand file and a congestion charging scheme set;
the macroscopic traffic simulation software is used for: acquiring a road network information file from a database management system, and creating a corresponding simulation road network; acquiring a traffic demand file and a congestion charging scheme from a database management system, performing simulation evaluation, and outputting the flow and the travel time of each road section;
the control platform is used for: constructing a dynamic congestion charging optimization model of the road by taking the total travel time as a target evaluation index; initializing N charging schemes xiI is 1,2, …, N forms a set S, and M traffic demand parameters e are obtained through initializationjJ ═ 1,2, …, M constitutes set E; controlling the combination of the macroscopic traffic simulation software for each billing plan and traffic demand parameter (x)i,ej) Carrying out simulation evaluation for a plurality of times, and solving an optimization model by adopting a distributed robust simulation optimization algorithm; updating the congestion charging scheme obtained by each iteration into a database management system, and finally obtaining the congestion charging scheme with optimal robustness after the algorithm is ended;
the road dynamic congestion charging optimization model comprises the following steps:
wherein, S is a congestion charging scheme set, and x is any congestion charging scheme in S; j is the number of times of simulationNumber index, NxTotal number of simulations for x; t is a time segment index, and T is the total number of time segments; o is an initial node set, D is an end node set, OD is an origin and belongs to an OD pair set of which the O end belongs to D, and OD is any OD pair in the OD pair set; a is a road section set in the simulation road network, and a is any road section in the road section set A;obey normal distribution for the traffic demand of od in the time segment t Andmean and variance of the traffic demand of od at time segment t, respectively;a traffic demand matrix for a time segment t, consisting ofForming;the random traffic demand matrix obeys normal distribution in the jth simulation; theta is the traffic demandThe set of fuzzy normal distributions that is satisfied,e is a parameter for controlling traffic demand distribution variance, referred to as traffic demand parameter for short, and represents the uncertainty of traffic demand, e belongs to [ e [ ]l,eu]Wherein e isl、euUpper and lower limits of e, respectively;representing macroscopic traffic simulation software xi according to simulation inputs x andcarrying out simulation evaluation on the simulation road network, and outputting the flow of the obtained road section a in the time slice tAnd travel time
8. An electronic device comprising a memory and a processor, the memory having stored therein a computer program, wherein the computer program, when executed by the processor, causes the processor to implement the method of any of claims 1-6.
9. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the method according to any one of claims 1 to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110494686.4A CN113177320B (en) | 2021-05-07 | 2021-05-07 | Simulation optimization method, system, equipment and medium for road dynamic congestion charging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110494686.4A CN113177320B (en) | 2021-05-07 | 2021-05-07 | Simulation optimization method, system, equipment and medium for road dynamic congestion charging |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113177320A CN113177320A (en) | 2021-07-27 |
CN113177320B true CN113177320B (en) | 2022-06-14 |
Family
ID=76928256
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110494686.4A Active CN113177320B (en) | 2021-05-07 | 2021-05-07 | Simulation optimization method, system, equipment and medium for road dynamic congestion charging |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113177320B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014187476A1 (en) * | 2013-05-21 | 2014-11-27 | Nec Europe Ltd. | Method and system for predicting mobility demand of users |
CN107293115A (en) * | 2017-05-09 | 2017-10-24 | 上海电科智能系统股份有限公司 | A kind of traffic flow forecasting method for microscopic simulation |
CN110380409A (en) * | 2019-07-16 | 2019-10-25 | 山东大学 | Consider the active distribution network distributed robust state estimation method and system of communication failure |
CN110661258A (en) * | 2019-09-29 | 2020-01-07 | 广东电网有限责任公司 | Flexible resource distributed robust optimization method for power system |
CN112347615A (en) * | 2020-10-20 | 2021-02-09 | 天津大学 | Power distribution network hybrid optimization scheduling method considering light storage and fast charging integrated station |
DE102020203199A1 (en) * | 2020-03-12 | 2021-04-15 | Vitesco Technologies GmbH | Method and system for charging and load sharing |
CN112701687A (en) * | 2021-01-25 | 2021-04-23 | 四川大学 | Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10382285B2 (en) * | 2011-08-25 | 2019-08-13 | Siemens Industry, Inc. | Smart grid communication assessment and co-simulation tool |
US10986516B2 (en) * | 2017-03-10 | 2021-04-20 | Huawei Technologies Co., Ltd. | System and method of network policy optimization |
US11010503B2 (en) * | 2018-05-15 | 2021-05-18 | Tata Consultancy Services Limited | Method and system providing temporal-spatial prediction of load demand |
-
2021
- 2021-05-07 CN CN202110494686.4A patent/CN113177320B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014187476A1 (en) * | 2013-05-21 | 2014-11-27 | Nec Europe Ltd. | Method and system for predicting mobility demand of users |
CN107293115A (en) * | 2017-05-09 | 2017-10-24 | 上海电科智能系统股份有限公司 | A kind of traffic flow forecasting method for microscopic simulation |
CN110380409A (en) * | 2019-07-16 | 2019-10-25 | 山东大学 | Consider the active distribution network distributed robust state estimation method and system of communication failure |
CN110661258A (en) * | 2019-09-29 | 2020-01-07 | 广东电网有限责任公司 | Flexible resource distributed robust optimization method for power system |
DE102020203199A1 (en) * | 2020-03-12 | 2021-04-15 | Vitesco Technologies GmbH | Method and system for charging and load sharing |
CN112347615A (en) * | 2020-10-20 | 2021-02-09 | 天津大学 | Power distribution network hybrid optimization scheduling method considering light storage and fast charging integrated station |
CN112701687A (en) * | 2021-01-25 | 2021-04-23 | 四川大学 | Robust optimization operation method of gas-electricity distribution network system considering price type combined demand response |
Non-Patent Citations (6)
Title |
---|
Detailed string stability analysis for bi-directional optimal velocity model;郑亮;《中南大学学报(英文版)》;20150422(第4期);第1563-1573页 * |
城市路网的分布式鲁棒预测控制;刘安东等;《浙江工业大学学报》;20161225(第06期);第633-638页 * |
基于NetLogo的地铁车站人员紧急疏散仿真研究;胡丽娟等;《铁道科学与工程学报》;20171215(第12期);第2730-2737页 * |
基于效率-安全评价的近距离错位交叉口信控优化;罗慧敏等;《交通信息与安全》;20180228(第01期);第65-73页 * |
考虑风电高阶不确定性的分布式鲁棒优化调度模型;夏鹏等;《电工技术学报》;20200110(第01期);第189-200页 * |
面向大型社会活动的快速路网控制策略仿真评价方法;杨珍珍等;《计算机应用研究》;20101215(第12期);第4473-4475页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113177320A (en) | 2021-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Task scheduling algorithm based on improved firework algorithm in fog computing | |
Liang et al. | An integrated reinforcement learning and centralized programming approach for online taxi dispatching | |
CN108021451B (en) | Self-adaptive container migration method in fog computing environment | |
Zhou et al. | Multi‐agent model‐based predictive control for large‐scale urban traffic networks using a serial scheme | |
CN106134136A (en) | Calculate the long-term dispatch transmitted for the data on wide area network | |
WO2024087512A1 (en) | Graph neural network compression method and apparatus, and electronic device and storage medium | |
CN109788489A (en) | A kind of base station planning method and device | |
CN111586146B (en) | Wireless internet of things resource allocation method based on probability transfer deep reinforcement learning | |
CN103957261A (en) | Cloud computing resource distributing method based on energy consumption optimization | |
CN111176784B (en) | Virtual machine integration method based on extreme learning machine and ant colony system | |
CN111694664A (en) | Calculation unloading distribution method of edge server | |
CN115860081B (en) | Core algorithm scheduling method, system, electronic equipment and storage medium | |
WO2023134403A1 (en) | Internet of things resource allocation method and system, terminal and storage medium | |
CN117094535B (en) | Artificial intelligence-based energy supply management method and system | |
CN108106624A (en) | A kind of more people's Dispatch by appointment paths planning methods and relevant apparatus | |
Li et al. | Resource usage prediction based on BiLSTM-GRU combination model | |
Tungom et al. | Hierarchical framework for demand prediction and iterative optimization of EV charging network infrastructure under uncertainty with cost and quality-of-service consideration | |
Xie et al. | A censored semi-bandit model for resource allocation in bike sharing systems | |
Cheng et al. | A novel task provisioning approach fusing reinforcement learning for big data | |
CN117112703B (en) | Space planning stock unit identification method based on multidimensional analysis | |
CN113177320B (en) | Simulation optimization method, system, equipment and medium for road dynamic congestion charging | |
Wang et al. | A study of situation awareness-based resource management scheme in cloud environment | |
Zhang et al. | Offloading demand prediction-driven latency-aware resource reservation in edge networks | |
CN115691140B (en) | Analysis and prediction method for space-time distribution of automobile charging demand | |
Zhang et al. | Energy efficient federated learning in internet of vehicles: A game theoretic scheme |
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 |