CN107368938B - 一种单体滑坡风险损失定量评估方法 - Google Patents

一种单体滑坡风险损失定量评估方法 Download PDF

Info

Publication number
CN107368938B
CN107368938B CN201710412562.0A CN201710412562A CN107368938B CN 107368938 B CN107368938 B CN 107368938B CN 201710412562 A CN201710412562 A CN 201710412562A CN 107368938 B CN107368938 B CN 107368938B
Authority
CN
China
Prior art keywords
landslide
risk
disaster
loss
range
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710412562.0A
Other languages
English (en)
Other versions
CN107368938A (zh
Inventor
宋强辉
李晓雷
徐洪
覃小华
孙从露
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing Institute of Geology and Mineral Resources
Logistical Engineering University of PLA
Original Assignee
Chongqing Institute of Geology and Mineral Resources
Logistical Engineering University of PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Chongqing Institute of Geology and Mineral Resources, Logistical Engineering University of PLA filed Critical Chongqing Institute of Geology and Mineral Resources
Priority to CN201710412562.0A priority Critical patent/CN107368938B/zh
Publication of CN107368938A publication Critical patent/CN107368938A/zh
Application granted granted Critical
Publication of CN107368938B publication Critical patent/CN107368938B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Educational Administration (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种单体滑坡(包括土质滑坡、岩质滑坡)风险损失定量评估方法,属于地质灾害防灾减灾与风险管理研究领域。该方法基于上、下限定理和区间分析法,从下限和上限两个方向来逼近滑坡灾害风险损失真值,克服了过去采用单一定值风险损失评估结果所带来的不确定性大、评估结果与滑坡灾害实际损失偏差较大等弊端,具有快速、动态、实用等特点,可服务于滑坡灾害定量风险评估,为制定滑坡灾害风险管理措施提供决策依据。

Description

一种单体滑坡风险损失定量评估方法
技术领域
本发明属于地质灾害防灾减灾与风险管理研究领域,涉及一种可用于单体滑坡灾害风险损失预测评估的方法。
背景技术
滑坡灾害是致灾体作用于承灾体的结果,同时也是滑坡的自然属性与社会经济属性共同作用的结果,应分别从滑坡发生的可能性和滑坡可能造成的损失,即风险损失(RiskLosses-RL)或期望损失(潜在损失)两方面来考虑。因此,滑坡灾害风险损失评估(亦称危害性评估)是滑坡灾害风险评估的两个重要基础工作之一。通过开展单体滑坡灾害风险损失定量评估方法研究,从而为滑坡灾害进行有效的风险管理奠定基础,以减少滑坡灾害造成的生命财产损失和节约防灾减灾工程费用。
由于滑坡风险损失评估研究具有重要的意义和工程价值,日益成为国内外众多专家学者和工程师的研究热点和难点。综观前人的研究成果可以发现,国内外自20世纪80年起开始研究滑坡灾害风险评估以来,已取得了较丰硕的成果,已基本建立了滑坡灾害风险评估和管理的理论框架,并有一定的工程应用。然而,尽管单体滑坡风险评估研究已经逐步深入和完善,但是却限于问题的复杂性和不确定因素太多,尚没有形成系统完善的理论和方法体系,尤其是理论研究与工程实践结合程度不高、风险评估方法的实用性、可推广性不强,有进一步发展与突破的可能和必要。特别是在滑坡风险损失进行定量评估方面,基本上将风险损失定量表示为单值形式,看似对滑坡可能造成的风险损失进行“精确的量化”,但滑坡灾害真正发生时,往往发现滑坡风险损失预测值与滑坡灾害实际损失(人员伤亡和直接经济损失)相差甚大,有时甚至是数倍之巨。这样的“误差”已经不是误差,而是错误了,显然是我们所不能接受的。
发明内容
该发明综合应用了上、下限定理和区间分析法,从下限和上限两个方向来逼近滑坡风险损失真值,克服了过去采用定值评估结果不确定性大,评估结果与实际损失偏差较大等问题,具有快速、动态、实用等特点,特别是倡导最大限度地利用现有的科学技术方法和信息有效开展单体滑坡风险损失定量评估,可为制定风险管理措施提供决策依据。
为实现本发明目的而采用的技术方案是这样的,一种单体滑坡灾害风险损失定量评估的方法,其特征在于按照以下步骤进行:
1)基于滑坡空间预测和运动模型,分析预测滑坡风险区范围的下限和上限:
首先采用本地区历史滑坡统计数据,获取本地区滑坡体规模—滑距(或风险区范围)关系,经过相似比分析和归一化处理,获得滑坡风险区范围(即滑坡危害范围)的下限和上限,有条件时还可以结合DAN3D软件、LS-RAPID软件等国际著名滑坡运动模拟分析软件对滑坡风险区范围进行定量数值模拟,结合数值模拟结果对前述经统计分析获得的风险区范围的下限与上限进行修正,若在统计分析风险区范围区间内,则以软件数值模拟结果作为新的风险区范围的下限和上限,综合确定滑坡风险区范围的下限和上限;
2)通过高分辨率遥感技术、无人机航拍及人工现场调查等方法,识别风险区内各种承灾体,并进行分类和承灾体价值量化,确定全部承灾体价值的下限Ei 和上限
Figure BDA0001312257370000021
(包括风险区内公众总人数和直接经济损失或财产损失的原始总价值,其中承灾体总数为n,下标i表示第i个承灾体,可为人或财产);
3)估算承灾体受灾(即遭受滑坡体冲击破坏或掩埋)的空间概率下限Psi 、上限
Figure BDA0001312257370000022
和时间概率下限Pti 、上限
Figure BDA0001312257370000023
Psti Psi ·Pti 计算承灾体受灾的时空概率下限Psti ,由
Figure BDA0001312257370000024
计算时空概率上限
Figure BDA0001312257370000025
4)根据财产损失率和生命/健康损害率的估计,估算所有承灾体易损性的下限Vi 和上限
Figure BDA0001312257370000026
5)采用公式
Figure BDA0001312257370000027
计算滑坡风险损失下限RL,采用公式
Figure BDA0001312257370000028
计算滑坡风险损失上限
Figure BDA0001312257370000029
6)计算滑坡风险损失极差
Figure BDA00013122573700000210
与滑坡风险损失均值
Figure BDA00013122573700000211
比较极差是否超过均值的30%,若
Figure BDA00013122573700000212
则进入步骤7;若
Figure BDA00013122573700000213
则进一步收集获取资料信息,进一步精细化滑坡空间预测和运动模型,重复步骤1)至6),直至满足
Figure BDA0001312257370000031
为止;
7)最终确定单体滑坡风险损失范围
Figure BDA0001312257370000032
值得说明的是,既然风险定义为损失的不确定性,那么风险损失应采用非单点值表示,应采用区间表达更和合理和准确。因此,本发明提出基于上、下限定理和区间分析法,从下限和上限两个方向来逼近滑坡风险损失真值,以克服过去采用定值评估结果不确定性大,评估结果与实际损失偏差较大等问题,具有快速、动态、实用等特点,可服务于滑坡灾害定量风险评估,为制定风险管理措施提供决策依据。
本发明提出的方法仅考虑滑坡可能造成的人员伤亡和财产损失等直接经济损失,暂未考虑对社会、资源与环境等间接经济损失,但所提出的模型方法对滑坡灾害其他损失评估同样具有一定的适用性,可供借鉴和参考。
附图说明
图1为本发明单体滑坡灾害风险损失定量评估方法的流程图
图2为滑坡要素关系示意图
图3为本地区滑坡体规模——滑距(或风险区范围)关系图
图4为滑坡风险区范围下限L
图5为滑坡风险区范围上限
Figure BDA0001312257370000033
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
在进行步骤1之前,应首先限定分析范围(滑坡要素关系示意图见图2),尽最大可能搜集各类相关资料,包括地形、地貌、工程地质、水文地质、气候降雨、地震活动、水位变化、人类工程活动等资料信息和相关数据以及历史文献资料等。
1)基于滑坡空间预测和运动模型,分析预测滑坡风险区范围下限和上限:首先采用本地区历史滑坡统计数据,获取本地区滑坡体规模——滑距(或危害范围)关系(如图3),经过相似比分析和归一化处理,获得风险区范围(即滑坡危害范围)下限(图3中斜线a)和上限(图3中斜线b),有条件时还可结合DAN3D软件、LS-RAPID软件等国际著名滑坡运动模拟分析软件对滑坡风险区范围进行定量数值模拟,结合数值模拟结果对前述经统计分析风险区范围的下限和上限进行修正,若在统计分析风险区范围区间内,则以软件数值模拟结果作为新的风险区范围的下限和上限,综合确定滑坡风险区范围下限L(如图4)和上限
Figure BDA0001312257370000048
(如图5)。
2)通过高分辨率遥感技术、无人机航拍及人工现场调查等方法,识别风险区内各种承灾体,并进行分类和承灾体价值量化,确定全部承灾体价值的下限Ei 和上限
Figure BDA0001312257370000047
(包括公众人数和直接经济损失或财产损失的价值,其中承灾体总数为n,下标i表示第i个承灾体,可为人或财产),经过风险识别,获得风险区承灾体为一栋3层砖混结构建筑物,以及建筑物内人员共计15人:
2-1)其中,60岁以上老人3人,下限与上限均为3人;20~60岁5人,下限与上限均为5人;10~20岁3人,下限与上限均为3人;10岁以下4人,下限与上限均为4人;
2-2)其中,建筑物总价值300万元,下限与上限均为300万元。
3)估算承灾体受灾(即遭受滑坡体冲击破坏或掩埋)的空间概率下限Psi 、上限
Figure BDA0001312257370000041
和时间概率下限Pti 、上限
Figure BDA0001312257370000042
Psti Psi ·Pti 计算承灾体受灾的时空概率下限Psti ,由
Figure BDA0001312257370000043
计算时空概率上限
Figure BDA0001312257370000044
由此可获得承灾体的时空概率如下表所示:
Figure BDA0001312257370000045
4)根据财产损失率和生命/健康损害率的估计,估算所有承灾体易损性下限Vi 和上限
Figure BDA0001312257370000046
如表所示:
Figure BDA0001312257370000051
5)采用公式
Figure BDA0001312257370000052
计算滑坡风险损失下限RL,采用公式
Figure BDA0001312257370000053
计算滑坡风险损失上限
Figure BDA0001312257370000054
求得结果如下表所示:
Figure BDA0001312257370000055
*备注:对生命损失(人)取整的原因是考虑到人的生命不可以为小数。
6)计算滑坡风险损失极差
Figure BDA0001312257370000056
与滑坡风险损失均值
Figure BDA0001312257370000057
比较极差是否超过均值的30%,若
Figure BDA0001312257370000058
则进入步骤7;若
Figure BDA0001312257370000059
则进一步收集获取资料信息,进一步精细化滑坡空间预测和运动模型,重复步骤1至6,直至满足
Figure BDA00013122573700000510
为止。
经过计算,滑坡风险损失极差:
6-1)对生命损失:
Figure BDA00013122573700000511
为4人-3人=1人,滑坡风险损失均值
Figure BDA00013122573700000512
满足
Figure BDA0001312257370000061
进入下一步。
*备注:对生命损失(人)取整的原因是考虑到人的生命不可以为小数。
6-2)对财产损失(直接经济损失):
Figure BDA0001312257370000062
Figure BDA0001312257370000063
滑坡风险损失均值
Figure BDA0001312257370000064
满足
Figure BDA0001312257370000065
进入下一步。
7)最终确定单体滑坡风险损失范围
Figure BDA0001312257370000066
7-1)对生命损失:
Figure BDA0001312257370000067
7-2)对财产损失(直接经济损失):
Figure BDA0001312257370000068

