随着中国经济的高速增长，国民经济各行业都得到了快速发展，但随之而来的环境问题也日益凸显。根据我国环保部2017年发布的《2016中国环境状况公报》，中国地表水和地下水水质较差的监测断面（监测井）分别占32.3%和60.1%，中国的水环境污染形势非常严峻。水环境污染根据污染源发生类型分为点源污染和非点源污染^[1]。随着工业源等点源污染控制技术的发展，点源污染得到有效控制，非点源污染成为全球范围的突出问题^[2]，且造成水源地污染的主要是农业非点源污染^[3]。2010年全国第一次污染源普查公报显示，农业非点源污染是中国地表水污染的主要来源，而畜禽养殖污染的COD、TN和TP分别占农业污染源的96%、38%和56%^[4]。如何在特定生态、社会、经济环境下发展畜养业，同时保证区域可持续发展，成为目前急需解决的重要问题^[5]。

近年来国内学者从粪便污染、耕地氮磷负荷、畜禽污染的时空特征等角度研究中国的畜禽养殖污染，研究结果表明：畜禽粪便氮磷污染风险较高的地区主要分布于中国的东部地区^[6]，如果不进行政策干预，全国畜禽粪便总污染量将大幅增至2020年的2.98亿t^[7]；从畜禽污染发展上看，我国畜禽养殖造成的污染较严峻，污染程度与人均GDP之间存在倒“U”型曲线关系，且已跨过曲线拐点^[8]；中国农田消纳畜禽粪便的潜力较大，通过畜禽粪便还田利用能够有效解决粪便污染问题^[9]。在相关学者的省域畜禽养殖污染研究中，全国畜禽粪便TN、TP产生量分别是380.56万t和63.89万t，按畜种生猪占41.2%和34.8%，按省域河南居全国第二（分别占13.7%和13.3%）^[9-10]。

河南省是农业大省，全省粮食产量约占全国的十分之一，畜产品产量近几年则一直位居首位^[11]。畜养业在河南农村经济和农民收入中占有重要地位，河南省已经形成以黄河滩区为主的黄河滩区绿色奶业示范带，以商丘、周口、南阳、驻马店、平顶山为主的中原肉牛、肉羊产业带，以许昌、漯河、驻马店等为主的京广铁路沿线生猪产业带，以新乡、安阳、鹤壁为主的豫北蛋肉鸡产业带^[12]。与此同时，畜养业废弃物污染已成为环境污染的重要来源，部分地市畜养污染已超过农田的消纳和承受能力^[13]。河南省生猪基本采用干清粪，堆积发酵后施肥；肉鸡多数为地面平养，有垫料，一次清出堆积发酵施肥；蛋鸡为笼养，干清粪；肉牛、奶牛多为散栏式饲养，干清粪，晾晒后用于双孢菇种植等^[12]。如何科学估算畜养污染物产生量、制定引导其资源化利用的管理和治理政策已成为河南亟待解决的问题。

然而，对接污染源普查畜种和普查污染物的畜禽养殖业污染物产生量（简称畜养产污，下文提到的畜养均包含畜和禽）核算研究不多，多数文献针对特定畜种^[14]、单污染物或较少污染物指标^[15]进行排污分析，河南省畜养产污研究更鲜有研究^[16]。农田消纳与承载负荷的研究文献主要结合具体研究区域实际^[10]，寻求潜在的规律性问题，而目前针对河南省域的畜养产污农田消纳和承载负荷问题也鲜有研究。鉴于此，本研究借助环保部《畜禽养殖业源产排污系数手册》的分畜种产污系数，从畜养产污角度分析河南近年来的畜养产污粪便、尿液、COD、TN、TP、Cu、Zn等污染物的规律，及河南省畜养氮磷的农田消纳与承载负荷，旨在探索一种对接污染源普查工作的省域畜养产污分析方法，同时为河南省畜养产污管理政策的制定提供依据。

本研究畜禽养殖产污的核算品种选取猪（生猪、能繁母猪）、牛（肉牛、奶牛）、鸡（肉鸡、蛋鸡），养殖量基础数据来源于2001—2015年河南统计年鉴、调查年鉴、河南农村统计年鉴等，其中生猪、肉牛、肉鸡为出栏量，能繁母猪、奶牛、蛋鸡为存栏量。在进行畜禽粪便的农田消纳量和农田承载负荷核算时，耕地面积应反映当年的种植情况，因此选用河南统计年鉴的农作物总播种面积指标。

畜禽养殖产污核算以第一次全国污染源普查畜禽养殖业源产排污系数手册的日产污系数为核算基数，核算粪便、尿液、TP、TN、COD、Cu、Zn等污染物的产生量，核算模型为：

(1)

式中：P_k为污染物产生量；L_i为第i种畜禽养殖量；θ_ijk为第i种畜禽第j养殖期第k种污染物的日产污系数；T_ij为第i种畜禽第j饲养周期。

饲养周期核算：生猪保育期70 d、生猪育肥期95 d、母猪365 d、奶牛育成期1.5 a、奶牛产奶期4.5 a、肉牛养殖期1 a、蛋鸡育雏育成期140 d、蛋鸡产蛋期1 a、肉鸡养殖期55 d。在畜禽养殖产污污染物核算过程中，将畜禽养殖周期按照统计年鉴数据统计周期（1 a）内的养殖天数进行折算。

借鉴杨世琦等^[10]的相关研究，引入畜禽粪便的农田消纳量（一定区域农田单位面积分摊的畜禽粪便产生数量或氮、磷产生数量）和农田畜禽粪便承载负荷（一定区域的畜禽粪便的农田容纳量与畜禽粪便安全消纳容量的比值）, 分析河南畜养产污粪便的农田消纳及承载负荷，相关计算模型如下：

(2)

式中：GAPF为农田消纳量；M为畜禽粪便量或氮磷产生量；S为耕地面积或播种面积。本研究中选取河南省作物总播种面积。

(3)

式中：MCL为承载负荷；SC为畜禽粪便安全消纳容量。当MCL取值为0、1和大于1分别代表零负荷、满负荷和超负荷。SC取值标准参考欧盟标准，氮170 kg·hm^-2，磷35 kg·hm^{-2 [17]}。

根据2000—2014年河南省畜禽养殖产污核算结果（表 1），河南省畜禽养殖污染物产生总量逐年增加，2000年为最低值，2014年为最高值，中间略有波动。为横向对比，本研究借助Z标准化方法消除量纲影响，将畜养产污核算值转换成均值为0，方差为1的数列，转换后的时间序列数据如图 1。通过2000—2014年河南省畜养产污Z标准化结果，整个时间段内河南畜禽养殖产污总量呈增加趋势，但明显以2006年为界分成两个阶段，即阶段Ⅰ（2000—2005）和阶段Ⅱ（2006—2014），与杨飞^[18]、付强^[19]等对全国畜养规律的研究一致。

表 1 2000—2014年河南畜养产污量核算（万t） Table 1 LPWG in Henan during 2000—2014

表选项

图 12000—2014年河南畜养产污量Z标准化 Figure 1Z-Score of LPWG in Henan during 2000—2014

图选项

2006年之前的畜牧业全面发展阶段（1996—2006年），畜牧业发展成为中国农业发展的重点，畜牧业养殖方式向规模化、专业化和集约化转变，依靠科技加快畜牧业结构优化和产业化经营步伐^[18]，期间中国规模化畜禽养殖处于慢速增长阶段^[19]。畜牧业现代化发展阶段（2007年—至今），畜牧业的综合生产能力和整体科技水平得到显著增强和提高，畜牧业向技术集约型、资源高效利用型、环境友好型转变，畜产品有效供给和质量安全得到保障，畜牧业发展的生态环境得到有效管理控制，总体上向现代畜牧业生产体系发展^[18]，期间中国规模化畜禽养殖处于快速增长阶段^[19]。

粪便产生量，在阶段Ⅰ，从4 565.11万t增加到6 046.98万t，增加约32.5%，年均增加6.5%。在阶段Ⅱ，从5 151.72万t增加到6 307.53万t，增加约22.4%，年均增加2.8%。从2000—2014年整个时间段看，河南省畜养粪便产生量增加约38.2%，年均增加约2.7%。

尿液产生量（表 1），在阶段Ⅰ，从4 427.25万t增加到5 967.17万t，增加约34.8%，年均增加7%。在阶段Ⅱ，从5 019.60万t增加到6 165.88万t，增加约22.8%，年均增加2.9%。从2000—2014年整个时间段看，河南省畜养产污尿液产生量增加约39.3%，年均增加约2.8%。

