楊 芳，吳永常，韋文珊

(1．中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與區(qū)劃研究所，北京 100081； 2.中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)經(jīng)濟(jì)與發(fā)展研究所，北京 100081)

0 引言

農(nóng)業(yè)可持續(xù)發(fā)展是持續(xù)發(fā)展戰(zhàn)略的重要組成部分，是生態(tài)文明建設(shè)的重要內(nèi)容，是實(shí)現(xiàn)“五位一體”戰(zhàn)略布局、建設(shè)美麗中國(guó)的必然選擇，是中國(guó)特色新型農(nóng)業(yè)現(xiàn)代化道路的內(nèi)在要求[1]。2015年5月28日，國(guó)家農(nóng)業(yè)部聯(lián)合其他7部委頒布了《全國(guó)農(nóng)業(yè)可持續(xù)發(fā)展規(guī)劃(2015-2030)》[2]。2016年12月國(guó)務(wù)院發(fā)布了《中國(guó)落實(shí)2030年可持續(xù)發(fā)展議程創(chuàng)新示范區(qū)建設(shè)方案》，方案提出將在現(xiàn)有的國(guó)家可持續(xù)發(fā)展實(shí)驗(yàn)區(qū)工作的基礎(chǔ)上，建設(shè)一批落實(shí)2030年可持續(xù)發(fā)展議程創(chuàng)新的示范區(qū)[3-4]，科學(xué)地監(jiān)測(cè)和評(píng)價(jià)農(nóng)業(yè)可持續(xù)發(fā)展水平是其中的重要內(nèi)容。

關(guān)于農(nóng)業(yè)可持續(xù)發(fā)展評(píng)價(jià)的方法有很多，大致包括壓力—狀態(tài)—反應(yīng)框架[5]，能值分析[6-8]、生態(tài)足跡[9]、層次分析法[10]、熵指法[11]、神經(jīng)網(wǎng)絡(luò)模型[12-13]、綜合指數(shù)評(píng)價(jià)和多元統(tǒng)計(jì)方法[14]、數(shù)據(jù)包絡(luò)分析方法等[15]。這些評(píng)價(jià)方法主要以指標(biāo)評(píng)價(jià)法的研究為主，由“主題—次主題—指標(biāo)—指標(biāo)群”構(gòu)成，再用某種方法對(duì)眾多指標(biāo)集成，以求對(duì)可持續(xù)發(fā)展進(jìn)行綜合的評(píng)價(jià)，這種評(píng)價(jià)體系雖能夠較好地反映出一個(gè)地區(qū)的可持續(xù)發(fā)展情況，但也存在著指標(biāo)多、操作困難、數(shù)據(jù)不易獲取等問(wèn)題，缺乏通用的標(biāo)準(zhǔn)。農(nóng)業(yè)可持續(xù)發(fā)展，涉及“農(nóng)村—農(nóng)業(yè)—產(chǎn)業(yè)—產(chǎn)品消費(fèi)”各個(gè)環(huán)節(jié)，目前以解決主要矛盾為主，重點(diǎn)是區(qū)域農(nóng)業(yè)生產(chǎn)方式的可持續(xù)問(wèn)題。文章從產(chǎn)業(yè)結(jié)構(gòu)、能量平衡關(guān)系、生產(chǎn)方式可持續(xù)性等角度，根據(jù)農(nóng)業(yè)生產(chǎn)全要素投入和產(chǎn)出的關(guān)系，基于能量分析的方法，探索構(gòu)建一種簡(jiǎn)便易行的評(píng)價(jià)方法：農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)(Agriculture Sustainable Development Index，ASDI)。

圖1 修正的農(nóng)業(yè)能量投入產(chǎn)出曲線(xiàn)

圖2 農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)評(píng)價(jià)方法

