改革开放以来，中国畜禽养殖业发展迅速，产业整体素质实现了质的飞越，现已成为世界上最大的肉、蛋生产国^[1]。伴随着畜禽养殖业的迅猛发展，养殖排放也随之增多，至2015年中国畜禽粪尿排放已达到31.58亿t^[2]，如何对这大量的畜禽粪尿进行有效处理已成为社会各界关注的热点。实践表明，将畜禽粪尿作为养分进行还田利用，是最为经济可行的处理方式^[3]。然而，其同时也是引起污染的重要因素，一旦耕地负荷量超过作物养分需求，则将造成地下水亚硝酸盐污染、江河水体富营养化、土壤矿物质损失等诸多危害^[4-6]。因此，从连接种植业与畜牧业的种养平衡视角出发，进行耕地畜禽粪尿负荷污染风险预警，以及分析污染治理模式效果，对于防治相关污染具有重要的现实意义。

已有大量学者基于种养结合视角，对各地区耕地畜禽粪尿负荷情况进行了分析。一些学者基于畜禽粪尿与作物养分需求之间相对大小这一实际情况，对各地区耕地畜禽粪尿负荷污染风险进行了评价，并由此提出进行总量控制的政策建议；养分的挥发、损耗、土壤供氮能力、灾害对作物养分需求的影响等情况相继被纳入到测算体系之中^[7-10]。另有一些学者在现状分析的基础上，基于不同有机无机配施比例这一假设，运用比较静态分析对各地区污染风险进行了预测，并认为在具有承载潜力的地区，可通过有机替代无机的方式扩大畜禽养殖规模，而养分盈余地区则需要移除多余部分^[11-12]。

上述研究对于摸清中国耕地畜禽粪尿负荷情况提供了坚实的基础，但仍存在以下不足：一是研究方法上，现有文献多是采用静态或比较静态分析，忽视了对中间过程的动态考量；二是研究对象上，现有研究主要围绕种植业、畜牧业两方面展开，忽视了社会、经济因素对于畜禽养殖排放污染存在的影响；三是研究内容上，现有研究多是对现状进行分析，或是在有机肥投入占总养分投入不同比例的假设下进行预测，缺乏对未来实际情况的预警，同时也缺乏对相关治理模式效果的评价。鉴于此，本研究基于种养结合视角，定义“畜禽粪尿污染风险指数”，同时将社会、经济、种植业、畜牧业纳入风险指数核算体系，运用系统动力学方法构建耕地畜禽粪尿负荷系统，动态分析系统各部分间存在的相互作用，预警2018—2025年中国耕地畜禽粪尿负荷污染风险，并从污染治理效果及对作物养分供给影响两个方面，对治理模式的效果进行评价，最终为制定畜禽粪尿治理政策措施提供参考。

1 数据来源与方法 1.1 数据来源

本研究所涉及的经济社会发展状况、畜禽养殖量和农作物产量等数据均来源于历年《中国统计年鉴》 《中国农业年鉴》《国民经济统计公报》《中国畜牧兽医年鉴》以及《中国农业产业发展报告》。

1.2 方法

预警模型实质上由评价模型和预测模型构成^[13]。本研究构建畜禽粪尿污染风险指数评价耕地畜禽粪尿负荷污染风险，采用系统动力学模型预测风险指数、仿真污染治理模式效果。

1.2.1 基于种养平衡的耕地畜禽粪尿负荷污染风险评价模型

（1）评价指标设定

(1)

(2)

(3)

式中:ε为畜禽粪尿污染风险指数，依据《畜禽粪污土地承载力测算技术指南》，以氮素含量进行测算；L为耕地畜禽粪尿氮素负荷量，t；C为耕地畜禽粪尿氮素承载力，t；S为畜禽粪尿氮素含量，t；θ为有机肥氮素损失率；μ为草地畜禽粪尿负荷占比，用以排除草地承载的畜禽粪尿量；D为作物氮肥需求量，t；γ为作物在种植过程中因自然灾害造成产量失真的氮肥需求校正系数；f为化肥氮施用量，t；δ为化肥氮素损失率；s为秸秆氮素含量，t；λ为秸秆养分还田率，p为农村生活氮素排放，t；ω为农村生活排放留存率；B为生物固氮量，t；E为经大气沉降、灌溉输入到耕地的氮素量，t。根据相关文献^[14-15]及《畜禽粪污土地承载力测算技术指南》将θ、δ分别设为36.5%、50%；以草地畜禽粪尿负荷量（根据草地最大载畜量测算^[16]）与总排放之比为草地负荷占比，以减产面积（受灾、成灾、绝收面积与各级灾害损失比例乘积之和）与播种面积之比为氮肥需求校正系数，取2008—2017年平均值，即令μ为8.50%、γ为7.66%；秸秆养分还田率随时间变化存在较大提升，依据《到2020年化肥使用量零增长行动方案》将2008—2017年λ设为35%，2018—2025年λ设为60%；依据杨淑静^[17]的研究将ω设为78.1%。S、D、s、B、p、E的计算公式分别如下：

