武春芳,徐明德,李 璐,景勝元 (太原理工大學(xué)環(huán)境科學(xué)與工程學(xué)院,山西 太原 030024)

太原市迎澤湖富營養(yǎng)化控制的模型研究

武春芳,徐明德*,李 璐,景勝元 (太原理工大學(xué)環(huán)境科學(xué)與工程學(xué)院,山西 太原 030024)

針對湖泊富營養(yǎng)化問題,對太原市迎澤湖進行資料收集和長期水體監(jiān)測,綜合水動力模型、水質(zhì)模型,融合參與湖泊富營養(yǎng)化的各種生命活動過程,建立湖泊富營養(yǎng)化耦合模型,用實測水質(zhì)數(shù)據(jù)進行參數(shù)率定與驗證,選取總氮、總磷、葉綠素a及透明度4項因子進行模擬,得出迎澤湖營養(yǎng)物質(zhì)輸移擴散及時空分布規(guī)律.結(jié)合迎澤湖實際情況,從補水方式、補水頻率以及改變湖泊柔性結(jié)構(gòu)三方面,提出改變水動力條件的方案并進行數(shù)值模擬,研究了各方案水動力、物質(zhì)輸移擴散的改進效果,結(jié)果表明,作為富營養(yǎng)化程度表征的葉綠素a在空間和時間分布上均存在一定的規(guī)律性,其濃度變化范圍在 0.035～0.105mg/L之間;藻類等浮游植物迅速增殖導(dǎo)致水質(zhì)惡化水華暴發(fā),而采取加大湖水水力循環(huán)和改善入口水質(zhì)的控制方案,可有效改善湖水水質(zhì).

迎澤湖；水動力模型；水質(zhì)模型；富營養(yǎng)化耦合模型；MIKE 21

城市景觀湖作為城市基礎(chǔ)設(shè)施的重要組成部分,具有排泄雨洪、供水、旅游觀光等綜合功能,它對于城市生態(tài)建設(shè)是十分重要的,因此,城市中小型淺水湖泊的富營養(yǎng)化問題不容忽視.

湖泊水動力和污染物輸移擴散研究[1]較早,但就富營養(yǎng)化研究而言,大都集中于對自然湖泊和中大型人工水庫的模擬研究[2-5],對中小型水域的模擬研究[6-8]較少,對干旱地區(qū)內(nèi)陸景觀湖的模擬研究[9-10]則更少,研究所使用的模型經(jīng)歷了從簡單的回歸模型、營養(yǎng)物平衡模塊到復(fù)雜的水動力、水質(zhì)綜合模型和生態(tài)結(jié)構(gòu)動力學(xué)模型的過程[11].在國內(nèi)外湖泊富營養(yǎng)化研究中,建立耦合模型是一個新的發(fā)展方向.本研究針對干旱地區(qū)城市景觀湖特點,以影響湖泊富營養(yǎng)化的有機、無機營養(yǎng)物及浮游生物生命活動過程[12]等為基礎(chǔ),建立富營養(yǎng)化耦合模型,以太原市迎澤湖為載體,用該耦合模型對湖水水質(zhì)演變和保障措施實施效果進行模擬預(yù)測,以期為景觀湖水質(zhì)管理提供有效的決策支持.

1 研究對象

1.1 監(jiān)測點布設(shè)

選取太原市迎澤湖為研究對象.迎澤湖為太原市區(qū)最大的人工景觀湖,位于太原市中心,根據(jù)湖體形狀將其分為南湖和北湖兩湖區(qū),南湖長1060m,東西最寬處為 631m,最窄南北湖區(qū)接連處為 54m,其水源主要來自降雨和汾河二庫的定期補水,二庫來水時進水口設(shè)在北湖區(qū)最北端,出水口設(shè)在南湖區(qū)最南端,平均水深 1.98m,全湖面積22.23hm2.

針對迎澤湖湖體形狀和水體流動特點,在北湖和南湖各設(shè)一個取樣點,分別記為A、B,見圖1. 1.2 監(jiān)測項目及方法

