寇 潔, 李向有, 張雅楠, 胡 婷, 郭熊果
(延安大學(xué) 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院, 延安 716000)
隨著延安市城市化規(guī)模的不斷擴(kuò)大和人口數(shù)量的增加,延安市的交通壓力越來(lái)越大,交通擁堵問(wèn)題日益突出,并且逐漸成為制約城市發(fā)展和居民生活質(zhì)量的重要因素[1]. 交通擁堵在一定程度上會(huì)增加市民的通勤成本和時(shí)間成本,增加交通事故的發(fā)生率,還會(huì)導(dǎo)致城市環(huán)境惡化. 這些問(wèn)題既影響了市民的出行,也對(duì)城市的經(jīng)濟(jì)發(fā)展和社會(huì)穩(wěn)定造成了一定的影響[2]. 因此為改善城市交通可持續(xù)發(fā)展和提高居民生活質(zhì)量,本研究將對(duì)緩解延安市交通擁堵問(wèn)題提出優(yōu)化策略.
本研究是調(diào)查延安市老城區(qū)主要路段在特定時(shí)段的交通擁堵問(wèn)題,選取延安市寶塔區(qū)楊家?guī)X路段交叉口為研究對(duì)象,調(diào)查其相關(guān)交通數(shù)據(jù)、十字交叉口交匯點(diǎn)的通行狀況、單雙行道的分布及綠信比等,分析造成交通擁堵的各種因素,并通過(guò)設(shè)計(jì)優(yōu)化策略,調(diào)整信號(hào)配時(shí)方案,緩解交叉口交通擁堵問(wèn)題. 最后利用VISSIM仿真軟件對(duì)比兩相位交叉口和多相位交叉口優(yōu)化前后的時(shí)延、隊(duì)列長(zhǎng)度等擁塞指標(biāo),展示方法優(yōu)化效果. 基于以上分析,研究延安市交通擁堵問(wèn)題的優(yōu)化策略,具有重要的現(xiàn)實(shí)意義和理論價(jià)值.
研究目的是緩解延安市交通擁堵問(wèn)題,研究對(duì)象是由楊家?guī)X隧道、楊家?guī)X大橋和延大路組成的交叉口處的交通情況,對(duì)交叉口的現(xiàn)狀畫(huà)CAD圖,如圖1所示.
圖1 楊家?guī)X路段交叉口原狀
楊家?guī)X大橋的交叉口進(jìn)口道為2條,有左轉(zhuǎn)和直行,無(wú)右轉(zhuǎn),出口道2條;楊家?guī)X隧道的交叉口進(jìn)口道3條,有右轉(zhuǎn)和直行,無(wú)左轉(zhuǎn),出口道3條;延大路由西北向東南方向的進(jìn)口道1條,為公交車專用道,只有直行,出口道2條;延大路由東南向西北方向的進(jìn)口道2條,有左轉(zhuǎn)、右轉(zhuǎn)和直行,出口道1條. 圣地路、楊家?guī)X大橋、楊家?guī)X隧道均無(wú)非機(jī)動(dòng)車道. 楊家?guī)X大橋的中央隔離設(shè)施為護(hù)欄;圣地路的中央隔離設(shè)施為黃實(shí)線隔離;楊家?guī)X隧道的中央隔離設(shè)施為綠化帶隔離.
楊家?guī)X隧道是連接新區(qū)和老區(qū)的主要道路,由于其特殊的地理位置,經(jīng)常會(huì)出現(xiàn)交通擁堵的情況. 具體表現(xiàn)為早、晚高峰時(shí)段車輛排隊(duì)等待進(jìn)出隧道,等待時(shí)間長(zhǎng),擁堵嚴(yán)重;隧道內(nèi)車流量大,車輛行駛緩慢,擁堵嚴(yán)重;道路交叉口信號(hào)燈配時(shí)不太合理,隧道兩側(cè)車輛需排隊(duì)等待進(jìn)入隧道,影響周邊道路的交通暢通. 擁堵情況在周五晚高峰時(shí)段和下雨天晚高峰時(shí)段頗為嚴(yán)重. 由于地形地勢(shì)所限,在短時(shí)間內(nèi)通過(guò)加寬路面和減少該路段的車流量不具有可行性,因此優(yōu)化該道路交叉口信號(hào)燈配時(shí),成為解決該交叉口擁堵的方案之一. 同時(shí),為驗(yàn)證本文提出優(yōu)化方法的適用性,本文分別采用楊家?guī)X路段交叉口和1個(gè)四相位交叉口為研究對(duì)象,證明該優(yōu)化方法在兩相位及多相位交叉口同樣具有優(yōu)化效果.
門(mén)控循環(huán)單元(Gated Recurrent Unit,GRU)為1種用于處理序列數(shù)據(jù)的循環(huán)神經(jīng)網(wǎng)絡(luò)(Recurrent Neural Network,RNN)的變體[1]. 于2014年由Kyunghyun Cho等提出,旨在克服傳統(tǒng)RNN面臨的梯度消失問(wèn)題和長(zhǎng)依賴問(wèn)題.
楊家?guī)X路段交叉口的交通流量是帶有時(shí)間序列的數(shù)據(jù),數(shù)據(jù)前后有較強(qiáng)的關(guān)聯(lián)性. 使用神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)時(shí)間序列數(shù)據(jù)表現(xiàn)出的性能較好,常用于時(shí)間序列預(yù)測(cè)的神經(jīng)網(wǎng)絡(luò)模型有RNN、LSTM、GRU. RNN循環(huán)神經(jīng)網(wǎng)絡(luò)模型可用于解決時(shí)間序列預(yù)測(cè)問(wèn)題,網(wǎng)絡(luò)的隱藏狀態(tài)可捕獲到當(dāng)前時(shí)間步長(zhǎng)的序列歷史,并且模型參量不隨時(shí)間步長(zhǎng)的增加而增長(zhǎng). 但如果預(yù)測(cè)較長(zhǎng)的時(shí)間序列數(shù)據(jù),RNN會(huì)產(chǎn)生梯度爆炸和梯度消失的問(wèn)題[4]. 為此,有學(xué)者提出將RNN隱藏層神經(jīng)元改造成記憶神經(jīng)元,引入門(mén)控循環(huán)網(wǎng)絡(luò),緩解RNN梯度消失帶來(lái)的問(wèn)題,控制信息流動(dòng),使模型更好地記住長(zhǎng)遠(yuǎn)時(shí)期的信息. 長(zhǎng)短期記憶網(wǎng)絡(luò)(LSTM)模型中有遺忘門(mén)、輸出門(mén)和輸入門(mén),可有效防止重要信息在數(shù)據(jù)訓(xùn)練時(shí)消失. GRU模型是改造了LSTM模型的記憶神經(jīng)元形成的,計(jì)算收斂速度更快,并且兩個(gè)算法預(yù)測(cè)精度相差不大,GRU模型較LSTM模型簡(jiǎn)單,其結(jié)構(gòu)如圖2所示[4]. 所以本研究選擇GRU模型來(lái)預(yù)測(cè)交通流量,并結(jié)合粒子群算法優(yōu)化信號(hào)配時(shí),從而建立基于多目標(biāo)粒子群算法-GRU循環(huán)神經(jīng)網(wǎng)絡(luò)的信號(hào)配時(shí)優(yōu)化模型.
圖2 GRU網(wǎng)絡(luò)模型的架構(gòu)
多目標(biāo)粒子群優(yōu)化算法(Multi-Objective Particle Swarm Optimization,MOPSO)為1種基于粒子群優(yōu)化算法(Particle Swarm Optimization,PSO)的智能優(yōu)化算法,廣泛應(yīng)用于解決具有多個(gè)優(yōu)化目標(biāo)的問(wèn)題[6-8]. MOPSO算法的思想來(lái)源于對(duì)鳥(niǎo)群捕食行為的研究,該算法結(jié)構(gòu)簡(jiǎn)單,操作起來(lái)容易實(shí)現(xiàn),并且具有良好的全局搜索能力. 然而,MOPSO算法也存在收斂速度慢、容易陷入局部最優(yōu)等問(wèn)題,并不能保證找到全局最優(yōu)解,因此在具體應(yīng)用中需要根據(jù)實(shí)際情況進(jìn)行參數(shù)設(shè)置和改進(jìn). 設(shè)計(jì)MOPSO算法的關(guān)鍵在于粒子的速度和位置的更替,進(jìn)而粒子根據(jù)個(gè)體經(jīng)驗(yàn)和群體經(jīng)驗(yàn)的權(quán)重對(duì)速度和位置進(jìn)行調(diào)整,使得粒子能在搜索空間中進(jìn)行探索和利用信息,從而逐步尋找到最優(yōu)解.
多目標(biāo)粒子群優(yōu)化算法的實(shí)現(xiàn)流程如下:
算法名稱:多目標(biāo)粒子群優(yōu)化算法
算法功能:通過(guò)模擬群體行為實(shí)現(xiàn)全局搜索和解空間探索,以找到最優(yōu)解.
輸入:參數(shù)w:0.5~0.8;c1,c2:0.1~2;vmax,xmax:取決于優(yōu)化函數(shù).
輸出:輸出最優(yōu)值.
過(guò)程:
步驟1:初始化粒子的速度和位置,設(shè)置粒子數(shù)量、迭代次數(shù)等參數(shù);
步驟2:計(jì)算個(gè)體適應(yīng)度;
步驟3:更替粒子速度→更替粒子位置;
步驟4:計(jì)算更新位置后的適應(yīng)度,若新位置的適應(yīng)度更高,則將該粒子的位置進(jìn)行更替,否則不更替.
步驟5:判斷是否滿足終止條件,如果滿足則退出,不滿足則返回步驟2.
3.1.1 目標(biāo)函數(shù)
本研究選取的配時(shí)優(yōu)化模型評(píng)價(jià)指標(biāo)有延誤(h),停車次數(shù),通行能力(pcu/h)[9-15]. 由于3個(gè)評(píng)價(jià)指標(biāo)的單位不同,需對(duì)單位進(jìn)行統(tǒng)一,這樣有助于簡(jiǎn)化計(jì)算和建立模型的過(guò)程. 進(jìn)而對(duì)優(yōu)化目標(biāo)分別賦予不一樣的權(quán)重系數(shù)X,Y,Z,建立式(1):
(1)
式中,D0為現(xiàn)狀配時(shí)方案下的車輛延誤;H0為現(xiàn)狀配時(shí)方案下的停車次數(shù);Q0為現(xiàn)狀配時(shí)方案下的通行能力;X,Y,Z代表不同優(yōu)化目標(biāo)的權(quán)重系數(shù),其中X+Y+Z=1,且X>0,Y>0,Z>0.系數(shù)越大,則對(duì)應(yīng)目標(biāo)需進(jìn)行優(yōu)化的程度也越大.
交通流量是個(gè)變量,影響交通流量的因素有很多,它隨時(shí)間的變化而變化.本研究通過(guò)計(jì)算各相位的飽和流量比α,設(shè)置不同權(quán)重系數(shù)的變量,意味著對(duì)于交叉口不同相位的車流量分配給予不同的重視程度或優(yōu)先級(jí).通過(guò)設(shè)置權(quán)重系數(shù),可調(diào)整各個(gè)相位在信號(hào)燈控制中的綠燈時(shí)間分配,以實(shí)現(xiàn)更合理的交通流調(diào)度.當(dāng)交叉口處于低飽和狀態(tài)時(shí),α值相對(duì)較小,在所需優(yōu)化目標(biāo)中,主要控制車輛延誤D和停車次數(shù)H最少;當(dāng)臨近飽和或飽和狀態(tài)時(shí),α值相對(duì)較大,在所需優(yōu)化目標(biāo)中,主要控制通行能力Q最大.因此,根據(jù)不同的交通優(yōu)化需求目標(biāo),調(diào)節(jié)3個(gè)優(yōu)化目標(biāo)不同的權(quán)重系數(shù).
因此,根據(jù)上述分析,將各個(gè)權(quán)重系數(shù)X,Y,Z分別設(shè)置為:
(2)
式中,C為交叉口的信號(hào)周期,α為各相位中最大飽和流量比.
則目標(biāo)函數(shù)建立式(3):
(3)
式中,gi為相位i的有效綠燈時(shí)間;C為交叉口的信號(hào)周期;L為損失時(shí)間.
3.1.2 約束條件
信號(hào)燈在交叉口的實(shí)際運(yùn)用中,信號(hào)周期C和綠信比λ都需要做一定的限制,不能設(shè)置的過(guò)大或者過(guò)小.一方面,信號(hào)周期時(shí)長(zhǎng)過(guò)長(zhǎng)時(shí),在較長(zhǎng)的綠燈時(shí)間的后期,車流較小,紅燈方向車輛等待時(shí)間過(guò)長(zhǎng),導(dǎo)致道路交叉口沖突區(qū)域使用效率較低,造成道路交叉口通行能力降低;另一方面,信號(hào)周期時(shí)長(zhǎng)較短時(shí),紅綠燈反復(fù)切換,導(dǎo)致車輛頻繁停車等待,損失時(shí)間增加,車輛多次啟動(dòng)、停車,致使通行能力降低.因此合理的周期時(shí)長(zhǎng)一般為24~180 s,因此最高值為180 s,最低值為24 s.最小周期時(shí)長(zhǎng)確定見(jiàn)式(4):
(4)
式中,Cmin為最短周期時(shí)長(zhǎng)(s);Y為全部相位的最大流量比之和;L為周期損失時(shí)間見(jiàn)(5):
(5)
式中,n為相位總數(shù);li為第i個(gè)相位的損失時(shí)間.
各相位綠信比(1個(gè)相位內(nèi)某一方向有效通行時(shí)間與周期長(zhǎng)度之比)之和為定值見(jiàn)式(6):
(6)
式中,C為周期時(shí)長(zhǎng);G為周期綠燈總時(shí)間;Ye為周期黃燈總時(shí)間;L為周期總損失時(shí)間.
建立MOPSO-GRU算法的目標(biāo)是為了設(shè)置較為合適的信號(hào)配時(shí),模型中的GRU負(fù)責(zé)預(yù)測(cè)交通流量. 本研究采用人為拍攝道路交叉口交通情況進(jìn)行記錄車流量數(shù)據(jù),將車流量按照5 min為間隔建立數(shù)據(jù)集,數(shù)據(jù)集分成訓(xùn)練集和測(cè)試集,使用GRU模型對(duì)交通流量進(jìn)行預(yù)測(cè),將均方根誤差(RMSE)、平均絕對(duì)誤差(MAE)、準(zhǔn)確率(ACC)作為模型預(yù)測(cè)精度評(píng)價(jià)指標(biāo)[14]. 結(jié)果如圖3所示.
圖3 GRU模型預(yù)測(cè)某道口4月7日車流量結(jié)果展示
為驗(yàn)證本文提出的信號(hào)配時(shí)優(yōu)化方法在兩相位交叉口和多相位交叉口均可提高交叉口運(yùn)行效率,本文選取了1個(gè)兩相位交叉口和1個(gè)四相位交叉口作為實(shí)驗(yàn)對(duì)象.
3.2.1 兩相位交叉口
選取GRU測(cè)試集預(yù)測(cè)值的高峰小時(shí)交通量作為粒子群算法調(diào)優(yōu)的輸入數(shù)據(jù),兩相位交叉口其數(shù)據(jù)展示如表1所示.
楊家?guī)X路段交叉口是兩相位交叉口,交叉口信號(hào)控制方案如圖4所示.
3.2.2 四相位交叉口
為驗(yàn)證本文所提出信號(hào)配時(shí)方法在多相位交叉口上的適用性,本文采用了1個(gè)四相位交叉口作為研究對(duì)象,該交叉口東西走向?yàn)槠桨猜?南北走向?yàn)槭兰o(jì)路,其數(shù)據(jù)展示如表2所示.
表1 兩相位交叉口高峰小時(shí)交通量
圖4 兩相位交叉口相位圖
圖5 四相位交叉口相位圖
平安路—世紀(jì)路交叉口是兩相位交叉口,交叉口信號(hào)控制方案如圖5所示.
3.2.3 實(shí)驗(yàn)結(jié)果
以兩相位交叉口為例,利用MOPSO算法進(jìn)行實(shí)驗(yàn),相關(guān)參數(shù)設(shè)置如下:黃燈時(shí)間及各相位損失時(shí)間都為3 s;粒子群初始數(shù)量為100,迭代次數(shù)設(shè)為100次,加速常數(shù)C1和C2設(shè)為1.5;最大有效綠燈時(shí)間、最小有效綠燈時(shí)間分別為90 s、10 s,最大周期時(shí)長(zhǎng)設(shè)為180 s.
粒子分布的初始狀態(tài)如圖6(a)所示,迭代20、50、100次下,粒子的分布狀態(tài)分別如6(b)~(d)所示.
圖6 粒子迭代示意圖
圖6(a)為初始狀態(tài)下粒子位置圖,圖中的粒子較為分散,紅色的點(diǎn)代表初始化時(shí)最優(yōu)解的位置. 1張圖中有4小幅圖,前3幅圖是粒子變化的二維圖,第4幅圖是粒子變化的三維圖. 根據(jù)圖 7粒子迭代過(guò)程,可看出:隨著迭代次數(shù)的增加,結(jié)果趨于穩(wěn)定狀態(tài). 對(duì)于楊家?guī)X路段交叉口的情形,預(yù)設(shè)信號(hào)配時(shí)優(yōu)化方法,迭代60次后,粒子群逐步向最優(yōu)解靠近. 3項(xiàng)評(píng)價(jià)指標(biāo)也逐漸趨于穩(wěn)定達(dá)到最優(yōu)解.
圖7 迭代過(guò)程
以兩相位交叉口為實(shí)驗(yàn)對(duì)象,程序運(yùn)行結(jié)果為:相位Ⅰ綠燈時(shí)間為76 s,相位Ⅱ綠燈時(shí)間為66 s,信號(hào)周期為148 s.
以四相位交叉口為實(shí)驗(yàn)對(duì)象,用同樣的方式優(yōu)化四相位交叉口,程序運(yùn)行結(jié)果為:相位Ⅰ綠燈時(shí)間為48 s,相位Ⅱ綠燈時(shí)間為18 s,相位Ⅲ綠燈時(shí)間為46 s,相位Ⅳ綠燈時(shí)間為23 s,信號(hào)周期為150 s.
本研究首先對(duì)楊家?guī)X路段交叉口的交通組織進(jìn)行優(yōu)化,其次優(yōu)化信號(hào)配時(shí),不優(yōu)化交叉口的相位,最后完成交叉口路網(wǎng)的繪制和相關(guān)參數(shù)設(shè)定,對(duì)優(yōu)化前后2個(gè)道路交叉口進(jìn)行仿真[17-18]. 仿真路網(wǎng)中車道數(shù)按楊家?guī)X路段交叉口的渠化條件設(shè)置,車道寬度均設(shè)置為3.5 m. 楊家?guī)X隧道至楊家?guī)X大橋的飽和交通流量為1 600 pcu/h,延大路東南方向至西北方向的飽和交通流量為1 464 pcu/h. 車輛構(gòu)成和車速等均采用VISSIM的默認(rèn)值. 黃燈時(shí)長(zhǎng)、各相位損失時(shí)間均設(shè)置為3 s. 信號(hào)配時(shí)設(shè)置時(shí)間由優(yōu)化前后配時(shí)方案的結(jié)果決定,圖8(a)為兩相位交叉口原信號(hào)配時(shí),圖8(b)為粒子群算法優(yōu)化后的兩相位交叉口信號(hào)配時(shí). 最后運(yùn)行仿真模型得到仿真后各評(píng)價(jià)指標(biāo)的值.
為驗(yàn)證本文所提出信號(hào)配時(shí)方法在多相位交叉口上的適用性,也對(duì)四相位交叉口進(jìn)行了仿真驗(yàn)證,圖9(a)為四相位交叉口原信號(hào)配時(shí),圖9(b)為粒子群算法優(yōu)化后的四相位交叉口信號(hào)配時(shí).
圖8 兩相位交叉口優(yōu)化前后信號(hào)配時(shí)對(duì)比
圖9 四相位交叉口優(yōu)化前后信號(hào)配時(shí)對(duì)比圖
對(duì)優(yōu)化前后兩相位交叉口和四相位交叉口分別進(jìn)行仿真,仿真持續(xù)3 600 s,仿真結(jié)果對(duì)比如圖 10所示.
圖10 優(yōu)化前后評(píng)價(jià)指標(biāo)結(jié)果對(duì)比
根據(jù)圖8顯示,楊家?guī)X路段交叉口在優(yōu)化前的信號(hào)周期時(shí)長(zhǎng)為110 s,而優(yōu)化后的周期時(shí)長(zhǎng)為140 s,導(dǎo)致周期時(shí)間增加的原因是:在優(yōu)化信號(hào)配時(shí)過(guò)程中,為了實(shí)現(xiàn)交通流的高效運(yùn)行,會(huì)增加某些進(jìn)口道的綠燈時(shí)間,以便更多車輛通過(guò). 這樣一來(lái),路口的總車流量會(huì)增加,從而導(dǎo)致信號(hào)周期時(shí)間增長(zhǎng).
楊家?guī)X路段交叉口在優(yōu)化前的綠燈信號(hào)時(shí)長(zhǎng)為70 s,優(yōu)化后的綠燈信號(hào)時(shí)長(zhǎng)為76 s,而優(yōu)化之后綠燈配時(shí)沒(méi)有發(fā)生很大變化,原因?yàn)?在優(yōu)化信號(hào)配時(shí)過(guò)程中,優(yōu)化目標(biāo)與綠燈時(shí)長(zhǎng)的變化之間存在一定的權(quán)衡關(guān)系,因此最終的配時(shí)結(jié)果可能在一定程度上保持穩(wěn)定.
根據(jù)優(yōu)化前后評(píng)價(jià)指標(biāo)結(jié)果對(duì)比可知,兩相位交叉口在新的配時(shí)方案下平均車輛延誤降低了19.83%、平均停車次數(shù)降低了27.16%、車輛通行能力提高了5.98%,四相位交叉口在新的配時(shí)方案下平均車輛延誤降低了17.29%、平均停車次數(shù)降低了32.04%、車輛通行能力提高了5.81%.
綜上所述,使用本研究提出的信號(hào)配時(shí)優(yōu)化模型,可顯著改善交叉口的交通狀況. 車輛的延誤減少,停車次數(shù)減少,通行能力提高,這意味著車輛能更加順暢地通過(guò)交叉口,交通效率得到了提升.
交通擁堵問(wèn)題是城市發(fā)展過(guò)程中普遍存在的問(wèn)題,緩解交通擁堵問(wèn)題需要綜合運(yùn)用各種措施和方法,找到最優(yōu)解決方案. 本研究首先分析了楊家?guī)X路段交叉口目前存在的問(wèn)題,繼而對(duì)交通組織結(jié)構(gòu)進(jìn)行優(yōu)化. 之后構(gòu)建了GRU模型對(duì)交通流量進(jìn)行預(yù)測(cè),預(yù)測(cè)結(jié)果接近楊家?guī)X路段交叉口車流量真實(shí)值,誤差較小. 然后引入MOPSO算法作為解決方案,本研究將MOPSO和GRU結(jié)合,應(yīng)用于信號(hào)配時(shí)優(yōu)化問(wèn)題,將交通流量預(yù)測(cè)數(shù)據(jù)的最大值作為粒子群算法調(diào)優(yōu)的輸入數(shù)據(jù),得到優(yōu)化后的信號(hào)配時(shí)方案. 最后搭建VISSIM仿真平臺(tái),分別使用了兩相位和四相位的交叉口為實(shí)驗(yàn)對(duì)象進(jìn)行驗(yàn)證,對(duì)比優(yōu)化前后的實(shí)驗(yàn)數(shù)據(jù),驗(yàn)證MOPSO-GRU模型的合理性. 實(shí)驗(yàn)結(jié)果表明,本研究所提出的基于MOPSO-GRU神經(jīng)網(wǎng)絡(luò)的信號(hào)配時(shí)優(yōu)化模型可較好提高交通流量的效率和整體交通系統(tǒng)的性能.