韓亞萍,郭亞男
(河南省水利勘測設計研究有限公司,鄭州 450016)
生態環境是人類生存和發展的根基,生態環境變化直接影響文明興衰演替。城市水生態環境破壞問題日益凸顯,嚴重制約著我國經濟社會的可持續發展。針對我國城市廣泛的水環境污染問題,各地區大力開展水環境綜合治理,著力解決突出的水環境問題,努力實現水清綠岸、魚翔淺底的美好愿景。本文以桃花塢水庫為例,利用灰水足跡理論對現狀的水污染程度進行評價,并結合桃花塢水庫污染源分析采取水環境綜合治理,通過灰水足跡理論與Mike 21水質模型的聯合運用對其治理效果進行綜合評價。
灰水足跡理論從水量的角度來評價水污染的程度,為定量評價水污染程度提供了新的思路[1]。但是,灰水足跡理論是基于污染物總量核算,缺少對水動力擴散、水體自凈能力的考慮[2]。而水質模型與水動力模型耦合后可以反映出水質的演變過程,將水質模型與灰水足跡理論相結合運用,能夠彌補灰水足跡未能考慮水動力擴散這一不足之處,可全面綜合地評價區域的水環境狀況。
桃花塢水庫位于河南省信陽市固始縣城,史河支流鄭小河、小史河上,是一座以防洪、城市工業用水為主,結合水產養殖、休閑旅游等綜合利用的小(Ⅰ)型水庫。鄭小河發源于根親文化園南入口的一處坑塘,向北作為根親文化園的主要水系橫穿文化園,最終匯入桃花塢水庫。桃花塢水庫上游為鄭小河,下游為小史河,承接梅山灌區南干渠下游退水作為調節水源。桃花塢水庫控制流域面積為5.82 km2,50年一遇洪水設計,1 000年一遇洪水校核,校核庫容178×104m3。
桃花塢水庫和上游鄭小河位于老城區,污水管網不健全,沿河居民產生的生活污水直排入河,河道淤積嚴重。桃花塢水庫沿岸口徑較大的污水排放口共有15處,沿庫周邊約2萬居民,桃花塢水庫上游與鄭小河水系不連通,造成桃花塢水庫缺乏清潔的補水水源。根據對桃花塢水庫的抽樣檢測,桃花塢水庫NH3-N、COD、BOD5此3項檢測值均大于地表水環境Ⅴ類標準值,因此桃花塢水庫水質屬于劣Ⅴ類。
灰水足跡的概念在2008年由Hoekstra和Chapagain提出,并將其定義為排放到水體的污染負荷除以最大容許濃度[3-4]。通過引入灰水足跡理論計算,從水量的角度評價水污染的程度,量化地分析區域水資源對水質的影響,利用灰水足跡可持續來評價區域水污染程度,可根據不同地區的水資源稀缺程度來量化分析同一污染物在不同區域水環境造成的污染影響。
灰水足跡的計算通常采用《水足跡評價手冊》的方法直接核算[5],采用將污染物稀釋至水質標準所需水量進行表示,其計算公式如下:
(1)
式中:WFgrey為灰水足跡,m3/a;L為區域內污染物排放總量,kg/a;Cmax為水體達到某一環境標準所容許的污染物最大濃度,kg/m3;Cnat為受納水體的自然本底濃度,kg/m3。
根據灰水足跡可持續評估的計算方法[5],流域內水污染程度(WPL)是評價區域內灰水足跡可持續的指標,WPL等于一個流域的灰水足跡(max[WFgrey(P)])與流域徑流量(Ract) 的比值。
(2)
流域水污染程度與流域內徑流成反相關,Ract越小,相對的WPL就越大,說明流域內徑流量越小,水污染程度越大,灰水足跡不可持續性的概率也越大。當WPL大于1,說明該流域內水體的納污能力消耗殆盡,能容納污染物的能力也基本消失,水體呈現污染且不可持續發展狀態。
2.2.1 參數選取
桃花塢水庫在城市區,匯入的污染物主要是以點源方式的直排生活污水,面源匯入的地表徑流污染相對比點源來說含量較小,可忽略不計。因此,本次核算主要考慮點源污染,污染物類型為桃花塢水庫抽樣檢測的超標因子COD和NH3-N。
污染物濃度參照當地污水處理廠進水濃度,COD和NH3-N濃度為271和25.6 mg/L,沿庫周邊約2萬人的生活污水不經處理進入庫區,參照《固始縣城市排水防澇綜合規劃》,污水量為用水量的0.8,日變化系數取1.3。桃花塢水庫需要達到地表水Ⅳ類標準,即Cmax(COD)為30 mg/L、Cmax(NH3-N)為1.5 mg/L。COD和NH3-N指標的自然本底值參照《地表水環境質量標準》(GB 3838-2002)中的Ⅰ類標準濃度限值,Cnat(COD)為15 mg/L、Cnat(NH3-N)為0.015 mg/L。
Ract為桃花塢水庫年均徑流總量186.2×104m3,桃花塢水庫及流域內無水文觀測站,無實測流量資料。參照河南省水文水資源局2007年8月編繪的《河南省水資源》,查算流域內多年平均徑流深320 mm左右,且年際變化劇烈,本文將多年平均徑流量作為流域徑流量。
2.2.2 計算結果
根據式(1),WFgrey(COD)為970.34×104m3/a,WFgrey(NH3-N)為1 016.44×104m3/a,選取較大的污染物數值作為該區域的灰水足跡,故桃花塢水庫的灰水足跡為1 016.44×104m3/a。根據式(2),WPL為5.46,水污染程度大于1,說明該區域灰水足跡不可持續。
桃花塢水庫有大量生活污水未經處理直接排入庫區,因此造成桃花塢水庫灰水足跡比較大;桃花塢水庫上游與鄭小河現狀水系不連通,上游徑流不能匯入該研究區域,流域內徑流量不足,遠小于區域承受的灰水足跡量,桃花塢水庫水污染程度大于1,灰水足跡呈現不可持續性。
根據灰水足跡計算結果,本方案采用“控源截污、內源治理、水系連通”的思路對該區域水環境進行治理,桃花塢水庫水污染程度降低可以通過減小灰水足跡和增大徑流量兩個方式入手,即分別是控源截污削減污染物排放量和水系連通進行補源,同時為了減少內源污染物的釋放,削減污染物,增加水體的自凈能力,在該區域進行內源清淤。
控源截污措施是為了防止外來的各種污染物直接或隨雨水排入水體,是水環境整治最直接有效的工程措施,也是采取其他技術措施的前提[6]??卦唇匚鄞胧┲饕窍一▔]水庫的點源污染,根據現場調查,桃花塢水庫周邊無污水管道,兩岸居民生活污水都是直排入庫。對成功大道-桃花塢路之間水庫沿庫周圍分別在東西向布置兩個截污管道,東西兩個管道系統各服務人口約1萬人,污水管網分別沿桃花塢水庫岸線東西兩側由南向北布置,最后與城市已有的污水干管連通,進入城市污水處理廠。
底泥疏浚能夠永久性去除富含污染物的底泥,從而在源頭上有效控制內源污染物的釋放,提升區域的自凈能力[7]。為了快速降低桃花塢水庫的內源污染負荷,需要對桃花塢水庫進行清淤。基于MIKE 21 FM搭建二維水動力模型,結合地形資料,對庫區進行合理概化。從桃花塢水庫的流場分布圖(圖1)可以看出,桃花塢水庫引水時,上中游的流速較大,流速為0.05~0.1 m/s,對庫區底部淤泥沖刷作用相對較強,下游流速偏小且流態穩定,流速為0.005~0.05 m/s,對庫區底部淤泥的擾動較??;同時庫底淤積層相對穩定,污染物在不擾動情況下不易釋放。根據引水時流場分布情況,治理工程樁號0+606(化肥廠家屬院東南側)以下庫區水體流速小,底泥擾動輕微,底泥對水體水質影響??;且庫底淤積層相對穩定,污染物在不擾動情況下不易釋放,該區域可不進行清淤。因此,確定清淤范圍為化肥廠家屬院以上區域(圖2),清淤深度0.5~1.5 m,清淤面積約7.5×104m2,清淤量約9.8×104m3。

