洪水灾害损失分析

洪水灾害损失的分析

问题和数据分析

1.洪水灾害指标

（1）洪水灾害自然特征指标。
洪水灾害发生的位置：洪水灾害发生的地理位置或区域，自然位置以经纬度或投影坐标表示，社会位置以所属行政单元表示。

洪水灾害影响的范围：直接过水或受淹地区，自然影响范围用淹没范围图表达，社会影响范围以洪水所影响的行政管辖范围表达。

洪水淹没深度：受淹地区积水深度，洪水淹没深度是度量洪水灾害严重程度的一个重要指标，是评价洪水灾害损失的一个重要因子。

（2）洪水灾害社会特征指标。
人口指标：包括受灾人口、死亡人口、受伤人口和影响人口。

淹没土地利用类型：淹没范围内的土地利用现状。

房屋：洪水淹没、冲垮和破坏的各种房屋。

农作物：洪水长时间淹没或冲毁农田而造成农作物减产、绝收的面积或产量损失。

传染病：因长期水灾引起的疾病。

（3）洪水灾害经济损失指标。
财产损失率：财产损失率是指洪水淹没区各类财产损失的价值与灾前原有价值或正常年份各类财产价值之比。显然，确定了各类财产的洪水灾害损失率，乘以灾前原有各类财产的价值，就可以得到遭受洪水灾害后各类财产的损失值。财产损失率是基于居民经验和相关科技人员调查、统计得来的。不同土地利用类型在不同淹没深度下的损失率是不同的。

面上综合经济损失描述指标：除了财产损失率直接用来描述经济损失外，常用面上综合经济损失描述指标来描述经济损失，主要有亩均损失值指标、单位面积损失值指标和人均损失值指标。

2.问题提出

洪水淹没有一个最高水位，因此可以根据高程分布数据，与地块类型多边形进行联合分析，通过条件检索得到小于洪水最高水位的淹没区内不同的地块类型，根据不同地块类型的估计财产及损失系数等参数计算财产损失，分析准则如下。

（1）洪水水位的相对高程为500。
（2）估计住宅用地被洪水淹没而造成的损失。损失大小和居民的财产、地基稳定性有关，计算公式为：估计损失 = 面积×地均财产×损失系数。

3.数据准备

针对洪水淹没损失的评估条件，现已获得该地区数字化后的地块边界线land和等高线height，土地利用类型表landuse和损失系数表damage。

该地区land要素类的地块编号如图8.5.1（a）所示，height要素类的高程分布如图8.5.1（b）所示。land要素类的属性如表8.5.1所示。



土地使用类型中，R1为一类住宅用地，R2为二类住宅用地，C为公共设施用地。
对每一类地基，根据其地基类型，确定其稳定性及房屋倒塌的可能性，称之为损失系数，不同地基类型对应的损失系数如表8.5.2所示。

数据预处理

由于在分析中需要对多边形数据进行操作，而获取的原始数据为线类型的数据，并且原始线数据中的属性值也不能支持灾害损失估计的分析，所以要先将边界线表达的地块和等高线表达的高程生成多边形数据。

1、处理land数据

加载数据
首先加载数据，将land线要素类、landuse表和damage表添加到内容列表和地图窗口中。

要素转面
打开工具箱，选择数据管理工具-要素-要素转面，打开要素转面对话框
输入要素选择land图层，输出要素命名为landA，选择存放路径
点击确定，完成要素转面，如图所示。

为landA添加属性字段
由边界线要素类land转换成的landA面此时还没有面积、编号等属性，为了后续的分析，需要添加这些属性。打开工具箱，依次点击数据管理工具-字段-添加字段。

在添加字段对话框中，输入选择landA，字段名设置为area，字段类型设置为float，其他默认，表示面积，点击确定。

用同样的方法为landA定义一个lnum的字段，类型设置为LONG，表示地块编号；添加一个名为avgp的字段，字段类型为FLOAT，表示地均财产。

添加完成之后，在landA图层右键，选择打开属性表，可以看见添加的三个字段，目前均为空值，因此我们要对它们进行赋值。

在进行赋值前，需要启动编辑器，在landA图层右键，选择编辑-开始编辑，选择之后，属性表下面会多一行。

注：也可以直接在编辑工具条上面下拉编辑器，启动编辑。

随后对area字段进行赋值，在area字段右键，字体会变蓝，选择计算几何，打开计算几何对话框，将对话框中的属性设置为面积，其他不变，点击计算即可。

计算完成之后，点击编辑工具条，下拉，选择保存编辑，并停止编辑。

估计损失需要的土地使用类型、地基类型、估计财产、损失系数等都需要添加到landA要素类的属性表中。本例中通过属性连接将这些属性添加到landA的属性表中。
在图层landA右键，依次选择连接和关联-连接，打开连接数据对话框。
如图设置对话框数据，单击确定，将landuse表的属性连接到landA图层，为每个地块赋予了土地使用、估计财产和地基类型。
注：在连接之前，需要对地块赋值，右键landA，打开属性表，启动编辑，根据前面提供的land地块属性表，依次赋值。

