基于彈性需求的高速鐵路列車開行方案優(yōu)化*

楊宇正，周文梁

(1.廣州鐵路(集團)公司，廣東廣州 510088;2中南大學(xué) 商學(xué)院，湖南長沙 410083)

列車開行方案確定了高速鐵路上各類開行列車的數(shù)量、運行區(qū)段以及沿途停站序列等，是高速鐵路列車運營組織的綜合計劃。列車開行方案的編制質(zhì)量直接決定了高速鐵路列車運營質(zhì)量與旅客出行服務(wù)水平，高質(zhì)量的列車開行方案不但能夠有效提高高速鐵路運營效益與效率，而且還能因提供旅客高水平的出行服務(wù)而吸引更多客流選擇高速鐵路出行，提高高速鐵路客運市場份額。由于列車開行方案質(zhì)量影響到旅客出行服務(wù)水平，而旅客出行服務(wù)水平在一定程度上決定了旅客出行需求量，當旅客列車開行方案改變而導(dǎo)致旅客出行服務(wù)水平變化時，旅客出行需求量將會作相應(yīng)的變化。目前，關(guān)于固定需求下旅客列車開行方案優(yōu)化問題的研究成果［1－6］較多，而考慮旅客列車開行方案質(zhì)量差別對客運需求的影響來優(yōu)化列車開行方案的研究相對較少，文獻［7－8］設(shè)計了基于彈性需求的鐵路旅客列車開行方案優(yōu)化方法，但這僅考慮了鐵路全天OD需求量與旅客出行服務(wù)水平的變化關(guān)系，但對于高速鐵路而言，由于全天不同時段客流需求的差異性，列車開行方案分時段進行編制，各時段列車開行方案質(zhì)量不僅決定了本時段客流出行需求量，而且還影響到其他時段客流需求量大小。換而言之，各時段客流出行需求量同時受到自身時段與其它時段內(nèi)列車開行方案的影響。本文考慮不同時段列車開行方案與客流出行需求量之間的相互影響關(guān)系，通過建立高速鐵路各時段客流出行需求量關(guān)于各時段旅客出行服務(wù)水平間的函數(shù)關(guān)系，將旅客是否選擇高速鐵路出行以及選擇出行列車行為描述為彈性需求下的用戶均衡分配問題。在此基礎(chǔ)上，構(gòu)建彈性需求下高速鐵路列車開行方案優(yōu)化的雙層規(guī)劃模型，并設(shè)計相應(yīng)的模擬退火求解算法。

1 高速鐵路旅客出行彈性需求分析

1.1 旅客出行需求影響因素與出行需求結(jié)構(gòu)分析

高速鐵路旅客出行需求受多方面因素影響，主要包括旅客出行心理因素、高速鐵路運輸產(chǎn)品性能因素以及外部環(huán)境因素。這里僅考慮列車開行質(zhì)量對旅客出行需求的影響，故只需分析與旅客列車開行質(zhì)量相關(guān)的旅客出行各項費用(如換乘費用、票價)因素，而假設(shè)其他影響因素不變。

根據(jù)旅客出行需求受出行費用的影響程度將其分為剛性出行需求、轉(zhuǎn)移出行需求及誘增出行需求3部分。其中，剛性出行需求指高速鐵路上不受旅客出行費用影響的出行需求量，為高速鐵路上出行需求的最小值，如圖1(a)所示;轉(zhuǎn)移出行需求指因旅客出行費用的改變在高速鐵路與其他客運方式之間轉(zhuǎn)移的出行需求，當旅客出行費用減少，則存在部分需求從其他客運方式轉(zhuǎn)移到高速鐵路上，反之，則存在部分需求從高速鐵路上轉(zhuǎn)移到其他客運方式，如圖1(b)所示;誘增出行需求指因旅客出行費用減少而誘增產(chǎn)生的部分出行需求，旅客出行費用越小，誘增的出行需求越多，如圖1(c)所示。

1.2 旅客出行彈性需求函數(shù)建立

旅客出行彈性需求函數(shù)反映了高速鐵路旅客出行需求量與旅客出行費用之間的變化關(guān)系。在彈性需求的概念下，可以視確定需求為彈性需求的均值，或者說，彈性需求可以由確定需求這一均值波動而來。這種波動關(guān)系與旅客出行費用相關(guān)，若旅客出行費用相對較低，則存在更多需求從其他客運方式中轉(zhuǎn)移到高速鐵路上，同時誘增更多的出行客流，從而導(dǎo)致高速鐵路出行需求量增加;反之，高速鐵路出行需求減少。以高速鐵路出行需求預(yù)測量作為這一確定需求量，結(jié)合預(yù)測出行需求時的出行費用水平便可建立高速鐵路旅客出行需求與出行費用間的彈性需求函數(shù)關(guān)系式如下:

圖1 旅客出行需求量與出行費用變化關(guān)系圖Fig.1 Varying relation between passenger travel demand and travel cost

其中:fijt為時段t起點i至終點j的旅客出行需求量;f0ijt和分別為起點i至終點j的剛性出行需求量、出行需求預(yù)測量;μijt為客流fijt的最小出行費用;為起點i至終點j的時段t之后時段的平均出行費用;為預(yù)測出行需求量ˉfijt時的出行費用水平;R1(μijt－)和R2(μijt－)分別為轉(zhuǎn)移出行需求、誘增出行需求波動率與μijt和兩出行費用之差以及后續(xù)時段平均出行費用的變化關(guān)系。

2 高速鐵路列車開行方案優(yōu)化的雙層規(guī)劃模型

高速鐵路列車開行方案的制定涉及到高速鐵路企業(yè)與旅客雙方面利益，對于高速鐵路企業(yè)，在一定的動車組運用數(shù)量、線路資源能力限制條件下，以降低列車開行成本、提高列車開行收益為目的制定列車開行方案;旅客只能在高速鐵路企業(yè)制定的列車開行方案的限制下，根據(jù)其出行廣義費用決定是否選擇高速鐵路出行以及選擇出行列車，而旅客的這種出行選擇結(jié)果將直接影響到高速鐵路列車開行效益，成為高速鐵路企業(yè)進一步完善列車開行方案依據(jù)。由此，可將高速鐵路列車開行方案制定問題描述為一個Leader－Follower問題，其中高速鐵路企業(yè)屬于指導(dǎo)者(Leader)，而旅客屬于跟隨者(Follower)。高速鐵路通過列車開行方案來改變旅客廣義出行費用，從而影響旅客出行選擇行為，但不能控制旅客出行選擇行為，以下將構(gòu)建一個雙層規(guī)劃來描述這種關(guān)系。

根據(jù)高速鐵路客流需求將高速鐵路全天運營時間劃分為t=1，2，…，h個時段，設(shè)高速鐵路網(wǎng)N=(S，E)上時段t內(nèi)的列車開行方案為Ωt={T=(Q，q)}，其中，對于旅客列車T∈Ωt，用QT和qT分別表示列車T的停站序列與開行數(shù)量(列車編組一般為固定的一組動車組)，S={1，2，…，n}為高速鐵路車站集，E為高速鐵路路段集，de表示路段e∈E的里程。設(shè)時段t起點i到終點j的客流fijt可選擇 σtij種換乘方案Hkijt，k=1，2，…，σtij，則相應(yīng)的客流量為fkijt。由此便可將上層規(guī)劃，即高速鐵路企業(yè)制定列車開行方案的決策過程描述如下:

式中:Pkijt為旅客選擇出行路徑Hkijt出行所需花費的票價;cL為列車行駛的單位公里費用，cH為動車組使用的單位小時費用;v為高速鐵路列車的旅行速度;cM為列車的開行組織費用。當車站i為旅客列車 T 的始發(fā)站時，η(i，T)=1，否則，η(i，T)=0;L(i)為車站i的最大發(fā)車能力;M為可使用的動車組總數(shù)。

