GIS洪水淹没模拟报告?水文分析步骤是?

GIS基础理论
Dr.GIS
wowwwai GIS研习社 · 工具流程与项目排障

很多同学在做课程设计、项目汇报或应急分析时,会直接问:GIS洪水淹没模拟报告?水文分析步骤是? 这类报告不是简单把一张淹没范围图贴上去,而是要说明数据来源、地形处理、汇流分析、水位或流量设定、淹没范围提取、结果核查和不确定性。

GIS洪水淹没模拟报告 水文分析步骤流程图
GIS洪水淹没模拟报告通常需要从DEM处理、水文分析、淹没范围计算到成果核查形成完整流程。

引言:GIS洪水淹没模拟报告要解决什么问题

GIS洪水淹没模拟报告的核心目标,是回答某一水位、降雨情景或洪水流量条件下,哪些区域可能被淹没、淹没深度大约是多少、受影响的道路、居民地、耕地或设施有哪些。

对GIS学习者来说,常见难点不在于“画一张水面覆盖图”,而在于报告是否能说明清楚以下问题:

  • 使用了什么地形数据,例如DEM、等高线或实测断面。
  • 采用的是静态水位淹没法,还是结合水文水动力模型的结果。
  • DEM是否经过填洼、裁剪、坐标统一和高程基准检查。
  • 淹没范围是否与河道、堤防、道路、建筑物等现实地物相符。
  • 输出成果是否包含淹没范围、淹没深度、影响统计和风险说明。

本文按一份实用报告的写法,梳理GIS洪水淹没模拟的水文分析步骤,适合使用ArcGIS Pro、QGIS、HEC-RAS、Arc Hydro、SAGA GIS或Python进行基础分析的读者参考。

背景:洪水淹没模拟常见应用场景

在GIS项目中,洪水淹没模拟常用于防洪规划、应急预案、城市内涝分析、库区回水影响分析、河道整治评估和灾害风险制图。不同场景对精度要求不同,报告写法也应有所区分。

1. 课程设计或实训报告

课程设计通常更重视流程完整性。报告需要展示DEM预处理、流向流量计算、河网提取、集水区划分、淹没范围生成和地图制图成果。即使没有实测水文数据,也要说明假设水位或模拟情景。

2. 初步风险筛查

如果目标是快速识别低洼易涝区,可以使用DEM和设定水位进行静态淹没分析。优点是操作快,适合宏观筛查;缺点是不能准确表达洪水传播过程、流速、堤防阻隔和排水设施影响。

3. 工程或应急分析

如果用于工程决策或应急调度,通常不能只依赖GIS静态叠加。需要结合水文模型、水动力模型、实测水位站、河道断面、糙率参数和边界条件,GIS主要负责数据准备、空间叠加、成果表达和影响统计。

原理:GIS洪水淹没模拟的基本逻辑

GIS洪水淹没模拟可以简单理解为:在统一坐标和高程基准下,把“可能达到的洪水水面”与“地表高程”进行比较。如果某处地表高程低于洪水水面,并且与河道或水体连通,就可能被判定为淹没区。

常见方法主要有三类。

1. 静态水位淹没法

静态水位淹没法是最容易入门的GIS方法。它假设研究区存在一个统一或分区的洪水水位,例如20米、25米或设计洪水位。将DEM中低于该水位的区域提取出来,即可得到初始淹没范围。

基本判断逻辑如下:

淹没深度 = 设定洪水水位 - DEM地面高程
当 淹没深度 > 0 时,像元可能被淹没

但要注意,“低于水位”不一定等于“会被淹”。如果低洼区与河道不连通,或者被堤防、道路高填方、山脊阻隔,直接提取会产生虚假淹没斑块。因此报告中应增加连通性检查。

2. 基于河道水位线的淹没分析

在河流型洪水中,不同河段水位通常不同。更合理的做法是根据水位站、断面或模型输出结果,沿河道建立水面高程面,再与DEM比较。这样比全区统一水位更接近实际。

这种方法常见于河道洪水影响范围分析,关键在于水位点、断面线和河道中心线的空间位置必须准确。

3. 水文水动力模型结合GIS

HEC-RAS、MIKE、InfoWorks、LISFLOOD-FP等模型可以模拟水位、流速和洪水传播过程。GIS在其中承担数据准备和成果制图工作,例如生成地形、土地利用、糙率分区、断面、边界条件,以及把模型输出的水深栅格叠加到道路、建筑物、人口和土地利用数据上。

步骤:GIS洪水淹没模拟报告的水文分析流程

