張 濱 呂潔華

(東北林業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院 哈爾濱 150040)

林業(yè)經(jīng)濟(jì)增長(zhǎng)是林業(yè)事業(yè)發(fā)展過(guò)程中的重要議題，如何有效推動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)與林業(yè)經(jīng)濟(jì)轉(zhuǎn)型升級(jí)一直是國(guó)有林區(qū)改革過(guò)程中所面臨的核心問(wèn)題。隨著“天保工程”等林業(yè)政策的實(shí)施，森林資源開(kāi)發(fā)利用受到了諸多限制，以木質(zhì)林產(chǎn)品為主的勞動(dòng)密集型林業(yè)經(jīng)濟(jì)發(fā)展模式幾近崩潰，以非木質(zhì)林產(chǎn)品為主的資本密集型林業(yè)產(chǎn)業(yè)布局尚未完成，國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)發(fā)展面臨著新的挑戰(zhàn)。在這種情況下，《林業(yè)產(chǎn)業(yè)發(fā)展“十三五”規(guī)劃》中提出“推進(jìn)林業(yè)產(chǎn)業(yè)集聚發(fā)展，通過(guò)優(yōu)化生產(chǎn)要素、調(diào)整產(chǎn)業(yè)結(jié)構(gòu)等途徑來(lái)完成經(jīng)濟(jì)轉(zhuǎn)型和升級(jí)”。為了響應(yīng)國(guó)家林業(yè)產(chǎn)業(yè)發(fā)展號(hào)召，黑龍江省重點(diǎn)國(guó)有林區(qū)也發(fā)布了《黑龍江省森工重點(diǎn)國(guó)有林區(qū)經(jīng)濟(jì)和社會(huì)發(fā)展“十三五”規(guī)劃綱要》，提出在林業(yè)總產(chǎn)值和增加值顯著提高的目標(biāo)下，推進(jìn)林業(yè)產(chǎn)業(yè)向生態(tài)產(chǎn)業(yè)、非木質(zhì)林產(chǎn)品產(chǎn)業(yè)和境外森林資源開(kāi)發(fā)3個(gè)方向轉(zhuǎn)型，形成林業(yè)產(chǎn)業(yè)轉(zhuǎn)型新格局。

實(shí)際上，自“天保工程”實(shí)施后，黑龍江省重點(diǎn)國(guó)有林區(qū)積極嘗試經(jīng)濟(jì)體制改革與經(jīng)濟(jì)結(jié)構(gòu)轉(zhuǎn)型，林區(qū)開(kāi)始出現(xiàn)經(jīng)濟(jì)穩(wěn)定增長(zhǎng)、結(jié)構(gòu)調(diào)整加速、轉(zhuǎn)型發(fā)展路徑日漸清晰等良好發(fā)展勢(shì)頭，但是，“天保工程”等林業(yè)政策所帶來(lái)的負(fù)面沖擊無(wú)法在短時(shí)期內(nèi)得到有效修復(fù)，林區(qū)林業(yè)經(jīng)濟(jì)發(fā)展仍存在很多問(wèn)題。一是全面停伐后，林業(yè)局收入減少，職工收入增長(zhǎng)緩慢，產(chǎn)生了大量剩余勞動(dòng)力； 二是“天保工程”資金投入使得林業(yè)經(jīng)濟(jì)發(fā)展產(chǎn)生了極大的外部依賴性，林業(yè)經(jīng)濟(jì)增長(zhǎng)嚴(yán)重缺乏內(nèi)生動(dòng)力； 三是森林資源利用模式尚未完全轉(zhuǎn)變，林區(qū)資源優(yōu)勢(shì)、生態(tài)優(yōu)勢(shì)依然無(wú)法有效轉(zhuǎn)化為經(jīng)濟(jì)優(yōu)勢(shì)； 四是林業(yè)產(chǎn)業(yè)布局不盡合理，主導(dǎo)產(chǎn)業(yè)不突出，產(chǎn)業(yè)規(guī)?；?、集約化的發(fā)展?jié)摿€沒(méi)有釋放出來(lái)。那么，如何解決上述問(wèn)題，推動(dòng)林業(yè)產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型升級(jí)，形成林業(yè)經(jīng)濟(jì)發(fā)展新格局？這一系列問(wèn)題都可以歸結(jié)為林業(yè)經(jīng)濟(jì)增長(zhǎng)問(wèn)題和林業(yè)經(jīng)濟(jì)轉(zhuǎn)型問(wèn)題，在眾多解決方案中，從空間效應(yīng)角度解析林業(yè)經(jīng)濟(jì)增長(zhǎng)問(wèn)題，為破解林業(yè)經(jīng)濟(jì)發(fā)展困境提供了新思路。本文從空間角度出發(fā)，探討要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)在空間層面對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)存在怎樣的影響效果，是否存在不同形式的空間效應(yīng)，以及產(chǎn)生空間效應(yīng)的內(nèi)在原因和表現(xiàn)形式。

圍繞林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程中的空間效應(yīng)與投入產(chǎn)出等問(wèn)題，現(xiàn)有文獻(xiàn)已經(jīng)呈現(xiàn)較為豐富的研究成果。大部分學(xué)者認(rèn)為資源稟賦、要素投入、空間區(qū)位、林業(yè)產(chǎn)業(yè)集聚、林業(yè)產(chǎn)業(yè)結(jié)構(gòu)等對(duì)推動(dòng)林業(yè)產(chǎn)出增長(zhǎng)具有顯著影響(孔凡斌等， 2013； 魏素豪等， 2017； 程鵬飛等， 2018)，勞動(dòng)力和資本仍是推動(dòng)林業(yè)產(chǎn)出增長(zhǎng)的主要影響因素，并且具有規(guī)模報(bào)酬遞增的特征(胡晨沛， 2017)，林業(yè)投資主要在短期內(nèi)有效拉動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)(才琪等， 2015)； 林業(yè)產(chǎn)業(yè)結(jié)構(gòu)變動(dòng)同樣是林業(yè)經(jīng)濟(jì)增長(zhǎng)的重要?jiǎng)恿?lái)源(羅小鋒等， 2017)，急需通過(guò)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與產(chǎn)業(yè)升級(jí)等方式推動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)(曹玉昆等， 2018； 劉鐵鐸等， 2018)，通過(guò)產(chǎn)品深加工等產(chǎn)業(yè)升級(jí)方式提高木材的附加價(jià)值(Sandbergetal.， 2014)，或者通過(guò)發(fā)展森林旅游業(yè)有效帶動(dòng)經(jīng)濟(jì)增長(zhǎng)(譚益民等， 2012)。林業(yè)資本和勞動(dòng)力投入的增加對(duì)我國(guó)林業(yè)產(chǎn)出增長(zhǎng)的貢獻(xiàn)率逐漸降低(李研等， 2018)。林業(yè)產(chǎn)業(yè)集聚不僅對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)存在顯著影響，也與很多因素存在潛在關(guān)系，林業(yè)產(chǎn)業(yè)集聚與就業(yè)趨勢(shì)存在密切聯(lián)系，投資補(bǔ)貼、基礎(chǔ)設(shè)施狀況等也對(duì)林業(yè)產(chǎn)業(yè)集聚有很大影響(Kleinetal.， 2009； Kiesetal.， 2010)，林業(yè)產(chǎn)業(yè)集聚對(duì)投資的影響最大(汪浩， 2011)。也有學(xué)者認(rèn)為林業(yè)資源稟賦對(duì)林業(yè)產(chǎn)業(yè)集聚有顯著的非空間正向影響，人力資源條件則有顯著的空間負(fù)向影響(魏肖杰等， 2018)，林業(yè)產(chǎn)業(yè)間的共生關(guān)系有助于提升林業(yè)系統(tǒng)穩(wěn)定性和經(jīng)濟(jì)效益(Karlssonetal.， 2008)。另外，相關(guān)學(xué)者對(duì)黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)發(fā)展的實(shí)證研究顯示，人均林業(yè)總產(chǎn)值和人均林業(yè)增加值均存在顯著的空間正自相關(guān)，但是空間集聚水平較低(趙曉光等， 2016； 2017)，黑龍江省的林業(yè)產(chǎn)業(yè)結(jié)構(gòu)變動(dòng)對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)的貢獻(xiàn)較大(鐘艷等， 2011)，勞動(dòng)力和資本是推動(dòng)黑龍江省國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)的主要?jiǎng)恿?潘啟英等， 2018)。

