伍海泉，印 慧，曹玉昆

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

據(jù)第九次全國森林資源清查（2014～2018年）結(jié)果，我國森林面積為2.2 億hm2，森林覆蓋率為22.96%，森林蓄積量為175.6 億m3，實(shí)現(xiàn)30年連續(xù)保持森林面積、森林蓄積量“雙增長”[1]。大氣二氧化碳（CO2）濃度升高是影響陸地生態(tài)系統(tǒng)可持續(xù)性的主要全球變化因子之一[2]，而森林則是陸地生態(tài)系統(tǒng)中重要的碳氮儲存庫。應(yīng)繼續(xù)加大力度推進(jìn)生態(tài)文明建設(shè)，堅(jiān)持節(jié)約自然資源與保護(hù)生態(tài)環(huán)境。20 世紀(jì)80年代以來，我國提出并實(shí)施多項(xiàng)生態(tài)計(jì)劃，包括退耕還林工程、長江和珠江流域及農(nóng)田防護(hù)林體系建設(shè)工程、天然林保護(hù)工程等。1998年長江流域和松花江、嫩江流域特大洪災(zāi)后，黨中央、國務(wù)院決定在云南等12 個省（區(qū)、市）重點(diǎn)國有林區(qū)開展天然林保護(hù)工程試點(diǎn)。2000年《長江上游、黃河上中游地區(qū)天然林資源保護(hù)工程實(shí)施方案》經(jīng)國務(wù)院批準(zhǔn)實(shí)施，《東北、內(nèi)蒙古等重點(diǎn)國有林區(qū)天然林資源保護(hù)工程實(shí)施方案》隨即出臺。至此，天然林保護(hù)工程在我國17 個?。▍^(qū)、市）重點(diǎn)林區(qū)實(shí)施，具體包括長江上游地區(qū)云南、四川、貴州、重慶、湖北、西藏6 ?。▍^(qū)、市），黃河上中游地區(qū)陜西、甘肅、青海、寧夏、內(nèi)蒙古、山西、河南7 ?。▍^(qū)），以及黑龍江、吉林、海南、新疆4?。▍^(qū)）。天然林保護(hù)工程實(shí)施是針對我國長期以來天然林資源過度消耗引起生態(tài)環(huán)境惡化，國家從社會經(jīng)濟(jì)可持續(xù)發(fā)展戰(zhàn)略高度做出的重大決策[3]。天然林保護(hù)工程作為全國最大林業(yè)生態(tài)保護(hù)工程之一，實(shí)施20年來，中央財(cái)政總投入已超過4 000 億元[4]，林區(qū)自然資源與社會面貌發(fā)生巨大變化。

目前，天然林保護(hù)一期工程已結(jié)束，二期工程進(jìn)入關(guān)鍵時段。如何客觀評價天然林保護(hù)工程實(shí)施效果，哪些因素對其影響顯著，是該項(xiàng)生態(tài)保護(hù)工程尚待解決的問題。以往研究主要集中在以下三個方面：一是天然林保護(hù)工程實(shí)踐綜述，研究內(nèi)容包括天然林保護(hù)工程實(shí)施手段和取得成績。陳麗榮等通過繪制可視化知識圖譜，展示天然林保護(hù)工程研究領(lǐng)域時間分布和研究熱點(diǎn)演進(jìn)，揭示我國天然林保護(hù)工程研究整體脈絡(luò)[5]；曹玉昆等通過對我國天然林保護(hù)工程效益評價研究文獻(xiàn)梳理，系統(tǒng)分析天然林保護(hù)工程效益評價主要特點(diǎn)和方法[6]。二是天然林保護(hù)工程績效和價值研究，例如秦利等通過構(gòu)建天然林保護(hù)工程政策實(shí)施績效綜合評價指標(biāo)體系，采用模糊綜合評價法對伊春林區(qū)天保工程二期政策實(shí)施績效進(jìn)行綜合評價[7]；崔方茹等利用靖邊縣天保工程相關(guān)資料和既有成果，構(gòu)建森林生態(tài)系統(tǒng)服務(wù)功能評估指標(biāo)體系，評估靖邊縣天保工程森林生態(tài)系統(tǒng)服務(wù)功能價值[8]。三是天保工程相關(guān)影響研究。如臧良震運(yùn)用西部11 個省級面板數(shù)據(jù)，采用面板向量自回歸和脈沖響應(yīng)分析研究工程建設(shè)對農(nóng)民收入產(chǎn)生影響[9]；翟緒軍等采用黑龍江省國有林區(qū)林業(yè)產(chǎn)業(yè)數(shù)據(jù)，分析實(shí)施天然林保護(hù)工程對國有林區(qū)林業(yè)產(chǎn)業(yè)總體結(jié)構(gòu)變動產(chǎn)生的影響[10]。已有報(bào)道多是通過文獻(xiàn)歸納整理對天然林保護(hù)工程定性評價分析為主，缺乏以計(jì)量核算為主的定量評價，且研究對象或面板數(shù)據(jù)多集中于西部、東北部等區(qū)域天然林保護(hù)工程，針對全國天然林保護(hù)工程政策實(shí)施效果評價尤其與林業(yè)全要素生產(chǎn)率相結(jié)合研究較少。因此，本文使用多期雙重差分法（Difference in difference，DID），運(yùn)用我國1997～2018年省級面板數(shù)據(jù)，基于林業(yè)全要素生產(chǎn)率視角，天然林保護(hù)工程項(xiàng)目實(shí)施效果作出評價，為我國森林資源經(jīng)營管理和生態(tài)資源保護(hù)提供參考。

1 模型設(shè)定

1.1 方法選取與模型設(shè)定

為研究天然林保護(hù)工程對林業(yè)全要素生產(chǎn)率政策效應(yīng)，本文選取雙重差分法（DID）對政策實(shí)施地區(qū)作驗(yàn)證。雙重差分模型適用基礎(chǔ)是先將觀測樣本準(zhǔn)確分為試驗(yàn)組和對照組，本文選擇實(shí)施天然林保護(hù)工程省份為“試驗(yàn)組”，其他省份作為相應(yīng)“對照組”，設(shè)置DID回歸方程如下：

上述公式被解釋變量中，設(shè)置政策影響變量為yit；觀測組別為d i，分別代表試驗(yàn)組和對照組；觀測年份為dt；μit是殘差項(xiàng)。在借鑒相關(guān)研究前提下，本文在原基礎(chǔ)模型中加入部分控制變量。拓展后DID方程如下:

方程中，下標(biāo)p代表省份，t代表觀測年份。在Treatpt賦值中，甘肅、貴州、海南、河南、黑龍江、湖北、吉林、內(nèi)蒙古、寧夏、青海、山西、陜西、四川、西藏、新疆、云南、重慶地區(qū)為試驗(yàn)組1，其他地區(qū)為對照組0。Yearpt為天保工程實(shí)施年份，1998年及以后賦值為1，其他為0；Xpt為控制變量；γp為地區(qū)固定效應(yīng)；ηt為年份固定效應(yīng)；μpt為隨機(jī)擾動項(xiàng)。β1為重點(diǎn)關(guān)注系數(shù)，Treatpt×Yearpt則是反映政策凈效應(yīng)交乘變量。其次，參數(shù)β2表示無天然林保護(hù)工程政策試驗(yàn)組和對照組之間目標(biāo)變量時間變化情況，參數(shù)β3表示試驗(yàn)組與對照組之間不隨時間變動差異，式（2）中包含γp地區(qū)虛擬變量和ηt年份虛擬變量。

1.2 指標(biāo)選取與數(shù)據(jù)來源

本文被解釋變量為通過全要素生產(chǎn)率動態(tài)變化的Malmquist指數(shù)DEA模型測度我國31個省級林業(yè)全要素生產(chǎn)率，Malmquist指數(shù)以t期為基期，公式如下：

上述公式中，MALM為從t到t+1期林業(yè)全要素生產(chǎn)率、PECH為從t到t+1 期純技術(shù)效率變動指數(shù)、SECH為從t到t+1 期規(guī)模效率變動指數(shù)、TECHCH為從t到t+1 期技術(shù)進(jìn)步變動指數(shù)、EF?FCH為從t到t+1 期技術(shù)效率變動指數(shù)。xt和xt+1分別是t到t+1期單位投入向量，yt和yt+1分別是t到t+1 期單位產(chǎn)出向量，Dt(xt,yt)、Dt+1(xt+1,yt+1)分別是以t和t+1 時期技術(shù)為參照距離函數(shù)。上述公式中下標(biāo)v是可變規(guī)模報(bào)酬假設(shè)，c 是不變規(guī)模報(bào)酬假設(shè)，計(jì)算結(jié)果如表1所示。

