楊　銳,牟瑞芳,閆海峰

(西南交通大學(xué) 交通運(yùn)輸與物流學(xué)院，四川 成都　610031)

列車運(yùn)行計(jì)劃及列車運(yùn)行調(diào)整計(jì)劃(以下簡(jiǎn)稱計(jì)劃)是鐵路組織日常運(yùn)輸生產(chǎn)的核心和依據(jù)，對(duì)其實(shí)施過(guò)程進(jìn)行仿真，既可經(jīng)濟(jì)、安全地檢驗(yàn)所編制計(jì)劃的可執(zhí)行性及合理性，快速、方便地對(duì)其動(dòng)態(tài)性能進(jìn)行評(píng)價(jià)，為多計(jì)劃比選提供參考，又可進(jìn)一步分析和優(yōu)化計(jì)劃。

現(xiàn)有針對(duì)列車運(yùn)行仿真的研究[1-13]，大多根據(jù)運(yùn)行中列車的受力情況、列車運(yùn)行控制的過(guò)程或設(shè)備作用原理等進(jìn)行建模和設(shè)計(jì)仿真系統(tǒng)，較少考慮包括人為因素在內(nèi)的多種隨機(jī)擾動(dòng)對(duì)列車運(yùn)行的影響。因此，本文從列車運(yùn)行實(shí)際效果出發(fā)，考慮主、客觀因素相互疊加造成的列車運(yùn)行偏移，同時(shí)簡(jiǎn)化列車運(yùn)行過(guò)程描述，建立能夠快速、完整模擬計(jì)劃實(shí)施過(guò)程的仿真模型，用于支撐列車運(yùn)行仿真平臺(tái)的開(kāi)發(fā)。

1　建模工具的選擇

計(jì)劃在實(shí)施過(guò)程中，由于通常會(huì)有多列車同時(shí)處于運(yùn)行狀態(tài)，因此具有典型的并發(fā)特性；計(jì)劃所涉路網(wǎng)，由于列車所經(jīng)車站和區(qū)間(本文特指區(qū)間的單條線路)具有相對(duì)穩(wěn)定的拓?fù)浣Y(jié)構(gòu)，因此具有網(wǎng)絡(luò)特性；列車在運(yùn)行過(guò)程為實(shí)現(xiàn)位移，會(huì)不斷申請(qǐng)、占用和釋放不同的運(yùn)行資源，因此具有明顯的離散動(dòng)態(tài)系統(tǒng)特性。

鑒于顏色Petri網(wǎng)是描述離散并發(fā)系統(tǒng)動(dòng)態(tài)行為過(guò)程的優(yōu)良工具，又是高級(jí)網(wǎng)系統(tǒng)，且使用元素較少；而賦時(shí)Petri網(wǎng)則是在傳統(tǒng)Petri網(wǎng)的基礎(chǔ)上增設(shè)了時(shí)間參數(shù)。因此，本文綜合顏色Petri網(wǎng)和賦時(shí)Petri網(wǎng)的優(yōu)勢(shì)，定義多元組∑(P,T,F,B,I,O,C,ξ,M0)模擬計(jì)劃的實(shí)施過(guò)程。其中，P為庫(kù)所的有限集，描述列車運(yùn)行環(huán)境；T為變遷的有限集，對(duì)應(yīng)與列車運(yùn)行過(guò)程有關(guān)的各事件；F為有向弧的集合，表示各變遷對(duì)應(yīng)事件發(fā)生時(shí)列車運(yùn)行或運(yùn)行信息流動(dòng)的方向；B為無(wú)向弧的集合，表示觸發(fā)各變遷對(duì)應(yīng)事件的數(shù)據(jù)來(lái)源；I和O均為有向弧的權(quán)函數(shù)集合，分別表示各庫(kù)所在相關(guān)變遷對(duì)應(yīng)事件發(fā)生時(shí)失去的和獲得的數(shù)據(jù)，共同構(gòu)成各變遷的變遷規(guī)則；C為顏色的有限集，表示列車及運(yùn)行環(huán)境的各抽象要素；ξ為全局時(shí)鐘，表明模擬的計(jì)劃實(shí)施時(shí)刻；M0為初始標(biāo)識(shí)，描述列車在仿真開(kāi)始前的分布狀況。

2　模型的建立

2.1　列車描述

2.2　運(yùn)行環(huán)境描述

列車運(yùn)行環(huán)境包括所有圖定列車經(jīng)過(guò)的車站和區(qū)間。對(duì)于容納列車的車站和區(qū)間，在Petri網(wǎng)中均可用能容納資源的庫(kù)所表示，但由于車站和區(qū)間的功能不同而使得抽象的要素和所需存儲(chǔ)的參數(shù)也不同，因此，需以不同形式的庫(kù)所分別表示車站和區(qū)間。

2.2.1車站

2.2.2區(qū)間

對(duì)于部分存在渡線或越行線的非單一串聯(lián)形式的區(qū)間，可將渡線道岔區(qū)域或越行線區(qū)域抽象為特殊車站，將單區(qū)間拆分為多個(gè)區(qū)間后，再分別進(jìn)行描述。

2.2.3列車源站和匯站

2.3　運(yùn)行過(guò)程描述

列車運(yùn)行過(guò)程可用一系列離散事件描述，同時(shí)，限制列車運(yùn)行行為的運(yùn)行規(guī)則所依賴的動(dòng)態(tài)運(yùn)行信息也可由不同事件驅(qū)動(dòng)改變。由于不同事件所造成的列車位置或數(shù)據(jù)變化以及觸發(fā)的條件均不相同，而在Petri網(wǎng)中任意物質(zhì)和信息的流動(dòng)或改變均靠變遷實(shí)現(xiàn)，因此列車運(yùn)行過(guò)程需采用多種變遷共同描述。

2.3.1列車出站

圖1　列車出站變遷

{ns[ns+1],yks|k∈Rsj∧k≠ksj}

(1)

(2)

(3)

其中，

K=ksj

a=gK

a′=gK+

式中：sK+表示列車K前進(jìn)方向上位于s站的前方車站；Rj表示區(qū)間j內(nèi)的列車集合，|Rj|則對(duì)應(yīng)區(qū)間j內(nèi)的列車數(shù)；而N(j)和β(j)均為自定義函數(shù)， 分別返回符合μj中參數(shù)所述的正態(tài)分布隨機(jī)數(shù)和β分布隨機(jī)數(shù)。

(3) (|Rj|=0)∧(tξ≥(dj+τs))。

2.3.2列車進(jìn)站

K′}

(4)

{ns[ns-1],yK′s[0]}-

(5)

(6)

(2)tξ≥dj+τs。

2.3.3到站狀態(tài)修改

(7)

2.3.4股道安排及調(diào)整

(8)

式中：K2表示在s站未安排作業(yè)股道的所有即將到站列車中預(yù)計(jì)最早到達(dá)的列車。

(9)

2.3.5列車產(chǎn)生

(10)

(11)

(12)

2.3.6列車消失

{ns[ns+1]}

(13)

(14)

式中：列車K4表示正在從s站進(jìn)入PD的列車。

2.3.7發(fā)車準(zhǔn)備

(15)

其中ksj的取值規(guī)則需參照車站其他連接方向上區(qū)間內(nèi)的列車運(yùn)行情況、車站內(nèi)的列車作業(yè)情況和被越行情況等，以及所采用的調(diào)度調(diào)整策略制定，并且要求在出現(xiàn)列車運(yùn)行偏離計(jì)劃較大時(shí)能根據(jù)規(guī)則適當(dāng)調(diào)整s站發(fā)往區(qū)間j的列車的順序。

2.4　仿真時(shí)間步選定

從以上列車運(yùn)行相關(guān)事件所造成列車位置和運(yùn)行信息的變化以及各事件的觸發(fā)條件可以看出，全局時(shí)鐘不僅表征系統(tǒng)狀態(tài)所對(duì)應(yīng)的時(shí)刻，也驅(qū)動(dòng)系統(tǒng)運(yùn)行。因此，為能按照列車運(yùn)行實(shí)際過(guò)程順序觸發(fā)各事件，應(yīng)將全局時(shí)鐘ξ定義成一個(gè)可不斷遞增的動(dòng)態(tài)變量，而變量的增幅即仿真時(shí)間步長(zhǎng)Δtξ因?yàn)榕c仿真效率呈正比關(guān)系、與仿真度呈反比關(guān)系，故而應(yīng)均衡考慮仿真效率及仿真度高低，在兼顧仿真對(duì)象及各相關(guān)時(shí)間參數(shù)表示精度的基礎(chǔ)上合理選定。考慮到計(jì)劃中有關(guān)列車在沿途各站到發(fā)或通過(guò)的時(shí)刻以及各技術(shù)作業(yè)時(shí)間標(biāo)準(zhǔn)基本上都只精確到分鐘，而各車站間隔時(shí)間卻通常精確到秒，因此，本文推薦將全局時(shí)鐘的增幅設(shè)置為秒，即Δtξ=1 s。

2.5　標(biāo)識(shí)初始化

(16)

式中參數(shù)可根據(jù)列車運(yùn)行線路細(xì)分為本線和跨線列車分別予以設(shè)定。

至此，針對(duì)計(jì)劃實(shí)施過(guò)程的仿真建模已全部完成，通過(guò)借助時(shí)間步推進(jìn)驅(qū)動(dòng)所創(chuàng)建的Petri網(wǎng)系統(tǒng)不斷運(yùn)行，即可實(shí)現(xiàn)對(duì)所選計(jì)劃在指定時(shí)間段內(nèi)的實(shí)施過(guò)程進(jìn)行仿真。模型雖然在確定列車實(shí)際區(qū)間運(yùn)行時(shí)分時(shí)沒(méi)有單獨(dú)考慮起停車附加時(shí)分，但是鑒于因停車而導(dǎo)致增加的起停車附加時(shí)分已在編制計(jì)劃時(shí)考慮，計(jì)劃外停車是隨機(jī)擾動(dòng)的結(jié)果或本身就是隨機(jī)擾動(dòng)的范疇，由此增加的起停車附加時(shí)分屬于運(yùn)行偏移的一部分，而區(qū)間運(yùn)行時(shí)分標(biāo)準(zhǔn)并非列車所能達(dá)到的極限，通過(guò)搶點(diǎn)提速，還能部分甚至全部抵消計(jì)劃外停車所增加的起停車附加時(shí)分，因此不會(huì)對(duì)最終的仿真結(jié)果造成太大偏差。

3　仿真實(shí)現(xiàn)及實(shí)例

雖然模型基于改進(jìn)Petri網(wǎng)，但由于用以描述與列車運(yùn)行有關(guān)各事件的變遷采用了自定義的發(fā)生規(guī)則及觸發(fā)條件，因此無(wú)法使用現(xiàn)有Petri網(wǎng)軟件所構(gòu)建的網(wǎng)系統(tǒng)。為此，本文先基于.NET平臺(tái)，以所述模型為核心，輔以基礎(chǔ)數(shù)據(jù)獲取、仿真設(shè)置以及仿真結(jié)果顯示等相關(guān)模塊，運(yùn)用C#語(yǔ)言自主開(kāi)發(fā)了列車運(yùn)行計(jì)劃仿真平臺(tái)；而后以武廣高鐵為背景，參考文獻(xiàn)[15]，假設(shè)單區(qū)間實(shí)際運(yùn)行時(shí)分偏離標(biāo)準(zhǔn)的范圍：針對(duì)任意區(qū)間，在仿真的過(guò)程中，如果列車正點(diǎn)或早點(diǎn)進(jìn)入該區(qū)間，則列車在該區(qū)間的實(shí)際運(yùn)行時(shí)分最少不得少于標(biāo)準(zhǔn)的77.5%，最多不得多于標(biāo)準(zhǔn)的142.5%；如果列車晚點(diǎn)進(jìn)入該區(qū)間，則列車在該區(qū)間的實(shí)際運(yùn)行時(shí)分最少不得少于標(biāo)準(zhǔn)的77.8%，最多不得多于標(biāo)準(zhǔn)的122.8%。單項(xiàng)技術(shù)作業(yè)完成時(shí)間偏離標(biāo)準(zhǔn)的范圍取值亦類似。

選用2011-11-30的列車運(yùn)行圖，基于Win7系統(tǒng)進(jìn)行了50次全天仿真，每次仿真均在29 s內(nèi)完成，其中10次仿真結(jié)果的檢驗(yàn)統(tǒng)計(jì)數(shù)據(jù)如表1所示，而由系統(tǒng)自動(dòng)生成的其中1次仿真結(jié)果與計(jì)劃對(duì)比的局部如圖2所示。

表1　仿真結(jié)果檢驗(yàn)統(tǒng)計(jì)表(a)

表1　仿真結(jié)果檢驗(yàn)統(tǒng)計(jì)表(b)

圖2　某次仿真與計(jì)劃對(duì)比局部圖

