1978—2017年中國農(nóng)業(yè)污染物影子價(jià)格及污染成本測算

鄒利林，劉彥隨，王永生

(1.華僑大學(xué)政治與公共管理學(xué)院/政治發(fā)展與公共治理研究中心，泉州362021；2.中國科學(xué)院地理科學(xué)與資源研究所，北京100101)

0 引 言

改革開放以來，中國糧食產(chǎn)量由1978年30 476.5萬t增長至2017年66 160.7萬t，同期畜禽肉類產(chǎn)量由856.3萬t增長至8 654.6萬t。與此同時(shí)，化肥施用量從1978年884.0萬t增長至2017年5 859.4萬t，農(nóng)用塑料薄膜與農(nóng)藥使用量分別從1990年48.2萬t與73.3萬t增長至2017年252.8萬t與165.5萬t。由此可見，糧食產(chǎn)量與畜禽肉產(chǎn)品的增長主要依賴于生產(chǎn)要素的密集投入[1-2]，而正是這種“高消耗、高污染、高投入”的生產(chǎn)方式造成了嚴(yán)重的農(nóng)業(yè)污染問題[3-4]。農(nóng)業(yè)污染主要指農(nóng)業(yè)生產(chǎn)過程中，溶解的或固體的污染物通過水循環(huán)引起水質(zhì)污染的過程，具有來源分散性、擴(kuò)散隨機(jī)性、污染異質(zhì)性和治理低效性等特點(diǎn)[5-6]。在全面建設(shè)資源節(jié)約型與環(huán)境友好型社會的背景下，如何降低農(nóng)業(yè)污染排放，核算農(nóng)業(yè)污染治理成本，建立農(nóng)業(yè)污染管理制度[7-9]，成為國家推進(jìn)鄉(xiāng)村振興戰(zhàn)略必須解決的科學(xué)難題。為此，在制定農(nóng)業(yè)污染減排政策及衡量農(nóng)村發(fā)展水平時(shí)有必要將農(nóng)業(yè)污染排放納入到考核體系，以便全面反映農(nóng)村社會經(jīng)濟(jì)發(fā)展所帶來的環(huán)境成本。

相較于工業(yè)污染排放的排污許可制度而言，農(nóng)業(yè)污染排放至今仍未建立相應(yīng)的交易市場，污染物價(jià)格也無法通過市場機(jī)制得以反映，導(dǎo)致農(nóng)業(yè)污染排放產(chǎn)生的環(huán)境成本難以估計(jì)[7,10]。影子價(jià)格作為供求市場完全均衡時(shí)污染物真實(shí)經(jīng)濟(jì)價(jià)值的體現(xiàn)，已成為測算資源環(huán)境產(chǎn)品交易價(jià)格的主要依據(jù)[11-13]。Shephard較早利用投入與產(chǎn)出距離函數(shù)，對價(jià)格信息缺失下的污染物影子價(jià)格進(jìn)行了估算[14]。由于該方法無法克服投入與產(chǎn)出同向收縮或擴(kuò)張產(chǎn)生的隨機(jī)誤差，F(xiàn)?re等采用二次型產(chǎn)出方向性距離函數(shù)和利潤函數(shù)的對偶性對該函數(shù)模型進(jìn)行了優(yōu)化[15]，并被廣泛用于C、CO2、SO2與NOx等污染物影子價(jià)格的測算[16-18]。為破解因農(nóng)業(yè)污染物價(jià)格信息缺失造成環(huán)境成本估算依據(jù)不足的困境，汪慧玲等采用該方法測算2000—2012年中國省級農(nóng)業(yè)生產(chǎn)過程中排放的總氮、總磷和化學(xué)需氧量的影子價(jià)格[19]，這為本研究奠定了良好的基礎(chǔ)。此外，這一指標(biāo)還可以衡量“金山銀山”與“綠水青山”之間的替代關(guān)系[20-21]，在農(nóng)業(yè)污染排放交易缺失的情況下，可為政府制定農(nóng)業(yè)污染減排政策提供決策參考。

近年來，由于農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)與農(nóng)村經(jīng)濟(jì)結(jié)構(gòu)的轉(zhuǎn)變，農(nóng)藥化肥的過度使用、畜禽養(yǎng)殖量的不斷增加、農(nóng)田秸稈的不合理利用以及水產(chǎn)養(yǎng)殖規(guī)模的持續(xù)擴(kuò)大，農(nóng)

業(yè)污染已經(jīng)超過工業(yè)污染成為水污染的主要來源[22-24]。但相應(yīng)的治理政策因?yàn)閮r(jià)格信息的缺失難以出臺，導(dǎo)致農(nóng)業(yè)污染治理無法像工業(yè)污染治理那樣實(shí)現(xiàn)從命令式調(diào)控向市場式調(diào)控轉(zhuǎn)變[20,25]。因此，測算農(nóng)業(yè)污染物的影子價(jià)格成為推進(jìn)農(nóng)業(yè)環(huán)境政策制定的關(guān)鍵[26-27]。盡管已有研究對不同時(shí)期與不同尺度下農(nóng)業(yè)污染物影子價(jià)格進(jìn)行了測算，但對全國尺度下長時(shí)間序列農(nóng)業(yè)污染物影子價(jià)格的關(guān)注仍然較少[7,19]，尤其是缺少從制度環(huán)境變革的角度對其長期演變特征展開深入分析。因此，本文采用F?re等提出的二次型方向性距離函數(shù)和收益函數(shù)的對偶關(guān)系[15]，對1978—2017年中國省級農(nóng)業(yè)排放污染物總氮（TN）、總磷（TP）和化學(xué)需氧量（COD）的影子價(jià)格進(jìn)行估算，并測算農(nóng)業(yè)生產(chǎn)過程中污染物的邊際減排成本，以及指導(dǎo)環(huán)境保護(hù)政策的制定和綠色增長考核體系的建立。

1 研究方法

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

隨著農(nóng)業(yè)生產(chǎn)要素投入的增加，農(nóng)業(yè)期望產(chǎn)出成倍增長，與此同時(shí)非期望產(chǎn)出也逐漸增加，由此引發(fā)的農(nóng)業(yè)污染成為社會必須關(guān)注的問題。農(nóng)業(yè)生產(chǎn)投入、期望產(chǎn)出與非期望產(chǎn)出之間的量化關(guān)系被稱為環(huán)境技術(shù)效率。環(huán)境技術(shù)將生產(chǎn)部門要素投入、期望產(chǎn)出、非期望產(chǎn) 出 的 集 合 分 別 定 義 為 x=(x1,x2…,xN)∈，y=(y1,y2…,yM)∈，b=(b1,b2…,bJ)∈。因此，環(huán)境技術(shù)的生產(chǎn)集P(x),x∈可表示為

方向性距離函數(shù)建立在環(huán)境技術(shù)假設(shè)的前提下，根據(jù)F?re等[28-29]對方向性距離的描述，引入方向性向量g=(gy,gb)，且g≠0，用β表征方向向量的大小，則方向性距離函數(shù)可以定義為

進(jìn)一步在一次性假設(shè)下，生產(chǎn)集P(x),x∈RN+可以從方向性距離函數(shù)中推導(dǎo)，即：因此，環(huán)境技術(shù)可以用P(x)表示，或者說等于基于這一定義，方向性距離函數(shù)具有如下的轉(zhuǎn)移特征：

式中α為方向向量前進(jìn)距離。

進(jìn)一步令方向向量為(y,-b)時(shí)，方向性距離函數(shù)可以表示為

Do(x,y,z)為Shephard距離函數(shù)，式（5）揭示了方向性距離函數(shù)和Shephard距離函數(shù)的關(guān)系。結(jié)合環(huán)境技術(shù)的性質(zhì)，方向性距離函數(shù)滿足以下約束條件，具體形

