華祖林,程浩淼

(1. 河海大學(xué)淺水湖泊綜合治理與資源開(kāi)發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 南京　210098；2. 河海大學(xué)環(huán)境學(xué)院，江蘇 南京　 210098；3. 河海大學(xué)水資源高效利用與工程安全國(guó)家工程研究中心，江蘇 南京　210098)

?

一種考慮容量沿程變化的水環(huán)境容量計(jì)算方法

華祖林1,2,3,程浩淼1,2

(1. 河海大學(xué)淺水湖泊綜合治理與資源開(kāi)發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，江蘇 南京210098；2. 河海大學(xué)環(huán)境學(xué)院，江蘇 南京 210098；3. 河海大學(xué)水資源高效利用與工程安全國(guó)家工程研究中心，江蘇 南京210098)

摘要：基于水環(huán)境容量的定義和本質(zhì)內(nèi)涵，厘清目前常用的計(jì)算中小型河流的零維和一維環(huán)境容量計(jì)算方法中存在的不足。針對(duì)現(xiàn)有計(jì)算方法的不足與弊端，考慮了河流環(huán)境容量是沿程變化的這一事實(shí)，提出一種中小型河流環(huán)境容量的計(jì)算新方法，該方法消除了現(xiàn)行方法中夸大的部分環(huán)境容量，更契合天然實(shí)用情形；在此基礎(chǔ)上，推求出河道存在環(huán)境容量的臨界判據(jù)，即最大允許排污濃度或河道最大允許超標(biāo)長(zhǎng)度比例，并進(jìn)行了算例的應(yīng)用與比較。

關(guān)鍵詞：水環(huán)境容量；沿程變化；計(jì)算方法；臨界判據(jù)

隨著當(dāng)代社會(huì)經(jīng)濟(jì)的高速發(fā)展，工業(yè)生活廢污水大量排入河流，河流的水環(huán)境質(zhì)量不斷惡化，給我國(guó)社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展帶來(lái)嚴(yán)峻的挑戰(zhàn)。因此各級(jí)政府在環(huán)境管理中，準(zhǔn)確確定河流的水環(huán)境容量具有重要的作用。精確計(jì)算出水體的環(huán)境容量值可為采取相關(guān)措施控污和改善水質(zhì)提供技術(shù)支撐。

環(huán)境容量的概念最先是由日本學(xué)者借喻物理學(xué)上的電容概念并根據(jù)總量控制原則提出的，歐美國(guó)家則更多地使用日最大總負(fù)荷[1-3](total maximum daily load, TMDL)、同化容量[4-5](assimilative capacity)、污染負(fù)荷分配[6-7](wasteload allocation ,WLA)等術(shù)語(yǔ)。我國(guó)也在環(huán)境容量理論方面開(kāi)展了相應(yīng)的研究[8-10]，但目前對(duì)環(huán)境容量仍沒(méi)有達(dá)成完整統(tǒng)一的定義，目前學(xué)術(shù)界普遍認(rèn)可的水環(huán)境容量的概念定義為“某一水環(huán)境單元在特定的環(huán)境目標(biāo)下所能容納污染物的量”。其概念上強(qiáng)調(diào)了水環(huán)境容量的本質(zhì)三要素：水質(zhì)標(biāo)準(zhǔn)、水體條件和污染物特性。因此水環(huán)境容量的內(nèi)涵主要應(yīng)該有：其一，水環(huán)境容量是環(huán)境系統(tǒng)本身的自然屬性，是由自凈容量和稀釋容量組成的；其二，水環(huán)境容量是與污染物對(duì)環(huán)境造成的影響息息相關(guān)的，其大小是由水環(huán)境系統(tǒng)功能決定的?？傊?，水環(huán)境容量的目的是在水環(huán)境功能區(qū)水質(zhì)達(dá)標(biāo)的前提下，為制定污染物總量控制方案提供技術(shù)依據(jù)。

