姜 微，劉俊昌，胡 皓

（1.北京林業(yè)大學(xué) 經(jīng)濟(jì)管理學(xué)院，北京 100083；2.中南林業(yè)科技大學(xué) 商學(xué)院，湖南 長(zhǎng)沙 410004）

十九大報(bào)告中提出堅(jiān)持人與自然和諧共生，“綠水青山就是金山銀山”的理念，為建設(shè)美麗中國(guó)創(chuàng)造良好的生活環(huán)境[1]。林業(yè)對(duì)于生態(tài)文明建設(shè)的重要意義不言而喻，是生態(tài)文明建設(shè)的重要抓手。林業(yè)對(duì)于改善生態(tài)環(huán)境以及促進(jìn)區(qū)域綠色和可持續(xù)發(fā)展有著不可替代的重要作用。林業(yè)不僅產(chǎn)生經(jīng)濟(jì)效益還兼具生態(tài)效益和社會(huì)效益，故如何有效地利用、使用林業(yè)資源是值得探討的問題。但是如何保證在林業(yè)經(jīng)濟(jì)高速發(fā)展的同時(shí)維持資源消耗和生態(tài)環(huán)境污染在合理范圍內(nèi)是林業(yè)經(jīng)濟(jì)發(fā)展中亟待解決的問題。

1 文獻(xiàn)綜述

國(guó)內(nèi)最早的林業(yè)生態(tài)效率研究是劉璨等人[2]運(yùn)用隨機(jī)前沿生產(chǎn)函數(shù)和超越對(duì)數(shù)生產(chǎn)函數(shù)，分析江蘇省淮安市的農(nóng)田林網(wǎng)和小片林對(duì)農(nóng)業(yè)總產(chǎn)值、種植業(yè)產(chǎn)值和畜牧業(yè)產(chǎn)值的生態(tài)貢獻(xiàn)程度，隨著研究深度和計(jì)量經(jīng)濟(jì)學(xué)等不斷發(fā)展，隨后對(duì)林業(yè)生態(tài)效率研究開展了一系列研究。楊瑋[3]運(yùn)用數(shù)據(jù)包絡(luò)分析法對(duì)中國(guó)各?。▍^(qū)、市）的營(yíng)林產(chǎn)業(yè)和森工產(chǎn)業(yè)的生產(chǎn)效率進(jìn)行了多方面和多層次的研究。白世秀[4]將全要素與偏要素兩者相結(jié)合分析黑龍江省13 個(gè)市的生態(tài)效率；鄭宇梅等[5]運(yùn)用規(guī)模報(bào)酬不變的超效率DEA 分析林業(yè)第二產(chǎn)業(yè)的生態(tài)效率，引入?yún)^(qū)位熵概念分析林業(yè)產(chǎn)業(yè)集聚、區(qū)域經(jīng)濟(jì)效益與生態(tài)環(huán)境間的關(guān)系。羅小鋒等[6]利用ESDA技術(shù)、基尼系數(shù)等方法分析林業(yè)生產(chǎn)效率的差異并分析引起林業(yè)生產(chǎn)效率變化的影響因素。徐瑋等[7]運(yùn)用Malmquist-DEA 模型測(cè)算中國(guó)林業(yè)生產(chǎn)效率并對(duì)不同區(qū)域的林業(yè)生產(chǎn)效率展開對(duì)比分析。史常亮等[8]運(yùn)用SFA（“一步法”隨機(jī)前沿模型）測(cè)算各省林業(yè)技術(shù)效率，然后結(jié)合Malmquist 生產(chǎn)率指數(shù)，分析4 大區(qū)域的林業(yè)TFP 時(shí)序增長(zhǎng)與空間分布特征。黃韶海等[9]采用數(shù)據(jù)包絡(luò)分析方法測(cè)算中國(guó)的林業(yè)生產(chǎn)效率，通過局部空間相關(guān)檢驗(yàn)對(duì)中國(guó)林業(yè)生產(chǎn)的空間集聚進(jìn)行分析，旨在為實(shí)現(xiàn)林業(yè)可持續(xù)發(fā)展提供借鑒。管時(shí)一[10]運(yùn)用數(shù)據(jù)包絡(luò)分析法和空間相關(guān)性分析法對(duì)區(qū)域林業(yè)全要素生產(chǎn)率進(jìn)行縱向發(fā)展分析和橫向比較分析，并研究其地區(qū)收斂性和空間集聚性。魏肖杰等[11]采用DEA-Malmquist 指數(shù)法測(cè)算我國(guó)林業(yè)全要素生產(chǎn)率，運(yùn)用莫蘭指數(shù)驗(yàn)證其空間相關(guān)特征；使用動(dòng)態(tài)空間計(jì)量模型對(duì)林業(yè)產(chǎn)業(yè)集聚影響林業(yè)全要素生產(chǎn)率的機(jī)制進(jìn)行實(shí)證分析。王兆峰[12]利用環(huán)境影響與產(chǎn)品或服務(wù)價(jià)值的比值作為生態(tài)效率的測(cè)算公式，探討湖南武陵山片區(qū)旅游產(chǎn)業(yè)生態(tài)效率的演化過程及其影響因素。

