CN101188022A - 一种面向大城市灾害展示的洪涝淹没分析方法 - Google Patents

一种面向大城市灾害展示的洪涝淹没分析方法 Download PDF

Info

Publication number
CN101188022A
CN101188022A CNA2007103005810A CN200710300581A CN101188022A CN 101188022 A CN101188022 A CN 101188022A CN A2007103005810 A CNA2007103005810 A CN A2007103005810A CN 200710300581 A CN200710300581 A CN 200710300581A CN 101188022 A CN101188022 A CN 101188022A
Authority
CN
China
Prior art keywords
flood
seed points
grid
grid model
flooding
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
Application number
CNA2007103005810A
Other languages
English (en)
Other versions
CN101188022B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2007103005810A priority Critical patent/CN101188022B/zh
Publication of CN101188022A publication Critical patent/CN101188022A/zh
Application granted granted Critical
Publication of CN101188022B publication Critical patent/CN101188022B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Alarm Systems (AREA)

Abstract

本发明涉及一种面向大城市灾害展示中的洪涝淹没分析方法,具体步骤如下:确定待研究的淹没区域:首先,圈定一个最大的可能淹没范围,并且将河流两岸分成左右两半分别进行处理,靠河流边的边界处理为淹没区的进水边界;其次,根据国家测绘局提供的1∶1万的DEM数据,对洪水最大可能淹没范围进行剪裁,得到的区域就是所需要进行淹没分析研究的范围;生成所需的格网模型:基于格网模型的洪水淹没连通算法采用有种子点的填充算法。本发明有益的效果是:洪水淹没的特性与格网的这种大小分布特性是一致的,采用这种格网能更好地模拟洪水的淹没特性;采用三角单元格网和任意多边形格网模型能够更好地模拟这种不规则的边界。

Description

