李茹雪，劉 瀾,2

(1.西南交通大學 交通運輸與物流學院，四川 成都 610031; 2.西南交通大學 綜合交通運輸智能化國家地方聯(lián)合工程實驗室，四川 成都 610031)

0 引言

隨著城鎮(zhèn)化進程的發(fā)展和人們生活水平的提高，私人小汽車的擁有量快速增長，交通擁擠成為制約大城市經(jīng)濟發(fā)展的一個重要問題。同時，人們也逐漸從世界各地解決交通擁堵的實踐中認識到：從交通供給層面新建或擴展道路，在緩解擁堵的同時也產(chǎn)生了新的需求，結果便是導致交通出行環(huán)境的進一步惡化。因此，擁擠收費作為一有效的交通需求管理手段被提出，通過在特定的時段或路段對車輛實行收費的方法，限制車輛的使用頻率，降低用戶出行需求的增長速度，從而維持系統(tǒng)合理的運行和服務水平。

在過去的幾十年里，擁擠收費在倫敦[1]、新加坡[2]、斯德哥爾摩[3]等多個國家或地區(qū)得到了實踐。這些擁擠收費策略大都采用了邊界收費的方法，即對劃定邊界區(qū)域實行收費從而達到目標函數(shù)最優(yōu)。大多數(shù)邊界收費研究采用收費相同的單一費率方法，如 Verhoef建立了仿真模型來研究邊界收費，并通過啟發(fā)式算法求解得到費率和收費點的最佳位置[4]。但隨著研究的深入，單一費率的收費模式逐漸顯現(xiàn)出不公平、造成資源浪費等弊端。2000年，May 和Milne首先提出了分別基于時間、擁擠度和距離的3種多費率邊界收費策略[5]。其中前兩種收費策略通常與其他收費方法結合使用，而基于距離的邊界收費策略因其公平高效的特點得到了眾多學者的關注。在May的基礎之上，早期的學者們對英國[6]和德國[7]進行了基于距離的收費策略研究，但將費用定義為行駛距離的線性函數(shù)的做法與現(xiàn)實不符。隨后，Lawphongpanich和Yin提出了一種非線性的收費策略，即利用分段線性函數(shù)表達非線性距離函數(shù)進而求得費用，但該分段線性方法只能將分段數(shù)量限制在兩段及以下[8]。為了滿足交通網(wǎng)絡對更多分段形式的要求，Meng研究了將非線性距離函數(shù)等分成n段線性函數(shù)來獲取最優(yōu)距離收費的方法，并提出最優(yōu)距離收費函數(shù)是關于距離的非遞減正函數(shù)[9]。國內的孫鑫[10]和程啟秀[11]分別研究了基于距離的靜態(tài)與動態(tài)收費策略，前者對求解非線性距離函數(shù)中基于路段的算法進行了改進，后者則考慮交通流的時變性，研究隨機逐日動態(tài)擁擠收費下的最優(yōu)費率。Zheng首先將宏觀基本圖(Macroscopic Fundamental Diagram)作為判斷最優(yōu)收費的工具來實施邊界收費策略[12]，證明了MFD相比其他方法更為直觀有效，但單一費率的邊界收費忽略了車輛在邊界區(qū)域中行駛的距離。因此為了提高擁擠收費的公平性與有效性，本研究提出將MFD與多費率邊界收費策略相結合的方法，即考慮車輛在邊界區(qū)域中行駛的距離對費率產(chǎn)生的影響，將MFD的相關特性應用于收費模型，得到使邊界收費區(qū)域的通行能力維持在最優(yōu)水平的收費策略。

智能交通技術的迅猛發(fā)展為該收費策略在數(shù)據(jù)的獲取上提供了技術保障。以電子警察系統(tǒng)、治安卡口系統(tǒng)和各類城市交通監(jiān)控系統(tǒng)為代表的智能交通系統(tǒng)建設，基本上實現(xiàn)了在我國大中城市主要路段和交叉口信息采集點位的完整覆蓋[13]；結合GPS和北斗衛(wèi)星導航定位系統(tǒng)、電子地圖和ETC等技術的應用，能夠實現(xiàn)車輛的精確定位，從而獲取車輛在邊界區(qū)域內的行駛距離，再根據(jù)收費函數(shù)即可計算出出行者應支付的費用?，F(xiàn)代信息技術提供的大數(shù)據(jù)環(huán)境，加之數(shù)據(jù)采集、存儲、傳輸和計算等方面軟硬件技術的飛速進步，為在實際交通管控業(yè)務中完備、精細地處理各種復雜因素提供了可能。

為此，本研究首先在第1部分介紹了MFD的基本概念和相關特性，在此基礎上第2部分定義了收費函數(shù)的非線性形式，并結合MFD特性建立了固定需求下的雙層規(guī)劃收費模型，第3部分則提出了一種基于非線性收費函數(shù)的改進FW算法，以求解該收費模型。為了驗證模型和算法的可行性，論文在第4部分中進行了算例分析，并在第5部分總結了研究工作的結論，提出了面向實際應用，進一步深化研究的方向。

1 基于MFD的邊界收費策略的提出

圖1 區(qū)域路網(wǎng)宏觀基本圖Fig.1 Macroscopic fundamental diagram of regional road network

(1)

式中，a和da分別為路段a和該路段的長度；qa為路段a上的流量。

從圖中可看出，當路網(wǎng)內的車輛數(shù)累積到一定數(shù)值時，加權流量到達臨界值，路網(wǎng)整體運行效益最優(yōu)[16]；若車輛數(shù)繼續(xù)增長，則加權流量會隨之降低，此時路網(wǎng)處于過擁擠狀態(tài)?；诖诵再|，許多學者將MFD應用于構建控制模型[17]或經(jīng)濟學模型[9]中。

因此，借鑒相關策略中通過MFD調節(jié)進入邊界的車輛數(shù)這一思想，在實施收費策略時，將MFD的相關參數(shù)與收費模型相結合，讓區(qū)域內的累計車輛數(shù)N盡量保持在臨界車輛數(shù)N*附近，從而使收費區(qū)域的運行效益最大，保證區(qū)域內交通運行暢通。

2 基于距離的邊界收費方案設計

2.1 距離收費函數(shù)

在以往將距離收費函數(shù)設置為線性形式的收費策略中，曲線斜率多是人為主觀設定，導致無法準確獲得系統(tǒng)最優(yōu)目標值。因此，本研究將基于距離的收費函數(shù)定義為非線性函數(shù)：

(2)

(3)

由于該收費函數(shù)是非線性的，造成廣義路徑出行時間不等于該路徑上各條路段出行時間之和，這被稱為路徑成本的不可加性[16]。在此采用分段線性逼近的方法來解決此非線性收費問題。首先用lmax和lmin分別為邊界區(qū)域中最大和最小的路徑長度，對應于收費費率ymin和ymax。將該范圍等分為n個相同的間隔，并假設收費費率是車輛在收費區(qū)域中行駛距離的非遞減函數(shù)，則非線性函數(shù)近似由n個線性函數(shù)組成，且每個線性函數(shù)由間隔的兩個邊界值唯一確定。如圖2所示，當車輛在邊界區(qū)域中的行駛距離值屬于某個間隔范圍內時，就使用該間隔的線性函數(shù)計算費率。

圖2 分段線性收費原理Fig.2 Piecewise-linear tolling principle

因此得到計算基于距離的收費費率公式為：

(4)

2.2 固定用戶平衡下的雙層規(guī)劃模型

用MFD判斷上層目標函數(shù)是否達到最優(yōu)的方法，既簡化了求解雙層規(guī)劃模型過程中的計算，又能根據(jù)收費區(qū)域的MFD圖像，較為直觀清晰地觀測出系統(tǒng)處于最優(yōu)運行狀態(tài)下的臨界值，從而找到使收費區(qū)域輸出流量最大的費率組合y=[y0,y1,…，yk,yk+1,…，yn]Τ。因此，以收費區(qū)域輸出流量最大為目標建立上層模型，目標函數(shù)如下：

(5)

(6)

式中，a和da分別為收費區(qū)域內路段a和該路段的長度;xa為收費區(qū)域內路段a上的流量。

假定交通網(wǎng)絡中用戶的路徑選擇行為服從Wardrop第一定理，即用戶知道每條可選擇路徑的交通時間，并能選擇其中的最短路，此時就達到了用戶平衡狀態(tài)(UE)。在實行收費方案后，用戶在交通網(wǎng)絡內進行路徑選擇時不僅要考慮路徑的出行時間，還要考慮收費所造成的出行成本的增加。利用出行時間價值這一參數(shù)，將收費所造成的出行成本換算為時間費用，得到廣義出行成本如下：

(7)

式中，β為用戶的出行時間價值，路段出行時間ta(xa)采用BPR函數(shù)計算，如下公式所示：

(8)

式中，Ca為路段a的通行能力。在給定的距離收費函數(shù)下，得到的下層用戶均衡模型如下：

(9)

(10)

(11)

(12)

3 基于路段的FW算法

Frank-Wolfe(FW)算法作為一種可行方向法，廣泛運用于求解含線性約束的非線性規(guī)劃問題。該算法主要可分為確定可行下降方向和確定步長兩個步驟，其中確定可行下降方向等價于求解線性規(guī)劃問題，進一步等價于求解最短路問題[9]。但由于非線性距離收費函數(shù)導致了路徑成本的不可加性，使用傳統(tǒng)的最短路算法不再可行。因此Meng等人提出了一種網(wǎng)絡轉換法(Network Reformation Technique)以解決路徑成本的不可加性問題[18]。該方法通過用虛擬路段(Dummy Link)替代收費區(qū)域中的路徑，從而將路徑不可加成本轉換為基于路段的可加成本。

對收費區(qū)域中的每一個OD對(m,n)∈ψ，至少有1條非循環(huán)路徑連接起點m和終點n，即內部路徑，收費區(qū)域中所有的內部路徑構成集合E。用虛擬路段替代OD對(m,n)間所有的內部路徑e∈E后，包含虛擬路段的新收費區(qū)域與外部網(wǎng)絡就組成了新的交通網(wǎng)絡G′=(N′,A′)。

圖3 網(wǎng)絡變換法示例Fig.3 Illustrative example of network transformation method

如圖3所示，原始收費區(qū)域由節(jié)點1，2，3和路段1，2，3構成，OD對1-3之間的兩條內部路徑分別1→2→3和1→3。經(jīng)過網(wǎng)絡變換后，內部路徑由1-3-1和1-3-2兩條虛擬路段代替。這就將原本基于距離的路徑不可加收費轉化為只存在于虛擬路段上的可加收費，每條路徑上的廣義出行成本就等于其上各路段的成本之和[10]。解決了路徑成本的不可加性問題后，基于虛擬路段的FW算法步驟如下：

4 算例分析

圖4 路網(wǎng)拓撲結構Fig.4 Road network topology

由表1中的數(shù)據(jù)可知，收費區(qū)域中的最短和最長距離分別為300 m和900 m。假設分段線性函數(shù)有3個間隔，每個間隔長度為200 m，得到4個邊界值l0=3，l1=5，l2=7，l3=9。參考新加坡的收費制度，將最小收費和最大收費分別定為y0=ymin=5，ymax=20。

表1 簡單網(wǎng)絡路段參數(shù)Tab.1 Link parameters of simple network

表2 內部路徑參數(shù)Tab.2 Parameters of internal routes

