基于能值分析的种植业生产系统可持续发展水平评价--以山东省郓城县为例.pdf

江西农业学报 2023，35（06）：200206ActaAgriculturaeJiangxiDOI:10.19386/ki.jxnyxb.2023.06.028基于能值分析的种植业生产系统可持续发展水平评价以山东省郓城县为例丁 豪1,2，李宗新3，张文君4，钱 欣3，樊庆琦2，李升东2*（1.青岛农业大学 农学院，山东 青岛 266109；2.山东省农业科学院 作物研究所，山东 济南 250100；3.山东省农业科学院 玉米研究所，山东 济南 250100；4.山东省农业科学院 国资处，山东 济南 250100）摘 要:为探究郓城县种植业生产系统的可持续发展水平，采用能值分析法，分析了郓城县19782019年种植业生产系统可持续发展指数的变化趋势，并具体分析了可持续发展指数对能值投入的敏感性以及单位面积上不同作物的可持续发展指数。结果表明：19782019年，郓城县种植业生产系统可持续发展指数呈现先上升后下降趋势，其中，2002年可持续发展指数最低，为1.21；2019年可持续发展指数为1.29，接近指数1的警戒线。在投入敏感性方面，不可更新工业辅助能对可持续发展指数影响最为显著，不可更新工业辅助能每增加15%，可持续发展指数平均下降约7%。从单位面积上不同作物的可持续发展指数来看，大豆单位面积的可持续发展指数最高。基于以上结果，建议通过减少化肥的不合理使用、改善耕作措施、升级农机装备、加强种植基础设施建设、合理调整种植结构来提高郓城县种植业的可持续发展水平。关键词:能值分析；可持续发展；种植业；敏感性；郓城县 中图分类号：F326.1 文献标志码：A 文章编号：1001-8581（2023）06-0200-07Evaluation of Sustainable Development Level of Planting Production System Based on Emergy Analysis:A Case Study of Yuncheng County,Shandong Province DINGHao1,2,LIZong-xin3,ZHANGWen-jun4,QIANXin3,FANQing-qi2,LISheng-dong2*(1.CollegeofAgronomy,QingdaoAgriculturalUniversity,Qingdao266109,China;2.CropResearchInstitute,ShandongAcademyofAgriculturalSciences,Jinan250100,China;3.MaizeResearchInstitute,ShandongAcademyofAgriculturalSciences,Jinan250100,China;4.ManagementOfficeofStateOwnedAssets,ShandongAcademyofAgriculturalSciences,Jinan250100,China)Abstract:InordertoexplorethesustainabledevelopmentleveloftheplantingproductionsysteminYunchengCounty,theemergyanalysismethodwasusedtoanalyzethechangetrendofthesustainabledevelopmentindexoftheplantingproductionsysteminYunchengCountyfrom1978to2019,andthesensitivityofthesustainableindextotheenergyinputandthesustainableindexofdifferentcropsinunitareawerespecificallyanalyzed.Theresultsshowedthatfrom1978to2019,thesustainabledevelopmentindexoftheplantingproductionsysteminYunchengCountyshowedatrendoffirstrisingandthendeclining.Thesustainableindexin2002wasthelowest,withavalueof1.21.Thesustainableindexin2019was1.29,whichwasclosetothewarninglineofindex1.Intermsofinputsensitivity,nonrenewableindustrialauxiliaryenergyhasthemostsignificantimpactonthesustainabledevelopmentindex.Thenonrenewableindustrialauxiliaryenergyincreasesby15%,andthesustainableindexdecreasesbyabout7%onaverage.Fromtheperspectiveofsustainabilityindexofdifferentcropsperunitarea,thesustainabilityindexofsoybeanperunitareaisthehighest.Basedontheaboveresults,itissuggestedtoimprovethesustainabledevelopmentlevelofplantingindustryinYunchengCountybyreducingtheunreasonableuseofchemicalfertilizer,improvingfarmingmeasures,upgradingagriculturalmachineryandequipment,strengtheningtheconstructionofplantinginfrastructure,andreasonablyadjustingtheplantingstructure.Key words:Emergyanalysis;Sustainabledevelopment;Plantingindustry;Sensitivity;YunchengCounty收稿日期：2023-04-25基金项目：国家现代农业产业技术体系(CARS-03-22)；山东省农业重大技术协同推广计划(SDNYXTTG-2022-10)；山东省重点研发计划(重大科技创新工程)(2020CXGC0108052)。作者简介：丁豪(1996)，男，山东莱阳人，硕士研究生，研究方向为小麦栽培生理。*通信作者：李升东。6 期201丁豪等：基于能值分析的种植业生产系统可持续发展水平评价以山东省郓城县为例农业是国民经济的基础，一个国家的农业发展水平在很大程度上决定着该国的经济发展水平。其中，种植业以其较大的规模和较高的可操作性而成为整个农业的重要组成部分之一。但随着科学技术的发展和社会的进步，人口、资源与环境的矛盾日益凸显，严重威胁着种植业的可持续发展1。因此，如何有效评价一个国家或地区的种植业可持续发展水平，对于提高其农业综合生产能力，保护资源环境，加快农业现代化等均具有重要意义2。由于种植业生产是自然因素与人为因素双重作用的结果3，因而无法直接对其可持续性进行比较和评价。因此，本文采用能值分析法，将种植业生产过程中的各个投入与产出因素统一为太阳能值，而后进行比较和分析。随着可持续发展概念的逐渐普及，能值分析法被用于分析与评价各个领域及系统的可持续能力，如方敏哲等4分析了内蒙古磴口县20032015年土地生态经济系统的可持续性情况；张永杰等5对2014年甘肃省定西市农田生态系统的环境资源条件、投入产出结构和可持续性进行了分析；钱欣等6研究了滨州市2019年种粮大户冬小麦夏玉米粮食生产系统的效率和可持续性；李洋等7分析了19992015年山东省种养系统可持续发展指数的空间格局及其演变趋势；刘新茂等8分析了广东省1993年种植业生态系统的可持续发展状况；Zhuang等9分析了在本地化战略条件下，中国三大主粮作物可持续能力的变化情况。综合来看，多数研究聚焦于大尺度，对于整个县域尺度的种植业可持续能力变化及其主要影响因素的研究较少。而县域单位又是连接城市与农村的桥梁，是中央农业政策的具体执行者10，县域地区种植业发展状况的好坏在一定程度上关系着国家能否更好地调整产业结构和推进农业现代化，因此，对于县域尺度的研究具有重要意义。基于此，本文选择县域级作为本研究的空间尺度，通过延长时间尺度来分析县域级种植业的可持续能力变化情况，并尝试分析其影响因素。本文选取山东省菏泽市郓城县作为研究对象，以郓城县19782019年种植业的相关数据为基础，采用能值分析法，计算其历年来种植业投入与产出，并得出相关的可持续发展指数，同时分析其可持续发展指数的变化情况，以期为提高该地区的可持续发展能力以及探索其种植业高质量发展道路提供参考。1 材料与方法1.1 研究区域概况郓城县位于山东省西南部，介于115 40 11608 E、35 19 35 52 N之间，总面积为1643km2。气候上属暖温带半湿润季风气候，降水集中，雨热同季，年均气温12.314.4，年均降水量550.6633.9mm。全县属于黄河冲积平原，耕地面积113232.1hm2，农作物播种面积为199656.8hm2，常年种植小麦、玉米、蔬菜等农作物，粮食产量连续多年突破10亿kg，是山东省的产粮大县。1.2 研究方法能值分析法由美国著名生态学家Odum于20世纪80年代创立11，其理论依据在于：地球上的一切能量大都来自太阳12，根据特定的太阳能值转换率，可将种植业生产系统中的不同物质以及能量统一转换为太阳能值，单位为太阳能焦耳(Solaremjoules,sej)13-24。本文将郓城县种植业生产系统的各项投入与产出作为研究对象，运用相应的能值转化率将种植业系统中的各项指标统一转化为太阳能值，其中包括可更新的环境资源R(包括太阳能、风能、雨水势能、雨水化学能等)、不可更新的环境资源N(包括表土损失能、灌溉水)、不可更新的工业辅助能FN(包括化肥、农药、农膜、柴油、机械、电力等)、可更新的有机能FR(主要包括劳力、种子和还田的秸秆)、系统产出能值Y(包括各类作物的籽粒、果实和未还田秸秆)，由此确定的系统边界如图1所示13。为方便计算，暂不考虑资金对系统能值的影响作用。1.3 数据来源本研究所采用的数据主要来自 郓城统计年鉴(19782019)菏泽统计年鉴(19942020)山东统计年鉴(19782020)，由于郓城县部分农作物秸秆直接还田，缺少记录数据，故在本研究中，参考李书田等25的研究成果，确定不同作物的草谷比；参考毕于运等26的研究成果，通过利用不同作物的草谷比，计算出不同作物的秸秆产量。还田部分的秸秆属于有机肥范畴，故将其能值归为可更新有机能，未还田部分的秸秆归类为系统产出能值。通过调查发现，郓城县主要秸秆还田农作物为小麦、玉米、稻谷、豆类、花生，平均秸秆还田率约为70%。此外，本研究采用的地球生物圈能值基准为最新值1.21025sej/a27，能量折算系数与能值江 西 农 业 学 报35 卷202转换率参照李洋等7,14,15,18,22的研究成果。1.4 郓城县种植业系统能值的评价指标为更好地研究和分析郓城县种植业可持续发展水平及相关驱动因素，在参考相关文献9,14-15,23,28的基础上，建立了3项郓城县种植业系统能值评价指标，分别为净能值产出率(EYR)、环境负载率(ELR)、能值可持续发展指数(ESI)，具体计算公式为：EYR=Y/(FN+FR)（1）ELR=(N+FN)/(R+FR)（2）ESI=EYR/ELR （3）式（1）式（3）中，EYR为净能值产出率，Y为能值总产出量，FN为不可更新工业辅助能值，FR为可更新有机能值，ELR为环境负载率，N为不可更新自然资源能值，R为可更新自然资源能值，ESI为能值可持续发展指数。2 结果与分析2.1 能值投入分析从能值投入情况来看(图2)：19782019年郓城县能值总投入整体呈波动式增长趋势，1978年郓城县能值总投入为2.921020sej，增长到2019年的1.191021sej，总体增加了307.5%，年均增加7.5%。能值总投入由可更新环境资源、不可更新环境资源、不可更新工业辅助能和可更新有机能构成，1978年这4种能值的比例分别为22.9%、4.1%、34.9%、38.1%，至2019年这4种能值的比例则分别为4.9%、1.5%、58.5%、35.1%。1978年郓城县不可更新工业辅助能与可更新有机能投入分别为1.021020、1.111020sej，2019年不可更新工业辅助能增长到6.941020sej，增长率为580.4%，可更新有机能增长到4.161020sej，增长率为274.8%。图1 郓城县种植业生产系统边界图2 郓城县能值投入的变化趋势2.2 能值产出分析从图3中可以看出，19782019年郓城县能值产出呈波动式增长趋势。2019年郓城县种植业生产系统总能值产出比1978年增加了1.7610216 期203丁豪等：基于能值分析的种植业生产系统可持续发展水平评价以山东省郓城县为例sej，增长了435.5%。能值产出类型可大致分为粮食作物、油料作物、棉麻作物、蔬菜作物、水果作物产出五大类。从几种作物产出能值的变化趋势来看，粮食作物呈现先上升后下降再上升的趋势，1978年粮食作物产出能值为3.201020sej，1999年上升到1.241021sej，而后逐年下降，2003年降至7.311020sej，2004年又开始逐渐上升，截至2019年郓城县粮食作物能值产出为1.871021sej。油料作物呈现先上升后下降的趋势，1978年油料作物能值产出为1.171018sej，而后逐年上升，2001年达到最大值1.621020sej，而后逐年下降，到2019年，郓城县油料作物产出能值为3.771019sej。棉麻类作物则呈现波动式增长的趋势，其中19831985、20022008年期间产出较高，19781981、19982000年期间产出较低，2019年郓城县棉麻作物产出能值为3.821019sej。蔬菜、水果类作物产出能值从1994年之后开始逐渐增加，1994年郓城县蔬菜作物与水果作物产出能值分别为1.611019和5.981019sej，2019年 分 别 增 加 到1.331020、7.931019sej，增长率分别为726.1%和32.6%。图3 郓城县能值产出的变化趋势2.3 系统能值评价2.3.1 环境负载率 环境负载率主要指可更新资源与不可更新资源之间的比值，其数值的大小代表着该地区种植业生产系统所承受环境压力的大小。由图4可看出，19782019年郓城县种植业生产系统的环境负载率整体呈现上升的趋势，由1978年的0.64上涨到2019年的1.50，年均上涨3.28%；20042015年，环境负载率有所下降，2016年之后，环境负载率则逐年增加，说明随着生产资料的投入，郓城县种植业生产系统的资源环境压力正在逐渐增大。图4 19782019年郓城县环境负载率的变化趋势2.3.2 净能值产出率 净能值产出率主要指作物产出能值与不可更新工业辅助能、可更新有机能之和的比值，可用来衡量种植业生产系统的投入产出效率。由图5可看出，19782019年期间，郓城县种植业生产系统的净能值产出率整体上呈现下降的趋势。其中，1983年的净能值产出率较高，为5.45；2018年数值最低，为1.90。整体来看，20092019年，郓城县净能值产出率较前几年有所下降，说明在种植业生产过程中，存在生产资料不合理使用的现象。图5 19782019年郓城县净能值产出率的变化趋势2.3.3 可持续发展指数 可持续发展指数为净能值产出率与环境负载率的比值，可用来衡量该地区的种植业生产系统的可持续发展能力。一般来说，可持续发展指数在110之间，说明此地区种江 西 农 业 学 报35 卷204植业生产系统具有发展潜力；可持续发展指数在10以上，说明此地区种植业生产系统的自然资源利用不足；可持续发展指数小于1，说明此地区种植业生产系统的资源环境压力较大24。由图6可看出，19782019年，郓城县种植业生产系统的可持续发展指数呈现先上升后下降的趋势。其中，19781988年10a间，郓城县种植业生产系统的可持续发展指数保持在较高水平，1984年达到极值9.81，其余年份的可持续发展指数基本保持在4.0以上。说明此阶段，郓城县作物产出有所增加，但资源利用不足。19892002年期间，郓城县种植业生产系统的可持续发展指数呈现下降的趋势，由1989年的6.77下降至2002年的最低值1.21，年均下降1.4036%。说明此阶段，随着种植业生产资料投入的增加，作物产量有所提升，但由于化石能源的过量使用，导致环境压力逐渐增加，使得2002年郓城县种植业生产系统的可持续发展指数降至最低点。20032019年，受种植业结构调整、农业生产措施转变等因素的影响，郓城县种植业生产系统的可持续发展指数有所回升，截至2019年，可持续发展指数为1.29。整体来看，郓城县种植业生产系统的可持续发展指数均在1.0以上，理论上具有发展潜力，但也应加以重视，在可持续发展指数接近1.0的警戒线时，应减少对资源的不合理使用，提高其环境承载能力，以保证郓城县种植业发展的可持续性。图6 19782019年郓城县可持续发展指数的变化趋势2.3.4 敏感性分析 为了更好地探究可持续发展指数对各能值投入的敏感性，在其他条件不变的情况下，分别将可更新环境资源、不可更新环境资源、不可更新工业辅助能、可更新有机能的其中之一增加或减少15%，观察能值投入指标变化对可持续发展指数的影响程度。由图7可知，在各能值投入指标中，可持续发展指数对不可更新工业辅助能的敏感程度最高，不可更新工业辅助能每增加15%，可持续发展指数平均下降约7%；其次是对可更新环境资源，而对不可更新环境资源的敏感程度最低，且可更新环境资源与可更新有机能的增加对可持续发展指数的影响是正向的，不可更新工业辅助能与不可更新环境资源的增加对可持续发展指数的影响是负向的。由于不可更新工业辅助能主要包含化肥、电力、机械、地膜、农药、柴油等能值投入，其中化肥为氮、磷、钾、复合肥投入能值总和。为具体探究各项投入对可持续发展指数的影响程度，每次只针对一项指标，分别上下调整5%、10%和15%，观察其可持续发展指数的变化幅度，具体结果如图8所示。由此可以看出，化肥投入对可持续发展指数的影响程度最高，其次为机械和电力投入。2.3.5 不同作物的可持续发展指数的比较 由于粮食作物大于经济作物的播种面积，从而使得粮食作物的产出能值远高于经济作物。为更好地比较不同作物的可持续性，本文选用1hm2面积上不同作物的投入与产出数据，通过进行能值换算，进而求出各个作物的可持续发展指数。由图9可以看出，不同作物的可持续发展指数由大到小依次为大豆棉花花生玉米小麦水果稻谷蔬菜，环境负载率由大到小依次为水果蔬菜稻谷棉花小麦花生大豆玉米，净能值产出率由大到小依次为棉花水果大豆小麦花生稻谷玉米蔬菜。3 结论本研究以19782019年为时间尺度，利用能值分析法，通过计算分析郓城县种植业生产系统的能值投入与产出，得到如下结论：(1)19782019年郓城县种植业生产系统总投入能值增加了307.5%，总产出能值增加了435.5%；(2)19782019年郓城县种植业生产系统的环境负载率整体呈上升的趋势，净能值产出率整体呈下降的趋势，可持续发展指数呈先上升后下降的趋势，2002年之后的可持续发展指数波动稳定，整体在1以上，具有可持续发展潜力，但也应注意保护环境，提高资源利用率，防止可持续发展指数降至1以下；(3)在能值投入中，可持续发展指数对不可更新工业辅助能的敏感程度最高，其中化肥、机械、电力为主要的影响因素；(4)从单位面积可持续发展指数来看，大豆的可持续发展指数最高。6 期205丁豪等：基于能值分析的种植业生产系统可持续发展水平评价以山东省郓城县为例图7 可持续发展指数对各能值投入的敏感性分析图8 不可更新工业辅助能敏感性分析数据根据 全国农产品成本收益汇编 整理和计算所得，均为近5a的平均值。图9 不同作物的环境负载率、净能值产出率 及可持续发展指数比较4 建议4.1 减少化肥的不合理使用本研究发现，化肥是对可持续发展指数影响最为显著的1个因素，过量地施用化肥不仅不能提高作物产量，反而会对环境造成污染。结合相关试验以及实地调研结果，可针对部分施肥较高的作物进行合理减肥。其中，小麦、玉米的氮肥投入较高，可在小麦、玉米周年轮作地块减少氮肥施用量。以小麦、玉米周年轮作纯氮施用量400kg/hm2为对照，在此基础上小麦、玉米分别减少氮肥15%，可显著提高氮肥利用率。同时，针对蔬菜、水果等作物，可适当减少复合肥施用量，并增施有机肥，以此来提高郓城县的可持续发展指数，并保障其产品品质。4.2 实行保护性耕作措施推广免耕、少耕、秸秆还田等保护性耕作措施，其优势在于，可在一定程度上促进农业生产中水资源以及土壤资源的可持续利用。如免耕播种技术可将农作物秸秆进行粉碎并覆盖在土地表面，以达到减少土地水分流失和减少耕作流程的目的29。4.3 更新农机设备，减少农机作业流程本研究发现，农机使用也是影响可持续发展指数的主要因素之一，这主要与农机作业过程中排放大量二氧化碳以及土壤粉尘有关。传统农业生产过程繁琐，包括收获、秸秆粉碎、旋耕、播种、施肥、镇压等一系列农事操作流程，繁琐的流程需要反复使用农机，由此降低了种植业的可持续发展水平。因此，可与农机公司合作，将整地、播种等流江 西 农 业 学 报35 卷206程融合到一台农机上，一次性完成整地、播种、镇压等流程，从而大大提高农机的使用效率，进而提高农业的可持续发展水平。4.4 加强农田基础设施建设本研究中，电力投入对郓城县种植业可持续发展指数的影响仅次于化肥和机械投入。调查发现，电力的投入主要是用于排水和灌溉。建议加强农田基础设施建设，修建排水渠和灌溉水井，做到农田“旱能灌，涝能排”，推广滴灌、渗灌等节水灌溉技术，在节省农业用电量的同时，还可以节约水资源的投入。4.5 合理调整种植业结构从单位面积上不同作物的可持续发展指数来看，大豆、花生、玉米等作物的可持续发展指数较高，因此可在部分地区推广“玉米+大豆”间作、“玉米+花生”间作等种植模式，在增加大豆供给量的同时，提高郓城县种植业生产的可持续能力。参考文献：1 姚瑶,李广,王钧,等.甘肃省种植业系统的能值动态J.草业科学,2021,38(4):785-795.2 姜力玮,李绍辉,陈强强.甘肃省农业生态系统可持续性评价 J.新疆农垦经济,2020(8):23-30.3 蒋桂芹,吴泽宁,胡彩虹.水资源对农业生产贡献率的能值计算方法及应用 C/第七届中国水论坛论文集.2009:395-400.4 方敏哲,岳德鹏,张启斌,等.基于能值分析的磴口县土地生态经济系统可持续性研究 J.西北林学院学报,2017,32(4):178-185,228.5 张永杰,王平.定西市农田生态系统能值空间分异特征研究 J.甘肃农业大学学报,2017,52(5):97-102,111.6 钱欣,王小龙,张慧,等.华北平原种粮大户小麦玉米轮作系统能值分析:以山东省滨州市为例 J.山东农业科学,2020,52(10):42-47.7 李洋,孙志刚,张旭博,等.种养系统可持续发展指数的空间格局及其演变趋势:以山东省为例 J.应用生态学报,2019,30(7):2371-2383.8 刘新茂,蓝盛芳,陈飞鹏.广东省种植业系统能值分析J.华南农业大学学报,1999,20(4):111-115.9 ZhuangMH,LiuYZ,YangY,etal.ThesustainabilityofstaplecropsinChinacanbesubstantiallyimprovedthro-ughlocalizedstrategies J.RenewableandSustainableEnergyReviews,2022,154:111893.10 戎承法.县域农业资源配置效率研究:以阜南县种植业为例 D.北京:中国农业大学,2004.11 OdumHT.Environmentalaccounting:Emergyandenvi-ronmentaldecisionmaking M.NewYork:JohnWiley&Sons,1996:182-303.12 王旭熙,彭立,苏春江,等.四川省农业生态系统能值分析与可持续性评估:基于五大经济区差异分析 J.地域研究与开发,2015,34(5):128-132,166.13 蓝盛芳,钦佩,陆宏芳.生态经济系统能值分析 M.北京:化学工业出版社,2002:3-95.14 杨增玲,于璐,王德睿,等.呼和浩特市种植业生产系统能值分析与评价 J.中国农业大学学报,2019,24(8):179-193.15 徐宏.江苏省农业生态可持续发展评价 J.中国农业资源与区划,2019,40(8):164-170.16 张思静.张掖市耕地可持续集约利用评价研究 D.兰州:西北师范大学,2021.17 任钇蒙.河南省种植业生态系统能值分析 D.长沙:中南林业科技大学,2016.18 邢开成,龚宇,王璞.华北平原集约农区种植业生态经济系统的能值分析:以河北沧州为例 J.生态环境,2007,16(2):592-597.19 孙雯.基于能值的江苏省耕作系统可持续性分析 D.扬州:扬州大学,2021.20 王红彦.基于生命周期评价的秸秆沼气集中供气工程能值分析 D.北京:中国农业科学院,2016.21 郭勇.郓城县农业有机废弃物资源调查和环保昆虫转化处理模式探究 D.泰安:山东农业大学,2020.22 朱玉林.基于能值的湖南农业生态系统可持续发展研究 D.长沙:中南林业科技大学,2010.23 钟颖,沙之敏,顾麦云,等.基于能值分析的稻蛙生态种养模式效益评价 J.中国生态农业学报:中英文,2021,29(3):572-580.24 李俊莉,曹明明.生态脆弱区资源型城市农业生态系统的能值分析:以榆林市为例 J.中国农业科学,2012,45(12):2552-2560.25 李书田,金继运.中国不同区域农田养分输入、输出与平衡 J.中国农业科学,2011,44(20):4207-4299.26 毕于运,高春雨,王亚静,等.中国秸秆资源数量估算J.农业工程学报,2009,25(12):211-217.27 BrownMT,UlgiatiS.Emergyassessmentofglobalrenewablesources J.EcologicalModelling,2016,339:148-156.28 张黎黎,李子君.黄河三角洲耕地利用集约度变化及其驱动因素 J.中国农业资源与区划,2022,43(1):124-134.29 于海波.保护性耕作技术发展现状及展望 J.农机使用与维修,2022(10):61-63.（责任编辑：曾小军）