陳 欽，黃巧龍，游玲娜，謝麗星，程秋旺，潘 輝

（1.福建農(nóng)林大學(xué) 經(jīng)濟(jì)學(xué)院，福州 350002；2福建農(nóng)林大學(xué) 管理學(xué)院，福州 350002；3.閩江學(xué)院，福州 350108）

隨著福建省鐵路、公路建設(shè)的不斷擴(kuò)大，水源保護(hù)區(qū)域的增加，生態(tài)保護(hù)的加強(qiáng)，有些商品林被政府劃為生態(tài)公益林，有些商品林被政府界定為重點(diǎn)生態(tài)區(qū)位商品林，禁止采伐。重點(diǎn)生態(tài)區(qū)位商品林雖然形式上屬于商品林，但是經(jīng)濟(jì)實(shí)質(zhì)上雷同于生態(tài)公益林。這會(huì)對(duì)林農(nóng)產(chǎn)生經(jīng)濟(jì)影響。本文從林農(nóng)的微觀角度出發(fā)，分析福建省生態(tài)公益林保護(hù)對(duì)林農(nóng)的經(jīng)濟(jì)影響。從理論研究來看，針對(duì)生態(tài)公益林保護(hù)對(duì)林農(nóng)經(jīng)濟(jì)影響的研究較少，未見報(bào)道采用大樣本和DID模型進(jìn)行研究；從實(shí)踐方面來看，本文在問卷調(diào)查的基礎(chǔ)上分析福建省生態(tài)公益林保護(hù)對(duì)林農(nóng)的經(jīng)濟(jì)影響以及影響金額，對(duì)政府制定福建省生態(tài)公益林生態(tài)補(bǔ)償價(jià)格，保障林農(nóng)的切身利益，促進(jìn)林業(yè)生態(tài)效益和經(jīng)濟(jì)效益的提高，進(jìn)而推進(jìn)生態(tài)公益林保護(hù)，具有重要的實(shí)踐意義。

1 模型設(shè)定

本文使用DID方法來評(píng)估生態(tài)公益林保護(hù)對(duì)林農(nóng)的經(jīng)濟(jì)影響，用林農(nóng)家庭林業(yè)收入作為經(jīng)濟(jì)影響的衡量指標(biāo)。DID模型的基本思想是把樣本分為控制組（Control）與處理組（Treated）(史桂芬,2012)[1]，其中，控制組為未受到生態(tài)公益林保護(hù)政策影響的樣本，處理組為受到生態(tài)公益林保護(hù)政策影響的樣本，將兩組進(jìn)行比較從而評(píng)估該政策對(duì)林農(nóng)林業(yè)收入的影響，首先收集兩個(gè)時(shí)期的數(shù)據(jù)，即2002年和2013年相關(guān)變量數(shù)據(jù)；其次將樣本分為兩組，2002年只有商品林，2013年擁有生態(tài)公益林的樣本為處理組；而2002年和2013年都只有商品林，沒有生態(tài)公益林的樣本為控制組。通過實(shí)證分析，分別得到處理組和控制組在2002年和2013年林業(yè)收入的變化值，這兩個(gè)變化值的差額就反映了生態(tài)公益林保護(hù)對(duì)林農(nóng)林業(yè)收入的凈影響(古曉,2013)[2]。

在沒有考慮其他控制變量的情況下，DID模型的表達(dá)式如下：

其中，i代表林農(nóng)，t代表時(shí)期，在t=0階段，所調(diào)查的林農(nóng)的森林是商品林；在t=1階段，上述一部分林農(nóng)的森林由商品林變成生態(tài)公益林，而另一部分林農(nóng)的森林仍然是商品林。關(guān)鍵解釋變量FPCit是一組二值虛擬變量，代表t時(shí)期i林農(nóng)情況，分別測量林農(nóng)的森林是否由商品林變成生態(tài)公益林，1表示是，0表示否。Yit是t時(shí)期i林農(nóng)的收入；t是時(shí)期虛擬變量，生態(tài)公益林保護(hù)實(shí)施前取值為0，實(shí)施后取值為1。ci是非觀測效應(yīng)，控制隨時(shí)間不變的不可觀測因素；μit是特異性擾動(dòng)項(xiàng)，代表因林農(nóng)因時(shí)而變且影響Yit的那些非觀測擾動(dòng)因素。α0、α1和δ是待估計(jì)參數(shù)，其中δ是研究最為關(guān)注的參數(shù)，衡量了生態(tài)公益林保護(hù)對(duì)變量Yit的經(jīng)濟(jì)影響。

