吳金卓 王卉 劉濱凡

(東北林業(yè)大學(xué)，哈爾濱，150040) (黑龍江省林業(yè)科學(xué)院)

黑龍江省近些年進(jìn)口木質(zhì)林產(chǎn)品數(shù)量顯著增加。根據(jù)黑龍江省統(tǒng)計(jì)年鑒，1999年黑龍江省進(jìn)口原木和鋸材合計(jì)183.9萬(wàn)m3、2014年進(jìn)口717.3萬(wàn)m3、2018年進(jìn)口1 089.2萬(wàn)m3(鋸材按照原木當(dāng)量系數(shù)1.3折算)[1-3]。而黑龍江省木材產(chǎn)量逐年降低，1999年為976.03萬(wàn)m3、2014年下降為212.8萬(wàn)m3、2018年僅為70.6萬(wàn)m3[4-6]。進(jìn)口木材在黑龍江省木材供給市場(chǎng)中占有相當(dāng)大的份額，且呈現(xiàn)逐年上升的趨勢(shì)。木材加工業(yè)是黑龍江省經(jīng)濟(jì)發(fā)展的優(yōu)勢(shì)產(chǎn)業(yè)之一，進(jìn)口木質(zhì)林產(chǎn)品對(duì)于黑龍江省經(jīng)濟(jì)發(fā)展具有重要的影響。

進(jìn)口木質(zhì)林產(chǎn)品受到多種因素的影響，主要?dú)w結(jié)為國(guó)內(nèi)對(duì)于木材的需求、相關(guān)林業(yè)政策的實(shí)施以及木制成品的出口需求[7-8]。目前，多數(shù)研究依據(jù)國(guó)家層面的數(shù)據(jù)對(duì)木材進(jìn)口影響因素進(jìn)行分析[9-12]，而針對(duì)黑龍江省木材進(jìn)口的研究，多集中在對(duì)俄木材進(jìn)口貿(mào)易的可持續(xù)性、對(duì)俄木材貿(mào)易與環(huán)境協(xié)調(diào)發(fā)展等方面[13-15]，較少對(duì)天然林保護(hù)工程實(shí)施后影響黑龍江省進(jìn)口木質(zhì)林產(chǎn)品的主要因素進(jìn)行研究。為此，本研究依據(jù)黑龍江省1999—2018年木質(zhì)林產(chǎn)品進(jìn)口量數(shù)據(jù)，建立經(jīng)驗(yàn)?zāi)B(tài)分解模型，對(duì)木質(zhì)林產(chǎn)品進(jìn)口量隨時(shí)間動(dòng)態(tài)變化進(jìn)行分解，得到不同時(shí)間尺度的木質(zhì)林產(chǎn)品進(jìn)口波動(dòng)量以及長(zhǎng)期變化趨勢(shì)，并構(gòu)建不同尺度進(jìn)口波動(dòng)量、長(zhǎng)期變化趨勢(shì)與不同影響因素的回歸模型，分析天然林保護(hù)工程實(shí)施后黑龍江省木質(zhì)林產(chǎn)品進(jìn)口量波動(dòng)的影響因素，旨在為黑龍江省制定相關(guān)的林業(yè)產(chǎn)業(yè)政策提供參考。

1 研究方法

1.1 經(jīng)驗(yàn)?zāi)B(tài)分解法