生產(chǎn)本質(zhì)上是一個(gè)能量與物質(zhì)的轉(zhuǎn)換過(guò)程，從資源經(jīng)濟(jì)學(xué)角度來(lái)看，農(nóng)業(yè)是一個(gè)典型的農(nóng)業(yè)經(jīng)濟(jì)生產(chǎn)系統(tǒng)，投入各種資源和要素，生產(chǎn)出各類(lèi)具有價(jià)值的農(nóng)產(chǎn)品[16]。農(nóng)業(yè)投入和產(chǎn)出有實(shí)物、價(jià)值和勞務(wù)消耗等不同的形式，在生態(tài)環(huán)境的分析中具有很大的局限性[17]。將生態(tài)經(jīng)濟(jì)系統(tǒng)內(nèi)的不同投入和不同產(chǎn)出轉(zhuǎn)換為統(tǒng)一標(biāo)準(zhǔn)的能量單位，應(yīng)用能量這一概念和度量標(biāo)準(zhǔn)，克服了經(jīng)濟(jì)分析的主觀(guān)性和片面性，實(shí)現(xiàn)農(nóng)業(yè)復(fù)雜系統(tǒng)評(píng)價(jià)的一致性標(biāo)準(zhǔn)[18]。之后，有很多文獻(xiàn)利用能量從投入產(chǎn)出的角度進(jìn)行分析[19-24]，對(duì)農(nóng)業(yè)生態(tài)系統(tǒng)的能量投入效率進(jìn)行研究[25-27]。很多研究發(fā)現(xiàn)，農(nóng)業(yè)生態(tài)系統(tǒng)的投入存在類(lèi)似于邊際報(bào)酬遞減的現(xiàn)象，生態(tài)系統(tǒng)投入與產(chǎn)出之間呈S型的關(guān)系。隨著能量投入逐漸增加達(dá)到一定程度，生產(chǎn)系統(tǒng)的能量產(chǎn)出就會(huì)出現(xiàn)下降[28-33]。尹鈞等[34]于1999年在山西的萬(wàn)榮縣和山西農(nóng)大的試驗(yàn)農(nóng)場(chǎng)，以小麥的投入能和產(chǎn)出能作為研究對(duì)象進(jìn)行農(nóng)田人工輔助能各項(xiàng)投入能的投入產(chǎn)出規(guī)律研究，得出各項(xiàng)人工輔助能的能量投入與能量產(chǎn)出都遵循一元二次拋物線(xiàn)的變化規(guī)律。2012年，尹鋼吉等[35]在對(duì)施肥水平與產(chǎn)量關(guān)系數(shù)學(xué)模型分析中發(fā)現(xiàn)，在現(xiàn)有施肥技術(shù)水平條件下，增量施肥的增產(chǎn)效應(yīng)逐漸減少，在施肥水平(420kg/hm2)時(shí)產(chǎn)量達(dá)到最高，其后隨施肥水平的增加反而產(chǎn)生了負(fù)面的減產(chǎn)效應(yīng)。類(lèi)似地，當(dāng)能量投入產(chǎn)出處于邊際報(bào)酬遞增階段反映生態(tài)系統(tǒng)還處于良好的運(yùn)行狀態(tài)，說(shuō)明生態(tài)系統(tǒng)的可持續(xù)性良好，繼續(xù)增加投入可以提高產(chǎn)出； 當(dāng)能量投入產(chǎn)出處于邊際遞減階段時(shí)，說(shuō)明生態(tài)系統(tǒng)的可持續(xù)性已經(jīng)較弱，繼續(xù)增加投入使得產(chǎn)出增長(zhǎng)的速率降低。

1 研究方法與數(shù)據(jù)來(lái)源

1.1 研究方法

經(jīng)濟(jì)學(xué)上的規(guī)模報(bào)酬遞減規(guī)律認(rèn)為，在技術(shù)水平和要素價(jià)格保持不變的條件下，隨著生產(chǎn)規(guī)模的擴(kuò)大，所有投入要素按同一比例不斷增加會(huì)引起產(chǎn)量增加，但到一定程度，產(chǎn)量增加比例會(huì)逐漸遞減。類(lèi)似于規(guī)模報(bào)酬遞減規(guī)律，把農(nóng)業(yè)作為一個(gè)能量投入產(chǎn)出的生產(chǎn)系統(tǒng)，在氣候、技術(shù)等條件不變的情況下，隨著能量不斷投入，每一能量投入帶來(lái)的產(chǎn)出也呈現(xiàn)先遞增后遞減的趨勢(shì)。因此，根據(jù)生態(tài)最大承載量的理論，當(dāng)投入超過(guò)區(qū)域生態(tài)最大承載量所需要的投入時(shí)，能量投入產(chǎn)出將處于邊際報(bào)酬遞減的加速階段，從而出現(xiàn)生態(tài)負(fù)效應(yīng)使得可持續(xù)性惡化(圖1)。

