黎育紅，賀石磊

（華中科技大學(xué)水電與數(shù)字化工程學(xué)院，武漢 430074）

淺水湖泊群連通與調(diào)水的二維水動(dòng)力－水質(zhì)耦合模型研究

黎育紅，賀石磊

（華中科技大學(xué)水電與數(shù)字化工程學(xué)院，武漢 430074）

將湖泊群連通并進(jìn)行引清調(diào)水的水環(huán)境污染治理模式在淺水湖泊分布比較集中的地區(qū)逐漸受到人們的關(guān)注。為了探究獨(dú)立湖泊群相互連通引水后各湖區(qū)的水動(dòng)力水質(zhì)變化情況，以武漢大東湖生態(tài)水網(wǎng)工程的建設(shè)為背景，建立了湖泊群二維水動(dòng)力－水質(zhì)耦合模型。采用適合復(fù)雜邊界的非結(jié)構(gòu)化網(wǎng)格，考慮湖底地形和氣候條件，考慮污染物的輸入、遷移和轉(zhuǎn)化，以及藍(lán)藻等浮游生物的生長(zhǎng)條件，分別利用東湖2006年6月和2007年6月的野外數(shù)據(jù)對(duì)所建模型進(jìn)行參數(shù)的率定和校核，并利用2012年6月實(shí)測(cè)數(shù)據(jù)對(duì)3種引水方案與3種連通模式組合情況下湖泊群的BOD5，TP，TN以及Chl-a等生化指標(biāo)的變化情況進(jìn)行模擬，對(duì)穩(wěn)定運(yùn)行30 d后的模擬結(jié)果進(jìn)行比較分析。結(jié)果表明：引水工程中，對(duì)于水域面積較大的湖泊，風(fēng)力作用依舊是影響流場(chǎng)的主要因素；同樣的引水流量下，不同的調(diào)水方案對(duì)水質(zhì)改善結(jié)果有較大差別；湖泊群中隔堤的存在會(huì)給湖泊污染治理帶來(lái)困難，而在湖泊群中適當(dāng)區(qū)域設(shè)置生態(tài)小島對(duì)水質(zhì)變化過(guò)程影響微弱。

淺水湖泊群；水動(dòng)力－水質(zhì)耦合模型；調(diào)水方案；風(fēng)生流；情景規(guī)劃方法

2015，32（01）：21－27，38

1 研究背景

近幾十年來(lái)，隨著社會(huì)經(jīng)濟(jì)的發(fā)展和人口增長(zhǎng)，湖泊水體水質(zhì)惡化問(wèn)題逐漸凸顯［1－2］。為了研究湖泊水質(zhì)惡化的過(guò)程與機(jī)制，國(guó)內(nèi)外學(xué)者進(jìn)行了眾多的水動(dòng)力水質(zhì)模擬方法及相應(yīng)模型的研究，Hahram Missaghi等［3］運(yùn)用三維水動(dòng)力模型（EFDC）耦合生態(tài)模型（CAEDYM）對(duì)具有復(fù)雜地形的Minnetonka淺水湖進(jìn)行水質(zhì)參數(shù)的模擬研究，并對(duì)各水質(zhì)參數(shù)與生態(tài)變化之間的聯(lián)系進(jìn)行了探討。Changsheng Chen等［4］構(gòu)建了美國(guó)Michigan湖水動(dòng)力與生態(tài)模型，模擬結(jié)果很好地描述了Michigan湖藻類及浮游微生物隨季節(jié)變化的生長(zhǎng)繁殖狀況。Louise C.Bruce等［5］通過(guò)耦合的水動(dòng)力模型（DYRESM）和水生態(tài)系統(tǒng)動(dòng)力學(xué)模型（CAEDYM）對(duì)位于以色列的Kinneret湖進(jìn)行營(yíng)養(yǎng)鹽、浮游動(dòng)物和浮游植物的季節(jié)性變化模擬，結(jié)果很好地顯示了浮游動(dòng)物在Kinneret湖N，P等營(yíng)養(yǎng)鹽循環(huán)中的作用。國(guó)內(nèi)龔春生［6］等對(duì)淺水湖泊平面二維水流－水質(zhì)－底泥污染模型研究，其數(shù)值模擬的科學(xué)性得到普遍認(rèn)可。韓龍喜［7］等根據(jù)大型內(nèi)陸湖泊水動(dòng)力、氣象、水環(huán)境特征，建立了博斯騰湖二維水動(dòng)力、礦化度數(shù)學(xué)模型，預(yù)測(cè)分析了不同風(fēng)場(chǎng)、水文條件、水資源調(diào)度方案對(duì)湖流、礦化度空間分布的影響。萬(wàn)金保［8］等詳細(xì)介紹了湖泊水質(zhì)模型國(guó)內(nèi)外研究動(dòng)態(tài)、類型、常用軟件和應(yīng)用實(shí)例，并綜觀湖泊水質(zhì)模型的研究歷史和應(yīng)用前景，系統(tǒng)分析了湖泊水質(zhì)模型研究的發(fā)展趨勢(shì)等等。但是通過(guò)總結(jié)可以發(fā)現(xiàn)，這些模型和方法基本上都是針對(duì)單個(gè)湖泊，而以連通的淺水湖泊群為對(duì)象的模型研究并不多見(jiàn)。對(duì)于湖泊分布比較集中的區(qū)域，湖泊間連通能夠彼此建立水力聯(lián)系，有利于污染物的輸移和轉(zhuǎn)化，實(shí)施調(diào)水沖污過(guò)程中則能改善整個(gè)湖泊群的水質(zhì)。不過(guò)，湖泊連通和引水對(duì)于湖泊群水質(zhì)改善的具體效果，還需要進(jìn)行科學(xué)的模型分析。因此建立對(duì)象為淺水湖泊群的連通與調(diào)水二維水動(dòng)力－水質(zhì)耦合模型，對(duì)于淺水湖泊群水質(zhì)水生態(tài)研究以及實(shí)際湖泊連通引水工程的設(shè)計(jì)，都具有十分重要的意義。

以武漢東湖為主的大東湖湖泊群落是國(guó)內(nèi)最大的城市淺水湖泊群，近年來(lái)由于人類活動(dòng)的強(qiáng)力干涉，使得各子湖水質(zhì)都存在不同程度的惡化［9］。為改善湖泊群水環(huán)境，武漢市政府正在實(shí)施“大東湖生態(tài)水網(wǎng)構(gòu)建工程”，即通過(guò)污染控制、水網(wǎng)連通和生態(tài)修復(fù)3大工程對(duì)湖泊群進(jìn)行綜合治理［10］。已經(jīng)實(shí)施的截污工程初見(jiàn)成效，2006年?yáng)|湖水質(zhì)實(shí)現(xiàn)從1990年以來(lái)的首次好轉(zhuǎn)，但依舊沒(méi)有達(dá)到水功能區(qū)劃要求。基于截污工程實(shí)施引起入湖水量減少以及湖水更新的實(shí)際需要考慮，應(yīng)向湖泊進(jìn)行生態(tài)補(bǔ)水［11］。生態(tài)水網(wǎng)工程充分利用東湖毗鄰長(zhǎng)江便于引水的地理優(yōu)勢(shì)，將沙湖、東湖、楊春湖、嚴(yán)西湖、嚴(yán)東湖、北湖6個(gè)湖泊貫通并與長(zhǎng)江相連，構(gòu)建江湖連通的生態(tài)水網(wǎng)［12］。期望扭轉(zhuǎn)大東湖區(qū)水污染及水生態(tài)持續(xù)惡化的趨勢(shì)，改善水環(huán)境，豐富生物多樣性，加快實(shí)現(xiàn)湖區(qū)的生態(tài)平衡［13］。