目前水環(huán)境容量的計(jì)算方法有零維、一維、二維、三維等方法，其中對(duì)大江大河、湖庫(kù)、近海等寬闊水域常采用二維、三維方法計(jì)算[11-13]；而針對(duì)中小型河流則常用零維、一維方法，并在中小河流環(huán)境容量的確定方法上也取得了一定的研究進(jìn)展。譬如：龔若愚等[14]采用了零維水環(huán)境容量計(jì)算公式，計(jì)算了柳江河新圩至窯埠段的COD和氨氮2項(xiàng)指標(biāo)的水環(huán)境容量值；王華等[15]基于零維方法對(duì)感潮河段進(jìn)行了環(huán)境容量計(jì)算；羅縉等[16]也是運(yùn)用零維計(jì)算方法建立了太湖流域平原河網(wǎng)區(qū)往復(fù)流河道的水環(huán)境容量模型；周孝德等[17]提出了計(jì)算河流水環(huán)境容量的一維方法，并應(yīng)用于渭河干流的環(huán)境容量計(jì)算；陳丁江等[18]在一維方法的基礎(chǔ)上建立了針對(duì)飲用水水源保護(hù)區(qū)內(nèi)河流環(huán)境容量的計(jì)算模型和方法；楊杰軍等[19]運(yùn)用河流環(huán)境容量的零維、一維計(jì)算方法核算了中國(guó)北方河道、攔河閘水體和感潮河段3種水體的水環(huán)境容量。但是中小型河流的水環(huán)境容量計(jì)算運(yùn)用零維和一維2種方法，仍存在以下明顯不足。

a. 零維方法。水環(huán)境容量計(jì)算公式為

(1)

式中:E——水環(huán)境容量，g/s；E1——稀釋容量，g/s；E2——自凈容量，g/s；Q——河段上游來(lái)水流量，m3/s；CS——水質(zhì)標(biāo)準(zhǔn)，mg/L；C0——上游來(lái)水本底濃度，mg/L；k——綜合降解系數(shù)，s-1；V——河段水體體積，m3。

零維方法計(jì)算水環(huán)境容量主要存在的問(wèn)題是：該方法認(rèn)為河道環(huán)境容量沿程分布是均勻的，沒(méi)有考慮沿程變化，同時(shí)也沒(méi)有考慮排污情況與排污位置等，因此是理論上的最大值，只具有一定的指導(dǎo)意義。而在實(shí)際天然河道中，河道容量沿程分布是沿程變化且不均勻的，同時(shí)考慮排污狀況在實(shí)際應(yīng)用時(shí)是很必要的，因此零維方法計(jì)算的水環(huán)境容量在實(shí)際天然河道的環(huán)境管理中直接應(yīng)用存在一定的偏差。

(2)

式中：q——排污口的平均排污流量，m3/s；l——河段長(zhǎng)度，m；u——河流平均流速，m/s。

一維方法較零維方法更符合實(shí)際情況，也是目前應(yīng)用最廣泛的方法。根據(jù)式(2)，一維方法可以理解和解釋為：目標(biāo)河段的水環(huán)境容量即為段首處污染物濃度與河段上游來(lái)水本底濃度之差所形成的容量，如圖1所示。

圖1　一維方法計(jì)算水環(huán)境容量示意圖Fig. 1　Sketch of one-dimensional method for calculating water environmental capacity

顯然，一維方法存在著明顯弊端。(a)該方法中排污口斷面到下游控制斷面全河段的污染物濃度都是超標(biāo)的，只在控制斷面段尾處水質(zhì)達(dá)標(biāo)，因此該河段水環(huán)境功能區(qū)全部不達(dá)標(biāo)，這是不符合水環(huán)境容量?jī)?nèi)涵的，也不利于污染物質(zhì)的削減與控制。(b)該方法夸大了目標(biāo)河段的自凈容量，即圖1中實(shí)際降解曲線上方部分的容量是不存在的；同時(shí)也不能真實(shí)地反映環(huán)境容量的沿程變化過(guò)程，應(yīng)用于生產(chǎn)實(shí)際通常會(huì)產(chǎn)生較大的偏差。(c)忽略了河流由于本身斷面不均勻產(chǎn)生的縱向分散效應(yīng)，但在實(shí)際河段中由于斷面不均勻?qū)е铝魉俨痪鶆蚴呛茏匀坏?，因此污染物的縱向分散作用常常是不能忽略的。

