一种余震影响范围的模拟方法技术

本发明专利技术公开了一种余震影响范围的模拟方法，包括：建立以余震点经度为X轴、余震点纬度为Y轴与余震点震级为Z轴的坐标系；以X轴和Y轴所在面为投影面，将坐标系变换到投影坐标系中；利用点格局分析的关联分析方法，拟合得到余震点间距离的关系模式，将其分解为不同的无标度区间，无标度区间内的数值小于相邻的下一个无标度区间内的数值；在无标度区间内进行点的核密度模拟分析，分别得到无标度区间对应的图层；判断相邻图层Li与Li+1之间的遮挡内容，去除Li+1的相应遮挡内容后，Li与Li+1进行组合，依次处理所有图层，输出组合后的图层Ri。本发明专利技术基于公开发布的地震余震的点数据，可快速评估地震影响范围。

【技术实现步骤摘要】

本专利技术涉及地震测量
，特别涉及一种余震影响范围的模拟方法。
技术介绍
地震影响范围模拟是地震灾害综合评估的基础工作之一。地震灾害突发性强，可预测性差，近年来频度较高，伴生的严重次生灾害对社会影响很大。在此背景下，研究基于公开、可获取的地震余震点序列数据，进行地震影响范围快速模拟，其现实价值是显而易见的。现有对地震灾害影响范围评估的方法，多以最大地震震级、地面峰值加速度等参数为主，实施地震灾害影响范围分区评估，这种方法需要较多的地震物理参数，因此过程繁琐，可操作性不强。基于公开发布的地震余震的点数据，快速估计地震影响范围，数据来源简单，成本更低；能够基于实时发布的数据量，进行多次评估和逼近，实时性及操作简易性更强。
技术实现思路
本专利技术旨在至少解决现有技术中存在的问题之一。本专利技术要达到的目的之一在于基于公开发布的地震余震的点数据，可快速实现估计地震影响范围。为解决上述技术问题，本专利技术提供一种余震影响范围的模拟方法，其步骤包括：建立以余震点经度为X轴、余震点纬度为Y轴与余震点震级为Z轴的坐标系；以所述X轴和Y轴所在面为投影面，将所述坐标系变换到投影坐标系中；利用点格局分析的关联分析方法，拟合得到余震点间距离的关系模式，将其分解为不同的无标度区间S1、S2、…、Sn，其中n为余震点个数且n大于等于2，其中无标度区间Si内的数值小于无标度区间Si+1内的数值，其中1≤i≤n-1；在所述无标度区间内进行点的核密度模拟分析，分别得到所述无标度区间S1、S2、…、Sn对应的图层L1、L2、…、Ln；判断相邻所述图层Li与Li+1之间的遮挡内容，去除Li+1的相应遮挡内容后，Li与Li+1进行组合，依次处理所有图层L1、L2、…、Ln，输出组合后的图层Ri。进一步的，其步骤还包括根据所述图层Ri识别余震影响范围。进一步的，用于制作地震影响范围的地图。本专利技术的有益效果在于提供了一种余震影响范围的模拟方法，本专利技术基于公开发布的地震余震的点数据，可快速评估地震影响范围，数据来源简单，成本更低；同时，本专利技术能够基于实时发布的数据量，进行多次评估和逼近，实时性及操作简易性更强。附图说明图1所示为本专利技术实施例一种余震影响范围的模拟方法的流程图。图2所示为本专利技术实施例单一尺度为320km的余震的核密度估计图。图3所示为本专利技术实施例多尺度融合的余震的核密度估计图。具体实施方式下文将结合具体实施例详细描述本专利技术。应当注意的是，下述实施例中描述的技术特征或者技术特征的组合不应当被认为是孤立的，它们可以被相互组合从而达到更好的技术效果。如图1所述，本专利技术提供了一种余震影响范围的模拟方法，其步骤包括：100：建立以余震点经度为X轴、余震点纬度为Y轴与余震点震级为Z轴的坐标系；200：以所述X轴和Y轴所在面为投影面，将所述坐标系变换到投影坐标系中；300：利用点格局分析的关联分析方法，拟合得到余震点间距离的关系模式，将其分解为不同的无标度区间S1、S2、…、Sn，其中n为余震点个数且n大于等于2，其中无标度区间Si内的数值小于无标度区间Si+1内的数值，其中1≤i≤n-1；400：在所述无标度区间内进行点的核密度模拟分析，分别得到所述无标度区间S1、S2、…、Sn对应的图层L1、L2、…、Ln；500：判断相邻所述图层Li与Li+1之间的遮挡内容，去除Li+1的相应遮挡内容后，Li与Li+1进行组合，依次处理所有图层L1、L2、…、Ln，输出组合后的图层Ri。600：根据所述图层Ri识别余震影响范围。本专利技术的一种余震影响范围的模拟方法还可以应用于制作地震影响范围的地图。本专利技术实施例以汶川地区为实施对象：1.在国家地震科学数据分享中心(国家地震科学数据分享中心)，下载汶川地震后的余震数据，录入到EXCEL中并进行整理，数据属性包括经度(X)、纬度(Y)和震级(Z)；在ArcGIS的ArcMap中，选择Toos下的AddXYdata，导入整理好的EXCEL数据，建立以余震点经度为X轴、余震点纬度为Y轴与余震点震级为Z轴的坐标系并进行投影转化，最后投影为公里网；加载汶川余震涉及市、县的矢量边界；2.在EXCEL中对余震点进行关联维分析，其分析公式为：C(r)＝(2Nr)/N(N-1)，(R<r)，其中R为某两对地震之间的距离，r为选取的用来衡量地震间距的标尺，Nr为距离R小于r的地震对数，N为选用的地震总对数；对于选定的一系列r值，做出双对数曲线lg[C(r)]-lg(r)，在曲线上找出近似直线的那一段，用最小二乘法进行拟合D2＝lg[C(r)]/lg(r)；r采用等间距递增方法，增量为500m，取值范围为1km—324km，得到最终拟合结果表明：余震点在13.5—20km、30—43km、66—82km、225—236km、317—321.5km五个区间内惯关联特征显著；3.在空间关联分析的基础上，选取无标度区间的中值16km，36km，74km，230km，320km作为核密度估计带宽；4.在ArcGIS的ArcMap中打开KernelDensity工具，inputpoint选择汶川余震点数据，populationfield选择汶川余震点数据的震级字段，searchradius首先设置为320km，输出图层320；5.在不影响实验结果，保证效果图美观的条件下，过滤掉320的最外层。SpatialAnalyst，RasterCalculator，320>图层最小值，得到calculation；6.在ArcGIS的ArcMap中打开ExtractbyMask工具，inputraster选择calculation，inputrasterorfeaturemaskdata选择320，输出c-320；7.SpatialAnalyst，RasterCalculator，[c-320]*[320]得到calculation2，int整理，得到calculation3，如图2所示；8.按照上述步骤分别计算核密度估计带宽为230km，74km，36km，16km时的余震影响范围；将最小带宽16km的结果分为两级，其余带宽结果为一级，如图3所示。本专利技术提供了一种余震影响范围的模拟方法，本专利技术基于公开发布的地震余震的点数据，可快速评估地震影响范围，数据来源简单，成本更低；同时，本专利技术能够基于实时发布的数据量，进行多次评估和逼近，实时性及操作简易性更强。本文虽然已经给出了本专利技术的一些实施例，但是本领域的技术人员应当理解，在不脱离本本文档来自技高网...