本文以武漢“大東湖生態(tài)水網(wǎng)工程”實(shí)施為工程應(yīng)用背景，借助對(duì)前人有關(guān)淺水湖泊群水環(huán)境模型研究的理論和經(jīng)驗(yàn)的分析，建立東湖水網(wǎng)二維水動(dòng)力－水質(zhì)耦合模型，旨在探究淺水湖泊群連通和調(diào)水對(duì)于湖泊水動(dòng)力水質(zhì)影響的大小和規(guī)律，為湖泊污染治理及具體工程的實(shí)施提供科學(xué)的參考和依據(jù)。模型搭建于DHI水動(dòng)力學(xué)模型MIKE21之上，綜合考慮風(fēng)場(chǎng)、氣溫、湖底地形、入湖水量水質(zhì)、物質(zhì)遷移轉(zhuǎn)化等因素，耦合水動(dòng)力模塊（Hydrodynamic）和富營(yíng)養(yǎng)化生態(tài)實(shí)驗(yàn)室模塊（ECO Lab Module），對(duì)湖泊群流場(chǎng)及BOD5，TP，TN，Chl-a等水質(zhì)、水生態(tài)指標(biāo)進(jìn)行模擬。為了探究不同引水和連通方案對(duì)N，P等營(yíng)養(yǎng)鹽濃度場(chǎng)分布的影響，以及由于水環(huán)境變化引起的藻類遷移和數(shù)量變化，文中應(yīng)用情景規(guī)劃方法將3種引水模式和3種連通方案組合成9種模擬場(chǎng)景。對(duì)不同引水連通方案下湖泊中污染物的時(shí)空變化結(jié)果進(jìn)行分析，得到各子湖區(qū)具體的水質(zhì)改善情況以及最佳的連通調(diào)度方案，為淺水湖泊群水環(huán)境恢復(fù)與重建模型的探究提供新的思路，也為大東湖生態(tài)水網(wǎng)連通工程的實(shí)施提供科學(xué)的決策與參考。

2 研究區(qū)概況

為了使模型更加符合實(shí)際情況，模型主要以已經(jīng)實(shí)現(xiàn)連通的東湖和沙湖為研究區(qū)域。東湖（114° 09′～114°39′E，30°22′～30°41′N(xiāo)）位于武漢市武昌區(qū)東北部，平均水深約2.21 m，最大水深4.5 m，水域面積約為27.899 km2，為典型的城市淺水湖泊。沙湖（114°18′46′～114°20′49″E，30°33′36″～30°34′43″N）位于武漢市武昌東北部，東鄰中北路，南至小龜山，西抵武昌至大冶的鐵路線，北達(dá)徐東路。清末修筑的粵漢鐵路穿湖而過(guò)，路西為小沙湖，又名內(nèi)沙湖，現(xiàn)已近乎湮沒(méi)；路東為大沙湖，又名外沙湖，即現(xiàn)在的沙湖。按武漢市水務(wù)局最新公布的測(cè)量數(shù)據(jù)，內(nèi)沙湖現(xiàn)實(shí)有面積0.134 km2，外沙湖現(xiàn)實(shí)有面積3.197 km2，是武漢市僅次于東湖的第二大“城中湖”，也是武漢市區(qū)內(nèi)環(huán)線內(nèi)唯一的湖泊。

圖1 研究區(qū)域Fig.1 Study area

東湖－沙湖區(qū)是武漢市水源污染最嚴(yán)重的地區(qū)之一，近幾十年來(lái)，由于湖區(qū)周?chē)丝诿芏燃ぴ觯由瞎I(yè)、養(yǎng)殖以及旅游業(yè)的迅猛發(fā)展，大量生活和工業(yè)污水排入湖區(qū)。外源污染源源不斷，湖內(nèi)污染日益加重，導(dǎo)致浮游植物在數(shù)量上大幅增加，群落結(jié)構(gòu)上發(fā)生明顯變化，以致水體濁度增大，透明度降低，水質(zhì)下降，富營(yíng)養(yǎng)化日趨嚴(yán)重［14］。針對(duì)東湖的水生態(tài)環(huán)境，武漢市提出了“大東湖生態(tài)水網(wǎng)工程”，首期工程便是將東湖和沙湖相連，從長(zhǎng)江引水改善沙湖和東湖的水質(zhì)和生態(tài)，為下期工程的開(kāi)展提供經(jīng)驗(yàn)和參考。

3 模型原理

研究湖區(qū)平均水深低于3 m，屬于典型的淺水湖泊，垂直方向不存在明顯的分層現(xiàn)象，又考慮到東湖邊界的復(fù)雜性，因此選擇在MIKE21模型之上搭建二維數(shù)值模型。MIKE21是丹麥水環(huán)境研究所（DHI）開(kāi)發(fā)的一個(gè)具有非結(jié)構(gòu)化網(wǎng)格的二維水質(zhì)模擬軟件，能夠擬合復(fù)雜的邊界條件［15］。本文選擇了MIKE21的水動(dòng)力模塊（Hydrodynamic）和富營(yíng)養(yǎng)化生態(tài)實(shí)驗(yàn)室模塊（ECO Lab Module）。水動(dòng)力模塊能夠模擬水體水位和流場(chǎng)變化，是其他模塊的基礎(chǔ)。富營(yíng)養(yǎng)化生態(tài)實(shí)驗(yàn)室模塊能模擬不同外界條件下湖泊中營(yíng)養(yǎng)物質(zhì)的擴(kuò)散、遷移及轉(zhuǎn)化過(guò)程，同時(shí)基于實(shí)時(shí)的水動(dòng)力和水質(zhì)條件，模擬Chl-a等富營(yíng)養(yǎng)化因子的濃度及分布變化。