環(huán)境規(guī)制對(duì)生態(tài)效率的影響是復(fù)雜、階段性波動(dòng)的。目前普遍存在兩種觀點(diǎn)：一是“波特假說”[13]，雖然環(huán)境規(guī)制對(duì)生產(chǎn)有抑制作用但是會(huì)反向刺激企業(yè)進(jìn)行創(chuàng)新，促使企業(yè)進(jìn)行技術(shù)升級(jí)和效率的提升；二是“遵循成本說”[14]，認(rèn)為由于環(huán)境規(guī)制的要求使得企業(yè)的生產(chǎn)成本增加進(jìn)而導(dǎo)致企業(yè)利潤(rùn)下降，會(huì)降低企業(yè)的競(jìng)爭(zhēng)力和生產(chǎn)績(jī)效。環(huán)境規(guī)制對(duì)工業(yè)生產(chǎn)效率、農(nóng)業(yè)生產(chǎn)效率的研究頗為豐富[15-16]，但是對(duì)于林業(yè)生態(tài)效率的影響卻少有研究，故環(huán)境規(guī)制對(duì)于林業(yè)生態(tài)效率不同區(qū)域的影響存在抑制作用還是會(huì)倒逼技術(shù)升級(jí)、提高經(jīng)濟(jì)效率還有待進(jìn)一步研究?；诖?，本研究以全國(guó)26 個(gè)?。▍^(qū)、市）為研究對(duì)象，首先利用非徑向方法距離函數(shù)測(cè)算出各區(qū)域的林業(yè)生態(tài)效率，對(duì)不同區(qū)域的林業(yè)生態(tài)效率基于時(shí)空和空間進(jìn)行異質(zhì)性分析，而后運(yùn)用門檻回歸模型對(duì)不同區(qū)域的環(huán)境規(guī)制的門檻效應(yīng)以及空間異質(zhì)性進(jìn)行分析對(duì)比，為后續(xù)林業(yè)經(jīng)濟(jì)發(fā)展中的生態(tài)環(huán)境保護(hù)提供參考信息。

2 研究方法

2.1 非徑向方向距離函數(shù)

非徑向方法距離函數(shù)與DEA 和SBM 效率方法區(qū)別在于，前者在投入要素中加入了能源效率，產(chǎn)出中不僅考慮期望產(chǎn)出，還考慮非期望產(chǎn)出的要素。不同于非期望產(chǎn)出效率測(cè)算中將松弛變量按照同比例縮減（或擴(kuò)張）進(jìn)行測(cè)算[17]，為避免產(chǎn)生“松弛誤偏”問題，特構(gòu)建松弛變量可按不同比例進(jìn)行增減的非期望產(chǎn)出的非徑向方法距離函數(shù)[18-19]。其基本形式如下：

式中：β表示投入、產(chǎn)出的松弛變量，ω是各變量的權(quán)重，g為方向向量，diag(β)表示對(duì)向量β進(jìn)行對(duì)角化處理。Y,Z,X分別表示為投入要素、期望產(chǎn)出和非期望產(chǎn)出。

本文假設(shè)投入要素、期望產(chǎn)出和非期望產(chǎn)出三者同樣重要。對(duì)于林業(yè)產(chǎn)業(yè)生態(tài)效率求解目標(biāo)是在既定的資本和勞動(dòng)要素投入下，期待實(shí)現(xiàn)期望產(chǎn)出最大化以及資源要素投入和環(huán)境要素非期望產(chǎn)出最小化。其實(shí)，投入要素包括勞動(dòng)力投入、資本投入以及能源投入，分別由K，L，E符號(hào)表示；非期望產(chǎn)出包括廢水排放量、廢氣排放量、固體廢棄物排放量，分別由S，Q，W符號(hào)表示。將ω，g代入公式（1），通過線性優(yōu)化求解如下：

求解公式（2）可以得到β的最優(yōu)解即依據(jù)β的最優(yōu)解可以求解出各要素效率，計(jì)算公式如下：

本文將能源效率和環(huán)境效率進(jìn)行加權(quán)，得出林業(yè)生態(tài)效率FEE：

