湖泊科学   2018, Vol. 30 Issue (1): 199-210.  DOI: 10.18307/2018.0120.
0

研究论文

引用本文 [复制中英文]

张素香, 许新发, 徐兴, Chu Vicnent H, 邬年华, 付莎莎, 李熙, 鄱阳湖水龄季节性变化特征. 湖泊科学, 2018, 30(1): 199-210. DOI: 10.18307/2018.0120.
[复制中文]
ZHANG Suxiang, XU Xinfa, XU Xing, CHU Vicnent H, WU Nianhua, FU Shasha, LI Xi. Seasonal variations of water age in Lake Poyang. Journal of Lake Sciences, 2018, 30(1): 199-210. DOI: 10.18307/2018.0120.
[复制英文]

基金项目

国家自然科学基金项目(51369011)、江西省水利科技项目(KT201405)和江苏省水利科技项目(2017039)联合资助

作者简介

张素香(1971~), 女, 讲师, 博士; E-mail:zhangsuxiang@nuist.edu.cn

通信作者

李熙; E-mail:xili@hhu.edu.cn

文章历史

2017-03-23 收稿
2017-05-22 收修改稿

码上扫一扫

鄱阳湖水龄季节性变化特征
张素香 1,4, 许新发 2, 徐兴 3, Chu Vicnent H 4, 邬年华 2, 付莎莎 2, 李熙 5     
(1: 南京信息工程大学气象灾害教育部重点实验室, 南京 210044)
(2: 江西省水利科学研究院, 南昌 330029)
(3: 江苏省水文水资源勘测局无锡分局, 无锡 214031)
(4: McGill University, Department of Civil Engineering and Applied Mechanics, Montreal, Quebec H3A 2K6, Canada)
(5: 河海大学港口海岸与近海工程学院, 南京 210098)
摘要:基于环境水动力学模型EFDC源程序,建立了染色剂模型和水龄模型,在将模型与航测水文数据验证吻合的基础上,分别计算了鄱阳湖自然条件下春、夏、秋、冬季的水龄和倒灌前后鄱阳湖染色剂和水龄分布的变化,以及五河水系各分支河流水龄.分季节的水龄计算表明鄱阳湖水体交换受季节性来水影响明显.夏、秋季的水龄相对较小,在多数年份又受到长江水倒灌的影响导致水龄有所增大;冬、春季水龄较大,亦无长江水倒灌现象,相较于夏、秋季,水域面积明显减少.分支流的水龄计算表明,西南湖区的水体交换主要受到赣江的影响,西北湖区水体交换主要受到修水和赣江的影响,南部湖区主要受到抚河与信江的影响,东部湖区主要受到饶河的影响,湖心区和入江水道则受到五河水系的综合影响.同时水龄的研究表明拟建的鄱阳湖水利枢纽工程"调枯不调洪"的原则是合理的,为鄱阳湖水利枢纽工程论证提供了重要的参考依据.
关键词鄱阳湖    鄱阳湖水利枢纽工程    EFDC    风场    染色剂    水龄    季节性变化    
Seasonal variations of water age in Lake Poyang
ZHANG Suxiang 1,4, XU Xinfa 2, XU Xing 3, CHU Vicnent H 4, WU Nianhua 2, FU Shasha 2, LI Xi 5     
(1: Key Laboratory of Meteorological Disaster of Ministry of Education, Nanjing University of Information Science & Technology, Nanjing 210044, P. R. China)
(2: Jiangxi Provincial Research Institute of Water Science, Nanchang 330029, P. R. China)
(3: Wuxi Branch of Jiangsu Province Hydrology and Water Resource Investigation Bureau, Wuxi 214031, P. R. China)
(4: Department of Civil Engineering and Applied Mechanics, McGill University, Montreal, Quebec H3A 2K6, Canada)
(5: College of Harbor, Coastal and Offshore Engineering, Hohai University, Nanjing 210098, P. R. China)
Abstract: Based on source code of Environmental Fluid Dynamics Model (EFDC), the water age model and conservative dye model were set up to study water age of seasons and individual rivers of Lake Poyang under natural conditions. Water age calculations indicate that seasonal variations of water age in Lake Poyang are significant. The water age of summer or autumn is smaller, which is increased during reversing process in flood season in most of consecutive years. The water age of winter or spring is larger without occurrence of reversing process. The water age of individual river in five river systems indicates that the exchange of water in the southwest lake zone is mainly influenced by Ganjiang River, whereas in the northwest zone influenced by Ganjiang River and Xinjiang River, in the southern lake zone by Ganjiang River and Xinjiang River, in the eastern lake zone by Raohe River, respectively. The central zone and channel to Yangtze River are influenced by the five rivers together. The proposed Lake Poyang Complex will be rational with principles "water level control during dry seasons not during flood seasons", and the model is proved of reference value for feasibility study of Lake Poyang Complex.
Keywords: Lake Poyang    Lake Poyang Complex    Environmental Fluid Dynamics Model (EFDC)    wind    dye    water age    seasonal variation    

