CN109920028A - 一种宽度平均立面二维模型的地形修正方法 - Google Patents
一种宽度平均立面二维模型的地形修正方法 Download PDFInfo
- Publication number
- CN109920028A CN109920028A CN201910182669.XA CN201910182669A CN109920028A CN 109920028 A CN109920028 A CN 109920028A CN 201910182669 A CN201910182669 A CN 201910182669A CN 109920028 A CN109920028 A CN 109920028A
- Authority
- CN
- China
- Prior art keywords
- section
- landform
- bank
- line
- ridge line
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种宽度平均立面二维模型的地形修正方法,包括:步骤1,对于拟设置网格的断面,垂直于沿河道走向做剖面,若该剖面左岸或右岸剖切到的高程最高的等高线低于模型计算时可能出现的最高水位,则进入步骤2;步骤2,以左岸山脊线与剖面的交点为左岸分界点;以右岸山脊线与剖面的交点为右岸分界点;左、右岸分界点高程以下,断面地形为剖面与地形曲面的交线;分界点上,断面地形为山脊线;步骤3,找到断面地形最低点作为划分点,其左侧为左侧地形线,右侧为右侧地形线;步骤4,计算断面垂向各网格中心高程的过流宽度。本发明还原了水流、地形和网格三者之间的关系,让立面二维模型内水流流向与地形的关系接近真实场景,模型精度高。
Description
技术领域
本发明属于工程领域,特别涉及一种宽度平均立面二维模型的地形修正方法。
背景技术
宽度平均立面二维模型是指计算网格带有宽度属性的立面二维模型,多用于水库等深窄型水体的水动力学、水质、的数值模拟计算,通常宽度平均立面二维水库模型的网格,在尺寸上具有高度和宽度两个属性。其中高度属性由用户给定,范围多为1m~5m。宽度属性则由地形计算得出。
对于常规的宽度平均立面二维水库模型,在根据地形建模时,多根据大断面数据直接换算各高程对应的网格宽度。即沿垂直于水流方向的平面剖切地形,得到与地形形态相对应的大断面数据。而后根据垂向网格划分和大断面数据,计算该断面各高程的断面宽度,作为沿垂向各网格的宽度。
在宽度平均立面二维水库模型内,水库水流沿水平方向的投影垂直于断面网格。而在库区,水流沿河道走向流动,用垂直于河道走向的平面剖切地形,得到计算网格并进行计算不会造成明显的计算误差。但是在近坝区,水流受电站进水口取水的影响,不再沿河道走向流动,三维效应明显,依然用平面来直接剖切近坝区地形,而后根据得到的断面形态确定近坝区的网格尺寸将导致较大的计算误差。
发明内容
本发明的目的在于,针对上述现有技术的不足,提供一种宽度平均立面二维模型的近坝区地形修正方法,针对近坝区水流的三维特征,根据水流流向计算各高程的过水通道宽度,并将该宽度赋值于对应高程网格,保证宽度平均立面二维模型在该处的水流过流通道接近实际情况,提高模型计算精度。
为解决上述技术问题,本发明所采用的技术方案是:
一种宽度平均立面二维模型的地形修正方法,其特点是包括以下步骤:
步骤1,对于拟设置网格的断面,在断面处垂直于沿河道走向做剖面,若该剖面左岸或右岸剖切到的高程最高的等高线低于模型计算时可能出现的最高水位,则判定该断面可能需进行地形修正,进入步骤2;
步骤2,若剖面在左、右岸剖切到的山脊线不存在,则无需进行地形修正并结束;
若剖面在左、右岸剖切到的山脊线存在,则找到剖面在左、右岸剖切到的山脊线,以左岸山脊线与剖面的交点作为左岸分界点;以右岸山脊线与剖面的交点作为右岸分界点;左、右岸分界点高程以下,断面地形为剖面与地形曲面的交线;分界点上,断面地形为山脊线;由此确定出的断面地形曲线,为一条三维曲线;
步骤3,根据步骤2确定的断面地形,找到该断面地形最低点作为左右岸的划分点;该划分点左侧为左侧地形线,右侧为右侧地形线;对于一个高于断面地形最低点的高程值z,在断面高度范围内,左侧地形线和右侧地形线各有一个与该高程值对应的节点;
步骤4,根据步骤3的结果,计算该断面垂向各网格中心高程的过流宽度,其值为该高程在左侧地形线和右侧地形线对应的节点沿断面地形在垂向投影线的长度。
与现有技术相比,本发明具有以下优点:
通过宽度平均立面二维模型内,水流流速在水平方向的投影垂直于断面的特性,提出了断面地形的修正方法,该方法可根据水流方向计算出以水流走向为法线的平面所剖切出的非平面断面,还原了水流、地形和网格三者之间的关系,让立面二维模型内水流流向与地形的关系尽量接近真实场景,有效降低了由于断面剖切带来的计算误差,提高了模型精度。
附图说明
图1是发明的方法流程示意图。
图2是断面划分示意图。
图3是地形曲线示意图。
图4是地形曲线大断面数据曲线。
图5是地形曲线垂向投影示意图。
图6是网格中心高程与网格宽度的关系曲线。
其中,
1——等高线
2——左岸山脊线
3——断面剖面
4——电站进水口
5——分界点
6——断面最低点
7——断面曲线
具体实施方式
本发明的计算流程如图1所示。白莲崖水库正常蓄水位208m,模拟时最高水位不超过正常蓄水位,模拟所用网格垂向间隔为10m。以白莲崖电站取水口上游侧断面为例,介绍本发明的实施方法,实施步骤如下:
步骤1,如图2所示,在电站上游侧做剖面。
步骤2,该水库模拟时最高水位为208m,步骤1所做剖面在右岸剖面剖切到的最高等高线高程为220m,大于208m,在左岸剖切到的最高的等高线为175m,低于208m,因此该断面左岸可能需进行地形修正,进入步骤3。
步骤3,绘制左岸地形山脊线,找到山脊线与剖面的交点,作为左岸分界点。合并分界点至208m高程的山脊线,以及分界点以右的,由步骤1剖切所得的地形线,即为该断面修正后的地形曲线(图3),该地形曲线的大断面数据见图4。
步骤4,如图5所示,找到该地形曲线最低点,将地形线分为左侧地形线和右侧地形线,在左侧和右侧地形线分别找到各网格中心高程的对应节点。
步骤5,计算各网格中心高程在左岸节点沿网格地形线沿垂向投影线至左岸节点的长度,即为该网格的过流宽度,该步骤计算结果见图6。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是局限性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护范围之内。
Claims (1)
1.一种宽度平均立面二维模型的地形修正方法,其特征在于,包括以下步骤:
步骤1,对于拟设置网格的断面,在断面处垂直于沿河道走向做剖面,若该剖面左岸或右岸剖切到的高程最高的等高线低于模型计算时可能出现的最高水位,则判定该断面可能需进行地形修正,进入步骤2;
步骤2,若剖面在左、右岸剖切到的山脊线不存在,则无需进行地形修正并结束;
若剖面在左、右岸剖切到的山脊线存在,则找到剖面在左、右岸剖切到的山脊线,以左岸山脊线与剖面的交点作为左岸分界点;以右岸山脊线与剖面的交点作为右岸分界点;左、右岸分界点高程以下,断面地形为剖面与地形曲面的交线;分界点上,断面地形为山脊线;由此确定出的断面地形曲线,为一条三维曲线;
步骤3,根据步骤2确定的断面地形,找到该断面地形最低点作为左右岸的划分点;该划分点左侧为左侧地形线,右侧为右侧地形线;对于一个高于断面地形最低点的高程值z,在断面高度范围内,左侧地形线和右侧地形线各有一个与该高程值对应的节点;
步骤4,根据步骤3的结果,计算该断面垂向各网格中心高程的过流宽度,其值为该高程在左侧地形线和右侧地形线对应的节点沿断面地形在垂向投影线的长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910182669.XA CN109920028B (zh) | 2019-03-12 | 2019-03-12 | 一种宽度平均立面二维模型的地形修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910182669.XA CN109920028B (zh) | 2019-03-12 | 2019-03-12 | 一种宽度平均立面二维模型的地形修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109920028A true CN109920028A (zh) | 2019-06-21 |
CN109920028B CN109920028B (zh) | 2022-03-11 |
Family
ID=66964279
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910182669.XA Active CN109920028B (zh) | 2019-03-12 | 2019-03-12 | 一种宽度平均立面二维模型的地形修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109920028B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489511A (zh) * | 2019-08-23 | 2019-11-22 | 西南大学 | 等高线接边高程错误修正方法、系统及电子设备和介质 |
CN111006645A (zh) * | 2019-12-23 | 2020-04-14 | 青岛黄海学院 | 基于运动与结构重建的无人机测绘方法 |
CN112507519A (zh) * | 2020-11-12 | 2021-03-16 | 中国电建集团中南勘测设计研究院有限公司 | 一种隔水幕墙全年荷载变化过程分析方法 |
CN114108659A (zh) * | 2021-11-30 | 2022-03-01 | 中交路桥建设有限公司 | 景观绿化堆坡施工方法 |
CN116226586A (zh) * | 2023-05-10 | 2023-06-06 | 长江三峡集团实业发展(北京)有限公司 | 一种模拟微囊藻群体迁移分布的理论解析方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070024623A1 (en) * | 2003-06-06 | 2007-02-01 | Peter Dugge | Method for generation of a three-dimensional terrain model |
CN105740464A (zh) * | 2016-03-03 | 2016-07-06 | 中国国土资源航空物探遥感中心 | 一种基于dem的河谷形态参数自动提取方法 |
CN108614950A (zh) * | 2018-05-18 | 2018-10-02 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于bim的河道类线性工程设计建模方法 |
CN109137816A (zh) * | 2018-08-23 | 2019-01-04 | 湖北省水利水电规划勘测设计院 | 一种基于主流摆型波震荡机制的河型划分方法 |
-
2019
- 2019-03-12 CN CN201910182669.XA patent/CN109920028B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070024623A1 (en) * | 2003-06-06 | 2007-02-01 | Peter Dugge | Method for generation of a three-dimensional terrain model |
CN105740464A (zh) * | 2016-03-03 | 2016-07-06 | 中国国土资源航空物探遥感中心 | 一种基于dem的河谷形态参数自动提取方法 |
CN108614950A (zh) * | 2018-05-18 | 2018-10-02 | 中国电建集团华东勘测设计研究院有限公司 | 一种基于bim的河道类线性工程设计建模方法 |
CN109137816A (zh) * | 2018-08-23 | 2019-01-04 | 湖北省水利水电规划勘测设计院 | 一种基于主流摆型波震荡机制的河型划分方法 |
Non-Patent Citations (2)
Title |
---|
FENG D 等: "Digital terrain model extraction in SUAS clearance survey using LiDar data", 《GEOSCIENCE & REMOTE SENSING SYMPOSIUM. IEEE》 * |
沈锐燕 等: "复杂地形河道的二维数值模拟及问题处理", 《河北工程大学学报(自然科学版)》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110489511A (zh) * | 2019-08-23 | 2019-11-22 | 西南大学 | 等高线接边高程错误修正方法、系统及电子设备和介质 |
CN110489511B (zh) * | 2019-08-23 | 2021-07-30 | 西南大学 | 等高线接边高程错误修正方法、系统及电子设备和介质 |
CN111006645A (zh) * | 2019-12-23 | 2020-04-14 | 青岛黄海学院 | 基于运动与结构重建的无人机测绘方法 |
CN112507519A (zh) * | 2020-11-12 | 2021-03-16 | 中国电建集团中南勘测设计研究院有限公司 | 一种隔水幕墙全年荷载变化过程分析方法 |
CN112507519B (zh) * | 2020-11-12 | 2024-04-16 | 中国电建集团中南勘测设计研究院有限公司 | 一种隔水幕墙全年荷载变化过程分析方法 |
CN114108659A (zh) * | 2021-11-30 | 2022-03-01 | 中交路桥建设有限公司 | 景观绿化堆坡施工方法 |
CN116226586A (zh) * | 2023-05-10 | 2023-06-06 | 长江三峡集团实业发展(北京)有限公司 | 一种模拟微囊藻群体迁移分布的理论解析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109920028B (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109920028A (zh) | 一种宽度平均立面二维模型的地形修正方法 | |
CN105931291B (zh) | 一种数字化完整牙颌建模方法 | |
CN105160658B (zh) | 一种基于子流域边界和流路特征的山脊线提取方法 | |
CN103020368B (zh) | 列车头优化设计方法 | |
CN109165476A (zh) | 一种模块化风场模型的建模方法及风场模拟方法 | |
CN106295056B (zh) | 一种复式断面河道河漫滩和主槽的自动识别方法 | |
CN109191566B (zh) | 基于tin的三维梯田建模方法 | |
CN105405164B (zh) | 一种复杂地形风电场地形改造方法 | |
CN108563872A (zh) | 网格参数化方法及基于该网格参数化方法的轴流涡轮气动优化设计方法 | |
CN112287615B (zh) | 一种复杂结构弹簧阀非定常流场多块结构网格划分方法 | |
CN106192935A (zh) | 一种鱼道进口处河道流场的调整方法和系统 | |
CN110472375A (zh) | 一种基于bim的非对称悬索桥施工管控方法 | |
CN105631101A (zh) | 一种地震作用下高耸升船机塔柱结构动态分布系数分析方法 | |
CN114059501B (zh) | 一种表深孔水舌空中碰撞后入水后流量分布的估算方法 | |
CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
CN105608287B (zh) | 高性能大流量泵站出水流道三维形体过流面设计方法 | |
CN117454724B (zh) | 一种通过cfd确定护城河式防洪堤坝几何参数的方法 | |
CN110990923A (zh) | 水利水电工程双曲拱坝的快速工程计算方法、系统及介质 | |
CN109657290A (zh) | 基于CATIA与Excel的心墙坝设计方法 | |
CN106202711A (zh) | 液力变矩器流道模型的参数化方法 | |
CN102567594A (zh) | 一种近岸岛礁型人工鱼礁群流场仿真建模的方法 | |
CN103413019A (zh) | 一种不规则轮廓双曲率外形面整体壁板离散方法 | |
CN108595857B (zh) | 铁路线路防护栅栏自动设计方法 | |
CN110223390A (zh) | 一种基于直线体素遍历原理的多段线嵌入tin算法 | |
KR101910274B1 (ko) | 지상공간 및 지하공간을 동시 고려한 침수 해석 방법 |
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 |