在Vissim中構建上述簡單網(wǎng)絡的仿真路網(wǎng)，將該仿真路網(wǎng)的仿真周期設置為3 600 s，每隔90 s采集1次數(shù)據(jù)，對出入路網(wǎng)的車輛數(shù)、路網(wǎng)內各路段流量、路網(wǎng)平均行程時間、路網(wǎng)平均延誤時間進行統(tǒng)計。由于擁擠收費策略的對象不包括行人、自行車和公交車，因此僅針對小汽車進行仿真。通過在路網(wǎng)進出路段布設數(shù)據(jù)采集點，得到仿真時間間隔路網(wǎng)內的累積車輛數(shù)，再對各路段流量進行距離加權，最后得到表示路網(wǎng)加權流量與累積車輛數(shù)的關系如圖5(a)所示。

圖5 模擬迭代過程Fig.5 Simulation and interative process

利用MATLAB對MFD散點圖進行擬合后，發(fā)現(xiàn)得到的擬合函數(shù)符合二次曲線形式：

qw=0.000 4N2+0.696 1N+56.582。

(13)

從MFD散點擬合圖和擬合函數(shù)中可以看出，當路網(wǎng)內累積車輛數(shù)處于870輛附近時，路網(wǎng)的加權流量達到其峰值359 veh/h，此時的流出量最大。而在本研究提出的收費模型中，上層目標函數(shù)正是要使路網(wǎng)的輸出流量最大，從而使系統(tǒng)的整體運行效益最優(yōu)。因此，這就確定了收費模型中上層模型的目標函數(shù)值，將離去率λ取為0.8，則路網(wǎng)最大輸出流量為287 veh/h。通過MATLAB求解文中提出的雙層規(guī)劃收費模型，在收費條件下得到目標函數(shù)值的迭代情況如圖5(b)所示。在算法迭代到33次時，輸出流量最接近于目標函數(shù)值，為286.79 veh/h。這一結果表明，此時在模型中應用的收費函數(shù)使路網(wǎng)整體運行效益最優(yōu)，即為最優(yōu)收費策略。由此得到的非線性最優(yōu)收費函數(shù)如圖6所示。

圖6 基于距離的最優(yōu)收費函數(shù)Fig.6 Distance-based optimal tolling function

從圖6中可以看出，該路網(wǎng)最低收費為5元，最高收費為16元。對在邊界內行駛距離較短的車輛，費用曲線增長平緩，收費力度較小，因此行駛中小距離的車輛可能還是傾向于選擇汽車出行。而對于邊界內的長距離行駛車輛，費用增長幅度較大，出行費用的明顯增加能在一定程度上鼓勵用戶縮短在邊界區(qū)域中的出行距離，同時也促進用戶的出行方式向公共交通轉變。

5 結論

本研究建立了一種將MFD與距離收費理論相結合的最優(yōu)定價模型，該模型中考慮固定需求下的用戶平衡，采用非線性距離收費函數(shù)影響用戶的出行選擇行為，從而使網(wǎng)絡的輸出流量最大。在求解過程中通過網(wǎng)絡變換的方式構建虛擬網(wǎng)絡，以解決非線性距離函數(shù)所引起的路徑成本不可加性問題，并基于虛擬網(wǎng)絡設計了基于路段的FW算法。通過算例驗證得到了使路網(wǎng)內車輛數(shù)維持在最優(yōu)水平附近的分段線性收費函數(shù)，證明了該收費模型和算法的可行性。值得指出的是，由于現(xiàn)實應用場景的多樣性，受制于基礎數(shù)據(jù)獲取的困難，各種影響收費的復雜因素往往難以完整考慮，需要根據(jù)不同的現(xiàn)場實際，選取適用因素并簡化模型應用。特別是實施收費后可能引發(fā)的各路段交通量變化和數(shù)據(jù)疊加效應，為制定收費方案增加了新的難度；同時，考慮到更好的公平性和接受度，后續(xù)的研究可在基于距離收費的基礎上加入時間或擁擠度等因素，結合MFD進一步探討多費率下的最優(yōu)定價及動態(tài)優(yōu)化問題，提出一種不僅精確并且具有廣泛實用性的改進模型及方法。