基于GIS技术的洪水淹没模拟及灾害评估

基于GIS技术的洪水淹没模拟及灾害评估基于GIS技术的洪水淹没模拟及灾害评估姓名：孙珂单位：郑州测绘学校地址：河南省郑州市大学南路1号邮编：450000摘要：针对洪水灾害评估的特点 ，将GIS技术与RS技术相结合，根据数字高程模型DEM提供的三维数据和遥感影象数据，运用GIS的空间分析功能来预测、模拟显示红水河流域洪水淹没场景，并结合该流域水文站降雨量数据和各乡镇人口密度数据以及其他辅助数据来进行洪水灾害评估。关键字：洪水灾害、地理信息系统、泰森多边形、加权叠加1前言洪水灾害是最频发的自然灾害，严重影响国民经济发展危害人民生命财产安全，破坏生态环境。随着现代经济的高速发展和水利工程的增加，洪水灾害对人类的危害仍在加重。因此 ，快速、准确、科学地模拟、预测洪水淹没范围，对防洪减灾具有重要意义。特别是对于一些重点防洪城市和行蓄洪区，如果能够预先获知洪水的淹没范围和水深的分布情况，对于预先转移受灾区的生命财产，减少损失具有非常重要的价值，而且对于洪水造成的灾害损失进行评估也是非常有用的。近几年来，将GIS技术与RS技术相结合，根据数字高程模型DEM提供的三维数据和遥感影象数据来预测、模拟显示洪水淹没场景，并进行洪水灾害评估，已成为GIS在洪水方面主要研究领域。本研究以数字高程模型DEM和RS影象为基础，运用GIS的空间分析功能，研究试验区洪水河流域的洪水淹没情况。2研究区域及数据简介21 研究区域地理概括红水河是珠江流域西江水系的中上游河段，发源于云南省沾益县马雄山，流经滇、黔、桂三省（区），上游主流称南盘江，流至庶香双江口与北盘江汇合后称红水河，到广西三江口与柳江相汇合后称黔江。红水河流域位于东经1022010930，北纬23042650之间，流域四周为群山环绕，整个地势自西北向东南倾斜，平均海拔高程1450m。本次实验重点研究范围为红水河流域中的整个龙滩流域及其六个子流域(甲板、平腊、八茂、蔗香、这洞、高车)。722 实验数据 本研究采用的基本数据分为空间数据和水文数据以及其他辅助数据。其中空间 图1：研究区域红水河流域图数据包括龙滩流域的DEM底图、modis遥感影象底图、省市县行政边界、城市分布图、站点分布图、河网、龙滩流域及其子流域分布图等。水文数据是各子流域水文站获取的降雨量数据，辅助数据是流域流经区域经济数据以及为实现真实三维场景所采集到的部分建筑物纹理数据。3研究的主要任务及思路 通过利用研究区域的数字高程模型（DEM）和遥感影象数据对其进行三维场景仿真，再通过给定洪水水位高程值，认为在水位以下则是淹没区域，反之则不属于淹没范围，并对淹没过程进行模拟，通过人工方法确定连通区域来计算洪水淹没面积。同时结合研究区域土地利用类型数据、降雨量情况以及辅助经济数据进行加权叠置生成洪水灾害风险图，为防洪抗灾做决策。4研究区域三维场景建立。由于红水河流域覆盖面积广阔，短时间内不可能将整个流域的三维场景进行模拟，本次研究只对红水河流域某小山谷进行三维建模。为了真实的反映地物的纹理特征，我们采用的原始数据是quickbird遥感影象以及对应的DEM数据。三维场景的制作是通过将数字地形模型 (DEM)和遥感图像数据进行叠加，再将通过实地考察利用三维建模软件google sketchup建立好的建筑物模型和树木模型导入场景中，生成具有三维可视的地貌景观图（如图3）。在此基础上可以进行红水河流域水资源的研究 、洪涝灾害快速监测与评估及周边地表起伏形态特征等。技术路线如图2所示。Arcmap平面底图数字化Google sketchup建立模型场景导出为.mdb数据库花草树木模型符号化外业数据采集（纹理）Phtoshop纹理处理DEM数据遥感影象图红水河流域三维地形浏览可视化在Arcscene 9.3中进行叠加利用ESRI sketchup插件关联sketchup 关联图2 三维场景建立技术路线图3 三维可视化场景图5洪水淹没模拟分析及淹没面积计算洪水淹没模型是以三维地形和不同水位来综合演示洪水淹没行为的时空变化及水体形态运用 ArcGIS9.2中的3D扩展模块 ，在研究区域数字高程模型的基础上 ，在 ArcScene中实现了该区域洪水淹没演进过程，直观的显示不同时刻和水位该区域洪水淹没情况，为防洪救灾作出重要的决策。洪水淹没由多种因素造成，降雨、上游来水都可以造成淹没。按照洪水淹没的成因，可将其分为两大类：一类是无源淹没，一类是有源淹没。本次操作主要是利用arcmap软件完成，通过已知红水河流域各水文站水文数据确定洪水位，由已知DEM数据和给定洪水水位来确定淹没区域并计算淹没面积。具体操作流程在arctoolbox中新建model builder模型，（如图4），只需给定水文数据，可以直观的观察到其淹没区域及面积（如图5和图6）。输入水位高程图4：淹没面积计算model builder模型51 无源淹没分析只考虑受淹区的高程与给定水位的高程情况，而不用考虑淹没区的连通问题，凡是高程低于给定水位的点都记入淹没区，算作被淹没的点，这种情形相当于整个区域大面积均匀降水，所有低洼处都可能积水成灾。其淹没面积计算比较简单，所有低于或等于预测水位高程的像元都将计人淹没区，经累加计算得出淹没面积52 有源淹没分析水流受到地表起伏特征的影响，在这种情况下，即使在低洼处，也可能由于地形的阻挡而不会被淹没。造成的淹没原因除了自然降水外，还包括上游来水、洼地溢出水等。面积计算稍微有点复杂，它是在无源淹没的基础上，考虑到连通要求的淹没面积的计算。本文介绍的是用人为的方法确定符合连通条件的区域，将其取出，并计算出它的面积。具体做法是把决堤口定为源头，在转换生成的Raster中寻找与源头连通的区域。将所有 Raster的值为 1的区域进行累加计算，得出淹没区的面积。 （a） 无源淹没区 （b） 有源淹没区图5：水位为650m时红水河流域淹没情况 （a） 无源淹没区 （b） 有源淹没区图6：水位为800m时红水河流域淹没情况 6洪水灾害风险图制作：洪水风险图是以图的形式直观反映洪水威胁区域发生某一频率洪水后，可能淹没的范围、水深等洪水要素以及不同量级洪水可能造成的灾害风险和对社会经济的损害程度的工具。根据该图并结合影响区域内社会经济发展状况，合理制定防洪指挥方案，合理评价各项防洪措施的经济效益，合理估计洪灾损失，为防洪保险提供依据。本文研究的风险是一种相对的概念，即不同区域之间风险大小的相对比较。在 GIS的支持下，利用洪水淹没计算得到的淹没范围、淹没水深，与各种专题图层(包括行政区划、土地利用、居民点分布、人口分布、重要地物分布、交通线路等)进行拓扑叠加和空间分析 (空间查询、网络分析等)，即可生成有关专题淹没图件，包括重要设施淹没图、抢险迁安路线图等。如进一步关联背景数据库中的社会经济信息，并结合相关数学模型，可统计由洪水淹没造成的灾害损失。表一：红水河龙滩流域水文站日降雨量数据站码站点名称2005年降雨量（毫升）2006年降雨量（毫升）2007年降雨量（毫升）5月6月7月5月6月7月5月6月7月LT060001天生桥167.5321136194.5369.5192.5111.5283.5277.5LT060011盘江桥71.53183858126380181315305.5LT060012这洞57190571165353.5157.513953735*LT060017平湖219.5196.547139.534718261.5239.5353.5LT060019高车56.5338.5168.5107.519622915534554395LT060051大田河5429191.5122.5282.5114152.5294.5406LT060078惠水15636997.5173265183217.5337.5390.5LT060134蔗香146.533721777197204.5 *LT060143甲板258284.582182.52882565752492075LT060144八茂17133782.5216.5288.5202.56752734055LT060146平腊71201.561162323171512085205561 评价因子及评价原则编制洪水风险图有三大要素，即地形地貌、淹没水位和社会经济数据。其中本次研究主要有4个评价因子，包括高程、土地利用类型图、各子流域降雨量及人口密度。其中高程值是由红水河流域DEM数据生成；土地利用类型图利用原始TM遥感影象通过监督分类后生成；各子流域水文站测取的2002-2007年日降雨量数据，由于各水文站全年平均日降雨量比较接近，再加上洪水灾害一般发生在5、6、7三月，本次研究中各水文站平均日降雨量数据只取这四个月数据（如表一）。通过将该流域11个水文站数据利用arcgis生成泰森多边形，并转换生成红水河面雨量栅格数据。人口密度数据，通过利用乡镇行政区划图以及搜集到的各乡镇面积和人口数，将人口数除以乡镇面积并转换得到人口密度图。62 制作基本流程首先对淹没致灾的各因子进行分析，分别计算各影响因素形成洪水淹没可能性大小的比例，制作成相应的栅格图层；在 GIS技术的支持下，对各栅格图层重分类后进行加权叠置分析（具体操作流程见图7 ），对各影响图层确定评价系数以及子类型进行打分，得到洪水风险等级（如图8）。重分类后叠加流域DEM数据流域modis 遥感影象水文站降雨量数据乡镇人口数据土地利用类型图DEM图降雨量分布图洪水风险图人口密度分布图 图7：洪水风险图制作流程 图8：红水河流域洪水灾害风险图7结语：本文通过将GIS空间分析方法应用在红水河流域洪水灾害模拟及损失评估方面，初步得出以下结论：（1）利用DEM数据以及遥感影象数据和建筑物属性数据可以真实的模拟该地区洪水淹没真实场景，对于确定洪水淹没区地形起着重要的作用。（2）依据洪水水位确定洪水淹没范围后，可以结合其土地利用类型以及经济数据，对可能受淹地区的林地、耕地、居民用地等以及人口数据等进行快速评估，并进一步预测洪水灾害损失。（3）充分利用该地区水文数据，初步预测近几年洪水水位以及降水量情况，结合该区域经济数据，制作洪水灾害风险图，估计洪水发生时影响范围，这样可以为防灾救灾作决策支持。本次研究主要是通过arcgis空间分析功能，整个实现过程无需编程，而且方法简单适用，值得广泛应用。【参考文献】1刘仁义,刘南.基于 GIS的复杂地形洪水淹没区计算方法J地理学报，2001年6月2袁红梅,邓彩琼.应用GIS技术制作洪水风险图.江西水利科技,2004年l2月3丁志雄,李纪人,李琳 .基于GIS格网模型的洪水淹没分析方法.水利学报，2004年6月4刘小生,黄玉生.基于Arc/Info 的洪水淹没面积的计算方法。测绘通报，2003年第6期5郭利华,龙毅.基于DEM的洪水淹没分析.测绘通报,2002年第11期6张彦召.基于三维可视化技术的流域仿真模拟研究,硕士论文7翟宜峰,殷峻暹.基于GISRS的洪水灾害评估模型,人民黄河,2003年4月8向素玉,陈 军.基于 GIS城市洪水淹没模拟分析,中国地质大学学报,l995年9月9冯丽丽,李天文等.基于ArcGIS的渭河下游洪水淹没面积的计算,干旱区地理,2007年11月10封殿波,赵 杰等.基于ARCGIS的太平湖可视化研究,遥感与航空摄影,2008年第4期11李红英.基于GIS的洪灾损失评估研究以黑河为例，西安理工大学，硕士论文