柴 干 高祥云

(東南大學(xué)智能運(yùn)輸系統(tǒng)研究中心,南京 210096)

隨著高速公路車輛保有量的增長(zhǎng),路段擁擠狀況不斷出現(xiàn).當(dāng)路段出現(xiàn)擁擠時(shí),會(huì)形成道路瓶頸,產(chǎn)生嚴(yán)重的車輛延誤現(xiàn)象[1].若未對(duì)擁擠路段進(jìn)行及時(shí)控制,上游車輛對(duì)其不斷沖擊將導(dǎo)致交通環(huán)境惡化,降低路網(wǎng)交通的整體效益.緩解交通擁擠,提高運(yùn)營(yíng)效率與安全性,已經(jīng)成為亟待解決的現(xiàn)實(shí)問(wèn)題[2].研究表明,匝道協(xié)調(diào)控制是緩解高速公路擁擠的最有效控制方式[3].開(kāi)展擁擠條件下的匝道協(xié)調(diào)控制方法研究,對(duì)緩解擁擠、提高通行能力具有理論與工程應(yīng)用價(jià)值.匝道協(xié)調(diào)控制主要包括最優(yōu)控制、多變量控制、模糊控制和啟發(fā)式控制.最優(yōu)控制通過(guò)優(yōu)化性能指標(biāo),使主線交通處于期望狀態(tài)[4],但算法復(fù)雜、計(jì)算量大,難以建立擁擠下的控制方案.多變量控制使交通狀態(tài)維持在設(shè)定值附近,含有反饋機(jī)制,可減少誤差并抵抗干擾,但擁擠時(shí)控制效果并不理想[5].模糊控制不需要抽象模型,利用專家經(jīng)驗(yàn)克服原始數(shù)據(jù)的誤差[6],但對(duì)于多個(gè)匝道的協(xié)調(diào)控制,模糊規(guī)則呈指數(shù)增加,難以應(yīng)用于實(shí)際中.啟發(fā)式控制無(wú)需確定控制目標(biāo),基于經(jīng)驗(yàn)和實(shí)測(cè)數(shù)據(jù),采取局部與協(xié)調(diào)結(jié)合的方法生成調(diào)節(jié)率,技術(shù)復(fù)雜性低、可實(shí)施性好[7].

啟發(fā)式協(xié)調(diào)控制方法已經(jīng)成功應(yīng)用于美國(guó)高速公路的匝道控制中[5].其中,Zone算法具有實(shí)時(shí)性與靈活性,但由于控制目標(biāo)不清楚且沒(méi)有充分考慮交通突變,難以應(yīng)對(duì)交通擁擠;Helper算法無(wú)需準(zhǔn)確的OD信息和交通流模型,但由于缺乏優(yōu)化目標(biāo),難以實(shí)現(xiàn)優(yōu)化控制,并且不具備預(yù)測(cè)機(jī)制,難以避免時(shí)滯引起的控制誤差.Swarm算法依據(jù)預(yù)測(cè)交通參數(shù)識(shí)別交通瓶頸,期望在萌芽階段消除交通擁擠;但是該方法需要準(zhǔn)確的預(yù)測(cè)模型和精確的OD數(shù)據(jù),才能實(shí)現(xiàn)調(diào)節(jié)率的合理選擇.Bottleneck控制方法具有實(shí)時(shí)性、靈活性和邏輯簡(jiǎn)單等特點(diǎn),可通過(guò)改進(jìn)控制策略的魯棒性和時(shí)滯性,提高解決路段交通擁擠的有效性.

為了充分利用Bottleneck控制方法的優(yōu)點(diǎn),克服其缺陷,以解決擁擠條件下的匝道控制問(wèn)題,本文依據(jù)擁擠持續(xù)時(shí)間內(nèi)交通流波陣面向上游移動(dòng)的距離,推導(dǎo)擁擠影響范圍,獲得需要調(diào)節(jié)的匝道數(shù).然后,通過(guò)變換控制變量,對(duì)局部調(diào)節(jié)率和協(xié)調(diào)調(diào)節(jié)率進(jìn)行優(yōu)化,獲得改進(jìn)的Bottleneck控制方法,克服原方法超調(diào)量過(guò)大、調(diào)節(jié)時(shí)間過(guò)長(zhǎng)的缺陷.

1 Bottleneck控制方法

1.1 局部調(diào)節(jié)率

通過(guò)測(cè)得匝道上、下游鄰近處的占有率,根據(jù)調(diào)節(jié)率-占有率曲線,采用在線查表法[8]確定該匝道的局部調(diào)節(jié)率,即