基于影響湖泊富營養(yǎng)化主要因素,選擇溶解氧DO、總氮(TN)、總磷(TP)作為監(jiān)測項目.監(jiān)測方法參照《水和廢水監(jiān)測分析方法》[13],分別選取碘量法、過硫酸鉀氧化紫外分光光度法、鉬酸銨分光光度法進行測量.

1.3 監(jiān)測時限及頻率

圖1 迎澤湖水質(zhì)監(jiān)測點示意Fig.1 The monitor position of Yingze lake

綜合考慮氣溫?水溫變化及藻類植物生長發(fā)展規(guī)律,根據(jù)湖泊夏季水質(zhì)較差?易于發(fā)生富營養(yǎng)化的特點,監(jiān)測時間選為 6～11月,頻率為每月采樣2次.

1.4 監(jiān)測結(jié)果

監(jiān)測點A?B的監(jiān)測結(jié)果見表1.

表1 監(jiān)測點A、B的溶解氧、TN、TP濃度(mg/L)Table 1 The concentration of DO、TN、TP at A、B points(mg/L)

2 湖泊富營養(yǎng)化耦合模型

2.1 模型建立

在二維水動力模型的基礎(chǔ)上,綜合考慮水中溶解氧狀態(tài)?營養(yǎng)物循環(huán)、浮游動植物、根系植被及大型藻類的生長活動等復(fù)雜的生物過程

[14-15],將影響湖泊富營養(yǎng)化的有機、無機營養(yǎng)物的輸移擴散過程組耦合到水質(zhì)模型中,形成富營養(yǎng)化模型,并與二維水動力模型進行耦合,最終形成淺水湖泊富營養(yǎng)化耦合模型.

耦合過程①:將營養(yǎng)物質(zhì)的輸移擴散過程融合到湖泊二維水質(zhì)模型中,形成富營養(yǎng)化模型,公式如下:

mg/L;h為水深m;Pc為富營養(yǎng)化過程組.

式中:ci為富營養(yǎng)化狀態(tài)變量濃度;n為狀態(tài)變量的過程數(shù);Pc變化率的單位有兩種:g/(m3?d)或mg/(L?d).

耦合過程②:將富營養(yǎng)化模型與湖泊二維水動力模型進行耦合,最終形成湖泊富營養(yǎng)化耦合模型,公式如下:

式中:f為科氏力系數(shù);Pa為大氣壓強;η為水面高度;us、vs為點源排入周圍水體的速度;ρ為水密度;ρ0為參考水密度;sχχ、sχy、syχ、syy為輻射應(yīng)力張量;為湖水表面 χ、y方向風(fēng)應(yīng)力;為湖水底部χ、y方向摩擦應(yīng)力為垂向平均速度;

水平應(yīng)力 Tij包括黏性摩擦、湍流摩擦、差動平流,是由建立在平均流速梯度基礎(chǔ)上的渦流黏度方程估算的:式中:A為水平渦黏系數(shù).

2.2 數(shù)值解法

本次計算利用MIKE21FM非結(jié)構(gòu)網(wǎng)格模型進行計算.非結(jié)構(gòu)網(wǎng)格模型[16]中采用的數(shù)值方法是單元中心的有限體積法.控制方程離散時,結(jié)果變量 U、V位于單元中心,跨邊界通量垂直于單元邊.有限體積法中法向通量通過在沿外法向建立單元水力模型并求解一維黎曼問題而得到.

2.3 水質(zhì)模擬邊界

本次選擇全湖水域進行模擬,利用遙感影像圖資料結(jié)合現(xiàn)場 GPS定位,確定湖岸點位坐標(biāo),并導(dǎo)入地理信息系統(tǒng)進行邊界數(shù)字化處理,最終形成迎澤湖水陸邊界.考慮到迎澤湖的面積?地形?水文特征等情況,采用三角網(wǎng)格系統(tǒng)將該湖概化為740個三角網(wǎng)格,最小允許角度為26°,輪廓由514個節(jié)點組成,垂直方向不設(shè)分層,迎澤湖模擬邊界見圖1.