林業(yè)生態(tài)效率指標(biāo)選取要體現(xiàn)林業(yè)產(chǎn)業(yè)的特色，故指標(biāo)選取側(cè)重于林業(yè)特有指標(biāo)，在參考現(xiàn)有文獻(xiàn)的基礎(chǔ)上，多投入多產(chǎn)出決定了林業(yè)生態(tài)效率測(cè)算的多樣性與復(fù)雜性，它涉及到經(jīng)濟(jì)、資源、環(huán)境等指標(biāo)的各個(gè)方面，且各指標(biāo)之間存在錯(cuò)綜復(fù)雜的關(guān)系。為保證林業(yè)生態(tài)效率測(cè)度的真實(shí)性與可信度，本文選取的相關(guān)指標(biāo)體系見表1。

2.2 門檻回歸函數(shù)

環(huán)境規(guī)制對(duì)林業(yè)生態(tài)效率的影響是多維度的，其影響可能會(huì)隨著環(huán)境規(guī)制強(qiáng)度水平處于不同區(qū)間而呈現(xiàn)不同特點(diǎn)，即變量間可能存在非線性關(guān)系。為了檢驗(yàn)變量間是否存在非線性關(guān)系，采用Hansen 提出的面板門檻回歸模型對(duì)上述非線性關(guān)系進(jìn)行檢驗(yàn)?！霸撃Ｐ偷膬?yōu)點(diǎn)在于門檻值的確定是由樣本數(shù)據(jù)內(nèi)生決定，并運(yùn)用Bootstrap 方法估計(jì)門檻值的顯著性，客觀性較強(qiáng)。門檻回歸”（Threshold regression）作為一種非線性計(jì)量經(jīng)濟(jì)學(xué)模型，其實(shí)質(zhì)就是在反映因果關(guān)系的變量中尋找門檻變量，門檻值根據(jù)樣本數(shù)據(jù)進(jìn)行估算，并檢驗(yàn)根據(jù)門檻值劃分的樣本組參數(shù)是否存在著顯著不同。針對(duì)本文的計(jì)量模型，設(shè)定的面板門檻回歸模型如下：

表1 林業(yè)生態(tài)效率投入產(chǎn)出變量及數(shù)據(jù)說明Table 1 Input-output variables and data description of forestry ecological efficiency

式中：FEEi,t為被解釋變量，表示為i地域的第t年的林業(yè)生態(tài)效率值。數(shù)值采用上一節(jié)測(cè)算結(jié)果。ERi,t為核心解釋變量，為環(huán)境規(guī)制強(qiáng)度。環(huán)境規(guī)制測(cè)度方法有多種，本文采用綜合指數(shù)法，選取工業(yè)煙（粉）塵去除率、工業(yè)二氧化硫去除率、工業(yè)固體廢物綜合利用率、生活垃圾無害化處理率、生活污水處理率等5 個(gè)單項(xiàng)指標(biāo)數(shù)據(jù)，利用歸一法和熵值法計(jì)算出環(huán)境規(guī)制強(qiáng)度。Xit為控制變量向量，包括投資規(guī)模（IL），本文用該地區(qū)的林業(yè)固定資產(chǎn)投資數(shù)值進(jìn)行衡量；經(jīng)濟(jì)發(fā)展水平（GDP），為了更好反映出當(dāng)?shù)亟?jīng)濟(jì)發(fā)展水平，本文用當(dāng)?shù)厝司鵊DP 水平進(jìn)行衡量；科技投入（R&D），本研究用林業(yè)研究、試驗(yàn)發(fā)展、信息化等經(jīng)費(fèi)支出表示科技投入；對(duì)外開放程度EO 用外商直接投資額。

β0為常數(shù)項(xiàng)，itε為隨機(jī)擾動(dòng)項(xiàng)。1(·)代表示性函數(shù)，取值規(guī)則如下：當(dāng)表達(dá)式為真時(shí)，取值為1；當(dāng)表達(dá)式為假時(shí)，取值為0。β1，β2分別代表當(dāng)q≤γ和q＞γ時(shí)環(huán)境規(guī)制強(qiáng)度對(duì)林業(yè)生態(tài)效率的彈性系數(shù)，若門檻效應(yīng)存在，則β1，β2的系數(shù)符合或估計(jì)值會(huì)出現(xiàn)顯著差異。

2.3 數(shù)據(jù)來源