筆者針對(duì)目前常用于中小型河流的零維、一維環(huán)境容量計(jì)算方法中存在的弊端，認(rèn)為河流環(huán)境容量應(yīng)是沿程變化的，同時(shí)考慮了河流縱向分散作用和河流自身可能存在固定污染源的情況，提出了一種水環(huán)境容量計(jì)算的新方法。

1河流水環(huán)境容量計(jì)算新方法

1.1　水環(huán)境容量計(jì)算新方法

河流污染物在恒定情況的一維輸運(yùn)基本方程為

(3)

其中

式中：c——斷面污染物平均濃度，mg/L；x——縱向距離，m；Ex——縱向分散系數(shù)，m2/s；B——河段平均河寬，m；g——重力加速度，m/s2；h——河段平均水深，m；I——河底比降。

可得到考慮河流自凈和縱向分散作用的一維污染物輸運(yùn)方程的解：

(4)

式中：C1——段首處的污染物濃度，mg/L。

根據(jù)定義和內(nèi)涵，水環(huán)境容量是由自凈容量和稀釋容量?jī)刹糠纸M成的，前人研究顧及了當(dāng)上游來(lái)水的污染物濃度超出水質(zhì)標(biāo)準(zhǔn)時(shí)，河段也是不具備稀釋容量的，應(yīng)該存在負(fù)稀釋容量。但是，在河道中存在固定污染源的情況下，河流的污染物濃度是沿程變化的，這樣就造成了在全河段中不同河段部分的自凈容量是變化的；實(shí)際上該水質(zhì)超標(biāo)段是不應(yīng)具備自凈容量的，換言之，即應(yīng)該存在負(fù)自凈容量，故河段的自凈容量理應(yīng)考慮其沿程的變化。進(jìn)而推出最終的環(huán)境容量值。推導(dǎo)過(guò)程如下：

如圖2所示，對(duì)于河長(zhǎng)為l的中小型河段，若其段首處存在一固定排污濃度的污染源，根據(jù)河道容量的沿程變化，結(jié)合式(4)可得出在時(shí)間段t內(nèi)目標(biāo)河段的稀釋總量與自凈總量之和W，見(jiàn)式(5)。

圖2　水環(huán)境容量計(jì)算新方法示意圖Fig. 2　Sketch of new method for calculating water environmental capacity

(5)

式中：W1——稀釋總量，g；W2——自凈總量，g；A——河流斷面面積，m2；V1——河段上游來(lái)水體積，m3；V2——污水排放體積，m3；其他符號(hào)含義同前。

最后，計(jì)算河道的環(huán)境容量E(即為單位時(shí)間內(nèi)的W值):

(6)

特別地，河段水質(zhì)在段首處、段尾處達(dá)標(biāo)時(shí)的2種極限情形如下：

a. 當(dāng)C1=CS時(shí)，目標(biāo)河段在段首處水質(zhì)達(dá)標(biāo)，由于水體的降解分散作用，河段全段水質(zhì)達(dá)標(biāo)。此時(shí)嚴(yán)格控制了河段的水質(zhì)不超標(biāo)，要求太嚴(yán)格了，如圖3所示。計(jì)算公式如下：

(7)

圖3　全段達(dá)標(biāo)時(shí)水環(huán)境容量計(jì)算示意圖Fig. 3　Sketch of water environmental capacity calculation when water quality indicators throughout reach are within standards

(8)

圖4　全段超標(biāo)時(shí)水環(huán)境容量計(jì)算示意圖Fig. 4　Sketch of water environmental capacity calculationwhen water quality indicators throughout reach exceed standards

因此，合理的做法應(yīng)首先判斷該河段是否存在環(huán)境容量，并確定其存在正環(huán)境容量時(shí)的臨界條件，再進(jìn)行河段水環(huán)境容量的計(jì)算。

1.2　河道存在正環(huán)境容量的臨界判據(jù)

1.2.1河道存在正環(huán)境容量的最大允許排污濃度

對(duì)于實(shí)際河段，由于水環(huán)境容量的稀釋部分和自凈部分都是可能存在正或負(fù)的情況，可令目標(biāo)河段總的環(huán)境容量為零，求解出河段存在正環(huán)境容量的排污口臨界排放濃度，即最大允許排污濃度值。

