博士学位论文答辩报告博士学位论文答辩报告 作者：张俊香作者：张俊香 导师：黄崇福导师：黄崇福 教授教授 专业：自然地理学专业：自然地理学 自然灾害与风险管理自然灾害与风险管理2005年5月26日北京师范大学资源学院College of Resources Science and Technology,Beijing Normal University论文题目新一代自然灾害风险区划新一代自然灾害风险区划原理与方法研究原理与方法研究以地震灾害为例以地震灾害为例选题原因选题原因v国内外的减灾实践表明，用风险的理念认识和管理灾害，才能在最大程度减轻灾害影响的同时，谋求社会经济的持续发展。v人们永远不可能准确估计任何自然灾害的风险。“风险值估计不准”是自然灾害风险区划的一大难题。如何在风险区划图中将“风险值估计不准”这一信息揭示出来，让用户在决策分析时清楚何处数据有质量问题及其严重程度，是新的风险区划图图型设计时应关注的问题。主要内容1风险区划研究现状与展望2基于信息扩散技术的模糊风险理论与方法综述3软风险区划图模型设计4应用实例5软风险区划图与概率风险区划图的比较6结论和讨论1 风险区划研究现状与展望1.1 研究现状1.2 存在问题1.3 发展趋势1.1 研究现状自然灾害区划: 自然灾害区划图，实际上是自然灾害分布图，以自然灾害分布的空间特点进行区划。自然灾害风险区划随着对自然灾害研究的深入，逐步将灾害的量级概念引入区划图，以某一地区可能发生的灾害最大等级划分区划等级.进一步的研究，又将概率概念引入区划图，以某一地区可能发生的灾害最大等级的概率或超越某一概率的灾害最大等级为基础编制灾害区划图。1.2 现有风险区划图及其不足图4 某个研究区 图5 极值风险分析及极值风险区划示意图 令 ; 是 的估计值。 极值风险区划不足： 极值风险区划对风险的描述采用确定性的 方法，没有顾及到灾害发生的不确定性。概率风险区划图6 超越概率风险分析及超越概率风险区划示意图不足：表达了风险的随机不确定性，而忽视了模糊不确定性（小样本情况下）。风险软计算模糊风险风险区划软风险区划1.3 发展趋势2 基于信息扩散技术的模糊风险理论与方法综述2.1 自然灾害风险系统中的不确定性2.2 提高风险估计精度的信息扩散技术2.3 自然灾害模糊风险2.4 内集外集模型2.1 自然灾害风险系统中的不确定性自然灾害系统中的不确定性信息扩散技术随机不确定性模糊不确定性概率论2.2 提高风险估计精度的信息扩散技术信息扩散技术原理：令X是一个给定的样本，假定用它可以估计一个在论域U 上的关系。当且仅当X不完备时，必定存在一个适当的扩散函数和相应的算子，使得扩散估计比非扩散估计更靠近真实关系。目的：挖掘出尽可能多的有用信息，提高系统识别精度。最简单的方法是信息分配方法，最简单的扩散函数是正态扩散函数。2.3 自然灾害模糊风险 设I=Ij是自然灾害论域，P=pk是概率论域，为Ij的概率是pk的可能性分布， 是一个可能性概率分布（Possibility-Probability Distribution, PPD），称作I的模糊风险（Fuzzy Risk）。 2.4 内集外集模型.ujuiukuluhDxj1 xj2 xj3 xk1 xk2 .xi1xi2xi3qj2,j1-qj1,j令 Qj =def (1-qj1,j), (1-qj2,j), (1-qj3,j) sort(Qj) = (1-qj2,j), (1-qj3,j), (1-qj1,j)_令 Qj =def qi1,j, qi2,j, qi3,j, qk1,j, . sort(Qj) = qi3,j, qk1,j, qi2,j, qi1,j, . +|Qj | = nj |Qj | = n - nj_+2.4 内集外集模型1st (smallest) element of p = p02nd element ofp = p1.last (largest) element ofp = pnj 11 p = pnj 21st (largest) element ofp = pnj + 132nd element of p = pnj + 24.5last element of p = pnQj+Qj_Qj_Qj_Qj+Qj+pj(p) =令 pk = k/n , k = 0, 1, ., n表4 可能性概率分布表达的模糊风险 模糊风险相对于概率风险的主要区别在于：对超越概率引入可能性分布，指明了风险事件以某概率出现的可能性，从而既体现了超越概率估算并不精确这一事实，又为模型容纳模糊信息提供了一条途径。3 软风险区划图模型设计3.1 不确定性可视化的相关技术分析3.2 可能性概率分布的模糊期望值3.3 软风险区划图的总思路和技术途径3.4 软风险区划图模型3.1 不确定性可视化的相关技术分析超越概率区划图技术GIS制图技术 模糊集截集技术3.2 可能性概率分布的模糊期望值图8 可能性概率分布图 ，然后将隶属度 的元素即 挑选出来。令称有限闭区间 是 关于 的 水平截集。 (4-2)(4-3)可能性概率分布的模糊集截集示意图 令 称 和 分别是 和 的归一化。 令称 是 关于 的 水平截集的期望区间。 (4-5)(4-6)(4-7)表5 某一地区的灾害模糊风险值 3.3 软风险区划图的总思路和技术途径总思路： 以自然灾害模糊风险为约束条件，综合模糊集截集技术、GIS技术和超越概率区划图技术在表达不确定性方面的优势，给出既能指明风险值可信度的信息，又包含传统风险图的信息的风险图图型模式。 软风险区划图编制技术方案历史灾情分析历史灾情分析建立历史灾害数据库建立历史灾害数据库灾情分布图灾情分布图灾害强度分析灾害强度分析孕灾环境分析孕灾环境分析风险区划基本单元风险区划基本单元编制软风险区划图编制软风险区划图用户需求用户需求空间叠加空间叠加模糊风险计算模糊风险计算设防标准研究设防标准研究现有风险图基础数据现有风险图基础数据区划图应用政策研究区划图应用政策研究选取研究区选取研究区图10 可能性概率风险分析与软风险区划图3.4 软风险区划图模型4 应用实例4.1 研究区概况4.2 地震模糊风险计算4.3 模糊期望值的计算4.4 软风险区划图的绘制4.1 研究区概况斜线表示的区域为新生代印度板块与欧亚板块碰撞构造应力场影响的主要区域 图11 云南所处的区域构造环境图12 云南区域的活动构造与块体划分图13 云南地区M5.0级地震的震中分布图14研究区示意图 资料来源基础资料 ：（1）历史地震目录资料（包括发生时间、空间（地表坐标、 震中位置）和强度（震级）三个基本参数 ）。（2）地质资料：活动构造分布图。来源： 中国历史强震目录（国家地震局震害防御司， 1995）、中国近代地震目录（中国地震局震害防御司，1999）和中国地震年鉴（19922001）（国家地震局，19922001）。（由云南省地震局提供） 图15 研究区M5.0地震震中分布与地震活动带（区）划分 表7 各地震带（区）历史强震震级记录样本 4.2 地震模糊风险计算由于我们不知道该地区强震震级所服从概率分布的函数形式（例如，是正态分布还是指数分布），而且各地震带（区）强震记录样本的容量都很小，分别为14、7、13、11、6。因此，依据所给样本进行的概率估计必然很不精确，我们将用可能性概率分布来表达该地区的强震风险。计算地震模糊风险的基本步骤地震区的样本容量仅为14，是小样本，我们考虑用可能性概率来表达该地区的强震风险。考虑区间 相应的离散论域为： 此处控制点步长为 =0.6。 第1步：由式（3-4）计算各样本点的分配信息，列于表8中。表8 各样本点的分配信息第2步：由式（3-6）计算游离信息量 ，由表9所示。对给定的j， 升序排列 ，得游离信息升序 ，如表10所示。 表9 游离信息表10 游离信息升序 第3步：由式（3-7）计算漂入信息量 ，由表11所示。对给定的j，降序排列 ，得漂入信息降序 ，如表12所示。 第4步：令： 将表10和表12的有关数值代入式（3-9），可得由可能性概率分布表达的地震模糊风险如表13所示。 表13 可能性概率分布表达的地震区的地震灾害模糊风险I1I2I3I400.330.1700.170.070.330.330.170.170.1410.330.330.170.210.50.330.510.290.330.510.330.360.3310.170.170.430.330.50.1700.500.170.330.1700.5700.33000.6400000.7100000.7900000.8600000.9300001.00000表14 地震带的地震灾害模糊风险 表15 地震区的地震灾害模糊风险 表16 地震带的地震灾害模糊风险 表17 地震区的地震灾害模糊风险 表18 各地震带（区）地震灾害风险的模糊期望值 4.3 模糊期望值的计算4.4 软风险区划图的绘制图16 昆明楚雄大理丽江地区地震保守风险图（ 0.25） 图17 昆明楚雄大理丽江地区地震冒险风险图（ 0.25）图18 昆明楚雄大理丽江地区地震保守风险图（ 0.5） 图19 昆明楚雄大理丽江地区地震冒险风险图（ 0.5） 5 软风险区划与概率风险区划的比较5.1 研究区的概率风险区划图5.2 研究区的地震动参数区划图5.3 软风险区划与概率风险区划的比较5.1 研究区的概率风险区划图直方图模型：称 是概率密度函数的直方图估计。表20 昆明楚雄大理丽江地区的概率风险（超越概率10） 图20昆明楚雄大理丽江地区概率风险区划图（100年超越概率10） 图21 云南省有效峰值加速度区划图（50年超越概率10） 5.2 地震动参数区划图5.3 软风险区划图与概率风险区划图的比较区划图的技术方法 概率风险区划图：软风险区划图概率风险计算理论区划图技术模糊风险计算理论区划图技术概率风险区划： 地质资料（断层分布）、地球物理场资料、工程地震资料、强震资料地震软风险区划 强震资料、地质资料（断层分布）采用资料 区划图的结果 地震动峰值加速度区划图（50年超越概率10）地震动反应谱特征周期区划图（50年超越概率10）概率风险区划软风险区划采用模糊风险分析方法；采用“层”结构的图型模式，是一种多值风险区划图；充分考虑到不同用户的需求，体现了“人性化”的服务趋向。软风险区划图的特点 6 结论和讨论6.1 结论6.2主要的创新点6.3讨论6.1 结论“风险值估不准”是现有超越概率风险区划图的不足。如何表达风险值估计不准是深层次风险信息可视化研究的问题。GIS技术、超越概率区划图技术、模糊集截集技术在表达不确定现象方面各有其长处，为“风险值估不准”条件下的自然灾害风险信息可视化研究提供了重要的技术和方法。 自然灾害模糊风险计算理论是自然灾害风险评估研究不断深入的必然，为新一代自然灾害风险区划理论与方法研究提供了科学依据。 结论云南省昆明楚雄大理丽江地区的地震灾害软风险区划图的实证分析说明，自然灾害软风险区划图是一种多值风险区划图。其中可能性水平代表了风险信息的可靠性水平，多值风险表达了风险值估计不准这一信息。自然灾害软风险区划是对自然灾害风险信息的可喜的图型展示，是对现有风险区划图的一次重大改进，能为风险决策者提供更多风险信息，提高风险决策的科学性和有用性，是未来灾害区划与风险区划的发展方向。 6.2 主要的创新点首次针对现有自然灾害风险区划图“风险值估不准”这一致命弱点,引入处理小样本的信息扩散技术和相关的风险分析模型用于自然灾害风险区划研究，建立风险值估不准条件下的自然灾害风险区划的理论和方法综合模糊集截集技术、GIS技术和超越概率区划图技术在表达不确定性方面的优势，给出了既包含传统风险区划图信息，又能指明风险值可靠度的软风险区划图模型。首次将软风险区划理论与方法应用到我国多震省之一的云南地区。6.3 讨论自然灾害软风险区划图研究尚处于初步探索阶段，目前很多问题还没有搞清楚。例如：如何选取可能性水平；软风险区划图使用年限的确定问题；自然灾害软风险区划在区域经济可持续发展以及防灾减灾等方面的潜在价值分析。 6.3 讨论本文只是自然灾害软风险区划图研究的一个开端，重在理论的创新和思路的完整，与现有成熟的地震风险区划理论和方法相比，尚待成长。如能将已有的研究成果纳入这一新领域（例如随机过程的认识，地质资料的使用，地球物理分析方法等），则软风险区划图将会有更强的生命力。 演讲完毕，谢谢观看！