文中所涉及解釋變量和被解釋變量的數(shù)據(jù)均來源于《中國(guó)林業(yè)統(tǒng)計(jì)年鑒（2005—2016）》《中國(guó)環(huán)境統(tǒng)計(jì)年鑒（2005—2017）》、各省統(tǒng)計(jì)年鑒以及各省年度統(tǒng)計(jì)報(bào)告。其中由于天津市、上海市歷年相關(guān)數(shù)據(jù)缺失過多，青海、寧夏、西藏的林業(yè)產(chǎn)值過少，所以該文對(duì)以上?。▍^(qū)、市）不予分析。所以該文最終選擇全國(guó)26 個(gè)省（區(qū)、市）作為研究樣本。

3 結(jié)果與分析

3.1 中國(guó)林業(yè)生態(tài)效率測(cè)度結(jié)果與分析

本文運(yùn)用DEAMAX 軟件，按照公式（1）～（5）求解出2005—2016年我國(guó)各省（區(qū)、市）的林業(yè)生態(tài)效率值，結(jié)果如表2所示。其中，林業(yè)生態(tài)效率均值最高為廣東（0.821 1），林業(yè)生態(tài)效率數(shù)值大于0.7 的省市有：上海（0.774 7）、北京（0.707 4）、海南（0.707 0），這些省市大都屬于東部沿海地區(qū)，當(dāng)?shù)亟?jīng)濟(jì)水平較高，以先進(jìn)的制造業(yè)為基礎(chǔ)的現(xiàn)代化工業(yè)體系以及環(huán)保意識(shí)不斷增強(qiáng)和注重環(huán)境保護(hù)，故這些地區(qū)的林業(yè)生態(tài)效率值較高。林業(yè)生態(tài)效率均值最低的為陜西（0.082 2），數(shù)值低于0.3 的?。▍^(qū)、市）有：云南（0.122 2）、黑龍江（0.189 9）、內(nèi)蒙古（0.214 1）、廣西（0.229 7）、四川（0.232 8）、吉林（0.270 4），屬于東北部地區(qū)以及西部地區(qū)。東北部地區(qū)森林資源豐富，但該區(qū)域林業(yè)生產(chǎn)技術(shù)水平較為落后，而又想追求高經(jīng)濟(jì)效益由此惡性循環(huán)導(dǎo)致資源粗放式的開發(fā)利用，造成了資源環(huán)境的嚴(yán)重破壞。云南、廣西地處西部地區(qū)，采取資源換經(jīng)濟(jì)的模式發(fā)展，導(dǎo)致其資源過度開發(fā)利用，對(duì)資源環(huán)境影響較大。

我國(guó)林業(yè)生態(tài)效率整體上呈現(xiàn)從西部到東南部慢慢升高的趨勢(shì)，以東部地區(qū)的林業(yè)生態(tài)效率最高，中部地區(qū)林業(yè)生態(tài)效率值居中；西部、東北部林業(yè)生態(tài)效率值較低。因此，各?。▍^(qū)、市）的林業(yè)生態(tài)效率值分別與各地區(qū)經(jīng)濟(jì)實(shí)力保持基本趨同趨勢(shì)。東部地區(qū)林業(yè)生態(tài)效率水平基本上都保持較高的水平，北京因?yàn)槠涮厥獾膮^(qū)位優(yōu)勢(shì)、政治經(jīng)濟(jì)中心，使得其也擁有較高的林業(yè)生態(tài)效率水平；河北省是個(gè)特例，承接京津地區(qū)高能耗和高污染的產(chǎn)業(yè)轉(zhuǎn)移的首選陣地，導(dǎo)致其林業(yè)生態(tài)效率較低。中部地區(qū)林業(yè)生態(tài)效率除了山西、湖南基本上均值在0.5 以上。山西省2016年森林覆蓋率為20.5%，林業(yè)資源相對(duì)較少，主要以煤炭能源生產(chǎn)為主，對(duì)于林業(yè)資源開發(fā)利用相對(duì)來說不那么重視，導(dǎo)致生產(chǎn)技術(shù)水平落后，從而導(dǎo)致林業(yè)生態(tài)效率偏低；湖南省林業(yè)資源豐富，以資源換經(jīng)濟(jì)發(fā)展的模式導(dǎo)致前期的林業(yè)生態(tài)效率水平較低，而后在承接產(chǎn)業(yè)轉(zhuǎn)型的過程中，積極調(diào)整戰(zhàn)略結(jié)構(gòu)進(jìn)行產(chǎn)業(yè)升級(jí)改造，2012年以后林業(yè)生態(tài)效率達(dá)到0.5 以上。西部地區(qū)整體林業(yè)生態(tài)效率值偏低，而重慶由于其在西部特殊地理位置以及經(jīng)濟(jì)水平，故其林業(yè)生態(tài)效率相比其他西部地區(qū)高很多。

3.2 我國(guó)林業(yè)生態(tài)效率的時(shí)空演變趨勢(shì)