經(jīng)驗(yàn)?zāi)B(tài)分解(EMD)方法，是由Huang et al.[16]于1998年提出的一種處理非線性、非平穩(wěn)數(shù)據(jù)的信號(hào)分析方法。該方法對(duì)非線性波動(dòng)的信號(hào)進(jìn)行不斷的高頻分量輸出，直到所有的頻率成分都被分離出來(lái)，得到不同時(shí)間尺度的波動(dòng)分量。國(guó)內(nèi)學(xué)者將其應(yīng)用到國(guó)際原油價(jià)格波動(dòng)、股票市場(chǎng)價(jià)格分析、故障特征提取等領(lǐng)域[17-19]。本研究利用經(jīng)驗(yàn)?zāi)B(tài)分解法對(duì)木質(zhì)林產(chǎn)品進(jìn)口量隨時(shí)間動(dòng)態(tài)變化進(jìn)行分解，得到各自的本征模函數(shù)模型與一個(gè)長(zhǎng)期變化趨勢(shì)，分析木質(zhì)林產(chǎn)品進(jìn)口量的波動(dòng)情況。本征模函數(shù)模型按照頻率的平均周期由小到大排列，分別代表木質(zhì)林產(chǎn)品進(jìn)口量在不同時(shí)間尺度時(shí)的波動(dòng)特征。長(zhǎng)期變化趨勢(shì)代表木質(zhì)林產(chǎn)品進(jìn)口量的長(zhǎng)期走勢(shì)，當(dāng)長(zhǎng)期變化趨勢(shì)為單調(diào)遞增函數(shù)時(shí)，木質(zhì)林產(chǎn)品進(jìn)口量長(zhǎng)期處于不斷增長(zhǎng)的狀態(tài)[20]，反之亦然。各本征模函數(shù)值與長(zhǎng)期變化趨勢(shì)量相加仍與林產(chǎn)品進(jìn)口量相等，實(shí)現(xiàn)數(shù)據(jù)重構(gòu)。本研究采用Pearson相關(guān)系數(shù)、Kendall相關(guān)系數(shù)分析不同木質(zhì)林產(chǎn)品進(jìn)口波動(dòng)量、趨勢(shì)量與木質(zhì)林產(chǎn)品進(jìn)口量之間的線性相關(guān)程度，使用快速傅里葉變換(FFT)計(jì)算各個(gè)木質(zhì)林產(chǎn)品進(jìn)口波動(dòng)量的平均周期，并根據(jù)不同木質(zhì)林產(chǎn)品進(jìn)口波動(dòng)量、趨勢(shì)量的方差得出貢獻(xiàn)率，從而構(gòu)成相關(guān)統(tǒng)計(jì)量，將分解結(jié)果進(jìn)一步進(jìn)行定量化。

1.2 多元線性回歸法

多元線性回歸模型，通常研究一個(gè)評(píng)價(jià)指標(biāo)與多個(gè)影響因素的變化關(guān)系，表達(dá)式：yi=β0+β1x1i+β2x2i+…+βjxji+ui，i=1、2、…、n。式中：j為影響因素個(gè)數(shù)；βj(j=0、1、…、k)為回歸系數(shù)；ui為隨機(jī)擾動(dòng)項(xiàng)；n為樣本數(shù)。為了有效地估計(jì)模型中的參數(shù)，上述模型需要滿足一些假定，如：隨機(jī)擾動(dòng)項(xiàng)的均值為零，方差相同，且隨機(jī)擾動(dòng)項(xiàng)之間無(wú)自相關(guān)性；隨機(jī)擾動(dòng)項(xiàng)與影響因素之間不相關(guān)以及隨機(jī)擾動(dòng)項(xiàng)服從正態(tài)分布等。多元線性回歸模型，已經(jīng)被廣泛地應(yīng)用于不同領(lǐng)域中的影響因素分析[21-23]。本研究將以需求因素、供給因素、政策因素、貿(mào)易因素為影響因素，以方差貢獻(xiàn)率較大的木質(zhì)林產(chǎn)品進(jìn)口波動(dòng)量、長(zhǎng)期趨勢(shì)量為評(píng)價(jià)指標(biāo)，建立多元線性回歸模型，分析影響黑龍江省木質(zhì)林產(chǎn)品進(jìn)口的主要因素。

1.3 評(píng)價(jià)指標(biāo)體系的構(gòu)建

根據(jù)2000—2019年《黑龍江省統(tǒng)計(jì)年鑒》中統(tǒng)計(jì)的海關(guān)主要進(jìn)口商品，黑龍江省進(jìn)口的木質(zhì)林產(chǎn)品主要集中在原木、鋸材、木漿、紙(包括紙板)4種產(chǎn)品，因此本研究以黑龍江省1999—2018年原木、鋸材、木漿、紙(包括紙板)的進(jìn)口量為研究對(duì)象，分析影響木質(zhì)林產(chǎn)品進(jìn)口的主要因素。大量研究表明，木質(zhì)林產(chǎn)品進(jìn)口量受多種因素影響，如需求因素、供給因素、政策因素、貿(mào)易因素以及其他因素等[24-27]。考慮到數(shù)據(jù)的可得性，本研究選取黑龍江省固定資產(chǎn)投資(房地產(chǎn)開發(fā))(FAI)、地區(qū)生產(chǎn)總值(GDP)、中國(guó)木材消費(fèi)量(WC)為需求方面的影響因素，選取黑龍江省森林蓄積量(FS)、木材產(chǎn)量(TY)為供給方面的影響因素，選取加入WTO(WTO)、金融危機(jī)爆發(fā)(FC)、俄羅斯木材出口關(guān)稅調(diào)整(RTP)為國(guó)際政策環(huán)境方面的影響因素，選取進(jìn)口原油價(jià)格(ICP)、人民幣實(shí)際有效匯率指數(shù)(RMR)為貿(mào)易方面的影響因素。

