中國鐵路快捷貨物運輸網絡復雜性分析

馬 亮,金福才,胡宸瀚

(1.西南交通大學信息科學與技術學院,成都 611756;2.國家鐵路智能運輸系統工程技術研究中心,北京 100844;3.中國鐵道科學研究院集團有限公司電子計算技術研究所,北京 100081)

0 引言

地震、風災、滑坡等自然災害的發生和線路、信號、供電等關鍵設備的故障可能會降低鐵路運輸效率和通過能力[1],甚至會發生大面積車流擠壓、列車晚點和運行取消等情況[2]。為降低鐵路事故的發生頻率和嚴重度,需要對鐵路網絡進行建模,并分析鐵路網絡的靜態拓撲特性和動態傳播特性,為鐵路的日常維護和應急管理提供支持。鐵路網絡是規模龐大、復雜度高、受外界影響大的巨型動態復雜網絡系統[3]?；趶碗s網絡建模和分析方法,可以從不同的邏輯角度分析鐵路網絡的靜態拓撲特征和動態演化機理,為解決鐵路網絡發展規劃、運營保障和運輸組織優化等問題提供理論依據。隨著高鐵快運、跨境聯運等多種新型快運業務的不斷增長,對鐵路快捷貨物運輸網絡進行復雜性分析有助于提高鐵路運營維護和應急管理水平,是鐵路運輸領域的關鍵科學問題。

針對鐵路復雜網絡建模和拓撲結構研究:胡宸瀚等[4]歸納了P,R,L等3種空間下的鐵路網絡建模方法和特性分析方法;馮佳等[5]在L空間構建了城市軌道交通流網絡模型,分析了小世界特性與交通流特性;WANG Li等[6]建立了鐵路物理網絡和鐵路車流網絡,基于網絡效率貢獻度指標對網絡中的關鍵站點進行了定量分析;王小榮[7]和詹斌等[8]分別構建了鐵路快速貨運物理網絡和L空間的無權網絡,分析了網絡的小世界、無標度和節點中心等拓撲特性。針對鐵路網絡魯棒性、脆弱性或者抗毀性分析研究:Yin等[9]提出了地震災害高發區的地理鐵路網的魯棒性評估方法;Gu等[10]提出了基于聚集社團指標的高速鐵路網絡魯棒性評價方法;于寶等[11]通過計算網絡效率和最大連通子圖得到高速鐵路網絡在隨機攻擊下和蓄意攻擊下均比較脆弱;劉朝陽等[12]仿真計算了城市軌道交通網絡在不同失效策略下的級聯失效抗毀性。上述研究取得了較好的成果,然而已有的復雜網絡模型中較少考慮貨運車流量對車站間連接緊密程度的刻畫,導致模型和分析的結論脫離實際?，F有無權網絡復雜性分析方法已經不適用于加權網絡,單一的重要度評價指標與魯棒性分析方法無法全面分析網絡的復雜性,需要繼續完善拓撲性能指標、節點的重要度評價指標和魯棒性分析方法。

本文選取中國鐵路快捷貨運系統網絡(Chinese Railway Express Freight Transportation Network,CREFTN)為研究對象,提出L空間下中國鐵路快捷貨運系統的無向加權復雜網絡模型(Undirected and Weighted Complex Network Model of CREFTN,UWCN-CREFTN)構建方法、節點重要度綜合評價方法和基于綜合重要度節點失效的網絡魯棒性分析方法。依據實際的中國鐵路快捷貨運時刻表和車站間日均交互車流量數據,客觀全面地定量分析網絡的靜態拓撲特性和節點重要性;基于加權網絡效率指標,定量對比分析節點隨機失效和蓄意破壞2個方面4種情形下網絡的魯棒性;最后,給出對網絡魯棒性影響較大的關鍵車站列表。

1 中國鐵路快捷貨運加權復雜網絡建模

鐵路物理網絡由車站和線路構成,鐵路快捷貨運系統是在鐵路物理網絡的基礎上,通過快捷貨運組織技術實現零散貨物、批量散貨、行包、高鐵、聯運等多種制式的快捷運輸,是鐵路物理網絡和快捷貨運列車之間耦合產生的邏輯網絡,包括P空間、R空間和L空間等3種邏輯復雜網絡模型[1]。為研究CREFTN結構對貨物快捷運輸效率的影響,在鐵路快捷貨運列車成對開行和車站間交互車流均衡的前提下,在L空間建立CREFTN的加權復雜網絡模型。

