肖新成 何丙輝 倪九派 謝德體

摘要：農(nóng)業(yè)面源污染排放效率反映農(nóng)業(yè)經(jīng)濟(jì)發(fā)展與環(huán)境的協(xié)調(diào)程度，農(nóng)業(yè)面源污染物的影子價(jià)格體現(xiàn)了污染物的邊際減排成本。提高農(nóng)業(yè)生產(chǎn)的環(huán)境效率，降低農(nóng)業(yè)面源污染的影子價(jià)格對(duì)保護(hù)農(nóng)村生態(tài)環(huán)境和增強(qiáng)流域水資源的安全性起到較好的促進(jìn)作用。利用三峽生態(tài)屏障區(qū)重慶段2000-2012年間的面板數(shù)據(jù)，采用參數(shù)化方向性距離函數(shù)，對(duì)農(nóng)業(yè)面源污染的排放效率和影子價(jià)格進(jìn)行了測(cè)算，結(jié)合面板數(shù)據(jù)隨機(jī)效應(yīng)Tobit模型分析了影響農(nóng)業(yè)面源污染的排放效率和影子價(jià)格。結(jié)果表明：2000-2012年期間，農(nóng)業(yè)面源污染平均排放效率為0.649 9，TN、TP、COD和NH3-N 4種污染物的平均影子價(jià)格分別為0.372 9，0.032 6，0.137 1，0.053 3萬(wàn)元/t；不同區(qū)域農(nóng)業(yè)面源污染排放效率和影子價(jià)格存在一定的差異，庫(kù)中平行嶺谷區(qū)、庫(kù)區(qū)腹地沿江山地丘陵區(qū)和庫(kù)區(qū)腹地外圍山地區(qū)農(nóng)業(yè)面源污染排放效率低，污染物的影子價(jià)格高；影響農(nóng)業(yè)面源污染的排放效率和影子價(jià)格因素主要包括農(nóng)業(yè)結(jié)構(gòu)中經(jīng)濟(jì)作物與糧食作物的比例、農(nóng)業(yè)生產(chǎn)設(shè)施條件、農(nóng)村居民受教育程度。因此，確保糧食作物安全的前提下，進(jìn)一步優(yōu)化農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，減少農(nóng)藥、化肥投入大的農(nóng)作物種植規(guī)模，走“兩型農(nóng)業(yè)”之路；進(jìn)一步加大對(duì)農(nóng)田水利設(shè)施、節(jié)水灌溉設(shè)施和農(nóng)業(yè)生態(tài)建設(shè)工程的投資，減少和控制農(nóng)業(yè)生產(chǎn)中的污染物排放對(duì)生態(tài)環(huán)境的破壞；農(nóng)業(yè)面源污染源頭減排離不開農(nóng)戶的積極參與，提高農(nóng)戶對(duì)農(nóng)業(yè)生態(tài)環(huán)境退化的認(rèn)知也是解決農(nóng)業(yè)面源污染的重要途徑。

關(guān)鍵詞三峽生態(tài)屏障區(qū)；農(nóng)業(yè)面源污染；排放效率；影子價(jià)格；方向性距離函數(shù)

中圖分類號(hào)F323

文獻(xiàn)標(biāo)識(shí)碼A

文章編號(hào)1002-2104（2014）11-0060-09

在我國(guó)農(nóng)業(yè)集約化程度較高的大江大河流域，化肥、農(nóng)藥的大量投入和畜禽養(yǎng)殖業(yè)規(guī)模的擴(kuò)大對(duì)區(qū)域農(nóng)村經(jīng)濟(jì)的發(fā)展起到了舉足輕重的作用，同時(shí)農(nóng)業(yè)化學(xué)投入品的過(guò)量使用和畜禽糞便直接排入河流已成為流域水體富營(yíng)養(yǎng)化的主要來(lái)源[1-2]。三峽生態(tài)屏障區(qū)獨(dú)特的地理位置與氣候，使攔截、消納和過(guò)濾農(nóng)業(yè)面源污染的功能受到了嚴(yán)重的影響，加重了長(zhǎng)江污染[3]，致使庫(kù)區(qū)已經(jīng)成為國(guó)家生態(tài)環(huán)境建設(shè)重點(diǎn)治理區(qū)。在土地資源和環(huán)境約束的條件下，農(nóng)業(yè)經(jīng)濟(jì)的發(fā)展必須提高農(nóng)業(yè)生產(chǎn)效率，降低庫(kù)區(qū)農(nóng)業(yè)面源污染排放量，削減農(nóng)業(yè)面源污染的控制成本。因此，在衡量農(nóng)業(yè)和農(nóng)村發(fā)展水平時(shí)有必要將農(nóng)業(yè)面源污染因素納入到農(nóng)業(yè)生產(chǎn)效率的范疇，一方面可以更加全面地反映區(qū)域農(nóng)業(yè)和農(nóng)村經(jīng)濟(jì)增長(zhǎng)的數(shù)量與質(zhì)量，另一方面也可以為農(nóng)業(yè)和環(huán)境保護(hù)部門提供發(fā)展農(nóng)業(yè)經(jīng)濟(jì)，保護(hù)農(nóng)村環(huán)境的建議。

