地理科学进展  2016 , 35 (2): 137-147 https://doi.org/10.18306/dlkxjz.2016.02.001

本刊特稿

山洪泥石流风险评估与风险管理理论与方法

崔鹏13, 邹强2

1. 中国科学院山地灾害与地表过程重点实验室/中国科学院成都山地灾害与环境研究所,成都 610041
2. 西南科技大学环境与资源学院,四川 绵阳 621010
3. 中国科学院青藏高原地球科学卓越创新中心,北京 100101

Theory and method of risk assessment and risk management of debris flows and flash floods

CUI Peng13, ZOU Qiang2

1. Key Laboratory of Mountain Hazards and Earth Surface Process/Institute of Mountain Hazards and Environment, CAS, Chengdu 610041, China
2. School of Environmental and Resource, Southwest University of Science and Technology, Mianyang 621010, Sichuan, China
3. Center for Excellence in Tibetan Plateau Earth Sciences, CAS, Beijing 100101, China

收稿日期: 2015-12-10

接受日期:  2016-01-10

网络出版日期:  2016-02-10

版权声明:  2016 地理科学进展 《地理科学进展》杂志 版权所有

基金资助:  国家自然科学基金国际合作与交流重点项目(41520104002)中国科学院重点部署项目(KZZD-EW-05-01)四川省科技计划项目(2014SZ0163)

作者简介:

作者简介:崔鹏(1957-),男,陕西西安人,中国科学院院士,主要从事山地灾害与水土保持研究,E-mail:pengcui@imde.ac.cn

展开

摘要

山洪泥石流是中国常见的自然灾害,充分认识其形成机制与潜在风险是防灾减灾的关键。本文阐述了山洪泥石流形成机理,以及风险分析与管理的方法和内容,系统认识了地表产流流量激增、土体破坏物质供给激增、沟道堵塞体级联溃决流量放大和动床侵蚀规模增大等4个山洪泥石流的形成过程,介绍了基于动力过程的山洪泥石流风险评估方法和承灾体易损性评估方法,构建了基于灾害动力过程的风险评估与风险制图方法。进而,基于风险评估结果,提出可用于具体灾害点减灾的风险管理内容和风险调控技术、灾害防治的工程与非工程措施与制技术方案。最后,重点讨论了包括灾害风险预测、临灾预案、灾害防治工程方案等内容的风险处置对策,并形成一套基于山洪泥石流动力过程的风险评估与风险管理理论与方法体系。

关键词: 泥石流 ; 山洪 ; 危险性 ; 易损性 ; 风险评估 ; 风险管理 ; 动力过程

Abstract

Debris flows and flash floods are widely distributed mountain hazards in China. Effective hazard mitigation and prevention require understanding of hazards formation mechanisms and their potential risks. This article elaborates on the formation mechanism, risk analysis, and risk management of debris flows and flash floods. Surface runoff and material supply volume incensement, hazard scale amplification due to outburst of multiple channel blockages and bed erosion as hazard formation mechanism are discussed. Base on the dynamic process of debris flows and flash floods as well as vulnerability assessment of elements at risk, methods of risk assessment and mapping are proposed. Comprehensive engineering and non-engineering measures for hazards control shall be guided by the result of risk analysis that identifies the hazardous level of debris flows and flash floods and incorporates the vulnerability of different elements at risk. Finally, this article discusses strategies when facing risk of these hazards and puts forward a risk management system that involves the participation of local communities.

Keywords: debris flow ; flash flood ; hazard analysis ; vulnerability ; risk assessment ; risk management ; dynamic process

0

