一种恶劣霾天气的城市区域交通紧急预案方法
【专利摘要】本发明公开一种恶劣霾天气的城市区域交通紧急预案方法。针对近期出现的雾霾天气，本发明建立了由数据库模块、中观仿真器模块、离线优化平台模块、择优预案模块组成的紧急预案系统。系统模拟了城市区域机动车污染物排放情况和污染物浓度扩散情况。基于这些情况，预案系统通过计算和优化，获得了合理有效的紧急预案方法，即限制机动车进入城市区域的比例和区域内信号灯配时的新方案。在恶劣天气下，本发明方法保证了城市区域的天气不会由于机动车污染物排放而进一步恶化，也避免了污染物威胁人类的身体健康。
【专利说明】-种恶劣霾天气的城市区域交通紧急预案方法
【技术领域】
[0001] 本发明属于交通环境工程【技术领域】，具体涉及一种恶劣霾天气的城市区域交通紧 急预案方法。
【背景技术】
[0002] 随着我国国民经济和汽车工业的不断发展，在北京、上海等大城市的PM2. 5污染 中，机动车排放的分担率占22 %?34 %。解决上述交通环境问题，首先需调研当地机动车 的交通流数据和交通能耗情况，其次定量（或定性）的估计该地区的污染物排放因子和相 应的道路区域分布情况，然后在此基础上通过实施各种交通管理与控制措施，比如信号灯 控制、限制机动车数量、变化禁左车道数、变限速等，分析其对交通状况的影响，最后通过微 观交通仿真软件选择一种较优的方案并予以实施，从而达到改善城市交通环境的目的。
[0003] 目前，我国关于交通环境控制的研究已初具成果，但仍有一些不足。比如，在实现 控制时没有考虑机动车的动态限行；在考虑环境污染指标时，没有将污染物浓度纳入交通 管理控制的指标；在择优预案时，没有对微观交通仿真软件进行参数校正，直接将其应用则 不符合当地实际交通流状况。
[0004] 本发明综合考虑以上不足，同时结合前人的研究成果，提出了一种恶劣霾天气的 城市区域交通紧急预案方法。
【发明内容】
[0005] 本发明旨在针对现有技术的不足和近年来出现的霾天气，提供一种恶劣霾天气的 城市区域交通紧急预案方法，通过提出合理的预案方法，限制部分车辆通行以及重新调整 限制区域的信号灯配时，保证了当地天气不会由于尾气排放而进一步恶化，避免有害物质 侵袭行人的身体健康。
[0006] 本发明方法应用于以下预案系统，该预案系统包括数据库模块、中观仿真器模块、 离线优化平台模块、择优预案模块；数据库模块将采集到的数据传送至中观仿真器模块、离 线优化平台模块、择优预案模块；中观仿真器模块通过交通流模型、机动车污染物排放模型 和污染物扩散模型，模拟得到机动车总旅行时间、实地区域污染物浓度两个数据，并将上 述两个数据传送至离线优化平台模块；离线优化平台模块根据实地区域污染物浓度确立应 急预案的集合，得到最优预案集，并将最优预案集传送至择优预案模块；择优预案模块通过 微观交通仿真软件平台选取实地区域总旅行时间最少的预案作为最优预案，并且在实地区 域实施。
[0007] 本发明方法包括以下步骤：
[0008] 步骤（1).对实地区域采集环境参数、机动车流量、信号灯配时、道路几何参数四 个数据，并建立数据库模块；
[0009] 1. 1)对实地区域采集环境参数；其中环境参数包括实地区域污染物浓度背景值 (即预测的PM2. 5浓度背景值）、污染物排放因子、实地区域路网的近似箱体体积；
[0010] 1. 1. 1)实地预测未来一天的实地区域污染物浓度背景值PQ :
[0011] 通过某浓度监测网站查询进行统计距离当天η天之内的所有空气质量指数心? Κ ；
[0012] 采用加权移动平均法对未来1天的预测实地空气质量（AQI)指数进行预测，计算 如下：
[0013] = w〇 · A0+W! · A2+··· +wn_! · An_!+wn · An (1)
[0014] 其中A_i为未来1天的空气质量指数，wn为距离当天n天的相应权重。
[0015] 可通过查询现有资料将未来1天的空气质量指数A_i转化为实地区域污染物浓度 背景值P〇 ;
[0016] 1. 1. 2)通过机动车车载测试对实地区域采集机动车污染物排放因子：
[0017] 所述的机动车污染物排放因子是指实地区域内，机动车在单位时间内以不同 工况行驶排放污染物的量，单位为g/s，其中不同工况为机动车的行驶状态，具体是怠速 (idling)、加速（accelerating)、减速（decelerating)以及巡航（cruise);
[0018] 所述的机动车车载测试内设置有专业汽车尾气排放测量仪、GPS定位系统；通过 GPS定位系统测得机动车逐秒的地理位置，对地理位置进行差分和二次差分计算即可得到 机动车逐秒的速度与加速度，根据速度和加速度的数值将行驶时间归类到不同工况中；通 过专业汽车尾气排放测量仪测得逐秒的污染物的排放量，相应地，也将污染物排放量也以 不同工况进行统计，污染物排放因子计算公式如下：
【权利要求】
1. 一种恶劣霾天气的城市区域交通紧急预案方法，应用于以下预案系统，该预案系统 包括数据库模块、中观仿真器模块、离线优化平台模块、择优预案模块；数据库模块将采集 到的数据传送至中观仿真器模块、离线优化平台模块、择优预案模块；中观仿真器模块通过 交通流模型、机动车污染物排放模型和污染物扩散模型，模拟得到机动车总旅行时间、实地 区域污染物浓度两个数据，并将上述两个数据传送至离线优化平台模块；离线优化平台模 块根据实地区域污染物浓度确立应急预案的集合，得到最优预案集，并将最优预案集传送 至择优预案模块；择优预案模块通过微观交通仿真软件平台选取区域总旅行时间最少的预 案作为最优预案，并且在实地区域实施；其特征在于该方法包括以下步骤： 步骤（1).对实地区域采集环境参数、机动车流量、信号灯配时、道路几何参数四个数 据，并建立数据库模块： 1. 1)对实地区域采集环境参数；其中环境参数包括实地区域污染物浓度背景值、污染 物排放因子、实地区域路网的近似箱体体积； 1. 1. 1)实地预测未来一天的实地区域污染物浓度背景值Ρο; 采用加权移动平均法对未来1天的预测实地空气质量（AQI)指数进行预测，见公式 (1)，得到未来1天的空气质量指数A_i ;然后将未来1天的空气质量指数A_i转化为实地区 域污染物浓度背景值Ρ〇; Α-! = W0 · A0+W! · Α2 + · ··+￥"_! · An_!+wn · Αη (1) 其中A_i为未来1天的空气质量指数，wn为距离当天η天的相应权重，Αη为距离当天η 天的空气质量指数； 1. 1. 2)通过机动车车载测试对实地区域采集机动车污染物排放因子； 所述的机动车污染物排放因子是指实地区域内，机动车在单位时间内以不同工况行驶 排放污染物的量，单位为g/s ;其中不同工况为机动车的行驶状态，具体是怠速、加速、减速 以及巡航； 污染物排放因子计算公式如下：
(2) 其中h为工况j时机动车污染物排放因子，4为工况j时污染物排放量总和，L为工 况j时时间统计； 1. 1. 3)通过电子地图计算实际实地区域的近似箱体体积Λ V，如下： AV = S · Η (3) 其中箱体底面积S为覆盖实际实地区域的最小面积，Η为高度； 1. 2)对实地区域采集机动车流量并进行统计； 所述的机动车流量包括所有交叉口进口路段的车流量Q，以及左、直、右转弯的车流量 转弯比例β ; 1. 3)通过实地区域采集获取信号灯配时； 所述的信号灯配时包括所有交叉口的信号灯周期Υ，各个相位的开始b时刻和持续时 间d ; 1. 4)通过实地区域采集获取路网的道路几何参数； 所述的道路几何参数为相邻交叉口路段的长度、车道数、交叉口渠化、车道连接以及车 道限制； 步骤（2).步骤（1)采集得到的数据通过中观仿真器模块，得到实地区域机动车总旅行 时间、实地区域污染物浓度： 所述的中观仿真器模块包括交通流模型、机动车污染物排放模型、污染物扩散模型； 2. 1)通过现有的交通流模型建立交通流仿真器模拟机动车车流的运转；交通流仿真 器根据步骤（1)采集得到的机动车流量、信号灯配时、道路几何参数三个数据进行仿真，得 到实地区域机动车总旅行时间J、各个工况时间统计T ; 2. 1. 1)通过现有的交通流模型建立交通流仿真器；交通流仿真器仿真过程中记录第k 个采样周期的路段剩余容量c(k)、车辆排队长度ω〇〇 ;仿真结束时，将路段车辆总容量C 减去第k个采样周期的路段剩余容量c (k)，差值即为在该路网上第k个采样周期的机动车 总数目C-c(k);然后第k个采样周期的机动车总数目C-c(k)乘以采样周期时长t s即可得 到第k个采样周期实地区域的机动车总旅行时间J (k)，如下式： J (k) = ts · [C-c (k) ] (6)； 2. 1. 2)各个工况时间统计分为机动车怠速、加速、减速、巡航四个工况的时间统计，计 算公式如下： Ti (k) = ts · ω (k) (7) Ta (k) = ts · a (k) (8) Td(k) = ts · d(k) (9) Tc (k) = J (k) -?\ (k) -Ta (k) -Td (k) (10) 通过在第k个采样周期内统计车辆排队的长度ω (k)，乘以采样周期时长ts即可获得 工况为怠速的时间?\ ;通过在第k个采样周期内统计排队车辆驶离停车线的数目a(k)，乘 以采样周期时长ts即可获得加速工况的时间T a ;通过在第k个采样周期内统计机动车加 入车辆排队队伍的数目d(k)，乘以采样周期时长ts即可获得减速工况的时间T d ;通过在第 k个采样周期内机动车总旅行时间J减去怠速、加速、减速的时间即可获得在巡航工况下的 时间T。； 2. 2)将各个工况时间统计作为机动车污染物排放模型的输入，使其输出机动车排放的 污染物量； 所述的机动车污染物排放模型将在数据库模块输出的不同行驶工况的机动车污染物 排放因子与各行驶工况的时间统计线性叠加，从而输出机动车排放的污染物量，如下式： E (k) = ξ 山(k) + ξ 2Ta (k) + ξ 3Td (k) + ξ 4TC (k) (11) 其中E(k)为在第k个采样周期时间段内污染物排放量，ξρ ξ2, ξ3, ξ4分别为怠速、 加速、减速、巡航上述四种不同行驶工况的机动车污染物排放因子； 2. 3)将上述步骤2. 2得到的机动车排放的污染物量E (k)，步骤1. 1. 1得到的预测的污 染物浓度背景值Ρ〇,以及步骤1. 1. 3得到的近似箱体体积Λ V输入到污染物扩散模型，即可 输出实地区域污染物浓度ps，其计算公式如下： 步骤（3).离线优化平台模块建立： (12)； 所述的离线优化平台模块包括两个子模块：放行比例优化子模块和信号灯配时优化子 模块；模块的输入为数据库模块的所有信息，输出是最优预案集，预案集合中每一个预案包 括已优化的放行比例和整个实地路网信号灯配时方案； 3. 1)放行比例优化子模块建立： 所述的放行比例是指为减小实地区域的污染物浓度，放行的机动车数目占实地区域机 动车总数量的比例；具体算法如下： a) 首先确定放行比例的上下界：下限为low = 0,上限为high = 1 ; b) 令放行比例为
判断low < high-δ，若是则直接获取最优放行比 例mid，算法结束；若否则跳转步骤b)，然后将数据库模块的机动车流量乘以mid，其余的数 据不变化，然后将环境参数、新机动车流量、信号灯配时、道路几何参数输入到中观仿真器 模块中，获取污染物浓度Ps;其中S为容许误差； c) 设定plimit为日均PM2. 5浓度的上界阈值，判断ps < plimit- ε，若是则将mid赋值给 low，返回到步骤b);若否则执行步骤d);其中ε为容许误差； d) 判断ps>plimit，若是则将mid赋值给high，返回到步骤b);若否则直接得到最优放 行比例mid，算法结束； 3. 2)信号灯配时优化子模块建立： 信号灯配时优化子模块分为三个部分：实地区域主干道扫描算法，主干道双向协调控 制优化算法和非主干道相位差优化算法；子模块的外部以数据库模块和放行比例优化子模 块作为输入，而将已优化的放行比例mid和整个路网信号灯配时方案作为最优预案输出； 3. 2. 1)实地区域利用主干道扫描算法建立主干道路径： 首先通过数据库模块中机动车流量数据确定车流量最大的交叉口作为关键路口；然后 在关键路口的几个出口方向，依据车流量最大和次大的出口方向作为主干道的伸展，并且 延伸到下一个交叉口；然后在下一个交叉口根据机动车流量数据确定车流量最大的出口方 向进行伸展，直到接触到实地区域的边界，得到主干道路径； 3. 2. 2)利用主干道双向协调控制优化算法优化主干道： 根据数据库模块的机动车流量数据、步骤3. 2. 1得到的主干道路径，使用现有的主干 道双向协调控制优化算法将该主干道进行信号灯配时优化，得到新的信号灯周期Y'以及 主干道路径所有交叉口中各相位的新开始时间b'、新持续时间d'，即主干道的新信号灯 配时；具体算法如下： a) 使用模糊控制确定主干道公共信号灯周期，将交叉口饱和度X与周期的增量ΛΥ 作为模糊变量；其次确定论域和隶属度函数；最后去模糊化后得到新的周期预案f = Υ+ΛΥ，且Y' e [Ymin，Ymax];其中Y'为新的信号灯周期，Y原先的信号灯周期，饱和度X计 算见公式（13)，fl 〇V为r'对应的实际车流量，为相位r的有效绿灯时间，s为饱和流 率，饱和度X为所有相位的饱和度最大值：
(13) b) 实际交通流流量确定主干道各相位新持续时间d'，持续时间的计算如式14所示： (14) 其中fl〇K为相位r对应的实际车流量，cT ^为相位r的持续时间； c)确定相位的持续时间之后，通过现有的智能算法确定主干道的相位的起始时间； 3. 2. 3)非主干道相位差优化算法建立： 为依次优化对非主干道的路口相位差，得到非主干交叉口的新信号灯配时，其算法流 程如下： a) 将步骤3. 2. 1得到的主干道路径上所有交叉口作为元素放入集合U中； b) 寻找实地区域中与集合U每一个元素相邻的交叉口，放入集合Z中； c) 在集合U中，每一个元素的信号灯配时优化已经完成；而在集合Z中，每一个元素的 信号灯配时优化并未完成；由于集合Z的某一个元素 z在集合U中至少有一个元素 u与其 相邻，所以使用单向协调相位差方法即可优化集合Z中元素 z，单向协调相位差方法计算公 式如下： Oz，u = mod(lz，u/vz，u,Y'） （17) 其中〇z，u为交叉口 z与交叉口 u的相位差，lz，u为交叉口 z与交叉口 u的距离，vz，u为 交叉口 z与交叉口 u的机动车车流速度，Y'为主干道的新信号灯周期； d) 交换集合U、Z中元素，交换之后将集合Z清除； e) 判断新集合U是否为空，若是则表明实地区域路网所有的交叉口均已优化完毕，结 束；若否则跳到步骤b)继续优化； 步骤3. 2. 4)由上述步骤得到的放行比例mid、主干道的新信号灯配时Y'以及非主干 道的新信号灯配时获得若干个可行预案，构成最优预案集； 步骤（4).择优预案模块建立： 在步骤（3)输出的最优预案集中，使用微观交通仿真软件，找出评价指标最高的预案 作为最终的最优预案；其中评价指标是指所有车辆在实地区域中的旅行时间的总和，称为 总旅行时间，符号为J，单位为h ;模块建立主要分为以下二步： 4. 1)微观交通仿真软件平台上搭建交通仿真路网： 根据数据库模块内的道路几何参数、信号灯配时、机动车流量数据，在微观交通仿真软 件上搭建交通仿真路网； 4. 2)微观交通仿真软件平台上仿真每一个预案，进行筛选； 根据最优预案集中每一预案的机动车流量放行比例mid以及实地区域新信号灯配时 在微观交通仿真软件平台上进行仿真； 所述的实地区域新信号灯配时包括上述步骤得到的主干道新信号灯配时、非主干道新 信号灯配时； 4. 3)仿真完毕后，选择最优预案集中总旅行时间最小J的预案作为最优预案； 步骤5).实地区域实施预案： 将获得的最优预案在实地区域路网中进行实施。
【文档编号】G08G1/00GK104112357SQ201410317894
【公开日】2014年10月22日 申请日期:2014年7月4日 优先权日:2014年7月4日
【发明者】王绍楠, 王慧 申请人:浙江大学