由于農(nóng)業(yè)面源污染排放物沒(méi)有相應(yīng)的交易市場(chǎng)，難以獲取相應(yīng)的市場(chǎng)價(jià)格[4]，因而農(nóng)業(yè)生產(chǎn)的污染物排放效率和污染物帶來(lái)的損失并不容易被估計(jì)。Shephard[5]提出的產(chǎn)出和投入距離函數(shù)，無(wú)需產(chǎn)出最大化、成本最小化假設(shè)和要素投入的價(jià)格信息，為學(xué)者們研究環(huán)境經(jīng)濟(jì)效率提供了很好的方法。Chung[6]引入方向性距離函數(shù)，提出了MalmquistLuenberger生產(chǎn)率指數(shù)，分析瑞典紙漿行業(yè)的全要素生產(chǎn)率。這種方法能夠較好地分析合意產(chǎn)出的增加和不合意產(chǎn)出的減少，并且不需要考慮生產(chǎn)要素的價(jià)格信息。Fare[7]使用二次方向性產(chǎn)出距離函數(shù)估算了美國(guó)農(nóng)業(yè)生產(chǎn)效率、影子價(jià)格和相關(guān)成本。Karimov[8]使用二次方向性投入距離函數(shù)對(duì)烏茲別克斯坦蔬菜種植農(nóng)場(chǎng)的經(jīng)濟(jì)效率和影子價(jià)格進(jìn)行了分析。近年來(lái)國(guó)內(nèi)學(xué)者也使用方向性距離函數(shù)對(duì)我國(guó)的農(nóng)業(yè)環(huán)境效率進(jìn)行了研究。楊俊[9]考慮了農(nóng)業(yè)非點(diǎn)源污染，運(yùn)用方向性距離函數(shù)構(gòu)建ML生產(chǎn)率指數(shù)對(duì)中國(guó)農(nóng)業(yè)全要素生產(chǎn)率進(jìn)行了分析。梁流濤[10]利用中國(guó)省際面板數(shù)據(jù)和方向性距離函數(shù)評(píng)價(jià)了我國(guó)各區(qū)域農(nóng)業(yè)發(fā)展與環(huán)境的協(xié)調(diào)程度。潘丹[11-12]通過(guò)使用SBM方向性距離函數(shù)，測(cè)算了環(huán)境污染約束下農(nóng)業(yè)生產(chǎn)率地區(qū)差距的動(dòng)態(tài)分布演進(jìn)。閔銳[13]采用超效率DEA及方向性距離函數(shù)對(duì)國(guó)內(nèi)糧食生產(chǎn)與資源環(huán)境協(xié)調(diào)性空間分異特征進(jìn)行了分析，發(fā)現(xiàn)環(huán)境污染因素的引入對(duì)我國(guó)整體及各省區(qū)糧食生產(chǎn)技術(shù)效率產(chǎn)生了一定影響。國(guó)內(nèi)外環(huán)境技術(shù)效率的研究基本上采用非參數(shù)化的DEA方法，這種方法從投入或產(chǎn)出的角度來(lái)分析，無(wú)需考慮利潤(rùn)最大化的問(wèn)題[14-15]，但是，他們很少研究多種污染物條件下的農(nóng)業(yè)面源污染排放效率和污染物的影子價(jià)格，在一定程度上會(huì)影響到農(nóng)業(yè)面源污染排放效率評(píng)價(jià)和影子價(jià)格的準(zhǔn)確性。本文采用參數(shù)化的方向性距離函數(shù)并結(jié)合分參數(shù)化DEA方法的優(yōu)點(diǎn)，分析多種農(nóng)業(yè)面源污染物存在下的排放效率和影子價(jià)格，主要是因?yàn)閰?shù)化方程考慮了合意產(chǎn)出與不合意產(chǎn)出后的利潤(rùn)最大化問(wèn)題，能夠計(jì)算不合意產(chǎn)出與合意產(chǎn)出的彈性之比，而且使用參數(shù)化方程可以測(cè)算出多種不合意產(chǎn)出條件下農(nóng)業(yè)面源污染物的影子價(jià)格和排放效率，相對(duì)其他方法更為簡(jiǎn)潔方便。

1研究方法

1.1變量界定與數(shù)據(jù)選取

農(nóng)業(yè)生產(chǎn)過(guò)程中的生產(chǎn)要素投入不僅會(huì)生產(chǎn)人們希望獲得的合意產(chǎn)出（糧食、肉類產(chǎn)品、GDP），同時(shí)也不可避免地伴隨著一些不合意產(chǎn)出（農(nóng)業(yè)面源污染物TN、TP、COD、NH3N）在內(nèi)的產(chǎn)出。因此，本文選取農(nóng)業(yè)生產(chǎn)資金投入、土地、勞動(dòng)力作為農(nóng)業(yè)生產(chǎn)投入要素，農(nóng)業(yè)GDP為合意產(chǎn)出，TN、TP、COD、NH3N為不合意產(chǎn)出，農(nóng)業(yè)生產(chǎn)資金投入、土地、勞動(dòng)力和農(nóng)業(yè)GDP的數(shù)據(jù)來(lái)源于庫(kù)區(qū)各區(qū)縣2000-2012年間的統(tǒng)計(jì)年鑒，綜合陳敏鵬[16]、賴斯蕓[17]的清單分析法和錢曉雍[18]的等標(biāo)污染負(fù)荷法，估算三峽生態(tài)屏障區(qū)重慶段2000-2012年TN、TP、COD、NH3N排放量[19-20]。

1.2方向性距離函數(shù)

衡量農(nóng)業(yè)生產(chǎn)效率的傳統(tǒng)方法只關(guān)注農(nóng)業(yè)投入與產(chǎn)出，而忽略了農(nóng)業(yè)生產(chǎn)也會(huì)造成環(huán)境污染，這就需要構(gòu)建農(nóng)業(yè)生產(chǎn)中既包含合意產(chǎn)出也包含不合意產(chǎn)出的生產(chǎn)可能性集合，測(cè)算農(nóng)業(yè)投入與合意產(chǎn)出和不合意產(chǎn)出之間的技術(shù)結(jié)構(gòu)關(guān)系即農(nóng)業(yè)污染物排放效率[21]。假定農(nóng)業(yè)生產(chǎn)使用N種投入x=（x1，x2，…，xN）∈RN+，能夠得到M種合意產(chǎn)出y=（y1，y2，…，yM）∈RM+和T種不合意產(chǎn)出z=（z1，z2，…，zT）∈RT+，其污染物排放效率的生產(chǎn)可行性集可以表示為：

P（x）={（y，z）：x可以生產(chǎn)（y，z）}，x∈RN+（1）

根據(jù)Chunge等和Fare等的研究[6-7]，在投入產(chǎn)出理論中，污染物排放生產(chǎn)的可行性集P（x）是一個(gè)凸的、有界的閉集合，在P（x）中有限的生產(chǎn)要素投入只能生產(chǎn)出有限的產(chǎn)出。因此，生產(chǎn)可行性集具有如下的性質(zhì)：①合意產(chǎn)出和不合意產(chǎn)出具有聯(lián)合弱可處置性。如果（y，z）∈P（x），并且0≤θ≤1，那么（θy，θz）∈P（x），這表明在投入一定的情況下，減少不合意產(chǎn)出是有成本的，減少不合意產(chǎn)出也會(huì)減少合意產(chǎn)出，因此要控制農(nóng)業(yè)面源污染必須付出一定的代價(jià)，也就是說(shuō)存在如果z′

Shephard產(chǎn)出距離函數(shù)尋求的是投入增加的情況下，合意產(chǎn)出與不合意產(chǎn)出的同時(shí)增加，而方向性距離函數(shù)尋求的是在投入增加的情況下，合意產(chǎn)出增加，不合意產(chǎn)出減少，因此，方向性距離函數(shù)是Shephard產(chǎn)出距離函數(shù)的一般化[22-23]。根據(jù)Fare[24]和Anders[25]對(duì)方向性距離的描述，引入方向性向量g=（gy，gz），且g≠0，則方向性距離產(chǎn)出函數(shù)就可以表示為：

Do（x，y，z；gy，gz）=max{β：（y+βy，z-βz）}∈P（x）（2）

該方向性距離函數(shù)表示在給定P（x）的條件下，最大限度范圍內(nèi)擴(kuò)充合意產(chǎn)出的同時(shí)減少不合意產(chǎn)出是可行的。

1.3農(nóng)業(yè)面源污染影子價(jià)格的推導(dǎo)