表1 1997～2018年我國31個省級林業(yè)全要素生產(chǎn)率Malmquist指數(shù)及分解Table1 Malmquist index and decomposition of average total factor productivity of forestry in 31 provinces in China from 1997 to 2018

本文參考相關(guān)研究成果后選取系列控制變量是在兼顧數(shù)據(jù)指標(biāo)可得性前提下，基于天然林保護(hù)工程實(shí)施制度、管理和環(huán)境3個維度，確定經(jīng)濟(jì)發(fā)展水平、城鎮(zhèn)化水平、對外開放水平、科技發(fā)展水平、產(chǎn)業(yè)現(xiàn)代化程度、政府干預(yù)程度、自然災(zāi)害情況、林業(yè)人力資本水平等變量。確定變量值為經(jīng)計(jì)算處理比率值，分別采用地區(qū)人均GDP對數(shù)代表經(jīng)濟(jì)發(fā)展水平，常住城鎮(zhèn)人口占總?cè)丝诒戎氐某擎?zhèn)化率衡量城鎮(zhèn)化水平，進(jìn)出口貿(mào)易額比重代表對外開放程度，研究與試驗(yàn)發(fā)展經(jīng)費(fèi)投入代表科技發(fā)展水平，林業(yè)第二、三產(chǎn)業(yè)占林業(yè)產(chǎn)業(yè)總產(chǎn)值比重代表產(chǎn)業(yè)結(jié)構(gòu)，中央政府投資占林業(yè)總投資比重代表政府干預(yù)程度，森林病蟲害發(fā)生面積代表自然災(zāi)害情況，在崗職工年平均工資代表林業(yè)人力資本投入水平。數(shù)據(jù)主要來源于1998～2018年林業(yè)統(tǒng)計(jì)年鑒和國家統(tǒng)計(jì)局?jǐn)?shù)據(jù)庫。本文選擇1998～2018年我國31 個省級地區(qū)，合計(jì)5 859 個樣本開展經(jīng)驗(yàn)驗(yàn)證，具體變量描述性統(tǒng)計(jì)如表2所示。

表2 1997～2018年樣本描述性統(tǒng)計(jì)Table 2 Descriptive statistics of samples from 1997 to 2018

續(xù)表

2 實(shí)證結(jié)果

2.1 實(shí)證結(jié)果與差分分析

本文在模型測算中逐步加入控制變量以研究天然林保護(hù)工程政策凈效應(yīng)，估計(jì)結(jié)果如表3 所示。在天然林保護(hù)工程實(shí)施制度、管理和環(huán)境等控制變量逐步加入后，解釋變量系數(shù)符號和顯著性并未發(fā)生本質(zhì)變化，說明估計(jì)結(jié)果穩(wěn)健性良好。由于實(shí)證結(jié)果為0.0866～0.1057，在5%置信水平顯著，說明天然林保護(hù)工程實(shí)施顯著提高山西、內(nèi)蒙古、吉林等17 個政策實(shí)施省級地區(qū)林業(yè)全要素生產(chǎn)效率。在影響天然林保護(hù)工程提高生產(chǎn)率控制變量中，產(chǎn)業(yè)結(jié)構(gòu)、人力資本、科技發(fā)展水平和經(jīng)濟(jì)發(fā)展水平影響顯著。其中，人力資本投入程度在5%置信水平顯著，產(chǎn)業(yè)結(jié)構(gòu)、科技和經(jīng)濟(jì)發(fā)展水平在10%置信水平顯著。分析結(jié)果顯示控制變量中人力資本投入程度最為顯著，主要是因隨著科技發(fā)展水平提高和天然林保護(hù)工程資金投入力度加大，人均生產(chǎn)效率得到顯著提高，直接推動林業(yè)全生產(chǎn)要素提高；另一方面，林區(qū)經(jīng)濟(jì)發(fā)展水平提高帶來工資收入水平上升，人均勞動力成本加大，也促使相關(guān)企業(yè)和經(jīng)營主體提高管理水平，實(shí)現(xiàn)減員增效。同時，產(chǎn)業(yè)結(jié)構(gòu)、科技和經(jīng)濟(jì)發(fā)展水平提升一定程度上對沖禁止商業(yè)性砍伐對林業(yè)企業(yè)收入影響，開辟新財(cái)源，為提升林業(yè)全要素生產(chǎn)率提供資金支持、物質(zhì)基礎(chǔ)和科技保障。與人力資本、產(chǎn)業(yè)結(jié)構(gòu)、科技發(fā)展水平和經(jīng)濟(jì)發(fā)展水平相比，城鎮(zhèn)化水平、對外開放程度、政府干預(yù)和自然災(zāi)害情況在林業(yè)全要素生產(chǎn)效率中影響不顯著。從區(qū)域經(jīng)濟(jì)結(jié)構(gòu)特點(diǎn)看，天然林保護(hù)工程實(shí)施加快第三產(chǎn)業(yè)占比提高，但由于林區(qū)氣候、土地等資源限制以及《森林法》等相關(guān)法律法規(guī)規(guī)制，歷史上第一產(chǎn)業(yè)所占比例較低，農(nóng)業(yè)人口比重相對較低。天然林保護(hù)工程實(shí)施引起農(nóng)業(yè)人口向城鎮(zhèn)人口轉(zhuǎn)化空間有限，并不足構(gòu)成顯著影響。同時，改革開發(fā)以來，我國對木材等林業(yè)產(chǎn)品進(jìn)口數(shù)量逐年提升，但林業(yè)相關(guān)資金、勞動力和科技等生產(chǎn)要素跨國界流動并不顯著。對外開放程度變化主要受林區(qū)內(nèi)生經(jīng)濟(jì)增長和外部貿(mào)易環(huán)境影響，與林業(yè)全要素生產(chǎn)率相關(guān)程度低。天然林保護(hù)工程實(shí)施區(qū)域大部分為國有林區(qū)，政府干預(yù)程度較高，政策實(shí)施并未顯著提高政府干預(yù)水平。天然林保護(hù)工程實(shí)施增加政府對林區(qū)發(fā)展資金投入，但資金相當(dāng)比例被用于消耗性和保障性支出，受林業(yè)管理水平提升影響不顯著。

表3 DID 回歸結(jié)果分析Table 3 Regression results

2.2 分組驗(yàn)證結(jié)果與差分分析

為進(jìn)一步研究天然林保護(hù)工程實(shí)施過程中資源配置效率和影響因素作用機(jī)制，本文將各變量再次分組。首先將我國31 個省域按照1997～2018年年均森林覆蓋率進(jìn)行排名，然后將森林覆蓋率排名前15且實(shí)施天然林保護(hù)工程省級地區(qū)，即海南、云南、黑龍江、吉林、陜西、湖北、四川、重慶8個地區(qū)分為試驗(yàn)組，其余年均森林覆蓋率未排全國前15 且實(shí)施天然林保護(hù)工程省級地區(qū)，即貴州、內(nèi)蒙古、河南、山西、西藏、寧夏、青海、新疆、甘肅9個地區(qū)為對照組。首先對森林覆蓋率排名前15 名8 個省級地區(qū)即試驗(yàn)組進(jìn)行異質(zhì)性分析（見表4），回歸結(jié)果說明，森林資源覆蓋率高與天然林保護(hù)工程對該地區(qū)林業(yè)全要素生產(chǎn)率影響有正向作用，政策實(shí)施效果影響更為顯著。

表4 異質(zhì)性分析和森林資源覆蓋率分組回歸結(jié)果Table 4 Heterogeneity analysis and grouping regression results of forest resource coverage

式中，minori為虛擬變量，若i省是森林資源排名前15，則minori=1，否則為0。

森林覆蓋率是反映森林資源豐富程度重要參考指標(biāo)，天然林保護(hù)工程對林業(yè)全要素生產(chǎn)率影響與森林覆蓋率顯著正相關(guān)，說明政策實(shí)施后所在地區(qū)森林蓄積、面積、覆蓋率顯著提高，森林結(jié)構(gòu)趨于合理、森林質(zhì)量明顯改善[11]。同時天然林保護(hù)工程投入巨額資金，顯著加強(qiáng)政策實(shí)施省域森林生態(tài)系統(tǒng)、濕地生態(tài)系統(tǒng)、荒漠生態(tài)系統(tǒng)建設(shè)和生物多樣性保護(hù)，促進(jìn)該地區(qū)森林資源和林業(yè)產(chǎn)業(yè)協(xié)調(diào)發(fā)展。

3 結(jié)論與政策啟示

3.1 結(jié)論

本文以我國31個省域林業(yè)全要素生產(chǎn)率（TFP）Malmquist 指數(shù)為被解釋變量，基于天然林保護(hù)工程實(shí)施政策、管理和環(huán)境等維度，將天然林保護(hù)工程這一重點(diǎn)生態(tài)工程實(shí)施視為一項(xiàng)政策沖擊，通過DID模型對政策實(shí)施省域和未實(shí)施省域分組測算分析，引入經(jīng)濟(jì)發(fā)展水平、城鎮(zhèn)化水平、對外開放水平、科技發(fā)展水平、產(chǎn)業(yè)現(xiàn)代化程度、政府干預(yù)程度、自然災(zāi)害情況、林業(yè)人力資本水平等其他變量因素。結(jié)果表明，天然林保護(hù)工程實(shí)施整體上對林業(yè)全要素生產(chǎn)率影響顯著。人力投入減少、產(chǎn)業(yè)結(jié)構(gòu)合理化、科技投入和經(jīng)濟(jì)發(fā)展水平提高，可有效提升要素資源在林業(yè)生產(chǎn)之間配置效率，城鎮(zhèn)化水平和對外開放程度提高、政府干預(yù)和自然災(zāi)害情況對林業(yè)全要素生產(chǎn)效率影響并不顯著。同時，森林覆蓋率高地區(qū)，天然林保護(hù)工程對其林業(yè)全要素生產(chǎn)率提升影響更加顯著。充分說明天然林保護(hù)工程政策舉措和資金投入，有效促進(jìn)政策實(shí)施省域森林結(jié)構(gòu)更合理、森林資源更豐富。

3.2 政策啟示

①全面總結(jié)天然林保護(hù)工程實(shí)施經(jīng)驗(yàn)，適度擴(kuò)大工程實(shí)施范圍。優(yōu)化天保工程資金支出結(jié)構(gòu)，同時積極填補(bǔ)森林撫育資金缺口，適度上調(diào)森林撫育補(bǔ)貼標(biāo)準(zhǔn)，通過擴(kuò)大森林撫育補(bǔ)貼規(guī)模和范圍，逐步提高森林資源質(zhì)量加快低質(zhì)森林資源改造速度。鼓勵開發(fā)林業(yè)生態(tài)產(chǎn)品，走綠色發(fā)展之路，加快經(jīng)濟(jì)發(fā)展[12]，保護(hù)林業(yè)生態(tài)多樣性，開展“山水林田湖草”整體開發(fā)，提升林區(qū)生態(tài)產(chǎn)品開發(fā)利用效率。

②積極培養(yǎng)壯大林業(yè)科研隊(duì)伍，加速林業(yè)科技成果轉(zhuǎn)化。建設(shè)培養(yǎng)專業(yè)化的林業(yè)技術(shù)應(yīng)用和推廣團(tuán)隊(duì)，應(yīng)用現(xiàn)代化信息技術(shù)和媒體渠道，幫助基層技術(shù)人員掌握多種技術(shù)和裝備。科學(xué)編制人才培養(yǎng)計(jì)劃，創(chuàng)新人才使用機(jī)制，培養(yǎng)高素質(zhì)林業(yè)工程管理人員、高層次科研人員、高水平技術(shù)人員、高技能產(chǎn)業(yè)工人。完善營林管護(hù)技術(shù)手段，支持利用衛(wèi)星遙感探測、無人機(jī)巡航等技術(shù)手段，擴(kuò)大人均管護(hù)面積，提高人均工作效率，全面保護(hù)天然林。

③加強(qiáng)森林資源生態(tài)保護(hù)，提升地區(qū)森林覆蓋率。要加大管護(hù)和執(zhí)法力度，采取有效手段保護(hù)天然林；加強(qiáng)人工林封育管理，加強(qiáng)退化林分改造，提高有林地郁閉度和生物多樣性。同時加大國有林區(qū)基礎(chǔ)設(shè)施建設(shè)支持力度，將電網(wǎng)與飲水安全以及管護(hù)站點(diǎn)用房等基礎(chǔ)設(shè)施建設(shè)納入上級政府建設(shè)規(guī)劃統(tǒng)籌安排，同時將國有林區(qū)道路按屬性納入相關(guān)公路網(wǎng)規(guī)劃。依法依規(guī)嚴(yán)管森林資源，將非林業(yè)部門管理林地納入森林資源統(tǒng)計(jì)管理，加強(qiáng)林地監(jiān)管，確保林地面積只增不減。