基于MIKE 21的圩区河网水动力调控方法研究.pdf

水利水电技术(中英文)第 54 卷 2023 年第 7 期Water Resources and Hydropower Engineering Vol.54 No.7徐存东,訾亚辉,黄嵩,等.基于 Mike21 的圩区河网水动力调控方法研究J.水利水电技术(中英文),2023,54(7):161-170.XU Cundong,ZI Yahui,HUANG Song,et al.Research on the hydrodynamic regulation method of river network in polder area based on MIKE 21J.Water Resources and Hydropower Engineering,2023,54(7):161-170.基于 MIKE 21 的圩区河网水动力调控方法研究徐存东1,2,3,訾亚辉1,2,黄 嵩1,2,赵志宏1,任子豪1,2,胡小萌1,王海若1(1.华北水利水电大学 水利学院,河南 郑州 450046;2.浙江省农村水利水电资源配置与调控关键技术重点实验室,浙江 杭州 310018;3.河南省水工结构安全工程技术研究中心,河南 郑州 450046)收稿日期:2022-11-26;修回日期:2023-03-04;录用日期:2023-03-07;网络出版日期:2023-05-10基金项目:国家自然科学基金项目(51579102);浙江省基础公益研究计划项目(LZJWD22E090001);浙 江省重 大科技 计划项 目(2021C03019);中原科技创新领军人才支持计划(204200510048);河南省科技攻关项目(212102310273);河南省高等学校重点科研项目计划(20A570006)作者简介:徐存东(1972),男,教授,博士研究生导师,博士,主要从事水利工程优化设计理论方法及平原河网调控研究。E-mail:xcun-dong 通信作者:訾亚辉(1998),男,硕士研究生,主要从事河网水动力方面研究。E-mail:1971685442 Editorial Department of Water Resources and Hydropower Engineering.This is an open access article under the CC BY-NC-ND license.摘 要:【目的】我国南方平原地区地势低缓,河网结构复杂,为快速解决该地区因河网水动力不足而引起的水环境问题,急需研究河网水体运行状况调控方案。【方法】以太湖流域南浔区沈庄漾“溇港圩田”河网作为研究对象,构建 Mike 21 FM 二维河网水动力模型,结合研究区河网改造后的水利工程布置和调度规则,提出定位增压的 7 种水动力调控方案,模拟分析河网在不同方案下流场分布规律,探究圩区河网水动力提升和流场改善的优化调度方法。【结果】模拟结果表明:活水和引排水方案虽能有效改善河网水动力,但对于中心水塘的水动力驱动效果不佳。在开启西南侧河道水泵和西部泵闸的情况下,河网内大部分区域的流速基本能维持在 0.05 m/s 左右,部分河道的流速能达到 0.3 m/s 以上,并且对中心水塘水体驱动效果最好。【结论】综合分析表明:所设计的 7 种调控方案均能有效地提升河网水动力,在开启西南侧河道 1 台水泵和西部泵闸时,流场分布最均匀且水动力整体提升效果最佳。研究结果可为南方平原河网的水利工程水动力调控提供参考。关键词:水动力模型;调控方案;圩区河网;流场分布DOI:10.13928/ki.wrahe.2023.07.014开放科学(资源服务)标志码(OSID):中图分类号:TV131.2文献标志码:A文章编号:1000-0860(2023)07-0161-10Research on the hydrodynamic regulation method of river network in polder area based on MIKE 21XU Cundong1,2,3,ZI Yahui1,2,HUANG Song1,2,ZHAO Zhihong1,REN Zihao1,2,HU Xiaomeng1,WANG Hairuo1(1.School of Water Conservancy,North China University of Water Resources and Electric Power,Zhengzhou 450046,Henan,China;2.Key Laboratory for Technology in Rural Water Management of Zhejiang Province,Hangzhou 310018,Zhejiang,China;3.Henan Provincial Hydraulic Structure Safety Engineering Research Center,Zhengzhou 450046,Henan,China)Abstract:ObjectiveThe southern plain area of China is low and gentle,and the river network structure is complex.In order to quickly solve the water environment problems caused by the lack of water power in the river network in this area,it is urgent to study the regulation scheme of the river network water operation.Methods This study took the“Lougang Polder”river 161徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期network in Shenzhuangyang in Nanxun District of the Taihu Lake Basin as the research object,constructed the hydrodynamic model of Mike 21FM two-dimensional river network,combined with the layout and scheduling rules of water conservancy projects after the transformation of the river network in the study area,proposed seven hydrodynamic regulation schemes for positioning and pressurization,simulated and analyzed the distribution of flow fields under different schemes of river net-work,and explored the optimal scheduling method of hydrodynamic improvement and flow field improvement of river network in the polder area.Results The simulation result show that although the live water and water diversion and drainage scheme can effectively improve the hydrodynamic force of the river network,the hydrodynamic driving effect of the central pond is not good.With the southwest river water pump and the west pump gate open,the flow velocity in most areas of the river network can basically maintain around 0.05 m/s,and the flow velocity in some rivers can reach more than 0.3 m/s,and the driving effect on the central pond water body is the best.ConclusionThe comprehensive analysis shows that the seven control schemes designed can effectively improve the hydrodynamic force of the river network.When one water pump and the west pump gate are opened in the southwest river channel,the flow field distribution is the most uniform and the overall hydrodynamic force improvement effect is the best.The research result can provide a reference for the hydrodynamic regulation of water conservancy projects in the southern plain river network.Keywords:hydrodynamic model;regulation scheme;polder river network;flow field distribution0 0 引引 言言 近年来,随着城市的快速发展,人类活动与平原河网水系间的相互作用愈发显著。流域内河网的水动力不足以及水质恶化等问题已经成为制约流域经济社会可持续发展的主要瓶颈1-4。天然河网水系的水体动力主要受限于水文、河床地形、汇流条件等结构性因素,在一定的尺度内,其水动力变化与河网水环境存在着一定的相关性5-6。我国南方平原地区作为我国社会经济最发达的区域,随着人类活动的增加和水资源的过度开发,导致该区域水资源承载力已趋近于极限,平 原 河 网 水 生 态、水 环 境 等 问 题 日 益 凸显7-8。由于南方平原河网水体自身流动性较差、空间较封闭的特点难以彻底改变,所以在天然河网上兴修闸、泵等工程设施,是我国南方地区流域水资源有效开发、兴利排涝的重要工程手段。太湖流域的南浔区流域河道密布,纵横交错,湖泊星罗棋布,是典型“江南水网”。南浔区河网水流平均流速普遍不足0.01 m/s,河道中的各种污水长期游荡、停滞在河道中,导致其不能顺利排出。天然的河网被人工改造,不同程度地改变了河网的自然结构,原本完整的河网系统被分割成多个不连续的片段,截断的河道使得河网整体水动力不足,在引起其水动力特性改变的同时,也使得天然河网水系的连通性大幅降低,部分河网甚至出现了死水区9。与此同时,国内外学者开展了水动力数值模拟。STELL-ING10基于浅水河网特性,对河网结构进行了系统的描述,并对平面浅水方程详细地作了研究,成功模拟浅水河网的水动力变化,进一步推动了二维模型的发展;BOUCHARD 等11利用水动 力 分 析 模 拟 程 序WASP 模拟纳米管在地表水运输中的效果,分析其水动力变化情况;AJIWIBOWO 等12以 CBL 河道为研究区,利用 MIKE 11 进行一维数值模拟,成功模拟了丰水期和枯水期疏浚调水对河道水动力的影响,研究结果证明了不同流量的引调与水动力的改善之间存在着一定的关联性,该结果为后续研究提供了理论支持;WALDMAN 等13将 MIKE 3 与 Delft 3D 进行有效的结合,成功开发了 Pentland Firth 和 Orkney 两个区域尺度的水动力模型,通过模型的使用,成功预测了潮 汐 涡 轮 机 的 运 行 对 两 个 水 域 水 动 力 的 影 响;MAWANDHA 等14以本加湾索罗河下游水系为研究对象,并基于 HEC-RAS 模型对该水系进行模拟,成功利用该模型对水利调度措施进行了有效评估,为未来调度方法的研究提供理论参考。对于国内研究,在水系连通方面,潘剑光等15以平原区城市为例,在河网水系连通的条件下,对泵站进行合理的设置,增强水体的流动性,通过利用 MIKE 11 水动力模型,进行模拟不同情景下河网水系的水动力条件,并为水环境的治理提供新思路。在闸泵联合调度方面,何健等16构建水动力-水质耦合模型,设置河网闸泵不同的运行方案,分析不同方案下河道水动力改善效果,表明泵站引水量可在不同程度上改善河网水动力条件;王思如等17建立了基于长江下游典型平原城市的河网水动力学模型,并进行数值模拟,具体分析工程设施的布置和调度方式等对河网水体流动的改善效果,研究结果为平原河网城市地区的水动力调控提供新技术;闫毓等18以上海中心城区蕰南水利控制片261徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期 图 1 研究区河网地理位置示意Fig.1 Geographic location map of river network in research area为研究对象,通过对研究区的闸泵进行联调联控调度,对河网进行水动力时空再造,依据现有的水利工程设施进行多种调控方案设计,并对其进行模拟计算,分析其活水改善效果并对其进行结果优选,研究结果大幅度增强了水体的流动,对水动力运行状态的维持具有重要意义。综上所述,随着数值模拟技术的不断发展,在水系连通和闸泵联合调度方面对水动力改善研究均取得了丰富的成果,但发现大多数平原河网的研究只关注于完整水系之间所存在的水动力问题,而面对河网自然结构的改变而引发的水动力不足现象,系统地开展在特定的封闭活水区域生态补水方案的研究尚存不足,导致无法验证调水、补水方案的普适性与有效性。本研究以太湖流域南浔区沈庄漾“溇港圩田”河网作为研究对象,通过分析研究区域水文特征和河网地形,选取具有代表性点位进行模型率定,建立相应的二维水动力模型,从河网内活水和外河网引排水两个方面,模拟分析不同工况下河网水动力调控效果,水体流动均得到了不同程度的提高,解决了平原河网水动力不足问题。研究成果对于修复南方平原河网,促进南方平原河网区可持续发展具有积极意义。1 1 研研究究区区概概况况1.1 研究范围 南浔区隶属于浙江省湖州市,位于杭州市东北部。其地势较低缓,河道交错,结构复杂,是典型的水网平原。当地地势较低,平均海拔在 5 m 以下,境内河流纵横,湖漾密布。本研究目标河网位于南浔区东南侧的沈庄漾圩区,北邻南长兴港,东倚沈庄漾,西靠阳安塘,南依博成桥港。圩区面积约 0.98 km2,水域面积占比约 15.3%,超过南浔区平均水域面率。研究区河网现状如图 1 所示。1.2 研究区现状 研究区河网共有 3 处与外河网相连,分别位于研究区河网的西、北、东侧,通常情况下,3 处河道进行水闸关闭,仅依靠南部的水渠向河网内供水,以满足日常使用需求。研究区河网分为外圈河道与中心水塘2 个部分,水位保持在2.5 m,并设有2 道拦水堰。为有效提取各处的水体流速情况,本研究在河网内设置了 8 个监测点,用于提取数据与分析,河网地形及各水利工程与监测点位位置如图 2 所示。361徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期图 2 研究区河网地形及水利工程与监测点示意Fig.2 Schematic diagram of river network topography,hydraulic engineering and monitoring points in the research area2 2 模模型型构构建建2.1 水动力模块控制方程 MIKE 21 FM 模拟模型的水动力模块基于三向不可压缩和 Reynolds 值均布的 Navier-Stokes 方程,并服从于 Boussinesq 假定和静水压力的假定19。二维非恒定平面浅水流方程组如下ht+hux+hvy=hS(1)hut+hu2x+huvy=fvh-ghx-h0Pax-gh220 x+sx0-bx0-10sxxx+sxyy()+(hTxx)x+(hTxy)y+husS(2)hvt+huvx+hv2y=-fuh-ghy-h0Pay-gh220y+sy0-by0-10syxx+syyy()+(hTxy)x+(hTyy)y+hvsS(3)式中,t为时间;x、y为笛卡尔坐标系坐标方向;为水位;d 为静止水深;h 为总水深,h=+d;u 为 x方向上的速度分量;v 为 y 方向上的速度分量;f 为哥氏力系数,f=2sin,代表地球自转角速度,代表当地纬度;g 为重力加速度;为水的密度;0为淡水的密度;Pa为当地大气压强;sxx、sxy、syx、syy为辐射应力分量;sx、sy为水面风应力张量;bx、by为河床床面应力张量;S 为源项;us、vs 为源项在x、y 方向的水流流速。2.2 初始条件与边界条件的确定 在水动力模块,设定输出结果为水位与流速。在初始条件设置上,圩区河网的初始水位设置为日常水位 2.5 m。在边界条件设置上,设置西、北、东部泵闸限制水位 2.5 m,流量根据实际工程进行确定。2.3 模型参数的设置 在构建模型时,除了河网地形、初始与边界条件外,还需要设置求解格式、干湿水深、河床糙率等20。求解格式选择低阶,干湿边界使用软件中存在的预设值,即 h干为 0.005 m,h淹没为 0.05 m,h湿为 0.1 m。涡粘系数设定为 0.28。本研究选择曼宁系数反应河床床底阻力,曼宁系数值初定为模型默认值32 m1/3/s。2.4 参数率定 本研究需要对河床糙率进行率定,即模型中的曼宁系数,采用 2021 年 11 月 3 日 8:0016:00 的水位实测数据对模型进行率定,模拟过程设置步长 60 s,选取水位验证点为 No.2 与 No.8。经过率定,确定河床曼宁系数 32 m1/3/s 能满足研究的精度要求。为保证研究的科学性,采用 2021 年 12 月 6 日 8:0016:00 的实测数据进行进一步验证,结果表明,水位点模拟值误差小于 0.001 m 的合格率较高,且高于461徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期90%,模型满足模拟要求。两处监测点位的率定与验证结果如图 3、图 4 所示。3 3 调调控控方方案案设设计计与与模模拟拟3.1 调控方案 为解决研究区河网水动力不足的问题,本研究基于圩区河网改造后的水利工程布置以及调度规则,设计了 7 种圩区河网水动力调控方案。方案编号为,各方案的泵闸调控规则如表 1 所列。3.2 方案模拟与分析3.2.1 活水方案、水动力情况由于外河网水位与内河网水位存在一定水位差,活水方案遵循研究区的河网调度规则,因此,水体由 西向东方向抽出。活水方案、河网内部初始水位均设置为 2.5 m,活水方案仅依靠河道水泵抽水运行调度,因此共同模拟分析,为更直观观察方案与方案工况下的流场情况,故将其做单独分析,河网流速云图与流场分布情况如图 5 所示。在方案和方案工况下,圩区河网的水动力情况均得到了明显的提升,外圈河道的水动力提升效果要强于中心水塘。这表明泵站对中心水塘内水体的驱动能力有限,随着水塘内水体远离与外圈河道相连的进出口河道,流速也随之快速减弱。相比之下,方案对于外圈河道水动力的提升效果更好,流速基本能维持在 0.080.24 m/s,外圈河道的流速几乎提升了1 倍。图 3 模型水位率定结果Fig.3 Model water level calibration results图 4 模型水位验证结果Fig.4 Verification results of model water level表 1 研究区河网水动力调控方案Table 1 Hydrodynamic reconstruction scheme of river network in research area方案编号调度方案类 型河道水泵抽活水(1 台,1 m3/s),西、北、东三处泵闸关闭河道水泵抽活水(2 台,2 m3/s),西、北、东三处泵闸关闭西部泵闸引水(1.24 m3/s),北部泵闸排水(控制闸门处水位 2.5 m)西部泵闸引水(1.24 m3/s),东部泵闸排水(控制闸门处水位 2.5 m)西部泵闸引水(1.24 m3/s),北部+东部泵闸排水(控制闸门处水位 2.5 m)西部泵闸引水(1.24m3/s)+河道水泵抽活水(1 台,1 m3/s)西部泵闸引水(1.24 m3/s)+河道水泵抽活水(2 台,2 m3/s)活水方案外河网引排水活水+引排水组合方案561徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期图 5 方案、工况下河网流速、流场分布Fig.5 Flow velocity and flow field distribution of river network under scheme and working conditions 分析比较图 5(b)与(d)的流场分布可知,方案和方案满足了方案对于外圈河道的水动力驱动目标,且方案和方案所示的外环河段的流场分布基本完全一致,通过西南部水泵将水体排出,顺时针流向河道,在 No.5 附近分支,并在 No.6 处重新汇集。方案与方案的中心水塘的流场分布规律也基本接近,流场由 No.5 所在河道进入水塘,推动中心水塘水体运动,在水塘内逆时针旋转 180后汇入No.6 所在河道,且中心水塘的流场分布均呈现“北密南疏”的特征。整个中心水塘流场流速最大的位置在 No.5 与 No.6 所在的河道内,水塘内部流场分布比较稀疏,说明随着水塘内水体远离进出口河道,其水动力随之快速减弱,这一点也验证了在方案与方案的工况下河网各监测点位的流速模拟结果,位于水塘内的 No.7 与 No.8 的流速明显小于河网内其他各点。虽然 2 个方案的中心水塘的流场分布比较接近,但活水方案工况下的水塘流场分布相比于方案更密集,且流场密集分布的范围更大,这也证实了方案条件下 No.7 的流速比方案条件下的流速大206%的原因。3.2.2 引排水方案、水动力情况对于引排水方案,设置引水流量为 1.24 m3/s,排水的边界条件遵循 MIKE 软件中的边界条件设置,包括陆地边界和进出水口边界,泵闸控制水位均设置为 2.5 m,排水流量为在给定闸门水位的基础上下泄河网的整体流量,同时,保持河网初始水位 2.5 m 不变。同样为更直观观察方案工况下的流场情况,将其做单独分析,方案工况下的河网流速云图与流场分布情况如图 6 所示。在引排水方案的作用下,外圈河道的流速得到661徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期 图 6 方案、工况下河网流速、流场分布Fig.6 The distribution of flow velocity and flow field of river network under scheme III,IV and V working conditions了显著提升,效果与活水方案接近。引排水方案工况下,外圈河道的流速相比方案更加均匀,整体流速处于0.030.12 m/s 范围内,存在部分点位流速较快,可达 0.120.20 m/s。但是,与方案不同的是,这些速度快的位置分布在外环河网的不同区域,而不是在某一个特定的河段上。方案作为方案与方案的结合,具有方案与方案的部分流速分布特征,外圈河道的整体流速基本在 0.030.12 m/s。但是,由于河道水体在流经北部与东部水闸时排出河网,导致北部水闸-No.2-东部水闸一线的河道流速761徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期 图 7 方案、工况下河网流速、流场分布Fig.7 The distribution of flow velocity and flow field of river network under scheme、working conditions相比方案与方案明显减小。通过对图 6(b)、(d)、(f)的流场分布分析,可以看出,在引水方案条件下,河网水体的运动规律与活水方案有很大的区别。水体自西闸门流入河网后,分成两股,一股沿顺时针向北方向流动,在No.5 与 No.6 所在河道与中央水塘水体交换,则另一条往南逆流,最后在北边的水闸汇合。方案与方案的水体运动趋势基本相似,不同的是,在方案条件下,东部闸门两侧的两个流场汇合点并不在闸门的中心,而是在闸门的右侧,这是由于闸门北侧的河道水体流速比南侧河道水体流速要快,再加上东侧水闸附近的河道形态共同造成的。在方案工况下,东部水闸处的流场,依旧分布在水闸右侧,东部水闸北段河道的水流流量和流速均都比方案小,但其流场分布与方案基本相似。结果表明,河网东侧闸门处的流场分布受河道地形的影响较大,而闸门两侧流速的大小对其影响甚微。通过对引排水方案中心水塘流场分布分析,方案与方案中心水塘的流场分布基本一致,但其流场分布的密度和范围要比方案小。通过对 No.7 与No.8 的流速分析,结果表明,方案对水塘内水体的驱动作用比活水方案要小,但对于水塘进出口河道附近的水体,其水动力驱动效果要优于活水方案。河网中心水塘在方案的工况下,流场分布要比方案工况下的水塘流场分布稀疏,与方案工况下的水塘流场分布相似。在靠近水塘进出口河道的位置,方案工况下的流场分布范围介于方案和方案范围之间,并且一旦远离进出水口河道,流场分布将快速衰减。3.2.3 活水+引排水组合方案、水动力情况方案开启了位于圩区河网西南侧(1 台)和西部水闸处的水泵,而方案增添了西南侧的 2 台水泵,该两种调控方案排水边界条件与外河网引排水方案边界条件设置一致,并且泵闸控制水位同样设置为2.5 m,方案、工况下的河网流速云图与流场分布情况如图 7 所示。861徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期由图 7(a)、(c)分析可知,方案、的河网流速情况基本接近,河网西部水闸处和外圈河道的西北部分流速相对较快,在0.100.50 m/s 之间,虽然方案比方案多开启 1 台水泵,在泵闸附近河道流速相比方案有所提升,但是位于两处水泵间的河道流速却并未加快,与河网大部分区域相似,流速在0.10 m/s 以下。由于西南河道开设水泵数量的增加,在水泵东侧河段,存在部分河道流速加快,但范围较小。通过对图 7(b)、(d)的流场分布分析可以发现,在该两种组合方案作用下的河网流场分布非常相似,相对于方案来说,并未因为水泵的增加流场分布而发生明显的变化,与之前的方案不同,外圈河网存在两个流场分布比较稀薄的区域,即北部水闸处以及两处水泵之间的河道,这与相对应的流速分布情况基本吻合。在设计的两种组合方案下,水塘流场分布相比于之前所有方案更密集,且流场分布的范围更广,表明该工况下单位时间内中心水塘水体流动范围更大,强度更高,对水动力提升效果更佳。3.3 调控效果综合分析 通过对 7 种方案模拟研究和流场分析表明,7 个调控方案均能有效提升河网的水动力状况,河道流速的提升在外圈河道表现的犹为明显,中心水塘内流速较快的位置主要集中在 No.5 与 No.6 附近。活水和引排水方案下的河网水体流场分布与流速规律相似,但是不同的方案使得外圈河道的流场分布略有不同,所有设计方案中,外圈河道的流场分布均要比中心水塘的密集。相比较而言,活水方案的水塘流场分布比引排水方案影响范围更大,而引排水工况下的流场分布比活水方案更为集中,两种方案对中心水塘的水动力驱动方面各有优劣。而在活水和引排水组合方案下,方案、对外圈河道流速的提升相对较快,两种组合方案对中心水塘的的流场分布表现的更密集,流场分布的范围更广,对中心水塘水体流动的驱动效果更好。由于方案比方案多增添了 1 台水泵,相对而言,考虑到经济性与有效性的协调发展,方案工况下的消耗更大,所以从整体改善效果与经济性方面定量分析得,方案流场分布最均匀,且范围分布更广并对河网整体水动力提升效果最优。4 4 结结 论论 (1)靠近闸泵的进出口处流速最大,并且呈现随着水塘内水体远离与外圈河道相连的进出口河道流速也随之快速减弱的负相关趋势,所设计调控方案中,流速值最大点位和最小点位均是 No.1 和 No.8,且外圈河道的流场分布均要比中心水塘密集。(2)开启西南侧河道的水泵和西部泵闸,能在较短时间内提升河网水体流速。组合方案不仅流场分布最均匀、范围最广,而且对中心水塘的驱动效果最佳,对河网整体水动力的提升效果最好。引排水方案除提升河网流速外,还可以对河网水体进行置换,若河网出现水体恶化时,则可以考虑使用引排水方案提升河网水动力,功能更为全面。(3)所建立的二维河网水动力模型精确度较高,模拟状况更加直观清晰,符合平原河网水体运动规律,充分体现了对平原河网水动力调控方法研究的科学性,为未来平原河网水体的治理提供了一定的参考价值。参考文献(References):1 陈建标,钱小娟,朱友银,等.南通市引江调水对河网水环境改善效果的模拟J.水资源保护,2014,30(1):38-42.CHEN Jianbiao,QIAN Xiaojuan,ZHU Youyin,et al.Simulation of improvement of water environment in river network of nantong city by water diversion fromyangtze riverJ.Water Resources Protection,2014,30(1):38-42.2 GOLDBERG E D,GRIFFIN J J,HODGE V,et al.Pollution history of thesavannah river estuaryJ.Environmental Science&Technolo-gy,2002,13(5):588-594.3 WOITKE P,WELLMITZ J,HELM D,et al.Analysis and assess-ment of heavy metal pollution in suspended solids and sediments of the river danubeJ.Chemosphere,2003,51(8):632-642.4 PENG Z,PANG Y,PAN H,et al.Factors contributing to hypoxia in the minjiang river estuary,southeast chinaJ.International Journal of Environmental Research and Public Health,2015,12(8):9357-9374.5 张力小,张鹏鹏,郝岩,等.城市食物-能源-水关联关系:概念框架与研究展望.生态学报,2019,39(4):1143-1153.ZHANG Lixiao,ZHANG Pengpeng,HAO Yan,et al.Urban food-energy-water(FEW)nexus:conceptual frameworks and prospects.Acta Ecologica Sinica,2019,39(4):1143-1153.6 夏玉宝,王华,何新辰,等.太湖流域典型滨湖河网水动力与水质时空异质性J.湖泊科学,2021,33(4):1100-1111.XIA Yubao,WANG Hua,HE Xinchen,et al.Spatiotemporal het-erogeneity of hydrodynamic forces and water quality in typical lakeside river networks in taihu basinJ.Journal of Lake Sciences,2021,33(4):1100-1111.7 管志军.水环境质量影响因素及水生态环境保护措施研究J.皮革制作与环保科技,2021,2(18):83-84.GUAN Zhijun.Study on influencing factors of water environment quality and protection measures of water ecological environmentJ.Leather Manufacture and Environmental Technology,2021,2(18):83-84.961徐存东,等/基于 MIKE 21 的圩区河网水动力调控方法研究水利水电技术(中英文)第 54 卷 2023 年第 7 期8 王若冰.小流域水环境承载力评估预警分析J.环境保护科学,2022,48(1):115-119.WANG Ruobing.Early warning analysis of water environment carry-ing capacity assessment in small watershedJ.Environmental Pro-tection Science,2022,48(1):115-119.9 王俊钗.强干扰区域水动力水环境数值模拟研究D.武汉:武汉大学,2017.WANG Junchai.Numerical simulation on hydrodynamics and water environment in strongly disturbed regionD.Wuhan:Wuhan Uni-versity,2017.10 王鑫东.昌图县浮桥工程对辽河水流流场的影响分析J.水利发展研究,2016,16(5):55-57.WANG Xindong.Analysis of the impact of Changtu County pontoon project on the flow field of Liaohe RiverJ.Water Resources Devel-opment Research,2016,16(5):55-57.11 BOUCHARD D,KNIGHTES C,CHANG X J,et al.Simulating multiwalled carbon nanotube transport in surface water systems using the water quality analysis simulation program(WASP)J.Environ-mental Science&Technology,2017,51(19):11174-11184.12 AJIWIBOWO H,PRATAMA M B.Hydrodynamic changes impacted by the waterway capital dredging in cika