一种地震导致滑坡所造成的人员死亡的快速评估方法技术

一种地震导致滑坡所造成的人员死亡的快速评估方法，其首先将研究区域按地质类别划分为不同的地质区域，然后获得研究区域的坡度数据；并计算各地质区域的安全系数Fs，又计算出各地质区域的极限加速度及地震动峰值加速度数据，之后计算可能造成各地质区域人口死亡的滑坡危险性参数并以县为单位，求得每个县滑坡危险性参数平均值；然后计算每个县的人口脆弱性参数，对于每个县所辖区域中地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域，计算折减系数，最后计算出地震引起的次生滑坡可能造成的人员死亡数量。本发明专利技术所提供的一种地震导致滑坡所造成的人员死亡的快速评估方法，其可快速评估由地震引发的次生滑坡灾害可能造成的人员死亡情况。

【技术实现步骤摘要】
一种地震导致滑坡所造成的人员死亡的快速评估方法
本专利技术涉及一种评估地震造成的损失的方法，特别是一种对地震导致的次生滑坡灾害可能造成的人员死亡的快速评估方法。
技术介绍
滑坡是斜坡上的岩土体，在重力和地表水，地下水等作用下，沿一定的滑坡面整体向下滑动的一种现象。强烈地震时，地震诱发的滑坡次生地质灾害，特别是在山地、丘陵地区，其经济损失和人员伤亡甚至比地震直接造成的还要大。在中国大陆地区，尤其是地形相对复杂的山区，地震造成的滑坡是最为常见且破坏力最强的次生地震灾害。然而，地震造成的滑坡具有发生突然、机理复杂、运动形态多样的特征，预测比较困难。各国政府为减轻地震造成的滑坡灾害，都采取了积极有效的措施。其中，在地震应急体系中，特别提到“防御次生灾害”，通过加强次生灾害监测预警，组织专家开展险情排查、评估和除险加固等方式，提高政府应对震后滑坡灾害的应急指挥和救援能力，减少人民的生命财产损失。在针对地震造成的滑坡应急指挥系统中，地震造成的滑坡灾害快速评估是应急决策系统的重要基础和科学依据，如果能够快速预估出地震造成的滑坡伤亡情况，就可以为应急资源需求提供参考数据；也可以根据人员伤亡情况，制定物资医疗需求和救灾人员需求方案，在最短的时间内使灾民得到营救和安置，从而使人员伤亡水平降到最低。我国在“九五”期间建立了地震灾害快速评估模型，此外，在“灾害学”杂志2010年10月第25卷增刊中的文章“地震人员伤亡快速评估模型研究”(高惠瑛、李清霞)中，作者提出了对地震中人口损失的评估方法。然而，关于地震造成的滑坡灾害人员伤亡快速评估方法研究却寥寥无几，在“岩石力学与工程学报”杂志2010年10月第10期第29卷的文章“地震滑坡灾害快速评估技术及对应急影响研究”(王秀英)中，作者也仅提出了一些研究思路，却并没有涉及具体的解决方法。可见，在现有评估方法中，尚无针对地震造成的滑坡次生灾害可能造成人员死亡的有效评估方法，特别是适于震后应急期间的快速评估方法。
技术实现思路
本专利技术要解决的技术问题是提供一种地震导致滑坡所造成的人员死亡的快速评估方法，以减少或避免前面所提到的问题。为解决上述技术问题，本专利技术提供了一种地震导致滑坡所造成的人员死亡的快速评估方法，其用于在地震灾害发生后，快速评估研究区域内由地震引起的次生滑坡可能造成的人员死亡数量C，其包括如下步骤，步骤A，在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域，并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数；步骤B，在计算机上根据DEM的30米分辨率数据，利用Arcgis软件获得所述研究区域的坡度数据；步骤C，在计算机上根据公式获得各所述地质区域的安全系数Fs，上式中β为坡度；步骤D，在计算机上根据公式ac＝(FS-1)gsinβ获得各所述地质区域的极限加速度ac；步骤E，将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据，利用公式求得平均值，然后进行插值计算，获得各所述地质区域的地震动峰值加速度amax数据，并利用Arcgis软件获得各所述地质区域amax-ac的数据；步骤F，在计算机上根据下列公式计算可能造成各地质区域人口死亡概率的滑坡危险性参数P1数据，当amax-ac≥0时，求得P1＝0.05-0.12(amax-ac-0.446)2，当上式P1＜0时，则P1取值为0；当amax-ac＜0时，则P1＝0，根据行政区划数据，以县为单位，利用Arcgis软件求得每个县滑坡危险性参数平均值P1i，i表示所述研究区域的第i县；步骤G，按照下列公式计算每个县的人口脆弱性参数，P2i＝(0.28-0.1Ia)5.4-0.55M，其中，上式中M为地震震级，Ca为该县15-59岁人口数量，Ci为该县人口总人数，i表示所述研究区域的第i县；步骤H，对于每个县所辖区域中地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域，按照下列公式计算折减系数S折：上式中，S总为所述县总面积，SVI为所述县地震烈度大于等于VI或地震动峰值加速度大于等于0.05g的区域的面积；对于所述研究区域的第i县，其折减系数即为S折i；步骤I，根据下列公式计算在所述研究区域由地震引起的次生滑坡可能造成的人员死亡数量C，本专利技术所提供的一种地震导致滑坡所造成的人员死亡的快速评估方法，其可在地震发生后，快速评估由地震引发的次生滑坡灾害可能发生的位置以及其可能造成的人员死亡情况，对选择启动应急预案的应急响应级别具有较高参考价值，为决策者指定灾民的营救和安置计划提供依据。附图说明以下附图仅旨在于对本专利技术做示意性说明和解释，并不限定本专利技术的范围。其中，图1显示的是根据本技术的一个具体实施例的一种地震导致滑坡所造成的人员死亡的快速评估方法的研究区域的地质区域划分图；图2为图1基础上的研究区域的坡度图；图3为图2基础上的安全系数数据图；图4为图3基础上的极限加速度图；图5为图1基础上的地震动峰值加速度图；图6为图4、5基础上的amax-ac差值图；图7为图6基础上的地震导致滑坡危险性参数图；图8为图7基础上的县级地震导致滑坡危险性参数平均值数据图。具体实施方式为了对本专利技术的技术特征、目的和效果有更加清楚的理解，现说明本专利技术的具体实施方式。本专利技术提供了一种地震导致滑坡所造成的人员死亡的快速评估方法，其用于在地震灾害发生后，快速评估研究区域内由地震引起的次生滑坡可能造成的人员死亡数量C，其包括如下步骤，步骤A，在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域，并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数；不同的地质类别对应的参数如下表所示；类别内聚力(MPa)土容重(kN/m3)内摩擦角(°)火山岩0.2227.540变质岩0.1726.533碎屑岩0.125.524碳酸盐岩0.0723.516第四系松散堆积物0.0521.513步骤B，在计算机上根据DEM的30米数据，利用Arcgis软件获得所述研究区域的坡度数据；基于30米的DEM数字高程模型(DigitalElevationModel)可通过“地理空间数据云”的ASTERGDEM获取，并利用Arcgis软件获得所述研究区域的坡度数据，步骤C，在计算机上根据公式获得各所述地质区域的安全系数Fs，上式中β为坡度；根据步骤A、B所获得的数据，可在Arcgis软件中根据上述公式计算出各所述地质区域的分辨率为30米的各基础地面数据单元的安全系数Fs，在计算之前，可先将内摩擦角和坡度β先换算为同一单位，例如，Arcgis软件在计算三角函数时需要以弧度为单位输入数据，则在计算前可先将内摩擦角和坡度β均换算为弧度进行输入。步骤D，在计算机上根据公式ac＝(FS-1)gsinβ获得各所述地质区域的极限加速度ac；在计算得到步骤C的安全系数Fs数据后，可在Arcgis软件中根据上述公式进一步计算出各所述地质区域的分辨率为30米的各基础地面数据单元的极限加速度ac。步骤E，将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据，利用公式求得平均值，然后进行插值计算，获得各所述地质区域的地震动峰值加速度amax数据，并利用Arcgis软件获得各所述地质区域amax-ac的数据；也就是说，可在A本文档来自技高网...