令鐵路物理網絡為G=(P,E),其中P={p1,p2,…,pN}為鐵路貨運車站集合,N為車站數量,E={eij|i,j=1,2,…,N}為邊集合,若車站pi和pj之間存在線路,則eij=1,否則eij=0。令鐵路快捷貨運列車集合為T={t1,t2,…,tM},M為列車數量,任意列車tk=(Pk,Ck),Pk∈P為列車tk停靠并辦理貨車中轉的車站集合,Ck為列車tk在車站辦理貨運中轉車輛數集合。通過鐵路物理網絡與鐵路快捷貨運列車耦合得到L空間UWCN-CREFTN模型記為三元組[1-2]:GW=(P,L,W)=G⊕T,耦合關系“⊕”表示:若G中的兩個車站至少是同一趟快捷貨運列車連續?？壳肄k理貨運中轉業務的車站,則這兩個車站之間存在一條邊。L={lij|i,j=1,2,…,N}為表征節點間耦合連接關系的鄰接矩陣,若車站pi和pj之間存在邊,則lij=1,否則lij=0。W={wij|i,j=1,2,…,N}為邊的權值集合,元素wij為車站之間日均交互車輛數。若任意兩個節點wij=lij時,則為無權網絡。如圖1所示為UWCN-CREFTN模型建立過程示意圖,模型中線條寬度與車站間交互的貨運量成正比。由于p1和p4,p3和p4,p4和p6等節點對之間不存在車流耦合關系,因此在模型中這些節點之間的邊被約減。由于至少有2趟列車連續?？寇囌緋2和p4,p5和p7,因此模型中車站p2和p4,p5和p7間的連線較寬。

圖1 鐵路物理網絡與鐵路快捷貨運列車耦合示意圖

通過上述建模方法,以“中國鐵路95306網站”上公開的161列快運貨物班列和362列直達貨物班列時刻表,以及車站間日均交互車流量數據為基礎,在L空間建立的UWCN-CREFTN模型如圖2所示,網絡包含了350個節點和468條邊,圖中紅色點表示承擔快捷貨運中轉業務的車站,藍色連線表示相鄰車站之間有快捷貨運列車停留并有車流中轉業務,連線的寬度與邊的權值成正比。由圖2可見鐵路承擔快捷貨運的車站較多分布在中國環渤海、長江中游、長三角、成渝、珠三角和中原等城市群或者經濟群,與這些區域的經濟活躍和鐵路貨運需求程度相關。

圖2 在L空間下UWCN-CREFTN模型結構圖

2 鐵路快捷貨運網絡復雜性分析方法

2.1 加權復雜網絡性能指標

2.2 加權復雜網絡綜合重要度指標

2.3 基于效率的加權復雜網絡魯棒性分析方法

1)節點失效或者遭受蓄意破壞是指車站喪失貨運車流中轉能力,保留列車接發和通過能力。節點失效之后只影響模型中相鄰節點間的車流交互,并不影響其他節點間的車流交互。若在節點被移除之后出現部分剩余節點之間路徑不連通時,那么不連通節點之間的效率記為0。

2)節點隨機失效下網絡的魯棒性:隨機從網絡中移除節點,并移除與該節點相連的所有邊,計算Ri,直至網絡中所有的節點被移除或者效率累積下降率達到100%。為降低隨機性,重復多次上述過程,計算節點失效率對應Ri的平均值。

3)從最大加權介數、最大強度和最大綜合重要度三種情況分析節點受到蓄意破壞時網絡的魯棒性:對網絡中節點按照加權介數(強度、綜合重要度)從大到小進行排序,依次移除網絡中加權介數(強度、綜合重要度)最大的節點,并移除與該節點相連的所有邊,重新計算網絡的效率Ri,直至網絡中所有的節點被移除或者效率下降率達到100%。

3 鐵路快捷貨運網絡復雜性分析

3.1 加權復雜網絡模型的拓撲特性分析

圖3 無權網絡和加權網絡距離概率分布曲線

通過計算得到:網絡中所有點的平均強度為105.906,表明每日平均有105.906輛貨車到達這些車站,圖4所示為UWCN-CREFTN的強度分布曲線圖,圖5所示為強度區間內車站數分布柱狀圖,從圖4和圖5可以看出只有少部分車站的強度較大,絕大部分車站的強度小于350。

圖4 點強度分布曲線圖

圖5 強度區間內車站數分布柱狀圖

圖6所示為節點累積強度分布曲線,經擬合得到具有典型的冪率分布R(s≥si)∝e-0.015 7si,表明:UWCN-CREFTN具有高度的異質性,只有少數車站在貨物快捷運輸中起到了關鍵作用,這些車站具有較強的車流處理能力。通過計算得到節點的平均強度與度之間的關系曲線如圖7所示,經曲線擬合得到節點的平均強度與節點度之間存在冪律關系:〈si〉∝ki2.372 7,表明:度數越大的節點強度越大,強度的增長快于度的增長,即車站服務能力的增長快于連通度的增長。由此可以說明:在當前鐵路網絡基礎設施條件下,可以通過提高冪指數的取值來提高鐵路快捷貨運的服務能力,但是物理網絡中區間通過能力和車站作業能力的限制使得冪指數的增大受限,需要通過運輸組織策略調整和路網結構優化共同來提高運輸效率,譬如:調整裝車流、優化列車編組計劃、優化運行徑路、提高瓶頸區段和車站的能力等,實現車流動態分布與路網運輸能力的最優匹配。

