CN106295056B - 一种复式断面河道河漫滩和主槽的自动识别方法 - Google Patents

一种复式断面河道河漫滩和主槽的自动识别方法 Download PDF

Info

Publication number
CN106295056B
CN106295056B CN201610700632.8A CN201610700632A CN106295056B CN 106295056 B CN106295056 B CN 106295056B CN 201610700632 A CN201610700632 A CN 201610700632A CN 106295056 B CN106295056 B CN 106295056B
Authority
CN
China
Prior art keywords
water level
section
river
hydraulic radius
major trough
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
Application number
CN201610700632.8A
Other languages
English (en)
Other versions
CN106295056A (zh
Inventor
吴腾
吴玲莉
秦杰
朱瑞虎
丁坚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN201610700632.8A priority Critical patent/CN106295056B/zh
Publication of CN106295056A publication Critical patent/CN106295056A/zh
Application granted granted Critical
Publication of CN106295056B publication Critical patent/CN106295056B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • 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)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Sewage (AREA)

Abstract

本发明公开了一种复式断面河道河漫滩和主槽的自动识别方法，该方法选取复式断面的最低点高程和最高点高程，将其作为断面过流的最低水位和最高水位，从最低水位开始，每次抬升0.05m直至最高水位，计算每一级水位对应的断面过流面积和湿周，确定对应的水力半径，建立水位与水力半径的关系曲线。针对该曲线的特点，按照水位由低到高，计算相临两水位对应的水力半径的差值，设定水力半径变幅的临界值，捕捉水位与水力半径关系曲线中的突变，将突变点作为滩与槽、滩与滩的分界点。该方法能定量自动识别复式断面的主槽和河漫滩，且简易可行，可应用于准二维水沙数值模拟中，也可用于复式断面河道分滩槽的演变分析。

Description

