太湖水龄分布特征及“引江济太”工程对其的影响

黄春琳,李 熙,孙永远

(1:河海大学海岸灾害及防护教育部重点实验室，南京 210098)(2：江苏省水文水资源勘测局，南京 210029)

太湖水龄分布特征及“引江济太”工程对其的影响

黄春琳1,李 熙1,孙永远2

(1:河海大学海岸灾害及防护教育部重点实验室，南京 210098)(2：江苏省水文水资源勘测局，南京 210029)

太湖作为典型的风生流湖泊，风场对水体运输和交换的过程起着重要的作用. 基于环境水动力学模型EFDC源程序建立了染色剂模型以及水龄模型，借助水龄研究太湖水体的长期输运过程和更新速率特征，从而为太湖的调水工程管理提供科学依据. 本文在研究春、夏、秋、冬季不同风场作用下太湖的水龄季节性分布特征的基础上，结合太湖实测风速、流量数据及“引江济太”工程调水运行的现状，着重分析了望虞河枢纽调水运行以及新沟河工程对太湖水龄分布的影响. 数值试验的模拟结果表明：太湖水体交换受季风影响明显，春、夏季的水龄相对较大、水体交换较差；靠近湖流入口处的地方水龄较小，远离入口的地方水龄较大；水流流向与风向一致时水龄减小，水循环加快，反之则减慢. 望虞河引水引工程能够减小贡湖及湖心区的水龄，加快贡湖湖区及湖心区的水循环；新沟河工程引水能够减小梅梁湖区的水龄，改善该湖区的水质. 引水工程的实施对加快整个太湖的水循环做出重要贡献.

太湖；风场；染色剂；水龄；EFDC模型

入湖流量、风场等在湖体的水体运动过程中起很重要的作用，三维环境流体动力学模型EFDC提供了一个高分辨率的工具来研究太湖的水体交换. 水龄可以定量反映水体的运动和交换程度以及滞留情况，对评估水质的变化具有重要意义. Deleersnijder等[1]介绍了基于对流-扩散输运方程的水龄理论，这个新的水龄定义提供了一个健全的方法来定量地研究欧拉坐标框架下的物质溶解的输运特征. Li等[2]运用EFDC三维水动力模型，借助水龄的定义，研究了“引江济太”引水量对太湖水循环的影响.

太湖作为中国的第三大淡水湖，由于富营养化而严重影响了饮用水的安全，政府部门近年来采取了许多措施来提升太湖水质，包括建更多的污水处理厂，禁止使用含磷清洁剂等. 此外，一个重要的生态工程——“引江济太”工程在2002年启动，“引江济太”工程将更优质的长江水通过望虞河引入太湖，最终从太浦河流出，以期通过此举来加快太湖的水循环，改善太湖的水质状况[3-5]. 此外，在2014年太湖北部流域内开始疏通新沟河、新孟河两条河道，太湖将通过新孟河、望虞河补水，借助走马塘、新沟河、太浦河向长江排水，形成“两进三出”的循环系统，通过让水体动起来改善水质. 本文基于EFDC三维水动力模型，建立水龄模型来研究太湖水体交换速率. 基于实测径流、风场数据、引水工程的调水运行数据，对太湖进行三维水动力模拟，研究风场作用下太湖水龄分布的季节性特征以及“引江济太”工程对太湖水龄分布特征的影响.

1 研究区域

太湖(30°55′40″～31°32′58″N，119°52′32″～120°36′10″E)横跨江、浙两省，西和西南侧为丘陵山地，东侧以平原和水网为主. 地处亚热带，气候温和湿润，属季风气候，夏季受热带海洋气团影响，盛行东南季风，温和多雨；冬季受北方高压气团控制，盛行偏北季风，寒冷干燥. 太湖是典型的风生流湖泊，且水深极浅，湖域面积几乎全年不变. 太湖水域面积为2338.1 km2，平均年出湖径流量75×108m3，蓄水量44×108m3，年降水量1100～1150 mm，年平均气温16～18℃, 平均水深1.89 m，最大水深2.6 m，年平均风速3.5～5.0 m/s. 为分析“引江济太”工程建设对太湖水龄的影响，环太湖共采集12个站点进行分析对比，其中1#、5#～7#站点位于引水口，1#～4#、12#站点位于取水口，2#、10#站点为水位验证站，其余站点为流量控制点(图1).

图1 太湖地形及引水工程示意(1：望亭；2：犊山闸；3：梅梁湖泵站；4：大渲河泵站；5：龚巷桥；6：雅浦桥；7：百渎口；8：浯溪桥；9：城东港；10：大浦口；11：长兜港；12：太浦)Fig.1 Bathymetry and sketch for water transfer layout of Lake Taihu