根據(jù)該模型，可以分別得到處理組和控制組在生態(tài)公益林保護(hù)實(shí)施前后林農(nóng)林業(yè)收入變動(dòng)的模型：當(dāng)FPCit=0，t=0時(shí)，得到控制組和處理組在生態(tài)公益林保實(shí)施前林業(yè)收入的變動(dòng)模型；當(dāng)FPCit=0，t=1時(shí)，得到控制組在生態(tài)公益林保實(shí)施后林業(yè)收入的變動(dòng)模型；當(dāng)FPCit=1，t=1時(shí)，得到處理組在生態(tài)公益林保實(shí)施后林業(yè)收入的變動(dòng)模型，具體如下所示：

根據(jù)式（2），可得到控制組在生態(tài)公益林保護(hù)實(shí)施前后林農(nóng)林業(yè)收入的變化值為：

根據(jù)式（3），可得到處理組在生態(tài)公益林保護(hù)實(shí)施前后林農(nóng)林業(yè)收入的變化值為：

根據(jù)式（4）和（5），生態(tài)公益林保護(hù)對(duì)林農(nóng)林業(yè)收入的凈影響為：

為了剔除其它因素對(duì)林業(yè)收入的影響，本文進(jìn)一步引入控制變量，采用非觀測效應(yīng)綜列數(shù)據(jù)模型如下：

式中，Zit是一組隨時(shí)間變化的可觀測的影響因變量Yit的控制變量；Xit是一組不隨時(shí)間變化，或隨時(shí)間同等變化的可觀測的影響因變量Yit的控制變量；β和γ是待估參數(shù)矩陣，衡量了控制變量對(duì)因變量Yit的作用。

2 理論分析

經(jīng)濟(jì)學(xué)認(rèn)為在經(jīng)濟(jì)活動(dòng)中，人們都是理性的，人們參加經(jīng)濟(jì)活動(dòng)時(shí)都持有預(yù)期收益大于預(yù)期成本的觀點(diǎn)，這是人們參與活動(dòng)的基本出發(fā)點(diǎn)和動(dòng)力。人們投入的主要目的是獲得收益，只有當(dāng)其收益大于成本時(shí)，才是經(jīng)濟(jì)合理的。一般而言，農(nóng)戶作為理性的經(jīng)濟(jì)主體，其生態(tài)公益林建設(shè)和保護(hù)是為了追求潛在的收益，即通過比較其成本和收益水平，在可接受的范圍內(nèi)，產(chǎn)生行為。由此可見，在生態(tài)公益林建設(shè)和保護(hù)過程中，成本和收益水平是農(nóng)戶決策的核心與關(guān)鍵。

2.1 生態(tài)公益林保護(hù)對(duì)林農(nóng)林業(yè)收入的影響

生態(tài)公益林保護(hù)對(duì)林農(nóng)林業(yè)收入的影響主要體現(xiàn)在商品林和生態(tài)公益林的收益差距上。為了保護(hù)森林資源，提高生態(tài)環(huán)境質(zhì)量，政府采用自上而下的辦法將部分商品林強(qiáng)制劃定為生態(tài)公益林，禁止對(duì)其進(jìn)行采伐，導(dǎo)致林業(yè)收入減少(沈潔,2014)[3]。為此，政府制定了生態(tài)補(bǔ)償標(biāo)準(zhǔn)，生態(tài)公益林的經(jīng)營者可以獲得一定的經(jīng)濟(jì)補(bǔ)償，但是現(xiàn)行的生態(tài)補(bǔ)償標(biāo)準(zhǔn)低，對(duì)于福建省絕大部分生態(tài)公益林而言明顯不夠，使得很多林農(nóng)獲得的補(bǔ)償價(jià)不能彌補(bǔ)營林成本。

