基于质量平衡的太湖氮、磷自净能力计算*

翟淑华，韩 涛，陈 方

(1:水利部太湖流域管理局，上海 200434)

(2:中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室，南京 210008)

太湖是我国第三大淡水湖泊，具有供水、防洪、抗旱、旅游、养殖、航运等多种功能，1990s 以来，伴随着太湖流域经济的高速发展，入湖河道水质恶化，加之污染物直接入湖量增加，导致太湖水质下降，富营养化问题突出.2007年太湖蓝藻暴发引发无锡供水危机后，太湖水环境综合治理的力度进一步加大，对太湖问题的研究成为国家层面、国际间合作的热点和焦点.对太湖水体自净能力以及通过水动力条件的改变对太湖自净能力的作用进行研究，对于富营养化湖泊治理具有重要的理论指导意义，引起了国内外不同层面的管理人员和专家、学者的关心、关注.国外对湖泊营养盐自净能力问题的研究较早且较深入[1]，美国环境保护局于1972年基于水体自净能力最先提出了实际环境能够承担的最大日负荷总量TMDL(total maximum daily loads，TMDL)的概念.污染排放量超过该负荷，水环境就可能恶化，进而对水域使用功能和生态系统产生损伤.Vagnetti 等[2]较为全面地从物理、化学和生物过程等方面，通过现场实测资料分析和历史资料的对比分析，对不同底质下河流对氮、磷、重金属等多种污染物的自净能力进行了定量分析，得出该河流对于不同污染物的自净量和自净能力以及主要作用因子等研究结论.对于太湖水体自净能力的研究，在“国家水体污染控制与治理科技重大专项”中曾列有专题[3]，以太湖南部水域为对象开展了探索性研究，研究南太湖水体自净能力季节变化的特征，得出影响南太湖地区小型浅水湖泊自净能力的3 个主成分分别为浮游植物因子(水温、pH、高锰酸盐指数和叶绿素a)、农田排水因子(pH、溶解氧和总氮)和营养因子(总氮和总磷).

国内专家也曾结合水环境容量研究进行过探索[4-5].影响太湖自净能力的因素多样且十分复杂[6]，从物理角度分析，水流流速、流向的不同将直接对污染物迁移、扩散方向和强度带来影响，从生物化学角度分析，水体本身的组分决定了生物和化学进程对水域自净能力作用的大小，如湖泊中挺水生物的多少和水质本底浓度的高低直接影响着湖体自净能力的大小.此外，风对水体的强烈扰动又会增加水体中的溶解氧，对增加水域的自净能力有益，在水域床底长期积累的底质污染，即内源污染积累对水域的自净能力也有不可忽视的间接影响.自净能力包括了稀释、吸附、挥发、沉淀等物理自净以及氧化还原、絮凝等化学自净和生物净化等能力.它不是一个固定值，而是随年度水情、水流、水质以及水生生物条件等自然环境和人类活动而变化.随着最严格水资源管理制度的实行，严格控制入河、湖排污总量已成为国家水资源管理的要求，这对水域净化能力的研究提出了更高的要求，有必要进一步深入研究太湖水体自净能力与净化能力，为从严核定水域纳污容量提供理论支撑.为此，本文选取资料较为全面的2010年，开展进出太湖的水量平衡计算和营养盐的质量平衡计算，通过对2010年太湖进出湖水量和氮、磷通量计算，以及在进出太湖质量平衡相关要素的资料收集和计算分析的基础上，对太湖氮、磷自净能力进行探索性量化研究.

1 材料与方法

1.1 水量平衡要素及计算

1.1.1 环太湖河流进出湖水量计算 太湖地处平原水网区，水系复杂，河道纵横交错，出入太湖溇港河道共计228 条，其中河道口门敞开的河道有62 条.为弄清太湖河道出入水量以及水质和掌握太湖水情及污染状况，国内学者展开了大量的进出太湖污染负荷通量的研究[8-19]，太湖流域水文部门在主要出入湖河道设置了多个控制断面，实施环太湖水文巡测以及资料整编，是环湖河道通量计算的主要数据基础.

2010年环湖出入湖河道水量数据，采用环太湖国家水文基准站经过整编后的“2010年环太湖水文巡测资料”.由于巡测线的实际工况不断变化，2010年湖州、苏州、无锡三市水文部门对环太湖水文巡测段进行了新的调整:由原来的“11 段12 站共计131 个进出水口门”调整为“10 段13 站共计130 个进出水口门”，其中湖州有2 段3 站共计21 个进出水口门;巡测线移至太湖边，撤销三里桥段，设幻溇段;新增城北水闸站.苏州有5 段2 站共60 个进出水口门;巡测线移至环太湖公路，撤销五福桥段和新通安桥段，设胥江大桥段和铜坑闸段.无锡有3 段8 站共49 个进出水口门;撤销漕桥、黄埝桥基点站，设浯溪桥基点站代替，漕桥+黄埝桥段改名为浯溪桥段;浯溪桥段新增小金桥、葛渎桥、师渎桥3 个口门;陈东港桥段新增庙渎桥、双桥、定跨桥3 个口门;另外有大港桥、雅浦桥、龚巷桥、湖山桥、大渲河泵站、犊山闸、梅梁湖泵站、五里湖闸等8 个单站.

