陶 玲，彭 亮，代梨梨，杨 镇，陈思媛，可 毅，李 谷
(中国水产科学研究院长江水产研究所，武汉 430223)
据《中国小龙虾产业发展报告(2021)》[1]，2020年我国小龙虾养殖总面积146万hm2，总产量239万t，其中小龙虾稻田养殖面积约为126 万hm2，养殖产量206多万t，分别占小龙虾养殖总面积和养殖总产量的86.61%和86.15%。江汉平原地区水土资源优越、气候适宜，稻田养殖小龙虾模式发展迅速，其中位于江汉平原腹地的湖北省作为小龙虾稻田养殖第一大省，养殖面积及其产值均居全国之首。虽然稻田养殖小龙虾模式目前被认为是绿色生态健康养殖模式，但该模式快速发展过程中也逐渐暴露出了水资源消耗大、耕地破碎化、土壤潜育化和面源污染潜在风险等问题[2-3]。未来稻田养殖小龙虾模式的可持续发展问题受到关注[4]。
目前，较多的研究表明，由于在稻田养殖小龙虾期间进行肥水培藻、投饲、用药导致残余的饲料、肥料以及粪便残留于稻田水体中，使小龙虾养殖期间稻田的田面水中氮、磷浓度显著提升[5-7]。稻虾轮作模式是传统的小龙虾稻田养殖模式，也是目前的主流模式，该模式普遍在5-6月将成熟小龙虾基本收获，并在极短时间内集中将稻田的田面水即小龙虾养殖阶段的养殖尾水从稻田内排放至周围沟渠管网，以进行整田和种植水稻，排出的水量远大于一般农田，易造成氮、磷等营养集中大量排放，且该模式面积较大，可能会加剧周边水域水体富营养化，破坏水体生态系统，需要引起农业和环境保护工作者的重视[8-9]。目前相关部门对池塘养殖、工厂化养殖等不同养殖模式不同养殖品种的排污系数进行了测算，先后发布了《第一次全国污染源普查水产养殖业污染源产排污系数手册》[10]和第二次全国污染源普查《排放源统计调查产排污核算方法和系数手册》[11]。进行稻虾轮作的稻田兼具农田和养殖塘的特点，既包含种植业也包含水产养殖业，在周年的农业生产中进行种植和养殖的转换，水稻种植过程中稻田面源污染排放系数已有估算[11]，但稻田小龙虾养殖阶段养殖尾水排放造成的面源污染程度不得而知，有必要对其排污系数进行估算。
经济效益高且对环境影响小的可持续发展养殖模式是未来水产养殖发展的方向和焦点。能值(emergy)理论及分析方法是美国著名生态学家、系统能量分析先驱Odum[12]于1980s创立的。依据能值这一新的科学概念和度量标准将生态经济系统内不同类别的物质和能量转换为统一标准的能值，进行定量分析研究，从而评价自然生产和人类经济活动，实现对生态系统和生产模式的环境效应和可持续性评价分析[13]。目前，能值分析方法已被应用于小龙虾养殖系统的可持续性评价。Chen等[4]利用能值分析方法评价了稻虾轮作、稻麦轮作和水稻单作3种模式的环境效应和可持续性。Hou等[14]也利用能值分析方法评价了稻虾轮作、稻虾共作和水稻单作3种模式的环境效应和可持续发展能力。环境污染排放是对系统能值产生的负向影响，会影响养殖模式的可持续发展性能，但以上研究对稻虾轮作模式的能值分析均未将其面源污染负值产出计算在内，不能反映系统面源污染物质排放的环境影响，得出的可持续发展指标不能体现其真实的可持续性。因此，有必要在获得稻虾轮作模式排污系数的基础上，对该模式进行包含面源污染排放的能值分析，定量评价面源污染排放对稻虾轮作模式可持续性发展的影响，以帮助政府、养殖户和环保部门进行养殖模式选择的决策分析。
因此，本研究在江汉平原稻虾轮作模式推广的核心区域，夏季小龙虾捕捞基本结束水稻种植前集中排水期间，通过野外采样获取江汉平原稻虾轮作模式小龙虾养殖排放的尾水中污染物浓度，并计算其排放系数，采用等标污染负荷法进行主要污染物解析。通过入户调研获取稻虾田排、换水以及养殖生产投入和产出情况等数据，运用能值分析法，计算得出包含面源污染负产出的能值评价指数，对江汉平原稻虾轮作模式的可持续发展能力进行定量评价，为最终摸清稻虾轮作模式面源污染排放量以及江汉平原农业面源污染的防治提供科学依据和决策参考，确保农业增效、农民增收的同时，保障稻虾轮作模式的可持续发展。
1.1 采样地点
采样和调查地点位于江汉平原，该区属北亚热带湿润季风气候，四季分明，雨量充沛，光照适宜，春秋季短，冬夏季长。多年平均气温15.9～16.6℃，无霜期一般为242～263 d。多年平均降水量1000～1300 mm，呈自北往南递增的趋势，东西相对均匀。区内主要河流有长江、汉江，次有沮漳河、西荆河、东荆河等，湖泊星罗棋布，渠系纵横，水系十分发达。
在湖北省潜江市和荆州市所辖的5个县(市、区)：包括洪湖市、公安县、江陵县、荆州区和沙市区，选择了能代表当地实际情况的家庭农场/合作社或农户的稻虾轮作田进行样品采集和调研(图1)。农户养殖规模为4669～46690 m2，家庭农场/合作社养殖规模为80040～800400 m2，共采集了8个农户和7个家庭农场/合作社的36块小龙虾养殖稻田的田面水，共计36个水样。样品采集的稻田均采用稻虾轮作模式：水稻一般于每年6月进行整地、插秧，10月收获。当年10-11月投放虾苗或种虾，至翌年3月幼虾在稻田中自然生长，3-5月投放饲料和补充虾苗投放，并进行成虾捕捞，6月下旬捕捞基本结束，再进行下一季的水稻种植，依次循环。需要说明的是，潜江市各采样点均采用了“虾稻轮作、共作一体”的养殖模式，即在稻虾轮作基础上增加了稻虾共作养殖，在水稻插秧前后继续投放小龙虾幼虾养殖或继续养殖留存在稻田中的小龙虾幼虾，7-8月捕捞第二季成虾，但此次仅采集了潜江各样点稻虾轮作阶段排放的养殖尾水，因此将潜江各采样点与其它县(市、区)采样点的小龙虾养殖排放尾水一起分析。
图1 采样点分布Fig.1 Distribution of sampling sites
通过采用入户访谈的方式获取了稻虾轮作模式生产的基本情况：包括种植/养殖面积、生产模式；
生产投入情况：稻种、虾苗、肥料、地笼和防逃设施、饲料、机耕/机收、雇工、水电等；
产出：水稻产量和小龙虾产量；
换水和补水方式、水量、时间和周期等。调研结果显示：各采样点农户采用自繁虾苗或自繁虾苗+补苗的方式养殖，主要靠经验进行田间水肥管理和投饲，投放饲料均为虾蟹专用料，小龙虾产量为1000～2700 kg/hm2(表1)。养殖水体主要来源于附近河道或沟渠，养殖期间基本不换水，仅进行补水，小龙虾养殖结束水稻种植前将稻田的田面水基本排干，直接排入附近河道或沟渠。
表1 稻虾轮作模式生产管理情况Tab.1 Survey of the rice-crayfish rotation mode
1.2 样品采集及分析方法
在2020年12月初，采集了稻田的田面水或进水沟渠渠水作为小龙虾养殖初期进水初始水样，2021年6月初即稻田养殖的商品小龙虾捕捞基本结束，水稻移栽前、小龙虾养殖尾水集中排放期间，采集稻田的田面水水样作为小龙虾养殖尾水排放水样。每个稻田采集5个点位的田面水进行混合后再取500 mL水样带回实验室进行测定，同时现场测定稻田的田面水深度。
水质测定指标包括总氮(TN)、总磷(TP)、化学需氧量(COD)和氨氮(NH3-N)4个指标，测定方法均参照《水和废水监测分析方法》[15]进行。
1.3 排污系数计算
根据《排放源统计调查产排污核算方法和系数手册》[11]，水产养殖排污系数指在正常养殖生产条件下，养殖生产1 t水产品所产生的污染物量中，经不同排放渠道直接排放到湖泊、河流及海洋等(不包括排放到农田及水产养殖再利用等部分)外部水体环境中的污染物量(kg/t)；
种植业氮磷排放(流失)系数指土壤和肥料中的氮磷在降雨或灌溉水作用下溶解或悬浮于径流水中，随径流迁移出田块而导致的农田氮磷流失的发生量(kg/hm2)。稻虾轮作模式既包含水稻种植又包含小龙虾养殖，由于本研究只采集了小龙虾养殖阶段排放尾水中氮磷数据，因此本文采用水产养殖排污系数方法计算小龙虾养殖阶段的排污系数。但为了比较不同稻田生产模式转换带来的面源污染排放差异，便于与江汉平原现有的稻麦轮作、稻油轮作和单季稻等稻田种植模式进行周年的面源污染负荷对比，还估算了稻虾轮作模式的单位面积排放系数。
因各采样点稻田养殖小龙虾过程中基本不进行换水，根据小龙虾捕捞结束、水稻种植前排放的养殖尾水中TN、TP、COD和NH3-N浓度值和排水量，以及初始补水浓度值和水量，参考高月香等[16]的方法，小龙虾单位养殖产量排污系数公式计算如下：
(1)
式中，Mi为第i种污染物的单位养殖产量排污系数，kg/t；
Q排j,i为养殖过程中的第j次排水所排放的水量，m3；
C排j,i为第j次排放到水体中第i种污染物的浓度，mg/L；
Q补j,i为养殖过程中第j次补水的水量，m3；
C补j,i为第j次补水水体中第i种污染物的浓度，mg/L；
ΔW为养殖生物净产量或质量增量，t。
单位面积TN、TP、COD 和NH3-N的排放系数(Pi)参考文献[17]，计算公式如下：
Pi=Q×ΔC×10-3
(2)
式中，Pi为第i种污染物的单位面积排放系数，kg/hm2；
ΔC为小龙虾养殖尾水排放浓度与初始补水浓度的差值，mg/L；
Q为水稻种植前每公顷稻田排放的小龙虾养殖尾水的排水量，m3。
1.4 等标污染负荷计算
等标污染负荷法用污染物的绝对排放量与排放标准进行比较，从而能将污染源污染物的排放量转化为“把污染物全部稀释到评价标准所需的介质量”中，使同一污染源所排放的污染物之间、不同污染源之间对环境的潜在影响进行比较成为可能[18]，其主要反映污染源本身潜在的污染水平，可用于评价区域内不同污染物排放量对环境的影响程度，确定主要污染源和污染物[19]。采用等标污染负荷法，对小龙虾养殖排放的尾水中TN、TP、COD和NH3-N等污染物进行数据标准化处理，计算得到相应的等标污染负荷值和等标污染负荷比，比较虾稻轮作模式小龙虾养殖尾水中不同污染物对环境的潜在影响大小。等标污染负荷相关公式计算如下：
Pi=Ci/C0i×Q
(3)
(4)
式中，Pi为污染物i的等标污染负荷，m3；
Ci为污染物的实测浓度值，mg/L；
Q为污染物的废水排放量，m3；
C0i为污染物i基于水环境功能分区的水质控制类别标准值(统一按照《地表水环境质量标准》(GB 3838-2002)[20]中Ⅲ类标准进行计算，COD 取值为20 mg/L，TN取值为1 mg/L，TP 取值为0.2 mg/L，NH3-N取值为1 mg/L)；
Ki为污染物i的等标污染负荷比。遵循污染源评价经验，按照污染物等标污染负荷比由大到小排列，分别计算其累计百分比，规定百分比累计到80%的污染物为主要污染物[18]。
1.5 包含面源污染排放的能值估算
1.5.1 数据来源 本研究中稻虾轮作模式小龙虾养殖结束排放的养殖尾水中面源污染排放数据来自本实验测定。模式经济投入与产出相关数据来自本次调研数据、文献和官方公开统计数据。所有数据均以对应于1 hm2农田进行折算，以一个完整的虾稻生产周期为界限。
1.5.2 能值分析 使用太阳能值转化率(UEV，sej/J或 sej/g)将各项物质、能量和资金数据转化为能值；
在能值核算中引入可更新比例系数(renewable fraction)，从而将投入资源分为可更新和不可更新两部分[21]。由于系统中投入的秸秆来源于稻田本身，在收获后又被归还于农田，为避免重复计算，秸秆投入在系统总能值投入中不计入[22]。采用2016 年最新修订的能值基准(12.00×1024sej/a)[23]作为计算基础。
1.5.3 能值评价指标 计算出各模式的本地可更新自然资源能值(LR)、本地不可更新自然资源能值(LN)、外部购买资源能值(F)、不可更新外部购买资源能值(FN)、可更新外部购买资源能值(FR)和系统的总投入能值(U)。根据 Odum[12]提出的能值评价指标体系，选取农业生态系统常用的能值指标：可更新资源投入占比(renewable fraction，%R)、能值产出率(emergy yield ratio，EYR)、环境负载率(environment loading ratio，ELR)和能值可持续指数(emergy sustainability index，ESI)作为评价系统可持续性的指标。相关计算公式及其相关解释意义如下：
(1)可更新资源投入占比(%R)
反映生产过程中所利用的可更新资源的比率。%R值越大，表明系统运作更多地依赖可更新资源的投入。
%R= 100 × (LR+FR)/U
(5)
(2)能值产出率(EYR)
为系统中投入的总能值与购买的外部资源能值总和的比值。EYR值越大，表明系统的生产效率越高，经济效益越好。由于面源污染排放对环境的影响表现为负面效应，参考钟珍梅[24]等，将面源污染排放能值设为负能值，公式优化为总投入能值减去面源污染能值之差与购买的外部资源能值总和的比值。
EYR=(U-P)/F
(6)
(3)环境负载率(ELR)
为系统不可更新资源投入总量与可更新资源投入总量之比，表示系统对环境造成的生态压力，值越大表示，对环境的压力越大[25]。由于面源污染排放对环境的影响是负面的，即对环境负载率的贡献为正值，因此指标的计算公式优化为不可更新资源投入能值和面源污染排放能值之和与可更新资源能值投入总量的比值。
ELR= (LN+FN+P)/(LR+FR)
(7)
(4)能值可持续指数(ESI)
ESI是一个广泛用于大范围长时期研究持续性的能值指标，该指标指示系统的可持续性。
ESI=EYR/ELR
(8)
1.6 数据处理及分析
所有数据采用SPSS 26.0统计分析软件进行One-way ANOVA方差分析，用Origin 8.5软件绘图。
2.1 小龙虾养殖排放尾水中TN、TP、COD和NH3-N浓度
夏季小龙虾基本收获后、水稻种植前的小龙虾养殖排放尾水中TN浓度范围为0.53～5.36 mg/L，平均值为2.03 mg/L，中位数为1.84 mg/L，最小值出现在沙市区，最大值出现在潜江市；
TP浓度范围为0.12～0.70 mg/L，平均值为0.24 mg/L，中位数为0.21 mg/L，最小值出现在洪湖市，最大值出现在沙市区；
COD浓度范围为6.60～78.39 mg/L，平均值为25.91 mg/L，中位数为23.35 mg/L，最小值出现在公安县，最大值出现在潜江市；
NH3-N浓度范围为0.34～1.75 mg/L，平均值为0.76 mg/L，中位数为0.62 mg/L，最小值出现在洪湖市，最大值出现在公安县。总体来说，江汉平原稻虾轮作模式小龙虾养殖排放的尾水中TN、TP和COD平均浓度均高于《地表水环境质量标准》(GB 3838-2002)[20]Ⅲ类水质标准，NH3-N 浓度高于Ⅱ类水质标准(NH3-N浓度≤0.5 mg/L)。按各采样点排放尾水浓度的达标率来说，依据《淡水池塘养殖水排放要求》(SC/T 9101-2007)[26]，TN一级达标率达到86.11%，二级达标率达97.22%，TP一级达标率达到97.22%，二级达标率达100%；
依据《地表水环境质量标准》(GB 3838-2002)Ⅲ类水质标准，TN达标率仅达到8.33%，TP达标率为50.00%，COD达标率达到38.89%，NH3-N达标率为80.56%(表2)。在采样的不同县(区、市)中，潜江市稻虾轮作模式小龙虾养殖排放尾水中TN、COD和NH3-N浓度均最高，TP浓度也仅次于沙市区，处于较高水平(图2)。
表2 小龙虾稻田养殖排放尾水中TN、TP、COD和NH3-N浓度Tab.2 Concentrations of TN、TP、COD and NH3-N in the discharge water from crayfish culture paddy field
图2 不同县(市、区)排放尾水中TN、TP、COD和NH3-N浓度(不同字母标注表示存在显著差异)Fig.2 The concentrations of TN、TP、COD and NH3-N in the discharge water from different counties (significant differences were marked with different letters above the column)
2.2 面源污染排放综合评价
稻虾轮作模式中，小龙虾养殖排放的尾水中各污染物质等标负荷和污染负荷比计算结果见表3。各污染物等标排放负荷量依次为TN>TP >COD>NH3-N，等标污染负荷比分别为38.08%、24.02%、22.69%和15.21%，不同县(市、区)稍有差别，其中除了沙市区的TP等标负荷最高，其他地区均为TN的等标负荷最高。这表明稻虾轮作模式小龙虾养殖排放尾水的主要污染物为TN，其次为TP和COD，TN是污染防治的关键污染物。
表3 各污染物等标污染负荷及污染负荷比Tab.3 Equal standard pollution load and load ratio
2.3 排污系数
稻虾轮作模式中，根据小龙虾养殖排放尾水中TN、TP、COD和NH3-N的初始进水浓度和排水浓度，再根据调研所得的稻虾轮作模式的平均产量增量1492 kg/hm2，平均排水量为4558 m3/hm2，根据排污系数计算公式，得到TN、TP、COD和NH3-N的单位产量排污系数分别为2.994、0.458、35.132和1.405 kg/t(表4)。可以看出，小龙虾养殖排放尾水中COD的排污系数最高，其次为TN和NH3-N，TP的排污系数最低。排污系数的得出可为进一步估算区域稻虾轮作模式面源污染排放量提供计算依据。稻虾轮作模式小龙虾养殖排放尾水中TN、TP、COD和NH3-N的单位面积排放系数分别为4.46、0.68、52.42和2.10 kg/hm2。由于本研究未对稻虾轮作模式水稻种植期间稻田的氮磷排放数据进行监测，参考陈玲等[27]的研究结果，稻虾轮作模式水稻种植过程中TN和TP排放系数与单季稻的流失量相当，将本研究得出的小龙虾养殖阶段尾水单位面积排放系数与单季稻的排放系数相加，估算得出稻虾轮作模式TN、TP和NH3-N的周年面源污染排放系数分别为16.31、1.19 和4.22 kg/hm2。与《全国农田面源污染排放系数手册》[28]中相同地理条件区域的单季稻模式、稻麦轮作模式和稻油轮作模式的单位面积氮磷排放系数进行对比可知，稻虾轮作模式的TN排放系数低于稻油轮作模式和稻麦轮作模式，高于单季稻模式，而TP和NH3-N的单位面积排放系数均高于其它3种模式(表5)。
表4 排污系数估算Tab.4 Pollutant discharge rate of the rice-crayfish rotation mode
表5 单位面积排放系数比较Tab.5 Comparison of pollutant discharge rate per hectare
2.4 基于面源污染的能值评估
2.4.1 能值投入与输出分析 江汉平原单季稻和稻虾轮作模式的能值分析结果见表6，各项能值投入比例在能值分析表中列出。两种模式投入的资源可分为两部分：本地自然资源和购买的外部资源，由氮磷面源污染排放产生的负能值单独列出。
表6 包含面源污染的能值分析*Tab.6 Emergy analysis with the integration of environmental pollution
由表6可以看出，稻虾轮作模式的能值总投入为2.15×1016sej/(hm2·a)，高于单季稻模式的9.97×1015sej/(hm2·a)，其主要原因是稻虾轮作模式增加了小龙虾养殖期间灌溉水、电力、饲料、虾苗等投入。从投入的资源结构来看，本地可更新资源年投入能值分别为1.75×1015和4.18×1015sej/(hm2·a)，分别占能值投入总和的17.51%和19.46%，其中灌溉河水能在可更新资源中所占比例最大，分别占总投入能值的13.04%和17.38%。由于稻虾轮作模式中小龙虾养殖消耗更多的灌溉水资源，使稻虾轮作模式的本地可更新资源能值高于单季稻模式。稻虾轮作模式的购买性外部资源能值为1.71×1016sej/(hm2·a)，而单季稻模式为8.02×1015sej/(hm2·a)，分别占总能值投入的80.40%和79.63%，增加了0.91×1016sej/(hm2·a)，表明两种模式均主要依靠外部购买资源，且稻虾轮作模式需要更多的外部资源维持系统。鉴于可获取数据的限制，主要计算了面源污染排放的氮和磷的污染能值，参考水产养殖业和种植业污染能值转换率[29]，计算得出两种模式氮磷面源污染物排放的负能值产出分别为3.74×1015和4.81×1015sej/(hm2·a)，且稻虾轮作模式的负能值产出更高。
2.4.2 能值评价指标分析 能值评价指标结果(表7)显示，包含面源污染负产出的条件下，江汉平原稻虾轮作模式能值分析所得的可更新资源投入占比为51.53%，高于单季稻模式的41.89%，表明稻虾轮作模式对可更新资源的利用效率更高。稻虾轮作模式的能值产出率也高于单季稻模式，表明稻虾轮作模式的生产效率更高。稻虾轮作模式的环境负载率低于单季稻模式，而能值可持续指数高于单季稻模式，表明稻虾轮作模式相较于单季稻模式对环境产生的压力更低，具有更高的可持续发展能力。
表7 系统能值评价指标Tab.7 Emergy indices for the two rice production modes
3.1 稻虾轮作模式小龙虾养殖尾水排放污染物质浓度
潜江市虾稻轮作模式小龙虾收获后排放的尾水中 TN、TP和COD平均浓度分别为 2.21、0.50和59.0 mg/L[7]。本研究中测定的排水中TN、TP、COD和NH3-N浓度平均值分别为2.03、0.24、25.91 和0.76 mg/L，浓度范围与该结果较一致，其中TN和TP的浓度值满足《淡水池塘养殖水排放要求》(SC/T 9101-2007)[26]一级排放标准，也低于报道的池塘养殖尾水的排放浓度，与河蟹养殖池塘尾水排放浓度相当[30]，符合对环境污染较小的客观现象。但排放水中TN、TP和COD平均浓度均高于《地表水环境质量标准》(GB 3838-2002)Ⅲ类水质标准[19]，而此时为满足稻田种稻的需求，这部分养虾尾水需集中排放，直接排入周围沟渠管网水系，可能会增加周边水体富营养化风险。另外，因其从稻田直接排放，也需要和稻田径流排水浓度进行比较，与湖北省农田地表径流排放水中 TN、TP和NH3-N的浓度(8.20、0.36、0.49 mg/L)[31]相比，小龙虾养殖排放的尾水中TN和TP浓度均更低，NH3-N的排放浓度略高，表明小龙虾养殖尾水排放的营养冲击负荷更小。
等标污染负荷比指标能够有效地分析不同污染物对受纳水体污染的相对重要性。本研究中采用等标污染负荷法分析结果表明，小龙虾养殖排放尾水的主要污染物为TN，其次为TP和COD，这与小龙虾养殖过程中产生的残饵、粪便、残留虾肥等在田面水中累积造成了氮磷浓度升高，以及水稻秸秆全量还田后有机质分解使田面水中COD浓度升高的生产实际相关。
3.2 稻虾轮作模式排污系数
如何对非点源污染进行定量计算，量化农业活动与污染负荷之间的关系，从根本上控制污染，一直是研究的热点与难点[32]。水产养殖排污系数测算的科学性和准确性将直接影响水产养殖污染负荷的结果和结论，从而影响养殖产业发展规划和产业政策的制订，对渔业产业的可持续发展和水域环境保护将产生重大影响[16]。水产养殖污染物排放总量的测算是开展总量控制工作的前提和重要参考，而单位产量排污强度是衡量水产养殖业污染物排放情况的最主要指标[33]。2017年全国单位水产品养殖产量的TN、TP、COD和NH3-N排污系数分别为2.02、0.33、13.6 和0.45 kg/t[34]。第二次全国污染源普查《排放源统计调查产排污核算方法和系数手册》[11]中湖北省水产养殖TN、TP、COD和NH3-N的排污系数分别为1.345、0.095、28.358 和0.441 kg/t。本研究中计算得出稻虾轮作模式小龙虾养殖尾水排放的排污系数均高于全国和湖北省，主要原因是稻虾轮作模式单位养殖面积小龙虾产量更低，导致单位产量的排污系数更高。
《第一次全国污染源普查水产养殖业污染源产排污系数手册》[10]中湖北省池塘养殖小龙虾的TN、TP和COD的排污系数分别为2.281、0.485和2.2135 kg/t，池塘养殖青虾的TN、TP和COD的排污系数分别为2.658、0.566和2.488 kg/t，池塘养殖河蟹的TN、TP和COD的排污系数分别为2.597、0.458和54.979 kg/t。太湖流域池塘养殖河蟹的TN和TP排污系数均值分别为2.49和0.701 kg/t[16]。表明本研究中稻虾轮作模式小龙虾养殖尾水TN和TP排污系数与池塘养殖虾、蟹的TN和TP排污系数相当。当某种水产品排污系数之和>100 g/kg 时，表明养殖该类水产品存在对周边河流、湖泊等外部水体环境的污染风险，应高度关注该类水产品[35]。本研究中稻虾轮作模式的TN、TP和COD的排污系数之和<50 kg/t，也表明稻田养殖小龙虾为低排污系数种类。通过更全面的稻虾轮作模式的排放污染监测，较精确地计算得出该模式的排污系数后，可以测算得出湖北省乃至江汉平原地区稻虾轮作模式小龙虾养殖尾水中TN、TP、COD和NH3-N的排污量，以进一步评估这些排放的营养物质对该区域湖泊富营养化程度的影响。目前，已估算和报道的虾蟹池塘的污染排放系数虽然较低，但也已纳入农业面源污染排放监督体系，因此，为全面估算农业面源污染排放量，稻虾轮作模式面源污染排放也应纳入农业面源污染源排放监测体系和监管当中。鉴于可获取数据的限制，由于稻虾轮作模式兼具种养特性，其面源污染既包含种植业污染，又包含养殖业污染，可能出现漏算或重复计算的问题，建议进行农业面源污染普查时，建立相应的计算方法和统计口径，准确得出相关数值。
佀国涵等[8]报道潜江地区虾稻模式夏季排水TN、TP和COD的单位面积排放系数分别为9.7、2.2 和258 kg/hm2。本研究中，稻虾轮作模式小龙虾养殖尾水排放的TN、TP和COD的单位面积排放系数均低于该报道，这可能是潜江市小龙虾稻田养殖年限更长，且大多数养殖户在种稻期间开展了虾稻共作，其更高的饵料投入量和多年的残饵累积可能导致更多的氮磷和有机质随排水排出，导致更高的排水浓度和排污强度。稻虾轮作模式周年的TN、TP和NH3-N排放系数均高于单季稻种植模式，其原因是与单季稻模式相比，稻虾轮作模式增加了小龙虾养殖周期和更多的投入品而使环境负荷增加。但与稻麦轮作和稻油轮作模式相比，稻虾轮作模式TN排放系数更低，TP和NH3-N排放系数更高，表明由稻麦轮作和稻油轮作等模式转换为稻虾轮作模式后，可能会造成TP和NH3-N排放量的增加，但TN排放量可能会降低。当前，稻虾田与耕地(水田和旱地)、水域、林地等利用类型之间产生较多转换，很多种植油菜、小麦等作物的旱地也改为虾稻田[2]。因此，对区域的农业污染负荷进行评估时，需精确区分稻麦轮作、稻油轮作、稻虾轮作等不同模式面积，以更加准确地估算区域农业面源污染总量。值得注意的是，稻虾轮作模式COD的排放系数最高，但其等标排放量却不是最高的，说明COD的排放量虽然很大，但是对于水环境恶化的影响不是最强的，结合等标排放负荷比和单位面积排放系数指标，稻虾轮作模式小龙虾养殖尾水中排放的TN最应受到重视，是农业面源污染防治的首要指标。
3.3 包含面源污染的模式能值评估
Hou等[14]运用能值分析方法对小龙虾稻田养殖模式的可持续性进行评价的结果表明，相比较于单季稻模式，虾稻模式可更新资源投入占比降低10%～14%，能值产出率降低9%，能值可持续指数降低23%～26%，环境负载率增加18%～23%。Chen等[4]的研究结果也表明单季稻模式的能值产出率最高，虾稻模式环境负载率最高，能值可持续指数最低。与稻麦轮作模式和单季稻等模式相比，稻虾模式面临更高的环境压力。另外，Hu等[36]应用生命周期法评价虾稻模式，结果表明虾稻模式比单季稻模式具有更高的环境富营养化潜力。不同地区的实际情况可能影响其可持续性，能解释得出相反结论的原因[37]，而且以上研究仅从资源利用角度计算的环境负载率可能不能全面阐述生态系统对环境的压力。本研究在包含面源污染的条件下，与单季稻模式相比，虽然稻虾轮作模式的氮磷污染排放系数更高，尾水排放造成环境污染负能值更高，且投入的外部资源中增加了饵料和苗种投入能值，但稻虾轮作模式中水资源消耗能值投入远远高于单季稻模式，更好地利用了水资源这一可更新资源，使稻虾轮作模式的环境负载率更低。值得注意的是当前灌溉水资源是可更新资源，目前还比较充裕，不需要购买，在水资源日趋紧张的形式下，水资源的可更新程度降低，进行能值计算时，稻虾轮作模式的环境负载率可能会更高。另外，尽管稻虾轮作模式的环境负载率低于单季稻模式，但两者的可持续指数均 <1，提示两种模式的可持续生产性均较差[38]，均需进一步优化，以提高模式的环境安全性。
3.4 存在的问题
本研究仅监测了小龙虾养殖向水稻种植转换期小龙虾养殖尾水的面源污染排放浓度，且仅监测了江汉平原6个县(市、区)的小龙虾养殖排放尾水，未覆盖整个江汉平原的小龙虾稻田养殖区域。另外，由于稻田排水呈间歇性，实际监测中难以采到正在排放的尾水，所以大多采集小龙虾养殖末期稻田的田面水样作为养殖尾水排水水样，且未监测小龙虾养殖过程中补水的水量和污染物质浓度，进行排污系数计算时仅减去初始补水污染物质的量，可能使计算的氮磷排放系数产生偏差。稻虾轮作模式中水稻种植期间未现场监测稻田氮磷排放数据，仅参考文献数据进行计算，覆盖全年的面源污染排放数据还需进一步监测分析。能值分析中，由于小龙虾养殖排放的氮磷面源污染的太阳能值转化率没有参考数值，本文采用鱼类养殖的太阳能值转化率，从而导致估算的模式能值可能有误差。稻油轮作和稻麦轮作两种种植模式的投入产出数据未进行调研，未进行能值评价比较。另外，本研究仅监测了一个稻虾轮作周期内小龙虾养殖阶段的面源污染排放情况，尚需进行多年连续监测，才能准确得出稻虾轮作模式的面源污染排强度。
1)采用现场监测方法得出江汉平原稻虾轮作模式由小龙虾养殖向水稻种植转换期，小龙虾养殖排放尾水中TN、TP、COD和NH3-N的浓度范围分别为0.53～5.36、0.12～0.70、6.60～78.39和0.34～1.75 mg/L，依据《淡水池塘养殖水排放要求》(SC/T 9101-2007)，TN一级达标率达到86.11%，二级达标率达到97.22%，TP一级达标率达到97.22%，二级达标率达100%，符合其对环境污染较小的客观现象。但依据《地表水环境质量标准》(GB 3838-2002)Ⅲ类水质标准，TN达标率仅达到8.33%，TP达标率为50.00%，COD达标率为38.89%，NH3-N达标率为80.56%， TN、TP和COD平均排放浓度高于《地表水环境质量标准》(GB 3838-2002)Ⅲ类水质标准，表明当前的稻虾轮作模式存在一定面源污染风险，其面源污染排放应受到重视，并纳入农业面源污染源排放监管体系和水产养殖业污染源普查当中。等标污染负荷法分析结果表明TN的等标负荷比最高，其次是COD和TP，TN是稻虾轮作模式面源污染控制的关键污染物。
2)稻虾轮作模式小龙虾养殖排放尾水中TN、TP、COD和NH3-N的单位面积排放系数分别为4.46、0.68、52.42和2.10 kg/hm2，单位产量排污系数分别为2.994、0.458、35.132和1.405 kg/t。
3)包含面源污染的能值分析结果表明，江汉平原稻虾轮作生态系统的可更新资源投入占比、能值产出率和能值可持续指数均高于单季稻模式，环境负载率也更低。表明稻虾轮作模式仍具有较强的可持续发展能力，但其能值可持续指数<1，还需进一步优化，提高模式的环境安全性，保障模式的稳步可持续发展。
4)研究结果为进一步开展稻虾轮作模式的面源污染负荷评估提供了基础数据和参考，建议继续开展尾水排放监测工作，以准确估算出稻虾轮作模式的排污系数和面源污染排放量，以期为农业面源污染防治和稻虾轮作模式可持续发展提供一定依据。
猜你喜欢能值稻虾面源基于国家粮食安全下的农业面源污染综合防治体系思考中国农业科学(2022年17期)2022-09-19农业面源污染的危害与治理商品与质量(2021年43期)2022-01-18澄江市农业面源污染成因及对策云南农业科技(2021年6期)2021-12-302019年第一季稻虾行情报告新鲜出炉，致敬奋战在一线的广大虾农们当代水产(2019年7期)2019-09-03安徽省农业生态经济系统能值分析*中国农业资源与区划(2019年12期)2019-04-14基于能值分析法的大庆石化企业生态效率研究经济研究导刊(2018年22期)2018-08-24发展稻虾产业，重振鱼米之乡当代水产(2018年12期)2018-05-16——以湖南新晃县（2006年～2015年）为例">生态经济系统的动态能值分析——以湖南新晃县（2006年～2015年）为例吉林农业(2017年21期)2017-11-02南县稻虾产业富民强县湖南农业(2017年1期)2017-03-20肉鸡日粮添加或不添加酶条件下不同木薯产品能值的评定广东饲料(2016年6期)2016-12-01