近年来，畜禽养殖业快速发展，而大多数畜禽粪污并没有得到有效处理或资源化利用，直接或间接排入周边环境，对土壤、地表水甚至地下水造成污染，致使耕地土壤中重金属超标，引发较高的生态环境风险。刘洋^[1]研究表明，畜禽粪便中含有的大量N、P等元素进入水体中后会促进藻类的生长，使水中溶解氧含量降低，导致水体富营养化；张丽君等^[2]研究表明畜禽粪便中含有大量钾盐和钠盐，过量的钾和钠排放到农田中会导致土壤发生板结；据Liu等^[3]估算，2017年我国通过畜禽粪便排入农田环境的重金属（Zn、Cu、Pb、Cd、Cr、Hg、As和Ni）总量为2.86×10⁵ t（干质量），其中猪粪贡献占比高达71.52%。因此，开展畜禽养殖污染耕地承载负荷估算及风险评价，对于科学规划种养规模和布局、合理推进畜禽粪污资源化利用、防止农田环境污染具有重要意义。

国内学者已经针对典型地区畜禽粪污耕地承载负荷开展了相关研究。易秀等^[4]以陕西省为例，通过日排污系数法测算畜禽粪便负荷量，以市域尺度对陕西省畜禽养殖环境进行风险评价；张英等^[5]为评估河南省畜禽养殖环境污染现状，利用日排污系数法等方法估算了河南省畜禽养殖粪尿负荷量，并评价了不同地区畜禽养殖现状的合理性及对环境的污染状况。以上研究均假定畜禽养殖粪污全部被耕地消纳，未考虑粪污综合利用率。根据国务院办公厅印发的《关于加快推进畜禽养殖废弃物资源化利用的意见》，到2020年，全国畜禽粪污综合利用率要达到75% 以上^[6]。惠州市2018年的粪污综合利用率低于75%，因此在进行畜禽养殖污染耕地承载负荷估算时必须考虑粪污综合利用率才能更为准确地反映当地实际情况。

本研究在全面分析惠州市各区（县）养殖现状的基础上，根据畜禽产污系数和粪污综合利用率估算了2018年惠州市各区（县）畜禽粪污及主要污染物（COD、NH₃-N、TN、TP）的产生量和排放量，根据各区（县）实际耕地面积估算粪污、TN和TP的耕地承载负荷，并进行风险评价，进而结合分析情况和估算结果针对性地提出对策建议，旨在为惠州市畜禽养殖和农业种植的科学规划、畜禽粪污的合理利用、畜禽养殖污染的有效防治、东江水环境保护等工作提供科学依据，同时也为其他地区开展类似研究和有关工作提供参考。

1 材料与方法 1.1 研究区概况

惠州市地处珠江三角洲，位于粤港澳大湾区东岸，下辖惠城区、惠阳区、惠东县、博罗县、龙门县以及大亚湾经济技术开发区和仲恺高新技术产业开发区（图 1）。流经惠州市的东江，作为珠江三角洲地区重要的饮用水源，为惠州、东莞、广州、深圳和香港3 000多万人口提供了必不可少的生产、生活和生态用水。惠州市畜牧业发达，其中博罗县为养殖大县。若区域内产生的畜禽粪污没有得到有效处理或资源化利用，将会对东江水环境造成较大威胁。近些年惠州市种养情况调整较大，各区（县）的畜禽养殖现状、畜禽粪污及其中污染物产排情况并不完全清楚。因此，有必要系统分析惠州市各区（县）畜禽养殖污染现状，并对各区（县）耕地承载负荷进行估算和风险评价。

1.2 数据来源

以惠州市惠城区、惠阳区、惠东县、博罗县、龙门县5个区（县）的生猪、奶牛、肉牛、蛋鸡、肉鸡为研究对象，以各区（县）为基本单位，根据最新的相关统计数据（含第二次全国污染源普查数据），全面分析评估各区（县）畜禽养殖现状及其耕地承载负荷。依据2019年《广东农村统计年鉴》，大亚湾经济技术开发区和仲恺高新技术产业开发区分别并入惠阳区和惠城区统计。

惠州市各区（县）2018年末耕地面积，2018年生猪、奶牛、肉牛、蛋鸡、肉鸡年末存栏数和年出栏数，均来源于2019年《广东农村统计年鉴》^[7]和2019年《惠州统计年鉴》^[8]（表 1）。考虑到各区（县）全年的畜禽养殖量是动态变化的，畜禽年出栏量或年末存栏量都难以准确反映全年养殖量，故为了最大限度地反映畜禽粪污对生态环境的影响，选取年出栏量或年末存栏量中较大的一项作为全年养殖量。经比较后，选取了奶牛、肉牛、蛋鸡年末存栏量和生猪、肉鸡年出栏量进行后续估算。估算时采用的畜禽粪污产污系数及其中主要污染物产污系数，来源于第二次全国污染源普查中南区畜禽养殖饲养期产污系数（广东省被划为中南区）。

1.3 计算方法 1.3.1 畜禽粪污产生量和排放量估算

畜禽粪污产生量计算公式：

(1)

式中：Q为畜禽粪污产生量，t；N为畜禽养殖量，头或羽；P为畜禽饲养期产污系数，kg·头^-1或kg·羽^-1；f为单位转换因子，取值为10^-3。

畜禽饲养期产污系数为畜禽在整个饲养期内产生的所有粪污量。由于第二次污染源普查将畜禽饲养期产污系数分为规模化和非规模化，因此所采用的畜禽饲养期产污系数根据惠州市各类畜禽规模化与非规模化养殖量占比进行折算，见表 2。

实际生产过程中，不是所有畜禽粪污都被综合利用，未利用的部分直接排入受纳环境或者进行其他处理。对惠州市第二次全国污染源普查结果进行统计分析，得出2018年惠州市畜禽粪污综合利用率为70.6%，具体如下：

(2)

式中：E为畜禽（生猪、奶牛、肉牛、蛋鸡、肉鸡）粪污排放量，t。

1.3.2 畜禽养殖主要污染物产生量和排放量估算

畜禽粪污中的污染物主要有COD、NH₃-N、TN和TP，畜禽粪污中各类污染物产生量估算公式如下：

(3)

式中：W为畜禽粪污中污染物（COD、NH₃-N、TN、TP）产生量，kg；C为畜禽粪污中某污染物的饲养期产污系数，kg·头^-1或kg·羽^-1。由于第二次污染源普查将畜禽主要污染物饲养期产污系数分为规模化和非规模化，因此所采用的畜禽主要污染物饲养期产污系数根据惠州市各类畜禽规模化与非规模化养殖量占比进行折算，见表 3。

同理，畜禽粪污排入受纳环境的污染物排放量也应扣除被综合利用的部分（70.6%），具体估算公式如下：

(4)

式中：V为畜禽（生猪、奶牛、肉牛、蛋鸡、肉鸡）粪污主要污染物（COD、NH₃-N、TN、TP）排放量，kg。

1.3.3 畜禽养殖污染耕地承载负荷评价

由于相关文献中只获得了耕地畜禽粪污、TN、TP最大理论承载量。故本研究只选择畜禽粪污、TN、TP为参照进行单位面积耕地承载负荷评价，估算公式如下：

(5)

式中：L为畜禽污染物（粪污、TN、TP）耕地承载负荷，t·hm^-2或kg·hm^-2；S为年末实际耕地面积，hm²。

采用畜禽养殖污染负荷警报值评估研究区域内畜禽养殖污染对环境造成的压力以及潜在环境风险，估算公式如下：

(6)

式中：R为耕地负荷警报值，分级方法参考文献[9]，具体见表 4；M为各类畜禽养殖污染最大理论承载量，参考文献[10]，本研究以30 t·hm^-2、170 kg·hm^-2、35 kg· hm^-2作为耕地对畜禽粪污及其中TN、TP的最大理论承载量。

2 结果与讨论 2.1 惠州市畜禽养殖现状分析

根据统计数据（表 1）可知，惠州市畜禽养殖量较大，主要养殖品种为生猪、奶牛、肉牛、蛋鸡、肉鸡。2018年全市年出栏生猪1.38×10⁶头，奶牛年末存栏量5 500头，肉牛年末存栏量3.73×10⁴头，蛋鸡年末存栏4.43×10⁶羽，年出栏肉鸡3.65×10⁷羽。按照2018年农业部发布的《畜禽粪污土地承载力测算技术指南》的相关系数（100头猪相当于15头奶牛、30头肉牛、2 500只家禽）^[11]，将畜禽养殖量统一转化为猪当量以方便比较各区（县）、各畜禽的养殖情况，结果见图 2。2018年全市养殖总量（猪当量）为3.18×10⁶头，但各畜禽之间及各区（县）之间的养殖量分布差异较大。生猪、奶牛、肉牛、蛋鸡、肉鸡对全市畜禽养殖总量的贡献占比分别为43.4%、1.1%、3.9%、5.6%、46.0%，可见，肉鸡养殖量最多，为1.46×10⁶头猪当量；养殖量最少的是奶牛，为3.7×10⁴头猪当量。从空间分布来看，博罗县畜禽养殖量最大，为1.46×10⁶头猪当量，占全市的45.9%，是名副其实的养殖大县；惠东县和惠城区的畜禽养殖量居中，分别为7.11×10⁵、5.96×10⁵头猪当量，占全市的22.4%、18.7%；龙门县和惠阳区畜禽养殖量较少，分别占全市的7.3% 和5.7%。一般情况下，各区（县）畜禽养殖情况主要与当地产业结构、生产习惯以及禁养区、限养区划定等密切相关。针对惠州市畜禽养殖空间分布不平衡的现状，应建立总体协调机制，合理规划惠州市各区（县）禁养区、限养区和适养区，鼓励大规模的现代化畜禽养殖，控制小规模畜禽养殖场的养殖量，同时加大资金投入帮助小规模畜禽养殖场完成转型升级。

2.2 惠州市畜禽粪污产生量和排放量

根据2018年惠州市各类畜禽的养殖量、饲养期产污系数和粪污综合利用率，估算出惠州市2018年畜禽粪污产生量和排放量，结果分别见图 3和表 5。由于统一采用70.6%作为畜禽粪污综合利用率，故各区（县）、各畜禽对全市畜禽粪污产生量与排放量的贡献排序及其占比是相同的。2018年惠州市全市畜禽粪污产生量达1.51×10⁶ t，排放量达4.44×10⁵ t，各区（县）畜禽粪污产生量和排放量分布差异较大，贡献占比排序依次为博罗县（41.9%）>惠东县（25.3%）>惠城区（18.9%）>龙门县（8.1%）>惠阳区（5.8%）。博罗县作为养殖大县，畜禽粪污排放量最大，为1.86×10⁵ t；惠阳区排放量最小，为2.56×10⁴ t。不同畜禽之间粪污产生量和排放量差异较大，贡献率依次为生猪（50.3%）>肉鸡（17.2%）>肉牛（14.3%）>蛋鸡（13.4%）>奶牛（4.8%），生猪的粪尿排放量最大，为2.23×10⁵ t；奶牛的粪尿排放量最小，为2.16×10⁴ t。为了减少惠州市畜禽粪污排放量，应减少化肥使用，加强对畜禽粪污处理及资源化利用的研究及工程示范，探索出适合惠州市畜禽养殖特点的粪污处理处置模式，使畜禽粪污得到有效的处理及合理的资源化利用，既增加了土壤肥力，有效提高农作物产量，也能提高畜禽粪污经济效益，降低粪污资源化利用成本，减少畜禽粪污对生态环境的污染。

2.3 惠州市畜禽养殖污染物产生量和排放量

根据表 1畜禽养殖量和表 3畜禽粪污主要污染物饲养期产污系数，计算得出2018年惠州市畜禽养殖污染物产生量及排放量，结果分别见图 4和表 6。2018年惠州市畜禽养殖4种主要污染物（COD、NH₃-N、TN和TP）的产生量为2.61×10⁸kg，各区（县）贡献占比排序依次为博罗县（42.9%）>惠东县（24.4%）>惠城区（18.8%）>龙门县（7.8%）>惠阳区（6.1%）；养殖大县博罗县的畜禽养殖污染物产生量和排放量均最大，分别为1.12×10⁸kg和3.29×10⁷ kg。全市畜禽养殖4种污染物排放量为7.66 ×10⁷ kg，其中COD排放占比为92.9%，而TN、NH₃-N和TP排放量占比较低。可见，各类畜禽粪污中有机污染物含量较高，产生和排放的COD负荷较大，为畜禽粪污中的主要污染因子，假如处理利用不当，可能会造成周边农田土壤和水环境污染。而惠州市降水多，河流众多，畜禽污染物可能已经对境内的珠江水系造成了一定的污染，亟需加强畜禽粪污处理及资源化利用监管，建立健全畜禽养殖污染防治有关奖惩制度。对于规模畜禽养殖场，必须严格遵守环境影响评价制度、“三同时”制度、排污收费制度等环保制度；对于积极响应国家政策、配套粪污处理设备、合理进行粪污资源化利用的畜禽养殖场，应给予一定的资金补贴；对于偷排乱排养殖场，给予相应资金处罚，督促养殖场配套资源化利用设施，情节严重的养殖场应予以关停或者拆除。