眾多仿真結(jié)果顯示：包括動(dòng)檢列車在內(nèi)的205列圖定列車均能沿指定徑路完成運(yùn)行，且跨線列車DJ5536/5均能在廣州北站正常上下線；各次仿真得到的列車平均旅速與計(jì)劃的255.723 km·h-1沒(méi)有較大偏差；各次仿真全過(guò)程均未出現(xiàn)因車站能力不足而導(dǎo)致列車計(jì)劃未完成的情況；所有列車到站、出站及追蹤全部滿足各自的間隔時(shí)間；部分仿真出現(xiàn)列車計(jì)劃內(nèi)停站的時(shí)間略小于標(biāo)準(zhǔn)，這可視為為趕點(diǎn)而組織旅客快速乘降或乘降旅客人數(shù)較少所致；相比于計(jì)劃運(yùn)行時(shí)分，所有區(qū)間運(yùn)行仿真時(shí)分的減少值均在設(shè)置范圍內(nèi)，而部分區(qū)間運(yùn)行時(shí)分增加值超出了設(shè)置范圍，這是由于高等級(jí)車連續(xù)不停站早點(diǎn)通過(guò)多個(gè)車站后，受追蹤或進(jìn)站限制而降速運(yùn)行所致；各次仿真中列車在各站到發(fā)或通過(guò)的時(shí)刻雖時(shí)常偏離計(jì)劃，但均在合理范圍內(nèi)波動(dòng)，并且大多通過(guò)冗余時(shí)間得到消解，另有小部分雖改變了列車在區(qū)間的運(yùn)行秩序，但從涉及的列車數(shù)也表明影響范圍有限，屬于正常的列車運(yùn)行調(diào)度范疇。因此，可以認(rèn)為仿真結(jié)果基本符合客觀實(shí)際，說(shuō)明所建模型可行、有效。

4　結(jié)　語(yǔ)

基于部件組合思想，在考慮區(qū)間雙向行車需求的同時(shí)，兼顧仿真需求的多樣性而建立的列車運(yùn)行仿真模型可根據(jù)路網(wǎng)結(jié)構(gòu)快速構(gòu)建仿真框架，在對(duì)各靜態(tài)參數(shù)根據(jù)相關(guān)數(shù)據(jù)進(jìn)行初始化后，可按設(shè)置的區(qū)域范圍及時(shí)段要求模擬列車運(yùn)行計(jì)劃。由于該模型屬于介觀模型，仿真速度較快。仿真實(shí)現(xiàn)的過(guò)程也表明，模型直觀，易于編程實(shí)現(xiàn)，仿真結(jié)果能夠較為完整地反映計(jì)劃實(shí)施的全過(guò)程。

[1]張波，馬大煒. 中高速列車共線運(yùn)行的仿真研究[J]. 中國(guó)鐵道科學(xué)，2003, 24(3): 119-124.

(ZHANG Bo, MA Dawei. Simulation of High-Speed and Medium-Speed Trains Sharing the Same Railway Line[J]. China Railway Science, 2003,24(3):119-124. in Chinese)

[2]丁國(guó)富，張衛(wèi)華，劉伯興，等. 列車運(yùn)行動(dòng)態(tài)仿真研究[J]. 系統(tǒng)仿真學(xué)報(bào)，2004, 16(8): 1697-1700.

(DING Guofu, ZHANG Weihua, LIU Boxing,et al. Study on Dynamical Simulation of Running Locomotive[J]. Journal of System Simulation, 2004,16(8):1697-1700. in Chinese)

[3]周磊，余祖俊，史紅梅. 列車運(yùn)行線路仿真系統(tǒng)的研究[J]. 系統(tǒng)仿真學(xué)報(bào)，2004, 16(7): 1463-1466.

(ZHOU Lei, YU Zujun, SHI Hongmei. Research on Train Running Line Simulation System[J]. Journal of System Simulation, 2004,16(7):1463-1466. in Chinese)

[4]劉云. 列車運(yùn)行仿真系統(tǒng)的建模與實(shí)現(xiàn)[J]. 鐵道學(xué)報(bào)，1995，17(增2): 20-26.

(LIU Yun.Modelling and Implementation of a Simulation System for Train Operation[J]. Journal of the China Railway Society,1995，17(Supplement 2):20-26. in Chinese)

[5]毛保華，何天鍵，袁振洲，等. 通用列車運(yùn)行模擬軟件系統(tǒng)研究[J]. 鐵道學(xué)報(bào)，2000, 22(1): 1-6.

(MAO Baohua, HE Tianjian, YUAN Zhenzhou, et al.A General-Purposed Simulation System on Train Movement[J]. Journal of the China Railway Society, 2000,22(1):1-6. in Chinese)

