三峽庫區漢豐湖流域水環境容量實證研究

何家軍 樊連生 張慧

摘要?采用二維水質模型對三峽庫區漢豐湖流域流量、水位、污染源排放、水動力和水質進行模擬，依據國家標準和規范計算了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水環境容量超負荷現象。