基于种养平衡的吉林省辽河流域农田畜禽粪便负荷研究

黄鑫，赵兴敏*，苏伟，郑玉昕，王鸿斌，王玉军

（1.吉林农业大学资源与环境学院，长春 130118；2.吉林省环境科学研究院，长春 130012）

改革开放40多年来，我国经济迅猛发展，已经成为世界第二大经济体，但随着我国经济由高速增长阶段转向高质量发展阶段，农业现代化和可持续性发展成为国民经济高质量发展的重要组成部分，绿色发展建设已成为当前我国亟待解决的重要问题。

东北地区是我国重要的粮食生产基地与畜牧业生产基地，根据《2018年中国统计年鉴》数据，东北粮食产量占全国总产量的21.0%，其中玉米产量占全国总产量的33.7%，大豆更是占比43.9%，这在保障国家粮食安全与农产品供需平衡中占有重要的战略地位。辽河平原作为东北三大平原之一，在当地农业生产中发挥着重要作用，而吉林省辽河流域作为辽河的上游地区，其生态环境形势却不容乐观。吉林省第二次全国污染源普查公报显示，2017年吉林省各大类污染源中，畜禽养殖业带来的污染问题比较突出，其排放的化学需氧量（COD）、总氮（TN）、总磷（TP）分别占总排放量的68.9%、47.3%、69.8%，远超工业污染源所带来的危害。随着我国工业现代化不断发展，逐步完善了工业减排管理体制，工业点源污染排放得到了有效控制。按照当前国内外农业面源污染排放的情况，畜禽粪便若处置不当将会给当地生态环境带来巨大风险。目前针对吉林省辽河流域的研究主要集中在水体环境和生态系统等方面，但对农田畜禽粪便污染的研究鲜有报道。

为此，本研究以吉林省辽河流域（四平市区、双辽市、梨树县、公主岭市、伊通县、辽源市区及东辽县）为研究区域，通过各市县2017年畜禽养殖基础数据，结合种植业年度统计数据，进行畜禽粪便负荷及环境风险分析研究，以为防治畜禽养殖污染，科学规划养殖业发展提供科学依据。

1 材料与方法

1.1 数据来源

1.1.1 畜禽粪便负荷核算数据

以吉林省辽河流域畜禽养殖业和种植业为研究对象，养殖业数据来源于各市县第二次全国污染源普查，具体情况如图1所示。

图1 2017年吉林省辽河流域畜禽养殖数量Figure 1 Amount of livestock husbandry in Liaohe River basin of Jilin Province in 2017

通过四平市和辽源市《统计年鉴》，获取2017年耕地面积。畜禽粪便的排泄系数是负荷核算的基础，考虑到不同地区饲养周期和养殖方式差异，根据前人研究的东北地区畜禽有关数据，并结合调研资料确定了吉林省辽河流域畜禽粪便产生量核算相关参数，见表1。

表1 畜禽粪便产生量核算相关参数Table 1 Parameters related to the accounting of livestock manure production

1.1.2畜禽养殖环境承载力数据

根据《畜禽粪污土地承载力测算技术指南》及吉林省农田耕地土壤养分资料，确定了不同作物生长的氮磷养分需求量（表2）。施肥供给氮养分占比为35%，磷占比为45%；粪肥占施肥比例为50%；粪肥当季氮养分利用率为27.5%，磷养分利用率为32.5%。各市县农作物生产情况来源于当地统计年鉴数据（表3）。

表2 吉林省辽河流域主要农作物单位经济产量养分需求量（kg·t-1）Table 2 Demand nutrients per ton economic yield in Liaohe River basin of Jilin Province（kg·t-1）

表3 2017年吉林省辽河流域主要农作物生产情况（t）Table 3 Production of main crops in Liaohe River basin of Jilin Province in 2017（t）

1.2 核算方法

1.2.1 畜禽粪便产生量及耕地负荷核算

畜禽粪便产生量核算采用国家环保总局推荐的测算方法，计算公式为：

式中：为各类畜禽粪便产生量，t·a；为出（存）栏量，头（羽）；为饲养周期，d；为各类粪便日排泄系数，kg·头（羽）·a。

先将各类畜禽粪便统一折算成猪粪当量，再进行畜禽粪便负荷量核算，计算公式为：

式中：为畜禽粪便负荷量，t·hm·a；为各类畜禽粪便猪粪当量总量，t·a；为有效耕地面积，hm；为各类畜禽粪便总量，t·a；为猪粪当量折算系数。

为了更全面反映各地区畜禽粪便环境污染威胁程度，引用畜禽粪便污染警戒值概念，具体警戒值分级情况见表4，其计算公式为：

表4 畜禽粪便污染警戒值分级Table 4 Classification of alarming values for loading amount of livestock feces