上述研究成果表明要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)等能夠有效解釋林業(yè)經(jīng)濟(jì)增長(zhǎng)，空間區(qū)位與林業(yè)經(jīng)濟(jì)增長(zhǎng)也存在密切聯(lián)系，大部分研究利用全局或局部自關(guān)聯(lián)Moran指數(shù)，缺乏考慮不同因素的交叉影響，空間溢出效應(yīng)的解析也沒(méi)有從根本上回答空間溢出效應(yīng)的本質(zhì)是什么。從林業(yè)經(jīng)濟(jì)增長(zhǎng)的機(jī)理來(lái)看，勞動(dòng)力、資本、森林資源等投入要素增加是推動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)的主要原因，同時(shí)，產(chǎn)業(yè)結(jié)構(gòu)的轉(zhuǎn)型升級(jí)也會(huì)優(yōu)化林產(chǎn)品產(chǎn)出結(jié)構(gòu)，從而帶動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)； 將空間因素引入林業(yè)經(jīng)濟(jì)增長(zhǎng)問(wèn)題時(shí)，由于存在空間依賴關(guān)系和空間異質(zhì)性，生產(chǎn)要素和相關(guān)信息具備空間流動(dòng)條件，有利于資源的合理配置，進(jìn)而衍生出林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程中的空間關(guān)聯(lián)效應(yīng)、空間集聚效應(yīng)和空間溢出效應(yīng)等內(nèi)容，空間關(guān)聯(lián)表明空間區(qū)位對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)存在影響，空間集聚會(huì)產(chǎn)生“中心-外圍”等問(wèn)題，從而產(chǎn)生以空間溢出效應(yīng)為主的外部經(jīng)濟(jì)。本文在上述研究和機(jī)理分析的基礎(chǔ)上，將交叉Moran指數(shù)引入空間關(guān)聯(lián)分析和空間集聚分析中，證明要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)存在空間交叉關(guān)聯(lián)和空間集聚影響； 進(jìn)而建立空間回歸模型，討論要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)在林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程中的空間溢出效應(yīng)及其表現(xiàn)形式，為推動(dòng)黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)與林業(yè)產(chǎn)業(yè)轉(zhuǎn)型升級(jí)提供一定的參考依據(jù)。

1)由于林業(yè)投資統(tǒng)計(jì)口徑變化，缺少2003—2010年林業(yè)投資完成額，利用2011—2017年40個(gè)林業(yè)局林業(yè)固定資產(chǎn)投資額和林業(yè)投資完成額的平均比例關(guān)系，對(duì)林業(yè)固定資產(chǎn)投資額乘以平均比例系數(shù)1.594 6以替代缺失數(shù)據(jù)。

1 模型、數(shù)據(jù)與變量說(shuō)明

1.1 Moran指數(shù)

Moran指數(shù)是觀測(cè)值與關(guān)聯(lián)指標(biāo)空間滯后值的相關(guān)系數(shù)(Moran， 1950)。考察全局空間關(guān)聯(lián)效應(yīng)時(shí)，可以計(jì)算全局自關(guān)聯(lián)Moran指數(shù)I和全局交叉關(guān)聯(lián)Moran指數(shù)Ixy，在Moran散點(diǎn)圖上根據(jù)散點(diǎn)分布和Moran指數(shù)判斷觀測(cè)值與關(guān)聯(lián)指標(biāo)之間是否存在顯著的空間關(guān)聯(lián)效應(yīng)。考察某一特定區(qū)域與鄰近區(qū)域之間的空間集聚特征時(shí)，可以計(jì)算局部自關(guān)聯(lián)Moran指數(shù)Ii和局部交叉關(guān)聯(lián)Moran指數(shù)Ii,xy，以判斷是否存在空間集聚效應(yīng)以及哪些指標(biāo)對(duì)空間集聚存在顯著影響(Anselin， 1995)。

Moran指數(shù)大于0時(shí)表示正相關(guān)，即高值被高值包圍或低值被低值包圍； 小于0表示負(fù)相關(guān)，即高值被低值包圍或低值被高值包圍； 接近0時(shí)表示隨機(jī)空間分布，不存在空間關(guān)系。

對(duì)Moran指數(shù)的可靠性進(jìn)行統(tǒng)計(jì)檢驗(yàn)時(shí)，標(biāo)準(zhǔn)化Moran指數(shù)為Z統(tǒng)計(jì)量(方差是n次隨機(jī)生成，伴隨概率是偽P值，n越大誤差越小，本文取n=999)，Z服從漸進(jìn)標(biāo)準(zhǔn)正態(tài)分布。

1.2 空間面板回歸模型

空間面板回歸模型是探討空間溢出效應(yīng)的有效方法(Elhorst， 2014)，一般形式如下：

式中：Y是因變量向量，X是自變量矩陣，W是空間權(quán)重矩陣，ξ是殘差向量，ε是白噪聲向量且服從均值為0、方差為σ2In的正態(tài)分布(0是零向量，In是單位矩陣)，β是非空間項(xiàng)系數(shù)，δ、ρ和λ表示空間項(xiàng)系數(shù)，u和v表示固定效應(yīng)和隨機(jī)效應(yīng)。

當(dāng)λ=0，ρ=0,δ=0時(shí)，為普通面板回歸模型(ordinary least square，OLS)；

當(dāng)λ=0，ρ≠0,δ=0時(shí)，為空間自回歸模型(spatial auto-regression model，SAR)；

當(dāng)λ≠0，ρ=0,δ=0時(shí)，為空間誤差模型(spatial error model，SEM)；

當(dāng)λ=0，ρ≠0,δ≠0時(shí)，為空間杜賓模型(spatial Durbin model，SDM)。

