孫良斌 謝枰飛 蔣桂容

[摘要]運用2006～2015年中國省際面板數(shù)據(jù)實證檢驗了中國農(nóng)業(yè)面源污的排放總量和排放總強度的經(jīng)濟機制，實證結(jié)果表明，中國的農(nóng)村經(jīng)濟總量、林業(yè)經(jīng)濟總量、農(nóng)村人口教育程度和政府防治政策投入對農(nóng)業(yè)面源污染的排放總量和排放總強度有顯著的負效應(yīng)，而農(nóng)業(yè)生產(chǎn)總量，農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)、特別是漁業(yè)和畜牧業(yè)生產(chǎn)總量對排放總量和排放總強度有顯著的正影響。這表明，中國要實施綠色發(fā)展的鄉(xiāng)村振興戰(zhàn)略，需強化農(nóng)村經(jīng)濟發(fā)展，特別是農(nóng)業(yè)面源污染嚴重的貧困地區(qū)的經(jīng)濟發(fā)展，使農(nóng)戶內(nèi)生出防治污染的需求動力；要把漁業(yè)和畜牧業(yè)發(fā)達的地區(qū)作為農(nóng)業(yè)面源污染防治的重點區(qū)域；加強農(nóng)村人口教育，特別是對農(nóng)業(yè)面源污染的認知教育，繼續(xù)加大農(nóng)業(yè)面源污染防治的政策投入，為農(nóng)業(yè)面源污染的攻堅克難和綠色發(fā)展鄉(xiāng)村振興戰(zhàn)略提供政策保障。

[關(guān)鍵詞]農(nóng)業(yè)面源污染；經(jīng)濟機制；面板數(shù)據(jù)

[中圖分類號]F323.22 [文獻標識碼]A

農(nóng)業(yè)面源污染是中國水體和土壤污染最重要的污染源，已經(jīng)成為水體污染最主要的原因（張維理等，2004）。從20世紀70年代以來，中國國家就開始對農(nóng)田水污染開始治理，當前，中國的農(nóng)業(yè)面源正處于治理的攻堅階段（金書秦，邢曉旭，2018），隨著中國實施生態(tài)文明建設(shè)、綠色發(fā)展戰(zhàn)略和鄉(xiāng)村振興戰(zhàn)略，農(nóng)業(yè)面源污染的治理尤其凸顯其重要性和緊迫性，國家和政府相關(guān)部門相繼出臺了一系列的政策措施：2015年農(nóng)業(yè)部發(fā)布《關(guān)于打好農(nóng)業(yè)面源污染防治攻堅戰(zhàn)的實施意見》 批準實施，標志著中國對農(nóng)業(yè)面源污染防控治理攻堅戰(zhàn)正式打響；2017年3月農(nóng)業(yè)部印發(fā) 《2017年農(nóng)業(yè)面源污染防治攻堅戰(zhàn)重點工作安排》，標志著中國對農(nóng)業(yè)面源污染防控治理攻堅戰(zhàn)的細化和持續(xù)化；2018年中央一號文件全面部署實施鄉(xiāng)村振興戰(zhàn)略，將生態(tài)宜居綠色發(fā)展作為鄉(xiāng)村振興發(fā)展的重要內(nèi)容；2018年2月27日，農(nóng)業(yè)部召開農(nóng)業(yè)面源污染防治推進工作組第五次會議，會議強調(diào)，要認真貫徹落實黨的十九大報告精神，圍繞實施鄉(xiāng)村振興戰(zhàn)略，堅決打贏農(nóng)業(yè)面源污染防治攻堅戰(zhàn)，促進農(nóng)業(yè)綠色發(fā)展。