【技术保护点】
一种地震导致滑坡所造成的人员死亡的快速评估方法，其特征在于，其用于在地震灾害发生后，快速评估研究区域内由地震引起的次生滑坡可能造成的人员死亡数量C，其包括如下步骤，步骤A，在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域，并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数；步骤B，在计算机上根据DEM的30米数据，利用Arcgis软件获得所述研究区域的坡度数据；步骤C，在计算机上根据公式获得各所述地质区域的安全系数Fs，上式中β为坡度；步骤D，在计算机上根据公式ac＝(Fs‑1)gsinβ获得各所述地质区域的极限加速度ac；步骤E，将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据，利用公式求得平均值，然后进行插值计算，获得各所述地质区域的地震动峰值加速度amax数据，并利用Arcgis软件获得各所述地质区域amax‑ac的数据；步骤F，在计算机上根据下列公式计算可能造成各地质区域人口死亡概率的滑坡危险性参数P1数据，当amax‑ac 0时，求得P1＝0.05‑0.12(amax‑ac‑0.446)2，当上式P1＜0时，则P1取值为0；当amax‑ac＜0时，则P1＝0，根据行政区划数据，以县为单位，利用Arcgis软件求得每个县滑坡危险性参数平均值P1i，i表示所述研究区域的第i县；步骤G，按照下列公式计算每个县的人口脆弱性参数，P2i＝(0.28‑0.1Ia)5.4‑0.55M，其中，上式中M为地震震级，Ca为该县15‑59岁人口数量，Ci为该县人口总人数，i表示所述研究区域的第i县；步骤H，对于每个县所辖区域中地震烈度大于等于VI(或地震动峰值加速度大于等于0.05g)的区域，按照下列公式 计算折减系数S折：上式中，S总为所述县总面积，SVI为所述县地震烈度大于等于VI(或地震动峰值加速度大于等于0.05g)的区域的面积；步骤I，根据下列公式计算在所述研究区域由地震引起的次生滑坡可能造成的人员死亡数量...

【技术特征摘要】
1.一种地震导致滑坡所造成的人员死亡的快速评估方法，其特征在于，其用于在地震灾害发生后，快速评估研究区域内由地震引起的次生滑坡可能造成的人员死亡数量C，其包括如下步骤，步骤A，在计算机上利用现有地质图将所述研究区域按地质类别划分为不同的地质区域，并获得各所述地质区域的内聚力c、土容重r及内摩擦角参数；步骤B，在计算机上根据DEM的30米数据，利用Arcgis软件获得所述研究区域的坡度数据；步骤C，在计算机上根据公式获得各所述地质区域的安全系数Fs，上式中β为坡度；步骤D，在计算机上根据公式ac＝(FS-1)gsinβ获得各所述地质区域的极限加速度ac；步骤E，将地震时监测台站采集到的南北向aEW和东西向aNS两个峰值加速度数据，利用公式求得平均值，然后进行插值计算，获得各所述地质区域的地震动峰值加速度amax数据，并利用Arcgis软件获得各所述地质区域amax-ac的数据；步骤F，在计算机上根据下列公式计算可能造成各地质区域人口死亡概率的滑坡危险...

【专利技术属性】
技术研发人员：刘吉夫，商璟璐，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：北京;11