令河段水環(huán)境容量E=0，可由式(6)推導(dǎo)出河段起始斷面濃度的最大值C1,max：

(9)

則可求得河段存在環(huán)境容量的排污口最大允許排放濃度：

(10)

1.2.2河道存在正環(huán)境容量的最大允許超標(biāo)長(zhǎng)度比例

當(dāng)河段總的水環(huán)境容量為零時(shí)，可計(jì)算出河道存在正環(huán)境容量的最大允許超標(biāo)長(zhǎng)度比例αmax：

(11)

1.3　計(jì)算步驟

采用本文計(jì)算方法的流程如圖5所示。

圖5　水環(huán)境容量計(jì)算流程Fig. 5　Flow chart of water environmental capacity calculation

2算 例 比 較

某河段長(zhǎng)l=25 km，流量Q=10 m3/s，流速u(mài)=0.2 m/s，本底濃度C0= 18 mg/L，其相應(yīng)的CODMn水質(zhì)目標(biāo)CS= 20 mg/L，CODMn綜合降解系數(shù)k=0.25/d。該河道段首處設(shè)有固定排污濃度q=0.8 m3/s的排污源，由經(jīng)驗(yàn)公式得到縱向分散系數(shù)Ex= 0.53 m2/s。分別設(shè)定3種不同的排污濃度背景值進(jìn)行算例比較，濃度值分別為80、150、200 mg/L。

根據(jù)本文提出的環(huán)境容量計(jì)算方法的步驟，首先依據(jù)上述基礎(chǔ)數(shù)據(jù)，由臨界判據(jù)公式(10)或(11)，得到當(dāng)存在正環(huán)境容量時(shí)的最大允許排污濃度為Cp,max=150 mg/L或該河道最大允許超標(biāo)長(zhǎng)度比例為αmax= 0.91。故判斷出排污濃度為80 mg/L時(shí)，該河段存在正值的環(huán)境容量；排污濃度為200 mg/L時(shí)，該河段已不具備環(huán)境容量。具體計(jì)算結(jié)果見(jiàn)表1。

表1　不同排污狀況下的環(huán)境容量計(jì)算值Table 1　Water environmental capacity calculation results under different drainage conditions

可以看出，不論排污濃度如何變化，對(duì)零維和一維方法的環(huán)境容量計(jì)算結(jié)果都是沒(méi)有影響的；而河流存在天然排污口的情況時(shí)，該河段的容量是沿程變化的，這也恰恰是決定河段是否存在環(huán)境容量的重要依據(jù)。對(duì)于算例3，排污濃度過(guò)大已經(jīng)導(dǎo)致全河段的水質(zhì)超標(biāo)，環(huán)境容量已為負(fù)值，根本就不具備環(huán)境容量了，因此零維、一維方法的結(jié)果是不符合天然實(shí)際情況的。本文新方法的計(jì)算結(jié)果可以反映實(shí)際河段水環(huán)境容量的沿程變化，將臨界排污濃度或最大允許超標(biāo)長(zhǎng)度比例作為評(píng)判是否存在河流容量的依據(jù)，消除了自凈容量的夸大部分，結(jié)果更接近環(huán)境容量的真實(shí)值，更具實(shí)用性和可操作性。

3結(jié)論

本文以水環(huán)境容量的定義和本質(zhì)內(nèi)涵為基石，剖析并指出了常用于中小型河流的零維、一維環(huán)境容量計(jì)算方法的不足和弊端；根據(jù)河流的天然情況和排污狀況，從水環(huán)境容量的內(nèi)涵出發(fā)，考慮了河流環(huán)境容量的沿程變化過(guò)程，推導(dǎo)出一種計(jì)算河道環(huán)境容量的新方法，使其與實(shí)際情況更為契合；在此基礎(chǔ)上，提出了河道存在環(huán)境容量時(shí)的臨界判據(jù)，分別得到排污口最大允許排放濃度、河道最大允許超標(biāo)長(zhǎng)度比例的值。通過(guò)算例和比較分析，證明了該計(jì)算新方法的結(jié)果更符合客觀實(shí)際情形，也易于操作應(yīng)用。該計(jì)算方法可為環(huán)境管理部門(mén)更準(zhǔn)確地確定環(huán)境容量和合理控制排污口源強(qiáng)提供可靠的技術(shù)支撐。

參考文獻(xiàn)：

[1] SAHOO G B, NOVER D M, REUTER G E, et al. Nutrient and particle load estimates to Lake Tahoe (CA-NV, USA) for total maximum daily load establishment[J]. Science of the Total Environment,2013, 444(2): 579-590.

[2] CHEN D, DAHLGRENC R A, SHEN Y A. Bayesian approach for calculating variable total maximum daily loads and uncertainty assessment[J]. Science of the Total Environment, 2012, 430(14): 59-67.

[3] KIM J, ENGELB B A, PARKA Y S, et al. Development of Web-based Load Duration Curve system for analysis of total maximum daily load and water quality characteristics in a water-body[J]. Environmental Management, 2012, 97(2): 46-55.

[4] SEHERER C R．On the efficient allocation of environmental assimilative capacity: the case of thermal emissions to large body of water[J]. Water Resource Research, 1975, 11(1): 180-181.

[5] GOYAL P, ANAND S, GERA B S. Assimilative capacity and pollutant dispersion studies for Gangtok city[J]. Atmospheric Environment, 2006, 40(9): 1671-1682.

[6] CHADDERTON R A, MILLER A C, MCDONNELL A J. Analysis of wasteload allocation procedures[J]. Water Resources Bulletin, l981, 17(5): 760-766.

[7] MAEDA S, KAWACHI T, UNAMI K, et al. Controlling wasteloads from point and nonpoint sources to river system by GIS-aided epsilon robust optimization model[J]. Journal of Hydro-environment Research, 2010, 4(1): 27-36.

[8] 張永良, 劉培哲. 水環(huán)境容量綜合手冊(cè)[M]. 北京: 清華大學(xué)出版社, 1991.

[9] 關(guān)卉. 湛江市水環(huán)境容量分析與應(yīng)用研究[M].北京: 中國(guó)環(huán)境科學(xué)出版社，2005.

[10] 張永良. 水環(huán)境容量基本概念的發(fā)展[J]. 環(huán)境科學(xué)研究, 1992, 5(3)：59-61.(ZHANG Yongliang. The development of basic concept of water environmental capacity[J]. Research of Environmental Sciences, 1992, 5(3)：59-61. (in Chinese))

[11] 宿俊英, 劉樹(shù)坤, 何少答, 等. 太湖水環(huán)境容量的研究[J]. 水利學(xué)報(bào), 1992, 23(11)：22-36.(SU Junying, LIU Shukun, HE Shaoda, et al. Water quality management in Tai Lake[J]. Journal of Hydraulic Engineering, 1992, 23(11): 22-36.(in Chinese))

[12] 陳益明.規(guī)劃環(huán)評(píng)中海灣水環(huán)境容量的計(jì)算方法[J]. 水資源保護(hù),2014,30(3):70-75(CHEN Yiming.A method for calculation of environmental capacity of bay water during planning environmental impact assessment[J].Water Resources Protection，2014,30(3):70-75.(in Chinese))

[13] 黃真理, 李玉粱, 李錦秀, 等. 三峽水庫(kù)水環(huán)境容量計(jì)算[J]. 水利學(xué)報(bào)，2004,35(3)：7-14.(HUANG Zhengli, LI Yuliang, LI Jingxiu, et al. Water environmental capacity for the reservoir of three gorges project[J]. Journal of Hydraulic Engineering, 2004,35(3)：7-14.(in Chinese))

[14] 龔若愚, 周源崗. 柳江柳州段水環(huán)境容量研究[J]. 水資源保護(hù),2001,17(1): 31-32.(GONG Ruoyu, ZHOU Yuangang. Study on capacity of water environment of Liuzhou Section of Liujiang River[J].Water Resources Protection，2001,16(1)：31-32.(in Chinese))