(1)

式中,R(j,k)為周期k內(nèi)匝道j的局部調(diào)節(jié)率;Pu,Pv分別為占有率-調(diào)節(jié)率曲線上控制點(diǎn)u,v的相應(yīng)占有率;Ru,Rv為Pu,Pv對(duì)應(yīng)的調(diào)節(jié)率;P(j,k)為周期k內(nèi)匝道j上游的主線占有率.

1.2 協(xié)調(diào)調(diào)節(jié)率

若路段i下游占有率超過(guò)閾值,出現(xiàn)交通瓶頸,則

P(i,k-1)≥Pc(i)

(2)

Qup(i,k-1)+Qon(i,k-1)≥Qoff(i,k-1)+

Qdown(i,k-1)

(3)

式中,P(i,k-1)為周期k-1內(nèi)路段i的下游占有率;Pc(i)為路段i的下游占有率閾值;Qup(i,k-1),Qdown(i,k-1)分別為周期k-1內(nèi)路段i的上、下游流進(jìn)交通量;Qon(i,k-1),Qoff(i,k-1)分別為周期k-1內(nèi)路段i內(nèi)的入、出口匝道流出交通量.

針對(duì)每一個(gè)瓶頸路段i,計(jì)算其上游各入口匝道調(diào)節(jié)率的降低總量,并按照經(jīng)驗(yàn)權(quán)重將降低總量分配至相關(guān)入口匝道.權(quán)重系數(shù)根據(jù)上游入口匝道到路段i的距離設(shè)定,即

Qreduction(i,k)=[Qup(i,k-1)+Qon(i,k-1)]-

[Qoff(i,k-1)+Qdown(i,k-1)]

(4)

(5)

式中,Qreduction(i,k)為周期k內(nèi)路段i上游關(guān)聯(lián)入口匝道流入的降低流量;Rreduction(j,i,k)為周期k內(nèi)路段i上游關(guān)聯(lián)入口匝道j流入的降低流量;n為需要協(xié)調(diào)的匝道數(shù);Wji為路段i上游關(guān)聯(lián)入口匝道j的權(quán)重系數(shù).各瓶頸路段的影響區(qū)域可能發(fā)生重疊,每個(gè)入口匝道常會(huì)處于不同瓶頸路段的影響區(qū)域中.因此,需要選取各瓶頸路段的入口匝道流入交通降低值中的最大值,作為各個(gè)匝道的協(xié)調(diào)調(diào)節(jié)率,即

(6)

式中,Qon(j,k-1)為周期k-1內(nèi)匝道j的實(shí)際流入量;r(j,k)為周期k內(nèi)匝道j的協(xié)調(diào)調(diào)節(jié)率.

1.3 最終調(diào)節(jié)率

將局部層面和協(xié)調(diào)層面上匝道j的調(diào)節(jié)率R(j,x)和r(j,k)進(jìn)行比較,取最小值作為匝道j的最終調(diào)節(jié)率Rend,即

Rend=min[R(j,k),r(j,k)]

(7)

Bottleneck控制方法的不足在于:① 協(xié)調(diào)層面上,缺乏需要調(diào)節(jié)的入口匝道數(shù)量和位置的準(zhǔn)確表達(dá)式.② 主線最大流量隨時(shí)間、天氣等因素的變化而變化,故將流量作為控制變量,難以確切反映主線的交通瓶頸;臨界密度則相對(duì)穩(wěn)定且不受天氣等因素的影響,故將密度作為控制變量,適于反映主線的交通擁擠狀態(tài).③ 頻繁對(duì)局部調(diào)節(jié)率和協(xié)調(diào)調(diào)節(jié)率取極小值,易形成主線交通運(yùn)量的欠載或過(guò)載現(xiàn)象.

2 改進(jìn)的Bottleneck控制方法

2.1 控制匝道數(shù)的識(shí)別

根據(jù)路段擁擠時(shí)交通波運(yùn)行機(jī)理,應(yīng)用Greenshields模型推導(dǎo)擁擠影響范圍,并結(jié)合匝道樁號(hào)確定所需調(diào)節(jié)的匝道數(shù)與位置.設(shè)ρi為路段i的堵塞密度,路段i有2個(gè)相鄰且不同密度的區(qū)域ρi1和ρi2,用垂直面S加以分割(稱S為波陣面);交通波的波速為ω,且按圖中箭頭x方向運(yùn)行(見(jiàn)圖1).

圖1 交通波示意圖

高速公路發(fā)生擁擠致使車輛難以通行時(shí),上游車輛將在波陣面S處開(kāi)始排隊(duì).根據(jù)交通波波速和Greenshields模型的車流速度-密度線性關(guān)系,求得擁擠點(diǎn)上游影響范圍L為

(8)

式中,μf為自由流車速；T為擁擠持續(xù)時(shí)間.

擁擠狀態(tài)下的匝道控制即是對(duì)擁擠點(diǎn)上游流入交通進(jìn)行控制.從探測(cè)到擁擠時(shí)刻開(kāi)始,就需要估計(jì)匝道控制范圍,確定需要進(jìn)行協(xié)調(diào)的匝道數(shù)n,其計(jì)算公式為

n={L,Xn}

(9)

式中,Xn為波陣面上游第n個(gè)匝道的樁號(hào).傳統(tǒng)的匝道協(xié)調(diào)方法沒(méi)有給出需調(diào)節(jié)的匝道數(shù)量與位置,故實(shí)際操作性不強(qiáng).此處定義XS為波陣面樁號(hào),在[XS,XS+L]內(nèi),得到用于辨識(shí)匝道數(shù)n及其位置的表達(dá)式為

(10)

n的確定還與高速公路具體的結(jié)構(gòu)有關(guān).當(dāng)L相同、匝道間距不同時(shí),n的取值一般不同.

2.2 匝道流入交通降低量的確定

控制變量是反映交通控制系統(tǒng)特性的關(guān)鍵參數(shù),其選取的合理性直接關(guān)系到控制效果.本節(jié)將依據(jù)流量與密度的關(guān)系,闡述應(yīng)用密度代替流量作為控制變量的原理,并采用比例積分控制方法,推導(dǎo)流入交通降低量的表達(dá)式.交通流量和密度的關(guān)系如圖2所示.圖中,ρm為最大流量qm對(duì)應(yīng)的密度.任一流量值q對(duì)應(yīng)2個(gè)密度值ρi1和ρi2,前者小于ρm(即位于不擁擠區(qū)域),后者大于ρm(即位于擁擠區(qū)域).交通擁擠時(shí),事發(fā)點(diǎn)上游密度變大,下游密度變小,而上、下游交通量的計(jì)算值可能相差并不大[9],故Qreduction可能很小.若按照經(jīng)驗(yàn)權(quán)重將Qreduction分配至各個(gè)關(guān)聯(lián)匝道,則難以實(shí)現(xiàn)擁擠條件下的控制效果,故需要對(duì)Qreduction的計(jì)算公式進(jìn)行改進(jìn).

圖2 交通流量和密度的關(guān)系圖

PID控制器采用反饋控制方式控制變量與期望值的誤差,其積分單元用于消除穩(wěn)態(tài)誤差,比例和微分單元用于加快控制反饋.由于微分項(xiàng)產(chǎn)生沖擊效應(yīng),PID控制器常用于反饋速度較慢的系統(tǒng)[10].將密度作為控制變量,則

Qreduction(i,k)=-[K1e(i,k)+K2(e(i,k)-e(i,k-1))]

(11)

式中,K1,K2為PID的控制參數(shù)[11],由實(shí)驗(yàn)獲得;e(i,k)為周期k內(nèi)的控制誤差.將e(i,k)=ρd-ρS(i,k)代入式(11)可得

Qreduction(i,k)=-(K1+K2)[ρd-ρS(i,k)]+

K2[ρd-ρS(i,k-1)]=-K1ρd+

(K1+K2)ρS(i,k)-K2ρS(i,k-1)

(12)

式中,ρd為擁擠點(diǎn)上游的期望密度;ρS(i,k)為周期k內(nèi)路段i中波陣面S的上游密度.

將式(12)代入式(5),則

Rreduction(j,i,k)=[-K1ρd+(K1+K2)ρS(i,k)-

(13)

2.3 最終調(diào)節(jié)率的優(yōu)化

路段擁擠時(shí),Bottleneck控制方法取局部調(diào)節(jié)率與協(xié)調(diào)調(diào)節(jié)率中的較小值,這樣雖然可以迅速減小路段密度,但當(dāng)路段沒(méi)有發(fā)生擁擠或交通密度小于期望密度時(shí),容易使超調(diào)量過(guò)大,并產(chǎn)生較長(zhǎng)的調(diào)節(jié)時(shí)間.對(duì)其進(jìn)行改進(jìn)的思路為:當(dāng)被調(diào)密度大于期望密度時(shí),按照最嚴(yán)格的方式調(diào)節(jié)(取局部調(diào)節(jié)率與協(xié)調(diào)調(diào)節(jié)率中的較小值);當(dāng)實(shí)際密度小于期望密度時(shí),對(duì)調(diào)節(jié)率進(jìn)行重新取值,最終調(diào)節(jié)率取局部調(diào)節(jié)率和協(xié)調(diào)調(diào)節(jié)率的加權(quán)組合,即

(14)

式中,ρtarget=35 pcu/(km·lane)為路段i的密度理想值[11];αj為匝道j的比例因子，且0≤αj≤1.

3 仿真示例與分析

選取寧杭高速溧陽(yáng)西至丁山6個(gè)路段作為示例路徑,路徑拓?fù)鋱D見(jiàn)圖3.

圖3 寧杭高速公路路徑的拓?fù)鋱D

3.1 參數(shù)初始化

交通擁擠發(fā)生時(shí),檢測(cè)器檢測(cè)到的各路段密度和調(diào)節(jié)率初值見(jiàn)表1.由表可知,第6路段(即丁山路段)的密度大于期望密度,發(fā)生擁擠.

表1 各個(gè)路段的初始密度與初始調(diào)節(jié)率

根據(jù)該路段交通運(yùn)行的基本狀況,設(shè)期望密度為35 pcu/(km·lane);上游主線進(jìn)入的流量平均值為6 000 pch/h;仿真步長(zhǎng)為20 s;仿真時(shí)間為2 000 s.仿真實(shí)驗(yàn)中,K1=97,K2=29時(shí),期望密度的動(dòng)態(tài)跟蹤性能較好.對(duì)寧杭高速上296起碰撞或拋錨等事件持續(xù)時(shí)間進(jìn)行統(tǒng)計(jì)分析,取事件持續(xù)時(shí)間T=60 min,代入式(8)得事件影響范圍為19.6 km.根據(jù)式(10)得事件影響范圍內(nèi)的匝道為丁山、宜興和西塢.協(xié)調(diào)層面的調(diào)節(jié)率分配權(quán)值可根據(jù)匝道與擁擠斷面的距離決定,離擁擠斷面越遠(yuǎn),匝道承擔(dān)的調(diào)節(jié)率降低值越小.

3.2 比例因子的敏感性分析

根據(jù)改進(jìn)的Bottleneck控制方法,應(yīng)用數(shù)字仿真技術(shù),分別對(duì)丁山、宜興、西塢3個(gè)入口匝道的比例因子αj(j=1,2,3)進(jìn)行敏感性分析.當(dāng)αj=0,0.2,0.4,0.6,0.8,1.0時(shí),相應(yīng)路段密度隨時(shí)間變化的曲線見(jiàn)圖4.

當(dāng)權(quán)重系數(shù)為1.0時(shí),密度曲線快速且平穩(wěn)地變化,超調(diào)量小,調(diào)節(jié)時(shí)間少,幾乎無(wú)振蕩,控制效果最為理想,可以使路段交通密度快速平穩(wěn)地達(dá)到所期望的目標(biāo)值.

由此可知,針對(duì)本示例,式(14)實(shí)際演變?yōu)?/p>

(15)

3.3 控制算法的比較

采用Bottleneck控制方法、改進(jìn)的Bottleneck控制方法以及未施加控制方法時(shí),對(duì)擁擠路段及其上游相鄰路段進(jìn)行仿真分析，路段車流密度隨時(shí)間的變化軌跡見(jiàn)圖5.由圖可知,未施加控制方案時(shí),3個(gè)路段的擁擠狀況隨時(shí)間延長(zhǎng)而加劇,并逐步影響上游各路段.在路段發(fā)生擁擠時(shí),利用Bottleneck控制方法,取局部調(diào)節(jié)率與協(xié)調(diào)調(diào)節(jié)率中的較小值可以迅速減小路段密度.當(dāng)路段沒(méi)有發(fā)生擁擠時(shí),仍然采用局部調(diào)節(jié)率與協(xié)調(diào)調(diào)節(jié)率中的較小值,將產(chǎn)生較大的超調(diào)量和較長(zhǎng)的調(diào)節(jié)時(shí)間.控制變量需要經(jīng)過(guò)1 500 s才能穩(wěn)定在理想密度上,且期望密度與路段密度的最大絕對(duì)誤差達(dá)到11 pcu/(km·lane).利用改進(jìn)的Bottleneck控制方法,能夠迅速使擁擠的第6路段(丁山段)在450 s時(shí)消散擁擠,達(dá)到平穩(wěn)交通流狀態(tài).由于采取了優(yōu)化的匝道調(diào)節(jié)率,有效避免了原Bottleneck控制方法的不足,使調(diào)節(jié)時(shí)間減少70%,超調(diào)量減小75%,無(wú)余差且?guī)缀鯚o(wú)振蕩.利用該方法不僅能迅速疏導(dǎo)路段擁擠,而且能夠使擁擠路段的上游路段密度快速趨向期望密度,保持路徑上交通流運(yùn)行的通暢與穩(wěn)定.

