湖泊科学   2021, Vol. 33 Issue (2): 449-461.  
0

研究论文——水环境保护与修复

引用本文 [复制中英文]

李金华, 李一平, 唐春燕, 程一鑫, 施媛媛, 程月, 伍欣, “引江济星”工程及风场环境影响广东肇庆星湖水龄分布特征. 湖泊科学, 2021, 33(2): 449-461. DOI: .
[复制中文]
Li Jinhua, Li Yiping, Tang Chunyan, Cheng Yixin, Shi Yuanyuan, Cheng Yue, Wu Xin. Water age distribution characteristics of Lake Star (Zhaoqing, Guangdong) influenced by the Water Diversion Project and wind field. Journal of Lake Sciences, 2021, 33(2): 449-461. DOI: .
[复制英文]

基金项目

国家重点研发计划项目(2017YFC0405203,2016YFC0401703)、中央高校建设世界一流大学(学科)和特色发展引导专项资金项目和国家自然科学基金项目(51779072,51809102)联合资助

通信作者

李一平, E-mail:liyiping_hhu@163.com

文章历史

2020-04-21 收稿
2020-07-01 收修改稿

码上扫一扫

“引江济星”工程及风场环境影响广东肇庆星湖水龄分布特征
李金华1,2 , 李一平1,2 , 唐春燕1,2 , 程一鑫1,2 , 施媛媛1,2 , 程月1,2 , 伍欣3     
(1: 河海大学浅水湖泊综合治理与资源开发教育部重点实验室, 南京 210098)
(2: 河海大学环境学院, 南京 210098)
(3: 肇庆市环境保护监测站, 肇庆市环境科学研究所, 肇庆 526000)
摘要:目前大量景观类湖泊在人工布局下呈现形状不规整的特点,容易造成湖湾水体滞留形成死水区,水质恶化问题日益突出,本文以星湖为研究区域,建立基于环境流体生态动力学(EFDC)模型的星湖水动力-风场耦合模型,以水龄作为衡量水体交换速率的指标,设计12种工况模拟计算不同出入湖通道布局及流量配比、"引江济星"工程和风场对星湖水龄时空分布的影响.结果表明:星湖水龄存在季节性和空间性分布差异,春季主导风向与湖体流动主方向互斥,抑制了波海湖水体交换,导致水龄在春季较长,而秋季南西南季风促进中心湖等水体交换,致使水龄偏短,也反映出风向的影响存在湖区空间异质性,风向为南东南时星湖整体水龄最小;入湖口通道布局及流量配比对湖区水体交换速率存在影响,以外坑为波海湖入湖口且波海湖、中心湖、仙女湖入流比为4:3:3时,星湖整体水龄最小,拟建的青莲湖出水口和仙女湖出水口能有效加快青莲湖东北部和仙女湖东南部的水体交换速率;综合考虑调水引流效益、经济投资、生态开发和占地的成本,方案二为最优调水引流方案.本研究强调入湖口布局、流量大小和比例、及风场等因素对城市内湖水体置换的影响,为人为干扰下的景观湖泊的治理和管控提供科学支撑.
关键词星湖    水龄    EFDC模型    “引江济星”工程    风场    
Water age distribution characteristics of Lake Star (Zhaoqing, Guangdong) influenced by the Water Diversion Project and wind field
Li Jinhua1,2 , Li Yiping1,2 , Tang Chunyan1,2 , Cheng Yixin1,2 , Shi Yuanyuan1,2 , Cheng Yue1,2 , Wu Xin3     
(1: Key Laboratory of Integrated Regulation and Resource Development on Shallow Lakes, Ministry of Education, Hohai University, Nanjing 210098, P. R. China)
(2: College of Environment, Hohai University, Nanjing 210098, P. R. China)
(3: Environmental Protection Monitoring Station, Institute of Environmental Science, Zhaoqing 526000, P. R. China)
Abstract: At present, a large number of Landscape lakes have the characteristics of irregular shape under the artificial layout, which is easy to cause bay water retention and form dead zones. The water quality deterioration is increasingly prominent. This paper takes Lake Star as the research area, based on EFDC (environmental fluid dynamics code) model that hydrodynamic-wind field coupling model was set up, the water age is taken as the index to measure the water exchange rate, 12 working conditions were designed to calculate the influence of different access points and flow ratio, Water Diversion Project from the West River to Lake Star and wind field on the spatial and temporal distribution of water age in Lake Star. The result shows that:there are seasonal and spatial distribution differences in the water age of Lake Star. The dominant wind direction in spring and the main flow direction of the lake are mutually exclusive, which restrained the water exchange with Lake Bohai, resulting in a longer water age in spring. The south-southwest monsoon in autumn promotes the water exchange of Lake Central and other water bodies, resulting in a shorter water age, which also reflects the influence of wind direction on the spatial heterogeneity of the lake area. When the wind direction is southeast, the water age of Lake Star is the smallest; the channel layout and flow ratio have influence on the water exchange rate in the lake area, the outer pit is the entrance of Lake Bohai and when the inflow ratio of Lake Bohai, Lake Central and Lake Xiannv is 4:3:3, the overall water age of Lake Star is the smallest, the proposed outlets of Lake Qinglian and Lake Xiannv can effectively accelerate the water exchange rate in the northeast of Lake Xiannv and southeast of Lake Qinglian; comprehensive study considering the benefit of water diversion, economic investment, cost of ecological development and land occupation, the second scheme is the optimal water diversion scheme. This study emphasizes the influence of the layout of the lake entrance, the size and proportion of the flow, and the wind field on the water replacement of the municipal lakes, which provides scientific support for the management and control of the landscape lakes under human disturbance.
Keywords: Lake Star    water age    Water Diversion Project from the West River to Lake Star    EFDC model    wind field    