Claims (1)

1.一种单体滑坡灾害风险损失定量评估的方法,其特征在于按照以下步骤进行:
1)搜集滑坡要素资料;所述滑坡要素资料包括地形、地貌、工程地质、水文地质、气候降雨、地震活动、水位变化、人类工程活动资料信息和相关数据以及历史文献资料;
2)基于滑坡空间预测和运动模型,分析预测滑坡风险区范围的下限和上限:
首先采用本地区历史滑坡统计数据,获取本地区滑坡体规模—滑距或风险区范围关系,经过相似比分析和归一化处理,获得滑坡风险区范围的下限和上限;结合DAN3D软件或LS-RAPID软件对滑坡风险区范围进行定量数值模拟;结合数值模拟结果对经统计分析获得的滑坡风险区范围的下限和上限进行修正,若在统计分析风险区范围区间内,则以软件数值模拟结果作为新的风险区范围的下限和上限,综合确定滑坡风险区范围的下限和上限;
3)通过高分辨率遥感技术、无人机航拍或人工现场调查方法,识别风险区内各种承灾体,并进行分类和承灾体价值量化,确定全部承灾体价值的下限Ei 和上限
Figure FDF00000126806400000112
其中,承灾体为人或财产;承灾体总数为n;承灾体价值包括风险区内公众总人数和直接经济损失或财产损失的原始总价值;下标i表示第i个承灾体;
4)估算承灾体受灾的空间概率下限Psi 、上限
Figure FDF0000012680640000011
和时间概率下限Pti 、上限
Figure FDF0000012680640000012
Psti Psi ·Pti 计算承灾体受灾的时空概率下限Psti ,由
Figure FDF0000012680640000013
计算时空概率上限
Figure FDF0000012680640000014
5)根据财产损失率和生命/健康损害率的估计,估算所有承灾体易损性的下限Vi 和上限
Figure FDF0000012680640000015
6)采用公式
Figure FDF0000012680640000016
计算滑坡风险损失下限RL,采用公式
Figure FDF0000012680640000017
计算滑坡风险损失上限
Figure FDF0000012680640000018
7)计算滑坡风险损失极差
Figure FDF0000012680640000019
与滑坡风险损失均值
Figure FDF00000126806400000110
比较极差是否超过均值的30%,若
Figure FDF00000126806400000111
则进入步骤8);若
Figure FDF0000012680640000021
则进一步收集获取资料信息,进一步精细化滑坡空间预测和运动模型,重复步骤1)至7),直至满足
Figure FDF0000012680640000022
为止;
8)最终确定单体滑坡风险损失范围
Figure FDF0000012680640000023
CN201710412562.0A 2017-06-03 2017-06-03 一种单体滑坡风险损失定量评估方法 Active CN107368938B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710412562.0A CN107368938B (zh) 2017-06-03 2017-06-03 一种单体滑坡风险损失定量评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710412562.0A CN107368938B (zh) 2017-06-03 2017-06-03 一种单体滑坡风险损失定量评估方法