為了探索各?。▍^(qū)、市）林業(yè)生態(tài)效率時(shí)空演變情況，選擇2005年、2008年、2011年、2014年和2016年為觀測(cè)點(diǎn)，運(yùn)用非參數(shù)的Kernel密度函數(shù)得出其不同時(shí)期的分別差異（圖1）。具體來看，相對(duì)于2005年，2008年密度函數(shù)中心明顯向左移動(dòng)，峰值明顯增大，變化區(qū)間縮小，這說明在此期間整體生態(tài)效率顯著降低，同時(shí)兩極分化現(xiàn)象有所改善，地區(qū)差距縮小。與2008年相比，2011年、2014年密度函數(shù)中心向右移動(dòng)，說明在此期間整體生態(tài)效率有所提高。2016年密度函數(shù)相比2014年峰值明顯下降，這表明在此期間整體生態(tài)效率繼續(xù)改善，同時(shí)地區(qū)差距有所縮小。觀察2005—2016 東部、東北部、中部、西部以及全國(guó)的2005—2016年林業(yè)生態(tài)效率的時(shí)空演變趨勢(shì)，結(jié)果如圖2所示。由圖2可以看出，西部、中部以及全國(guó)的林業(yè)生態(tài)效率都是呈現(xiàn)先下降后上升再波動(dòng)后穩(wěn)定的趨勢(shì)。2005—2009 西部、中部及中國(guó)都是年呈下降趨勢(shì)，而且2009年都是整個(gè)研究期的最低點(diǎn)，而后從2010年開始上升，到2012年后就開始波動(dòng)以致保持穩(wěn)定趨勢(shì)。2005—2010 東部地區(qū)年林業(yè)生態(tài)效率值在0.52～0.55 的范圍波動(dòng)，2010年后就直線上升，2012年后就保持較為穩(wěn)定的趨勢(shì)，數(shù)值接近0.7。東北部地區(qū)林業(yè)生態(tài)效率開始呈現(xiàn)下降的趨勢(shì)，也是在2009年到達(dá)最低點(diǎn)，而后開始上升，但是2012年以后又呈下降的趨勢(shì)。2009年西部和東北部林業(yè)生態(tài)效率數(shù)值最低，低于0.2，全國(guó)林業(yè)生態(tài)效率值高于中部地區(qū)，東部地區(qū)最高。其中，2010年以后中部地區(qū)林業(yè)生態(tài)效率值均高于全國(guó)林業(yè)生態(tài)效率值；2012年中部地區(qū)林業(yè)生態(tài)效率值呈接近東部地區(qū)數(shù)值；2014年以前東北部地區(qū)林業(yè)生態(tài)效率值高于西部地區(qū)，而2014年以后均低于西部地區(qū)；東部地區(qū)林業(yè)生態(tài)效率一直保持較為明顯的領(lǐng)先優(yōu)勢(shì)。

將圖1～2的圖形結(jié)合起來看，各年份的核密度函數(shù)圖與各區(qū)域的林業(yè)生態(tài)效率變化趨勢(shì)有一致性。開始林業(yè)生態(tài)效率整體呈下降趨勢(shì)，而后保持上升態(tài)勢(shì)，兩者直接可以相互輔證變化趨勢(shì)。在2009年前，各地區(qū)都是以經(jīng)濟(jì)發(fā)展為主要目標(biāo)，以資源過度消耗為代價(jià)發(fā)展經(jīng)濟(jì)，與此同時(shí)對(duì)于生態(tài)環(huán)境缺乏相關(guān)意識(shí)和保護(hù)，雖然經(jīng)濟(jì)發(fā)展了但是給生態(tài)環(huán)境造成了嚴(yán)重后果，如2008年的特大洪災(zāi)、冰凍災(zāi)害還讓人記憶猶新。而且由于技術(shù)水平低，資源過度消耗，環(huán)境不斷惡化，從而導(dǎo)致林業(yè)生態(tài)效率持續(xù)走低。自2008年提出的集體林權(quán)制度改革，生態(tài)文明理念的提出和不斷推進(jìn)，使得各地區(qū)逐漸重視生態(tài)環(huán)境，積極引進(jìn)先進(jìn)的生產(chǎn)技術(shù)和提高生產(chǎn)效率以及降低污染能耗，減輕生態(tài)環(huán)境的壓力。東部地區(qū)利用其地域優(yōu)勢(shì)，經(jīng)濟(jì)較為發(fā)達(dá)，生產(chǎn)條件、技術(shù)水平相比其他地區(qū)優(yōu)勢(shì)更為明顯，調(diào)整了發(fā)展模式由粗放型向集約型的轉(zhuǎn)變，因此2012年后可以長(zhǎng)期保持較為穩(wěn)定的態(tài)勢(shì)。中部地區(qū)在承接?xùn)|部地區(qū)的產(chǎn)業(yè)轉(zhuǎn)移過程中，引進(jìn)東部地區(qū)的先進(jìn)工藝和產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化升級(jí)，但由于還處于轉(zhuǎn)型期其尚未完全擺脫傳統(tǒng)模式的影響，故林業(yè)生態(tài)效率值出現(xiàn)一定程度的下降。西部地區(qū)由于地域因素造成的交通不便，經(jīng)濟(jì)較為落后，導(dǎo)致仍以資源的過度消耗為主發(fā)展經(jīng)濟(jì)，生產(chǎn)技術(shù)水平較為低下，從而導(dǎo)致林業(yè)生態(tài)效率水平較低。