以中國(guó)耕作制度區(qū)劃為區(qū)域劃分的生態(tài)最大承載力下，其投入產(chǎn)出的“S”型曲線(xiàn)共分3個(gè)區(qū)，最左側(cè)區(qū)段表示健康發(fā)展，投入產(chǎn)出報(bào)酬遞增，到P點(diǎn)達(dá)到邊際效益PX最大； 中間區(qū)段表示進(jìn)入預(yù)警期，報(bào)酬遞減，到區(qū)域生態(tài)最大承載量時(shí)可獲得生物產(chǎn)出能最大時(shí)所需要的投入，即投能臨界點(diǎn)E，投入能和產(chǎn)出能平衡； 最右側(cè)區(qū)段表示不可持續(xù)，生態(tài)負(fù)效應(yīng)遞增。

根據(jù)上述規(guī)律，該文通過(guò)傳統(tǒng)的Logistic曲線(xiàn)來(lái)求出產(chǎn)投效益的拐點(diǎn)P劃分邊際報(bào)酬遞增和遞減的區(qū)間，通過(guò)拋物線(xiàn)的一階導(dǎo)數(shù)來(lái)作為投能臨界點(diǎn)E確定生態(tài)負(fù)效應(yīng)遞增區(qū)間(圖2)。

具體計(jì)算方法為：

首先，根據(jù)邊際生產(chǎn)力、平均生產(chǎn)力和生產(chǎn)彈性系數(shù)找到產(chǎn)投效益的拐點(diǎn)P、相應(yīng)的PX，將Logistic曲線(xiàn)劃分成兩個(gè)階段。

(1)

其中，生物量y即種植業(yè)生物產(chǎn)出能； 農(nóng)業(yè)生產(chǎn)方式x即各種農(nóng)業(yè)生產(chǎn)模式下的要素投入能； 最大生態(tài)最大承載量K即在各生產(chǎn)要素均為合理?xiàng)l件下的農(nóng)牧業(yè)能夠達(dá)到的最大產(chǎn)出能；a、b分別為估算系數(shù)。

通過(guò)式(1)求出式(2)～(4)，當(dāng)邊際生產(chǎn)力等于平均生產(chǎn)力，即當(dāng)EEI=1的時(shí)候，拐點(diǎn)產(chǎn)投效益拐點(diǎn)值P出現(xiàn)，因此可以找到對(duì)應(yīng)的PX。

(2)

(3)

(4)

其次，通過(guò)拋物線(xiàn)方程求出投能臨界點(diǎn)E，尋找負(fù)效應(yīng)指數(shù)出現(xiàn)的第三階段。通過(guò)對(duì)拋物線(xiàn)方程進(jìn)行求導(dǎo)，令之等于0求得能投產(chǎn)出最大的點(diǎn)所需要的投入能，即投能臨界點(diǎn)E。

最后，根據(jù)PX和E，將農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)(ASDI)劃分的3個(gè)區(qū)間，將取值范圍設(shè)定在[-1， 1]之間。第一區(qū)間：當(dāng)XE時(shí)，指數(shù)為(E-X)/X。

(5)

如圖2所示：第一區(qū)間ASDI[0， 1]表示健康發(fā)展，投入產(chǎn)出報(bào)酬遞增，到P點(diǎn)達(dá)到邊際效益PX最大； 第二區(qū)間ASDI[1， 0]進(jìn)入預(yù)警期，報(bào)酬遞減，到投能臨界點(diǎn)E達(dá)到區(qū)域生態(tài)最大承載量，投入能和產(chǎn)出能平衡； 第三區(qū)間ASDI[0，-1]表示不可持續(xù)，生態(tài)負(fù)效應(yīng)遞增。

1.2 數(shù)據(jù)來(lái)源