農(nóng)業(yè)生產(chǎn)過(guò)程中，會(huì)伴隨農(nóng)業(yè)面源污染物TN、TP、COD、NH3N的產(chǎn)生，這些不合意產(chǎn)出通常沒(méi)有市場(chǎng)價(jià)格，需要引入影子價(jià)格的概念，污染物的影子價(jià)格實(shí)際反映了污染削減的邊際成本[26-27]。對(duì)不合意產(chǎn)出影子價(jià)格的估計(jì)可以從方向性距離函數(shù)與利潤(rùn)函數(shù)之間的關(guān)系中推導(dǎo)出來(lái)。本文在參照Chung和Fare提出的環(huán)境污染影子價(jià)格方向性距離函數(shù)的基礎(chǔ)上，利用超越對(duì)數(shù)生產(chǎn)函數(shù)來(lái)推導(dǎo)農(nóng)業(yè)面源污染的影子價(jià)格的參數(shù)化形式。為便于分析，將農(nóng)業(yè)生產(chǎn)中1種合意產(chǎn)出與4種不合意產(chǎn)出的超越對(duì)數(shù)產(chǎn)出函數(shù)定義為：

lnD0=α0+∑3i=1αilnxi+β1lny+∑4i=1γilnzi

+12∑3i=1∑3j=1αijlnxilnxj+12β11（lny）2

+12∑4i=1∑4j=1γijlnzilnzj+∑3i=1θilnxilny

+∑3i=1∑4j=1εijlnxilnzj+∑4i=1jlnylnzj

（3）

其中，x1，x2，x3為農(nóng)業(yè)生產(chǎn)性資金投入、耕地面積、農(nóng)業(yè)勞動(dòng)力三種要素的投入量，y為合意產(chǎn)出GDP，z1，z2，z3，z4為不合意產(chǎn)品TN、TP、COD、NH3N的產(chǎn)出數(shù)量。D0表示含有農(nóng)業(yè)面源污染物的產(chǎn)出距離函數(shù)，且0

β1-（∑4i=1γi）=-1，β11=（∑4j=1γij）=∑4j=1φj（i=1，2，3，4），（∑3j=1θj）=∑4j=1εij（i=1，2，3）

αij=αji，γij=γji，lnD0lnxi≤0，lnD0lny≥0lnD0lnzi≤0，

假設(shè)R為農(nóng)業(yè)的產(chǎn)出收益，其收益函數(shù)應(yīng)為：

R=p[y+（1-D0）gy]-{pz1[z1-（1-D0）gz1]

+pz2[z2-（1-D0）gz2]+pz3[z3-（1-D0）gz3]

+pz4[z4-（1-D0）gz4]}（4）

其中，p為合意產(chǎn)出的價(jià)格，pz1，pz2，pz3，pz4為4種不合意產(chǎn)出的價(jià)格，（1-D0）表示產(chǎn)出的無(wú)效率，則gy是增加合意產(chǎn)出的固定量，且gz1，gz2，gz3，gz4分別為減少4種不合意產(chǎn)出的固定量，那么，農(nóng)業(yè)生產(chǎn)的成本函數(shù)為：

C=ω1x1+ω2x2+ω3x3（5）

其中，ω1，ω2，ω3分別為三種要素投入的價(jià)格，x1，x2，x3為三種要素投入的數(shù)量。則農(nóng)業(yè)生產(chǎn)的利潤(rùn)函數(shù)為：

π=R-C=p[y+（1-D0）gy]

-{pz1[z1-（1-D0）gz1]+pz2[z2-（1-D0）gz2]

+pz3[z3-（1-D0）gz3]+pz4[z4-（1-D0）gz4]}

-（ω1x1+ω2x2+ω3x3）（6）

按照農(nóng)業(yè)生產(chǎn)實(shí)現(xiàn)利潤(rùn)最大化的原則，分別對(duì)式（6）的y，z1，z2，z3，z4求偏導(dǎo)數(shù)并令它們等于0。

πy=p+（pgy+pz1gz1+pz2gz2+pz3gz3+pz4gz4）

×（-D0y）=0（7）

πz1=-pz1+（pgy+pz1gz1+pz2gz2+pz3gz3+pz4gz4）

×（-D0z1）=0（8）

πz2=-pz1+（pgy+pz1gz1+pz2gz2+pz3gz3+pz4gz4）

×（-D0z2）=0（9）

πz3=-pz1+（pgy+pz1gz1+pz2gz2+pz3gz3+pz4gz4）

×（-D0z3）=0（10）

πz4=-pz1+（pgy+pz1gz1+pz2gz2+pz3gz3+pz4gz4）

×（-D0z4）=0（11）

將式（7）分別代入式（8）、（9）、（10）、（11）式得：

pz1=（-D0z1/D0z）p，

pz2=（-D0z2/D0z）p

pz3=（-D0z3/D0z）p，

pz4=（-D0z4/D0z）p（12）

由于式（3）是關(guān)于合意產(chǎn)出與不合意產(chǎn)出的一階齊次方程，為了計(jì)算不合意產(chǎn)出的影子價(jià)格，需要將將式（1）兩邊用合意產(chǎn)出（或不合意產(chǎn)出）正規(guī)化，結(jié)合約束條件，運(yùn)用Lingo9.0軟件編程計(jì)算農(nóng)業(yè)面源污染的排放效率。并計(jì)算出lnD0lny、lnD0lnz1、lnD0lnz2、lnD0lnz3、lnD0lnz4，那么式（12）中4種不合意產(chǎn)出的影子價(jià)格又可以表示為：pz1=-yz1（lnD0lnz1/lnD0lny）p，pz2=-yz2（-lnD0lnz2/lnD0lny）p，pz3=-yz3（-lnD0lnz3/lnD0lny）p，pz4=-yz4（-lnD0lnz4/lnD0lny）p，其中l(wèi)nD0lny、lnD0lnz1、-lnD0lnz2、-lnD0lnz3、-lnD0lnz4分別為農(nóng)業(yè)面源污染排放效率對(duì)合意產(chǎn)出和4種不合意產(chǎn)出的彈性，其經(jīng)濟(jì)學(xué)的含義是：合意產(chǎn)出每增加1%引起農(nóng)業(yè)面源污染排放效率的增長(zhǎng)率，反之亦然；不合意產(chǎn)出每減少1%引起農(nóng)業(yè)面源污染排放效率的增長(zhǎng)率，反之亦然。通過(guò)lnD0（x1，x2，x3，y，z1，z2，z3，z4）=-μ，μ≥0，可以直接求出農(nóng)業(yè)面源污染的排放效率為：D0=exp（-μ）。

2研究結(jié)果

2.1農(nóng)業(yè)面源污染物的影子價(jià)格分析

2.1.1估計(jì)結(jié)果分析與檢驗(yàn)

通過(guò)對(duì)2000-2012年期間三峽生態(tài)屏障區(qū)重慶段21個(gè)區(qū)縣的投入、產(chǎn)出變量進(jìn)行規(guī)范化處理后，采用上述約束條件對(duì)方向性產(chǎn)出距離函數(shù)進(jìn)行參數(shù)估計(jì)，具體結(jié)果見表1。根據(jù)參數(shù)估計(jì)值，可以對(duì)參數(shù)化超越對(duì)數(shù)方向性距離函數(shù)進(jìn)行相關(guān)檢驗(yàn)。首先，檢驗(yàn)合意產(chǎn)出與不合意產(chǎn)出的零結(jié)合性假定條件能否得到了較好的適配。在上述的模型分析中，可得知（y，z）∈P（x）與D0（x，y，z；gy，gz）具有一定的等價(jià)性。因此，可以分析y>0的情況下，計(jì)算D0（x，y，0；gy，gz）的值來(lái)判斷零結(jié)合性成立，即沒(méi)有不合意產(chǎn)出也就不會(huì)有合意