在對(duì)比不同模型時(shí)，除了R2、T值等常見(jiàn)統(tǒng)計(jì)量，還包括對(duì)數(shù)極大似然值(log likelihood)、判斷使用SAR還是SEM的拉格朗日乘數(shù)(Lagrange multiplier，LM)和穩(wěn)健拉格朗日乘數(shù)(robust Lagrange multiplier，R-LM)； 似然比(likelihood ratio，LR)和沃爾德檢驗(yàn)(Wald test)則是判斷SDM是否優(yōu)于SAR和SEM； 豪斯曼檢驗(yàn)判斷使用固定效應(yīng)還是隨機(jī)效應(yīng)； 進(jìn)一步基于LeSage等(2009)的方法可以將變量的邊際效應(yīng)分解為本地效應(yīng)和溢出效應(yīng)。

1.3 數(shù)據(jù)來(lái)源與變量說(shuō)明

數(shù)據(jù)來(lái)源于2003—2017年《中國(guó)林業(yè)統(tǒng)計(jì)年鑒》和龍江森工集團(tuán)內(nèi)部統(tǒng)計(jì)資料，樣本量600個(gè)(15年40個(gè)林業(yè)局)。林業(yè)經(jīng)濟(jì)發(fā)展水平用林業(yè)總產(chǎn)值(Y)表示； 要素投入包括勞動(dòng)力、投資和森林資源稟賦，其中林業(yè)勞動(dòng)力用在崗職工人數(shù)(L)表示，林業(yè)投資用林業(yè)固定資產(chǎn)投資額和林業(yè)投資完成額(K)表示1)，森林資源稟賦用森林管護(hù)面積(F)表示1)； 產(chǎn)業(yè)結(jié)構(gòu)升級(jí)(industrial structure upgrading，ISU)則通過(guò)木質(zhì)-非木質(zhì)產(chǎn)業(yè)結(jié)構(gòu)表示，借鑒徐德云(2008)設(shè)計(jì)的產(chǎn)業(yè)結(jié)構(gòu)升級(jí)指標(biāo)，將產(chǎn)業(yè)結(jié)構(gòu)升級(jí)定義為ISU=木質(zhì)產(chǎn)業(yè)產(chǎn)值比重+非木質(zhì)產(chǎn)業(yè)產(chǎn)值比重*2，ISU∈[1,2]。對(duì)林業(yè)總產(chǎn)值、林業(yè)投資分別剔除物價(jià)影響(以2003年黑龍江省CPI為基準(zhǔn))。另外，空間權(quán)重矩陣使用一階“后”相鄰矩陣(W)，“后”相鄰是指2個(gè)區(qū)域有共同邊界或頂點(diǎn)(相鄰取1，不相鄰取0)。描述性統(tǒng)計(jì)見(jiàn)表1。

表1 各變量的描述性統(tǒng)計(jì)結(jié)果Tab.1 Descriptive statistical results of variables

在計(jì)算空間關(guān)聯(lián)效應(yīng)和空間集聚效應(yīng)時(shí)，林業(yè)總產(chǎn)值(ΔY/Y，單位：%，下同)、勞動(dòng)力投入(ΔL/L)和林業(yè)投資(ΔK/K)采用年增長(zhǎng)率，森林管護(hù)面積取自然對(duì)數(shù)(lnF)； 在計(jì)算回歸模型時(shí)，對(duì)主要投入產(chǎn)出指標(biāo)取自然對(duì)數(shù)； 產(chǎn)業(yè)結(jié)構(gòu)升級(jí)均使用ISU。另外，自變量間的相關(guān)系數(shù)均為中低度相關(guān)，內(nèi)生性問(wèn)題對(duì)結(jié)果影響較小。

1)森林管護(hù)面積是指森林資源管理、保護(hù)、培育和利用等活動(dòng)所對(duì)應(yīng)的森林面積，其中包括人工培育、養(yǎng)殖等林下多種經(jīng)營(yíng)活動(dòng)所對(duì)應(yīng)的森林面積，因此，以森林管護(hù)面積衡量森林資源稟賦水平明顯優(yōu)于森林面積。另外，森林管護(hù)面積數(shù)據(jù)連續(xù)性和穩(wěn)定性較好，受投資、勞動(dòng)力的干擾較小，具有良好的替代性。

2 結(jié)果與分析

2.1 空間關(guān)聯(lián)效應(yīng)

計(jì)算2017年林業(yè)總產(chǎn)值增長(zhǎng)率的全局自關(guān)聯(lián)Moran指數(shù)(Geoda軟件)，圖1(a)顯示40個(gè)林業(yè)局存在3個(gè)明顯的離群值(葦河、山河屯、鐵力)，導(dǎo)致Moran指數(shù)為-0.024 1(偽P=0.448 0)； 從圖1(a)來(lái)看大部分林業(yè)局都比較集中，剔除離群值后，如圖1(b)所示，修正后的全局自關(guān)聯(lián)Moran指數(shù)為0.271 8(偽P=0.001 0)，即2017年林業(yè)總產(chǎn)值增長(zhǎng)整體上存在顯著的空間關(guān)聯(lián)效應(yīng)。

圖1 2017年林業(yè)總產(chǎn)值增長(zhǎng)率(ΔY/Y)全局自關(guān)聯(lián)Moran散點(diǎn)分布Fig.1 Global auto-correlation Moran scatter distribution of forestry gross output growth rate(ΔY/Y) in 2017

同理，計(jì)算全部Moran指數(shù)，當(dāng)存在圖1中情況時(shí)也剔除相應(yīng)離群值(表2)。從林業(yè)總產(chǎn)值增長(zhǎng)的全局自關(guān)聯(lián)Moran指數(shù)來(lái)看， 2004—2017年林業(yè)局林業(yè)總產(chǎn)值增長(zhǎng)整體表現(xiàn)為顯著的正向空間自關(guān)聯(lián)； 從全局交叉關(guān)聯(lián)Moran指數(shù)來(lái)看，大部分在0.10顯著性水平上成立，其中，林業(yè)投資增長(zhǎng)、森林管護(hù)面積、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)與林業(yè)總產(chǎn)值增長(zhǎng)之間整體上為正的空間交叉關(guān)聯(lián)，勞動(dòng)力投入增長(zhǎng)與林業(yè)總產(chǎn)值增長(zhǎng)之間整體上為負(fù)的空間交叉關(guān)聯(lián)。剔除離群值后，林業(yè)總產(chǎn)值增長(zhǎng)整體存在顯著空間關(guān)聯(lián)效應(yīng)，林業(yè)總產(chǎn)值增長(zhǎng)不僅顯著受周?chē)謽I(yè)局林業(yè)經(jīng)濟(jì)增長(zhǎng)影響，還暗示對(duì)投資、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)和森林管護(hù)面積有較強(qiáng)的依賴性，與勞動(dòng)力增長(zhǎng)的負(fù)向交叉關(guān)聯(lián)效應(yīng)反映出勞動(dòng)力投入存在冗余。空間全局Moran指數(shù)結(jié)果意味著林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程中存在空間集聚和空間溢出等問(wèn)題。

表2 全局自關(guān)聯(lián)和交叉關(guān)聯(lián)Moran指數(shù)(修正后)①Tab.2 Global auto-correlation and cross-correlation Moran indexes (Adjusted)

2.2 空間集聚效應(yīng)

為驗(yàn)證林業(yè)經(jīng)濟(jì)增長(zhǎng)是否存在顯著的空間集聚效應(yīng)，計(jì)算2017年40個(gè)林業(yè)局局部Moran指數(shù)(表3)。從局部Moran指數(shù)值來(lái)看，林口、鐵力、美溪、紅星、葦河5個(gè)林業(yè)局均存在顯著的空間相鄰特征(Moran指數(shù)大于0.10或者小于-0.10)，其他林業(yè)局也在不同程度上存在空間集聚或空間分散特征，其中，ΔY/Y-ΔY/Y、ΔY/Y-ISU的高增長(zhǎng)空間集聚均以綏棱為中心，ΔY/Y-ΔL/L主要以綏棱、沾河、湯旺河、雙鴨山為中心，ΔY/Y-ΔK/K主要以綏陽(yáng)、大海林、穆棱為中心，ΔY/Y-lnF主要以通北、綏棱、八面通為中心(0.10顯著性水平)。

表3 局部Moran指數(shù)及其分類(lèi)①Tab.3 Local Moran indexes and classifications

從40個(gè)林業(yè)局Moran指數(shù)的正、負(fù)值來(lái)看，ΔY/Y-ΔY/Y的空間相鄰特征中空間集聚和空間分散各有20個(gè)林業(yè)局，林業(yè)經(jīng)濟(jì)增長(zhǎng)存在局部空間集聚特征； ΔY/Y-ΔL/L主要是空間分散特征，ΔY/Y-ΔK/K、ΔY/Y-lnF、ΔY/Y-ISU主要是空間集聚特征(正值個(gè)數(shù)大于負(fù)值個(gè)數(shù))。另外，在0.05、0.10、0.20和0.30顯著性水平上對(duì)集聚特征進(jìn)行分類(lèi)，結(jié)果表明： 除ΔY/Y-ΔL/L外，其他Moran指數(shù)在 0.05和0.10顯著性水平上空間集聚特征整體不顯著，在0.20和0.30顯著性水平上空間集聚特征逐漸清晰。結(jié)合40個(gè)林業(yè)局的空間地理位置，在0.30顯著性水平上考察空間集聚趨勢(shì)，林業(yè)總產(chǎn)值增長(zhǎng)呈現(xiàn)較為顯著的“高增長(zhǎng)”包圍“低增長(zhǎng)”的空間集聚特征，即40個(gè)林業(yè)局內(nèi)部地區(qū)林業(yè)總產(chǎn)值增長(zhǎng)緩慢，外圍地區(qū)增長(zhǎng)較快； 勞動(dòng)力投入所帶來(lái)的空間分散效應(yīng)最為明顯(a+bc+d)。