[15] 王華, 逄勇, 丁玲. 濱江水體水環(huán)境容量計(jì)算研究[J].環(huán)境科學(xué)學(xué)報(bào), 2007, 27(12): 2067-2073.(WANG Hua, PANG Yong, DING Ling, Calculation of the water environment capacity for a water front body[J]. Acta Scientiae Circumstantiae, 2007, 27(12): 2067-2073.(in Chinese))

[16] 羅縉, 逄勇, 羅清吉. 太湖流域平原河網(wǎng)區(qū)往復(fù)流河道水環(huán)境容量研究[J].河海大學(xué)學(xué)報(bào)：自然科學(xué)版, 2004, 32(2): 144-146.(LUO Jin, PANG Yong, LUO Qingji, et al. Study on water environment capacity for reversing current channels in plain river network region in Taihu Lake Basin[J]. Journal of Hohai University：Natural Sciences, 2004, 32(2): 144-146.(in Chinese))

[17] 周孝德, 郭瑾瓏, 程文, 等. 水環(huán)境容量計(jì)算方法研究[J].西安理工大學(xué)學(xué)報(bào)，1999,15(3): 1-6.(ZHOU Xiaode, CHENG Jinong, CHENG Wen, et al. The comparison of the environmental capacity calculation methods[J].Journal of Xi’an University of Technology, 1999,15(3): 1-6.(in Chinese))

[18] 陳丁江, 呂軍, 沈曄娜, 等. 飲用水水源保護(hù)區(qū)河流水環(huán)境容量計(jì)算模型[J].環(huán)境科學(xué)，2008，29(9): 2437-2440.(CHEN Dingjiang , LYU Jun, SHEN Yena, et al. Water environmental capacity calculation model for the rivers in drinking water source conservation area[J]. Environmental Science, 2008，29(9): 2437-2440.(in Chinese))

[19] 楊杰軍, 王琳, 王成見(jiàn), 等. 中國(guó)北方河流環(huán)境容量核算方法研究[J].水利學(xué)報(bào)，2009, 40(2)：194-200.(YANG Jiejun, WANG Ling, WANG Chengjian, et al. Environmental capacity of rivers in North China[J]. Journal of Hydraulic Engineering, 2009, 40(2)：194-200.(in Chinese))

A method for estimating water environmental capacity

based on its variation along river

HUA Zulin1, 2, 3, CHENG Haomiao1, 2

(1.KeyLaboratoryofIntegratedRegulationandResourceDevelopmentonShallowLakesofMinistryofEducation,

HohaiUniversity,Nanjing210098,China;

2.CollegeofEnvironment,HohaiUniversity,Nanjing210098,China;

3.NationalEngineeringResearchCenterofWaterResourcesEfficientUtilizationandEngineeringSafety,

HohaiUniversity,Nanjing210098,China)

Abstract:Based on the definition and connotation of water environmental capacity, the shortcomings of currently prevalent zero-dimensional and one-dimensional methods for estimating the water environmental capacity of medium- and small-sized rivers are pointed out. According to these shortcomings, a new method considering the variation of water environmental capacity along a river is put forward. This method eliminates the exaggerated part of the environmental capacity estimation generated with the previous methods, and it is more practical. Moreover, the critical judgment criteria for estimating the existence of water environmental capacity are derived, including the maximum allowable concentrations of discharged pollutants and the maximum ratio of the length of the river reach with water quality indicators exceeding the standards to the total length of the river. Finally, practical cases are provided in order to compare the proposed method with previous methods.

Key words:water environmental capacity; variation along river; calculation method; critical judgment criteria

中圖分類(lèi)號(hào)：X824

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-1980(2015)06-0505-06

作者簡(jiǎn)介：華祖林(1965—)，男，江蘇江陰人，教授，主要從事水環(huán)境模擬與污染物輸運(yùn)機(jī)理研究。E-mail:zulinhua@hhu.edu.cn

基金項(xiàng)目：水體污染控制與治理科技重大專(zhuān)項(xiàng)(2012ZX07103-005)；“十二五”國(guó)家科技支撐計(jì)劃(2012BAB03B04)；國(guó)家自然科學(xué)基金(51179052)；上海市水務(wù)局科研項(xiàng)目(2014-10)

收稿日期：2015-06-04

DOI：10.3876/j.issn.1000-1980.2015.06.001