產(chǎn)出，不合意產(chǎn)出是生產(chǎn)合意產(chǎn)出過(guò)程中不可避免的附屬

產(chǎn)品。對(duì)21個(gè)區(qū)縣的273個(gè)樣本（21×13）的檢驗(yàn)，結(jié)果得到有249個(gè)樣本符合零結(jié)合性，占樣本的比例達(dá)到了93.77%，絕大部分樣本滿足該條件，說(shuō)明合意產(chǎn)出與不合意產(chǎn)出的零結(jié)合性條件成立。其次，檢驗(yàn)0

2.1.2農(nóng)業(yè)面源污染排放效率分析

根據(jù)方向性距離函數(shù)的理論，采用參數(shù)化的超越對(duì)數(shù)生產(chǎn)函數(shù)，測(cè)算出2000-2012年三峽生態(tài)屏障區(qū)重慶段21個(gè)區(qū)縣的農(nóng)業(yè)面源污染排放的效率，根據(jù)農(nóng)業(yè)面源污染排放效率取值的大小衡量農(nóng)業(yè)經(jīng)濟(jì)發(fā)展與環(huán)境的相對(duì)協(xié)調(diào)程度。2000-2012年期間農(nóng)業(yè)面源污染平均排放效率為0.649 9，并且呈現(xiàn)“V”形波動(dòng)狀態(tài)（見表2）。為了進(jìn)一步分析協(xié)調(diào)性的動(dòng)態(tài)變化，將2000-2012年分為三個(gè)階段（2000-2004年、2005-2008年和2009-2012年），這三個(gè)階段協(xié)調(diào)度的平均值分別為0.645 8、0.620 5、0.684 4。尤其是第三個(gè)階段農(nóng)業(yè)經(jīng)濟(jì)發(fā)展與環(huán)境協(xié)調(diào)性指數(shù)有了較大的提高，主要原因在于區(qū)域大規(guī)模的農(nóng)業(yè)生產(chǎn)基礎(chǔ)設(shè)施投資、大力推行農(nóng)業(yè)清潔生產(chǎn)、重視農(nóng)業(yè)面源污染的危害性、大力加強(qiáng)農(nóng)業(yè)和農(nóng)村環(huán)境的整治等措施，提升了農(nóng)業(yè)經(jīng)濟(jì)發(fā)展與農(nóng)村環(huán)境的協(xié)調(diào)性。

根據(jù)三峽生態(tài)屏障區(qū)重慶段地理位置，將其劃分為庫(kù)尾都市核心區(qū)（渝中、大渡口、江北、沙坪壩、九龍坡和南岸）、庫(kù)尾都市外圍區(qū)（北碚、渝北和巴南）、庫(kù)尾低山丘陵區(qū)（江津）、庫(kù)中平行嶺谷區(qū)（長(zhǎng)壽、涪陵、武隆、豐都、石柱和忠縣）、腹地沿江山地丘陵區(qū)（萬(wàn)州、云陽(yáng)和奉節(jié)）和腹地外圍山地（開縣、巫山和巫溪），不同階段各區(qū)域的農(nóng)業(yè)面源污染平均排放效率如表3。農(nóng)業(yè)面源污染排放效率最高的是庫(kù)尾都市核心區(qū)，其次是庫(kù)尾都市外圍區(qū)和庫(kù)尾低山丘陵區(qū)，效率最低的是庫(kù)中平行嶺谷區(qū)。在經(jīng)濟(jì)發(fā)展水平較高的都市核心區(qū)、都市外圍區(qū)，人們比較關(guān)注農(nóng)產(chǎn)品質(zhì)量和環(huán)境質(zhì)量，區(qū)域產(chǎn)業(yè)結(jié)構(gòu)的及時(shí)調(diào)整，使得農(nóng)業(yè)污染物排放效率得到了較大的提高。

2.2農(nóng)業(yè)面源污染影子價(jià)格測(cè)算

在獲取了方向性距離函數(shù)的參數(shù)后，可以計(jì)算農(nóng)業(yè)面源污染4種污染物的影子價(jià)格。由于選取的合意產(chǎn)出為

農(nóng)業(yè)生產(chǎn)總值，將合意產(chǎn)出的價(jià)格定為1，估算出的2000-2012年4種污染物的影子價(jià)格如表4。農(nóng)業(yè)面源污染物的影子價(jià)格反映的是減少1單位的污染物排放量，將會(huì)產(chǎn)生多少合意產(chǎn)出（農(nóng)業(yè)產(chǎn)值）的損失。從表中可以看出，在2000-2012年期間，TN、TP、COD、NH3N4種污染物的平均影子價(jià)格分別為0.372 9萬(wàn)元/t、0.032 6萬(wàn)元/t、0.137 1萬(wàn)元/t、0.053 3萬(wàn)元/t，即減少1tTN、TP、COD、NH3N帶來(lái)的農(nóng)業(yè)總產(chǎn)值的損失分別為0.372 9萬(wàn)元、0.032 6萬(wàn)元、0.137 1萬(wàn)元、0.053 3萬(wàn)元。4種農(nóng)業(yè)面源污染物之間影子價(jià)格相差比較大，因而處理不同的農(nóng)業(yè)面源污染物面臨的成本壓力有所差別。從4種污染物的平均影子價(jià)格看，對(duì)TN的處理難度較大，對(duì)TP的處理成本最小，4種污染物在庫(kù)區(qū)排放規(guī)模巨大，污染物總成本仍然相當(dāng)高。從污染物影子價(jià)格的年份來(lái)看，2000-2012年期間污染物的平均影子價(jià)格有所波動(dòng)，但總體上波動(dòng)幅度不大，比較平穩(wěn)。2000-2006年間各種污染物的影子價(jià)格

處于上升趨勢(shì)，2006-2012年處于下降的趨勢(shì)。