2010年环湖出入湖河道水量计算依据调整后的环太湖水文巡测线，即10 段13 基点站共计130 个进出水口门(图1)的逐日进出湖流量.2010年各基点站与巡测单站水量测验频次为每天1 ～2 次，在水情存在明显变化时进行加密监测.根据单站每日水量测验实测资料以及单站与所代表巡测段总流量的相关关系，推求出巡测段逐日进出太湖的水量，据此计算得到环湖出入湖河道逐月进出湖水量和年进出湖水量.2010年环太湖河流入湖水量为118.81×108m3，其中包括望虞河引江济太的水量，出湖水量为110.06×108m3.

1.1.2 降雨量、蒸发量计算 2010年太湖湖面降水深为1102.1 mm[7]，按全年平均水位3.37 m，对应的太湖水面面积为2332.75 km2，折算为降雨量为25.71×108m3;湖面蒸发量为884.23 mm(湖面蒸发量采用杭长桥、西山、瓜泾口、枫桥4 站平均值)，按全年平均水位3.37 m，对应的太湖水面面积为2332.75 km2，折算为蒸发量为20.63×108m3.

1.1.3 太湖调蓄水量计算 根据“2010年度太湖流域及东南诸河水资源公报”，2010年年初太湖水位为3.26 m，蓄水量为50.83×108m3;年末水位为3.10 m，蓄水量为47.17×108m3;太湖蓄变量(蓄变量为年末蓄水量减年初蓄水量)为-3.66×108m3.

1.1.4 其他水量要素计算 根据“2010年度太湖流域及东南诸河水资源公报”，直接从太湖取水水量(包括自备水源、自来水厂等取水户)约为15.0×108m3.根据国务院批复的《太湖流域水资源综合规划》数据，太湖巡测线内陆地产水模数为28.18×104m3/km2，环太湖水文巡测线所围面积为3454 km2，扣除农业灌溉耗水量3.40×108m3，估算得到巡测线内陆地产水入湖量约为0.82×108m3.

图1 2010年环太湖进出湖河流水量巡测线示意图Fig.1 Sketch map of hydrometric gauging line of rivers flowing in and out of Lake Taihu in 2010

1.1.5 水量平衡计算 考虑上述各水量平衡要素，开展对各太湖水量平衡项进行统计计算，其计算公式为:

式中，WI为环太湖河流入湖水量，根据“2010年环太湖出水文巡测资料”逐日统计得到;W0为环太湖河流出湖水量，根据“2010年环太湖出水文巡测资料”逐日统计得到;P 为湖面降水量，为太湖年降水量与2010年太湖平均水位对应的水面面积的乘积;E 为湖面蒸发量，为太湖年蒸发量与2010年太湖平均水位对应的水面面积的乘积;D 为陆地产水入湖量，为环太湖水文巡测线范围内的陆地区域的产水入湖量;WC为取水户直接取水量，直接采用当年统计资料;ΔV 为太湖蓄变量，为太湖年初和年末蓄水量的变化量，当年初蓄水量大于年末蓄水量时ΔV 取负值，当年初蓄水量小于年末蓄水量时ΔV 取正值;ΔE1为水量平衡计算的绝对误差.

水量平衡计算结果显示:2010年太湖水量平衡入湖项 WI、P、D 分别为 118.81×108、25.71×108、0.82×108m3，出湖项 W0、E、WC分别为110.06×108、20.63×108、15×108m3，蓄变量 ΔV 为-3.66×108m3，水量平衡计算的绝对误差为-4.01×108m3，相对误差为-2.8%，在规范要求的误差范围以内.

1.2 水质资料及通量计算

1.2.1 环太湖河流水质与出入湖通量计算 水质监测资料采用太湖流域水资源保护局发布的太湖和环太湖主要进出湖河流每月1 次的常规水质监测资料，相关水质监测的断面布置、监测方法、监测频次以及环太湖河道出入湖氮、磷负荷通量计算采用翟淑华等[8]的方法，即采用与水量巡测段相匹配的河流水质实测资料，计算得到2010年环太湖河流入湖和出湖水量分别为118.81×108和110.06×108m3，带入太湖的总磷和总氮负荷通量分别为0.28×104和5.64×104t，带出太湖的总磷和总氮负荷通量分别为0.068×104和2.54×104t.