PDF (1089KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

崔鹏, 邹强. 山洪泥石流风险评估与风险管理理论与方法[J]. , 2016, 35(2): 137-147 https://doi.org/10.18306/dlkxjz.2016.02.001

CUI Peng, ZOU Qiang. Theory and method of risk assessment and risk management of debris flows and flash floods[J]. 地理科学进展, 2016, 35(2): 137-147 https://doi.org/10.18306/dlkxjz.2016.02.001

1 引言

山洪和泥石流是山区小流域水土物质快速输移过程,暴发突然,破坏力强,往往造成毁灭性灾害。活跃的地质构造、巨大的地形高差、季风气候带来的丰沛降雨、密集的人口分布和高强度的人类活动影响等因素的组合与叠加,使得山洪和泥石流广泛发育于约占中国陆地面积2/3的山地,频繁成灾,损失巨大,严重威胁广大山区人民生命财产和重大工程安全,制约社会经济发展。

根据国家防洪抗旱总指挥部统计,全国2058个县(市)有山洪和泥石流分布,分布区面积达487万km2,受威胁人口达到5.7亿人。2010年8月8日甘肃省舟曲山洪泥石流就造成1765人死亡和失踪,4321 间房屋毁坏,22667人无家可归;2013年7月10日前后的降雨过程导致都江堰、汶川、北川等10个县(市)暴发大范围群发性山洪和泥石流,导致2.1×104 km2范围内269个乡镇、118.7×104人受灾,毁坏房屋6975间,直接经济损失114亿元,共造成100余人死亡;2013年8月16日辽宁清原特大型山洪灾害造成77人死亡,直接经济损失76.34亿元。据涂勇等人统计(2013),2013年全国24个省(自治区、直辖市)共发生了181起造成人员死亡的山洪灾害事件,死亡774人,山洪灾害死亡人数占洪涝灾害死亡人数比例为72.4%;近十多年来,全国因山洪和泥石流灾害造成的死亡人数占洪涝灾害死亡人数的比例呈增长趋势,占比高达70%~90%(图1),成为造成人员伤亡的主要灾种。

图1   2000-2014年山洪灾害死亡人数及其占洪涝灾害总死亡人数比例(数据来源:国家防汛抗旱总指挥部等, 2013; 张志彤, 2015)

Fig.1   Casualties from flash flood disasters and their proportions in annual total casualties due to floods, 2000-2014(Data source:The State Flood Control and Drought Relief Headquarters et al, 2013; Zhang, 2014)

山洪泥石流灾情表明,预防和减轻山洪泥石流是国家减灾战略的重要部分。从国际上看,联合国减灾署2015年3月在日本仙台召开的第三届世界减灾大会上,通过了《2015-2030年仙台减少灾害风险框架》,力求在未来15年内大幅减少在生命、生计方面以及在经济、实物、社会、文化、环境资产方面的灾害风险和损失,采取包括经济、法律、社会、卫生、文化、环境、技术、政治和体制等综合措施,减少受灾体暴露程度和脆弱性,提高抗灾和灾后恢复能力(UN/ISDR, 2015)。新的国际减灾行动计划必将有力地促进山洪、泥石流灾害风险研究与减灾工作。

灾害风险评估是减灾的基础和依据。自然灾害风险评估旨在预测灾害期望损失值和评估灾害损失程度,不仅是减轻自然灾害损失的重要途径(UNDP, 2004; UN/ISDR, 2005),也为土地利用规划提供科学依据(AGS, 2007)。山洪和泥石流风险评估分为区域评估与灾点评估。区域评估主要依据灾害成因和承灾体的数量、质量与空间位置等信息,采用统计分析方法,分别确定灾害的危险度与成灾体的易损度,进而计算出风险值。区域评估的结果可以用于不同级别行政单元甚至全国和全球的灾害风险管理,服务于不同级别的减灾规划和减灾方针政策制定。这方面的研究已经取得了突出进展,在国家层面减灾中发挥了重要作用(Huggel et al, 2003; Liu, 2003; Chau, 2004; Bisson et al, 2005; 史培军, 2005; Pradhan, 2010)。近年来考虑到承灾体的暴露性,对风险评估方法进行了新的改进,使得风险评估结果更加符合实际(葛全胜等, 2008; Fuchs, 2012, Cui, Zou et al, 2013; 邹强, 2014)。

由于区域风险评估是基于统计方法的因素分析,对具体灾点的详细信息反映不足,评估结果难以应用到一个具体灾害点的临灾预案和减灾工程设计。对服务于具体灾点减灾工作的风险评估,需要进一步依据具体灾害的物理力学性质,考虑灾害的动力过程,开展灾害危险性分析与易损性分析,在此基础上计算该灾点的风险值并进行风险等级划分。我们把这种基于灾害体动力过程和灾害体与承灾结构体动力响应的风险分析方法,称为基于动力过程的风险评估。本文将介绍基于动力过程的风险评估方法以及在此基础上提出的社区灾害风险管理。

2 基于动力过程的风险评估

基于动力过程的灾害风险评估的核心科学问题,是基于灾害动力过程认识进行危险性分析和易损性分析。

2.1 山洪泥石流动力过程

山洪和泥石流都是水沙混合体在沟道中的快速运动过程。其形成与发展主要包括以下四个过程:①降雨在集水区产流并迅速汇聚形成沟道洪水;②水流冲蚀裹挟沟道内的松散固体物质形成泥石流;③在其沿着具有滑坡崩塌体堵塞段沟道的流动过程中,存在物质和能量暂时蓄积与瞬间释放过程,使得洪水或泥石流规模增大(Cui, Zhou et al, 2013a);④山洪和泥石流在较陡的沟道中高速运动,具有强烈的侵蚀冲刷能力和泥沙携带能力,侵蚀沟床并携带松散物质,不断补给物质而增大规模,随着物质和能量不断积累,最后形成具有较大破坏能力的洪流。山洪和泥石流形成发展的四个关键过程为:地表产流流量激增过程,土体破坏物质供给激增过程,沟道堵塞体级联溃决流量放大过程和动床侵蚀规模增大过程(图2)。

图2   山洪泥石流形成过程示意图

Fig.2   Formation process for torrential floods and debris flows

山洪和泥石流的破坏动力作用主要表现在高速运动的冲击作用、快速增大水位(流深)的洪水淹没和泥沙淤埋作用,还有对构筑物基础的冲刷和淘蚀作用。其中,冲击和淤埋(淹没)为最常见的破坏方式。

2.2 山洪泥石流数值模拟

流速是确定山洪和泥石流冲击破坏作用的关键参数,流深反映其淤埋或淹没程度。依据山洪和泥石流运动方程,通过动力学数值模拟方法,确定山洪与泥石流的危害范围、以及致险范围内流速和流深的时空分布,从而可依据动力参数确定其危险性。

山洪一般可视为清水或含沙水流,浅水波方程可以较好描述其运动过程。目前,基于运动波和扩散波近似的分布式水文模型能够预测山洪区域的内部水位、流速过程。但由于忽略了流动的力学因子,难以准确描述山洪的形成与演化过程。因此有必要发展完整的二维水动力学模型。传统数值计算格式难以准确计算激波等水流现象,曹志先团队等建立了新的多时间尺度山洪水沙数值模型,开辟了山洪数值模拟的新途径(王鑫等, 2009)。该方法控制方程采用二维浅水波方程,包括质量守恒方程和动量守恒方程,考虑降雨和土壤入渗,构成如下守恒形式:

Ut+F(U)x+G(U)y=S(U)(1)

式中:U为守恒变量;F(U)、G(U)为通量向量;S(U)为源项向量。方程中各项的定义如下:

U=hhuhv,F=huhu2+gh2/2huv,G=hvhuvhv2+gh2/2,S=r-fgh(Sbx-Sfx)gh(Sby-Sfy)(2)

式中:h为水深(m);(u, υ)为(x, y)方向的流速分量(m/s);g为重力加速度(m/s2);r为降雨强度(mm/h)或入流量(m3);f为入渗强度(mm/h)或出流量(m3);Sbx=-∂z/∂x,Sby=-∂z/∂y,分别表示xy方向的坡面比降(‰);z为地表高程(m);SfxSfy分别为xy方向的阻力坡度(°),采用明渠水流阻力公式;入渗强度采用Green-Ampt公式计算(Smith et al, 2002):

f=dFdt=KS(1+(θS+θI)HcF)(3)

式中:Ks为土壤饱和入渗率(mm/h);Hc为土壤吸力(m); θSθI分别为饱和、初始含水量(%);F为累积入渗深度(m)。

泥石流为多相流体,其动力学过程比较复杂,需要通过假设简化,建立动力学方程,实现数值模拟。通常采用基于Saint-Venant方程的泥石流运动方程,该方程只包含三个重要变量——泥深、x方向的速度、y方向的速度,可以忽略密度的变化和压力的作用。在泥石流数值模拟中,一般只需要考虑泥石流体的动量守恒方程。在Lagrange坐标系中,泥石流运动方程可写成:

DuDt=gSsx-gSfxDvDt=gSsy-gSfy(4)

式中:SsxSsy分别是xy方向泥石流堆积区底面坡降(‰);SfxSfy采用O'Brien等(1993)的摩阻坡降公式:

Sfx=τBγmhsgn(u)+2μBuγmh2+kcuu2+v2ghSfy=τBγmhsgn(v)+2μBvγmh2+kcvu2+v2gh(5)

式中:τB为屈服应力(N/m2);γm为泥石流密度(kg/m3);h为泥石流淤埋深度(m);μB为泥石流粘滞系数(N·s/m2);kc为糙率系数。右边第1项反映了屈服应力对阻力的贡献,第2项反映了粘性应力对阻力的贡献,第3项反映了碰撞应力对阻力的贡献。

模拟计算时把整个区域划分成正方形的网格,对每个交叉点(网格的角点)计算流深、速度值及高程值,作为危险性的评价因子。

2.3 基于动力过程的危险性分析

山洪致灾作用主要有:水位陡涨导致的淹没危害、携带大量泥沙产生的淤埋危害、水流对构筑物基础的冲蚀。泥石流致灾作用主要有对承灾体的冲击、淤埋和冲蚀作用。冲蚀作用是山洪和泥石流作用在具体建筑物上的危害方式,在对具体工程结构体的安全性分析时需要精确计算,在进行不针对具体结构体安全性计算而主要描述灾害的危害能力时,为了简化计算,可以暂时不考虑冲刷作用,只考虑冲击作用、淤埋作用和淹没作用。

(1) 泥石流淤埋深度的确定方法

用泥石流运动过程中的最大流深来表征泥石流的淤积危害。在泥石流运动模拟计算中,每个计算网格(i,j)的流深由网格内的所有颗粒体积除以网格面积得到(Cui et al, 2011),计算公式如下:

Hd=maxt>0(Ni,jΔV)A(6)

式中:Hd为泥石流淤埋深度(m),其值越大,泥石流危险性越大;Ni, j等于以点(i, j)为中心的控制格网内的颗粒数; ΔV为颗粒的体积(m3);A是数值模拟采用的网格面积(m2);t为数值模拟时间(s)。

(2) 泥石流冲击破坏能力的确定方法

用泥石流运动过程中的最大动能来表征泥石流的冲击破坏能力。采用每个网格在泥石流运动全过程中的最大动能值,反映各网格最大冲击作用产生的危险性,计算公式如下:

Ed=Amaxt>0(u2+v2)](7)

式中:Ed为泥石流最大动能;A为数值模拟采用的网格的面积(m2);uv分别是xy方向的速度(m/s);h为泥石流泥深(m);ρ为泥石流密度(kg/m3);t为数值模拟时间(s)。

(3) 洪水淹没破坏能力的确定方法

用洪水运动过程中山洪的最大动量来表示山洪的淹没破坏能力。基于山洪运动过程的数值模拟,获得计算网格单元的山洪流速与水深值,通过式(8)即可计算洪水的最大动量。

Hf=Amaxt>0ρvh(8)

式中:Hf为任一计算网格的动量(m2/s);A为数值模拟采用的网格的面积(m2);ρ为山洪密度(kg/m3);v为任一计算网格洪水流速(m/s);h为任一计算网格洪水淹没深度(m);t为数值模拟时间(s)。

根据承灾体与流深(泥石流淤埋高度和洪水水位)的相对高度,可以确定是否受到山洪或泥石流的威胁及其危险程度。泥石流的危害作用除了冲击和淤埋危害之外,大规模泥石流堵江后还会导致上游淹没和下游溃决洪水灾害,这种灾害链的复合致灾作用主要表现为水的淹没。反映这种淹没的危险性可用山洪的危险度计算方法确定,但在危险性评估时需要与泥石流本身的冲击与淤埋危险性相结合。

基于山洪和泥石流动力过程的数值模拟,可求算泥石流冲毁、淤埋以及洪水淹没的动力学参数量值;结合山洪、泥石流危险性分析方法,利用GIS技术的空间分析技术,进一步确定泥石流冲击破坏能力、泥石流淤埋深度以及洪水淹没破坏能力的空间分布,从而进行危险性分区。考虑灾害的复合致灾作用,一般上游堰塞湖溃决在泥石流堆积之后,在计算溃决洪水淹没范围时,首先利用胡凯衡等(2005)和Cui等(2011)提出的方法,计算泥石流的堆积范围和堆积厚度,再把被泥石流堆积改变以后的地形作为输入条件,计算上游堰塞湖溃决洪水的淹没范围和淹没深度。在进行危险性分区时,对于同时受到泥石流淤埋和洪水淹没的区域,则选取危害最重者作为该区的危险性评估依据。

