三峡库区汉丰湖流域水环境容量实证研究

何家军 樊连生 张慧

摘要 采用二维水质模型对三峡库区汉丰湖流域流量、水位、污染源排放、水动力和水质进行模拟，依据国家标准和规范计算了2017年的三峡库区汉丰湖流域水环境容量。结果表明，三峡库区汉丰湖流域COD、氨氮（NH3-N）、总氮（TN）、总磷（TP）的水环境容量分别为30 596.93、2 795.97、2 882.35、144.12 t/a。目前三峡库区汉丰湖流域水体的COD、TN、TP已经达到水环境容量限值，属于饱和运行状态;NH3-N污染负荷已基本接近维持现有水质功能下的水环境容量阈值。

关键词 水环境容量;污染负荷;二维模型;汉丰湖流域

中图分类号 X 26文献标识码 A文章编号 0517-6611（2020）01-0062-03

doi：10.3969/j.issn.0517-6611.2020.01.020

开放科学（资源服务）标识码（OSID）：

Empirical Study on Water Environmental Capacity of Hanfeng Lake Basin in Three Gorges Reservoir Region

HE Jia jun1，FAN Lian sheng2，ZHANG Hui3

（1. School of Management， Wuhan University of Engineering， Wuhan， Hubei 430205; 2. Yangtze River Engineering Supervision Consulting Co.， Ltd.， Wuhan， Hubei 430072; 3. School of Chemistry and Environmental Engineering， Wuhan University of Technology， Wuhan， Hubei 430205）

Abstract A two dimensional water quality model was used to simulate the discharge， water level， pollutant discharge， hydrodynamics and water quality of Hanfeng Lake Basin in the Three Gorges Reservoir Region. The water environmental capacity of Hanfeng Lake Basin in Three Gorges Reservoir Region in 2017 was calculated according to the national standards and standards. The results showed that the water environmental capacity of COD， NH3 N， TN and TP in the Hanfeng Lake Basin of Three Gorges Reservoir Region were 30 596.93 t/a， 2 795.97 t/a， 2 882.35 t/a and 144.12 t/a， respectively. At present， the COD， TN and TP of Hanfeng Lake Basin in Three Gorges Reservoir Region were close to the water environmental capacity， which belongs to the saturated operation state; the NH3 N pollution load had basically approached the water environmental capacity threshold under the existing water quality function.

Key words Water environmental capacity;Pollution load;Two dimensional model; Hanfeng Lake Basin

环境容量是环境目标管理的基本依据，是环境规划的主要环境约束条件[1]。研究水环境容量可以为水环境规划与管理提供技术支持，对于水域污染物排放的总量控制，保护和改善目前的水生态环境具有重要现实意义。水环境容量是指在满足水环境质量标准的前提下，水体最大允许污染负荷量，又称水体的纳污能力[2-3]，反映了污染物在环境中的迁移、转化和赋存规律，也反映了水环境在满足特定功能条件下对污染物的承受能力[4]，其容量大小与水体特征、水質目标及污染物特性有关。水环境容量计算方法主要以水动力学水质模型为主[5]，通常首先设定目标水质和一定保证率下的设计流量，再建立一维或二维的水环境容量模型，并在模型参数的率定基础上进行估算;也有研究通过随机规划、概率模型和未确知数学法等方法探讨水环境容量计算[6-7]。

三峡库区作为长江上游生态涵养区[8]，库区重要支流的水资源保护和水环境治理直接关系到整个长江流域水安全和生态安全[9]，对长江流域经济社会的可持续发展具有重要意义。汉丰湖是三峡库区内最大的库中湖，同时也是生态环境问题集中反映区域[10]。研究三峡库区汉丰湖流域水环境容量对长江流域社会经济发展具有重要意义。该研究根据2017年环境监测数据，采用二维水环境容量模型对三峡库区汉丰湖流域的COD、氨氮（NH3-N）、总磷（TP）、总氮（TN）的水环境容量进行计算，以期为该流域的水资源管理、水污染防治和总量控制规划提供科学依据。

1 资料与方法

1.1 研究区域概况

汉丰湖流域位于小江乌杨桥水位调节坝以上，包括东河、南河2条河流。东起乌杨桥水位调节坝，西南至开州区县境巫山镇;北到东河源头白泉乡。流域位于107°55′48″～108°53′00″E、30°52′31″～31°30′00″N，涉及云枫街道、汉丰街道、文峰街道、镇东街道、丰乐街道、白鹤街道、大德乡、临江镇、竹溪镇、镇安镇、九龙山镇、郭家镇、温泉镇、和谦镇、谭家镇、关面乡、白泉乡、河堰镇、大进镇、满月乡、铁桥镇、巫山镇、中和镇、三汇口乡、义和镇、南雅镇、麻柳乡、敦好镇、天和镇、高桥镇、紫水乡，共31个乡镇街道办。总面积2 534.70 km2、人口数97.93 万人。汉丰湖流域行政区划及水系见图1。

1.2 研究方法

1.2.1 因子选择。

根据汉丰湖流域水污染现状调查和污染物特征分析，选择影响汉丰湖水环境的最为敏感的指标因子COD、NH3-N、TN、TP。

