黃春琳,李 熙,孫永遠(yuǎn)

(1:河海大學(xué)海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室，南京 210098)(2：江蘇省水文水資源勘測(cè)局，南京 210029)

太湖水齡分布特征及“引江濟(jì)太”工程對(duì)其的影響

黃春琳1,李 熙1,孫永遠(yuǎn)2

(1:河海大學(xué)海岸災(zāi)害及防護(hù)教育部重點(diǎn)實(shí)驗(yàn)室，南京 210098)(2：江蘇省水文水資源勘測(cè)局，南京 210029)

太湖作為典型的風(fēng)生流湖泊，風(fēng)場(chǎng)對(duì)水體運(yùn)輸和交換的過(guò)程起著重要的作用. 基于環(huán)境水動(dòng)力學(xué)模型EFDC源程序建立了染色劑模型以及水齡模型，借助水齡研究太湖水體的長(zhǎng)期輸運(yùn)過(guò)程和更新速率特征，從而為太湖的調(diào)水工程管理提供科學(xué)依據(jù). 本文在研究春、夏、秋、冬季不同風(fēng)場(chǎng)作用下太湖的水齡季節(jié)性分布特征的基礎(chǔ)上，結(jié)合太湖實(shí)測(cè)風(fēng)速、流量數(shù)據(jù)及“引江濟(jì)太”工程調(diào)水運(yùn)行的現(xiàn)狀，著重分析了望虞河樞紐調(diào)水運(yùn)行以及新溝河工程對(duì)太湖水齡分布的影響. 數(shù)值試驗(yàn)的模擬結(jié)果表明：太湖水體交換受季風(fēng)影響明顯，春、夏季的水齡相對(duì)較大、水體交換較差；靠近湖流入口處的地方水齡較小，遠(yuǎn)離入口的地方水齡較大；水流流向與風(fēng)向一致時(shí)水齡減小，水循環(huán)加快，反之則減慢. 望虞河引水引工程能夠減小貢湖及湖心區(qū)的水齡，加快貢湖湖區(qū)及湖心區(qū)的水循環(huán)；新溝河工程引水能夠減小梅梁湖區(qū)的水齡，改善該湖區(qū)的水質(zhì). 引水工程的實(shí)施對(duì)加快整個(gè)太湖的水循環(huán)做出重要貢獻(xiàn).

太湖；風(fēng)場(chǎng)；染色劑；水齡；EFDC模型

入湖流量、風(fēng)場(chǎng)等在湖體的水體運(yùn)動(dòng)過(guò)程中起很重要的作用，三維環(huán)境流體動(dòng)力學(xué)模型EFDC提供了一個(gè)高分辨率的工具來(lái)研究太湖的水體交換. 水齡可以定量反映水體的運(yùn)動(dòng)和交換程度以及滯留情況，對(duì)評(píng)估水質(zhì)的變化具有重要意義. Deleersnijder等[1]介紹了基于對(duì)流-擴(kuò)散輸運(yùn)方程的水齡理論，這個(gè)新的水齡定義提供了一個(gè)健全的方法來(lái)定量地研究歐拉坐標(biāo)框架下的物質(zhì)溶解的輸運(yùn)特征. Li等[2]運(yùn)用EFDC三維水動(dòng)力模型，借助水齡的定義，研究了“引江濟(jì)太”引水量對(duì)太湖水循環(huán)的影響.

