中國和巴西港口擁塞指數多重分形研究

林國龍，馬理博，張辰彥，何紅弟

(上海海事大學 物流研究中心，上海 201306)

中國和巴西港口擁塞指數多重分形研究

林國龍，馬理博，張辰彥，何紅弟

(上海海事大學 物流研究中心，上海 201306)

運用多重分形消除趨勢波動分析法(MF-DFA)與多重分形消除趨勢波動交叉相關分析法(MF-DCCA)對全球港口擁塞指數的中國港口ETA(Estimated Time of Arrival)海岬型船與巴西港口ETA海岬型船港口擁塞指數進行研究。研究表明：中國港口ETA海岬型船與巴西港口ETA海岬型船時間序列廣義Hurst指數及MF-DCCA廣義Hurst指數，均單調遞減且值域在(0.5,1)之間，標度指數均表現凹向橫軸的函數，分形譜密集鐘形分布；兩組時間序列均具有明顯的多重分形特征，同時具有長期記憶性與多標度性，且二者多重分形交叉相關，表現長期相互影響。

交通運輸工程；港口擁塞指數；多重分形；消除趨勢波動分析；交叉相關分析

0 引 言

由全球港口〔G-ports (UK) Ltd〕發布的全球港口擁塞指數報告(The Global Ports Congestion Index，GPCI)是衡量全球港口擁堵狀況的有效參考。GPCI以周報的形式發布，詳細及時地統計反映全球滯留主要港口的煤炭、礦石等干散貨船數量，分析由于船舶壓港、港口擁塞對干散貨市場供求關系以及散貨船租金所帶來的影響，及時追蹤全球主要港口煤炭、礦石等干散貨船滯期情況，覆蓋了全球不同國家80個煤炭、礦石碼頭的擁塞情況，包括澳大利亞、巴西、中國、印度以及南非等國。GPCI能夠反映全球港口擁堵現象及程度。報告中港口擁塞指數是非平穩的時間序列。通過研究港口擁塞指數時間序列的變化情況，探索港口擁塞指數間的內在相關性，發現港口擁塞指數的發展特征，從而能夠預測港口擁堵現象的發展方向。

港口擁堵問題備受全球關注，源于2002年美西碼頭工人罷工導致洛杉磯、長灘碼頭擁堵。近年來進出口貿易往來頻繁對國際海運航運行業刺激加劇，港口擁堵問題也隨之加劇且成為全球性問題。尤其是集裝箱貿易運輸，運價運力皆不斷攀升。張榮忠[1]在全球性港口擁堵現象透析文中就港口擁堵現狀的原因及解決辦法給出過介紹。尹銘鈺[2-3]依據全球港口擁塞指數報告，詳細分析全球干散貨港口擁堵現狀，定性指出滯港船舶情況呈現區域性不均衡、巴西、澳大利亞港口擁堵嚴重、中國港口擁堵集中于鐵礦石運輸港的港口擁堵問題；隨后針對巴西港口擁堵情況分析得出，2014年港口擁堵情況略有好轉，但整體而言，巴西港口的擁堵現狀明顯要嚴峻于世界平均水平。徐劍華[4]指出美西港口擁堵的深層原因包括船舶大型化、班輪公司超級聯盟和承運人為節約成本出售自有底盤車隊。 P.T.LEACH[5-6]就美西港口擁堵問題仍對航運市場造成的影響，提出重新定義航運業務旺季。

分形(fractal)是20世紀70年代Mandelbort表征復雜圖形和過程時首先引入于科學領域，表示不規則的物體。C.K.PENG等[7]在研究DNA序列時首次提出消除趨勢波動分析方法(detrended fluctuation analysis，DFA)，該方法能有效計算一個非平穩時間序列地長程相關性，即非平穩時間序列的內在自相似性。J.W.KANTELHARDT等[8]提出檢驗一個非平穩時間序列是否具有多重分形特征的方法，多重分形消除趨勢波動分析法(multifractal detrended fluctuation analysis，MF-DFA)是在DFA基礎上更加準確有效的判斷一個時間序列長程相關性的定量方法。B.PODOBNIK等[9]首次提出定量分析兩列非平穩時間序列相關性的科學分析方法消除趨勢波動交叉相關分析法(Detrended Cross-Correlation Analysis，DCCA)。L.Y.HE等[10]在檢驗中美兩國農產品期貨市場的量價關系時，首次采用多分形消除趨勢波動交叉相關分析法(Multifractal Detrended Cross-Correlation Analysis，MF-DCCA)，并證實其多重分形特征及交叉相關性。近年來，DFA，MF-DFA，DCCA和MF-DCCA被應用于多個領域時間序列研究，包括生命自然科學[11]、證券黃金期貨市場[12-14]、原油集裝箱市場[15-16]、水文地理地質[17-18]、設備故障診斷[19]、能源[20]、交通[21]等多方領域。如此廣泛地應用足以說明消除趨勢波動分析法與消除趨勢波動交叉相關分析法在研究復雜非平穩時間序列內在自相關和互相關的價值與作用。

