CN109710959A - 一种水资源引清调度水体置换效果数值模拟方法 - Google Patents
一种水资源引清调度水体置换效果数值模拟方法 Download PDFInfo
- Publication number
- CN109710959A CN109710959A CN201810968490.2A CN201810968490A CN109710959A CN 109710959 A CN109710959 A CN 109710959A CN 201810968490 A CN201810968490 A CN 201810968490A CN 109710959 A CN109710959 A CN 109710959A
- Authority
- CN
- China
- Prior art keywords
- water
- concentration
- equation
- mathematical model
- water substitution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开的一种水资源引清调度水体置换效果数值模拟方法,包括如下步骤:1)基于引清调度的以新换旧基本原理,遴选并定义表征水体置换效果的水体置换率;2)通过选用质量百分比浓度来表达水体浓度,并经合理假定与简化,建立水体浓度与水体置换率的数值上的等同关系;3)建立包含连续方程、运动方程、对流扩散方程在内的流体力学基本方程;4)分析引清调度关注区域所处的河网片区的特点,确定数学模型的建模范围及维度,搭建水动力与对流扩散耦合数学模型,设置对流扩散模块的初始浓度和边界浓度;5)运行计算水动力与对流扩散耦合数学模型,分析研究区域内不同时间、不同位置的水体置换情况。
Description
技术领域
本发明涉及水力学技术领域,尤其涉及一种水资源引清调度水体 置换效果数值模拟方法。
背景技术
水资源引清调度是利用外海潮汐动力与清洁水源,通过泵闸等水 工建筑物的调度,使得河网内主要河道水体定向有序流动,从而达到 加快水体更新、改善内河水质的目的。通过水资源引清调度改善内河 水质的方式自上海于20世纪80年代中期开创先例以来,已陆续被国 内诸多地区所采纳,成为提高现有河网水环境承载能力的一种重要手 段。在水资源引清调度工作中,水质改善效果历来是各方关注的焦点 问题,目前主要通过水质数学模型模拟或现场水质监测等手段分析研 究,但是这些方法存在所需要的水质基础数据众多、工作量大、耗时 长、费用高等缺点,因此有必要提出一种更为简便易行的研究方法来满足实际工作的需要。
发明内容
本发明的目的在于:提出一种水资源引清调度水体置换效果数值 模拟方法。
为了实现上述目的,本发明采用的技术方案如下:
一种水资源引清调度水体置换效果数值模拟方法,包括如下步骤:
步骤1,基于引清调度的以新换旧基本原理,遴选并定义表征水 体置换效果的水体置换率;
步骤2,通过选用质量百分比浓度来表达水体浓度,并经合理假 定与简化,建立水体浓度与水体置换率的数值上的等同关系;
步骤3,建立包含连续方程、运动方程、对流扩散方程在内的流 体力学基本方程;
步骤4,分析引清调度关注区域所处的河网片区的特点,确定数 学模型的建模范围及维度,并通过步骤3建立的流体力学基本方程搭 建水动力与对流扩散耦合数学模型,根据步骤2对水体浓度表达方式 的选用设置对流扩散模块的初始浓度和边界浓度;
步骤5,运行计算步骤4所搭建的水动力与对流扩散耦合数学模 型,根据步骤2建立的等同关系分析研究区域内不同时间、不同位置 的水体置换情况。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)基于水资源引清调度的以新换旧基本原理,推荐了“水体 置换率”这一简便指标来表征引清调度水体置换效果,简单明了且具 有代表性。
(2)通过选用“质量百分比浓度”来表达水体浓度,建立水体 浓度与水体置换率数值上的等同关系,基于水动力模型与对流扩散数 学模型的耦合计算实现了水体置换率数值的自动模拟,简便易操作。
附图说明
图1是本发明的流程示意图。
图2是应用案例位置及一维河网数学模型计算范围的示意图。
图3是三次引排水后研究区域及附近河道水体置换率计算结果 分布的示意图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
参见图1,图中给出的是一种水资源引清调度水体置换效果数值 模拟方法,包括如下步骤:
步骤1,基于引清调度的以新换旧基本原理,遴选并定义表征水 体置换效果的水体置换率。具体如下:
水资源引清调度的基本原理是“以新换旧”,即利用外部“新” 的清洁水体替换河网内“旧”的水体。基于该原理,推荐采用“水体 置换率”作为表征水体置换效果的简便指标,将其定义为“新水体的 占比”,即某处来自外部水源的新水体占总水体的比例。水体置换率 越接近于1则说明水体置换效果越好。水体置换率基本公式如下:
式(1)中:f为水体置换率;V新为某处来自“新”水体的体积; V旧为某处“旧”水体的体积。
步骤2,通过选用质量百分比浓度来表达水体浓度,并经合理假 定与简化,建立水体浓度与水体置换率的数值上的等同关系。具体如 下:
一般而言,溶液浓度可分为质量百分比浓度、体积浓度(如摩尔 浓度、当量浓度)和质量-体积浓度三类。在水力学领域,常用的是 质量-体积浓度,如悬沙浓度、水质指标等。鉴于水资源引清调度中 “新”、“旧”水体均是液体,提出选用“质量百分比浓度”来表达水 体浓度。通过选用“质量百分比浓度”来表达 水体浓度,并经合理假定与简化,就建立了水体浓度与水体置换率数 值上的等同关系。
步骤3,建立包含连续方程、运动方程、对流扩散方程在内的流 体力学基本方程。具体如下:
(1)一维模型基本方程:
①连续方程
式(2)中:A为过水断面面积;Q为过水断面流量;x为沿程坐 标;t为时间坐标。
②水流运动方程:
式(3)中:h为断面水位;C为谢才系数;R为水力半径。
③对流扩散方程:
式(4)中:C为浓度;D为扩散系数;K为线性降解系数;C2为 源/汇项浓度;q为侧向入流。
(2)二维模型基本方程
①水流连续方程
式(5)中:x、y为直角坐标系坐标;t为时间;h为平均水深; ζ为相对于平均海平面的潮位;Ux、Uy分别为x、y方向上的垂线平 均速度。
②水流运动方程:
式(6)中:ρ为水体密度;g为重力加速度;Nx、Ny分别为x、 y方向的水平紊动粘性系数;f为科氏参数;τx、τy分别为波流共同 作用下床面剪切应力在x、y方向的分量。
③对流扩散方程:
式(7)中:c为垂线平均浓度;Kx、Ky分别为x、y方向的扩 散系数。
步骤4,分析引清调度关注区域所处的河网片区的特点,确定数 学模型的建模范围及维度,并通过步骤3建立的流体力学基本方程搭 建水动力与对流扩散耦合数学模型,根据步骤2对水体浓度表达方式 的选用设置对流扩散模块的初始浓度和边界浓度。具体如下:
河网地区水系相互连通,数学模型建模范围需要足够大以消除开 边界对模拟结果的影响,一般应选取相对独立的水利片区作为建模范 围。根据经验,水资源引清调度研究中最常用的数值模拟维度为一维 或二维,当研究区域内河道众多且泵闸等水工建筑物运行调度方式相 对复杂时建议选用一维数值模拟手段,反之则可选用二维,也可根据 实际需要选用一维与二维相耦合的方式。以MIKE11或MIKE21为工具, 基于步骤3的基本方程搭建一维或二维水动力与对流扩散相耦合的 数学模型。建模过程中需要注意以下几个问题:
(1)需采用合适的空间步长或网格尺寸以便能较为准确的刻画 模型范围内的河道特征。
(2)需要对时间步长、糙率、扩散系数等关键参数进行合理赋 值。
(3)需要选择模拟引清调度的合理时间段,水动力模块需要设 置模拟时间段内的开边界条件、初始条件(其中初始条件应根据引清 调度规则合理赋值)。
(4)需要根据引清调度规则给泵闸等水工建筑物设置合理的调 度运行规则。
(5)需要设置对流扩散模块开边界入流浓度、出流方式及模型 范围内初始浓度。根据步骤2对水体浓度表达方式的选用,可将开边 界处“新”水体的入流质量百分比浓度设置为100%,河网内“旧” 水体的初始质量百分比浓度设置为0,开边界处出流设置为自由出流 方式。
步骤5,运行计算步骤4所搭建的水动力与对流扩散耦合数学模 型,根据步骤2建立的等同关系分析研究区域内不同时间、不同位置 的水体置换情况。具体如下:
运行计算步骤4搭建的一维或二维水动力与对流扩散耦合数学 模型。根据步骤2建立的等同关系,对流扩散模型计算得到的水体浓 度结果在数值上即等于水体置换率,因此可以直接根据对流扩散模型 计算得到的浓度分析研究区域内不同时间、不同位置的水体置换情况。
针对上海市崇明岛绿华镇水资源引清调度水体置换效果的数值 模拟对本发明方法进行分析说明。
步骤1,基于引清调度“以新换旧”的基本原理,遴选并定义“水 体置换率”作为表征水体置换效果的简洁指标,如公式(1)所示。
步骤2,通过选用“质量百分比浓度”来表达水体浓度,并经合 理假定与简化,建立水体浓度与水体置换率数值上的等同关系。
步骤3,建立包含连续方程、运动方程、对流扩散方程在内的一 维流体力学基本方程,如公式(2)至公式(4)或公式(5)至公式 (7)所示。
步骤4,分析引清调度关注区域所处的河网片区的特点,确定数 学模型的建模范围及维度,搭建水动力与对流扩散耦合数学模型,设 置对流扩散模块的初始浓度和边界浓度。由于案例绿华镇位于上海市 崇明岛水利片,建议模型范围涵盖整个崇明岛。参见图2,崇明岛河 网密布、引排水口门众多,其中仅绿华镇就涉及市级河道1条(环岛 运河)、区级河道1条(新建港)、镇级河道27条,另有2个距离较 近的引水口门(崇西水闸和新建水闸)、1个距离较近的排水口门(跃 进水闸),研究手段选用一维数值模拟技术。以MIKE11为工具,搭建 崇明岛包含所有镇级河道在内的一维水动力与对流扩散耦合数学模 型,设置边界条件与初始条件,根据实测资料率定糙率、扩散系数等 关键参数。模型中各水动力开边界设置为2009年9月3日~9月5 日的潮位过程。根据崇明岛引清调度中河道平均控制水位为2.60-3.00m的要求以及沿江水闸“南引北排”的调度运行规则,河 网初始水位设置为2.60m,南支沿线水闸调度规则按“能引则引”、 当控制河道水位高于3m时关闸控制;北支沿线水闸按“能排则排” 控制。南支沿线开边界入流质量百分比浓度赋值为100%、北支沿线开边界出流设置为自由出流方式,模型范围内河网水体的初始质量百 分比浓度赋值为0。
步骤5,运行计算步骤4搭建的崇明岛一维水动力与对流扩散耦 合数学模型,并根据对流扩散模型计算得到的浓度分析研究区域内不 同时间、不同位置的水体置换情况,如图3所示。
Claims (1)
1.一种水资源引清调度水体置换效果数值模拟方法,其特征在于,包括如下步骤:
步骤1,基于引清调度的以新换旧基本原理,遴选并定义表征水体置换效果的水体置换率;
步骤2,通过选用质量百分比浓度来表达水体浓度,并经合理假定与简化,建立水体浓度与水体置换率的数值上的等同关系;
步骤3,建立包含连续方程、运动方程、对流扩散方程在内的流体力学基本方程;
步骤4,分析引清调度关注区域所处的河网片区的特点,确定数学模型的建模范围及维度,并通过步骤3建立的流体力学基本方程搭建水动力与对流扩散耦合数学模型,根据步骤2对水体浓度表达方式的选用设置对流扩散模块的初始浓度和边界浓度;
步骤5,运行计算步骤4所搭建的水动力与对流扩散耦合数学模型,根据步骤2建立的等同关系分析研究区域内不同时间、不同位置的水体置换情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810968490.2A CN109710959B (zh) | 2018-08-23 | 2018-08-23 | 一种水资源引清调度水体置换效果数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810968490.2A CN109710959B (zh) | 2018-08-23 | 2018-08-23 | 一种水资源引清调度水体置换效果数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109710959A true CN109710959A (zh) | 2019-05-03 |
CN109710959B CN109710959B (zh) | 2023-08-18 |
Family
ID=66253671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810968490.2A Active CN109710959B (zh) | 2018-08-23 | 2018-08-23 | 一种水资源引清调度水体置换效果数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109710959B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539153A (zh) * | 2020-04-07 | 2020-08-14 | 长江水利委员会长江科学院 | 一种基于预构泥沙信息库的水沙联合优化调度方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0946431A1 (en) * | 1996-11-15 | 1999-10-06 | K-V Associates, Inc. | Groundwater remediation with microporous diffusor |
CN106202618A (zh) * | 2016-06-24 | 2016-12-07 | 珠江水利委员会珠江水利科学研究院 | 一种工程调度与感潮河网污染物输移过程耦合的数值模拟方法 |
CN106817904A (zh) * | 2014-02-27 | 2017-06-09 | 博德研究所 | T细胞平衡基因表达、物质组合物及其使用方法 |
CN107516172A (zh) * | 2017-09-04 | 2017-12-26 | 珠江水利委员会珠江水利科学研究院 | 咸潮影响地区抢淡蓄淡应急供水调度方法 |
-
2018
- 2018-08-23 CN CN201810968490.2A patent/CN109710959B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0946431A1 (en) * | 1996-11-15 | 1999-10-06 | K-V Associates, Inc. | Groundwater remediation with microporous diffusor |
CN106817904A (zh) * | 2014-02-27 | 2017-06-09 | 博德研究所 | T细胞平衡基因表达、物质组合物及其使用方法 |
CN106202618A (zh) * | 2016-06-24 | 2016-12-07 | 珠江水利委员会珠江水利科学研究院 | 一种工程调度与感潮河网污染物输移过程耦合的数值模拟方法 |
CN107516172A (zh) * | 2017-09-04 | 2017-12-26 | 珠江水利委员会珠江水利科学研究院 | 咸潮影响地区抢淡蓄淡应急供水调度方法 |
Non-Patent Citations (1)
Title |
---|
哈佳;顾霜妹;孟潇;薛海;: "《水动力-水质模型在人工湖优化设计中的应用》" * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111539153A (zh) * | 2020-04-07 | 2020-08-14 | 长江水利委员会长江科学院 | 一种基于预构泥沙信息库的水沙联合优化调度方法 |
CN111539153B (zh) * | 2020-04-07 | 2023-12-08 | 长江水利委员会长江科学院 | 一种基于预构泥沙信息库的水沙联合优化调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109710959B (zh) | 2023-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101403929B (zh) | 一种自动控制闸门泄流量的方法和系统 | |
CN108009349B (zh) | 一种河流二维水质数值模拟模型计算网格优化绘制方法 | |
CN103886135B (zh) | 基于二维非恒定流数值模型的电力工程选址方法 | |
CN109271672B (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN111950211A (zh) | 基于海流设计流速的海床基础局部冲刷深度确定方法及系统 | |
CN103149600B (zh) | 一种基于优化控制点确定水文地质参数的自动配线方法 | |
CN105631529A (zh) | 一种边水气藏见水时间预测方法 | |
CN104933268B (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN114218840A (zh) | 河口航道水沙运动及其地形演变整体建模及可视化系统 | |
CN115408887B (zh) | 一种近海电厂取水卷吸机制数值模拟与评估方法 | |
CN103774605B (zh) | 一种提高环抱式港池水体交换能力的设计方法 | |
Tung et al. | Morphological modeling of tidal inlet migration and closure | |
CN115689058A (zh) | 一种确定大型水库工程生态调度的方法 | |
CN114117952A (zh) | 水动力与基质相耦合的植被消长模型构建方法及装置 | |
CN109710959A (zh) | 一种水资源引清调度水体置换效果数值模拟方法 | |
CN109271661A (zh) | 开阔水域疏浚吹填泥沙成型的数值模拟方法 | |
CN103412991B (zh) | 一种网箱对水流影响的模拟方法 | |
CN110442997B (zh) | 一种用于指导河口生境改善工程的河口水域范围划定方法 | |
Paudel | Role of sediment in the design and management of irrigation canals: Sunsari Morang Irrigation Scheme, Nepal | |
Yi et al. | The bedform morphology of Chinese sturgeon spawning sites in the Yangtze River | |
CN116204956A (zh) | 竖缝式鱼道池室隔板结构优化方法和装置 | |
CN116187001A (zh) | 一种地浸采铀有效溶浸范围的自动提取方法 | |
Xie et al. | Mathematical model for flood routing in Jingjiang River and Dongting Lake network | |
CN114357903A (zh) | 基于河床形态参数的冲积河流阻力计算、模拟方法及装置 | |
Yang et al. | Potential causes of salinity trend in seawater overlaying a mudflat inside a semi-enclosed bay |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |