王海兵，杜延帥，王文憲

WANG Hai-bing, DU Yan-shuai, WANG Wen-xian

（西南交通大學交通運輸與物流學院，四川成都610031）

（School of Transportation and Logistics， Southwest Jiaotong University， Chengdu 610031， Sichuan， China）

基于集對分析的高速列車開行方案比選

王海兵，杜延帥，王文憲

WANG Hai-bing, DU Yan-shuai, WANG Wen-xian

（西南交通大學交通運輸與物流學院，四川成都610031）

（School of Transportation and Logistics， Southwest Jiaotong University， Chengdu 610031， Sichuan， China）

在闡述集對分析多方案評價方法的基礎上，針對高速列車開行方案比選情況構(gòu)建因素層、指標層和方案層3層評價指標體系，根據(jù)集對分析原理選取20個評價指標，并且給出各指標具體定義和指標值確定方法。最后以某高速鐵路為例，通過對3組備選方案20個指標值進行綜合評價，構(gòu)造決策矩陣，再根據(jù)同一度原則計算得到同一度決策矩陣，采用層次分析法確定決策矩陣中各指標所占權(quán)重進一步計算得到帶權(quán)同一度矩陣，對帶權(quán)同一度矩陣中各備選方案相應值進行比較，通過比選得到高速鐵路最優(yōu)實施方案。

開行方案；綜合評價指標；集對分析；帶權(quán)同一度矩陣

鐵路旅客運輸組織質(zhì)量與旅客列車開行方案質(zhì)量密切相關(guān)。針對各旅客列車開行方案建立一套完整的評價體系，結(jié)合經(jīng)濟和服務等影響因素對整個系統(tǒng)體系進行評估，不僅可以用于評價旅客列車開行方案質(zhì)量，還可以作為基礎來優(yōu)化實際選定的開行方案[1]。旅客列車開行方案通常會受列車徑路、列車種類、?？空拘蛄屑傲熊嚾臻_行對數(shù)等多種路線方案因素的影響。雖然目前國內(nèi)外對旅客列車開行方案比選進行了大量研究，但都側(cè)重于各種方案的技術(shù)經(jīng)濟比較[2-5]，未考慮旅客服務質(zhì)量和舒適性等特性，從而使使用范圍受限。在研究基礎上，從鐵路運營部門的效益、能力、技術(shù)特性及旅客服務特性出發(fā)，構(gòu)建旅客列車開行方案評價體系模型，并且在集對分析原理的基礎上，綜合考慮實際工程中的確定和不確定因素，最后結(jié)合某高速鐵路實例，對 3 種旅客列車開行方案進行方案比選，得到最優(yōu)實施方案。

1 集對分析多方案評價方法

集對分析 (Set Pair Analysis，SPA)[6]理論核心是辯證分析系統(tǒng)的確定性和不確定性，通過假設所有事物的本質(zhì)屬性為不確定性，綜合考察一個系統(tǒng)的確定性與不確定性，對其各子系統(tǒng)采用數(shù)學理論及模型進行處理，再結(jié)合構(gòu)建的同一度決策矩陣對所選方案進行全面評價，選出旅客列車開行方案中最符合實際的方案組合。

1.1 構(gòu)造同一度決策矩陣

假設某高速鐵路全線中有 A1，A2，…，Am共m 個待選開行方案，而每個方案又包含 f1，f2，…，fn共 n 個評價指標，進一步可以得到多方案評價的決策矩陣 H。

從 m 個方案的各指標中選取最優(yōu)值 (效益型指標選取最大值，成本型指標選取最小值) 組成新的理想方案 A0，并且用 f0j表示新理想方案中第 j ( j = 1，2，…，n) 個指標的值；計算所選理想方案中各指標值 fij；計算所選理想開行方案 A0中各對應的指標值 f0j，求得二者的同一度 aij，由此得到基于集對分析的理想方案指標和被評價方案指標不帶權(quán)的同一度矩陣 Q 如下。

1.2 評價指標權(quán)重

使用層次分析法確定評價指標權(quán)重的計算步驟如下[7]。

（1）構(gòu)建判斷矩陣。結(jié)合已經(jīng)建立的遞階層次結(jié)構(gòu)模型，對模型中各因素采用 1—9 及其倒數(shù)標度的方法進行相互間對比分析，得到判斷矩陣 B = {bij}m×n。