2.4 惠州市畜禽养殖污染耕地承载负荷评价

目前，还田利用是畜禽粪污的主要处理方式之一，耕地是承载畜禽粪污的主要受纳场所。耕地环境具有一定的自净能力，可消纳一定负荷的畜禽粪污，耕地承载负荷这一指标可以直接反映区域耕地对畜禽养殖污染的消纳能力。为了科学评估惠州市耕地对畜禽粪污的消纳能力以及养殖污染对周边环境造成的影响，本研究估算了2018年惠州市各区（县）耕地对畜禽粪污、TN和TP的承载负荷，并根据相应的警报值进行环境风险评价，结果见表 7和图 5。

从畜禽粪污来看，全市畜禽粪污耕地承载负荷为7.69 t·hm^-2，警报值为0.26，属于Ⅰ级，对周边环境无威胁。其中惠城区的畜禽粪污耕地承载负荷最高，为9.31 t·hm^-2，博罗县次之，为8.97 t·hm^-2，但所有区（县）的畜禽粪污警报级别均为Ⅰ级，无潜在环境风险。对于TN来说，全市TN耕地承载负荷为65.38 kg· hm^-2，警报值为0.38，属于Ⅰ级，对周边环境暂无威胁。而惠城区和博罗县TN耕地承载负荷较高，警报级别为Ⅱ级，对周围环境稍有威胁。对于TP来说，全市TP耕地承载负荷为15.63 kg·hm^-2，警报级别属于Ⅱ级，潜在环境风险低，对周边环境稍有威胁。除了惠阳区和龙门县，其余三个区（县）的TP承载负荷较高（15.30~19.62 kg·hm^-2），警报级别均为Ⅱ级，潜在环境风险低，对环境稍有威胁。

综上，目前惠州市畜禽粪污及主要污染物排放量较小，粪污、TN、TP耕地承载负荷较低，对周围环境威胁小，潜在环境风险低。与其他区（县）相比，博罗县、惠城区耕地承载负荷较高，是因为近年来深圳市、东莞市以及惠阳区部分养殖场迁入。而全市耕地承载负荷并不高，这与全市开展禁养区划定、拆除不合格养殖场，同时严格要求畜禽业配套环保措施有关。从长远来看，惠州市应制定全市养殖业和种植业发展规划，实现区域种养平衡，并根据各区（县）对畜禽粪污的最大土地消纳能力，合理地调整各区（县）畜禽养殖规模和作物种植面积。对于耕地承载负荷较高的区域，应减少畜禽养殖量，扩大作物种植面积，加强粪污资源化利用，实现区域种植业和养殖业在空间上和时间上有机结合、平衡发展，防止畜禽粪污对当地生态环境造成污染和破坏，保障东江饮用水安全。

3 结论

（1）2018年惠州全市养殖总量（猪当量）为3.18× 106头，各区（县）占比排序为博罗县（45.9%）>惠东县（22.4%）>惠城区（18.7%）>龙门县（7.3%）>惠阳区（5.7%）。

（2）2018年惠州市畜禽粪污产生量达1.51×10⁶ t，排放量达4.44×10⁵ t，各畜禽贡献占比排序为生猪（50.3%）>肉鸡（17.2%）>肉牛（14.3%）>蛋鸡（13.4%）>奶牛（4.8%）。

（3）2018年惠州市畜禽养殖污染物产生量为2.61×10⁸kg，排放量7.66×10⁷ kg，各污染物贡献占比排序为COD（92.9%）>TN（4.9%）>TP（1.2%）>NH₃-N（1.0%）。

（4）2018年惠州市畜禽粪污和TN耕地承载负荷分别为7.69 t·hm^-2和65.38 kg·hm^-2，警报级别均为Ⅰ级，对周围环境无威胁；而TP耕地承载负荷为15.63 kg·hm^-2，警报级别为Ⅱ级，对周边环境稍有威胁。

（5）综合畜禽粪污、TP、TN耕地承载负荷情况，目前惠州市畜禽养殖量对周围耕地环境压力较小，可以通过制定种养平衡规划提升畜禽粪污资源化利用水平。