CN113421341A - 一种水库冲刷淤积计算方法 - Google Patents
一种水库冲刷淤积计算方法 Download PDFInfo
- Publication number
- CN113421341A CN113421341A CN202110711302.XA CN202110711302A CN113421341A CN 113421341 A CN113421341 A CN 113421341A CN 202110711302 A CN202110711302 A CN 202110711302A CN 113421341 A CN113421341 A CN 113421341A
- Authority
- CN
- China
- Prior art keywords
- volume
- unit
- reservoir
- calculating
- scouring
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computer Graphics (AREA)
- Operations Research (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Evolutionary Computation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种水库冲刷淤积计算方法,属于工程测量数据应用技术领域,具体是一种水库冲刷淤积计算方法。本发明的方法,采用多期实测水库地形数据成果创建规则格网DEM数据,然后对DEM格网单元进行三角剖分并利用不同类别的三角柱模型计算各体积单元的冲刷体积和淤积体积,再进行累加得到总体的冲刷体积和淤积体积,实现了水库冲刷淤积的准确计算。本发明的方法设计合理,具有自动化程度高、效率高、准确性好、成本低等特点。
Description
技术领域
本发明属于工程测量数据应用技术领域,具体是一种水库冲刷淤积计算方法。
背景技术
在水库运行中,水库的冲刷淤积计算对于了解水库的运行状态具有重要作用,水库运维管理人员可依据该信息把握水库运行的总体情况,必要时可采取库区清淤工作以保证水库的运行安全。
目前对于水库冲刷淤积的计算方式有多种,主要包括规则网格镶嵌法、断面法、DTM法、三角柱法等。其中,三角柱法是较为精确的一种计算模型,但目前的三角柱计算方法中,只考虑了常规形式下的三角柱模型,没有顾及到非常规的情况,另外常用的基于不规则三角网的DEM构建方式不能保证两期DEM数据节点互相对应,给三角柱法的应用带来不便,因此存在一定不足。
发明内容
本发明的目的是解决目前水库冲刷淤积计算模型中存在的不足的问题,提出一种冲刷淤积计算方法。
本发明采用如下技术方案实现。
一种水库冲刷淤积计算方法,所述的水库冲刷淤积计算方法包括以下步骤:
步骤1,利用不同期的实测水库地形数据成果,创建至少两期规则格网形式的数字高程模型:DEM;
步骤2,将至少两期规则格网DEM的每个格网单元以一条对角线进行划分,得到冲刷淤积计算基本单元;所述的冲刷淤积计算基本单元即体积单元;
步骤3,对每个体积单元计算冲刷体积Vei和淤积体积Vsi。
本发明所述的水库冲刷淤积计算方法包括步骤4,将每个体积单元的冲刷体积和淤积体积累加,得到两期数据总共的冲刷体积Ve和淤积体积Vs。
本发明所述的步骤1具体为:(1)收集同一个水库的至少两期实测地形数据成果,检查数据的完整性和有效性;(2)设定规则格网的起止坐标及格网大小d;(3)利用多期实测地形数据,按照设定好的规则格网,分别创建多期规则格网DEM数据。
本发明所述的步骤2具体为:将两期规则格网DEM数据的各个格网单元按照对角线进行剖分,则每一个格网单元可以划分为两个体积单元,构建体积单元集合V。
本发明所述的步骤3具体为:(1)从体积单元集合V中任取一个体积单元Vi,(2)根据各个节点两期高程的大小判断该体积单元属于哪一个类别;(3)根据各类别的体积计算公式分别计算该体积单元中的冲刷体积Vei和淤积体积Vsi;(4)重复上述过程,知道所有的体积单元遍历完毕,则得到了所有体积单元的冲刷体积和淤积体积。
本发明所述的步骤4具体为:将每个体积单元的冲刷体积Vei和淤积体积Vsi进行累加,得到两期数据总共的冲刷体积Ve和淤积体积Vs。
本发明所述的步骤2中,冲刷淤积计算基本单元包含类别1:纯冲刷,具体为,体积单元中后一期节点的高程全部不大于前一期节点高程,则该体积单元内全部为冲刷;冲刷体积为三角柱体A1B1C1ABC的体积,设点A、B、C的平面坐标分别为(x1,y1)、(x2,y2)、(x3,y3),h1、h2、h3分别为AA1、BB1、CC1的高程差,则体积计算公式如下:
本发明所述的步骤2中,冲刷淤积计算基本单元包含类别2:纯淤积,具体为,体积单元中后一期节点的高程全部不小于前一期节点高程,则该体积单元内全部为淤积。淤积体积为三角柱体A1B1C1ABC的体积,设点A、B、C的平面坐标分别为(x1,y1)、(x2,y2)、(x3,y3),h1、h2、h3分别为AA1、BB1、CC1的高程差,则体积计算公式如下:
本发明所述的步骤2中,冲刷淤积计算基本单元包含类别3:半冲刷,具体为,体积单元中后一期节点有且只有一个点高程大于前一期对应节点高程,则该体积单元内部分为冲刷。冲刷部分为五面体EFAA1B1B,其体积等于四棱锥C1AA1B1B的体积减去四棱锥C1EFAB的体积,记四边形AA1B1B的面积为S1,C1到平面AA1B1B的距离为h1,四边形EFAB的面积为S2,C1到平面EFAB的距离为h2;淤积部分为三棱锥CC1EF,记三角形C1EF的面积为S3,C到平面C1EF的距离为h3,则体积计算公式如下:
本发明所述的步骤2中,冲刷淤积计算基本单元包含类别4:半淤积,具体为,体积单元中后一期节点有且只有一个点高程小于前一期对应节点高程,则该体积单元内部分为淤积。淤积部分为五面体EFAA1B1B,其体积等于四棱锥C1AA1B1B的体积减去四棱锥C1EFAB的体积,记四边形AA1B1B的面积为S1,C1到平面AA1B1B的距离为h1,四边形EFAB的面积为S2,C1到平面EFAB的距离为h2;冲刷部分为三棱锥CC1EF,记三角形C1EF的面积为S3,C到平面C1EF的距离为h3,则体积计算公式如下:
本发明设计合理,具有自动化程度高、效率高、准确性好、成本低等特点。具体体现为:
1、本发明利用实测地形数据创建规则格网DEM,能够较好的保持地形的完整性,能够比较真实的反应实际地形变化,根据两期DEM节点高程的大小关系运用不同类别的体积计算公式,较好的反应了地形的变化情况,能够得到更为精确的冲刷淤积计算结果。
2、本发明采用的冲刷淤积计算模型能够较好的反应出实际地形的变化情况,计算结果准确。
3、本发明自动化程度高,冲刷淤积计算无需人工干预,可以直接得到计算结果。
附图说明
图1是本发明的总体处理流程图。
图2是步骤1的处理流程图。
图3是步骤2的示意图。
图4是步骤3的处理流程图。
图5是不同类别体积计算单元的示意图。
图6是两期规则格网DEM(部分)示例图。
图7是两期规则格网DEM三角剖分示例图。
图8是体积单元内计算冲刷淤积体积(灰色为冲刷,黑色为淤积)示例图。
图9是体积累加得到总共的冲刷体积和淤积体积示例图。
具体实施方式
见图1,一种水库冲刷淤积计算方法,包括以下步骤:
步骤1,利用不同期的实测水库地形数据成果,创建两期规则格网形式的数字高程模型(以下简称DEM);具体步骤为(如图2):
(1)收集同一个水库的至少两期实测地形数据成果,检查数据的完整性和有效性;
(2)设定规则格网的起止坐标及格网大小d;
(3)利用多期实测地形数据,按照设定好的规则格网,分别创建多期规则格网DEM数据。
步骤2,将两期规则格网DEM的每个格网单元以一条对角线进行划分,得到冲刷淤积计算基本单元(以下简称体积单元,如图3);具体处理过程为:
将两期规则格网DEM数据的各个格网单元按照对角线进行剖分,则每一个格网单元可以划分为两个体积单元,构建体积单元集合V。
步骤3,对每个体积单元计算冲刷体积Vei和淤积体积Vsi;具体过程为(如图4):
(1)从体积单元集合V中任取一个体积单元Vi,
(2)根据各个节点两期高程的大小判断该体积单元属于哪一个类别;
(3)根据各类别的体积计算公式分别计算该体积单元中的冲刷体积Vei和淤积体积Vsi;
(4)重复上述过程,知道所有的体积单元遍历完毕,则得到了所有体积单元的冲刷体积和淤积体积。
步骤4,将每个体积单元的冲刷体积和淤积体积累加,得到两期数据总共的冲刷体积Ve和淤积体积Vs。具体处理过程为:
将每个体积单元的冲刷体积Vei和淤积体积Vsi进行累加,得到两期数据总共的冲刷体积Ve和淤积体积Vs。
本发明的具体过程中涉及到的体积计算基本单元类别及相应的体积计算公式如下:
类别1:纯冲刷(图5中类别1,A、B、C表示后一期节点,A1、B1、C1表示前一期节点)。体积单元中后一期节点的高程全部不大于前一期节点高程,则该体积单元内全部为冲刷。冲刷体积为三角柱体A1B1C1ABC的体积,设点A、B、C的平面坐标分别为(x1,y1)、(x2,y2)、(x3,y3),h1、h2、h3分别为AA1、BB1、CC1的高程差,则体积计算公式如下:
类别2:纯淤积(图5中类别2,A、B、C表示后一期节点,A1、B1、C1表示前一期节点)。体积单元中后一期节点的高程全部不小于前一期节点高程,则该体积单元内全部为淤积。淤积体积为三角柱体A1B1C1ABC的体积,设点A、B、C的平面坐标分别为(x1,y1)、(x2,y2)、(x3,y3),h1、h2、h3分别为AA1、BB1、CC1的高程差,则体积计算公式如下:
类别3:半冲刷(图5中类别3,A、B、C表示后一期节点,A1、B1、C1表示前一期节点)。体积单元中后一期节点有且只有一个点高程大于前一期对应节点高程,则该体积单元内部分为冲刷。冲刷部分为五面体EFAA1B1B,其体积等于四棱锥C1AA1B1B的体积减去四棱锥C1EFAB的体积,记四边形AA1B1B的面积为S1,C1到平面AA1B1B的距离为h1,四边形EFAB的面积为S2,C1到平面EFAB的距离为h2;淤积部分为三棱锥CC1EF,记三角形C1EF的面积为S3,C到平面C1EF的距离为h3,则体积计算公式如下:
类别4:半淤积(图5中类别4,A、B、C表示后一期节点,A1、B1、C1表示前一期节点)。体积单元中后一期节点有且只有一个点高程小于前一期对应节点高程,则该体积单元内部分为淤积。淤积部分为五面体EFAA1B1B,其体积等于四棱锥C1AA1B1B的体积减去四棱锥C1EFAB的体积,记四边形AA1B1B的面积为S1,C1到平面AA1B1B的距离为h1,四边形EFAB的面积为S2,C1到平面EFAB的距离为h2;冲刷部分为三棱锥CC1EF,记三角形C1EF的面积为S3,C到平面C1EF的距离为h3,则体积计算公式如下:
如图6是由某水库的两期实测地形数据成果根据步骤1所创建的两期规则格网DEM数据(部分),格网大小设置为1m。
如图7是对两期DEM数据的每个格网单元以一条对角线进行划分的示意图,得到了体积单元集合。
如图8是对各个体积单元内进行冲刷淤积计算的示意,通过步骤3可以得到每个体积单元内的冲刷体积Vei和淤积体积Vsi。
如图9是各个体积单元的冲刷体积Vei和淤积体积Vsi进行累加得到的最终结果,得到了两期数据总共的冲刷体积Ve为5211539.963m3和淤积体积Vs为8538880.850m3,可以得出该水库在两次测量期间总体上呈淤积状态,总体淤积体积为3327340.888m3。
Claims (10)
1.一种水库冲刷淤积计算方法,其特征在于,所述的水库冲刷淤积计算方法包括以下步骤:
步骤1,利用不同期的实测水库地形数据成果,创建至少两期规则格网形式的数字高程模型:DEM;
步骤2,将至少两期规则格网DEM的每个格网单元以一条对角线进行划分,得到冲刷淤积计算基本单元;所述的冲刷淤积计算基本单元即体积单元;
步骤3,对每个体积单元计算冲刷体积Vei和淤积体积Vsi。
2.根据权利要求1所述的水库冲刷淤积计算方法,其特征在于,所述的水库冲刷淤积计算方法包括步骤4,将每个体积单元的冲刷体积和淤积体积累加,得到两期数据总共的冲刷体积Ve和淤积体积Vs。
3.根据权利要求1所述的水库冲刷淤积计算方法,其特征在于,所述的步骤1具体为:(1)收集同一个水库的至少两期实测地形数据成果;(2)设定规则格网的起止坐标及格网大小d;(3)利用多期实测地形数据,按照设定好的规则格网,分别创建多期规则格网DEM数据。
4.根据权利要求1所述的水库冲刷淤积计算方法,其特征在于,所述的步骤2具体为:将两期规则格网DEM数据的各个格网单元按照对角线进行剖分,则每一个格网单元可以划分为两个体积单元,构建体积单元集合V。
5.根据权利要求1所述的水库冲刷淤积计算方法,其特征在于,所述的步骤3具体为:(1)从体积单元集合V中任取一个体积单元Vi,(2)根据各个节点两期高程的大小判断该体积单元属于哪一个类别;(3)根据各类别的体积计算公式分别计算该体积单元中的冲刷体积Vei和淤积体积Vsi;(4)重复上述过程,知道所有的体积单元遍历完毕,则得到了所有体积单元的冲刷体积和淤积体积。
6.根据权利要求1所述的水库冲刷淤积计算方法,其特征在于,所述的步骤4具体为:将每个体积单元的冲刷体积Vei和淤积体积Vsi进行累加,得到两期数据总共的冲刷体积Ve和淤积体积Vs。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110711302.XA CN113421341B (zh) | 2021-06-25 | 2021-06-25 | 一种水库冲刷淤积计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110711302.XA CN113421341B (zh) | 2021-06-25 | 2021-06-25 | 一种水库冲刷淤积计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113421341A true CN113421341A (zh) | 2021-09-21 |
CN113421341B CN113421341B (zh) | 2023-04-21 |
Family
ID=77716715
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110711302.XA Active CN113421341B (zh) | 2021-06-25 | 2021-06-25 | 一种水库冲刷淤积计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113421341B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114595596A (zh) * | 2022-05-07 | 2022-06-07 | 长江水利委员会水文局 | 一种水库淤积物密实沉降速率的模拟方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278115A (zh) * | 2012-12-27 | 2013-09-04 | 北京地拓科技发展有限公司 | 一种基于dem计算淤地坝淤积量的方法及系统 |
CN106446910A (zh) * | 2016-09-12 | 2017-02-22 | 电子科技大学 | 一种复杂地质曲面特征提取与重构方法 |
-
2021
- 2021-06-25 CN CN202110711302.XA patent/CN113421341B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278115A (zh) * | 2012-12-27 | 2013-09-04 | 北京地拓科技发展有限公司 | 一种基于dem计算淤地坝淤积量的方法及系统 |
CN106446910A (zh) * | 2016-09-12 | 2017-02-22 | 电子科技大学 | 一种复杂地质曲面特征提取与重构方法 |
Non-Patent Citations (3)
Title |
---|
燕荔: "基于GIS技术暨泥沙水下休止角的水库淤积量测算方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
袁伟: "DEM在水库冲刷淤积分析中的应用", 《云南水力发电》 * |
赵吉祥: "数字河床洲滩演变空间建模分析", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114595596A (zh) * | 2022-05-07 | 2022-06-07 | 长江水利委员会水文局 | 一种水库淤积物密实沉降速率的模拟方法 |
CN114595596B (zh) * | 2022-05-07 | 2022-08-30 | 长江水利委员会水文局 | 一种水库淤积物密实沉降速率的模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113421341B (zh) | 2023-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110298076B (zh) | 一种基于gis和swmm的城市内涝智能建模及分析方法 | |
US10147057B2 (en) | Systems and methods for flood zone modeling | |
CN107239657B (zh) | 一种面向对象的水动力学建模要素管理方法 | |
CN111369059A (zh) | 基于内涝快速模拟耦合模型的城市内涝预测方法及系统 | |
CN110232737B (zh) | 一种城市汇水区划分方法 | |
CN114218840B (zh) | 河口航道水沙运动及其地形演变整体建模及可视化系统 | |
CN115167212B (zh) | 基于监测平台的基坑动态施工控制系统及方法 | |
CN114647881B (zh) | 顾及建筑物微观水文过程的城市内涝建模方法 | |
CN113807008B (zh) | 一种基于深度学习的城市暴雨内涝模拟方法 | |
CN111046567B (zh) | 一种基于Godunov格式的城市排水管网水流数值模拟方法 | |
CN113421341A (zh) | 一种水库冲刷淤积计算方法 | |
CN112417562B (zh) | 一种基于Dynamo的土方施工建模方法和装置 | |
CN109271662A (zh) | 一种基于实时监测数据处理的桥梁桩基健康状态评估方法 | |
CN106126911B (zh) | 一种考虑潮河水位顶托影响的城市内涝风险评估方法 | |
JP6103083B1 (ja) | 情報処理装置、情報処理システム、情報処理方法及びプログラム | |
CN115062389A (zh) | 一种面向坝前泥沙减淤的水库闸门调度方法、装置和设备 | |
CN112268595A (zh) | 一种区域平均地下水水位监测方法 | |
CN114648617A (zh) | 一种基于数字高程模型dem的水系提取方法 | |
CN113627817A (zh) | 一种基于城市内涝条件下道路连通性分析方法 | |
CN116361973B (zh) | 考虑排水管网节点水位自动修正的城市洪涝过程模拟方法 | |
CN111581603A (zh) | 载荷谱处理方法及装置 | |
CN107180287B (zh) | 一种基于汶川强震区的泥石流物源转换为泥石流的转换率计算方法 | |
CN113947628B (zh) | 一种基于不同比尺河道地形图计算冲淤量的方法 | |
CN115828786A (zh) | 一种次洪过程中淤地坝系减沙量的计算方法 | |
CN113609683B (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 |