（2）求解矩陣 B 的特征向量 W 和最大特征值λmax。特征向量是指某一層與上一層各因素間的重要程度經(jīng)排序后得到的權(quán)重值其中 bi=并且依次得到特征向量 W = (w1，w2，…，wn)，最大特征根 λmax=

（3）計算一致性指標。對判斷矩陣進行驗證，計算一致性指標 CI = (λmax-n) / (n-1)，得到其平均隨機一致性指標值如表1 所示。

表1 平均隨機一致性指標

如果隨機一致性比率滿足條件CR = CI/RI＜0.10，即可認為所得的排序結(jié)果滿足上述條件，從而得到合理分配的權(quán)重系數(shù)；如果不能滿足條件，則需要重新調(diào)整矩陣 B 中數(shù)值。

1.3 綜合評價模型

對被評價方案 Ai和理想方案 A0的帶權(quán)同一度矩陣 R 進行求解，計算公式為

i = 1，2，…，m

ai是第 i 個被評價方案與理想方案的帶權(quán)同一度之和，計算公式為

求解各被評價方案的權(quán)重，次序可以通過帶權(quán)同一度矩陣 R 中 ai值得到，ai最大值所對應的方案即為最佳。

2 高速列車開行方案評價體系

2.1 評價指標體系的構(gòu)建

高速鐵路是一個集多層次和多功能于一體的綜合服務系統(tǒng)，而旅客列車開行方案的制訂不僅會影響到鐵路部門的經(jīng)濟效益和技術(shù)指標，還會對旅客服務質(zhì)量造成不同程度影響[8]。通過綜合考慮整個鐵路系統(tǒng)內(nèi)部的影響因子，在考慮不同旅客需求情況下對不同列車開行方案進行比選，再結(jié)合列車開行方案特征，通過構(gòu)建因素層、指標層和方案層 3層評價指標體系模型來對整個方案進行綜合評價。針對高速鐵路旅客列車開行方案比選特征，構(gòu)建高速鐵路列車開行方案評價指標體系如圖1 所示。

圖1 高速鐵路列車開行方案評價指標體系

2.2 評價指標的量化

根據(jù)集對分析原理對所選的 20 個指標計算同一度值，從而得到所選方案中各指標值，進一步將初始不確定性指標值轉(zhuǎn)化成可衡量指標值。

（1）效益評價指標 U1。①列車客票收入 U11。該指標指評價范圍內(nèi)列車的開行票價總收入，即其中，fr為第 r 類旅客列車的頻率；Qr為所選旅客列車的固定乘客數(shù)；λr為旅客列車的上座率；Pr為旅客列車的票價率；Lr為第 r 類旅客列車運行的實際距離。②列車開行成本 U12。該指標指評價范圍內(nèi)列車的開行總成本即 C成本=其中，C1為固定成本，元/車組；C2為變動成本，元/ (車組·km)；nr為各個子方案的動車組數(shù)，組；Kr為列車總運行距離，km。

（2）能力評價指標 U2。①列車發(fā)送能力 U21。該指標指評價范圍內(nèi)列車能夠提供的座位總數(shù)，②列車運送能力 U22，人。該指標指選定方案中全部列車運送旅客量 Q運，采用旅客座位數(shù)與客座復用率的乘積計算。③客座周轉(zhuǎn)量U23，人公里。該指標指評價范圍內(nèi)列車開行時及列車滿員條件下所有旅客的位移量：④旅客周轉(zhuǎn)量 U24，人公里。該指標指評價范圍內(nèi)列車開行時旅客的位移量 M旅客，可以根據(jù)客流量統(tǒng)計結(jié)果進行分配。⑤列車平均運距 U25。該指標指所有列車平均行駛的距離，即所有列車行駛總距離與所有列車對數(shù)的比值⑥旅客平均運距 U26。該指標指評價范圍內(nèi)平均運送每位旅客的距離，km，L旅客= M旅客/Q旅客，其中 Q旅客為實際運送旅客數(shù)，人。⑦客座利用率 U27。該指標指平均每一客座公里的人公里數(shù)，即旅客周轉(zhuǎn)量與客座周轉(zhuǎn)量之比，L客座= M旅客/M客座×100%。⑧平均上座率 U28。該指標指列車運行時的總客流量與列車能提供的總座位數(shù)之比，而平均上座率指在選定方案的評價范圍內(nèi)全部區(qū)段上座率的平均值其中，εk為選定方案中第 k 個區(qū)段的客流密度；δr是0-1變量，當列車在 k 區(qū)段運行時 δr= 1，反之 δr= 0。