2.4 基于动力过程的易损性分析

山洪和泥石流灾害易损性分析是对承灾体抗灾能力的综合分析。承灾体易损度大小取决于承灾体对山洪、泥石流作用的敏感程度,通常用承灾体的价值(或数量)及其脆弱性指数来表示。承灾体的易损度可以描述为:

V=V(u)×C(9)

式中:V为承灾体的易损度;V(u)为承灾体的综合价值;C为承灾体的脆弱性指数值。

承灾体的综合价值可以通过承灾体的平均单价P和受灾的实际数量N计算获得。承灾体脆弱性是描述承灾体受到灾害侵袭时损毁的难易程度,可用介于0和1之间的数表示,数值越大越易损毁。承灾体脆弱性的定量描述比较复杂,主要受承灾体自身的结构强度与山洪、泥石流破坏形式的影响。如承灾体相对泥石流位置不同,受损的方式不同,在泥石流主流线上的建筑物和构筑物主要受到冲击危害;位于泥石流边缘或扇形地前缘者,则主要受到淤埋危害。对于被淤埋的承灾体,多数难以再利用,即使结构完好,由于地形改变也很难再发挥原来设计的功效。为此,以房屋和构筑物为例说明承灾体的脆弱性计算方法。

对于受淤埋危害的承灾体,脆弱性指数Cd定义为其几何高度与淤埋深度的比值,即

Cd=HdHc(10)

式中:Hd为泥石流淤埋深度(m);Hc为建筑物或构筑物有效高度(m),如桥面、房屋等的高度。当 HdHc1时,表示建筑物或构筑物已经完全被泥石流淤埋,Cd值取1。

对于受到冲击的承灾体(建筑物或构筑物),其脆弱性指数Ce可通过泥石流最大动能Ed与建筑物或构筑物的结构强度Ec之比来确定,即:

Ce=EdEc(11)

EdEc1时,表示建筑物或构筑物已经完全被泥石流冲毁,Ce值取1。

2.5 风险评估与风险制图

(1) 风险度计算

山洪和泥石流风险评估是对风险的定性或定量计算,其结果可为风险管理提供依据,用于指导灾害风险处置与决策等。在灾害危险性和易损性分析结果的基础上,根据灾害风险的定义(UNDHA, 1992)计算每个评估单元的风险值:

R=H×V(12)

式中:R为评估单元的风险性指标值;H为评估单元的危险性指标值;V为评估单元的易损性指标值。在计算以前,式(12)中各类指数均需经过归一化处理。归一化处理方法如下:

H'=(H-Hmin)(Hmax-Hmin)V'=(V-Vmin)(Vmax-Vmin)(13)

式中: H'是危险度的归一化值,H是危险度指标值,HmaxHmin分别为最大危险度值和最小危险度值; V'为易损度的归一化值,V为易损度指标值,VmaxVmin分别为最大易损度值和最小易损度值。

(2) 风险分级

计算出每个评估单元的风险值以后,还需要对风险度进行等级划分。危险度、易损度和风险度计算值的分级方法通常有等间距法、分位数法、标准差法和自然断点法。其中标准差法适用于符合正态分布规律的数值;而等间距法、分位数法和自然断点法均具有较好制图效果,得到广泛应用。风险分级时,依据灾害风险特征,选取适合的分级方法寻找数据集的自然转折点和特征点对评估结果进行等级划分。结合风险管理的目标,一般将风险度划分为高度风险、中度风险、低度风险和微度风险四个等级,各个等级风险特征描述如表1所示。

表1   山洪泥石流风险分级及特征描述

Tab.1   Risk classification and characteristics of debris flows and flash floods

风险编号风险等级空间分布与危害特征
I微度风险山洪、泥石流分布极少,规模很小,其危险度与易损度均很低,受危害导致损毁的风险很小,一般不影响正常生产和生活。
II低度风险山洪、泥石流分布少,规模不大,承灾体易损度较低,遭受轻度危害,破坏小,灾害综合风险值较低,需要采取一定预防措施。
III中度风险山洪、泥石流分布较广泛,规模较大,承灾体易损度较高,灾害危害较重,灾害风险水平为中等,需布设灾害防治措施以保障生产和生活安全。
IV高度风险山洪、泥石流分布广泛,规模大,破坏力强,危险性与易损度均较高,风险很高,频繁受灾,严重影响生产和生活,应实施综合减灾工程,加强风险管理。

新窗口打开

(3) 风险制图

山洪和泥石流风险制图过程中,风险度分布、风险等级、风险区位置是风险制图中的关键因素。以危险性与易损性的分析结果为基础,应用GIS的地图代数功能进行栅格数值计算,获得研究区风险度分布,用不同颜色表示各风险等级,合并相同风险级别的栅格单元,结合ArcGIS多边形构面方法,完成山洪、泥石流灾害风险分区图。风险制图的流程见图3

图3   山洪泥石流风险制图流程

Fig.3   Flow chart of risk assessment for debris flows and flash floods

2.6 研究实例——清平场镇泥石流风险评估

清平场镇位于四川省绵竹市西北部山区。2010年“8.13”泥石流在清平场镇形成了长3.5 km,宽400~500 m,平均厚度约5 m(最大厚度超13 m),总土石方量约600×104 m3的淤积区(图4)。泥石流造成7人遇难,5人失踪,33人受伤;379户房屋严重受损,占总户数的20.9%,其中,205户房屋全部淹没,75户房屋部分淹没,99户房屋进水。

图4   2010年四川省清平场镇遭受“8·13”特大泥石流危害影像图(影像来源:四川省国土资源厅)

Fig.4   Debris flow disaster in Qingping Township in Sichuan Province on 13 August 2010 (Image from: Land and Resources Department of Sichuan Province)

根据野外勘测与调查,分析泥石流堆积物样品,确定出这些沟谷的泥石流基本物理力学参数,如底床坡度、沟床糙率、流动时间、流团大小、粘滞系数、屈服应力等,结合清平场镇的地形数据,模拟2010年“8·13”泥石流在堆积扇区的运动过程。通过对泥石流运动过程的模拟,可以获取堆积区每个模拟计算网格的流速、流深数据。经过ArcGIS栅格-矢量数据转换即可获取泥石流堆积范围。比较数值模拟的堆积范围与场镇受灾分布,泥石流对场镇实际危害的面积为148.2345×104 m2,而模拟的受泥石流影响范围的面积为143.9488×104 m2,模拟值较实际值小2.9%。

在清平场镇危险性与易损性分析结果的基础上,基于式(13)对危险度和易损度进行归一化处理,利用式(12)的计算方法,通过泥石流危险性归一化数据与承灾体易损性归一化数据进行乘积运算,获取研究区的泥石流风险度。依据风险分级,并用不同图斑表示各风险等级的区域,合并相同风险级别的栅格单元,编制清平镇泥石流灾害风险分区图(图5)。结果表明:泥石流危险区内,高风险区面积约83.7×104 m2,占研究区总面积的33.4%;中风险区面积79.2×104 m2,占研究区面积的31.6%;低风险区面积87.5×104 m2,占研究区总面积35.0%。经野外调查验证,风险分区结果与清平场镇泥石流承灾体的分布情况及危害严重程度比较吻合,评价结果与实际情况基本一致(Cui, Zou et al, 2013b)。

图5   清平场镇风险分区结果

Fig.5   Risk map of Qingping Township

3 风险管理

3.1 山洪泥石流风险管理内容

灾害风险管理(disaster risk management)是在对灾害风险形成过程与机理(致灾体成因、力学机理、运动规律、成灾机制)系统认识的基础上,根据不同类型灾害的物理特性及其成灾特点,评估可能产生的灾害风险,采用适宜的风险调控措施(预案、预防、预警与治理措施),最大限度地减少灾害风险和可能损失的手段。它贯穿于灾害发生发展的全过程,涉及防灾、备灾、预警、响应(应急和救助)、处置(治理)和恢复/重建等内容,是在以往灾害防治理念基础上发展的一项自然规律与人类认知、工程技术与社会管理融合的系统工程。

山洪和泥石流风险管理的目的是认识灾害规律,协调人与山地环境的关系,在适度风险范围内规范山洪与泥石流灾害调控行为,增强人类适应与抗灾能力(resilience),合理利用灾害资源以满足社会经济可持续发展需要。

风险管理是当前山洪和泥石流灾害管理研究的重要课题。基于灾害风险理论及国内外相关研究成果(Blaikie et al, 2004; UN/ISDR, 2004; 史培军等, 2005; 崔鹏, 2014),山洪和泥石流风险管理主要包括灾害风险辨识、风险分析与评估、风险决策与处理等三方面的内容。

(1) 灾害风险辨识。是指对尚未发生灾害或具有潜在灾害形成客观条件的集水区,进行系统调查,分析影响山洪和泥石流灾害风险的各种因素,确定灾害形成的基本条件,辨识风险源,推断和预测灾害风险关键因素与可能风险的程度。

(2) 风险分析与评估。包括山洪和泥石流灾害危险性分析、承灾体易损性分析和风险评估三个方面。危险性分析主要针对孕灾环境和致灾因子,分析山洪和泥石流灾害形成因素,研究其动力过程,分析可能受灾害影响的范围、强度和频率(概率)。承灾体易损性分析主要研究承灾体遭受不同类型和强度山洪、泥石流袭击的动力响应,并综合分析承灾体抗灾能力,包括暴露要素分析、脆弱性分析、易损曲线分析等,确定易损度。风险评估是在危险性分析和易损性分析的基础上,估算不同强度和频率山洪、泥石流可能造成的风险损失,并进行风险等级划分,编制风险图。

(3) 风险决策与处置。根据山洪和泥石流灾害风险管理的目标,在灾害风险辨识和风险评估的基础上,确定可以接受的风险度,进行可能风险与可接受风险水平的权衡,做出应对风险的决策。综合考虑投资—收益与实施条件,从可以采取的减灾技术与途径(避难搬迁、监测预警、预防措施、治理工程、居民点规划与重建等)中选择适宜的措施与措施组合,综合制定风险处置方案,包括临预案与灾害防治工程方案,有效控制灾害风险。进而,通过政策、教育和社会管理等手段,提高受灾区域居民的减灾知识与防灾意识,形成有效的鼓励、激励和保障机制,构建高效的组织管理模式,使得政府减灾部门、减灾专家和社区居民形成有机配合的协调一致并可持续的减灾管理体系,保障风险决策和风险处置方案发挥最好的减灾作用,取得理想的风险调控成效。

在风险管理内容中,基于动力过程的灾害风险评估可以直接服务于居民点建设与防灾规划,居民参与式社区灾害风险管理是新的研究课题。

3.2 风险处置对策

本文重点讨论基于山洪、泥石流动力过程灾害风险评估结果的风险处置对策,主要包括:基于动力过程的灾害风险预测、科学合理的临灾预案和灾害防治工程方案。

(1) 基于动力过程的灾害风险预测

通过山洪、泥石流的动力学模拟方法,计算山洪和泥石流的致险范围内流速和流深的时空分布,确定其危险性。通过危险区内不同位置、不同类型承灾体的综合抗灾能力调查分析,确定易损度。在危险性分析和易损性分析的基础上,估算不同强度和频率山洪、泥石流可能危害的区域与潜在风险损失,并进行风险等级划分,编制风险图。根据调查的山洪、泥石流堆(冲)积扇形态与危害特征,以山洪和泥石流主流线为沟口堆积扇中轴线,可通过极坐标的形式来表达堆(冲)积扇上承灾体的空间分布(图6)(曾超, 2014)及其可能的风险,预测山洪和泥石流在不同等级风险区内灾害损失,为制定科学合理的土地利用规划与防灾方案提供科学依据。

图6   山洪和泥石流堆(冲)积扇上承灾体的分布与风险区

Fig.6   Distribution of exposure units in at-risk area of debris flows and flash floods

(2) 制定合理的临灾预案

山洪和泥石流灾害临灾预案要求:一要调查确定山洪和泥石流灾害的基本情况,包括灾害危险区、风险等级等;二要针对不同风险级别区,制定灾害的应急措施,明确不同危险等级灾害点监测内容、巡查要求、预警方式、临灾抢险救灾措施、疏散撤离路线等;三要明确组织指挥机构和抢险救灾队伍,公布责任人员名单和联系电话;四要明确灾害应急队伍人员分工与职责。临灾预案需通过宣传、教育、培训、综合演练等方式不断提高公众对临灾预案的知晓度,要采用居民喜闻乐见的方式(挂图、影像、卡片等)提高群众的防灾减灾意识,充分发动群众参与社区防灾减灾,认真贯彻落实临灾预案,加强群测群防和群专结合的减灾管理,落实联防机制。

(3) 建立灾害综合防治工程方案

山洪和泥石流灾害防治是对一个流域内水—土物质的综合调控。根据流域灾害发育规律和运动特征,提出了“坡面防护—支沟的稳坡固沟—主沟的拦挡调节—沟口的排导防护”的基于物质与能量过程调控的小流域山洪、泥石流灾害综合治理技术体系(图7)。具体工程措施主要包括:通过改善地表产汇流条件,减少降雨过多渗入土层,控制坡面产流和汇流过程,消减沟道洪峰,从而控制流水侵蚀动力作用;通过稳定坡脚和坡体,防治滑坡崩塌产生,减少固体物质供给;通过沟道拦挡和排导,减小泥石流的规模和破坏作用,减轻灾害风险。

图7   山洪和泥石流灾害综合治理技术体系

Fig.7   Engineering and biological control measures for debris flows and flash floods

3.3 社区风险管理

目前,社区参与式灾害风险管理是一种有效的灾害管理模式(图8),这种“自下而上”的灾害管理体制,重在提高社区居民的防灾认知、防灾减灾能力建设和灾害风险管理水平。《兵库行动框架协议》指出了社区减灾的重要作用,强调要利用知识、创新和教育建立一个安全的、在各个层面对灾害具有适应力的文化(UN/ISDR, 2005)。社区灾害风险管理的挑战之一是保持社区层面防灾减灾的持久化(Mano, 2011)。因此,如何建立长效的社会参与激励机制,营造居民自觉参与社区灾害风险管理全过程的文化氛围,探索政府主导的“自上而下”型和发挥风险社区居民自觉性的“自下而上”型相结合的社区参与式灾害风险管理模式,是值得进一步探索的重要课题。

图8   参与式社区风险管理模式(崔鹏, 2015)

Fig.8   Participatory community-based disaster risk management model (Cui, 2015)

4 结论与讨论

4.1 结论

(1) 认识山洪和泥石流形成过程与机理对防灾减灾具有重要意义。本文分析了山洪和泥石流的演化过程,阐述了四个关键形成演化过程:即地表产流流量激增过程、土体破坏物质供给激增过程、沟道堵塞体级联溃决流量放大过程和动床侵蚀规模增大过程。

(2) 山洪和泥石流风险评估分为区域评估与灾点评估。区域风险评估对具体灾点的详细信息反映不足,评估结果难以满足一个具体灾害点的临灾预案和工程设计需求。针对灾点减灾需求,依据具体灾害的物理力学性质,考虑灾害的动力过程,分析灾害体动力过程和灾害体与承灾结构体动力响应,介绍了灾害危险性分析与承灾体易损性分析方法,在此基础上计算灾点的风险值并进行风险等级划分,建立了基于动力过程的山洪、泥石流风险评估方法。

(3) 阐述了山洪、泥石流灾害风险管理的定义与主要内容,重点讨论了包括灾害风险预测、临灾预案、灾害防治工程方案等内容的基于山洪、泥石流动力过程的灾害风险处置对策,提出基于社区的参与式灾害风险管理模式,形成一套基于山洪、泥石流动力过程的风险评估与风险管理理论与方法。

4.2 讨论

(1) 本文在分析山洪和泥石流动力过程的基础上,结合风险管理目标,建立了基于动力过程的山洪、泥石流风险评估方法,可以直接服务于居民点建设与防灾规划,是一种值得推广的方法。然而,山地灾害风险分析和风险管理近期才得到重视,需要进一步深化和完善基于动力过程的灾害风险分析和风险管理理论和方法,更好地服务于防灾减灾。

(2) 科学合理的灾害风险管理可以较好地协调人与山地环境的关系,增强人类适应与抗灾能力,合理规划风险区并利用灾害资源满足社会经济可持续发展需要。中国目前还存在山区经济与教育水平不高、社会参与减灾积极性不够、社区减灾综合协调机制有待健全等问题。随着高强度人类活动而出现的灾害危险区人口密度与经济密度的增加,加强山洪、泥石流灾害风险评估和风险管理理论方法的研究,提高灾害风险管理水平,建立健全灾害风险管理体制和机制,有效降低特大灾害所造成的风险,是中国山洪和泥石流减灾的重要课题。因此,可操作的基于灾害动力过程评估与风险管理方法将会在社区减灾中发挥重要作用。