一种面向大城市灾害展示的洪涝淹没分析方法
技术领域
本发明涉及灾害信息管理领域,主要是一种面向大城市灾害展示的洪涝淹没分析方法。
背景技术
应用于大城市灾害展示中的洪涝淹没分析方法的基础建立在GIS技术之上。在给定洪水水位和洪量两种条件下,对基于三角形格网模型和任意多边形格网模型,采用平面模拟方法进行洪水淹没范围和水深分布的模拟,分别得出洪水淹没范围和深度结果。该方法通过已知的洪水水位和洪量,能够得出比较精确的洪水淹没范围和水深分布,为灾害展示提供精确的数据支持。
自然灾害在我国发生非常频繁,其中以洪涝灾害尤为严重。近十几年来,自然资源的开发利用不断扩大,城乡经济建设飞速发展,洪水出现的频率及其造成的损失也不断地增加。因此,快速、准确、科学地模拟、预测洪水淹没范围,对防洪减灾具有重要意义。特别是对于一些重点防洪城市和行蓄洪区,如果能够预先获知洪水的淹没范围和水深的分布情况,对于预先转移受灾区的生命财产,减少损失具有非常重要的价值,而且对于洪水造成的灾害损失进行评估也是非常有用的。
洪涝淹没分析方法是一种基于GIS的分析方法,目的是在提供洪水淹没范围和水深分布的服务的前提下,保持数据的准确性。其中,基于格网模型的洪水淹没算法是至关重要的。当前的洪水淹没分析方法都只进行了简单平面模拟,抛弃了复杂的地形因素的影响,输出的淹没范围和水深分布相当粗糙。本发明在分析洪水淹没范围和水深分布中考虑了复杂的地形因素,在必要的时候通过GIS技术的格网模型,从而增加了输出数据的准确性。
发明内容
有鉴于此,本发明提出了一种面向大城市灾害展示中的洪涝淹没分析方法,本方法的目的在于通过现有的洪水水位和洪量,结合地形因素,提供准确的洪水淹没范围和水深分布。
本发明解决其技术问题所采用的技术方案:这种面向大城市灾害展示中的洪涝淹没分析
方法,具体步骤如下:
1.确定待研究的淹没区域:
首先,圈定一个最大的可能淹没范围,并且将河流两岸分成左右两半分别进行处理,靠河流边的边界处理为淹没区的进水边界。这样处理对于防洪减灾来说是合理的,一般在防洪区域,河流两岸堤防建设的洪水保证率是不一样的,有重点地保护一些地区和放弃一些地区,所以需要将两岸分开处理。
其次,根据国家测绘局提供的1∶1万的DEM数据(地形数据),对洪水最大可能淹没范围进行剪裁,得到的区域就是所需要进行淹没分析研究的范围。
2.生成所需的格网模型:
(1)三角单元格网模型。将研究区的DEM转换为TIN模型,提取三角单元格网,并对每个三角网格赋高程值,高程值按3个顶点从GRID上取得的高程值取平均求得。生成的三角格网就是要进行洪水淹没分析的三角单元格网模型,如图1所示。
(2)任意多边形格网模型。将研究区的DEM转换为几何特征图层,在处理时自动将具有相同高程并且相邻的单元合并为一个多边形,多边形的高程自动取GRID的高程值,这样保证DEM的高程数据原始精度完全不损失。生成的多边形格网就是要进行洪水淹没分析的任意多边形格网模型,如图2所示。
3.基于格网模型的洪水淹没连通算法:
对洪水淹没区域连通性的考虑,在一些淹没分析软件中,仅考虑高程平铺的问题,即在任何地势低洼的区域都同时进水,实际上从洪水本身淹没的角度来说这是不准确的,洪水首先是从洪水源处开始向外扩散淹没,只有水位高程达到一定程度之后,洪水才能越过某一地势较高的区域到达另一个洼地。洪水淹没的连通性算法可以采用有种子点的填充算法。
种子点的填充算法步骤:
我们总是先给定区域内一点,称为种子点。对种子点赋予指定的颜色,并将这种颜色扩充到区域内的所有点。区域填充的扫描线算法相对于前面的递归算法是一种较好的区域填充扫描线算法。它很好地避免了递归算法由于多层递归、反复进出堆栈操作而造成的费时费内存的缺点。这一借助于堆栈技术实现的种子点入堆栈的区域填充扫描线算法可以由如下几个步骤实现:
步骤1:(初始化)设置一空的堆栈,将预先给定的种子点(x,y)压入堆栈。
步骤2:(种子点出栈)弹出栈顶元素(x,y)作为当前种子点。
步骤3:(向左填充)沿当前种子点所在的扫描线从种子点开始向左侧逐个像素点进行填色,直到区域边界,像素点颜色值置为new color,填充区段左端点记为xl。
步骤4:(向右填充)类似步骤3的操作向右填充,相应地获得一右端点值记为xr。
步骤5:(新种子点入栈)分别确定与当前扫描线相邻的上下两条扫描线上位于给定的区段[xl;xr]内含于区域的区段。如果有这样的区段,则取区段的右端点为新的种子点压入堆栈(这样的区段可不止一个)。如果区段内的像素点值为newcolor或者边界像素点颜色值为boundary color,则不产生新的种子点入栈。
步骤6:(结束判断)如果种子了点堆栈非空则转步骤2;否则算法结束。
本发明有益的效果是:
(1)洪水淹没的特性与格网的这种大小分布特性是一致的,即在平坦的地区淹没面积大,在陡峭的区域淹没面积小,所以采用这种格网能更好地模拟洪水的淹没特性;
(2)洪水的淹没边界和河流边界都是非常不规则的,采用三角单元格网和任意多边形格网模型比规则的四边形格网模型,能够更好地模拟这种不规则的边界。
(3)网格大小疏密变化不一致,既能满足模型物理意义上的需求,也能节省计算机的存储空间,提高计算速度。
附图说明
图1为三角形格网模型图。
图2为任意多边形格网模型图。
图3为种子点填充算法的流程图。
具体实施方式
下面结合附图和实施例对本发明作进一步介绍:
图1所示为生成的三角格网,也就是要进行洪水淹没分析的三角单元格网模型。
图2所示为生成的三角格网,也就是要进行洪水淹没分析的三角单元格网模型。
图3所示为种子点填充算法流程。流程共分六个步骤:
步骤1,初始化一个空堆栈,将给定的种子点(x,y)压入堆栈。
步骤2,从堆栈中弹出一个种子,作为当前的种子点。
步骤3,沿当前种子点所在的扫描线从种子点开始向左侧逐个像素点进行填色,直到区域边界,像素点颜色值置为new color,填充区段左端点记为xl。
步骤4,类似步骤3的操作向右填充,相应地获得一右端点值记为xr。
步骤5,分别确定与当前扫描线相邻的上下两条扫描线上位于给定的区段[xl;xr]内含于区域的区段。如果有这样的区段,则取区段的右端点为新的种子点压入堆栈(这样的区段可不止一个)。如果区段内的像素点值为new color或者边界像素点颜色值为boundary color,则不产生新的种子点入栈。
步骤6,如果种子了点堆栈非空则转步骤2;否则算法结束。
上述实施例用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。

