陳曉龍,胡志坤,雷 霆,谷 豐,劉光勇,王文祥
(1.中南大學 物理與電子學院,湖南 長沙410083;2.株洲南車時代電氣股份有限公司 智能交通事業(yè)部,湖南 株洲412001)
交通堵塞嚴重已成為阻礙城市發(fā)展的重要因素。對此,設計一套合理的區(qū)域協(xié)調(diào)控制方案是解決道路交通堵塞、優(yōu)化交通分布最主要的方式[1]。
目前,國外較為成功且應用廣泛的區(qū)域控制系統(tǒng)有SCATS、TRANSYT、SCOOT。其中,SCATS系統(tǒng)中區(qū)域協(xié)調(diào)控制方案是一種選擇式的自適應方案,以綜合流量和飽和度為依據(jù)對信號周期、綠信比、相位差分別進行獨立地優(yōu)選,然而該方案未使用交通模型,因此方案的確定具有一定的主觀性[2,3]。TRANSYT系統(tǒng)中區(qū)域協(xié)調(diào)控制方案是一種脫機定時控制方案,以延誤時間和停車次數(shù)為評價指標對綠信比和相位差進行優(yōu)化,但是這種方案計算量太大、需要大量的網(wǎng)絡幾何尺寸和交通流數(shù)據(jù)[4]。SCOOT系統(tǒng)中區(qū)域協(xié)調(diào)控制方案是一種在線配時優(yōu)化方案,以實時測量得到的交通流數(shù)據(jù)為依據(jù)對周期、相位差和綠信比進行頻繁且適量的調(diào)整。同樣這種方案也需要大量的路網(wǎng)幾何尺寸和交通流數(shù)據(jù),并且小步長地調(diào)整不足以響應每個周期的交通需求[5]。更關鍵的一點,我國的路網(wǎng)復雜、交通流獨特,因此國外的方案很難直接適用于我國的應用環(huán)境。而國內(nèi)文獻 [6,7]分別提出基于分層并行災變粒子群算法和雙層規(guī)劃模型的區(qū)域協(xié)調(diào)控制方案,都以平均延誤時間或平均排隊長度等為評價指標對周期、相位差和綠信比進行優(yōu)化,但是這些方案都呈現(xiàn)出算法復雜、收斂速度緩慢的缺點。
針對以上論述的問題,為了降低算法的復雜度、減少系統(tǒng)所需的輸入量,本文提出了一套干支線混合區(qū)域協(xié)調(diào)方案,以子區(qū)為基礎,動態(tài)地對子區(qū)進行合并和拆分操作,重點對子區(qū)內(nèi)部干線進行綠波協(xié)調(diào)控制以及對所有交叉口進行綠信比優(yōu)化,確定出子區(qū)的最佳周期、干線上所有交叉口的最佳相位差和所有交叉口下一周期最佳綠信比方案,解決了城市某些干支線混合區(qū)域的交通信號同時需要進行區(qū)域協(xié)調(diào)、干線控制和單個交叉口綠信比優(yōu)化的問題。通過實際數(shù)據(jù)驗證了本方案具有良好的控制效果。
由于路網(wǎng)交通流的特性總是以區(qū)域為單位,即一個區(qū)域內(nèi)相鄰交叉口的交通流特性會在某一個時段內(nèi)具有相似性[8]。因此,對干支線混合的城市道路進行區(qū)域協(xié)調(diào)控制能有效緩解堵塞、降低車輛行駛延誤、減少紅燈停車次數(shù)、縮短車輛在路網(wǎng)內(nèi)的行駛時間、提高路網(wǎng)整體通行能力。以某個包含干線、支線的混合區(qū)域為例,如圖1所示。
圖1 支干道混合區(qū)域
設干線和支線的交叉口個數(shù)總計為n;交叉口序號為i(i=1,2…n);各交叉口相位數(shù)為τi,相位序號為j(j=1,2…τi);根據(jù)路況測得所有交叉口各相位關鍵進口道的飽和流量為qmaxi,j,為常量;系統(tǒng)周期損失時間為l,一般取l=10s[9];車輛行駛速度為v。
根據(jù)檢測器測得所有交叉口中各相位當前飽和度和其關鍵進口道的當前流量分別為di,j、qi,j,并計算出所有交叉口的飽和度
區(qū)域協(xié)調(diào)的具體步驟如下:
步驟1 劃分子區(qū)
為了維持道路的完整性,將圖1中的某條干線所有交叉口劃成一個子區(qū),進而將整個區(qū)域劃分成m個子區(qū),如圖2所示?;谠撟訁^(qū)劃分原則,沿干線方向把整個區(qū)域劃分成3個子區(qū)A,B,C。
步驟2 計算各子區(qū)的初始周期
在子區(qū)的內(nèi)部,首先按單個交叉口信號控制的周期配時方法計算各交叉口的信號周期,然后從中選出最大的信號周期作為該子區(qū)的公共初始周期,對應的最大信號周期的交叉路口為關鍵交叉口[10]。單個交叉口初始周期計算如下式
圖2 劃分子區(qū)
式中:c——該子區(qū)下一周期采用的初始周期;ci——交叉口i計算的理論周期;ysumi——交叉口i總流率比;yi,j——交叉口i相位j流率比。
步驟3 子區(qū)的合并和拆分
子區(qū)的合并和拆分主要基于周期時長、流量這兩個基本原則[11]。由于子區(qū)的流量存在不均勻性,為了簡化模型,只用周期原則作為本文子區(qū)合并與劃分的唯一原則:設有兩個相鄰的子區(qū)a、b,周期分別為ca、cb。當時,ξ為一常量,一般取ξ<10,將子區(qū)a,b合并成新的子區(qū)ab,其周期為max(ca,cb);當>ξ時,如果子區(qū)a、b已經(jīng)合并,則將合并的子區(qū)ab重新拆分成兩個子區(qū)a、b,拆分后的子區(qū)a、b采取各自的周期ca、cb。
在步驟2的基礎上,根據(jù)上述周期原則對子區(qū)進行拆分和合并操作。
步驟4 計算子區(qū)的最佳周期和各交叉口的最佳相位差
在步驟3的基礎上,以劃分好的子區(qū)為基礎,對子區(qū)內(nèi)部干線的交叉口進行綠波協(xié)調(diào)控制得到子區(qū)的最佳周期和干線上所有交叉口的最佳相位差。
步驟5 確定子區(qū)內(nèi)各交叉口下一周期的綠信比方案
對子區(qū)內(nèi)部各個交叉口獨立地進行綠信比優(yōu)化確定出下個周期各個交叉口的最優(yōu)綠信比方案。
綠波協(xié)調(diào)控制是將某條干線上多個交叉口以一定方式聯(lián)結起來作為研究對象,對各個交叉口進行相位差設置,使盡可能多的車流沿這條干線行駛中以規(guī)定車輛行駛速度v行駛能連續(xù)得到一個接一個的綠燈放行時間,暢通無阻地通過沿途所有交叉口。
在理想狀態(tài)下,干線上相鄰交叉口的間距相等時,車輛通過各交叉口的時間相等,那么通過設置相同的相位差使車輛能夠按照規(guī)定行程速度連續(xù)通過所有交叉口。然而在實際情況中,干線上相鄰交叉口的間距大小相異,通過尋找最佳的理想交叉口使得實際交叉口集中在理想交叉口附近使得綠信比損失最小,那么子區(qū)內(nèi)部綠波協(xié)調(diào)能得到最大綠波帶。
初始條件
如圖1所示,子區(qū)內(nèi)某干線的交叉口有g1,g2…gu,交叉口序號為ζ(ζ=1,2…u);子系統(tǒng)的初始周期為cβ;相鄰交叉口間距為sa(α=1,2…u-1);各交叉口沿干道方向相位的綠信比為aζ(ζ=1,2…u)。
計算流程
(1)首先求出車輛以系統(tǒng)速度在1/2周期內(nèi)可運行的距離L,不失一般性,以10米為單位計。然后以L為中心,在它前后各取10個數(shù)值,即21個數(shù)值,作為理想交叉口間距mμ(μ=1,2…21),據(jù)此確定出理想交叉口最佳間距的范圍
(2)假定第一個理想交叉口與第一個實際交叉口重合,則其余的實際交叉口會與理想交叉口有一定的挪移量,如圖3所示。計算出理想交叉口的間距為mμ下u個實際交叉口距離最近理想交叉口的挪移量
圖3 實際交叉口與理想交叉口挪移圖示
式中:Δαmμ,i代表在理想交叉口間距為mμ下交叉口ζ距離最近理想交叉口的挪移量。
根據(jù)式 (7)、(8),可以求得21個理想交叉口間距值下的各個交叉口距離最近理想交叉口的挪移量,整理得到挪移量矩陣Δα21×u
(3)為了計算各交叉口與理想交叉口的挪移量差值,首先將挪移量矩陣Δα21×u中每一行中的值按從小到大順序排列得到矩陣β21×u
然后將矩陣β21×u每一行相鄰兩個數(shù)相減,得到挪移量差值矩陣Δβ21×(u-1)
(4)找出矩陣Δβ21×(u-1)中最大值Δβmax=Δβmopt,ζ,Δβmax所對應的理想交叉口的間距為mopt。根據(jù) “挪移差值最大原則”,即挪移差值越大,實際交叉口位置越集中在理想交叉口附近。因此理想交叉口的最佳間距為mopt。
(5)調(diào)整初始周期,使周期時長與最佳間距相對應。根據(jù)步驟④的最佳間距mopt計算子區(qū)的最佳周期
式中:copt——子區(qū)最佳周期。
(6)根據(jù)最佳理想間距mopt以及最大挪移量Δβmax,確定實際交叉口與理想交叉口之間的最大偏移量為δ,以及最大偏移量所對應的實際交叉口序號為λ
(7)以交叉口λ為基點向前推進δ即為一個理想交叉口,然后根據(jù)實際交叉口間距與理想交叉口間距遞推求出其他交叉口與理想交叉口的偏移量,并確定實際交叉口臨近的是奇數(shù)理想交叉口還是偶數(shù)理想交叉口
式中:lζ——實際交叉口與理想交叉口的偏移量;tζ——交叉口奇偶標識符。
(8)根據(jù)步驟 (7),臨近奇數(shù)理想交叉口的一組交叉口采用同步協(xié)調(diào),另一組交叉采用同步協(xié)調(diào),兩組交叉口之間采用交互協(xié)調(diào)。根據(jù)協(xié)調(diào)方式確定各交叉口的最佳相位差。
式中:ζ代表各交叉口的相位差。
在單點綠信比優(yōu)化中,交叉口獨立實行控制,信號機根據(jù)交叉口當前的飽和度、當前綠信比方案、以及信號機內(nèi)儲存的可供選擇綠信比方案基于飽和度估計來確定該交叉口下一周期所執(zhí)行的綠信比方案。
目前綠信比方案的確定采用 Webster配時算法[12]。各相位的綠信比是根據(jù)該相位流率比之間的比例進行分配
式中:aj代表該交叉口相位序號為j的綠信比。信號機根據(jù)該交叉口的歷史流量存儲一定數(shù)量的綠信比方案。
綠信比優(yōu)化算法如下:
步驟1 根據(jù)當前各相位的飽和度和當前綠信比方案估計下一周期采用系統(tǒng)中N套綠信比方案時該交叉口各相位的飽和度
式中:t代表當前周期時刻;t+1代表下一周期時刻;dt+1i,k,j代表該交叉口下一周期采用第k套方案j相位預計飽和度;dti,j代表當前該交叉口相位j的飽和度;aj代表當前綠信比方案中相位j的綠信比;ak.j代表第k套方案中相位j的綠信比。
步驟2 估計下一周期使用N套方案后該交叉口的飽和度
步驟3 從步驟2選出的所有方案中交叉口預計飽和度的最小值,則下一周期信號機采用這個值對應的方案號后可以使得該交叉口飽和度最低,運行效果最好。
已知某城市區(qū)域交叉口以及相鄰交叉口間距如圖4所示,將該區(qū)域沿干線方向劃分成3個子區(qū)A、B、C。規(guī)定車輛行駛速度為v=40km/h。
圖4 實際交叉口與理想交叉口挪移
某一時刻由各自子區(qū)的關鍵路口確定其初始周期,如表1所示。
表1 各子區(qū)的初始周期
根據(jù)周期原則,將子區(qū)A與子區(qū)B合并成新的子區(qū)AB,新系統(tǒng)周期為cAB=65s,對子區(qū)AB、子區(qū)C進行綠波協(xié)調(diào)控制。首先按照綠波協(xié)調(diào)控制步驟①計算子區(qū)AB的理想交叉口的最佳間距區(qū)間mμ為 [25,45],然后通過步驟②③算出子區(qū)AB在理想交叉口間距mμ在25~45的條件下交叉口最大挪移量差值Δβmax=23,對應理想交叉口的間距為mμ=45,因此子區(qū)AB交叉口最佳間距為mopt=45,即450m,最佳周期為:c==81s,接著根據(jù)步驟⑥計算距離最近理想交叉口的最大偏移量為δ==11,對應的交叉口號為B3,根據(jù)步驟⑦計算出其他交叉口偏移量lζ,最后根據(jù)步驟⑧計算出每個交叉口的相位差ζ。同理按上述方法計算出子區(qū)C的最佳周期和各個交叉口的相位差。子區(qū)AB和子區(qū)C的綠波效果如表2,表3所示,其中損失綠信比、有效綠信比、綠波帶寬度計算公式分別為
式中:Δaζ——相位ζ損失綠信比,a′ζ——相位ζ有效綠信比,D——子區(qū)的綠波帶寬度。
表2 子系統(tǒng)AB綠波協(xié)調(diào)表
表3 子系統(tǒng)C綠波協(xié)調(diào)表
上述數(shù)據(jù)表明:對合并后的子系統(tǒng)AB和子系統(tǒng)C進行綠波協(xié)調(diào)控制,綠波帶寬度分別為0.365、0.388,大于文獻 [13]中32.5%,綠波效果良好。
對子區(qū)C的關鍵交叉口C5進行單點綠信比優(yōu)化,建立交通流量模型,模型中交通流量數(shù)據(jù)來自國內(nèi)某三線城市某個關鍵交叉口的高峰期數(shù)據(jù),并且按照該交叉口的歷史交通流量數(shù)據(jù)制定了4套綠信比方案,以交叉口每個周期的停車次數(shù)為評價標準算出4套綠信比方案以及優(yōu)化綠信比方案運行60個周期的停車次數(shù)。仿真結果如圖5和表4所示??梢杂^察出相對于傳統(tǒng)區(qū)域控制的方案中全程采用某個單套方案的數(shù)據(jù),采用優(yōu)化方案后總停車次數(shù)以及平均停車次數(shù)大大降低,其中平均停車次數(shù)至少減少65%。
表4 平均停車次數(shù)
圖5 綠信比優(yōu)化算法仿真數(shù)據(jù)
本文首次提出在子區(qū)的基礎上,實時動態(tài)地對子區(qū)進行拆分合并、對子區(qū)內(nèi)部的干線進行綠波協(xié)調(diào)控制以及對各交叉口進行獨立地綠信比優(yōu)化,得到適應當前交通流的最佳周期、最佳相位差、最優(yōu)綠信比方案,成功將點控、線控、面控有機的結合在一起,達到將整個區(qū)域進行協(xié)調(diào)控制的目的,并且方案中系統(tǒng)所需的輸入量少、算法簡單。結果表明,使用該方案后子區(qū)內(nèi)部的干線具有較寬的綠波帶寬度,單個交叉口的平均停車次數(shù)大大減少,具有良好的控制效果,有一定的實際意義。
[1]LI Qunzu,XIA Qingguo,BA Mingchun,et al.Status and development of urban traffic signal control system [J].Science Technology and Engineering,2009,24 (9):7438-7435 (in Chinese).[李群祖,夏清國,巴明春,等.城市交通信號控制系統(tǒng)現(xiàn)狀與發(fā)展[J].科學技術與工程,2009,24 (9):7436-7442.]
[2]Jatmiko M,Azurat A,Herry A.Self-organizing urban traffic control architecture with swarm-self organizing map in Jakarta:Signal control system and simulator [J].International journal on smart sensing and intelligent systems,2010,3 (3):444-465.
[3]YANG Zhaosheng,LIN Ciyun,GONG Bowen,et al.Generation technique of adaptive control strategy for SCATS combined with probe vehicle [J].Journal of Jilin University (Engineering and Technology Edition),2010,40 (1):35-41 (in Chinese).[楊兆升,林賜云,龔勃文,等.結合浮動車的SCATS自適應控制策略生成技術 [J].吉林大學學報 (工學版),2010,40 (1):35-41.]
[4]Bus O L,Errnst D,De S B.Cross-entropy optimization of control polices with adaptive basis functions [J].IEEE Transaction on Automatic Control,Man and Cybernetics,Part B,2011,41(1):196-209.
[5]Blanc N,Kroening D,Sharygina N.Scoot:A tool for the analysis of System C models [C]//Tools and Algorithms for Construction and Analysis of System,2008:467-470.
[6]XU Jianmin,SHOU Yanfang,LU Kai.Area coordinated control of traffic signals based on bi-level programming model [J].Journal of South China University of Technology (Natural Science Edition),2011,39 (3):95-106 (in Chinese).[徐建閩,首艷芳,盧凱.基于雙層規(guī)劃模型的交通信號區(qū)域協(xié)調(diào)控制 [J].華南理工大學學報 (自然科學版),2011,39 (3):95-106.]
[7]SUN Youxin,QIAN Yongsheng,WANG Chunlei,et al.Traffic signal control based on hierarchic parallel catastrophic particle swarm optimization algorithm [J].Computer Engineering and Design,2008,29 (6):1497-1500 (in Chinese).[孫有信,錢勇生,王春雷,等.基于分層并行災變粒子群算法的交通控制[J].計算機工程與設計,2008,29 (6):1497-1500.]
[8]LU Kai,XU Jianmin,ZHENG Shujian,et al.Research on fast dynamic division method of coordinated control subarea [J].Acta Automatica Sinica,2012,38 (2):279-286 (in Chinese).[盧凱,徐建閩,鄭淑建,等.協(xié)調(diào)控制子區(qū)快速動態(tài)劃分方法研究 [J].自動化學報,2012,38 (2):279-286.]
[9]YANG Xiaoguang,ZHAO Jing,WANG Tao.Optimal cycle calculation method of signal control at roundabout [J].China Journal of Highway and Transport,2008,21 (6):90-94 (in Chinese).[楊曉光,趙靖,汪濤.環(huán)形交叉口信號控制最佳周期計算方法 [J].中國公路學報,2008,21 (6):90-94.]
[10]SHEN Guojing,XU Weiming.Study on traffic trunk dynamic two-direction green wave control technique [J].Journal of Zhejiang University (Engineering Science),2008,42 (9):1625-1630(in Chinese).[沈國江,許衛(wèi)明.交通干線動態(tài)雙向綠波帶控制術研究 [J].浙江大學學報 (工學版),2008,42 (9):1625-1630.]
[11]SU Bing,F(xiàn)ANG Lin,YANG Qian.How to finding apath to zoning a traffic network efficiently [J].Information Technology Journal,2012,11 (4):551-553.
[12]TAO Zeming,JIANG Xiancai,PEI Yulong.An improve green ratio model based on attemperment priority in oversaturated traffic [C]//IEEE International Conference on Intelligent Computation Technology and Automation,2008:563-567.
[13]WANG Dianhai,YANG Xirui,SONG Xianming.Improvement of classical numerical method for arterial road signal coordinate control [J].Journal of Jilin University (Engineering and Technology Edition),2011,41 (1):29-34 (in Chinese).[王殿海,楊希瑞,宋現(xiàn)敏.交通信號干線協(xié)調(diào)控制經(jīng)典數(shù)值計算法的改進 [J].吉林大學學報 (工學版),2011,41 (1):29-34.]