式如下：

式中x′，y′分別為x與y的條件值，θ為約束條件。

1.2 影子價(jià)格模型構(gòu)建

本文從方向性距離函數(shù)與收益函數(shù)之間的對偶關(guān)系中推導(dǎo)出污染物的影子價(jià)格。為便于分析，假設(shè)p=(p1,p2…,pM)∈RM+與q=(q1,q2…,qJ)∈ RJ+分別為期望產(chǎn)出和非期望產(chǎn)出的價(jià)格。由于非期望產(chǎn)出是農(nóng)業(yè)污染物，其對收益的貢獻(xiàn)為負(fù)，因此收益函數(shù)被定義為式（7）：

R(x,p,q)=maxy,b{ }py-qb:(y,b)∈ P(x) （7）

收益函數(shù)R(x,p,q)給出了當(dāng)期望產(chǎn)出價(jià)格為p和非期望產(chǎn)出價(jià)格為q時(shí)，投入集x可能估算出的最大收益。利用式（3），最大的收益函數(shù)等同于

如果農(nóng)業(yè)部門產(chǎn)出集合(y,b)∈P(x)，那么存在：

也就是說，如果產(chǎn)出向量(y,b)是可行的，那么通過向 方向移動來消除與該產(chǎn)出向量相關(guān)聯(lián)的任何無效性也是可行的。因此，給定g=(gy,gb)，可以得出

由此可以利用收益函數(shù)將式（2）中的方向性距離函數(shù)推導(dǎo)出來：

對式（11）運(yùn)用兩次包絡(luò)定理，就可以得到影子價(jià)格模型式（12）和（13）。此時(shí)，假設(shè)知道第m種期望產(chǎn)出的價(jià)格pm，則可以通過式（14）求出第j種非期望產(chǎn)出的價(jià)格qj。

1.3 二次型函數(shù)

由于參數(shù)距離函數(shù)具有更好的微分性質(zhì)，因此本文采用二次參數(shù)形式的距離函數(shù)求解農(nóng)業(yè)污染物的影子價(jià)格。假設(shè)方向向量g=(1,1)，假設(shè)在t=1,…,T時(shí)期有k=1,…,K生產(chǎn)部門，那么k生產(chǎn)部門在t時(shí)期的二次型方向性距離函數(shù)形式為

上式中，αn，βm，γj，αnn′，βmm′，γjj′，δnm，ηnj，μmj為回歸系數(shù)，代表各二次型單元的偏效應(yīng)。交叉產(chǎn)出與投入效果 的 對 稱 性 要 求 αnn′= αn′n,n≠n′； βmm′

= βm′m,m ≠ m′； γjj′

為滿足生產(chǎn)部門每一個(gè)時(shí)期的生產(chǎn)都盡可能在生產(chǎn)前沿上，本文在求解影子價(jià)格時(shí)選取式（15）中參數(shù)使得目標(biāo)函數(shù)最小。該目標(biāo)函數(shù)必須滿足：

條件（i）表示各生產(chǎn)部門在不同時(shí)期的投入—產(chǎn)出向量是可行的；條件（ⅱ）和（ⅲ）是對式（16）和（17）的單調(diào)性約束；條件（ⅳ）是對式（4）轉(zhuǎn)移性質(zhì)的約束；條件（v）體現(xiàn)了投入和產(chǎn)出向量之間的對稱性。

1.4 農(nóng)業(yè)污染成本測算

污染物影子價(jià)格可以衡量農(nóng)業(yè)生產(chǎn)的污染成本，是指導(dǎo)環(huán)境保護(hù)政策和綠色增長核算的重要參考。借助已測算出的3種污染物影子價(jià)格，可以對農(nóng)業(yè)生產(chǎn)過程的污染成本進(jìn)行分析。計(jì)算公式如下：