城市景观湖泊通常兼具景观生态、接纳雨污、储蓄调洪等多种功能[1],由于地处人口相对密集的区域,周边分布着复杂的小区、工厂、商户,城市景观湖泊的水质往往更容易受到人为排污等活动的影响,生态环境十分脆弱,加之此类湖泊水系连通性差,缺少补水源水,出入流河道、闸门较少,湖泊相对封闭,导致湖泊换水周期较长.城市景观湖泊在人为布局下往往呈现形状不规整的特点,容易使湖湾水体滞留形成死水区,在夏季高温时期,藻类极易滋生堆积形成水华,湖泊水体健康面临严重危机[2].星湖位于广东省肇庆市,是国家级重点风景名胜区,近年来星湖富营养化问题日益加剧,为改善部分湖区水体流通性差、水质不达标的现状,肇庆市政府提出“一滴污水不进星湖”,并对星湖周边排污口、排水沟实行严格的控源截污,同时为缩短星湖换水周期,改善星湖水质,提升星湖水体透明度,提出了“引江济星”工程,即引西江水和北岭山山水进入星湖.

对于大型湖泊而言,入湖流量和风场在驱动湖泊水体运动和水动力交换过程中发挥着重要的作用,水龄作为衡量水体交换速率的重要指标,已有大量关于水龄研究的成果. Li等[3]引入水龄的概念,分析“引江济太”工程引水量对太湖水动力调控效果的影响,研究表明太湖水体的输移过程受出入流水动力条件和风场的强烈影响;Gao等[4]研究了引水工程对人为布局下的城市人工湖水龄的影响,揭示了引水量和湖区布局与水体交换速率间的关系,并基于EFDC模型模拟了风场及水位对七里海泻湖水体更新能力的影响[5],发现东北风显著提升了七里海泻湖的水体更新能力;王钟等[6]分析“引江济淮”工程和风应力对蔡子湖水龄分布的影响,结果表明:风应力对蔡子湖水龄分布有重要影响;黄春琳等[7]基于EFDC模型研究“引江济太”工程对太湖水龄分布的影响;张素香等[8]构建鄱阳湖染色剂和水龄模型,计算不同入流河道对湖区的影响,论证了鄱阳湖水利枢纽工程调水方案的合理性;戚文等[9]利用EFDC模型建立天镜湖的三维非稳态水量、水龄、污染物颗粒追踪数学模型.此外,通过建立水龄和水质的相关关系,可便于明晰水体中营养盐含量与水龄的联系;Zhang等[10]构建三维水动力水质耦合模型,分析牛栏山引水工程对滇池水龄和水质的影响,发现滇池水龄的大小与总磷浓度呈正相关,与总氮、叶绿素a浓度无明显相关性.以上研究虽然都取得了很多有意义的成果,但大多是针对引调水工程、风场或季节变化对大型湖泊水体输移速率方面的研究,涵概因素单一,特别是针对南方季风气候区城市景观湖泊,缺少全面考虑进出水通道布局、流量配比及风场对工程效益的影响等方面的研究.

近年来,随着引调水工程在西湖[11]、太湖、滇池的成功运用,调水引流已经作为提升湖泊水体交换速率、改善湖泊水质的重要措施之一[12],本研究以南方季风气候地区典型城市景观湖泊——星湖为例,以EFDC模型为工具,基于实测风场、水位和地形数据以及基于湖泊水量平衡换算的流量数据构建星湖三维水动力—风场耦合模型,综合考虑季节变化、风场、调水引流工程对星湖水龄时空分布的影响,并结合实际情况,对进水线路和进水比例进行优化,全面考虑引调水工程的效益、经济投资、生态占地和生态开发等各个方面,遴选最优引调水方案、最佳调水季节和风向,从水动力的角度模拟星湖在不同调水引流情景下以及不同风向下星湖水龄的时空分布情况,以提供星湖平均水龄较低的最优调水引流方案,为“引江济星”工程的实施管理提供科学依据.