1.4 數(shù)據(jù)來(lái)源與處理

本研究中的各個(gè)影響因素?cái)?shù)據(jù)主要由2000—2019年《中國(guó)統(tǒng)計(jì)年鑒》、《中國(guó)林業(yè)發(fā)展報(bào)告》及《黑龍江省統(tǒng)計(jì)年鑒》整理得到；人民幣實(shí)際有效匯率指數(shù)來(lái)源于國(guó)際清算銀行數(shù)據(jù)(以2010年的匯率進(jìn)行了指數(shù)化)；原油價(jià)格由聯(lián)合國(guó)數(shù)據(jù)庫(kù)中進(jìn)口貿(mào)易金額與進(jìn)口貿(mào)易數(shù)量計(jì)算而得；俄羅斯木材出口關(guān)稅調(diào)整(RTP)，根據(jù)俄羅斯關(guān)稅稅則等資料整理，其他國(guó)際政策影響因素以發(fā)生時(shí)間為轉(zhuǎn)折點(diǎn)，在發(fā)生到結(jié)束期間取值為1，其余為0。

2 結(jié)果與分析

2.1 木質(zhì)林產(chǎn)品進(jìn)口量的動(dòng)態(tài)變化

木質(zhì)林產(chǎn)品進(jìn)口量的變化受多種因素的共同影響，進(jìn)口量含有多種時(shí)間尺度的復(fù)雜變化和局部化特征。運(yùn)用矩陣實(shí)驗(yàn)室(MATLAB)軟件分別對(duì)1999—2018年黑龍江省原木、鋸材、木漿、紙(包括紙板)的進(jìn)口量隨時(shí)間動(dòng)態(tài)變化進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解及線性化和平穩(wěn)化處理，分別得到反映進(jìn)口量的多層次時(shí)間尺度和局部特征的進(jìn)口量變化，即：4種木質(zhì)林產(chǎn)品進(jìn)口量的短期波動(dòng)情況、中長(zhǎng)期波動(dòng)情況、20年內(nèi)的長(zhǎng)期變化趨勢(shì)。受短期、中長(zhǎng)期因素的影響，進(jìn)口量波動(dòng)曲線呈正負(fù)幅度差變化(見(jiàn)圖1)。

圖1 黑龍江省木質(zhì)林產(chǎn)品進(jìn)口量不同時(shí)間尺度的動(dòng)態(tài)變化