式中Cp為農(nóng)業(yè)污染總成本，萬元；CTN、CTP、CCOD分別為TN、TP、COD的污染成本，萬元；qTN、qTP、qCOD分別為TN、TP、COD的影子價(jià)格，萬元/t；bTN、bTP、bCOD分別為TN、TP、COD的排放量，t。

2 指標(biāo)體系與數(shù)據(jù)來源

2.1 指標(biāo)體系

農(nóng)業(yè)生產(chǎn)過程中由于生產(chǎn)要素的投入不僅會產(chǎn)生期望產(chǎn)出，同時(shí)也會排放各類污染物，即產(chǎn)生非期望產(chǎn)出。為測算中國農(nóng)業(yè)生產(chǎn)過程中污染物TN、TP、COD的影子價(jià)格，在構(gòu)建農(nóng)業(yè)生產(chǎn)投入產(chǎn)出指標(biāo)體系時(shí)應(yīng)當(dāng)綜合考慮資源節(jié)約、環(huán)境友好和農(nóng)業(yè)經(jīng)濟(jì)增長三者間的統(tǒng)籌協(xié)調(diào)發(fā)展[30]。由于本研究關(guān)注的是農(nóng)產(chǎn)品直接生產(chǎn)過程，而非加工過程，因此本研究選取農(nóng)林牧漁業(yè)總產(chǎn)值為期望產(chǎn)出指標(biāo)，選擇農(nóng)田化肥、畜禽養(yǎng)殖、水產(chǎn)養(yǎng)殖、農(nóng)田秸稈4類非點(diǎn)源農(nóng)業(yè)污染的TN、TP、COD排放量為非期望產(chǎn)出指標(biāo)。參照已有文獻(xiàn)及依據(jù)典型性、科學(xué)性、可獲取性等原則，選擇土地、勞動力、資本為投入要素指標(biāo)，分別用農(nóng)作物播種面積（hm2）、農(nóng)林牧漁業(yè)從業(yè)人員（萬人）、農(nóng)業(yè)資本存量（萬元）表示。

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

本研究選取除香港、臺灣、澳門及西藏、青海（由于西藏、青海的農(nóng)業(yè)生產(chǎn)模式與其他省市差異較大，暫不分析）以外的29個(gè)省市自治區(qū)1978—2017年期間1 131條農(nóng)業(yè)投入產(chǎn)出數(shù)據(jù)為研究樣本（1978—1996年重慶市數(shù)據(jù)包含在四川省內(nèi)，1978—1987年海南省數(shù)據(jù)包含在廣東省內(nèi)）。農(nóng)林牧漁業(yè)產(chǎn)值主要來源于《中國農(nóng)村統(tǒng)計(jì)年鑒》（1985—2018年）及各省《統(tǒng)計(jì)年鑒》（1979—1984年），該值以1978年為基期進(jìn)行價(jià)格平減。由于現(xiàn)有統(tǒng)計(jì)年鑒中缺少農(nóng)業(yè)污染物排放總量統(tǒng)計(jì)數(shù)據(jù)，本研究從《中國農(nóng)村統(tǒng)計(jì)年鑒》、《改革開放三十年農(nóng)業(yè)統(tǒng)計(jì)資料匯編》、《新中國農(nóng)業(yè)60年統(tǒng)計(jì)資料》、《中國統(tǒng)計(jì)年鑒》、《中國漁業(yè)統(tǒng)計(jì)年鑒》、《中國海洋統(tǒng)計(jì)年鑒》、《中國農(nóng)業(yè)統(tǒng)計(jì)資料匯編》等統(tǒng)計(jì)資料中獲取化肥折純施用量、畜禽出欄數(shù)與存欄數(shù)、水產(chǎn)品產(chǎn)量、農(nóng)作物產(chǎn)量等數(shù)值，采用賴斯蕓提出的清單分析法計(jì)算得到農(nóng)業(yè)污染物TN、TP、COD的排放量[31]。投入要素指標(biāo)數(shù)據(jù)也主要來自上述統(tǒng)計(jì)資料，其中農(nóng)業(yè)資本存量數(shù)據(jù)首先通過徐現(xiàn)祥等[32]的縮減指數(shù)構(gòu)造方法計(jì)算得出農(nóng)業(yè)投資數(shù)據(jù)的價(jià)格平減指數(shù)，然后借鑒宗振利等[33]計(jì)算的分省折舊率和基期資本存量數(shù)據(jù)，根據(jù)永續(xù)盤存法計(jì)算得出各省數(shù)值。

3 結(jié)果與分析

3.1 中國農(nóng)業(yè)污染物影子價(jià)格的變化趨勢

為檢驗(yàn)估計(jì)參數(shù)結(jié)果是否具有良好的適應(yīng)性，根據(jù)參數(shù)估計(jì)結(jié)果驗(yàn)證樣本是否滿足零點(diǎn)關(guān)聯(lián)性，即驗(yàn)證非期望產(chǎn)出與期望產(chǎn)出是否同步產(chǎn)生[15]。經(jīng)檢驗(yàn)發(fā)現(xiàn)，在所研究的1 131個(gè)樣本中，一共有1 078個(gè)樣本滿足零點(diǎn)關(guān)聯(lián)性，占比達(dá)到95.31%，說明絕大多數(shù)的樣本都滿足模型的假設(shè)。從估算結(jié)果來看，1978—2017年中國農(nóng)業(yè)污染物的影子價(jià)格總體呈下降趨勢，TN與TP的影子價(jià)格隨著農(nóng)業(yè)污染治理力度加強(qiáng)而呈現(xiàn)平緩下降趨勢，而COD影子價(jià)格由于畜禽養(yǎng)殖規(guī)模的變化而呈現(xiàn)較大的波動性。具體來看，TN、TP、COD的影子價(jià)格分別從7 941.5、8 990.1、2 087.6元/t，下降到5 811.9、5 759.5、841.9元/t，分別下降了26.8%、35.9%、59.7%，表明中國對農(nóng)業(yè)污染問題日益重視，農(nóng)業(yè)生產(chǎn)的機(jī)會成本不斷降低。從影子價(jià)格的變化趨勢來看，具有較為明顯的階段性，并大致經(jīng)歷了4個(gè)階段（圖1）。1978—1985年為自由發(fā)展階段，盡管國家實(shí)行改革開放及推進(jìn)家庭聯(lián)產(chǎn)承包責(zé)任制，但由于市場機(jī)制與國家政策的缺失，農(nóng)業(yè)發(fā)展緩慢，農(nóng)業(yè)污染物的影子價(jià)格基本保持不變。1986—1995年為改革推動階段，國家將主要農(nóng)副產(chǎn)品的統(tǒng)購統(tǒng)派制度逐步改革為以計(jì)劃為主與市場調(diào)節(jié)為輔的制度，充分調(diào)動了廣大農(nóng)業(yè)生產(chǎn)者的積極性，農(nóng)業(yè)進(jìn)入快速發(fā)展時(shí)期，在農(nóng)產(chǎn)品市場化改革的推動下農(nóng)業(yè)污染物的影子價(jià)格大幅度降低。1996—2006年為市場調(diào)控階段，農(nóng)業(yè)政策以保護(hù)農(nóng)業(yè)生產(chǎn)、支持農(nóng)民增收、促進(jìn)農(nóng)村發(fā)展及“取消農(nóng)業(yè)稅、工業(yè)反哺農(nóng)業(yè)”為主要特征，在市場機(jī)制的自我調(diào)控下農(nóng)業(yè)污染物的影子價(jià)格輕微變動。2007—2017年為政策激勵階段，中央一號文連續(xù)14 a聚焦“三農(nóng)”問題，特別是增加農(nóng)業(yè)補(bǔ)貼、實(shí)現(xiàn)農(nóng)業(yè)現(xiàn)代化、推進(jìn)農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革等系列農(nóng)業(yè)新政的出臺，激活了農(nóng)村與農(nóng)業(yè)發(fā)展的內(nèi)在活力，農(nóng)業(yè)污染物的影子價(jià)格再次緩慢下降。