采用MF-DFA與MF-DCCA對全球港口擁塞指數的中國港口ETA(Estimated Time of Arrival)海岬型船(100000mts+)與巴西港口ETA海岬型船港口擁塞指數進行研究，MF-DFA結果發現中國港口ETA海岬型船與巴西港口ETA海岬型船港口擁塞指數時間序列具有明顯的多重分形特征，廣義Hurst指數均大于0.5小于1，尺度函數均為凹向橫軸的函數，多重分形譜呈密集鐘形分布。MF-DCCA結果發現中國港口ETA海岬型船與巴西港口ETA海岬型船港口擁塞指數時間序列具有明顯的多重分形交叉相關特征, 廣義Hurst指數大于0.5小于1，尺度函數為凹向橫軸的函數，多重分形譜呈密集鐘形分布。

1 MF-DFA與MF-DCCA基本原理

多重分形消除趨勢波動交叉相關分析法(MF-DCCA)是計算兩個非平穩時間序列內在交叉相關性的一種定量方法，理論基礎是多重分形消除趨勢波動分析法(MF-DFA)，MF-DFA是在DFA[22]的基礎上檢驗一個非平穩時間序列多重分型特性的方法。對長度為N的時間序列{xi},i=1,2,…,N和等長度時間序列{yi},i=1,2,…,N，MF-DCCA的具體步驟：

1)計算兩個時間序列的累積離差：

(1)

(2)

2)將累積離差X(t),Y(t)分別等分成Ns個互不相關時間長度為s的等時距區間。由于時間序列總長N不一定是長度s的整數倍，為保證時間序列信息完整，將時間序列反向再計算一次累積離差。因此求得兩組2Ns個長度相同子區間,記每個子區間為v=(1,2,…,2Ns)。

3)對兩組時間序列每個子區間的時間數據采用最小二乘法分別進行多項式回歸擬合，得擬合多項式xv(t)和yv(t)，求局部協方差函數。

當v=1,2,…,2Ns時：

{Y[(v-1)s+i]-yv(i)}

(3)

當v=Ns+1,…,2Ns時：

(4)

4)把所有子區間的局部協方差取均值，計算所有區間的平均值，得到時間序列的q階波動函數

(5)

當q=0時，時間序列的波動函數為

(6)

當q=2時，對一列時間序列的MF-DFA退化成標準DFA。

5)若時間序列長程冪律交叉相關，Fq(s)與s應滿足冪律關系式：

Fq(s)～sHxy(q)

(7)

將式(7)兩邊同時取對數得

ln[Fq(s)]～lns

(8)

將時間序列在對數坐標中的散點圖用最小二乘法進行擬合，擬合函數斜率即為q階廣義Hurst指數Hxy(q)。

6)經過MF-DCCA的計算得到對應不同q下的廣義Hurst指數Hxy(q)，我們也稱其為標度指數，這個指數是衡量時間序列是否長程冪律交叉相關的參數，標度指數的變化特征描述為：①Hxy(q)是常數時，時間序列具有單分形特征；②是q的單減函數，時間序列具有多重分型特征；③Hxy(q)值域為(0,1)時，可判斷時間序列是否長程冪律交叉相關；④Hxy(q)值域為(0,0.5)時，時間序列具有反長程冪律交叉相關性；⑤Hxy(q)為常數0.5時，時間序列隨機游走；⑥Hxy(q)值域為(0.5,1)時，時間序列具有正長程冪律交叉相關性。

7)廣義Hurst指數Hxy(q)與質量指數τ(q)的關系式為

τ(q)=qHxy(q)-1

(9)