分區(qū)域來(lái)看，同一污染物在不同區(qū)域之間的影子價(jià)格也存在較大的差異，各區(qū)域2000-2012年間的平均影子價(jià)格如表5。庫(kù)區(qū)庫(kù)中平行嶺谷區(qū)4種污染物的平均影子價(jià)格最高，是因?yàn)樵搮^(qū)是重慶的糧倉(cāng)，耕地分布最為集中，農(nóng)業(yè)活動(dòng)較為活躍，這與地形相對(duì)平緩、農(nóng)業(yè)產(chǎn)業(yè)化程度較高、農(nóng)墾歷史悠久等有很大關(guān)系，同時(shí)該區(qū)域也具有畜禽養(yǎng)殖所需要的資源、場(chǎng)地、人力等條件，是重慶市主要的畜禽養(yǎng)殖基地分布區(qū)；庫(kù)區(qū)腹地沿江山地丘陵區(qū)和庫(kù)區(qū)腹地外圍山地區(qū)的影子價(jià)格次之，因地勢(shì)起伏較大，農(nóng)業(yè)生產(chǎn)便捷性程度較平行嶺谷區(qū)低，但是兩區(qū)獨(dú)特的立體氣候條件逐漸被利用，發(fā)展獨(dú)具山地特色的名、特、優(yōu)和反季節(jié)高附加值產(chǎn)品，化肥施用量逐年增加，山地特色畜禽養(yǎng)殖較為發(fā)達(dá)，為庫(kù)區(qū)第二大畜禽養(yǎng)殖分布區(qū)；庫(kù)區(qū)庫(kù)尾都市核心區(qū)的污染物平均影子價(jià)格較低，是由于該區(qū)域農(nóng)業(yè)活動(dòng)逐漸朝都市休閑、觀光和旅游方向轉(zhuǎn)變，體現(xiàn)景觀、生態(tài)、服務(wù)等觀賞價(jià)值。4種農(nóng)業(yè)面源污染物各區(qū)域的平均影子價(jià)格也存在一定差異，這意味著不同區(qū)域在農(nóng)業(yè)面源污染減排成本方面具有差異，在這種情況下，完全依靠命令控制型措施來(lái)控制污染不能有效地削減污染成本，也不能起到良好的減排效果，調(diào)整農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)、實(shí)施農(nóng)業(yè)清潔生產(chǎn)、合理引導(dǎo)農(nóng)戶從事環(huán)境友好型農(nóng)業(yè)行為和提高農(nóng)戶環(huán)境保護(hù)意識(shí)等措施可能更為有效。

3農(nóng)業(yè)面源污染排放效率和影子價(jià)格影響因素的實(shí)證分析

采用參數(shù)化方向性距離函數(shù)估計(jì)得到的三峽生態(tài)屏障區(qū)重慶段的農(nóng)業(yè)面源污染排放效率值明顯低于梁流濤[28]對(duì)重慶結(jié)果估計(jì)。梁流濤估計(jì)的1997-2009重慶市農(nóng)業(yè)污染物排放效率的平均值為0.712 0，而我們采用參數(shù)化方向性距離函數(shù)估計(jì)的2000-2012年期間三峽生態(tài)屏障區(qū)重慶段的農(nóng)業(yè)污染物排放效率的平均值為0.649 9。導(dǎo)致結(jié)果差異的原因除了方法上和時(shí)間段的差異外，在此，實(shí)證選取的是庫(kù)區(qū)的21個(gè)區(qū)縣，其它區(qū)縣因

不屬于庫(kù)區(qū)未予以考慮，同時(shí)我們選取的不合意產(chǎn)出不僅僅是TN、TP、COD，還包括了NH3N，我們采用的數(shù)據(jù)是縣域數(shù)據(jù)，相比梁流濤采用的省域數(shù)據(jù)更為微觀，數(shù)據(jù)越微觀，結(jié)果越接近真實(shí)的農(nóng)業(yè)面源污染排放效率。

農(nóng)業(yè)面源污染排放效率和影子價(jià)格在不同的區(qū)域存在較大的差別，同一種污染物的影子價(jià)格在不同區(qū)縣也存在巨大差異。然而，農(nóng)業(yè)面源污染的發(fā)生來(lái)自兩部分，其一是污染物遷移與轉(zhuǎn)化過(guò)程，主要是降雨徑流的遷移過(guò)程；其二是污染物的產(chǎn)生過(guò)程，主要是人為污染排放等[29]。因此，自然與社會(huì)經(jīng)濟(jì)因素是影響農(nóng)業(yè)面源污染排放效率和影子價(jià)格的主要變量。由于農(nóng)業(yè)生產(chǎn)具有規(guī)模小和分散經(jīng)營(yíng)的特點(diǎn)，農(nóng)戶生產(chǎn)經(jīng)營(yíng)更多地受家庭特征的影響，有關(guān)農(nóng)業(yè)面源污染排放效率和影子價(jià)格的決定因素主要從農(nóng)村家庭特征和降雨徑流的角度考慮。首先是經(jīng)濟(jì)因素，包括農(nóng)戶收入水平和農(nóng)產(chǎn)品價(jià)格指數(shù)。農(nóng)產(chǎn)品價(jià)格的提高會(huì)使農(nóng)戶收入增加，一方面可能會(huì)使農(nóng)戶通過(guò)加大農(nóng)藥、化肥等生產(chǎn)要素的投入獲取更多的農(nóng)業(yè)產(chǎn)出，但不可避免地造成污染物排放量的增加[30]，影響排放效率和影子價(jià)格；另一方面也有可能會(huì)使農(nóng)民投入高質(zhì)量的生產(chǎn)資料（控釋肥），這有利于農(nóng)作物吸收利用，從而有利于排放效率的提高和影子價(jià)格的降低；其次是農(nóng)戶的受教育水平與農(nóng)業(yè)技術(shù)培訓(xùn)。農(nóng)戶教育水平提高與農(nóng)業(yè)技術(shù)培訓(xùn)讓農(nóng)戶掌握了科學(xué)的施肥方法，農(nóng)戶化肥、農(nóng)藥的施用會(huì)考慮土壤作物和施肥制度等因素，能夠在生產(chǎn)中根據(jù)自身的情況合理施肥用藥，更有效地發(fā)揮化肥、農(nóng)藥的增產(chǎn)作用。農(nóng)戶教育水平提高與農(nóng)業(yè)技術(shù)培訓(xùn)也能使農(nóng)戶在畜禽養(yǎng)殖上做到更科學(xué)合理地循環(huán)利用畜禽糞便。再次是農(nóng)戶的耕地規(guī)模、產(chǎn)業(yè)結(jié)構(gòu)比重和農(nóng)村基礎(chǔ)設(shè)施。三峽庫(kù)區(qū)坡耕地約占耕地面積的95.3%，坡耕地面積比重較大，且主要是緩坡耕地[29]，考慮人均耕地面積變量時(shí)對(duì)耕地的地勢(shì)不再作細(xì)分。農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中畜禽養(yǎng)殖業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值比重和經(jīng)濟(jì)作物與糧食作物的比例都是影響農(nóng)業(yè)面源污染的重要因素。加強(qiáng)以農(nóng)田水利為重點(diǎn)的農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)和大力發(fā)展節(jié)水灌溉，完善農(nóng)田水利建設(shè)管護(hù)機(jī)制，有利于減緩農(nóng)業(yè)面源污染；最后降雨徑流是影響農(nóng)業(yè)面源污染發(fā)生的重要自然因素，一般來(lái)說(shuō)，降水量越大，水土流失就越嚴(yán)重，污染程度越高。這樣，函數(shù)所包含的變量有：區(qū)域農(nóng)業(yè)經(jīng)濟(jì)發(fā)展水平（農(nóng)業(yè)人均GDP），農(nóng)產(chǎn)品價(jià)格指數(shù)（FPI），農(nóng)村居民受教育程度用平均受教育年限（EDU），農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)分別用畜禽養(yǎng)殖業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值比重（Animal）、經(jīng)濟(jì)作物與糧食作物的比例（PS）、人均耕地面積（Land），農(nóng)業(yè)生產(chǎn)設(shè)施條件用農(nóng)業(yè)基礎(chǔ)設(shè)施的投資（Inv）和降雨徑流（RR）等方面為解釋變量，以農(nóng)業(yè)面源污染排放效率（ETE）和污染物的影子價(jià)格（SP）為被解釋變量，使用三峽生態(tài)屏障區(qū)重慶段21個(gè)區(qū)縣2000-2012年數(shù)據(jù)，分別構(gòu)建面板數(shù)據(jù)隨機(jī)效應(yīng)Tobit模型：

