一种道路动态拥挤收费的分布式鲁棒仿真优化方法、系统、设备及介质技术方案

本发明专利技术公开了一种道路动态拥挤收费的分布式鲁棒仿真优化方法、系统、设备及介质，其方法为：步骤1，对于待进行动态拥挤收费方案优化的路网，在宏观交通仿真软件中创建仿真路网；步骤2，以总旅行时间为目标评价指标，构建道路动态拥挤收费优化模型；步骤3，初始化收费方案集和交通需求参数集，使用宏观交通仿真软件对创建的仿真路网进行仿真评价，采用分布式鲁棒仿真优化算法求解步骤2建立的优化模型，得到鲁棒最优的拥挤收费方案。本发明专利技术可以解决交通需求不确定性条件下的道路动态拥挤收费仿真评价与优化问题。

【技术实现步骤摘要】
一种道路动态拥挤收费的分布式鲁棒仿真优化方法、系统、设备及介质
本专利技术涉及道路动态拥挤收费领域，特别涉及一种道路动态拥挤收费的分布式鲁棒仿真优化方法、系统、设备及介质。
技术介绍
交通拥挤收费，也即在特定时间段内对交通出行施加额外出行费用的管控措施，通过增加出行者使用收费区域交通设施的成本，诱导出行者重新做出出行时间和路径选择，实现部分出行者从拥堵区域到非拥堵区域的转移，或者选择在非收费时段出行，从而在一定程度上降低收费区域的交通需求，缓解交通供给与交通需求之间的矛盾，缓解交通拥挤，使路网资源得到更加的充分的利用，对城市健康发展具有积极意义。交通仿真技术即构建模型来模拟已有的或者正在规划中的交通系统，便于用户对研究的交通网络系统进行优化，或者预测未来的交通状态。然而，交通仿真技术也存在黑箱性与耗时长等问题。面向这些问题，仿真优化算法应运而生，大大提高了计算速度。主流的仿真优化方法主要有基于梯度的搜索法(Gradientbasedsearchmethods)、随机近似法(Stochasticapproximationmethods)、启发算法(Heuristicsearchmethods)和响应面法(Responsesurfacemethods)。响应面法即代理模型法具有节省仿真成本和拟合复杂未知的系统输入输出关系的优势，这种优势在仿真实验耗时的情况下更为突出。在复杂的现实工程问题中，即使忽略数据的不确定性影响，按照假设的条件求得的最优解，通常很可能是不稳定甚至不可行的。因此，相较于全局最优解，在实际问题中求解对扰动敏感较小的解，即分布式鲁棒解，将更具现实意义。
技术实现思路
本专利技术提供一种道路动态拥挤收费的分布式鲁棒仿真优化方法、系统、设备及介质，可以解决交通需求不确定性条件下的道路动态拥挤收费仿真与评价问题。为实现上述技术目的，本专利技术采用如下技术方案：一种道路动态拥挤收费的分布式鲁棒仿真优化方法，包括以下步骤：步骤1，对于待进行动态拥挤收费方案优化的路网，在宏观交通仿真软件中创建仿真路网；步骤2，以总旅行时间为目标评价指标，构建以下道路动态拥挤收费优化模型：式中，S为拥挤收费方案集，x为S中的任意拥挤收费方案；j为仿真次数指标，Nx为x的仿真总次数；t为时间片段指标，T为时间片段总数；O为起始节点集合，D为终点节点集合，OD为起点属于O终点属于D的OD对集合，od为OD对集合OD中的任意OD对；A为仿真路网中的路段集合，a为路段集合A中的任意路段；为od在时间片段t的交通需求，服从正态分布和分别为od在时间片段t的交通需求的均值和方差；为时间片段t的交通需求矩阵，由构成；是第j次仿真时服从正态分布的随机交通需求矩阵；Θ为交通需求满足的模糊正态分布集，e是控制交通需求分布方差的参数(简称交通需求参数)，表示交通需求不确定的程度，e∈[el，eu]，其中el、eu分别为e的上下限；表示宏观交通仿真软件Ξ根据仿真输入x和对仿真路网进行仿真评价，输出得到路段a在时间片段t的流量和旅行时间步骤3，初始化获得N个收费方案xi，i＝1，2，…，N构成集合S，初始化获得M个交通需求参数ej，j＝1，2，…，M构成集合E；每个收费方案和交通需求参数的组合(xi，ej)，均使用宏观交通仿真软件进行若干次仿真评价，采用分布式鲁棒仿真优化算法求解步骤2建立的优化模型，得到鲁棒最优的拥挤收费方案。在更优的技术方案中，步骤3中采用的分布式鲁棒仿真优化算法具体过程为：步骤3.1，对每个收费方案与交通需求参数的组合(xi，ej)，先使用宏观交通仿真软件进行N0/M次仿真评价，N0为预设的初始仿真评价次数；然后根据得到的流量和旅行时间计算每次仿真对应的目标函数值，再计算所有仿真评价的目标函数值的均值与方差步骤3.2，令迭代次数Num的初始值为0，最大迭代次数Nummax为(B-N·N0)/N0，其中B为预设的仿真评价次数上限，分配阶段的仿真次数的初始值Na(0)＝0，设置算法收敛指标的初始值；当算法收敛指标未达到收敛阈值，且还有仿真剩余资源时，执行以下步骤：步骤3.2.1，对于每个收费方案xi，根据其与交通需求参数组合对应的目标函数值的均值与方差，构建近似目标函数值与交通需求参数之间响应关系的随机代理模型，作为收费方案xi的第一随机代理模型；步骤3.2.2，基于第一随机代理模型，搜索目标函数值最大时所对应的交通需求参数，记为worst-case的交通需求参数eopt，i；步骤3.2.3，根据收费方案集S和eopt，i对应的目标函数值，构建近似收费方案与worst-case目标函数值之间响应关系的随机代理模型，作为第二随机代理模型；步骤3.2.4，基于第二随机代理模型，搜索新的收费方案样本点x*，更新拥挤收费方案集S←SUx*，N←N+1；步骤3.2.5，采用与步骤3.1相同的方法对x*与交通需求参数的组合进行Ns(Num+1)/M次仿真评价，Ns(Num+1)＝N0-Na(Num+1)，对应得到x*与交通需求参数的组合的目标函数值的均值与方差；步骤3.2.6，采用仿真资源分配方法，在第Num+1次迭代中将分配阶段总共的仿真次数Na(Num+1)分配给所有收费方案与交通需求参数的组合；步骤3.2.7，按照步骤3.2.6分配的额外仿真次数，对每个组合(xi，ej)执行额外的仿真评价，并统计相应的目标函数值的均值和方差；步骤3.2.8，若达到收敛条件，且还有仿真剩余资源时，更新迭代次数Num←Num+1，返回执行步骤3.2.1；否则，输出第二随机代理模型评估的已有样本集中最小worst-case目标函数值对应的方案作为鲁棒最优的收费方案。在更优的技术方案中，采用改进的期望改善准则(简称MEI)，并通过最大化MEI来搜索x*，所述MEI的表达式为：式中，和是由第二随机代理模型预测的在xi处的均值与方差，Φ是标准正态分布函数，Ψ是标准正态概率密度函数；是已有收费方案中第二随机代理模型评估的最小目标值。在更优的技术方案中，步骤3.2.7在第Num+1次迭代中可供分配的总共仿真次数Na(Num+1)为：其中，Nmin为预设的最小仿真次数。在更优的技术方案中，步骤3.2.1采用随机克里金模型对于每一个收费方案近似交通需求参数与目标函数值之间的响应关系，和/或，步骤3.2.3采用随机克里金模型近似收费方案与worst-case目标函数值之间的响应关系。在更优的技术方案中，初始化获得的N个收费方案和M个交通需求参数，均是通过空间均匀抽样方法例如拉丁超立方抽样获得。在更优的技术方案中，所述宏观交通仿真软件采用VISUM。一种道路动态拥挤收费的分布式鲁棒仿真优化平台，包括宏观交通仿真软件、控制平台和数据库管理系统；所述数据库管理系统用于：存储路网信息文件、交通需求文件以及拥挤收费方案集；所述宏观交通仿真软件用于：从数本文档来自技高网...