2.4 模型率定及驗證

選取A、B兩點溶解氧濃度、TN濃度值、TP濃度值實測數(shù)據(jù)作為率定與驗證資料,經(jīng)過反復(fù)的模型參數(shù)調(diào)整,確定最終計算參數(shù),最終模型率定與驗證結(jié)果比較見圖2.

通過上述結(jié)果對比分析,可知不同變量模擬值與實測值相對誤差較小,變化規(guī)律基本一致,以DO為例,模擬出的溶解氧值接近實測值,兩者濃度值擬合程度較好.根據(jù)偏差統(tǒng)計法,其中最大絕對值偏差為0.83mg/L,相對誤差為10.36%,因此,應(yīng)用本文所構(gòu)建的富營養(yǎng)化耦合模型對迎澤湖進行水動力、富營養(yǎng)化模擬是有效的,進一步設(shè)定的生態(tài)補水方案也是可行的.

2.5 參數(shù)的確定

根據(jù)相關(guān)參考文獻[17]及監(jiān)測資料率定驗證結(jié)果選取參數(shù)如下:

(1) 水動力學(xué)參數(shù):

①糙率系數(shù).通過模型的率定和驗證來確定,并參照一定的經(jīng)驗,取值范圍為0.028～0.032.

圖2 模型率定與驗證結(jié)果對比Fig.2 The model calibration and validation results contrast

②渦黏系數(shù).根據(jù)Smagorinsky公式確定:式中:U、V為X、Y方向垂線平均流速;Δ為網(wǎng)格間距;Cs為計算參數(shù),一般選0.25＜Cs＜1.0.

③動邊界處理.為避免過強淺水效應(yīng),采用“干濕點判別法”來處理,即退水時水深<0.005m時,視為“干點”,不作為水域處理,水域流速為0;當(dāng)增水時水深>0.1m時,視為“濕點”,作為水域處理,流速取其計算值.

表2 主要敏感參數(shù)率定結(jié)果Table 2 Main sensitive parameter calibration results

④源匯項.分別將迎澤湖的進出水流作為源匯項,補水中總磷、總氮平均濃度分別為0.090、2.084mg/L.由于現(xiàn)有資料有限,不考慮降雨和蒸發(fā)對湖水量的影響.

(2) 水質(zhì)參數(shù):

富營養(yǎng)化耦合模型主要敏感參數(shù)率定結(jié)果見表2.

2.6 水動力學(xué)模擬

考慮到太原市主導(dǎo)風(fēng)向及頻率,湖泊主要動力為風(fēng)作用力,設(shè)置3種不同情景見表3,進行迎澤湖水動力模擬分析.限于篇幅,此處僅給出迎澤湖情景二下自由擴散狀態(tài)的典型流場例圖,見圖3.

表3 情景設(shè)置分類Table 3 The classifications of scenario

圖3 情景二補水期與非補水期流場示意Fig.3 Flow field with inflow and no inflow in scenario II

由圖 3分析可得出以下結(jié)論:在有風(fēng)情況下,補水期與非補水期間總體流場分布變化不大.由于補水期間,補水量與湖水總量相比很小,因此補水對湖流的形態(tài)影響遠小于風(fēng)場作用,風(fēng)生流為迎澤湖水流流場的主要流動形式;湖中流場分布不均勻,邊界區(qū)域流速較大,湖中心區(qū)域流速相對較小.北湖和南湖都明顯存在著環(huán)流區(qū),環(huán)流中心區(qū)域湖水流速緩慢.南湖北部存在較多狹窄區(qū),這些區(qū)域水流速幾乎為零,對流擴散過程緩慢,形成“死水區(qū)”;分析迎澤湖補水時期湖水流場可以看出,北湖入流補水時,由于南北湖通道狹窄,致使初期南湖湖水流動受進水水流影響不明顯,對流擴散作用很大程度還是依靠風(fēng)生流,這就造成入湖補水?dāng)U散至南湖歷時較長.