太湖作為中國(guó)的第三大淡水湖，由于富營(yíng)養(yǎng)化而嚴(yán)重影響了飲用水的安全，政府部門(mén)近年來(lái)采取了許多措施來(lái)提升太湖水質(zhì)，包括建更多的污水處理廠，禁止使用含磷清潔劑等. 此外，一個(gè)重要的生態(tài)工程——“引江濟(jì)太”工程在2002年啟動(dòng)，“引江濟(jì)太”工程將更優(yōu)質(zhì)的長(zhǎng)江水通過(guò)望虞河引入太湖，最終從太浦河流出，以期通過(guò)此舉來(lái)加快太湖的水循環(huán)，改善太湖的水質(zhì)狀況[3-5]. 此外，在2014年太湖北部流域內(nèi)開(kāi)始疏通新溝河、新孟河兩條河道，太湖將通過(guò)新孟河、望虞河補(bǔ)水，借助走馬塘、新溝河、太浦河向長(zhǎng)江排水，形成“兩進(jìn)三出”的循環(huán)系統(tǒng)，通過(guò)讓水體動(dòng)起來(lái)改善水質(zhì). 本文基于EFDC三維水動(dòng)力模型，建立水齡模型來(lái)研究太湖水體交換速率. 基于實(shí)測(cè)徑流、風(fēng)場(chǎng)數(shù)據(jù)、引水工程的調(diào)水運(yùn)行數(shù)據(jù)，對(duì)太湖進(jìn)行三維水動(dòng)力模擬，研究風(fēng)場(chǎng)作用下太湖水齡分布的季節(jié)性特征以及“引江濟(jì)太”工程對(duì)太湖水齡分布特征的影響.

1 研究區(qū)域

太湖(30°55′40″～31°32′58″N，119°52′32″～120°36′10″E)橫跨江、浙兩省，西和西南側(cè)為丘陵山地，東側(cè)以平原和水網(wǎng)為主. 地處亞熱帶，氣候溫和濕潤(rùn)，屬季風(fēng)氣候，夏季受熱帶海洋氣團(tuán)影響，盛行東南季風(fēng)，溫和多雨；冬季受北方高壓氣團(tuán)控制，盛行偏北季風(fēng)，寒冷干燥. 太湖是典型的風(fēng)生流湖泊，且水深極淺，湖域面積幾乎全年不變. 太湖水域面積為2338.1 km2，平均年出湖徑流量75×108m3，蓄水量44×108m3，年降水量1100～1150 mm，年平均氣溫16～18℃, 平均水深1.89 m，最大水深2.6 m，年平均風(fēng)速3.5～5.0 m/s. 為分析“引江濟(jì)太”工程建設(shè)對(duì)太湖水齡的影響，環(huán)太湖共采集12個(gè)站點(diǎn)進(jìn)行分析對(duì)比，其中1#、5#～7#站點(diǎn)位于引水口，1#～4#、12#站點(diǎn)位于取水口，2#、10#站點(diǎn)為水位驗(yàn)證站，其余站點(diǎn)為流量控制點(diǎn)(圖1).

圖1 太湖地形及引水工程示意(1：望亭；2：犢山閘；3：梅梁湖泵站；4：大渲河泵站；5：龔巷橋；6：雅浦橋；7：百瀆口；8：浯溪橋；9：城東港；10：大浦口；11：長(zhǎng)兜港；12：太浦)Fig.1 Bathymetry and sketch for water transfer layout of Lake Taihu

2 數(shù)值模型

2.1 EFDC三維水動(dòng)力模型

2.1.1 三維水動(dòng)力模型的建立 EFDC模型是一個(gè)多參數(shù)的有限差分模型，可用于河流、湖泊、水庫(kù)、濕地系統(tǒng)、河口和海洋等水體的水動(dòng)力學(xué)和水質(zhì)模擬[6]. 本文將其應(yīng)用于太湖的研究. 控制方程在水平方向上采用笛卡爾直角坐標(biāo)系或者曲線正交坐標(biāo)系，垂向采用Sigma坐標(biāo). 在Sigma坐標(biāo)系下EFDC模型的連續(xù)性方程和動(dòng)量方程為[7-9]：

?t(mH)+?x(myHu)+?y(mxHv)+?z(mω)=QH

(1)

(2)

(3)

(τxz,τyz)=AvH-1?z(u,v)

(4)

(5)

(6)