The authors have declared that no competing interests exist.


参考文献

[1] 崔鹏. 2014.

中国山地灾害研究进展与未来应关注的科学问题

[J]. 地理科学进展, 33(2): 145-152.

https://doi.org/10.11820/dlkxjz.2014.02.001      Magsci      [本文引用: 1]      摘要

本文首先简要回顾了山地灾害研究与防治方面的新进展:认识了山地灾害的空间分布规律,建立了山洪、泥石流、滑坡危险性评价方法;发展了滑坡稳定性分析的原理和计算方法,建立了泥石流流体应力本构关系、泥石流流速流量和冲击力计算公式、粘性泥石流起动模型,提出了山洪和泥石流规模放大效应;基于降雨和地面成灾环境要素耦合分析,发展了山地灾害气象预报方法;基于对灾害物理特性的认识,研发了一系列灾害监测预警仪器、数字流域平台与智能手机网络相结合的山洪预警系统;发展了灾害治理工程技术,形成了适合欠发达地区特点的灾害治理技术体系。在此基础上,分析了在灾害形成、运动、预测预报、防治技术和风险管理等方面还需要进一步深化研究的问题,提出山地灾害学科今后面临的任务。最后,针对国家减灾需求和学科发展目标,提出灾害对生态的响应机制、气候变化对山地灾害的影响与巨灾预测、水—土耦合的细观结构力学、灾害风险的理论与方法、基于灾害形成理论的机理预报模式、灾害防治技术规程的健全等未来应该关注的科学技术问题。

[Cui P.2014.

Progress and prospects in research on mountain hazards in China

[J]. Progress in Geography, 33(2): 145-152.]

https://doi.org/10.11820/dlkxjz.2014.02.001      Magsci      [本文引用: 1]      摘要

本文首先简要回顾了山地灾害研究与防治方面的新进展:认识了山地灾害的空间分布规律,建立了山洪、泥石流、滑坡危险性评价方法;发展了滑坡稳定性分析的原理和计算方法,建立了泥石流流体应力本构关系、泥石流流速流量和冲击力计算公式、粘性泥石流起动模型,提出了山洪和泥石流规模放大效应;基于降雨和地面成灾环境要素耦合分析,发展了山地灾害气象预报方法;基于对灾害物理特性的认识,研发了一系列灾害监测预警仪器、数字流域平台与智能手机网络相结合的山洪预警系统;发展了灾害治理工程技术,形成了适合欠发达地区特点的灾害治理技术体系。在此基础上,分析了在灾害形成、运动、预测预报、防治技术和风险管理等方面还需要进一步深化研究的问题,提出山地灾害学科今后面临的任务。最后,针对国家减灾需求和学科发展目标,提出灾害对生态的响应机制、气候变化对山地灾害的影响与巨灾预测、水—土耦合的细观结构力学、灾害风险的理论与方法、基于灾害形成理论的机理预报模式、灾害防治技术规程的健全等未来应该关注的科学技术问题。
[2] 崔鹏. 2015. 长江上游山地灾害与水土流失地图集[M]. 北京: 科学出版社.

[本文引用: 2]     

[Cui P.2015. Atlas of mountain hazards and soil erosion in the upper Yangtze[M]. Beijing: Science Press.]

[本文引用: 2]     

[3] 葛全胜, 邹铭, 郑景云. 2008. 中国自然灾害风险综合评估初步研究[M]. 北京: 科学出版社.

[本文引用: 1]     

[Ge Q S, Zou M, Zheng J Y.2008. Integrated assessment of natural disaster risk in China[M]. Beijing: Science Press.]

[本文引用: 1]     

[4] 国家防汛抗旱总指挥部, 中华人民共和国水利部. 2013. 2013中国水旱灾害公报[M]. 北京: 中国水利水电出版社.

[The State Flood Control and Drought Relief Headquarters, the Ministry of Water Resources of the People’s Republic of China. 2013. Bulletin of flood and drought disasters in China, 2013[M]. Beijing: China: Water & Power Press.]

[5] 胡凯衡, 韦方强. 2005.

基于数值模拟的泥石流危险性分区方法

[J]. 自然灾害学报, 14(1): 10-14.

https://doi.org/10.3969/j.issn.1004-4574.2005.01.002      URL      [本文引用: 1]      摘要

由于泥石流运动方程研究和计算 机技术的迅速发展使得基于数值模拟的泥石流危险性分区方法已成为泥石流危险性分区的主要方法。该方法也是这一领域最为有效的定量方法,并被大量应用于泥石 流灾害评估、土地利用规划和灾害制图中。在分析泥石流运动数值模拟方法和根据模拟结果进行分区研究现状的基础上,指出了该方法在数值模型的建立、分区指标 的确定以及指标分级中存在的问题,并针对这些问题提出了解决方法。

[Hu K H, Wei F Q.2005.

Numerical-simulation-based debris flow risk zoning

[J]. Journal of Natural Disasters, 14(1): 10-14.]

https://doi.org/10.3969/j.issn.1004-4574.2005.01.002      URL      [本文引用: 1]      摘要

由于泥石流运动方程研究和计算 机技术的迅速发展使得基于数值模拟的泥石流危险性分区方法已成为泥石流危险性分区的主要方法。该方法也是这一领域最为有效的定量方法,并被大量应用于泥石 流灾害评估、土地利用规划和灾害制图中。在分析泥石流运动数值模拟方法和根据模拟结果进行分区研究现状的基础上,指出了该方法在数值模型的建立、分区指标 的确定以及指标分级中存在的问题,并针对这些问题提出了解决方法。
[6] 史培军. 2005.

四论灾害系统研究的理论与实践

[J]. 自然灾害学报, 14(6): 1-7.

https://doi.org/10.3969/j.issn.1004-4574.2005.06.001      URL      [本文引用: 2]      摘要

作者曾分别于 1991,1995和2001年发表了关于灾害研究理论与实践的3篇文章,对灾害系统的性质、动力学机制等进行了探讨。本文在前3篇文章的基础上,从综合 灾害风险管理的角度,完善了灾害系统的结构与功能体系,论证了灾情形成过程中恢复力的作用机制,分析了区域开发与安全建设的互馈关系,构建了区域综合减灾 的行政管理体系,提出了由政府、企业与社区构成的区域综合减灾范式。研究结果表明,区域灾情形成过程中,脆弱性与恢复力有着明显的区别,脆弱性是区域灾害 系统中致灾因子、承灾体和孕灾环境综合作用过程的状态量,它主要取决于区域的经济发达程度与社区安全建设水平;恢复力则是灾害发生后,区域恢复、重建及安 全建设与区域发展相互作用的动态量,它主要取决于区域综合灾害风险行政管理能力、政府与企业投入和社会援助水平。区域安全水平与土地利用的时空格局和产业 结构关系密切,通过划定区域高风险“红线区”的办法,调整土地利用时空格局和产业结构,有利于建立区域可持续发展的综合减灾范式。针对区域自然灾害系统存 在着相互作用、互为因果的灾害链规律,以及灾害系统所具有的结构与功能特征,完善由纵向、横向和政策协调共同组成的一个“三维矩阵式”的区域综合减灾行政 管理体系,构建以政府为主导、企业为主体、社区全面参与的区域综合减灾范式。以此促进在发展中提高区域减灾能力,并在一定安全水平下,建设区域可持续发展 模式。

[Shi P J.2005.

Theory and practice on disaster system research in a fourth time

[J]. Journal of Natural Disasters, 14(6): 1-7.]

https://doi.org/10.3969/j.issn.1004-4574.2005.06.001      URL      [本文引用: 2]      摘要

