CN110046469B - 多约束条件下水电站坝前河床冲淤变形的计算方法 - Google Patents
多约束条件下水电站坝前河床冲淤变形的计算方法 Download PDFInfo
- Publication number
- CN110046469B CN110046469B CN201910392663.5A CN201910392663A CN110046469B CN 110046469 B CN110046469 B CN 110046469B CN 201910392663 A CN201910392663 A CN 201910392663A CN 110046469 B CN110046469 B CN 110046469B
- Authority
- CN
- China
- Prior art keywords
- water
- flow
- riverbed
- water intake
- sand
- 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
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B1/00—Equipment or apparatus for, or methods of, general hydraulic engineering, e.g. protection of constructions against ice-strains
-
- 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
Abstract
本发明公开了多约束条件下水电站坝前河床冲淤变形的计算方法,属于水工结构和河流模拟技术领域,该方法通过针对水电站坝前的多控制边界条件水沙数学模型问题,提出了分步迭代的求解方法。第一步,不考虑发电站取水口等约束条件的影响,根据水电站的调度方式给定出口断面的水位控制条件,进口边界给定流量过程;第二步,根据第一步的条件计算得到水位作为模型上游边界条件,水电站取水口的流量作为下游边界条件;计算完成后对坝前的平均流速和水流挟沙力进行修正,然后代入河床变形的计算式,即可得到多约束条件下水电站坝前河床冲淤变形的计算量。本发明能准确模拟复杂条件下的流场,修正后的结果与实际情况更加符合。
Description
技术领域
本发明属于水工结构和河流模拟技术领域,具体涉及多约束条件下水电站坝前河床冲淤变形的计算方法。
背景技术
水利枢纽工程坝前水工建筑物众多,主要包括:发电引水建筑物进水口、岸边溢洪道、生态放水孔、泄洪排沙洞、放空洞等。发电引水建筑物主要功能为发电供水,按照分层取水要求,通过机组发电后进入河道向下游灌溉、生态供水。岸边溢洪道:作为主要泄水建筑物,具有超泄能力,承担泄放校核洪水、设计洪水和常遇洪水。生态放水孔:初期蓄水时向下游泄放生态流量,初期运行时具备紧急降低和放空水库的功能。永久运行时,与电站、放空洞联合承担下游灌溉及生态放水的任务。泄洪排沙洞:作为次要泄水建筑物,辅助溢洪道泄放校核洪水、设计洪水和常遇洪水。放空洞:在紧急状况下降低水位、放空水库,以便检修大坝等建筑物。与电站及生态放水孔联合承担下游灌溉及生态放水的功能。
水利枢纽坝前流态关系到水工建筑物尺寸的设计和布置点的选取,因此,水流枢纽设计过程中需要对坝前流态有明确的认识。实际工程中常采用水工物理模型试验进行研究,在模型试验中开展不同条件的水流试验,得到坝前的流态。该方法较为准确,但研究周期较长。
随着计算机的发展,许多学者采用数值模拟的方法来模拟水工建筑物前的流态。常规的数值模型只有两个控制条件,上游边界采用水位或流量,下游边界也采用水位或流量,但水电站取水口前除了上游的流量条件、下游的水位条件外,在下游边界还增加了孔口出流的条件,增加了约束控制,按照常规的数学模型难以充分利用已有的控制条件,模拟坝前的流态和冲淤变化。
发明内容
发明目的:本发明的目的是提供多约束条件下水电站坝前河床冲淤变形的计算方法,该方法通过针对水电站坝前的多控制边界条件水沙数学模型问题,提出了分步迭代的求解方法。
多约束条件下水电站坝前河床冲淤变形的计算方法,包括如下步骤:
步骤一,确定包含水工建筑物的研究区域,对研究区域进行网格剖分,采用实测地形高程数据进行插值,得到网格节点上的高程;
步骤二,分类列出步骤一确定的研究区域中所有的边界条件,进行水文资料分析,选取典型的计算条件;
步骤三,建立包含研究区域的平面二维数学模型并进行验证,率定得到研究区域的糙率;
步骤四,不考虑发电站取水口的影响,采用步骤二中的边界条件,模型进口断面给定流量边界条件,出口断面采用开边界条件,即流速梯度为0,并根据水电站的调度方式给出出口断面的水位控制条件;
步骤五,采用步骤三建立的平面二维数学模型,按照步骤四中的边界条件,计算得到上游边界的水位;
步骤六,将步骤五中计算得到的水位作为模型的上游进口边界条件,水电站取水口的流量作为下游边界条件;确定取水口中心点至水面的距离H1、中心点至河床的距离H2以及取水口的过水高度H3,通过取水口的流量和相关水深,可以计算得到取水口的初始流速;
步骤七,采用步骤六中的边界条件进行数值模拟计算,记录计算得到的研究区域网格节点上的平均流速U和平均含沙量S;
步骤八,对取水口前10H3距离内的平均流速进行修正,修正值U0=λ1U,λ1为流速修正系数;
步骤九,对取水口前10H3距离内的含沙量进行修正,修正值S0=λ2S,λ2为含沙量修正系数;
步骤十,根据河床变形方程,任意一时段坝前网格节点上的冲淤变形数值为:其中,α为恢复饱和系数,当河床淤积时取0.25,当河床冲刷时取1.0;ρ′为淤积泥沙容重;S*为水流挟沙力;S0为平均含沙量;将水流挟沙力和平均含沙量的表达式代入上式,即可得到取水口附近的河床变形值:
有益效果:与现有技术相比,本发明的多约束条件下水电站坝前河床冲淤变形的计算方法,克服了传统数学模型不能充分考虑边界条件的不足,能准确模拟复杂条件下的流场。在取水口附近,该方法修正了传统二维数学模型计算得到的平均含沙量和平均流速,修正后的结果与实际情况更加符合。
附图说明
图1为多约束条件下水电站坝前河床冲淤变形的计算方法流程图;
图2为研究区域网格划分示意图;
图3为水电站研究区域示意图;
图4为一步计算后研究区域水位分布图;
图5为第二步计算后水工建筑物前流态;
图6为水库运行50年后计算得到的坝前局部三维地形图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步的说明。
多约束条件下水电站坝前河床冲淤变形的计算方法,包括如下步骤:
步骤一,确定包含水工建筑物的研究区域,对研究区域进行网格剖分,采用实测地形高程数据进行插值,得到网格节点上的高程;
步骤二,分类列出步骤一确定的研究区域中所有的边界条件,进行水文资料分析,选取典型的计算条件;
步骤三,建立包含研究区域的平面二维数学模型并进行验证,率定得到研究区域的糙率;
步骤四,不考虑发电站取水口的影响,采用步骤二中的边界条件,模型进口断面给定流量边界条件,出口断面采用开边界条件,即流速梯度为0,并根据水电站的调度方式给出出口断面的水位控制条件;
步骤五,采用步骤三建立的平面二维数学模型,按照步骤四中的边界条件,计算得到上游边界的水位;
步骤六,将步骤五中计算得到的水位作为模型的上游进口边界条件,水电站取水口的流量作为下游边界条件;确定取水口中心点至水面的距离H1、中心点至河床的距离H2以及取水口的过水高度H3,通过取水口的流量和相关水深,可以计算得到取水口的初始流速;
步骤七,采用步骤六中的边界条件进行数值模拟计算,记录计算得到的研究区域网格节点上的平均流速U和平均含沙量S;
步骤八,对取水口前10H3距离内的平均流速进行修正,修正值U0=λ1U,λ1为流速修正系数;
步骤九,对取水口前10H3距离内的含沙量进行修正,修正值S0=λ2S,λ2为含沙量修正系数;
步骤十,根据河床变形方程,任意一时段坝前网格节点上的冲淤变形数值为:其中,α为恢复饱和系数,当河床淤积时取0.25,当河床冲刷时取1.0;ρ′为淤积泥沙容重;S*为水流挟沙力;S0为平均含沙量;将水流挟沙力和平均含沙量的表达式代入上式,即可得到取水口附近的河床变形值:
步骤六中,取水口处的初始流速为其中Q为取水口的流量,单位m3/s,V为取水口处的平均流速,π为圆周率,取3.14。
以新疆大石峡水利枢纽工程为例进行计算。
(1)计算范围的选取
二维数模计算范围:根据设计单位提供的坝区地形资料,选取建模范围为坝前1.1km至大坝坝址。上游断面为河道进口,下游断面布置水工建筑物,主要包括混凝土面板砂砾石坝、溢洪道、泄洪洞、电站机组、放空洞。
(2)网格划分和地形插值
为详细反映河道形态、枢纽布置,采用贴体网格划分计算区域,对溢洪道、泄洪洞、电站机组、放空洞等建筑物区域局部加密,布置8482个网格,计算网格如图1所示。图2为模拟河段插值后的三维地形。
(3)确定计算条件
上游为条件为来水来沙量,下游有水位条件和溢洪道、泄洪洞、电站机组、放空洞的流量控制条件。根据水文资料分析和水电站的运用方式确定100年一遇的水文条件,如表1所示。
表1泄水建筑物的泄洪运行组合表
(4)糙率的率定
采用大石峡水库建库前天然实测水文资料对二维数学模型进行了验证。表2为天然条件下实测水位与计算水位对比。由表可知,模型计算水位与实测较为接近,最大差值为7cm,符合相关要求,说明模型选取的糙率合理。
表2天然条件下实测与计算水位对比
(5)第一步计算得到进口断面的水位
进口断面采用流量条件2931m3/s,出口断面采用1700m水位控制,模拟计算得到进口断面的水位1700.1m,图2为研究区域水位分布图。
(6)更换模型边界条件,重新计算研究区域的流场。
进口断面采用水位边界条件1700.1m,下游采用流量边界条件控制,其中溢洪道下泄流量1609.8m3/s,泄洪洞下泄流量1025m3/s,电站机组下泄流量296.2m3/s,总共下泄流量2931m3/s。图3为计算得到的水工建筑物前流态,水流进入溢洪道的水流主要有两股,一是大坝前沿水流;二是溢洪道上游水流。受泄洪洞的影响,水电站前的水流均偏向泄洪洞。泄洪洞前沿流速最大;由于坝前控制水位较高,坝区整体水深较大,坝区内的流速均较小。该计算结果与实际相符。
(8)按照公式计算某一时段坝前河床的冲淤变化,并多个时段累加,得到多约束条件下水电站坝前河床冲淤变形结果。图4为一步计算后研究区域水位分布图;图5为第二步计算后水工建筑物前流态;图6为水库运行50年后计算得到的坝前局部三维地形图。
本发明的描述和应用实例是说明性的,发明使用范围并非限制在上述实施例中。本发明专利使用人员及本领域技术人员须知晓并执行的是,如对本发明技术方案的细节或形式进行的合理修改或替换,均落入本发明的专利保护范围内。
Claims (5)
1.多约束条件下水电站坝前河床冲淤变形的计算方法,其特征在于:包括如下步骤:
步骤一,确定包含水工建筑物的研究区域,对研究区域进行网格剖分,采用实测地形高程数据进行插值,得到网格节点上的高程;
步骤二,分类列出步骤一确定的研究区域中所有的边界条件,进行水文资料分析,选取典型的计算条件;
步骤三,建立包含研究区域的平面二维数学模型并进行验证,率定得到研究区域的糙率;
步骤四,不考虑发电站取水口的影响,采用步骤二中的边界条件,模型进口断面给定流量边界条件,出口断面采用开边界条件,流速梯度为0,并根据水电站的调度方式给出出口断面的水位控制条件;
步骤五,采用步骤三建立的平面二维数学模型,按照步骤四中的边界条件,计算得到上游边界的水位;
步骤六,将步骤五中计算得到的水位作为模型的上游进口边界条件,水电站取水口的流量作为下游边界条件;确定取水口中心点至水面的距离H1、中心点至河床的距离H2以及取水口的过水高度H3,通过取水口的流量和相关水深,可以计算得到取水口的初始流速;
步骤七,采用步骤六中的边界条件进行数值模拟计算,记录计算得到的研究区域网格节点上的平均流速U和平均含沙量S;
步骤八,对取水口前10H3距离内的平均流速进行修正,修正值U0=λ1U,λ1为流速修正系数;
步骤九,对取水口前10H3距离内的含沙量进行修正,修正值S0=λ2S,λ2为含沙量修正系数;
步骤十,根据河床变形方程,任意一时段坝前网格节点上的冲淤变形数值为:其中,α为恢复饱和系数,当河床淤积时取0.25,当河床冲刷时取1.0;ρ′为淤积泥沙容重;S*为水流挟沙力;ω为沉速;S0为平均含沙量;将水流挟沙力和平均含沙量的表达式代入上式,即可得到取水口附近的河床变形值:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910392663.5A CN110046469B (zh) | 2019-05-13 | 2019-05-13 | 多约束条件下水电站坝前河床冲淤变形的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910392663.5A CN110046469B (zh) | 2019-05-13 | 2019-05-13 | 多约束条件下水电站坝前河床冲淤变形的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110046469A CN110046469A (zh) | 2019-07-23 |
CN110046469B true CN110046469B (zh) | 2020-02-11 |
Family
ID=67281684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910392663.5A Active CN110046469B (zh) | 2019-05-13 | 2019-05-13 | 多约束条件下水电站坝前河床冲淤变形的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110046469B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110837692B (zh) * | 2019-10-18 | 2022-08-19 | 河海大学 | 一种基于挟沙力变化的河漫滩平均高程确定方法 |
CN110795792B (zh) * | 2019-11-13 | 2020-06-23 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种由工程建设引发的河道强紊动区的河床变形预测方法 |
CN111008459B (zh) * | 2019-11-14 | 2021-09-28 | 天津大学 | 一种融合人工神经网络与水沙传输机制的河道水沙冲淤演变模拟方法 |
CN111241467B (zh) * | 2020-01-07 | 2022-06-03 | 河海大学 | 一种多沙河流水库排沙比的计算方法 |
CN112834471B (zh) * | 2021-01-04 | 2022-05-20 | 大连理工大学 | 一种基于物质浓度的时空梯度反演均匀扁平微流通道内平均流速的优化方法 |
CN112685918A (zh) * | 2021-02-04 | 2021-04-20 | 中建一局华江建设有限公司 | 一种基于无人机倾斜测量技术的河道清淤方法 |
CN113297753B (zh) * | 2021-06-24 | 2023-02-28 | 四川大学 | 山区河流透过性坝沿程水深预测方法 |
CN113821862B (zh) * | 2021-10-27 | 2023-08-08 | 杭州京杭运河二通道建设投资有限公司 | 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法 |
CN114578088B (zh) * | 2022-02-18 | 2024-04-12 | 长江水利委员会水文局长江三峡水文水资源勘测局 | 一种强约束河断面或龙口平均流速测量方法 |
CN115034468B (zh) * | 2022-06-07 | 2023-04-07 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种淤泥质海岸丁坝或防波堤拆除后海床冲淤变化预测方法 |
CN116205634A (zh) * | 2023-03-17 | 2023-06-02 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种水库大坝维修信息提取方法 |
CN117408173B (zh) * | 2023-12-16 | 2024-03-01 | 长江水利委员会水文局长江中游水文水资源勘测局 | 一种基于机器学习的水文流量整编智能模型构建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271672A (zh) * | 2018-08-21 | 2019-01-25 | 长江水资源保护科学研究所 | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105045954B (zh) * | 2015-06-09 | 2017-12-19 | 北京交通大学 | 一种陡坎河床洪水冲刷演变模拟方法及系统 |
US20170140842A1 (en) * | 2015-11-12 | 2017-05-18 | Westinghouse Electric Company Llc | Subcritical Reactivity Monitor Utilizing Prompt Self-Powered Incore Detectors |
CN105714727B (zh) * | 2016-01-21 | 2017-06-23 | 河海大学 | 一种弯曲河道河床冲淤变形的估算方法 |
CN108256177B (zh) * | 2017-12-28 | 2021-02-05 | 中国水利水电科学研究院 | 一种河流水沙模型的参数优化方法和系统 |
CN108334660B (zh) * | 2017-12-28 | 2021-07-20 | 中国水利水电科学研究院 | 一种基于数据同化的强冲积河流的水沙预测方法和系统 |
CN108647449B (zh) * | 2018-05-15 | 2022-01-28 | 长江水利委员会长江科学院 | 一种基于絮凝动力学的粘性泥沙运动数值模拟方法 |
-
2019
- 2019-05-13 CN CN201910392663.5A patent/CN110046469B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109271672A (zh) * | 2018-08-21 | 2019-01-25 | 长江水资源保护科学研究所 | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110046469A (zh) | 2019-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110046469B (zh) | 多约束条件下水电站坝前河床冲淤变形的计算方法 | |
CN110362925B (zh) | 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 | |
CN113723024B (zh) | 一种适用于滨海地区的“溪流”-“河道”-“河口”分布式洪水过程模拟方法 | |
CN111046574B (zh) | 一种平原湖区湖河闸泵系统防洪排涝演算方法 | |
CN107885958A (zh) | 一种平原感潮河网区纳污能力计算方法 | |
CN108629055B (zh) | 一种基于饱和输沙原理的沙质内河航道回淤量预报方法 | |
CN104933268B (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN111581828B (zh) | 一种感潮河段闸下水位流量关系的计算方法 | |
CN110046400A (zh) | 一种基于h-p型有限元法对混凝土重力坝进行静力特性分析的高精度数值模拟计算方法 | |
CN111428401A (zh) | 一种堰塞湖溃决过程模拟方法 | |
CN108532532A (zh) | 感潮河道防洪防潮岸线制定方法 | |
CN110532641A (zh) | 一种地表网格分层建模方法及系统 | |
CN112084671A (zh) | 城市时变增益降雨-径流过程模拟计算方法 | |
CN113742820B (zh) | 一种分洪隧洞出口消力池消能效果分析的数值模拟方法 | |
Pinho et al. | Mathematical modelling of salt water intrusion in a northern Portuguese estuary | |
CN114232557B (zh) | 跟踪截流模型水下地形的模拟确定方法和试验方法 | |
Zeng et al. | Evaluating the hydrological effects of the Three Gorges Reservoir based on a large-scale coupled hydrological-hydrodynamic-dam operation model | |
CN115994396A (zh) | 一种水系连通数学模型的构建方法 | |
CN104047255A (zh) | 人为引发河床演变对河流水位变化敏感性分析方法 | |
CN108959816A (zh) | 高拱坝施工初-中期导流风险率检测方法 | |
CN111340649B (zh) | 一种水系结构连通性的量测方法 | |
Tallón-Ballesteros | Operation Dispatching Feature Analysis Based on TELEMAC-2D | |
CN108491634A (zh) | 一种滩地防浪林种植不利于河道行洪的定量分析方法 | |
CN114459571B (zh) | 一种基于初态水位曲线簇河道初始水位确定方法 | |
CN113627050B (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 |