3.1 水動(dòng)力方程

水動(dòng)力模塊的控制方程是沿水深平均的二維淺水流動(dòng)質(zhì)量和動(dòng)量守恒控制方程組，其連續(xù)性方程和動(dòng)量方程可分別表示為

式中：ε為自由水面水位（m）；h為水深（m）；u，v分別為x，y方向上的垂線平均流速（m/s）；p＝hu，q＝hv分別為單寬流量在x，y方向上的分量［m·（m/s）］；g為重力加速度（m/s2）；f為阻力系數(shù)，f＝f1＋f2，f1為河床底部摩阻系數(shù)，f1＝g，C為反映河床底部對(duì)水流阻礙作用的謝才系數(shù)，（m1/2/s；f2為局部摩阻系數(shù)，f2＝，ξ為局部阻力系數(shù)；Ω＝ 2ωsinφ為科氏力系數(shù)，反映了地球自轉(zhuǎn)偏心力的作用，ω為地球自轉(zhuǎn)角速度，為2π/（24×3 600）rad/s，φ為所在點(diǎn)的緯度；v為紊流渦黏性系數(shù)；λ為風(fēng)應(yīng)力系數(shù)；ρa(bǔ)為空氣密度；ρw為水體密度；w，ua，va分別為風(fēng)速及其在x，y方向的分量（m/s）。

3.2 污染物擴(kuò)散方程

可溶污染物的擴(kuò)散控制方程式為

式中：h為水深（m）；c為污染物濃度；u，v分別為x，y方向上的水平流速（m/s）；F為線性衰減系數(shù)（s－1）；λx，λy為x、y方向上的擴(kuò)散系數(shù)（m2/s）；S＝Qs·（cs－c），其中Qs為源匯項(xiàng)流量［m3/（s·m2）］，cs為源匯項(xiàng)處污染物相對(duì)濃度。

4 二維水動(dòng)力－水質(zhì)耦合模型構(gòu)建

4.1 模型區(qū)域和地形

東湖地形邊界復(fù)雜，且湖區(qū)之間存在較多隔堤，綜合考慮模型的精確度和適用性，在保留主要隔堤的前提下，根據(jù)東湖－沙湖區(qū)的地形矢量數(shù)據(jù)生成模型邊界，然后利用MIKE21模型的網(wǎng)格生成器生成內(nèi)部非結(jié)構(gòu)化網(wǎng)格，模型網(wǎng)格由3 231個(gè)三角形組成，共計(jì)節(jié)點(diǎn)2 123個(gè)（圖2a）。最后利用2012年6月采集的水深數(shù)據(jù)進(jìn)行插值，得到模型的地形文件（圖2b）。

圖2 東湖－沙湖地形網(wǎng)格和水下地形Fig.2 Topography grids and terrain of East Lake and Sand Lake

4.2 初始及邊界條件的設(shè)定

4.2.1 區(qū)域劃分及初始條件設(shè)定

為了研究方便，將東湖－沙湖區(qū)按地理位置和污染狀況劃分為6塊子區(qū)域：沙湖區(qū)（包括沙湖、水果湖及郭鄭湖西灣）、郭鄭湖區(qū)（包括郭鄭湖、筲箕湖以及菱角湖）、湯菱湖區(qū)（包括湯菱湖和小潭湖）、團(tuán)湖區(qū)（包括團(tuán)湖和后湖）、官橋湖區(qū)（官橋湖）以及喻家湖區(qū)（喻家湖）。每個(gè)子區(qū)域設(shè)置3個(gè)測(cè)點(diǎn)（喻家湖除外），以3測(cè)點(diǎn)的污染物平均實(shí)測(cè)濃度值作為子區(qū)域的初始濃度值，由于喻家湖面積較小，僅設(shè)1個(gè)測(cè)點(diǎn)。模型采用武漢市6月多年平均風(fēng)速3.5 m/s，風(fēng)向105°；溫度為各模擬時(shí)段實(shí)測(cè)值；鹽度95‰；源匯項(xiàng)共設(shè)2個(gè)引水口（青山港、曾家港），3個(gè)排水口（新溝渠、九峰渠以及沙湖港），7個(gè)污染排放口（圖3）。

圖3 區(qū)域監(jiān)測(cè)點(diǎn)及排污口位置Fig.3 Location of sam p ling sites and drainage outlets

4.2.2 邊界條件及調(diào)水方案設(shè)置

東湖－沙湖的污染主要來(lái)自于點(diǎn)源污染、面源污染以及底泥釋放，其中點(diǎn)源占絕對(duì)的主導(dǎo)地位，所以模型中主要考慮點(diǎn)源污染排放量。東湖－沙湖周邊的排污口眾多，無(wú)法一一統(tǒng)計(jì)，模型中把污染口進(jìn)行整合，最后整合成7大排污口，如表1所示。

表1 排污口位置及排放量Table 1 Location and discharge of drainage outlets 103t/d

2006年?yáng)|湖－沙湖減排后日排約16萬(wàn)t污水，相當(dāng)于流量1.85 m3/s。按照甘義群《武漢東湖富營(yíng)養(yǎng)化現(xiàn)狀分析及治理對(duì)策》［16］中各污染物質(zhì)的總量折算，污水水質(zhì)大致如下：BOD 34.22 mg/L，TN 7.63 mg/L，TP 0.752 mg/L，NH3-N 7.58 mg/L，Chl-a 0.003 mg/L，DO 6 mg/L。

根據(jù)湖北省武漢市大東湖生態(tài)水網(wǎng)構(gòu)建水網(wǎng)連通工程（近期）可行性研究報(bào)告，確定模型調(diào)水方案如表2。長(zhǎng)江水質(zhì)狀況為：BOD 3.21 mg/L，TN 0.72 mg/L，TP 0.05 mg/L，NH3-N 0.67 mg/L，Chl-a 0.006 7 mg/L，DO 7.6 mg/L。

表2 引水調(diào)度模式Table 2 Scenarios of water diversion into the lake group 103m/s

4.3 模型的校核與驗(yàn)證

依據(jù)東湖－沙湖未連通引水前的實(shí)際情形，利用2006年6月16個(gè)監(jiān)測(cè)點(diǎn)的實(shí)測(cè)數(shù)據(jù)對(duì)模型率定，模型率定結(jié)果為床底摩擦系數(shù)（Manning number）為50，渦黏系數(shù)（Smagorinsky formulation）為0.31。再利用2007年6月的相關(guān)數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證，驗(yàn)證結(jié)果如圖4所示水深、TP、TN及BOD的平均誤差分別為7.16%，19.59%，25.51%，29.50%。

5 結(jié)果與討論

5.1 方案設(shè)定

在東湖生態(tài)水網(wǎng)構(gòu)建中，東湖中的隔堤對(duì)于東湖水動(dòng)力及水質(zhì)過(guò)程影響較大。所以，本研究在東湖－沙湖連通水動(dòng)力水質(zhì)模型的基礎(chǔ)上，考慮保留隔堤、拆除隔堤以及拆除隔堤將土方堆成湖心島（減少運(yùn)輸費(fèi)用，同時(shí)湖心島可用于旅游開(kāi)發(fā)）3種情況下各湖區(qū)的水質(zhì)改善情況，然后與上文中提出的3種引水調(diào)度模式相結(jié)合得到9種調(diào)水－隔堤方案，如表3所示。

圖4 模型水深，TP，TN和BOD率定及驗(yàn)證結(jié)果Fig.4 Calibration and validation results of water depth，TP，TN and BOD of themodel

表3 研究方案Table 3 Research scenarios

5.2 模擬結(jié)果