由圖1(a)可見(jiàn)：原木進(jìn)口量短期波動(dòng)較頻繁，尤其在2002—2007年呈現(xiàn)出明顯的“波峰-波谷-波峰”，變化幅度較大；原木進(jìn)口量中長(zhǎng)期波動(dòng)頻率放緩，在1999—2012年呈現(xiàn)出顯著的先增加后減少的趨勢(shì)，之后雖有上升，但上升幅度小于下降幅度；原木進(jìn)口量的總體趨勢(shì)在1999—2007年呈現(xiàn)上升狀態(tài)，在2007年后下降趨勢(shì)較為明顯。由圖1(b)可見(jiàn)：鋸材進(jìn)口量短期波動(dòng)明顯，尤其在2010—2014年呈現(xiàn)較頻繁的波動(dòng)；鋸材進(jìn)口量中長(zhǎng)期波動(dòng)比較平緩，僅在2011—2014年有所下降，之后持續(xù)回升，上升幅度大于之前的下降幅度；鋸材進(jìn)口量總體呈持續(xù)性上升的態(tài)勢(shì)。由圖1(c)可見(jiàn)：木漿進(jìn)口量短期內(nèi)先呈現(xiàn)大幅度波動(dòng)，之后在2006—2018年連續(xù)呈現(xiàn)小幅鋸齒狀波動(dòng)；木漿進(jìn)口量中長(zhǎng)期波動(dòng)比較均勻，類似于正弦函數(shù)波動(dòng)，變化幅度相對(duì)穩(wěn)定；木漿進(jìn)口量總體趨勢(shì)在1999—2010年持續(xù)下降，在2010年后緩慢上升，但上升幅度小于下降幅度。由圖1(d)可見(jiàn)：紙(包括紙板)進(jìn)口量短期波動(dòng)劇烈，在2004—2011年呈現(xiàn)出顯著的鋸齒狀波動(dòng)；紙(包括紙板)進(jìn)口量中長(zhǎng)期波動(dòng)比較均勻，呈現(xiàn)兩次波峰，變化幅度比較穩(wěn)定；紙(包括紙板)進(jìn)口量總體趨勢(shì)呈持續(xù)性緩慢下降的態(tài)勢(shì)。

由表1可見(jiàn)：原木、鋸材、木漿、紙(包括紙板)的進(jìn)口量波動(dòng)均包含2個(gè)周期，但是周期長(zhǎng)度存在差異。方差貢獻(xiàn)率可以判定各個(gè)波動(dòng)狀態(tài)和長(zhǎng)期趨勢(shì)在木質(zhì)林產(chǎn)品進(jìn)口量中作用大小的程度，比較各自的方差貢獻(xiàn)率可知，原木進(jìn)口量波動(dòng)主要以中長(zhǎng)期波動(dòng)為主，波動(dòng)呈現(xiàn)先上升然后下降再上升的趨勢(shì)；鋸材進(jìn)口量長(zhǎng)期趨勢(shì)量的方差貢獻(xiàn)率達(dá)到了94.11%，說(shuō)明鋸材進(jìn)口量波動(dòng)依賴于長(zhǎng)期趨勢(shì)變化，總體趨勢(shì)以上升為主；木漿進(jìn)口量波動(dòng)主要以短期波動(dòng)為主，波動(dòng)經(jīng)過(guò)短期的大幅上升和下降后維持小幅波動(dòng)的狀態(tài)；紙(包括紙板)進(jìn)口量波動(dòng)主要以長(zhǎng)期趨勢(shì)變化為主，總體呈現(xiàn)持續(xù)性下降的態(tài)勢(shì)。從Pearson相關(guān)系數(shù)和Kendall相關(guān)系數(shù)也可看出，中長(zhǎng)期波動(dòng)與原木進(jìn)口量波動(dòng)的相關(guān)性最強(qiáng)；長(zhǎng)期趨勢(shì)量與鋸材、紙(包括紙板)進(jìn)口量波動(dòng)的相關(guān)性最強(qiáng)；短期波動(dòng)與木漿進(jìn)口量波動(dòng)的相關(guān)性最強(qiáng)。因此，原木進(jìn)口量的波動(dòng)情況，主要來(lái)源于中長(zhǎng)期波動(dòng)；鋸材進(jìn)口量的波動(dòng)，主要來(lái)源于長(zhǎng)期變化趨勢(shì)；木漿進(jìn)口量的波動(dòng)，主要來(lái)源于短期波動(dòng)；紙(包括紙板)進(jìn)口量的波動(dòng)，主要來(lái)源于長(zhǎng)期變化趨勢(shì)。

表1 黑龍江省不同木質(zhì)林產(chǎn)品進(jìn)口量與波動(dòng)狀態(tài)、趨勢(shì)的相關(guān)性

2.2 木質(zhì)林產(chǎn)品進(jìn)口量逐步回歸分析

為了進(jìn)一步定量地分析不同因素對(duì)木質(zhì)林產(chǎn)品進(jìn)口量波動(dòng)的影響程度，本研究對(duì)不同影響因素與木質(zhì)林產(chǎn)品進(jìn)口量方差貢獻(xiàn)率較大的木質(zhì)林產(chǎn)品進(jìn)口波動(dòng)量、長(zhǎng)期趨勢(shì)量進(jìn)行逐步回歸分析(見(jiàn)表2)。