圖1 1978—2017年中國農(nóng)業(yè)污染物影子價(jià)格的階段性變化趨勢Fig.1 Periodic change trends of shadow prices of agricultural pollutants in China from 1978 to 2017

3.2 中國省級農(nóng)業(yè)污染物影子價(jià)格的差異分析

為探索農(nóng)業(yè)污染物影子價(jià)格的區(qū)域差異，本文分區(qū)方案采用《中國自然地理》提出的七大分區(qū)，并參照劉彥隨等[34]提出的現(xiàn)代農(nóng)業(yè)區(qū)劃方案適當(dāng)調(diào)整。從地理分區(qū)尺度來看，農(nóng)業(yè)污染物的影子價(jià)格存在明顯的區(qū)域差異。圖2顯示，東北地區(qū)TN的影子價(jià)格最低為5 511.9元/t，而西南地區(qū)TN的影子價(jià)格最高為8 430.4元/t；華東地區(qū)TP的影子價(jià)格最低為6 406.0元/t，而西南地區(qū)TP的影子價(jià)格最高為8 821.4元/t；西南地區(qū)COD的影子價(jià)格最低為570.8元/t，而華東地區(qū)COD的影子價(jià)格最高為1 727.2元/t。究其原因在于，東北地區(qū)作為國家糧食主產(chǎn)區(qū)，優(yōu)越的農(nóng)業(yè)資源稟賦使得農(nóng)業(yè)污染排放效率較高，農(nóng)業(yè)機(jī)械化水平也相對更高，因此農(nóng)業(yè)污染物的影子價(jià)格低于全國平均水平。華北與華中地區(qū)地勢平坦、水土條件優(yōu)沃且勞動力充足，農(nóng)業(yè)生產(chǎn)水平與農(nóng)業(yè)現(xiàn)代化進(jìn)程基本同步，農(nóng)業(yè)污染物的影子價(jià)格與全國平均水平基本一致。華東地區(qū)經(jīng)濟(jì)較發(fā)達(dá)，農(nóng)業(yè)種植規(guī)模較小，農(nóng)業(yè)投入強(qiáng)度較低，TN與TP影子價(jià)格較低，但由于畜禽養(yǎng)殖產(chǎn)值占農(nóng)業(yè)產(chǎn)值的比例較大，COD的影子價(jià)格較高。華南地區(qū)的農(nóng)業(yè)氣候資源及社會經(jīng)濟(jì)環(huán)境差異較大，導(dǎo)致區(qū)域農(nóng)業(yè)污染物的影子價(jià)格也存在較大差異，但總體低于全國平均水平。西北地區(qū)農(nóng)業(yè)生產(chǎn)的自然條件與要素投入水平較差，農(nóng)業(yè)生產(chǎn)方式較粗放，農(nóng)業(yè)污染物的影子價(jià)格較高。西南地區(qū)由于人均耕地規(guī)模較少，土地細(xì)碎化程度較高，農(nóng)業(yè)生產(chǎn)效率相對較低，TN與TP的影子價(jià)格較高，但由于其畜禽養(yǎng)殖產(chǎn)值占農(nóng)業(yè)產(chǎn)值的比重較小，COD的影子價(jià)格較低。

圖2 1978—2017年中國七大農(nóng)業(yè)分區(qū)的農(nóng)業(yè)污染物平均影子價(jià)格Fig.2 Average shadow prices of agricultural pollutants in seven agricultural regions of China from 1978 to 2017

從省級尺度來看，由于農(nóng)業(yè)生產(chǎn)條件、農(nóng)業(yè)種植強(qiáng)度、農(nóng)業(yè)經(jīng)濟(jì)構(gòu)成、產(chǎn)業(yè)結(jié)構(gòu)等因素的差異性，使得不同省份的農(nóng)業(yè)污染物影子價(jià)格差異明顯。廣西、四川、云南、貴州等以喀斯特地貌為主的省份水土流失較嚴(yán)重，土層較薄且滲水性較強(qiáng)，農(nóng)業(yè)種植污染治理難度較大，TN與TP的影子價(jià)格較高，而北京、上海、廣東、浙江等省市農(nóng)業(yè)經(jīng)濟(jì)在整個(gè)國民經(jīng)濟(jì)中占比較低，農(nóng)業(yè)污染減排措施較嚴(yán)，TN與TP影子價(jià)格也較低；河南、山東、河北等農(nóng)業(yè)大省由于復(fù)種、間作、規(guī)?；笄蒺B(yǎng)殖、過量施用肥料等集約化農(nóng)業(yè)措施的實(shí)施，農(nóng)業(yè)污染排放較高，COD影子價(jià)格較高；而云南、內(nèi)蒙古、廣西、重慶、貴州等省市為促進(jìn)旅游業(yè)發(fā)展而強(qiáng)化生態(tài)環(huán)境的保護(hù)，畜禽養(yǎng)殖規(guī)模較小且排放管制較嚴(yán)，COD影子價(jià)格較低（表1）。

表1 1978—2017年中國省級農(nóng)業(yè)污染物的平均影子價(jià)格Table 1 Average shadow prices of agricultural pollutants at provincial level in China from 1978 to 2017

3.3 中國省級農(nóng)業(yè)污染成本分析

1978—2017年中國農(nóng)業(yè)生產(chǎn)年均污染成本為760.6億元，占全國年均農(nóng)業(yè)總產(chǎn)值的10.8%，幾乎相當(dāng)于東北三省的農(nóng)業(yè)總產(chǎn)值之和，表明農(nóng)業(yè)污染成本巨大。從省市層面來看，海南、重慶、黑龍江、新疆的污染成本占各省年均農(nóng)業(yè)總產(chǎn)值的比例低于6%，而寧夏、貴州、河北和山東的污染成本占各省年均農(nóng)業(yè)總產(chǎn)值的比例高于14%。特別是山東作為全國農(nóng)業(yè)規(guī)模和農(nóng)業(yè)產(chǎn)值最大的省份，其污染成本占比最高為16.5%（圖3）。這表明，盡管農(nóng)業(yè)生產(chǎn)規(guī)模是決定污染成本的重要因素（如山東、河南），但不同的地理空間格局、農(nóng)業(yè)資源稟賦、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)與農(nóng)業(yè)管理措施對農(nóng)業(yè)污染成本也具有重要影響。因此，國家可以通過制定一些政策將與農(nóng)業(yè)生產(chǎn)有關(guān)的外部成本內(nèi)部化，以平衡農(nóng)業(yè)污染治理的省際差異。譬如說，在污染成本較高的省份對化肥、農(nóng)藥、地膜等農(nóng)業(yè)投入品加征稅收，以激勵農(nóng)業(yè)部門采用新技術(shù)、新產(chǎn)品減少污染排放，同時(shí)對污染成本較低省份的農(nóng)業(yè)部門給予補(bǔ)貼，以彌補(bǔ)這些地方因農(nóng)業(yè)污染減排而造成的損失。

圖3 1978—2017年中國省級農(nóng)業(yè)污染成本占農(nóng)業(yè)總產(chǎn)值的比例Fig.3 Proportion of provincial agricultural pollution cost in total agricultural output value in China from 1978 to 2017

4 結(jié)論與討論

本文以1978—2017年中國農(nóng)業(yè)部門的省級面板數(shù)據(jù)為研究單元，利用二次型方向性距離函數(shù)和收益函數(shù)的對偶關(guān)系估算了農(nóng)業(yè)污染物的影子價(jià)格，并測算了農(nóng)業(yè)污染成本。結(jié)論如下：

1）1978—2017年，中國農(nóng)業(yè)污染物TN、TP、COD的影子價(jià)格總體呈下降趨勢，分別下降了26.8%、35.9%、59.7%，并大致經(jīng)歷了自由發(fā)展、改革推動、市場調(diào)控與政策激勵4個(gè)階段。這既說明中國對農(nóng)業(yè)污染問題日益重視，農(nóng)業(yè)生產(chǎn)的機(jī)會成本不斷降低，也表明制度變革與政策演替對農(nóng)業(yè)生產(chǎn)效率及污染物影子價(jià)格具有重要的激勵與約束作用。

2）農(nóng)業(yè)污染物的影子價(jià)格存在明顯的區(qū)域差異，東北地區(qū)農(nóng)業(yè)污染物的影子價(jià)格明顯低于國家平均水平；西南地區(qū)TN與TP的影子價(jià)格明顯高于國家平均水平，而COD的影子價(jià)格明顯低于國家平均水平；華東地區(qū)TP的影子價(jià)格明顯低于國家平均水平，而COD的影子價(jià)格明顯高于國家平均水平。農(nóng)業(yè)污染物影子價(jià)格的區(qū)域差異表明，農(nóng)業(yè)既受到光、熱、水、土等自然要素的直接影響，又受到不同經(jīng)濟(jì)社會發(fā)展條件下農(nóng)業(yè)生產(chǎn)力水平、產(chǎn)業(yè)政策、區(qū)域經(jīng)濟(jì)水平與農(nóng)業(yè)管理水平的綜合影響。

3）1978—2017年，中國農(nóng)業(yè)生產(chǎn)年均污染成本為760.6億元，占全國年均農(nóng)業(yè)總產(chǎn)值的10.8%，表明農(nóng)業(yè)污染成本巨大。從省市層面來看，海南、重慶、黑龍江、新疆的污染成本占各省年均農(nóng)業(yè)總產(chǎn)值的比例低于6%，而寧夏、貴州、河北、山東高于14%。這說明盡管農(nóng)業(yè)生產(chǎn)規(guī)模是決定污染成本的主要因素，但不同的地理空間格局、農(nóng)業(yè)資源稟賦、農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)與農(nóng)業(yè)管理措施對農(nóng)業(yè)污染成本也具有重要影響。為此，國家可以制定一些差異化的稅收政策來平衡農(nóng)業(yè)污染治理的省際差異。

中國農(nóng)業(yè)污染物影子價(jià)格的長期演變趨勢及其污染成本說明，農(nóng)業(yè)生產(chǎn)是一個(gè)受自然、社會、經(jīng)濟(jì)與政策等因素綜合影響的復(fù)雜問題。農(nóng)業(yè)政策的制定要按照“兩型農(nóng)業(yè)”發(fā)展的要求，通過農(nóng)業(yè)供給側(cè)結(jié)構(gòu)性改革提高農(nóng)業(yè)供給質(zhì)量和效率，并避免掉入“逐利政策”陷阱。農(nóng)業(yè)發(fā)展方式的轉(zhuǎn)變要以提高土地產(chǎn)出率、資源利用率與勞動生產(chǎn)率為核心，降低對石化農(nóng)業(yè)的依賴與減少農(nóng)業(yè)面源污染物的排放，大力培養(yǎng)農(nóng)業(yè)主體的資源節(jié)約和環(huán)境保護(hù)觀念，探索循環(huán)農(nóng)業(yè)、生態(tài)農(nóng)業(yè)、集約農(nóng)業(yè)等新理念與新模式。