式中，x、y為水平笛卡爾坐標(biāo)系；z為豎向Sigma坐標(biāo)；H為水深，ζ是水面高程，η為相對(duì)水位，定義ζ=H+η；u、v、w分別為x、y、z方向的流速；QH為體積源項(xiàng)，包括降雨、蒸發(fā)、滲流；Qu、Qv為x、y方向的動(dòng)量源項(xiàng)；p為附加靜水壓力；f為柯氏參數(shù)；τxz、τyz為x、y方向的垂向剪切應(yīng)力；Av為垂向湍流動(dòng)量擴(kuò)散系數(shù)；ρ和ρ0分別為實(shí)際和參考密度，b為浮力. mx、my為水平坐標(biāo)因子，且 m=mxmy.Uw和Vw為水平面10 m高度處風(fēng)速x和y方向的分量，風(fēng)應(yīng)力系數(shù)為：

(7)

式中，風(fēng)速單位為m/s，ρa(bǔ)和ρw分別代表空氣和水的密度.

2.1.2 模型主要計(jì)算參數(shù)以及模型的率定 本文采用2014年太湖的流量、氣象、水位以及地形資料，基于EFDC建立三維水動(dòng)力模型，EFDC模型的應(yīng)用已十分成熟，時(shí)間步長(zhǎng)通過(guò)CFL(科朗數(shù))穩(wěn)定條件確定. 糙率與水深、床面形態(tài)等因素有關(guān)，根據(jù)計(jì)算率定. 紊動(dòng)粘性系數(shù)采用Smagorinsky公式計(jì)，Cs為Smagorinsky系數(shù)，通常取0.28[8]. 數(shù)學(xué)模型的初始條件一般分熱啟動(dòng)和冷啟動(dòng)兩種. 前者一般通過(guò)計(jì)算得到計(jì)算區(qū)域內(nèi)相對(duì)準(zhǔn)確的水位、流速場(chǎng)，后者則假定一種均勻分布，與實(shí)際情況不一定相符，但通過(guò)邊界信息的傳入，誤差會(huì)逐漸消失；本文采用了冷啟動(dòng). 數(shù)學(xué)模型的邊界條件包括開(kāi)邊界(水邊界)和固定邊界(岸邊界). 因入湖河流的斷面尺度相對(duì)于湖域面積很小，故本文的計(jì)算將太湖視為一個(gè)封閉的區(qū)域，入湖河流的流量通過(guò)源匯項(xiàng)輸入. 模型網(wǎng)格采用正交曲線網(wǎng)格，最大網(wǎng)格尺度1000 m，最小網(wǎng)格800 m，網(wǎng)格數(shù)15625個(gè)，垂向分為4層，每層的厚度由湖底地形以及水面高程來(lái)決定. 模型主要受徑流、表層風(fēng)應(yīng)力的驅(qū)動(dòng)，溫度在模型中視為常數(shù). 模型的率定通過(guò)驗(yàn)證站的實(shí)測(cè)水位進(jìn)行，選取犢山閘和大浦口站的水位進(jìn)行模型的率定，模擬結(jié)果與實(shí)測(cè)值吻合較好，表明模型精度滿(mǎn)足要求(圖2). 由于實(shí)際情況下太湖的換水周期為300 d左右[3]，本文采用的流量邊界基本接近實(shí)際情況，因此模擬周期為365 d.

圖2 犢山閘(a)和大浦口(b)計(jì)算水位和實(shí)測(cè)水位的對(duì)比驗(yàn)證Fig.2 Comparison of the measured and simulated water levels in Dushan Gate (a) and Dapukou (b) sites

2.2 水齡的定義

水齡被定義為某個(gè)固定的水質(zhì)點(diǎn)從入口傳輸?shù)侥硞€(gè)指定點(diǎn)所用的時(shí)間，通常入口處的水齡設(shè)為0. 水齡可以通過(guò)系統(tǒng)中示蹤劑的濃度分布以及一個(gè)虛擬的水齡濃度來(lái)計(jì)算. 假設(shè)系統(tǒng)中只有一種示蹤劑，那么示蹤劑濃度和水齡濃度滿(mǎn)足以下兩個(gè)控制方程[10-16]：

(8)

(9)

(10)