在保留隔堤的情況下，以校驗(yàn)?zāi)Ｐ椭蟹€(wěn)定后污染物濃度分布作為初值（東湖－沙湖未連通引水），3種引水模式連續(xù)引水30 d后湖區(qū)流場(chǎng)與污染物分布狀況分別如圖5、圖6所示。由圖5可知，3種調(diào)水模式下在研究區(qū)域的主湖區(qū)郭鄭湖都形成了2個(gè)較大的環(huán)流，靠北的環(huán)流呈逆時(shí)針?lè)较颍磕系沫h(huán)流呈順時(shí)間方向，團(tuán)湖區(qū)則在南部形成1個(gè)順時(shí)針環(huán)流，說(shuō)明風(fēng)力作用是影響東湖－沙湖流場(chǎng)的主要因素。在相對(duì)封閉、水域面較小且有引水注入的湖區(qū)，比如沙湖、小潭湖、筲箕湖等則主要受引水流量的影響，因此這些區(qū)域在不同的調(diào)水方案下呈現(xiàn)出各自獨(dú)有的流場(chǎng)分布情況。湯菱湖區(qū)在3種方案下流速和流向都有差別，說(shuō)明其受引水流量和風(fēng)力作用共同影響。官橋湖和喻家湖由于隔堤阻隔，且面積小，受引水流量和風(fēng)力作用的影響都較微弱，流速緩慢，與主湖區(qū)水力交換十分平緩。

圖5 保留隔堤條件下3種調(diào)水方案流場(chǎng)分布Fig.5 Flow field in the presence of dividing dike in three diversion schemes

圖6 調(diào)水前后的TP濃度空間分布Fig.6 Spatial distribution of TP concentration before and after water diversion

東湖的主要污染因子為T(mén)P，由于篇幅所限，文中主要以TP為代表對(duì)各方案進(jìn)行比較分析。如圖6所示，除了官橋湖和喻家湖外，3種引水模式都能夠有效地降低TP濃度，但是不同的引水方式對(duì)各子湖區(qū)的水質(zhì)改善大小稍有差別。方案3對(duì)沙湖區(qū)TP濃度改善效果比方案1好，但在郭鄭湖區(qū)情況則恰好相反，余下湖區(qū)改善效果較為接近。方案2在郭鄭湖區(qū)TP改善效果最好，但是在其他區(qū)域效果不如方案1，3顯著。由上文中流場(chǎng)分布情況可知，官橋湖和喻家湖與主湖區(qū)水力交換緩慢，同時(shí)自身還伴有穩(wěn)定的污水排入，3種調(diào)水方案下TP濃度變化幅度都很小，污染物濃度始終維持在較高水平，成為東湖最難治理的區(qū)域。通過(guò)對(duì)比分析可知，方案1在總體上對(duì)TP的改善效果最佳，分布也最均勻。

在大東湖生態(tài)水網(wǎng)工程長(zhǎng)遠(yuǎn)規(guī)劃中，將逐步對(duì)東湖中的隔堤進(jìn)行拆除。文中考慮東湖主要隔堤拆除后（圖7a）3種調(diào)水模式下各湖區(qū)的水質(zhì)改善情況，并進(jìn)一步創(chuàng)造性地在拆除隔堤后模型中設(shè)置虛擬小島模型（圖7b），小島假設(shè)由隔堤土方堆成，旨在解決土方輸運(yùn)及堆積方面難題。3種隔堤存在形式及3種調(diào)水方案構(gòu)成9類組合，各組合方案下各湖區(qū)的水質(zhì)改善效果如表4。

由表4中統(tǒng)計(jì)數(shù)據(jù)可知，保留隔堤時(shí)，方案1即第1種調(diào)水模式下水質(zhì)改善效果最佳，與上文中分析結(jié)果一致。拆除隔堤后，無(wú)論運(yùn)用哪種調(diào)水模式各湖區(qū)水質(zhì)改善效果都有明顯提高，而3種調(diào)水模式中依舊是第1種調(diào)水模式呈現(xiàn)最佳的引水效果?？臻g上比較，隔堤拆除后，官橋湖和喻家湖的改善效果最為明顯，這是由于隔堤拆除后改變了兩處的水動(dòng)力條件，使其與主湖區(qū)的水力聯(lián)系更加頻繁。假設(shè)將拆除隔堤的土方就近堆放于東湖偏北位置，形成湖心島，模擬結(jié)果顯示該措施對(duì)水質(zhì)改善過(guò)程基本上沒(méi)有影響，所以在實(shí)際工程中，可以考慮通過(guò)此措施來(lái)降低土方運(yùn)輸成本。

圖7 東湖－沙湖隔堤拆除和堆成小島地形網(wǎng)格Fig.7 Topography grids of East Lake and Sand Lake when dividing dike is rem oved and when earthwork is piled up into small island

表4 9種組合方案水質(zhì)改善效果Table 4 Improvements of lake water quality in the presence of different combinatorial schemes