該方法是基于農(nóng)業(yè)統(tǒng)計(jì)年鑒及相關(guān)資料，進(jìn)行數(shù)據(jù)的收集、整理與分析。經(jīng)過(guò)各類(lèi)數(shù)據(jù)的折能系數(shù)進(jìn)行折能轉(zhuǎn)換，計(jì)算出農(nóng)業(yè)生態(tài)系統(tǒng)單位面積上的能量投入和產(chǎn)出。

其中，投入能主要為人工輔助能的投入，包括有機(jī)能和無(wú)機(jī)能。有機(jī)能包括勞動(dòng)力和有機(jī)肥； 無(wú)機(jī)能包括農(nóng)機(jī)械總動(dòng)力、農(nóng)村用電量、化肥、農(nóng)藥、農(nóng)膜的投入； 其他投入為耕地。產(chǎn)出能則包括各類(lèi)種植作物的第一性產(chǎn)出，包括經(jīng)濟(jì)產(chǎn)出和秸稈。各類(lèi)投入要素和產(chǎn)出要素的折算系數(shù)主要參照劉巽浩[36]、吳湘淦[37]、駱世明[39]等的研究，折算標(biāo)準(zhǔn)如表1和表2。

表1 農(nóng)業(yè)生產(chǎn)方式投入要素能量折算標(biāo)準(zhǔn)

能量類(lèi)型103MJ/hm2來(lái)源折算方法無(wú)機(jī)能農(nóng)用機(jī)械210吳湘淦將農(nóng)業(yè)機(jī)具的馬力數(shù)或臺(tái)數(shù)折算成公斤數(shù)，乘以0 1(10年折舊期)燃油45國(guó)家標(biāo)準(zhǔn)標(biāo)準(zhǔn)柴油農(nóng)田用電12 5國(guó)家標(biāo)準(zhǔn)氮肥91吳湘淦磷肥13 3吳湘淦鉀肥9 0吳湘淦農(nóng)藥102吳湘淦塑料薄膜107李波有機(jī)能勞力3500劉巽浩等1個(gè)勞動(dòng)力1年按300d計(jì)，一天需要投入食物1255萬(wàn)J，但扣除糞便部分的能；農(nóng)村勞動(dòng)力-外出務(wù)工人員有機(jī)肥13 5劉巽浩等以畜牧業(yè)和人的糞便中的有機(jī)質(zhì)計(jì)算其實(shí)有量，但需要因地區(qū)考慮丟失量；秸稈還田量的有機(jī)質(zhì)

表2 各類(lèi)作物經(jīng)濟(jì)系數(shù)和能量折算標(biāo)準(zhǔn)

籽 粒秸 稈種類(lèi)103MJ/hm2來(lái)源收獲指數(shù)103MJ/hm2來(lái)源糧食產(chǎn)量 稻谷15 1劉巽浩0 35～0 614 1駱世明 小麥16 3劉巽浩0 3～0 414 6駱世明 玉米16 3劉巽浩0 25～0 414 6駱世明 高粱15 9劉巽浩0 3915 1駱世明 谷子16 1劉巽浩0 3915 1駱世明 大豆20 9劉巽浩0 12～0 2415 1駱世明 薯類(lèi)0 42劉巽浩0 60～0 7814 6駱世明經(jīng)濟(jì)作物 棉花16 7劉巽浩0 11～0 2215 1駱世明 花生23 0劉巽浩0 41～0 5115 1駱世明 蔬菜20 4×0 004186劉巽浩0 7315 1駱世明 向日葵27 86劉巽浩0 21～0 4015 1駱世明 油菜籽26 31劉巽浩0 25114 3駱世明 黑芝麻22 2劉巽浩0 3415 1駱世明 麻類(lèi)16 3劉巽浩0 315 1駱世明 甘蔗2 68劉巽浩0 715 1駱世明 甜菜2 5劉巽浩0 60～0 8515 1駱世明 煙葉16 128劉巽浩0 52～0 6715 1駱世明 茶葉17 3駱世明0 516 7駱世明 蘋(píng)果2 18劉巽浩0 316 7駱世明 柑桔1 80劉巽浩0 316 7駱世明 梨園1 84劉巽浩0 316 7駱世明 葡萄1 80劉巽浩0 316 7駱世明 棗子4 66駱世明0 316 7駱世明 柿子3 45駱世明0 316 7駱世明 核桃31 2駱世明0 316 7駱世明 板栗8 7駱世明0 316 7駱世明 木材16 7劉巽浩 注：該文所用的能量折算方法主要參考劉巽浩[36]，吳湘淦[37]，駱世明[39]和農(nóng)業(yè)技術(shù)經(jīng)濟(jì)手冊(cè)等