Claims (2)

1.一种面向大城市灾害展示中的洪涝淹没分析方法,其特征在于:步骤如下:
(1)、确定待研究的淹没区域:首先,圈定一个最大的可能淹没范围,并且将河流两岸分成左右两半分别进行处理,靠河流边的边界处理为淹没区的进水边界;其次,根据国家测绘局提供的1∶1万的DEM数据,对洪水最大可能淹没范围进行剪裁,得到的区域就是所需要进行淹没分析研究的范围;
(2)、生成所需的格网模型:
(2.1)、三角单元格网模型,将研究区的DEM转换为TIN模型,提取三角单元格网,并对每个三角网格赋高程值,高程值按3个顶点从GRID上取得的高程值取平均求得,生成的三角格网就是要进行洪水淹没分析的三角单元格网模型;
(2.2)、任意多边形格网模型,将研究区的DEM转换为几何特征图层,在处理时自动将具有相同高程并且相邻的单元合并为一个多边形,多边形的高程自动取GRID的高程值,生成的多边形格网就是要进行洪水淹没分析的任意多边形格网模型,;
(3)、基于格网模型的洪水淹没连通算法:对洪水淹没区域连通性,洪水淹没的连通性算法采用有种子点的填充算法。
2.根据权利要求1所述的面向大城市灾害展示中的洪涝淹没分析方法,其特征在于:种子点的填充算法步骤:先给定区域内一点,称为种子点,对种子点赋予指定的颜色,并将这种颜色扩充到区域内的所有点,区域填充的扫描线算法由如下几个步骤实现:
步骤1:设置一空的堆栈,将预先给定的种子点(x,y)压入堆栈;
步骤2:弹出栈顶元素(x,y)作为当前种子点;
步骤3:沿当前种子点所在的扫描线从种子点开始向左侧逐个像素点进行填色,直到区域边界,像素点颜色值置为new color,填充区段左端点记为xl;
步骤4:类似步骤3的操作向右填充,相应地获得一右端点值记为xr;
步骤5:分别确定与当前扫描线相邻的上下两条扫描线上位于给定的区段[xl;xr]内含于区域的区段,如果有这样的区段,则取区段的右端点为新的种子点压入堆栈,如果区段内的像素点值为new color或者边界像素点颜色值为boundary color,则不产生新的种子点入栈;
步骤6:如果种子了点堆栈非空则转步骤2;否则算法结束。
CN2007103005810A 2007-12-20 2007-12-20 一种面向大城市灾害展示的洪涝淹没分析方法 Expired - Fee Related CN101188022B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2007103005810A CN101188022B (zh) 2007-12-20 2007-12-20 一种面向大城市灾害展示的洪涝淹没分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2007103005810A CN101188022B (zh) 2007-12-20 2007-12-20 一种面向大城市灾害展示的洪涝淹没分析方法

Publications (2)

Publication Number Publication Date
CN101188022A true CN101188022A (zh) 2008-05-28
CN101188022B CN101188022B (zh) 2010-09-29

Family

ID=39480395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2007103005810A Expired - Fee Related CN101188022B (zh) 2007-12-20 2007-12-20 一种面向大城市灾害展示的洪涝淹没分析方法

Country Status (1)