作者曾分别于 1991,1995和2001年发表了关于灾害研究理论与实践的3篇文章,对灾害系统的性质、动力学机制等进行了探讨。本文在前3篇文章的基础上,从综合 灾害风险管理的角度,完善了灾害系统的结构与功能体系,论证了灾情形成过程中恢复力的作用机制,分析了区域开发与安全建设的互馈关系,构建了区域综合减灾 的行政管理体系,提出了由政府、企业与社区构成的区域综合减灾范式。研究结果表明,区域灾情形成过程中,脆弱性与恢复力有着明显的区别,脆弱性是区域灾害 系统中致灾因子、承灾体和孕灾环境综合作用过程的状态量,它主要取决于区域的经济发达程度与社区安全建设水平;恢复力则是灾害发生后,区域恢复、重建及安 全建设与区域发展相互作用的动态量,它主要取决于区域综合灾害风险行政管理能力、政府与企业投入和社会援助水平。区域安全水平与土地利用的时空格局和产业 结构关系密切,通过划定区域高风险“红线区”的办法,调整土地利用时空格局和产业结构,有利于建立区域可持续发展的综合减灾范式。针对区域自然灾害系统存 在着相互作用、互为因果的灾害链规律,以及灾害系统所具有的结构与功能特征,完善由纵向、横向和政策协调共同组成的一个“三维矩阵式”的区域综合减灾行政 管理体系,构建以政府为主导、企业为主体、社区全面参与的区域综合减灾范式。以此促进在发展中提高区域减灾能力,并在一定安全水平下,建设区域可持续发展 模式。
[7] 史培军, 邹铭, 李保俊, . 2005.

从区域安全建设到风险管理体系的形成: 从第一届世界风险大会看灾害与风险研究的现状与发展趋向

[J]. 地球科学进展, 20(2): 173-179.

Magsci      摘要

<p>综述了第一届世界风险大会关于灾害与风险科学研究的进展,分析了其发展趋势,并对我国加强这一领域的研究提出了对策。灾害与风险研究的目标正在趋向一致,都集中在为可持续发展模式提供科学依据;基于灾害研究的区域安全建设正在向基于风险研究的区域风险管理体系转变,由此认为可持续发展模式就是接受一定风险水平条件下的区域发展模式,一定风险水平依赖于由区域经济水平与社会文化体系相关联的区域安全建设水平;加强我国灾害与风险科学研究的学科建设、国家实验室建设、重大项目的组织与实施已势在必行;重点开展灾害脆弱性评价体系与评价模型、风险评价体系与风险评价模型、典型区域灾害与风险管理范式、灾害与风险应急控制预案(情景分析)、灾害与风险区划等方面的研究。为此把灾害科学和风险科学,与全球变化科学和可持续发展科学统一纳入到地球系统科学体系之中,以加强与完善我国在这些领域的综合多学科交叉研究能力。</p>

[Shi P J, Zou M, Li B J, et al.2005.

Regional safety construction and risk management system: the actuality and trend of the study of disaster and risk based on the world congress on risk

[J]. Advances in Earth Science, 20(2): 173-179.]

Magsci      摘要

<p>综述了第一届世界风险大会关于灾害与风险科学研究的进展,分析了其发展趋势,并对我国加强这一领域的研究提出了对策。灾害与风险研究的目标正在趋向一致,都集中在为可持续发展模式提供科学依据;基于灾害研究的区域安全建设正在向基于风险研究的区域风险管理体系转变,由此认为可持续发展模式就是接受一定风险水平条件下的区域发展模式,一定风险水平依赖于由区域经济水平与社会文化体系相关联的区域安全建设水平;加强我国灾害与风险科学研究的学科建设、国家实验室建设、重大项目的组织与实施已势在必行;重点开展灾害脆弱性评价体系与评价模型、风险评价体系与风险评价模型、典型区域灾害与风险管理范式、灾害与风险应急控制预案(情景分析)、灾害与风险区划等方面的研究。为此把灾害科学和风险科学,与全球变化科学和可持续发展科学统一纳入到地球系统科学体系之中,以加强与完善我国在这些领域的综合多学科交叉研究能力。</p>
[8] 涂勇, 何秉顺, 褚明华, . 2014.

2013年全国山洪灾害特征分析

[J]. 中国水利, (18): 18-22.

https://doi.org/10.3969/j.issn.1000-1123.2014.18.005      URL      摘要

利用国家防办2012年、2013年洪涝灾情统计数据和各省上报 的灾情信息,对2013年有人死亡的山洪灾害事件进行了整理,通过统计分析和现场调研的方法,从区域分布、历史对比、灾害类型、灾害等级、发生时间、人员 构成等方面对2013年全国山洪灾害特征进行了分析,绘制了全国2013年山洪灾害等级分布图和山洪灾害类型分布图,分析了2000年以来山洪灾害死亡人 数及其占洪涝灾害总死亡人数比例情况以及与2012年山洪灾害灾情的对比.结果表明:2013年共发生了181起造成人员死亡的山洪灾害事件,共有560 人死亡、221人失踪.2013年全国因洪涝灾害死亡774人,山洪灾害死亡人数占洪涝灾害死亡人数比例为72.4%,山洪灾害仍然是洪涝灾害死亡的主 因,溪沟洪水仍是山洪灾害中主要的灾害类型.全国山洪灾害事件持续时间长,主要集中在7、8月.2013年山洪灾害灾情区域性差异显著,局部省份山洪灾害 较重.多地外来务工、务农、旅游人员频遭山洪袭击,老年人和未成年人占总死亡人数比例较大.极端强降雨、流动人口山洪灾害防御意识淡薄、预警信息失效是造 成2013年山洪灾害频发并造成较大伤亡的主要原因.

[Tu Y, He B S, Chu M H, et al.2014.

Characteristics of mountain flood disasters in 2013

[J]. China Water Resources, (18): 18-22.]

https://doi.org/10.3969/j.issn.1000-1123.2014.18.005      URL      摘要

利用国家防办2012年、2013年洪涝灾情统计数据和各省上报 的灾情信息,对2013年有人死亡的山洪灾害事件进行了整理,通过统计分析和现场调研的方法,从区域分布、历史对比、灾害类型、灾害等级、发生时间、人员 构成等方面对2013年全国山洪灾害特征进行了分析,绘制了全国2013年山洪灾害等级分布图和山洪灾害类型分布图,分析了2000年以来山洪灾害死亡人 数及其占洪涝灾害总死亡人数比例情况以及与2012年山洪灾害灾情的对比.结果表明:2013年共发生了181起造成人员死亡的山洪灾害事件,共有560 人死亡、221人失踪.2013年全国因洪涝灾害死亡774人,山洪灾害死亡人数占洪涝灾害死亡人数比例为72.4%,山洪灾害仍然是洪涝灾害死亡的主 因,溪沟洪水仍是山洪灾害中主要的灾害类型.全国山洪灾害事件持续时间长,主要集中在7、8月.2013年山洪灾害灾情区域性差异显著,局部省份山洪灾害 较重.多地外来务工、务农、旅游人员频遭山洪袭击,老年人和未成年人占总死亡人数比例较大.极端强降雨、流动人口山洪灾害防御意识淡薄、预警信息失效是造 成2013年山洪灾害频发并造成较大伤亡的主要原因.
[9] 王鑫, 曹志先, 岳志远. 2009.

强不规则地形上浅水二维流动的数值计算研究

[J]. 水动力学研究与进展, 24(1): 56-62.

URL      [本文引用: 1]      摘要

浅水流动数值模拟技术取得了显著进步,已能够有效地捕捉激波,但是应用于强不规则地形时计算 困难,其原因在于:陡峭坡面上水流重力作用对数值稳定性有重要影响。采用简单的算子分裂法,即使是时间步长满足CFL稳定性准则,也可能致使计算不稳定, 或者要求远比CFL稳定性准则更加苛刻的时间步长,从而严重影响计算效率。本文建立对强不规则地形和一般地形条件均能保持高效率的浅水二维数学模型:基本 控制方程采用二维浅水方程;在空间上采用能捕捉激波的WAFTVD二阶精度格式;在时间上采用三阶TVD龙格库塔法,并采用自适应步长满足数值稳定性要 求。将该模型应用于冰湖溃决洪水的模拟,试验结果表明:该文模型从本质上解除了强不规则地形对时间步长的约束,与现有数学模型相比,计算效率达到数量级上 的提高,为强不规则地形上浅水二维流动数值模型的广泛应用(如高原与山丘区溃坝洪水数值预报等)提供了新的技术基础。

[Wang X, Cao Z X, Yue Z Y.2009.

Numerical modeling of shallow flows over irregular topography

[J]. Journal of Hydrodynamics (Ser. A), 24(1): 56-62.]

URL      [本文引用: 1]      摘要

浅水流动数值模拟技术取得了显著进步,已能够有效地捕捉激波,但是应用于强不规则地形时计算 困难,其原因在于:陡峭坡面上水流重力作用对数值稳定性有重要影响。采用简单的算子分裂法,即使是时间步长满足CFL稳定性准则,也可能致使计算不稳定, 或者要求远比CFL稳定性准则更加苛刻的时间步长,从而严重影响计算效率。本文建立对强不规则地形和一般地形条件均能保持高效率的浅水二维数学模型:基本 控制方程采用二维浅水方程;在空间上采用能捕捉激波的WAFTVD二阶精度格式;在时间上采用三阶TVD龙格库塔法,并采用自适应步长满足数值稳定性要 求。将该模型应用于冰湖溃决洪水的模拟,试验结果表明:该文模型从本质上解除了强不规则地形对时间步长的约束,与现有数学模型相比,计算效率达到数量级上 的提高,为强不规则地形上浅水二维流动数值模型的广泛应用(如高原与山丘区溃坝洪水数值预报等)提供了新的技术基础。
[10] 曾超. 2014.