但是，現(xiàn)有文獻研究表明，對農(nóng)業(yè)面源污染形成的經(jīng)濟機制實證文獻欠缺，代表成果有梁流濤，馮淑怡，曲福田（2010）和葛繼紅，周曙東（2011）等，雖較好實證了農(nóng)業(yè)面源的經(jīng)濟形成機制，但尚未考慮到農(nóng)業(yè)塑料地膜的農(nóng)業(yè)污染情況；梁流濤等（2010）較早地從農(nóng)業(yè)面源污染排放總量和強度進行了實證研究，但使用的面板數(shù)據(jù)是2007年前的數(shù)據(jù)，未能體現(xiàn)2007年后中國農(nóng)村和農(nóng)業(yè)轉(zhuǎn)型給農(nóng)業(yè)面源污染的帶來的影響，也未能反映2007年后中央政府農(nóng)業(yè)面源污染治理政策對農(nóng)業(yè)面源污染的影響。因此，本文利用中國2006～2015年31個省的省際面板數(shù)據(jù)從農(nóng)業(yè)面源污染排放總量和排放總強度來探討農(nóng)業(yè)面源形成的經(jīng)濟機制，為中國當前農(nóng)業(yè)面源治理攻堅克難和綠色發(fā)展的鄉(xiāng)村振興戰(zhàn)略提供政策參考的理論依據(jù)，具有重要的現(xiàn)實意義。

1 中國農(nóng)業(yè)面源污染的來源及核算

1.1 中國農(nóng)業(yè)面源污染的主要來源

中國農(nóng)業(yè)面源污染主要來源禽畜牧業(yè)、漁業(yè)，農(nóng)業(yè)塑料地薄，農(nóng)藥化肥施用，以及農(nóng)村生產(chǎn)生活的各種垃圾、廢水等。因此，本文農(nóng)業(yè)面源污染的核算主要從農(nóng)業(yè)化肥農(nóng)藥使用量、農(nóng)業(yè)塑料地膜使用量和農(nóng)業(yè)養(yǎng)殖業(yè)肉豬出欄量來進行核算。數(shù)據(jù)來源于全國以及各省市其他數(shù)據(jù)主要來自歷年的《中國統(tǒng)計年鑒》 、《中國農(nóng)村統(tǒng)計年鑒》 、《中國農(nóng)業(yè)年鑒》 、《環(huán)境統(tǒng)計數(shù)據(jù)》（ 2007～2016年），以及各省市統(tǒng)計年鑒和農(nóng)業(yè)統(tǒng)計資料 。2006～2015年中國農(nóng)業(yè)面源污染源變化情況，見表1。

1.2 中國農(nóng)業(yè)面源污染的核算

由于農(nóng)業(yè)面源污染的污染源分布十分廣泛，具有分散性和隨機性等特點，因此，不同于點源污染，難以通過觀察發(fā)現(xiàn)規(guī)律和有效監(jiān)測，精準核算具有一定的難點，目前中國在農(nóng)業(yè)面源污染核算中采用比較普遍的科學(xué)方法是清單分析法（陳敏鵬，陳吉寧等，2006），清單分析方法是綜合運用多種分析方法來核算的，基本思路是：特定的農(nóng)業(yè)活動對應(yīng)特定的農(nóng)業(yè)污染排放量和排放強度，中國的農(nóng)業(yè)面源污染主要來源農(nóng)村生活垃圾，農(nóng)藥化肥排放量和畜牧養(yǎng)殖等這幾種所造成的污染，其主要的污染排放物是氮（N），磷（P），化學(xué)需氧量（COD）的排放量，因此，根據(jù)特定農(nóng)業(yè)面源污染源的氮（N），磷（P），化學(xué)需氧量（COD）的排放量就可以核算出其農(nóng)業(yè)面源污染的排放總量和排放總強度。本文主要參考計陳敏鵬，陳吉寧等（2014）的清單分析方法進行核算。計算公式如下：

其中，TE為農(nóng)業(yè)面源污染物的排放總量，這里主要是農(nóng)業(yè)面源污染物為N，P，COD的排放總量：氮排放總量（TNq），磷排放總量（TPq），化學(xué)需氧量排放總量（TCODq）；SUi為i污染單元污染物產(chǎn)生系數(shù)，這里面包括農(nóng)林牧漁產(chǎn)量，鄉(xiāng)村人口數(shù)及化肥折純量；ρi為單元i污染物的產(chǎn)污強度系數(shù)；LCi代表的是考慮資源綜合應(yīng)用和管理因素時i污染單元污染物的產(chǎn)污排放系數(shù)；TEl為農(nóng)業(yè)面源污染物的排放總強度： 氮排放總強度（TNs），磷排放總量（TPs），化學(xué)需氧量排放總量（TCODs）；中國農(nóng)業(yè)面源污染物的排放總量和排放強度的核算結(jié)果，見表2。

根據(jù)表2做出圖1和圖2。由圖1和圖2可以看出，2006～2015年間全國農(nóng)業(yè)面源TN，TP和COD的排放量，在整體上仍然呈現(xiàn)緩慢增長的趨勢；在排放強度上，氮和磷的排放強度顯著上升，化學(xué)需氧量的排放強度也先減后增 ，由此可見，農(nóng)業(yè)面源污染雖然經(jīng)過十多年的治理，但形勢依然十分嚴峻。

2 中國農(nóng)業(yè)面源污染形成經(jīng)濟機制的實證分析

2.1 理論分析框架和實證模型構(gòu)建

農(nóng)業(yè)面源污染的形成的經(jīng)濟機理是與農(nóng)村生產(chǎn)生活方式、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)發(fā)展，農(nóng)村人口數(shù)量，農(nóng)村人口對農(nóng)業(yè)生態(tài)環(huán)境的認知程度，以及政府的農(nóng)業(yè)面源污染治理政策等有關(guān)，具體包括以下六個方面。

2.1.1 農(nóng)村經(jīng)濟總量。根據(jù)環(huán)境庫茲涅茨曲線（EKC）（Dinda，S.2004；孟祥海等，2015）可知，人均GDP如果達到一定程度，人們對生活的環(huán)境質(zhì)量有更高需求，對環(huán)境保護有較強意愿，因此農(nóng)村經(jīng)濟總量越高，更注重農(nóng)業(yè)面源污染的防治，對農(nóng)業(yè)面源污染有負效應(yīng)作用。

2.1.2 農(nóng)業(yè)生產(chǎn)總值。一個地區(qū)農(nóng)業(yè)生產(chǎn)總值越高，即農(nóng)業(yè)越發(fā)達，在中國廣大農(nóng)村生產(chǎn)技術(shù)沒有發(fā)生顯著提高的條件下，高投入高產(chǎn)出是常見的辦法，因此農(nóng)業(yè)生產(chǎn)總值越高，越可能造成農(nóng)業(yè)面源污染，我農(nóng)業(yè)面源污染的重點區(qū)域。

2.1.3 農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)。廣義農(nóng)業(yè)包括農(nóng)林牧副漁，從理論上講，林業(yè)的發(fā)展有利于農(nóng)業(yè)面源污染的防治，而畜牧業(yè)和漁業(yè)的發(fā)展更有可能惡化農(nóng)業(yè)面源污染，是重點防治區(qū)域。

2.1.4 農(nóng)村人口規(guī)模數(shù)。農(nóng)村人口越多，產(chǎn)生的垃圾、廢水越多，也越可能造成的農(nóng)業(yè)面源污染。

2.1.5 農(nóng)村人口受教育程度。農(nóng)村人口受教育程度越高，可能越容易掌握先進的農(nóng)業(yè)綠色發(fā)展技術(shù)，對農(nóng)業(yè)面源污染的認知水平和意愿也可能越高，因此有利于農(nóng)業(yè)面源污染的防治。

2.1.6 農(nóng)業(yè)面源污染防治政策投入。政府投入到農(nóng)業(yè)面源污染治理的資金越多，可能越有利于防治農(nóng)業(yè)面源污染。

因此，本文根據(jù)以上理論分析，建立中國農(nóng)業(yè)面源污染形成的經(jīng)濟機理分析框架，如圖3所示。

從圖3的理論分析框架，構(gòu)建中國農(nóng)業(yè)面源污染經(jīng)濟影響因素實證分析模型：