模型的估計(jì)結(jié)果如表6所示。農(nóng)業(yè)面源污染排放效率和影子價(jià)格影響因素的回歸結(jié)果可以發(fā)現(xiàn)，兩個(gè)方程均在1%的水平上顯著。從農(nóng)業(yè)面源污染排放效率的影響因素看，農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中經(jīng)濟(jì)作物所占的比重越大效率越低，說(shuō)明提高農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中經(jīng)濟(jì)作物與糧食作物的比例在一定程度上加重了農(nóng)業(yè)面源污染；農(nóng)業(yè)生產(chǎn)設(shè)施條件的改善能有效地提高排放效率，農(nóng)業(yè)基礎(chǔ)設(shè)施投入增加，對(duì)改善農(nóng)業(yè)生產(chǎn)條件，提高農(nóng)業(yè)綜合生產(chǎn)能力，減少農(nóng)業(yè)生產(chǎn)中的污染物排放，促進(jìn)農(nóng)業(yè)生態(tài)環(huán)境的改善能夠起到積極的作用；農(nóng)戶受教育水平程度能影響農(nóng)業(yè)面源污染排放效率，農(nóng)戶教育水平的提高，將方便他們獲取農(nóng)業(yè)方面的信息，減少農(nóng)業(yè)生產(chǎn)資料的投入，科學(xué)合理地施肥；農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中畜禽養(yǎng)殖業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值比重越大排放效率越低，主要是由于在庫(kù)區(qū)，畜禽養(yǎng)殖的糞便沒(méi)有經(jīng)過(guò)處理，直接排放導(dǎo)致農(nóng)業(yè)面源污染排放壓力增大；降雨徑流越大，農(nóng)業(yè)面源污染排放效率越低，其它變量對(duì)農(nóng)業(yè)面源污染排放效率的作用不很明顯。從農(nóng)業(yè)面源污染的影子價(jià)格影響因素看，農(nóng)業(yè)結(jié)構(gòu)中糧食作物與經(jīng)濟(jì)作物的比例、農(nóng)業(yè)生產(chǎn)設(shè)施投資、農(nóng)業(yè)人均GDP、農(nóng)村居民受教育程度與農(nóng)業(yè)面源污染的影子價(jià)格負(fù)相關(guān)，這些變量會(huì)降低農(nóng)業(yè)面源污染的影子價(jià)格即邊際減排成本越低；農(nóng)產(chǎn)品價(jià)格指數(shù)、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)畜禽養(yǎng)殖業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值比重、人均耕地面積和降雨徑流的提高都會(huì)使影子價(jià)格提高即邊際減排成本增加。無(wú)論從提高農(nóng)業(yè)面源污染的排放效率還是從降低影子價(jià)格來(lái)看，需要采取調(diào)整農(nóng)業(yè)結(jié)構(gòu)中經(jīng)濟(jì)作物與糧食作物的比例、改善農(nóng)業(yè)生產(chǎn)設(shè)施條件和提高農(nóng)村居民受教育程度等措施，從源頭調(diào)控農(nóng)業(yè)面源污染。

4結(jié)論

提高農(nóng)業(yè)生產(chǎn)的環(huán)境效率，降低農(nóng)業(yè)面源污染的影子價(jià)格是控制農(nóng)業(yè)面源污染的重要組成部分，這對(duì)保護(hù)農(nóng)村生態(tài)環(huán)境和增強(qiáng)流域水資源的安全性起到較好的促進(jìn)作用。通過(guò)采用參數(shù)化的方向性距離函數(shù)，對(duì)2000-2012年三峽生態(tài)屏障區(qū)重慶段21個(gè)區(qū)縣的農(nóng)業(yè)面源污染的排放效率和影子價(jià)格進(jìn)行了估計(jì)，并對(duì)影響農(nóng)業(yè)面源污染排放效率和農(nóng)業(yè)面源污染影子價(jià)格的影響因素進(jìn)行了分析。主要結(jié)論如下：

2000-2012年期間，三峽生態(tài)屏障區(qū)重慶段農(nóng)業(yè)面源污染平均排放效率為0.649 9，并且呈現(xiàn)“V”形波動(dòng)狀態(tài)。農(nóng)業(yè)面源污染排放效率最高的是庫(kù)尾都市核心區(qū)，其次是庫(kù)尾都市外圍區(qū)和庫(kù)尾低山丘陵區(qū)，效率最低的是庫(kù)中平行嶺谷區(qū)。TN、TP、COD、NH3N4種污染物的平均影子價(jià)格分別為0.372 9萬(wàn)元/t、0.032 6萬(wàn)元/t、0.137 1萬(wàn)元/t、0.053 3萬(wàn)元/t。庫(kù)區(qū)庫(kù)中平行嶺谷區(qū)4種污染物的平均影子價(jià)格最高，庫(kù)區(qū)腹地沿江山地丘陵區(qū)和庫(kù)區(qū)腹地外圍山地區(qū)的影子價(jià)格次之，庫(kù)區(qū)庫(kù)尾都市核心區(qū)和污染物的平均影子價(jià)格最低。從農(nóng)業(yè)面源污染排放效率和影子價(jià)格來(lái)看，庫(kù)中平行嶺谷區(qū)、庫(kù)區(qū)腹地沿江山地丘陵區(qū)和庫(kù)區(qū)腹地外圍山地區(qū)農(nóng)業(yè)面源污染排放的效率低，污染物的影子價(jià)格高，它們是庫(kù)區(qū)農(nóng)業(yè)面源污染治理的重點(diǎn)地區(qū)。