1 研究区域及方法 1.1 研究区域概况

星湖(23°3′27″~23°5′15″N, 112°26′38″~112°29′20″E)位于广东省肇庆七星岩风景区内,北倚北岭山,南接端州城区,地处南亚热带,年均温度21.2℃,气候温暖湿润,属于季风气候,且降水主要集中在4-9月,年平均降水量约1650 mm,2018年全年主导风向为西南风,年平均风速约1.25 m/s,星湖为城市内景观湖泊,平均水深为2.34 m,最大水深为3.65 m,总面积约为5.82×106 m2(相应水位为4.91 m),总容积约为13×106 m3,目前星湖水生植物较少,常见水生植物约15种,其中沉水植物仅有4种,群落结构简单[13].

星湖由波海湖、中心湖、青莲湖、仙女湖、里湖五大湖区组成(图 1),其中波海湖、仙女湖、中心湖和青莲湖通过桥涵互相连通,水可以自由流动,补水源水由北部及周边集水区的涝水和北岭截洪渠的山洪溢流提供,里湖相对封闭,主要由湖面降雨和中心湖补给水量,且与其他湖区存在约1 m的水位落差,本文以波海湖、中心湖、青莲湖、仙女湖4个子湖区为主要研究对象,由于自然条件下里湖与其它4个子湖区互不存在水体交换,故暂不考虑引水工程及风场对里湖水龄分布的影响.星湖的主要入湖口为外坑A1、盘古坑A2、石牌坑A3以及概化的雨水口B1、排污口B2、排污口B3,主要出水口为中心湖南部的牌坊换水涵C1、青莲湖东南部的芹田水闸C2以及仙女湖北部的东岗节制闸C5,青莲湖出水口C3和仙女湖出水口C4为规划拟建的出水口.

图 1 研究区域概况 Fig.1 Map of the study area
1.2 模型的构建及率定验证 1.2.1 模型构建

EFDC模型最早由美国美国弗吉尼亚大学和海洋科学研究所开发,应用笛卡尔坐标或曲线正交坐标系统[14],在空间上运用二阶精度有限差分格式,水平扩散方程采用时间显示、空间隐式格式,水平输运方程采用Blumberg-Mellor模型的中心差分格式或者正定迎风差分格式[15].被广泛应用于河流[16]、湖泊[17]、湿地[18]、河口[19]以及水库[20]的水动力、水质和水生态模拟中,能够动态展现三维水动力及污染物迁移扩散的轨迹.本文将EFDC模型应用于星湖的水动力模拟研究中,构建三维水动力模型,在水平上采用笛卡尔直角坐标系,垂向上采用Sigma坐标,水平网格数共计2212个,网格空间分辨率为52 m,垂向分为3层,每层厚度由水面高度和湖底地形决定. EFDC模型以风速、风向、入湖流量、出湖水位为水动力模块的边界条件,其中入湖边界为3条排水坑(外坑、盘古坑、石牌坑)、6个排污口和1个雨水口,出湖边界包括3个出水闸门(牌坊换水涵、芹田水闸和东岗节制闸).

由4215个实测点位高程数据内插形成星湖底部地形,初始水位为2018年1月1日星湖水位监测站实测水位(4.88 m),四季风速(图 2)数据采用2018年实测值(数据来源:http://data.cma.cn/),模拟时间2018年1月1日—12月31日,动态时间步长取0.5 s.为了适应水位波动,尤其是相对浅水区域,模型设置临界干水深0.05 m.

图 2 星湖2018年风玫瑰图 Fig.2 Wind rose of Lake Star in 2018
1.2.2 模型率定验证

模型中的参数取值会对模拟结果产生重要的影响,本研究重点率定风场参数及底部粗糙系数,最终取糙率参数为0.02,风拖曳系数为3×10-3,风遮挡系数为1.通过对星湖水位监测点位2018年全年的模拟值与实测值比较可知(图 3),星湖水位的绝对误差为0.007 m,均方根误差为4.4 %,选取2017年星湖水位监测站实测水位及基于水位平衡换算的出入流量数据作为模型的验证数据,模拟值与实测值比较可知(图 3),2017年星湖水位绝对误差为0.024 m,均方根误差为3.4 %,模拟值与实测值误差较小且吻合较好,说明模型的模拟结果能较好地反映星湖的水位变化情况.

图 3 星湖水动力模型率定(2018年)和验证(2017年)结果(珠基高程) Fig.3 Calibration(2018) and verification(2017) results of the hydrodynamic model of Lake Star (Pearl River elevation reference)
1.2.3 水龄的定义

本研究利用水龄和拉格朗日颗粒物追踪的概念来描述湖体交换快慢及交换程度,获取最佳调水流量、调水时机,提出多种改善湖体水龄、加速水体交换的有效组合.

水龄模型是基于对流扩散模型发展起来,用于衡量湖泊和水库污染物迁移的有效参数[21].水龄定义为:对单个粒子来说,粒子从进入水体到达指定位置所需的时间.针对不同的研究问题,水龄可以定义为不同的形式.例如研究水体垂向交换时,水龄定义为水体离开水体表面,到达垂向某一位置所需的时间;当研究水体与外部水体交换时可定义为颗粒物从入口传输到指定点的时间(往往入口的水龄设为零).水龄能够反映出水质点或者污染物从边界传输到研究区域内任何一点的时间,定量反映出水体交换能力的时间和空间异质性,定性反映污染物传输特征和水体富营养化程度.水龄越小,说明水体交换程度越强,反之亦然.水龄计算公式如下:

$ \begin{array}{c} \frac{\partial c(t, \vec{x})}{\partial t}+\nabla(u \cdot c(t, \vec{x})-K \nabla c(t, \vec{x}))=0 \end{array} $ (1)
$ \begin{array}{c} \frac{\partial \alpha(t, \vec{x})}{\partial t}+\nabla(u \cdot \alpha(t, \vec{x})-K \nabla \alpha(t, \vec{x}))=c(t, \vec{x}) \end{array} $ (2)

式中,c为示踪剂浓度(mg/L),α为水龄初始设置值(d),u为流速(m/s),K为扩散张量,t为时间(s),$\overrightarrow x $为坐标.水龄a可根据下式计算:

$ a(t, \vec{x})=\frac{\alpha(t, \vec{x})}{c(t, \vec{x})} $ (3)
1.2.4 计算方案

为确定入湖口通道布局、入湖流量分配比例、出入湖口数量以及风场对星湖水龄的影响,本文共设计4个方案共计12种工况(表 1)进行计算.星湖是典型的城景观湖泊,湖区面积较小,工况1和工况10分别计算现状情景下有风和无风时星湖的水龄分布,并基于实测气象数据模拟星湖水龄的季节性变化规律,方案一、方案二和方案三分别计算了不同引调水方案下星湖的水龄空间分布,其中工况3、4、5、7、8、9计算不同入湖口位置分布和入湖流量比例情景下星湖水龄空间分布,工况6模拟拟建的青莲湖出水口和仙女湖出水口对星湖水龄的影响,工况61~16则模拟恒定风速,16种不同风向下星湖水龄的时空分布特征.

表 1 模型计算方案* Tab. 1 Model calculation scenarios
2 结果与讨论 2.1 星湖水龄时空分布特征及风场对水龄分布的影响

现状情景下工况1和工况10分别计算了有风和无风情景下星湖水龄的季节性分布情况(图 4),在无风情景下水龄的季节性变化从小到大依次为:秋季(123.95 d)<夏季(131.89 d)<冬季(155.21 d)<春季(182.23 d),而有风情景下四季水龄从小大依次为:秋季(96.63 d)<夏季(114.78 d)<冬季(142.42 d)<春季(179.77 d),较无风情景而言,有风情景下星湖春、夏、秋、冬季整体平均水龄分别降低了2.46、17.11、27.31、12.78 d,主要原因是秋季入流流量大,且秋季主导风向与湖流运动方向相近,利于水体交换,而春、冬季节出入流流量小,水体交换速率较弱;由于湖区边界形状差异,水龄的季节变化在湖区空间分布上也有显著区别,对比无风情景,有风情景下中心湖、青莲湖、仙女湖水龄分别降低了52.37、11.65、43.95 d,波海湖水龄反而增加了21.6 d,主要是因为波海湖南、西、北三面为岸线形成“小湖湾”,A1、A2两个入流口与风向相悖,唯一的出流口仅为波海湖与中心湖中间的相通桥涵,因此主导风抑制了波海湖的水体交换,而中心湖、青莲湖和仙女湖的水体交换受到风场正向驱动.夏、秋季风(西南风)作用提升了星湖的水体交换速率,而春季季风(南西南风)对星湖的水龄分布影响较小.

图 4 星湖水龄季节性差异(a)和空间分布差异(b) Fig.4 Seasonal (a) and spatial (b) variations of water age in Lake Star

对于大型浅水湖泊,风应力是湖泊形成环流的主要驱动力之一,分析不同风向下湖泊水龄的时空分布特征可以更好地论证调水引流工程对湖泊水体交换速率的改善效益,同时,水龄的时空分布也可以作为评估调水引流对湖泊水质改善效益的重要参考指标.

研究中模拟了恒定风速2 m/s时16种不同风向下、以及强风向(南东南风)下不同风速时星湖湖体水龄的时空分布,计算了波海湖、中心湖、青莲湖以及仙女湖在不同风向下的水龄,计算结果见图 5.模拟结果表明:在相同风速、不同风向下,南东南风对星湖整体水龄的改善最为显著(湖区整体平均水龄为19.95 d),北东北风对水体交换的影响最弱(湖区整体平均水龄为27.76 d),不同子湖区的影响结果也有很大的差异,青莲湖是水龄最大的湖区,仙女湖是受风场影响最小的湖区(16个风向下水龄在13.52~21.54 d之间),这表明靠近引水点区域的水体最先被交换,其它区域在不同风向情景下呈现不同的差异,波海湖平均水龄范围为6.53 d(南风)~21.54 d(西西北风);中心湖平均水龄范围为16.89 d(南西南风)~27.93 d(东北风);青莲湖平均水龄范围为39.86 d(西风)~63.46 d(北风);仙女湖平均水龄范围为13.52 d(西南风)~21.54 d(北风).

图 5 各湖区水龄受风向影响情况 Fig.5 Spatial distribution of water age in each lake district

在南东南风(强风向)0、2、4、6 m/s风速下,星湖平均水龄分别为18.44、19.95、18.18和18.35 d,对于同一湖区,由风向差异引起的水龄的最大差值为34.38 d(波海湖),最小为8.03 d(仙女湖),不同湖区在对同一风向(南东南风)不同风速下水龄的表现也有所差异,最小为2.50 d(波海湖),最大为26.31 d(青莲湖);不同风向对星湖整体水龄的影响在0.02~7.81 d之间,因此风场对星湖水龄时空分布影响较小.研究结果表明:不同风向对星湖各个湖区的改善效果不同,南风有助于波海湖的水体交换,南西南风有助于中心湖的水体交换,西西南风有助于青莲湖的水体交换,东风有助于仙女湖的水体交换,而对于星湖整体而言,风向为南东南风时水龄较小,表明南东南风是“引江济星”工程提升星湖换水速率和改善水质最有效的风向.

2.2 不同出入湖点位设置及入流配比情景下水龄分布特征

为了进一步明晰不同入湖通道布局以及配水比例下“引江济星”工程对星湖水龄的影响,结合星湖现有进水通道和“引江济星”工程拟规划设计3个进水口,对于波海湖入湖口的选择,一是由外坑直接入波海湖,二是新建输水管道后由盘古坑入波海湖,其入湖点需进一步评估,其余入中心湖和仙女湖的点位相对明确;补水进入星湖后,通过芹田水闸、牌坊换水涵和东岗节制闸3处向外排水,星湖形成“三进三出”的水体交换格局,让星湖各个湖区的水体得到充分交换,减少死水区面积,通过增强水体循环来改善水质,由模拟结果可知,经外坑入波海湖的配水点位设置要略优于经盘古坑入湖的路线,其中波海湖的改善效果尤为显著,波海湖水龄由68 d降至43 d,水体交换速率提高了37 %,但对其它湖区影响甚微.

此外,进水口的配水比例也会对湖泊水体交换产生一定的影响,本文模拟了4种配水比例工况下“引江济星”工程对星湖水动力的影响,结果表明:工况5为最优方案,即波海湖入湖口、中心湖入湖口、仙女湖入湖口的配水比例为4 :3 :3时星湖整体水龄最低,此时波海湖、中心湖、青莲湖、仙女湖以及星湖整体平均水龄分别为42.79、33.20、73.28、28.25和37.40 d(表 2).

表 2 不同工况下各湖区水龄统计 Tab. 2 Water age statistics of each lake area under different working conditions
2.3 不同引调水方案下星湖水龄分布特征

为评估不同工况下“引江济星”工程对星湖水动力的改善效果,基于最优入湖通道布局以及最优入湖流量比例,模拟3种调水引流方案下星湖的水龄分布情况,即方案一:北岭山来水(q=0.988 m3/s),方案二:西江引水(q=6 m3/s),方案三:北岭山来水+西江引水(q=0.988 m3/s+5.8 m3/s),方案一充分考虑北岭山山洪和降雨资源,通过新修截洪渠和集水池收集水资源,并利用现有进水坑道(盘古坑和石牌坑)对星湖实施补水;方案二从西江引水,通过修建引配水管道达到合理分配引水资源、盘活星湖水体和提升星湖换水速率的目的;方案三综合考虑利用北岭山山洪、降雨资源和西江引水资源.星湖现阶段主要通过牌坊换水涵、芹田水闸和东岗节制闸向外排水,另有两处规划建设排水口,即青莲湖出水口和仙女湖出水口,本研究在评估“引江济星”工程对星湖水动力改善的过程中,假定:1)忽略风场的影响;2)各个工况除入、出湖边界条件(边界位置及入流流量)等不同之外,其它条件及参数设定均相同,模型模拟时长为365 d.

基于水动力-风场耦合模型得出的不同引调水工况下星湖的水龄空间分布(图 6~图 9),颜色为蓝色代表水龄越小,颜色为红色则说明该点水龄越大,模拟结果表明,入湖流量与水龄呈反比,即随着流量的提升,湖区整体水龄水龄越小,且随着流量的增加,4、15、30、60 d等水龄线分布范围不断增加(表 3).其中方案三为降低湖区水龄值的最优调水引流方案(工况10),方案二次之(工况5),方案一对星湖的水龄改善效果最差(工况1),三种工况下星湖水龄均值分别为30.62、33.30、106.70 d,表明“引江济星”工程可以加快星湖的水体交换速率,减小星湖水龄,特别针对波海湖、中心湖和仙女湖中北部,水体交换速率有很大的提升,但青莲湖仍存在水体交换速率较差的区域,而青莲湖东部恰是藻类暴发的高发区.

图 6 现状情景水龄时空变化 Fig.6 Current situation water age spatiotemporal variation
图 7 方案一(工况1)水龄时空变化 Fig.7 Scheme one (working condition 1) water age spatiotemporal variation
图 8 方案二(工况5)水龄时空变化 Fig.8 Scheme two (working condition 5) water age spatiotemporal variation
图 9 方案三(工况10)水龄时空变化 Fig.9 Scheme three (working condition 10) water age spatiotemporal variation
表 3 不同工况下水龄区间所占湖区面积 Tab. 3 The area covered by different working conditions of water age

为减少青莲湖东部和仙女湖东部“死水区”的分布范围,降低水华的暴发风险,“引江济星”工程规划拟建青莲湖出水口和仙女湖出水口,其建设的效果和必要性仍需进一步评估,工况6模拟了新增青莲湖出水口和仙女湖出水口对青莲湖和仙女湖的水龄改善,结果表明,工况6情景下青莲湖和仙女湖水龄分别为30.45和12.92 d,工况5情景下青莲湖和仙女湖水龄分别为73.70和33.30 d,较工况5而言,工况6情景下青莲湖和仙女湖水体交换速率分别提升了58.68 %和61.20 %,说明青莲湖出水口和仙女湖出水口能有效增加青莲湖和仙女湖的水体交换速率.

2.4 调水引流工程及湖泊水龄改善的讨论

基于EFDC构建的三维水动力模型对“引江济星”工程的评估是以水龄为定量评价指标,同时考虑风场对湖流运动的影响,能够准确反映不同引调水方案实施后对星湖水体交换速率的提升效果,为工程方案选取、实施及管理提供科学依据,并已成功应用到蔡子湖[6]、鄱阳湖[8]、滇池[10]、太湖[22]等引调水工程评估中,但由于湖泊类型不同,所受风场的影响也有所差异,导致各研究结果均有所差异.太湖是大型浅水湖泊,其出流入边界明确,湖边带相对平滑,风生流在湖流运动中起着至关重要的作用,在相同点位由风向引起的水龄差异大于100 d,在不同点位的水龄差异则超过150 d[22].星湖是小型景观湖泊,其在人为布局的影响下呈现形状不规整的特点,各湖区水体连通性较弱,在湖湾地区易发生水体滞留形成死水区,其水动力交换过程主要受到入湖通道布局(包括流量配比、点位设置、出入湖道数量)的影响,受风场影响较小,其由风向引起的水龄最大差值为34.38 d,由风向引起的水龄最小差值为26.18 d.引调水工程对湖泊的影响表现在水质、水动力、透明度、水生生物群落等各个方面,为得出更为科学客观的引调水方案,水龄值不能作为选定调水引流工程的唯一判据,引调水工程的实施需综合考虑以上指标并结合当地水资源的配置、引水工程的经济投资[23]、对生态环境的开发及影响[24]等进行综合评定.在论证调水引流方案可达性、合理性、可行性的过程中,本文以减小水龄为主要目标的同时,综合了调水引流方案在实施过程中的经济投资、生态占地和开发,最终拟定方案二为最优调水引流方案,方案二(西江引水)和方案三(西江引水+北岭山来水)情景下星湖平均水龄分别为33.30和30.62 d,方案三水龄略小于方案二,但方案三需新建更多的箱涵、集水池、输水管道,其中涉及对北岭山森林资源的开发,相较于方案二,方案三需多投资近四千万(其中不含征地、拆迁费用).

3 结论

“引江济星”工程对星湖的影响可通过水龄的时空分布来体现,本文基于EFDC构建了星湖三维水动力—风场耦合模型,研究了星湖水龄的季节性分布特征、不同入湖通道布局和不同流量配比下星湖的水龄时空分布及“引江济星”工程和风场对湖泊水龄时空分布的影响.模拟结果表明,“引江济星”工程对星湖水龄时空分布有十分重要的影响,风场对星湖水龄的时空分布影响较小.

1) 星湖水龄分布存在季节和空间差异.季风作用有效地降低了星湖平均水龄,秋季季风对水体交换影响最显著;风场促进了中心湖、青莲湖、仙女湖的水体交换,抑制了波海湖的水体交换,对于星湖整体而言,由风向为南东南风时水龄最小.

2) 波海湖的最优入湖线路河道为外坑,且波海湖入流流量:中心湖入流流量:仙女湖入流流量最优分配比为4 :3 :3.综合比较“引江济星”工程的不同方案的平均水龄,结合考虑调水引流综合效益、经济投资、生态开发和和占地等成本,方案二(工况5,引西江水由外坑、雨水口和石牌坑入星湖)为最优调水引流方案.规划拟建的青莲湖出水口和仙女湖出水口能有效加快青莲湖东北部和仙女湖东南部的水体交换速率.

3) 针对形状不规则的城市内湖,在湖心大面积敞风区,风场是驱动湖体环流的主导因素,在岸线复杂的湖湾区,出入湖口布局及风场风向的一致性是影响湖区水龄的主控因素.

4 参考文献

[1]
Gong R, Xu J, Xu LG et al. Study on hydrodynamics in urban landscape lake based on EFDC. Environmental Engineering, 2015, 33(4): 58-62, 91. [龚然, 徐进, 徐力刚等. 基于EFDC城市景观湖泊水动力模拟研究. 环境工程, 2015, 33(4): 58-62, 91. DOI:10.13205/j.hjgc.201504013]
[2]
Wu X, He J. The application of EFDC model in the design of ecological water supplement project for urban landscape water-case study of thegreenlake in Kunming. Environmental Protection Science, 2019, 45(1): 73-78. [吴雪, 何佳. EFDC模型在城市景观水体生态补水工程设计中的应用——以昆明市翠湖为例. 环境保护科学, 2019, 45(1): 73-78.]
[3]
Li YP, Acharya K, Yu ZB. Modeling impacts of Yangtze River water transfer on water ages in Lake Taihu, China. Ecological Engineering, 2011, 37(2): 325-334. DOI:10.1016/j.ecoleng.2010.11.024
[4]
Gao XP, Xu LP, Zhang C. Estimating renewal timescales with residence time and connectivity in an urban man-made lake in China. Environmental Science and Pollution Research, 2016, 23(14): 13973-13983. DOI:10.1007/s11356-016-6569-5
[5]
Gao XP, Chen YY, Zhang C. Water renewal timescales in an ecological reconstructed lagoon in China. Journal of Hydroinformatics, 2013, 15(3): 991-1001. DOI:10.2166/hydro.2013.136
[6]
Wang Z, Fan ZY, Yang ZY et al. Effects of water diversion project from the Yangtze River to Huaihe river on the water age distribution of Lake Caizi, Anhui Province. J Lake Sci, 2018, 30(6): 1576-1586. [王钟, 范中亚, 杨忠勇等. "引江济淮"工程对安徽菜子湖水龄分布的影响. 湖泊科学, 2018, 30(6): 1576-1586. DOI:10.18307/2018.0609]
[7]
Huang CL, Li X, Sun YY. Water age distribution of the Lake Taihu and impact of the Yangtze River to Lake Taihu Water Transfer Project on the water age. J Lake Sci, 2017, 29(1): 22-31. [黄春琳, 李熙, 孙永远. 太湖水龄分布特征及"引江济太"工程对其的影响. 湖泊科学, 2017, 29(1): 22-31. DOI:10.18307/2017.0103]
[8]
Zhang SX, Xu XF, Xu X et al. Seasonal variations of water age in Lake Poyang. J Lake Sci, 2018, 30(1): 199-210. [张素香, 许新发, 徐兴等. 鄱阳湖水龄季节性变化特征. 湖泊科学, 2018, 30(1): 199-210. DOI:10.18307/2018.0120]
[9]
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, 33(1): 47-51. [戚文, 李一平, 王莹等. 基于EFDC模型的天镜湖水动力优化调控方案研究. 水电能源科学, 2015, 33(1): 47-51.]
[10]
Zhang XL, Zou R, Wang YL et al. Is water age a reliable indicator for evaluating water quality effectiveness of water diversion projects in eutrophic lakes?. Journal of Hydrology, 2016, 542: 281-291. DOI:10.1016/j.jhydrol.2016.09.002
[11]
You AJ, Wu ZY, Han ZC et al. Spatial and temporal distributions and variations of nutrients in the West Lake, Hang-zhou, after the implementation of integrated water management program (1985-2013). J Lake Sci, 2015, 27(3): 371-377. [尤爱菊, 吴芝瑛, 韩曾萃等. 引水等综合整治后杭州西湖氮、磷营养盐时空变化(1985-2013年). 湖泊科学, 2015, 27(3): 371-377. DOI:10.18307/2015.0301]
[12]
Gao XP, Xu LP, Zhang C. Modelling the effect of water diversion projects on renewal capacity in an urban artificial lake in China. Journal of Hydroinformatics, 2015, 17(6): 990-1002. DOI:10.2166/hydro.2015.004
[13]
Zhou LP, Huang JC, Guo CH et al. Study on the plant disposition of Xinghu national wetland park. Journal of Zhaoqing University, 2015, 36(2): 49-53. [周丽萍, 黄建昌, 郭春华等. 星湖国家湿地公园植物配置研究. 肇庆学院学报, 2015, 36(2): 49-53. DOI:10.3969/j.issn.1009-8445.2015.02.011]
[14]
Zhang YF, Wang YL, Wang L. EFDC overview and application analysis. Environmental Impact Assessment, 2015, 37(3): 70-72, 92. [张以飞, 王玉琳, 汪靓. EFDC模型概述与应用分析. 环境影响评价, 2015, 37(3): 70-72, 92. DOI:10.14068/j.ceia.2015.03.018]
[15]
Liu XM, Li JQ, Dou XM et al. The application and advance of environmental fluid dynamics code(EFDC) in estuarine water environment. Environmental Science & Technology, 2011, 34(S1): 136-140, 360. [刘夏明, 李俊清, 豆小敏等. EFDC模型在河口水环境模拟中的应用及进展. 环境科学与技术, 2011, 34(S1): 136-140, 360.]
[16]
Song WW, Xu Q, Fu XQ et al. EFDC simulation of fishway in the diversion dahaerteng river to Danghe reservoir, China. Ecological Indicators, 2019, 102: 704-715. DOI:10.1016/j.ecolind.2019.03.025
[17]
Li YP, Tang CY, Wang C et al. Assessing and modeling impacts of different inter-basin water transfer routes on Lake Taihu and the Yangtze River, China. Ecological Engineering, 2013, 60: 399-413. DOI:10.1016/j.ecoleng.2013.09.067
[18]
Jiang CB, Zhang MW, Yang XL. Water quality evaluation for the Hengshui Lake wetland in Northern China. Journal of Tsinghua University:Science and Technology, 2010, 50(6): 848-851. [江春波, 张明武, 杨晓蕾. 华北衡水湖湿地的水质评价. 清华大学学报:自然科学版, 2010, 50(6): 848-851.]
[19]
Dai J, Liu S, Han JF et al. Study on river water quality response relationship under multi scenario change of pollution load. China Environmental Science, 2018, 38(2): 776-783. [戴君, 刘硕, 韩金凤等. 污染负荷多情景变化下河流水质响应关系研究. 中国环境科学, 2018, 38(2): 776-783. DOI:10.3969/j.issn.1000-6923.2018.02.043]
[20]
Yang Q. Water environment dynamics simulation research and emergency water pollution incidents based on EFDC-Miyun Resorvior as an example[Dissertation]. Beijing: China University of Geosciences(Beijing), 2014. [杨倩.基于EFDC的密云水库水环境及应急处理模型研究[学位论文].北京: 中国地质大学(北京), 2014. ]
[21]
Deleersnijder E, Campin JM, Delhez EJM. The concept of age in marine modelling:I. Theory and preliminary model results. Journal of Marine Systems, 2001, 28(3/4): 229-267. DOI:10.1016/S0924-7963(01)00026-4
[22]
Xu RC, Pang Y, Hu ZB et al. Relationship between water age and hydraulic residence time in Taihu Lake and parameter sensitivity. Water Resources Protection, 2020, 36(3): 34-39. [胥瑞晨, 逄勇, 胡祉冰等. 太湖水龄与水力停留时间关系及参数敏感性. 水资源保护, 2020, 36(3): 34-39.]
[23]
Hao WB, Tang CY, Hua L et al. Effects of water diversion from Yangtze River to Taihu Lake on hydrodynamic regulation of Taihu Lake. Journal of Hohai University:Natural Sciences, 2012, 40(2): 129-133. [郝文彬, 唐春燕, 滑磊等. 引江济太调水工程对太湖水动力的调控效果. 河海大学学报:自然科学版, 2012, 40(2): 129-133.]
[24]
Wang ZQ. Project of water diversion to north of Zhejiang province[Dissertation]. Hangzhou: Zhejiang University, 2006. [王贞琴.浙北引水工程方案选优综合评价[学位论文].杭州: 浙江大学, 2006. ]
[25]
Lu HM, Zou Y, Feng HL. Analysis and enlightenment of eco-environmental impact of typical water diversion and transfer projects at home and abroad. Water Resources Planning and Design, 2018(12): 88-92, 166. [陆海明, 邹鹰, 丰华丽. 国内外典型引调水工程生态环境影响分析及启示. 水利规划与设计, 2018(12): 88-92, 166.]