[发明专利]一种快速道路孤立瓶颈路段的可变限速控制策略优化方法在审

专利信息
申请号: 201610316223.8 申请日: 2016-05-13
公开(公告)号: CN105931459A 公开(公告)日: 2016-09-07
发明(设计)人: 李志斌;刘攀;王炜;徐铖铖 申请(专利权)人: 东南大学
主分类号: G08G1/01 分类号: G08G1/01;G08G1/09
代理公司: 江苏永衡昭辉律师事务所 32250 代理人: 王斌
地址: 210096*** 国省代码: 江苏;32
权利要求书: 查看更多 说明书: 查看更多
摘要: 一种快速道路孤立瓶颈路段的可变限速控制策略优化方法。在快速道路孤立瓶颈路段上下游设置检测器和可变限速指示牌,采集实时交通流数据和发布可变限速控制策略,基于遗传算法挖掘控制效果最优的可变限速控制参数值,依据事故预测模型计算各路段内实时追尾事故风险,当路段事故风险超过启动阈值时启动可变限速控制,逐步调整限速值至目标值,同时用相邻路段最大限速差值控制限速值的空间波动幅度;当事故风险小于阈值时,可变限速值逐步恢复到默认值。本发明弥补了可变限速控制中核心参数取值的随意性,使可变限速控制的限速值逐步连续变化,降低了限速值在时空上的大幅波动,缓和驾驶员的速度变化,有效降低限速值突变导致的交通流波动与紊乱。
搜索关键词: 一种 快速 道路 孤立 瓶颈 路段 可变 限速 控制 策略 优化 方法
【主权项】:
一种快速道路孤立瓶颈路段的可变限速控制策略优化方法,其特征是包括以下步骤:1)确定孤立瓶颈路段范围及交通流检测器设置合理间距,在快速道路路段内合理位置配套设置交通流检测器和可变限速控制指示牌。通过交通流检测器周期性检测孤立瓶颈区域及其上下游交通流数据;2)基于事故预测模型及瓶颈处交通流检测器检测到的交通流数据,计算实时追尾事故风险,具体步骤包括:201)计算追尾事故风险指数RCRI,计算公式如下:<mrow><mi>R</mi><mi>C</mi><mi>R</mi><mi>I</mi><mo>=</mo><mfrac><mrow><mo>&lsqb;</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>U</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>T</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>V</mi><mo>&OverBar;</mo></mover><mi>D</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>T</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>&CenterDot;</mo><msub><mover><mi>O</mi><mo>&OverBar;</mo></mover><mi>U</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>T</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><msub><mover><mi>O</mi><mo>&OverBar;</mo></mover><mi>U</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>&Delta;</mi><mi>T</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow></mfrac></mrow>其中,为时间Δt内上游检测器位置平均速度,为时间Δt内下游检测器位置平均速度,为时间Δt内上游检测器位置平均占有率,分别为第j个时间Δt内上下游检测器位置车道m平均速度和上游检测器位置车道m的平均占有率;J为一个时间段内集计交通流数据个数(J=ΔT/Δt,Δt=30s);M为路段断面车道数;202)计算瓶颈所处路段i内当前时刻t的实时追尾事故风险Ri(t),计算公式如下:<mrow><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mrow><mo>(</mo><mi>Y</mi><mo>=</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mn>3.095</mn><mo>+</mo><mn>0.191</mn><mo>&CenterDot;</mo><mi>R</mi><mi>C</mi><mi>R</mi><mi>I</mi><mo>+</mo><mn>0.178</mn><mo>&CenterDot;</mo><mi>&sigma;</mi><mo>(</mo><msub><mi>O</mi><mi>U</mi></msub><mo>)</mo><mo>+</mo><mn>0.172</mn><mo>&CenterDot;</mo><mi>&sigma;</mi><mo>(</mo><msub><mi>O</mi><mi>D</mi></msub><mo>)</mo><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mn>3.095</mn><mo>+</mo><mn>0.191</mn><mo>&CenterDot;</mo><mi>R</mi><mi>C</mi><mi>R</mi><mi>I</mi><mo>+</mo><mn>0.178</mn><mo>&CenterDot;</mo><mi>&sigma;</mi><mo>(</mo><msub><mi>O</mi><mi>U</mi></msub><mo>)</mo><mo>+</mo><mn>0.172</mn><mo>&CenterDot;</mo><mi>&sigma;</mi><mo>(</mo><msub><mi>O</mi><mi>D</mi></msub><mo>)</mo><mo>)</mo></mrow></mrow></mfrac></mrow>其中,P(Y=1)为追尾事故发生概率;σ(OU)为上游检测器占有率标准差,σ(OD)为下游检测器占有率标准差,3)判断是否启动可变限速控制策略并计算当前时刻各位置可变限速控制指示牌显示的限速值,具体步骤包括:301)比较步骤2)中计算得到的当前时刻t瓶颈所处路段i的事故风险值Ri(t)和可变限速控制启动阈值RTH,若该路段事故风险超过启动阈值则启动瓶颈处上游可变限速标志的可变限速控制,逐步将限速值调整至目标安全限速值;反之,若该路段事故风险小于启动阈值,则可变限速值逐步恢复到默认值。具体公式如下:<mrow><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mi>d</mi><mi>e</mi><mi>f</mi><mi>a</mi><mi>u</mi><mi>l</mi><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&Delta;V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&gt;</mo><msub><mi>R</mi><mrow><mi>T</mi><mi>H</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>min</mi><mo>{</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mi>d</mi><mi>e</mi><mi>f</mi><mi>a</mi><mi>u</mi><mi>l</mi><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>}</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>a</mi><mi>l</mi><mi>l</mi><mi> </mi><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi> </mi><mi>c</mi><mi>a</mi><mi>s</mi><mi>e</mi><mi>s</mi></mrow></mtd></mtr></mtable></mfenced></mrow>302)依据下式计算可变限速控制启动后位置xi的可变限速标志在t时刻的变化步长值:<mrow><msub><mi>&Delta;V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><mi>&Delta;</mi><mi>V</mi><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>T</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi></mrow><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow><mo>-</mo><mi>&Delta;</mi><mi>V</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow><mo>-</mo><mi>&Delta;</mi><mi>V</mi><mo>&le;</mo><msub><mi>T</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi></mrow><mo>)</mo></mrow><mo>&le;</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow><mo>+</mo><mi>&Delta;</mi><mi>V</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>&Delta;</mi><mi>V</mi><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi></mrow></mtd><mtd><mrow><msub><mi>T</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi></mrow><mo>)</mo></mrow><mo>&gt;</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi></mrow><mo>)</mo></mrow><mo>+</mo><mi>&Delta;</mi><mi>V</mi></mrow></mtd></mtr></mtable></mfenced></mrow>其中,ΔV为路段限速值变化步长;TVSL为目标安全限速值;303)将步骤302)中计算得到的位置xi的可变限速标志在t时刻的变化步长值带入可变限速控制限速值的计算公式,计算得到路段i内位置xi的可变限速标志在t时刻的可变限速控制状态计算其在t+Δt时刻的限速值,计算公式如下:<mrow><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mi>d</mi><mi>e</mi><mi>f</mi><mi>a</mi><mi>u</mi><mi>l</mi><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>&gt;</mo><msub><mi>R</mi><mrow><mi>T</mi><mi>H</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mi>min</mi><mo>{</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><mi>d</mi><mi>e</mi><mi>f</mi><mi>a</mi><mi>u</mi><mi>l</mi><mi>t</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&Delta;V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>}</mo><mo>,</mo><mi>a</mi><mi>l</mi><mi>l</mi><mi> </mi><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi> </mi><mi>c</mi><mi>a</mi><mi>s</mi><mi>e</mi><mi>s</mi></mrow></mtd></mtr></mtable></mfenced></mrow>其中,RTH为可变限速控制启动阈值;Ri(t)为当前时刻t瓶颈所处路段i的事故风险值;VSL(xi,t+Δt)为路段i内位置xi的可变限速标志在t+Δt时刻的限速值;VSL(default)为路段默认限速值;ΔVSL(xi,t)为位置xi的可变限速标志在t时刻的变化步长;4)在步骤3)中确定了各路段内限速值后,计算位置xi的可变限速标志与相邻上游可变限速标志的限速值差,若差值大于预设的相邻路段最大允许限速值差,则启动相邻上游可变限速标志的可变限速控制,依据下式计算限速值改变步长并更新发布相邻上游可变限速标志的限速值:<mrow><msub><mi>&Delta;V</mi><mrow><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>-</mo><msup><mi>&Delta;V</mi><mo>&prime;</mo></msup><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mi>f</mi><mi> </mi><msub><mi>T</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>&le;</mo><msub><mi>T</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>,</mo><mi>t</mi><mo>+</mo><mi>&Delta;</mi><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msup><mi>&Delta;V</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>a</mi><mi>l</mi><mi>l</mi><mi> </mi><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi> </mi><mi>c</mi><mi>a</mi><mi>s</mi><mi>e</mi><mi>s</mi></mrow></mtd></mtr></mtable></mfenced></mrow>其中,ΔV′为相邻路段最大限速值差;5)确定可变限速控制策略中的四个关键参数的取值范围及遗传算法中的优化目标函数,具体步骤包括:501)可变限速控制策略包含的四个核心控制参数分别为启动可变限速的事故风险阈值RTH、目标安全限速值RTH、限速值改变步长ΔV和相邻路段最大限速值差ΔV′,启动可变限速控制的事故风险阈值的取值范围为最大事故风险的10%‑30%,以5%的步长递进;目标安全限速值的取值范围从拥堵中交通流速度25mph到50mph,以5mph的步长递进;限速值改变步长从每2分钟改变5mph到每30秒改变15mph;相邻路段最大限速值差为从5mph到25mph;502)遗传算法中的优化目标函数的计算公式如下:<mrow><mi>F</mi><mi>i</mi><mi>t</mi><mi>n</mi><mi>e</mi><mi>s</mi><mi>s</mi><mo>=</mo><mo>-</mo><mo>&lsqb;</mo><mi>&alpha;</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mrow><msub><mi>R</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>N</mi><mi>o</mi></mrow></msub></mrow><msub><mi>R</mi><mrow><mi>N</mi><mi>o</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mrow><msub><mi>T</mi><mrow><mi>V</mi><mi>S</mi><mi>L</mi></mrow></msub><mo>-</mo><msub><mi>T</mi><mrow><mi>N</mi><mi>o</mi></mrow></msub></mrow><msub><mi>T</mi><mrow><mi>N</mi><mi>o</mi></mrow></msub></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>其中,Fitness为适应度;a为权重系数,其目的为在事故风险下降和出行时间增加之间进行权衡,缺省取值为0.5;RVSL和RNO分别为可变限速控制下和无控制下事故风险;TVSL和TNO分别为可变限速控制下和无控制下总通行时间;6)基于遗传算法和交通流仿真模型对步骤5)中的四个关键参数在其取值范围内对其取值进行优化;在单目标遗传算法中四个关键参数的取值如下:种群大小M为20,、最大代数T为100、交叉概率Pc为0.9、变异概率Pm为0.1;初始化设置进化代数计数器为t=0,设置最大进化代数T,设置遗传算法中交叉概率Pc和变异概率Pm;随机生成M个个体作为初始群体P(0),第m个个体中包含各可变限速控制参数取值信息编码m(x1,x2,x3,x4);将群体P(t)中的个体参数值输入交通流仿真模型,依据仿真输出结果计算每个个体的适应度;从群体中选择优胜的个体,个体被选择的概率为个体适应度在整个群体适应度综合中所占比例,被选择后的个体随机组成交配对;以一定概率在个体串中随机设定一个交叉点,实行交叉时该点前或后的两个个体部分结构进行互换,并生成两个新个体;按照变异概率选取群体中部分个体改动其部分参数值;群体P(t)经过选择、交叉、变异运算之后,得到下一代群体P(t+1);当迭代次数达到最大代数T时,遗传算法终止;否则,重复初始化以后的各步骤直至终止条件满足;将具有最大适应度的个体对应的最优解(x1,x2,x3,x4)翻译成可变限速控制的四个核心参数的取值;7)将步骤6)中得到的事故风险阈值RTH、目标安全限速值RTH、限速值改变步长ΔV和相邻路段最大限速值差ΔV′的优化取值带入步骤3)和步骤4)中作为各关键参数的推荐取值,然后依据步骤1)到步骤4)的控制策略对快速道路鼓励瓶颈路段进行可变限速控制。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于东南大学,未经东南大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201610316223.8/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top