(4)

式中：n_i为i种畜禽猪当量换算系数，X_i为i种畜禽存栏量，头；u为单位猪当量氮素排泄量，kg。依据《畜禽粪污土地承载力测算技术指南》，按存栏量折算：100头猪相当于15头奶牛、30头肉牛、250只羊、2500只家禽；1个猪当量的氮素排泄量为11 kg。

(5)

(6)

(7)

式中：m_i、g_i、h_i、k_i分为i种作物氮肥需求系数（t·100 t^-1）、谷草比、秸秆含氮量（kg·t^-1）、固氮系数（kg·hm^-2）；Y_i、Z_i分别为i种作物产量（t）、播种面积或水田、旱地面积（hm²）。收集国内公开发表的文章数据，确定作物相关系数（表 1），此外旱地非共生固氮系数为15，水田非共生固氮系数为44.8^{[14-15, 18]}。

(8)

式中：r为常住人口城镇化率；i为人均粪尿氮素排放量，kg；H为总人口，人；依据秦钟等^[12]的研究将i设为4.4。

(9)

式中：j、k分别为大气降氮系数及灌溉供氮系数，kg·hm^-2；W、Z分别为耕地面积和有效灌溉面积，hm²。根据李书田等^[14]的研究，将j、k分别设为20.2、12.1。

（2）预警值级设定

沈根祥等^[7]较早进行耕地畜禽粪尿负荷污染风险的警报分级研究，其依据种植过程中存在的有机无机配施比例在1:5到1:1之间这一实际情况，将预警级别设为5级，认为耕地畜禽粪尿负荷占到作物养分需求的0.4以上就对环境轻微污染，达到0.7、1、1.5以上分别造成稍有、较严重、很严重环境污染，这一分级标准被诸多学者采用^[19-20]。冯爱萍等^[21]依据农作物实际消纳情况及欧盟农业政策，将预警级别设为4级，认为每公顷耕地畜禽粪尿负荷量超过100 kg即对环境构成污染，超过170、225 kg分别构成较高污染和高污染。陈天宝等^[10]综合上述研究，将预警级别设为5级，认为耕地粪尿负荷量占到作物养分需求的0.45以上就对环境构成威胁，达到0.6、0.8、1以上分别构成明显、严重、重大污染威胁。通过对比上述3种分级方式可知，虽然各分级点绝对数值有所不同，但之间的相对比值十分接近，不同之处在于陈天宝等^[10]对较低级别的污染作了更为详细划分。结合本研究需要，借鉴陈天宝等^[10]研究，设置预警值级标准（表 2）。

1.2.2 基于系统动力学的畜禽粪尿污染风险指数预测模型

（1）系统框架及因果关系图

系统动力学（System dynamics，SD）以反馈控制理论为基础，通过建立子系统间的因果关系，实现对具有多元性、动态性等特征的复杂系统的刻画，以及实现对系统发展中行为和趋势的研究，是一种研究地区环境承载及污染情况所常用的仿真方法^[22]。本研究运用系统动力学构建耕地畜禽粪尿负荷系统对2018—2025年中国畜禽粪尿污染风险指数进行预警，并仿真污染治理模式效果。

环境承载系统研究一般由社会、经济、生态3个子系统构成，本研究根据需要，按照各子系统功能将耕地畜禽粪尿负荷系统具体划分为社会、经济、种植业和畜禽养殖业4个子系统（图 1）。其中，种植业子系统主要反映耕地畜禽粪尿氮素承载力，畜禽养殖子系统主要反映耕地畜禽粪尿氮素负荷量，二者直接决定畜禽粪尿污染风险指数。社会子系统、经济子系统分别通过畜禽、作物需求量、农村生活氮素排放和畜禽、作物贸易量、工业净化量对种植业子系统、畜禽养殖子系统产生影响，是驱动畜禽粪尿污染风险指数发生变化的根本力量。

（2）存量流图

根据各子系统间的相互作用，以及子系统内部变量间的逻辑关系对因果关系图进行适当拓展，同时结合畜禽粪尿污染风险指数相关公式，运用Vensim软件建立存量流图（图 2）。

各子系统变量、系统方程及参数设置如下：社会子系统以总人口为状态变量，包含新增人口流速变量、常住人口城镇化率、各畜禽产品人均消费量7个表函数；同时包括计算各畜禽产品需求量、农村常住人口和农村生活氮素排放（公式7）的8个辅助方程，各畜禽人均需求量直接写入辅助方程。经济子系统由各畜禽贸易量6个表函数、治理单价、工业净化投入占比2个参数以及计算畜牧业总产值、环保治理投入、工业净化量的3个辅助方程构成，各肉类价格6个常数直接写入辅助方程中；根据吕志奎^[23]的调研数据，以含氮量计算，将工业净化的治理单价设为39 000元·t^-1，由于当前畜禽养殖场较少采用工业净化的方式对畜禽粪尿进行处理^[3]，故假定当前工业净化投入占比为0，各肉类价格以2017年不变价格计算。种植业子系统包括化肥氮施用量1个状态变量，化肥氮增量1个流速变量，水稻、玉米等作物预计产量以及根据公式5~7测算的预计氮肥需求量、秸秆供氮量、生物固氮量、沉降灌溉供氮量共15个表函数，同时还包括各作物氮肥产出弹性、其他养分需求占比、秸秆养分还田率等共15个参数，以及计算各作物产量、氮肥需求量、氮肥盈余情况等变量的24个辅助方程，各作物消纳系数共11个常数写入辅助方程；粮食作物氮肥产出弹性来源于史常亮等^[24]测算结果，其他作物氮肥产出弹性根据各自“氮肥用量-作物产量关系函数”测算得出^[25-30]，其他养分需求占比指烟草茶叶、豆类等其他作物养分需求占总养分需求的比例，取2008—2017年平均值表示，将其设定为5.9%，依据魏静等^[31]的调研数据，农村生活氮素排放还田率设定为70%，依据《全国农业可持续发展规划（2015—2030年）》将2008年畜禽养殖排放还田率设定为32.5%，2008—2017年间每年提升3.5个百分点。畜禽养殖子系统，由其他畜禽养殖排放占比1个参数以及计算各肉类产量、各畜禽存栏量、各畜禽粪尿氮素含量、畜禽粪尿氮素含量（公式4）的19个辅助方程构成，以每5年平均值表示的各单位畜禽产肉量、出栏率10个分段函数及排放系数5个常数写入辅助方程；依据相关研究^[2]，将驴等其他畜禽排放占比设为2.5%；以2008—2017年平均值表示产肉量、出栏率等变量。此外，存量流图中还包括计算畜禽粪尿污染风险指数的3个辅助方程（公式1~公式3）。

（3）历史性检验

从系统中选取生猪粪尿氮素含量、玉米氮肥需求量和氮肥需求量为检验变量，将2008—2017年的数据代入模型进行模型历史性检验，结果见表 3。

由表 3可知，检验变量实际值与仿真值相关误差均小于10%，表明模型运行结果与实际数据之间拟合度较高，模型能较为准确地描述研究系统的基本现状，具有良好预测效果，可进行下一阶段的模拟仿真。

1.3 情景设置 1.3.1 惯性趋势设置

2018—2025年预测数据，如新增人口、人均畜禽消费量、畜禽贸易量、作物预计产量等变量均来源于中国农科院“中国农业产业模型”（China Agricultural Sector Model，CASM）预测结果^[32]；根据作物预计产量测算出预计氮肥需求量，并结合《到2020年化肥使用量零增长行动方案》，假定氮肥施用量以每年0.3%的速度增长。其他变量预测值依据相关发展规划进行假设，如依据《国家人口发展规划（2016—2030年）》假设常住人口城镇化率每年增长1%；依据《全国农业可持续发展规划（2015—2030年）》，设定2018—2025年畜禽粪尿还田率每年提升1.5个百分点，同时假定农村生活氮素排放按照畜禽养殖排放还田率相同增速变化。

1.3.2 治理情景设置

（1）干预点选择

减量化、资源化、达标排放是应对畜禽粪尿污染的3类主要办法^[33]。“减量化模式”主要通过限制地区养殖规模，从源头降低畜禽粪尿排放总量，最终缓解由此造成的环境污染^[34]；主要涉及系统中各畜禽人均消费量、贸易量2个变量，考虑到进行污染治理的同时应当保障畜禽产品充足供给，满足中国社会饮食需求，故只选取贸易量作为“减量化模式”的干预点。“资源化模式”主要通过提升还田率，并等量调减化肥用量，以“有机替代无机”的方式提升耕地畜禽粪尿氮素承载力，最终避免排放污染^[35]；主要涉及系统中化肥供氮量、畜禽粪尿还田率2个干预点。“达标排放模式”主要强调畜禽养殖场运用工业净化等方式，对养殖排放进行深度处理，实现达标排放；主要涉及系统中工业净化投入占比这一干预点。

（2）情景因子参数设置

在干预点选择基础上，设置“贸易因子”“还田因子”和“环保投入因子”，并结合相关规划及现实约束，调整情景因子，改变相应干预点取值，以形成不同治理情景。具体而言，“贸易因子”主要影响禽肉产品自给率的调整，故以2025年我国畜禽产品自给率（总贸易量/总需求量）为参照，分别设置自给率80%、75%、70% 3个情景，即相当于调整贸易因子，使所有畜禽产品贸易量均提升1.64、2.05、2.46倍。“还田因子”主要影响畜禽粪尿排放还田率的调整，设置还田率在相关规划基础上分别提升10%、20%、30% 3个情景，与之对应的2025年畜禽粪尿还田率分别为74.25%、81.00%、87.75%。“环保投入因子”主要影响工业净化投入占总产值比重的调整，依据《全国城市生态保护与建设规划（2015—2020年）》中对环保投入占GDP比重的要求，以及畜禽养殖企业平均利润率，分别设置环保投入占比2.5%、3.5%、4.5% 3个情景。依据干预点选择及情景因子设置，进行治理情景设置（表 4）。

2 结果与分析 2.1 中国耕地畜禽粪尿负荷惯性趋势预警

仿真结果显示（表 5），随着畜禽需求的持续增长，中国畜禽粪尿氮素含量呈现持续上升态势，到2025年，畜禽粪尿氮素含量将达到1.381 38×10⁷ t，较2017年增长11.59%，其中肉鸡、肉羊、奶牛粪尿氮素含量增幅较大，分别达到35.76%、16.87%、14.88%；绝对数值方面，虽然生猪粪尿氮素含量增幅最小，仅为2.44%，但仍然是主要的氮素来源，占到氮素总量的35.98%，较2017年降低3.23%，其次是肉牛、蛋鸡，粪尿氮素含量分别占到氮素总量的20.17%、13.17%，其他畜禽粪尿氮素含量占30.68%。2025年，化肥氮施用量、秸秆氮素含量、农村生活氮素排放将分别达到2.275 71×10⁷、7.748 2×10⁶、2.174 0×10⁶ t，较2017年分别增长2.43%、4.89%、-14.31%。随着社会农作物需求及氮肥供给的增长，作物产量持续上涨，表现为氮肥需求量持续增长，至2025年，氮肥需求量将达到2.774 30×10⁷ t，约增长3.10%，其中粮食作物氮肥需求上升约2.80%，达到1.670 10×10⁷ t，占氮肥总需求的60.20%，饲料需求上涨带动玉米产量大幅增长是粮食作物氮肥需求上涨的主要推动力量，玉米氮肥需求上升7.16%；蔬菜/瓜果氮肥需求上升4.32%，达到4.647 5×10⁶ t，占氮肥总需求的16.75%；油料、糖类氮肥需求上升3.38%，达到2.239 8×10⁶ t，占氮肥总需求的8.07%。2018—2025年畜禽粪尿污染风险指数持续下降，但始终高于1.83，表明预警区间内耕地畜禽粪尿负荷将持续对环境构成明显污染风险。

2.2 治理情景仿真

在对畜禽粪尿污染进行治理的过程中，须同时保障作物生长所需养分的充足供给，所以本研究将畜禽粪尿污染风险指数与氮肥盈余情况同时设为治理情景的输出结果。

从表 6仿真结果可知，通过提升“减量化模式”治理强度，可有效降低污染风险指数，但即便进行高强度治理，使畜禽产品自给率降低至70%，耕地畜禽粪尿负荷仍将对环境构成明显污染，同时将导致畜牧业总产值减少3693亿元，这表明由于我国畜禽粪尿排放总量巨大，远超过耕地承载力，若主要依靠“减量化模式”，需大幅进口国外畜禽产品，将极大降低我国畜禽产品的自给率，同时畜禽养殖是我国农民的重要收入来源，简单地进行源头减量也不利于我国农民脱贫、增收。“资源化模式”能够在维持氮肥盈余不变的情况下，有效降低风险指数，但即便进行高强度治理，将2025年还田率提升到87.75%，耕地畜禽粪尿负荷仍将对环境构成轻度污染，意味着畜禽粪尿氮素排放总量过大，超过了耕地的承载力，单凭提升资源化强度难以彻底消除污染风险。“达标排放模式”对于降低畜禽粪尿污染风险具有明显效果，进行低强度治理就可达到“减量化模式”高强度治理效果，此时污染治理支出仅为754.05亿元。但此模式同时会减少氮肥供给，当达到中等强度，氮肥供给将减少1.25×10⁶ t，此时氮肥供给将低于预计作物氮肥需求量，这表明，此模式难以在满足作物氮肥需求的情况下，消除污染风险。同时意味着在化肥用量受到控制的情况下，畜禽粪尿将成为重要的氮肥供给来源，所以应当主要将其进行资源化利用以保障氮肥的充足供给，即“达标排放模式”只能作为畜禽粪尿污染治理的一种补充模式。