【技术保护点】
一种余震影响范围的模拟方法，其特征在于，其步骤包括：建立以余震点经度为X轴、余震点纬度为Y轴与余震点震级为Z轴的坐标系；以所述X轴和Y轴所在面为投影面，将所述坐标系变换到投影坐标系中；利用点格局分析的关联分析方法，拟合得到余震点间距离的关系模式，将其分解为不同的无标度区间S1、S2、…、Sn，其中n为余震点个数且n大于等于2，其中无标度区间Si内的数值小于无标度区间Si+1内的数值，其中1≤i≤n‑1；在所述无标度区间内进行点的核密度模拟分析，分别得到所述无标度区间S1、S2、…、Sn对应的图层L1、L2、…、Ln；判断相邻所述图层Li与Li+1之间的遮挡内容，去除Li+1的相应遮挡内容后，Li与Li+1进行组合，依次处理所有图层L1、L2、…、Ln，输出组合后的图层Ri。

【技术特征摘要】
1.一种余震影响范围的模拟方法，其特征在于，其步骤包括：
建立以余震点经度为X轴、余震点纬度为Y轴与余震点震级为Z轴的坐标系；
以所述X轴和Y轴所在面为投影面，将所述坐标系变换到投影坐标系中；
利用点格局分析的关联分析方法，拟合得到余震点间距离的关系模式，将其分解为不同
的无标度区间S1、S2、…、Sn，其中n为余震点个数且n大于等于2，其中无标度区间Si内
的数值小于无标度区间Si+1内的数值，其中1≤i≤n-1；
在所述无标度区间内进行点的核密度模拟分析...

【专利技术属性】
技术研发人员：罗明良，董丞妍，白雷超，谭壮，昌小莉，邹宝裕，张斌，
申请(专利权)人：西华师范大学，
类型：发明
国别省市：四川;51