黃帥勛，李琦琦

（1.西南交通大學(xué) 交通運(yùn)輸與物流學(xué)院，四川 成都 611731；2.中國鐵路濟(jì)南局集團(tuán)有限公司車務(wù)段，山東 濟(jì)南 250000）

1 引言

自2013年6月鐵路貨運(yùn)市場(chǎng)化改革以來，政府逐步放松對(duì)鐵路定價(jià)權(quán)的管制，鐵路公司可以對(duì)部分鐵路在政府規(guī)定的運(yùn)價(jià)上限內(nèi)自由定價(jià)。鐵路運(yùn)價(jià)作為調(diào)節(jié)運(yùn)輸市場(chǎng)供給與需求的手段，直接關(guān)系到鐵路公司的盈利能力與運(yùn)輸生產(chǎn)安排。當(dāng)價(jià)格定得過高，運(yùn)輸單位貨物的利潤會(huì)提高，但是貨運(yùn)量會(huì)降低；當(dāng)價(jià)格定得過低，貨運(yùn)量會(huì)提升，但是運(yùn)輸單位貨物的利潤也會(huì)降低，所以必然存在使得鐵路收益最大的最優(yōu)價(jià)格。如圖1所示，當(dāng)鐵路公司制定出未來一段時(shí)間內(nèi)鐵路網(wǎng)絡(luò)中各站點(diǎn)之間的運(yùn)價(jià)后，運(yùn)輸市場(chǎng)在此運(yùn)價(jià)的吸引下會(huì)產(chǎn)生運(yùn)輸需求，鐵路公司根據(jù)該運(yùn)輸需求將會(huì)進(jìn)行車流徑路的規(guī)劃。所以鐵路公司科學(xué)地統(tǒng)籌制定鐵路運(yùn)價(jià)與規(guī)劃車流徑路，將是鐵路市場(chǎng)化改革的重要舉措。

圖1 貨運(yùn)定價(jià)與車流徑路的關(guān)系

在制定鐵路運(yùn)輸價(jià)格策略的研究中，李煜[1]將需求價(jià)格反應(yīng)函數(shù)引入到馬爾科夫動(dòng)態(tài)定價(jià)模型中，并運(yùn)用最大凹向包絡(luò)優(yōu)化方法進(jìn)行優(yōu)化。張小強(qiáng)等[2]充分考慮了多種運(yùn)輸方式競(jìng)爭(zhēng)下的鐵路貨運(yùn)定價(jià)，以客戶廣義費(fèi)用最小為目標(biāo)函數(shù)建立了雙層規(guī)劃模型。在車流徑路研究方面，田亞明等[3]以多商品流概念為基礎(chǔ)，建立了車流分配的點(diǎn)-弧模型和弧-路模型。將鐵路定價(jià)與運(yùn)輸生產(chǎn)結(jié)合起來的研究較少，劉丹[4]研究了競(jìng)爭(zhēng)環(huán)境下集裝箱的開行方案與動(dòng)態(tài)定價(jià)問題，但是研究重點(diǎn)是點(diǎn)對(duì)點(diǎn)單線的定價(jià)與開行方案綜合研究，網(wǎng)絡(luò)條件下的定價(jià)涉及較少。張旭[5]將定價(jià)問題與貨流分配問題結(jié)合起來研究，得出了鐵路與公路的貨流分配方案以及鐵路與公路的最優(yōu)運(yùn)價(jià)，但是未能將貨流轉(zhuǎn)化成列流。本文將在上述研究的基礎(chǔ)上，針對(duì)集裝箱運(yùn)輸?shù)亩▋r(jià)與車流徑路進(jìn)行綜合優(yōu)化，制定出鐵路運(yùn)輸網(wǎng)絡(luò)中各OD之間的集裝箱運(yùn)輸價(jià)格以及各OD之間的車流徑路，使得鐵路公司的總利潤最大。

2 模型建立

2.1 模型假設(shè)

假設(shè)1：在決策時(shí)間內(nèi)，運(yùn)輸市場(chǎng)能夠?qū)﹁F路運(yùn)價(jià)充分反應(yīng)；

假設(shè)2：在決策時(shí)間內(nèi)，運(yùn)輸市場(chǎng)沒有突發(fā)事件導(dǎo)致運(yùn)輸需求大幅度變化；

假設(shè)3：價(jià)格需求反應(yīng)函數(shù)在較小的價(jià)格波動(dòng)區(qū)間內(nèi)為線性函數(shù)。

2.2 符號(hào)說明

表1 符號(hào)說明

2.3 目標(biāo)函數(shù)建立

（1）運(yùn)輸重箱的收入

（2）運(yùn)輸重箱的成本

（3）箱流不平衡導(dǎo)致的空箱調(diào)運(yùn)費(fèi)用

式（3）中除以2的原因是消除正方向和反方向的箱流不平衡量重復(fù)計(jì)算。

綜合式（1）-式（3），目標(biāo)函數(shù)為：

2.4 約束條件建立

（1）決策期內(nèi)各站點(diǎn)的實(shí)際發(fā)送量小于各站點(diǎn)的集裝箱發(fā)送作業(yè)能力。

（2）決策期內(nèi)各站點(diǎn)的實(shí)際運(yùn)到量小于各站點(diǎn)的集裝箱到達(dá)作業(yè)能力。

（3）根據(jù)車流不可分原則，各站點(diǎn)之間的車流只能選擇唯一的路徑。

（4）決策期內(nèi)各個(gè)區(qū)段內(nèi)的列車開行數(shù)量小于區(qū)段通過能力。

（5）市場(chǎng)規(guī)律約束，反應(yīng)運(yùn)輸需求與運(yùn)輸定價(jià)之間的規(guī)律，本文采用線性的價(jià)格反應(yīng)函數(shù)。

（6）各站點(diǎn)間的實(shí)際運(yùn)輸量小于市場(chǎng)需求量。

（7）價(jià)格調(diào)整區(qū)間約束，一般最低價(jià)格不低于成本，最高價(jià)格不高于政府規(guī)定的上限。

（8）決策變量基本性質(zhì)約束。

3 模型求解

3.1 模型分析

通過分析式（4）-（12），可以看出該模型為非線性混合0-1規(guī)劃模型，決策變量有0-1變量和實(shí)數(shù)變量兩種類型。0-1變量的數(shù)量與鐵路運(yùn)輸網(wǎng)絡(luò)拓?fù)鋱D密切相關(guān)，但是可以確定至少有|個(gè)；實(shí)數(shù)變量有個(gè)，因此該模型中決策變量的數(shù)量與鐵路網(wǎng)絡(luò)中站點(diǎn)的數(shù)量呈二次冪函數(shù)增長。

目標(biāo)函數(shù)第一項(xiàng)為兩個(gè)實(shí)變量的乘積，第二項(xiàng)為實(shí)變量與0-1變量的乘積。約束中除了式（8）外，其余皆為線性約束，所以本模型目前沒有成熟的精確算法能夠求解。再者由于定價(jià)與車流徑路問題屬于戰(zhàn)術(shù)層面的決策問題，決策結(jié)果需要在未來較長時(shí)間內(nèi)執(zhí)行，對(duì)于運(yùn)輸市場(chǎng)的預(yù)測(cè)不會(huì)非常精確，因此對(duì)于決策精度的要求并不是太高，所以本文決定采用遺傳算法進(jìn)行求解。

3.2 遺傳算法求解過程

遺傳算法屬于啟發(fā)式算法，是模仿生物界種群的進(jìn)化過程，其核心內(nèi)容包括基因型與表現(xiàn)型的對(duì)應(yīng)關(guān)系、自然選擇、個(gè)體交配、基因變異等，其流程如圖2所示，具體求解步驟如下所示：

圖2 遺傳算法流程