1.2.2 模型选择。

污染物进入水体后，在一定范围内经过平流输移、纵向离散和横向混合后达到充分混合[11]，相对于污染物浓度仅在纵向上发生变化的一维模型，当水中污染物浓度仅在一个方向是均匀的，而在其余2个方向变化的条件下，一维水质模型不再适用，这时排入水中的污染物的分布用二维水质模型来计算[12]。

二维水质模型通常假定污染物浓度在水深方向是均匀的，而在纵向和横向上是变化的。二维水质模型的控制偏严，计算量大，但是较之零维、一维水质模型，计算结果更为精确。因此，在汉丰湖水环境容量计算时，考虑到汉丰湖的实际污染分布状况，同时为了保证计算的精确性，采用了二维水质模型。基本方程组为：

式中，Z为水位（m）;h为水深（m），h=Z-ZB;ZB为湖底高程（m）;u为x方向分速度（m/s）;v为y方向分速度（m/s）;q为湖面降雨、蒸发及湖底渗漏等水量源汇项（m/s）;f为柯氏加速度，f=2ωsinφ，φ为纬度，ω为地球自转速度（ad/s）;Sfx、Sfy分别为x方向和y方向的切应力，有床面阻力和水面风切应力两部分组成。

水质计算方程如下：

（hC）t+（hUC）x+（hVC）y=

x（hExCx）+y（hEyCy）+hS86 400+Sw（4）

式中，h为水深（m）;C为某种水质指标的浓度（mg/L）;U为x向沿垂向的平均流速（m/s）;V为y方向沿垂向的平均流速（m/s）;t为时间（s）;Ex为x方向扩散系数（m2/s）;Ey为y方向扩散系数（m2/s）;S为某种水质指标的生化反应项，g/（m3·d）;Sw为某种水质指标的外部源汇项（g/s）。

1.2.3 模型概化。

根据汉丰湖流域长系列水文统计资料和风况资料，确定汉丰湖流域设计水文条件和设计风况，利用二维风生湖流数学模型，模拟太湖设计流场分布，作为汉丰湖纳污能力计算的设计水动力条件[13]。采用模型概化湖泊型水（环境）功能区纳污能力计算采用湖（库）均匀混合模型，COD、NH3-N采用计算公式如下：

W=Q0×（Cs-C0）×0.086 4+K×V×Cs×10-6（5）

式中，W为水环境功能区纳污能力（t/d）;Q0为入湖（库）流量（m3/s）;Cs为水质目标浓度（mg/L）;C0为入湖（库）污染物浓度（mg/L）;K为污染物综合降解系数（d-1）;V为湖（库）容积（m3）。

对于较难降解的TN和TP，采用Dillon模型。依据物质平衡原理，考虑TN、TP负荷、滞留系数、湖水替换率、水库平均深度等因素在稳定条件下与水库中磷浓度关系，其计算方程可以概化为：

Ws=Cs×ρw×Z×A1-R（6）

式中，Ws为湖泊水体中氮、磷的环境容量（g/a）;Cs为湖泊水体中氮、磷的允许浓度（mg/L）;Z为湖泊水体平均水深（m）;ρw为水力冲涮系数（a-1）;A为湖泊面積（m2）;R为滞留系数。

1.2.4 参数率定。

根据《GB/T 25173—2010水域纳污能力计算规程》，通过资料收集情况及现场调查情况分析，采用汉丰湖入湖近30年月均流量（33.4 m3/s）为设计水文条件。

根据水功能区区划，确定各计算单元COD、NH3-N、TP、TN目标浓度值CS，根据地表水环境质量标准（GB 3838—2002）Ⅲ类水水质标准。C0根据汉丰2017—2018年湖乌杨大坝、巫山乡、津关、赵家大桥4个监测断面污染物监测的平均值确定。

降解系数，通过类比国内外部分河流降解系数的研究成果，以此来确定汉丰湖流域COD、NH3-N、TP和TN的降解系数。

2 结果与分析

2.1 水环境容量

根据水环境计算模型，该研究选择影响汉丰湖水环境的最为敏感的指标因子COD、NH3-N、TN、TP估算其环境容量。

计算结果显示，三峡库区汉丰湖流域COD、NH3-N、TN、TP的水环境容量分别为30 596.93、2 795.97、2 882.35、144.12 t/a。根据相关规划设计报告，三峡库区汉丰湖流域COD、TN、TP污染源入河量分别为41 743.15、3 710.16、344.48 t/a。模型计算结果表明，目前三峡库区汉丰湖流域水体的COD、TN、TP已达到水环境容量限值，属于饱和运行状态。NH3-N的水环境容量为2 795.97 t/a，除去调查已知的污染源入河量1 880.58 t/a，目前最大允许排放量仅为915.39 t/a。

但就NH3- N的水环境容量而言，目前三峡库区汉丰湖流域NH3-N污染负荷已接近维持现有水质功能下的水环境容量阈值，如不能及时对NH3-N污染排放源进行控制和削减，则三峡库区汉丰湖流域水体也极易发生NH3-N水环境容量超负荷现象。