圖1 我國(guó)林業(yè)生態(tài)效率的Kernel 密度估計(jì)Fig.1 Kernel density estimation of China’s forestry ecological efficiency

圖2 我國(guó)林業(yè)生態(tài)效率的時(shí)空演變趨勢(shì)Fig.2 Temporal and spatial evolution trend of China’s forestry ecological efficiency

3.3 門檻回歸結(jié)果

本文利用上節(jié)測(cè)度出來的2005—2016年我國(guó)林業(yè)生態(tài)效率數(shù)據(jù)為被解釋變量，環(huán)境規(guī)制為核心解釋變量，借鑒Hansen（1999）的“自助法”（bootstrap），利用STATA12.0 軟件，通過反復(fù)抽樣1 000 次測(cè)算出P值，以判斷是否存在門檻值、存在幾門檻值以及門檻效應(yīng)如何，檢驗(yàn)結(jié)果見表3。

由表3可知，F(xiàn)統(tǒng)計(jì)量無論是在單門檻、雙門檻還是三門檻模型中，都通過10%的顯著性檢驗(yàn)，因此說明環(huán)境規(guī)制對(duì)林業(yè)生態(tài)效率的影響存在基于環(huán)境規(guī)制為門檻變量的三門檻效應(yīng)。因此，本文對(duì)該三門檻模型進(jìn)行進(jìn)一步檢驗(yàn)和估計(jì)（表4）。

設(shè)定環(huán)境規(guī)制為門檻變量，數(shù)值的不同對(duì)林業(yè)生態(tài)效率影響有較大差異。當(dāng)環(huán)境規(guī)制條件寬松時(shí)（ER ≤0.599 5），對(duì)林業(yè)生態(tài)效率有正向影響，影響系數(shù)為0.324 1。因?yàn)榄h(huán)境規(guī)制條件寬松時(shí)，對(duì)于污染排放要求較低，此時(shí)企業(yè)利潤(rùn)最大化的產(chǎn)量要高于最優(yōu)化產(chǎn)量，促進(jìn)經(jīng)濟(jì)效益的增長(zhǎng)。但由于污染排放對(duì)于環(huán)境具有典型的負(fù)外部效應(yīng)，導(dǎo)致無法實(shí)現(xiàn)資源的最優(yōu)化配置，產(chǎn)生資源過度消耗以及污染過度排放等現(xiàn)象。當(dāng)環(huán)境規(guī)制條件較為寬松時(shí)（0.599 5＜ER ≤0.709 7），對(duì)林業(yè)生態(tài)效率正向影響，影響系數(shù)為0.213 5；當(dāng)環(huán)境規(guī)制條件較為嚴(yán)格時(shí)（0.709 7＜ER ≤0.816 5），對(duì)林業(yè)生態(tài)效率正向影響，影響系數(shù)為0.149 3。由于在0.599 5 和0.816 5 區(qū)間水平里，企業(yè)會(huì)根據(jù)環(huán)境規(guī)制要求在允許的條件下進(jìn)行技術(shù)改造升級(jí)，以適應(yīng)產(chǎn)業(yè)發(fā)展的要求，故這是環(huán)境規(guī)制對(duì)林業(yè)生態(tài)效率是促進(jìn)作用，但是隨著環(huán)境規(guī)制的強(qiáng)度越高，其產(chǎn)生的成本越來越高，導(dǎo)致影響系數(shù)的降低。當(dāng)環(huán)境規(guī)制條件嚴(yán)格時(shí)（ER＞0.816 5），對(duì)林業(yè)生態(tài)效率負(fù)向影響，影響系數(shù)為-0.05 4。因?yàn)殡S著環(huán)境規(guī)制條件越來越嚴(yán)苛，其對(duì)于林業(yè)生態(tài)效率影響為有利變?yōu)椴焕?。環(huán)境規(guī)制越嚴(yán)苛，由此會(huì)增加企業(yè)的生產(chǎn)成本，導(dǎo)致企業(yè)會(huì)調(diào)整產(chǎn)量以適應(yīng)環(huán)境規(guī)制要求，甚至有企業(yè)為此而停產(chǎn)，由此環(huán)境規(guī)制對(duì)生態(tài)效率產(chǎn)生負(fù)面影響，這也是為什么在環(huán)境規(guī)制水平大于0.816 5 時(shí)，會(huì)導(dǎo)致影響為負(fù)的原因。