式中：为畜禽粪便负荷量警戒值；为有机肥最大适宜施肥量，t·hm·a。由于各地区的地貌条件、土壤类型及气候等原因值的选取不尽相同，基于吉林省辽河流域种植业实际情况考虑，本研究选用30 t·hm·a。

1.2.2 畜禽养殖环境承载力计算

畜禽粪便的总量大、运输成本高，决定了其无法进行大范围、跨地区移动，因此畜禽粪便必须在一定区域内自我消化。本研究基于所有的畜禽粪便就地还田，通过种养平衡核算畜禽养殖环境承载力。首先计算流域内农作物生长对氮、磷养分的需求量，计算公式为：

式中：为农作物生长对氮、磷养分的需求量，kg；G为第类农产品的年总产量，t；H为第类农产品单位产量所需的氮、磷养分量，kg·t。

农作物生长所需的养分主要由土壤、化肥和畜禽粪便还田等提供，粪肥养分按照理论推荐值为50%，所以农作物粪肥养分需求量计算公式为：

式中：为域农作物粪肥养分需求量，kg；为域农作物养分需求量，kg；为施肥供给养分占比，%；为粪肥占施肥比例，%；为粪肥当季利用率，%。

综合考虑畜禽粪便养分在收集、处理和贮存过程中的损失，单位猪当量的氮养分供给量为7 kg，磷为1.2 kg，计算氮（磷）畜禽养殖环境承载力公式为：

式中：为氮（磷）畜禽养殖环境承载力，头猪当量·hm；为域农作物粪肥氮（磷）养分需求量，kg；为单位猪当量粪肥氮（磷）养分供给量，kg；为耕地面积，hm。

为评估当前实际养殖产生量相较于理论环境承载力的饱和程度，研究各地畜禽粪便产生量是否会对农田带来环境污染风险，本研究引入畜禽粪便污染风险指数，并将其进行污染威胁程度分级（表5）。

表5 畜禽养殖环境污染风险分级Table 5 Classification of environmental pollution risk for livestock husbandry

式中：为畜禽粪便污染风险指数；为单位耕地面积实际承载畜禽数量，头猪当量·hm；为氮（磷）畜禽养殖环境承载力，头猪当量·hm。

2 结果与分析

2.1 吉林省辽河流域畜禽粪便产生量

吉林省辽河流域2017年畜禽粪便构成比例和产生量如图2所示，粪便总量为537.59万t。总体来看，畜禽粪便产生量大小依次为生猪>肉牛>蛋鸡>奶牛>肉鸡，分别为229.52万、199.96万、60.26万、36.29万、11.56万t，占比依次为42.69%、37.20%、11.21%、6.75%、2.15%。畜禽粪便主要来源地也不尽相同，其中奶牛和肉牛粪便主要集中在公主岭市、东辽县、双辽市，这3个市县的奶牛和肉牛粪便占牛粪产生总量的68.63%；梨树县、伊通县、公主岭市的猪粪产生量均超过了当地畜禽粪便量的50%，同时这3个市县的猪粪产生量也占据了猪粪总量的65.42%；蛋鸡和肉鸡粪便主要分布在东辽县、梨树县、公主岭市，这3个市县的蛋鸡和肉鸡粪便占鸡粪产生总量的72.91%。

图2 2017年吉林省辽河流域畜禽粪便产生量与耕地面积Figure 2 The production of livestock manure and area of farmland in Liaohe River basin of Jilin Province in 2017

对市县畜禽粪便构成比例进行分析，双辽市、四平市区和东辽县的牛粪占当地畜禽粪便总量的比例分别为71.54%、61.95%和54.94%，梨树县、伊通县、公主岭市和辽源市区猪粪占比分别达到65.61%、52.84%、50.73%和41.62%。由此可以看出，2017年吉林省辽河流域最主要的粪便污染是由牛产生的，其次是猪，这两者的产生量占总量的86.65%。

从各市县的畜禽粪便的猪粪当量看，公主岭市的产生量最高，达159.91万t，占流域畜禽粪便总量的29.39%；其次是梨树县和东辽县，产生量分别为114.48万、109.81万t，各占比为21.04%、20.18%；再次是伊通县和双辽市，这两个市县的畜禽粪便产生量均达到了45万t；而产生量较小的是四平市区和辽源市区，分别为30.03万t和12.13万t。从各市县耕地面积进行分析，公主岭市的耕地面积最大为31.57万hm，其次是梨树县，耕地面积为24.18万hm，这两个市县的耕地面积占总耕地面积的55.69%。四平市区和辽源市区耕地面积少，仅有2.09万hm和1.00万hm，仅占比4.09%，其余市县耕地面积相差较小，但是均达到10万hm以上。由此可以看出，畜禽粪便产生量大的市县往往耕地面积也更多，这为消纳畜禽粪便提供了足够的空间，从而减轻畜禽粪便对当地农田的污染风险。

2.2 吉林省辽河流域畜禽粪便负荷量

目前，吉林省畜禽粪便的主要处理方式是直接还田，将各市县的耕地作为畜禽粪便实际负载区域，核算出各地区的畜禽粪便负荷量，如图3所示。结果表明，该流域的平均畜禽粪便负荷量为5.43 t·hm·a，其中辽源市区的负荷量最高，为12.09 t·hm·a，其次为东辽县、四平市区，负荷量分别为10.50、9.72 t·hm·a，另外伊通县也超过平均畜禽粪便负荷量，为5.81 t·hm·a。而梨树县、公主岭市和双辽市位于流域中下游，农业种植面积相对较多，其负荷量均在平均值以下，故这些区域消纳畜禽粪便的潜能大，对周边环境影响小。

图3 吉林省辽河流域畜禽粪便负荷量Figure 3 The manure load of livestock in Liaohe River basin of Jilin Province

2.3 吉林省辽河流域畜禽粪便负荷预警分析

根据畜禽粪便负荷风险分级，吉林省辽河流域平均警戒值为0.18，对环境暂无明显影响。但是按照县区划分进行核算发现，辽源市区的警戒值超过0.4，对环境已经构成一定程度的威胁，东辽县、四平市区的警戒值处于Ⅱ级预警的边缘，如果畜禽粪便处置不当也会对周边农田产生威胁。其余市县的畜禽粪便污染警戒值较小，对当地环境基本无影响。结合GIS地理信息技术分析发现（图4），辽源市区、东辽县和四平市区位于吉林省辽河流域的中上游，而畜禽养殖场大多靠近水体，这将会对流域水质环境产生较为严重的威胁。

图4 吉林省辽河流域畜禽粪便负荷风险等级Figure 4 The risk level of manure load for livestock in Liaohe River basin of Jilin Province

2.4 畜禽养殖环境承载力

基于种养平衡理念进行核算，环境承载力可以反映当地耕地对畜禽粪便的消纳能力，该数值在一定程度上决定了区域养殖发展的空间。结合欧盟国家耕地消纳粪肥限量值（以氮、磷计分别为170、35 kg·hm，相当于24.29、29.17头猪当量·hm），根据畜禽养殖环境承载力的高低（图5），将吉林省辽河流域承载力分为3个层次：梨树县和公主岭市属于高承载力区，这两个市县基于氮的环境承载力分别为152.54 kg·hm和143.36 kg·hm，基于磷的环境承载力分别是23.63 kg·hm和21.99 kg·hm；伊通县和双辽市为一般承载力区，基于氮的环境承载力处于负荷上限标准的65%~75%，磷的环境承载力分别是19.08 kg·hm和20.49 kg·hm；辽源市区、东辽县、四平市区为承载力偏低区，畜禽养殖环境的氮、磷承载力均小于限量标准的60%，其中四平市区最低，仅有78.64 kg·hm和11.81 kg·hm。吉林省辽河流域平均氮、磷环境承载力分别为131.39 kg·hm和20.65 kg·hm，梨树县和公主岭市超过平均值，畜禽养殖发展潜力较大，而四平市区和辽源市区的环境承载力相对较低，具有畜禽养殖环境污染的风险。

图5 吉林省辽河流域畜禽养殖环境承载力Figure 5 The environmental carrying capacity of livestock husbandry in Liaohe River basin of Jilin Province

2.5 畜禽养殖环境污染风险分析

根据畜禽养殖环境风险指数核算，如图6所示，四平市区和辽源市区的氮污染风险指数分别为1.25、1.07，说明这两个地区没有足够的耕地消纳禽畜粪便，这将会对农田产生一定程度的危害。辽源市区的磷污染风险指数为1.08，属于磷环境中等污染风险，四平市区和东辽县的磷污染风险指数分别为0.97、0.63，属于磷畜禽养殖环境低污染风险。综合氮和磷畜禽养殖环境风险评估发现，辽源市区、东辽县和四平市区属于畜禽养殖风险区，且具有污染地区成块分布、氮和磷养分同时失衡的特征，而其余市县属于畜禽养殖环境安全区，对现有的养殖规模具有一定的承载潜力。

图6 吉林省辽河流域畜禽养殖环境污染风险评估Figure 6 The risk assessment for pollution of livestock husbandry in Liaohe River basin of Jilin Province

3 讨论

3.1 关于畜禽粪便负荷量的核算

吉林省辽河流域近5年来养殖规模不断扩大，畜禽数量相较于2011年增加了近20%，但是畜禽养殖种类构成基本不变。目前流域内主要养殖种类是鸡，占养殖总量的88.45%，牛的养殖量最少，仅占总量的0.78%。不同畜禽粪便的排泄系数差异显著，以奶牛和蛋鸡为例，两者相差186倍。畜禽粪便排泄系数的差异，导致牛的养殖量虽然远不及鸡和生猪，但是前者的粪便产生量却大于后者。

公主岭市、梨树县和双辽市区虽然畜禽粪便量很高，猪粪当量分别是159.9万、114.5万、45.7万t，但是这3个市县的农田负荷风险属于安全等级，结合地区耕地面积进行分析发现，公主岭市、梨树县和双辽市区拥有吉林省辽河流域31.5%、24.2%和17.4%的耕地面积，这使得畜禽粪便拥有更多消纳的空间，从而降低了畜禽粪便负荷污染风险；以四平市区、辽源市区为例，这两个区畜禽猪粪当量最低，但是由于这两个区的耕地面积小，仅占总耕地面积的3.08%和1.01%，从而导致了这两个区农田负荷量大，分别达到9.72 t·hm·a及12.09 t·hm·a。

为了进一步明确畜禽养殖环境风险范围，将辽源市区、四平市区按照行政区域划分核算，结果如表6所示。研究表明龙山区和铁西区污染警戒值分别为0.43和0.52，而西安区和铁东区为0.36和0.22，说明龙山区和铁西区为畜禽粪便污染Ⅱ级风险区，而西安区和铁东区则属于安全区，这进一步缩小了畜禽污染风险范围，为划分畜禽养殖安全区提供了依据。

表6 2017年四平及辽源市区畜禽粪便污染警戒值与风险等级Table 6 Warning valuesand risk levels of livestock manure pollution in Siping City and Liaoyuan City in 2017

从整体上看，吉林省辽河流域大多数市县畜禽粪便负荷风险等级处于Ⅰ级，对周围环境暂无威胁，仅龙山区和铁西区对环境稍有威胁。通过对比流域畜禽粪便相关研究发现，畜禽粪便负荷风险区域划分结果与本文一致，但是负荷量核算结果却相差较大，其平均畜禽粪便负荷量达到了61.16 t·hm·a，超出目前公认的30 t·hm·a的畜禽粪便还田限值，远高于同一时期吉林省平均畜禽粪便负荷量15~18 t·hm·a。早期研究通常采用全国统一的畜禽产排污系数，核算数值较大，而不同省份和地区的养殖模式、密度和粪便处理方式等因素有很大的差异，本研究采用东北地区系数进行核算，结果则较小，造成了本研究畜禽粪便负荷量核算结果与其他学者存在一定的差异。

3.2 畜禽养殖环境污染及风险控制

对比吉林省辽河流域畜禽粪便负荷和环境承载力发现，虽然四平市区和东辽县属于畜禽粪便负荷安全区，但是区域内的氮、磷畜禽养殖环境属于有污染风险状态，耕地土壤的氮磷含量有超标风险，这些市县需要加紧对畜禽粪便的资源化利用，控制对环境的污染影响。

从种养平衡角度进行分析，辽源市区、四平市区和东辽县的实际养殖数量超过或接近该地区环境承载能力的上限，可能引起氮磷供求失衡，造成较大的畜禽环境污染风险，需要减少畜禽养殖规模或者扩大高标准农田建设面积，以使地区种养保持平衡。从这3个市县的作物产量看，玉米产量最大，占区域作物总产量的81%，由于玉米单位产量养分需求量相较于大豆、油料作物要小很多，这些区域若要减少畜禽养殖环境风险，需要改变种植作物类型，进一步增加大豆、油料作物的播种面积，扩大种植业对氮、磷养分的需求量，力求作物养分供求平衡，避免种养结合的盲目性和环境污染问题。

4 结论

（1）2017年吉林省辽河流域畜禽粪便产生量达537.59万t，其中公主岭市产生量最多，为171.33万t，其次为梨树县，为109.84万t。

（2）吉林省辽河流域的畜禽粪便负荷量平均为5.43 t·hm·a，畜禽粪便污染警戒平均值为0.18，辽源市区的粪便负荷量最大，为12.09 t·hm·a，达到Ⅱ级环境污染风险，对当地农田环境产生了威胁。其余市县粪便负荷量都处于安全范围，对周围环境基本无影响。

（3）畜禽养殖氮、磷环境承载力平均值分别为131.39 kg·hm和20.65 kg·hm，四平市区的氮环境承载力最弱，为78.64 kg·hm，辽源市区的磷污染风险指数最大为1.08。根据畜禽养殖环境风险划分，辽源市区属于氮和磷中等污染风险区，东辽县属于氮和磷低污染风险区，而四平市区分别是氮中等风险区和磷低风险区。其余市县对畜禽粪便具有一定的消纳潜力，属于畜禽养殖环境安全区域。