水环境容量计算方法研究及应用赵君(河海大学，江苏 南京 210098)E-mail:zsmzyq@摘要：一维稳态条件下计算水环境容量的3种方法，即段首控制方法、段尾控制方法和功能区段尾控制方法。
本文通过分析比较各方法的优劣及其相互联系，针对曹娥江支流--长乐河的具体情况,采用段首控制对其水环境容量进行计算，系统地将各方法的物理含义及其适用奈件推广到实际中。
计算结果证明了方法的可靠性。
关键词:水环境容量；段首控制；段尾控制；功能区段末控制1 计算方法1．1基本概念和方程水环境容量是在给定水域范围和水文条件，规定排污方式和水质目标的前提下，单位时间内该水域最大允许纳污量，称作水环境容量。
水环境容量具有资源性、区域性、系统性、发展需要性四个基本特征,其大小主要与水域特性、环境功能要求、污染物质以及排污方式有关，这些因素直接影响入流污染物的稀释能力以及污染物质在水体中的时空分布。
由于河流具有对污染物质的稀释、输移、降解能力，因此河流环境容量可分为以下三个组成部分：输移容量：污染物在水体中随水流的对流运动产生的输移量,它只与水力要素和水质目标有关，因此输移容量是有限的不可再生的。
较大的输移容量并不代表较大的允许排放量。
对保守物质来说，河段总的环境容量只由输移容量组成。
稀释容量：当水体本底水质浓度低于水质标准时，由于对流及扩散作用，使排入的污染物逐步均匀分布到整个水体，其浓度达到标准浓度的限值时，水体所增加的污染物容量。
稀释容量在数量上等于标准浓度时的输移容量与本底浓度时输移容量的差值，也称差值容量。
自净容量：由于沉降、生化、吸附等物理、化学和生物作用，给定水域达到水质目标所能自净的污染物量称为自净容量。
自净容量是反映水体对污染物的自净能力，也称同化容量。
自净容量是水环境容量中最重要的组成部分，河流水环境容量的计算关键在于自净容量的计算。
它是可不断再生的量。
河流是我国最常见、最基本的纳污水域。
河流的水环境容量占在我国的很大的比重。
污染物进入河流后，在一定范围内经过平流输移、纵向离散和横向混合后达到充分混合，或者根据水质管理的精度要求，允许不考虑混合过程而假定在排污口断面瞬时完成均匀混合，可按一维问题概化计算条件，建立水质模型。
考虑到国内，尤其是长乐河的具体治理能力及有关水环境保护要求[2],同时考虑到稳态模型发展较成熟，本文采用一维稳态模型进行容量计算。
其次，当前国内河流的污染都以有机物污染为主，曹娥江支流--长乐河的有机物污染相当严重，所以本文以COD,NH3-N为代表对有机物的环境容量计算方法进行研究。
第三，为了简化计算，本文只对点源污染（即排污口）进行计算，同时假定各排污口连续，均匀排污。
一维稳态情况下有机物的降解过程可用托马斯（Thomas）模型表述（忽略弥散）[3]:13(CUk k )C x∂=−+∂ (1) 在C(0)=C 0的初值条件下得到托马斯模型的积分：130()()()exp[]k k C x f x C x u+==−（2） 式中，C(x )、C 0为x ＝x、x＝0处河水各水质指标浓度（mg/L）；x为到排污口（x ＝0）的河水流动距离（m）；u为河水平均流速（m/s）；k 1，k 3分别是各水质指标衰减系数和沉浮系数（d －1）。
1．2段首控制段首控制中的段是指沿河任何两个排污口断面之间的河段，而段首则是指各段的上游第一个排污口断面。
某一功能区内各段的划分情况如图1所示。
段首控制就是控制上游断面（也即段首）的水质达到功能区段的要求，那么由于有机物的降解，则在该段内的水质处处达到或高于功能区段的控制指标。
段首控制严格控制了功能区段的水质不超标。
图1 功能区段内各小段的划分下面就段首控制的水环境容量计算方法进行分析。
在功能区段的段首，由于来水的COD,NH3-N 浓度与功能区段水质要求的差别，为来水提供了稀释容量：000(s E Q C C )=− （3）式中，E 0为功能区段段首的稀释容量（t/d）；C s 为功能区段水质标准（mg/L）； Q 0、C 0分别代表来水的流量（m 3/s）及浓度（mg/L）。
功能区段内任意一段的容量计算如图2所示。
由于控制各段段首为水质标准，那么经过一段降解（如图2中A、B 所示）后，到达段末时的降解量即为该断面处的环境容量（如图中C、D所示）。
A.第i-1段浓度衰减曲线；B.第i 段浓度衰减曲线；C.第i-1断面的环境容量；D.第i 断面的环境容量图2 段首控制容量计算示意图第i 个断面处的环境容量为：（4）1()(i i i s i s i i E Q q C Q C f x x −=+−−)1)]]式中，E i 为第i个断面处的环境容量（t/d）；q i 为第i个断面处的排污流量（m 3/s）；Q i 为混合后干流流量；其余符号意义同上。
则功能区段内所具有的总环境容量为：00011()[()(nni s i i s i s i i i i E E E Q C C Q q C Q C f x x −===+=−++−−∑∑ （5）化简式（5）得：0011()[(1())ns s i i i i i E Q C C C Q f x x q −==−+−−+∑（6） 1.3段尾控制段尾控制中段的划分与段首控制相似，只是段的控制断面在下游的排污口断面，亦即段尾。
段尾控制的目的在于让水质在各段末达到功能区段水质标准，那么我们可以反推出该段段首处的环境容量。
在段尾控制水环境容量计算中，功能区段全段水质低于水质要求，但考虑到降解能力很低，而且各小段较短，超标不会太高，因此水质超标很小。
功能区段内小段的划分如图1所示。
下面就段尾控制的水环境容量计算进行分析。
类似于段首控制，功能区段段首的稀释容量为：000()S E Q C C =−功能区段内任意一段的容量计算如图3所示。
由于控制断面在段末，那么可由段末按降解曲线（如图3中A 、B 所示）发推到段首，即可求得段首处的环境容量，如图3中的C 、D 所示。
11111()()i s i i i i i s E C Q q Q C f x x ++++=+−−（7） 式中，Q i+1为第i+1个断面前的干流流量（m 3/s）；q i+1是第i+1个断面处的排污流量（m 3/s）；其余各符号意义同上。
A.第i-1段浓度衰减曲线；B.第i 段浓度衰减曲线；C.第i-1断面的环境容量；D.第i 断面的环境容量图3 段尾控制容量计算示意图功能区段内所具有的总环境容量为：100011111()[()]()n n i s s i i i i i i i E E E Q C C C Q q Q C f x x −+++==+=+=−++−−∑∑1s （8）化简式（8）得：1001011()[(1)()n s s i i i i E Q C C C Q q f x x −1]i ++=+=−+−+−∑ （9）1.4功能区段段末控制功能区段段末控制就是在功能区段的最末断面控制水质。
这里的段末常常是特定功能区段的段首。
方法的实质是控制功能区段最终断面，而不是考虑段内水质变化是否超标。
类似于段首控制，功能区段段首的稀释容量为：000()S E Q C C =−功能区段内各排污口的浓度变化如图4所示，具体公式表述如下：1(()i i i i i i i i iQ C f x x q c c Q q −)−+=+ （10）同理可求得功能区段最终控制断面的浓度为：1(()n n n n n n n n nQ C f x x q c c Q q −)−+=+ （11）式中，q i 、c i 为断面处排污口的排污流量（m 3/s）及浓度（mg/L）；Q n 、 C n 为最终控制断面的干流流量（m 3/s）及浓度（mg/L）；q n 、 c n 为最终控制断面的排污流量（m 3/s）及浓度（mg/L）。
A.第i-1段浓度衰减曲线；B.第i 段浓度衰减曲线；图4 功能区段段末控制容量计算示意图本论文第二部分提供的计算结果是按排污总量比例消减后得到的数据。
所谓的按总量比例消减是指如果最终断面水质低于控制标准（c n >c s ）,则按各排污口的排污总量进行消减;反之,则相应让各排污口增排放量(c n <c s ),以计算环境容量。
若c n =c s ,则各排污口排放总量即为该段在现行情况下的环境容量。
经校核,消减后的水质在最终控制断面可以达标。
2 计算实例分析2.1水环境容量计算方法分析段首控制方法,因其控制非常严格,所以适用于对水质要求高、经济发达污染治理能力强的地区,或是水质较好的源头地区,该方法的适用对象应为污染较轻或者旨在改善水质条件的地区。
段尾控制方法对水质的要求略低于标准,所以适用于污染较严重、经济不太发达的地区,使用该方法的目的在于控制污染,而非改善水质。
功能区段段末控制方法适用于对特定区段水质有很高要求而对其它区段不限制的地区,在特定段段首处严格限制水质,其它段任意排放,只要满足特定段段首要求即可,特定段内不许排污。
该方法面向污染十分严重而又急需用水的地区.从长远角度看,该方法不可取,但在我国目前情况下还是很实用的。
2.2实例分析基于理论方法对不同地区适用性的考虑,本文采用段首控制的方法对长乐江进行水环境容量计算。
长乐江是曹娥江的支流之一,长度32.79km, 水质目标Ⅲ。
本次计算,水质监测资料采用2004年1月至12月长乐江的实地监测资料,范围选取的是22、23、24、18号断面。
COD Mn 与COD Cr 的转换系数, 即C(COD Cr )=KC(COD Mn ),通常同一水体中,K的取值为3.0左右,本次算取3.2。
另外河流流量、流速及水质浓度等在一年中的不同水期中变化很大,导致了河流的环境容量在一年中变化很大,所以本文以月为时段来计算每个月的环境容量。
对于段首控制:lim [1()](1)lk uz s s n E nQC f x nQC e−→∞=−=− (12)利用倒代换,令t=1/n,而后再利用罗必达法则求极限得: Z s lE k QC u= (13)表1 地表水环境质量标准 mg/L等级 ρ(DO) ρ(COD) ρ(BOD 5)ρ(NH3-N)I ≥7.5 ≤2 ≤3 ≤0.15 Ⅱ ≥6 ≤4 ≤3 ≤0.50 Ⅲ ≥5 ≤6 ≤4 ≤1.00 Ⅳ ≥3 ≤10 ≤6 ≤1.50 Ⅴ≥2 ≤15 ≤10 ≤2.00表2 长乐江COD cr 水环境容量(t/月)月份 12 3 456 总计 测算 结果 56.12171.53 110.68 69.23-24.4986.15 月份 7 89 10 11 12 测算 结果 182.45 24.46 577.58 146.72 41.1729.232370.92表3 长乐江NH3-N 水环境容量(t/月)月份 1 2 3 4 5 6 总计 测算 结果 22.14 -29.27 39.10 53.34176.07 79.07 月份 789 10 11 12 测算 结果29.2445.17933.53 38.231.4612.531400.61段尾控制与段首控制方法类似,其功能区段段首的稀释容量是相等的，区别是:段首控制的计算使全段水质高于标准，而段尾控制使全段略低于标准,故只对自净容量E Z 进行比较即可。