王昊+仇蕾

摘 要：傳統(tǒng)農(nóng)業(yè)生產(chǎn)過程中對(duì)能源的不合理使用以及大量廢棄物的排放加劇了農(nóng)村生態(tài)環(huán)境的惡化。通過建立系統(tǒng)動(dòng)力學(xué)模型來分析研究綠色農(nóng)業(yè)的未來發(fā)展趨勢(shì)，以此來分析解決相關(guān)環(huán)境問題。該模型的主要作用是在江蘇省農(nóng)村相關(guān)數(shù)據(jù)的基礎(chǔ)上，分析不同政策對(duì)綠色農(nóng)業(yè)的影響。通過對(duì)比三種不同的政策組合的模擬結(jié)果，找出最適宜當(dāng)?shù)貙?shí)際情況的政策，并給出相應(yīng)的建議。

關(guān)鍵詞：綠色農(nóng)業(yè)；系統(tǒng)動(dòng)力學(xué)；政策模擬

中圖分類號(hào)： F127 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1673-291X（2017）12-0012-07

引言

盡管作為主要產(chǎn)業(yè)，農(nóng)業(yè)在中國一直扮演著重要的角色，但高速發(fā)展的農(nóng)業(yè)經(jīng)濟(jì)與日益惡化的生態(tài)環(huán)境之間的矛盾正威脅著我國長遠(yuǎn)發(fā)展。因此，采取措施緩和這種矛盾勢(shì)在必行。

近年來，綠色農(nóng)業(yè)這個(gè)概念不斷進(jìn)入公眾視線，愈來愈多的人已經(jīng)意識(shí)到綠色農(nóng)業(yè)在改善生態(tài)環(huán)境和公眾健康方面的重要意義。在我國，綠色農(nóng)業(yè)這個(gè)概念是20世紀(jì)90年代提出，當(dāng)時(shí)是指利用先進(jìn)的技術(shù)、設(shè)備和經(jīng)驗(yàn)來改善農(nóng)業(yè)生產(chǎn)與生態(tài)環(huán)境之間的關(guān)系[1]。鑒于此概念的廣泛性與復(fù)雜性，本文將建立相關(guān)指標(biāo)分析未來數(shù)十年中綠色農(nóng)業(yè)發(fā)展水平。

顯然，用于分析的相關(guān)指標(biāo)受到不同復(fù)雜且動(dòng)態(tài)的因素影響，并且該指標(biāo)將在一定程度上隨相關(guān)政策因子變化而改變。因此，使用系統(tǒng)動(dòng)力學(xué)建立相關(guān)模型分析生態(tài)環(huán)境與綠色農(nóng)業(yè)之間的關(guān)系是十分適合的。

一、研究方法

（一）系統(tǒng)動(dòng)力學(xué)的起源與概念

系統(tǒng)動(dòng)力學(xué)理論是Forrester教授于1958年在麻省理工學(xué)院提出。起初，該理論是用于研究諸如生產(chǎn)與雇員狀態(tài)之間的波動(dòng)以及股票市場(chǎng)的不穩(wěn)定性這類問題的工具[2～3]。1968年，F(xiàn)orrester教授在其著作《系統(tǒng)的原則》一書中闡釋了動(dòng)力學(xué)行為的基本原則以及系統(tǒng)結(jié)構(gòu)的相關(guān)概念[4]。該書中所討論的相關(guān)原則如今被廣泛且普遍地應(yīng)用于系統(tǒng)預(yù)測(cè)與政策制定方面。

系統(tǒng)動(dòng)力學(xué)是系統(tǒng)科學(xué)和管理科學(xué)的一條重要分支。它將計(jì)算機(jī)仿真技術(shù)與系統(tǒng)科學(xué)理論相結(jié)合，主要研究系統(tǒng)行為與反饋的結(jié)構(gòu)[5]。通過建立系統(tǒng)動(dòng)力學(xué)模型，人們不僅可以繪制出反映復(fù)雜變量之間關(guān)系的圖像，還可以預(yù)測(cè)關(guān)鍵變量在不同假設(shè)條件下的未來趨勢(shì)[6～8]。此外，系統(tǒng)動(dòng)力學(xué)對(duì)數(shù)據(jù)的準(zhǔn)確性要求偏低，使得人們可以利用有限的資源在不同變量之間建立復(fù)雜、動(dòng)態(tài)且非線性的關(guān)系[9～10]。

（二）系統(tǒng)動(dòng)力學(xué)應(yīng)用現(xiàn)狀

在系統(tǒng)動(dòng)力學(xué)理論被提出后的數(shù)年間，其被廣泛應(yīng)用于多個(gè)領(lǐng)域。Forrester教授的學(xué)生Meadows教授提出了一個(gè)重要結(jié)論：如果像人口、污染、工業(yè)化、農(nóng)作物產(chǎn)量以及資源消耗這些因素按照當(dāng)時(shí)的趨勢(shì)不受控制地增長下去，那么增長的極限就將在一百年內(nèi)發(fā)生[11]。這個(gè)觀點(diǎn)的提出背景就是羅馬俱樂部的建立，該非政府組織主要召集科研人員使用系統(tǒng)動(dòng)力學(xué)解決人口過剩和資源枯竭等問題。

如今，系統(tǒng)動(dòng)力學(xué)已經(jīng)成為一個(gè)方便的數(shù)理仿真工具，應(yīng)用于環(huán)境工程、商業(yè)管理、政策評(píng)估等領(lǐng)域中。Geum和Lee通過將技術(shù)路線圖與系統(tǒng)動(dòng)力學(xué)方法結(jié)合，找到了一種分析比較三種汽車產(chǎn)業(yè)分紅的方法[12]。Zawadzki將系統(tǒng)動(dòng)力學(xué)應(yīng)用到IT工程管理，并通過系統(tǒng)動(dòng)力學(xué)仿真找到了有關(guān)IT管理的一個(gè)最佳方案[13]。通過將系統(tǒng)動(dòng)力學(xué)與蒙特卡羅方法相結(jié)合，Jeon與Shin提出了一個(gè)新的長期技術(shù)評(píng)估方法用于評(píng)估可再生能源技術(shù)[14]。此外，其他領(lǐng)域如人口管理和可持續(xù)發(fā)展的研究者也成功地將系統(tǒng)動(dòng)力學(xué)模型與生態(tài)系統(tǒng)評(píng)估相結(jié)合[15～16]。

（三）研究范圍

本文主要研究對(duì)象是江蘇省農(nóng)村的生態(tài)環(huán)境與綠色農(nóng)業(yè)。江蘇省作為經(jīng)濟(jì)強(qiáng)省，其綠色農(nóng)業(yè)的發(fā)展水平在全國處于領(lǐng)先地位，其農(nóng)業(yè)產(chǎn)值僅次于山東省與河南省。研究的時(shí)間跨度為2010—2050年。

二、建模與測(cè)試

（一）SD模型設(shè)計(jì)

GAE-SD的完整形式（如圖1所示）。原始數(shù)據(jù)來源于《中國環(huán)境年鑒》《中國農(nóng)村統(tǒng)計(jì)年鑒》《中國農(nóng)業(yè)年鑒》和中國農(nóng)業(yè)統(tǒng)計(jì)資料。

（二）模型測(cè)試

軟件Vensim本身提供了非常方便的檢驗(yàn)?zāi)Ｐ偷墓ぞ?。點(diǎn)擊菜單項(xiàng)“檢驗(yàn)?zāi)Ｐ汀焙汀皢挝粰z驗(yàn)”，屏幕上顯示“模型無問題”和“單位無問題”的信息，這表明該模型通過的軟件檢驗(yàn)。

總之，GAE-SD模型是完全可以用于分析當(dāng)前政策背景下綠色農(nóng)業(yè)生態(tài)系統(tǒng)的影響程度以及預(yù)測(cè)不同政策背景下的未來趨勢(shì)。

三、模擬與討論

上述測(cè)試已證明GAE-SD模型可以用于模擬當(dāng)前和未來的綠色農(nóng)業(yè)指數(shù)以及其他變量的趨勢(shì)。首先，基于當(dāng)前政策的趨勢(shì)將被模擬并用于分析當(dāng)前綠色農(nóng)業(yè)水平。其次，幾個(gè)政策方案將會(huì)被分別模擬以找出改善綠色農(nóng)業(yè)現(xiàn)狀的理想模式。

（一）基于當(dāng)前政策的模擬

圖2b曲線1在2012年劇烈上升后就保持輕微的波動(dòng)，直到2044年后才開始增長。由于這條曲線始終處于原點(diǎn)之上，圖2a所代表的綠色農(nóng)業(yè)指數(shù)才能維持穩(wěn)定的增長。這說明了在當(dāng)前政策的影響下，綠色農(nóng)業(yè)的發(fā)展水平將隨之時(shí)間而上升。

圖2c中唯一的曲線表示了綠色農(nóng)業(yè)對(duì)當(dāng)?shù)剞r(nóng)業(yè)產(chǎn)出的貢獻(xiàn)程度。這條曲線之所以幾乎沒有上升，是因?yàn)榫G色農(nóng)業(yè)指數(shù)的增速?zèng)]有當(dāng)?shù)剞r(nóng)業(yè)產(chǎn)出的增速快。這意味著在當(dāng)?shù)剞r(nóng)業(yè)經(jīng)濟(jì)中，綠色農(nóng)業(yè)的地位并不足以高到引起政府投入的重視程度。

簡而言之，當(dāng)前綠色農(nóng)業(yè)的相關(guān)政策會(huì)取得一定程度的效果，但是當(dāng)?shù)卣矐?yīng)該對(duì)綠色農(nóng)業(yè)中不成熟的部分有足夠的重視。

（二）基于不同政策方案的模擬

三個(gè)模塊的政策變量將會(huì)被分別調(diào)整，其他兩個(gè)模塊的變量保持不變。在模擬開始前，為了方便，8個(gè)政策變量將被重新命名為政策變量1—8。原始的政策變量和新的政策變量人（如下表所示）。

方案1將代表清潔能源模塊的政策變量3—5從0.02、0.02、0.02調(diào)整至0.05、0.05、0.05。相關(guān)變量的模擬趨勢(shì)（如圖3a所示）。從圖中可以看出，代表清潔能源模塊的三個(gè)變量在方案1作用下明顯比原始方案增長的要快。特別是曲線2遠(yuǎn)遠(yuǎn)高于曲線1，其代表的清潔能源發(fā)電量幾乎是原始方案下的3倍，即使是之前表現(xiàn)不佳的秸稈固化量也有了明顯的增加。

方案2將代表污染管制模塊的政策變量6—8從0.03、0.03、0.03調(diào)整至0.06、0.06、0.06。圖3b顯示，在方案2作用下，2010-2046年間，化肥使用量發(fā)生了明顯的減少，表現(xiàn)在圖像上就是曲線1遠(yuǎn)遠(yuǎn)低于曲線2，并且在2046年后接近底部。這種明顯的減少，意味著綠色農(nóng)業(yè)將更加依賴有機(jī)肥而非對(duì)生態(tài)環(huán)境會(huì)造成危害的化肥。由于農(nóng)藥使用量和COD本身就處于低水平，所以曲線3、4和曲線5、6這兩組曲線間并沒有多少明顯的差別。

宏觀政策變量1和變量2從0.01、0.01調(diào)整至0.03、0.03。同預(yù)期一樣，圖3c的曲線1、3明顯高于曲線2、4，這表明在耕地保護(hù)和教育上的小小投資就可以收獲客觀的效益。

新政策方案使得綠色農(nóng)業(yè)指數(shù)占每單位農(nóng)業(yè)產(chǎn)出的數(shù)值要好于原始方案，具體狀況（見圖3c）。值得注意的是，代表方案2的曲線2在2030年后開始下降，原因是方案2使得綠色農(nóng)業(yè)指數(shù)的增速不如原始方案下的增速，以至于當(dāng)?shù)剞r(nóng)業(yè)產(chǎn)出的增速高于綠色農(nóng)業(yè)指數(shù)的增速。2034年后，代表方案1的曲線1有了一個(gè)明顯的上升趨勢(shì)，這也說明了方案1對(duì)于長期而言，是個(gè)最佳方案。

參考文獻(xiàn)：

[1] 王德勝.綠色農(nóng)業(yè)的發(fā)展現(xiàn)狀與未來展望[J].中國農(nóng)業(yè)資源與區(qū)劃，2016，（2）：226-230.

[2] Forrester J.W.Industrial Dynamics：A Major Breakthrough for Decision Makers[J].Harvard Business Review，1958，（4）.

[3] Forrester J.W.Industrial Dynamics[J].Journal of the Operational Research Society，1997，（10）：1037-1041.

[4] Forrester J.W.Principles of systems[J].1990.

[5] 鐘永光，賈曉菁，錢穎.系統(tǒng)動(dòng)力學(xué).第2版[M].北京：科學(xué)出版社，2013.

[6] 王其藩.系統(tǒng)動(dòng)力學(xué)：修訂版[M].上海：上海財(cái)經(jīng)大學(xué)出版社，2009.

[7] Chunyang H.E.，Shi P.，Chen J.，et al.Developing land use scenario dynamics model by the integration of system dynamics model and cellular automata model[J].中國科學(xué)：地球科學(xué)，2005，（11）：1979-1989.

[8] Teng J.，Wang P.，Wu X，et al.Decision-making tools for evaluation the impact on the eco-footprint and eco-environmental quality of green building development policy[J].Sustainable Cities & Society，2016，（23）：50-58.

[9] 陳雪萍.基于系統(tǒng)動(dòng)力學(xué)的雅安市土地綜合承載力的可持續(xù)研究[D].成都：四川師范大學(xué)，2013.

[10] 徐曉燕.太原循環(huán)經(jīng)濟(jì)的系統(tǒng)動(dòng)力學(xué)分析[D].太原：山西財(cái)經(jīng)大學(xué)，2014.

[11] Meadows D.H.，et al.The limits to growth：a report for the club of Romes project on the predicament of mankind[J].New York：U-niverse Books，1972.

[12] Geum Y.，Lee S.，Park Y.Combining technology roadmap and system dynamics simulation to support scenario-planning：A case of car-sharing service[J].Computers & Industrial Engineering，2014，（1）：37-49.

[13] Zawadzki T.Applying system dynamics modeling to IT project management[C]// Methods and Models in Automation and Robotics.2009：152-157.

[14] Jeon C.，Shin J.Long-term renewable energy technology valuation using system dynamics and Monte Carlo simulation：Photovoltaic technology case[J].Energy，2014，（4）：447-457.

[15] Vafa-Arani H.，Jahani S.，Dashti H.，et al.A system dynamics modeling for urban air pollution：A case study of Tehran，Iran[J].Trans-portation Research Part D Transport & Environment，2014，（31）：21-36.

[16] Nabavi E.，Daniell K.A.，Najafi H.Boundary matters：the potential of system dynamics to support sustainability？[J].Journal of Cleaner Production，2017，（140）：312-323.

[責(zé)任編輯 劉嬌嬌]