圖4 情景二Chla濃度分布Fig.4 Chla distribution in scenario II

3 湖泊富營養(yǎng)化耦合模型應(yīng)用

3.1 湖泊富營養(yǎng)化模擬

富營養(yǎng)化耦合模型中的富營養(yǎng)化模擬部分是在水動力模擬的基礎(chǔ)上進行的,因此富營養(yǎng)化模擬的數(shù)值模擬條件與所處風(fēng)場水動力模擬完全一致.根據(jù)國外湖泊富營養(yǎng)化分類設(shè)計的物理、化學(xué)和生物參數(shù)以及我國湖泊水質(zhì)參數(shù)與營養(yǎng)狀態(tài)之間關(guān)系,采用總氮?總磷、葉綠素a和透明度作為分析因子,對補水時段和非補水期自由擴散階段迎澤湖在各個主導(dǎo)風(fēng)場作用下湖泊富營養(yǎng)化耦合數(shù)值模擬進行研究.

由于所使用的是耦合模型,依據(jù)相應(yīng)的水動力情況,設(shè)置 3種不同的情景對迎澤湖富營養(yǎng)化狀態(tài)進行模擬分析,具體設(shè)置分類見表 3.限于篇幅,此處僅給出情景二下富營養(yǎng)化的表征因子葉綠素a濃度分布,見圖4.

由圖 4可以看出,低濃度的補水對湖泊水環(huán)境的改善起到了一定作用,補水進入湖區(qū)后,沿湖區(qū)邊界線向南擴散,中央?yún)^(qū)域的擴散較快,而岸邊的污染物擴散就相對較慢,導(dǎo)致岸邊水質(zhì)比中心水質(zhì)較差.在物質(zhì)擴散過程中,由于南北湖通道水動力不足,新鮮補水不能順利的到達全部湖區(qū),使南湖的營養(yǎng)物質(zhì)擴散輸移過程受到一定限制,水質(zhì)改善過程較為漫長,同時南湖結(jié)構(gòu)上存在“死角”和狹窄通道,湖流流速很低,形成“死水區(qū)”,營養(yǎng)物質(zhì)濃度高度集中,水質(zhì)較差.國際上一般以總氮濃度為0.2mg/L或總磷濃度0.02mg/L為湖泊富營養(yǎng)化發(fā)生濃度,迎澤湖中整湖區(qū)總氮?總磷值都遠遠超過此標(biāo)準(zhǔn),迎澤湖已達到富營養(yǎng)化發(fā)生所需要的營養(yǎng)狀態(tài).迎澤湖葉綠素 a的濃度都較高,變化范圍在0.035～0.105mg/L之間,北湖湖區(qū)葉綠素 a濃度偏低,南湖湖區(qū)相對偏高,分析其原因為北湖比較背光,水體溫度相對較低,而南湖隨著溫度的升高,光照的增強,藻類加速繁殖導(dǎo)致葉綠素a濃度增大,故南湖湖區(qū)葉綠素a濃度高于北湖. 3.2 水動力條件改進方案研究

由于迎澤湖水體已呈現(xiàn)中度富營養(yǎng)化狀態(tài)

[18],為保證迎澤湖的景觀效果,抑制藻類過度繁殖產(chǎn)生水華等不良現(xiàn)象,擬采用加大水力循環(huán)和對入水進行處理的方式保障水質(zhì).為考察方案的施行效果,本文對方案施行后分 3種情景進行了水質(zhì)模擬預(yù)測,方案列表見表4.

表4 改進方案Table 4 Improved proposals

圖5 新增補水口位置Fig.5 The new replenishment ports

圖6 連接通道布置Fig.6 The connected channels

圖7 方案一東東南風(fēng)作用下Chla濃度分布Fig.7 Chla distribution under ESE wind in scheme I

圖8 方案二東東南風(fēng)作用下Chla濃度分布Fig.8 Chla distribution under ESE wind in scheme II

限于篇幅,此處給出方案一的葉綠素 a濃度分布見圖 7,方案二葉綠素 a最不利情況濃度分布見圖8,方案三湖水流場分布見圖9.

對比圖4?圖7可得,增設(shè)1#?2#出水口后,新鮮補水可以在較短時間內(nèi)到達大部分湖區(qū),尤其是南湖區(qū)湖水更新速度加快,使這片區(qū)域水質(zhì)得到明顯改善,由于湖道狹長,水流速較快,浮游植物停留時間短,葉綠素a濃度值并不高.

圖9 方案三補水期與非補水期流場Fig.9 Flow field with inflow and no inflow in scenario III

對比圖 4?圖 8,由于大量的低濃度補水進入湖區(qū),一方面改善了湖區(qū)水動力條件,加快了營養(yǎng)物質(zhì)輸移擴散速度,另一方面,大量新鮮的補水縮短了湖水更新周期.在改進補水頻率后的最不利情況下,湖泊中各營養(yǎng)物質(zhì)濃度值都很低,透明度值較高,湖水水質(zhì)情況有明顯改善.

對比圖3、圖9可得,北湖與南湖連接通道流速加快,從而加快了北湖的新鮮補水向南湖的擴散速度.南湖中部及其西部“死角”湖區(qū)的水流速也明顯加快,水動力提升,從而避免了“死水區(qū)”的形成,為整個湖泊營養(yǎng)物質(zhì)輸移擴散提供了良好的動力環(huán)境.可見增加湖泊內(nèi)循環(huán)后,湖泊水動力環(huán)境明顯改善.

4 結(jié)論

4.1 綜合水動力模型、水質(zhì)模型和復(fù)雜的生物過程建立的富營養(yǎng)化耦合模型,可有效模擬迎澤湖的水流狀態(tài)及水質(zhì)演變情況,得出風(fēng)生流為迎澤湖水流流場的主要流動形式,為其富營養(yǎng)化規(guī)劃治理提供了理論參考和決策支持.

4.2 該類小尺度湖體的富營養(yǎng)化耦合數(shù)值模擬,參數(shù)的率定是難點,采用手動試錯法實現(xiàn)了參數(shù)的率定,但由于該方法具有主觀性,難以定量測量參數(shù)優(yōu)化水平,有待尋找更好的參數(shù)率定方法. 4.3 迎澤湖水體呈現(xiàn)中度富營養(yǎng)化狀態(tài),湖體中葉綠素a的濃度較高,濃度變化范圍在0.035～ 0.105mg/L之間,通過對不同方案富營養(yǎng)化模擬結(jié)果的分析,3種方案的實施,改變了迎澤湖的水動力條件,促進營養(yǎng)物質(zhì)的輸移擴散,改善了湖泊水質(zhì).另外,也可根據(jù)實際情況選擇多種方案同時進行,對湖泊水環(huán)境的改善會更為明顯.

[1] Vollenweider R A. Input-output models with special reference to the phosphorus loading concept in limnology [J]. Schweizerische Zeitschrift Hyrol, 1975,37:53-84.

[2] 蔡慶華,劉建康,Lorenz King.評價湖泊富營養(yǎng)化的一個綜合模型 [J]. 應(yīng)用生態(tài)學(xué)報, 2002,13(12):1674-1678.

[3] 韓 菲,陳永燦,劉昭偉.湖泊及水庫富營養(yǎng)化模型研究綜述 [J].水科學(xué)進展, 2003,14(6):785-791.

[4] 李錦秀,禹雪中,幸治國.三峽庫區(qū)支流富營養(yǎng)化模型開發(fā)研究[J]. 水科學(xué)進展, 2005,16(6):777-783.

[5] 王建平,蘇保林,賈海峰,等.密云水庫及其流域營養(yǎng)物集成模擬的模型體系研究 [J]. 環(huán)境科學(xué), 2006,27(7):1286-1291.

[6] 黃 莉,何 強,彭述娟.小型景觀湖水質(zhì)模擬及應(yīng)用研究 [J].環(huán)境工程學(xué)報, 2011,5(11):2453-2458.