下面给出一套适合写入报告的方法流程。实际项目可以根据数据条件删减,但建议不要跳过坐标、高程、DEM质量和结果核查这几个关键环节。

步骤1:明确研究范围和模拟情景

报告开头应先说明研究对象,而不是直接进入软件操作。

  • 研究区:例如某河段两岸、某城区、某水库下游或某乡镇范围。
  • 模拟目标:识别淹没范围、计算淹没深度、统计受影响要素。
  • 洪水情景:例如10年一遇、20年一遇、设计洪水位、历史最高水位或假定水位。
  • 输出成果:淹没范围图、淹没深度图、风险分级图、面积统计表、影响对象清单。

如果是课程报告,可以写明“本次分析采用静态水位法,以DEM为基础,在假定洪水水位条件下提取淹没范围”。这样能避免把简化模型说成精确预报。

步骤2:准备基础GIS数据

洪水淹没模拟对数据质量非常敏感。建议至少准备以下数据:

数据类型 用途 注意事项
DEM 计算地面高程、坡度、流向、淹没深度 分辨率越高越好,需检查高程单位和垂直基准
河流或水系 确定洪水来源和连通区域 河道位置要与影像或DEM沟谷一致
行政区或研究区边界 裁剪数据和统计成果 坐标系需统一
土地利用、建筑、道路 进行影响分析 最好使用最新数据
水位或流量资料 设定模拟情景 需说明来源、年份和代表性

如果使用ArcGIS Pro,可用“投影栅格”“裁剪栅格”“填洼”“流向”“流量”等工具;如果使用QGIS,可结合Processing工具箱中的GDAL、SAGA和GRASS工具完成类似处理。

步骤3:统一坐标系和高程基准

这是很多水文分析步骤中最容易被忽视的一步。水平坐标系影响面积和距离,高程基准影响淹没深度判断。如果DEM是米单位,而水位资料是另一套基准或单位,结果会整体偏高或偏低。

建议检查:

  • 所有矢量和栅格是否在同一投影坐标系下。
  • DEM高程单位是否为米。
  • 水位高程与DEM是否使用同一高程基准。
  • 研究区是否适合使用经纬度坐标进行面积统计,通常应转换到投影坐标系。

报告中可以写明:“为保证面积和高程计算一致,本文将所有数据统一至同一投影坐标系,并在分析前检查DEM高程单位与水位资料单位。”

步骤4:DEM预处理

DEM预处理直接决定后续水文分析结果。常见处理包括裁剪、异常值检查、填洼和平滑处理。

  1. 按研究区边界裁剪DEM,减少计算量。
  2. 检查DEM是否存在无效值、异常高程或空洞。
  3. 执行填洼处理,消除因数据误差形成的局部洼地。
  4. 根据需要进行河道刻蚀或水系校正,使DEM流向更符合实际河网。

填洼不是越强越好。如果研究目标是城市内涝或真实低洼地识别,过度填洼可能抹掉实际洼地。因此报告要说明填洼目的:是为了保证水文汇流计算连续,还是为了保留真实洼地进行内涝识别。

步骤5:进行流向和流量分析

水文分析的基础步骤通常包括流向和流量计算。流向表示每个栅格像元的水往哪个相邻方向流,流量累积表示上游有多少像元汇入该位置。

常见流程如下:

  1. 输入填洼后的DEM。
  2. 计算流向栅格。
  3. 基于流向计算流量累积栅格。
  4. 设置流量阈值提取河网。
  5. 将栅格河网转换为矢量河流。

在报告中,流量阈值必须说明。阈值太小,提取出的沟道过密;阈值太大,支流会被遗漏。可以通过与遥感影像、已有水系图或实地资料对比来调整。

步骤6:划分集水区或汇水范围

如果分析对象是某个断面、桥梁、闸口、水库或城市排水出口,需要划分其上游集水区。集水区用于说明洪水来源范围,也可作为后续流量估算和降雨径流分析的空间基础。

常用步骤包括:

  1. 确定出水口或控制断面位置。
  2. 将出水口捕捉到高流量累积河道像元上。
  3. 基于流向数据计算上游集水区。
  4. 统计集水区面积、平均坡度、土地利用类型和不透水面比例。

如果出水口没有捕捉到正确河道,划分出来的集水区可能明显偏小或偏移,这是报告核查时必须重点检查的问题。

步骤7:设定洪水水位或导入模型结果

淹没模拟必须有明确的水位或水深来源。常见输入方式有三种:

  • 使用历史洪水位,例如某水位站最高水位。
  • 使用设计洪水位,例如10年一遇、50年一遇对应水位。
  • 导入水动力模型输出的水面高程或水深栅格。

如果使用静态水位法,可以设置一个或多个水位情景,例如:

情景 设定水位 用途
低水位情景 18 m 识别低洼敏感区
中水位情景 20 m 模拟一般洪水影响
高水位情景 22 m 评估极端洪水风险

如果不同河段水位不同,建议根据断面水位点插值生成水面高程栅格,而不是全区使用同一个常数水位。

步骤8:计算淹没深度和淹没范围

在GIS中,淹没深度通常通过栅格计算器完成。核心公式为:

淹没深度栅格 = 洪水水面高程栅格 - DEM高程栅格

如果使用统一水位,公式可以写成:

淹没深度 = 设定水位 - DEM

然后提取淹没深度大于0的区域:

当 淹没深度 > 0,判定为初始淹没区
当 淹没深度 <= 0,判定为未淹没区

在QGIS中,可以使用“栅格计算器”生成水深栅格,再用“重分类”或“按表达式提取”得到淹没范围。在ArcGIS Pro中,可以使用“Raster Calculator”“Con”“Reclassify”等工具。

步骤9:进行连通性修正

仅凭高程低于水位会产生很多孤立低洼斑块。对于河流洪水,应检查这些区域是否与河道或主淹没区连通。

可采用以下方法:

  • 将初始淹没区转为面,删除面积很小且与河道不连通的斑块。
  • 使用区域分组工具识别连通区域,只保留与河道相交的区域。
  • 结合堤防、道路、铁路等线状阻隔物进行人工核查。
  • 对明显不合理的淹没孤岛在报告中说明处理原则。

这一步是让GIS洪水淹没模拟报告更可信的关键。否则结果图上会出现许多看似被淹、实际与洪水来源无关的低洼区域。

步骤10:分级制图与影响统计

淹没范围只是基础成果,报告还应统计影响对象。常见统计内容包括:

  • 总淹没面积。
  • 不同水深等级面积,例如0到0.5米、0.5到1米、1到2米、大于2米。
  • 受影响土地利用类型面积。
  • 受影响道路长度。
  • 受影响建筑物数量或居民点数量。
  • 重点设施是否位于淹没区内。

水深分级可参考实际需要,不建议盲目套用固定等级。对于人员通行风险,0.5米水深已经可能造成明显影响;对于农田或生态分析,分级标准可能不同。

步骤11:结果验证和合理性检查

任何洪水淹没结果都需要验证。即使是课程作业,也应加入基本的合理性检查。

  • 与历史洪水范围、灾情点或遥感水体范围对比。
  • 检查淹没区是否沿河道、低洼地和谷地分布。
  • 检查高地、山脊和明显高程区域是否被错误淹没。
  • 检查河道两岸是否存在不合理断裂。
  • 检查面积统计是否受投影坐标系影响。

如果没有实测资料,也可以使用高分影像、地形阴影图、等高线和已有水系图做逻辑核查,并在报告中说明验证条件有限。

常见坑:GIS洪水淹没模拟报告最容易出错的地方

1. 坐标系不统一导致面积错误

在经纬度坐标系下直接计算面积,结果往往不适合用于报告。应将数据转换到适合研究区的投影坐标系后再统计面积和长度。

2. 高程基准不一致导致水深整体偏差

水位和DEM必须在同一高程基准下比较。如果一个是1985国家高程基准,另一个来自未知来源DEM,可能产生系统误差。报告中应说明数据来源和可能误差。

3. 把静态淹没图当作洪水过程模拟

静态水位法不能表达洪水到达时间、流速、涨落过程和排水过程。它适合快速识别潜在淹没范围,不适合作为精细工程调度依据。

4. 忽略堤防和道路阻隔

DEM分辨率较低时,堤防、道路高填方、挡墙等小尺度阻隔物可能无法体现。结果可能显示水越过堤防进入保护区。若研究区有堤防工程,应单独加入堤防高程或人工修正。

5. DEM分辨率不足

30米DEM适合区域尺度分析,但不适合精细到街区、建筑物或道路级别的城市内涝模拟。城市洪水分析最好使用更高分辨率地形数据,并考虑雨水管网和地表排水条件。

6. 没有说明水位来源

报告中只写“设定水位为20米”是不够的。应说明该水位来自历史记录、设计标准、模型输出,还是教学假设。来源不明会削弱报告可信度。

方法比较:静态水位法、GIS水文分析和水动力模型怎么选

方法 适用场景 优点 局限
静态水位淹没法 课程实训、快速筛查、缺少水文资料的初步分析 操作简单,GIS即可完成 不能表达洪水传播过程和流速
GIS水文分析 流域识别、河网提取、汇水区划分、辅助风险判断 能解释地形汇流关系 不能单独替代水动力模拟
HEC-RAS等水动力模型结合GIS 河道洪水、工程评估、应急预案 可输出水位、水深、流速和过程结果 需要断面、糙率、边界条件和较高建模经验
遥感水体提取结合GIS 历史洪水范围验证、灾后评估 可基于真实影像识别已发生淹没 受云雨、影像时间和水体识别精度影响

对于初学者,建议先掌握静态水位法和基础水文分析步骤,理解DEM、水位、流向、流量和连通性的关系;如果要做更严谨的河道洪水评估,再学习HEC-RAS与GIS的联合建模。

检查清单:提交GIS洪水淹没模拟报告前逐项核对

  • 是否明确写出研究区、模拟目标和洪水情景。
  • 是否列出DEM、水系、边界、水位、土地利用等数据来源。
  • 是否统一坐标系,并说明面积统计使用投影坐标系。
  • 是否检查DEM高程单位和水位高程基准。
  • 是否完成DEM裁剪、填洼和异常值检查。
  • 是否说明流向、流量、河网提取和集水区划分方法。
  • 是否说明洪水水位或水面高程来源。
  • 是否计算淹没深度,而不仅仅输出淹没范围。
  • 是否对孤立低洼斑块做连通性检查。
  • 是否统计淹没面积、不同水深等级和受影响对象。
  • 是否对结果进行影像、历史洪水或地形逻辑验证。
  • 是否在结论中说明方法局限和不确定性。

FAQ:GIS洪水淹没模拟报告常见问题

1. GIS洪水淹没模拟报告必须使用HEC-RAS吗?

不一定。如果是课程练习、基础风险识别或缺少水文断面资料,可以用DEM和设定水位做静态淹没分析。但如果报告用于工程评估、河道治理或应急调度,建议结合HEC-RAS等水动力模型。

2. 水文分析步骤和洪水淹没模拟有什么区别?

水文分析步骤更关注地形汇流关系,例如填洼、流向、流量、河网和集水区。洪水淹没模拟则更关注在特定水位或水动力结果下,哪些地表区域被水覆盖。两者可以结合,但不是同一个概念。

3. 只有DEM,没有水位资料,能做洪水淹没模拟吗?

可以做情景假设分析,例如设置不同水位进行敏感性测试。但报告必须明确说明这是“假定水位情景”,不能把结果写成真实洪水预报。

4. 为什么低于水位的地方没有全部作为淹没区?

因为洪水需要与河道或水源连通。某些低洼地虽然高程低于设定水位,但可能被山脊、道路、堤防或其他高地阻隔。报告中应进行连通性修正,避免虚假淹没。

5. QGIS可以完成GIS洪水淹没模拟报告吗?

可以。QGIS可以完成DEM裁剪、填洼、流向流量分析、栅格计算、淹没范围提取、矢量叠加和地图制图。复杂水动力模拟仍建议结合专门模型软件。

6. ArcGIS Pro中常用哪些工具?

常用工具包括Fill、Flow Direction、Flow Accumulation、Watershed、Raster Calculator、Con、Reclassify、Raster to Polygon、Intersect、Tabulate Area等。具体工具名称会随扩展模块和版本有所不同,使用前应确认已启用相应分析工具箱。

7. 淹没深度图应该如何分级?

常见做法是按0到0.5米、0.5到1米、1到2米、大于2米分级。但分级应根据研究目的调整,例如道路通行、居民风险、农田受灾或设施安全的关注阈值并不相同。

8. 报告结论应该写什么?

结论应概括主要淹没范围、最大或主要水深区域、受影响对象、风险集中区域,并说明数据精度、模型假设和结果适用范围。不要只写“完成了洪水淹没模拟”。

结论:一份合格报告要讲清数据、方法、结果和局限

GIS洪水淹没模拟报告的重点不是把软件步骤堆在一起,而是建立清晰的分析链条:数据是否可靠,水位从哪里来,DEM如何处理,淹没范围如何提取,结果如何验证,局限在哪里。

如果是入门学习,可以从静态水位法开始,掌握DEM预处理、流向流量分析、栅格计算和空间叠加统计。若用于更严肃的防洪决策,则应把GIS作为数据处理和成果表达平台,结合水文水动力模型、实测水位和工程资料共同完成分析。

最终,好的洪水淹没模拟报告应让读者看完后能够判断:在什么条件下、哪些区域可能被淹、淹多深、影响什么,以及这个结果能用到什么程度。