CN116955889A - 一种基于曼宁公式的河道枯水流量计算方法 - Google Patents
一种基于曼宁公式的河道枯水流量计算方法 Download PDFInfo
- Publication number
- CN116955889A CN116955889A CN202310893339.8A CN202310893339A CN116955889A CN 116955889 A CN116955889 A CN 116955889A CN 202310893339 A CN202310893339 A CN 202310893339A CN 116955889 A CN116955889 A CN 116955889A
- Authority
- CN
- China
- Prior art keywords
- cell
- water
- cells
- water flow
- calculating
- 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.)
- Pending
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 248
- 238000004364 calculation method Methods 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 claims abstract description 72
- 238000013480 data collection Methods 0.000 claims abstract description 4
- 230000008595 infiltration Effects 0.000 claims description 29
- 238000001764 infiltration Methods 0.000 claims description 29
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 18
- 238000012937 correction Methods 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 238000012876 topography Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000009826 distribution Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000004836 empirical method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000032823 cell division Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于曼宁公式的河道枯水流量计算方法,包括如下步骤,S1、数据收集及格栅化处理;S2、单元格水流流向计算;S3、单元格水流流速计算;S4、单元格水量交换量和单元格水位计算;S5、单元格水量修正;S6、单元格水位信息更新。优点是:提出了一种能解决河道枯水流量演算不稳定问题,同时能实现空间二维度演算的水流计算方法,该方法保持了一维度水动力学方法的简单实用,在保证水流计算效率和计算精度的基础上,解决了以往水动力学方法在小水深、干河床情况下的不适用或不稳定问题,实现了枯水期河道的水流演算,为枯水期河道补水调度等工程实践提供技术支持。
Description
技术领域
本发明涉及河道水动力技术领域,尤其涉及一种基于曼宁公式的河道枯水流量计算方法。
背景技术
随着人类活动和全球气候变化的影响,枯水期补水调度成为了流域生态管理最为重要的工作之一。最大可能节约用水、科学有效调度有限的水资源以精准满足流域生态用水需求具有现实意义,而这项工作离不开河道枯水流量的计算。因此,高效、准确模拟出河道枯水流量的演进过程是流域实施水量调度和生态保护等工作的前提。
河道水流模拟计算目前主要有经验方法、水文学方法和水动力学方法三大类。其中,经验方法依据统计相关分析,构建以线性系统为主的经验性模型;水文学方法从水量平衡原理和河道槽蓄角度出发,建立河道水流演进的半经验性模型;水动力学方法则主要通过求解圣维南方程组实现非恒定水流数值计算。其中经验方法类似于黑箱模型,对资料要求低,计算精度也相对较低,一定程度上可解释为运动波的理论求解,方法主要运用于水文模型中的产汇流计算,在流域枯水、补水调度等水流计算中很少应用;水文学方法是水动力学方法的简化,方法没有对地形资料的要求,是运动波方程的差分形式,计算远没有动力波求解复杂,但模型结果在精度上不占优势,在河道洪水演算中应用较为广泛;而水动力学方法通过求解完整的圣维南方程组计算得到流量水位时空分布信息,能综合反映各水力要素之间的动态联系,演算结果更能满足枯水期水量调度的精度需求。尽管如此,枯水期水流小、水位低、甚至有干河床情形,采用水动力学方法常出现数值不稳定、迭代不收敛等问题,即便采用一些激波捕捉技术或者窄缝方法进行处理也无法解决频繁的干湿交替问题。
发明内容
本发明的目的在于提供一种基于曼宁公式的河道枯水流量计算方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明采用的技术方案如下:
一种基于曼宁公式的河道枯水流量计算方法,包括如下步骤,
S1、数据收集及格栅化处理:
对收集或获取的河道地形数据、河床下渗数据和模型参数进行格栅化处理,以将各河道划分为多个单元格,并为各单元格进行数据输入和参数赋值;
S2、单元格水流流向计算:
基于D8算法,通过计算当前单元格与其周围单元格之间的坡降,将单元格水流方向与相邻八个单元格之间坡降最大值代表的方向作为该单元格的水流流向;
S3、单元格水流流速计算:
基于单元格的水流流向及对应坡降,利用曼宁公式计算单元格的水流流速,并对流速大的急流情况进行修正;
S4、单元格水量交换量和单元格水位计算:
分别计算单元格水量交换量所包含的单元格流入流出量和单元格渗漏水量,基于获取的单元格水量交换量,利用水量平衡原理,计算各单元格水位;
S5、单元格水量修正:
所有单元格水位计算完毕后,利用各单元格流入扣除流出和下渗水量得到的单位面积不平衡水量对所有单元格逐个进行水量平衡修正;
S6、单元格水位信息更新:
更新所有单元格水位信息,进入下一水流计算时段。
优选的,步骤S1具体包括如下内容,
S11、收集河道地形数据和河床下渗数据,并根据河道河床情况确定模型参数;
S12、将河道地形数据、河床下渗数据和模型参数进行格栅化处理,按一定分辨率将河道划分为多个正方形的单元格,并对所划分的每一个单元格进行数据输入和参数赋值。
优选的,步骤S12中,单元格的最大行数由河道沿程最大宽度决定,单元格的最大列数由河道长度决定;
M=BMAX/DS (1)
N=Len/DS (2)
其中,M为河网划分单元格的最大行数;N为河网划分单元格的最大列数;BMAX为沿着水流方向最大河道宽度;Len为河道长度;DS为河道网格划分分辨率,即所划分的单元格的边长。
优选的,步骤S2具体包括如下内容,
S21、以中心单元格为计算对象,依次将其周围的单元格进行编码;
S22、根据上一时段计算的每个单元格水量结果,计算当前第i行、第j列的单元格(i,j)与周围单元格之间的坡降;
Slope(i,j)=dZ/DS (3)
其中,Z为单元格水位,dZ为中心单元格和相邻八个单元格的水位差;单元格(i,j)与相邻单元格之间的距离与方向有关,对于位于单元格(i,j)正上、正左、正右、正下方向上的相邻单元格,单元格(i,j)与相邻单元格之间的距离为DS,坡降计算公式采用方程式(3);对于位于单元格(i,j)左上、右上、左下、右下方向上的相邻单元格,单元格(i,j)与相邻单元格之间的距离为坡降计算公式采用方程式(4);
S23、根据单元格(i,j)坡降计算结果,假设每个计算时间步长内水流方向唯一,单元格(i,j)水流方向为与该单元格与相邻八个单元格间坡度最大值代表的方向;
Aspect(i,j)=f(max Slope(i,j)) (5)
其中,Aspect(i,j)为计算单元格与周围单元格的水流交换方向代号;f()为最大坡度对应单元格指示的水流交换方向代号。
优选的,步骤S23中最大坡度对应单元格指示的水流交换方向代号f()采用如下规则确定,
若maxSlope(i,j)≤0,则赋以Aspect(i,j)为当前单元格代号值,说明当前单元格水流不流出;
若maxSlope(i,j)>0,且只有一个相邻单元格对应坡降最大,则赋以Aspect(i,j)为对应相邻单元格的编码值;
若maxSlope(i,j)>0,且有两个及以上相邻单元格对应坡降最大值,则按照顺时针方向赋以Aspect(i,j)编码号。
优选的,步骤S3具体包括如下内容,
S31、根据单元格水流流向及对应坡降,由曼宁公式计算单元格水流流速,
其中,v为单元格水流流速;ns为曼宁糙率系数;R为水力半径,近似为单元水深;
S32、根据弗洛德数判定流态的方法,对速度限定上限阀值,以对流速大的急流情况进行修正;
其中,g为重力加速度;h(i,j)为单元格(i,j)的水深。
优选的,步骤S4具体包括如下内容,
S41、基于单元格水流流速,计算单元格流出水量;
dw(i,j)=v·Δt·DS·h(i,j) (8)
其中,Δt为计算时段长;dw(i,j)为单元格(i,j)流出水量,当流向为正上、正左、正右、正下时,dw(i,j)由公式(8)计算,当流向为左上、右上、左下、右下时,dw(i,j)由公式(9)计算;
S42、基于霍顿产流公式,计算单元格渗漏水量;
Fp=[fc+(f0-fc)e-βt]h(i,j) (10)
其中,Fp为单元格渗漏水量;f0为干河床时最大下渗率,随着入渗过程进行,减小到近似常数的入渗率fc;参数β控制着入渗能力的下降速度;f0、fc和β根据实测入渗资料率定;
S43、基于单元格流出水量和单元格渗漏水量,根据水量平衡原理,计算单元格水位;
其中,Z(t)和Z(t+Δt)分别为计算时段初和时段末的单元格水位;∑dw′(i,j)为相邻单元格流入单元格(i,j)的水量。
优选的,步骤S5具体为,计算完毕所有单元格水位后,利用如下公式,对所有单元格逐个进行水量平衡修正,
Z‘(t+Δt)=Z(t+Δt)+ΔZ (12)
其中,ΔZ为每个单元格流入扣除流出和下渗水量得到的单位面积不平衡水量;Z‘(t+Δt)为平衡修正后的水位。
优选的,步骤S6具体为,完成当前时段所有单元格的水流流速、水位计算后,将相关计算结果作为下一时段的初始值,返回步骤S2。
本发明的有益效果是:1、本发明提供的河道枯水流量计算方法,解决了河道枯水流量小水深或干河床情况下的演算不稳定问题。一方面,方法保持了一维度水动力学方法的简单实用,在保证水流计算效率和计算精度的基础上,解决了以往水动力学方法在小水深、干河床情况下的不适用或不稳定问题,另一方面,方法能实现河道水流的空间二维演算,可更为准确地刻画出天然河道小水深情况下横向流速、水位的空间差异。2、本发明方法实现了枯水期河道水流的稳定演算,可为枯水河道补水调度和生态保护等工程实践提供技术支持。3、本发明方法提供了一种权衡计算效率和计算精度的水流演算新思路,可推广于相关类似工程实践中,尤其是枯水期补水调度、生态调度、水资源管理等问题。
附图说明
图1是本发明实施例中方法的技术路线图;
图2是本发明实施例中模型输入数据和参数格栅化处理示意图;
图3是本发明实施例中D8栅格编码示意图;
图4是本发明实施例中河网网格单元划分示意图;
图5是本发明实施例中初始无水状态下河道48h模拟水位分布(m);
图6是本发明实施例中初始无水状态下河道2h流速大小分布(m/s);
图7是本发明实施例中初始无水状态下河道48h流速大小分布(m/s);
图8是本发明实施例中初始水深0.02m下的模拟不同断面流量过程;
图9是本发明实施例中初始水深0.05m下的模拟不同断面流量过程;
图10是本发明实施例中点源入流且初始无水状态下河道2h流速大小分布(m/s)。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
实施例一
如图1所示,本实施例中,提供了一种基于曼宁公式的河道枯水流量计算方法,包括如下步骤,
一、数据收集及格栅化处理:
对收集或获取的河道地形数据、河床下渗数据和模型参数进行格栅化处理,以将各河道划分为多个单元格,并为各单元格进行数据输入和参数赋值。
具体包括如下内容,
1、收集河道地形数据、河床下渗数据等,并根据河道河床情况确定模型参数(糙率系数)。
2、将河道地形数据、河床下渗数据和模型参数进行格栅化处理,按一定分辨率将河道划分为多个正方形的单元格(即正方形网格,如图2所示),并对所划分的每一个单元格进行数据输入和参数赋值。
其中,单元格的最大行数由河道沿程最大宽度决定,单元格的最大列数由河道长度决定;
M=BMAX/DS (1)
N=Len/DS (2)
其中,M为河网划分单元格的最大行数;N为河网划分单元格的最大列数;BMAX为沿着水流方向最大河道宽度;Len为河道长度;DS为河道网格划分分辨率,即所划分的单元格的边长。
二、单元格水流流向计算:
基于D8算法,通过计算当前单元格与其周围单元格之间的坡降,将单元格水流方向与相邻八个单元格之间坡降最大值代表的方向作为该单元格的水流流向。
具体包括如下内容,
1、以中心单元格为计算对象,依次将其周围的单元格进行编码(如图3所示)。
2、根据上一时段计算的每个单元格水量结果,计算当前第i行、第j列的单元格(i,j)与周围单元格之间的坡降;
Slope(i,j)=dZ/DS (3)
其中,Z为单元格水位,dZ为中心单元格和相邻八个单元格的水位差;单元格(i,j)与相邻单元格之间的距离与方向有关,对于位于单元格(i,j)正上、正左、正右、正下方向上的相邻单元格(即图3中编号为2、4、6、8的相邻单元格),单元格(i,j)与相邻单元格之间的距离为DS,坡降计算公式采用方程式(3);对于位于单元格(i,j)左上、右上、左下、右下方向上的相邻单元格(即图3中编号为1、3、7、9的相邻单元格),单元格(i,j)与相邻单元格之间的距离为坡降计算公式采用方程式(4)。
3、根据单元格(i,j)坡降计算结果,假设每个计算时间步长内水流方向唯一,单元格(i,j)水流方向为与该单元格与相邻八个单元格间坡度最大值代表的方向;
Aspect(i,j)=f(max Slope(i,j)) (5)
其中,Aspect(i,j)为计算单元格与周围单元格的水流交换方向代号(取值如附图3所示);f()为最大坡度对应单元格指示的水流交换方向代号。
其中,最大坡度对应单元格指示的水流交换方向代号f()采用如下规则确定,
(i)若maxSlope(i,j)≤0,则赋以Aspect(i,j)为当前单元格代号值(即图3中5),说明当前单元格水流不流出;
(ii)若maxSlope(i,j)>0,且只有一个相邻单元格对应坡降最大,则赋以Aspect(i,j)为对应相邻单元格的编码值(附图3所示);
(iii)若maxSlope(i,j)>0,且有两个及以上相邻单元格对应坡降最大值,则按照顺时针方向(即1-2-3-6-9-8-7-4)赋以Aspect(i,j)编码号;例如编码为2、9、4的单元格都对应最大坡降,则Aspect(i,j)赋值为2。
三、单元格水流流速计算:
基于单元格的水流流向及对应坡降,利用曼宁公式计算单元格的水流流速,并对流速大的急流情况进行修正。
具体包括如下内容,
1、根据单元格水流流向及对应坡降,由曼宁公式计算单元格水流流速,
其中,v为单元格水流流速;ns为曼宁糙率系数(简称糙率),使用曼宁系数表可查到不同河床类型糙率值;R为水力半径,近似为单元水深。
2、对于上一步计算出的水流流速,计算公式基于均匀流动得到,曼宁公式对一般枯水流量有效,但对流速较大的急流情况则需要修正。
因此,根据弗洛德数判定流态的方法,对速度限定上限阀值,以对流速大的急流情况进行修正;
其中,g为重力加速度;h(i,j)为单元格(i,j)的水深。
四、单元格水量交换量和单元格水位计算:
分别计算单元格水量交换量所包含的单元格流入流出量和单元格渗漏水量,基于获取的单元格水量交换量,利用水量平衡原理,计算各单元格水位。
单元格水量交换量包括单元格流入流出量和单元格渗漏水量。
具体包括如下内容,
1、基于单元格水流流速,计算单元格流出水量;
dw(i,j)=v·Δt·DS·h(i,j) (8)
其中,Δt为计算时段长;dw(i,j)为单元格(i,j)流出水量,当流向为正上、正左、正右、正下时(即流向编号为2、4、6、8时),dw(i,j)由公式(8)计算,当流向为左上、右上、左下、右下时(即流向编号为1、3、7、9时),dw(i,j)由公式(9)计算。
2、基于霍顿产流公式,计算单元格渗漏水量;
Fp=[fc+(f0-fc)e-βt]h(i,j) (10)
其中,Ep为单元格渗漏水量;f0为干河床时最大下渗率,随着入渗过程进行,减小到近似常数的入渗率fc;参数β控制着入渗能力的下降速度;f0、fc和β根据实测入渗资料率定。
3、基于单元格流出水量和单元格渗漏水量,根据水量平衡原理,计算单元格水位;
其中,Z(t)和Z(t+Δt)分别为计算时段初和时段末的单元格水位;∑dw′(i,j)为相邻单元格流入单元格(i,j)的水量。
五、单元格水量修正:
所有单元格水位计算完毕后,利用各单元格流入扣除流出和下渗水量得到的单位面积不平衡水量对所有单元格逐个进行水量平衡修正。
具体地:计算完毕所有单元格水位后,利用如下公式,对所有单元格逐个进行水量平衡修正,
Z‘(t+Δt)=Z(t+Δt)+ΔZ (12)
其中,ΔZ为每个单元格流入扣除流出和下渗水量得到的单位面积不平衡水量;Z‘(t+Δt)为平衡修正后的水位。
六、单元格水位信息更新:
更新所有单元格水位信息,进入下一水流计算时段。
具体地:完成当前时段所有单元格的水流流速、水位计算后,将相关计算结果作为下一时段的初始值,返回步骤S2。
实施例二
本实施例中,以一个长10km,宽100m的矩形河道为例对本发明方法进行实施,如附图4所示。所选实施例河道从上游有一非恒定小流量序列输入模型(线源输入,即边界处所有单元格流速相同),考虑三种不同情形的初始水深,包括干河床、初始0.02m水深和初始0.05m水深。本实施例步骤和效益分析如下:
1、根据河道地形数据,插值生成分辨率为10m的河道地形网格单元,如附图4所示。赋值所有单元格相同糙率0.02、相同初始水深(分别0,0.02m和0.05m),相同下渗系数(暂不考虑下渗,三个系数皆取值为0)。
2、采用Fortran编程实现水流数值计算,并输出流速、水位信息,如附图5和6所示。从图中可看出,所提出方法可计算出干河床条件下的流场信息,解决了以往河道水动力模型初始为干河床不能计算的问题。
3、将本实施例计算结果与传统一维水动力学方法结果进行比较,对方法进行效益分析,如附图7和8所示。从图中可看出,所提出方法可解决以往方法在枯水条件下计算不稳定的问题。当初始水深越小时,本发明所提出方法优势越明显。此外,根据附图对比结果还可看出,本发明所提出方法在精度上并不逊于一维水动力学方法。
4、本实施例中,由于地形和边界条件都满足河宽方向的均匀性,因此,模拟的结果也呈现出一维特征。尽管如此,方法可用于分析二维空间结果,如附图9所示的在点源入流情况下的模拟流速平面分布情况。根据附图9结果,本发明所提出方法可反映出在点源入流情况下的流场空间差异性。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:
本发明提供了一种基于曼宁公式的河道枯水流量计算方法,本发明提供的河道枯水流量计算方法,解决了河道枯水流量小水深或干河床情况下的演算不稳定问题。一方面,方法保持了一维度水动力学方法的简单实用,在保证水流计算效率和计算精度的基础上,解决了以往水动力学方法在小水深、干河床情况下的不适用或不稳定问题,另一方面,方法能实现河道水流的空间二维演算,可更为准确地刻画出天然河道小水深情况下横向流速、水位的空间差异。本发明方法实现了枯水期河道水流的稳定演算,可为枯水河道补水调度和生态保护等工程实践提供技术支持。本发明方法提供了一种权衡计算效率和计算精度的水流演算新思路,可推广于相关类似工程实践中,尤其是枯水期补水调度、生态调度、水资源管理等问题。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (9)
1.一种基于曼宁公式的河道枯水流量计算方法,其特征在于:包括如下步骤,
S1、数据收集及格栅化处理:
对收集或获取的河道地形数据、河床下渗数据和模型参数进行格栅化处理,以将各河道划分为多个单元格,并为各单元格进行数据输入和参数赋值;
S2、单元格水流流向计算:
基于D8算法,通过计算当前单元格与其周围单元格之间的坡降,将单元格水流方向与相邻八个单元格之间坡降最大值代表的方向作为该单元格的水流流向;
S3、单元格水流流速计算:
基于单元格的水流流向及对应坡降,利用曼宁公式计算单元格的水流流速,并对流速大的急流情况进行修正;
S4、单元格水量交换量和单元格水位计算:
分别计算单元格水量交换量所包含的单元格流入流出量和单元格渗漏水量,基于获取的单元格水量交换量,利用水量平衡原理,计算各单元格水位;
S5、单元格水量修正:
所有单元格水位计算完毕后,利用各单元格流入扣除流出和下渗水量得到的单位面积不平衡水量对所有单元格逐个进行水量平衡修正;
S6、单元格水位信息更新:
更新所有单元格水位信息,进入下一水流计算时段。
2.根据权利要求1所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S1具体包括如下内容,
S11、收集河道地形数据和河床下渗数据,并根据河道河床情况确定模型参数;
S12、将河道地形数据、河床下渗数据和模型参数进行格栅化处理,按一定分辨率将河道划分为多个正方形的单元格,并对所划分的每一个单元格进行数据输入和参数赋值。
3.根据权利要求2所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S12中,单元格的最大行数由河道沿程最大宽度决定,单元格的最大列数由河道长度决定;
M=BMAX/DS (1)
N=Len/DS (2)
其中,M为河网划分单元格的最大行数;N为河网划分单元格的最大列数;BMAX为沿着水流方向最大河道宽度;Len为河道长度;DS为河道网格划分分辨率,即所划分的单元格的边长。
4.根据权利要求3所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S2具体包括如下内容,
S21、以中心单元格为计算对象,依次将其周围的单元格进行编码;
S22、根据上一时段计算的每个单元格水量结果,计算当前第i行、第j列的单元格(i,j)与周围单元格之间的坡降;
Slope(i,j)=dZ/DS (3)
其中,Z为单元格水位,dZ为中心单元格和相邻八个单元格的水位差;单元格(i,j)与相邻单元格之间的距离与方向有关,对于位于单元格(i,j)正上、正左、正右、正下方向上的相邻单元格,单元格(i,j)与相邻单元格之间的距离为DS,坡降计算公式采用方程式(3);对于位于单元格(i,j)左上、右上、左下、右下方向上的相邻单元格,单元格(i,j)与相邻单元格之间的距离为坡降计算公式采用方程式(4);
S23、根据单元格(i,j)坡降计算结果,假设每个计算时间步长内水流方向唯一,单元格(i,j)水流方向为与该单元格与相邻八个单元格间坡度最大值代表的方向;
Aspect(i,j)=f(maxSlope(i,j)) (5)
其中,Aspect(i,j)为计算单元格与周围单元格的水流交换方向代号;f()为最大坡度对应单元格指示的水流交换方向代号。
5.根据权利要求4所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S23中最大坡度对应单元格指示的水流交换方向代号f()采用如下规则确定,
若maxSlope(i,j)≤0,则赋以Aspect(i,j)为当前单元格代号值,说明当前单元格水流不流出;
若maxSlope(i,j)>0,且只有一个相邻单元格对应坡降最大,则赋以Aspect(i,j)为对应相邻单元格的编码值;
若maxSlope(i,j)>0,且有两个及以上相邻单元格对应坡降最大值,则按照顺时针方向赋以Aspect(i,j)编码号。
6.根据权利要求5所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S3具体包括如下内容,
S31、根据单元格水流流向及对应坡降,由曼宁公式计算单元格水流流速,
其中,v为单元格水流流速;ns为曼宁糙率系数;R为水力半径,近似为单元水深;
S32、根据弗洛德数判定流态的方法,对速度限定上限阀值,以对流速大的急流情况进行修正;
其中,g为重力加速度;h(i,j)为单元格(i,j)的水深。
7.根据权利要求6所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S4具体包括如下内容,
S41、基于单元格水流流速,计算单元格流出水量;
dw(i,j)=v·Δt·DS·h(i,j) (8)
其中,Δt为计算时段长;dw(i,j)为单元格(i,j)流出水量,当流向为正上、正左、正右、正下时,dw(i,j)由公式(8)计算,当流向为左上、右上、左下、右下时,dw(i,j)由公式(9)计算;
S42、基于霍顿产流公式,计算单元格渗漏水量;
Fp=[fc+(f0-fc)e-βt]h(i,j) (10)
其中,Fp为单元格渗漏水量;f0为干河床时最大下渗率,随着入渗过程进行,减小到近似常数的入渗率fc;参数β控制着入渗能力的下降速度;f0、fc和β根据实测入渗资料率定;
S43、基于单元格流出水量和单元格渗漏水量,根据水量平衡原理,计算单元格水位;
其中,Z(t)和Z(t+Δt)分别为计算时段初和时段末的单元格水位;∑dw′(i,j)为相邻单元格流入单元格(i,j)的水量。
8.根据权利要求7所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S5具体为,计算完毕所有单元格水位后,利用如下公式,对所有单元格逐个进行水量平衡修正,
Z‘(t+Δt)=Z(t+Δt)+ΔZ (12)
其中,ΔZ为每个单元格流入扣除流出和下渗水量得到的单位面积不平衡水量;Z‘(t+Δt)为平衡修正后的水位。
9.根据权利要求1所述的基于曼宁公式的河道枯水流量计算方法,其特征在于:步骤S6具体为,完成当前时段所有单元格的水流流速、水位计算后,将相关计算结果作为下一时段的初始值,返回步骤S2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310893339.8A CN116955889A (zh) | 2023-07-20 | 2023-07-20 | 一种基于曼宁公式的河道枯水流量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310893339.8A CN116955889A (zh) | 2023-07-20 | 2023-07-20 | 一种基于曼宁公式的河道枯水流量计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116955889A true CN116955889A (zh) | 2023-10-27 |
Family
ID=88452351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310893339.8A Pending CN116955889A (zh) | 2023-07-20 | 2023-07-20 | 一种基于曼宁公式的河道枯水流量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116955889A (zh) |
-
2023
- 2023-07-20 CN CN202310893339.8A patent/CN116955889A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kirkby | A two-dimensional simulation model for slope and stream evolution | |
Wolock et al. | Comparison of single and multiple flow direction algorithms for computing topographic parameters in TOPMODEL | |
CN110457651B (zh) | 一种利用数据同化改进的滑坡位移预测方法 | |
CN111080157B (zh) | 一种梯级水电站排磷量的调度方法和系统 | |
CN112464584B (zh) | 自由表面流的水位和流量推求方法 | |
CN104268645B (zh) | 水资源配置中考虑供水优先序的行业供水量确定方法 | |
CN110598290A (zh) | 考虑气候变化的流域未来水电发电能力预测方法和系统 | |
CN109871621A (zh) | 城市暴雨内涝汇水区分析方法 | |
CN103093114A (zh) | 一种基于地形和土壤特性的分布式流域缺水量测算方法 | |
CN106320256B (zh) | 一种多沙河流抽水蓄能电站过机泥沙通量的确定方法 | |
CN104156570A (zh) | 一种土壤侵蚀类面源污染物入河量的计算方法及装置 | |
CN108763798B (zh) | 一种湖泊与地下水非稳定流作用模拟方法 | |
CN113486556B (zh) | 一种改进的油气藏高效自动历史拟合方法 | |
CN103870699B (zh) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 | |
CN116955889A (zh) | 一种基于曼宁公式的河道枯水流量计算方法 | |
CN112231913A (zh) | 一种城市内涝黑点的内涝模拟方法及装置 | |
CN113449982A (zh) | 基于控制性生态因子比尺的湖泊生态水文节律确定方法 | |
CN110837684B (zh) | 基于通用计算的流域水质实时模拟方法 | |
CN115983158B (zh) | 一种地下水模型与二维水动力模型松散耦合的方法 | |
CN114742299B (zh) | 基于径流侵蚀功率的流域泥沙预报模型构建方法 | |
CN112084643A (zh) | 基于数字高程和土壤参数的流域提取方法 | |
CN115564238B (zh) | 一种基于水动力机制的小水库洪水预报调洪分析方法 | |
CN113435631B (zh) | 一种洪水预报方法、系统、可读存储介质及计算设备 | |
CN110578317A (zh) | 一种水文模型水库泄水量模拟方法 | |
CN114428989A (zh) | 一种游荡分汊河道洪峰流量分配方法及系统 |
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 |