（1）編碼。編碼是模擬遺傳過程中基因型與表現(xiàn)型的關(guān)系。分析本文的模型，其解空間由puv和λluv組成，其中puv為實(shí)數(shù)變量，其數(shù)量有個(gè)，為0-1變量，其數(shù)量至少有個(gè)，并且與鐵路網(wǎng)絡(luò)拓?fù)鋱D密切相關(guān)。本文采用二進(jìn)制編碼，但是puv的含義是價(jià)格，所以其表現(xiàn)型為十進(jìn)制，必須將其轉(zhuǎn)化為二進(jìn)制的基因型。其中二進(jìn)制編碼位數(shù)與puv的上下限和精度要求密切相關(guān)，從模型中可以看出puv的上下限即為和。λluv的含義為OD是否選擇某條路徑，其表現(xiàn)型本身為二進(jìn)制，所以其二進(jìn)制編碼位數(shù)為1位，具體的編碼次序排列如圖3所示。

圖3 決策變量編碼

（2）產(chǎn)生初始種群。初始種群通過計(jì)算機(jī)隨機(jī)產(chǎn)生。每個(gè)基因的每個(gè)位置從0和1中隨機(jī)選取一個(gè)數(shù)字組成一個(gè)染色體，稱為種群的一個(gè)個(gè)體。通過觀察發(fā)現(xiàn)約束（3），即式（7）表明了變量λluv之間存在約束關(guān)系（車流不分離原則，一個(gè)車流只能選擇一個(gè)徑路），也就意味著通過完全隨機(jī)產(chǎn)生的初始λluv不在解空間中，這勢(shì)必會(huì)拖慢收斂速度。為了避免這種情況出現(xiàn)，本文將改進(jìn)原始的遺傳算法，通過帶有約束的隨機(jī)數(shù)產(chǎn)生方法，保證產(chǎn)生的初始種群全部滿足約束（3），這樣將大大提高遺傳算法的收斂速度。

（3）適應(yīng)度函數(shù)。適應(yīng)度函數(shù)是模仿進(jìn)化過程中的自然選擇，當(dāng)個(gè)體的適應(yīng)度函數(shù)越大，就越能被自然選擇中繁殖下一代。本文模型的適應(yīng)度函數(shù)為：

（4）新種群的產(chǎn)生。根據(jù)產(chǎn)生的初始種群或者上一代種群的個(gè)體適應(yīng)度函數(shù)來決定哪些個(gè)體能夠生產(chǎn)后代。每個(gè)個(gè)體的選取過程作為一個(gè)獨(dú)立時(shí)間，并且個(gè)體被選擇的概率與其適應(yīng)度函數(shù)的大小呈正比。在實(shí)際的操作過程中，本文用“輪盤賭”的方式進(jìn)行模擬選擇，對(duì)所有個(gè)體的適應(yīng)度函數(shù)值求和；然后將每個(gè)個(gè)體的適應(yīng)度函數(shù)值除以所有個(gè)體適應(yīng)度函數(shù)之和的商作為每個(gè)個(gè)體所對(duì)應(yīng)扇形面積的比例；最后，轉(zhuǎn)動(dòng)輪盤確定出具有交配權(quán)的個(gè)體。通過這種方式可以保證較優(yōu)的解被選擇的概率大于次優(yōu)解的概率，從而使得新解可以不斷向最優(yōu)值收斂。

（5）交叉。交叉是模擬種群進(jìn)化過程中的交配繁殖過程。生物的遺傳繁殖是通過父母的染色體配對(duì)產(chǎn)生新的染色體的過程，在這個(gè)過程中子代染色體會(huì)將父親和母親的基因各繼承一半。本文模型由于變量具有明顯的約束關(guān)系（式（7）），所以在進(jìn)行交叉時(shí)，為了避免產(chǎn)生不滿足式（7）的新個(gè)體，本文將改進(jìn)傳統(tǒng)遺傳算法的交叉操作，只對(duì)控制puv的基因位進(jìn)行交叉操作，而對(duì)控制λluv的基因進(jìn)行特殊的交叉操作，保證產(chǎn)生的新個(gè)體的控制λluv的基因滿足式（7）。

（6）變異。變異是模擬生物進(jìn)化過程中的基因突變。在模擬的過程中，首先產(chǎn)生一個(gè)0-1的隨機(jī)數(shù)，如果該值小于設(shè)定的變異概率則個(gè)體發(fā)生變異，將該個(gè)體基因從隨機(jī)位置開始0變異為1，1變異為0.與交叉過程的特殊操作類似，由于變量具有明顯的約束關(guān)系（式（7）），變異僅對(duì)控制puv的基因位進(jìn)行變異操作，對(duì)控制λluv的基因重新隨機(jī)生成滿足式（7）的變異操作。

（7）終止條件。遺傳算法的終止條件為當(dāng)進(jìn)化的最優(yōu)解收斂到特定的值，并且很難再有更優(yōu)的解出現(xiàn)。一般通過限制種群的進(jìn)化代數(shù)來實(shí)現(xiàn)，當(dāng)種群的進(jìn)化代數(shù)設(shè)置過大時(shí)，雖然會(huì)有更大的概率搜索出更優(yōu)的解，但是計(jì)算量會(huì)相應(yīng)地增加；當(dāng)進(jìn)化代數(shù)設(shè)置過小時(shí)，雖然計(jì)算量較少，但是有可能會(huì)出現(xiàn)結(jié)果不收斂的情況。

4 案例分析

4.1 案例介紹

本文選取廣州鐵路（集團(tuán)）公司管內(nèi)的6個(gè)集裝箱站作為案例，如圖4所示，鐵路運(yùn)輸網(wǎng)絡(luò)中的6個(gè)集裝箱站用方框表示，從1到6依次為：郴州站、湘潭東站、衡陽站、長沙東站、岳陽北站和益陽東站。需要注意的，本文研究的是技術(shù)站之間的車流徑路問題，本案例為了突出研究重點(diǎn)，未將各技術(shù)站之間的中間站標(biāo)出，中間站之間的車流較簡(jiǎn)單，屬于管內(nèi)區(qū)段列車的研究范疇。各站點(diǎn)之間的連線表示鐵路區(qū)段，圖中的區(qū)段都為雙向的區(qū)段，區(qū)段編號(hào)從1-12，其中上行區(qū)段用白圓圈表示，有1，2，3，4，5，11；下行區(qū)段用黑圓圈表示，有6，7，8，9，10，12。表2為各站點(diǎn)的相關(guān)參數(shù)，表3為鐵路區(qū)段的相關(guān)參數(shù)。

圖4 鐵路運(yùn)輸網(wǎng)絡(luò)拓?fù)鋱D

表2 各站點(diǎn)集裝箱發(fā)送與到達(dá)作業(yè)能力（箱/天）

表3 鐵路區(qū)段相關(guān)參數(shù)

根據(jù)廣州鐵路公司市場(chǎng)調(diào)查的結(jié)果，各OD之間的價(jià)格需求反應(yīng)函數(shù)參數(shù)見表4和表5，優(yōu)化前各OD之間的集裝運(yùn)輸價(jià)格見表6，政府規(guī)定的運(yùn)價(jià)上限見表7。廣州鐵路公司需要確定各OD之間的集裝箱運(yùn)輸價(jià)格以及貨流徑路，使得鐵路公司的總收益最大化。

表4 各OD之間集裝箱運(yùn)輸價(jià)格需求反應(yīng)函數(shù)參數(shù)a

表5 各OD之間集裝箱運(yùn)輸價(jià)格需求反應(yīng)函數(shù)參數(shù)b

表6 優(yōu)化前各OD之間的集裝箱運(yùn)輸價(jià)格（元/箱）

表7 各OD之間政府規(guī)定的運(yùn)價(jià)上限（元/箱）

4.2 求解結(jié)果