[6]唐金金，周磊山，佟路，等. 單列高速列車運(yùn)行仿真模型與算法[J]. 中國(guó)鐵道科學(xué)，2012, 33(3): 109-115.

(TANG Jinjin, ZHOU Leishan, TONG Lu, et al. Simulation Model and Algorithm for Single High-Speed Train Operation[J]. China Railway Science, 2012,33(3):109-115. in Chinese)

[7]葉陽(yáng)東，程少芬，王旭，等. 基于一種混合Petri網(wǎng)的列車運(yùn)行系統(tǒng)的建模與分析[J]. 鐵道學(xué)報(bào)，2009, 31(5): 42-49.

(YE Yangdong, CHENG Shaofen, WANG Xu, et al. Modeling and Analyzing of Train Operation Systems Based on a Kind of Hybrid Petri Net[J]. Journal of the China Railway Society, 2009,31(5):42-49. in Chinese)

[8]白紫熙，周磊山，唐金金. 基于運(yùn)行圖驅(qū)動(dòng)的列車運(yùn)行控制半實(shí)物仿真系統(tǒng)[J]. 中國(guó)鐵道科學(xué)，2015, 36(4): 139-144.

(BAI Zixi, ZHOU Leishan, TANG Jinjin.Hardware-in-the-Loop Simulation System for Train Operation Control Based on Working Diagram Driven[J]. China Railway Science,2015,36(4):139-144. in Chinese)

[9]周妍，周磊山. 高速鐵路行車調(diào)度指揮一體化仿真實(shí)驗(yàn)平臺(tái)設(shè)計(jì)與研究[J]. 鐵道學(xué)報(bào)，2012，34(6): 1-7.

(ZHOU Yan, ZHOU Leishan. Study on Simulation and Experiment Platform of Integrated High-Speed Railway Traffic Control[J]. Journal of the China Railway Society,2012,34(6):1-7. in Chinese)

[10]王宏剛. 復(fù)線列車運(yùn)行調(diào)度的層次Petri模型及仿真[J]. 系統(tǒng)仿真學(xué)報(bào)，2011, 23(12): 2793-2798, 2804.

(WANG Honggang. Hierarchical Petri Net Model and Simulation for Train Traffic Control on Double Railway[J]. Journal of System Simulation, 2011,23(12):2793-2798, 2804. in Chinese)

[11]劉海東，毛保華，何天健，等. 不同閉塞方式下城軌列車追蹤運(yùn)行過(guò)程及其仿真系統(tǒng)的研究[J]. 鐵道學(xué)報(bào)， 2005, 27(2): 120-125.

(LIU Haidong, MAO Baohua,He Tianjian, et al. Study on Tracking Operations between Trains of Different Block Modes and Simulation System[J]. Journal of the China Railway Society,2005,27(2):120-125. in Chinese)

[12]孟令云，蔣熙，楊肇夏. 基于Agent的客運(yùn)專線列車運(yùn)行仿真研究[J]. 系統(tǒng)仿真學(xué)報(bào)，2009, 21(6): 1538-1542.

(MENG Lingyun, JIANG Xi, YANG Zhaoxia. Study on Train Running Simulation of Railway Passenger Dedicated Line Based on Agent[J]. Journal of System Simulation, 2009,21(6):1538-1542. in Chinese)

[13]郭彤城，慕春棣. ATP鐵路運(yùn)行仿真系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J]. 系統(tǒng)仿真學(xué)報(bào)，1999, 11(4): 293-296.

(GUO Tongcheng, MU Chundi. The Design and Implementation of a Railway Running Simulation System[J]. Journal of System Simulation, 1999,11(4):293-296. in Chinese)

[14]張星臣，胡安洲，楊浩.一種基于有控隨機(jī)與隨機(jī)有控約束的列車運(yùn)行仿真模型[J]. 鐵道學(xué)報(bào),1995,17(3)：21-27.

(ZHANG Xingchen, HU Anzhou, YANG Hao. A Train Operating Simulation Model Based on Constraints of Control in Random Origin and Control in Random Process[J]. Journal of the China Railway Society,1995,17(3):21-27. in Chinese)

[15]張星臣，胡安洲. 列車區(qū)間運(yùn)行時(shí)分偏離的β函數(shù)分布分析[J]. 鐵道學(xué)報(bào),1996，18(3): 1-6.

(ZHANG Xingchen, HU Anzhou. Analysis of β-Function Distribution for Deviation of Train Running Time in Section[J]. Journal of the China Railway Society, 1996,18(3):1-6. in Chinese)