基于EFDC的源程序，增加了水齡的數(shù)值模塊. 為了率定風(fēng)場(chǎng)作用下的水齡，與Li等[2]的計(jì)算結(jié)果進(jìn)行對(duì)比率定，在SE、N風(fēng)向下假設(shè)單入口流量下各測(cè)站水齡計(jì)算結(jié)果擬合結(jié)果良好(表1)，因此本文計(jì)算的風(fēng)場(chǎng)下水齡的精度滿(mǎn)足要求.

表1 本文水齡結(jié)果與Li等2011年水齡計(jì)算結(jié)果[2]對(duì)比

Tab.1 Comparison of water ages between present calculation and Lietal(2011)[2]

監(jiān)測(cè)點(diǎn)(風(fēng)向)Li等2011年計(jì)算的水齡/d本文計(jì)算的水齡/dB(N)310308B(SE)300295C(N)300300C(SE)200194

2.3 計(jì)算方案

本文的計(jì)算將太湖視為一個(gè)封閉的區(qū)域，太湖的環(huán)湖河道眾多，但其流量均較小，因此模型在設(shè)置流量邊界時(shí)，考慮上游主要入湖河流有從西苕溪水系流入的城東港、大渲河泵站、百瀆口和長(zhǎng)兜港，下游的主要出湖河流有位于東南部的太浦河，太浦河的出水量占太湖總出水量的65%左右[17]. 由于缺乏所有站點(diǎn)的實(shí)測(cè)資料，根據(jù)水量平衡原理，由2014年城東港、百瀆口、大渲河、浯溪橋、梅梁湖泵站、望亭站的實(shí)測(cè)資料以及2014年的出湖總流量(約為104×108m3)，推算出長(zhǎng)兜港以及太浦河的流量. 太湖是典型的風(fēng)生流湖泊，風(fēng)場(chǎng)在水體循環(huán)中起很重要的作用，為研究不同風(fēng)向?qū)μw交換的影響，對(duì)2014年NCEP再分析風(fēng)場(chǎng)資料(圖3)進(jìn)行分析后，得出春、夏、秋、冬的常風(fēng)向分別為SE、ESE、NNE和N，采用1994-2014年的多年平均風(fēng)速.在太湖湖流入口處持續(xù)釋放濃度為1 mg/L的示蹤劑，在不改變流量的情況下分別對(duì)四季不同風(fēng)向作用下的水齡進(jìn)行模擬計(jì)算，計(jì)算方案參見(jiàn)表2.

表2 計(jì)算方案

Tab.2 Calculation setup

計(jì)算方案風(fēng)向風(fēng)速/(m/s)染色劑入口(測(cè)站)出口(無(wú)染色劑測(cè)站)1SE(春季)3．442ESE(夏季)2．753NNE(秋季)3．074N(冬季)3．64城東港、大渲河、浯溪橋、百瀆口、長(zhǎng)兜港(各方案相同)梅梁湖、太浦河(各方案相同)

圖3 太湖風(fēng)速矢量圖 Fig.3 Wind velocity vector diagram of Lake Taihu

3 水齡的計(jì)算以及模擬結(jié)果分析

3.1 染色劑和流場(chǎng)的季節(jié)性分布特征

梅梁湖口的流速實(shí)測(cè)值在0.02～0.06 m/s之間[17]，本模型計(jì)算結(jié)果顯示流速在0.03～0.07 m/s之間，與實(shí)測(cè)值吻合較好，表明模型的精度滿(mǎn)足要求. 染色劑和流場(chǎng)分布參見(jiàn)圖4，由模擬結(jié)果可知，湖流方向隨風(fēng)向的改變而改變，染色劑從太湖入口處放入，隨著湖流方向開(kāi)始擴(kuò)散. 春、夏偏南向季風(fēng)作用下，竺山湖、西北湖區(qū)、西南湖區(qū)的染色劑濃度較大；秋、冬偏北向季風(fēng)作用下，竺山湖、湖心區(qū)以及東部湖區(qū)的染色劑濃度較大.