COD产生量，在阶段Ⅰ，从915.37万t增加到1 231.86万t，增加约34.6%，年均增加6.9%。在阶段Ⅱ，从1 055.15万t增加到1 315.67万t，增加约24.7%，年均增加3.1%。从2000—2014年整个时间段看，河南省畜养产污COD产生量增加约43.7%，年均增加约3.1%。

总氮（TN）产生量，在阶段Ⅰ，从50.77万t增加到71.31万t，增加约40.5%，年均增加8.1%。在阶段Ⅱ，从62.57万t增加到80.67万t，增加约28.9%，年均增加3.6%。从2000—2014年整个时间段看，河南省畜养产污TN产生量增加约58.9%，年均增加约4.2%。

总磷（TP）产生量，在阶段Ⅰ，从7.93万t增加到11.16万t，增加约40.7%，年均增加8.1%。在阶段Ⅱ，从9.61万t增加到12.60万t，增加约31.2%，年均增加3.9%。从2000—2014年整个时间段看，河南省畜养产污TP产生量增加约58.8%，年均增加约4.2%。

Cu产生量，在阶段Ⅰ，从0.10万t增加到0.14万t，增加约40.8%，年均增加8.2%。在阶段Ⅱ，从0.13万t增加到0.16万t，增加约25.4%，年均增加3.2%。从2000—2014年整个时间段看，河南省畜养产污Cu产生量增加约56.3%，年均增加约4.0%。

Zn产生量，在阶段Ⅰ，从0.30万t增加到0.43万t，增加约41.7%，年均增加8.3%。在阶段Ⅱ，从0.37万t增加到0.48万t，增加约29.0%，年均增加3.6%。从2000—2014年整个时间段看，河南省畜养产污Zn产生量增加约60.5%，年均增加约4.3%。

2000年以来河南畜养产污粪便的农田消纳及承载负荷核算结果如表 2。

表 2 2000—2014年河南畜养粪便的农田消纳及承载负荷 Table 2 GAPF and MCL of LPWG in Henan during 2000—2014

表选项

在阶段Ⅰ，河南省畜养产污TN的农田消纳值从2000年的38.65 kg·hm^-2提升到2005年的51.22 kg·hm^-2，增加约32.5%，年均增加约6.5%。相应的，其承载负荷从0.23提高到0.30。在阶段Ⅱ，河南省畜养产污TN的农田消纳值从2006年的44.71 kg·hm^-2提升到2014年的56.11 kg·hm^-2，增加约25.5%，年均增加约3.2%。相应的，其承载负荷从0.26提高到0.33。从2000—2014年整个时间段看，河南省畜养产污TN的农田消纳增加约45.2%，年均增加约3.2%，其承载负荷从0.23提高到0.33。

在阶段Ⅱ，河南省畜养产污TP的农田消纳值从2000年的6.04 kg·hm^-2提升到2005年的8.02 kg·hm^-2，增加约32.7%，年均增加约6.5%。相应的，其承载负荷从0.17提高到0.23。在阶段Ⅱ，河南省畜养产污TP的农田消纳值从2006年的6.86 kg·hm^-2提升到2014年的8.76 kg·hm^-2，增加约27.7%，年均增加约3.5%。相应的，其承载负荷从0.20提高到0.25。从2000—2014年整个时间段看，河南省畜养产污TP的农田消纳增加约45.1%，年均增加约3.2%，其承载负荷从0.17提高到0.25。

根据河南省畜养产污、农田消纳及承载负荷的计算结果，将阶段Ⅰ、阶段Ⅱ及全阶段各污染物增长率整理如图 2，有如下规律：

图 22000—2014年河南省畜养产污增长特征 Figure 2Growth characteristics of LPWG in Henan during 2000—2014

图选项

第一，各污染物及TN、TP农田消纳在阶段Ⅰ的增长率普遍高于各自在阶段Ⅱ的增长率，这种情况在年均增长率上更加显著（阶段Ⅰ年均7%左右，阶段Ⅱ年均3%左右）。阶段Ⅰ处于中国畜牧业全面发展阶段，规模化畜禽养殖慢速增长，在注重规模增长的同时忽视了专业化和污染物的处理，导致该阶段畜养产污的增长率偏高。阶段Ⅱ处于中国畜牧业现代化发展阶段，科技因素的引入带来了规模化、专业化的成熟发展，该阶段畜养产污得到一定治理，增长率较阶段Ⅰ有所下降。