在以上上層規(guī)劃中，式(2)表示為高速鐵路開行列車效益(開行列車票價收入減去各項開行成本)最大化，式(3)表示表示車站發(fā)車能力限制，式(4)表示動車組數(shù)量限制。

上層規(guī)劃中各條出行路徑的客流量fkijt是在給定的列車開行方案Ωt條件下旅客決定是否出行及選擇換乘方案所形成的路徑流量，這是一個基于彈性需求的用戶平衡分配問題。它構(gòu)成求解f={fkijt}的下層規(guī)劃如下:

其中:Ukijt為旅客選擇出行路徑Hkijt出行所需花費的出行時間價值;yT(x)為旅客列車T乘載客流量x的擁擠費用函數(shù)。利用旅客列車T的定員xT定義擁擠費用函數(shù)如下:

3 模擬退火算法設(shè)計

將采用模擬退火算法求解以上所構(gòu)建的列車開行方案優(yōu)化的雙層規(guī)劃模型。首先，以“按流開車”為原則生成初始列車開行方案，進而通過構(gòu)造模擬退火算法鄰域系不斷產(chǎn)生新解進行反復(fù)迭代。模擬退火算法的設(shè)計必須解決好以下3個問題。

(1)初始列車開行方案生成。根據(jù)高速鐵路線路區(qū)間里程尋找各OD對之間的里程最短路，為提高列車旅行速度，按照路徑里程由大到小順序逐個判斷，只要最短路上客流量達到開行列車要求便可開行，并盡可能地吸引經(jīng)過此路徑的其他OD客流量，直至將所有客流量安排到列車上為止獲得初始列車開行方案。初始列車開行方案中列車停站采用站站停方式，在列車開行方案調(diào)整階段將逐漸調(diào)整停站使之滿足要求。

(2)解的鄰域系的構(gòu)造。構(gòu)造開行方案Ω的鄰域系包括剔除和添加列車兩個步驟。剔除列車指從旅客列車開行方案Ω中剔除部分列車;添加列車指向旅客列車開行方案Ω中添加新的列車，而新列車采用“按流開車”的原則［4］產(chǎn)生。

(3)模擬退火計劃表確定。模擬退火計劃表構(gòu)造為:初始溫度Ψ0由數(shù)值方法確定;溫度下降規(guī)則采用比例下降法，比例值θ隨著溫度下降而逐步提高;同一溫度的迭代次數(shù)采用迭代次數(shù)下限、上限和接受頻率來綜合控制;算法終止規(guī)則采用溫度下降至接近0或目標函數(shù)在λ次迭代無改變來綜合控制。

在生成初始列車開行方案、構(gòu)造解的鄰域系、安排模擬退火計劃表的基礎(chǔ)上，便可設(shè)計求解高速鐵路列車開行方案雙層規(guī)劃模型的模擬退火算法如下。

第0步:初始化。求取初始溫度值Ψ0，設(shè)置溫度下降比例θ，終止溫度Ψend，當前溫度Ψ =Ψ0，迭代次數(shù)Κ =1，轉(zhuǎn)第1步。

第1步:生成初始解Ω0。按以上方法生成初始列車開行方案Ω0，并置Ω =Ω0，轉(zhuǎn)第2步。

第2步:用戶平衡配流。基于列車開行方案Ω求解下層規(guī)劃得到符合客流均衡條件的彈性需求量f，計算目標函數(shù)值G，轉(zhuǎn)第3步。

第3步:生成鄰域解。計算列車開行方案Ω中每列列車的開行效益GT，根據(jù)GT值以概率形式從列車開行方案Ω中剔除一定比例的列車，相應(yīng)地獲得這部分旅客列車上的客流，將這部分客流按“按流開車”的原則產(chǎn)生部分新列車，并將它們添加到列車開行方案Ω中去得到鄰域解Ω'。轉(zhuǎn)第4步。

第4步:用戶平衡配流。基于列車開行方案Ω'求解下層規(guī)劃得到符合客流均衡條件的彈性需求量f'，計算目標函數(shù)值G'，轉(zhuǎn)第5步。

第5步:比較最優(yōu)解與鄰域解。若G'＞G，或G'≤ G 且 Random(0，1)＜exp［(G － G')/Ψ ］，則接受新解，令Ω =Ω'，G=G';否則，保持原有最優(yōu)解不變。轉(zhuǎn)第6步。

第6步:迭代次數(shù)檢驗。令Κ =Κ+1，若Κ大于溫度Ψ下的迭代次數(shù)下限LΨ且較劣解的接受頻率大于給定的接受概率，或若Κ小于溫度Ψ時的迭代次數(shù)上限UΨ，則令Ψ =Ψ×θ后，轉(zhuǎn)第7步，否則返回第3步。

第7步:算法終止檢驗。若Ψ ＜Ψend，或最優(yōu)解Ω已經(jīng)在λ次迭代中沒有發(fā)生改變，則算法終止，輸出最優(yōu)解Ω;否則，重置迭代次數(shù)Κ =1，并返回第3步。

4 算例分析

為了驗證模型與算法的有效性，構(gòu)造如圖2所示的由13個客運站和20個復(fù)線運行區(qū)間組成的高速鐵路網(wǎng)。該網(wǎng)絡(luò)上各車站之間存在高時間價值和低時間價值兩類出行旅客，其中，高時間價值旅客的時間價值β為40元/h，低時間價值旅客的時間價值為20元/h。高時間價值旅客所采用的費用權(quán)重系數(shù)α1和α2分別為0.5和0.5，而低時間價值旅客所采用的費用權(quán)重系數(shù)分別為0.7和0.3。考慮開行高、中速兩類列車，其中，高速列車的開行組織費用為4 200元，列車公里費用為59.2元/km，車輛公里費用為5.5元/km，車輛的平均載客量為70人，列車上的人公里票價率為0.35元/(人·km)，列車平均運行速度為325 km/h;而中速列車的開行組織費用為3 700元，列車公里費用為45.6元/km，車輛公里費用為3.8元/km，車輛的平均載客量為70人，列車上的人公里票價率為0.30元/(人·km)，列車平均運行速度為250 km/h。兩類列車上旅客擁擠費用計算參數(shù)a，b分別為10.0和2.5，列車開行所需要的最低載客量為800人/列。車站內(nèi)任意2列車之間的平均換乘時間取為0.5 h。

圖2 高速鐵路網(wǎng)絡(luò)結(jié)構(gòu)示意圖Fig.2 The structure diagram of high －speed railway network

為便于對比，首先將旅客運輸需求函數(shù)中需求波動強度系數(shù)設(shè)為r=0/元，此時各OD對之間旅客出行需求量，故可將此情形下確定的旅客列車開行方案看作為固定需求下優(yōu)化得到的旅客列車開行方案，其中，高速旅客列車204對，中速旅客列車159對，鐵路運輸企業(yè)總收益為19 065 441.5元。在其他數(shù)據(jù)保持不變的條件下，將旅客運輸需求函數(shù)中最大需求比例系數(shù)設(shè)為η=1.2，需求波動強度系數(shù)設(shè)為r=0.05/元，優(yōu)化計算得到彈性需求條件下旅客列車開行方案，其中，高速旅客列車191對，中速旅客列車158對，鐵路運輸企業(yè)總收益為19 967 276.3元，各項列車開行費用如表1中彈性需求部分所示。

表1 旅客列車組織收益表Table 1 The benefit and cost of train organization

由表1可知:雖然固定需求條件下旅客列車開行方案的票價收入比彈性需求條件下旅客列車開行方案的票價收入高，但由于其開行的旅客列車較多，相應(yīng)的列車開行費用增多，從而導(dǎo)致其總效益反而比彈性需求條件下總效益少。這主要是因為固定需求條件時沒有考慮旅客出行需求量受其廣義出行費用的影響，廣義出行費用相對低的OD對出行需求量沒有增加;而彈性需求條件下則充分考慮了旅客需求量與出行費用的關(guān)系，產(chǎn)生了更好的社會與經(jīng)濟效益。

5 結(jié)論

高速鐵路各時段列車開行方案質(zhì)量不僅決定了本時段客流出行需求量，而且還影響到其它時段客流需求量。本文考慮不同時段列車開行方案與客流出行需求量之間的相互影響關(guān)系，在建立高速鐵路各時段客流出行需求量關(guān)于各時段旅客出行服務(wù)水平間的函數(shù)關(guān)系的基礎(chǔ)上，結(jié)合彈性需求下旅客出行乘車選擇行為，構(gòu)建彈性需求下高速鐵路列車開行方案優(yōu)化的雙層規(guī)劃模型，并設(shè)計相應(yīng)的模擬退火求解算法。由此確定的高速鐵路列車開行方案充分體現(xiàn)了高速鐵路企業(yè)與旅客出行選擇之間的內(nèi)在制約關(guān)系。

［1］史 峰，周文梁，陳 彥.基于彈性需求的旅客列車開行方案優(yōu)化研究［J］.鐵道學(xué)報，2008，30(3):1 －6.SHI Feng，ZHOU Wen-liang，CHEN Yan.Optimization study on passenger train plans with elastic demands［J］.Journal of the China Railway Society，2008，30(3):1 －6.

［2］史 峰，鄧連波.旅客換乘網(wǎng)絡(luò)優(yōu)化設(shè)計［J］.鐵道科學(xué)與工程學(xué)報，2004，1(1):78 －84.SHI Feng，DENGLian-bo.Optimal design of passenger transfer network［J］.Journal of Railway Science and Engineering，2004，1(1):78 －82.

［3］周文梁，史 峰，陳 彥.客運專線網(wǎng)絡(luò)列車開行方案與運行圖綜合優(yōu)化方法［J］.鐵道學(xué)報，2011，32(1):1 －7.ZHOUWen-liang，SHIFeng，CHENYan.Integrated optimization method of train plan and diagram on dedicated passenger network［J］.Journal of the China Railway Society，2011，32(1):1－7.

［4］史 峰，鄧連波，霍 亮.旅客列車開行方案的雙層規(guī)劃模型和算法［J］.中國鐵道科學(xué)，2007，28(3):110 －116.SHI Feng，DENG Lian-bo，HUO Liang.Bi－ level programming model and algorithm of passenger train operation plan［J］.China Railway Science，2007，28(3):110 －116.

［5］何宇強，毛保華.高速鐵路旅客列車開行方案的多目標雙層規(guī)劃模型［J］.鐵道學(xué)報，2006，28(5):6 －10.HE Yu-qiang，MAO Bao-hua.Multiobjective bi－ level programming model of making train working plan for passenger－ only line［J］.Journal of the China Railway Society，2006，28(5):6－10.

［6］查偉雄，符 卓.直通旅客列車開行方案優(yōu)化方法研究［J］.鐵道學(xué)報，2000，22(5):1 －5.ZHA Wei-xiong，F(xiàn)U Zhuo.Research on the optimization method of through passenger train plan ［J］.Journal of the China Railway Society，2000，22(5):1 －5.

［7］史 峰，周文梁，陳 彥.基于彈性需求的旅客列車開行方案優(yōu)化研究［J］.鐵道學(xué)報，2008，30(3):1 －6.SHI Feng，ZHOU Wen-liang，CHEN Yan.Optimization study on passenger train plans with elastic demands［J］.Journal of the China Railway Society，2008，30(3):1 －6.

［8］周文梁.客運專線網(wǎng)絡(luò)列車開行方案與運行圖綜合優(yōu)化模型及算法［D］.長沙:中南大學(xué)，2010.ZHOU Wen－liang.Intergrated optimization model and algorithm of train plan and train schedule on dedicated passenger network［D］.Changsha:Central South University，2010.