基于GIS的天水市滑坡危险性评价.doc

 甘 肃 地 质 2008 第 17 卷 第 4 期 2008 Vol.17 No.4 GANSU GEOLOGY 文 章 编 号 ：1004-4116(2008)04-0072-0004 基于GIS的天水市滑坡危险性评价 白晓华1，马金辉2，冯乐涛1，韩庆宪1 （1.甘肃省科学院地质自然灾害防治研究所 ，甘 肃 兰 州 730000；2.兰州大学资源环境学院 ，甘 肃 兰 州 730000） 摘 要 ：本文通过统计学方法提取可能影响该地区滑坡的因子 ，并依据二元逻辑回归结果 ，利 用GIS的 空 间 分 析 和 建 模 功 能 ，对研究区滑坡进行建模 ，再由相关性等级分析方法 进一步获取相关影响因子对滑坡的影响范围 ，最 终 得 出该区域滑坡危险性评价图 。 而在不同分辨率尺度上 ，对滑坡产生影响的因子有所相同 ，当 分 辨 率 由1 000m~60m时 ， 影响因子明显增多 。以逻辑回归模型自身精度 、滑坡实际发生比和坡度法为模型判定标准 ，发 现120m分辨率下的滑坡 发生概率模型能有效表达研究区滑坡发生状况 。 关 键 词 ：滑 坡 ；GIS；危 险 性 评 价 中 图 分 类 号 ：P642.22 文 献 标 识 码 ：B 通 过 ArcGIS 9 空 间 分析模块中的地形栅格转换插 值（Topo to Raster）方法得到的[2]。 2．1．3 滑坡影响因子 本文在现有资料基础之上 ， 主要针对影响滑坡 的内在因素 ， 选取出了与该区域滑坡可能相关的地 质岩性和地形因子 ，各影响因素选取如下 ： 1 引言 天水市位于陇西黄土高原东南部 ， 渭 河 及 其 支 流 河 谷 地 带 ，市 区 东 西 长 30km，南 北 宽 2~3km，呈 条 带状展布 ， 南北两侧黄土山梁连绵起伏 ， 是 我 国 滑 坡、泥石流等地质灾害多发的城市之一 ，也 是 我 省 地 质灾害重灾区之一[1]。 长期以来 ，许多科研机构 都 不 断引进新技术和新方法来拓展滑坡研究 手 段 。 我 们 利用统计学方法 ， 确定出该区域滑坡最 为 相 关 的 影 响因子 ，并采用GIS空间分析和环境建模 方 法 对 天 水 地区滑坡危险性进行分区评价 。 表 1 滑坡影响因子表 Table 1 The landslide factors 地 质 因 子 地 质 类 型 、岩 土 类 型 高 程 、坡 向 、坡 度 、坡 长 曲 率 （平 面 曲 率 、坡面曲率和下流域曲率 ） 河 流 等 级 、流 长 、上 集 水 区 面 积 、径 流 能 量 地 形 指 数 、地 表 粗 糙 度 、地 形 湿 度 地 形 因 子 （由 DEM 数 据 派 生 而 成 ） 2 资料与方法 2．2 研究尺度的选择 为了探讨滑坡影响因子在不同尺度下的差异性 和寻求适合于该区域的最佳研究尺度 ， 本 文 基 于 研 究区滑坡发生的规模 、 现有的数据资料状况和区域 环 境 特 性 ， 采 用60m、120m、250m、500m和1 000m 5 个分辨率尺度 ，对研究区域进行滑坡建模 。 由于滑坡数据最初在 GIS中 是 以 点 数 据 格 式 所 表 示 ， 因此在处理过程中须将其分别转化为上述 5 2．1 数据准备 2．1．1 滑坡数据 滑坡数据来自于天水市滑坡普查数 据 。 为 了 符 合滑坡建模的数据需 求 ， 在ArcGIS中将其进行面域 化处理 ，最终得到研究区的滑坡栅格数据 。 2．1．2 DEM数据 研究区的DEM数据是由1：5万地形图数字化后 ， 收 稿 日 期:2008-09-25 基 金 项 目:甘肃省科学事业费科研项目 （QS031－C32－02）资 助 。 作 者 简 介:白 晓 华 （1965~），女 ，1988年毕业于西北大学地理系 ，副 研 究 员 ，主 要 从 事GIS与地质灾害防治研究 。 个分辨率尺度下的ESRI GRID格 式 ， 同 时 进 行 重 采 样 ，对 有 滑坡的区域采样为 1，无滑坡的区域采样为 0。 3 结果 3．1 采用逻辑回归计算滑坡概率 对研究区滑坡数据而言 ， 在一组自变量作用下 所发生的结果用指示变量Y表 示 ，对 于Y而 言 其 赋 值 规则为 2.3 数据处理总体流程 本文以天水市滑坡区域为研究区 ，以GIS空 间 分 析和环境建模功能为技术支持 ， 首先提 取 与 该 区 域 滑坡可能相关的多个影响因子 ， 然后利 用 统 计 学 方 法 确 定 出 对该区域滑坡最具影响力的影响因子 ，并 在此基础上对该区域滑坡进行建模 ， 最 终 得 到 该 区 域滑坡危险