CN109271672B - 一种河-湖-泵站相互影响作用下的河道水面线计算方法 - Google Patents

一种河-湖-泵站相互影响作用下的河道水面线计算方法 Download PDF

Info

Publication number
CN109271672B
CN109271672B CN201810955378.5A CN201810955378A CN109271672B CN 109271672 B CN109271672 B CN 109271672B CN 201810955378 A CN201810955378 A CN 201810955378A CN 109271672 B CN109271672 B CN 109271672B
Authority
CN
China
Prior art keywords
pump station
flow
river
lake
pumping
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.)
Active
Application number
CN201810955378.5A
Other languages
English (en)
Other versions
CN109271672A (zh
Inventor
白凤朋
辛小康
贾海燕
李建
王超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
YANGTZE RIVER WATER RESOURCES PROTECTION SCIENCE RESEARCH INSTITUTE
China Yangtze Power Co Ltd
Original Assignee
YANGTZE RIVER WATER RESOURCES PROTECTION SCIENCE RESEARCH INSTITUTE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by YANGTZE RIVER WATER RESOURCES PROTECTION SCIENCE RESEARCH INSTITUTE filed Critical YANGTZE RIVER WATER RESOURCES PROTECTION SCIENCE RESEARCH INSTITUTE
Priority to CN201810955378.5A priority Critical patent/CN109271672B/zh
Publication of CN109271672A publication Critical patent/CN109271672A/zh
Application granted granted Critical
Publication of CN109271672B publication Critical patent/CN109271672B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Non-Positive-Displacement Pumps (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种河‑湖‑泵站相互影响作用下的河道水面线计算方法，是将泵站作为河道和湖泊的内边界条件，其中内边界条件大小由泵站的约束条件决定，利用有限体积方法数值求解河道一维水流流动和湖泊二维浅水水流流动，并将计算的水流参数同泵站的约束条件进行对比，动态修改泵站的内边界条件，成功的解决了如何求解河‑湖‑泵站三者相互影响作用下的水面线的问题。通过本发明获得的结果此方法有明确的物理意义，且能够较快速精确地计算河‑湖‑泵站三者相互影响作用下的水面线。

Description

一种河-湖-泵站相互影响作用下的河道水面线计算方法
技术领域
本发明属于计算水力学领域，尤其涉及河-湖-泵站复杂系统下的河道水面线计算方法。
背景技术
人类为了方便自己而对河流的改造会对河流形态和流态产生巨大的影响，从而在经历洪水时使河道水面线发生巨大的变化。水面线的计算是计算水力学的重要内容，是河道整治的基础工作,其推算结果直接影响到河道断面的规划设计,进而影响到河道整治的工程量和工程造价。水面线的计算结果偏低会导致堤防建设防水效果较差，洪水来临时无法起到堤防应有的作用；计算结果较高会导致堤防建设偏高，增加建设成本。总之，正确的求解河道水面线在河流治理和岸边防护等方面有着不可忽视的意义。
近年来，我国建设了多项长距离、跨流域、跨地区的长距离调水工程，例如南水北调东线工程采用梯级泵站提水方式将长江流域的水运往华北地区。受泵站约束条件和配水量等因素的变化，梯级泵站输水过程中水位、流量一直处于动态变化中；水位、流量的变化反过来又会约束泵站的抽水能力，影响河道的水流流态。因此，准确、科学的求解梯级泵站系统作用下的水面线是一个难题。目前，常用的方法是静态水量方法，不考虑泵站抽水过程中河道中水体的动态变化，没有将河-湖-泵站复杂系统耦合起来计算，计算的结果精度很低，可靠性较低。
发明内容
本发明针对现有技术的不足，提供一种考虑河-湖-泵站相互影响作用下的河道水面线计算方法。河道水流模拟采用采用Godunov格式的有限体积法求解守恒形式的一维圣维南方程组，对于不规则的河道断面，提出了采用复化Simpson公式计算天然河道下任意不规则断面的静水压力的方法。湖泊水流流动模拟采用适用于复杂地形条件的二维浅水水流高精度数学模型，根据半隐式格式的构造方法，推导了摩阻源项的半隐式离散格式，增加了数值模型的稳定性。泵站布置在河道和湖泊之间，成梯级分布；泵站的抽水量受到其约束条件限制，包括水位约束条件、扬程约束条件和流量约束条件。泵站的抽水量作为河道一维数学模型和湖泊二维浅水数学模型的内边界条件；同时，河道和湖泊数学模型计算的泵站前后水位和流量值又会约束泵站的抽水量，河-湖-泵站三者相互约束。此方法有明确的物理意义，且能够较快速精确地计算河-湖-泵站三者相互影响作用下的水面线。
本发明的技术方案为一种河-湖-泵站相互影响作用下的河道水面线计算方法，其特征在于，包含以下步骤：
步骤1，收集资料：收集河道、湖泊的地形资料和水文数据，包括水深、水位、流速、河道流量；根据实际需求确定泵站机组最大数量限制条件；根据泵站的性能确定泵站约束条件，包括泵站特征水位变化范围约束条件、泵站特征扬程变化范围约束条件；根据泵站约束条件和已有实测流量、泵站扬程资料，拟合出泵站扬程-流量变化曲线；
步骤2，数值计算：泵站将河道分为若干段，将泵站的抽水流量作为河道一维数学模型和湖泊二维浅水数学模型的内边界条件，泵站上游河道将泵站抽水流量作为下边界条件，泵站下游河道将泵站抽水流量作为上边界条件；采用数值离散方法数值模拟河道水流流动和湖泊水流流动，得到河道和湖泊的水位、流速、流量三个水流参数；
步骤3，判断泵站约束条件：根据步骤2数值模拟得到的泵站上游和下游水位计算泵站扬程，判断扬程是否满足泵站约束条件，若满足约束条件，根据扬程-流量变化曲线计算得到泵站最新的抽水流量；若不满足，根据泵站扬程约束条件的最大值计算泵站抽水流量；
步骤4，再次数值计算：根据步骤3得到的泵站抽水量作为数学模型的新的内边界条件，返回步骤2，重新采用数值离散方法求解河道水流流动和湖泊水流流动，得到新的水位、水深、流速三个水流参数；
步骤5，数据输出和处理数据步骤。
在上述的一种河-湖-泵站相互影响作用下的河道水面线计算方法，步骤2中数值离散的具体操作方法如下：
步骤2.1、将计算区域进行网格剖分；
步骤2.2、输入相关计算参数，包括计算区域长度、计算结束时间、时间步长、河床曼宁系数；
步骤2.3、初始化：给定计算区域初始值，区域初始值包括水深、流速；
步骤2.4、确定河道和湖泊入流、出流边界条件；
步骤2.5、泵站的约束条件作为模型计算的边界条件，泵站的抽水流量作为河道的内边界条件，泵站将河道分为若干个河道，泵站上游河道将泵站抽水流量作为下边界条件，泵站下游河道将泵站抽水流量作为上边界条件；
步骤2.6、采用一维数值模型对每个河道分别进行数值求解，湖泊采用二维数学模型进行数值求解，对任意网格控制体采用有限体积Godunov格式对控制方程进行积分，运用格林公式将对流通量梯度项由控制体的面积分转化为沿其边界的线积分；
河道一维圣维南控制方程组为：
Figure BDA0001772552270000041
其中，t表示时间变量，B表示河道水面宽度，z表示水位，x表示沿着河道深泓线的坐标，A表示过水断面面积；Q表示流量；qin表示旁侧入流单宽流量，若不考虑旁侧入流，可设置为零；g表示重力加速度；Sb表示河床斜率，zb表示河床高程；Sf表示沿程阻力损失，n表示曼宁糙率系数，R＝A/P表示水力半径，P表示湿周，gI1代表静水压力；
一维数值模型更新水位、流速、流量的公式如下：
Figure BDA0001772552270000042
式中，i表示控制体编号，n表示计算时间层，Δt表示计算时间步长，Δx表示控制体间距；当控制体编号i位于泵站所在位置时，需采用泵站的抽水流量作为边界条件；
湖泊二维浅水数学模型的控制方程为：
Figure BDA0001772552270000051
其中，x表示沿着河道深泓线的坐标，y表示沿河道横向方向，g表示重力加速度，zb表示河床高程；n表示曼宁糙率系数，h表示水深，u，v分别代表x、y方向上的流速，水位z＝h+zb；
二维数值模型更新水位、流速、流量的公式如下：
Figure BDA0001772552270000052
式中，i表示x方向控制体编号，j表示y方向控制体编号，n表示计算时间层，Δt表示计算时间步长，Δx表示x方向控制体间距，Δy表示y方向控制体间距,流量
Figure BDA0001772552270000053
步骤2.7、根据步骤2.3的初始条件、步骤2.4和2.5的边界条件以及2.6的控制方程，采用有限体积方法进行数值求解，计算得到河道和湖泊水深、水位、流速三个水流参数；
步骤2.8循环步骤2.4至2.7，直至满足计算时间和精度要求，输出和处理数据。
在上述的一种河-湖-泵站相互影响作用下的河道水面线计算方法，步骤1中的泵站约束条件往往包括机组数量约束、泵站特征水位变化范围约束、泵站特征扬程变化范围约束和泵站扬程-流量变化曲线约束；受到经济性、安全性和工程管理等不同因素的影响，不同的泵站对不同约束条件的响应程度不同。
在上述的一种河-湖-泵站相互影响作用下的河道水面线计算方法，步骤1中的泵站扬程-流量变化曲线可以根据泵站约束条件和泵站之前运行得到的数据进行拟合得到。
因此，本发明具有如下优点：1、成功的将泵站抽水系统嵌套至河道、湖泊数学模型中；2、根据每个时间的水面线数值计算结果，可以实时判断水流条件是否满足泵站约束条件。
附图说明
图1是本发明实施例的南水北调东线东平湖和南四湖上级湖之间的河-湖-泵站复杂系统下的水面线求解的流程图。
图2是本发明实施例的南水北调东线东平湖和南四湖上级湖之间的河-湖-泵站复杂系统示意图。
图3a是河道最后一个断面高程图。
图3b是河道第一个断面高程图。
图4是本发明实施例的东平湖和南四湖上级湖之间的实测水面线图。
图5是本发明实施例的水面线计算水面线与实测水面线对比图。
图6是本发明实施例的泵站扬程计算水面线图。
具体实施方式
本发明将泵站作为河道和湖泊的内边界条件，其中内边界条件大小由泵站的约束条件决定，利用有限体积方法数值求解河道一维水流流动和湖泊二维浅水水流流动，并将计算的水流参数同泵站的约束条件进行对比，动态修改泵站的内边界条件，成功的解决了如何求解河-湖-泵站三者相互影响作用下的水面线的问题。通过本发明获得的结果更加方便和准确。
一、首先介绍本发明的具体方法。
本发明包含以下步骤：
步骤1，给定河道和湖泊的地形资料和初始水流状态，包括水深、流速等；
步骤2，收集河道、湖泊的地形资料和水文数据，包括水深、水位、流速、河道流量；根据实际需求确定泵站机组最大数量限制条件；根据泵站的性能确定泵站约束条件，包括泵站特征水位变化范围约束条件、泵站特征扬程变化范围约束条件；根据泵站约束条件和已有实测流量、泵站扬程资料，拟合出泵站扬程-流量变化曲线；
步骤3，泵站将河道分为若干段，将泵站的抽水流量作为河道一维数学模型和湖泊二维浅水数学模型的内边界条件，泵站上游河道将泵站抽水流量作为下边界条件，泵站下游河道将泵站抽水流量作为上边界条件；采用一维数值模型对每个河道分别进行数值求解，湖泊采用二维数学模型进行数值求解，对任意网格控制体采用有限体积Godunov格式对控制方程进行积分，运用格林公式将对流通量梯度项由控制体的面积分转化为沿其边界的线积分。
河道一维圣维南控制方程组为：
Figure BDA0001772552270000071
其中，t表示时间变量，B表示河道水面宽度，z表示水位，x表示沿着河道深泓线的坐标，A表示过水断面面积；Q表示流量；qin表示旁侧入流单宽流量，若不考虑旁侧入流，可设置为零；g表示重力加速度；Sb表示河床斜率，zb表示河床高程；Sf表示沿程阻力损失，n表示曼宁糙率系数，R＝A/P表示水力半径，P表示湿周，gI1代表静水压力。
一维数值模型更新水位、流速、流量的公式如下：
Figure BDA0001772552270000081
式中，i表示控制体编号，n表示计算时间层，Δt表示计算时间步长，Δx表示控制体间距。当控制体编号i位于泵站所在位置时，需采用泵站的抽水流量作为边界条件。
湖泊二维浅水数学模型的控制方程为：
Figure BDA0001772552270000082
其中，x表示沿着河道深泓线的坐标，y表示沿河道横向方向，g表示重力加速度，zb表示河床高程；n表示曼宁糙率系数，h表示水深，u，v分别代表x、y方向上的流速，水位z＝h+zb。
二维数值模型更新水位、流速、流量的公式如下：
Figure BDA0001772552270000091
式中，i表示x方向控制体编号，j表示y方向控制体编号，n表示计算时间层，Δt表示计算时间步长，Δx表示x方向控制体间距，Δy表示y方向控制体间距,流量
Figure BDA0001772552270000092
步骤4，根据步骤3数值模拟得到的泵站上游和下游水位计算泵站扬程，判断扬程是否满足泵站约束条件，若满足约束条件，根据扬程-流量变化曲线计算得到泵站最新的抽水流量；若不满足，根据泵站扬程约束条件的最大值计算泵站抽水流量；
步骤5，根据新确定的泵站内部边界条件，重复步骤3，求解得到新的水位、水深、流速三个水流参数；
步骤6，重复步骤2至步骤5，直至满足计算时间和精度要求；
步骤7，输出和处理计算数据。
二、本发明提供的方法能够用计算机编程技术实现流程，参见图1。以南水北调东线东平湖和南四湖上级湖之间的河-湖-泵站复杂系统下的水面线求解为实施案例对本发明的流程进行一个具体的阐述。东平湖和南四湖上级湖之间存在八里湾泵站、邓楼泵站、长沟泵站三个泵站，参见图2。具体实施方式如下：
步骤1，收集东平湖、南四湖上级湖地形资料以及两者之间河道的断面地形，参见图3a和图3b；
步骤2，给定东平湖、南四湖上级湖以及两者之间河道的初始水流状态，确定三个泵站前后的水位值，参见图4；
步骤3，确定三个泵站的约束条件，包括水位约束条件、机组不同角度下扬程-流量关系曲线和流量约束条件，三个泵站水位和扬程约束条件如下表：
泵站水位和扬程约束条件
Figure BDA0001772552270000101
步骤4，根据三个泵站的扬程-流量变化曲线确定每个计算时间段的泵站抽水流量，并将其作为一维河道数学模型和二维湖泊数学模型的内边界条件，上游河道将泵站抽水流量作为下边界条件，下游河道将泵站抽水作为上边界条件；
步骤5，根据东平湖的入流边界条件、南四湖上级湖的下游边界条件和泵站的内边界条件，采用有限体积方法数值求解一维河道数学模型和二维湖泊数学模型，得到河道的水流参数，包括水位、流量、流速等；
步骤6，将泵站前后计算的水位、流量值同步骤3的泵站约束条件进行对比，根据计算值和泵站约束条件重新确定泵站的抽水流量，并将新的泵站抽水流量作为一维河道数学模型和二维湖泊数学模型的内边界条件；
步骤7，重复步骤5，采用有限体积方法数值求解一维河道数学模型和二维湖泊数学模型，得到河道新的水流参数；
步骤8，重复步骤3至步骤7直至满足计算时间条件；
步骤9，输出和处理计算数据，计算结果参见图5和图6。
具体的泵站数量在不同的工程中会各不相同，本领域技术人员可以实际情况输入模型当中。基于经济性、安全性等因素的考量，泵站的约束条件在不同时间也存在差异，本领域技术人员可以需要根据具体情况设定，需要注意的是泵站受到约束的时间需要与模型计算的时间一致。
本文中所描述的具体实施例仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施案例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (4)

1.一种河-湖-泵站相互影响作用下的河道水面线计算方法，其特征在于，包含以下步骤：
步骤1，收集资料：收集河道、湖泊的地形资料和水文数据，包括水深、水位、流速、河道流量；根据实际需求确定泵站机组最大数量限制条件；根据泵站的性能确定泵站约束条件，包括泵站特征水位变化范围约束条件、泵站特征扬程变化范围约束条件；根据泵站约束条件和已有实测流量、泵站扬程资料，拟合出泵站扬程-流量变化曲线；
步骤2，数值计算：泵站将河道分为若干段，将泵站的抽水流量作为河道一维数学模型和湖泊二维浅水数学模型的内边界条件，泵站上游河道将泵站抽水流量作为下边界条件，泵站下游河道将泵站抽水流量作为上边界条件；采用数值离散方法数值模拟河道水流流动和湖泊水流流动，得到河道和湖泊的水位、流速、流量三个水流参数；
步骤3，判断泵站约束条件：根据步骤2数值模拟得到的泵站上游和下游水位计算泵站扬程，判断扬程是否满足泵站约束条件，若满足约束条件，根据扬程-流量变化曲线计算得到泵站最新的抽水流量；若不满足，根据泵站扬程约束条件的最大值计算泵站抽水流量；
步骤4，再次数值计算：根据步骤3得到的泵站抽水量作为数学模型的新的内边界条件，返回步骤2，重新采用数值离散方法求解河道水流流动和湖泊水流流动，得到新的水位、水深、流速三个水流参数；
步骤5，数据输出和处理数据步骤。
2.根据权利要求1所述的一种河-湖-泵站相互影响作用下的河道水面线计算方法，其特征在于，步骤2中数值离散的具体操作方法如下：
步骤2.1、将计算区域进行网格剖分；
步骤2.2、输入相关计算参数，包括计算区域长度、计算结束时间、时间步长、河床曼宁系数；
步骤2.3、初始化：给定计算区域初始值，区域初始值包括水深、流速；
步骤2.4、确定河道和湖泊入流、出流边界条件；
步骤2.5、泵站的约束条件作为模型计算的边界条件，泵站的抽水流量作为河道的内边界条件，泵站将河道分为若干个河道，泵站上游河道将泵站抽水流量作为下边界条件，泵站下游河道将泵站抽水流量作为上边界条件；
步骤2.6、采用一维数值模型对每个河道分别进行数值求解，湖泊采用二维数学模型进行数值求解，对任意网格控制体采用有限体积Godunov格式对控制方程进行积分，运用格林公式将对流通量梯度项由控制体的面积分转化为沿其边界的线积分；
河道一维圣维南控制方程组为：
Figure FDA0001772552260000021
Figure FDA0001772552260000022
Figure FDA0001772552260000023
其中，t表示时间变量，B表示河道水面宽度，z表示水位，x表示沿着河道深泓线的坐标，A表示过水断面面积；Q表示流量；qin表示旁侧入流单宽流量，若不考虑旁侧入流，可设置为零；g表示重力加速度；Sb表示河床斜率，zb表示河床高程；Sf表示沿程阻力损失，n表示曼宁糙率系数，R＝A/P表示水力半径，P表示湿周，gI1代表静水压力；
一维数值模型更新水位、流速、流量的公式如下：
Figure FDA0001772552260000031
Figure FDA0001772552260000032
Figure FDA0001772552260000033
式中，i表示控制体编号，n表示计算时间层，Δt表示计算时间步长，Δx表示控制体间距；当控制体编号i位于泵站所在位置时，需采用泵站的抽水流量作为边界条件；
湖泊二维浅水数学模型的控制方程为：
Figure FDA0001772552260000034
Figure FDA0001772552260000035
Figure FDA0001772552260000036
其中，x表示沿着河道深泓线的坐标，y表示沿河道横向方向，g表示重力加速度，zb表示河床高程；n表示曼宁糙率系数，h表示水深，u，v分别代表x、y方向上的流速，水位z＝h+zb；
二维数值模型更新水位、流速、流量的公式如下：
Figure FDA0001772552260000037
Figure FDA0001772552260000038
Figure FDA0001772552260000039
式中，i表示x方向控制体编号，j表示y方向控制体编号，n表示计算时间层，Δt表示计算时间步长，Δx表示x方向控制体间距，Δy表示y方向控制体间距,流量
Figure FDA0001772552260000041
步骤2.7、根据步骤2.3的初始条件、步骤2.4和2.5的边界条件以及2.6的控制方程，采用有限体积方法进行数值求解，计算得到河道和湖泊水深、水位、流速三个水流参数；
步骤2.8循环步骤2.4至2.7，直至满足计算时间和精度要求，输出和处理数据。
3.根据权利要求1所述的一种河-湖-泵站相互影响作用下的河道水面线计算方法，其特征在于，步骤1中的泵站约束条件往往包括机组数量约束、泵站特征水位变化范围约束、泵站特征扬程变化范围约束和泵站扬程-流量变化曲线约束；受到经济性、安全性和工程管理等不同因素的影响，不同的泵站对不同约束条件的响应程度不同。
4.根据权利要求1所述的一种河-湖-泵站相互影响作用下的河道水面线计算方法，其特征在于，步骤1中的泵站扬程-流量变化曲线可以根据泵站约束条件和泵站之前运行得到的数据进行拟合得到。
CN201810955378.5A 2018-08-21 2018-08-21 一种河-湖-泵站相互影响作用下的河道水面线计算方法 Active CN109271672B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810955378.5A CN109271672B (zh) 2018-08-21 2018-08-21 一种河-湖-泵站相互影响作用下的河道水面线计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810955378.5A CN109271672B (zh) 2018-08-21 2018-08-21 一种河-湖-泵站相互影响作用下的河道水面线计算方法

Publications (2)

Publication Number Publication Date
CN109271672A CN109271672A (zh) 2019-01-25
CN109271672B true CN109271672B (zh) 2022-11-29

Family

ID=65154120

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810955378.5A Active CN109271672B (zh) 2018-08-21 2018-08-21 一种河-湖-泵站相互影响作用下的河道水面线计算方法

Country Status (1)

Country Link
CN (1) CN109271672B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109918821B (zh) * 2019-03-15 2020-01-14 中国水利水电科学研究院 一种迎风守恒型河道漫溢出流数值模拟方法
CN110046469B (zh) * 2019-05-13 2020-02-11 水利部交通运输部国家能源局南京水利科学研究院 多约束条件下水电站坝前河床冲淤变形的计算方法
CN110984088A (zh) * 2019-11-26 2020-04-10 昆明理工大学 一种竖直消力井消能率的实验计算方法
CN111259607B (zh) * 2020-01-08 2022-02-18 北京师范大学 一种河湖过渡区水文边界界定方法
CN111767684B (zh) * 2020-06-30 2024-04-19 西安理工大学 一种优化的摩阻源项隐式格式二维浅水方程建模方法
CN112597667B (zh) * 2021-01-08 2022-02-01 中国水利水电科学研究院 一种梯级泵站系统调水工程泵站优化调度方法
CN112784505B (zh) * 2021-01-29 2023-03-10 三峡大学 基于数值-解析联合求解的河道三维流场数据处理方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722611A (zh) * 2012-05-29 2012-10-10 清华大学 一种含梯级水电站河流水动力条件并行化数值模拟方法
CN104750972A (zh) * 2015-01-30 2015-07-01 长江水利委员会水文局 复杂水网地区楔形水体容积的计算方法
CN106815448A (zh) * 2017-02-07 2017-06-09 长江水资源保护科学研究所 一种河流衰减型污染物模拟方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722611A (zh) * 2012-05-29 2012-10-10 清华大学 一种含梯级水电站河流水动力条件并行化数值模拟方法
CN104750972A (zh) * 2015-01-30 2015-07-01 长江水利委员会水文局 复杂水网地区楔形水体容积的计算方法
CN106815448A (zh) * 2017-02-07 2017-06-09 长江水资源保护科学研究所 一种河流衰减型污染物模拟方法

Also Published As

Publication number Publication date
CN109271672A (zh) 2019-01-25

Similar Documents

Publication Publication Date Title
CN109271672B (zh) 一种河-湖-泵站相互影响作用下的河道水面线计算方法
CN110046469B (zh) 多约束条件下水电站坝前河床冲淤变形的计算方法
CN108009349B (zh) 一种河流二维水质数值模拟模型计算网格优化绘制方法
CN107045568B (zh) 基于动态规划逐次逼近法的河道糙率反演方法
CN104933268B (zh) 一种基于一维非恒定流数值模型的洪水分析方法
CN110135069B (zh) 输水隧洞输水时的泥沙特征获取方法、装置、计算机设备
CN104091065A (zh) 一种求解浅水问题模拟间断水流数值的方法
CN112464584A (zh) 自由表面流的水位和流量推求方法
CN106295056A (zh) 一种复式断面河道河漫滩和主槽的自动识别方法
CN107256312B (zh) 一种基于潮流环境下海湾纳潮变化量计算方法
CN106844963B (zh) 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法
Wang et al. Three-dimensional flood routing of a dam break based on a high-precision digital model of a dense urban area
CN102567594A (zh) 一种近岸岛礁型人工鱼礁群流场仿真建模的方法
Kazemzadeh-Parsi Numerical flow simulation in gated hydraulic structures using smoothed fixed grid finite element method
CN103412991B (zh) 一种网箱对水流影响的模拟方法
KR20110072551A (ko) 텐서형 와점성계수를 가진 2차원 하천흐름모형을 이용하여 천수흐름을 해석하는 방법
CN107169227A (zh) 一种分段压裂水平井的粗网格模拟方法及系统
CN114896909A (zh) 一种基于水位高度的明渠流量计算方法
Sharma Mathematical modelling and braid indicators
CN115994396A (zh) 一种水系连通数学模型的构建方法
Xie et al. Mathematical model for flood routing in Jingjiang River and Dongting Lake network
Khorrami et al. Improving multi-block sigma-coordinate for 3D simulation of sediment transport and steep slope bed evolution
CN115796077B (zh) 潮排潮灌区灌溉用水配置方法、计算机装置及存储介质
CN206616513U (zh) 兼顾床面稳定和快速排水的河工模型排水系统
Wang et al. Research and application of hydrodynamics modeling of channel in reservoir area-Case of Feilaixia station to Qingyuan station section

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
TR01 Transfer of patent right

Effective date of registration: 20230831

Address after: 430051 no.515 Qintai Avenue, Hanyang District, Wuhan City, Hubei Province

Patentee after: YANGTZE RIVER WATER RESOURCES PROTECTION SCIENCE Research Institute

Patentee after: CHINA YANGTZE POWER Co.,Ltd.

Address before: 430051 no.515 Qintai Avenue, Hanyang District, Wuhan City, Hubei Province

Patentee before: YANGTZE RIVER WATER RESOURCES PROTECTION SCIENCE Research Institute

TR01 Transfer of patent right