若實驗時間序列表現多重分形特征，則質量指數為凹向橫軸的函數，質量指數可進一步檢驗時間序列是否具有長程冪律交叉相關的特性。

8)經過Legendere變換，可得到多標度時間序列的多重分形譜f(α)：

(10)

f(α)=q[α-Hxy(q)]+1

(11)

式中：α為時間序列區間上的奇異指數，α值越小表示時間區間奇異性越大；f(α)為多重分形譜，多重分形譜能說明實驗時間序列發展趨勢，分形譜越寬則實驗時間序列的波動越大，分形譜曲線極大值為1，若左端高于右端說明實驗時間序列出現高位置概率更大，反之，若右端高于左端則說明實驗時間序列出現低位置概率更大，而且圖像左右高度距離差值越大，時間序列出現高低位置概率相差也就越大。

2 港口擁塞指數多重分形研究

筆者用全球港口擁塞指數報告中時間序列數據進行MF-DFA，MF-DCCA研究，發現中國港口ETA海岬型船與巴西港口ETA海岬型船港口擁塞指數時間序列具有較為明顯的多重分形特征，且二者明顯多重分形交叉相關。

2.1 數據來源

選擇地理位置上無明顯相關關系的中國和巴西港口，選用由全球港口公司發布的全球港口擁塞指數報告的中國港口ETA海岬型船和巴西港口ETA海岬型船港口擁塞指數，2008年第2周至2015年第17周共380個周數據進行研究，其原始數據分布如圖1。

圖1 中國和巴西港口ETA海岬型船港口擁塞指數波動Fig.1 Congestion index fluctuations of Chinese and Brazil ETA Capesize port

2.2 研究分析

運用MF-DCCA對港口擁塞指數時間序列進行實驗，將實驗數據導入已編碼的MF-DCCA MATLAB R2013a程序并運行，提取運行數據用OriginPro8.0作圖，見圖2。

圖2 Fq(s)與s的對數關系圖Fig.2 The logarithmic relation diagram of Fq(s) and s

由圖2可見，對應式(1)～式(6)，表示q變化時波動函數Fq(s)與s呈雙對數關系，圖中從上到下3條線表示q=5,0，-5時時間序列的擬合狀況，隨著q的增大擬合函數圖像單調遞增且斜率減小，因為擬合函數的斜率即是為廣義Hurst指數h(q)，因此實驗結果能夠初步說明廣義Hurst指數h(q)隨q的增大而減小。

圖3為對應式(7)和式(8)給出MF-DFA與MF-DCCA廣義Hurst指數。

圖3 港口擁塞指數時間序列廣義Hurst指數Fig.3 Generalized Hurst index of port congestion index in time series

由圖3可見，q增大h(q)單調遞減，h(q)值域在(0.5,1)內，充分說明巴西港口ETA海岬型船與中國港口ETA海岬型船港口擁塞指數時間序列均具有多重分形特征，均表現正的長程冪律相關。巴西港口ETA海岬型船港口擁塞指數時間序列廣義Hurst指數全部位于中國港口ETA海岬型船港口擁塞指數時間序列廣義Hurst指數下方，且值域更大，說明巴西港口ETA海岬型船港口擁塞指數序多重分形特征更為明顯，其函數圖像的斜率下降更快。MF-DCCA分析結果表明，廣義Hurst指數Hxy(q)位于MF-DFA Hurst指數h(q)結果之間且單調遞減。值域在(0.5,1)內，說明巴西港口ETA海岬型船與中國港口ETA海岬型船港口擁塞指數時間序列之間具有多重分形交叉相關特性，且呈正的冪律相關關系。

為了進一步驗證時間序列多重分形特征，同時給出港口擁塞指數時間序列尺度函數如圖4。

圖4 港口擁塞指數時間序列尺度函數Fig.4 Port congestion index time series scale function chart

由圖4可見，巴西港口ETA海岬型船與中國港口ETA海岬型船港口擁塞指數尺度函數τ(q)均為單調遞增的凹向橫軸的函數，證明實驗時間序列是多重分形的。MF-DCCA對比結果顯示，三者為基本重合且凹向橫軸的單調遞增曲線，更進一步說明二者之間的多重分形交叉相關特性。

此外，根據式(10)和式(11)計算給出港口擁塞指數時間序列分形譜，如圖5。

圖5 港口擁塞指數時間序列分形譜Fig.5 Fractal spectrum of port congestion index in time series