空間集聚分析結(jié)果表明黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)的空間集聚特征相對(duì)較弱，即空間集聚效應(yīng)正處于開(kāi)始形成的階段。從增長(zhǎng)集聚形成動(dòng)因來(lái)看，勞動(dòng)力的流動(dòng)(負(fù)向交叉關(guān)聯(lián)和空間分散)使得人力資源得到合理配置，提高了人均勞動(dòng)生產(chǎn)率，對(duì)林業(yè)總產(chǎn)值增長(zhǎng)集聚有顯著的帶動(dòng)作用(0.10顯著性水平)； 林業(yè)投資和森林資源的空間流動(dòng)性較差，存在非生產(chǎn)性林業(yè)投資和純粹的森林管護(hù)活動(dòng)，從而對(duì)林業(yè)總產(chǎn)值增長(zhǎng)集聚的推動(dòng)作用較弱； 產(chǎn)業(yè)結(jié)構(gòu)升級(jí)反映林業(yè)生產(chǎn)結(jié)構(gòu)變化，產(chǎn)業(yè)結(jié)構(gòu)升級(jí)尚未完成，無(wú)法形成有效的規(guī)模經(jīng)濟(jì)以帶動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)的空間集聚。

2.3 空間溢出效應(yīng)

林業(yè)經(jīng)濟(jì)增長(zhǎng)的普通面板回歸模型(OLS)、空間自回歸模型(SAR)、空間誤差模型(SEM)和空間杜賓模型(SDM)如表4所示。從林業(yè)經(jīng)濟(jì)增長(zhǎng)的普通面板回歸模型來(lái)看，應(yīng)建立固定效應(yīng)模型，回歸系數(shù)依次為： 林業(yè)產(chǎn)業(yè)結(jié)構(gòu)升級(jí)(1.161 2)、森林管護(hù)面積增加(0.954 1)、林業(yè)投資增長(zhǎng)(0.251 3)、勞動(dòng)力投入增長(zhǎng)(-0.301 3)，與之對(duì)應(yīng)， 2003—2017年40個(gè)林業(yè)局產(chǎn)業(yè)結(jié)構(gòu)升級(jí)、森林管護(hù)面積、林業(yè)投資、勞動(dòng)力投入的年平均變化量分別為0.023 8(單位： 1)、0.488 2%、8.967 6%、-3.421 7%，乘以回歸系數(shù)得到對(duì)應(yīng)變量年平均貢獻(xiàn)率依次為0.027 6%、0.465 8%、2.253 6%和1.031 0%，在不考慮空間效應(yīng)的前提下，現(xiàn)階段驅(qū)動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)的核心動(dòng)力仍是林業(yè)投資和勞動(dòng)力投入，森林管護(hù)與產(chǎn)業(yè)結(jié)構(gòu)升級(jí)的貢獻(xiàn)率相對(duì)較小。

表4 林業(yè)經(jīng)濟(jì)增長(zhǎng)的面板回歸結(jié)果①Tab.4 Panel regression results of forestry economic growth

將空間因素引入面板回歸模型中，建立空間自回歸模型(SAR)和空間誤差模型(SEM)以探討林業(yè)總產(chǎn)值增長(zhǎng)是否存在空間溢出效應(yīng)。SAR假定空間溢出效應(yīng)僅存在于因變量，SEM假定空間溢出效應(yīng)通過(guò)誤差形式表現(xiàn)出來(lái)。LM值和R-LM值顯示SAR優(yōu)于SEM，相關(guān)統(tǒng)計(jì)指標(biāo)也表明固定效應(yīng)最佳。此時(shí)，空間自回歸固定效應(yīng)模型中林業(yè)總產(chǎn)值空間項(xiàng)的系數(shù)為0.530 0，即外圍林業(yè)局林業(yè)總產(chǎn)值每增長(zhǎng)1%，將平均帶動(dòng)本地林業(yè)局林業(yè)總產(chǎn)值增長(zhǎng)0.530 0%，林業(yè)總產(chǎn)值增長(zhǎng)存在顯著的空間溢出效應(yīng)，有必要進(jìn)一步探討因變量的空間溢出效應(yīng)是否源于自變量的空間溢出效應(yīng)。