泥石流作用下建筑物易损性评价方法[D]

. 成都: 中国科学院大学.

[本文引用: 1]     

[Zeng C.2014.

Vulnerability assessment of building to debris flow hazard[D].

Chengdu: University of Chinese Academy of Sciences.]

[本文引用: 1]     

[11] 张志彤. 2015.

关于2014年防汛抗旱工作的报告

[J]. 中国防汛抗旱, 25(1): 8-13.

URL      摘要

正(2015年1月20日)根据会议安排,我就2014年全国防汛抗旱工作和2015年安排意见报告如下。1.1汛情2014年全国平均降雨量622 mm,与常年基本持平。华东、西南等地部分地区降雨量较常年偏多1~2成。华西秋1 2014年汛情、旱情和灾情

[Zhang Z T.2015.

Guanyu 2014 nian fangxun kanghan gongzuo de baogao

, 25(1): 8-13.]

URL      摘要

正(2015年1月20日)根据会议安排,我就2014年全国防汛抗旱工作和2015年安排意见报告如下。1.1汛情2014年全国平均降雨量622 mm,与常年基本持平。华东、西南等地部分地区降雨量较常年偏多1~2成。华西秋1 2014年汛情、旱情和灾情
[12] 邹强. 2014.

川藏公路泥石流灾害风险分析[D]

. 成都: 中国科学院大学.

[本文引用: 1]     

[Zou Q.2014.

Risk analysis of debris flows along Sichuan-Tibet highways[D].

Chengdu, China: University of Chinese Academy of Sciences.]

[本文引用: 1]     

[13] AGS.2007.

Commentary on guideline for landslide susceptibility, hazard and risk zoning for land use management[R].

Australian: Australian Geomechanics Society Landslide Task force Landslide Zoning Working Group, 42(1): 37-62.

[本文引用: 1]     

[14] Bisson M, Favalli M, Fornaciai A, et al.2005.

A rapid method to assess fire-related debris flow hazard in the Mediterranean region: an example from Sicily (southern Italy)

[J]. International Journal of Applied Earth Observation and Geoinformation, 7(3): 217-231.

https://doi.org/10.1016/j.jag.2005.04.003      URL      [本文引用: 1]      摘要

Increased atmospheric temperatures during the high-pressure which characterise the Mediterranean climate in the dry summer time, coupled with an increase in the intensity of storms in the following wet season over recent years, increase the risk of summer fires as well as debris flows and floods in the autumn and/or in the following years. In addition, the diffuse urbanization of Italy requires a rapid and reliable tool be available in order to obtain preliminary information, at the end of the summer season, that identifies newly fired areas that present a significant hazard to human populations. In such burned zones, soil instability may be more severe favouring debris flows which may impact on populated zones. Thus, in this paper we discuss a rapid methodology to: (i) identify burned areas using band ratio's using multitemporal LANDSAT ETM images; (ii) evaluate the potential of the burned areas as the source of debris flows based on morphometric parameters (slope and hill slope curvature); (iii) evaluate the structures, such as houses and roads, exposed to potential damage by debris flows. Hazardous areas were evaluated using a stochastical model coupled with an empirical relationship which accounts for the mobility of the debris flows. The methodology provides a classification of the most “dangerous” burned areas and the potentially maximum inundated downslope areas. This has been applied to Sicily for the period autumn 2001–autumn 2002. The total burned area was 76.37km 2 . According to the classification proposed 6.4% of the burned areas were consider of very high to high hazard potential, 54.4% of medium hazard and 43.2% of low hazard potential.
[15] Blaikie P, Cannon T, Davis I, et al.2004. At risk: natural hazards, people's vulnerability and disasters[M]. New York: Roublisher Press.

[本文引用: 1]     

[16] Chau K T, Lo K H.2004.

Hazard assessment of debris flows for Leung King Estate of Hong Kong by incorporating GIS with numerical simulations

[J]. Natural Hazards and Earth System Science, 4(1): 103-116.

https://doi.org/10.5194/nhess-4-103-2004      URL      [本文引用: 1]      摘要

As over seventy percent of the land of Hong Kong is mountainous, rainfall-induced debris flows are not uncommon in Hong Kong. The objective of this study is to incorporate numerical simulations of debris flows with GIS to identify potential debris flow hazard areas. To illustrate this approach, the proposed methodology is applied to Leung King Estate in Tuen Mun. A Digital Elevation Model (DEM) of the terrain and the potential debris-flow sources were generated by using GIS to provide the required terrain and flow source data for the numerical simulations. A theoretical model by Takahashi et al. (1992) improved by incorporating a new erosion initiation criterion was used for simulating the runout distances of debris flows. The well-documented 1990 Tsing Shan debris flow, which occurred not too far from Leung King Estate, was used to calibrate most of the flow parameters needed for computer simulations. Based on the simulation results, a potential hazard zone was identified and presented by using GIS. Our proposed hazard map was thus determined by flow dynamics and a deposition mechanism through computer simulations without using any so- called expert opinions, which are bounded to be subjective and biased.
[17] Cui P, Hu K H, Zhang J Q, et al.2011.

Prediction of the debris-flow danger area by combining hydrological and inundation simulation methods

[J]. Journal of Mountain Science, 8(1): 1-9.

https://doi.org/10.1007/s11629-011-2040-8      URL      [本文引用: 1]     

[18] Cui P, Zhou G G D, Zhu X H, et al.2013.

Scale amplification of natural debris flows caused by cascading landslide dam failures

[J]. Geomorphology, 182: 173-189.

https://doi.org/10.1016/j.geomorph.2012.11.009      URL      Magsci      [本文引用: 2]      摘要

Debris flows are typically caused by natural terrain landslides triggered by intense rainfalls. If an incoming mountain torrent collapses a series of landslide dams, large debris flows can form in a very short period. Moreover, the torrent can amplify the scale of the debris flow in the flow direction. The catastrophic debris flows that occurred in Zhouqu, China, on 8 August 2010 were caused by intense rainfall and the upstream cascading failure of landslide dams along the gullies. In the wake of the incident, a field study was conducted to better understand the process of cascading landslide dam failures and the formation of debris flows. This paper looks at the geomorphic properties of the debris-flow gullies, estimates the peak flow discharges at different locations using three different methods, and analyzes the key modes (i.e., different landslide dam types and their combinations) of cascading landslide dam failures and their effect on the scale amplification of debris flows. The results show that five key modes in Luojiayu gully and two modes in Sanyanyu gully accounted for the scale amplification of downstream debris flows in the Zhouqu event. This study illustrates how the hazardous process of natural debris flows can begin several kilometers upstream as a complex cascade of geomorphic events (failure of landslide dams and erosion of the sloping bed) can scale to become catastrophic discharges. Neglecting recognition of these hazardous geomorphic and hydrodynamic processes may result in a high cost. (C) 2012 Elsevier B.V. All rights reserved.
[19] Cui P, Zou Q, Xiang L Z, et al.2013.

Risk assessment of simultaneous debris flows in mountain townships

[J]. Progress in Physical Geography, 37(4): 516-542.

https://doi.org/10.1177/0309133313491445      Magsci      [本文引用: 1]      摘要

Many mountain towns in China are located on the joint alluvial fans of multiple and adjacent past debris flows, making them vulnerable to large, multiple, and simultaneous debris flows during heavy rainfall. Without emergency management planning, such flows, often appearing with interconnecting and chain-reaction processes, can lead to extensive loss of life and property. In the Wenchuan earthquake-affected area, such disasters are common. We analyzed the compound effects of simultaneous debris flow events, and proposed three quantitative methods of debris risk assessment based on kinetic energy, flow depth, and inundation depth. Validated using a field study of actual debris flow disasters, these analyses are useful in determining the type, quantity, distribution, economic worth, and susceptibility of hazard-affected objects in a region. Subsequently, we established a method to determine the vulnerability of different hazard-affected objects, particularly concerning the susceptibility indexes of buildings or structures. By analyzing the elements underlying hazard formation conditions, damage potential, and the socio-economic conditions of mountain townships, we proposed a systematic and quantitative method for risk analysis of mountain townships. Finally, the proposed method was applied to a case study of Qingping Township, which was affected by 21 simultaneous debris flows triggered by a 50-year return period precipitation event. The proposed method analyzed the superposition and chain-reaction effects of disasters and divided the affected area of the township into three risk zones. The analysis indicated that the calculated risk zones coincide with the actual distribution and severity of damage in the debris flow event, which suggests that the risk assessment is consistent with results from the actual disaster.
[20] Fuchs S, Ornetsmüller C, Totschnig R.2012.

Spatial scan statistics in vulnerability assessment: an application to mountain hazards

[J]. Natural Hazards, 64(3): 2129-2151.

https://doi.org/10.1007/s11069-011-0081-5      [本文引用: 1]      摘要

Abstract<br/>In the European Alps, the concept of risk has increasingly been applied in order to reduce the susceptibility of society to mountain hazards. Risk is defined as a function of the magnitude and frequency of a hazard process times consequences; the latter being quantified by the value of elements at risk exposed and their vulnerability. Vulnerability is defined by the degree of loss to a given element at risk resulting from the impact of a natural hazard. Recent empirical studies suggested a dependency of the degree of loss on the hazard impact, and respective vulnerability (or damage-loss) functions were developed. However, until now, only little information is available on the spatial characteristics of vulnerability on a local scale; considerable ranges in the loss ratio for medium process intensities only provide a hint that there might be mutual reasons for lower or higher loss rates. In this paper, we therefore focus on the spatial dimension of vulnerability by searching for spatial clusters in the damage ratio of elements at risk exposed. By using the software SaTScan, we applied an ordinal data model and a normal data model in order to detect spatial distribution patterns of five individual torrent events in Austria. For both models, we detected some significant clusters of high damage ratios, and consequently high vulnerability. Moreover, secondary clusters of high and low values were found. Based on our results, the assumption that lower process intensities result in lower damage ratios, and therefore in lower vulnerability, and vice versa, has to be partly rejected. The spatial distribution of vulnerability is not only dependent on the process intensities but also on the overall land use pattern and the individual constructive characteristics of the buildings exposed. Generally, we suggest the use of a normal data model for test sites exceeding a minimum of 30 elements at risk exposed. As such, the study enhanced our understanding of spatial vulnerability patterns on a local scale.<br/>
[21] Huggel C, Kääb A, Haeberli W, et al.2003.

Regional-scale GIS-models for assessment of hazards from glacier lake outbursts: evaluation and application in the Swiss Alps

[J]. Natural Hazards and Earth System Science, 3(6): 647-662.

https://doi.org/10.5194/nhess-3-647-2003      URL      [本文引用: 1]      摘要

Debris flows triggered by glacier lake outbursts have repeatedly caused disasters in various high-mountain regions of the world. Accelerated change of glacial and periglacial environments due to atmospheric warming and increased anthropogenic development in most of these areas raise the need for an adequate hazard assessment and corresponding modelling. The purpose of this paper is to pro-vide a modelling approach which takes into account the current evolution of the glacial environment and satisfies a robust first-order assessment of hazards from glacier-lake outbursts. Two topography-based GIS-models simulating debris flows related to outbursts from glacier lakes are presented and applied for two lake outburst events in the southern Swiss Alps. The models are based on information about glacier lakes derived from remote sensing data, and on digital elevation models (DEM). Hydrological flow routing is used to simulate the debris flow resulting from the lake outburst. Thereby, a multiple- and a single-flow-direction approach are applied. Debris-flow propagation is given in probability-related values indicating the hazard potential of a certain location. The debris flow runout distance is calculated on the basis of empirical data on average slope trajectory. The results show that the multiple-flow-direction approach generally yields a more detailed propagation. The single-flow-direction approach, however, is more robust against DEM artifacts and, hence, more suited for process automation. The model is tested with three differently generated DEMs (including aero-photogrammetry- and satellite image-derived). Potential application of the respective DEMs is discussed with a special focus on satellite-derived DEMs for use in remote high-mountain areas.
[22] Liu X L, Lei J Z.2003.

A method for assessing regional debris flow risk: an application in Zhaotong of Yunnan province (SW China)

[J]. Geomorphology, 52(3-4): 181-191.

https://doi.org/10.1016/S0169-555X(02)00242-8      URL      [本文引用: 1]      摘要

中国科学院机构知识库(中国科学院机构知识库网格(CAS IR GRID))以发展机构知识能力和知识管理能力为目标,快速实现对本机构知识资产的收集、长期保存、合理传播利用,积极建设对知识内容进行捕获、转化、传播、利用和审计的能力,逐步建设包括知识内容分析、关系分析和能力审计在内的知识服务能力,开展综合知识管理。
[23] Mano T.2011.

Community-based disaster management and public awareness[R]. Disaster Risk Vulnerablity Conference. New Delhi, India:

Mahatma Gandhi University.

[本文引用: 1]     

[24] O'Brien J S, Julien P Y, Fullerton W T.1993.

Two-dimensional water flood and mudflow simulation

[J]. Journal of Hydraulic Engineering, 119(2): 244-261.

[25] Pradhan B.2010.

Remote sensing and GIS-based landslide hazard analysis and cross-validation using multivariate logistic regression model on three test areas in Malaysia

[J]. Advances in Space Research, 45(10): 1244-1256.

https://doi.org/10.1016/j.asr.2010.01.006      [本文引用: 1]      摘要

<h2 class="secHeading" id="section_abstract">Abstract</h2><p id="">This paper presents the results of the cross-validation of a multivariate logistic regression model using remote sensing data and GIS for landslide hazard analysis on the Penang, Cameron, and Selangor areas in Malaysia. Landslide locations in the study areas were identified by interpreting aerial photographs and satellite images, supported by field surveys. SPOT 5 and Landsat TM satellite imagery were used to map landcover and vegetation index, respectively. Maps of topography, soil type, lineaments and land cover were constructed from the spatial datasets. Ten factors which influence landslide occurrence, i.e., slope, aspect, curvature, distance from drainage, lithology, distance from lineaments, soil type, landcover, rainfall precipitation, and normalized difference vegetation index (ndvi), were extracted from the spatial database and the logistic regression coefficient of each factor was computed. Then the landslide hazard was analysed using the multivariate logistic regression coefficients derived not only from the data for the respective area but also using the logistic regression coefficients calculated from each of the other two areas (nine hazard maps in all) as a cross-validation of the model. For verification of the model, the results of the analyses were then compared with the field-verified landslide locations. Among the three cases of the application of logistic regression coefficient in the same study area, the case of Selangor based on the Selangor logistic regression coefficients showed the highest accuracy (94%), where as Penang based on the Penang coefficients showed the lowest accuracy (86%). Similarly, among the six cases from the cross application of logistic regression coefficient in other two areas, the case of Selangor based on logistic coefficient of Cameron showed highest (90%) prediction accuracy where as the case of Penang based on the Selangor logistic regression coefficients showed the lowest accuracy (79%). Qualitatively, the cross application model yields reasonable results which can be used for preliminary landslide hazard mapping.</p>
[26] Smith R E, Smettem K R J, Broadbridge P, et al.2002.

Infiltration theory for hydrologic applications

[M]. Washington DC: American Geophysical Union.

[本文引用: 1]     

[27] UN/ISDR (United National International Strategy for Disaster Reduction). 2004.

Living with risk: a global review of disaster reduction initiatives

[M]. Geneva: United Nations Publications.

[本文引用: 1]     

[28] UN/ISDR (The United Nations Office for Disaster Risk Reduction). 2005.

Hyogo framework for action 2005-2015: building the resilience of nations and communities to disasters

[R]. Kobe, Hyogo, Japan: UN World Conference on Disaster Risk Reduction.

[本文引用: 2]     

[29] UN/ISDR (The United Nations Office for Disaster Risk Reduction). 2015.

Sendai framework for action 2015-2030: building the resilience of nations and communities to disasters

[R]. Sendai, Japan: Third UN World Conference on Disaster Risk Reduction.

[本文引用: 1]     

[30] UNDHA (United Nations, Department of Humanitarian Affairs). 1992.

Internationally Agreed Glossary of Basic Terms Related to Disaster Management

[R]. Geneva: United Nations Department of Humanitarian Affairs.

[本文引用: 1]     

[31] UNDP (United Nations Development Program). 2004.

A global report reducing disaster risk: a challenge for development[R].

New York: UNDP, l-14.

[本文引用: 1]     

/