表3 門檻效應(yīng)顯著性檢驗(yàn)結(jié)果?Table 3 Significant test results of threshold effect

表4 三門檻模型估計(jì)結(jié)果?Table 4 Estimation results of three threshold models

4 結(jié)論與建議

本文運(yùn)用非徑向方法距離函數(shù)求得綠色發(fā)展理念的我國(guó)26 個(gè)?。▍^(qū)、市）2005—2016年的林業(yè)生態(tài)效率值，然后利用非參數(shù)的Kernel 密度函數(shù)分析林業(yè)生態(tài)效率的時(shí)空演變趨勢(shì)，再利用面板門檻回歸模型分析環(huán)境規(guī)制強(qiáng)度對(duì)林業(yè)生態(tài)效率的影響程度。

1）我國(guó)林業(yè)生態(tài)效率整體上呈現(xiàn)從西部到東南部慢慢升高的趨勢(shì)，以東部地區(qū)的林業(yè)生態(tài)效率最高，中部地區(qū)林業(yè)生態(tài)效率值居中；西部、東北部林業(yè)生態(tài)效率值較低。因此，我國(guó)各省（區(qū)、市）的林業(yè)生態(tài)效率值分別與各省份經(jīng)濟(jì)實(shí)力保持基本趨同趨勢(shì)。其中，林業(yè)生態(tài)效率最高值為廣東省，最低為陜西省。

2）分析林業(yè)生態(tài)效率整體時(shí)空演變，2005—2009年間密度函數(shù)中心先向左移動(dòng)，表明整體生態(tài)效率顯著降低，2009年后密度函數(shù)中心向右移動(dòng)，表明整體生態(tài)效率有所提升，2016年密度函數(shù)相比2014年峰值明顯下降，表明整體生態(tài)效率繼續(xù)改善，同時(shí)地區(qū)差距有所縮小；分區(qū)域分析林業(yè)生態(tài)效率時(shí)空演變，西部、中部以及全國(guó)的林業(yè)生態(tài)效率都是呈現(xiàn)先下降而后上升再波動(dòng)后穩(wěn)定的趨勢(shì)，東部地區(qū)林業(yè)生態(tài)效率值先固定范圍波動(dòng)，而后直線上升，再保持較為穩(wěn)定的趨勢(shì)，東北部地區(qū)林業(yè)生態(tài)效率先下降再上升而后下降的趨勢(shì)。

3）環(huán)境規(guī)制對(duì)于林業(yè)生態(tài)效率影響表明：當(dāng)環(huán)境規(guī)制條件小于0.599 5 時(shí)，對(duì)林業(yè)生態(tài)效率正向影響；當(dāng)環(huán)境規(guī)制條件在0.599 5 和0.816 5 區(qū)間水平里，環(huán)境規(guī)制對(duì)林業(yè)生態(tài)效率是促進(jìn)作用，但是隨著環(huán)境規(guī)制的強(qiáng)度越高，其產(chǎn)生的成本越來越高，導(dǎo)致影響系數(shù)的降低。環(huán)境規(guī)制條件大于0.816 5 時(shí)，對(duì)生態(tài)效率產(chǎn)生負(fù)面影響。

基于以上實(shí)證分析提出以下4 點(diǎn)建議：

1）完善政策引導(dǎo)扶持方向，大力發(fā)展林業(yè)產(chǎn)業(yè)。開辟產(chǎn)業(yè)融合發(fā)展新路徑，實(shí)現(xiàn)林業(yè)從單純的生產(chǎn)功能向生態(tài)、文化功能的全面拓展，特色林業(yè)與傳統(tǒng)文化、現(xiàn)代科技、文化旅游、創(chuàng)意體驗(yàn)、互聯(lián)網(wǎng)以及金融等有機(jī)融合，推動(dòng)區(qū)域產(chǎn)業(yè)升級(jí)。由此可見，應(yīng)在做強(qiáng)做大主導(dǎo)產(chǎn)業(yè)的同時(shí)，加大研發(fā)力度，挖掘林產(chǎn)品潛在價(jià)值、生態(tài)價(jià)值、文化價(jià)值，提升價(jià)值鏈。促進(jìn)林業(yè)與旅游、文化、教育、健康養(yǎng)老等深度融合發(fā)展，延伸產(chǎn)業(yè)鏈，形成種養(yǎng)、加工、科技、營(yíng)銷等全產(chǎn)業(yè)鏈格局，實(shí)現(xiàn)三大產(chǎn)業(yè)融合發(fā)展。