圖1 桃花塢庫區流場分布圖

圖2 桃花塢水庫清淤范圍圖
水系連通對河道水質改善作用明顯,能夠增加水動力條件,縮短水體水力停留時間,改善水系水質[8]。 本方案通過將梅山南干渠和鄭小河連通,增加桃花塢水庫的生態補水量,提高水體的流動性,增加水體的水環境容量,桃花塢水庫生態需水量為317.2×104m3。在梅山灌區南干渠九支渠下游60 m處建引調水建筑物(此處上游20 m處建有節制閘,可以控制引調水不會沿渠道向渠道上游回水),將鲇魚山東干渠水源引調至梅山灌區南干渠,將引調水引至鄭小河中并進入桃花塢水庫。
根據灰水足跡計算方法及工程措施,工程實施后該部分污水被截留進入污水處理廠,出水經處理后達到一級A排放標準進入水體,即污染物入河濃度COD和NH3-N分別為50和5 mg/L;Cmax和Cnat參數不變;Ract為桃花塢水庫連通后通過梅山南干渠可補給的水量317.2×104m3。
根據式(1),WFgrey(COD)為179.69×104m3/a,WFgrey(NH3-N)為181.51×104m3/a,選取較大的污染物數值作為該區域的灰水足跡,故桃花塢水庫的灰水足跡為181.51×104m3/a。根據式(2),WPL為0.57,該數值小于1,說明該區域灰水足跡持續。
通過控源截污、內源治理的工程措施后,污染物排放量得到削減,桃花塢水庫的灰水足跡從1 016.44×104m3/a降低至181.51×104m3/a。通過水系連通工程措施后,桃花塢水庫徑流量增加,能夠用來稀釋污染物的水源變大,桃花塢水庫的水污染程度由5.46降低至0.57,實現了區域灰水足跡的可持續性。
4.2.1 水質濃度變化分析評價
根據MIKE 21二維生態水質的模擬,將引水、清淤、截污等工程措施添加到模型的外邊界條件,模型中考慮降雨蒸發及大氣沉降的影響,并根據換水運行工況來設置入流的水質邊界,引水水質采用梅山南干渠實測水質資料。在每年的3、6、8和11月初進行補水,補水流量1.5 m3/s,全部換水一次需6 d。該工況主要水質指標值見表1,模擬計算得到TP、TN及Chla濃度的年變化情況見圖3。

表1 主要水質指標模擬結果

圖3 模擬水質濃度年變化曲線
藻類生長受溫度影響較大,溫度越高藻類生長越快,5-10月份,隨著溫度的升高藻類生長速率加快,在水體內累積的營養鹽也相應增加。而在補水換水的時候,由于稀釋作用營養鹽濃度有一定的降低,10月份以后溫度降低藻類不適宜生存而開始衰亡,TN、TP等營養鹽濃度也迅速降低。
根據《地表水環境質量標準》(GB 3838-2002),水質模型模擬出的TN、TP指標年度變化值最大值分別是1.45 和0.07 mg/L,均在地表水Ⅳ類標準以內。根據《地表水資質量評價技術規程》(SL 395-2007),水質模型模擬出的Chla指標年度變化值最大值是0.025 mg/L,小于湖泊富營養狀態時葉綠素濃度要求。
4.2.2 水華風險分析
ECOLab是DHI開發的生態數值模擬軟件,主要回來描述水環境中的物理、化學、生物及生態過程,與Mike11或者Mike21模型耦合一起解決水質問題,對湖泊和水庫的富營養化問題進行模擬[9]。
由于水溫從年初逐漸升高,到夏季7、8和9月份水溫達到25℃左右,處于藻類生長適宜溫度范圍內,隨后溫度下降。本文對桃花塢水庫進行模型構建,模擬高峰期8月份的Chla濃度變化情況,通過在連續入流條件下的模擬,得到高溫、光照等自然條件及營養鹽濃度較充足情況下的Chla濃度月變化曲線(圖4)。Chla濃度持續升高,經過約25 d后Chla濃度為0.025 mg/L,小于湖泊富營養狀態時葉綠素濃度0.026 mg/L的要求,可以認為庫區的水華風險程度為低風險。

圖4 Chla濃度8月份變化曲線
1) 利用灰水足跡理論和水質模型對研究區域水環境措施效果進行評價,工程實施后灰水足跡降低至原來的17.8%,水污染程度從5.46降低至0.57,灰水足跡由不可持續變為可持續;根據水質模型結果,TN和TP指標的濃度全年都在地表水Ⅳ類標準內,水華風險處于低風險。
2) 利用水質模型來彌補單一的灰水足跡計算未考慮水體自凈能力這一因素,水質模型驗證結果與灰水足跡理論計算結果一致,說明本次治理措施能夠有效提升改善桃花塢水庫的生態環境。
3) 水污染程度與徑流量關系密切,年際變化下或者年內豐、平、枯水期條件變化下水污染程度均不相同,本文取多年平均徑流量核算一個多年平均的水污染程度。