2 實(shí)例驗(yàn)證

印江縣位于貴州省東北部之銅仁市中部，屬云貴高原向湘西丘陵和四川盆地過(guò)渡的大斜坡地帶，地處喀斯特化較為嚴(yán)重的亞熱帶高原山區(qū)，全縣總面積19.69萬(wàn)hm2，喀斯特巖溶面積達(dá)10.147萬(wàn)hm2，占全縣國(guó)土面積的51.74%，各種喀斯特形態(tài)并存，是典型的喀斯特地貌，石漠化比較嚴(yán)重。海拔在600～1000m之間，年均氣溫16.8℃，日照時(shí)間長(zhǎng)達(dá)1296h，無(wú)霜期近300d，年降雨量1100mm左右，冬無(wú)嚴(yán)寒，夏無(wú)酷暑?！鞍松揭凰环痔铩?，是貴州省典型的山區(qū)農(nóng)業(yè)縣，地球北回歸線(xiàn)上唯一的綠洲。全縣總?cè)丝?2萬(wàn)，土家族苗族占總?cè)丝诘?0%以上。

該文選擇收集了1980～2014年印江縣的農(nóng)業(yè)統(tǒng)計(jì)年鑒及相關(guān)資料，進(jìn)行審核、整理與分析。經(jīng)過(guò)各類(lèi)數(shù)據(jù)的折能系數(shù)進(jìn)行折能轉(zhuǎn)換，計(jì)算出該縣農(nóng)業(yè)生態(tài)系統(tǒng)單位耕地面積上的能量投入和產(chǎn)出。

2.1 模型擬合結(jié)果

將整理出的單位面積上農(nóng)業(yè)能量投入和產(chǎn)出數(shù)據(jù)根據(jù)式(1)輸入SPSS Statistics 20統(tǒng)計(jì)軟件，進(jìn)行S曲線(xiàn)方程估計(jì)，估計(jì)結(jié)果如(6)所示。根據(jù)估計(jì)結(jié)果，得出最大產(chǎn)出能K值為276.87MJ/hm2，參數(shù)a、b分別為2.316、-0.023。模型可決系數(shù)R2為0.587，各參數(shù)都通過(guò)了顯著性檢驗(yàn)，模型擬合程度較好。

(6)

把擬合得到的a、b、K值帶入式(2)、式(3)、式(4)等，求出產(chǎn)投效益拐點(diǎn)值P，即得出對(duì)應(yīng)的PX，即PX=17.196萬(wàn)MJ/hm2。如圖3所示：

通過(guò)S曲線(xiàn)和拋物線(xiàn)的參數(shù)估計(jì)情況，得出拋物線(xiàn)參數(shù)a、b、c分別為： 2.571、-0.005、-67.04，對(duì)拋物線(xiàn)求導(dǎo)并令之等于0求得拋物線(xiàn)頂點(diǎn)對(duì)應(yīng)橫坐標(biāo)，即投能臨界值E，E=21.642萬(wàn)MJ/hm2。

最后將計(jì)算得出的 PX=17.196萬(wàn)MJ/hm2和E=21.642萬(wàn)MJ/hm2，根據(jù)式(5)計(jì)算出印江縣的農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)。

圖3 印江縣產(chǎn)投效益拐點(diǎn)值

圖4 拋物線(xiàn)的參數(shù)估計(jì)情況

圖5 印江縣農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)

2.2 農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)分析

通過(guò)圖5可以看出，印江縣的農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)具有波動(dòng)性； 1980～2001年呈現(xiàn)緩慢遞增的趨勢(shì)， 2002年后波動(dòng)下降至2007年出現(xiàn)了負(fù)效應(yīng)指數(shù)，而后迅速回升至較高可持續(xù)的發(fā)展水平。其中2007年下降的主要原因?yàn)榛实氖┯昧坑?006年的6757t增加到2.3383萬(wàn)t。氮肥由上一年的3484t增加到7761t，磷肥由2006年的939t增加到5425t。