2）大力發(fā)展綠色制造業(yè)。一是加快傳統(tǒng)產(chǎn)業(yè)轉(zhuǎn)型升級(jí)。按照高端化、綠色化、集約化要求，繼續(xù)淘汰落后產(chǎn)能，推動(dòng)傳統(tǒng)林業(yè)產(chǎn)業(yè)采用先進(jìn)技術(shù)，不斷優(yōu)化產(chǎn)品結(jié)構(gòu)，逐步向低能耗、綠色環(huán)保、高附加值的新領(lǐng)域發(fā)展。二是運(yùn)用大數(shù)據(jù)提供優(yōu)質(zhì)生態(tài)產(chǎn)品。積極發(fā)展生態(tài)產(chǎn)品全產(chǎn)業(yè)鏈的數(shù)據(jù)存儲(chǔ)、采集、處理、教育、安全等核心業(yè)務(wù)，培育生態(tài)相關(guān)產(chǎn)業(yè)智能終端產(chǎn)品及服務(wù)的相關(guān)業(yè)態(tài)，不斷豐富生態(tài)相關(guān)產(chǎn)業(yè)中的智能制造、智慧健康、旅游、物流等衍生業(yè)態(tài)發(fā)展，實(shí)現(xiàn)大數(shù)據(jù)生態(tài)產(chǎn)業(yè)的快速增長(zhǎng)。三是鼓勵(lì)金融機(jī)構(gòu)進(jìn)一步加大對(duì)林業(yè)發(fā)展的支持力度，創(chuàng)新林業(yè)金融產(chǎn)品和服務(wù)，開發(fā)林業(yè)信貸產(chǎn)品，支持技術(shù)改造、綠色林產(chǎn)品加工、循環(huán)產(chǎn)業(yè)、清潔能源、新材料和大數(shù)據(jù)等產(chǎn)業(yè)，特別是支持實(shí)體經(jīng)濟(jì)綠色發(fā)展。如林權(quán)抵押+政府增信、林業(yè)APP、龍頭企業(yè)+基地+新型經(jīng)營(yíng)主體+林農(nóng)等融資模式，推廣以自身現(xiàn)金流還款的模型，多渠道為林業(yè)企業(yè)或林農(nóng)建立融資模式。

3）發(fā)展新旅游業(yè)態(tài)，提高旅游經(jīng)濟(jì)效益。發(fā)揮林業(yè)生態(tài)資源優(yōu)勢(shì)，特別是地域特色、國(guó)家森林公園、自然保護(hù)區(qū)等有著自身特點(diǎn)的旅游景區(qū)，開發(fā)多樣式的旅游模式如自駕車營(yíng)地、房車露營(yíng)、戶外運(yùn)動(dòng)、低空飛行等新業(yè)態(tài)旅游產(chǎn)品，以適應(yīng)人們的旅游需求。

4）利用大數(shù)據(jù)更好地保護(hù)生態(tài)環(huán)境。推動(dòng)生態(tài)環(huán)境保護(hù)大數(shù)據(jù)應(yīng)用，運(yùn)用遙感、大數(shù)據(jù)、云計(jì)算等技術(shù)，更加精準(zhǔn)有效地推動(dòng)生態(tài)建設(shè)和污染防治，提高自然生態(tài)修復(fù)能力。以提升環(huán)境管理信息化和智能化水平為目標(biāo)，建立以大數(shù)據(jù)為技術(shù)手段的環(huán)境在線監(jiān)控體系，加大行業(yè)重點(diǎn)污染源自動(dòng)監(jiān)控設(shè)施建設(shè)力度。建設(shè)涵蓋大氣、水、土壤、噪音等要素的環(huán)境質(zhì)量監(jiān)測(cè)網(wǎng)絡(luò)，完善生態(tài)環(huán)境監(jiān)測(cè)數(shù)據(jù)集成共享機(jī)制，推動(dòng)環(huán)境質(zhì)量、重點(diǎn)污染源、生態(tài)狀況大數(shù)據(jù)監(jiān)測(cè)全覆蓋，實(shí)現(xiàn)多部門數(shù)據(jù)融合集成，互通共享，推動(dòng)生態(tài)環(huán)境治理體系和治理能力現(xiàn)代化。