【技术保护点】
1.一种道路动态拥挤收费的分布式鲁棒仿真优化方法，其特征在于，包括以下步骤：/n步骤1，对于待进行动态拥挤收费方案优化的路网，在宏观交通仿真软件中创建仿真路网；/n步骤2，以总旅行时间为目标评价指标，构建以下道路动态拥挤收费优化模型：/n

【技术特征摘要】
1.一种道路动态拥挤收费的分布式鲁棒仿真优化方法，其特征在于，包括以下步骤：
步骤1，对于待进行动态拥挤收费方案优化的路网，在宏观交通仿真软件中创建仿真路网；
步骤2，以总旅行时间为目标评价指标，构建以下道路动态拥挤收费优化模型：






式中，S为拥挤收费方案集，x为S中的任意拥挤收费方案；j为仿真次数指标，Nx为x的仿真总次数；t为时间片段指标，T为时间片段总数；O为起始节点集合，D为终点节点集合，OD为起点属于O终点属于D的OD对集合，od为OD对集合OD中的任意OD对；A为仿真路网中的路段集合，a为路段集合A中的任意路段；为od在时间片段t的交通需求，服从正态分布和分别为od在时间片段t的交通需求的均值和方差；为时间片段t的交通需求矩阵，由构成；是第j次仿真时服从正态分布的随机交通需求矩阵；Θ为交通需求满足的模糊正态分布集，e为交通需求参数，表示交通需求不确定的程度，e∈[el,eu]，其中el、eu分别为e的上下限；表示宏观交通仿真软件Ξ在仿真路网中对仿真输入x和进行仿真评价，输出得到路段a在时间片段t的流量和旅行时间
步骤3，初始化获得N个收费方案xi,i＝1,2,…,N构成集合S，初始化获得M个交通需求参数ej,j＝1,2,…,M构成集合E；每个收费方案和交通需求参数的组合(xi,ej)，均使用宏观交通仿真软件进行若干次仿真评价，采用分布式鲁棒仿真优化算法求解步骤2建立的优化模型，得到鲁棒最优的拥挤收费方案。