3.2 水齡季節(jié)性分布特征

基于染色劑模型得出太湖水齡分布(圖4). 結(jié)果表明，染色劑濃度較大的地方水齡較??；相反，染色劑濃度較小的地方水齡較大. 湖體水齡分布統(tǒng)計(jì)參見(jiàn)表3. 春、夏、秋、冬季的湖體的平均水齡分別為220.45、235.62、217.47、211.80 d. 秋、冬季的平均水齡小于春、夏季. 春、夏在偏南向季風(fēng)作用下，竺山湖、西北湖區(qū)以及西南湖區(qū)的水齡較小，200 d等水齡線分布范圍較?。磺?、冬季在偏北風(fēng)作用下，竺山湖、梅梁湖、西南湖區(qū)以及湖心區(qū)的水齡較小，東部湖區(qū)、東太湖湖區(qū)的水齡較春、夏季有所減小，200 d等水齡線分布范圍較大. 靠近湖流入口處的地方水齡較小，遠(yuǎn)離入口的地方水齡較大；水流流向與風(fēng)向一致時(shí)水齡減小，水循環(huán)加快，反之則減慢. 竺山湖是湖體水齡最小的地區(qū)，貢湖則是湖體水齡最大的地區(qū).

表3 不同水齡區(qū)間(天數(shù))所占湖區(qū)面積

Tab.3 The area covered by different water ages(days)

季節(jié)面積/km20～100d100～200d200～300d300d以上平均水齡/d春季422．36366．02714．10575．00220．45夏季309．06400．33598．28769．81235．62秋季433．69450．68554．53638．57217．47冬季417．85505．44627．87526．21211．80

3.3 “引江濟(jì)太”工程對(duì)水齡分布的影響

每年太湖都有低水位的“枯水期”，往年都是通過(guò)望虞河“引江濟(jì)太”，確保太湖水位保持在3 m以上. 為了更好地盤(pán)活太湖水體，在從長(zhǎng)江引水進(jìn)入太湖的同時(shí)，太湖流域內(nèi)在2014年開(kāi)始疏通新溝河、新孟河兩條河道，太湖將通過(guò)新孟河、望虞河補(bǔ)水，借助走馬塘、新溝河、太浦河向長(zhǎng)江排水，太湖將形成“兩進(jìn)三出”的循環(huán)系統(tǒng)，通過(guò)讓水體動(dòng)起來(lái)，改善水質(zhì). 為防止汛期太湖水位過(guò)高，新溝河將承擔(dān)太湖北部區(qū)域向長(zhǎng)江排水的功能. 同時(shí)，應(yīng)急情況下也能通過(guò)新溝河向太湖調(diào)水，改善梅梁湖水質(zhì).

據(jù)實(shí)測(cè)資料分析，2014年新溝河和望虞河的引水時(shí)間并不在同一區(qū)間，2014年新溝河在6-8月進(jìn)行引水20 d，共引水1.1522×108m3，而望虞河2014年1-3月、11-12月共引水130 d(圖5)，引水量10.5759×108m3. 據(jù)此本文對(duì)2014年望虞河引水和新溝河引水分別進(jìn)行了計(jì)算對(duì)比，在引水時(shí)放入1 mg/L的染色劑，其余湖流入口采用實(shí)測(cè)流量，不放入染色劑，通過(guò)染色劑的面積可以得出進(jìn)入湖體的引水水體面積，進(jìn)一步分析對(duì)比從而評(píng)估引水的效益，模擬結(jié)果參見(jiàn)表4和圖6. 望虞河的引水量大于新溝河，因此染色劑擴(kuò)散面積較大，染色劑濃度大于0.01 mg/L的面積約為新溝河的5倍；染色劑大于0.001 mg/L的面積約為太湖總面積的1/4，因此可以看出望虞河引水效果良好，表明引水量越大，引水效果越好. 但是望虞河引水量約為新溝河的10倍，染色劑面積未達(dá)到10倍，是由于新溝河入口附近水體交換較快.

圖4 季風(fēng)影響下太湖4個(gè)季節(jié)染色劑、流場(chǎng)和水齡分布Fig.4 Distributions of dye, flow field and water age of Lake Taihu in each season

圖5 2014年望亭站流量Fig.5 Discharge of Wangting site in 2014

圖6 2014年新溝河和望虞河引水運(yùn)行條件下染色劑及水齡分布(其余湖流不放入染色劑)Fig.6 Distributions of dye and water age under operation of Xingou River and Wangyu River transfer in 2014 (other inlets with 0 mg/L dye)

為進(jìn)一步評(píng)估引水方案對(duì)水齡的影響，在對(duì)實(shí)際風(fēng)場(chǎng)、流場(chǎng)下太湖的水齡進(jìn)行計(jì)算的前提下，分別在望虞河和新溝河引水時(shí)也放入染色劑，其余湖流入口采用實(shí)測(cè)流量，均放入1 mg/L染色劑，模擬結(jié)果參見(jiàn)圖7～8. 結(jié)果表明望虞河引水工程可以加快貢湖的水體交換，減小該地區(qū)的水齡；從新溝河引水可以加快梅梁湖地區(qū)的水體交換，改善梅梁湖的水質(zhì)，但是由于引水量太小，因此對(duì)水齡的影響并不明顯.

表4 染色劑分布面積

Tab.4 Distribution area of dye

引水河流引水量/(×108m3)染色劑濃度大于0．01mg/L的面積/km2染色劑濃度大于0．001mg/L的面積/km2新溝河1．152267．0484．35望虞河10．5759324．16625．67

圖7 實(shí)際風(fēng)場(chǎng)下染色劑和水齡分布(其余湖流放入染色劑) Fig.7 Distributions of dye and water age in actual wind field (other inlets with 1 mg/L dye)

圖8 2014年望虞河和新溝河引水運(yùn)行條件下染色劑及水齡分布(其余湖流放入染色劑)Fig.8 Distributions of dye and water age under operation of Wangyu River and Xingou River transfer in 2014 (other inlets with 1 mg/L dye)

4 結(jié)論

對(duì)于湖泊的水體交換評(píng)估，水齡是一個(gè)非常重要的指標(biāo). 它反映了水質(zhì)點(diǎn)從進(jìn)入湖體到離開(kāi)湖體的平均時(shí)間，并受到流量、風(fēng)場(chǎng)等因素的影響，因此可以用來(lái)評(píng)估水體的更新速率. 本文建立了基于EFDC的三維水動(dòng)力模型來(lái)研究太湖水齡的分布以及“引江濟(jì)太”工程對(duì)其的影響. 春、夏季在偏南向季風(fēng)作用下，竺山湖、西北湖區(qū)以及西南湖區(qū)的水齡較小，200 d等水齡線分布范圍較小；秋、冬季在偏北風(fēng)作用下，竺山湖、梅梁湖、西南湖區(qū)以及湖心區(qū)的水齡較小，東部湖區(qū)、東太湖湖區(qū)的水齡較春、夏季有所減小，200 d等水齡線分布范圍較大. 靠近入口處的地方水齡較小，遠(yuǎn)離入口的地方水齡較大；水流流向與風(fēng)向一致時(shí)，水齡減小，水循環(huán)加快，反之則減慢. 望虞河引水工程能夠減小貢湖湖區(qū)的水齡，加快該湖區(qū)及部分湖心區(qū)的水循環(huán)，盤(pán)活太湖的東部地區(qū)，對(duì)太湖水污染的治理具有重要意義；而2014年新溝河工程采用的引水方案可以減小梅梁湖區(qū)的水齡，加快梅梁湖區(qū)的水循環(huán)，改善梅梁湖水質(zhì)，并且通過(guò)流水帶動(dòng)，改善入湖口直湖港、武進(jìn)港的水質(zhì). 在需要引水時(shí)，新孟河、新溝河、望虞河均可引水，可以組合使用，以實(shí)現(xiàn)最優(yōu)引水方案；需要排水時(shí)太浦河、新溝河、走馬塘均可排水，太浦河是主要通道. 引水工程的實(shí)施對(duì)加快整個(gè)太湖水體循環(huán)做出了重要貢獻(xiàn).