6 結(jié) 論

（1）在引水過(guò)程中，除了局部區(qū)域主要受引水形式影響以外，風(fēng)力作用依舊是影響東湖－沙湖流場(chǎng)的主要因素?，F(xiàn)今（存在隔堤）條件下，官橋湖和喻家湖受湖流影響很小，所以在引水工程中受益也最小。

（2）同一種隔堤存在形式下，以第1種調(diào)水模式引水，即青山港引水流量40 m3/s，九峰渠出流量20 m3/s，新溝渠、沙湖港出流量10m3/s時(shí)水質(zhì)改善效果最佳。

（3）與拆除隔堤前比較，隔堤拆除后能夠明顯地提高各湖區(qū)的水質(zhì)改善效果，在官橋湖和喻家湖這2個(gè)最難治理的區(qū)域表現(xiàn)得尤為明顯。在東湖偏北位置添加小島對(duì)水質(zhì)變化過(guò)程影響不大。

［1］CHAU KW.A Three-dimensional Eutrophication Modeling in Tolo Harbour［J］.Applied Mathematical Modelling，2004，28（9）：849－861.

［2］萬(wàn) 蕾，我國(guó)湖泊富營(yíng)養(yǎng)化問(wèn)題與治理現(xiàn)狀［J］.生態(tài)環(huán)境，2012，（1）：378－381，389.（WAN Lei.Present Situation of Lake Eutrophication and Control in China［J］.Ecological Economy，2012，（1）：378－381，389.（in Chinese））

［3］MISSAGHIS，HONDZO M.Evaluation and Application of a Three-dimensionalWater Quality Model in a Shallow Lake with Complex Morphometry［J］.Ecological Modelling，2010，221（11）：1512－1525.

［4］CHEN C，JIR，SCHWAB D J，et al.A Model Study of the Coupled Biological and Physical Dynamics in Lake Michigan［J］.Ecological Modelling，2002，152（2/3）：145－168.

［5］BRUCE L C，HAMILTON D，IMBERGER J，et al.A Numerical Simulation of the Role of Zooplankton in C，N and P Cycling in Lake Kinneret，Israel［J］.Ecological Modelling，2006，193（3/4）：412－436.

［6］龔春生，姚 琪，趙棣華.淺水湖泊平面二維水流-水質(zhì)-底泥污染模型研究［J］.水科學(xué)進(jìn)展，2006，17（4）：496－501.（GONG Chun-sheng，YAO Qi，ZHAO Dihua.Plane 2D Flow-water Quality-Sediment Pollution Couple Model in Shallow Lake［J］.Advances in Water Science，2006，17（4）：496－501.（in Chinese））

［7］韓龍喜，張防修，劉協(xié)亭，等.博斯騰湖湖流及礦化度分布研究［J］.水利學(xué)報(bào)，2004，10（1）：100－105.（HAN Long-xi，ZHANG Fang-xiu，LIU Xie-ting，et al.Flow Field and Salinity Distribution of Large Inland Lake［J］.Journal of Hydraulic Engineering，2004，10（1）：100－105.（in Chinese））

［8］萬(wàn)金保，李媛媛.湖泊水質(zhì)模型研究進(jìn)展［J］.長(zhǎng)江流域資源與環(huán)境，2007，16（6）：805－809.（WAN Jinbao，LI Yuan-yuan.Progress in the Research of Lake Water Quality Model［J］.Resources and Environment in the Yangtze Basin，2007，16（6）：805－809.（in Chinese））

［9］XIE P，LIU J.Practical Success of Biomanipulation U-sing Filter-feeding Fish to Control Cyanobacteria Bloom［J］.The Scientific World，2001，（1）：337－356.

［10］張軍花，項(xiàng)久華，邱漢明.大東湖生態(tài)水網(wǎng)構(gòu)建工程總體構(gòu)架設(shè)計(jì)［J］.城市道橋與防洪，2008，（10）：24－27.（ZHANG Jun-ha，XIANG Jiu-hua，QIU Han-ming.Design of Overall Frame of the Lake East Ecology-oriented Water Network Project［J］.Urban Roads Bridges and Flood Control，2008，（10）：24－27.（in Chinese））

［11］康 玲，郭曉明，王學(xué)立.大型城市湖泊群引水調(diào)度模式研究［J］.水利發(fā)電學(xué)報(bào)，2012，31（3）：65－69.（KANG Ling，GUO Xiao-ming，WANGXue-li.Study on Water Diversion Schemes of Large Urban Lake Group［J］.Journal of Hydroelectric Engineering，2012，31（3）：65－69.（in Chinese））