2.根据权利要求1所述的方法，其特征在于，步骤3中采用的分布式鲁棒仿真优化算法具体过程为：
步骤3.1，对每个收费方案与交通需求参数的组合(xi,ej)，先使用宏观交通仿真软件进行N0/M次仿真评价，N0为预设的初始仿真评价次数；根据仿真评价得到的流量和旅行时间计算每次仿真对应的目标函数值，再计算所有目标函数值的均值与方差
步骤3.2，令迭代次数Num的初始值为0，最大迭代次数Nummax为(B-N·N0)/N0，其中B为预设的仿真评价次数上限，分配阶段的仿真次数Na(Num)，设置算法收敛指标的初始值；当算法收敛指标未达到收敛阈值，且还有仿真剩余资源时，执行以下步骤：
步骤3.2.1，对于每个收费方案xi，根据其与交通需求参数组合对应的目标函数值的均值与方差，构建近似目标函数值与交通需求参数之间响应关系的随机代理模型，作为收费方案xi的第一随机代理模型；
步骤3.2.2，基于第一随机代理模型，搜索目标函数值最大时所对应的交通需求参数，记为worst-case的交通需求参数eopt,i；
步骤3.2.3，根据收费方案集S和eopt,i对应的目标函数值，构建近似收费方案与worst-case目标函数值之间响应关系的随机代理模型，作为第二随机代理模型；
步骤3.2.4，基于第二随机代理模型，搜索新的收费方案x*，更新拥挤收费方案集S←S∪x*,N←N+1；
步骤3.2.5，采用与步骤3.1相同的方法对x*与交通需求参数的组合进行Ns(Num+1)/M次仿真评价，Ns(Num+1)＝N0-Na(Num+1)，对应得到x*与交通需求参数组合的目标函数值均值与方差；
步骤3.2.6，采用仿真资源分配方法，在第Num+1次迭代中将分配阶段总共的仿真次数Na(Num+1)分配给所有收费方案与交通需求参数的组合；
步骤3.2.7，按照步骤3.2.6分配的额外仿真次数，对每个组合(xi，ej)执行额外仿真评价，并统计相应的目标函数值的均值和方差；
步骤3.2.8，若未达到收敛条件，且还有仿真剩余资源时，更新迭代次数Num←Num+1，返回执行步骤3.2.1；否则，输出第二随机代理模型评估的已有样本集中...

【专利技术属性】
技术研发人员：郑亮，冯敏，黄会敏，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南;43