邱 菲
(蘭州鐵路局 客貨運輸統(tǒng)計所,甘肅 蘭州 730000)
基于灰色馬爾可夫鏈組合模型的蘭州鐵路局集裝箱運量預測
邱 菲
(蘭州鐵路局 客貨運輸統(tǒng)計所,甘肅 蘭州 730000)
闡述灰色馬爾可夫鏈組合模型理論,在分析蘭州鐵路局集裝箱運輸現(xiàn)狀及特點的基礎上,以蘭州鐵路局集裝箱運量現(xiàn)狀及歷史數(shù)據(jù)為依據(jù),運用灰色馬爾可夫鏈組合模型對蘭州鐵路局集裝箱運量進行預測分析,為鐵路部門制定物流營銷政策、科學組織貨物運輸提供參考。
鐵路;集裝箱;運量預測;灰色預測;馬爾可夫模型
集裝箱運輸作為我國鐵路運輸重點發(fā)展方向,具有標準統(tǒng)一、方便運輸、環(huán)境友好等優(yōu)勢,目前全路正在大力建設集裝箱中心站,積極將集裝箱運輸作為鐵路運輸發(fā)展的主要方向之一。選用基于灰色馬爾可夫鏈組合模型進行鐵路集裝箱運量預測,為科學分析和預測集裝箱運量提供理論支持。
1.1 灰色馬爾可夫鏈組合模型概述
隨著自然科學的發(fā)展和計算機應用技術的不斷推廣,多種交通量預測模型產生,預測的準確性也隨之提高[1]。常用的預測方法可以分為定性預測和定量預測。定性預測方法主要以專家為索取信息的對象,組織相關專家,通過對過去和現(xiàn)在發(fā)生的問題進行綜合分析,從中找出規(guī)律,對未來作出判斷[2]。定性預測方法包括集體意見法、頭腦風暴法、德爾菲法、情景分析法等。定量預測方法是用定量分析研究運量的發(fā)展趨勢,它以歷史統(tǒng)計資料和有關信息為依據(jù),運用各種數(shù)學方法預測未來貨運市場需求情況[3]。定量預測方法包括時間序列預測法、因果分析預測法、灰色預測法、組合預測法等。
鐵路集裝箱運量預測受多種因素影響,主要包括社會經濟因素、箱源基礎、設施條件和服務水平、政策及其他因素。在鐵路集裝箱運量的影響因素中,有的因素可以量化,有的因素只能做定性分析。隨著經濟社會的快速發(fā)展,以及鐵路推進現(xiàn)代物流轉型發(fā)展、實施鐵路供給側改革、轉變服務觀念,在這種背景下,影響鐵路貨物運量的因素更加復雜,而這些因素部分是確定的,部分是不確定的?;疑到y(tǒng)預測避開系統(tǒng)復雜的相互關系,著眼于系統(tǒng)本身的灰色信
息,通過處理原始數(shù)據(jù)和建立灰色模型,掌握系統(tǒng)發(fā)展規(guī)律,對系統(tǒng)的未來狀態(tài)做出科學的定量預測,達到使灰色系統(tǒng)白化的目的[4]。GM (1,1) 模型是灰色系統(tǒng)預測理論的基本模型和具體應用。GM 代表灰色模型(Grey Model)[5],第 1 個數(shù)字“l(fā)”表示微分方程的階數(shù),第 2 個數(shù)字“1”表示只含 1 個變量?;疑到y(tǒng)預測理論的基本思路是將已知的數(shù)據(jù)序列按照某種規(guī)則構成動態(tài)或非動態(tài)的白色模塊,再按照某種變換或解法來求解未來的灰色模型,在灰色模塊中再按照某種準則,逐步提高白度,從而確定系統(tǒng)在未來發(fā)展變化的趨勢,為事物的規(guī)劃決策、系統(tǒng)的控制與狀態(tài)的評估提供依據(jù)[6]?;疑到y(tǒng)預測具有預測精準度高、無需大量有規(guī)律樣本、計算工作量小等優(yōu)點[7],但其預測結果精確度受原始數(shù)據(jù)變化幅度的影響較大,可以通過結合馬爾可夫鏈修正的灰色預測來體現(xiàn)GM (1,1) 預測數(shù)據(jù)時間序列宏觀上的變化規(guī)律,以及馬爾可夫鏈預測數(shù)據(jù)時間序列的微觀波動規(guī)律,使預測結果得到較大改善,其精度明顯優(yōu)于傳統(tǒng)的灰色馬爾可夫模型或單一的灰色預測模型[8]。因此,選用基于灰色馬爾可夫鏈組合模型對集裝箱運量進行預測,并且進行算例驗證。
1.2 構建灰色馬爾可夫鏈組合模型
(1)原始數(shù)列生成累加數(shù)列[9]。記原始數(shù)列為x(0)= (x(0)(1),x(0)(2),…,x(0)(k),…,x(0)(n)),x(0)表示 x(0)(k) 中第 k 個原始數(shù)據(jù),n 表示原始數(shù)據(jù)的個數(shù)。對原始數(shù)列經過一次累加后生成累加數(shù)列為x(1)= (x(1)(1),x(1)(2),…,x(1)(k),…,x(1)(n)),其中,。
(2)構造均值數(shù)列。對 x(1)作緊臨均值生成均值數(shù)列 z(1)= (z(1)(1),z(1)(2),z(1)(3),…,z(1)(k),…,z(1)(n)),其中,z(1)(k) = 0.5x(1)(k) + 0.5x(1)(k-1),k = 2,3,…,n。
(3)求解灰色參數(shù)。將緊臨均值數(shù)據(jù)代入灰色微分方程 x(0)(k) + ax(1)(k) = b 中,k = 2,3,…,n,a,b為灰色參數(shù),記
用最小二乘法求解灰色參數(shù),得
(4)建立灰色模型預測公式。
(5)殘差檢驗。模型相對殘差可以表示為
如果 |δ (k+ 1)| < 16%,則模型達到一般要求;如果 |δ (k+ 1)| < 8%,則模型達到較高要求。
(6)對預測結果進行馬爾可夫鏈修正。
①預測效果狀態(tài)劃分。根據(jù)馬爾可夫鏈隨機過程分析方法的應用經驗及相對殘差幅度分布情況,對預測效果進行狀態(tài)劃分,如表 1 所示。
表 1 預測效果狀態(tài)劃分標準
②生成預測結果轉移概率矩陣。記 aij表示狀態(tài) i轉移到狀態(tài) j 的次數(shù)之和,ai表示狀態(tài) i 轉移到其他狀態(tài)的次數(shù)之和,從狀態(tài)劃分情況可得出狀態(tài)轉移情況,如表 2 所示。
表 2 預測結果轉移
根據(jù)表 2,可得預測結果轉移概率矩陣為
式中:P 為 aij占 ai的比率,即狀態(tài) i 轉移到狀態(tài) j 的次數(shù)之和占狀態(tài) i 轉移到其他狀態(tài)的次數(shù)之和的比率。pij的計算公式為
③生成預測狀態(tài)向量。根據(jù)馬爾可夫鏈預測模型,得到預測狀態(tài)向量
式中:k=0,1,…,n,A(0)為初始時刻的狀態(tài)概率向量,,為初始時刻的狀態(tài) i 對應的取值,當初始時刻的狀態(tài)為狀態(tài)i 時,,否則為 k + 1時刻的預測狀態(tài)概率向量,為 k + 1 時刻的狀態(tài) i 對應的取值。
④生成修正后的預測值。對公式 ⑵ 得出的預測值進行修正,修正表達式為
式中:Qi表示狀態(tài)i 的相對殘差幅度中間值。
2.1 蘭州鐵路局集裝箱運輸現(xiàn)狀及特點
蘭州鐵路局位于全國鐵路網西北部,跨越甘肅省和寧夏回族自治區(qū),處于亞歐大陸橋在中國境內的重要區(qū)段,東連西安,西通烏魯木齊,南接西寧,北望銀川,是西北交通運輸和經濟建設的大動脈。蘭州鐵路局貨源結構以資源輸出為主,產成品輸出為輔。自2005 年起,蘭州鐵路局貨物發(fā)送量穩(wěn)步攀升,2012年達到歷史高點,2013 年同比開始下降,但鐵路集裝箱運量除2013 年有所回落外,整體呈線性上漲的趨勢,自2005 年的 113 萬 t 至 2015 年的 445 萬 t,累計上漲 332 萬 t。2015 年以來,蘭州鐵路局實施了現(xiàn)代物流轉型發(fā)展戰(zhàn)略,貨物運輸水平和服務水平得到顯著提升,其中集裝箱運輸對全局鐵路貨運量增長貢獻突出。蘭州鐵路綜合貨場項目是全路 18 個區(qū)域性集裝箱中心站之一,該項目預計 2017 年初建成,投入運營后,將成為西北地區(qū)集裝箱裝卸、運輸?shù)闹行暮蜆屑~,通過鐵路實現(xiàn)與全國其他中心站的聯(lián)通,蘭州鐵路局集裝箱運量將呈現(xiàn)快速增長態(tài)勢。蘭州鐵路局2005—2015 年集裝箱發(fā)送量趨勢如圖 1 所示。
圖 1 蘭州鐵路局 2005—2015 年集裝箱發(fā)送量趨勢圖
蘭州鐵路局集裝箱運輸有以下特點:一是集裝箱發(fā)送量增幅明顯,占比逐年增加。由 2006 年的 2.3%增加到 2015 年的 5.9%,10 年時間集裝箱占比翻了近2 番,可見集裝箱運輸越來越得到市場的認可,也是未來蘭州鐵路局貨運增長的一個亮點。二是主要集裝箱辦理站發(fā)送量增幅明顯,發(fā)貨人較為集中。三是集裝箱裝運的品名眾多,主要品類占比較大,如 2015 年集裝箱鋼鐵運量占集裝箱發(fā)送量的近 80%。四是集裝箱運輸多為長距離運輸,運距集中在 1 000~3 500 km。五是鐵路集裝箱運費較低。
2.2 基于灰色馬爾可夫鏈組合模型的集裝箱運量預測
(1)以蘭州鐵路局 2005—2015 年集裝箱運量為參考數(shù)據(jù),生成原始數(shù)列 x(0),并且對原始數(shù)列進行一次累加生成累加數(shù)列 x(1),如表 3 所示。
表 3 2005—2015 年蘭州鐵路局鐵路集裝箱發(fā)送量原始數(shù)列及累加數(shù)列
(2)對 x(1)作緊臨均值生成新的數(shù)列 z(1),具體數(shù)值如表 4 所示。
表 4 x(1)的緊臨均值生成數(shù)列z(1)
(3)將表 4 中 z(1)(k) 數(shù)據(jù)代入公式⑴,得
(4)根據(jù)公式 ⑵,蘭州鐵路局集裝箱貨物發(fā)送量灰色預測模型可表示為
(5)檢驗殘差。根據(jù)集裝箱發(fā)送量的灰色模型計算 2006—2015 年集裝箱運量預測值,預測結果如表 5 所示。
表 5 蘭州鐵路局集裝箱發(fā)送量的 GM (1,1) 預測結果
由表 5 中預測數(shù)據(jù),可得蘭州鐵路局集裝箱發(fā)送量預測結果平均殘差為 7.7%,在可允許和接受的范圍內,可以用來預測,但準確度仍然有待進一步提高。
(6)對預測結果進行馬爾可夫鏈修正。根據(jù)馬爾可夫鏈隨機過程分析方法的應用經驗及相對殘差幅度分布情況,對其進行狀態(tài)劃分。根據(jù)表 5 中相對殘差幅度值及表 1 的預測效果狀態(tài)劃分標準可得,狀態(tài)2 的弱上升年份為 2010—2012 年、2015 年,狀態(tài) 3 的正常年份為 2007 年,狀態(tài) 4 的弱下降年份為 2006年、2008 年、2009 年、2014 年,狀態(tài) 5 的強下降年份為 2013 年。蘭州鐵路局 2006—2015 年集裝箱運量預測結果狀態(tài)轉移情況如表 6 所示。
表 6 集裝箱運量預測結果轉移
根據(jù)表 6 可得蘭州鐵路局集裝箱運量轉移概率矩陣 P。 2015 年相對殘差為 -4.56%,其預測結果狀態(tài)劃分至狀態(tài) 2,其狀態(tài)概率向量可表示為:A(0)= (1,
0,0,0)。根據(jù)公式 ⑹,得到蘭州鐵路局 2016—
2020 年集裝箱運量預測狀態(tài)向量,如表 7 所示。
表 7 蘭州鐵路局 2016—2020 年集裝箱運量預測狀態(tài)向量
根據(jù)公式 ⑺,經馬爾可夫鏈修正后得到的 2016—2020 年的預測結果如表 8 所示。
通過馬爾可夫鏈對集裝箱運量預測值進行修正后,平均相對殘差為 3.05%,較之前灰色 GM(1,
1) 預測模型的平均相對殘差 7.7% 相比,精度明顯提高,預測準確度得到較大改善,體現(xiàn)了灰色馬爾可夫鏈組合模型的優(yōu)點。在使用灰色 GM(1,1) 進行集裝箱運量預測的過程中,發(fā)現(xiàn)突變性因素無法在灰色預測模型中反映,如鐵路集裝箱運量在 2010—2012 年上升過快,而在 2013 年有所下降等。能夠明顯看到2008 年、2011 年、2012年、2013 年和 2014 年的擬合數(shù)據(jù)與實際數(shù)據(jù)存在較大的差距,相對殘差分別為9.99%、-9.55%、-9.44%、11.83% 和 9.7%,遠遠高于平均相對殘差 7.7%。一般的模型無法處理突變因素,甚至受突異點的影響,將誤差帶到預測中,而馬爾可夫狀態(tài)轉移將歷史的波動信息考慮進來,使得預測結果準確度提高[10]。
表 8 基于灰色馬爾可夫鏈組合模型的蘭州鐵路局集裝箱運量預測結果
集裝化運輸由于中途運輸換裝次數(shù)少、機械化程度和裝卸效率高、運費比較優(yōu)惠,能夠提升運輸貨品的安全性,降低貨品的損毀率,在市場競爭中具有優(yōu)勢。大力發(fā)展集裝箱運輸,突出集裝箱運輸優(yōu)勢,擴大鐵路集裝箱運輸占比,是實現(xiàn)鐵路向現(xiàn)代物流轉型的手段之一。利用灰色馬爾可夫鏈組合模型進行預測,能夠發(fā)揮灰色系統(tǒng)預測模型的優(yōu)點,并在考慮原始數(shù)據(jù)影響的基礎上,考慮歷年狀態(tài)轉移情況,避免了灰色系統(tǒng)預測模型受原始數(shù)據(jù)變化幅度影響大的缺點?;诨疑R爾可夫鏈組合模型的集裝箱運量預測,為鐵路科學確定考核指標、合理安排車輛運用計劃提供重要的依據(jù),其分析結果對于鐵路部門制定物流營銷政策、科學組織貨物運輸具有重要作用。
[1] 宋光平. 鐵路貨運量預測方法研究[D]. 北京:北京交通大學,2007.
[2] 羅安全. 鐵路現(xiàn)代物流中心的構建與運作模式研究[D].長沙:中南大學,2009.
[3] 程禮芬. 區(qū)域貨運量預測方法的比較[J]. 重慶交通學院學報,2005,24(4):114-116. CHENG Li-fen. Compare of Regional Freight Volume Forecast Methods[J]. Journal of Chongqing Jiaotong University,2005,24(4):114-116.
[4] 郭 晶,姜秀山. 鐵路集裝箱運量發(fā)展預測[J]. 物流科技,2006,29(1):31-34. GUO Jing,JIANG Xiu-shan. Application of Grey System to the Forecast of the Railroad Container Volume[J]. Logistics Technology,2006,29(1):31-34.
[5] 吳 璇,王 烈. 組合預測模型在鐵路集裝箱運量預測中的應用[J]. 鐵道運輸與經濟,2010,32(10):90-94. WU Xuan,WANG Lie. Application of Combined Forecasting Model in Railway Container Prediction[J]. Railway Transport and Economy,2010,32(10):90-94.
[6] 張 誠,周湘峰. 基于灰色預測-馬爾可夫鏈-定性分析的鐵路貨運量預測[J]. 鐵道學報,2007,29(5):15-21. ZHANG Cheng,ZHOU Xiang-feng. Railway Freight Volume Forecasting based on Grey Prediction-Markov-Qualitative Analysis[J]. Railway Transaction,2007,29(5):15-21.
[7] 張 玥,帥 斌. 基于改進灰色模型的東北地區(qū)鐵路貨運量預測[J]. 鐵道科學與工程學報,2012,9(5):125-128. ZHANG Yue,SHUAI Bin. Railway Freight Volume Forecast in Northeast Region based on the Improved Gray
Model[J]. Journal of Railway Science and Engineering,2012,9(5):125-128.
[8] 謝建文,張元標,王志偉. 基于無偏灰色模糊馬爾可夫鏈法的鐵路貨運量預測研究[J]. 鐵道學報,2009,31(1):1-7. XIE Jian-wen,ZHANG Yuan-biao,WANG Zhi-wei. Railway Freight Volume Forecasting based on Unbiased Grey -Fuzzy-Markov Chain Method[J]. Journal of the China Railway Society,2009,31(1):1-7.
[9] 于 洋,楊學斌,杜 文. 基于灰色-馬爾可夫鏈改進方法的預測方法[J]. 統(tǒng)計與決策,2013(13):53-55.
[10] 李金龍,張紅亮. 我國鐵路集裝箱運輸發(fā)展對策的思考[J]. 鐵道貨運,2013,31(1):47-50. LI Jin-long,ZHANG Hong-liang. Thoughts on the Development of Railway Container Transportation in China[J]. Railway Freight Transport,2013,31(1):47-50.
(責任編輯 王 靜)
Container Transportation Volume Prediction of Lanzhou Railway Administration based on Gray Markov Combination Chain Model
QIU Fei
(Passenger and Freight Transportation Statistical Institute, Lanzhou Railway Administration, Lanzhou 730000, Gansu, China)
Through elaborating on the Gray Markov combination chain model theory, this paper analyses Present situation and characteristics of container transportation of Lanzhou Railway Administration, and applies Markov combination chain model for prediction and analysis for container transportation volume of Lanzhou Railway Administration based on current status and historical data, which provides
for railway department to formulate logistic marketing policy and scientifically organize freight transportation.
Railway; Container; Freight Volume Prediction; Gray Prediction, Markov Model
1004-2024(2016)09-0024-06
U294.3
B
10.16669/j.cnki.issn.1004-2024.2016.09.06
2016-07-25
邱菲(1982—),女,天津人,大學本科。