基于gis的洪水灾害淹没损失计算模型.docxVIP

基于gis的洪水灾害淹没损失计算模型 信息系统是与地理学、计算机科学、遥感等学科和技术一起交换、处理和分析数据的结果。如果说第一个语言是文献，第二个语言是地图，地理数学是第二个语言，信息系统是第二个语言。根据数字描述空间的本质。信息系统是资源和环境研究、城市规划管理、国土开发调整、洪水风险分析、防洪救灾等服务的重要工具。在美国、加拿大、日本和其他发达国家，评估洪水灾害损失和整体管理洪水是应用信息的一个重要方面。 1 洪水灾害损失计算流程 洪水灾害淹没损失计算涉及众多的影响因子,这些因子包括区域环境背景因子(如地形、坡度、土地利用)、洪水特征(如流量、水位、洪水重现期)和社会经济因子(如人口、农业、工业)等,其中土地利用类型的分布,洪水淹没水深、历时与范围是洪灾损失评估的重要因子. 在GIS支持下,利用洪水水流模型计算出的淹没范围、淹没水深、淹没历时,各种专题图层,以及DEM等来计算灾害损失情况和再现洪水场景的详细流程,如图1所示. 具体步骤如下: (1)把洪水水流模型计算的淹没范围,转换为Arc/Info的淹没面积图层. (2)根据淹没水深估算灾情的程度,按水深分级,然后将此数据导入Arc/Info,并用不同的颜色来表示水深的分布情况. (3)使用淹没面积图层与人口分布图层、耕地面积图层、道路分布图层、乡镇分布图层和土地利用图层叠加分析,得到淹没区的人口数量、耕地面积、淹没道路等数据. (4)利用DEM与各时段的淹没范围、水深,在Arc/Info的GRID模块中,用GIS的图形处理和显示功能再现洪水淹没场景. 当然在求算淹没面积时,一定要扣除淹没前的水面面积,即原水面面积. 洪水灾害损失计算需要建立洪水水流模拟模型,以便再现或预测洪水灾害过程,确定洪水灾害发生的频率、量级、危害区域以及对象的变化规律与趋势,是洪水灾害损失计算的基础. 洪水淹没范围信息可以从遥感影像中提取洪水边界,但是它只能用于已经发生的洪水,对于预报洪水无能为力,而洪水水流模拟模型无论对已发生的洪水还是预报洪水均适用. 2 洪水流量模拟模型 2.1 控制方程的建立 由于洪泛区是大面积的浅层短期蓄水,其洪水波是浅水长波,水中压强服从静水压强分布. 洪水演进过程是湍流的非定常运动过程,可以采用笛卡尔坐标的圣维南方程为控制方程: 式中:Z是水位;H是水深;U,V分别是x,y方向的平均流速;Z=Zb+H,Zb为河底高程. 式中:f1为河床底部摩阻系数;f2为局部摩阻系数;可令Q=HU,R=HV,分别为x,y方向的单宽流量;C为谢才系数;n为河床糙率;Ω为哥氏力系数;γ为紊流涡旋性系数;λ为风应力系数;Ua,Va为分别是x,y方向的风速. 2.2 适体坐标变换计算 为适应分滞洪区的地形变化,采用将不规则区域变换为矩形区域的适体坐标变换方法,计算中考虑地形起伏使网格变形. 在洪水水流计算中,使用坐标变换、交替方向、交错网格、显隐格式差分求解的方法. 具体算法如下: 2.2.1 流圣维南方程 为适应河流、水库、分洪区不规则的边界形状,采用曲线坐标(δ,η)代替笛卡尔坐标(x,y),变换关系为 令变换系数αi,j的转置阵为 此为坐标变换的雅可比矩阵,其逆矩阵为 式中:J为α的行列式,J=xδyη-xηyδ;M是以α的代数余子式为元素组成的矩阵. 经过变换后的一般曲线坐标下的二维非恒定流圣维南方程为 式中:U,V为原x,y轴方向的速度分量,而不是δ,η轴方向的切线分量;Zδ,Zη,Uδ,Uη分别是Z对δ,Z对η,U对δ,U对η的导数. 在(11)、(12)、(13)式中左边第一项为时间变化项,是非恒定项,若没有这一项,流动就是恒定的;动量方程(12)、(13)中,左边第二项为对流项,具有非线性性质;左边第三项为反映重力作用的水位变化项;左边第四项为底部摩阻项;左边第五项为哥氏力项;右边第一项为涡粘性应力项,右边第二项为风应力项. 用交替方向进行计算,不计涡粘性及风应力,沿δ方向用下面一组方程: 沿η方向用下面一组方程: 为了保持计算的稳定,用交错网格,其中Z,U,V,H,X,Y分别处于不同的网格点上. 2.2.2 隐格式离散求解 前半步在δ方向求Z,U用隐格式差分离散,随后用显格式求V. 解(14a)、(14b): 其中: 整理式(16a),(16b),结合边界条件,可以求得Zn+1/2,Un+1/2,Vn+1/2. 同样方法,在η方向,进行后半步,由n+1/2到n+1步,对式(15a),(15b)进行隐格式离散求解,结合边界条件先解得Z,V,再显格式求得U. 3 基于dem的洪水模拟 洪水场景模拟是利用计算机软、硬件系统,模拟或再现洪水淹没的发展过程,给人以动态视觉效果. 系统使用Arc/Info软件的GRID和TIN模块来实现. 首先把防洪区的DEM转换为GRID,显示研