林國龍， 王鋮, 丁一

(上海海事大學(xué) 科學(xué)研究院， 上海 201306)

0 引 言

干散貨航運市場是國際貿(mào)易的派生市場，具有明顯的周期波動特性.一方面由于世界經(jīng)濟的變化導(dǎo)致市場對煤炭、谷物等干散貨需求的變化，另一方面由于超靈便型船舶市場運力供需的變化,使作為干散貨航運市場組成部分的超靈便型船舶市場也產(chǎn)生周期性波動.波羅的海干散貨運價指數(shù)(Baltic Dry Index, BDI)是反映世界經(jīng)濟走勢的晴雨表，而超靈便型船運價指數(shù)(Baltic Supramax Index, BSI)是BDI的重要組成部分.超靈便型船既能適應(yīng)船舶大型化的趨勢又具有自身的靈便性.它可以根據(jù)航運市場的走勢調(diào)節(jié)貨源，既可以運輸巴拿馬型船的煤炭谷物又可以運輸小靈便型船的小宗散貨等.干散貨航運市場方面，學(xué)者對BDI指數(shù)的波動特性等進行過大量分析和研究.XU等[1]通過建立AR-GARCH模型并進行GMM回歸分析得出船隊規(guī)模增長與干散貨市場運價波動的關(guān)系.DURU[2]運用多變量模糊邏輯預(yù)測模型(M-FILF)和與模糊時間序列相關(guān)的理論對干散貨市場的波動性進行研究并作出短期預(yù)測.GEMAN等[3]運用不變方差彈性模型(CEV)對1988—2010年間的BDI的隨機波動特性進行分析，得出現(xiàn)貨價格模型不僅對未來運費率期權(quán)的增長、而且對谷物生產(chǎn)者和商人都是非常有價值的工具.聶金龍等[4]針對BDI指數(shù)收益序列存在長期記憶性引入ARFIMA模型來改變傳統(tǒng)的預(yù)測模型.李耀鼎等[5]運用ADF檢驗和ARCH LM檢驗得出BDI指數(shù)序列存在一階單整性和高階ARCH效應(yīng),并運用ARCH(1,1)消除殘差序列的條件異方差性.陸克從等[6]通過GARCH(1,1)模型分析出現(xiàn)今的BDI收益率序列有明顯波動集聚效應(yīng),并結(jié)合實際分析特別波動的原因,同時利用EGARCH模型刻畫不同時期外部信息對國際干散貨市場造成的影響.張頁等[7]對小靈便型船運價指數(shù)(BHI)進行研究并建立ARMA模型.在濾波分析方法方面，林國龍等[8]運用CF濾波分析及譜分析方法分析國際干散貨航運市場的周期，得出國際干散貨航運市場存在3.5 a左右的周期,BDI指數(shù)是新造船市場、二手船市場和拆船市場的領(lǐng)先指標.王磊等[9]運用譜分析方法研究海岬型船平均收入與新造船價格之間的周期波動關(guān)系.韓荊[10]運用Census X12季節(jié)調(diào)整法和HP濾波法，結(jié)合吉林省人民銀行的調(diào)研數(shù)據(jù),對現(xiàn)金投放的特點進行分析并提出參考建議.鄭衛(wèi)國等[11]運用HP濾波法分析我國經(jīng)濟周期的波動特征，總結(jié)出我國經(jīng)濟周期波動的典型化事實.徐國祥等[12]采用加窗平均周期圖譜分析和多次分辨法分析我國房地產(chǎn)市場周期波動特性，并根據(jù)波動特性提出相應(yīng)的對策建議.

BSI指數(shù)自2005年7月正式發(fā)布以來對其研究很少.本文將Census X12季節(jié)調(diào)整法、HP濾波和BP濾波方法相結(jié)合，將帶有季節(jié)因素的BSI月度值序列進行分解，分別找出其中的趨勢要素、季節(jié)要素和循環(huán)要素,以便更好地發(fā)現(xiàn)BSI序列的本質(zhì)波動規(guī)律.

1 研究方法

1.1 Census X12季節(jié)調(diào)整法