由圖5可知，巴西港口ETA海岬型船港口擁塞指數時間序列比中國港口ETA海岬型船港口擁塞指數時間序列分形譜更寬，且左端明顯高于右端，幅度更大。說明巴西港口ETA海岬型船港口擁塞指數時間序列波動比中國港口ETA海岬型船港口擁塞指數時間序列波動較大，兩者港口擁塞指數均表現為出現高位置時概率大，而巴西港口此特征更為明顯。MF-DCCA分析結果表示，兩組實驗時間序列多重分形交叉相關，因此港口擁塞指數時間序列受自身歷史數據影響，同時受另外一組港口擁塞指數實驗時間序列歷史數據影響。

3 結 語

通過對中國港口ETA海岬型船與巴西港口ETA海岬型船港口擁塞指數進行多重分形分析，結果發現中國港口和巴西港口擁塞指數的廣義Hurst指數均單調遞減且值域在(0.5,1)之間，并且巴西港口的擁塞指數的分形譜明顯大于中國港口的多重分形譜。這說明兩個國家的港口擁塞指數都具有多重分形的特征和長程相關的特征，并且巴西港口的多重分析特征相對更明顯。此外，還對兩個國家港口擁塞指數的互相關行為進行研究，結果發現兩者之間具有明顯的長程互相關特征。這些結果對研究不同港口擁塞指數的相互行為具有很好的參考價值，從而為研究港口擁堵問題提供依據。

[1] 張榮忠. 全球性港口擁堵現象透析[J]. 交通企業管理,2005,20(5):31-32. ZHANG Rongzhong. Analysis of the global port congestion phenomenon[J].TransportationEnterpriseManagement,2005,20(5):31-32.

[2] 尹銘鈺. 全球干散貨港口擁堵現狀分析[J]. 中國港口,2013(11):4-7. YIN Mingyu. Analysis of the current situation of global dry bulk port congestion[J].ChinaPorts,2013(11):4-7.

[3] 尹銘鈺. 巴西港口在探索變革中成長[J]. 中國港口,2014(8):61-62. YIN Mingyu. The development of Brazil port in the exploration of the changes[J].ChinaPorts,2014(8):61-62.

[4] 徐劍華. 美西港口擁堵迷局[J]. 中國船檢,2015(3):18-22. XU Jianhua. West Port congestion puzzles[J].ChineseShipSurvey,2015(3):18-22.

[5] LEACH P T. 美西港口擁堵余波未平[J]. 中國遠洋航務,2015(4):60-61. LEACH P T. West Port congestion are still repercussions[J].MaritimeChina,2015 (4):60-61.

[6] LEACH P T. 重新定義航運旺季[J]. 中國遠洋航務,2015(8):48-49. LEACH P T. Redefining shipping season[J].MaritimeChina,2015(8):48-49.

[7] PENG C K, BULDYREV S V, HAVLIN S,et al. Mosaic organization of DNA nucleotides[J].PhysicalReviewE,1994,49(2):1685-1689.

[8] KANTELHARDT J W, ZSCHIEGNER S A, KOSCIELNY B E, et al. Multifractal detrended fluctuation analysis of nonstationary time series[J].PhysicaAStatisticalMechanics&ItsApplications,2002,316(1/2/3/4):87-114.

[9] PODOBNIK B, STANLEY H E. Detrended cross-correlation analysis: a new method for analyzing two nonstationary time series[J].PhysicalReviewLetters,2007,100(8):38-71.

[10] HE L Y, CHEN S P. Nonlinear bivariate dependency of price-volume relationships in agricultural commodity futures markets: A perspective from Multifractal Detrended Cross-Correlation Analysis[J].PhysicaAStatisticalMechanics&ItsApplications,2011,390(2):297-308.

[11] WANG Jun, ZHAO Daqing. Detrended cross-correlation analysis of electroencephalogram[J].ChinesePhysicsB,2012,21(2):577-580.

[12] 曾志堅, 張倩倩. 基于MF-DCCA方法的證券市場間交叉相關性研究[J]. 財經理論與實踐,2013,34(6):45-49. ZENG Zhijian, ZHANG Qianqian. Study on the cross correlation of securities market based on MF-DCCA method[J].TheTheoryandPracticeofFinanceandEconomics,2013,34(6):45-49.

[13] 曹建軍,顧榮寶. 基于MF-DFA方法的黃金市場多重分形分析[J]. 黃金,2011,32(3):5-7. CAO Jianjun, GU Rongbao. Multi fractal analysis of gold market based on MF-DFA method[J].Gold,2011,32 (3):5-7.

[14] 郭麗娜,王宏勇. 大宗商品期貨市場收益率的多重分形分析 [J]. 南京財經大學學報,2013,20(1):43-50. GUO Li’na, WANG Hongyong. Return rate of commodity futures market multifractal analysis[J].JournalofNanjingUniversityofFinanceandEconomics,2013,20(1):43-50.

[15] 張德園,王雪飛. 基于多重分形理論的原油市場耦合關系研究 [J]. 成都理工大學學報(社會科學版),2014,22(5):58-64. ZHANG Deyuan, WANG Xuefei. Research on the coupling relationship of crude oil markets based on multifractal theory[J].JournalofChengduUniversityofTechnology(SocialSciences) ,2014,22(5):58-64.

[16] 劉俊超,陳秀乾. 上海集裝箱運價衍生品市場多重分形特征研究[J]. 航海,2014(1):47-51. LIU Junchao, CHEN Xiuqian. Research on multi fractal characteristics of container transport market in Shanghai[J].Navigation,2014(1):47-51.

[17] SEO S K, KIM K, CHANG K H, et al. Determination of the dynamical behavior of rainfalls by using a multifractal detrended fluctuation analysis[J].JournaloftheKoreanPhysicalSociety,2012,61(4):658-661.

[18] LIM G, KIM K, PARK J K, et al. Dynamical analyses of the time series for the temperature and the humidity[J].JournaloftheKoreanPhysicalSociety,2013,62(1):193-196.

[19] 胡江,蘇懷智,馬福恒,等. MF-DFA在大壩安全監測序列分析和整體性態識別中的應用[J]. 水利水電科技進展,2014,34(3):50-55. HU Jiang, SU Huaizhi, MA Fuheng, et al. The application of MF-DFA in time series analysis and global state recognition of dam safety monitoring[J].AdvancesinScienceandTechnologyofWaterResources,2014,34(3):50-55.

[20] WANG F, LIAO G P, ZHOU X Y, et al. Multifractal detrended cross-correlation analysis for power markets[J].NonlinearDynamics,2013,72(1/2):353-363.

[21] XU N, SHANG P, KAMAE S. Modeling traffic flow correlation using DFA and DCCA[J].NonlinearDynamics,2010,61(1/2):207-216.

[22] 魏宇,黃登仕. 中國股票市場波動持久性特征的DFA分析[J]. 中國管理科學,2004,12(4):13-20. WEI Yu, HUANG Dengshi. A DFA study on the persistence of fluctuations in china stock market[J].ChineseJournalofManagementScience,2004,12(4):13-20.

Multifractal Research on Congestion Index of Chinese and Brazil Port

LIN Guolong, MA Libo, ZHANG Chenyan, HE Hongdi

(Logistics Research Center, Shanghai Maritime University, Shanghai 201306, P.R.China)

The multifractal detrended fluctuation analysis method (MF-DFA) and multifractal detrended fluctuation cross correlation analysis method (MF-DCCA) were used to study the congestion index of ETA (Etimated Ttime of Arrival) capesize in Chinese ports and in Brazil ports, which were components of GPCI (the global port congestion index). It is found that generalized Hurst index of ETA capesize in Chinese ports and in Brazil ports in time series and MF-DCCA generalized Hurst index are both monotonically decreasing and the range of their values is between 0.5 and 1; scaling index shows a concave function and fractal spectrum shows a intensive bell shaped distribution. Research results show that the two groups of experimental time series both have obvious multifractal characteristics and also have a long memory and multi-marked degree; moreover, the two multifractals are cross-correlated, and their performances are mutually influenced in a long-term.

traffic and transportation engineering; port congestion index; multifractal; detrended fluctuation analysis; cross correlation analysis

2015-08-21；

2016-03-15

林國龍(1951—)，男，浙江象山人，教授，博士生導師，主要從事國際航運、自由貿易區、綜合物流環境與供應鏈一體化管理方面的研究。E-mail：LingLzm@163.com。

10.3969/j.issn.1674-0696.2016.06.30

U691+.32

A

1674-0696(2016)06-148-05