由表2可知，需求因素對(duì)黑龍江省木質(zhì)林產(chǎn)品進(jìn)口量的影響最為顯著。其中，黑龍江省地區(qū)生產(chǎn)總值，與原木進(jìn)口量的中長(zhǎng)期波動(dòng)量、鋸材進(jìn)口量的長(zhǎng)期趨勢(shì)量呈顯著正相關(guān)；中國(guó)木材消費(fèi)量，與木漿進(jìn)口量的短期波動(dòng)量、紙(包括紙板)進(jìn)口量的長(zhǎng)期趨勢(shì)量呈顯著負(fù)相關(guān)。國(guó)際政策環(huán)境因素對(duì)木質(zhì)林產(chǎn)品進(jìn)口量的變化有較顯著的影響。2001年中國(guó)加入WTO，與原木進(jìn)口量的中長(zhǎng)期波動(dòng)量、木漿進(jìn)口量的短期波動(dòng)量呈顯著正相關(guān)，與鋸材進(jìn)口量的長(zhǎng)期趨勢(shì)量呈顯著負(fù)相關(guān)。2007年開始的俄羅斯木材出口關(guān)稅調(diào)整，與原木進(jìn)口量的中長(zhǎng)期波動(dòng)量呈顯著負(fù)相關(guān)。2008年國(guó)際經(jīng)濟(jì)危機(jī)的爆發(fā)，與鋸材、紙(包括紙板)的進(jìn)口量長(zhǎng)期趨勢(shì)量呈顯著負(fù)相關(guān)。貿(mào)易因素也對(duì)木質(zhì)林產(chǎn)品進(jìn)口量的變化有重要的影響。人民幣實(shí)際有效匯率，與鋸材進(jìn)口量的長(zhǎng)期趨勢(shì)量、木漿進(jìn)口量的短期波動(dòng)量呈顯著正相關(guān)，與原木進(jìn)口量的中長(zhǎng)期波動(dòng)量呈顯著負(fù)相關(guān)。進(jìn)口原油價(jià)格，與鋸材進(jìn)口量的長(zhǎng)期趨勢(shì)量呈顯著正相關(guān)。另外，黑龍江省森林蓄積量和木材產(chǎn)量對(duì)黑龍江省木質(zhì)林產(chǎn)品進(jìn)口量的變化沒(méi)有顯著影響，所以這兩個(gè)因素沒(méi)有進(jìn)入回歸方程。綜上分析，黑龍江省原木進(jìn)口量中長(zhǎng)期波動(dòng)，主要受到黑龍江省地區(qū)生產(chǎn)總值、國(guó)際政策環(huán)境和人民幣實(shí)際有效匯率的共同影響；地區(qū)生產(chǎn)總值的增加和人民幣升值，是鋸材進(jìn)口量增長(zhǎng)的主要原因；加入世界貿(mào)易組織，是黑龍江省木漿進(jìn)口量短期快速上升的主要原因；中國(guó)木材消費(fèi)量中，造紙用材林的培育快速增加，因此對(duì)木漿、紙(包括紙板)的進(jìn)口起到抑制作用。

表2 各因素對(duì)木質(zhì)林產(chǎn)品進(jìn)口量波動(dòng)/趨勢(shì)影響的回歸系數(shù)

3 結(jié)論與討論

應(yīng)用經(jīng)驗(yàn)?zāi)B(tài)分解方法，將1999—2018年黑龍江省原木、鋸材、木漿、紙(包括紙板)的進(jìn)口量各自分解為短期波動(dòng)量、中長(zhǎng)期波動(dòng)量和長(zhǎng)期趨勢(shì)量。從各分解量與木質(zhì)林產(chǎn)品進(jìn)口量之間的相關(guān)統(tǒng)計(jì)結(jié)果可知：原木進(jìn)口量波動(dòng)主要以中長(zhǎng)期波動(dòng)為主，波動(dòng)持續(xù)影響時(shí)間較長(zhǎng)，約為10 a；鋸材進(jìn)口量波動(dòng)主要以長(zhǎng)期趨勢(shì)為主，整體趨勢(shì)上升，并伴隨著一定的中短期波動(dòng)；木漿進(jìn)口量波動(dòng)主要以短期波動(dòng)為主，平均周期為5 a；紙(包括紙板)進(jìn)口量波動(dòng)主要以長(zhǎng)期變化為主，整體呈持續(xù)下降的趨勢(shì)。