（3）技術(shù)評價指標 U3。①平均技術(shù)速度 U31。該指標指評價范圍內(nèi)的中途平均速度 (計算時間不包括中間站停車時間)，其中，為所有列車總運行時間，min；為所有列車總的起停車附加時間，min。②平均旅行速度 U32。該指標指選定方案中評價范圍內(nèi)所有列車運行速度的平均值，即。其中，為所有列車在沿途中間站的停站總時間，min。③車底需求數(shù) U33。該指標指評價體系下所有開行列車總的運用車底數(shù)，即總周轉(zhuǎn)時間與每列車每晝夜平均運營時間之比，即其中，為第 r 類列車在折返站總停留時間；為第 r 類列車的平均旅行速度；σ 為平均每個車底每晝夜的運營時間，h，參考值取 σ = 13。④車底走行公里 U34。該指標指在選定的一組方案中所有列車平均行駛的距離，km，即選定方案中所有列車行駛的總距離與所有列車底組數(shù)之比⑤車底周轉(zhuǎn)時間 U35。該指標指在選定的一組方案中每進行 1 次周轉(zhuǎn)而消耗時間的平均值，其中，為估算總運行時間，min；為估算總停站時間，min。

（4）服務評價指標 U4。①列車開行對數(shù) U41。該指標指開行方案中的旅客列車對數(shù)，由開行方案直接統(tǒng)計得到。②平均停站距離 U42。該指標指在行駛途中，所有列車行駛的總距離與總停站次數(shù)的比值，其中，為各停站列車(不包括一站直達式列車) 的總停站次數(shù)。③旅行時間損失 U43。該指標指在旅途中由于列車停靠造成的損失時間總和，根據(jù)客流量的統(tǒng)計結(jié)果計算得到，…，(N-1)。其中，N 為停站數(shù)；Wi為列車在車站 i 的停站時間，min，參考值取 3；Pir為列車在車站 i 停站時的車上旅客數(shù)，人。④不停車直達旅客比率 U44。該指標指途中不停站直達目的地而不需要在中途站點停止的旅客人數(shù)與整個旅途中總客流OD 的比例，結(jié)合各個客流量的分配結(jié)果得到。⑤不換乘直達旅客比率 U45。該指標指直達目的地而不需要換乘的旅客數(shù)與整個旅途中總旅客數(shù)之比，可以根據(jù)客流分配結(jié)果統(tǒng)計得到。

3 算例

以某高速鐵路為研究對象，分別以開行成本、旅客出行便捷程度及兩者線性加權(quán)為目標，以滿足客流出行需求為基本約束，根據(jù)實際客流數(shù)據(jù)和成本參數(shù)經(jīng)過分析計算得到 3 組旅客列車開行方案的備選計劃，再應用集對分析多方案評價法對 3 組備選方案進行綜合評價，得到最適宜的一組方案。各備選方案的 20 個評價指標計算結(jié)果如表2 所示。

由表2 構(gòu)造決策矩陣 H 如下。

表2 各備選方案評價指標

通過文獻檢索和問卷調(diào)查確定效益型指標包括列車客票收入、列車發(fā)送能力、列車運送能力、客座周轉(zhuǎn)量、旅客周轉(zhuǎn)量、列車平均運距、旅客平均運距、客座利用率、平均上座率、平均技術(shù)速度、平均旅行速度、車底走行公里、不停車直達旅客比率及不換乘直達旅客比率等指標，成本型指標包括列車開行成本、車底需求數(shù)、車底周轉(zhuǎn)時間、列車開行對數(shù)、平均停站距離及旅行時間損失等指標，從而得到理想方案矩陣 A0如下。

A0= [11 484.5 3 286.4 340 902 … 0.83]

計算得到同一度矩陣 Q 為

根據(jù)高速鐵路旅客列車開行方案的運作特性進行綜合調(diào)查分析，并且向有關(guān)專家進行詳細咨詢，利用應用層次分析法得到一組可以判斷矩陣的特征向量 W = [0.063 0.071 0.010 … 0.075]。計算特征根 λmax= 11.243，隨機一致性比率 CR = 0.016 8＜0.1，均滿足一致性要求。最后得到備選方案的帶權(quán)同一度矩陣 R = [0.914 0.923 0.888]。綜上所述，由帶權(quán)同一度矩陣得到不同的備選方案，判斷比較方案的優(yōu)越性：方案Ⅱ＞方案Ⅰ＞方案Ⅲ，即方案Ⅱ為最佳方案。

