ArcGIS DEM水文分析与山洪风险评估

关键词：ArcGIS | DEM | 水文分析 | 流域划分 | 山洪风险

一、研究方法：DEM水文分析与Google Earth人工解译

Hamdy等（Geosciences 2016）以埃及阿斯旺省Abouelreesh村为案例——该村位于红海山脉西麓，曾于2010年遭受山洪导致25座房屋完全损毁、50余座部分毁坏。研究目标是通过免费公开数据实现山洪风险区识别与城市增长趋势评估，为缺乏水文实测资料的欠发达地区提供可复用的风险评估方法。

研究分为两条并行主线：(1) DEM水文分析——使用SRTM 30m分辨率DEM数据，通过ArcGIS Spatial Analyst水文工具箱提取流域边界和地表水系网络（Table 1列出了完整11步工具链：Fill填洼→Flow Direction流向计算→Flow Accumulation流量累计→Stream Link水系链接→Stream Order分级→Stream to Feature栅格转要素→Flow Length流长→Basin流域划分→Snap Pour Point汇流点捕捉→Watershed子流域提取→Raster to Polygon栅格转面）。提取出三个子流域覆盖案例区全境，流域边界内部即为山洪风险区。(2) Google Earth历史影像人工目视解译——在Google Earth桌面应用中调取2001/2005/2009/2013四个年份的高分辨率历史影像，通过人工判读逐栋标注建筑轮廓和建筑年代，生成四个时段的建筑矢量图斑。该方法仅需照片解译培训和Internet连接，成本极低但可靠性高——多篇文献证实人工目视解译在精细尺度城市用地分类中优于半自动分类方法。

图1 DEM高程数据与ArcGIS水文分析提取的水系网络：(a)案例区DEM文件（像元值5~458m）；(b)水系网络提取结果（Hamdy et al., 2016）

图2 Google Earth历史影像人工目视解译的Abouelreesh村建筑年代分布图：2001年前（蓝色）、2001-2005年（绿色）、2005-2009年（黄色）、2009-2013年（橙色）

二、结果分析：城市增长持续向风险区倾斜

将ArcGIS提取的流域边界与Google Earth解译的建筑图斑做空间叠加（Intersect），落在流域边界内部的建筑标记为「风险区」、外部标记为「安全区」。四个年份的叠加统计数据（Table 2）揭示了清晰的趋势：风险区建筑占比从2001年的53.64%持续攀升至2005年71.35%、2009年74.69%、2013年82.13%。12年间总建筑面积从118.02ha增至184.66ha，其中风险区增长51.79ha（占地83%的新增面积），安全区仅增长14.85ha。

2009年是最大增长峰值年——新增24.56ha建筑全部位于风险区。年增长率数据（Figure 10）显示风险区年增长率从2005年的5.19%飙升到2009年的36.88%，随后回落至2013年的26.26%。图9指出78.65%的安全区建筑建成于2001年以前，而风险区仅55.01%，说明新增开发更倾向择址于流域覆盖区域。

建筑尺寸演化进一步证实了这一趋势（Table 4）：风险区房屋平均面积从2001年的239m²增至2013年的365m²（+53%），安全区仅从235m²增至272m²（+15%）。到2013年，风险区住宅面积已达到安全区对应数值的130%以上。论文指出这与地形直接相关——安全区等高线坡度更陡，居民为追求更大的居住空间而向地势更平缓的风险区扩张。

图3 ArcGIS叠加分析结果：DEM提取的流域边界（山洪风险区，彩框）与安全区（无框）划分。三个子流域覆盖全村，风险区占比从2001年53.64%升至2013年82.13%

图4 2001-2013年风险区与安全区城市建筑面积对比柱状图：12年间总建成区从118.02ha增至184.66ha，其中风险区从63.31ha增至115.10ha（Hamdy et al., 2016）

三、关键发现与规划启示

论文最后对政府2020版总体规划进行了风险审计（Table 5）：在划定的97.40ha拟建用地中，59.63ha（61.22%）位于DEM提取的流域边界内部，仅37.78ha（38.78%）位于安全区。这一比例与现状建成区几乎相同（62.33%风险 vs 37.67%安全），说明官方规划层面并未主动引导城市增长远离风险区——政府在风险区内规划了供水供电等公共服务，客观上鼓励了居民继续向危险区域扩张。

三个子流域中，Watershed 2面积最大（72.62ha，占比63.07%），是人口聚集密度最高的风险区。论文结论指出：山洪风险区城市增长并非仅由居民对大居住空间的追求驱动，也与政府规划决策密切相关。未来研究需要关注如何将城市增长引导至安全备选区域，并为风险区现有建筑制定工程防护方案。该方法完全基于免费DEM数据和Google Earth影像，无需实地水文站或付费遥感数据，特别适合数据匮乏的发展中国家村镇级风险评估。

我们团队在GIS空间分析与自然灾害风险评估领域积累了丰富的工程实践经验，可提供以下服务：①基于ArcGIS/ArcGIS Pro的DEM水文分析全流程——从SRTM/ASTER DEM下载、Fill填洼预处理到Flow Direction/Accumulation/Watershed流域边界提取，交付完整操作手册与中间栅格。②多时相城市增长与灾害风险叠加分析——Google Earth/Sentinel-2/Landsat影像建筑图斑目视解译、ArcGIS空间叠置与风险区定量评估、时空趋势制图。③流域山洪风险专题评估——SWAT/HEC-HMS水文模型建模、不同暴雨重现期淹没区制图、风险情景对比与防护策略建议。④城市总体规划风险审计与备选方案——将官方规划矢量数据与DEM提取的风险区边界叠加分析，输出规划调整建议与安全备选用地选址。⑤CA-Markov/SLEUTH城市增长模拟——结合DEM地形约束的未来城市扩张方向预测，支撑国土空间规划与灾害韧性决策。

图5 2001-2013年平均住宅建筑面积对比柱状图：风险区房屋从239m2增至365m2（+53%），安全区仅从235m2增至272m2（+15%）。2013年风险区房屋面积超过安全区130%

图6 政府2020总体规划拟建用地分布图：拟建区97.40ha，其中59.63ha（61.22%）位于流域边界内，37.78ha（38.78%）位于安全区（Hamdy et al., 2016）

如有需要，欢迎通过公众号"320科技工作室"联系我们。