Publications (2)

Publication Number Publication Date
CN107368938A CN107368938A (zh) 2017-11-21
CN107368938B true CN107368938B (zh) 2021-09-21

Family

ID=60305296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710412562.0A Active CN107368938B (zh) 2017-06-03 2017-06-03 一种单体滑坡风险损失定量评估方法

Country Status (1)

Country Link
CN (1) CN107368938B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110046454A (zh) * 2019-04-25 2019-07-23 中国地震局地质研究所 概率地震经济损失计算方法及系统
CN110298489A (zh) * 2019-05-31 2019-10-01 成都理工大学 一种隧道大变形灾害的风险量化及支护决策方法
CN112183874B (zh) * 2020-10-12 2021-06-18 兰州大学 一种滑坡滑距预测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002288386A (ja) * 2001-03-28 2002-10-04 National Maritime Research Institute 船舶災害のリスク評価方法
CN103605897A (zh) * 2013-11-29 2014-02-26 民政部国家减灾中心 一种区域滑坡灾害风险估计方法
CN105654414A (zh) * 2015-12-25 2016-06-08 浙江大学城市学院 基于开源体系架构和建筑物空间数据库的城市多灾害风险损失评测系统及方法
CN105765582A (zh) * 2014-10-30 2016-07-13 伊迈格创新技术学院 用于评估由地下采矿产生的高能大地突裂风险的方法及系统
CN112686522A (zh) * 2020-12-25 2021-04-20 中南大学 一种Newmark修正模型地震滑坡危险性评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002288386A (ja) * 2001-03-28 2002-10-04 National Maritime Research Institute 船舶災害のリスク評価方法
CN103605897A (zh) * 2013-11-29 2014-02-26 民政部国家减灾中心 一种区域滑坡灾害风险估计方法
CN105765582A (zh) * 2014-10-30 2016-07-13 伊迈格创新技术学院 用于评估由地下采矿产生的高能大地突裂风险的方法及系统
CN105654414A (zh) * 2015-12-25 2016-06-08 浙江大学城市学院 基于开源体系架构和建筑物空间数据库的城市多灾害风险损失评测系统及方法
CN112686522A (zh) * 2020-12-25 2021-04-20 中南大学 一种Newmark修正模型地震滑坡危险性评估方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
单体滑坡灾害财产风险定量评估模型;吴越 等;《岩土力学》;20101130;第31卷(第S2期);第342-348页 *

Also Published As

Publication number Publication date
CN107368938A (zh) 2017-11-21

Similar Documents

Publication Publication Date Title
Pullan et al. Spatial parasite ecology and epidemiology: a review of methods and applications
Lee et al. Assessment of shallow landslide susceptibility using the transient infiltration flow model and GIS-based probabilistic approach
CN107368938B (zh) 一种单体滑坡风险损失定量评估方法
Wang et al. Probabilistic seismic landslide hazard maps including epistemic uncertainty
Chen et al. Resilience assessment of regional areas against earthquakes using multi-source information fusion
CN104616093A (zh) 一种地震灾害救援指挥调度系统及方法
CN114201922A (zh) 基于InSAR技术的动态滑坡敏感性预测方法及系统
Ma et al. Spatial prediction strategy for landslides triggered by large earthquakes oriented to emergency response, mid-term resettlement and later reconstruction
Kolathayar et al. Comprehensive probabilistic seismic hazard analysis of the Andaman–Nicobar regions
Oliver Geostatistics, rare disease and the environment
Haining The special nature of spatial data
CN113487105A (zh) 一种滑坡地质灾害的空间预测方法、装置和计算机设备
Zhang et al. Developing GIS-based earthquake loss model: a case study of Baqiao District, China
Vishnu et al. Challenges of modeling rainfall triggered landslides in a data-sparse region: A case study from the Western Ghats, India
Li et al. Seismic hazard models for typical urban masonry structures considering optimized regression algorithms
Wei et al. Risk assessment of people trapped in earthquake disasters based on a single building: A case study in Xichang city, Sichuan Province, China
Pailoplee et al. CU-PSHA: A Matlab software for probabilistic seismic hazard analysis
CN112200399A (zh) 一种地震成灾风险评估和经济损失预测方法
Wang et al. Identifying anomalous signals in GPS data using HMMs: An increased likelihood of earthquakes?
CN112259248A (zh) 一种covid-19境外输入风险评估预测方法
Subramanian et al. Constructing and validating geographically refined HAZUS-MH4 hurricane wind risk models: A machine learning approach
CN114925576A (zh) 一种基于灾害动力过程的地震滑坡危险性评估方法与系统
Hwang et al. Automated extraction of community mobility measures from GPS stream data using temporal DBSCAN
Khodashahi et al. Earthquake-induced landslides hazard zonation of Rudbar-Manjil using CAMEL model
Liu et al. Arias intensity attenuation relationship in Sichuan–Yunnan region, China

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant