CN110110455B - 给定流量下调节闸门开度的方法 - Google Patents
给定流量下调节闸门开度的方法 Download PDFInfo
- Publication number
- CN110110455B CN110110455B CN201910395474.3A CN201910395474A CN110110455B CN 110110455 B CN110110455 B CN 110110455B CN 201910395474 A CN201910395474 A CN 201910395474A CN 110110455 B CN110110455 B CN 110110455B
- Authority
- CN
- China
- Prior art keywords
- water level
- gate
- level value
- flow
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种给定流量下调节闸门开度的方法。其包括根据闸门的孔的数量以及给定流量Q,利用Saint‑Venant方程组,获取渠道的横断面尺寸;对渠道俯视平面用非结构化网格进行空间离散;并以实线交点定义水位信息,虚线交点定义流量信息,且闸门位于虚线上;采用全隐式标量耗散有限体积法求解Saint‑Venant方程组,并根据Saint‑Venant方程组的解建立过闸水运动实时并行计算模型;以给定流量Q为过闸水运动实时并行计算模型的内边界条件,并以小于等于闸门最高点处水位值为上游水位边界条件,运算过闸水运动实时并行计算模型,获得恒定流态下的模拟上游水位值和下游水位值;调节渠道的闸门开度,直至实测的上游水位值和下游水位值与唯一恒定流态下的模拟上游水位值和下游水位值吻合。
Description
技术领域
本发明涉及水利工程领域,具体涉及一种给定流量下调节闸门开度的方法。
背景技术
给定渠道控制节点处的流量后,如何设定闸门开度,在流域或灌区地表水网调控过程中极其重要。传统的确定闸门开度的方法是比例分配法,即基于最大开度对应的过流量和无开度时零过流量,开等比例的分配流量与开度。该方法简洁易用,但完全未考虑水运动过程的势能与动能梯度驱动的非线性特征,故精度不高。
发明内容
针对现有技术中的上述不足,本发明旨在提供一种更加精准的给定流量下调节闸门开度的方法。
为了达到上述发明创造的目的,本发明采用的技术方案为:
提供一种给定流量下调节闸门开度的方法,其包括:
S1、根据闸门的孔的数量以及给定流量Q,利用Saint-Venant方程组,获取渠道的横断面尺寸;
S2、对渠道俯视平面用非结构化网格进行空间离散;
S3、在空间离散后的非结构化网格中,以实线交点定义水位信息,虚线交点定义流量信息,且闸门位于虚线上;
S4、根据步骤S3中的流量和水位的定义,采用全隐式标量耗散有限体积法求解Saint-Venant方程组,并根据Saint-Venant方程组的解建立过闸水运动实时并行计算模型;
S5、以给定流量Q为过闸水运动实时并行计算模型的内边界条件,并以小于等于闸门最高点处水位值为上游水位边界条件,运算过闸水运动实时并行计算模型,获得过闸水流运动过程唯一恒定流态下的模拟上游水位值和下游水位值;
S6、调节渠道的闸门开度,直至实测的上游水位值和下游水位值与唯一恒定流态下的模拟上游水位值和下游水位值吻合,此时的闸门开度即给定流量Q下的闸门开度。
进一步地,当闸门的孔的数量为1时,Saint-Venant方程组为:
其中:地表水位ζ=h+zb,h为地表水深(m)、zb是地表相对高程(m);t为时间坐标(s);s为沿渠道方向的距离坐标(m);Q为给定流量(m3/s·m);为渠道的过流速度(m/s);A为渠道的过流断面面积(m2);R为渠道的水力半径(m);g=9.8为重力加速度(m/s2);n为地表糙率(s/m1/3);
当闸门的孔的数量大于1时,Saint-Venant方程组为:
其中,x和y为以闸门方向为X轴,渠道方向为Y轴的坐标系中的坐标(m);u和v分别为水流沿x和y坐标向的垂向均布流速(m/s);qx和qy分别为水流沿x和y坐标向的单宽流量(m2/s)。
进一步地,步骤S4中,以闸门最高点处水位值为上游水位边界条件,若过闸水流运动过程无法达到恒定流态,则逐渐减小上游水位边界条件后运算过闸水运动实时并行计算模型,直至获得唯一恒定流态下的模拟上游水位值和下游水位值。
进一步地,当实测的上游水位值和下游水位值与唯一恒定流态下的上游水位值和下游水位值的误差值小于预设误差值时,即判定为吻合。
本发明的有益效果为:
本发明以上游水位和下游水位为调控目标,通过Saint-Venant方程组完全考虑水运动过程的势能与动能梯度驱动机制的非线性特征来建立流量、水位和闸门开度这三者之间的定量映射关系,进而更加精准地确定给定流量下的闸门开度。并且本发明无需闸门具体的几何参数,即能确定单孔或多孔闸门的开度。
附图说明
图1为本方法的局部流程图;
图2为闸门的孔的数量大于1时,空间离散后的非结构化网格中水位、流量及闸门定位位置示意图;
图3为河北冶河灌区的鹿泉大李庄试验点的闸门和渠道底部几何尺寸图。
具体实施方式
下面结合附图,对本发明的具体实施方式做详细说明,以便于本技术领域的技术人员理解本发明。但应该清楚,下文所描述的实施例仅仅是本发明的一部分实施例,而不是全部实施例。在不脱离所附的权利要求限定和确定的本发明的精神和范围内,本领域普通技术人员在没有做出任何创造性劳动所获得的所有其他实施例,都属于本发明的保护范围。
如图1所示,该给定流量下调节闸门开度的方法,包括:
S1、根据闸门的孔的数量以及给定流量Q,利用Saint-Venant方程组,获取渠道的横断面尺寸;
S2、对渠道俯视平面用非结构化网格进行空间离散;
S3、在空间离散后的非结构化网格中,以实线交点定义水位信息,虚线交点定义流量信息,且闸门位于虚线上;
S4、根据步骤S3中的流量和水位的定义,采用全隐式标量耗散有限体积法求解Saint-Venant方程组,并根据Saint-Venant方程组的解建立过闸水运动实时并行计算模型;
S5、以给定流量Q为所述过闸水运动实时并行计算模型的内边界条件,并以小于等于闸门最高点处水位值为上游水位边界条件,运算过闸水运动实时并行计算模型,获得过闸水流运动过程唯一恒定流态下的模拟上游水位值和下游水位值;
S6、调节渠道的闸门开度,直至实测的上游水位值和下游水位值与唯一恒定流态下的模拟上游水位值和下游水位值吻合,此时的闸门开度即给定流量Q下的闸门开度。
实施时,本方案优选步骤S1中,当闸门的孔的数量为1时,Saint-Venant方程组为:
其中:地表水位ζ=h+zb,h为地表水深(m)、zb是地表相对高程(m);t为时间坐标(s);s为沿渠道方向的距离坐标(m);Q为给定流量(m3/s·m);为渠道的过流速度(m/s);A为渠道的过流断面面积(m2);R为渠道的水力半径(m);g=9.8为重力加速度(m/s2);为地表糙率(s/m1/3);
当闸门的孔的数量大于1时,Saint-Venant方程组为:
其中,x和y为以闸门方向为X轴,渠道方向为Y轴的坐标系中的坐标(m);u和v分别为水流沿x和y坐标向的垂向均布流速(m/s);qx和qy分别为水流沿x和y坐标向的单宽流量(m2/s),此时Q=hu。
并且,在步骤S4中,以闸门最高点处水位值为上游水位边界条件,若过闸水流运动过程无法达到恒定流态,则逐渐减小上游水位边界条件后运算过闸水运动实时并行计算模型,直至获得唯一恒定流态下的模拟上游水位值和下游水位值。
在步骤S5中,当实测的上游水位值和下游水位值与唯一恒定流态下的上游水位值和下游水位值的误差值小于预设误差值时,即判定为吻合。
下面以河北冶河灌区的鹿泉大李庄试验点为例,应用上述方法。
该试验点的闸门和渠底几何尺寸见图2,其中,渠道两侧的坡度为10°。基于这些渠道几何参数,在给定流量Q分别为16.85(m3/s)、17.36(m3/s)和18.21(m3/s)的条件下,上下游水位的实测值和唯一恒定流态下的模拟值以及流量的给定值和实测值与闸门开度的对应表如下:
以上述河北冶河灌区的鹿泉大李庄试验点为例,采用传统的比例分配法,在给定流量Q分别为16.85(m3/s)、17.36(m3/s)和18.21(m3/s)的条件下,确定的闸门开度及流量对应表如下:
结果表明当将上下游水位实测值与上下游水位的模拟值吻合时,流量的实测值与给定值之间具有极好的拟合度,由此确定了更加精准的闸门开度。
Claims (3)
1.给定流量下调节闸门开度的方法,其特征在于,包括:
S1、根据闸门的孔的数量以及给定流量Q,利用Saint-Venant方程组,获取渠道的横断面尺寸;
当闸门的孔的数量为1时,Saint-Venant方程组为:
其中:地表水位ζ=h+zb,h为地表水深(m)、zb是地表相对高程(m);t为时间坐标(s);s为沿渠道方向的距离坐标(m);Q为给定流量(m3/s·m);为渠道的过流速度(m/s);A为渠道的过流断面面积(m2);R为渠道的水力半径(m);g=9.8为重力加速度(m/s2);n为地表糙率(s/m1/3);
当闸门的孔的数量大于1时,Saint-Venant方程组为:
其中,x和y为以闸门方向为X轴,渠道方向为Y轴的坐标系中的坐标(m);u和v分别为水流沿x和y坐标向的垂向均布流速(m/s);qx和qy分别为水流沿x和y坐标向的单宽流量(m2/s);
S2、对渠道俯视平面用非结构化网格进行空间离散;
S3、在空间离散后的非结构化网格中,以实线交点定义水位信息,虚线交点定义流量信息,且闸门位于虚线上;
S4、根据步骤S3中的流量和水位的定义,采用全隐式标量耗散有限体积法求解Saint-Venant方程组,并根据Saint-Venant方程组的解建立过闸水运动实时并行计算模型;
S5、以给定流量Q为所述过闸水运动实时并行计算模型的内边界条件,并以小于等于闸门最高点处水位值为上游水位边界条件,运算过闸水运动实时并行计算模型,获得过闸水流运动过程唯一恒定流态下的模拟上游水位值和下游水位值;
S6、调节渠道的闸门开度,直至实测的上游水位值和下游水位值与唯一恒定流态下的模拟上游水位值和下游水位值吻合,此时的闸门开度即给定流量Q下的闸门开度。
2.根据权利要求1所述的方法,其特征在于,步骤S4中,以闸门最高点处水位值为上游水位边界条件,若过闸水流运动过程无法达到恒定流态,则逐渐减小上游水位边界条件后运算过闸水运动实时并行计算模型,直至获得唯一恒定流态下的模拟上游水位值和下游水位值。
3.根据权利要求2所述的方法,其特征在于,步骤S5中,当实测的上游水位值和下游水位值与唯一恒定流态下的上游水位值和下游水位值的误差值小于预设误差值时,即判定为吻合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910395474.3A CN110110455B (zh) | 2019-05-13 | 2019-05-13 | 给定流量下调节闸门开度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910395474.3A CN110110455B (zh) | 2019-05-13 | 2019-05-13 | 给定流量下调节闸门开度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110110455A CN110110455A (zh) | 2019-08-09 |
CN110110455B true CN110110455B (zh) | 2020-01-07 |
Family
ID=67489773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910395474.3A Active CN110110455B (zh) | 2019-05-13 | 2019-05-13 | 给定流量下调节闸门开度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110110455B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115423346B (zh) * | 2022-09-20 | 2023-04-07 | 中国水利水电科学研究院 | 一种明渠调水工程串联闸群汛期可调节能力实时测算方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04195082A (ja) * | 1990-11-28 | 1992-07-15 | Mitsubishi Electric Corp | 配管系統の流量圧力模擬方式 |
CN101403929B (zh) * | 2008-11-06 | 2010-07-07 | 清华大学 | 一种自动控制闸门泄流量的方法和系统 |
CN103017834A (zh) * | 2012-12-29 | 2013-04-03 | 大连理工大学 | 一种利用水电站泄水弧形闸门开度准确测定泄流量的方法 |
CN105780717B (zh) * | 2016-03-24 | 2018-03-06 | 华北水利水电大学 | 一种用于多沙渠道的闸门量水方法 |
CN106874622B (zh) * | 2017-03-13 | 2018-05-11 | 中国水利水电科学研究院 | 一种串联多渠池闸门过闸流量系数率定方法 |
-
2019
- 2019-05-13 CN CN201910395474.3A patent/CN110110455B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110110455A (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110287571B (zh) | 一种河流险工冲刷安全分析与岸坡稳定性判定方法 | |
WO2021129039A1 (zh) | 一种跨水库流域河道径流演算方法和系统 | |
CN103970989B (zh) | 一种基于流固界面一致条件的浸入边界流场计算方法 | |
WO2017201812A1 (zh) | 一种渠道流量的获取方法 | |
CN108021780B (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN109271672B (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN106498896B (zh) | 一种与泥沙冲淤相耦合的水库调度过程模拟方法 | |
CN104933268B (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
Oskuyi et al. | Vertical sluice gate discharge coefficient | |
CN110135056B (zh) | 一种河道内桥墩有效阻水比的快速自动分析方法 | |
CN112464584B (zh) | 自由表面流的水位和流量推求方法 | |
CN110110455B (zh) | 给定流量下调节闸门开度的方法 | |
CN111400974B (zh) | 一种估算矩形河渠壁面和床面切应力的方法 | |
CN107256312B (zh) | 一种基于潮流环境下海湾纳潮变化量计算方法 | |
CN110147646B (zh) | 一种数值模拟框架下线性挡水构筑物的过流处理方法 | |
WO2017201813A1 (zh) | 地面灌溉地表水流运动过程模拟方法 | |
CN105205215B (zh) | 一种计算不同状态挺水植被水流阻力的方法 | |
Bleninger et al. | Near-and far-field model coupling methodology for wastewater discharges | |
Zheng et al. | Parameter identification for discharge formulas of radial gates based on measured data | |
CN106651054A (zh) | 基于遗传算法的长距离调水工程输水渠道糙率辨识方法 | |
CN105956376A (zh) | 一种深埋隧洞施工开挖数值模拟分析的模型范围取值方法 | |
Fang et al. | Diagonal Cartesian method for the numerical simulation of flow and suspended sediment transport over complex boundaries | |
CN114896909A (zh) | 一种基于水位高度的明渠流量计算方法 | |
Hwang et al. | Hydrodynamic and water quality modeling for gate operation: A case study for the Seonakdong River basin in Korea | |
CN104239730B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |