孫浩然,劉明浩,姜志偉,楊盼盼
(長安大學 電子與控制工程學院,西安 710064)
在冰雪天氣下,高速公路道路環(huán)境、路面狀態(tài)和車輛性能均會發(fā)生一定程度的變化,駕駛?cè)说囊暰€、情緒受到負面干擾,車輛制動性能受限,發(fā)生交通事故的幾率增大[1-2],如何保障冰雪天氣下的行車安全,是高速公路運營管理部門面臨的一項重要挑戰(zhàn)。
限速控制是冰雪天氣下保障高速公路行車安全的常用措施。傳統(tǒng)限速方法限速值確定后一般不會改變,屬于靜態(tài)限速方法[3]。然而,高速公路是一個由人、車、路和環(huán)境組成的復雜動態(tài)系統(tǒng),傳統(tǒng)靜態(tài)限速方法無法在變化的道路交通環(huán)境下保障管控措施的持續(xù)適用,有一定的局限性[4]。靈活可變的限速值能夠使駕駛員行駛速度始終保持在當前道路交通環(huán)境下的安全范圍內(nèi),從而提升冰雪天氣下高速公路行車安全性[5]。
為滿足冰雪天氣下高速公路安全行車需求,研究人員對可變限速管控策略進行了持續(xù)探索。Chen等[6]通過模擬不同天氣條件下的跟車場景,對不同天氣條件和道路類型下的駕駛員感知能力進行了定量評估。馬超等[7]對包括冰雪天氣在內(nèi)的不同惡劣天氣類型與等級條件下的安全行車車速給出建議值,建立了惡劣天氣條件下高速公路安全限速評估體系。姚冬冬等[8]建立了考慮道路線型和路面摩擦系數(shù)的安全限速模型,并提出了一種冰雪條件下基于模糊控制的高速公路可變限速策略。Guo等[9]分析了惡劣天氣因素和道路線型與坡道通行能力折減對通行的影響,提出惡劣天氣下的可變限速控制策略。
然而,目前惡劣天氣下高速可變限速的研究多集中在通過限速值的選取保障交通安全上,而針對冰雪天氣條件下限速管控的同時有效保障通行效率的研究相對較少[10]。鑒于此,本文研究冰雪天氣下兼顧行車安全與通行效率的高速公路可變限速策略。從冰雪天氣高速公路交通特性出發(fā),分析冰雪天氣條件下的安全限速模型,探索一種根據(jù)冰雪天氣變化進行動態(tài)路況預測的可變限速控制方法,以改善冰雪天氣對高速公路行車效率與安全造成的負面影響。最后通過對比仿真實驗,驗證了所提出方法的有效性。
高速公路作為一個復雜的動態(tài)系統(tǒng),交通流各參數(shù)及其之間的關(guān)系會隨著交通流演變以及路況環(huán)境的改變發(fā)生變化,可變限速控制通過對交通流運行狀態(tài)進行人工干預,使高速公路交通流保持在較為穩(wěn)定的運行狀態(tài),實現(xiàn)交通安全與效率的均衡。本文將宏觀交通流預測應用于限速策略的反饋調(diào)節(jié),整體管控流程如圖1所示。
圖1 冰雪天氣下高速路可變限速管控流程圖
從圖1可以看出,冰雪天氣下高速路可變限速管控流程為:每個控制周期內(nèi),通過檢測器對交通流數(shù)據(jù)和氣象數(shù)據(jù)以及路況信息進行采集并進行最大安全限速值更新;以交通流當前時刻運行狀態(tài)作為輸入,結(jié)合交通流預測模型與優(yōu)化算法對控制目標函數(shù)進行求解,得出相應的限速策略,通過宏觀交通流預測模型對限速調(diào)控效果進行反饋校正,最終得出合理可行的限速策略作用于系統(tǒng)。隨著時間的推進,周期性地重復上述過程。
為保障冰雪天氣下管控路段的行車安全,首先應保證停車視距滿足要求。車輛的停車視距指從駕駛員發(fā)現(xiàn)前方障礙物開始,做出制動動作,使車輛從初始車速逐漸減速直至完全停止時所行駛的距離。本文將制動減速度作為關(guān)鍵參數(shù),從車輛運動學角度構(gòu)建停車視距計算模型,使停車視距計算更加合理。車輛的制動過程分析如圖2所示。
圖2 車輛制動過程分析圖
整個過程分為駕駛員反應時間、制動系統(tǒng)協(xié)調(diào)時間以及全制動時間3個階段,首先,當前方緊急情況后,駕駛員經(jīng)過視覺感知、危險識別以及判斷決策進行制動操作,此為駕駛員反應階段;受車輛制動結(jié)構(gòu)以及駕駛?cè)瞬忍ぐ逅俣扔绊?,因此?jīng)制動系統(tǒng)響應后才會形成地面制動力,制動減速度逐漸由零增加至最大值,此為制動系統(tǒng)協(xié)調(diào)階段;此后制動減速度保持最大值不變,直至車輛完全停止,此為全制動階段。
將制動力提升階段近似視作線性變化過程[11],可得車輛制動過程的最小停車視距計算公式如下:
(1)
式中,V0為車輛的初始速度;τ1為駕駛員反應時間;V為制動力提升階段任意時刻t的車輛速度;τ2為制動系統(tǒng)協(xié)調(diào)響應所需時間,τ3為制動力提升至最大值所需時間;amax為車輛所能達到的最大制動減速度;V1為車輛剛達到最大減速度amax時的行駛速度。
考慮到冰雪天氣的環(huán)境復雜性以及路況的不確定性,取最小安全距離LS=10 m[12],構(gòu)建冰雪天氣下車輛安全停車視距模型:
Ssafe=Smin+LS=
(2)
通常,車輛最大減速度amax與路面摩擦系數(shù)有直接關(guān)系[13]。冰雪路面是冰雪天氣對行車安全造成影響的最主要因素。在冰雪路面條件下,由于路面覆蓋有積雪結(jié)冰層,車輪無法像正常天氣下與路面形成良好接觸,路面附著系數(shù)會減小,車輛制動性能受到制約。文獻[14]在封閉車道環(huán)境下對典型瀝青混凝土路面處于不同的冰層厚度以及不同的冰層表面溫度條件下進行路面附著系數(shù)測定,并通過非線性回歸分析對兩參數(shù)與附著系數(shù)的定量關(guān)系進行擬合,得出冰雪天氣下路面摩擦系數(shù)與冰層表面溫度和冰層表面厚度的關(guān)系:
φ=[(0.006 1-0.003 4h)T+0.042 1h-0.751]N+
[(0.014 4T+0.061 3)h-0.008 7+0.494 2]Ce
(3)
式中,T為冰層表面溫度;h為冰層厚度;N為實驗測量周期次數(shù);Ce為環(huán)境因素折算系數(shù)。
注:實際中,對車輛的制動過程視作產(chǎn)生一次相對摩擦,取N=1;取折算因子Ce=0.37[15]。
駕駛員反應時間主要取決于前方障礙物的外形、顏色、司機的視力和機警程度以及氣象能見度等因素。冰雪天氣條件下,駕駛員的生理及心理特征都會受到不同程度的影響,使駕駛員產(chǎn)生緊張的情緒,導致其反應遲鈍,判斷決策能力下降。如果駕駛員在駕駛車輛的時候長時間處于緊張狀態(tài),就會使駕駛員產(chǎn)生嚴重的視覺疲勞感,有關(guān)研究顯示,大部分駕駛員在進入低能見度環(huán)境上時都會產(chǎn)生緊張不安的情緒;而在低能見度條件下駕駛疲勞的速度會顯著升高。此時能見度與降雪強度并非強線性對應關(guān)系,而是在一定范圍內(nèi)的動態(tài)變化,并且會受到雪粒子的大小、速度、密度、光譜特性等特征的影響[16],容易產(chǎn)生錯誤判斷及視覺疲勞,增加了行車危險性。鑒于冰雪天氣下行車環(huán)境以灰白背景色為主,本文取能見度與可視距離比值為2.5[17-18]。綜合上述冰雪天氣下駕駛員特性的相關(guān)研究,并結(jié)合式(2)~(3),由摩擦附著系數(shù)和安全停車視距反算得出安全行駛車速,得出不同附著系數(shù)和能見度條件下安全限速值計算公式為:
(4)
式中,Vsafe為車輛安全限速值;L為當前路段的氣象能見度;φ為路面摩擦附著系數(shù)。
準確預測交通流的演變趨勢是使用可變限速控制進行科學管控的先決條件。METANET 模型的基本思想是將交通流進行時空維度上的離散化處理,得到表示流量、速度、密度三參數(shù)間關(guān)系的模型。其對于交通流的描述采用宏觀的二階方程,便于嵌入可變限速控制參數(shù),參數(shù)標定過程直觀且準確,可以模擬交通流參數(shù)間的非線性關(guān)系,對于可變限速控制的作用機理的描述更加準確,可從時間和空間維度對交通流演變情況進行周期性預測,進一步滿足可變限速管控系統(tǒng)的準確性與實時性需求,為可變限速管控提供決策參考。
在此,以METANET交通流模型對限速策略下的交通流演化情況進行預測,作為可變限速控制系統(tǒng)的交通流預測反饋。METANET基本交通流模型描述如下:
(5)
其中:第一項為流量速度密度方程,k時刻路段i的交通流量等于同時段和路段的平均速度、密度和車道數(shù)相乘;第二項中,k+1時刻的第i路段的密度是由k時刻的密度,以及來自于i-1路段的車流輸入和離開i路段的車流輸出計算得出;第三項意義是i路段在k+1時刻的空間平均速度分別受本路段期望速度、上游速度、下游密度三項因素的影響;第四項中的V(ρi(k))是駕駛員在路段i行駛時的期望速度。vf,i為i路段的自由流速度,ρcr,i為i路段的關(guān)鍵密度,τ、am、η、θ為模型常量參數(shù),由具體的道路模型進行擬合標定。
在可變限速控制下,期望速度V(ρi(k))不但受到當前路段的車流密度影響,也受到系統(tǒng)發(fā)布的限速值限制。由此,V(ρi(k))的計算公式更正為[19]:
(6)
式中,Vctrl,i為當前時刻i路段由控制系統(tǒng)確定的限速值。
通過元胞自動機模型輸入初始交通流得到管控路段在各控制周期內(nèi)的交通流數(shù)據(jù),同時對預測模型輸入相同的初始交通流,并建立相應的目標函數(shù):
(7)
通過尋優(yōu)算法,求解目標函數(shù)JM的最小值,求解得到各參數(shù)組成的參數(shù)組即為與微觀模型相匹配的METANET模型參數(shù)。
針對多目標優(yōu)化問題,可以通過分別給出各個控制目標在最終優(yōu)化結(jié)果中的優(yōu)先級,通過權(quán)重分配的方式,將多個目標優(yōu)化函數(shù)通過線性加權(quán)求和的方式轉(zhuǎn)化成較簡單的單目標優(yōu)化計算,從而給出在考慮各控制目標權(quán)重情況下的最優(yōu)解。冰雪天氣條件下的可變限速管控策略主要針對行車安全和通行效率兩個控制目標進行優(yōu)化,可通過分別設置關(guān)于行車安全的優(yōu)化目標函數(shù)和關(guān)于通行效率的優(yōu)化目標函數(shù),并進行合理的權(quán)重分配對各優(yōu)化目標函數(shù)進行線性加權(quán)整合,得出作用于系統(tǒng)的優(yōu)化目標函數(shù)。
當前,許多研究表明提升道路通行能力等同于提升交通效率。根據(jù)對交通流通行能力突變現(xiàn)象的研究可知,將管控路段的交通密度維持在關(guān)鍵密度以下,是提升通行能力的關(guān)鍵所在[20]。由此,提出提升交通效率的目標函數(shù):
(8)
式中,ρi(k)和ρi(k+1)分別為第k時刻第i段管控道路的車流平均密度,ρcr為路段關(guān)鍵密度,N為管控路段數(shù),σ為常量參數(shù)。
在元胞自動機模型中,交通流速度的離散性表現(xiàn)在相鄰元胞之間的速度變化,引申到宏觀層面,減小相鄰路段平均速度在時間及空間維度的離散程度,可進一步減小碰撞追尾事故發(fā)生概率,據(jù)此構(gòu)建提升交通安全性的目標函數(shù)[21]:
(9)
式中,vi(k)為k時刻第i段管控道路的車流平均速度,Δvi,i-1(k)為k時刻第i段與第i-1管控道路的車流平均速度絕對差,Ci,i-1(k)為k時刻第i段與第i-1管控道路的通行流量之比。
將減小相鄰路段平均速度的離散性作為可變限速控制的安全性控制目標,將交通密度維持在關(guān)鍵密度附近作為可變限速控制的效率性控制目標,以各限速區(qū)段的限速值作為決策變量,建立冰雪天氣下高速公路可變限速控制多目標優(yōu)化模型,構(gòu)建如下優(yōu)化控制目標函數(shù):
(10)
式中,α、β分別為效率目標函數(shù)和安全目標函數(shù)的權(quán)重系數(shù)。
為保障限速策略的可行性與合理性,需針對道路交通環(huán)境并結(jié)合駕駛員特性,對優(yōu)化目標函數(shù)給出相關(guān)的約束條件。
2.2.1 安全車速約束
管控路段i在第n個管控周期內(nèi)的限速值應當不大于當前路段的固有設計限速值以及當前路況下基于安全視距所確定的安全限速值,且不小于冰雪天氣高速公路最小行車限速:
Vc(i,n)≤min{Vsafe(i,n),Vd}
(11)
式中,Vc(i,n)為第n個管控周期內(nèi)第i個路段的給定限速值;Vsafe(i,n)為第n個管控周期內(nèi)第i路段的最大安全限速值;Vd為當前高速公路設計限速值。
2.2.2 交通流穩(wěn)定性約束
通過對交通流穩(wěn)定性的相關(guān)約束,可保障交通流的穩(wěn)定性以及車輛行駛的安全性,抑制交通運動波的大量產(chǎn)生。一般情況下,相鄰限速段的限速值差值不宜過大,對同一管控路段而言,相鄰限速周期的限速值差值也應當控制在一定的范圍內(nèi),以滿足駕駛員對限速值變化跨度的心理接受范圍:
(12)
2.2.3 限速值合理性約束
考慮到駕駛員的視覺識別特性以及遵從度特征,限速值選取應作相應的取整處理,本文將限速值取為5的整倍數(shù):
Vc(i,n)=Vc(i,n)-Vc(i,n)mod5
(13)
式中,mod為取整函數(shù)。
本文使用粒子群優(yōu)化算法對可變限速控制策略的優(yōu)化目標函數(shù)進行求解。粒子群優(yōu)化算法 (PSO,particle swarm optimization)是由Eberhart和Kennedy于1995年提出的一種群體智能優(yōu)化算法[22],主要研究具有適應性的主體在與環(huán)境和其他主體的交互中進行學習行為,從而改變自身結(jié)構(gòu)或行為的現(xiàn)象。其靈感來自于自然界中的鳥群集體覓食行為,種群中的鳥類個體抽象為沒有質(zhì)量的粒子,分別代表潛在的可行解,而食物的位置則被認為是全局最優(yōu)解。通過模擬鳥群覓食的過程指引粒子在復雜的數(shù)據(jù)空間中逐步迭代搜索最優(yōu)解。與傳統(tǒng)的遺傳變異算法相比,粒子群優(yōu)化算法不需要使用遺傳變異算法的操作符,即交叉、變異和選擇等。此外,粒子群優(yōu)化算法可以處理一些傳統(tǒng)方法不能處理的不可導節(jié)點傳遞函數(shù)的極值點搜索問題,在解決實際工程應用問題中展現(xiàn)了其優(yōu)越性。粒子群優(yōu)化算法的收斂速度較快,算法參數(shù)設置靈活,算法結(jié)構(gòu)簡單清晰,易于工程化實現(xiàn),目前已廣泛應用于諸多的連續(xù)多目標優(yōu)化模型問題的求解[23]。
初始階段,先對預先設置好種群規(guī)模的粒子種群中的每個粒子在搜索空間中賦予隨機的初始速度和初始位置,粒子的位置代表了所求問題的可行解,通過粒子個體的迭代更新行為以及種群內(nèi)的信息交互實現(xiàn)優(yōu)化問題的求解。關(guān)于種群規(guī)模的選擇方面,當粒子數(shù)越多時,迭代收斂至最優(yōu)解位置附近的幾率越高,但與此同時每一步迭代的計算成本增高,容易造成求解速度的降低,因此需考慮迭代效率與求解速度的均衡,根據(jù)求解問題的復雜程度與實際需求確定合適的種群規(guī)模。粒子在搜索更新的過程中,依據(jù)以下3個方面的因素更新自己的速度:首先為“慣性”,代表粒子具有保持自身速度的趨勢;其次為“經(jīng)驗”,代表粒子搜索過程中具有飛往自身搜索到的最佳位置的傾向;最后為“共享”,代表粒子搜索過程中會與群體中其他個體互通信息,具有飛往種群搜索到的最佳位置的傾向,之后再根據(jù)更新后的速度信息對各粒子的位置信息進行逐個更新,隨著迭代次數(shù)的增加,整個粒子群從無序狀態(tài)向有序狀態(tài)過渡,從而實現(xiàn)向最優(yōu)解位置的靠近甚至成功搜索到最優(yōu)解位置[24]。
在n維搜索空間中,每個粒子擁有兩個屬性,即個體當前位置Xi=[xi1,xi2,…,xin]和粒子個體速度Vi=[vi1,vi2,…,vin],同時可獲取記錄中的個體最佳位置pbest和群體最佳位置gbest。其中個體當前位置是指每個粒子在本次迭代完成前所處的位置,個體最佳位置是指該粒子在搜索迭代過程中所達到的歷史極值,群體最佳位置是所有粒子在搜索迭代過程中所達到的歷史極值,在粒子群優(yōu)化算法中一般以最佳適應度表示。粒子個體速度則是根據(jù)上述3個屬性對粒子每次迭代所行進的方向和距離進行調(diào)節(jié)的屬性。在每一次的迭代中,粒子通過跟蹤粒子自身所找到的最優(yōu)解位置pbest以及整個種群目前找到的最優(yōu)解位置gbest,并在迭代中以上述兩極值在向量空間的位置為導向?qū)ψ陨淼奈恢脁i以及自身速度vi進行逐次更新,速度和位置信息更新公式分別如式(14)和式(15)所示:
(14)
xid(t+1)=xid(t)+vid(t+1)
(15)
由于基本粒子群算法容易陷入局部最優(yōu)的困境,學者們對粒子群算法進行了各種優(yōu)化改進,以期實現(xiàn)各種實際優(yōu)化問題的求解。對于基本粒子群優(yōu)化算法的優(yōu)化改進主要集中于對最大速度、慣性權(quán)重參數(shù)、學習因子進行改進、優(yōu)化變異策略以及優(yōu)化編碼方式等方面[25],本文的控制目標為行車安全以及通行效率的優(yōu)化,對算法輸出的準確性和實時性有著一定需求。通過限制粒子的最大速度來限制粒子群的最大運動范圍,防止粒子群因速度過快從最優(yōu)區(qū)域飛出,實現(xiàn)對粒子群優(yōu)化算法搜索準確性的優(yōu)化。另一方面,通過設置變化的慣性權(quán)重參數(shù)可以控制粒子速度的更新快慢,使粒子群優(yōu)化算法在運行初期能夠?qū)θ诌M行充分探索,后期能夠滿足局部重點區(qū)域搜索的需求,從而實現(xiàn)對粒子群算法搜索效率的優(yōu)化。
當慣性權(quán)重ω以及粒子最大速度限制vmax隨著迭代次數(shù)作動態(tài)更新時,粒子群優(yōu)化算法的性能會得到一定的提高。為求解得到使目標函數(shù)達到最優(yōu)值且滿足相應約束條件的決策變量,選擇針對粒子群算法中粒子的最大速度進行動態(tài)限制以及設置遞減的慣性權(quán)重參數(shù)的方法分別對粒子群算法的搜索準確性以及搜索效率進行優(yōu)化改進。粒子最大速度限制vmax和慣性權(quán)重參數(shù)ω的計算公式如式(16)和式(17)所示:
(16)
(17)
改進后的粒子群優(yōu)化算法的尋優(yōu)迭代性能可較好地滿足可變限速控制系統(tǒng)的準確性需求。在可變限速控制問題的場景中,各空間維度分別代表各限速路段的限速值,每一個粒子坐標代表一個候選的限速控制策略,根據(jù)優(yōu)化目標函數(shù)通過多次迭代確定最優(yōu)種群適應度,從而尋得最優(yōu)限速控制策略作用于可變限速控制系統(tǒng)。由于可變限速管控系統(tǒng)的實時性需求,需將粒子群優(yōu)化算法的運算時間成本控制在一定的時間范圍內(nèi),本文通過粒子群優(yōu)化算法優(yōu)化目標函數(shù)在可變限速控制場景下最優(yōu)適應度的實際收斂情況來確定優(yōu)化算法的最大迭代次數(shù),粒子群優(yōu)化算法最優(yōu)適應度迭代過程如圖3所示。
圖3 粒子群算法最優(yōu)適應度迭代過程示意圖
由圖3可知,當?shù)螖?shù)超過100次后最優(yōu)適應度的變化趨于穩(wěn)定,為簡化計算并滿足可變限速控制系統(tǒng)的實時性方面需求,將最大迭代次數(shù)設置為100。
基于粒子群優(yōu)化的可變限速控制算法流程如表1所示。
為驗證冰雪天氣高速公路可變限速控制策略的可行性與優(yōu)勢,分別將全線固定限速和分段安全限速兩種控制方法與本文方法進行對比仿真研究。
選取十天高速公路G7011段作為可變限速控制路段,如圖4所示。整個路段道路曲率較小,仿真中近似視作直線路段,自左向右分別將可變限速管控路段劃分為A、B、C、D四段,最左側(cè)為緩沖段,均為等長的1.2 km路段。
圖4 仿真路段示意圖
根據(jù)尋優(yōu)算法,求解目標函數(shù)JM的最小值,此時得到的參數(shù)組即為與微觀模型匹配性最佳的模型參數(shù)。仿真路段宏觀交通流預測模型參數(shù)標定值如表2所示。
表2 交通流預測模型標定參數(shù)組
仿真時長為120 min,其中前20 min為仿真預熱期,后100 min路況狀態(tài)設定列于表3,限速值更新周期為20 min,車輛速度隨機慢化概率設為0.3,粒子種群規(guī)模設置為40,粒子個體學習因子c1和粒子種群學習因子分別設置為0.8和0.9,初始粒子最大速度限制設為2.0,優(yōu)化目標函數(shù)權(quán)重系數(shù)α、β分別選為3.0和1.0。
表3 仿真路段路況狀態(tài)
根據(jù)仿真條件與上述模型計算所得可變限速方案列于表4,仿真結(jié)果如圖5~7所示。
表4 可變限速控制方案
圖5 不同限速方案下整體路段通行能力
圖6 不同限速方案下整體路段交通密度
圖7 不同限速方案下B路段車輛速度標準差
由圖5~6可知,對比其它兩種限速方案,本文提出的可變限速控制方案提升了道路通行能力,并有效降低了車流密度。由圖7可知,相比分段安全限速方案,本文提出的可變限速控制方案降低了車輛速度標準差,優(yōu)化了整體交通流的穩(wěn)定性,使行車安全性得到提升。
此外,采用總行程時間TTS(total time spend)、總行程延誤時間TTD(total time delay)以及車輛沖突暴露時間ETC(exposed time to collision)三項指標,分別從交通效率與安全兩個方面對該可變控制策略的優(yōu)化效果進行對比分析,三項交通指標的仿真結(jié)果如表5所示。
由表5可知,相比于全線固定限速方法,可變限速控制策略下的總行程時間與總行程延誤時間分別減少了20.5%和30.1%,車輛沖突暴露時間降低8.6%;相比于分段安全限速方法,可變限速控制策略下的總行程延誤時間減少25.5%,車輛沖突暴露時間降低了12.6%。
表5 不同限速控制策略的交通指標對比
由以上仿真結(jié)果可以看出,本文所提出的可變限速方案在提升道路通行效率的同時,通過對交通流的宏觀調(diào)控,有效降低了車輛碰撞風險,改善了行車安全性。
針對冰雪條件下高速公路車輛通行效率及安全性降低的問題,在建立冰雪天氣的高速公路安全限速模型以及交通流預測模型的基礎(chǔ)上,設計了兼顧通行效率與行車安全的優(yōu)化目標函數(shù)與相關(guān)約束條件,并通過粒子群優(yōu)化算法對可變限速值進行求解,從而提出一種兼顧行車安全與通行效率的高速公路可變限速策略。仿真結(jié)果表明,所提出的可變限速方法與固定限速方法以及分段安全限速方法相比,能有效地提升行車安全和通行效率。