上式中，ANP表示面源污染物排放總量和排放總強度：其中，排放總量包括：氮排放總量（TNq），磷排放總量（TPq），化學(xué)需氧量排放總量（TCODq）；排放總強度包括：氮的排放總強度（TNs），磷的排放總強度（TPs）和化學(xué)需氧量排放總強度（TCODs）；gdpi表示i地區(qū)農(nóng)村經(jīng)濟總量；ngpi表示i地區(qū)農(nóng)業(yè)生產(chǎn)總值；frgpi表示i地區(qū)林業(yè)生產(chǎn)總值；stgpi表示i地區(qū)畜牧業(yè)生產(chǎn)總值；figpi表示i地區(qū)漁業(yè)生產(chǎn)總值；edui表示i地區(qū)受教育水平程度；rgpi表示i地區(qū)農(nóng)村人口規(guī)模；esi表示i地區(qū)農(nóng)業(yè)面源污染防治政策投入；μ表示隨機擾動項，β1表示常數(shù)項。

2.2 實證結(jié)果分析

本文利用STATA12.0軟件，運用2006～2015年中國31個省的省際面板數(shù)據(jù)，對中國農(nóng)業(yè)面源污染形成的經(jīng)濟影響因素進行了實證檢驗，由于使用面板數(shù)據(jù)進行實證分析，省際數(shù)據(jù)可能會存在解釋變量和個體效應(yīng)相關(guān)的現(xiàn)象，因此，本文先進行了豪斯曼檢驗（Hausman）檢驗，檢驗結(jié)果的P值為0.0000，拒絕原假設(shè)，即拒絕隨機效應(yīng)模型。因此，本模型采取固定效應(yīng)模型對中國農(nóng)業(yè)面源污染排放總量和排放總強度進行了實證分析，估計結(jié)果見表3和表4。

由表3和表4可以看出：

（1）農(nóng)村的經(jīng)濟總量gdp對氮的排放總量（TNq）、氮的排放總強度（TNs），磷的排放總量（TPq）、磷的排放總強度（TPs）和化學(xué)需氧量排放總量（TCODq）、化學(xué)需氧量排放總強度（TCODs）在1%的顯著水平下，有顯著的負效應(yīng)，這說明農(nóng)村經(jīng)濟發(fā)達，農(nóng)村居民有較強的環(huán)境保護意愿和高環(huán)境質(zhì)量需求，這也符合庫茲涅佐曲線。

（2）農(nóng)業(yè)經(jīng)濟總量雖然只有在10%的顯著水平下通過檢驗，但在一定程度上說明了農(nóng)業(yè)生產(chǎn)總量越大，對農(nóng)業(yè)面源污染有正影響，農(nóng)業(yè)的發(fā)展還是粗放型的發(fā)展方式。

（3）農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)對農(nóng)業(yè)面源污染的影響。從估計結(jié)果來看，林業(yè)生產(chǎn)總值對氮的排放總量（TNq）、磷的排放總量（TPq）在10%的顯著水平下通過檢驗，有顯著的負效應(yīng)，對化學(xué)需氧量排放總量（TCODq）在1%的顯著水平下，有顯著的負效應(yīng)，在對排放總強度上，在10%的顯著水平下通過檢驗，有顯著的負效應(yīng)；漁業(yè)，特別是畜牧業(yè)的生產(chǎn)總值對農(nóng)業(yè)面源污染的排放總量和排放總強度，都產(chǎn)生顯著的正影響。

（4）農(nóng)村人口受教育程度和農(nóng)業(yè)面源污染防治政策投入對農(nóng)業(yè)面源污染的排放總量和排放總強度，都產(chǎn)生了顯著的負影響。

3 研究結(jié)論與政策建議

3.1 研究結(jié)論

從以上對農(nóng)業(yè)面源污染的排放總量和排放總強度的估計結(jié)果分析來看，本文可以得到以下結(jié)論。

（1）中國的農(nóng)業(yè)面源污染符合庫茲涅佐環(huán)境曲線，當農(nóng)村經(jīng)濟總量達到一定程度后，農(nóng)村居民有較強的環(huán)境質(zhì)量需求和較強的環(huán)境環(huán)境保護意愿，有利于農(nóng)業(yè)面源污染的防治。