同样的步骤，单击确定按钮，系统可能会弹出一个警告对话框提示未建立索引。单击是按钮建立索引，也可以单击否按钮不建立索引。是否建立索引不影响本例后续的分析。
即以地基类型为关键字将damage表的属性连接到landarea图层，为每个地块赋予了损失系数

计算avgp字段的属性值，avgp为地均财产，表示地块中单位面积上的财产值。
同理，打开landA的属性表，启动开始编辑，右键单击avgp字段名，在弹出菜单中单击字段计算器，打开字段计算器对话框。
Step2：在字段计算器中输入计算表达式：[landuse.估计财产] / [landA.area]，单击确定按钮，完成avgp的计算。

处理完之后的landA属性如图

2、处理height等高线数据

首先进行等高线转面，将height加载到地图中，在工具箱中依次点击数据管理工具-要素-要素转面，要素转面对话框中将输入要素设置为height，输出要素类命名为ht，同上方法一样。

由线要素类height转换成的ht面要素类此时还没有高程属性，为了后续的分析，需要添加高程属性。首先添加一个字段，在工具箱中打开数据管理工具-字段-添加字段，在添加字段对话框中将输入表设置为ht，字段名设置为height，字段类型设置为FLOAT，其他选项保持默认设置，单击确定按钮，完成height字段的添加。

右键ht图层，打开属性表，启动编辑器，按照前述的图片完成height的属性值添加，完成后点即确定并停止编辑。

此时数据处理完成。

受灾区域分析

根据分析条件，高程低于或等于500的居民地是受灾区域。在土地使用类型中属性值为R1和R2的两种用地类型为住宅用地。

首先启动编辑，使图层处于编辑状态。

1.检索出土地使用为住宅的区域

在ArcMap主菜单上依次单击“选择→按属性选择”，打开按属性选择对话框。

在按属性选择对话框中将图层设置为landA，方法设置为创建新选择内容，在属性文本框单击"landuse.土地使用"属性，再单击获取唯一值按钮，在选择表达式编辑框中输入"landuse.土地使用" = ‘R1’ OR “landuse.土地使用” = ‘R2’，

2.检索出高程低于或等于500的区域

在ArcMap主菜单上依次单击“选择→按属性选择”，打开按属性选择对话框。在按属性选择对话框中将图层设置为ht，方法设置为创建新选择内容，在属性文本框单击"height"属性，再单击获取唯一值按钮，在表达式编辑框中输入"height" ＜= 500，单击确定按钮，选出高程低于或等于500的区域

3.分析高程低于或等于500的住宅区域

在ArcMap主菜单上依次单击“地理处理→相交”，打开相交对话框。
在相交对话框中将输入要素设置为landA和ht，输出要素类命名为lh_inter。

处理完之后的图像。

灾害损失估计

灾害损失估计的计算公式为：估计损失 = 面积地均财产损失系数。式中，面积是被淹没的土地面积，lh_inter要素类中目前还没有表示面积的字段，需要添加。

1、添加并计算area字段

在ArcToolbox中依次单击“数据管理工具→字段-添加字段”，打开添加字段对话框。
在添加字段对话框中为输入表lh_inter添加一个字段名为area的字段，字段类型为FLOAT。

在lh属性表中右键单击字段名area，在弹出菜单中单击计算几何，打开计算几何对话框。在计算几何对话框中，将属性设置为面积，单击确定按钮，完成area字段的计算

2、添加并计算damage字段

步骤同上。在计算字段对话框中将输入表设置为lh，字段名设置为damage。根据损失估计的计算方法：估计损失 = 面积地均财产损失系数，输入表达式[area] * [landA_avgp] *[damage_损]，其他选项保持默认设置

3、将damage字段的值作为标注显示在地图上

注意：在ArcMap的内容列表中右键单击lh图层名，在弹出菜单中单击属性，打开图层属性对话框，在图层属性对话框中单击标注属性页。在标注属性页中，将方法设置为以相同方式为所有要素加标注，标注字段设置为damage

在ArcMap的内容列表中右键单击lh图层名，在弹出菜单中单击标注要素，将damage字段的值标注在地图上，

4、 将lh图层改变显示方式

在ArcMap的内容列表中右键单击lh图层名，在弹出菜单中单击属性，打开图层属性对话框。在图层属性对话框中单击符号系统属性页。在符号系统属性页中将显示设置为“类别→唯一值”，将值字段设置为damage，选择任意色带，单击添加所有值按钮

5、对损失进行统计

在ArcMap的内容列表中右键单击lh图层名，在弹出菜单中单击打开属性表，打开lh要素类的属性表。
在lh要素类属性表中右键单击damage字段名，在弹出菜单中单击统计，打开统计数据lh对话框，如图8.5.32所示。统计结果显示了damage字段，也就是洪水淹没估计损失的最大值、最小值、总和、平均值、标准差及频数分布等。

即计算完成。