一种地震成灾风险评估和经济损失预测方法技术

本发明专利技术提供一种地震成灾风险评估和经济损失预测方法，包括：收集在记录历史上地震灾害等级大于5级的成灾和不成灾数据集P{i}，其中包含地震发生地点的经纬度坐标(X

【技术实现步骤摘要】
一种地震成灾风险评估和经济损失预测方法
本专利技术属于地震成灾风险评估和地震经济损失预测
，具体涉及一种地震成灾风险评估和经济损失预测方法。
技术介绍
地震灾害是一种极为常见的自然灾害，它具有突发性强，续发性高和破坏力巨大等特点，其也被称之为“群灾之首”。尤其在人口密集的地方，地震灾害一旦发生，将极有可能给人类的生产，生活造成难以估量的损失。地震的成灾性是指地震发生是否会造成直接经济损失，或者人员伤亡事故。因此一旦地震灾害发生，如何对其成灾风险进行有效评估，以及对成灾地震灾害进行合理的经济损失预测，这对于地震灾害的定量化分析和之后的地震减灾工作的开展具有十分重要的意义。目前常见的地震灾害损失评估有两类方法，一种是基于建筑物，基础设施结构分类的清单易损性分析法，这种方法相对比较简便，但是该方法往往需要收集统计不同地区的大范围数据资料，将及其消耗人力和物力，这对于地震应急、救灾等的宏观决策并不经济。另一种方法是基于GDP的宏观易损性分析方法，该方法包括了建筑物和基础设施在内的全部社会财富遭受地震灾害的总损失。这种方法简便易行，而在资料数据上，国家每年都会更新各个地区的GDP数值，对于数据更新上更为容易。现如今已经广泛应用于地震灾害经济损失预测的实际工作中。例如，专利[CN110046454A]中在宏观易损性分析方法上，结合大数定理和全概率公式，来确定某一位置发生地震损失概率进而进一步估计经济损失。然而这种方法往往忽略了地震的成灾性分析。而通过上述方法进行分析评估，不仅会消耗一定的人力成本，并且也可能造成评估数据不准确的现象。如果在地震发生时，能够先对地震成灾风险进行评估，确定其成灾的可能性，然后再进行地震灾害的损失预测，这将极大提升工作效率以及预测的准确性。针对上述所存在的问题，需要设计出一种地震成灾风险评估和经济损失预测方法，通过该方法来满足当前对地震成灾风险评估和地震灾害损失预测的需求。
技术实现思路
有鉴于此，本专利技术提出一种地震成灾风险评估和经济损失预测方法。根据历史上的中国大陆地震灾害损失述评，并对地震数据进行收集整理，构建了一个快速、准确的地震成灾风险评估和经济损失预测方法，工作人员只需要根据新发生地震的地点的经纬度，震级大小，作为地震成灾风险评估模型的输入，如果判断其为成灾事件，则获取该地区当年的GPD值，通过宏观经济易损模型进行损失性预测其可能造成的直接经济损失。通过采用该方法能够极大的提升预测的准确性，工作人员的工作效率，并且降低了人工和硬件成本。为实现上述目的，本专利技术提供如下技术方案：本专利技术一种地震成灾风险评估和经济损失预测方法，其包括1、历史地震数据收集模块，主要为了提取历史数据中的5级以上成灾和不成灾的地震记录，其中包括了地震地点的经纬度，地震的震级，人员伤亡数，受灾人口数，地震烈度以及直接经济损失；2、地震成灾风险评估模型构造模块，该模块是为了通过历史数据判断当前地震灾害是否成灾；3、地震灾害数据采集模块，其目的是为了收集当前地震灾害的实时数据，其中包括发生地震的经纬度，震级，已经发生地当年的GDP，并对其进行分析；4、烈度回归模型构造模块，通过历史的地震级数对地震烈度进行线性回归。5、宏观经济易损性模型构造模块，该模块是为了根据历史的烈度数据和经济损失比拟合经济易损性模型；6、直接经济损失预测模块，目的是为了预测地震成灾的直接经济损失。一种地震成灾风险评估和经济损失预测方法，大致包括以下步骤：步骤(1)：收集在记录历史上地震灾害等级大于5级的成灾和不成灾数据集P{i}，其中包含地震发生地点的经纬度坐标(Xi，Yi)，地震震级大小Mi，地震烈度等级Di，地震的直接经济损失Li以及地震受灾人口数量Napi。步骤(2)：获取目标地震位置T的经纬度坐标(XT，YT)，地震等级MT，该地区当年的GDPT。步骤(3)：通过如下公式(1),(2)筛选出目标地震位置与历史数据P{i}中距离小于200千米的成灾次数N1和不成灾次数N2。其中R表示地球的半径，一般取6378.137千米，D(x)表示地震目标位置与历史地震位置之间的距离，如果距离小于等于200千米，则D(x)＝1,否则D(x)＝0。N1表示在所有历史地震数据集合P{i}中，地震成灾，并且与目标地震T的距离小于200千米的次数。N2表示在所有历史地震数据集合P{i}中，地震不成灾，并且与目标地震T的距离小于200千米的次数。步骤(4)：判断N1和N2的比例，如果N1/(N1+N2)>0.5，则说明目标地震是可能成灾则进行下一步，否则说明该地震成灾可能性不大，则返回步骤(2)继续获取相应地震数据。步骤(5)：如果判断目标地震T可能成灾，则继续筛选集合P{i}中距离目标地震T距离小于200千米成灾地震数据中的震级集合Ma和地震烈度集合Da，其中a<＝i。由于地震的烈度与震级往往都是正相关，所以这里可以采用一元线性回归的方法对其进行线性拟合，可获得D＝k×M+b。其中拟合系数b和k可通过公式(3)和公式(4)计算得到。步骤(6)：根据上述步骤(2)中所获取的目标地震的震级MT，通过步骤(5)中拟合方程可以计算目标地震位置的烈度MD，其中MD＝k×MT+b。步骤(7)：通过历史地震数据集P{i}，构造关于烈度与GDP的损失率公式如下：F(M,GDP)＝C×A×MB(5)其中，F(I,GDP)表示GDP损失率，M为地震烈度，A、B为易损性统计参数，C为修正系数。步骤(8)：根据步骤(2)中所获取的目标地震地区当年的GDP值MG，带入公式(5)中，即可预测出该地区的直接经济损失值LOSS：LOSS＝C×A×MDB×MG(6)与现有技术相比，本专利技术的有益效果为：通过现有分类清单易损性分析法进行地震灾害直接经济损失预测，该方法主要对于地震危险性的分析工作比较多，但是对于地震灾害的危害性分析比较少。此外，该方法要求需要及为详尽的建筑物、地质、社会经济方面的大量资料，往往这些资料的获取是十分困难的。而普通的宏观经济易损模型通常会忽略地震灾害的风险评估分析，通常情况下只要有地震灾害的发现该方法就会认定其是风险灾害，就会对其进行经济损失的预测，这往往缺乏一定的合理性。为此，本专利技术方法结合历史地震数据进行分析建模，通过对目标地震位置以及历史相邻位置的地震信息采用概率统计的方法进行风险评估，如果判断其为成灾地震，再通过宏观经济模型进行直接经济损失的预测。本专利技术方法能够更为准确的对目标地震的风险和损失进行有效的评估，对于地震灾害的抗御对策研究和保障社会公共安全都有很高的应用需求和前景。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。本文档来自技高网...

【技术保护点】
1.一种地震成灾风险评估和经济损失预测方法，其特征在于，包括以下步骤：/n步骤(1)：收集在记录历史上地震灾害等级大于5级的成灾和不成灾数据集P{i}，其中包含地震发生地点的经纬度坐标(X

【技术特征摘要】
1.一种地震成灾风险评估和经济损失预测方法，其特征在于，包括以下步骤：
步骤(1)：收集在记录历史上地震灾害等级大于5级的成灾和不成灾数据集P{i}，其中包含地震发生地点的经纬度坐标(Xi，Yi)，地震震级集合Mi，地震烈度等级集合Di，地震的直接经济损失Li以及地震受灾人口数量Napi；
步骤(2)：获取目标地震位置T的经纬度坐标(XT，YT)，地震震级集合MT，该地区当年的GDPT；
步骤(3)：通过如下公式(1)、公式(2)筛选出目标地震位置与历史数据P{i}中距离小于200千米的成灾次数N1和不成灾次数N2；



如果P{i}成灾，D(x)＝1(1)



如果P{i}不成灾，D(x)＝1(2)
其中R表示地球的半径，D(x)表示地震目标位置与历史地震位置之间的距离，如果距离小于等于200千米，则D(x)＝1,否则D(x)＝0；
N1表示在所有历史地震数据集合P{i}中，地震成灾，并且与目标地震位置T的距离小于200千米的次数；
N2表示在所有历史地震数据集合P{i}中，地震不成灾，并且与目标地震位置T的距离小于200千米的次数；
步骤(4)：判断N1和N2的比例，如果N1/(N1+...

【专利技术属性】
技术研发人员：俞一奇，邱彦林，陈尚武，
申请(专利权)人：杭州叙简科技股份有限公司，
类型：发明
国别省市：浙江;33