CN110263428B - 一种基于流量加权平均流线长度指标的河床演变分析方法 - Google Patents
一种基于流量加权平均流线长度指标的河床演变分析方法 Download PDFInfo
- Publication number
- CN110263428B CN110263428B CN201910527635.XA CN201910527635A CN110263428B CN 110263428 B CN110263428 B CN 110263428B CN 201910527635 A CN201910527635 A CN 201910527635A CN 110263428 B CN110263428 B CN 110263428B
- Authority
- CN
- China
- Prior art keywords
- flow
- streamline
- weighted average
- length
- len
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mechanical Engineering (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于流量加权平均流线长度指标的河床演变分析方法,其根据对象河段的二维流场和位于二维流场出口断面处的流线终点分别计算得出该对象河段在N个时期内的流量加权平均流线长度E_LEN,对比N个时期内的流量加权平均流线长度E_LEN实现对象河段的定量分析。本发明能够实现对河道主流、深泓摆动情况进行定量刻画。
Description
技术领域
本发明属于河床演变分析方法,是一种用于定量分析河道形态与水动力条件平面变化的方法,适用于对平面形态复杂、冲淤频繁的平面河流的主流、深泓进行定量刻画。
背景技术
以长江中下游平原河流为例,目前主要通过河道长度、河道内等高线(等深线)、深泓线等分析平面形态变化;以主流区、主流线等分析水动力条件平面变化。受限于河道形态以及水动力条件的复杂性,上述指标均有一定的局限性。比如堤岸受控情况下,河道长度是基本固定的,无法反映河道内滩槽冲淤调整变化;等高线(等深线)需结合具体的河段进行针对性选择,欠缺一般性,平面形态描述上也难以定量;严格定义上的深泓线、主流区、主流线等指标在天然河流平滑程度差,甚至不连续,往往需要进行经验性的判断处理,精确度与严谨性较为欠缺。
受限于现有指标体系的欠缺,河床演变中的“取直”、“坐弯”等概念定量刻画不足。尤其是在大型枢纽修建以后下游河床的剧烈演变期,不同河型的平面调整缺乏一致性描述,同样以长江中下游为例,三峡水库蓄水运行以后,弯曲河段凸冲凹淤、汊道段短汊发育等现象,仍只能定性阐述,无法对其进行定量分析。
发明内容
本发明针对现有技术的不足,提出了一种基于流量加权平均流线长度指标的河床演变分析方法,能够实现对河道主流、深泓摆动情况进行定量刻画。
本发明公开了一种基于流量加权平均流线长度指标的河床演变分析方法,其根据对象河段的二维流场和位于二维流场出口断面处的流线终点分别计算得出该对象河段在N个时期内的流量加权平均流线长度E_LEN,对比N个时期内的流量加权平均流线长度E_LEN实现对象河段的定量分析。
在本发明的一种优选实施方案中,获取每个时期内的流量加权平均流线长度E_LEN的具体步骤:包括步骤一,对象河段的二维流场建模;步骤二,在二维流场的出口位置合理布置多个流线终点;步骤三,基于二维流场数据反算每个流线终点所对应的流线长度LENi;步骤四,监测各流线间的过流流量ΔQi;步骤五,根据LENi和ΔQi计算得到流量加权平均流线长度E_LEN,
在本发明的一种优选实施方案中,步骤一中,建立分析对象河段的平面二维水流数学模型,在完成水位、流速验证的基础上,模拟不同测次地形下洪、中、枯不同流量级的二维流场。
在本发明的一种优选实施方案中,步骤二中,流线终点总数不少于10个。
在本发明的一种优选实施方案中,步骤二中,任意两个相邻流线终点之间的间距大于二维网格流场的节点间距。
在本发明的一种优选实施方案中,步骤三中,以每个流线终点为起点、ΔL为步长,每一步的溯点均取自计算网格的结点,并保证ΔL的前进方向与流速矢量平行,最大限度地逼近经过各流线终点的完整流线。
在本发明的一种优选实施方案中,步骤四中,ΔQi=ΔQi+1=恒定值。
在本发明的一种优选实施方案中,若前一时期的流量加权平均流线长度与后一时期的流量加权平均流线长度之间的差值为ΔE_LEN;当ΔE_LEN>0时,则流量加权平均流线长度减小,对应的河道深泓、主流撇弯取直,流程缩短;当ΔE_LEN<0时,则流量加权平均流线长度增加,对应的河道深泓、主流坐弯,流程增加。
本发明的有益效果时:本发明通过对对象河段的二维流场、流线终点的分析得出流量加权平均流线长度从而实现了对对象河段的定量分析,提高了对象河段的深泓线、主流区、主流线等指标的精确分析。
附图说明
图1本发明一种基于流量加权平均流线长度指标的河床演变分析方法的流程图。
图2本发明一种基于流量加权平均流线长度指标的河床演变分析方法的二维流场示意图。
图3本发明实施例马家咀河段2018年24000m3/s流量级下出口位置的流线终点总数为30个。
图4本发明实施例熊家洲至城陵矶河段2017年24000m3/s流量级下出口位置的流线终点总数为25个。
具体实施方式
具体实施方式如下:
本发明对于分析对象河段,利用平面二维水流数学模型模拟不同测次地形下的二维流场,即水流的二维欧拉场,根据二维流场,计算出二维流线分布。由于流线之间流量恒定,可对所有流线按各流线邻域内流量取流线加权平均值,作为整个河道某一测次地形条件下的流量加权平均流线长度。该指标融合了地形与水动力条件因素,较河道中轴线长度、深泓线长度等指标有更为清晰的物理意义。加权平均流线长度减小,对应了河道深泓、主流撇弯取直,流程缩短;加权平均流线长度增大,对应了河道深泓、主流坐弯,流程增加。以长江中游微弯分汊型河段马家咀河段、典型蜿蜒型河段熊家洲至城陵矶河段为例,本发明具体的实现步骤如下:
(1)建立分析对象河段的平面二维水流数学模型,在完成水位、流速验证的基础上,模拟不同测次地形下洪、中、枯不同流量级的二维流场。进行年际变化分析时,各测次地形应尽量对应每年的同一水文时期。马家咀河段分别模拟了2010、2012、2014、2018共四个年份枯水期测图的各级流场,熊家洲至城陵矶河段分别模拟了2010、2012、2014、2017共四个年份枯水期测图的各级流场。
(2)在二维流场出口位置结合水深、流速情况沿断面布置流线终点。水深、流速较大的位置相对密集布置,较小的位置相对稀疏布置;流线终点总数不少于10个,但点距不宜小于二维网格流场的节点间距。
比如马家咀河段2018年24000m3/s流量级下出口位置的流线终点总数为30个;熊家洲至城陵矶河段2017年24000m3/s流量级下出口位置的流线终点总数为25个。
(3)根据二维流场计算结果,由选定的流线终点溯源而上。以ΔL为计算步长,每一步的溯点均取自计算网格的结点,并保证每一段ΔL的前进方向与流速矢量平行,最大限度地逼近经过各布置终点的完整流线。而后通过将各步长点的坐标(Xj,Yj)导出,计算第j条流线的步长采取分段相加的方法计算第i条流线的总长LENi,以此提取所有布设的流线长度。
马家咀河段2018年24000m3/s流量级下各流线LENi计算结果如表1;熊家洲至城陵矶河段2017年30000m3/s流量级下各流线LENi计算结果如表2。
表1马家咀河段2018年24000m3/s流量级下流线长度计算过程
表2熊家洲至城陵矶河段2017年30000m3/s流量级下流线长度计算过程
(4)流线之间流量不变,计算流线li、li+1之间的流量Qi,再计算流量加权平均流线长度E_LEN。马家咀河段2018年24000m3/s流量级下E-LEN计算结果如表1;熊家洲至城陵矶河段2017年30000m3/s流量级下E_LEN计算结果如表2。
(5)对比不同测次地形条件下各流量级的E_LEN,若E_LEN减小,即加权平均流线长度减小,对应了河道深泓、主流撇弯取直,流程缩短。
以长江中游典型蜿蜒型河段熊家洲至城陵矶河段为例,30000m3/s流量的加权平均流线长度从2010年的33.66km减小至2017年的33.21km,这定量刻画了该弯曲河段凸岸冲刷、凹岸淤积的切滩撇弯取直过程,若E_LEN增大,即加权平均流线长度增大,对应了河道深泓、主流坐弯,流程增加。
以长江中游微弯分汊型河段马家咀河段为例,24000m3/s流量的加权平均流线长度从2010年的15.65km增大至2018年的15.86km,这定量刻画了该分汊河段左汊(较短的汊道)因工程控制逐渐淤塞,而深泓主流逐渐向右汊(较长的汊道)集中,进而造成主流、深泓坐弯的现象。
由上述实例可以看出,因变化幅度相对于总长度偏小,若无加权平均流向长度这一严格的统一处理指标,人为勾勒主流线、深泓线,引起的误差容易掩盖河道真实的变化过程。
Claims (6)
1.一种基于流量加权平均流线长度指标的河床演变分析方法,其特征在于:根据对象河段的二维流场和位于二维流场出口断面处的流线终点分别计算得出该对象河段在N个时期内的流量加权平均流线长度E_LEN,对比N个时期内的流量加权平均流线长度E_LEN实现对象河段的定量分析;
获取每个时期内的流量加权平均流线长度E_LEN的具体步骤包括
步骤一,对象河段的二维流场建模;
步骤二,在二维流场的出口位置合理布置多个流线终点,流线终点总数不少于10个;
步骤三,基于二维流场数据反算每个流线终点所对应的流线长度LENi;
步骤四,监测各流线间的过流流量ΔQi;
步骤五,根据LENi和ΔQi计算得到流量加权平均流线长度E_LEN,
若前一时期的流量加权平均流线长度与后一时期的流量加权平均流线长度之间的差值为ΔE_LEN;当ΔE_LEN>0时,则流量加权平均流线长度减小,对应的河道深泓、主流撇弯取直,流程缩短;当ΔE_LEN<0时,则流量加权平均流线长度增加,对应的河道深泓、主流坐弯,流程增加。
2.根据权利要求1所述的一种基于流量加权平均流线长度指标的河床演变分析方法,其特征在于:步骤一中,建立分析对象河段的平面二维水流数学模型,在完成水位、流速验证的基础上,模拟不同测次地形下洪、中、枯不同流量级的二维流场。
3.根据权利要求1所述的一种基于流量加权平均流线长度指标的河床演变分析方法,其特征在于:步骤二中,任意两个相邻流线终点之间的间距大于二维网格流场的节点间距。
4.根据权利要求1所述的一种基于流量加权平均流线长度指标的河床演变分析方法,其特征在于:步骤三中,以每个流线终点为起点、ΔL为步长,每一步的溯点均取自计算网格的结点,并保证ΔL的前进方向与流速矢量平行,最大限度地逼近经过各流线终点的完整流线。
6.根据权利要求1所述的一种基于流量加权平均流线长度指标的河床演变分析方法,其特征在于:步骤四中,ΔQi=ΔQi+1=恒定值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910527635.XA CN110263428B (zh) | 2019-06-18 | 2019-06-18 | 一种基于流量加权平均流线长度指标的河床演变分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910527635.XA CN110263428B (zh) | 2019-06-18 | 2019-06-18 | 一种基于流量加权平均流线长度指标的河床演变分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110263428A CN110263428A (zh) | 2019-09-20 |
CN110263428B true CN110263428B (zh) | 2023-01-24 |
Family
ID=67919078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910527635.XA Active CN110263428B (zh) | 2019-06-18 | 2019-06-18 | 一种基于流量加权平均流线长度指标的河床演变分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110263428B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112802069B (zh) * | 2020-12-30 | 2023-07-04 | 北京师范大学 | 一种浅水湖泊阻水围堤围埝的判定方法 |
CN112925040B (zh) * | 2021-01-14 | 2022-11-01 | 中国石油天然气股份有限公司 | 砂岩储层渗透率各向异性参数预测方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09185702A (ja) * | 1996-01-04 | 1997-07-15 | Hitachi Ltd | 河床形状補間方法および三次元情報補間方法および中線算出方法および外分線算出方法 |
KR20100044371A (ko) * | 2008-10-22 | 2010-04-30 | 연세대학교 산학협력단 | 자연하천의 하상 변동 수치해석 방법 및 그 장치 |
CN107908888A (zh) * | 2017-11-22 | 2018-04-13 | 河海大学 | 一种感潮河段生态环境需水量计算方法 |
CN108038574A (zh) * | 2017-12-15 | 2018-05-15 | 武汉大学 | 冲积河道深弘摆动预测及崩岸预警的方法 |
CN109137816A (zh) * | 2018-08-23 | 2019-01-04 | 湖北省水利水电规划勘测设计院 | 一种基于主流摆型波震荡机制的河型划分方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8348553B2 (en) * | 2010-06-01 | 2013-01-08 | Applied University Research, Inc. | Bridge pier and abutment scour preventing apparatus with vortex generators |
CN106759063B (zh) * | 2016-12-05 | 2018-08-28 | 河海大学 | 一种感潮河段造床流量的计算方法 |
CN106682317B (zh) * | 2016-12-29 | 2019-10-25 | 武汉大学 | 上下游边界共同控制下河段平滩流量预测模型的构建方法 |
-
2019
- 2019-06-18 CN CN201910527635.XA patent/CN110263428B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09185702A (ja) * | 1996-01-04 | 1997-07-15 | Hitachi Ltd | 河床形状補間方法および三次元情報補間方法および中線算出方法および外分線算出方法 |
KR20100044371A (ko) * | 2008-10-22 | 2010-04-30 | 연세대학교 산학협력단 | 자연하천의 하상 변동 수치해석 방법 및 그 장치 |
CN107908888A (zh) * | 2017-11-22 | 2018-04-13 | 河海大学 | 一种感潮河段生态环境需水量计算方法 |
CN108038574A (zh) * | 2017-12-15 | 2018-05-15 | 武汉大学 | 冲积河道深弘摆动预测及崩岸预警的方法 |
CN109137816A (zh) * | 2018-08-23 | 2019-01-04 | 湖北省水利水电规划勘测设计院 | 一种基于主流摆型波震荡机制的河型划分方法 |
Non-Patent Citations (4)
Title |
---|
Response of reach-scale bankfull channel geometry in the Lower Yellow River to the altered flow and sediment regime;Xia Junqiang;《Geomorphology》;20141231;第255-265页 * |
荆江石首河段近50年河床演变分析;周美蓉等;《泥沙研究》;20170215(第01期);第43-49页 * |
长江中下游砂卵石分汊河段演变机理;刘林双等;《武汉大学学报(工学版)》;20170601(第03期);第28-35、76页 * |
黄河下游游荡段河道平面形态与河势变化趋势预测;陆中臣等;《地理学报》;20000623(第06期);第91-98页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110263428A (zh) | 2019-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Extension of a GIS procedure for calculating the RUSLE equation LS factor | |
CN110906992B (zh) | 基于水平adcp施测垂线流速分布的河流流量测量方法 | |
CN110263428B (zh) | 一种基于流量加权平均流线长度指标的河床演变分析方法 | |
CN108021780B (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN111241758B (zh) | 基于可溶性污染物在水环境中输移扩散模型的评估方法 | |
CN109271672B (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN110135056B (zh) | 一种河道内桥墩有效阻水比的快速自动分析方法 | |
CN112464584B (zh) | 自由表面流的水位和流量推求方法 | |
Niazkar et al. | Assessment of artificial intelligence models for estimating lengths of gradually varied flow profiles | |
Hodskinson | Computational fluid dynamics as a tool for investigating separated flow in river bends | |
CN115455867B (zh) | 基于回归分析的坝区流态的推求方法 | |
CN110847112B (zh) | 一种基于水力学模拟的河道泄洪预警方法 | |
CN117648878A (zh) | 一种基于1d-cnn算法的洪水快速演进及淹没模拟方法 | |
CN113327323B (zh) | 基于散点数据的水体环境地形构建方法 | |
CN113639805B (zh) | 一种基于渠道断面流速场重建的流量测量方法 | |
Thakur et al. | Exploring CCHE2D and its sediment modelling capabilities | |
CN205373998U (zh) | 一种基于压力传感器的薄层水流滚波测量装置 | |
CN114757049A (zh) | 流域污水处理厂提标改造的必要性分析验证方法 | |
CN108763732B (zh) | 一种分流砂坝型三角洲前缘训练图像建立方法 | |
CN114297953B (zh) | 一种山区大风灾害风险区划方法 | |
CN111291457B (zh) | 滩地有淹没植被的弯曲复式河道水深平均二维水流流向预测方法 | |
He et al. | Modeling bedload transport trajectories along a Sine-generated channel | |
Abed et al. | Detect the Water Flow Characteristic of the Tigris River by Using HFC-RAS Program | |
CN108345738A (zh) | 一种中小流域暴雨洪水汇流产流模型参数自我率定方法 | |
CN114547869B (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 |