1.2.2 降雨带入与人工带出量估算 降雨带入量为太湖湖面年降水量与其浓度之积计算得到，其中总磷和总氮浓度采用黄漪平[9]的值.计算得到太湖2010年降雨带入太湖的总氮和总磷量分别约为2160 和51 t.人工带出太湖的总磷和总氮量主要通过取水、蓝藻打捞等生物收获资料估算，其中水厂直接取水带出的污染负荷量按照当年的取水量和原水水质浓度估算，带出的总磷和总氮量分别为76 和3075 t;2010年通过蓝藻打捞带出太湖的总磷和总氮量分别为67 和335 t[7]，通过水生生物带出的总磷和总氮量参考王鸿勇等[10]的值，约为116 和3028 t;2010年底泥疏浚量为418×104方[7]，依据江苏省水利厅底泥疏浚资料，2008-2011年底泥清淤工程年均疏浚562×104t，带出总磷和总氮分别为725 和2775 t，相应推求2010年底泥疏浚带出的总磷和总氮负荷量分别约为539 和2064 t;太湖地区污染源排放的负荷量大部分主要通过河道带入太湖，2010年环太湖巡测线已基本调整到太湖岸边，未调整的部分河道由于“引江济太”以出湖为主，因此在太湖污染负荷量平衡计算中，暂不考虑巡测线内污染源直接入湖量.

1.2.3 太湖年初原有量和年末残留量计算 太湖年初原有量和年末残留量计算，分别考虑太湖蓄水量、分区水面积和分区水质，其中蓄水量和分区水面积按照年初和年末水位条件下水位—库容曲线和水位面积差分计算.为提高湖区水质监测点的代表性，采用泰森多边形法对太湖进行分区，按照水质监测站点将太湖划分为33 个分区(图2).

图2 太湖湖区划分及水质代表站分布Fig.2 Sub zones and water quality monitoring sites in Lake Taihu

采用面积加权计算得到太湖污染负荷年初总磷和总氮原有量分别为477.802 和11487.58 t，年末残留量分别为297.171 和5566.06 t.

综上所述，得到2010年进出太湖氮、磷负荷通量计算结果，详见表1.

2 结果与讨论

本文利用实测资料，在上述水量平衡、氮磷负荷通量计算的基础上，以物质平衡原理为基础，在宏观层面将太湖水体对污染物的自净过程看作一个黑箱，不论其通过物理自净还是通过化学和生物自净，认为进入太湖的河流水质经过太湖水体的自净过程后，其自净量值即为在考虑水量平衡误差的条件下的入湖项扣除出湖项与水体本身含有的总量后的值.自净量与入湖量的比值即为自净能力.以2010年总磷、总氮负荷通量计算为例，对太湖总磷和总氮的自净量和自净能力计算结果进行分析及讨论.

表1 2010年进出太湖总磷和总氮负荷通量Tab.1 Total phosphorus and total nitrogen flux flows in and out of Lake Taihu in 2010

2.1 基于进出湖负荷通量的太湖自净能力

采用考虑出入湖负荷量的方法进行自净能力的计算，在不考虑太湖水体本身年初和年末调蓄量变化的条件下，太湖氮、磷自净量的计算公式为:

式中，W自净为太湖水体自净作用去除的污染物质量，包括污染物降解、生物转化等;W入为通过河流直接带入太湖的负荷量;W降雨为通过降雨直接带入太湖的负荷量;W出为通过河流直接带出太湖的负荷量;W人工为通过水厂取水、蓝藻打捞和太湖底泥疏浚等途径带出太湖的负荷量;ΔW 为质量平衡计算的绝对误差.

相应地，太湖自净能力计算公式[3]为:

表2 仅考虑出入湖负荷量计算的太湖总磷和总氮的自净能力Tab.2 Self-purification capacity of total phosphorus and total nitrogen with the nutrient flows in and out of Lake Taihu

式中，η 为自净能力;Qi为入湖污染物量，包括入湖河流、降雨等自然输入量等，即Qi=W入+W降雨;Qo包括出湖河流、直接取水带出量、蓝藻打捞、底泥疏浚等，即Qo=W出+W人工.

按照上述计算公式，在不考虑误差的条件下，即认为 ΔW=0 的条件下，计算2010年太湖总磷和总氮的自净率分别约为48%和42%(表2).

2.2 基于质量平衡的自净能力计算

由于湖泊调蓄水量和水体水质处于变化之中，因此，湖体本身的氮、磷负荷量的变化在质量平衡中应该予以考虑，即在自净能力计算中考虑太湖本身水情和水质变化而产生的氮、磷年初原有量和年末残留量[9]的不同.采用质量平衡方程，计算太湖自净能力公式为:

式中，W年初为年初太湖水体中的原有量;W年末为年末太湖水体中的残留量.该方法同陈小锋等[11]采用质量平衡的方法原理一致，但增加了平衡误差项，2010年进出太湖的水量平衡误差为-2.8%.