致謝：感謝江蘇省水文水資源勘測(cè)局無(wú)錫分局的各位同仁在野外采樣工作中給予的無(wú)私幫助.

[1] Deleersnijder E, Campin JM, Delhez EJM. The concept of age in marine modelling I. Theory and preliminary model results.JournalofMarineSystems, 2001, 28(3/4): 229-267．

[2] Li YP, Kumud A, Yu ZB. Modeling impacts of Yangtze River water transfer on water ages in Lake Taihu, China.EcologicalEngineering, 2011, 37(2): 325-334.

[3] Gao Yi, Mao Xinwei, Xu Weidongetal. Analysis of the influence on the Taihu Lake and the area around： Diversion from the Yangtze River to the Taihu Lake.JournalofChinaHydrology, 2006, 26(1): 92-94. [高怡, 毛新偉, 徐衛(wèi)東等. “引江濟(jì)太”工程對(duì)太湖及周邊地區(qū)的影響分析. 水文, 2006, 26(1): 92-94.]

[4] Wu Haoyun. Study on key technologies of Yangtze-Taihu Water Diversion.ChinaWaterResources, 2008, (1): 6-8. [吳浩云. 引江濟(jì)太調(diào)水試驗(yàn)關(guān)鍵技術(shù)研究和應(yīng)用. 中國(guó)水利, 2008, (1): 6-8.]

[5] Wu Haoyun, Zhou Danping, He Jiaetal. Integrated benefit assessment of the project water diversion from Yangtze River to Lake Taihu and discussion on the methodology.JLakeSci, 2008, 20(5): 639-647. DOI:10.18307/2008.0513. [吳浩云, 周丹平, 何佳等. 引江濟(jì)太工程綜合效益的評(píng)估及方法探討. 湖泊科學(xué), 2008, 20(5): 639-647.]

[6] Zhang Yifei, Wang Yulin, Wang Liang. EFDC overview and application analysis.EnvironmentalImpactAssessment, 2015, (3): 70-72, 92. [張以飛,王玉琳,汪靚. EFDC模型概述與應(yīng)用分析. 環(huán)境影響評(píng)價(jià), 2015, (3): 70-72, 92.]

[7] Iglesias C, Mazzeo N, Meerhoff Metal. High predation is of key importance for dominance of small-bodied zooplankton in warm shallow lakes: Evidence from lakes, fish exclosures and surface sediments.Hydrobiologia, 2011, 667: 133-147.

[8] Ji ZG, Hu G, Shen Jetal. Three-dimensional modeling of hydrodynamic processes in the St. Lucie Estuary.Estuarine,CoastalandShelfScience, 2007, 73: 188-200.

[9] Lu Hui, Ning Yawei, Yuan Yongling. EFDC overview and application analysis.WaterResourcesandPower, 2013, (4): 100-102. [盧慧, 寧亞偉, 袁永齡. 基于EFDC模型的人工湖生態(tài)換水優(yōu)化計(jì)算. 水電能源科學(xué), 2013, (4): 100-102.]

[10] Li XW, Shen YM. Numerical simulation of the impacts of water level variation on water age in Dahuofang Reservoir.FrontiersinEarthScience, 2015, 9(2): 209-224.

[11] Deleersnijder E, Campin JM, Delhez EJM. The concept of age in marine modelling I. Theory and preliminary model results.JournalofMarineSystems, 2001, 28(3/4): 229-267.

[12] Xiong Xiong. Hydrodynamic research of Poyang Lake based on the EFDC[Dissertation]. Nanchang: Nanchang University, 2014. [熊雄. 基于EFDC的鄱陽(yáng)湖水動(dòng)力研究[學(xué)位論文]. 南昌: 南昌大學(xué), 2014.]

[13] Shao Junrong, Wu Shiqiang, Zhou Jieetal. An age model for water transfer.AdvancesinWaterScience, 2014, 25(5): 695-703.[邵軍榮, 吳時(shí)強(qiáng), 周杰等. 水體交換年齡模型研究. 水科學(xué)進(jìn)展, 2014, 25(5): 695-703.]

[14] Liu Lüliu, Yang Zhifeng, Shen Zhenyaoetal. Temporal characteristics of exchange and transit of water bodies.JournalofNaturalResources, 2003, 18(1): 87-93. [劉綠柳, 楊志峰, 沈珍瑤等. 水體交換與傳輸?shù)臅r(shí)間維特征. 自然資源學(xué)報(bào), 2003, 18(1): 87-93.]

[15] Qi Wen, Li Yiping, Wang Yingetal. Control and optimization scheme of TianJing Lake hydrodynamics based on EFDC Model.WaterResourcesandPower, 2015, (1): 47-51. [戚文, 李一平, 王瑩等. 基于EFDC模型的天鏡湖水動(dòng)力優(yōu)化調(diào)控方案研究. 水電能源科學(xué), 2015, (1): 47-51.]

[16] Gong Ran, Xu Jin, Xu Ligangetal. Study on hydrodynamics in urban landscape lake based on EFDC.EnvironmentalEngineering, 2015, 33(4):58-62.[龔然, 徐進(jìn), 徐力剛等. 基于EFDC城市景觀湖泊水動(dòng)力模擬研究. 環(huán)境工程, 2015, 33(4): 58-62.]

[17] Qin Boqiang ed. Lake Taihu, China, dynamics and environmental change. Springer, 2008: 125.

Water age distribution of the Lake Taihu and impact of the Yangtze River to Lake Taihu Water Transfer Project on the water age

HUANG Chunlin1, LI Xi1**& SUN Yongyuan2

(1:KeyLaboratoryofCoastalDisasterandDefence,MinstryofEducation,HohaiUniversity,Nanjing210098,P.R.China)(2:JiangsuProvinceHydrologyandWaterResourceInvestigationBureau,Nangjing210029,P.R.China)

Lake Taihu is a typical wind-driven lake, and wind plays an important role in process of water exchange and concentration transportation. Based on Environmental Fluid Dynamics Model (EFDC), water age model and conservative dye model were set up to study characteristics of seasonal water age distribution and to further reveal the impacts of the Wangyu River Water Transfer Project and the Xingou River to the Lake Taihu Transfer Project on water age distribution on basis of winds, inflow and outflow, operations of existing water transfer projects. It is concluded that by the numerical test results, water exchange of Lake Taihu is mainly influenced by the monsoons. Water age is comparably larger in the spring and summer, which indicates the water exchange rate in the spring and summer is lower than that in the autumn and winter. Water ages near inflow locations are comparably smaller, and the locations far from the inflows are larger. Water age is smaller when the flow direction is the same as wind direction, while the exchange of water is increased at the same time and on the contrary water exchange is reduced. During operation of the Wangyu River Water Transfer Project, the water ages of the Gonghu Bay and Middle Lake Taihu are basically reduced, and water exchange is increased in the areas. During operation of the Xingou River Water Transfer Project, the water ages of the Meiliang Bay are basically reduced, and water quality is improved in the area. The water exchange of Lake Taihu will be significantly improved with continuous or occasional transfer operations.

Lake Taihu; wind; dye; water age; Environmental Fluid Dynamics Model

*江蘇省水利科技項(xiàng)目(2015054)和江西省水利科技項(xiàng)目(KT201405)聯(lián)合資助. 2016-06-06收稿；2016-07-03收修改稿. 黃春琳(1992～)，女，碩士研究生；E-mail：805282052@qq.com.

*通信作者；E-mail：xili@hhu.edu.cn.

J.LakeSci.(湖泊科學(xué))， 2017， 29(1)： 22-31

DOI 10.18307/2017.0103

?2017 byJournalofLakeSciences