［12］周 俊，王焰新，蔡鶴生.構(gòu)建“長(zhǎng)江-東湖-長(zhǎng)江”水循環(huán)及其在東湖污染中的意義［J］.武漢大學(xué)學(xué)報(bào)（工學(xué)版），2002，35（4）：64－67.（ZHOU Jun，WANG Yanxin，CAIHe-sheng.Construction of Water Exchange between Yangtze River and East Lake and Its Significance to Pollution Control of East Lake［J］.Journal of Wuhan University of Hydraulic and Electric Engineering，2002，35（4）：64－67.（in Chinese））

［13］嚴(yán)江涌，黎南關(guān).武漢市大東湖水網(wǎng)連通治理工程淺析［J］.人民長(zhǎng)江，2010，41（11）：82－84.（YAN Jiang-yong，LI Nan-guan.Water Network Connection Project of Donghu Lake in Wuhan City［J］.Yangtze River，2010，41（11）：82－84.（in Chinese））

［14］余 成，任憲友，班 璇，等.二維水質(zhì)模型在武漢東湖引水工程中的應(yīng)用［J］.湖泊科學(xué)，2012，24（1）：43－50.（YU Cheng，REN Xian-you，BAN Xuan，et al.Application of Two-dimensionalWater Quality Model in the Project of the Water Diversion in East Lake，Wuhan［J］.Journal of Lake Sciences，2012，24（1）：43－50.（in Chinese））

［15］Danish Hydraulic Institute（DHI）.MIKE 21 Flow Model：Hydrodynamic Module Scientific Documentation［M］.Denmark：DHI，2007.

［16］甘義群，郭永龍.武漢東湖富營(yíng)養(yǎng)化現(xiàn)狀分析及治理對(duì)策［J］.長(zhǎng)江流域資源與環(huán)境，2004，13（3）：277－281.（GAN Yi-qun，GUO Yong-long.Evaluation Analysis Strategy for Eutrophication in Wuhan Lake Donghu［J］.Resources and Environment in the Yangtze Basin，2004，13（3）：277－281.（in Chinese） ）

（編輯：趙衛(wèi)兵）

A 2-D M odel Coupling Hydrodynam ics and W ater Quality for the W ater Connection and Diversion of Shallow Lake Group

LIYu-hong，HE Shi-lei
（School of Hydropower and Information Engineering，Huazhong University of Science and Technology，Wuhan 430074，China）

Themethod of treatingwater pollution by reestablishing hydraulic connection between lakes and diverting water from a nearby river to the lake group which is seriously polluted is gettingwider attention in shallow lake areas.To investigate the variation of hydrodynamics and water quality after lake connection and diversion，we built a 2-D model coupling hydrodynamics and water quality.The Projectof East Lake Network in Wuhan was taken as research background.In thismodel，unstructured grid was adopted to fit the complex boundary，and lake topography and climate conditionswere adequately reflected.In addition，the input，transport and transformation of nutrients and pollutants aswell as the growing conditions of some algae such as cyanobacteria were taken into account.Moreover，field measurement data in June 2006 and June 2007 were used for the calibration and validation of themodel

parameters.Then the variation of indexes including BOD5，TP，TN and Chl-a in three diversion conditions and three connection conditions after 30 days of stable operation were simulated based on the measured data in June 2012.Results suggest thatwind force is still themajor factor which induces the change of lake flow field for larger lakes in awater diversion project；in the presence of the same diversion flow，different diversion conditions has different effects on thewater quality improvement.Dividing dikes in shallow lake group bring difficulties to the control of lake pollution，whereas setting up artificial island at suitable position of the shallow lake has no great impact on the water quality variation.

shallow lake group；hydrodynamic-water quality coupling model；scheme of water diversion；winddriven current；Scenario Planning Method

X524

A

1001－5485（2015）01－0021－07

10.3969/j.issn.1001－5485.2015.01.005

2013－09－12；

2013－11－22

黎育紅（1972－），男，湖北孝感人，副教授，博士，研究方向?yàn)閺?fù)雜系統(tǒng)分析及其工程應(yīng)用，（電話）13886101077（電子信箱）liyuhong＠hust.edu.cn。