相应地，太湖自净能力计算公式为:

按照上述计算公式，得到2010年太湖总磷和总氮的自净率分别约为52%和49%(表3).

表3 采用质量平衡法计算的太湖总磷和总氮的自净能力Tab.3 Self-purification capacity of total phosphorus and total nitrogen in Lake Taihu on the basis of mass balance

3 结论

进出太湖的河道水质对太湖水质的影响密切相关.本文利用已有的实测资料和相关研究数据，以质量平衡理论为基础，开展了太湖氮、磷自净能力计算，初步得出2010年太湖总磷、总氮的自净能力分别约为48% ～52%和42% ～49%.质量平衡方法建立在水量平衡和质量平衡的基础上，理论上较为严谨，也便于计算，可作为宏观上把握太湖自净能力的方法.但由于环太湖水文巡测线不是同步监测的，在计算水量平衡时可能会带来误差，特别是巡测期间有降雨或人工引调水存在时影响会更大.另外将每月1 ～2 次的资料作为月平均流量也会带来误差.为使计算结果更为准确，需要加强进出太湖水量和营养盐监测的频度和力度，此外要考虑对当年资料获取的完整性.

本文提出的太湖自净能力研究方法偏于宏观，尤其是太湖底泥带出量的计算受制于实际资料的获取，且难以得到底泥沉降和释放过程的实际资料，是影响自净能力估算结果精度的主要因素，故研究太湖自净能力作用机制更为科学的方法是采用太湖生态模型计算方法，即利用现有的资料和对太湖水质变化的追踪分析，结合实际观测资料，采用太湖生态模型对水质变化参数进行计算，分析铵态氮、硝态氮、亚硝态氮、底泥可交换态氮、正磷酸盐磷、底泥可交换态磷、底泥间隙水溶解性磷等不同形态营养盐转化关系，通过太湖水动力—生态模型探索物理、化学、生物过程对营养盐的净化作用，定量评估在不同污染负荷压力条件下，太湖水体不同的过程对总磷与总氮的净化能力，这将是太湖自净能力下一步研究的工作方向.

[1]Benoit RJ.Self-purification in natural waters.In:Ciaccio LL ed.Water and water pollution handbook，vol.1.New York:Dekker，1971.

[2]Vagnetti R，Miana P，Fabris M et al.Self-purification ability of a resurgence stream.Chemosphere，2003，52(10):1781-1795.

[3]许 磊，李 华，陈英旭等.南太湖地区小型浅水湖泊自净能力季节变化研究.环境科学，2010，31(4):924-930.

[4]范丽丽，沙海飞，逄 勇.太湖湖体水环境容量计算.湖泊科学，2012，24(5):693-697.

[5]胡维平.平原水网地区湖泊的水环境容量及允许负荷量.海洋湖沼通报，1992，(1):37-40.

[6]何少苓，彭 静.论提高水域纳污能力与自净能力的水动力潜力.中国水利，2003，(1):71-73.

[7]水利部太湖流域水资源管理局.2010年太湖健康状况白皮书，2010.

[8]翟淑华，张红举.环太湖河流进出湖水量及污染负荷(2000 2002年).湖泊科学，2006，18(3):225-230.

[9]黄漪平.太湖水环境及其污染控制.北京:科学出版社，2001.

[10]王鸿勇，张海泉，朱 喜等.太湖蓝藻治理创新与实践.北京:中国水利水电出版社，2012.

[11]陈小锋，揣小明，曾 巾等.太湖氮素出入湖通量与自净能力研究.环境科学，2012，33(7):2309-2314.

[12]燕姝雯，余 辉，张璐璐等.2009年环太湖入出湖河流水量及污染负荷通量.湖泊科学，2011，23(6):855-862.

[13]郑 一，王学军，江耀慈.环太湖河道水质分析与入湖污染物负荷量估算.地理学与国土研究，2001，17(1):40-44.

[14]沈国华，陈 颖.太湖进出水量变化分析.江苏水利，2008，(2):27-28.

[15]戈礼宾，张泉荣，黄利亚.太湖进出水量变化及水量平衡初步分析.江苏水利，2007，(11):33-34.

[16]罗 缙，逢 勇，林 颖等.太湖流域主要人湖河道污染物通量研究.河海大学学报:自然科学版，2005，33(2):131-135.

[17]马 倩，刘俊杰，高明远.江苏省入太湖污染量分析(1998 2007年).湖泊科学，2010，22(1):29-34.

[18]余 辉，燕妹雯，徐 军.太湖出人湖河流水质多元统计分析.长江流域资源与环境，2010，19(6):696-702.

[19]许 梅，任瑞丽，刘茂松.太湖人湖河流水质指标的年变化规律.南京林业人学学报:自然科学版，2007，31(6):121-124.