第二，TN、TP、Cu、Zn在阶段Ⅰ和整个时段的增长率高于其他污染物，阶段ⅡTN、TP与Zn的增长率高于其他污染物。

第三，在整个时段，TN、TP、Cu、Zn增长在60%左右，粪便、尿液、COD、TN与TP的农田消纳增长40%左右。TN、TP年均增长约4%，其余年均增长约3%。

方法层面，本研究在进行污染物产生量核算时采用的核算系数，区分了不同畜种的养殖周期、每个畜种养殖过程的不同阶段，加上不同畜种在研究时段内的结构化变动，各个污染物按照系数单独核算，是一个综合的矩阵核算方法。规律二和规律三的不同污染物变化的差异性体现了这种综合计算的优势。

针对以上规律，在制定河南省域畜养产污应对策略时，可做以下考虑：

第一，河南省畜养产污阶段性明显，高增长的阶段Ⅰ已经过去，目前处于相对稳定增长的阶段Ⅱ，且持续时间较长。管理、治理政策的制定也应适应其当前阶段增长规律，持续一定时期。

第二，在治理畜养污染物时，有的放矢，主要瞄准增长率较高的TN、TP、Cu、Zn，重点是TN与TP的处理。氮磷是作物生长所需的重要养分，经过堆肥处理的畜禽粪便是高效的天然有机肥料。根据袁彩凤等^[12]的研究，2007年河南省粪便综合利用率为82%，堆沤还田利用占48%，双孢菇种植利用占25%，沼气化利用占9%。在现有畜禽粪便消纳方式的基础上，引导果园、菜地施用畜禽粪便有机肥，引导生态农场优先施用，形成大田为主、设施农业和生态农场为辅，其他方式补充的农田消纳模式。

第三，养殖量与作物种植的农田面积均有不同程度的变化，作物种植面积的增加能够在一定程度上缓解TN、TP的排放。当MCL取值为1时为满负荷，而按照目前的状况，河南省畜养粪便TN的最大MCL值为0.33，TP的最大MCL值为0.25。因此，只要河南省三分之一的作物种植面积用于畜养产污TN的消纳，四分之一的农田面积用于TP的消纳就足够了。总体上，引导三分之一的作物种植面积用于畜养产污的农田消纳可满足畜养产污资源利用需要。

需要说明的是，本研究是根据核算结果的理论化解释，未考虑到畜禽粪便清扫、堆积过程中的氮损失率，若考虑相关学者研究中给出的清扫和堆积中的氮损失率^[20]，河南省畜养产污的氮消纳会进一步降低。

经过河南省畜养产污核算、农田消纳及承载负荷计算、增长规律和对策分析等的研究，本研究探索了对接污染源普查畜种和污染物的省域畜养产污核算方法，以河南省为例的省域畜养产污研究也得到了一些可供相关研究人员和管理人员参考的研究结论，总结如下：

第一，本研究中的省域畜养产污核算方法，基于第一次全国污染源普查畜禽养殖业源的产排污系数，以不同畜种、不同养殖周期的天数为基准核算产污量，核算结果的基本规律与相关学者前期研究一致，能够为省域畜养产污尤其是河南省畜养产污规律的相关研究提供借鉴。

第二，河南省畜养产污核算结果表明，2000—2014年河南省畜养产污存在阶段性特点。2006年之前的阶段Ⅰ，各污染物的增长率普遍较快；2006年之后的阶段Ⅱ，各污染物的增长率相对较慢。

第三，本研究核算方法综合了各畜种养殖周期、养殖阶段和在研究时段内的结构化变动，能够探索出不同畜养污染物在研究期内的变化差异性，结果表明河南省畜养产污增长速度存在着差异性。从污染物增长率来看，TN、TP、Cu和Zn增长最快，在2000—2014年整个时段均超过其他污染物。

第四，现阶段，河南省畜养产污TN和TP的农田消纳和承载负荷压力不大，在不考虑土壤环境自身的N、P含量情况下，以全部畜养污染物产生量为农田消纳的核算基准，则引导三分之一的作物种植面积用于畜养产污的农田消纳就能够满足畜养产污资源利用的需要。考虑到土壤养分情况、交通运输成本等因素，河南省畜养污染物的消纳需要政策引导和空间配额优化。