現(xiàn)實的時間序列經(jīng)濟數(shù)據(jù)會存在月度或者季度的循環(huán)變動.出現(xiàn)季度變動的時間序列數(shù)據(jù)可能會掩蓋其本質(zhì)規(guī)律，解決辦法是采用季節(jié)調(diào)整法從時間序列中去除季節(jié)變動要素，得到序列潛在的趨勢循環(huán)分量.趨勢循環(huán)分量在沒有季節(jié)要素干擾的情況下,能夠如實反映經(jīng)濟數(shù)據(jù)時間序列運動的客觀規(guī)律.本文采用Census X12[13]季節(jié)調(diào)整法進行季節(jié)調(diào)整.Census X12是在Census X11的基礎(chǔ)上發(fā)展而來的，包括4種季節(jié)調(diào)整的分解形式：加法、乘法、對數(shù)加法和偽加法模型.各模型將時間序列分為長期循環(huán)趨勢要素(TC)、季節(jié)變動要素(S)和不規(guī)則變動要素(I).

(1)加法模型的一般形式為

(1)

加法模型中季節(jié)要素和趨勢循環(huán)要素的影響采用絕對量表示，這樣分析起來比較直觀，但是因各經(jīng)濟變量之間的計量單位不同而缺乏可比性.

(2)乘法模型的一般形式為

(2)

乘法模型的主要特點在于以相對數(shù)表示季節(jié)要素，可以避免計量單位的影響，增加不同經(jīng)濟變量間的可比性，但是在直觀性方面不如加法模型.

(3)對數(shù)加法模型的一般形式為

(3)

對數(shù)加法模型是通過對乘法模型取自然對數(shù)而得到的特殊形式的加法模型.

(4)偽加法模型的一般形式為

(4)

偽加法模型主要針對一些非負的時間序列進行調(diào)整(如農(nóng)產(chǎn)品產(chǎn)量).通過偽加法模型進行季節(jié)調(diào)整,使調(diào)整后的結(jié)果更接近于序列趨勢的估計值.

季節(jié)調(diào)整法中運用比較成熟的方法是加法模型和乘法模型.當時間序列隨著時間等寬度推進或季節(jié)變動大致相等時一般采用加法模型處理；當時間序列的季節(jié)波動與長期趨勢大致相同時，即序列圖呈放射狀或者喇叭狀時,一般采用乘法模型處理.

1.2 HP濾波法

(5)

(6)

折算因子λ對于濾波效果非常重要.隨著λ的增大，估計的趨勢越光滑.根據(jù)一般經(jīng)驗值,月度數(shù)據(jù)時λ=14 400.

1.3 BP濾波法

BP濾波[13]的基本思想是把時間序列看作是互不相關(guān)的頻率分量疊加，通過研究和比較各分量的周期變化，充分揭示時間序列的頻域結(jié)構(gòu)，掌握其主要的波動特性.

頻譜濾波分析主要通過相關(guān)的頻率響應(yīng)函數(shù)對序列中存在的某些波動成分進行剔除分析.BP濾波的頻率響應(yīng)函數(shù)就是兩個低通濾波的響應(yīng)函數(shù)之差.利用BP濾波分析，將時間序列中緩慢變動和劇烈變動的成分剔除.這里，頻率為f，切斷頻率為fc，截斷點為n的低通濾波(low-pass filters)的頻率響應(yīng)函數(shù)：

(7)

HP濾波通過濾波分析分離趨勢成分和循環(huán)成分，通過分離出的趨勢曲線及循環(huán)曲線說明序列的相關(guān)特征;BP濾波為譜分析范疇，通過分析由HP濾波分離出的循環(huán)成分的頻率響應(yīng)函數(shù)的頻率值確定循環(huán)周期.通過HP濾波和BP濾波的結(jié)合分析確定序列的趨勢特征與周期循環(huán)特征.

2 實證研究

本文研究對象為BSI，數(shù)據(jù)選取自2005年7月至2012年7月波羅的海航交所發(fā)布的BSI指數(shù)的月度值，共85個數(shù)據(jù).BSI指數(shù)月度走勢見圖1.

2.1 季節(jié)波動實證分析

季節(jié)波動主要是由自然氣候及社會影響等因素造成的.船舶在不同的季節(jié)航行所面臨的自然條件迥異，導(dǎo)致航運成本在不同季節(jié)的差別.另外,季節(jié)波動還受到世界總體經(jīng)濟形勢的影響.國際航運市場是國際貿(mào)易市場的派生市場，國際航運市場運價隨著國際貿(mào)易需求量的變化而變化.整個航運市場的繁榮情況和四季氣候變化對煤炭、糧谷和一些小宗散貨的生產(chǎn)運輸產(chǎn)生影響.這些貨物的需求量和船舶供給量的波動引起B(yǎng)SI的季節(jié)性波動.由圖1可以看出序列圖存在喇叭狀放射、波動寬度不等的特性，因此對月度BSI數(shù)據(jù)采用基于乘法模型的Census X12季節(jié)調(diào)整法,并結(jié)合EViews[14]分析BSI序列的季節(jié)波動性，得出BSI序列的季節(jié)因素圖,見圖2.

圖1 BSI月度走勢

圖2 BSI序列季節(jié)因素圖

由圖2可知，4月至10月BSI季節(jié)因素百分比大于100%，表示這些月份是超靈便型船舶市場的旺季;12月至次年3月BSI季節(jié)因素百分比小于100%，說明這些月份為超靈便型船舶市場的淡季.整體而言，超靈便型船舶市場的季節(jié)波動高位點位于夏季，低位點位于冬季，表現(xiàn)出“夏高冬低”的特性.為了進一步驗證這一特性，取12月至次年1月的BSI月度數(shù)據(jù)均值代表冬季數(shù)據(jù)，6月至8月的BSI月度數(shù)據(jù)均值代表夏季數(shù)據(jù)，通過計算得出BSI冬季均值為2 142點，夏季均值為2 621點，從而通過超靈便型船舶市場的實際運價指數(shù)驗證夏季市場較冬季活躍的特征.超靈便型船舶的這一季節(jié)波動特性與干散貨航運市場“夏低冬高”的季節(jié)波動剛好相反，需要引起超靈便型船舶經(jīng)營者的注意.

2.2 HP濾波實證分析

經(jīng)過Census X12季節(jié)調(diào)整法去除序列中可能存在的長期趨勢和季節(jié)變化，得到經(jīng)過季節(jié)調(diào)整的序列Seasons Adjusted BSI(SA-BSI).對SA-BSI數(shù)據(jù)進行HP濾波分析得出循環(huán)(Cycle)序列和趨勢(Trend)序列.圖3為經(jīng)過HP濾波分析的序列圖.

圖3 季節(jié)調(diào)整序列(SA-BSI)HP濾波分析序列圖

由圖3可知，BSI的長期趨勢可以分解為2部分:第一部分為2008年經(jīng)濟危機爆發(fā)之前，運價指數(shù)表現(xiàn)為增強趨勢；第二部分為2008年經(jīng)濟危機爆發(fā)之后，作為國際貿(mào)易派生市場的航運市場受其影響，運價指數(shù)表現(xiàn)為減弱趨勢,該減弱趨勢一直延續(xù)至今.從循環(huán)序列中可以分析出BSI月度序列的波動特性：2007年7月至2009年8月波動異常，波動率超出-10%～10%的區(qū)間，其中2007年9月至2008年2月、2008年10月至2009年1月波動尤為異常，波動幅度超過20%，2007年10月波動率達到30%.BSI序列巨幅波動的根本原因在于經(jīng)濟危機的爆發(fā).經(jīng)濟危機爆發(fā)前夕，航運市場活動火爆，干散貨船舶運力劇增，2007年10月30日BSI指數(shù)達到歷史最高峰6 956點;2008年經(jīng)濟危機爆發(fā)，BSI指數(shù)開始大幅跳水下跌，2008年10月22日跌破千點.由于航運市場火爆時船舶所有人紛紛建造新船、延長船舶使用年限等原因造成運力在近年來劇增，而世界經(jīng)濟不景氣造成干散貨貿(mào)易量降低，從而使市場上的船舶運力供過于求，運費率一直在低位徘徊.

2.3 BP濾波實證分析

周期循環(huán)波動是以數(shù)年或數(shù)月為周期的一種變動，它可能是一種景氣變動也可能是經(jīng)濟變動或者其他周期變動.對HP濾波分離出的周期項進行BP濾波分析.周期項序列即為隱含周期的平穩(wěn)時間循環(huán)序列，其平穩(wěn)性檢驗結(jié)果見表1.其p值為0.001 4，小于0.05，因此循環(huán)序列為平穩(wěn)序列，可以對其進行BP濾波分析.

用EViews對循環(huán)項進行BP濾波分析，按照截斷點取值的一般原則通過比較分析，取截斷點n為28，得到循環(huán)序列的頻率響應(yīng)函數(shù)圖見圖4.

表1 循環(huán)序列ADF檢驗結(jié)果

圖4 BSI序列周期項(循環(huán))的頻率響應(yīng)函數(shù)圖

由圖4可見，BSI月度值序列低頻率處(f約為0.012 5)功率譜值較高，頻率響應(yīng)函數(shù)在大于等于0.012 5的頻率區(qū)間取值應(yīng)為1，說明長周期的變動所占比例較高，即BSI月度序列以中長周期波動為主.周期為頻率的倒數(shù)，即BSI月度序列存在80個月左右的波動周期.超靈便型船6.7 a的周期變動,與相關(guān)學(xué)者的受全球經(jīng)濟不景氣影響航運周期少于7 a的結(jié)論相呼應(yīng).

3 結(jié) 論

通過對BSI月度序列的研究可以得出以下結(jié)論與建議：

(1)超靈便型船舶市場存在“夏高冬低”的季節(jié)波動特性，該特性與干散貨航運市場“夏低冬高”的季節(jié)特性剛好相反，需要引起超靈便型船舶經(jīng)營者的注意.

(2)超靈便型船舶的長期趨勢表現(xiàn)為以2008年經(jīng)濟危機為轉(zhuǎn)折點的“鐘形”趨勢，2008年附近月份BSI波動異常.

(3)超靈便型船舶市場的周期波動以中長期波動為主，其運價指數(shù)存在80個月左右的周期性波動.

(4)超靈便型船舶經(jīng)營者可以根據(jù)船舶航運市場“夏高冬低”的特性并結(jié)合當年航運市場的總體行情，選擇性地增減投入市場的運力，比如通過租入或買進二手船等措施增加夏季航運市場上的運力投入以獲益；租出或賣出二手船等措施減少冬季市場上運力的投入以減少損失.另外，船舶經(jīng)營者可以根據(jù)航運周期合理選擇經(jīng)營策略，規(guī)避市場風(fēng)險.

參考文獻：

[1] XU J J, YIP T L, MARLOE P B. The dynamics between freight volatility and fleet size growth in dry bulk shipping markets[J]. Transportation Res: Part E, 2011, 47(6): 983-991.

[2] DURU O. A multivariate model of fuzzy integrated logical forecasting method(M-FILF) and multiplicative time series clusting: a model of time-varying volatility for dry cargo freight market[J]. Expert Systems with Applications, 2012, 39(4): 4135-4142.

[3] GEMAN H, SMITU W. Shipping markets and freight rates: an analysis of the Baltic dry index[J]. J Alternative Investments, 2012, 15(1): 98-109.

[4] 聶金龍, 李序穎. 波羅的海干散貨運價指數(shù)的ARFIMA模型研究[J]. 中國水運, 2009(4): 57-58.

[5] 李耀鼎, 宗蓓華. 波羅的海運價指數(shù)波動研究[J]. 上海海事大學(xué)學(xué)報, 2006, 27(12): 84-87.

[6] 陸克從, 趙剛, 胡佳驊. ARCH族模型在干散貨航運指數(shù)分析中的應(yīng)用[J]. 系統(tǒng)工程, 2008, 26(9): 50-56.

[7] 張頁, 王穎. 小靈便型船運價指數(shù)的分析與預(yù)測[J]. 中國航海, 2011, 34(1): 96-99.

[8] 林國龍, 陳言誠. 基于CF濾波的國際干散貨航運市場周期性分析[J]. 上海海事大學(xué)學(xué)報, 2012, 33(3): 69-74.

[9] 王磊, 李序穎. 海岬型船平均收益與新造船價格波動相關(guān)性的互譜分析[J]. 上海海事大學(xué)學(xué)報, 2010, 31(3): 74-77.

[10] 韓荊. 基于Census X12與HP濾波下的現(xiàn)金投放分析[J]. 工業(yè)技術(shù)經(jīng)濟, 2009, 28(10): 143-144.

[11] 鄭衛(wèi)國, 陳萍. 我國經(jīng)濟周期波動的實證研究——一個基于HP濾波的實證研究[J]. 海南金融, 2008(3): 8-12.

[12] 徐國祥, 王芳. 我國房地產(chǎn)市場周期波動譜分析及其實證研究[J]. 統(tǒng)計研究, 2010, 27(10): 18-24

[13] 高鐵梅. 計量經(jīng)濟方法與建模[M]. 北京: 清華大學(xué)出版社, 2009.

[14] 李敏, 陳勝可. EViews統(tǒng)計分析與應(yīng)用[M]. 北京: 電子工業(yè)出版社, 2010.