多閘壩河流水質(zhì)模擬

卓全錄，許 洋

(1.濮陽市環(huán)境保護科學研究所，河南 濮陽 457000；2.許昌市生態(tài)環(huán)境局，河南 許昌 461000)

1 引言

近年來水環(huán)境數(shù)學模型得到了長足的發(fā)展，并在此基礎上開發(fā)出了諸如MIKE[1]、Delft3D[2]、QUAL2K[3]和EFDC[4]等軟件，EFDC因其廉價、開源、高效被廣泛應用，雖在國內(nèi)的應用推廣較晚，但不乏成功的案例。相關學者借助EFDC模型對丹江口水庫、密云水庫、洱海、太湖、長江及長江口等進行了卓有成效的研究[5～10]。

利用EFDC的dye模塊對河流水質(zhì)進行一維模擬是一種簡單有效的水質(zhì)模擬方法，但該模式下模擬河流只能定義一個統(tǒng)一的衰減系數(shù)，而現(xiàn)實河流中往往閘壩較多，在閘壩密集的河段和無閘壩河段，河流污染物的衰減系數(shù)相差較大，因此往往會影響EFDC一維模擬的精準度。

修建閘壩等水利工程是人類開發(fā)和利用水資源的有效手段，閘壩通過對水量的調(diào)控實現(xiàn)防洪、供水、灌溉、發(fā)電、航運等功效的同時，污染物質(zhì)輸移過程也相應發(fā)生變化，從而對水質(zhì)產(chǎn)生一定影響。第一次全國水利普查成果顯示：截至2011年，中國共有水庫98002座、水電站46758座、過閘流量1m3/s及以上水閘268476座[11]。關于閘壩的界定，沒有權(quán)威的文獻進行論述。我國學者索麗生[12]、夏軍[13]、林巍[14]、張永勇[15,16]、蔣艷[17]、左其亭[18,19]等在多個文獻中引用“閘壩”一詞。本文繼續(xù)沿用“閘壩”一詞，并界定為：具有防洪、排澇、拒咸蓄淡、灌概供水、通航養(yǎng)殖、景觀娛樂、生態(tài)保護等功能的平板閘、液壓升降閘、橡皮壩等水工構(gòu)筑物。從近期研究看，陳煉鋼[20、21]、呂菲菲[22]、劉子輝[23]、楊洵[24]、米慶彬[25]等學者均對多閘壩河流的水質(zhì)模擬進行了研究，但使用MIKE和自研模型較多，采用EFDC的較少。基于EFDC探索一種簡單的多閘壩河流水質(zhì)模擬方法對于河流水質(zhì)模擬具有重要意義，本文針對清潩河，采取分段模擬的方法解決了多閘壩河流EFDC一維模擬不準確的問題，對相關研究具有指導意義。

2 清潩河流域概況

清潩河屬淮河流域沙潁河水系，發(fā)源于新鄭市，流經(jīng)長葛市、許昌縣、魏都區(qū)(許昌市區(qū))、臨潁縣、鄢陵縣等縣(市、區(qū))，至鄢陵縣陶城閘下匯入潁河，全長149 km，流域面積2362 km2，其中約67%在許昌境內(nèi)，主要擔負沿河區(qū)域防洪、排澇、納污及許昌市城區(qū)景觀水的重要功能。

清潩河干流上自增福廟水庫以下至匯入潁河口共有8閘2壩，基本分布于長葛市、許昌市及許昌縣段，10個閘壩中有9個均位于此范圍內(nèi)。在用(處于調(diào)控狀態(tài))的閘壩為關莊閘、小黃橋閘、滹沱閘、建安大道橡膠壩、前進路橡膠壩、祖師廟閘、陶城閘。閘壩的主要功能為調(diào)蓄清潩河許昌市區(qū)段水質(zhì)和水量。

3 模型概化方法

3.1 河流概化

結(jié)合清潩河流域河流水系情況和污染源分布情況，對清潩河流域進行概化。現(xiàn)狀情況下，清潩河干流長葛段無水，自杜村寺閘下游才有杜村寺濕地排水，且杜村寺閘-祖師廟閘段河流基本處于靜穩(wěn)狀態(tài)，河水較深，與祖師廟下游段區(qū)別明顯，在進行清潩河水動力和水質(zhì)模擬的過程中，針對全段進行模擬較為困難，且不符合客觀實際。因此，將清潩河干流分成兩段進行模擬，即上段(杜村寺閘-祖師廟閘段)和下段(祖師廟閘-匯入潁河口段)，將上段模擬得出的河流出口處流量和水質(zhì)作為下段進口處的邊界數(shù)據(jù)。

在清潩河干流由上游到下游選取點位進行河底高程監(jiān)測。將對應點位高程信息導入EFDC并進行插值，得出清潩河干流的高程分布。

河道糙率選用季振剛的《水動力學和水質(zhì)——河流、湖泊及河口數(shù)值模擬》一書中的河道糙率經(jīng)驗值進行設置[26]。清潩河上段大部分區(qū)域已經(jīng)過整治，河道情況已非自然狀態(tài)，介于普通的混凝土涂層和最優(yōu)條件下的泥土河道之間，其糙率取0.015。清潩河下段基本屬于天然狀態(tài)，具有平原地區(qū)天然河流的一般性特征，其糙率取為0.03。

3.2 閘壩設置

清潩河模擬河段內(nèi)對清潩河水流狀態(tài)有明顯影響的閘壩為小黃橋閘、建安大道橡膠壩、前進路橡膠壩、祖師廟閘和陶城閘。小黃橋閘、建安大道橡膠壩、前進路橡膠壩、祖師廟閘均為堰流，陶城閘為閘孔出流。在閘壩處需設置水工構(gòu)筑物邊界，將該水工構(gòu)筑物的過水流量導入模型。堰流和閘孔出流的流量計算方法如下：

3.2.1 堰流流量計算

堰流流量計算公式為：

(1)

式(1)中，μ為淹沒系數(shù)，取1；ε為側(cè)收縮系數(shù)，取1；m為堰的流量系數(shù)(取0.36)；b為寬度，m；g為重力加速度，m/s2；H0為堰頂總水頭，m，當堰前流速很小的情況下，H0近似為堰前水深和堰高的差值。

3.2.2 閘孔出流流量計算

閘孔出流流量計算公式為：

(2)

(3)

式(2)、(3)中，μ為流量系數(shù)；n為閘門開孔數(shù)；b為每孔的寬度，m；e為閘門開啟高度，m；H0為閘前全水頭，m，當閘前流速很小的情況下，H0近似為閘上水深；H為閘上水深，m。

4 水動力模擬

用2017年全年的數(shù)據(jù)進行模型校準。清潩河上段沒有流量、流速等監(jiān)測數(shù)據(jù)，實際調(diào)研中發(fā)現(xiàn)，清潩河上段由于閘壩眾多、蓄水量大，基本沒有流速，處于靜穩(wěn)狀態(tài)，所以清潩河上段水動力模型校準主要以水深對比為依據(jù)。清潩河下段臨潁高村橋斷面有流量監(jiān)測數(shù)據(jù)，因此清潩河下段主要以流量對比為依據(jù)進行水動力模型校準。

利用2017年數(shù)據(jù)模擬得出清潩河上段沿程水深變化。小黃橋閘上游由于沒有穩(wěn)定的匯流，導致在小黃橋閘處河水無法溢流；建安大道橡膠壩和前進路橡膠壩處由于橡膠壩的截流作用，上游河水被雍高約5 m深；祖師廟閘處由于閘板的作用上游水位雍高約3.6 m，閘板上有水溢流。以上模擬情況與河流實際情況吻合，清潩河上段水動力模型較為準確。

利用2017年數(shù)據(jù)對清潩河下段水動力進行模擬，將臨潁高村橋斷面模擬流量與實測流量進行對比，如圖1所示。實測流量與模擬流量較為吻合，平均誤差約為9%。可見清潩河下段的水動力模型較為準確。

4 水質(zhì)模擬

采用2017年的數(shù)據(jù)對清潩河進行模擬，經(jīng)過調(diào)整，將清潩河上段COD的衰減系數(shù)設置為0/d，下段COD的衰減系數(shù)設置為0.1/d。模擬得出清潩河107國道橋、鄢陵陶城閘2個斷面全年的水質(zhì)濃度變化，將其與實測值對比，如圖2～3所示。2個斷面的平均相對誤差分別為12.8%、4.5%。誤差較小，COD水質(zhì)模型較為準確。

圖1 清潩河臨潁高村橋斷面流量對比

圖2 清潩河107國道橋斷面COD濃度對比

圖3 清潩河鄢陵陶城閘斷面COD濃度對比

5 結(jié)論

(1)在分段模擬的情況下，清潩河多閘壩的上段衰減系數(shù)設置為0/d時，水質(zhì)模型較為準確，而下段衰減系數(shù)設置為0.1/d時，水質(zhì)模型較為準確。

(2)同時可以看出，因景觀用水需求清潩河許昌市區(qū)段設置大量閘壩，河水較深，流動遲緩，處于靜穩(wěn)狀態(tài)，加之河底和邊坡硬化導致河流失去自凈能力，而清潩河下段更為接近天然河流，對污染物仍保持一定的自凈能力。

本研究針對多閘壩河流，采用分段模擬的方法，解決了采用EFDC的dye模塊對河流的一維模擬中不能分段設置衰減系數(shù)的問題。該方法簡單可行，結(jié)果可靠，對相似河流的模擬具有指導意義。