Country Link
CN (1) CN101188022B (zh)

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101794495B (zh) * 2010-02-05 2011-08-24 国网电力科学研究院 洪水预报系统中实时校正模型的优选方法
CN102902893A (zh) * 2012-10-29 2013-01-30 南京信息工程大学 一种基于dem的汇水区降雨积水深度的计算方法
CN102968806A (zh) * 2012-10-26 2013-03-13 中国石油大学(华东) 一种不规则边界区域内自适应网格的生成方法
CN103116695A (zh) * 2013-01-16 2013-05-22 吴立新 一种基于cd-tin的城区内涝淹没过程模拟方法
CN103256914A (zh) * 2012-12-27 2013-08-21 北京地拓科技发展有限公司 一种基于dem计算淤地坝淹没面积的方法及系统
CN104537718A (zh) * 2014-12-31 2015-04-22 河南理工大学 基于栅格dem的大河流域洼地填充预处理方法
CN106022953A (zh) * 2016-05-13 2016-10-12 国网浙江省电力公司台州供电公司 电网设施暴雨风险评估方法
CN106845074A (zh) * 2016-12-19 2017-06-13 中国人民解放军信息工程大学 建立六边形格网模型的方法、洪水推演模拟方法及其系统
CN110348324A (zh) * 2019-06-20 2019-10-18 武汉大学 一种基于遥感大数据的洪涝实时淹没分析方法及系统
CN110956412A (zh) * 2019-12-16 2020-04-03 珠江水利委员会珠江水利科学研究院 基于实景模型的洪灾动态评估方法、装置、介质和设备
CN111260714A (zh) * 2020-01-17 2020-06-09 成都理工大学 一种洪涝受灾评估方法、装置、设备及计算机存储介质
CN111581756A (zh) * 2020-04-28 2020-08-25 国网通用航空有限公司 洪水淹没范围的确定方法及装置
CN111612908A (zh) * 2020-05-26 2020-09-01 江西省水利规划设计研究院 水利工程淹没范围展示方法、装置、设备和存储介质
CN111652777A (zh) * 2020-05-15 2020-09-11 长江勘测规划设计研究有限责任公司 洪灾应急避险方法
CN113884051A (zh) * 2021-09-24 2022-01-04 重庆市勘测院 用于获取建筑物被淹没面积的方法及装置、电子设备

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1246803C (zh) * 2003-11-24 2006-03-22 浙江大学 一种与视点相关的用于地形网格生成的误差度量方法
CN101009024A (zh) * 2007-02-09 2007-08-01 肖盛燮 一种滑坡灾害可视化的实现方法

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101794495B (zh) * 2010-02-05 2011-08-24 国网电力科学研究院 洪水预报系统中实时校正模型的优选方法
CN102968806A (zh) * 2012-10-26 2013-03-13 中国石油大学(华东) 一种不规则边界区域内自适应网格的生成方法
CN102902893A (zh) * 2012-10-29 2013-01-30 南京信息工程大学 一种基于dem的汇水区降雨积水深度的计算方法
CN102902893B (zh) * 2012-10-29 2015-03-11 南京信息工程大学 一种基于dem的汇水区降雨积水深度的计算方法
CN103256914B (zh) * 2012-12-27 2015-11-18 北京地拓科技发展有限公司 一种基于dem计算淤地坝淹没面积的方法及系统
CN103256914A (zh) * 2012-12-27 2013-08-21 北京地拓科技发展有限公司 一种基于dem计算淤地坝淹没面积的方法及系统
CN103116695A (zh) * 2013-01-16 2013-05-22 吴立新 一种基于cd-tin的城区内涝淹没过程模拟方法
CN104537718B (zh) * 2014-12-31 2017-06-16 河南理工大学 基于栅格dem的大河流域洼地填充预处理方法
CN104537718A (zh) * 2014-12-31 2015-04-22 河南理工大学 基于栅格dem的大河流域洼地填充预处理方法
CN106022953A (zh) * 2016-05-13 2016-10-12 国网浙江省电力公司台州供电公司 电网设施暴雨风险评估方法
CN106845074A (zh) * 2016-12-19 2017-06-13 中国人民解放军信息工程大学 建立六边形格网模型的方法、洪水推演模拟方法及其系统
CN106845074B (zh) * 2016-12-19 2019-03-19 中国人民解放军信息工程大学 建立六边形格网模型的方法、洪水推演模拟方法及其系统
CN110348324A (zh) * 2019-06-20 2019-10-18 武汉大学 一种基于遥感大数据的洪涝实时淹没分析方法及系统
CN110348324B (zh) * 2019-06-20 2022-06-14 武汉大学 一种基于遥感大数据的洪涝实时淹没分析方法及系统
CN110956412B (zh) * 2019-12-16 2020-11-13 珠江水利委员会珠江水利科学研究院 基于实景模型的洪灾动态评估方法、装置、介质和设备
CN110956412A (zh) * 2019-12-16 2020-04-03 珠江水利委员会珠江水利科学研究院 基于实景模型的洪灾动态评估方法、装置、介质和设备
CN111260714A (zh) * 2020-01-17 2020-06-09 成都理工大学 一种洪涝受灾评估方法、装置、设备及计算机存储介质
CN111260714B (zh) * 2020-01-17 2023-07-04 成都理工大学 一种洪涝受灾评估方法、装置、设备及计算机存储介质
CN111581756A (zh) * 2020-04-28 2020-08-25 国网通用航空有限公司 洪水淹没范围的确定方法及装置
CN111652777A (zh) * 2020-05-15 2020-09-11 长江勘测规划设计研究有限责任公司 洪灾应急避险方法
CN111612908A (zh) * 2020-05-26 2020-09-01 江西省水利规划设计研究院 水利工程淹没范围展示方法、装置、设备和存储介质
CN113884051A (zh) * 2021-09-24 2022-01-04 重庆市勘测院 用于获取建筑物被淹没面积的方法及装置、电子设备
CN113884051B (zh) * 2021-09-24 2023-12-05 重庆市勘测院 用于获取建筑物被淹没面积的方法及装置、电子设备

Also Published As

Publication number Publication date
CN101188022B (zh) 2010-09-29

Similar Documents

Publication Publication Date Title
CN101188022B (zh) 一种面向大城市灾害展示的洪涝淹没分析方法
CN102612682B (zh) 用于储层建模和模拟的方法和设备
CN105184867B (zh) 基于规则格网和角点网格技术的三维地层建模方法
CN103454686B (zh) 基于地层切片的小尺度沉积相进行储层预测的方法及系统
Steffler et al. Two-dimensional depth averaged model of river hydrodynamics and fish habitat
CN102156779A (zh) 地下水流仿真与预测分析方法
CN103279986A (zh) 三维水平地质剖面图制作方法及其用途
CN105093277B (zh) 地震建模中的浅中深层速度融合方法
Bonacci et al. Karst rivers’ particularity: an example from Dinaric karst (Croatia/Bosnia and Herzegovina)
CN102254349A (zh) 一种使用钻孔数据构建沉积地层系统三维实体模型的方法
CN102129712A (zh) 基于多地层及三维的土石方数量的三角网模型构建方法
CN102609982B (zh) 空间地质数据非结构化模式的拓扑发现方法
CN104573705A (zh) 一种建筑物激光扫描点云数据的聚类方法
CN102982566B (zh) 一种基于最小曲率法插值的含断层的等值线图自动生成方法
CN102819631B (zh) 一种三维地质勘察分析和交互方法
CN106355640A (zh) 一种地上地下三维一体化空间的处理方法和系统
CN108053458A (zh) 一种基于地理信息系统的流域概化图制作及展示方法
Pasternack et al. Hierarchically nested river landform sequences. Part 2: Bankfull channel morphodynamics governed by valley nesting structure
CN111754623A (zh) 一种基于多点地质统计学的离散裂缝建模方法
CN108957535A (zh) 油藏地震储层的预测方法及装置
Cornette et al. Refining discrete fracture networks with surface microseismic mechanism inversion and mechanism-driven event location
CN102968470B (zh) 一种准确计算具有水涮窝的冲沟沟头沟壁形态变化的方法
Jardine et al. Carbonate reservoir description
ZHU et al. An approach to computer modeling of geological faults in 3D and an application
CN102467754B (zh) 煤矿地质三维模型的构建方法和系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20080528

Assignee: Enjoyor Co., Ltd.

Assignor: Zhejiang University

Contract record no.: 2013330000105

Denomination of invention: A flood submerging analysis method oriented to a large city disaster demonstration

Granted publication date: 20100929

License type: Common License

Record date: 20130425

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100929

Termination date: 20171220