对于水体输运时间尺度,不同学者使用不同的模型和定义来计算分析. Deleersnijder等[1]介绍了基于对流扩散输运方程的水龄理论,提供了定量地研究保守染色剂输运特征的理论和方法. Shen等[2]运用CH3D模型(Chesapeake Bay model in 3 Dimensions)来研究Chesapeake湾的水龄,结果表明,密度流使得水体更易向下游运动,水体的垂直和侧向混合作用则会减弱密度流. Shen等[3]基于FVCOM模型,研究了径流量、风场和地表热通量对大伙房水库水龄和停留时间的影响,并将滞留时间和水龄进行了对比分析. Li等[4]选取了国内外常用的水动力模型(EFDC)和典型的浅水湖泊(太湖),以水龄作为指标,利用三维数值模型来衡量从长江调水到太湖产生的稀释和转移污染物的效果.龚然等[5]利用EFDC模型对天印湖的水动力进行了模拟,计算结果显示,天印湖的水龄主要随降水量和径流量的改变而改变,降雨量和径流量越大,水龄越小,而风场作用对水龄的影响非常小. Ribbe等[6]使用COHERENS模型来研究Hervey湾的水体更新速率,研究结果显示,Hervey湾内85 %的水体可以在50~80 d内完成更新. Huang等[7]使用淡水留分法研究了流量对Apalachicola湾水体停留时间的影响,结果表明,当径流量从177 m3/s增大到4561 m3/s时,停留时间从10 d降到3 d,这表示增大流量将会导致停留时间减少. Kennedy等[8]基于FLUENT软件,利用粒子跟踪法计算曼彻斯特Wachusett水库水体平均停留时间为3~4 d. Dabrowski等[9]使用ECOMSED模型计算爱尔兰海水体停留时间和空间滞留时间,结果表明,一些地区的冬季平均停留时间比夏季高出28 %.综合以上分析,水龄主要用于水体内部物质交换与传输时间的计算,反映保守染色剂的动态时间变化;滞留时间主要用于反映系统内水团流出前在系统内的停留时间;而更替周期则用于描述封闭水体与外界交换的总体时间特征.鄱阳湖作为典型过水型浅水湖泊,其湖泊动力特征主要受五河(修水、赣江、抚河、信江、饶河)水系来水的影响,具有流速大、水域面积季节性变化剧烈的特点;同时鄱阳湖与长江直接相连接,依据水位变化和长江水体发生季节性交换,所以也具有吞吐型湖泊的特点.因此,针对鄱阳湖典型过水型和吞吐型浅水湖泊的特征,研究鄱阳湖水体交换的季节性变化更适合使用水龄的计算方法.

迄今,水龄的相关模拟研究以采用二维模型居多[10-11],水龄以保守染色剂的浓度分布变化来反映,通常染色剂空间浓度分布具有典型的三维特征,理论上与水流的紊动切应力呈正相关,因此构建三维水动力模型框架下的水龄模型是目前区域性湖泊河网水体交换和水质模拟的重要发展方向. EFDC三维水动力模型提供了一个高分辨率的框架性数值模型工具来研究鄱阳湖的水体交换过程,建立并嵌入三维EFDC模型中的水龄数值模型,提供了一个健全的方法来定量研究染色剂的输运特征[12].

1 研究区域概况

鄱阳湖作为可以调节长江洪水的吞吐型和过水型湖泊,其水量和水质的变化关系到鄱阳湖流域和长江中下游地区的用水安全.鄱阳湖地处江西省北部,长江中下游南岸(图 1),以松门山为界,分为南北两部分,北面为入江水道(湖口水道),长40 km;南面为主湖体,长133 km.因为和长江直接相连,洪季长江水倒灌入鄱阳湖,使鄱阳湖有着独一无二的水动力特性.鄱阳湖作为中国最大的淡水湖泊,有着丰富的湿地资源,生物种类繁多.然而,近年来土地利用、水资源开发利用以及气候变暖使鄱阳湖盆地湿地区域面积急剧减少.自2000年后,鄱阳湖已经连续多年出现旱情,根据水文数据监测,鄱阳湖如今在夏季丰水期水域面积最大可能达4000 km2以上,到冬季枯水期却能减少到几百平方公里,尤其是三峡水库运行后,枯水情况日益加剧,因此加强鄱阳湖的治理和保护刻不容缓.本文基于EFDC三维水动力模型[13-16],建立染色剂模型和水龄模型,分别计算了春、夏、秋、冬季鄱阳湖的水龄以及倒灌前后鄱阳湖染色剂分布的变化;基于水龄评估了鄱阳湖的水体交换特征以及长江水倒灌对水体交换的影响;通过水龄模型对拟建的鄱阳湖枢纽工程“调枯不调洪”的运行原则进行研究和评价.

图 1 鄱阳湖地形及分区 Fig.1 Bathymetry and zones of Lake Poyang
2 鄱阳湖水龄数值模型设置 2.1 计算方案

由于鄱阳湖在洪季和枯季表现出不同的水动力特性,因此本文在计算水龄时采用分季节计算的方法,计算方案参见表 1.参照图 2,对应方案1~4,每个季节的起始日在鄱阳湖湖流主要入口(修水,赣江西、中、南支,抚河,信江,乐安河,昌江,饶河)各处均持续释放浓度为1(任意单位)的示踪剂,设置其衰减速率为0 d-1;采用2007年NECP风场(美国国家环境预报中心)统计的鄱阳湖风场(图 3);方案1~4的示踪剂标志自然条件下注入鄱阳湖的五河来水的水龄.对应方案5,参照图 4,为研究长江水倒灌过程对湖体水交换的影响,仅在湖口持续释放浓度为1(任意单位)的示踪剂,设置其衰减速率为0 d-1,此方案的示踪剂标志倒灌进入鄱阳湖的长江水的水龄.为分析五河水系进入鄱阳湖的分支流水龄,对应方案6,设置了五河水系各个河流的洪季和枯季水龄计算方案,在洪、枯两个季节的起始日在各水系持续释放浓度为1(任意单位)的示踪剂,风场仍采用2007年用NECP风场,此方案的示踪剂标志五河分支流的季节性水龄变化.鉴于鄱阳湖的水域面积变化剧烈,在方案6中为便于比较,采用定义1 %的染色剂浓度前锋所到达湖口时的水龄作为各河流水龄最大值.

表 1 计算方案 Tab.1 Calculations setup
图 2 实测五河入湖流量及出湖流量(2007春季至2008年春季, 倒灌为负值) Fig.2 Measured discharge into and out of the lake(from spring of 2007 to spring of 2008, the reverse flow is negative)
图 3 2007年3月-2008年3月鄱阳湖NCEP风场资料 Fig.3 Diagram of wind velocity vector above Lake Poyang from March of 2007 to March of 2008
图 4 实测洪季出湖流量(2007年湖口日流量,倒灌为负值) Fig.4 Measured discharge during flood season (Hukou daily discharge in 2007, the reverse flow is negative)

模型在水平向采用网格数为185×265的正交曲线网格,其中水域计算网格有9385个,网格空间步长采用52~655 m.动边界的处理根据两级水深阈值采用3种单元:完全湿单元适用于本单元及其相邻单元的水深均大于预定的第二级水深阈值(此处为0.05 m);完全干单元适用于本单元水深小于预定的第一级阈值(此处为0.01 m),且单元所有边界均为陆地边界,不进行任何水力计算;半干半湿单元适用于当本单元水深大于第一级水深阈值但小于第二级阈值,或本单元水深小于第一级水深阈值但单元的边界不全是陆地边界,洪水漫滩过程按照动边界处理采用半干半湿单元.垂向分为4层,每层的厚度由湖底地形以及水面高程来决定.模型主要受径流、表层风应力的驱动.鄱阳湖具有浅水湖泊的特点,平面尺度相对较大,因此采用垂线平均的染色剂浓度来计算水龄.风场对鄱阳湖的湖流影响主要是发生在南北方向上,即与鄱阳湖水域的走向一致,因而在季风影响上NECP再分析风场与实际的风场相似.在极端气象条件下或湖域非均匀分布风场情况下,实际风场与NECP风场是有差异的,产生差异的主要原因是NECP是再分析风场资料且网格尺度偏大.考虑到基于2007年NECP风场资料风场对湖流流速的贡献在0.1 m/s以内,此外受到缺乏实测资料的限制,这种与实际风场的差异及修正不在本文的研究范围之内.

2.2 数值模型的率定和验证

采用星子站水位和流速对模型进行精度的验证,因缺乏全年的同步水位流速实测资料,水位验证采用了2007春季至2008年春季的水位资料.江西省鄱阳湖水文局于2011年11月13-18日和2012年5月11-14日以及22-24日在鄱阳湖和长江实施了3个航次的水文测验[17],为模型的流速验证提供了可靠的数据.进行水位和流速验证前,通常需要调节底部粗糙厚度(Z0).因此对底部粗糙厚度进行了敏感性分析与测试,选定Z0为0.019 m(通常情况下水动力模型中的设定值).结果表明,星子站实测水位和模型计算水位数据吻合较好(图 5).对水位进行误差分析,模型的精度由平均绝对偏差(MAE)和均方根误差(RMSE)来衡量:

图 5 星子站水位验证(2007春季至2008年春季,85黄海基面) Fig.5 Verification of water level in Xingzi Station(from spring of 2007 to spring of 2008, National 85 Yellow Sea Datum)
$ MAE=\frac{\sum |{{X}_{\rm{obs}}}-{{X}_{\rm{sim}}}|}{N} $ (1)
$ RMSE=\sqrt{\frac{\sum {{({{X}_{\rm{obs}}}-{{X}_{\rm{sim}}})}^{2}}}{N}} $ (2)

式中,XobsXsim分别为实测和模型计算的数值.

通过MAERMSE可以直观地反映模型精度.偏差最大时间出现在4月,其余的时间模型偏差很小,总的来说MAERMSE均小于0.4 m(相对误差小于5 %),表明模型的精度满足要求(表 2).将星子站的实测垂向流速分布和本模型分层计算的垂向流速分布做对比并进行误差分析,可以看出计算的流速分布和实测值吻合(图 6);从表 3可以看出计算值与实测值吻合较好(相对误差均小于5 %).同时不同时刻垂向流速分布曲线的表层和底层差异不大,表明鄱阳湖湖流具有典型的浅水湖泊特征,因此下文以染色剂垂线平均分布来表示三维染色剂和水龄计算结果输出.

表 2 星子站实测水位与模拟水位MAERMSE Tab.2 MAE and RMSE between measured and simulated water level in Xingzi Station
图 6 星子站垂向流速分布验证 Fig.6 Verifications of vertical velocity distribution in Xingzi Station
表 3 星子站垂向平均流速和流向的实测值与模拟值比较 Tab.3 Comparisons between measured and simulated vertical averaged velocity and direction in Xingzi Station
3 鄱阳湖水龄数值模型试验

染色剂是指可以溶解于水中成为分子状态的着色剂,只能以溶解态存在于水中,通过水流来输运,只与水动力有关.在湖泊研究中,染色剂可以作为所有溶解于水且不吸附于泥沙上的物质(比如部分污染物)的代称,它通过各种物理、生物、化学过程发生衰减,其衰减的快慢通过衰减速率来确定.染色剂模型是水龄模型的基础,水龄被定义为在一个复杂的水体里的时间和空间上的变量.结合鄱阳湖的实际情况,基于环境水动力学模型EFDC源程序,建立了染色剂模型和水龄模型并纳入EFDC的三维水动力模型,通过水龄的计算,可以定量地评价五河水系水体的运动和交换程度以及滞留情况,对评估鄱阳湖水体的交换以及水质的变化具有重要意义.

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

本文采用鄱阳湖2007年春季至2008年冬季流量、水位、风场以及地形数据,建立EFDC三维水动力模型及染色剂模型.染色剂模型建立在水动力模型的基础上,水动力模型经过准确计算并率定分析后再建立染色剂模型.由模拟结果(图 7)可知,湖流流速的变化主要受径流影响,受风场的影响较小.春季星子站的流速变化范围为0.28~0.71 m/s;夏季星子站的流速变化范围为0.40~0.78 m/s(倒灌期间流速在0~0.50 m/s);秋季星子站的流速变化范围为0.26~0.72 m/s;冬季星子站的流速变化范围为0.19~0.69 m/s;夏季的流速变幅最小,而秋季流速变幅最大.湖流方向随径流的改变而变化,染色剂从五河水系入口处放入,随着湖流方向开始扩散.示踪剂的分布反映了注入鄱阳湖的五河来源的水体季节性变化,总体上,东北部水体交换缓慢,西南部水体交换迅速.浅水湖泊的湖流常常表现为急变流(弗劳德数大于1),譬如:鄱阳湖行洪漫滩过程或拟建枢纽工程水闸运行时的闸下出流,急变流模拟在数值上需要解决的关键问题是控制亚网格伪物理耗散效应.此外,急变流条件下染色剂的分布与缓变流有较大差别,受到表面波的影响.目前仍较多地依赖物理模型试验和现场实际测量.

图 7 季风影响下入湖染色剂和垂向平均流速流场分布 Fig.7 Dye and averaged vertical velocity distributions under monsoon
3.2 鄱阳湖水龄的季节性变化

基于染色剂模型对水龄进行计算,结果(图 8)表明,春季出湖口的水龄为45 d,夏季出湖口的水龄为35 d(在东北湖湾处的水龄为65 d),秋季出湖口的水龄为27 d,冬季出湖口的水龄为52 d.四季出湖水龄从小到大依次为秋季、夏季、春季、冬季,湖心区的水龄又小于东北湖湾处的水龄.水龄越大,水体交换的速率越慢.从湖区水龄的空间分布来看,湖心区的水龄远远小于东北湖湾处水龄,表明自然状况下东北湖湾处水体交换极为缓慢.水龄分布显示了鄱阳湖过水型湖泊的特点,鄱阳湖主湖区水体交换受季节性来水影响明显,夏、秋季的水龄相对小,在多数年份又受到长江倒灌的影响使水龄有所增大,总的来说夏、秋季水循环较好且水质由于长江水倒灌得到改善.冬、春季鄱阳湖湖流显著减弱,相较于夏、秋季水龄较大,尤其是冬季水龄最大,表明冬季水循环在自然条件下最为缓慢且水质较差.鄱阳湖夏季盛行偏南季风,冬季盛行偏北季风,其他两季存在过度期.鄱阳湖季节性水龄变化与季风也有相关性,在夏季偏南向季风的作用下,南北向的水域地形保持了较大的风区和风生流的量值,并且风生流的方向与入湖径流方向一致,偏南季风有减少水龄的趋势.与夏季偏南季风的作用相反,在冬季偏北季风作用下风生流的方向与入湖径流方向相反,冬季偏北季风有增大水龄的趋势.季风对水龄的影响与模型的计算结果相一致,由此得出季风加剧了鄱阳湖水龄季节性差异的结论.相较于洪季(4-9月),枯季(10月至次年3月)的水位较低、流速小,存在水域面积不足的问题.拟建的鄱阳湖大孔径水闸通过调节水位能够有效地调度湖区水量,是改善枯季水域面积不足的有效措施.此外,染色剂作为标志水团来源及其运动位置的一个指标,较为精确地反映了不同来源的水团在紊动水流作用下的对流和扩散过程.

图 8 季风影响下鄱阳湖水龄的季节分布 Fig.8 Water age distribution in Lake Poyang under monsoon
4 长江水倒灌对鄱阳湖水龄的影响

为研究长江水倒灌对鄱阳湖的影响[18-20],不仅要关注倒灌的天数更要关注倒灌水量的大小.在1951-2007年连续57个统计年中共发生倒灌48次,累计率为50 %的倒灌水量是25亿m3,2007年的倒灌水量是45.2亿m3,处于累计率约80 %.典型年宜选择倒灌量中等偏大的年份,因此以2007年作为典型年是合理的.根据历年的统计分析倒灌主要发生在洪季,但在1974年10月、1983年10月和1996年11月也都发生过倒灌,2007年鄱阳湖共发生两次倒灌(表 4),因为两次倒灌的时间相隔较近,因此本文选取两次倒灌连续计算模拟(22 d),在湖口(长江水进入鄱阳湖)释放任意单位为1的染色剂,染色剂不为0的地方即为长江水到达的地方.模拟结果(图 9)表明,2007-2008年的枯季,湖流方向从南向北;在2007年的汛期倒灌期间,湖口的湖流方向从北向南.这个倒灌过程的模拟还原了长江对鄱阳湖的淡水补充过程,水龄的模拟显示这种交换对鄱阳湖水质提升很重要,倒灌入湖的水体扩散和输运的范围越远对湖区水质的提升越有利,还为鱼类迁移和生物多样性提供了一个重要通道[21-23].洪季倒灌入湖的优质水能够在湖心对流扩散缓慢的区域保持到枯季,在鄱阳湖水利枢纽大孔闸蓄水运行的条件下,枯季的水质也能够得到提升.因此根据水位的变化设置拟建鄱阳湖枢纽工程的水闸“调枯不调洪”的调蓄原则,对于保持鄱阳湖和长江水的水体交换和湖流变化以及满足生物多样性需求有重要意义.

表 4 2007年鄱阳湖倒灌面积 Tab.4 Reverse flow area of Lake Poyang in 2007
图 9 2007年洪季倒灌期间流速染色剂和水龄分布 Fig.9 Distribution of velocity, dye and water age during reverse flow in flood season of 2007
5 鄱阳湖五河水系各河流水龄

为进一步分析水资源的分布,针对各入湖河流分别进行了水龄计算,选择适合鄱阳湖分支流水龄计算的1 %浓度前锋重新定义了水龄.各支流水龄进行计算对鄱阳湖河网与受水区水量分配和水体交换的分区有重要意义.修水洪季出湖口的水龄为32 d,枯季出湖口的水龄为55 d;赣江洪季出湖口的水龄为22 d,枯季出湖口的水龄为51 d;抚河洪季出湖口的水龄为107 d,枯季出湖口的水龄为87 d;信江洪季出湖口的水龄为82 d,枯季出湖口的水龄为69 d;饶河洪季出湖口的水龄为98 d,枯季出湖口的水龄为128 d(图 10图 11).从分支流的水龄计算可知,修水、饶河和赣江洪季水龄小于枯季水龄,而抚河、信江洪季水龄大于枯季水龄,表明在鄱阳湖水体交换上不同河流具有不同的特点,这与自然条件下的观测结果也是相吻合的,因此需要根据各支流的情况设计不同的调度运行方案.赣江的出湖水龄在洪季和枯季都相对较小,表明赣江的流量在鄱阳湖水体交换过程中起主导作用.对于西南湖区,水龄主要受到赣江的影响;对于西北部湖区,水龄受到赣江和修水的共同影响;对于南部湖区,受抚河和信江的影响较大;东部湖区则受饶河的影响较大;湖心区和入江水道的水龄则受到五河水系的综合影响.上述分区水龄的分析表明在入江水道选择建闸是合理的,各水系的水龄特征直接影响着枢纽工程水闸运行方案的制定及调度方案的具体实施,同时为各支流的水闸系统对湖区水资源进行预先调节提供了参考.

图 10 洪季各分支河流水龄分布(2007年4-9月) Fig.10 Water age distributions of individual river in flood season (from April to September of 2007)
图 11 枯季各分支河流水龄分布(2007年10月-2008年3月) Fig.11 Water age distributions of individual river in dry season (from October of 2007 to March of 2008)
6 结论

本文建立了基于EFDC的三维水动力水龄模型来研究鄱阳湖水龄的季节性变化以及倒灌过程对其的影响.鄱阳湖的水龄模型既反映了来源于五河水系的径流在湖区的更新速率,也反映了长江水倒灌入湖水体与湖体的交换.鄱阳湖的水体交换与更新速率主要受到季节性的流量变化影响,风场等因素对水龄的影响较小.对于2007-2008年来源于五河水系的水体,春季的出湖水龄约为45 d,夏季的出湖水龄约为35 d,秋季的出湖水龄约为27 d,冬季的出湖水龄约为52 d;倒灌入湖的长江水水体的水龄约为32 d,水龄计算显示夏、秋季的水龄小于冬、春季,表明冬、春季的水体交换速率相对缓慢,而夏、秋季水体交换速率相对较快.通常发生在夏、秋季的长江水倒灌过程使水龄增大,染色剂和水龄模型结果表明长江水可以进入湖心区并在湖区停留较长时间,相应地鄱阳湖湖区水质也因长江水的注入而改善.对于拟建中的鄱阳湖水利枢纽工程,其“调枯不调洪”的运行原则是基本符合鄱阳湖水体交换的自然规律的,能够改善湖区的枯季蓄水量和调节湖区的水体交换.分支流的水龄计算表明:赣江的洪季出湖水龄最小,对出湖水龄的主要影响来自赣江.西北部湖区水体交换主要受到赣江和修水的影响;西南部湖区水体交换主要受赣江的影响;南部湖区水体交换受抚河和信江的影响较大;东部湖区水体交换则受饶河的影响较大;湖心区和入江水道的水体交换则受到五河水系的综合影响.因此,在入江水道上建设水利枢纽工程是合理的选择,利用入湖河流的水闸系统进行湖区蓄水量的预调节将有利于鄱阳湖枢纽调度方案的制定和实施.此外,鄱阳湖水利枢纽的建设无疑会增大枯季的水龄,即减缓水体交换的速率,因此选择保持一定的枯季生态流量,在枯季选择持续的蓄水或间断的蓄水是制定合理的枢纽工程水闸运行规则的重要方面,而水龄的计算为鄱阳湖枢纽工程的可行性研究提供了重要的参照.水龄的分布和保守示踪剂的浓度计算针对水团的对流-耗散过程,同时为反应-耗散过程提供了泥沙水质等环境变量的背景浓度场,因此水龄的概念虽然源自水动力计算同时也是水质和生态环境变量的重要研究基础.实测资料的分析表明2007年是枯水年,从理论发展和实际应用的角度,结合不同水文情势分析不同代表年情景下水龄季节变化的变幅及差异也是今后的重要研究方向.此外,结合物模试验成果和现场观测资料,进一步开展高分辨率的急变流模式研究,是解决拟建鄱阳湖枢纽工程水闸运行方式及其影响的关键,也是今后浅水湖泊湖流及水体交换混合的重要研究方向.

致谢: 感谢江西省水利科学研究院的各位同仁在野外采样工作中给予的无私帮助,感谢论文完成过程中南昌大学胡振鹏教授、江西省水利厅谭国良高工给予的宝贵建议.
7 参考文献

[1]
Deleersnijder E, Campin JM, Delhez EJM. The concept of age in marine modelling Ⅰ. Theory and preliminary model results. Journal of Marine Systems, 2001, 28(3/4): 229-267.
[2]
Shen J, Wang HV. Determining the age of water and long-term transport timescale of the Chesapeake Bay. Estuarine, Coastal and Shelf Science, 2007, 74(4): 585-598. DOI:10.1016/j.ecss.2007.05.017
[3]
Shen YM, Wang JH, Zheng BH et al. Modeling study of residence time and water age in Dahuofang Reservoir in China. Science China-Physics, Mechanics & Astronomy, 2011, 54(1): 127-142.
[4]
Li YP, Acharya K, Chen D et al. Modeling water ages and thermal structure of Lake Mead under changing water levels. Lake and Reservoir Management, 2011, 26(4): 258-272.
[5]
Gong R, Xu J, Xu LG et al. Study on hydrodynamics in urban landscape lake based on EFDC. Environmental Engineering, 2015(4): 58-62. [龚然, 徐进, 徐力刚等. 基于EFDC城市景观湖泊水动力模拟研究. 环境工程, 2015(4): 58-62.]
[6]
Ribbe J, Wolff JO, Staneva J et al. Assessing water renewal time scales for marine environments from three-dimensional modelling:A case study for Hervey Bay, Australia. Environmental Modelling & Software, 2008, 23(10/11): 1217-1228.
[7]
Huang WR, Spaulding M. Modelling residence-time response to freshwater input in Apalachicola Bay, Florida, USA. Hydrological Processes, 2002, 16(15): 3051-3064. DOI:10.1002/(ISSN)1099-1085
[8]
Kennedy MG, Ahlfeld DP, Schmidt DP et al. Three dimensional modeling for estimation of hydraulic retention time in a reservoir. Journal of Environmental Engineering-ASCE, 2006, 132(9): 976-984. DOI:10.1061/(ASCE)0733-9372(2006)132:9(976)
[9]
Dabrowski T, Hartnett M, Olbert AI. Determination of flushing characteristics of the Irish Sea:A spatial approach. Computers & Geosciences, 2012, 45: 250-260.
[10]
Arega F, Armstrong S, Badr AW. Modeling of residence time in the East Scott Creek Estuary, South Carolina, USA. Journal of Hydro-environment Research, 2008, 2(2): 99-108. DOI:10.1016/j.jher.2008.07.003
[11]
Shao JR, Wu SQ, Zhou J et al. An age model for water transfer. Advances in Water Science, 2014, 25(5): 695-703. [邵军荣, 吴时强, 周杰等. 水体交换年龄模型研究. 水科学进展, 2014, 25(5): 695-703.]
[12]
Huang CL, Li X, Sun YY. Water age distribution of the Taihu Lake and impact of the Yangtze River to Taihu Lake Water Transfer Project on water age. J Lake Sci, 2017, 29(1): 22-31. [黄春琳, 李熙, 孙永远. 太湖水龄分布特征及"引江济太"工程对其的影响. 湖泊科学, 2017, 29(1): 22-31. DOI:10.18307/2017.0103]
[13]
Zhang XL, Li YL, Yu G et al. Runoff simulation and response to climate changes for Poyang Basin during the past 1000 years. J Lake Sci, 2016, 28(4): 887-898. [张小琳, 李云良, 于革等. 鄱阳湖流域过去1000a径流模拟以及对气候变化响应研究. 湖泊科学, 2016, 28(4): 887-898. DOI:10.18307/2016.0423]
[14]
Li X, Wang YG, Zhang SX. Numerical simulation of water quality in Yangtze Estuary. Water Science and Engineering, 2010, 2(4): 40-42.
[15]
Zhang SX, Mei M, Xu XF et al. Numerical modeling of water level and current in very shallow lake Poyang. 12th International conference on fuzzy systems and knowledge discovery, FSKD, 2015:2544-2548.
[16]
Qi W, Li YP, Wang Y et al. Control and optimization scheme of TianJing Lake hydrodynamics based on EFDC model. Water Resources and Power, 2015(1): 47-51. [戚文, 李一平, 王莹等. 基于EFDC模型的天镜湖水动力优化调控方案研究. 水电能源科学, 2015(1): 47-51.]
[17]
Zhang C, Min Q, Cao Y et al. The hydrology report of field investigations for verifications of similarity rule in Poyang Lake physical model. Poyanghu Hydrology Bureau of Jiangxi Province, 2012. [张纯, 闵骞, 曹美等. 鄱阳湖实体模型相似验证原型水文测验报告. 江西省鄱阳湖水文局, 2012. ]
[18]
Tang CX, Xiong X, Wu NH et al. Simulation of the impact of the reverse flow from Yangtze River on the hydrodynamic process of Lake Poyang. J Lake Sci, 2015, 27(4): 700-710. [唐昌新, 熊雄, 邬年华等. 长江倒灌对鄱阳湖水动力特征影响的数值模拟. 湖泊科学, 2015, 27(4): 700-710. DOI:10.18307/2015.0419]
[19]
Hu Q, Feng S, Guo H et al. Interactions of the Yangtze river flow and hydrologic processes of the Poyang Lake, China. Journal of Hydrology, 2007, 347(1/2): 90-100.
[20]
Guo H, Hu Q, Zhang Q et al. Effects of the Three Gorges Dam on Yangtze River flow and river interaction with Poyang Lake, China:2003-2008. Journal of Hydrology, 2012, s416/417(2): 19-27.
[21]
Xie DM, Deng HB, Wang DY et al. Zoning of ecological function importance for Lake Poyang wetlands. J Lake Sci, 2011, 23(1): 136-142. [谢冬明, 邓红兵, 王丹寅等. 鄱阳湖湿地生态功能重要性分区. 湖泊科学, 2011, 23(1): 136-142. DOI:10.18307/2011.0120]
[22]
Zhang QJ, Yu XB, Hu BH. Research on the characteristics of plant communities in the Poyang Nanji wetlands, China. Resources Science, 2013, 35(10): 42-49. [张全军, 于秀波, 胡斌华. 鄱阳湖南矶湿地植物群落分布特征研究. 资源科学, 2013, 35(10): 42-49.]
[23]
Dong YY ed. Contingent valuation of Yangtze finless porpoisesin Poyang Lake, China. Springer, 2013.