從影響農(nóng)業(yè)面源污染排放效率的因素看，農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中經(jīng)濟(jì)作物所占的比重越低、農(nóng)業(yè)生產(chǎn)設(shè)施條件的改善、農(nóng)戶受教育水平程度的提升能夠提高農(nóng)業(yè)面源污染排放效率；農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中畜禽養(yǎng)殖業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值比重越大和降雨徑流值越高，排放效率越低，其他變量對(duì)農(nóng)業(yè)面源污染排放效率的作用不很明顯。從影響農(nóng)業(yè)面源污染的影子價(jià)格因素看，農(nóng)業(yè)生產(chǎn)設(shè)施條件、農(nóng)業(yè)人均GDP增加、農(nóng)村居民受教育程度與農(nóng)業(yè)面源污染的影子價(jià)格負(fù)相關(guān)即這些變量值越大，農(nóng)業(yè)面源污染物的影子價(jià)格越低；而農(nóng)業(yè)結(jié)構(gòu)中經(jīng)濟(jì)作物與糧食作物的比例、農(nóng)產(chǎn)品價(jià)格指數(shù)、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)畜禽養(yǎng)殖業(yè)產(chǎn)值占農(nóng)業(yè)總產(chǎn)值比重、人均耕地面積和降雨徑流值越高都會(huì)使污染物的影子價(jià)格提高。

減少農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)中糧食作物與經(jīng)濟(jì)作物的比例、加大農(nóng)業(yè)生產(chǎn)設(shè)施的投資力度、提高農(nóng)村居民受教育程度3個(gè)因素對(duì)提高農(nóng)業(yè)面源污染排放效率和降低面源污染物的影子價(jià)格都能起到很好的作用。因此，進(jìn)一步優(yōu)化農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)，在統(tǒng)籌考慮區(qū)域農(nóng)業(yè)資源稟賦特征，確保糧食作物安全的前提下，減少農(nóng)藥、化肥投入大的農(nóng)作物種植規(guī)模，走“兩型農(nóng)業(yè)”之路。進(jìn)一步加大對(duì)農(nóng)田水利設(shè)施、節(jié)水灌溉設(shè)施和農(nóng)業(yè)生態(tài)建設(shè)工程的投資，減少和控制農(nóng)業(yè)生產(chǎn)中的污染物排放對(duì)生態(tài)環(huán)境的破壞。農(nóng)戶是農(nóng)業(yè)面源污染的主體，農(nóng)業(yè)面源污染源頭減排離不開農(nóng)戶的積極參與，提高農(nóng)戶對(duì)農(nóng)業(yè)生態(tài)環(huán)境退化的認(rèn)知也是提高農(nóng)業(yè)面源污染的排放效率和降低農(nóng)業(yè)面源污染物影子價(jià)格的重要途徑。

（編輯：劉照勝）

參考文獻(xiàn)（References）

[1]金書秦，沈貴銀.中國(guó)農(nóng)業(yè)面源污染的困境擺脫與綠色轉(zhuǎn)型[J].改革，2013，231（5）：79-87.[Jin Shuqin， Shen Guiyin. To Get Rid of the Predicament and Green Transformation of Agricultural Nonpoint Source Pollution [J].Reform， 2013， 231（5）：79-87.]

[2]嚴(yán)素定.黃石市農(nóng)業(yè)面源污染的解析及其空間異質(zhì)性研究[J].農(nóng)業(yè)工程學(xué)報(bào)，2008，24（9）：225-228.[Yan Suding. Source Apportionment and Spatial heterogeneity of Agricultural Nonpoint Source Pollution in Huangshi， Hubei Province [J].Transactions of the CSAE， 2008，24（9）：225-228.]

[3]熊偉，王久臣，湯文志，等.三峽環(huán)庫(kù)多業(yè)共生耦合循環(huán)農(nóng)業(yè)生態(tài)系統(tǒng)初構(gòu)[J].農(nóng)業(yè)工程學(xué)報(bào)，2013，29（14）：203-209.[Xiong Wei， Wang Jiuchen， Tang Wenzhi， et al. Establishment of Integrative Circular Agroecology System for Multiple Agricultural Industries in Three Gorges Reservoir Area [J].Transactions of the Chinese Society of Agricultural Engineering，2013，29（14）：203-209.]

[4]唐學(xué)玉，張海鵬，李世平.農(nóng)業(yè)面源污染防空的經(jīng)濟(jì)價(jià)值：基于安全農(nóng)產(chǎn)品生產(chǎn)戶視角的支付意愿[J].中國(guó)農(nóng)村經(jīng)濟(jì)，2012，（3）：53-67.[Tang Xueyu， Zhang Haipeng， Li Shiping. The Economic Value of Agricultural Nonpoint Source Pollution： An Analysis Based on the Perspective of Safe Agricultural Products Production Households [J].Chinese Rural Economy， 2012，（3）：53-67.]

[5]Shephard R W. Theory of Cost and Production Functions [M].New Jersey： Princeton University Press， 1970.

[6]Chung Y H， Fre R， Grosskopf S. Productivity and Undesirable Outputs：A Directional Distance Approach [J].Journal of Environmental Management，1997，51（3）：229-24.

[7]Rolf F， Shawna G， William L W. Shadow Prices and Pollution Costs in US Agriculture[J].Ecological Economics，2006，56（1）：89-10.

[8]Karimov A. Economic Inefficiency and Shadow Prices of Inputs： The Case of Vegetable Growing Farms in Uzbekistan [J]. Procedia Economics and Finance，2013，（5）：403-412.

[9]楊俊，陳怡.基于環(huán)境因素的中國(guó)農(nóng)業(yè)生產(chǎn)率增長(zhǎng)研究[J].中國(guó)人口·資源與環(huán)境，2011，21（6）：153-157.[Yang Jun， Chen Yi. Empirical Study on Chinas Agricultural Production Growth under the Binding of Environment [J].China Population， Resources and Environment，2011，21（6）：153-157.]

[10]梁流濤.農(nóng)業(yè)發(fā)展與環(huán)境協(xié)調(diào)性評(píng)價(jià)及其影響因素分析[J].中國(guó)環(huán)境科學(xué)，2012，32（9）：1702-1708.[Liang Liutao. Coordination of Agricultural Growth with Environment [J].China Environmental Science， 2012， 32（9）：1702-1708.]

[11]潘丹，應(yīng)瑞瑤.環(huán)境污染約束下農(nóng)業(yè)生產(chǎn)率增長(zhǎng)地區(qū)差異及其動(dòng)態(tài)分布演進(jìn)[J].中國(guó)科技論壇，2013，（5）：60-67.[Pan Dan， Ying Ruiyao. Regional Disparity and Dynamic Distribution Evolution of Agricultural Productivity Growth in China under Environmental Pollution Constraint [J].Forum on Science and Technology in China， 2013， （5）：60-67.]

[12]潘丹，應(yīng)瑞瑤.中國(guó)農(nóng)業(yè)生態(tài)效率評(píng)價(jià)方法與實(shí)證：基于非期望產(chǎn)出的SBM模型分析[J].生態(tài)學(xué)報(bào)，2013，33（12）：3837-3845.[Pan Dan， Ying Ruiyao. Agricultural Ecoefficiency Evaluation in China Based on SBM Mode [J].Acta Ecologica Sinica， 2013， 33（12）：3837-3845.]

[13]閔銳，李谷成.“兩型”視角下我國(guó)糧食生產(chǎn)技術(shù)效率的空間分異[J].經(jīng)濟(jì)地理，2013，33（3）：144-149.[Min Rui， Li Gucheng. Grains Production Efficiency and Its Spatial Distribution in China from the “Twooriented” Perspective [J].Economic Geography，2013，33（3）：144-149.]

[14]Falavigna G， Manello A， Pavone S. Environmental Efficiency， Productivity and Public Funds： The Case of the Italian Agricultural Industry [J].Agricultural Systems，2013，121：73–80.

[15]李谷成，范麗霞，閔銳.資源、環(huán)境與農(nóng)業(yè)發(fā)展的協(xié)調(diào)性：基于環(huán)境規(guī)制的省級(jí)農(nóng)業(yè)環(huán)境效率排名[J].數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究，2011，（10）：21-35.[Li Gucheng， Fan Lixia， Min Rui.The Coordination of Agricultural Development with Environment and Resource [J].Journal of Quantitative & Technical Economics，2011，（10）：21-35.]

[16]陳敏鵬，陳吉寧，賴斯蕓.中國(guó)農(nóng)業(yè)和農(nóng)村污染的清單分析與空間特征識(shí)別[J].中國(guó)環(huán)境科學(xué)，2006，26（6）：751-755.[Chen Minpeng， Chen Jining， Lai Siyun. Inventory Analysis and Spatial Distribution of Chinese Agricultural and Rural Pollution [J].China Environmental Science， 2006， 26（6）：751-755.]

[17]賴斯蕓，杜鵬飛，陳吉寧.基于單元分析的非點(diǎn)源污染調(diào)查評(píng)估方法[J].清華大學(xué)學(xué)報(bào)：自然科學(xué)版，2004，44（9）：1184-1187.[Lai Siyun， Du Pengfei， Chen Jining. Evaluation of Nonpoint Source Pollution Based on Unit Analysis [J].Journal of Tsinghua University：Science and Technology Edition，2004，44（9）：1184-1187.]

[18]錢曉雍，沈根祥，郭春霞，等.基于水環(huán)境功能區(qū)劃的農(nóng)業(yè)面源污染源解析及其空間異質(zhì)性[J].農(nóng)業(yè)工程學(xué)報(bào)，2011，27（2）：103-108.[Qian Xiaoyong， Shen Genxiang， Guo Chunxia， et al. Source Apportionment and Spatial Heterogeneity of Agricultural Nonpoint Source Pollution Based on Water Environmental Function Zoning [J].Transactions of the CSAE， 2011， 27（2）：103-108.]

[19]肖新成，何丙輝，倪九派，等.農(nóng)業(yè)面源污染視角下的三峽庫(kù)區(qū)重慶段水資源的安全性評(píng)價(jià)：基于DPSIR框架的分析[J].環(huán)境科學(xué)學(xué)報(bào)，2013，33（8）：2324-2331.[Xiao Xincheng， He Binghui， Ni Jiupai， et al. Safety Assessment of Water Resources in Chongqing Section of the Three Gorges Reservoir Area：Based on DPSIR Model from the Perspective of Agricultural Nonpoint Pollution Source [J].Acta Scientiae Circumstantiae，2013，33（8）：2324-2331.]

[20]肖新成，謝德體，倪九派. 面源污染減排增匯措施下的農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)耦合狀態(tài)分析：以三峽庫(kù)區(qū)忠縣為例[J].

中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2014，22（1）：111-119. [Xiao Xincheng，Xie Deti， Ni Jiupai. Coupling State of Agricultural Ecoeconomic System Under Emission Mitigation and Sink Enhancement of Nonpoint Source Pollution：A Case Study of Zhong County in the Three Gorges Reservoir Region [J]. Chinese Journal of Ecoagriculture， 2014， 22（1）：111-119.]

[21]國(guó)涓，劉豐，王維國(guó).中國(guó)區(qū)域環(huán)境績(jī)效動(dòng)態(tài)差異及影響因素：考慮可變規(guī)模報(bào)酬和技術(shù)異質(zhì)性的研究[J].資源科學(xué)，2013，35（12）：2444-2456.[Guo Juan， Liu Feng， Wang Weiguo. Dynamic Differences and Factors Affecting Regional Environmental Performance in China： Variable Return to Scale and Technology Heterogeneity [J].Resources Science， 2013， 35（12）：2444-2456.]

[22]石風(fēng)光.中國(guó)地區(qū)環(huán)境技術(shù)效率的測(cè)算及隨機(jī)收斂性檢驗(yàn)[J].系統(tǒng)工程，2013，31（1）：61-67.[Shi Fengguang. Calculation and Stochastic Convergence Test on the Regional Environmental Technical Efficiency of China [J].Systems Engineering， 2013， 31（1）：61-67.]

[23]Lozano S，Gutierrez E. Nonparametric Frontier Approach to Modeling the Relationships among Population GDP， Energy Consumption and CO2 Emissions [J].Ecol.Econ，2008，66（4）：687-699.

[24]Fare， Grosskopf. A Comment on Weak Disposability in Nonparametric Production Analysis [J].American Journal of Agricultural Economics， 2009， 91（2）：535-538.

[25]Andre， Herrero， Riesgo. A Modified DEA Model to Estimate the Importance of Objectives with an Application to Agricultural Economics [J].Omega，2010，38：371-3822.

[26]Murty M N， Kumar S， Dhavala K K. Measuring Environmental Efficiency of Industry： A Case Study of Thermal Power Generation in India [J].Environmental & Resource Economics， 2007， 38（1）：31-50.

[27]袁鵬，程施.我國(guó)工業(yè)污染物的影子價(jià)格估計(jì)[J].統(tǒng)計(jì)研究，2011，28（9）：66-73.[ Yuan Peng， Cheng Shi. Estimating Shadow Pricing of Industrial Pollutions in China [J].Statistical Research，2011，28（9）：66-73.]

[28]梁流濤，曲福田，馮淑怡.基于環(huán)境污染約束視角的農(nóng)業(yè)技術(shù)效率測(cè)度[J].自然資源學(xué)報(bào)，2012，27（9）：1580-1588.[Liang Liutao， Qu Futian， Feng Shuyi. Agricultural Technical Efficiency Measurement under the Environmental Constraints [J].Journal of Natural Resources， 2012， 27（9）：1580-1588.]

[29]石嫣，程存旺，朱藝，等.中國(guó)農(nóng)業(yè)源污染防治的制度創(chuàng)新與組織創(chuàng)新：兼析第一次全國(guó)污染源普查公報(bào)[J].農(nóng)業(yè)經(jīng)濟(jì)與管理，2011，12（2）：27-31.[Shi Yan， Cheng Cunwang， Zhu Yi， et al. Institutional and Organizational Innovation on China Agricultural Nonpoint Pollution Prevention： Analysis on the 1th National Survey of Pollution Sources Bulletin [J].Agricultural Economics and Management， 2011， 23（2）：27-31.]

[30]李曉紅，韓勇，鄭陽(yáng)華.三峽庫(kù)區(qū)坡耕地土壤侵蝕治理效益分析[J].