3 讨论

大量学者基于种养平衡视角，从畜禽存栏及作物产量或耕地面积的现状出发，对某一区域或全国范围的耕地畜禽粪尿负荷情况进行了测算，对评价畜禽粪尿环境污染风险具有重要意义。本研究则从中国禽肉产品、主要农产品消费需求及贸易量的预测数据出发，对中国畜禽粪尿环境污染风险进行预警，还通过情景模拟，分析了不同污染治理模式的实际效果，对未来相关政策制定提供参考。

由于一旦畜禽粪尿排放总量超过当地承载力，就将对环境构成威胁^[3]，所以与前人的研究相同，本研究在计算耕地畜禽粪尿负荷量时，只考虑了畜禽粪尿排放量及氮素损失率^{[9, 11]}。而与前人的研究不同，本研究还同时测算了氮肥盈余情况，并通过调整还田率，分析了“资源化模式”的治理效果。

本研究使用统一的畜禽粪尿排放系数测算了全国畜禽粪尿排放量，但不同地区畜禽粪尿排放系数存在一定差异，因此后续研究应分地区对畜禽产品需求等数据进行预测，并结合各地区畜禽粪尿排放系数，计算畜禽粪尿污染风险指数，从而使研究结果更具科学性和指导意义。

4 结论与政策启示 4.1 结论

本文基于种养平衡视角，定义畜禽粪尿污染风险指数，并应用系统动力学方法进行中国耕地畜禽粪尿负荷污染风险预警研究，以之为基础，通过情景分析，模拟不同治理模式的污染治理效果及对作物养分供给情况的影响。研究结果显示：

（1）2018—2025年，惯性趋势下我国耕地畜禽粪尿负荷始终将对环境构成明显污染。

（2）在满足社会畜禽产品需求这一前提下，“减量化模式”主要通过增加畜禽产品进口量来减少国内养殖规模，最终降低耕地畜禽粪尿负荷污染风险，但将导致我国畜禽产品自给率及总产值大幅下降。

（3）“资源化模式”通过有机替代无机的方式，能够在维持氮肥盈余不变的情况下，有效降低污染风险，但由于氮肥供给过量，单凭提升“资源化”强度难以消除耕地畜禽粪尿负荷污染风险。

（4）“达标排放模式”主要采用工业净化等方式，对畜禽粪尿进行深度处理，对于降低污染风险具有明显效果，但进行相应处理的同时将减少氮肥供给，难以在满足作物氮肥需求的情况下，消除耕地畜禽粪尿负荷污染风险。

4.2 政策启示

（1）进一步采取措施进行污染治理。虽然在《到2020年化肥使用量零增长行动方案》作用下，2018—2025年我国“耕地畜禽粪尿负荷污染风险指数”将逐年降低，但耕地畜禽粪尿负荷始终将对环境构成明显污染，有必要进一步采取措施进行污染治理。

（2）我国耕地畜禽粪尿负荷污染治理应以“资源化模式”为主，“达标排放模式”作为补充。促进畜禽粪污资源化利用，能在缓解污染的同时，减少化肥用量，实现有机无机肥配施，提升作物产量、质量，故应作为首选模式。在受到客观条件制约，无法对畜禽粪尿进行资源化利用或超过当地承载力的地方，采用“达标排放模式”进行治理。

（3）出台激励政策进一步促进畜禽粪尿养分还田。“资源化模式”顺利开展的关键在于畜禽粪尿还田率的提升。首先，应当依据地区耕地承载力优化畜禽养殖业布局，为畜禽粪尿还田创造客观条件。其次，应通过经济激励等方式，鼓励养殖户开展种养结合，实现畜禽粪尿还田。第三，畜禽粪尿制成有机肥更便于长距离运输、使用，故应当对有机肥的生产、流通进行补贴，加强有机肥市场建设，提升畜禽粪尿养分还田率。