4 結(jié)束語

以某高速鐵路旅客列車開行方案比選為例，根據(jù)多方案評價方法，采用集對分析理論構(gòu)建一套列車開行計劃的綜合評價指標體系，確定各指標的權(quán)重與取值的計算方法，構(gòu)建同一度決策矩陣模型，將定量指標與定性指標量綱化，最后進行統(tǒng)一分析以降低衡量主觀意見偏差，在確保模型簡潔、準確性的同時提高研究結(jié)果可靠性，使之更加符合高速鐵路實際情況。

[1] Chang Y，Yeh C，Shen C. A Multi-objective Model for Passenger Train Services Planning：Application to Taiwan’s High-speed Rail Line [J]. Transportation Research Part B，2000(34)：91-106.

[2] 陶思宇，查偉雄. 旅客列車開行方案經(jīng)濟效益的評價方法[J]. 中國鐵路，2005(8)：42-45. TAO Si-yu，ZHA Wei-xiong. The Evaluation Method of Economic Benefit for Passenger Trains’ Operation Plan [J]. China Railway，2005(8)：42-45.

[3] 李 劍，劉會林. 鐵路旅客列車開行方案經(jīng)濟效益評價系統(tǒng)的開發(fā)[J]. 鐵道運輸與經(jīng)濟，2008，30(5)：26-28. LI Jian，LIU Hui-lin. The Development of Economic Benefit Evaluation System for Railway Passenger Trains’ Operation Plan [J]. Railway Transport and Economy，2008，30(5)：26-28.

[4] 霍黎明，田志強，保魯昆. 旅客列車綜合效益評價方法的研究[J]. 鐵道運輸與經(jīng)濟，2009，31(10)：15-18. HUO Li-ming，TIAN Zhi-qiang，BAO Lu-kun. A Research about the Method of Comprehensive Benefits Evaluation for Passenger Trains[J]. Railway Transport and Economy，2009，31(10)：15-18.

[5] 馮 鑫，羅金寶. 一種基于經(jīng)濟效益的高速鐵路旅客列車開行方案研究[J]. 蘭州大學學報：社會科學版，2014(2)：139-143. FENG Xin，LUO Jin-bao. A Study of Operation Plan about the High-speed Passenger Trains based on an Economic Benefits[J]. Lanzhou University：Social Sciences，2014(2)：139-143.

[6] 趙克勤. 集對分析及其初步應用[M]. 杭州：浙江科學技術(shù)出版社，2000. ZHAO Ke-qin. Set Pair Analysis and Its Preliminary Application [M]. Hangzhou：Zhejiang Science & Technology Press，2000.

[7] 付 強，王慈光，王 丹. 快捷貨物列車開行方案評價指標體系與方法[J]. 鐵道貨運，2009，27(4)：12-15. FU Qiang，WANG Ci-guang，WANG Dan. The Evaluation Index System and Method of Express Freight Train Plan[J]. Railway Freight，2009，27(4)：12-15.

[8] 鄧連波，史 峰. 旅客列車開行方案評價指標體系[J]. 中國鐵道科學，2006，27(3)：106-110. DENG Lian-bo，SHI Feng. The Evaluation Index System of Passenger Train Plan[J]. China Railway Science, 2006，27(3)：106-110.

責任編輯：吳文娟

Comparison on High-speed Railway Train Operating Schemes based on Set Pair Analysis

Based on expounding the multi-scheme evaluation method of set pair analysis (SPA), targeting with the 3-layers evaluation index system of comparison status of high-speed railway train operating scheme including establishment factor layer, index layer and scheme layer, according to SPA principle, 20 evaluation indices are selected, and the detail concept of each index and determination method of the index value are provided. In the end, taking certain high-speed railway as an example, through taking comprehensive evaluation on 20 indices in 3 sets of alternative schemes, the decision matrix is established, then according to the principle of identity degree, the decision matrix of identity degree is calculated and achieved, the occupied weight of each index in the decision matrix is determined by using AHP, and the weighted identity matrix is further calculated, and then, through comparing the corresponding value of each alternative scheme in the weighted identity matrix, the optimized implementation scheme of high-speed railway is achieved.

Operating Scheme; Comprehensive Evaluation Index; SPA; Weighted Identity Matrix

1003-1421(2015)09-0013-05

U292.3+5

A

10.16668/j.cnki.issn.1003-1421.2015.09.03

2015-05-25

中國鐵路總公司科技研究計劃項目(2014X004-D，2015X008-B)