圖4 比例因子的敏感性分析

圖5 控制方法的比較結(jié)果

4 結(jié)語(yǔ)

基于交通波理論推導(dǎo)的交通擁擠影響范圍,能夠辨識(shí)影響范圍內(nèi)需要協(xié)調(diào)控制的匝道數(shù)量及位置,為改進(jìn)Bottleneck控制方法提供關(guān)聯(lián)入口匝道權(quán)重系數(shù)的計(jì)算依據(jù).針對(duì)Bottleneck控制方法使用流量累積判別動(dòng)態(tài)交通狀況,難以檢測(cè)主線交通擁擠狀態(tài)的缺陷,以密度代替流量作為控制變量,設(shè)計(jì)PID控制器,構(gòu)建改進(jìn)的Bottleneck控制方法,從而將流量囤積調(diào)節(jié)轉(zhuǎn)變?yōu)槊芏葥頂D控制.調(diào)節(jié)路段流量累積，使控制變量接近期望值時(shí)會(huì)產(chǎn)生調(diào)節(jié)過(guò)度現(xiàn)象,通過(guò)對(duì)局部調(diào)節(jié)率和協(xié)調(diào)調(diào)節(jié)率進(jìn)行組合,確定優(yōu)化的最終調(diào)節(jié)率,則可解決調(diào)節(jié)過(guò)度的問(wèn)題.

)

[1]Papamichail I,Papageorgiou M.Heuristic ramp-metering coordination strategy implemented at monash freeway[J].TransportationResearchRecord,2010,2179: 10-20.

[2]張海軍,張鈺,楊曉光.異常事件下高速道路交通狀態(tài)的分析與仿真[J].交通運(yùn)輸工程學(xué)報(bào),2008,8(2):116-121.

Zhang Haijun,Zhang Jue,Yang Xiaoguang.Analysis and simulation of traffic state on expressway during incident[J].JournalofTrafficandTransportationEngineering,2008,8(2): 116-121.(in Chinese)

[3]Chai Gan,Xie Yongli.A genetic-fuzzy control application to ramp-mainstream joint coordination control of the freeway[J].AdvancedScienceLetters,2012,7: 369-374.

[4]Kotsialos A,Papageorgiou M,Mangeas A.Coordinated and integrated control of motorway networks via nonlinear optimal control[J].TransportationResearchPartC:EmergingTechnologies,2002,10(1): 65-84.

[5]Papageorgiou M,Kotsialos A.Freeway ramp metering: an overview[J].IEEETransactionsonIntelligentTransportationSystems,2002,3(4): 171-181.

[6]溫凱歌,曲仕茹,張玉梅.基于模糊邏輯的高速公路入口匝道控制方法[J].中國(guó)公路學(xué)報(bào),2007,20(6): 100-104.

Wen Kaige,Qu Shiru,Zhang Yumei.Method for freeway on-ramp control based on fuzzy logic[J].ChinaJournalofHighwayandTransport,2007,20(6): 100-104.(in Chinese)

[7]Papamichail I,Papageorgiou M.Traffic-responsive linked ramp-metering control[J].IEEETransactionsonIntelligentTransportationSystems,2008,9(1): 111-121.

[8]Jacobsen L,Henry K.Real-time metering algorithm for centralized control[J].TransportationResearchRecord,1989,1232: 17-26.

[9]Zhang H M.Coordinated Traffic responsive ramp control via nonlinear state feedback[J].TransportationResearchPartC:EmergingTechnologies,2001,9(5): 337-352.

[10]Sheu J B.Stochastic modeling and real-time prediction of vehicular lane-changing behavior[J].TransportationResearchPartB:Methodological,2001,35(7): 695-716.

[11]梁新榮,劉志勇,徐建閩.高速公路多匝道協(xié)調(diào)控制系統(tǒng)設(shè)計(jì)與仿真[J].信息與控制,2005,34(3):308-311.

Liang Xinrong,Liu Zhiyong,Xu Jianmin.Design and simulation of coordinated ramp control system on freeway[J].InformationandControl,2005,34(3): 308-311.(in Chinese)