（2）農(nóng)業(yè)生產(chǎn)總規(guī)模，特別是農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)對農(nóng)業(yè)面源污染有顯著影響。隨著中國經(jīng)濟的發(fā)展，人們對魚肉產(chǎn)品需求的加大，漁業(yè)和畜牧產(chǎn)業(yè)的發(fā)展嚴重影響了農(nóng)業(yè)面源污染的防治，因此，漁業(yè)和畜牧業(yè)是農(nóng)業(yè)面源污染防治的重點產(chǎn)業(yè)和重點區(qū)域。

（3）農(nóng)村人口的教育程度對農(nóng)業(yè)面源污染有負效應(yīng)，說明提高農(nóng)村人口教育程度有利于防治農(nóng)村垃圾廢水等方面帶來的農(nóng)業(yè)面源污染。

（4）國家的農(nóng)業(yè)面源污染防治政策投入顯著地影響了農(nóng)業(yè)面源污染，說明國家的政策是有效的，對農(nóng)業(yè)面源的防治起到了很大的治理作用。

3.2 政策建議

根據(jù)以上研究結(jié)論，本文提出下列政策建議。

（1）發(fā)展農(nóng)村經(jīng)濟，特別是發(fā)展農(nóng)業(yè)面源污染嚴重的貧困地區(qū)的經(jīng)濟，才能使農(nóng)戶從內(nèi)生需求對農(nóng)業(yè)面源污染防治的內(nèi)生動力，這也是解決農(nóng)業(yè)面源污染的長效機制，否則，國家的防治政策可能只能起到短期的政策效應(yīng)。

（2）在當前農(nóng)業(yè)結(jié)構(gòu)調(diào)整和發(fā)展趨勢下，要把漁業(yè)和畜牧養(yǎng)殖業(yè)發(fā)展規(guī)模比較大的地區(qū)作為農(nóng)業(yè)面源污染防治的重點區(qū)域。

（3）加強農(nóng)村人口的教育水平和對農(nóng)業(yè)面源污染防治的認知教育，有利于農(nóng)業(yè)面源污染的防治，特別是生活產(chǎn)生的垃圾、廢水之類的農(nóng)業(yè)面源污染。

（4）中國的農(nóng)業(yè)面源污染防治政策發(fā)揮有效作用，政府應(yīng)繼續(xù)加大政策投入，繼續(xù)實施美麗鄉(xiāng)村建設(shè)的垃圾分類處理，廁所變革，農(nóng)業(yè)地膜回收，綠色化肥農(nóng)藥等項目，為實現(xiàn)中國美麗鄉(xiāng)村建設(shè)，生態(tài)文明建設(shè)和鄉(xiāng)村振興戰(zhàn)略提供政策保障。

[參考文獻]

[1] 張維理，武淑霞，冀宏杰，等.中國農(nóng)業(yè)面源污染形勢估計及控制對I.21世紀初期中國農(nóng)業(yè)面源污染的形勢估計[J].中國農(nóng)業(yè)科學(xué)，2004（7）：1008-1017.

[2] 金書秦，邢曉旭.農(nóng)業(yè)面源污染的趨勢研判、政策評述和對策建議[J].中國農(nóng)業(yè)科學(xué)，2018，51（3）：593-600.

[3] 梁流濤，馮淑怡，曲福田.農(nóng)業(yè)面源污染形成機制：理論與實證[J].中國人口資源與環(huán)境，2010（4）：74-80.

[4] 葛繼紅，周曙東.農(nóng)業(yè)面源污染的經(jīng)濟影響因素分析[J].中國農(nóng)村經(jīng)濟，2010（5）：72-80.

[5] 陳敏鵬，陳吉寧，等.中國農(nóng)業(yè)和農(nóng)村污染的清單分析與空間特征識別[J].中國環(huán)境科學(xué)， 2006，26（6）：751-755.

[6] Dinda，S.Environmental Kuznets Curve Hypothesis：A Survey[J].Ecological Economics，2004（48）：431-455.

[7] 孟祥海，周海川，張俊飚.中國畜禽污染時空特征分析與環(huán)境庫茲涅茨曲線驗證[J].干旱區(qū)資源與環(huán)境，2015，29（11）：104-108.