基于广义线性模型的地质灾害气象预警模型研究.pdfVIP

基于广义线性模型的地质灾害气象预警模型 徐 晶 国家气象中心天气与环境预报室 北京 100081 摘要 利用1949年-2003年的地质灾害 气象资料及遥感地质数据 基于以往的降水诱发地质灾害的统计模 型研究的基础上 本文尝试将传统的降水观测与遥感地质数据结合起来 运用统计学中已获得广泛应用的 广义线性模型建立起地质灾害模型 试验表明将气象因素与地质环境因素综合考虑纳入模型是可行的 为 提高地质灾害气象预警水平提供的有效的途径 而广义线性模型则是实现这一构想的有力工具 同时在2004 年的业务试验也证明该模型具有一定的实用价值 关键词 广义线性模型 地质灾害 降水 1 前言 降雨是滑坡 泥石流等地质灾害的主要诱发因素 由局地暴雨引发的滑坡 泥石流等灾 害占这类灾害总数的90%和95%以上 因此 根据降雨的观测进行滑坡 泥石流等地质灾害 的预报 成为地质灾害预报的主要方法[1]-[4] 2003年 11月 19日国务院第29次常务会议通 过了 地质灾害防治条例 以下简称 条例 2003年 11月24日温家宝总理签发国务院 第394号令 予以公布 自2004年3月 1日起施行 国土资源部和中国气象局根据 条例 的要求 决定联合开展地质灾害气象预报警报工作 并委托中央气象台和中国地质环境监测 院联合攻关 研究开发全国地质灾害气象预报模式 初步建立了国家级地质灾害气象预报预 警系统 并于2003年6月 1日正式运行 从此 我国地质灾害防治工作走上了法制化 规 范化的轨道 用降雨量预报泥石流等地质灾害的基本思路是 通过对雨量资料进行统计分析 确定滑 坡 泥石流等地质灾害的临界雨量或触发雨量 然后根据降雨量制作灾害发生可能性的预报 [5]-[13] 然而 滑坡 泥石流等灾害的发生除了直接受降水因素触发之外 还受地形地貌 土壤 植被等多种地质和生态环境的影响 以往的预报模型中 由于缺乏必要的地质环境数 据 往往只考虑降水因素 因子选择过于单一 不能全面反映诱发灾害的地质环境条件 另 外 建立模型时不够注重模型的统计检验 缺乏理论上的严格性 近年来 随着遥感技术的 应用 高分辨率的地质环境卫星遥感数据已经可以获得 统计科学中一些先进的方法也已在 地球科学的各领域得到广泛应用 因此 本文尝试将传统的降水观测与遥感地质数据结合起 来 运用统计学中已获得广泛应用的广义线性模型 建立更为全面 严格的地质灾害气象预 警模型 2 资料来源 本文使用的资料主要包括地质灾害资料 降雨量资料以及卫星遥感地质数据 地质灾害资料由国土资源部中国地质环境监测院提供 该资料的时间跨度从 1951年到 2002年 共 13400多条地质灾害资料 这部分资料分为基本历史资料和分县调查历史资料 另外 我们也从全国气象业务服务信息系统上收集到了2004年发生的滑坡 泥石流灾 害资料 用于预报效果 灾害在地质灾害气象预报方法研究中只使用 1951-2002年的资料 我们还使用了相应年份的降雨量资料 该资料来自国家气象中心降雨量历史资料库 地 质遥感数据由气象中心遥感业务室提供 包括地形坡度 构造因子 土壤类型 植被覆盖等 多个因子的分辨率为1公里的数据 3 广义线性模型 [14] 本文采用广义线性模型 中的Logistic模型来建立地质灾害的发生概率模型 传统的 统计预报采用的线性回归模型为 设Y (Y , , Y ) 为一n 1 随机向量 每个分量依赖于 1 n p 个预测因子 构成一n p 矩阵X 其元素X 对应于第i 个随机变量Y 的第j 个预测因 ij i 子 假设Y 遵从某种概率分布 其均值向量为ì (m , , m ) 则有 1 n