CN110909490B - 基于简化圣维南方程组的库区冰坝壅水模拟方法 - Google Patents
基于简化圣维南方程组的库区冰坝壅水模拟方法 Download PDFInfo
- Publication number
- CN110909490B CN110909490B CN201911257146.3A CN201911257146A CN110909490B CN 110909490 B CN110909490 B CN 110909490B CN 201911257146 A CN201911257146 A CN 201911257146A CN 110909490 B CN110909490 B CN 110909490B
- Authority
- CN
- China
- Prior art keywords
- ice
- dam
- ice dam
- water
- section
- 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
- E02B1/02—Hydraulic models
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Control Of Non-Electrical Variables (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于简化圣维南方程组的库区冰坝壅水模拟方法,S1、确定冰坝形成时期;S2、确定所述冰坝发生位置;S3、确定冰坝壅水计算所需库容曲线;S4、计算入库冰量及几何形态;S5、确定冰坝壅水计算所需泄流曲线;对不同流量级,考虑冰坝的阻水面积,得到该流量下的冰水混合面积,并根据断面形态反算水位,得到水位或冰位的流量关系,即为该断面的所述泄流曲线;S6、基于简化圣维南方程组进行冰坝调洪演算。本发明在有限资料条件下,利用库区坝前水位、库尾冰坝发生位置的关联关系及判别指标解决了冰坝发生位置依靠实测资料的难题。利用断面水力要素分析冰坝头部断面泄流能力,解决了断面过凌能力确定的难题。
Description
技术领域
本发明涉及水库防洪减灾技术领域,尤其是涉及基于简化圣维南方程组的库区冰坝壅水模拟方法。
背景技术
在北方河道流域地区,冰凌洪水频繁发生,往往造成十分严重的灾害,危及生命财产安全,影响社会经济发展。其中,冰坝作为一种特殊的冰凌灾害表现形式,广泛存在于冬季寒冷河道中,并显著壅高上游水位,成为冬季凌汛隐患和防凌重点。目前,在实际防凌工作中,常常需要通过数学模拟方法估算冰坝壅水情况,为冰凌洪水的前期防控提供技术支撑。
文献调研表明,由于冰坝突发性强,观测资料较少。目前工程上主要采用经验公式法、类比分析法、凌洪调洪演算法等计算冰坝壅水位。苏联P.B.多钦科等学者根据实验,认为冰坝壅高水位是冰或冰花堆积上游边缘河深的函数,并由此确定了冰坝壅水经验公式,通过该冰坝壅水经验公式可以计算无观测资料河段的冰坝水位,但此方法对不同河段的适应性较差。类比分析法,方法简单,但无法考虑水库运用前后水力要素变化的影响,精度较低。凌洪调洪演算法过于依靠实测资料来确定相关计算条件,导致实际应用存在瓶颈。除上述方法以外,对于库区的冰坝壅水模拟计算方法,特别是理论方法研究较少。因此,如何利用有限的资料条件,在水库库区的地形条件下,估算可能发生的冰坝壅水情况,及时调控,避免造成壅水淹没灾害,是本领域研究的重点和难点。
发明内容
本发明目的在于提供一种基于简化圣维南方程组的库区冰坝壅水模拟方法,用以解决北方河流防凌减灾和冰坝壅水防控问题,提高凌汛灾害应对能力。
为实现上述目的,本发明采取下述技术方案:
本发明所述基于简化圣维南方程组的库区冰坝壅水模拟方法,包括下述步骤:
S1、确定冰坝形成时期,即开河期的入库流量过程线,即:采用入库水文站开河期实测的最大洪峰流量作为最大值,采用典型年开河期流量过程作为典型洪峰流量过程,经同比例缩放得到所述开河期冰坝入库流量过程;
S2、确定所述冰坝发生位置,即:按照凌汛期开河流量,采用不同蓄水位,推算水库水面线,计算比降减小为整个河段比降的1/3所在的断面位置,获得不同蓄水位条件下的冰坝头部断面;
S3、确定冰坝壅水计算所需库容曲线,即:根据实测断面成果或设计淤积断面成果,将同级水位下断面面积和断面间距乘积,求得冰坝发生位置以上的库容曲线;
S4、计算入库冰量及几何形态,即:采用公式计算入库冰量,即:封河期蓄存于河道中的冰层在开河以后随着水流进入水库的冰块量;冰坝的几何形态采用概化的楔形形状,并假设冰坝在各断面分布为矩形,且在冰坝头部断面矩形面积达到最大;
S5、确定冰坝壅水计算所需泄流曲线,即:根据水库坝前水位,采用伯努利方程,选取适宜的流量级,向推算敞流条件下的回水水面线,获得库区断面的水位-面积-流量关系曲线;对不同流量级,考虑冰坝的阻水面积,得到该流量下的冰水混合面积,并根据断面形态反算水位,得到水位或冰位的流量关系,即为该断面的所述泄流曲线;
S6、基于简化圣维南方程组进行冰坝调洪演算,即:以冰坝最大尺寸作为调洪演算过程的初始条件,起调水位采用洪水过程起始流量对应的冰坝含冰水位,根据流量过程,冰坝头部断面对应的库容曲线和冰坝头部断面泄流曲线,进行冰坝调洪计算,得到最高壅水位。
本发明在有限资料条件下,利用库区坝前水位、库尾冰坝发生位置的关联关系及判别指标解决了冰坝发生位置依靠实测资料的难题。利用断面水力要素分析冰坝头部断面泄流能力,解决了断面过凌能力确定的难题。通过对冰坝入库流量过程、泄流曲线、库容曲线、入库冰量等的确定,基于简化圣维南方程组进行调洪演算,为凌汛多发河流水库的冰坝模拟提供了一种半经验半理论的计算方法。
本发明以位于黄河干流具有防凌、发电等综合利用要求的海勃湾水利枢纽为例,采用2016年实测库区断面、历史实测入库站石嘴山站流量等资料,对冰坝发生位置、冰坝入库冰量、入库流量过程、冰坝库容曲线等进行计算和分析,最后通过调洪计算,推算出冰坝壅水水面线,模拟精度达90%以上,为水库库区防凌提供技术支撑。
附图说明
图1是本发明实施例的方法流程图。
图2是本发明对冰坝几何形态的概化图。
图3是本发明计算冰坝壅水水位与冰坝断面面积示意图。
图4是本发明实施例海勃湾水库开河期入库流量过程线(石嘴山站)。
图5是本发明实施例冰坝头部断面以上库容曲线图。
图6是本发明实施例冰坝头部断面水位泄量关系曲线图。
图7是本发明实施例冰坝入流出流过程线图。
图8是本发明实施例冰坝壅水水位过程线图。
图9是本发明实施例冰坝壅水回水水面线图。
具体实施方式
下面结合附图对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述实施例。
如图1-3所示,本发明所述基于简化圣维南方程组的库区冰坝壅水模拟方法,现以黄河干流海勃湾水库为例作进一步详细说明:
S1、确定冰坝形成时期(开河期)入库流量过程线,即:采用海勃湾入库水文站石嘴山站开河期实测的最大洪峰流量1190m3/s(1995年2月18日)作为最大值,采用1975年开河期流量过程作为典型洪峰流量过程,经同比例缩放得到开河期冰坝入库流量过程;
S2、确定冰坝发生位置,即:按照凌汛期开河流量1190m3/s,基于海勃湾水库库区2016年实测断面条件,推求水库坝前防凌运用水位1073.5m时的水库回水水面线,计算比降减小为整个河段比降的1/3所在的断面位置,获得1073.5m水位下的冰坝头部断面;
S3、确定冰坝壅水计算所需库容曲线,即:根据海勃湾水库库区2016年实测断面成果,按照0.5m的间隔,将坝前水位范围1065.5m~1090m内不同水位条件下断面面积和断面间距乘积,求得冰坝发生位置以上的库容曲线;
S4、计算入库冰量及几何形态,即:采用公式计算海勃湾水库冰坝形成期的入库冰量;冰坝的几何形态采用概化的楔形形状,并假设冰坝在各断面分布为矩形,在冰坝头部断面矩形面积达到最大。其中,入库冰块量的计算公式如下:
式中:W为入库冰量,m3;
B为河段平均封河宽度,m;
H为河段内平均冰厚,m;
L为河段长度,km;
η为冰块沿程融化率;
S5、确定冰坝壅水计算所需泄流曲线,即:根据海勃湾水库坝前防凌运用水位1073.5m,采用伯努利方程,按照100m3/s的间隔,对设计流量范围100m3/s~2000m3/s内不同流量级,推算敞流条件下的回水水面线,根据冰坝头部断面位置,获得冰坝头部的水位~面积~流量关系曲线;对不同流量级,考虑冰坝的阻水面积,得到该流量下的冰水混合面积,并根据断面形态反算水位,得到水(冰)位~流量关系,即为冰坝头部泄流曲线。其中,断面冰水总面积的计算公式如下:
S6、基于简化圣维南方程组进行冰坝调洪演算,即:以冰坝最大尺寸作为调洪演算过程的初始条件,起调水位采用洪水过程起始流量665m3/s对应的冰坝头部含冰水位1074.86m,根据入库流量过程线、冰坝头部库容曲线、冰坝头部泄流曲线等条件进行冰坝调洪计算,并得到最高壅水位。
最后通过调洪计算,推算出冰坝壅水水面线,最高位置壅水2.3m。通过与经验公式等其他算法对比,模拟精度可达90%以上,可为海勃湾水库库区防凌提供技术支撑。
本实施例分析如下:
1、入库流量过程线
采用石嘴山水文站开河期实测的最大洪峰流量1190m3/s(1995年2月18日)作为最大值,采用1975年开河期流量过程作为典型洪峰流量过程,经同比例缩放得到入库流量过程线,见图4。
2、冰坝头部断面位置及库容曲线
采用开河期设计流量1200m3/s(接近实测最大洪峰流量1190m3/s),按2016年海勃湾水库库区实测断面,推求水库防凌运用水位1073.5m时的水库回水水面线,见表1。
表1 海勃湾水库水面线成果表
海勃湾库区比降约为0.04%,由表1成果可以看出,水面比降从上游向下游逐渐减小,D20断面处比降减小至0.013%以下,即总比降的1/3,因此可假定D20断面为冰坝可能发生的断面。根据D20断面及其上游的实测断面形态,计算其库容曲线,见表2和图5;
表2
3、冰坝头部断面泄流曲线
根据上述分析假定,现状条件下在水位为1073.5m时冰坝头部断面位置约在D20号断面,需计算该断面水位泄量关系。
本模型拟采用水力学法进行分析计算:首先采用伯努利方程计算各断面的水位-面积-流量关系,在不同流量级下,考虑冰坝的阻水面积,得到该流量下的冰水混合面积,并根据断面形态反算水位,得到水(冰)位~流量关系,即为该断面冰坝条件下的水位流量关系。
1)、断面水位流量关系(无冰)
采用伯努利方程,分别计算不同起始水位条件下,各个流量级在冰坝头部断面的水库回水水位,见表3。
表3 各流量级回水水位及断面水利要素
2)、冰坝水位流量关系
表4 冰坝计算过程表
4、起调水位
以冰坝最大尺寸为调洪过程的初始条件,根据洪水过程线起始流量665m3/s,查冰坝头部断面D20断面泄流曲线,确定调洪的起调水位1074.86m。
5、调洪成果
凌洪调洪计算过程见表5,成果见表6及图7-图8。由调洪结果可以看出,现状条件下的水库运用初期,当水库坝前水位为防凌运用水位1073.5m时,发生流量为1190m3/s的凌峰,可能会在D20号断面(距坝里程17km处)发生冰坝,并向上游延伸,冰坝阻水,造成最高冰坝壅水位1075.71m,水位较敞流水面壅高1.84m。
表5
续表5 凌洪调洪计算过程表
表6
6、冰坝水面线
根据冰坝壅水位分析,水库运行初期,当入库流量为1190m3/s,水库坝前水位为1073.5m时,冰坝头部断面在D20断面(铁路桥附近),壅水水位1075.71m,略低于正常蓄水位。海勃湾水利枢纽冰坝壅水回水水面线见表7,图9。
表7 冰坝壅水回水水面线
7、合理性分析
计算经验公式法,进行验证。根据苏联P.B.多钦科等研究,由经验公式计算得到的冰坝壅水高度1.5m,本发明方法计算成果为1.8m,二者相差不大,但本发明方法可以提供冰坝发展的全过程模拟,结果更加可靠。
Claims (1)
1.一种基于简化圣维南方程组的库区冰坝壅水模拟方法,其特征在于:包括下述步骤:
S1、确定冰坝形成时期,即开河期的入库流量过程线,即:采用入库水文站开河期实测的最大洪峰流量作为最大值,采用典型年开河期流量过程作为典型洪峰流量过程,经同比例缩放得到所述开河期冰坝入库流量过程;
S2、确定所述冰坝发生位置,即:按照凌汛期开河流量,采用不同蓄水位,推算水库水面线,计算比降减小为整个河段比降的1/3所在的断面位置,获得不同蓄水位条件下的冰坝头部断面;
S3、确定冰坝壅水计算所需库容曲线,即:根据实测断面成果或设计淤积断面成果,将同级水位下断面面积和断面间距乘积,求得冰坝发生位置以上的库容曲线;
S4、计算入库冰量及几何形态,即:采用公式计算入库冰量,即:封河期蓄存于河道中的冰层在开河以后随着水流进入水库的冰块量;冰坝的几何形态采用概化的楔形形状,并假设冰坝在各断面分布为矩形,且在冰坝头部断面矩形面积达到最大;
S5、确定冰坝壅水计算所需泄流曲线,即:根据水库坝前水位,采用伯努利方程,选取适宜的流量级,向推算敞流条件下的回水水面线,获得库区断面的水位-面积-流量关系曲线;对不同流量级,考虑冰坝的阻水面积,得到该流量下的冰水混合面积,并根据断面形态反算水位,得到水位或冰位的流量关系,即为该断面的所述泄流曲线;
S6、基于简化圣维南方程组进行冰坝调洪演算,即:以冰坝最大尺寸作为调洪演算过程的初始条件,起调水位采用洪水过程起始流量对应的冰坝含冰水位,根据流量过程,冰坝头部断面对应的库容曲线和冰坝头部断面泄流曲线,进行冰坝调洪计算,得到最高壅水位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911257146.3A CN110909490B (zh) | 2019-12-10 | 2019-12-10 | 基于简化圣维南方程组的库区冰坝壅水模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911257146.3A CN110909490B (zh) | 2019-12-10 | 2019-12-10 | 基于简化圣维南方程组的库区冰坝壅水模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110909490A CN110909490A (zh) | 2020-03-24 |
CN110909490B true CN110909490B (zh) | 2023-03-24 |
Family
ID=69824107
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911257146.3A Active CN110909490B (zh) | 2019-12-10 | 2019-12-10 | 基于简化圣维南方程组的库区冰坝壅水模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110909490B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112287528B (zh) * | 2020-10-16 | 2021-07-20 | 黄河勘测规划设计研究院有限公司 | 一种基于高滩高槽的多沙河流水库防洪高水位确定方法 |
CN112418491B (zh) * | 2020-11-06 | 2021-08-03 | 黄河勘测规划设计研究院有限公司 | 一种水库剩余拦沙库容动态配置方法 |
CN113139230B (zh) * | 2021-04-30 | 2023-05-16 | 黄河勘测规划设计研究院有限公司 | 城市中小型偏流河道潜坝减淤设计方法 |
CN114169047B (zh) * | 2021-11-23 | 2024-01-23 | 中国长江三峡集团有限公司 | 一种计算水库回水范围的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2304659C1 (ru) * | 2006-02-21 | 2007-08-20 | Федеральное государственное унитарное предприятие "Дальневосточный научно-исследовательский институт гидротехники и мелиорации" (ФГУП ДальНИИГиМ) | Способ регулирования речного стока для предупреждения наводнений |
CN101619972A (zh) * | 2009-08-12 | 2010-01-06 | 国网电力科学研究院 | 一种水库蓄水量的测算方法 |
CN101956382A (zh) * | 2010-04-08 | 2011-01-26 | 潜江市知识产权协会 | 用槽蓄增量关系求解圣维南方程组模拟洪水演进的方法 |
CN107761654A (zh) * | 2017-10-31 | 2018-03-06 | 中国电建集团成都勘测设计研究院有限公司 | 梯级建设条件下的高坝施工中期防洪度汛方案决策方法 |
CN108256780A (zh) * | 2018-02-05 | 2018-07-06 | 贵州乌江水电开发有限责任公司 | 一种基于迭代库区回水计算的梯级水库防洪调度计算方法 |
CN108399312A (zh) * | 2018-04-16 | 2018-08-14 | 江苏省水利科学研究院 | 基于桥梁上下游水位差的桥墩绕流阻力等效糙率概化方法 |
-
2019
- 2019-12-10 CN CN201911257146.3A patent/CN110909490B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2304659C1 (ru) * | 2006-02-21 | 2007-08-20 | Федеральное государственное унитарное предприятие "Дальневосточный научно-исследовательский институт гидротехники и мелиорации" (ФГУП ДальНИИГиМ) | Способ регулирования речного стока для предупреждения наводнений |
CN101619972A (zh) * | 2009-08-12 | 2010-01-06 | 国网电力科学研究院 | 一种水库蓄水量的测算方法 |
CN101956382A (zh) * | 2010-04-08 | 2011-01-26 | 潜江市知识产权协会 | 用槽蓄增量关系求解圣维南方程组模拟洪水演进的方法 |
CN107761654A (zh) * | 2017-10-31 | 2018-03-06 | 中国电建集团成都勘测设计研究院有限公司 | 梯级建设条件下的高坝施工中期防洪度汛方案决策方法 |
CN108256780A (zh) * | 2018-02-05 | 2018-07-06 | 贵州乌江水电开发有限责任公司 | 一种基于迭代库区回水计算的梯级水库防洪调度计算方法 |
CN108399312A (zh) * | 2018-04-16 | 2018-08-14 | 江苏省水利科学研究院 | 基于桥梁上下游水位差的桥墩绕流阻力等效糙率概化方法 |
Non-Patent Citations (2)
Title |
---|
万家寨水库开河期冰坝壅水计算及影响因素分析;钱云平等;《人民黄河》;20020330(第03期);全文 * |
河道凌期冰塞壅水计算;范中亮等;《中国水运(下半月)》;20150915(第09期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110909490A (zh) | 2020-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110909490B (zh) | 基于简化圣维南方程组的库区冰坝壅水模拟方法 | |
Wu et al. | Effect of altered flow regime on bankfull area of the Lower Yellow River, China | |
Han et al. | Interacting effects of multiple factors on the morphological evolution of the meandering reaches downstream the Three Gorges Dam | |
CN113792975B (zh) | 水工程联合调度防洪减灾效益评估方法 | |
CN107386200A (zh) | 一种适用于分汊河段航道整治的丁坝尺寸确定方法 | |
Tate et al. | Simulating reservoir management under the threat of sedimentation: the case of Tarbela dam on the river Indus | |
CN111881592A (zh) | 一种考虑河床锚冰的河冰数值模拟方法 | |
Barbé et al. | Effects on dredging due to diversions from the lower Mississippi River | |
CN113836477B (zh) | 基于平滩流量的稳封期河道防凌安全过流量计算方法 | |
Ikramov et al. | The height of a damless water intake structure threshold | |
Wiejaczka et al. | Water temperature dynamics in a complex of reservoirs and its effect on the temperature patterns of a mountain river | |
Zeng et al. | Application of Computational Fluid Dynamics Modeling to Study Maximum Allowable Gate Openings at Spillways | |
Yean et al. | Sedimentation in the lower reaches of the Yellow River and its basic laws | |
Widiarti et al. | The hydraulics analysis of the USBR stilling basin type of a dam under multiple flow discharge scenarios (A case study of Temef Dam, East Nusa Tenggara, Indonesia) | |
Lagasse | River response to dredging | |
CN112836450B (zh) | 流道断面体型及其设计方法、存储介质和计算机设备 | |
Altawash et al. | Velocity Patterns inside the Proposed Makhool Dam Reservoir with Different Operation Plans | |
CN111931276B (zh) | 一种快速优选跌坎消力池体型的方法 | |
Verelst et al. | Erosion and sedimentation near the renovated weir lock complex of Asper in Belgium | |
Cassidy | Fluid mechanics and design of hydraulic structures | |
Ji et al. | Branch Channel Evolution and Its Effects on Tail Channel in the Yellow River Estuary | |
Raza et al. | Exploring sediment management options of Mangla Reservoir using RESSASS | |
Yapa et al. | Effect of ice cover on hydropower production | |
Majewski | The lower vistula and its ice problems | |
Shuiping | Study on the Determination Method of Water Surface Ratio for Transfer Flooding and Drainage in Dyke Type River Network Area |
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 |