基于MIKE21数值模拟的滆湖换水能力评估分析研究.pdf

Journal of Water Resources Research 水资源研究水资源研究,2023,12(4),334-346 Published Online August 2023 in Hans.https:/www.hanspub.org/journal/jwrr https:/doi.org/10.12677/jwrr.2023.124038 文章引用文章引用:李承龙,郭生练,高勋,冯智敏,詹慧婕.基于 MIKE21 数值模拟的滆湖换水能力评估分析研究J.水资源研究,2023,12(4):334-346.DOI:10.12677/jwrr.2023.124038 基于基于MIKE21数值模拟的滆湖换水能力评估数值模拟的滆湖换水能力评估 分析研究分析研究 李承龙李承龙1，郭生练，郭生练1，高，高 勋勋2，冯智敏，冯智敏3，詹慧婕，詹慧婕4 1武汉大学水资源工程与调度国家重点实验室，湖北 武汉 2浙江省水利水电勘测设计院有限责任公司，浙江 杭州 3黄河水利委员会宁蒙水文水资源局，内蒙古 包头 4河海大学水文水资源学院，江苏 南京 收稿日期：2023年7月17日；录用日期：2023年8月15日；发布日期：2023年8月29日 摘摘 要要 湖泊换水能力直接影响着污染物迁移和水体自净能力，为探究不同引水流量入湖对滆湖换水能力的影响，文中湖泊换水能力直接影响着污染物迁移和水体自净能力，为探究不同引水流量入湖对滆湖换水能力的影响，文中使用使用MIKE21建立了滆湖水动力水质模型，模拟并对比分析多种工况及工况组合对湖区整体和局部的水体交换建立了滆湖水动力水质模型，模拟并对比分析多种工况及工况组合对湖区整体和局部的水体交换能力的影响；提出了使用水体更新时间和水体交换率两个指标，定量分析交叉验证水体交换能力。结果表明：能力的影响；提出了使用水体更新时间和水体交换率两个指标，定量分析交叉验证水体交换能力。结果表明：引水流量为引水流量为100、150、200、250 m3/s时整个湖泊水体交换率到达时整个湖泊水体交换率到达80%以上分别需要以上分别需要60 d、45 d、32 d、30 d。本文的网格划分方案为湖泊水动力水质模型建模提供了新思路，换水能力两个指标能够从时间和空间上综合全本文的网格划分方案为湖泊水动力水质模型建模提供了新思路，换水能力两个指标能够从时间和空间上综合全面的评估湖泊换水能力，研究成果可为滆湖水质改善工程提供重要科学参考。面的评估湖泊换水能力，研究成果可为滆湖水质改善工程提供重要科学参考。关键词关键词 滆湖，水体更新时间，水体交换率，换水能力，滆湖，水体更新时间，水体交换率，换水能力，MIKE21 Evaluation and Analysis of Water Exchange Capacity of Gehu Lake Based on MIKE21 Numerical Simulation Chenglong Li1,Shenglian Guo1,Xun Gao2,Zhimin Feng3,Huijie Zhan4 1State Key Laboratory of Water Resources and Hydropower Engineering Science,Wuhan University,Wuhan Hubei 2Zhejiang Design Institute of Water Conservancy and Hydro-Electric Power Co.,Ltd.,Hangzhou Zhejiang 3Ningmeng Hydrology and Water Resources Bureau of Yellow River Conservancy Commission,Baotou Inner Mongolia 4College of Hydrology and Water Resources,Hohai University,Nanjing Jiangsu Received:Jul.17th,2023;accepted:Aug.15th,2023;published:Aug.29th,2023 作者简介：李承龙，男，研究生，研究方向：水资源与水环境，Email: 基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 335 水资源研究 Abstract The water exchange capacity of lakes directly affects the migration of pollutants and the self-purification capacity of water bodies.To explore the influence of different water diversion flows on the water ex-change capacity of the lake,MIKE21 was used to establish a hydrodynamic water quality model of the lake,and the influence of various combinations of working conditions on the overall and local water ex-change capacity of the lake area was simulated and compared.The quantitative analysis of water renew-al time and water exchange rate was proposed to cross-verify the water exchange capacity.Considering the investment amount and effect,the results show that when the water diversion flow is 100,150,200 and 250 m3/s,the water exchange rate reaches over 80%of the whole lake needs 60 d,45 d,32 d and 30 d,respectively.The meshing scheme provides a new idea for the modeling of lake hydrodynamic water quality model,and the dual indicators of water exchange capacity can comprehensively evaluate the lake water change capacity both from time and space,which can provide an important scientific reference for the water quality improvement project of Gehu Lake.Keywords Gehu Lake,Water Renewal Time,Water Exchange Rate,Water Exchange Capacity,MIKE21 Copyright 2023 by author(s)and Wuhan University.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 湖泊在水文循环过程中扮演着重要的角色，是地球上大气圈、岩石圈、生物圈、水圈等各个圈层相互作用的重要节点。近几十年来城市建设进程迅猛，很多城市湖泊受到建设活动和周边用地的影响，出现了面积萎缩、水动力条件不足、水质恶化等一系列问题1。现有研究通常使用引水入湖方式来提升水体水动力条件，增强湖泊自我修复能力，改善湖泊水质水环境状况2。换水能力并没有一种统一的指标来定义，一般的评价指标有水体交换率、水体更新时间、水龄3、滞留时间、半交换时间4、净化时间、水体存留时间5等。黄春琳等6使用水龄分布特征来表征太湖水体交换能力，黄爱平等7建立鄱阳湖水龄模型描述其换水能力，唐继张8等以换水周期为指标观察昆明湖换水能力。当前关于湖泊换水能力的研究通常都是单一指标分析，或从空间上分析湖泊各局部区域的换水特性，或从时间上分析湖泊整体浓度变化，难以兼顾湖泊换水能力在时间和空间、局部和整体的变化特点。为探究引水入湖对滆湖水动力条件的影响，定量分析引水流量变化对滆湖换水能力的改善程度，本文提出水体更新时间和水体交换率两个指标来评估湖泊换水能力，既能兼顾其时空变化和局部与整体之间的关系，又能使两个指标的结果形成交叉验证提高结果的可信度。水体更新时间是指水体内物质量减少到原有物质总量的37%时所需的时间，常用于评估近岸海域水交换的研究，近几年被广泛应用到湖泊、河流等水体的运动特性和物质输移的研究，能够表现出水体各个区域水体更新的空间分布特征9。水体交换率可表征整个水体中某种物质含量的变化，体现整个水体中物质平均浓度随时间的变化情况。使用水体更新时间和水体交换率两个指标来描述滆湖这样较为封闭湖泊的换水能力，同时研究其局部水体更新特征和整体水体交换能力是更加全面的。目前湖泊水动力模型包括 MIKE21、EFDC、Delft3D、太湖模型10 11 12等国内外知名的水质水动力模拟软件。本文采用 MIKE21 软件，建立滆湖水动力水质模型来研究其换水能力。采用对流扩散方法，以水体更Open AccessOpen Access基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 336 水资源研究 新时间和水体交换率两个指标来定量评价换水能力，探讨引水入湖情况下滆湖的水动力条件和换水能力的变化。研究结果对水质改善工程规划具有的实际参考价值，为湖泊水体更新研究提供理论依据。2.研究区域概况研究区域概况 滆湖位于江苏省常州市西南方向，横跨常州市武进区和宜兴市，属于浅水湖泊，平均水深 1.3 m，最大水深4 m，湖泊长度 22 km，最大宽度 9 km，现状水域总面积约为 146 km2，总库容 2.1 亿 m3 13 14。如图 1 湖底地形 3D 图所示，湖泊形状西岸岸线较为平滑，东岸岸线较为曲折，相比湖泊西岸，其死水区域较多。湖泊上设有 5 个观测点(观测点位置见图 1)观测 1#断面处引水流量为 100 m3/s，2#、3#断面处水位为 0.2 m 时，5 个观测点坐标及各点处稳定流速度和水深见表 1 所示。此次模拟设置 1 个规划进水口(见图 1，断面 1#)和 2 个规划出水口(见图 1，断面 2#、3#)。Figure 1.3D map of bottom terrain in the Gehu Lake 图图 1.滆湖湖底地形 3D 图 Table 1.Coordinates,stable flow velocity and water depth of five observation points in the Gehu Lake 表表 1.滆湖 5 个观测点坐标和稳定流速及水深 观测点序号 经纬度 流速/cms1 水深/m 1 E 119.8484752,N 31.6820494 0.103 1.12 2 E 119.8804448,N 31.6396883 0.039 1.107 3 E 119.8467778,N 31.5834785 0.112 1.364 4 E 119.8239562,N 31.5389336 0.036 0.976 5 E 119.7661855,N 31.5051893 0.692 2.659 基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 337 水资源研究 3.研究方法研究方法 3.1.水水动力模型构建动力模型构建 1)数据预处理与模型建立 滆湖水动力学模型采用丹麦水利研究所(DHI)开发的 MIKE21 数学模型软件进行建模，MIKE21 模型是二维平面水流数学模型的一种。该模型在国内已经得到了广泛的运用，常用于河流、水库、湖泊、海湾等水体区域的流场、水质和泥沙运输的模拟。其详细的原理与数学求解方法在本领域很多研究中都有极为详尽的介绍15 16，本文就不再过多介绍。数据预处理过程中，利用 CAD 和 ArcGIS 软件将测绘得到的湖泊边界及高程散点数据的 dxf 文件转换成 xyz文件并导入 MIKE21 中生成模型。将湖泊 5 个观测点的经纬度坐标转换到相应的投影坐标下，利用 MIKE21 自带的数据转换工具进行批量转换并生成观测点坐标文件导入模型。2)网格划分及网格划分方案对比 为适应湖泊复杂的湖底地形和曲折的岸线情况，本文在建模过程中绘制了 4 种不同特点的网格划分方案(如图 2、图 3 所示)，来对比不同网格划分方法对模型模拟精度和计算效率的影响，并比选出最合适的网格划分方案。考虑到湖泊引水出入口处流态更加复杂，所以在 B、C、D 方案中在出入口处不同程度的增加了计算网格的 Figure 2.Mesh diagram of A and B mesh division schemes 图图 2.A、B 网格划分方案网格图 基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 338 水资源研究 Figure 3.Mesh diagram of C,D meshing schemes 图图 3.C、D 网格划分方案网格图 Table 2.Attribute tables for four meshing schemes 表表 2.四种网格划分方案的属性表 网格划 分方案 网格形态及大小 三角网 最小角度 网格 数目 节点数 入口(1#断面)左出口(2#断面)右出口(3#断面)中心湖泊 A 方案 三角网格，单位网格最大面积为 120,000 m2 三角网格，单位网格最大面积为 120,000 m2 三角网格，单位网格最大面积为 120,000 m2 三角网格，单位网格最大面积 为 120,000 m2 33 3379 1862 B 方案 三角网格，单位网格最大面积为 10,000 m2 三角网格，单位网格最大面积为 10,000 m2 三角网格，单位网格最大面积为 10,000 m2 三角网格，单位网格最大面积 为 100,000 m2 33 5695 3011 C 方案 三角网格，单位网格最大面积为 90,000 m2 边长为 200 m 的 矩形网格 三角网格，单位网格最大面积为 90,000 m2 三角网格，单位网格最大面积 为 200,000 m2 33 1461 832 D 方案 边长为 200 m 的 矩形网格 边长为 200 m 的 矩形网格 边长为 200 m 的 矩形网格 三角网格，单位网格最大面积 为 100,000 m2 34 2997 1713 基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 339 水资源研究 密度；同时考虑到引水出入口处流向更加确定，所以 C 方案中在左出口(2#断面)处设置了矩形网格，D 方案中出入口(1#、2#、3#断面)处均设置了矩形网格。4 个方案网格划分主要特点如下：A、B 方案均采用三角网格划分，A 方案网格划分在全湖区域内较为均匀，而 B 方案在入水口和出水口处网格划分进行了有明显的细化处理；C、D 方案采用湖泊中心使用三角网格、出湖口的河道区域使用矩形网格，D 方案相比 C 方案的网格更加密集。各方案的网格形态、网格参数设置见表 2。3)模型主要参数设置 模型重要参数及模型边界设置信息见表 3。Table 3.Construction and parameter setting of MIKE21 hydrodynamic model 表表 3.MIKE21 水动力模型的构建和参数设定 模型构建 参数设置 模型适用性 正压模式，深度平均(2D)1#断面流量边界 100 m3/s 2#、3#断面水位边界 0.2 m 空间计算域 146.18 km2 离散分辨率 网格尺寸变化范围 时间步长 3600 s 计算时长 90 d CFL 条件数 0.8 曼宁糙率系数 空间分布，变化范围 0.0200.091 3.2.模型率定与验证模型率定与验证 1)最优网格方案选取 分别对 4 种网格方案进行多次率定，得到每个网格划分方案下的最优糙率空间分布图，将相应的最优糙率分布图导入对应网格方案模拟计算，得到各方案观测点流速和水深的模拟值，进行误差分析和比较。由分析表4 和表 5 可知，A 方案 5 个观测点水深平均相对误差为 3.4%，相对误差最大值为 7.00%；B 方案水深平均相对误差为 0.44%，相对误差最大值为 0.75%；C 方案水深平均相对误差为 3.70%，最大相对误差为 9.4%；D 方案水深平均相对误差为 1.75%，最大相对误差为 3.90%。4 个方案水深模拟结果最大相对误差均在 10%以内，5 个观测点平均误差均在 4%以内，模拟结果均在合理范围之内。Table 4.Relative errors of water depth simulation values for each scheme 表表 4.各方案水深模拟值相对误差 观测点 实测值/m A 方案 B 方案 C 方案 D 方案 模拟值 相对误差 模拟值 相对误差 模拟值 相对误差 模拟值 相对误差 1 1.12 1.1284 0.7%1.1116 0.75%1.115 0.4%1.1345 1.30%2 1.107 1.1281 1.9%1.1091 0.19%1.089 1.6%1.1082 0.11%3 1.364 1.3367 2.0%1.364 0.00%1.326 2.8%1.3331 2.26%4 0.976 1.044 7.0%0.969 0.72%0.885 9.4%0.9647 1.16%5 2.659 2.5163 5.4%2.6448 0.53%2.546 4.3%2.5554 3.90%平均相对误差 3.40%0.44%3.70%1.75%基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 340 水资源研究 而流速的相对误差各方案之间差距较大，A、B 方案平均误差分别为 25.66%、90.76%，误差过大，与实测值差距过大；D 方案平均误差 16.75%，但最大误差为 47.19%，与实测值仍有较大差距；C 方案平均误差 12.26%，最大误差不超过 25.0%，流速模拟精度与实测值相差不大，可认为 C 方案模拟得到湖泊流场与实际情况基本相符。Table 5.Relative error of flow velocity simulation value of each scheme 表表 5.各方案流速模拟值相对误差 观测点 实测值 A 方案 B 方案 C 方案 D 方案 模拟值 相对误差 模拟值 相对误差 模拟值 相对误差 模拟值 相对误差 1 0.103 0.0715 30.6%0.066 35.9%0.096 6.5%0.1339 29.96%2 0.039 0.0515 32.1%0.104 166.7%0.042 6.4%0.0396 1.42%3 0.112 0.0941 16.0%0.288 157.1%0.096 14.0%0.1125 0.45%4 0.036 0.0412 14.4%0.066 83.3%0.045 25.0%0.0253 47.19%5 0.692 0.4482 35.2%0.767 10.8%0.757 9.4%0.7248 4.73%平均相对误差 25.66%90.76%12.26%16.75%在相同计算机配置、相同计算时长和步长情况下，各方案的运行时长见表 6 所示。从表 6 可见，在相同计算机配置和相同计算时长、计算步长情况下 C 方案的运行时间最短，说明 C 方案的计算效率最高，同时也说明计算网格数量越少计算效率越高。综合各网格划分方案的水深、流速模拟精度、计算效率，可认为四个方案中 C 方案为最优的模型网格划分方案，并使用 C 方案作为最终模型方案。Table 6.The running time of each program 表表 6.各方案的运行时长 网格划分方案 计算时长/d 计算步长/s 运行时间/min A 方案 90 3600 60 B 方案 90 3600 155 C 方案 90 3600 12 D 方案 90 3600 26 2)最优方案的结果 通过不断地对模拟期水动力参数进行率定和调整，得到 C 方案的最优糙率分布图如图 4 所示。湖泊糙率分布设置靠近湖中心位置糙率较小，岸边区域靠近死水区、淤塞区域糙率较大。进行水动力模拟得到滆湖流场矢量图如图 5 所示，模拟结果与实际流场基本一致。从图 6、图 7 可见，模型验证结果表明，经过一定时间的模拟模型的水深、流速均收敛到稳定值，并在小幅度范围内波动与实际情况相符。在模型运行过程中，水深在模型模拟 9 天便达到稳定，各观测点水深分别收敛至 1.115 m、1.089 m、1.326 m、0.885 m、2.546 m，各观测点水深与对应实测稳定水深平均相对误差在 3.70%以内。流速在运行过程中波动较大，在图 6 中已剔除前 25 天流速波动过大的数据。流速在 45 天左右达到稳定，各观测点流速分别收敛到 0.096 cm/s、0.042 cm/s、0.096 cm/s、0.045 cm/s、0.757 cm/s，各观测点流速与对应实测稳定流速平均相对误差在 12.26%以内。模型验证结果表明，该模型能够较为准确的模拟滆湖水动力情况，可以用于水体更新时间和水体交换率的研究。基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 341 水资源研究 Figure 4.Optimal roughness distribution map of scheme C Figure 5.Illustration of flow field vector in the Gehu Lake 图图 4.方案 C 最优糙率分布图 图图 5.滆湖流场矢量图 Figure 6.Variation curve of water depth at each observation points Figure 7.Variation curve of flow velocity at each observation points 图图 6.各观测点水深变化曲线 图图 7.各观测点流速变化曲线 3.3.换水能力两个指标换水能力两个指标定义定义 1)水体更新时间 水体更新时间是基于水箱模型的一种水交换定义指标5，是用于评估湖泊水体内纯粹由于对流扩散作用导致湖泊物质交换速率的一个常用指标。可以理解为湖泊内均匀分布某种稳定不可降解的溶解态示踪剂。本文设置这种示踪剂浓度为 1(无量纲)，从湖泊入口处引入清水且示踪剂浓度为 0。湖泊水体引水过程中，受对流扩散作用影响水体中各网格内示踪剂浓度降低到初始浓度的 37%时(即示踪剂浓度降低至 0.37 时)，所消耗的时间为水体更新时间。水体更新时间表达式如下：()()037%renewalTT C tC=(1)式中：renewalT表示水体更新时间，()C t表示 t 时刻示踪剂浓度，0C表示示踪剂初始浓度，()()037%T C tC表示示踪剂浓度降低到初始浓度 37%以下时的时间。基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 342 水资源研究 2)水体交换率 水体交换率是指湖泊水体内部浓度为 1 的与上述相同的稳定示踪剂，在经过湖泊对流扩散作用一段时间后，输移至水体外部的物质总量与初始物质总量的比值17。水体交换率反映湖泊内水体整体的平均换水程度。水体交换率的表达式如下：()()0000NnnCShCtS nhECSh=(2)式中：E表示水体交换率，0C、0S、h分别表示初始浓度、湖泊总面积、湖泊平均水深，N表示湖泊网格总数为 1461 个，n表示网格序号，()nCt第n个网格时刻的示踪剂浓度，()S n第n个网格的面积。3.4.模拟方案与模型工况设置模拟方案与模型工况设置 1)模拟方案 本文以水体更新时间和水体交换率作为指标来评估研究滆湖水体交换能力。模拟方案为：如图 1 所示的 1#断面作为湖泊引水入口，2#、3#断面为出水口，在引水入流的水力驱动下改善湖泊水动力条件，增加湖泊水体对流扩散能力。为方便进行水体更新时间和水体交换率的计算，模型假设湖泊水体充分混合，并存在某种不可降解的稳定溶解态示踪剂浓度为 1(无量纲)，引水清水入流示踪剂浓度为 0。2)工况设置 湖泊当前原工况引水流量为100 m3/s，相较于引水前湖泊的水动力条件和水质已经得到了较为明显的改善。但由于湖泊死水区域较多，特别是由于右岸区域岸线曲折导致右岸区域水体更新依然较慢。因此，为进一步提高湖泊水动力条件、对比和量化不同引水流量对湖泊水体交换能力的影响，使模拟结果能更好的服务于实际引水工程改造，本文设置了 3 种加大引水流量的工况。工况 1：加大引水流量至 150 m3/s；工况 2：加大引水流量至 200 m3/s；工况 3：加大引水流量至 250 m3/s。4.结果与分析结果与分析 4.1.水动力特征水动力特征 引水流量为 100 m3/s 时，根据图 5 流场矢量图显示，滆湖河口处流速较大且流速方向较为一致，在右侧出口处岸线凹槽处流向较为紊乱。在湖体中央流速相较湖边流速更大且方向更加有序，说明引水入湖能够有效的提高湖泊的流动性，增强水动力条件。两岸流速相较于湖中心流速较慢，且流向相对更加无序，右岸部分区域形成涡流，说明虽然湖泊的水动力条件增强但对湖泊的死水区域影响较小，仍需进一步加大流量或采取其他措施来增强右岸区域的水体流动性。流场整体状况与实际情况相近，说明模拟结果具有可靠性。4.2.水体更新时间水体更新时间 根据滆湖水动力模型模拟结果绘制的各工况下水体更新时间图 9 可见，湖泊的水体更新时间大致存在以下规律：越靠近出入水口，水体更新时间越短，反之则越长，且湖泊东岸存在较多死水区；湖泊中央区域较左右岸边区域的水体更新时间更短。根据图 8 可见，引水流量加大后湖泊中心区域、左岸的水体更新时间均缩短，右岸水体更新时间过长的区域面积有明显变小，水体更新时间小于 10 天的区域明显变大。根据 UNEP 流域综合管理指南建议湖泊水体更新时间一般不超过 30 天。当入水口流量达到 200 m3/s 时，湖泊由中心向两岸扩散的区域整体水体更新时间基本下降至 30 天以内。而引水流量由 200 m3/s 增加至 250 m3/s 时，工况 3 相比于工况 2 水体更新时间在 30 天以内的区域面积并没有很明显的增加，在右岸区域水体更新时间大于 60 天的区域虽然稍有减少，但依然存在死水区域 基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 343 水资源研究 Figure 8.Water renewal time chart of water body under various working conditions 图图 8.各工况水体更新时间图 水体更新时间过长的问题。加大入流流量能够缩短该湖泊的水体更新时间，但该湖泊右岸部分死水区仍然存在基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 344 水资源研究 换水周期较长，污染物滞留的现象。4.3.水体交换率水体交换率 根据滆湖水动力水质模型各工况下的示踪剂变化模拟结果，计算得到滆湖 15 d、30 d、45 d、60 d、75 d、90 d 的水体交换率。滆湖水体交换率计算结果表见表 7，水体交换率变化曲线见图 9。Table 7.Water exchange rate under various working conditions 表表 7.各工况下水体交换率 工况 引水流量/(m3/s)水体交换率/%15 d 30 d 45 d 60 d 75 d 90 d 原工况 100 24.67 51.25 71.07 81.92 88.43 92.31 工况 1 150 36.86 69.00 82.42 89.30 93.13 95.41 工况 2 200 49.52 78.10 88.31 93.26 95.90 97.39 工况 3 250 60.01 82.47 90.61 94.58 96.71 97.93 Figure 9.Variation diagram of water exchange rate under various working conditions 图图 9.各工况水体交换率变化图 由图 9 可见，各工况的水体交换率随引水时间的延长而增大，且均呈现前期水体交换率的增长较快，后期水体交换率的增长较慢的趋势。说明引水时间越长湖泊内原有的污染物因对流扩散作用运移到湖泊之外的污染物越多，且越到后期因湖泊内污染物浓度的降低水体交换率也随之变化得越慢。将四个工况作比较发现，引水流量越大水体交换率的增长越快，相同时间内水体交换率变化也越大，说明加大引水流量是加快降低湖泊整体污染物浓度的一个有效办法。从表 7 中展示的数据可以看到，经过 15 天的引水入湖，原工况和工况 1、2、3 的水体交换率分别达到了24.67%、36.86%、49.52%和 60.01%。说明在短时间内，加大引水流量能够快速明显提高水体交换率，当引水流量达到 200 m3/s 以上时，滆湖水体能够在 15 天就能将湖泊污染物整体浓度降低到原来一半的效果。进一步分析基于 MIKE21 数值模拟的滆湖换水能力评估分析研究 DOI:10.12677/jwrr.2023.124038 345 水资源研究 可以看到，水体交换率达到 80%以上，原工况需要 60 天、工况 1 需要 45 天、工况 2 需要 32 天、工况 3 需要30 天，说明引水流量越大整体浓度降低 80%以上所需的时间越短，与水体更新时间计算结果形成了相互验证。从湖泊水动力条件来看，湖泊流场与实际情况相符合，湖中心流速较大两岸流速较小，水动力特征明显。从水体更新时间来看，工况 2 引水流量为 200 m3/s 时湖泊中心区域大部分水体更新时间均已经达到 30 天以内，湖泊右岸局部区域水体更新时间过长；且工况 3 继续增大引水流量未见特别明显改善。从水体交换率来看，湖泊整体的水体交换率达到 80%以上工况 2 需要 30 天左右，工况 3 需要 30 天。5.结论结论 1)在建模划分网格时，加密出入水口处网格的密度能够降低水深模拟的相对误差，但对流速模拟的相对误差没有明显提升作用，甚至可能会增大误差。在二维湖泊水动力学模型当中，河道区域使用矩形网格或者三角网格对模型模拟精度没有明显影响。模型网格的数量能够明显影响到模型运行速度，所以在保证精度的情况下可以尽量减少网格数量以提高模型运行效率。2)模型模拟结果水深和流速均能够收敛到稳定值，且模拟值与实测值相对误差较小、结果相吻合，说明所构建的模型具有合理性。3)增大湖泊引水流量能够明显提高湖泊的水动力条件。当工况 2 流量增加到 200 m3/s 时湖泊大部分区域水体更新时间均小于 30 天，湖泊水体 15 天交换率可达到近 50%、30 天湖泊水体交换率达到近 80%。因考虑到工况 3 继续加大补水流量带来的边际效益递减，所以本文建议在实际工程中，引水流量加大至 200 m3/s 即可。既可以缩短湖泊各区域水体更新时间，又能提高湖泊整体的水体交换率，保持湖泊水质指标的稳定。4)同时使用水体更新时间和水体交换率两个指标来评估湖泊的换水能力是更加系统和全面的。水体更新时间可以更加清晰直观的分析湖泊换水能力空间上的分布特征，能看到湖泊内局部区域乃至每个网格的水体更新时间。水体交换率可帮助分析湖泊换水能力的时间变化特点，能够分析每天乃至每个小时湖泊整体的污染物减少量。两个指标结合，既能看到湖泊换水能力的空间特征，又能掌握湖泊换水能力的时间变化特点，同时观测局部变化又能把握整体情况。基金项目基金项目 国家自然科学基金区域创新发展联合基金(U20A20317)。参考文献参考文献 1 杨桂山,马荣华,张路,等.中国湖泊现状及面临的重大问题与保护策略J.湖泊科学,2010,22(6):799-810.YANG Guishan,MA Ronghua,ZHANG Lu,et al.Current status of lakes in China,major problems and protection strategies.Lake Science,2010,22(6):799-810.(in Chinese)2 王健,徐望朋,马方凯,等.城市湖泊治理思考与建议J.人民长江,2022,53(2):41-47.WANG Jian,XU Wangpeng,MA Fangkai,et al.Thoughts and suggestions on urban lake governance.Yangtze River,2022,53(2):41-47.(in Chinese)3 ZHANG,P.,LI,L.,WANG,Y.,et al.Influence of riverbed incision and hydrological evolution on water quality and water age based on numerical simulation:A case study of the Minjiang Estuary.International Journal of Environmental Research and Public Health,2021,18(11):6138-6147.https:/doi.org/10.3390/ijerph18116138 4 王平,孙家文,董祥科,等.生态修复后普兰店湾水体交换数值研究J.海洋环境科学,2021,40(6):937-946.WANG Ping,SUN Jiawen,DONG Xiangke,et al.Numerical study on water exchange in Pulandian Bay after ecological resto-ration.Marine Environmental Science,2021,40(6):937-946.(in Chinese)5 魏皓,田恬,周锋,等.渤海水交换的数值研究水质模型对半交换时间的模拟J.青岛海洋大学学报(自然科学版),2002(4):519-525.WEI Hao,TIAN Tian,ZHOU Feng,et al.Numerical study on the exchange of sea water in Bohai SeaSimulation of half ex-change time by water quality mode