CN113642274B - 一种基于流场模型的河流流量计算方法 - Google Patents
一种基于流场模型的河流流量计算方法 Download PDFInfo
- Publication number
- CN113642274B CN113642274B CN202111198670.5A CN202111198670A CN113642274B CN 113642274 B CN113642274 B CN 113642274B CN 202111198670 A CN202111198670 A CN 202111198670A CN 113642274 B CN113642274 B CN 113642274B
- Authority
- CN
- China
- Prior art keywords
- river
- flow
- section
- flow velocity
- river channel
- 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
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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/40—Protecting water resources
- Y02A20/402—River restoration
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于流场模型的河流流量计算方法,包括河道断面测绘:对河道的边界条件进行测绘,结合数值积分算法,对各个断面网格进行重构拟合,构建拟合断面;河流流速CFD仿真:利用流体力学对河道断面流速进行仿真,得到河道断面各点的流速分布;重构流速场:采用径向神经网络的方法重构建立流速场;修正流速场:采用基于修正函数的有限元分析方法对重构的流速场进行修正;河道流量测算:通过修正后的流速场和拟合断面,两者共同构成了流场模型,结合现场监测点数据和流场模型,反演出整个河道断面的流速分布及断面面积。本发明相较于现有技术,通过测绘、仿真、重构、反演能够因地制宜的实现河流流量的动态监测,极大节省了人力成本。
Description
技术领域
本发明涉及水资源管理技术领域,尤其涉及一种基于流场模型的河流流量计算方法。
背景技术
流量是河流的重要水文特征,是反映水资源的基本资料,也是小水电站可行性研究的必要的资料。
河流流量测量是水文工作者的重要任务之一,传统的河流流量测量方法包括人工船测、桥测、缆道测量和涉水测量。其基本原理是在测流断面上布设多条垂线,在每条垂线处测量水深并用流速仪测量一至几个点的流速从而得到垂线平均流速,这种传统方法费工费时,效率低。
在河流流量测试过程中,一般来讲,测试人员多是运用转子式的流速仪进行河流流速以及水深的测量,并利用断面面积将测量结果转化为流量的真实数据。在运用上述两种河流流量测量方法时,测量人员应先行修建缆道,盖房观测,牵高、低压线等,整个过程下来耗费时间较长,工作量较大,成本投入较高,且难以实现在线观测。针对这种情况,作为水资源发达的身份,防洪预警工作较为艰难。
实际的河流断面流量测量中,由于河道断面的不规整,水深随着断面位置和河岸距离变化,且流速的变化也更加复杂,这都使得河流流量测量十分棘手。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于流场模型的河流流量计算方法。
本发明的目的是通过以下技术方案来实现的:
一种基于流场模型的河流流量计算方法,包括:
河道断面测绘:将河道断面划分为多个断面网格,对河道的边界条件进行测绘,结合数值积分算法,对各个断面网格进行重构拟合,构建拟合断面;
河流流速CFD仿真:结合河流流速、水位和河道边界条件的测绘结果,利用流体力学对河道断面流速进行仿真,得到河道断面各点的流速分布;
重构流速场:在利用流体力学仿真得到的流速场数据的基础之上,采用径向神经网络的方法重构建立流速场;
修正流速场:以测绘阶段积累的流速数据以及河流的历史数据为依据,采用基于修正函数的有限元分析方法对重构的流速场进行修正;
河道流量测算:通过修正后的流速场和拟合断面,两者共同构成了流场模型,结合现场监测点数据和流场模型,反演出整个河道断面的流速分布及断面面积。
所述划分河道断面具体包括:根据河道宽度,将河道等间距划分,并在划分处设立垂线,使得河道横截面被划分成多个网格。
所述边界条件测绘的目标是得到河流水质情况、河道坡度和糙率。
所述数值积分算法包括梯形公式、辛普森公式、中矩形公式和高斯公式。
所述河道断面流速仿真采用流体力学中的3D湍流计算模型。
所述河道流量测算具体步骤为:沿河道横向布置k个监测点,监测此处的河流水位及其表面流速数据,i=1、2、...、k;通过河流水位、表面流速和测绘阶段划分的n个网格断面,利用建立的流场模型得到各网格断面的平均流速,网格断面面积,从而采用面积包围法计算断面流量:
本发明的有益效果:
首先相较于传统的人工测绘的流速面积法,本发明在传统基础上只需要一次的测绘,后续采用CFD和径向神经网络相结合的方式可以实现流量监测。
其次相较于现有技术采用的固定流速推演公式,本发明在利用CFD仿真的基础上采用径向神经网络进行重构,可以适应多种工况和环境,且后续还可以根据实际数据进行动态修正。
最后,相较于图像测流的方式而言,每一次测量都需要布置测流粒子来获取图像进行测流分析,且需要购买特定的摄像设备,本发明的优势在于和传统的测流方式相结合,成本降低,同时也可实现流量监测。
本发明相较于现有技术,通过测绘、仿真、重构、反演能够因地制宜的实现河流流量的动态监测,极大的节省了人力成本。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图示出的结构获得其他的附图。
图1是本发明的流程框图;
图2是本发明的河道断面网格划分图。
具体实施方式
应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,一种基于流场模型的河流流量计算方法,包括:
河道断面测绘:人为设置垂线将河道断面划分为多个断面网格,对河道的边界条件工况进行测绘,结合数值积分算法,对各个断面网格进行重构拟合,构建一个目标河道的及时拟合断面;
在实际的河流断面流量测量中,区别于传统的明渠测流,河道断面往往是不规则的断面。基于此,利用测量装置从河道的左岸开始每隔一定的距离设立垂线进行此垂线河流流速和水位的测量一直到河道的右岸;这些垂线将河道断面划分为了多个网格。最后对河道的边界条件工况进行测绘。
测绘过程依据河道的宽度施测n条横向水平测线,并记录其到河道一岸的距离;得到n个水位数据,并在每条水平测线出摘取m条垂线流速数据,其中i=1、2、3...n,j=1、2、3...m。边界条件需要测绘得到河流水质情况,河道的大致坡度、糙率等。(距离和水位单位为m,流速单位m/s)
垂线将河道断面划分为了多个断面网格,基于垂线的水位数据,结合适应的数值积分算法:梯形公式、梯形公式、辛普森公式、中矩形公式、高斯法等,可以对各断面网格分别进行重构拟合。
由于这些断面网格代表了河道的整个断面情况,从而可以构建出一个目标河道的近似拟合断面。
网格断面的重构拟合要依据实际的工况确定:梯形公式和中矩形公式适合于断面网格较为规整的情况;辛普森公式适合于断面网格不规整、水深呈现线性变化的情况;高斯公式则很大程度取决于权重系数和测点情况的选择。
河流流速CFD仿真:结合测绘阶段得到的流速、水位数据以及河道边界工况的测绘,利用流体力学对河道断面流速进行仿真,得到河道断面各点的流速分布;
对于规则断面的水力学问题的处理比较容易, 而天然河流断面形状很复杂, 水深随断面位置和与河岸的距离而变化, 流速变化就更复杂,它是一种三维动态函数, 其大小和方向都随位置和时间变化, 为了解决这种复杂水流的流场结构问题, 最好是采用流体力学中的3D湍流计算模型。
基于测绘阶段得到的流速和水位数据,以及河道边界工况的测绘,我们采用计算流体力学(CFD)的方法对河道的断面流速分布进行仿真以此来得到河道断面各点的流速分布。
重构流速场:利用流体力学仿真得到的流速场数据的基础之上,采用径向神经网络的方法重构建立流速场;
天然河道中由于流速变化的复杂性,单一的指数分布或对数分布规律不足以描述整体的流速分布情况,在利用流体力学仿真得到的流速场数据的基础之上,采用径向神经网络的方法重构建立流速场。
修正流速场:以测绘阶段积累的流速数据以及河流的历史数据为依据,采用基于修正函数的有限元分析方法对重构的流速场进行修正;
以测绘阶段积累的流速数据以及河流的历史数据为依据,采用基于修正函数的有限元分析方法对重构的流速场进行修正,使其能够更加的接近于实际的流速场。
河道流量测算:通过修正的流速场和拟合得到的河道断面,两者共同构成了流场模型,结合现场监测点数据和建立的流场模型,反演出整个河道断面的流速分布及断面面积。
如图2所示,所述划分河道断面具体包括:根据河道宽度,将河道等间距划分,并在划分处设立垂线,使得河道横截面被划分成多个网格。
所述边界条件测绘的目标是得到河流水质情况、河道坡度和糙率。
所述数值积分算法包括梯形公式、辛普森公式、中矩形公式和高斯公式。
所述河道断面流速仿真采用流体力学中的3D湍流计算模型。
在经历了上述的测绘-仿真-重构-修正四个阶段,可以得到一个修正的流速场和拟合得到的河道断面,两者共同构成了流场模型。沿河道横向布置k个监测点(监测点的选取应具有代表性,如河道中部和两侧),监测此处的河流水位及其表面流速数据,i=1、2、...、k;通过河流水位、表面流速和测绘阶段划分的n个网格断面,利用建立的流场模型得到各网格断面的平均流速,网格断面面积,从而采用面积包围法计算断面流量:
首先相较于传统的人工测绘的流速面积法,本专利在其基础上只需要一次的测绘,后续采用CFD和径向神经网络相结合的方式可以实现流量监测。
其次相较于部分专利采用的固定流速推演公式,本专利在利用CFD仿真的基础上采用径向神经网络进行重构,可以适应多种工况和环境,且后续还可以根据实际数据进行动态修正。
最后,相较于图像测流的方式而言,每一次测量都需要布置测流粒子来获取图像进行测流分析,且需要购买特定的摄像设备,本专利的优势在于和传统的测流方式相结合,成本降低,同时也可实现流量监测
本发明相较于现有技术,通过测绘、仿真、重构、反演能够因地制宜的实现河流流量的动态监测,极大的节省了人力成本。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护的范围由所附的权利要求书及其等效物界定。
Claims (6)
1.一种基于流场模型的河流流量计算方法,其特征在于,包括:
河道断面测绘:将河道断面划分为多个断面网格,对河道的边界条件进行测绘,结合数值积分算法,对各个断面网格进行重构拟合,构建拟合断面;
河流流速CFD仿真:结合河流流速、水位和河道边界条件的测绘结果,利用流体力学对河道断面流速进行仿真,得到河道断面各点的流速分布;
重构流速场:在利用流体力学仿真得到的流速场数据的基础之上,采用径向神经网络的方法重构建立流速场;
修正流速场:以测绘阶段积累的流速数据以及河流的历史数据为依据,采用基于修正函数的有限元分析方法对重构的流速场进行修正;
河道流量测算:通过修正后的流速场和拟合断面,两者共同构成了流场模型,结合现场监测点数据和流场模型,反演出整个河道断面的流速分布及断面面积。
2.根据权利要求1所述的一种基于流场模型的河流流量计算方法,其特征在于,划分所述河道断面具体包括:根据河道宽度,将河道等间距划分,并在划分处设立垂线,使得河道横截面被划分成多个网格。
3.根据权利要求1所述的一种基于流场模型的河流流量计算方法,其特征在于,测绘所述河道的边界条件的目标是得到河流水质情况、河道坡度和糙率。
4.根据权利要求1所述的一种基于流场模型的河流流量计算方法,其特征在于,所述数值积分算法包括梯形公式、辛普森公式、中矩形公式和高斯公式。
5.根据权利要求1所述的一种基于流场模型的河流流量计算方法,其特征在于,仿真所述河道断面流速采用流体力学中的3D湍流计算模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111198670.5A CN113642274B (zh) | 2021-10-14 | 2021-10-14 | 一种基于流场模型的河流流量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111198670.5A CN113642274B (zh) | 2021-10-14 | 2021-10-14 | 一种基于流场模型的河流流量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113642274A CN113642274A (zh) | 2021-11-12 |
CN113642274B true CN113642274B (zh) | 2022-01-14 |
Family
ID=78426880
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111198670.5A Active CN113642274B (zh) | 2021-10-14 | 2021-10-14 | 一种基于流场模型的河流流量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113642274B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114896909A (zh) * | 2022-05-23 | 2022-08-12 | 四川大学 | 一种基于水位高度的明渠流量计算方法 |
CN115600044B (zh) * | 2022-11-28 | 2023-05-05 | 湖南大学 | 一种河流断面流量计算方法、装置、设备及存储介质 |
CN116432551B (zh) * | 2023-03-30 | 2023-12-15 | 西北农林科技大学 | 基于卫星图像和cfd的河流交汇区水质分布仿真方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102087130A (zh) * | 2010-11-19 | 2011-06-08 | 北京工业大学 | 基于cfd技术的多声路超声流量计弯管安装声路的优化方法 |
CN107421853A (zh) * | 2017-04-26 | 2017-12-01 | 青岛理工大学 | 一种河流横向扩散系数的等浓度线快速观测方法 |
CN111581828A (zh) * | 2020-05-11 | 2020-08-25 | 福建省水利水电勘测设计研究院 | 一种感潮河段闸下水位流量关系的计算方法 |
CN113124941A (zh) * | 2021-04-01 | 2021-07-16 | 长江水利委员会长江科学院 | 一种河道流量非接触式测量及精确计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7457735B2 (en) * | 2001-11-14 | 2008-11-25 | Bentley Systems, Incorporated | Method and system for automatic water distribution model calibration |
-
2021
- 2021-10-14 CN CN202111198670.5A patent/CN113642274B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102087130A (zh) * | 2010-11-19 | 2011-06-08 | 北京工业大学 | 基于cfd技术的多声路超声流量计弯管安装声路的优化方法 |
CN107421853A (zh) * | 2017-04-26 | 2017-12-01 | 青岛理工大学 | 一种河流横向扩散系数的等浓度线快速观测方法 |
CN111581828A (zh) * | 2020-05-11 | 2020-08-25 | 福建省水利水电勘测设计研究院 | 一种感潮河段闸下水位流量关系的计算方法 |
CN113124941A (zh) * | 2021-04-01 | 2021-07-16 | 长江水利委员会长江科学院 | 一种河道流量非接触式测量及精确计算方法 |
Non-Patent Citations (4)
Title |
---|
北部湾经济区河流环境流量计算;陈毅 等;《中国农村水利水电》;20120131(第01期);8-17 * |
梯形断面明渠流速分布的研究;胡云进 等;《浙江大学学报(工学版)》;20090630;第43卷(第06期);1102-1106 * |
河流流场数值计算及其仿真研究;潘宇璐;《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》;20181115(第11期);C036-30 * |
黄河宁夏典型河段及水库水动力与水质数值模拟研究;黄凌霄;《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》;20200215(第02期);A004-16 * |
Also Published As
Publication number | Publication date |
---|---|
CN113642274A (zh) | 2021-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113642274B (zh) | 一种基于流场模型的河流流量计算方法 | |
EP4216117B1 (en) | Method and apparatus for training water-and-sediment prediction model for reservoir and method and apparatus for predicting water-and-sediment in reservoir | |
CN114997541B (zh) | 一种基于数字孪生技术的城市内涝预测方法及预警平台 | |
Zhong et al. | Tidal energy fluxes and dissipation in the Chesapeake Bay | |
CN108334660B (zh) | 一种基于数据同化的强冲积河流的水沙预测方法和系统 | |
JP6736731B2 (ja) | エネルギー過程の類似性に基づく水文モデル流域スケール決定方法 | |
CN108021780B (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN110646867A (zh) | 一种城市排水监测预警方法和系统 | |
CN112257352A (zh) | 一维水动力模型和二维水动力模型的耦合方法及系统 | |
CN112270115B (zh) | 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 | |
CN112464584A (zh) | 自由表面流的水位和流量推求方法 | |
CN115496015B (zh) | 一种基于流量梯度变化的水动力分析决策方法 | |
CN111428942A (zh) | 基于变网格技术提取微地形因子的线路覆冰厚度预测方法 | |
CN111241757A (zh) | 基于计算流体力学的铀尾矿库溃坝三维数值模拟方法 | |
CN115115262A (zh) | 一种洪水风险灾害评估的方法 | |
CN114580283B (zh) | 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法 | |
CN115290148A (zh) | 一种复杂水流条件下的流量测量方法及系统 | |
CN108196951B (zh) | Gpu流域径流模拟分布式调度系统及方法 | |
CN114254526A (zh) | 一种湖泊水-气界面二氧化碳交换量的评估方法及系统 | |
CN113642259B (zh) | 一种基于多层流速测量的河流流场重构方法 | |
CN117195603B (zh) | 基于高分辨率遥感要素的洪涝灾害推演方法、设备及介质 | |
CN114492984A (zh) | 粉尘浓度的时空分布预测方法、装置、设备和存储介质 | |
CN110847112B (zh) | 一种基于水力学模拟的河道泄洪预警方法 | |
CN113869804B (zh) | 一种洪涝灾害下的电网设备风险预警方法及系统 | |
Othata et al. | Irrigation water management strategies for salinity control in the Chao Phraya River using Sualyev finite difference method with Lagrange interpolation technique |
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 |