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

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

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

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

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

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

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

1 研究對(duì)象

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

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

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

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

1.3 監(jiān)測(cè)時(shí)限及頻率

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

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

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

監(jiān)測(cè)點(diǎn)A?B的監(jiān)測(cè)結(jié)果見(jiàn)表1.

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

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

2.1 模型建立

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

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

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

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

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

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

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

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

2.2 數(shù)值解法

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

2.3 水質(zhì)模擬邊界

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

2.4 模型率定及驗(yàn)證

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

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

2.5 參數(shù)的確定

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

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

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

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

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

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

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

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

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

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

2.6 水動(dòng)力學(xué)模擬

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

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

圖3 情景二補(bǔ)水期與非補(bǔ)水期流場(chǎng)示意Fig.3 Flow field with inflow and no inflow in scenario II

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

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

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

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

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

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

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

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

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

表4 改進(jìn)方案Table 4 Improved proposals

圖5 新增補(bǔ)水口位置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濃度分布見(jiàn)圖 7,方案二葉綠素 a最不利情況濃度分布見(jiàn)圖8,方案三湖水流場(chǎng)分布見(jiàn)圖9.

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

圖9 方案三補(bǔ)水期與非補(bǔ)水期流場(chǎng)Fig.9 Flow field with inflow and no inflow in scenario III

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

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

4 結(jié)論

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

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

[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.評(píng)價(jià)湖泊富營(yíng)養(yǎng)化的一個(gè)綜合模型 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2002,13(12):1674-1678.

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

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

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

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

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

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

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

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

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

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

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

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

[15] 李 璐.城市景觀湖富營(yíng)養(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].中國(guó)環(huán)境科學(xué), 2001,21(2):177-179.

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

分析 [J]. 太原理工大學(xué)學(xué)報(bào), 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