考虑交通冲突指标分布问题的交通仿真模型标定方法与流程

本发明涉及一种考虑交通冲突指标分布多约束优化问题的交通仿真模型标定方法，属于智能交通
技术领域：
。
背景技术：
：微观交通仿真模型运用大量的独立参数来描述交通系统运行、交通流特征以及驾驶员行为等，而模型的参数取值对仿真结果有较大的影响。微观交通仿真模型标定是指通过对模型参数的调整来提高模型拟合现实交通运行状况及车辆驾驶行为的精确度。模型标定所要达到的目标是使得仿真所得到的评价现实交通运行状况的指标值与实测值间的相对误差达到最小，以此再现现实交通运行状况。在微观交通仿真系统中，采取参数默认值而仿真得到的目标值与实际值的误差甚至达35％左右。而且在微观交通仿真模型中许多参数的取值无法通过实际调查获得，因而需要对模型参数进行标定。文章《MicroscopicSimulationModelCalibrationandValidation:ACaseStudyofVISSIMforaCoordinatedActuatedSignalSystem》中提出了微观交通仿真模型标定流程的九个步骤，即选取有效的性能评价指标数据采集确定标定参数实验设计仿真运行建立评价函数设置遗传算法标定结果评价选用新数据进行模型验证。该文以交通仿真模型为例，对这九个步骤进行了详细的说明。文中选取行程时间作为评价指标，采用最大排队长度作为验证指标。文章《基于遗传算法的交通仿真模型参数校正方法研究》提出了基于遗传算法的交通仿真模型参数校正方法。并根据实测数据，以仿真平台为例，对模型进行参数校正。研究给出了仿真模型的校正参数及参数的取值范围。文章《微观交通仿真系统参数校正研究》在之前研究的基础上，以合肥市大东门区域VISSIM仿真系统模型为例，建立了仿真模型参数校正的遗传模拟退火启发式算法，实现了对VISSIM的仿真参数的自动化校正。目前已有的仿真模型标定方法主要有两类，一类是基于效率指标的仿真模型标定方法，这种方法较为传统，一般选择效率指标如速度、行程时间、排队长度等进行参数标定；另一类是近年来提出的基于效率指标以及交通冲突的两阶段模型标定方法。这种方法主要是首先按照传统方法的标定思路对仿真模型的相关参数进行标定，接着以交通冲突的数量作为第二阶段的标定目标，通过标定相关的仿真参数，使得仿真冲突与实际观测冲突的数量接近，从而完成仿真模型的标定并用于交通安全分析。现有标定方法的主要缺点：1.传统的基于效率指标的仿真模型标定方法没有考虑安全性指标；2.考虑交通冲突数量标定的交通仿真模型并不一定能很好地反映真实的安全状况。交通冲突包含很多的不同参数，这些参数当中有一些已经被证明与交通安全具有非常强的关联性，例如碰撞时间TTC,后侵占时间PET。这些参数应该在标定过程中被仔细考虑，标定后的模型才能最大程度地反映真实的安全状况。技术实现要素：本发明为克服现有技术中存在的技术缺陷，提出一种考虑交通冲突指标分布问题的交通仿真模型标定方法，该方法针对交通冲突的重要指标，使仿真模型能够更好地服务于交通安全分析与评价。为解决上述技术问题，本发明提供了考虑交通冲突指标分布问题的交通仿真模型标定方法，选择交叉口效率指标作为交通仿真建模第一轮标定目标以及第二轮标定目标的约束条件，选择交叉口安全指标作为交通仿真建模第二轮标定目标；通过计算仿真效率指标的标准差以及选择置信区间对建立的仿真模型的最小所需仿真次数进行估计，筛选出效率指标和安全指标的仿真参数；建立多约束的标定优化问题，通过遗传算法对模型参数进行标定，选择适应度参数和模型约束条件，通过对于交通冲突参数分布进行模型标定。本发明具体包括以下步骤：S1.收集交通仿真建模所需数据，包括道路几何特征数据、交通控制数据、交通管理数据、交通流量数据和交通冲突指标；S2.根据步骤S1收集到的数据，进行交通仿真建模；S3.选择交叉口效率指标作为第一轮标定目标以及第二轮标定目标的约束条件，该效率指标包括交叉口平均排队长度、控制延误和两点间车辆行程时间；S4.选择交叉口安全指标作为第二轮标定目标，该安全指标包括交通冲突的避免碰撞时间TTC指标和后侵占时间PET指标；S5.对建立的仿真模型的最小所需仿真次数进行估计，先预设仿真次数，提取仿真输出的效率指标和安全指标；S6.计算效率指标中平均排队长度的标准差；其中，x为平均排队长度，N为仿真次数；S7.选择置信区间计算该仿真模型所需的最小仿真次数：其中，C为1-置信度，t为对应置信度的分位点，s为平均排队长度的标准差；S8.缺省参数的适用性检验：在不调整缺省参数的情况下，运行仿真模型输出仿真的效率指标和安全指标，与实际效率指标值和实际安全指标值进行对比，若单一指标，实际指标值落在仿真指标范围域的90％内，认为有效；若为多指标，则各个指标均同时落在其对应仿真指标90％范围域内为有效，否则，进行参数标定；S9.采用拉丁超方设计实验方法对模型参数进行进一步标定，将仿真参数集合按照其缺省的参数范围域进行随机采样；S10.将实际指标与通过拉丁超方设计实验得到的仿真指标范围域进行对比，若实际指标值落在仿真指标的90％范围域内，则认为目前的参数集合的范围域设定合理，如果并不落在范围域内，则重新进行拉丁超方设计随机抽样；S11.使用ANOVA检验或相关性检验，筛选出影响效率指标和安全指标的仿真参数；其中，nj为某仿真参数在j水平值下的样本个数，为该参数的在第j个水平值下的效率或安全指标均值，为效率或安全指标的均值，K为该参数的水平值个数，N为总的仿真次数；S13.针对选择的效率指标，选择相对误差值作为适应度函数，使用遗传算法进行参数的第一轮标定；S14.针对选择的安全指标，选择适应度函数为JS－divergence距离，其中DKL(P||Q)为P到Q的距离，p(x),q(x)分别为仿真冲突与实际冲突的后侵占时间PET分布概率密度函数；S15.定义约束条件一：仿真交通冲突的数量与实际交通冲突数量的差异性不能超过10％，并构造惩罚函数H1：H1:0.9-NCorNC-1.1其中NC为仿真冲突数量与实际冲突数量比值；S16.定义约束条件二：仿真得到的效率指标与实际的效率指标的差异性不能超过10％，并构造相应的惩罚函数H2：H2：0.9-EMorEM-1.1其中EM为仿真效率指标与实际效率指标的比值；S17.以步骤S15和S16定义的条件为约束，添加动态惩罚函数，定义优化目标函数，通过遗传算法进行标定：F＝JSD+(C*t)a(H1b+H2b)(5)其中：C，a,b是常数；S18.将标定后的模型与未标定的模型进行模型输出结果的对比，验证标定模型的优越性；S19.使用新的验证集，用标定的仿真模型的输出指标与实际指标进行对比，判断仿真模型的普适性。其中，步骤S13使用遗传算法进行参数的第一轮标定具体包括以下步骤：S131.随机选取初始种群，包含若干个体进行编码；S132.计算每个个体所对应的适应度函数，采用2人竞争算子法，选择出“优质”个体；S133：选择均匀交叉策略，对保留下来的优质个体进行“交叉”操作，产生新的个体；S134：选择均匀变异策略，对个体进行变异；迭代步骤S131-S134的过程，直到目标适应度函数收敛。其中，步骤S17通过遗传算法进行标定具体包括以下步骤：S171.随机选取初始种群，包含若干个体进行编码；S172.计算每个个体所对应的适应度函数，采用2人竞争算子法，选择出“优质”个体；S173.选择均匀交叉策略，对保留下来的优质个体进行“交叉”操作，产生新的个体；S174.选择均匀变异策略，对个体进行变异；重复步骤S171-S174，直至算法收敛，获得最优仿真参数集。其中，步骤S11和S13之间还包括步骤S12，所述步骤S12包括使用3Dcontour图对参数进行辅助筛选。其中，步骤S11的拒绝参数显著性的判断阈值为0.05。有益效果：本发明与现有技术相比，本发明提供的考虑交通冲突指标分布问题的交通仿真模型标定方法，区别于传统的基于效率指标的仿真模型标定方法，同时考虑了安全指标和效率指标，使仿真模型输出结果和真实情况更加接近。附图说明图1为本发明的流程图；图2为本发明的拉丁超方设计实验方法示意图；图3为本发明将实际指标与通过拉丁超方设计实验得到的仿真指标范围域进行对比示意图；图4为一个交叉路口基于无人机航拍与辅助视频分析的交通冲突提取示意图；图5为本实施例中标定后的分布示意图；图6为本实施例验证的分布比较示意图。具体实施方式下面结合附图对本发明作更进一步的说明。本发明提供一种考虑交通冲突指标分布多约束优化问题的交通仿真模型标定方法，同时考虑效率指标和安全指标，选择合适的交叉口效率指标作为交通仿真建模第一轮标定目标以及第二轮标定目标的约束条件，选择合适的交叉口安全指标作为交通仿真建模第二轮标定目标；通过计算仿真指标的标准差以及选择合适的置信区间对建立的仿真模型的最小所需仿真次数进行估计，有效减少参数校正的仿真试验次数；使用ANOVA检验或相关性检验，筛选出显著影响效率指标和安全指标的仿真参数；建立多约束的标定优化问题，提出使用遗传算法对模型参数进行标定的具体策略和流程；提出适应度参数选取方案；提出模型约束条件选择方案；通过对于交通冲突参数分布进行模型标定。参考图1，具体包括以下步骤：S1.收集交通仿真建模所有必要的数据，包括道路几何特征数据、交通控制数据、交通管理数据、交通流量数据和交通冲突指标；该道路几何特性数据包括CAD或商业地图数据，交通控制数据包括信号灯配时数据，交通管理数据包括标志标线、路权管理和专用车道，交通流量数据包括各个路口的交通流、车流比、转向比、车辆速度分布和车头间距分布，交通冲突指标包括避免碰撞时间TTC指标、后侵占时间PET指标和数量；S2.根据步骤S1收集到的数据，选择适合的交通仿真软件，比如VISSIM或Paramics，进行交通仿真的建模；S3.选择合适的交叉口效率指标作为第一轮标定目标以及第二轮标定目标的约束条件，该交叉口效率指标包括交叉口平均排队长度、控制延误、两点间车辆行程时间等；S4.选择合适的交叉口安全指标作为第二轮标定目标，可选用的安全指标包括交通冲突的TTC指标，PET指标等；S5.对建立的仿真模型的最小所需仿真次数进行估计，首先进行数次仿真(如50次)，提取仿真输出的效率指标和安全指标；S6.计算仿真效率指标的平均排队长度标准差；其中，x为平均排队长度，N为仿真次数。S7.选择置信区间计算该仿真模型所需的最小仿真次数，一般置信区间可选95％或90％；其中，C为1-置信度，t为对应置信度的分位点，s为标准差。S8.缺省参数的适用性检验(feasibilitytest)：在不调整缺省参数的情况下，运行仿真模型输出仿真的效率和安全指标，与实际的指标进行对比，如果仿真指标与实际指标的差异不大，则认为仿真模型缺省参数有效。如果差异过大，则需要进行参数标定。如果是单一指标，实际指标值落在仿真指标范围域的90％内，认为有效，若为多指标，则实际各个指标均同时落在其对应仿真指标90％范围域内为有效；S9.采用拉丁超方设计实验方法对模型参数进行进一步标定，将仿真参数集合按照其缺省的参数范围域进行随机采样，最终使得各仿真参数之间的相关性降到很低的水平，这样做的目的是使得参数集的取样尽可能覆盖整个参数空间；S10.将实际指标与通过拉丁超方设计实验得到的仿真指标范围域进行对比，如果实际指标值落在仿真指标的90％范围域内，则认为目前的参数集合的范围域设定合理。如果并不落在范围域内，则重新进行拉丁超方设计随机抽样；S11.使用ANOVA检验或相关性检验，筛选出显著影响效率指标和安全指标的仿真参数。拒绝参数显著性的判断阈值为0.05；其中，nj为某仿真参数在j水平值下的样本个数，等于在该具体参数值下的仿真次数；为该参数的在第j个水平值下的效率(安全)指标均值；为效率(安全)指标的均值；K为该参数的水平值个数，例如共10个不同的水平值，则k＝4；N为总的仿真次数。S12.使用3Dcontour图辅助筛选与显著因素存在共同显著效应(combinedeffect)但本身并不显著的参数，目的是检验根据先验信息判断可能存在共同效应的参数；S13.针对选择的效率指标，选择相对误差值作为适应度函数，使用遗传算法进行参数的第一轮标定，具体包括以下步骤：S131.随机选取初始种群，包含若干个体，该个体为与效率指标显著相关的参数值集合，进行编码；S132.计算每个个体所对应的适应度函数，采用2人竞争算子法，选择出“优质”个体；S133：选择均匀交叉策略，概率初始可设为0.5，根据拟合情况进行调整，对保留下来的优质个体进行“交叉”操作，产生新的个体；S134：选择均匀变异策略，概率可设为0.1，并根据拟合情况进行调整，对个体进行变异；迭代步骤S131-S134的过程，直到目标适应度函数收敛；S14.针对选择的安全指标，选择适应度函数为JS－divergence距离，其中DKL(P||Q)为P到Q的距离，p(x),q(x)分别为仿真冲突与实际冲突的后侵占时间PET分布概率密度函数。S15.定义约束条件一：仿真交通冲突的数量与实际交通冲突数量的差异性不能超过10％。并构造惩罚函数H1：H1:0.9-NCorNC-1.1其中NC为仿真冲突数量与实际冲突数量比值；S16.定义约束条件二，仿真得到的效率指标与实际的效率指标的差异性不能超过10％，并构造相应的惩罚函数H2：H2：0.9-EMorEM-1.1其中EM为仿真效率指标与实际效率指标的比值；S17.以步骤S15和S16定义的条件为约束，添加动态惩罚函数，定义优化目标函数。通过遗传算法进行标定。F＝JSD+(C*t)a(H1b+H2b)(5)其中：C，a,b是常数，取C＝0.5，a＝b＝2；具体包括以下步骤：S171.随机选取初始种群，包含若干个体，该个体为与安全指标显著相关的参数值集合，进行编码；S172.计算每个个体所对应的适应度函数，采用2人竞争算子法，选择出“优质”个体；S173.选择均匀交叉策略，概率初始可设为0.5，根据拟合情况进行调整，对保留下来的优质个体进行“交叉”操作，产生新的个体；S174.选择均匀变异策略，概率可设为0.1，并根据拟合情况进行调整，对个体进行变异；重复步骤S171-S174，直至算法收敛，获得最优仿真参数集；S18.将标定后的模型与为标定的模型进行模型输出结果的对比，验证标定模型的优越性；S19.使用全新的验证集，用标定的仿真模型的输出指标与实际指标进行对比，判断仿真模型的普适性。下面结合具体实施例对本发明作更进一步的说明。本实施例选取了上海市奉贤区望园-南奉路，通过无人机航拍的方法对每一个路口进行了4个小时的拍摄，将拍摄的视频结合人工与视频分析的方法，提取观测冲突数据，图4示例了一个交叉路口基于无人机航拍与辅助视频分析的交通冲突提取。交通仿真软件使用VISSIM，交通仿真冲突通过SSAM软件提取仿真输出轨迹获得，其中3个小时的数据用于标定，一个小时的数据用于验证，标定后的参数详见表1。表1标定后的参数标定参数缺省值标定值平均的停车间距21.6安全距离的附加部分22.2安全距离的倍数部分32.8可视距离100m89m前方最小空挡0.50.46TTC阈值1.51.6安全距离折减系数0.60.68参考图5和图6可知，分布基本吻合，均值和方差无显著性差异，k-s检验表示无显著差异。当前第1页1 2 3