將4.1節(jié)中的數(shù)據(jù)代入到第2節(jié)的模型中，根據(jù)第3節(jié)的改進(jìn)遺傳算法，利用MATLAB編程進(jìn)行求解。將種群規(guī)模設(shè)置為200個(gè)，進(jìn)化代數(shù)設(shè)置為1 000代，交叉概率設(shè)置為0.9，變異概率設(shè)置為0.01，用處理器為Intel(R)Core(TM)i3-2130M CPU@2.10GHz，內(nèi)存為4G的計(jì)算機(jī)進(jìn)行求解。目標(biāo)函數(shù)值隨著進(jìn)化代數(shù)的變化曲線如圖5所示，經(jīng)過400代開始收斂，目標(biāo)函數(shù)最優(yōu)值為38 489 692.85元。各OD之間的最優(yōu)集裝箱運(yùn)輸價(jià)格見表8，各OD之間的最優(yōu)徑路選擇見表9，各OD之間的最優(yōu)車流量見表10，各區(qū)段車流量與飽和度見表11。

圖5 改進(jìn)的遺傳算法求解過程

表8 各OD之間的最優(yōu)集裝箱運(yùn)輸價(jià)格（元/箱）

表9 各OD之間的最優(yōu)車流路徑

表10 各OD之間的最優(yōu)車流量（箱）

表11 各區(qū)段車流量以及飽和度

4.3 結(jié)果分析

（1）優(yōu)化前后結(jié)果對(duì)比。將4.2節(jié)的求解結(jié)果代入到式（1）、式（2）和式（3）中，得到最優(yōu)情況下的重箱運(yùn)輸收入為130 200 576.19元，重箱運(yùn)輸成本為90 683 705.27元，箱流不平衡導(dǎo)致的空箱調(diào)運(yùn)費(fèi)用為：1 027 178.07元，總利潤為：38 489 692.85元。將其與優(yōu)化前的利潤進(jìn)行對(duì)比，見表12。

表12 優(yōu)化前后結(jié)果對(duì)比

從表12可以看出雖然優(yōu)化前的總收入較高，但是優(yōu)化前的運(yùn)輸成本與空箱調(diào)運(yùn)費(fèi)用也較高，優(yōu)化后比優(yōu)化前的利潤提高了8 420 419.88元，增長28.00%。從優(yōu)化前與優(yōu)化后的價(jià)格看，除了湘潭東到衡陽站、衡陽站到郴州站、長沙東站到益陽東站之外的OD的運(yùn)價(jià)都有不同程度的上升。箱流不平衡量也減少了844箱，相應(yīng)的空箱調(diào)運(yùn)費(fèi)用也減少了84 379元。通過以上結(jié)果分析，可以看出本文模型的優(yōu)越性。

（2）與傳統(tǒng)的遺傳算法求解結(jié)果對(duì)比。由于本文改進(jìn)了傳統(tǒng)的遺傳算法，在生成初始解以及交叉、變異的過程中，采用特殊的策略，始終使產(chǎn)生的新的種群滿足式（7）。如圖6所示，將本文改進(jìn)的遺傳算法與傳統(tǒng)的遺傳算法比較，可以看出本文改進(jìn)的遺傳算法具有收斂快、結(jié)果更優(yōu)的特點(diǎn)。

5 結(jié)語

本文首先提出了鐵路集裝箱運(yùn)輸定價(jià)與車輛路徑問題的背景，接著以鐵路公司的總利潤為目標(biāo)函數(shù)建立了非線性混合整數(shù)規(guī)劃模型，然后設(shè)計(jì)了改進(jìn)遺傳算法對(duì)其進(jìn)行求解。最后以廣州鐵路公司的6個(gè)集裝箱站點(diǎn)作為案例，通過計(jì)算得出本文的模型和算法將使得鐵路公司的利潤提高了28.00%，并且比較了本文改進(jìn)的遺傳算法與傳統(tǒng)的遺傳算法，得出本文改進(jìn)的遺傳算法具有收斂速度快的優(yōu)點(diǎn)。