[7] 姜 毅,張龍濤,呂永濤,等.基于 WASP模型城市人工景觀湖水

體水質(zhì)的模擬 [J]. 西安科技大學(xué)學(xué)報, 2011,31(3):306-310. [8] 余 成,任憲友,班 璇,等.二維水質(zhì)模型在武漢東湖引水工程中的應(yīng)用 [J]. 湖泊科學(xué), 2012,24(1):43-50.

[9] 韓龍喜,張防修,劉協(xié)亭,等.博斯騰湖湖流及礦化度分布研究[J]. 水力學(xué)報, 2004,10:100-105.

[10] 李 穎,曹文志,張玉珍,等.九龍江流域上游淺水湖泊富營養(yǎng)化機制 [J]. 中國環(huán)境科學(xué), 2012,32(5):906-911.

[11] 劉 永,郭懷成,范英英,等.湖泊生態(tài)系統(tǒng)動力學(xué)模型研究進展[J]. 應(yīng)用生態(tài)學(xué)報, 2005,16(6):1169-1175.

[12] 劉 坤,楊正宇.MIKE軟件在水體富營養(yǎng)化研究中的應(yīng)用 [J].給水排水, 2009,35:456-458.

[13] 國家環(huán)境保護總局.水和廢水監(jiān)測分析方法 [M]. 4版.北京:中國環(huán)境科學(xué)出版社, 2002.

[14] 何金平.金山湖水體運行方案及物質(zhì)輸運數(shù)值模擬研究 [D].南京:河海大學(xué), 2007.

[15] 李 璐.城市景觀湖富營養(yǎng)化模擬與研究 [D]. 太原:太原理工大學(xué), 2012.

[16] DHI Water & Environment.Mike21&MIKE3FLOW MODEL FM Mud Transport Module Scientific Documentation [M]. Denmark, 2009.

[17] 劉瑞民,王學(xué)軍.湖泊水質(zhì)參數(shù)空間優(yōu)化估算的原理與方法 [J].中國環(huán)境科學(xué), 2001,21(2):177-179.

[18] 徐明德,李 璐,白金玲,等.太原市迎澤湖富營養(yǎng)化灰色關(guān)聯(lián)法

分析 [J]. 太原理工大學(xué)學(xué)報, 2011,42(6):650-653.

Model study on eutrophication control in Yingze Lake of Taiyuan city.

WU Chun-fang, XU Ming-de*, LI Lu, JING

Sheng-yuan (College of Environmental Science and Engineering, Taiyuan University of Technology, Taiyuan 030024, China). China Environmental Science, 2014,34(2)：485～491

Aiming at Yingze Lake eutrophication, based on gathering data and long-time water quality monitoring, hydromatical model and water quality model, combined with a variety of biological processes were applied to established lake eutrophication coupled model, calibrated and verified by measurement of water quality data, the total nitrogen, total phosphorus, chlorophyll-a and transparency were selected to simulate, and concluded transport diffusion and spatial and temporal distribution of nutrients. Combined with the practical situation, from replenishment methods, replenishment frequency and changing lakes flexible structure, schemes on hydrodynamic condition and numerical simulation were put forward, and the modified effects of hydromatical condition and material transport diffusion of different schemes were studied. The results showed obvious spatial and temporal distribution of chlorophyll-a concentration which was between 0.035～0.105mg/L. Phytoplankton such as algae grow quickly,which was the peak period of algal bloom and water quality deterioration．Strategies such as increasing the hydraulic circulation and improving the inflow water quality would achieve good effects in water quality management.

Yingze Lake；hydrodynamics model；water quality model；eutrophication coupled model；MIKE 21

X524

：A

：1000-6923(2014)02-0485-07

武春芳(1987-),女,山西朔州人,太原理工大學(xué)碩士研究生,主要從事污染控制系統(tǒng)工程研究.

2013-06-20

* 責(zé)任作者, 教授, mingdexu@126.com