建立空間杜賓模型(SDM)以探討要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)是否存在空間溢出效應(yīng)。LR值和Wald檢驗(yàn)顯示SDM優(yōu)于SAR和SEM，同樣固定效應(yīng)最佳。從回歸系數(shù)來(lái)看，除勞動(dòng)力的非空間項(xiàng)外，其他系數(shù)均通過(guò)T檢驗(yàn)，要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)均存在本地效應(yīng)和溢出效應(yīng)。另外，SDM同樣表明林業(yè)總產(chǎn)值存在空間溢出效應(yīng)。

對(duì)比上述模型結(jié)果，可以發(fā)現(xiàn)隨著空間項(xiàng)的引入，空間面板回歸模型中自變量對(duì)因變量的解釋能力逐漸增強(qiáng)，空間杜賓模型表明要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)和林業(yè)經(jīng)濟(jì)增長(zhǎng)均存在空間溢出效應(yīng)。由于空間杜賓模型中回歸系數(shù)包含滯后效應(yīng)、反饋效應(yīng)等，所以本地效應(yīng)、溢出效應(yīng)與非空間項(xiàng)、空間項(xiàng)系數(shù)略有不同。因此，需要對(duì)空間杜賓固定效應(yīng)模型的邊際效應(yīng)進(jìn)行分解重構(gòu)，將邊際效應(yīng)轉(zhuǎn)化為直接效應(yīng)和間接效應(yīng)，以準(zhǔn)確衡量空間溢出效應(yīng)(表5)。一般認(rèn)為，直接效應(yīng)是本地效應(yīng)，即要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)對(duì)林業(yè)局自身林業(yè)經(jīng)濟(jì)增長(zhǎng)的邊際效應(yīng)，間接效應(yīng)是溢出效應(yīng)，是受外圍林業(yè)局的邊際效應(yīng)。

從表5來(lái)看，空間杜賓模型中自變量的溢出效應(yīng)均明顯大于本地效應(yīng)，即要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)主要通過(guò)空間溢出效應(yīng)發(fā)揮作用； 通過(guò)對(duì)比OLS固定效應(yīng)系數(shù)(表4)與總效應(yīng)系數(shù)，表明空間項(xiàng)的引入重新解析了要素投入與產(chǎn)業(yè)結(jié)構(gòu)升級(jí)的邊際效應(yīng)； 另外， 2003—2017年勞動(dòng)力、林業(yè)投資、森林管護(hù)面積、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)的年平均貢獻(xiàn)率分別為2.097 2%、2.231 1%、0.889 8%和0.021 1%，與OLS相比，勞動(dòng)力投入、森林管護(hù)面積的年平均貢獻(xiàn)率顯著增大，林業(yè)投資增長(zhǎng)、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)的年平均貢獻(xiàn)率略微降低，平均貢獻(xiàn)率總和的變大在于模型解釋能力增強(qiáng)，即空間杜賓模型能更好地解釋林業(yè)總產(chǎn)值增長(zhǎng)，突顯空間溢出效應(yīng)的重要性?？臻g溢出效應(yīng)的存在表明林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程不是封閉的而是開(kāi)放的，林業(yè)局之間的“聯(lián)系”對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)產(chǎn)生了很大影響，林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程中存在顯著的外部經(jīng)濟(jì)。

表5 空間杜賓模型的邊際效應(yīng)分解①Tab.5 Decomposition of marginal effects in spatial Durbin model

一般情況下，空間溢出效應(yīng)主要通過(guò)要素流動(dòng)和知識(shí)溢出表現(xiàn)出來(lái)。國(guó)有林區(qū)經(jīng)營(yíng)管理體制的特殊性導(dǎo)致生產(chǎn)要素在林業(yè)局之間的流動(dòng)性很小，即要素流動(dòng)不是空間溢出的主要來(lái)源，黑龍江省重點(diǎn)國(guó)有林區(qū)要素投入與產(chǎn)業(yè)結(jié)構(gòu)升級(jí)的空間溢出主要屬于知識(shí)溢出，即林業(yè)局之間相互學(xué)習(xí)、模仿、競(jìng)爭(zhēng)，既有為滿足現(xiàn)實(shí)發(fā)展需求的主動(dòng)獲取，也有源于林業(yè)政策實(shí)施的被動(dòng)接受。從生產(chǎn)內(nèi)容與經(jīng)營(yíng)模式來(lái)看，隨著林業(yè)保護(hù)政策的實(shí)施，林業(yè)生產(chǎn)內(nèi)容整體上已經(jīng)調(diào)整為以非木質(zhì)林產(chǎn)品生產(chǎn)為主、木質(zhì)林產(chǎn)品加工為輔，各個(gè)林業(yè)局生產(chǎn)內(nèi)容與經(jīng)營(yíng)模式趨于一致； 從勞動(dòng)力-資本結(jié)構(gòu)來(lái)看， 2003年40個(gè)林業(yè)局勞動(dòng)力-投資比值的標(biāo)準(zhǔn)差為0.421 9(人·萬(wàn)元-1)， 2017年為0.093 1，勞動(dòng)力-投資結(jié)構(gòu)差異顯著降低； 同樣，木質(zhì)-非木質(zhì)產(chǎn)業(yè)結(jié)構(gòu)的標(biāo)準(zhǔn)差由0.127 2降為0.050 4(單位： 1)，產(chǎn)出結(jié)構(gòu)差異顯著降低，產(chǎn)業(yè)結(jié)構(gòu)升級(jí)得到同步發(fā)展。因此，黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)的空間溢出效應(yīng)本質(zhì)上是知識(shí)溢出，具體表現(xiàn)在3個(gè)方面： 生產(chǎn)內(nèi)容與經(jīng)營(yíng)模式的同化、要素配置結(jié)構(gòu)的趨同、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)的同步，并且在知識(shí)溢出過(guò)程中伴隨著林業(yè)經(jīng)濟(jì)增長(zhǎng)的空間關(guān)聯(lián)效應(yīng)與空間集聚效應(yīng)。

3 結(jié)論

黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)存在顯著的空間關(guān)聯(lián)效應(yīng)，并開(kāi)始形成高增長(zhǎng)包圍低增長(zhǎng)的空間集聚效應(yīng)； 要素投入、產(chǎn)業(yè)結(jié)構(gòu)升級(jí)和林業(yè)總產(chǎn)值增長(zhǎng)均存在顯著的空間溢出效應(yīng)，并且主要通過(guò)知識(shí)溢出的形式表現(xiàn)出來(lái)，空間溢出效應(yīng)作為間接效應(yīng)，相比本地效應(yīng)，對(duì)推動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)具有更強(qiáng)的驅(qū)動(dòng)效果； 森林管護(hù)面積和產(chǎn)業(yè)結(jié)構(gòu)升級(jí)對(duì)林業(yè)經(jīng)濟(jì)增長(zhǎng)有較為顯著的邊際作用，但是勞動(dòng)力和資本投入仍是推動(dòng)林業(yè)經(jīng)濟(jì)增長(zhǎng)的核心動(dòng)力。另外，林業(yè)勞動(dòng)力存在投入冗余，亟需通過(guò)勞動(dòng)力轉(zhuǎn)移與產(chǎn)業(yè)結(jié)構(gòu)升級(jí)吸收剩余勞動(dòng)力?？紤]到要素投入結(jié)構(gòu)與產(chǎn)業(yè)結(jié)構(gòu)的實(shí)際變動(dòng)，現(xiàn)階段黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)仍屬于典型的要素投入驅(qū)動(dòng)型，并且開(kāi)始從“勞動(dòng)密集型”向“資本密集型”過(guò)渡。

黑龍江省重點(diǎn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)過(guò)程中的空間溢出效應(yīng)屬于知識(shí)溢出為主的空間趨同，這種外部經(jīng)濟(jì)必須借助信息傳遞才能實(shí)現(xiàn)，因此，提升林業(yè)經(jīng)濟(jì)增長(zhǎng)的空間效應(yīng)可以通過(guò)林業(yè)局之間的知識(shí)傳遞、交流學(xué)習(xí)帶動(dòng)林業(yè)優(yōu)勢(shì)產(chǎn)業(yè)集聚收斂，形成規(guī)模經(jīng)濟(jì)和新的林業(yè)經(jīng)濟(jì)增長(zhǎng)極。同時(shí)，通過(guò)優(yōu)化要素投入結(jié)構(gòu)，引導(dǎo)林區(qū)職工創(chuàng)新就業(yè)，將生態(tài)建設(shè)與林業(yè)多種經(jīng)營(yíng)緊密結(jié)合，充分發(fā)揮國(guó)有林區(qū)森林資源稟賦優(yōu)勢(shì)，實(shí)現(xiàn)非木質(zhì)林產(chǎn)品產(chǎn)業(yè)的規(guī)?；?、集約化，深化產(chǎn)業(yè)轉(zhuǎn)型升級(jí)，進(jìn)一步增強(qiáng)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)的內(nèi)生增長(zhǎng)動(dòng)力，從要素配置、產(chǎn)業(yè)結(jié)構(gòu)、空間區(qū)位等多個(gè)方面共同促進(jìn)國(guó)有林區(qū)林業(yè)經(jīng)濟(jì)增長(zhǎng)。