2.2 收入影響因素的研究綜述

農(nóng)業(yè)收入會(huì)受到農(nóng)業(yè)生產(chǎn)條件、生產(chǎn)經(jīng)營費(fèi)用、勞動(dòng)力等因素的影響；同時(shí)，根據(jù)公共財(cái)政理論，財(cái)政支農(nóng)也會(huì)影響農(nóng)業(yè)收入(楊靜,2014)[4]。高陽、趙正等(2014)運(yùn)用實(shí)證方法分析林農(nóng)林業(yè)收入的影響因素，結(jié)果表明：林農(nóng)林業(yè)收入的影響因素包括受教育程度、林下經(jīng)營項(xiàng)目、林種等[5]。林和平(2009)分別分析了林地資源數(shù)量、林地資源經(jīng)營方式與林業(yè)收入之間的關(guān)系，得出：耕地資源、林地資源、人力資源以及樹種等都對(duì)林業(yè)收入有顯著的影響[6]。劉偉平(2009)認(rèn)為，經(jīng)濟(jì)環(huán)境和價(jià)格等經(jīng)濟(jì)因素、勞動(dòng)力數(shù)量、林業(yè)政策等都會(huì)影響林農(nóng)的林業(yè)收入[7]。石康橋、蘇時(shí)鵬等(2014)通過比較集體林權(quán)制度改革前后林農(nóng)林業(yè)收入的影響因素，得出：林改前林農(nóng)林業(yè)收入的主要影響因素包括林地投入、經(jīng)營類型等；而林改后林農(nóng)林業(yè)收入的主要影響因素為受教育程度、林地面積等[8]。姚林香、舒成(2010)通過對(duì)江西7縣農(nóng)戶樣本數(shù)據(jù)的分析，發(fā)現(xiàn)：農(nóng)業(yè)人口數(shù)、耕地面積等對(duì)農(nóng)民收入具有顯著影響[9]。

2.3 變量的選擇及說明

本文關(guān)注的是生態(tài)公益林保護(hù)對(duì)林農(nóng)的經(jīng)濟(jì)影響，選取林業(yè)收入作為模型的因變量。模型選取的自變量包括虛擬變量、隨時(shí)間變化的可觀測的影響林業(yè)收入的控制變量以及不隨時(shí)間變化，或隨時(shí)間同等變化的可觀測的影響林業(yè)收入的控制變量。綜合有關(guān)學(xué)者的研究成果以及理論分析，首先，本文選取時(shí)期、商品林是否變成生態(tài)公益林作為虛擬變量。如果由于高速公路、高鐵、水庫等建設(shè)，林農(nóng)的商品林被劃為生態(tài)公益林，禁止采伐，林農(nóng)就沒有木材采伐收入，為此林農(nóng)可以獲得森林生態(tài)效益補(bǔ)償收入、生態(tài)公益林林下經(jīng)營補(bǔ)貼收入等，可見，商品林是否變成生態(tài)公益林對(duì)林農(nóng)林業(yè)收入產(chǎn)生影響，本文試圖研究這個(gè)關(guān)鍵變量的影響大小。其次，選取家庭勞動(dòng)力、人均GDP、職務(wù)、受教育程度、耕地面積、林地面積和樹種作為控制變量，具體說明如下：①家庭勞動(dòng)力。家庭勞動(dòng)力越多，林業(yè)生產(chǎn)要素中的人力資源數(shù)量越多，所以林業(yè)收入越多，假設(shè)家庭勞動(dòng)力與林業(yè)收入正相關(guān)；②人均GDP。人均GDP越高，表明當(dāng)?shù)亟?jīng)濟(jì)發(fā)展水平越高，木材和林產(chǎn)品的銷售量越大，所以林農(nóng)的林業(yè)收入越多，假設(shè)人均GDP與林業(yè)收入正相關(guān)；③林農(nóng)職務(wù)。擔(dān)任職務(wù)（如村干部）的林農(nóng)，一般社會(huì)資本更多，更容易獲得木材采伐指標(biāo)和更快了解到林業(yè)信息，林業(yè)收入也會(huì)更高，假設(shè)林農(nóng)擔(dān)任的職務(wù)與林業(yè)收入正相關(guān)；④受教育程度。受教育程度越高，有利于林農(nóng)及時(shí)地對(duì)政策和市場作出反應(yīng)，林農(nóng)林業(yè)經(jīng)營能力越強(qiáng)，林業(yè)收入越高，假設(shè)受教育程度與林業(yè)收入正相關(guān)；⑤耕地面積。耕地面積越大，對(duì)耕地投入越多，農(nóng)業(yè)收入越多，對(duì)林業(yè)的經(jīng)濟(jì)依賴越小，對(duì)林業(yè)的投入隨之減少，導(dǎo)致林業(yè)收入減少，假設(shè)耕地面積與林業(yè)收入負(fù)相關(guān)；⑥家庭林地面積，一般情況下林地面積越多的林農(nóng)，林業(yè)收入越高，假設(shè)林地面積與林業(yè)收入正相關(guān)；⑦樹種。不同樹種的經(jīng)濟(jì)價(jià)值不一樣，一個(gè)樹種經(jīng)濟(jì)價(jià)值越大，其經(jīng)濟(jì)收益越高，假設(shè)樹種的經(jīng)濟(jì)價(jià)值與林業(yè)收入正相關(guān)。變量的選擇及其賦值說明具體如表1所示。