圖6 節點累積強度分布曲線

圖7 點度數與平均強度關系曲線

通過計算得到:無權網絡平均聚集系數C為0.035 3,加權網絡聚集系數CW為0.034 7,盡管遠小于規則網絡聚集系數0.5,但是比隨機網絡聚集系數0.002 85要大得多,介于規則和隨機之間,表明UWCN-CREFTN具有小世界特征[1],即網絡中任意兩個節點間都可以找到一條相對較短的路徑,且節點之間存在較強的聚集性。從圖8中聚集系數與度數關系的曲線可得:絕大多數C(k)>CW(k),也就是網絡的拓撲聚集取決于權低的邊,意味著具有較強車流處理能力的車站偏好與相關車站之間有更多的車流交互,這可以解釋為:比較重要的車站形成穩定的集團就能處理更多的車流。

圖8 網絡平均聚集系數與度數關系曲線

通過計算得到均等性Yi與點度數ki之間關系曲線如圖9所示,經曲線擬合得到存在冪律關系:Yi∝ki-0.864 37,接近于ki-1。說明:隨著度數的增大,點的不等性下降較快,邊的權值對節點的強度貢獻度比較均衡,即車站與其相鄰車站之間交互車流數比較均衡。

圖9 點度數與點權不等性關系曲線

通過計算得到最近鄰點平均度與最近鄰點加權平均度的關系曲線如圖10所示,從圖10可以得到:絕大多數強度較大的節點偏好與度數較小的節點相連,節點之間相互作用屬于加權非同類匹配性,即絕大多數節點偏好于與度數較小節點連接。通過計算得到節點強度與鄰接節點平均強度之間的關系曲線如圖11所示,從圖11可以得到:不同距離的鄰接節點平均強度與節點強度分布趨勢比較一致,節點間存在強度-強度負相關性,強度較小的節點偏好與強度較大的鄰接節點相連,強度較大的節點偏好與強度較小的鄰接節點相連。

圖10 節點間同類匹配關系曲線

圖11 節點強度與鄰接節點平均強度之間關系曲線

3.2 加權復雜網絡模型的節點重要性分析

表1 綜合重要度較大的20個車站表

圖12 前50個車站對應重要度對比分析曲線圖

圖13 綜合總要度較大的前100個車站的分布圖

從強度、加權介數和綜合重要度3個指標對比分析得到:1)比較重要的車站基本上是大型的技術站和邊境口岸站,且集中分布在中國中東部地區,這與相關區域的經濟活躍程度和貨運需求旺盛正相關,而部分出入境口岸站的重要性大與辦理大量的國際聯運業務相關。2)絕大多數車站的重要度比較小,少數比較重要的車站成為“中心節點”,部分比較重要的車站之間相互吸引形成集團,這些高等級車站的集中設置,為客戶辦理貨運業務提供了更多的選擇。

3.3 加權復雜網絡模型的魯棒性對比分析

圖14 基于隨機失效的網絡魯棒性變化曲線

圖15 三種蓄意破壞下網絡魯棒性變化曲線

從圖15對比分析得到:網絡效率下降率與3種破壞的車站比例近似呈指數分布;最大加權介數和最大綜合重要度破壞的車站比例為7%,最大強度破壞的車站比例為9.5%時,網絡效率下降達到90%;最大強度和最大綜合重要度破壞的車站比例達到38%,最大加權介數破壞的車站比例達到99.4%時,網絡接近完全崩潰。因此,最大加權介數破壞和最大綜合重要度破壞比最大強度破壞使網絡變得更脆弱,最大加權介數破壞比最大綜合重要度破壞更甚,最大強度破壞和最大綜合重要度破壞比最大加權介數破壞更快地造成網絡崩潰。

排序前20的車站依次受到最大強度、最大加權介數和最大綜合重要度破壞之后,網絡的累積效率下降了69.4%,85.1%和85.0%。表2為綜合重要度最大的車站依次受破壞后網絡效率累積下降率表,這些車站需要重點保護,一旦遭受破壞或者故障會對快捷貨運網絡的運輸效率產生較大的影響。

表2 對網絡魯棒性影響較大的關鍵車站

4 結論

本文提出了中國鐵路快捷貨運無向加權復雜網絡建模方法、節點重要度綜合評價方法和基于節點失效的網絡魯棒性分析方法。通過實例驗證表明:與無權網絡建模與分析方法相比,本文方法能夠更全面、客觀和準確地分析鐵路快捷貨運系統網絡的復雜性。在本文研究的基礎上,以提高網絡的魯棒性、連通性和效率為目標,通過對比分析新建改造鐵路基礎設施和調整鐵路貨運組織方案等措施對網絡性能的影響,進而為中國鐵路貨運系統網絡規劃和組織優化提供決策支持。