當(dāng)?shù)卦?007年以前以畜牧業(yè)為主，羊和豬的存欄量較高，因此得出有機(jī)能投入處于較高水平。其中， 2007年后由于產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，播種面積的增加，種植業(yè)主要向蔬菜和水果產(chǎn)業(yè)發(fā)展，因此導(dǎo)致化肥施用量增加，從而導(dǎo)致無(wú)機(jī)肥的大幅度增加，使得總投能超過(guò)了臨界點(diǎn)E=21.642萬(wàn)MJ/hm2，處于負(fù)效應(yīng)指數(shù)出現(xiàn)的投能階段。同時(shí)， 2007年后由于外出務(wù)工人員增加，從事農(nóng)業(yè)活動(dòng)的大都為老年人，再加上產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，畜牧業(yè)逐漸減少，種植業(yè)逐漸增加，因此2007年后有機(jī)能大幅度的降低，使得農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)回升至健康發(fā)展的階段。

3 結(jié)論與討論

農(nóng)業(yè)可持續(xù)發(fā)展涉及方面眾多，生態(tài)可持續(xù)是其中的核心內(nèi)容。當(dāng)前，我國(guó)農(nóng)業(yè)可持續(xù)發(fā)展面臨轉(zhuǎn)變生產(chǎn)方式的主要矛盾，該文采用能量投入產(chǎn)出的方法，進(jìn)行農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)的研究，以下是該文幾個(gè)方面的結(jié)論和討論。

(1)基于能量投入產(chǎn)出法，構(gòu)建農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)來(lái)確定農(nóng)業(yè)生產(chǎn)狀態(tài)的方法是可行的，通過(guò)S曲線(xiàn)和拋物線(xiàn)兩個(gè)方程來(lái)確定農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)的拐點(diǎn)和臨界點(diǎn)，能夠得出農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)的規(guī)律。說(shuō)明評(píng)價(jià)農(nóng)業(yè)可持續(xù)發(fā)展簡(jiǎn)便可行的方法是存在的，下一步考慮是否可以嘗試?yán)脝畏匠虂?lái)進(jìn)行農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)的研究。

(2)ASDI作為農(nóng)業(yè)可持續(xù)發(fā)展的研究，或?qū)檗r(nóng)業(yè)生態(tài)安全預(yù)警提供了一個(gè)可供參考的指標(biāo)。生態(tài)安全是可持續(xù)發(fā)展的前提和基礎(chǔ)，當(dāng)一個(gè)國(guó)家或地區(qū)所處的自然生態(tài)環(huán)境狀況能維系其經(jīng)濟(jì)社會(huì)的可持續(xù)發(fā)展時(shí)，其生態(tài)環(huán)境是安全的，反之則不安全。當(dāng)前通用的農(nóng)業(yè)生態(tài)安全預(yù)警模型及軟件尚不多，缺乏定量或半定量衡量特定區(qū)域農(nóng)業(yè)生態(tài)質(zhì)量的專(zhuān)業(yè)工具?；谀芰客度氘a(chǎn)出構(gòu)建的農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)，或可為農(nóng)業(yè)生態(tài)安全預(yù)警提供一個(gè)可供參考的指標(biāo)。

(3)然而該文也存在一定的不足之處，農(nóng)業(yè)可持續(xù)發(fā)展的研究涉及方面較多，該文僅從生態(tài)學(xué)的視角進(jìn)行分析，采用能量方法從農(nóng)業(yè)生產(chǎn)方式的角度來(lái)研究農(nóng)業(yè)可持續(xù)發(fā)展指數(shù)，對(duì)全面分析農(nóng)業(yè)可持續(xù)發(fā)展情況存在不足。下一步希望建立多步驟的評(píng)價(jià)體系，將經(jīng)濟(jì)、社會(huì)等其他因素考慮進(jìn)去以完善評(píng)價(jià)方法，以達(dá)到對(duì)農(nóng)業(yè)可持續(xù)發(fā)展的正確診斷和預(yù)測(cè)。

[1] 夏光.建立系統(tǒng)完整的生態(tài)文明制度體系——關(guān)于中國(guó)共產(chǎn)黨十八屆三中全會(huì)加強(qiáng)生態(tài)文明建設(shè)的思考.環(huán)境與可持續(xù)發(fā)展， 2014， 39(2)： 9～11

[2] 農(nóng)業(yè)部. 全國(guó)農(nóng)業(yè)可持續(xù)發(fā)展規(guī)劃(2015-2030).2015

[3] 國(guó)務(wù)院. 中國(guó)落實(shí)2030年可持續(xù)發(fā)展議程創(chuàng)新示范區(qū)建設(shè)方案.2016

[4] 孫新章. 中國(guó)參與2030年可持續(xù)發(fā)展議程的戰(zhàn)略思考.中國(guó)人口·資源與環(huán)境.2016, 26(1): 1～7

[5] 周海林. 農(nóng)業(yè)可持續(xù)發(fā)展?fàn)顟B(tài)評(píng)價(jià)指標(biāo)(體系)框架及其分析.農(nóng)村生態(tài)環(huán)境, 1999，(3)： 6～10

[6] Smith L G，Williams A G，Pearce B D.The energy efficiency of organic agriculture：a review.Renewable Agriculture &Food Systems， 2014, 30(3): 280～301

[7] Rosa A D L，Siracusa G，Cavallaro R.Emergy evaluation of Sicilian red orange production.A comparison between organic and conventional farming.Journal of Cleaner Production， 2008, 16(17): 1907～1914

[8] Laurence G.Smith The energy efficiency of organic agriculture：A review.Renewable Agriculture and Food Systems， 2014， 30(3)； 280～301

[9] 劉巽浩， 徐文修，李增嘉，等.農(nóng)田生態(tài)系統(tǒng)碳足跡法：誤區(qū)、改進(jìn)與應(yīng)用——兼析中國(guó)集約農(nóng)作碳效率.中國(guó)農(nóng)業(yè)資源與區(qū)劃， 2013, 34(6): 1～11

[10]黃益強(qiáng)， 張穩(wěn)成.基于AHP法的廣東省農(nóng)業(yè)可持續(xù)發(fā)展評(píng)價(jià).現(xiàn)代農(nóng)業(yè)科技, 2012，(19)： 308～310

[11]戶(hù)艷領(lǐng)， 陳志國(guó)，劉振國(guó)，等.基于熵值法的河北省農(nóng)業(yè)用水利用效率研究.中國(guó)農(nóng)業(yè)資源與區(qū)劃， 2015, 36(3): 136～142

[12]吳炎， 杜棟.基于BP神經(jīng)網(wǎng)絡(luò)的江蘇省農(nóng)業(yè)可持續(xù)發(fā)展能力評(píng)價(jià).中國(guó)系統(tǒng)工程學(xué)會(huì)決策科學(xué)專(zhuān)業(yè)委員會(huì)學(xué)術(shù)年會(huì).2009:113～119

[13]Khoshnevisan B，Rafiee S，Omid M，et al.Application of artificial neural networks for prediction of output energy and GHG emissions in potato production in Iran.Agricultural Systems， 2014, 123(3): 120～127

[14]虎陳霞， 傅伯杰，陳利頂.可持續(xù)農(nóng)業(yè)評(píng)價(jià)研究進(jìn)展.中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2005，(2)： 11～14

[15]劉麗英. 基于PCA和DEA方法的北京市可持續(xù)發(fā)展能力的評(píng)價(jià)研究.數(shù)理統(tǒng)計(jì)與管理.2013, 32(2): 202～210

[16]王韓偉. 投入產(chǎn)出分析在農(nóng)業(yè)經(jīng)濟(jì)系統(tǒng)中的應(yīng)用.農(nóng)機(jī)化研究, 2005，(3)： 36～38

[17]劉巽浩. 能量投入產(chǎn)出研究在農(nóng)業(yè)上的應(yīng)用.農(nóng)業(yè)現(xiàn)代化研究, 1984，(4)： 15～20

