CN110110455A - 给定流量下调节闸门开度的方法 - Google Patents
给定流量下调节闸门开度的方法 Download PDFInfo
- Publication number
- CN110110455A CN110110455A CN201910395474.3A CN201910395474A CN110110455A CN 110110455 A CN110110455 A CN 110110455A CN 201910395474 A CN201910395474 A CN 201910395474A CN 110110455 A CN110110455 A CN 110110455A
- Authority
- CN
- China
- Prior art keywords
- water
- place value
- flow
- channel
- gate
- 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
-
- 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);n为地表糙率(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 (4)
1.给定流量下调节闸门开度的方法,其特征在于,包括:
S1、根据闸门的孔的数量以及给定流量Q,利用Saint-Venant方程组,获取渠道的横断面尺寸;
S2、对渠道俯视平面用非结构化网格进行空间离散;
S3、在空间离散后的非结构化网格网格中,以实线交点定义水位信息,虚线交点定义流量信息,且闸门位于虚线上;
S4、根据步骤S3中的流量和水位的定义,采用全隐式标量耗散有限体积法求解Saint-Venant方程组,并根据Saint-Venant方程组的解建立过闸水运动实时并行计算模型;
S5、以给定流量Q为所述过闸水运动实时并行计算模型的内边界条件,并以小于等于闸门最高点处水位值为上游水位边界条件,运算过闸水运动实时并行计算模型,获得过闸水流运动过程唯一恒定流态下的模拟上游水位值和下游水位值;
S6、调节渠道的闸门开度,直至实测的上游水位值和下游水位值与唯一恒定流态下的模拟上游水位值和下游水位值吻合,此时的闸门开度即给定流量Q下的闸门开度。
2.根据权利要求1所述的方法,其特征在于,步骤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)。
3.根据权利要求1或2所述的方法,其特征在于,步骤S4中,以闸门最高点处水位值为上游水位边界条件,若过闸水流运动过程无法达到恒定流态,则逐渐减小上游水位边界条件后运算过闸水运动实时并行计算模型,直至获得唯一恒定流态下的模拟上游水位值和下游水位值。
4.根据权利要求3所述的方法,其特征在于,步骤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 true CN110110455A (zh) | 2019-08-09 |
CN110110455B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115423346A (zh) * | 2022-09-20 | 2022-12-02 | 中国水利水电科学研究院 | 一种明渠调水工程串联闸群汛期可调节能力实时测算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04195082A (ja) * | 1990-11-28 | 1992-07-15 | Mitsubishi Electric Corp | 配管系統の流量圧力模擬方式 |
CN101403929A (zh) * | 2008-11-06 | 2009-04-08 | 清华大学 | 一种自动控制闸门泄流量的方法和系统 |
CN103017834A (zh) * | 2012-12-29 | 2013-04-03 | 大连理工大学 | 一种利用水电站泄水弧形闸门开度准确测定泄流量的方法 |
CN105780717A (zh) * | 2016-03-24 | 2016-07-20 | 华北水利水电大学 | 一种多沙渠道引水闸门的水工模型 |
CN106874622A (zh) * | 2017-03-13 | 2017-06-20 | 中国水利水电科学研究院 | 一种串联多渠池闸门过闸流量系数率定方法 |
-
2019
- 2019-05-13 CN CN201910395474.3A patent/CN110110455B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04195082A (ja) * | 1990-11-28 | 1992-07-15 | Mitsubishi Electric Corp | 配管系統の流量圧力模擬方式 |
CN101403929A (zh) * | 2008-11-06 | 2009-04-08 | 清华大学 | 一种自动控制闸门泄流量的方法和系统 |
CN103017834A (zh) * | 2012-12-29 | 2013-04-03 | 大连理工大学 | 一种利用水电站泄水弧形闸门开度准确测定泄流量的方法 |
CN105780717A (zh) * | 2016-03-24 | 2016-07-20 | 华北水利水电大学 | 一种多沙渠道引水闸门的水工模型 |
CN106874622A (zh) * | 2017-03-13 | 2017-06-20 | 中国水利水电科学研究院 | 一种串联多渠池闸门过闸流量系数率定方法 |
Non-Patent Citations (1)
Title |
---|
章晋雄: "流量对事故闸门门体水力荷载影响数值模拟研究", 《水利水电技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115423346A (zh) * | 2022-09-20 | 2022-12-02 | 中国水利水电科学研究院 | 一种明渠调水工程串联闸群汛期可调节能力实时测算方法 |
CN115423346B (zh) * | 2022-09-20 | 2023-04-07 | 中国水利水电科学研究院 | 一种明渠调水工程串联闸群汛期可调节能力实时测算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110110455B (zh) | 2020-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109657418B (zh) | 一种基于mike21的湖泊水环境容量计算方法 | |
WO2021129039A1 (zh) | 一种跨水库流域河道径流演算方法和系统 | |
CN104933268B (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
WO2017201812A1 (zh) | 一种渠道流量的获取方法 | |
CN109271672B (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN107045568B (zh) | 基于动态规划逐次逼近法的河道糙率反演方法 | |
CN110705171A (zh) | 一种基于mike模型的感潮河网水环境治理方法 | |
CN112257352A (zh) | 一维水动力模型和二维水动力模型的耦合方法及系统 | |
CN112464584B (zh) | 自由表面流的水位和流量推求方法 | |
CN112227292A (zh) | 基于一维水动力模型的空腹拱桥壅水智能分析方法 | |
CN112163381B (zh) | 一种适用于复杂地形风场流动数值模拟的侧向边界条件设置方法 | |
CN113626923A (zh) | 一种平原河网区域水资源调度精准率计算与提升方法 | |
CN107256312B (zh) | 一种基于潮流环境下海湾纳潮变化量计算方法 | |
CN111400974B (zh) | 一种估算矩形河渠壁面和床面切应力的方法 | |
CN110990926B (zh) | 一种基于面积修正率的城市地表建筑水动力学仿真方法 | |
CN110147646B (zh) | 一种数值模拟框架下线性挡水构筑物的过流处理方法 | |
CN108897940A (zh) | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 | |
CN110110455A (zh) | 给定流量下调节闸门开度的方法 | |
CN107451399A (zh) | 一种基于实域总势能计算渗流溢出点的方法 | |
CN111311734A (zh) | 铀尾矿库排洪系统水力特性三维数值模拟方法 | |
CN108729405A (zh) | 一种河工模型试验方法 | |
CN110847112B (zh) | 一种基于水力学模拟的河道泄洪预警方法 | |
CN116562191A (zh) | 一种融合遥感影像的湖泊水动力模型水位参数校正方法 | |
CN114896909A (zh) | 一种基于水位高度的明渠流量计算方法 | |
Xie et al. | Mathematical model for flood routing in Jingjiang River and Dongting Lake network |
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 |