3 計(jì)量結(jié)果及其分析

考慮到可能存在影響林業(yè)收入的其他因素，為了剔除這些因素對(duì)林業(yè)收入的影響，本文選擇人均GDP、家庭勞動(dòng)力、林農(nóng)職務(wù)、受教育年限、家庭耕地面積、林地面積、樹種作為控制變量，分析福建省生態(tài)公益林保護(hù)對(duì)林農(nóng)林業(yè)收入的凈影響。運(yùn)用Stata12.0軟件，將上述控制變量放入模型中進(jìn)行估計(jì)，得出引入控制變量后DID模型的計(jì)量結(jié)果，如表2和表3所示，模型的Prob＞F=0.000 0，說明模型整體顯著；R方為0.631 7，調(diào)整后R方為0.624 0，說明模型擬合度較好。

表2 DID模型計(jì)量結(jié)果一

從表2可以看出，DID值的系數(shù)為-8 700，且在10%的統(tǒng)計(jì)水平上顯著，說明在控制其他因素的情況下，福建省生態(tài)公益林保護(hù)使樣本林農(nóng)戶均林業(yè)收入降低了8 700元，也就是說福建省生態(tài)公益林保護(hù)對(duì)林農(nóng)的經(jīng)濟(jì)影響為-8 700元。政府確定生態(tài)公益林補(bǔ)償價(jià)格可以參考這個(gè)數(shù)字，并且據(jù)此計(jì)算每年每畝生態(tài)公益林補(bǔ)償標(biāo)準(zhǔn)，從而為政府修訂生態(tài)公益林補(bǔ)償政策提供重要依據(jù)?？梢?，政府應(yīng)重視生態(tài)公益林補(bǔ)償政策的完善和實(shí)施，彌補(bǔ)生態(tài)公益林保護(hù)給林農(nóng)帶來的經(jīng)濟(jì)損失，從而促進(jìn)生態(tài)公益林保護(hù)工作的有效進(jìn)行。

從表3可知，時(shí)期P值為0.000＜0.01，說明2013年與2002年相比，樣本林農(nóng)林業(yè)收入發(fā)生了顯著變化；人均GDP對(duì)林業(yè)收入影響的回歸系數(shù)為正，且P值為0.002＜0.01，說明人均GDP對(duì)林農(nóng)林業(yè)收入在1%水平上具有正向的顯著影響；家庭耕地面積的回歸系數(shù)為負(fù)，且P值為0.071＜0.1，說明耕地面積對(duì)林業(yè)收入在10%水平上具有負(fù)向的顯著影響；家庭林地面積的回歸系數(shù)為正，且P值為0.000＜0.01，說明林地面積對(duì)林業(yè)收入在1%水平上具有正向的顯著影響。

表3 DID模型計(jì)量結(jié)果二