[18]王宏廣， 劉巽浩.黃淮海平原幾個(gè)典型縣農(nóng)田生態(tài)系統(tǒng)能量分析.北京農(nóng)業(yè)大學(xué)學(xué)報(bào), 1987，(4)： 475～483

[19]Klzllaslan N.Energy use and input-output energy analysis for apple production in Turkey.Journal of Food Agriculture & Environment， 2009, 7(2): 419～423

[20]Abdi R，Taki M，Akbarpour M.An analysis of energy input-output and emissions of greenhouse gases from agricultural productions.International Journal of Natural & Engineering Sciences， 2012： 73～79

[21]Karkacier O，Goktolga Z G.Input-output analysis of energy use in agriculture.Energy Conversion & Management， 2005， 46(9-10)： 1513～1521

[22]Ozkan B，Akcaoz H，F(xiàn)ert C.Energy input-output analysis in Turkish agriculture.Renewable Energy， 2004, 29(1): 39～51

[23]Tipi，T，?etin，B，Vardar，A.An analysis of energy use and input costs for wheat production in Turkey.Journal of Food Agriculture & Environment， 2009, 7(2): 352～356

[24]Blancard S，Martin E.Energy efficiency measurement in agriculture with imprecise energy content information.Energy Policy， 2014, 66(66): 198～208

[25]Panesar B S，F(xiàn)luck R C.Energy productivity of a production system：Analysis and measurement.Agricultural Systems， 1993, 43(4): 415～437

[26]Gurunathan S，Palanisami K.Energy Use and its Efficiency in Tamil Nadu Agriculture：A Case Study of Different Groundwater Development Regions in Coimbatore District of Tamil Nadu.Agricultural Economics Research Review， 2008, 21(1): 48～53

[27]Heap B，Kasemo B.Panel discussion on energy efficiency.Ambio， 2010， 39： 22

[28]Leach G.Slesser M.Energy Equivalents of Network Inputs to Food Production Process.Glasgow：University of Strathclyde.1973： 56～78

[29]Ma Z Y，Edwards-Jones G.Optimizing the external energy input into farmland ecosystems：a case study from Ningxia，China.Agricultural Systems， 1997， 53(2-3)： 269～283

[30]Wit C T D.Resource use efficiency in agriculture.Agricultural Systems， 1992， 40(1)： 125～151

[31]Bonny S.Is agriculture using more and more energy?A French case study.Agricultural Systems， 1993, 43(1): 51～66

[32]葉旭君， 李全勝.農(nóng)田生態(tài)系統(tǒng)能量投入產(chǎn)出模型及其應(yīng)用研究.生物數(shù)學(xué)學(xué)報(bào), 2002，(1)： 104～109

[33]魏小紅， 藺海明，胡恒覺(jué).甘肅河西地區(qū)農(nóng)田生態(tài)系統(tǒng)能量投入產(chǎn)出及效率分析.草業(yè)學(xué)報(bào), 2001，(3)： 79～84

[34]尹鈞， 曹衛(wèi)星，周乃健，等.農(nóng)田能量投入產(chǎn)出規(guī)律的研究.干旱地區(qū)農(nóng)業(yè)研究, 1999，(3)： 97～102

[35]尹鋼吉， 孫雪梅，孫艷玲，等.控制灌溉條件下水稻施肥水平與產(chǎn)量關(guān)系研究.現(xiàn)代農(nóng)業(yè)科技, 2012，(17)： 9～10

[36]劉巽浩. 能量投入產(chǎn)出研究在農(nóng)業(yè)上的應(yīng)用.農(nóng)業(yè)現(xiàn)代化研究, 1984，(4)： 15～20

[37]吳湘淦. 我國(guó)能源結(jié)構(gòu)分析與農(nóng)業(yè)現(xiàn)代化(續(xù)).農(nóng)業(yè)工程, 1980，(4)： 9～15

[38]朱希剛. 農(nóng)業(yè)技術(shù)經(jīng)濟(jì)分析方法及應(yīng)用.北京：中國(guó)農(nóng)業(yè)出版社.1997

[39]駱世明. 農(nóng)業(yè)生態(tài)學(xué).北京：中國(guó)農(nóng)業(yè)出版社， 2008