山东省干旱灾害风险评估与区划研究_李臻.pdf

第 卷 第 期 年 月人 民 长 江 ，收稿日期：基金项目：山东省重大科技创新工程项目（）作者简介：李 臻，女，硕士研究生，研究方向为自然灾害风险评估。：通信作者：季 民，男，教授，硕士生导师，博士，研究方向为空间数据组织和 系统应用开发。：文章编号：（）引用本文：李臻，季民，李云岭，等 山东省干旱灾害风险评估与区划研究 人民长江，（）：山东省干旱灾害风险评估与区划研究李臻，季民，李 云 岭，王 映 然（山东科技大学 测绘与空间信息学院，山东 青岛）摘要：为了提高山东省干旱灾害的风险评估和应急管理能力，减少旱灾损失，根据自然灾害风险理论，考虑干旱灾害的影响因素，从危险性、敏感性、易损性、防灾减灾能力 个方面选取评价指标，采用 熵权法确定指标权重，构建干旱灾害风险指数模型，并绘制出山东省干旱灾害风险区划图。风险区划结果表明：致灾因子的高危险区位于鲁西北地区，孕灾环境敏感性高值区主要位于中部及东部地区，承灾体易损性从阳信县、汶上县、陵城区、曹县、嘉祥县、宁阳县依次递减，防灾减灾能力从高到低依次为福山区、龙口市、荣成市、黄岛区、垦利区；从干旱风险区划图中可以看出，干旱综合风险大体呈现出由东向西逐渐递增的现象。关 键 词：干旱灾害；风险评估；风险区划；熵权法；山东省中图法分类号：文献标志码：引 言联合国政府间气候变化专门委员会（，）在其一系列报告中指出，干旱是世界上最常见的自然灾害。根据中国气象灾害大典山东卷灾情统计，年，除了 年未发生旱灾外，其他年份山东省均发生了不同程度的旱灾。全省平均每年农作物旱灾受灾面积为 万，成灾面积为 万。其中，年农作物的受旱面积达到最大，为 万。年，连续 大旱导致粮食短缺、企业停工、城镇供水严重短缺。年，山东省同样出现了严重的旱灾。年，山东省再次遭遇严重干旱，部分城镇供水受到影响，月，南四湖的蓄水量仅为常年同期的一半，水量亏缺导致水质下降，生态环境受到影响。由此可见，山东省常年发生旱灾，限制了农业生产和经济发展。近年来，灾害风险引起了国内外学者的高度关注，许多学者对干旱灾害风险评估与区划方面进行了广泛而深入的研究。等提出了基于干旱发生概率和潜在影响的干旱动态风险评估概括模型；等应用灰色物元分析方法对河南省农业旱灾进行风险评估，评估结果与旱灾发生的实际情况基本一致；等利用二维和三维 函数在中亚地区构建多维联合分布，并基于联合概率和重现期研究干旱风险特征；董涛等以安徽省淮河流域为研究区域，引入正态云模型，并基于该模型建立了农业旱灾风险评估模型；根据干旱风险指数（）危 险 指 数（）脆 弱 性 指 数（）基本模式，对土耳其各省的干旱综合风险进行评 估；刘 小 艳、王 志 春、杨 丰 政、尚 小宁、刘航等基于自然灾害风险理论，选取评价指标并建立了干旱灾害风险指数模型，分别对陕西省、赤峰市、徐水县、泾阳县以及淮河流域进行了干 第 期 李 臻，等：山东省干旱灾害风险评估与区划研究旱灾害风险评估。以往关于山东省干旱的研究多集中于干旱特征和时间序列分析方面等，但关于该省干旱灾害风险区划和评估方面的研究较少。鉴于此，本文基于自然灾害风险理论，根据山东省的实际情况，从自然灾害形成机制的角度出发，在参照前人研究的基础上，考虑并选择干旱风险四要素中的若干评价指标，对所选指标间的相关性进行分析，采用 熵权法确定相应的权重，并采用加权综合评价法构建了山东省旱灾风险评估模型。利用该风险评估模型，以 个地市的 个县区作为评价单元，开展了山东省干旱风险评估和区划研究。研究成果可为相关单位防灾减灾政策的制定提供一定的参考依据。研究区概况与数据来源 研究区概况山东省（，）位于中国东部沿海，地处黄河下游，全省面积近 万。如图 所示，山东省地形特征总体呈中部山地突起，西南、西北低洼平坦，东部缓丘起伏的态势。地貌类型以平原和丘陵为主，平原、丘陵和山地分别约占全省总面积的，和。土地利用类型主要包括耕地、林地、草地、水域、城乡居民用地等，耕地是山东省面积比例最大的地类。该省气候温和，降水集中，四季分明，属于暖温带季风气候。各地年平均气温在 之间，自鲁西南向半岛递减；年平均日照时数在 之间，自鲁南向鲁北递增；年平均降水量在 之间，自鲁西北向鲁东南递增，且全年 以上的降水量主要集中在 月，其他季节降雨量较少，因此该省干旱发生的频率高、持续时间长、范围广。图 山东省高程分布 数据来源气象资料来源于中国科学院资源环境科学与数据中心（：？），主要包括 年山东省 个气象站的逐日气温、降水量、日照时数、相对湿度以及平均风速等数据，剔除建站时间较晚以及缺失年份数据较多的站点，最后选取 个站点的气象资料。社会经济数据来源于 年各市的统计年鉴，主要包括山东省 个县（区）的年末总人口、国民生产总值、粮食播种面积、粮食产量、有效灌溉面积、大牲畜年末存栏数等。土地利用类型数据来源于中国科学院资源环境科学与数据中心。河网数据源于国家基础地理信息中心（：？）。分辨率数字高程模型数据来自地理空间数据云网站（：）。研究方法 指标规范化法为了消除灾害风险评价中各要素下指标之间的量纲和数值大小的影响，使评价具有可比性和一致性，需要对指标进行标准化处理，使各指标的数值分布在，范围内。需要注意的是，防灾减灾指数的数据范围需要置于，范围内，具体的计算公式如下：（）式中：为因子 的第 个指标的规范化值；和为第 个指标值的最小值和最大值；为因子 的第 个指标值。用该公式进行防灾减灾指数的指标处理时，不需要加 。确定指标权重方法 层次分析法层次分析法（，）是把与决策有关的元素分解成目标、准则和方案等层级，是一种定性和定量相结合的简便、快捷计算决策权重的研究方法，层次分析法求解权重的步骤如图 所示。本文利用 软件求解指标权重，最后检验判断矩阵的一致性比例是否小于 ，若小于，则认为判断矩阵通过一致性检验。熵值法熵值法是利用信息熵根据各指标的变异程度计算 人 民 长 江 年图 法求解权重的主要步骤 出各指标的熵值，然后用熵值修正各指标的权重，从而得出比较客观的指标权重。确定权重的计算方法及步骤如下：（）如果有 个评估指标和 个评估单元，则第 个指标的第 评价单元的值为，则指标原始数据矩阵 ()，其中 ，；，。（）对 的指标数据进行归一化处理，得到归一化矩阵为 ()。（）计算第 个指标下第 项评估对象的比重。（）（）计算第 个评价指标的熵值。（）（）计算第 个评价指标的权重系数。（）组合权重法组合权重法引入了距离函数的概念，通过计算层次分析法和熵权法两组权重之间的距离，使其等同于分配系数之间的差异程度，最后采用线性组合法获得干旱灾害风险评估体系中评价指标的组合权重。确定组合权重的表达式为 （）式中：为组合权重；、为权重的分配系数，且、需要满足 ；为层次分析法得到的权重；为熵权法得到的权重。根据上述内容，构造的方程组如下：，()()，()()（）通过求解方程组（）可以得到分配系数 和，将分配系数代入公式（）即可得出组合权重。加权综合评价法采用加权综合评价法计算灾害风险评价中 个因子（危险性、敏感性、易损性、防灾减灾能力）的指数。该方法将各评价指标的权重与指标值相乘，得到最终的综合值。具体公式为 （）式中：为 个因子加权求和后的值；为指标 的规范化值；为指标 的权重值（），通过组合权重法得出，之和为；为评价指标的个数。干旱灾害风险模型干旱灾害风险系统由 部分构成，包括致灾因子危险性（）、孕灾环境敏感性（）、承灾体易损性（）和防灾减灾能力（），其表达式如下：（）式中：为干旱灾害风险指数，表示旱灾风险大小；、是 个因子的权重指数。从公式中可以看出，高防灾能力会降低干旱风险，而其他因子会增加干旱风险。干燥度指数干燥度指数是指潜在蒸发量与降水量的比值，反映了当地的干燥程度，干燥度指数越大，说明该地区越容易发生干旱，其公式如下：（）式中：为干燥度指数；为潜在蒸散量，；为降水量，。潜在蒸散量计算常用的方法主要有 方法和 方法。本文中各气象站点的潜在蒸散量根据联合国粮农组织（）推荐的 模型计算，该模型包含的气象要素较多，计算的误差较小，其公式如下：()()()（）式中：为潜在蒸散量，；为地表净辐射，（）；为土壤热通量，（）；为干湿表常数，；为日平均温度，；为 高处风速，；为饱和水气压，；为实际水气压，；为饱和水汽压曲线斜率，。第 期 李 臻，等：山东省干旱灾害风险评估与区划研究 指标体系构建及权重确定 指标体系构建由于干旱形成的具体原因和过程十分复杂，涉及各种气象因素、社会经济因素、地形因素等方面，所以在干旱指标的选取上也会受到时间和空间尺度的诸多限制。参照国内外干旱灾害的研究情况，同时考虑数据的可获取性和可用性，以及数据是否易于处理和表达，本文筛选出构成风险要素的若干评价指标，最终确定了干燥度指数、气温、日照等 个指标，图 为建立的山东省干旱灾害风险评估指标体系。图 山东省干旱灾害风险评价指标体系 基于波段集统计方法的指标间一致性分析 利用波段集统计工具对选取的评估指标间的一致性进行分析，得到指标间的相关矩阵，表 中的参数名称如图 所示。从表 中可以看出，大部分参数之间的相关性低于 ，相关性较弱，只有气温与日照数、高程和坡度、人口密度与经济密度之间有较强的相关性，这表明所选参数间的总体相关性较弱，能够较为全面地反映干旱灾害风险评估所需要表达的指标。评价指标权重的确定采用层次分析法和熵权法确定 个影响因子及其子指标的权重，层次分析法构造的判断矩阵均通过了一致性检验，再利用线性组合法得出各指标的组合权重，具体指标和权重如表 所列，其中权重值的大小表征评价因子对干旱灾害的影响程度，权重值越大影响越大。干旱风险区划与分析 致灾因子危险性区划致灾因子是引起灾害发生的最重要的影响因素。本文致灾因子危险性选取气温、日照以及干燥度指数等指标，这些因素控制着干旱灾害的分布，决定了干旱灾害的变化趋势。根据山东省 个气象站 年的逐日气象数据，利用彭曼公式（），采用 软件计算潜在蒸散量，然后计算各站点的干燥度指数。采用反距离权重（）插值法对气象站点进行空间插值，得到气温、日照、降水、干燥度指数等指标的多年平均空间分布栅格图（见图）。山东省年平均降水量由东南地区向西北地区减少，其中鲁南地区降水量可达 ，是降雨量最多的地区，而干燥度指数的空间分布特征与降水量的分布相反，鲁北日照时间长，鲁西南气温高，气温在 之间。对各指标进行规范化处理后，根据公式（）建立风险分析模型，并利用自然断点分级法将致灾因子危险性指数（）划分为 个等级：为低危险区，表 干旱灾害风险评估参数间相关性 相关性 人 民 长 江 年表 干旱灾害风险评估指标体系 目标层因子指标权重层次分析法熵权法组合法干旱灾害风险指数 致灾因子危险性（）干燥度指数 气温 日照 孕灾环境敏感性（）高程 坡度 水系 土地利用类型 承灾体易损性（）人口密度 经济密度 耕地率 粮食单产 农林牧渔业产值密度 地均大牲畜 防灾减灾能力（）人均 有效灌溉率 图 致灾因子指标分布 为次低危险区，为中等危险区，为次高危险区，为高危险区。从图 可以看出，致灾因子危险性大体趋势是西北高、东南低。高危险区和次高危险区主要位于山东省西北地区，包括平原县、禹城市、宁津县、利津县、沾化区等地方，低危险区和次低危险区主要分布在威海市、日照市、临沂市、青岛市南部、枣庄市和烟台市南部，其他地方为中等危险区。主要由于鲁西北降水量少、蒸发量大、气温较高，成为高危险区，鲁南及东部沿海地区降水量多、蒸发量小、气候适宜，危险性低。孕灾环境敏感性区划孕灾环境包括当地的地形条件、植被状态、水文条件以及土壤条件等。在本文中，孕灾环境指标重点考虑地形、水系、土地利用类型等 个方面对干旱灾害所产生的综合影响。地形因素主要包括高程和坡度，坡图 山东省干旱灾害致灾因子危险性分布 度基于 提取。地势较高、坡度越陡的地带不容易蓄积足够的水分，容易发生旱灾。水系由河网密度表示，河网密度可以反映流域内一定面积范围内的河流密集程度，本研究中，利用 中的 模块计算了河网密度。河流、支流多且密集，水资源丰富，河网密度越高，干旱发生的概率就越小。不同的土地利用类