一种复式断面河道河漫滩和主槽的自动识别方法
技术领域
本发明涉及一种复式断面河道准二维数学模型中河漫滩和主槽的自动识别方法，属于河床演变分析中的滩槽划分技术领域。
背景技术
泥沙数学模型是研究泥沙问题的重要手段之一，具有投资少、周期短等优点，泥沙数学模型按照维数可分为：一维数学模型、二维数学模型和三维数学模型。随着计算机技术的快速发展，计算的速度越来越快，这三种数学模型在实际工程中均有应用，具体而言，一维数学模型多用于长距离、长时段的计算；二维数学模型和三维数学模型适用于流态较为复杂的局部区域或者较短距离的计算。对于水库，常采用一维数学模型进行长时段的计算分析。北方多沙河流，尤其是黄河中下游河道多为主槽和河漫滩组成的复式河槽，滩和槽上的水深不同，水力要素的数值也相差较多，为此常采用准二维数学模型进行求解。
准二维数学模型的一个重要步骤是划分滩槽，常用的方法是在进行模拟计算前根据断面图确定滩槽的位置和高程，划分过程中具有一定的人为性，并且划分完后在模拟计算中滩槽的位置和宽度保持不变，只有滩槽的高程发生变化，经过长时间的计算后断面形态可能发生较大变化，如果仍然沿用初始划分已难以符合实际情况，相应的糙率等参数都要进行相应的调整，因此，在计算中需要根据地形的变化实时调整滩槽的位置、宽度以及高程。此外，长河段的计算中，上下游的河型并不相同，上游可能是单一的河槽，到了中游和下游则是复式河槽，无论是用单一河槽的一维数学模型还是准二维数学模型进行模拟都不太适合。如果在准二维数学模型中加入滩槽自动划分的功能则可以将两种模型统一起来，可以在上游只划分一个子断面，而河道的下游可划分多个子断面，适应计算需求。有学者根据河道断面平均流速随水深的变化特征来定量判别主槽和河漫滩的位置，并应用于准二维数学模型中，但流速的计算过于复杂，需要率定的参数也较多，不易应用到准二维数学模型中。此外，也有学者通过流量模K(K＝AR2/3)来定量划分主边槽，由于流量模为断面面积A与水力半径R的组合，实际河道断面面积远大于水力半径，两者的乘积主要体现断面面积A的特征，断面面积A随高程递增，当河漫滩较宽时，流量模识别滩槽的能力极大降低，甚至难以识别。
发明内容
本发明所要解决的技术问题是：提供一种复式断面河道河漫滩和主槽的自动识别方法，该方法能定量自动识别复式断面的主槽和河漫滩，计算简单，可操作性强。
本发明为解决上述技术问题采用以下技术方案：
一种复式断面河道河漫滩和主槽的自动识别方法，包括如下步骤：
步骤1，获取河道断面的起点距和对应的高程数据，并查找最低点高程和最高点高程，将其分别作为计算河道断面水力半径的最低水位和最高水位；
步骤2，从最低水位开始，每次抬高水位0.05m，若抬高后的水位大于或等于最高水位时，将最高水位作为最后一个水位，依次记录所有的水位数据；
步骤3，对于步骤2记录的每个水位，计算河道断面在该水位的过水面积以及河道断面的湿周，根据过水面积和湿周计算河道断面在该水位的水力半径；
步骤4，水位由低至高，计算相邻两个水位对应的水力半径的差值ΔR0(i)＝R0(i+1)-R0(i)，i＝1,2,…,n-1；当ΔR0(i)＜0时，记录RR(p)＝R0(i)，Z1(p)＝Zlevel(i)，Z2(p)＝Zlevel(i+1)；当连续出现m个ΔR0(i)＜0且m≥2时，计算水力半径的变幅：同时记录Z1(p)＝Zlevel(i)，Z2(p)＝Zlevel(i+m)；p＝p+1，p＝1,2,…；其中，R0(·)表示水力半径，Zlevel(·)表示第·个水位的水位值，n表示所有水位的个数；
步骤5，设定河漫滩引起的河段水力半径的临界最小变幅ξ，当RR(p)≥ξ时，认为RR(p)对应的水位处出现河漫滩，低于最低河漫滩的河道断面部分为主槽。
作为本发明的一种优选方案，步骤2所述抬高水位的计算公式为：
Zlevel(i)＝Zmin+(i-1)*0.05
其中，Zlevel(i)表示第i个水位的水位值，Zmin表示最低水位的水位值，水位单位为：m，i＝1,2,…,n，n表示所有水位的个数。
作为本发明的一种优选方案，步骤3所述过水面积的计算方法为：对于任意水位，根据该水位河道断面的地形数据点，将河道断面分成若干梯形和三角形，所有梯形和三角形的面积之和即为过水面积。
作为本发明的一种优选方案，步骤5所述河漫滩引起的河段水力半径的临界最小变幅ξ＝0.1。
作为本发明的一种优选方案，步骤5所述河漫滩的代表高程为根据代表高程及河道断面地形得到该高程对应的河道宽度。
本发明采用以上技术方案与现有技术相比，具有以下技术效果：
本发明复式断面河道河漫滩和主槽的自动识别方法，采用水力半径定量自动识别复式断面的主槽和河漫滩，物理概念清晰，计算简单，可操作性强，能直接应用于准二维水沙数学模型的模拟计算，也能用于复式断面河道的河床演变分析，克服了人为划分滩槽带来的误差。
附图说明
图1是本发明复式断面河道河漫滩和主槽的自动识别方法的流程图。
图2是本发明实施例黄河花园口河道断面形态图。
图3是本发明实施例黄河花园口河道断面过水面积及湿周计算示意图。
图4是本发明实施例黄河花园口河道断面水力半径随水位变化图。
图5是本发明实施例黄河花园口河道断面滩槽识别结果图。
具体实施方式
下面详细描述本发明的实施方式，所述实施方式的示例在附图中示出。下面通过参考附图描述的实施方式是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
单一河槽断面水力半径随水位变化曲线是连续的，但复式河槽在漫滩瞬间水力半径迅速变小然后随水深逐渐增大，根据这种变化特性可以确定河漫滩的位置和高程。如图1所示，本发明复式断面河道河漫滩和主槽的自动识别方法的流程如下：
(1)获取河道断面地形点数据，河道断面左侧第一点为起始点，从左到右逐渐排列，记录每一点的起点距Di(至左侧第一点的距离)和河床高程Zi；
(2)查找河道断面河床高程最低的点高程Zmin和最高点高程Zmax，分别作为计算断面过水面积的最低水位和最高水位；
(3)从最低水位开始，作为初始计算水位，每次抬高水位0.05m，计算水位为：
Zlevel(i)＝Z min+(i-1)*0.05
其中，Zlevel(i)表示第i个水位的水位值，Z min表示最低水位的水位值，i＝1,2,…,n，n表示所有水位的个数，当Z(i)≥Zmax时，Zlevel(i)＝Zmax，该点为最后一个水位点；
(4)对于任意计算水位Zlevel(i)，根据河道断面的地形数据点，将河道断面分成若干梯形和三角形，计算河道断面在该水位下的过水面积A0(i)，统计河道断面的湿周χ0(i)，计算得到该水位下的断面水力半径R0(i)；
(5)计算所有水位下的水力半径R0(i)，i＝1、2、…、n；
(6)计算相临两水力半径的差值ΔR0(i)＝R0(i+1)-R0(i)，i＝1、2、…、n-1。统计ΔR0(i)＜0的结果，当ΔR0(i)＜0时，记录RR(p)＝R0(i)，Z1(p)＝Zlevel(i)，Z2(p)＝Zlevel(i+1)；当出现连续m个ΔR0(i)＜0且m≥2时，计算水力半径的变幅：记录Z1(p)＝Zlevel(i)，Z2(p)＝Zlevel(i+m)；p的取值为：p＝1、2、3、…；
(7)设ξ为河漫滩引起的该河段水力半径的最小变幅，当RR(p)≥ξ时，认为该水位处出现河漫滩，该河漫滩的高程处于Z1(p)～Z2(p)之间，河漫滩的代表高程可取二者的均值；
(8)根据(7)中得到的河漫滩高程范围Z1(p)～Z2(p)，结合断面地形容易得到该级河漫滩的宽度B(p)。
下面以以黄河花园口河道断面为例，应用本发明提出的复式断面河道河漫滩和主槽的自动识别方法，对滩槽进行自动识别，实施步骤如下：
(1)收集花园口断面的起点距和高程数据，断面形态如图2所示，从左至右对每个数据点进行编号。
(2)通过程序自动识别该断面的最低点90.36m，最高点100.02m。由最低点90.36m开始，每次增加5cm，确定计算水位。
(3)如图3所示，当水位为92.06m时，河道断面的过水区域可划分为三角形122'、梯形233'2'、梯形344'3'、梯形455'4'和三角形565'，可以分别求出每部分的面积，累加即可得到该水位下的过水面积；对于该水位下的湿周这可累加线段12'、线段2'3'、线段3'4'、线段4'5'、线段5'6的距离即可。采用该方法可以求出所有水位下的该断面过水面积、湿周，并计算对应的水力半径，结果如表1所示，图4为水力半径随水位的变化特征。
表1不同水位下断面特征参数
(4)计算相临的水位对应的水力半径的差值，当差值小于0时，记录累加连续的差值，根据河道断面的特性取ξ＝0.1，可以得到该断面存在3级河漫滩：第一阶高程为92.1m，宽434m；第二阶高程为93.3m，宽340m；第三阶高程为94.1m，宽2337m，划分形式如图5所示。
以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。

Claims (5)

1.一种复式断面河道河漫滩和主槽的自动识别方法，其特征在于，包括如下步骤：
步骤1，获取河道断面的起点距和对应的高程数据，并查找最低点高程和最高点高程，将其分别作为计算河道断面水力半径的最低水位和最高水位；
步骤2，从最低水位开始，每次抬高水位0.05m，若抬高后的水位大于或等于最高水位时，将最高水位作为最后一个水位，依次记录所有的水位数据；
步骤3，对于步骤2记录的每个水位，计算河道断面在该水位的过水面积以及河道断面的湿周，根据过水面积和湿周计算河道断面在该水位的水力半径；
步骤4，水位由低至高，计算相邻两个水位对应的水力半径的差值ΔR0(i)＝R0(i+1)-R0(i)，i＝1,2,…,n-1；当ΔR0(i)＜0时，记录RR(p)＝R0(i)，Z1(p)＝Zlevel(i)，Z2(p)＝Zlevel(i+1)；当连续出现m个ΔR0(i)＜0且m≥2时，计算水力半径的变幅：同时记录Z1(p)＝Zlevel(i)，Z2(p)＝Zlevel(i+m)；p＝p+1，p＝1,2,…；其中，R0(·)表示水力半径，Zlevel(·)表示第·个水位的水位值，n表示所有水位的个数；
步骤5，设定河漫滩引起的河段水力半径的临界最小变幅ξ，当RR(p)≥ξ时，认为RR(p)对应的水位处出现河漫滩，低于最低河漫滩的河道断面部分为主槽。
2.根据权利要求1所述复式断面河道河漫滩和主槽的自动识别方法，其特征在于，步骤2所述抬高水位的计算公式为：
Zlevel(i)＝Z min+(i-1)*0.05
其中，Zlevel(i)表示第i个水位的水位值，Zmin表示最低水位的水位值，水位单位为：m，i＝1,2,…,n，n表示所有水位的个数。
3.根据权利要求1所述复式断面河道河漫滩和主槽的自动识别方法，其特征在于，步骤3所述过水面积的计算方法为：对于任意水位，根据该水位河道断面的地形数据点，将河道断面分成若干梯形和三角形，所有梯形和三角形的面积之和即为过水面积。
4.根据权利要求1所述复式断面河道河漫滩和主槽的自动识别方法，其特征在于，步骤5所述河漫滩引起的河段水力半径的临界最小变幅ξ＝0.1。
5.根据权利要求1所述复式断面河道河漫滩和主槽的自动识别方法，其特征在于，步骤5所述河漫滩的代表高程为根据代表高程及河道断面地形得到该高程对应的河道宽度。
CN201610700632.8A 2016-08-22 2016-08-22 一种复式断面河道河漫滩和主槽的自动识别方法 Active CN106295056B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610700632.8A CN106295056B (zh) 2016-08-22 2016-08-22 一种复式断面河道河漫滩和主槽的自动识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610700632.8A CN106295056B (zh) 2016-08-22 2016-08-22 一种复式断面河道河漫滩和主槽的自动识别方法

Publications (2)

Publication Number Publication Date
CN106295056A CN106295056A (zh) 2017-01-04
CN106295056B true CN106295056B (zh) 2019-08-06

Family

ID=57661111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610700632.8A Active CN106295056B (zh) 2016-08-22 2016-08-22 一种复式断面河道河漫滩和主槽的自动识别方法

Country Status (1)

Country Link
CN (1) CN106295056B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109472109B (zh) * 2018-11-23 2019-10-01 中国水利水电科学研究院 一种用于分布式水文模型的复式河道断面概化及计算方法
CN109736259B (zh) * 2019-03-04 2019-12-31 四川大学 山区河流复式滩地利用范围及防洪安全位置的水力学确定方法
CN110837692B (zh) * 2019-10-18 2022-08-19 河海大学 一种基于挟沙力变化的河漫滩平均高程确定方法
CN111104746B (zh) * 2019-12-20 2022-09-23 河海大学 一种基于小波分析的河漫滩高程确定方法
CN113947628B (zh) * 2021-10-12 2024-06-11 杭州京杭运河二通道建设投资有限公司 一种基于不同比尺河道地形图计算冲淤量的方法
CN114997085B (zh) * 2022-08-02 2022-10-21 中国长江三峡集团有限公司 一维水动力模型构建方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103590363A (zh) * 2013-11-26 2014-02-19 武汉大学 下切性河段水位降幅的控制节点定位防护方法
CN105631168A (zh) * 2016-03-25 2016-06-01 中国水利水电科学研究院 一种实时高效的流域洪水演进可视化仿真方法
CN105844709A (zh) * 2016-03-25 2016-08-10 中国水利水电科学研究院 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8768663B2 (en) * 2010-10-08 2014-07-01 The United States Of America, As Represented By The Secretary Of The Navy Automated method and system for predicting high resolution tidal heights and currents in coastal zones

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103590363A (zh) * 2013-11-26 2014-02-19 武汉大学 下切性河段水位降幅的控制节点定位防护方法
CN105631168A (zh) * 2016-03-25 2016-06-01 中国水利水电科学研究院 一种实时高效的流域洪水演进可视化仿真方法
CN105844709A (zh) * 2016-03-25 2016-08-10 中国水利水电科学研究院 复杂河道地形流域洪水演进虚拟仿真的淹没线追踪方法

Also Published As

Publication number Publication date
CN106295056A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
CN106295056B (zh) 一种复式断面河道河漫滩和主槽的自动识别方法
Evangelista et al. Simulation of dam-break waves on movable beds using a multi-stage centered scheme
Arnold A new approach for dealing with depressions in digital elevation models when calculating flow accumulation values
CN115062389B (zh) 一种面向坝前泥沙减淤的水库闸门调度方法、装置和设备
CN107239657A (zh) 一种面向对象的水动力学建模要素管理方法
CN109271672B (zh) 一种河-湖-泵站相互影响作用下的河道水面线计算方法
CN103559368B (zh) 化学机械抛光模拟方法及其去除率计算方法
CN104933268B (zh) 一种基于一维非恒定流数值模型的洪水分析方法
CN111104746B (zh) 一种基于小波分析的河漫滩高程确定方法
CN109285219B (zh) 一种基于dem的网格型水文模型网格演算次序编码方法
CN102567594A (zh) 一种近岸岛礁型人工鱼礁群流场仿真建模的方法
CN106320255A (zh) 一种漫滩水流滩槽过流量的计算方法
CN103823936B (zh) 一种确定等流时线的方法
CN112163382B (zh) 一种平原圩区河网水动力模拟方法
Palalane et al. Analytical model of sand spit evolution
Yossef et al. Modelling large scale and long term morphological response to engineering interventions at river bifurcation
CN107704850B (zh) 一种潮沟系统地形识别处理方法
Busnelli et al. Hydro-morphological management to improve navigation and ecological functions on the ‘Canal del Dique’Colombia
Abad et al. Hydraulic model study of canoe chute and fish passage for the Chicago River north branch dam
Dissanayake et al. Process-based approach on tidal inlet evolution–Part 1
Harb et al. Numerical analysis of sediment transport processes during a flushing event of an Alpine reservoir
CN114103127B (zh) 流域地形3d打印方法及装置、电子设备、存储介质
CN114155350B (zh) 一种航道的水下地形冲淤演变分析方法及装置
CN110717247B (zh) 一种考虑梯田影响的坡面产汇流过程模拟方法
Zhang et al. Flow dynamic modeling of a large-scaled vegetated constructed wetland

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant