自改革开放以来，随着我国人民生活水平的提高和对动物蛋白摄入的需求加大，居民对肉类的消费量大幅增加，我国畜禽养殖业也因此得到迅速发展^[1-3]。自1991年以来，我国肉、禽、蛋总产量已连续多年保持世界第一^[4]，并在2016年达到13 523万t，其中肉、禽和蛋类总产量分别为8540万、1888万、3095万t^[5]。与此同时，未经有效处理的畜禽排泄物进入到环境中带来的水体污染问题也越来越严重^[6-10]。畜禽粪便、尿液中含有大量的氮、磷等营养元素，不恰当的处理方式会导致相当一部分营养元素进入水体，造成水环境污染^[11]。据陈瑶等^[12]统计，在2007年，全国畜禽养殖业粪便和尿液产生量分别达到2.43亿t和1.63亿t，由此带来的化学需氧量、氮元素和磷元素年排放量分别占全国污染物排放总量的42%、22%和38%；Gu等^[13]估算，在2010年，全国由畜禽养殖排放到水体中的氮元素通量达到300万t，是所有农业源排放的重要组成；孙良媛等^[3]估算出2011年畜禽养殖的磷元素入水通量约为40万t。由此可以看出，庞大的畜禽养殖已经给我国水环境带来了巨大的压力^[14-15]。

近十年来，我国畜禽养殖在数量上稍有增加，如2006年我国生猪年末出栏量为68 050万头，在2012年增加到69 790万头^[16]。同时，随着农业生产方式的改变，我国畜禽养殖模式也逐步由家庭散养向规模化养殖过渡^{[3, 17-19]}。2010年我国生猪出栏量500头以上的小型养殖户的数量占全国总生猪养殖户的0.4%，而其生猪出栏量却占到了34.5%^[16]。根据畜禽养殖统计数据，在2006—2010年间，全国畜禽养殖中猪、肉牛、奶牛、羊、蛋鸡和肉鸡的规模化养殖数量分别提高了28%、10%、26%、6%、17%和23%，其中猪和肉鸡的规模化养殖比例已经超过了50%。但是需要指出的是，畜禽养殖模式的转变在我国不同区域内的进展程度仍然存在较大差异。比如在2010年，北京和上海地区的蛋鸡规模化养殖比例已经分别达到67%和65%；而西部地区规模化进程则较慢，如内蒙古自治区蛋鸡规模化养殖比例仅为14%^[16]。

当前国内外学术界关于畜禽养殖环境影响问题的研究成果颇丰，国外学者多运用自然环境物理模型模拟和预测污染物的产生、迁移和危害等^{[3, 9]}。而当前国内对于畜禽养殖污染物排污量计算的方法主要采用排放因子法，现有研究对污染物排放因子多着重于定性描述^{[3-5, 13]}，通过选取不同畜禽产排污系数和粪便污染物含量参数，估算不同畜禽类别的污染物产生总量，再通过污染物的水体流失率估算畜禽养殖业的污染物入水通量，从而对畜禽粪便产生的环境效应进行分析。从已有研究来看，对畜禽养殖业氮元素污水通量的计算方法差别不大，但大部分研究只集中于个别年份，对全国各个省份的氮排放通量进行计算^{[3-5, 13]}，没有在时间和空间尺度上分析我国畜禽养殖业氮元素入水通量的变化特征，也没有将畜禽规模化养殖和散养排污分别进行分析。因此，为了评估近十年来我国畜禽养殖氮元素排放通量的变化以及养殖模式转变带来的影响，本文利用2006—2016年中国畜牧业统计数据、畜禽粪便含氮量、水体流失率等统计资料或文献数据，以畜禽养殖业中最主要的猪、羊、肉牛、奶牛、蛋鸡和肉鸡为调查对象，估算了分别由散养和规模化养殖向水环境中排放的氮元素通量，着重比较了不同地区的空间分布特征和不同畜禽间的差异。本研究的开展可以为评估畜禽养殖对我国水环境的影响和今后相应环境管理政策制定提供科学依据。

1 材料与方法 1.1 数据来源

本文涉及的主要数据包括：各类畜禽产排污系数、排泄物水体流失率、2006—2016年间我国31个省份和不同类型畜禽（包括羊、猪、肉牛、奶牛、蛋鸡和肉鸡）出栏量、存栏量、规模化养殖比例等^{[3, 20-22]}，养殖规模较小的畜禽不予以考虑。规模化养殖的畜禽粪便和尿液中氮元素水体流失率采用国家环保总局南京环科所测定出的数据，如表 1所示，各类规模化养殖畜禽粪便的氮元素水体流失率大约在5%~9%，而尿液氮元素流失率则达到50%，由于现有研究主要集中于规模化畜禽养殖的污染物流失率，散养畜禽排泄物的氮元素流失率统一采用30%^{[13, 21]}。利用《第一次全国污染源普查畜禽养殖业产排污系数与排污系数手册》，将主要畜禽的产污系数划分为五个区域取值，这五个地区具体包括华北地区（北京、天津、河北、山西、内蒙古、辽宁、吉林、黑龙江）、华东地区（上海、江苏、浙江、安徽、福建、江西、山东）、中南地区（河南、湖北、湖南、广东、广西、海南）、西南地区（重庆、四川、贵州、云南、西藏）和西北地区（山西、甘肃、青海、宁夏、新疆）。

1.2 畜禽规模化养殖界定以及氮元素入水通量计算

对于规模化养殖的具体标准，不同统计资料或行业标准间存在一定差异，综合《畜禽养殖业污染治理工程技术规范》《全国农产品成本收益资料汇编》、《中国畜牧业统计年鉴》信息^{[16, 23-24]}，本文采用的规模化标准如表 2所示。

畜禽排泄物产生的氮元素入水排放总量计算公式为：

(1)

(2)

(3)

式中，Q_1i、Q_2i和Q分别为第i种畜禽规模化养殖、散养以及所有畜禽总的氮元素入水通量（kg·a^-1）；N_i存、N_i出为第i种畜禽的年存栏量和出栏量；P_i为规模化养殖比例；E_i、N_Ei、L_Ei分别为畜禽的单位粪便量（kg·只^-1·d^-1）、粪便含氮量（表 1）、粪便水体流失率；U_i、N_Ui、L_Ui分别为畜禽的单位尿液量（kg·只^-1·d^-1）、尿液含氮量、尿液水体流失率；L_N为散养的粪便水体流失率。由于2011年以后畜禽规模化养殖数据缺失，本文利用2006—2010年的统计数据，采用线性回归方法来估算不同畜禽在各地区的规模化养殖比例。为了纠正出栏动物不再产生粪污这一事实，我们用因子0.542校正了一年中被杀死动物的数量，表示出栏动物是年饲养活畜的一部分，以此来估算出栏畜禽年实际排污量^[25-26]。

2 结果与讨论 2.1 2006—2016年全国畜禽养殖业氮元素入水通量

如图 1所示，在2006—2016年间，全国由畜禽养殖排入水体中的氮元素平均通量为4384 Gg，这11年中氮元素入水通量的最大和最小值分别为4832 Gg（2006年）和4010 Gg（2008年），而2006至2008年也是畜禽养殖业氮元素入水通量下降最为显著的时期。2016年全国畜禽养殖氮元素排放通量达到4426 Gg，相比于2006年下降了8%，但在2012年以后的氮元素入水通量几乎保持不变。Gu等^[13]和孙良媛等^[3]分别估算出的2010年和2011年畜禽养殖氮元素入水通量为3000 Gg和2140 Gg，本文的估算结果明显高于前两者，由于本文将畜禽规模化养殖和散养的粪便排污量分别进行计算，规模化养殖畜禽粪便的氮元素水体流失率远低于散养畜禽粪便氮元素水体流失率，而在这11年中，我国有相当一部分地区的某些畜禽仍然以散养为主，因此，畜禽散养和规模化的粪便水体流失率差异是造成本文结果明显高于前者计算结果的主要原因，但本文估算结果与《全国规模化畜禽养殖业污染情况调查技术报告》 ^[5]的估算结果（4071 Gg）相近。总体来看，尽管近十年中我国畜禽规模化养殖比例提升显著，但是入水氮元素通量下降并不显著。分析原因可能有两点：一是随着人口增加以及生活水平提高，国内肉类消费量逐年上升，同时活畜禽和鲜冻肉类的出口量也较大^{[16, 27]}，因此畜禽养殖量仍处于较高规模；二是规模化养殖也可能造成污染物的集中排放^[28]。规模化养殖对粪便集中收集处理效率较高，但是对畜禽尿液的处理并不理想^[29]。规模化养殖场的清粪方式主要有三种：水冲清粪、干清粪和垫草垫料清粪，其中水冲清粪的排污系数约为干清粪的3倍^[20]，垫草垫料清粪最为清洁，但水冲清粪造价低，操作方便，因此国内很多地方仍采用该方法来清理畜禽粪便^{[19, 29]}，如祝其丽等^[29]调查了2011年全国144处不同规模猪场的清粪方式，结果表明小型养猪场（500~ 2999头）水冲清粪的比例高达87%；刘永丰等^[19]对2014年海南省各地区畜禽规模化养殖场的调查结果显示全省各地区至少有60%的养猪场还在使用水冲清粪工艺。综合各省份畜禽养殖业氮元素入水通量，各省份间排放通量差异显著。2006年畜禽养殖业氮元素入水通量最多和最少的省市分别是河南（516 Gg）和上海（148 Gg），分别占该年氮元素总入水通量的11%和0.3%，到2016年，畜禽养殖业氮元素入水通量最多和最少的仍是河南（435 Gg）和上海（118 Gg），氮元素入水通量与10年前相比均有所下降，各自占比分别为10%和0.3%。除了河南省之外，山东省、四川省对于畜禽养殖业氮元素入水通量的贡献也较大，2016年的氮元素入水通量分别是359 Gg和298 Gg，占比分别为8%和7%。

在2006—2016年间，不同畜禽养殖类型的氮元素入水通量变化不大（图 2a）：在2006年和2016年，羊、猪、肉牛、奶牛、蛋鸡、肉鸡养殖的氮元素入水通量分别为31、163、177、41、28、42 Gg和22、164、138、59、22、36 Gg；从比例上看（图 2b），2006年和2016年上述各种畜禽的氮元素入水通量分别为7%、34%、36%、8%、6%、9%和5%、37%、31%、13%、5%、8%。总体来看，氮元素入水通量按猪>肉牛>奶牛>肉鸡>羊>蛋鸡的顺序依次降低。因此可以看出，我国畜禽养殖业中猪和肉牛是导致氮元素进入水体的主要排放源。有研究表明^[28]，在我国当前的畜禽养殖模式下，育肥猪贡献了约有一半的粪便污染物（包括氮元素、磷元素以及COD等）；另外，肉牛的单位产污量最大，约为猪肉的2~3倍，鸡肉的5~20倍^[28]。本文的结果也表明，若要有效控制我国氮元素入水通量，应密切关注畜禽养殖过程中的氮排放，而在畜禽养殖中猪和牛则是首要的控制因素^{[28, 30]}。

2.2 不同养殖方式的氮元素入水通量变化估算

由图 3所示，在2006—2016年间，全国由畜禽规模化养殖排入水中的氮元素通量呈现显著的上升趋势。在2006年，规模化养殖氮元素入水通量的比例占总入水通量的23%，在2016年，该比例上升到53%，而这种增加也可以显著反映我国规模化养殖比例的迅速增加。在我国不同的地理区域，畜禽散养和规模化养殖的氮元素入水通量的分布存在较大差异。如图 4所示，我们选择了2006、2010年和2016年为代表，比较了畜禽散养和规模化养殖的氮元素入水通量在我国空间分布上的差异。总体来看，中南和华北地区是我国的畜禽养殖氮元素排放的主要地区，华东和西南次之，西北地区最少。2016年华北、中南、华东、西南和西北地区的畜禽氮元素流失通量比例分别为27%、26%、20%、17%、10%。在所有主要区域中，在中南地区，我国畜禽散养氮元素入水通量下降趋势最为明显，由2006年的930 Gg下降至2016年的430 Gg，华北和华东地区次之，下降量分别为437 Gg和426 Gg。对于西南和西北地区，其畜禽散养带来的氮元素入水通量下降趋势并不明显，到2016年仍占相当高的比例，西南地区超过55%，西北地区超过70%。11年来，我国畜禽规模化养殖氮元素入水通量上升趋势都较为明显，中南地区最多，由2006年的340 Gg上升到2016年的744 Gg，华北、华东、西南、西北地区的增长量分别为234、326、154、111 Gg，到2016年，华东、中南地区的畜禽养殖业氮元素排放主要由规模化养殖引起，分别占各自地区氮元素入水通量的68%、63%，而华北地区的规模化养殖氮元素入水通量占50%，并有继续上升的趋势（图 3b）。出现这种差异的原因可能是我国西南和西北地区地广人稀，传统畜禽散养模式向规模化养殖的转换要远远慢于经济发展快、人口密度大的东南部地区。

2.3 畜禽类型和区域差异对氮元素入水通量的估算

不同畜禽类型（包括羊、猪、肉牛、奶牛、蛋鸡和肉鸡）通过散养和规模化养殖的氮元素入水通量变化如图 5所示。不同畜禽类型通过散养和规模化养殖带来的氮元素排放通量差异显著。具体来说，羊、肉牛和蛋鸡由散养向水中排放的氮元素量较大，在2006年和2016年的排放量分别为302、1499、270 Gg和201、930、172 Gg。在2016年，其排放在总排放中的比例分别高达91%、76%、67%。显著不同的是，猪和肉鸡由散养排放氮元素通量的比例下降明显，猪散养的氮元素入水通量在2006年和2016年分别为1042 Gg和347 Gg，占其总通量的64%和21%；肉鸡散养的氮元素入水通量在2006年和2016年分别为311 Gg和136 Gg，占其总通量的74%和38%。奶牛散养的氮元素入水通量在过去11年内几乎保持不变（均在295~297 Gg范围内），但是该部分的排放比例下降却较为显著，由2006的72%降至2016年的50%。这种差异也体现了不同畜禽养殖类型在规模化养殖比例中的差异。

不同区域的不同畜禽类型通过散养和规模化养殖的氮元素入水通量也有差异，如2016年（图 6），华北、华东、中南地区氮素入水量主要由猪规模化养殖排放，2016年猪规模化养殖排放的比例分别为36%、36%、47%，这些区域应该加强规模化猪养殖粪便污染配套设施建设。西南和西北地区氮素入水量主要由肉牛散养排放，2016年两区域的肉牛散养排放比例分别为39%和27%，因此应该加强散养肉牛的粪便污染控制。总体来看，猪规模化养殖和肉牛散养排放的氮元素在各个地区均占有较大比例。另外，由羊散养进入水体中的氮元素在西北地区的比例较其他地区高，该地区也应该加强散养羊的粪便污染控制。

3 结论

（1）2006—2016年，全国由畜禽养殖排入水体中的氮元素平均通量为4384 Gg，相比于2006年，2016年氮元素入水通量略有下降。规模化养殖氮元素入水通量所占的比例由2006年的23%上升到2016年的53%，因此，规模化养殖带来的污染物排放在未来水环境管理中需要优先考虑。

（2）过去11年内，我国中南、华北、华东、西南、西北地区规模化畜禽养殖的氮元素入水通量分别增长了404、234、326、154、111 Gg；在2016年，羊、肉牛、蛋鸡、奶牛、肉鸡和猪散养的氮元素入水通量分别占比91%、76%、67%、50%、38%、21%。

（3）对于不同畜禽养殖类型，氮元素入水通量按猪>肉牛>奶牛>肉鸡>羊>蛋鸡的顺序依次降低。2016年，华北、华东、中南地区氮素入水量主要由猪规模化养殖排放，排放比例分别为36%、36%、47%，西南和西北地区氮素入水量主要由肉牛散养排放，排放比例分别为39%和27%，总体来看，肉牛和猪成为我国畜禽养殖氮元素入水通量的主要排放源。