2 数值模型

2.1 EFDC三维水动力模型

2.1.1 三维水动力模型的建立 EFDC模型是一个多参数的有限差分模型，可用于河流、湖泊、水库、湿地系统、河口和海洋等水体的水动力学和水质模拟[6]. 本文将其应用于太湖的研究. 控制方程在水平方向上采用笛卡尔直角坐标系或者曲线正交坐标系，垂向采用Sigma坐标. 在Sigma坐标系下EFDC模型的连续性方程和动量方程为[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为水平笛卡尔坐标系；z为竖向Sigma坐标；H为水深，ζ是水面高程，η为相对水位，定义ζ=H+η；u、v、w分别为x、y、z方向的流速；QH为体积源项，包括降雨、蒸发、渗流；Qu、Qv为x、y方向的动量源项；p为附加静水压力；f为柯氏参数；τxz、τyz为x、y方向的垂向剪切应力；Av为垂向湍流动量扩散系数；ρ和ρ0分别为实际和参考密度，b为浮力. mx、my为水平坐标因子，且 m=mxmy.Uw和Vw为水平面10 m高度处风速x和y方向的分量，风应力系数为：

(7)

式中，风速单位为m/s，ρa和ρw分别代表空气和水的密度.

2.1.2 模型主要计算参数以及模型的率定 本文采用2014年太湖的流量、气象、水位以及地形资料，基于EFDC建立三维水动力模型，EFDC模型的应用已十分成熟，时间步长通过CFL(科朗数)稳定条件确定. 糙率与水深、床面形态等因素有关，根据计算率定. 紊动粘性系数采用Smagorinsky公式计，Cs为Smagorinsky系数，通常取0.28[8]. 数学模型的初始条件一般分热启动和冷启动两种. 前者一般通过计算得到计算区域内相对准确的水位、流速场，后者则假定一种均匀分布，与实际情况不一定相符，但通过边界信息的传入，误差会逐渐消失；本文采用了冷启动. 数学模型的边界条件包括开边界(水边界)和固定边界(岸边界). 因入湖河流的断面尺度相对于湖域面积很小，故本文的计算将太湖视为一个封闭的区域，入湖河流的流量通过源汇项输入. 模型网格采用正交曲线网格，最大网格尺度1000 m，最小网格800 m，网格数15625个，垂向分为4层，每层的厚度由湖底地形以及水面高程来决定. 模型主要受径流、表层风应力的驱动，温度在模型中视为常数. 模型的率定通过验证站的实测水位进行，选取犊山闸和大浦口站的水位进行模型的率定，模拟结果与实测值吻合较好，表明模型精度满足要求(图2). 由于实际情况下太湖的换水周期为300 d左右[3]，本文采用的流量边界基本接近实际情况，因此模拟周期为365 d.

图2 犊山闸(a)和大浦口(b)计算水位和实测水位的对比验证Fig.2 Comparison of the measured and simulated water levels in Dushan Gate (a) and Dapukou (b) sites

2.2 水龄的定义

水龄被定义为某个固定的水质点从入口传输到某个指定点所用的时间，通常入口处的水龄设为0. 水龄可以通过系统中示踪剂的浓度分布以及一个虚拟的水龄浓度来计算. 假设系统中只有一种示踪剂，那么示踪剂浓度和水龄浓度满足以下两个控制方程[10-16]：

(8)

(9)

(10)

基于EFDC的源程序，增加了水龄的数值模块. 为了率定风场作用下的水龄，与Li等[2]的计算结果进行对比率定，在SE、N风向下假设单入口流量下各测站水龄计算结果拟合结果良好(表1)，因此本文计算的风场下水龄的精度满足要求.

表1 本文水龄结果与Li等2011年水龄计算结果[2]对比

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

监测点(风向)Li等2011年计算的水龄/d本文计算的水龄/dB(N)310308B(SE)300295C(N)300300C(SE)200194

2.3 计算方案

本文的计算将太湖视为一个封闭的区域，太湖的环湖河道众多，但其流量均较小，因此模型在设置流量边界时，考虑上游主要入湖河流有从西苕溪水系流入的城东港、大渲河泵站、百渎口和长兜港，下游的主要出湖河流有位于东南部的太浦河，太浦河的出水量占太湖总出水量的65%左右[17]. 由于缺乏所有站点的实测资料，根据水量平衡原理，由2014年城东港、百渎口、大渲河、浯溪桥、梅梁湖泵站、望亭站的实测资料以及2014年的出湖总流量(约为104×108m3)，推算出长兜港以及太浦河的流量. 太湖是典型的风生流湖泊，风场在水体循环中起很重要的作用，为研究不同风向对太湖水体交换的影响，对2014年NCEP再分析风场资料(图3)进行分析后，得出春、夏、秋、冬的常风向分别为SE、ESE、NNE和N，采用1994-2014年的多年平均风速.在太湖湖流入口处持续释放浓度为1 mg/L的示踪剂，在不改变流量的情况下分别对四季不同风向作用下的水龄进行模拟计算，计算方案参见表2.

表2 计算方案

Tab.2 Calculation setup

计算方案风向风速/(m/s)染色剂入口(测站)出口(无染色剂测站)1SE(春季)3．442ESE(夏季)2．753NNE(秋季)3．074N(冬季)3．64城东港、大渲河、浯溪桥、百渎口、长兜港(各方案相同)梅梁湖、太浦河(各方案相同)

图3 太湖风速矢量图 Fig.3 Wind velocity vector diagram of Lake Taihu

3 水龄的计算以及模拟结果分析

3.1 染色剂和流场的季节性分布特征

梅梁湖口的流速实测值在0.02～0.06 m/s之间[17]，本模型计算结果显示流速在0.03～0.07 m/s之间，与实测值吻合较好，表明模型的精度满足要求. 染色剂和流场分布参见图4，由模拟结果可知，湖流方向随风向的改变而改变，染色剂从太湖入口处放入，随着湖流方向开始扩散. 春、夏偏南向季风作用下，竺山湖、西北湖区、西南湖区的染色剂浓度较大；秋、冬偏北向季风作用下，竺山湖、湖心区以及东部湖区的染色剂浓度较大.

3.2 水龄季节性分布特征

基于染色剂模型得出太湖水龄分布(图4). 结果表明，染色剂浓度较大的地方水龄较小；相反，染色剂浓度较小的地方水龄较大. 湖体水龄分布统计参见表3. 春、夏、秋、冬季的湖体的平均水龄分别为220.45、235.62、217.47、211.80 d. 秋、冬季的平均水龄小于春、夏季. 春、夏在偏南向季风作用下，竺山湖、西北湖区以及西南湖区的水龄较小，200 d等水龄线分布范围较小；秋、冬季在偏北风作用下，竺山湖、梅梁湖、西南湖区以及湖心区的水龄较小，东部湖区、东太湖湖区的水龄较春、夏季有所减小，200 d等水龄线分布范围较大. 靠近湖流入口处的地方水龄较小，远离入口的地方水龄较大；水流流向与风向一致时水龄减小，水循环加快，反之则减慢. 竺山湖是湖体水龄最小的地区，贡湖则是湖体水龄最大的地区.

表3 不同水龄区间(天数)所占湖区面积

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

季节面积/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 “引江济太”工程对水龄分布的影响

每年太湖都有低水位的“枯水期”，往年都是通过望虞河“引江济太”，确保太湖水位保持在3 m以上. 为了更好地盘活太湖水体，在从长江引水进入太湖的同时，太湖流域内在2014年开始疏通新沟河、新孟河两条河道，太湖将通过新孟河、望虞河补水，借助走马塘、新沟河、太浦河向长江排水，太湖将形成“两进三出”的循环系统，通过让水体动起来，改善水质. 为防止汛期太湖水位过高，新沟河将承担太湖北部区域向长江排水的功能. 同时，应急情况下也能通过新沟河向太湖调水，改善梅梁湖水质.

据实测资料分析，2014年新沟河和望虞河的引水时间并不在同一区间，2014年新沟河在6-8月进行引水20 d，共引水1.1522×108m3，而望虞河2014年1-3月、11-12月共引水130 d(图5)，引水量10.5759×108m3. 据此本文对2014年望虞河引水和新沟河引水分别进行了计算对比，在引水时放入1 mg/L的染色剂，其余湖流入口采用实测流量，不放入染色剂，通过染色剂的面积可以得出进入湖体的引水水体面积，进一步分析对比从而评估引水的效益，模拟结果参见表4和图6. 望虞河的引水量大于新沟河，因此染色剂扩散面积较大，染色剂浓度大于0.01 mg/L的面积约为新沟河的5倍；染色剂大于0.001 mg/L的面积约为太湖总面积的1/4，因此可以看出望虞河引水效果良好，表明引水量越大，引水效果越好. 但是望虞河引水量约为新沟河的10倍，染色剂面积未达到10倍，是由于新沟河入口附近水体交换较快.

图4 季风影响下太湖4个季节染色剂、流场和水龄分布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年新沟河和望虞河引水运行条件下染色剂及水龄分布(其余湖流不放入染色剂)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)

为进一步评估引水方案对水龄的影响，在对实际风场、流场下太湖的水龄进行计算的前提下，分别在望虞河和新沟河引水时也放入染色剂，其余湖流入口采用实测流量，均放入1 mg/L染色剂，模拟结果参见图7～8. 结果表明望虞河引水工程可以加快贡湖的水体交换，减小该地区的水龄；从新沟河引水可以加快梅梁湖地区的水体交换，改善梅梁湖的水质，但是由于引水量太小，因此对水龄的影响并不明显.

表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 实际风场下染色剂和水龄分布(其余湖流放入染色剂) Fig.7 Distributions of dye and water age in actual wind field (other inlets with 1 mg/L dye)

图8 2014年望虞河和新沟河引水运行条件下染色剂及水龄分布(其余湖流放入染色剂)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 结论

对于湖泊的水体交换评估，水龄是一个非常重要的指标. 它反映了水质点从进入湖体到离开湖体的平均时间，并受到流量、风场等因素的影响，因此可以用来评估水体的更新速率. 本文建立了基于EFDC的三维水动力模型来研究太湖水龄的分布以及“引江济太”工程对其的影响. 春、夏季在偏南向季风作用下，竺山湖、西北湖区以及西南湖区的水龄较小，200 d等水龄线分布范围较小；秋、冬季在偏北风作用下，竺山湖、梅梁湖、西南湖区以及湖心区的水龄较小，东部湖区、东太湖湖区的水龄较春、夏季有所减小，200 d等水龄线分布范围较大. 靠近入口处的地方水龄较小，远离入口的地方水龄较大；水流流向与风向一致时，水龄减小，水循环加快，反之则减慢. 望虞河引水工程能够减小贡湖湖区的水龄，加快该湖区及部分湖心区的水循环，盘活太湖的东部地区，对太湖水污染的治理具有重要意义；而2014年新沟河工程采用的引水方案可以减小梅梁湖区的水龄，加快梅梁湖区的水循环，改善梅梁湖水质，并且通过流水带动，改善入湖口直湖港、武进港的水质. 在需要引水时，新孟河、新沟河、望虞河均可引水，可以组合使用，以实现最优引水方案；需要排水时太浦河、新沟河、走马塘均可排水，太浦河是主要通道. 引水工程的实施对加快整个太湖水体循环做出了重要贡献.

致谢：感谢江苏省水文水资源勘测局无锡分局的各位同仁在野外采样工作中给予的无私帮助.

[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. [高怡, 毛新伟, 徐卫东等. “引江济太”工程对太湖及周边地区的影响分析. 水文, 2006, 26(1): 92-94.]

[4] Wu Haoyun. Study on key technologies of Yangtze-Taihu Water Diversion.ChinaWaterResources, 2008, (1): 6-8. [吴浩云. 引江济太调水试验关键技术研究和应用. 中国水利, 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. [吴浩云, 周丹平, 何佳等. 引江济太工程综合效益的评估及方法探讨. 湖泊科学, 2008, 20(5): 639-647.]

[6] Zhang Yifei, Wang Yulin, Wang Liang. EFDC overview and application analysis.EnvironmentalImpactAssessment, 2015, (3): 70-72, 92. [张以飞,王玉琳,汪靓. EFDC模型概述与应用分析. 环境影响评价, 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模型的人工湖生态换水优化计算. 水电能源科学, 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的鄱阳湖水动力研究[学位论文]. 南昌: 南昌大学, 2014.]

[13] Shao Junrong, Wu Shiqiang, Zhou Jieetal. An age model for water transfer.AdvancesinWaterScience, 2014, 25(5): 695-703.[邵军荣, 吴时强, 周杰等. 水体交换年龄模型研究. 水科学进展, 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. [刘绿柳, 杨志峰, 沈珍瑶等. 水体交换与传输的时间维特征. 自然资源学报, 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模型的天镜湖水动力优化调控方案研究. 水电能源科学, 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.[龚然, 徐进, 徐力刚等. 基于EFDC城市景观湖泊水动力模拟研究. 环境工程, 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

*江苏省水利科技项目(2015054)和江西省水利科技项目(KT201405)联合资助. 2016-06-06收稿；2016-07-03收修改稿. 黄春琳(1992～)，女，硕士研究生；E-mail：805282052@qq.com.

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

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

DOI 10.18307/2017.0103

©2017 byJournalofLakeSciences