黑龍省木質(zhì)林產(chǎn)品進(jìn)口量受到國(guó)內(nèi)木材需求、國(guó)際政策環(huán)境以及國(guó)際貿(mào)易等多種因素的影響，其中，需求因素的影響最為顯著，黑龍江省地區(qū)生產(chǎn)總值(GDP)的增加是原木進(jìn)口量中長(zhǎng)期波動(dòng)和鋸材進(jìn)口量長(zhǎng)期增長(zhǎng)的主要原因。這與魏旸艷等[9]通過(guò)回歸分析得出需求因素是影響我國(guó)木材進(jìn)口的最重要因素的結(jié)論相一致。蘇興國(guó)等[27]在構(gòu)建VAR模型分析木材進(jìn)出口量和木材進(jìn)口價(jià)格之間的雙向影響機(jī)制時(shí)，也表明國(guó)內(nèi)旺盛的木材需求是影響木材進(jìn)口的主要原因。國(guó)際政策環(huán)境和貿(mào)易因素，對(duì)木質(zhì)林產(chǎn)品進(jìn)口影響也較為顯著。2001年中國(guó)加入WTO，對(duì)進(jìn)口木材實(shí)行“零關(guān)稅”政策，促進(jìn)了黑龍江省原木和木漿進(jìn)口量的上升；由于原木進(jìn)口量增加，鋸材進(jìn)口量相對(duì)減少。從2007年開始，俄羅斯逐步提高木材的出口關(guān)稅，導(dǎo)致黑龍江省原木進(jìn)口量顯著下降。2008年全球爆發(fā)金融危機(jī)，導(dǎo)致黑龍江省2009年GDP增長(zhǎng)放緩，這也對(duì)鋸材、紙(包括紙板)的進(jìn)口產(chǎn)生了一定的消極影響。人民幣實(shí)際有效匯率對(duì)木質(zhì)林產(chǎn)品進(jìn)口影響顯著。人民幣升值有利于降低購(gòu)買成本，促進(jìn)了黑龍江省的鋸材和木漿的進(jìn)口，這也與林妙萍[28]在探究中國(guó)紙產(chǎn)品貿(mào)易條件及影響因素中，人民幣實(shí)際有效匯率和貿(mào)易條件之間存在較強(qiáng)的正向相關(guān)效應(yīng)的的觀點(diǎn)吻合。與其他方面的因素相比，木材供給因素對(duì)黑龍江省木質(zhì)林產(chǎn)品進(jìn)口的影響并不顯著。實(shí)際上，黑龍江省從1998年開始有計(jì)劃地對(duì)天然林實(shí)行禁伐和限伐，木材產(chǎn)量持續(xù)下降。隨著2012年黑龍江省開始嚴(yán)厲打擊非法砍伐、2014年對(duì)于龍江森工集團(tuán)和大興安嶺林業(yè)集團(tuán)全面停止木材商業(yè)性采伐政策的落實(shí)以及2018年天保工程二期建設(shè)的結(jié)束，森林植被逐步恢復(fù)，森林蓄積量趨于穩(wěn)定性上升，木材供給潛力逐漸恢復(fù)，但是還無(wú)法轉(zhuǎn)化為有效的木材供給能力，因此森林蓄積量對(duì)木質(zhì)林產(chǎn)品進(jìn)口量的影響較小。

本研究以天然林保護(hù)工程為分析背景，對(duì)黑龍江省木質(zhì)林產(chǎn)品進(jìn)口的波動(dòng)及其影響因素進(jìn)行分析，但是并未考慮林業(yè)碳匯項(xiàng)目對(duì)木質(zhì)林產(chǎn)品進(jìn)口變化的影響。另外，本研究只分析了供給、需求、國(guó)際政策環(huán)境和貿(mào)易方面的常用指標(biāo)對(duì)于木質(zhì)林產(chǎn)品進(jìn)口的影響變化，在未來(lái)的研究中，可以將林業(yè)碳匯項(xiàng)目及其他方面的因素考慮其中進(jìn)一步深入探討和分析。