基于深度震相自动匹配的远震起震点深度精定位技术的制作方法

1.本发明属于远震起震点定位技术领域，具体涉及一种基于深度震相自动匹配的远震起震点深度精定位技术。


背景技术：

2.世界测震权威机构通常采用震相走时方法确定地震起震点位置，进行地震速报和生成地震事件目录。例如美国地质调查局的国家地震信息中心、国际地震中心和国际数据中心，其主要通过短周期p波也即地震纵波来进行起震点定位。当地震台网密集且震中距较小时，起震点发震时刻以及经纬度和深度坐标均可通过基于走时的定位方法进行确定。但随着震中距增大，初至波走时对起震点深度变化开始变得不敏感。在只有远震台站数据可用的情况下，由于震中距远大于地震埋深，通过直达波走时确定的起震点深度存在很大的不确定性。
3.针对这个问题，现有技术通常使用深度震相约束的方法来解决。但是这些方法一般都无法实现远震的震源破裂起始点深度的快速全自动地重定位。而且有的因为需要人工介入深度震相筛选过程而导致地震目录更新周期非常长。
4.因此，如何准确快速自动地确定出远震起震点深度，是本领域一个亟需解决的技术问题。


技术实现要素：

5.本发明的目的是为了解决现有技术难以快速全自动定位远震起震点深度的技术问题，提出了一种基于深度震相自动匹配的远震起震点深度精定位技术。
6.本发明的技术方案包括以下步骤：
7.s1、根据gcmt目录的远震事件获取可用的远震台站三分量波形数据；
8.s2、对每个远震台站的三分量波形数据进行预处理；
9.s3、对预处理后的三分量波形数据通过自动拾取技术确定直达波到时，并截取直达波震相，然后通过傅里叶变换与相位旋转对所述直达波震相处理形成一个深度震相模板集；
10.s4、将每个远震台站每个分量直达波之后的波形与对应的深度震相模板集进行比对，若相关系数不小于预设阈值，则将对应的深度震相模板作为观测到的深度震相，且将相关系数最大值对应的时刻作为该深度震相的初至时刻；
11.s5、在已知震中位置由浅到深设置多个假设起震点，通过初至正演计算出每个假设起震点的理论直达波震相到时与深度震相到时，并计算两者的走时差；
12.s6、从浅到深依次将各个假设起震点的理论走时差跟观测到的深度震相与直达波震相之间的走时差进行比对，基于设定的误差标准确定出对应起震点成功匹配的深度震相个数与走时差总和，将成功匹配个数大于预设匹配个数且走时差总和最小的对应假设起震点作为确定出的起震点。
13.进一步地，所述步骤s2的预处理具体为将每个远震台站的三个分量波形数据进行去均值、去趋势、消除仪器响应和滤波。
14.进一步地，所述步骤s3具体包括以下分步骤：
15.s31、根据震中和远震台站的位置关系，将每个台站的三分量波形数据旋转变换到垂直、径向和水平三个方向；
16.s32、通过改良的峰态函数震相自动拾取技术截取每个地震台站每个分量的直达波震相波形数据；
17.s33、将每个分量的直达波震相波形数据通过傅里叶正变换变换到频率域，然后从
‑
180
°
到170
°
，以10
°
的间隔，将相位谱做旋转处理，并将所得的36个旋转结果经傅里叶逆变换变换回时间域，构成该分量的一组深度震相模板，对全部地震台站的各个分量进行相同处理，构成一个深度震相模板集。
18.进一步地，所述步骤s4具体包括以下分步骤：
19.s41、将每个远震台站对应的深度震相模板集与直达波后续波形进行比对；
20.s42、确定出每个地震台站对应的深度震相模板与直达波后续波形的相关系数，若相关系数不小于预设阈值，则将对应的深度震相模板作为观测到的深度震相，且将相关系数最大值对应的时刻作为该深度震相的初至时刻，若相关系数小于预设阈值，则继续向后选取扫描波形段并确定出对应相关系数。
21.进一步地，所述步骤s6具体包括以下分步骤：
22.s61、将各个假设起震点到达各个远震台站上的理论深度震相与观测到的深度震相进行比对，得到各假设起震点成功匹配的深度震相个数及走时差绝对值总和；
23.s62、在对所有假设起震点完成上述统计的基础上，将成功匹配个数大于预设匹配个数且走时差总和最小的对应假设起震点作为确定出的起震点。
24.与现有技术相比，本发明具备以下有益效果：
25.(1)本发明通过获取远震发生后可用台站的三个分量波形数据，并对这些台站波形数据进行预处理；然后通过傅里叶变换对直达波进行相位改造，构建深度震相模板；接着在直达波后面的波形搜寻与模板相似的波形段，将相关系数最大值所对应的时刻作为深度震相的初至时刻，计算深度震相与直达波的走时差；在已知震中位置由浅到深给定一系列假设起震点，通过初至正演计算每一假设起震点深度所有的理论深度震相到时，并计算这些深度震相与直达波的走时差；然后将理论走时差与观测走时差进行比对；最后根据地震台站走时差总和最小值对应的假设起震点作为最佳解。上述过程能够在起震点深度定位过程中大幅度节省人力物力，且可避免人为误差，为板块活动研究和地震灾害评价提供重要参考。
26.(2)本发明能够通过自动识别远震台站数据上的深度震相，并通过识别得到的深度震相自动确定出高精度的起震点深度，本发明在识别深度震相和确定起震点深度的过程中无需人工干预，能够完全自动识别深度震相和起震点深度定位。
附图说明
27.图1所示为本发明实施例提供的基于深度震相自动匹配的远震起震点深度精定位技术的流程示意图；
28.图2所示为本技术实施例中多个假设起震点深度示意图。
具体实施方式
29.下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
30.本技术提出了一种基于深度震相自动匹配的远震起震点深度精定位技术，如图1所示为该技术的流程示意图，包括以下步骤：
31.步骤s1、根据gcmt目录的远震事件获取可用的远震台站三分量波形数据。
32.步骤s2、对每个远震台站的三分量波形数据进行预处理。
33.步骤s3、对预处理后的三分量波形数据通过自动拾取技术确定直达波到时，并截取直达波震相，然后通过傅里叶变换与相位旋转对所述直达波震相处理形成一个深度震相模板集。
34.在本技术实施例中，所述步骤s3具体包括：
35.s31、根据震中和远震台站的位置关系，将每个台站的三分量波形数据旋转变换到垂直、径向和水平三个方向；
36.s32、通过改良的峰态函数震相自动拾取技术截取每个地震台站每个分量的直达波震相波形数据；
37.s33、将每个分量的直达波震相波形数据通过傅里叶正变换变换到频率域，然后从
‑
180
°
到170
°
，以10
°
的间隔，将相位谱做旋转处理，并将所得的36个旋转结果经傅里叶逆变换变换回时间域，构成该分量的一组深度震相模板，对全部地震台站的各个分量进行相同处理，构成一个深度震相模板集。
38.在本技术实施例中，所述步骤s2的预处理具体为将每个远震台站的三个分量波形数据进行去均值、去趋势、消除仪器响应和滤波。
39.步骤s4、将每个远震台站每个分量直达波之后的波形与对应的深度震相模板集进行比对，若相关系数不小于预设阈值，则将对应的深度震相模板作为观测到的深度震相，且将相关系数最大值对应的时刻作为该深度震相的初至时刻。
40.在本技术实施例中，所述步骤s4具体包括：
41.s41、将每个远震台站对应的深度震相模板集与直达波后续波形进行比对；
42.s42、确定出每个地震台站对应的深度震相模板与直达波后续波形的相关系数，若相关系数不小于预设阈值，则将对应的深度震相模板作为观测到的深度震相，且将相关系数最大值对应的时刻作为该深度震相的初至时刻，若相关系数小于预设阈值，则继续向后选取扫描波形段并确定出对应相关系数。
43.步骤s5、在已知震中位置由浅到深设置多个假设起震点(如图2所示)，通过初至正演计算出每个假设起震点的理论直达波震相到时与深度震相到时，并计算两者的走时差。
44.在本技术实施例中，所述步骤s5是由浅到深设置多个假设起震点计算对应理论深度震相集及通过初至正演计算具体走时。
45.步骤s6、从浅到深依次将各个假设起震点的理论走时差跟观测到的深度震相与直
达波震相之间的走时差进行比对，基于设定的误差标准确定出对应起震点成功匹配的深度震相个数与走时差总和，将成功匹配个数大于预设匹配个数且走时差总和最小的对应假设起震点作为确定出的起震点。
46.在本技术实施例中，所述步骤s6具体包括以下分步骤：
47.s61、将各个假设起震点到达各个远震台站上的理论深度震相与观测到的深度震相进行比对，得到各假设起震点成功匹配的深度震相个数及走时差绝对值总和；
48.s62、在对所有假设起震点完成上述统计的基础上，将成功匹配个数大于预设匹配个数且走时差总和最小的对应假设起震点作为确定出的起震点。
49.在具体应用场景中，从观测波形匹配出所有可能的深度震相及其初至时刻之后，在已知震中位置基础上由浅到深设置多个假设起震点，如图2所示为本技术实施例中多个假设起震点的示意图。然后通过初至正演计算出每个扫描深度所有的理论深度震相到时，并计算这些深度震相与直达波的理论走时差；然后将这些理论走时差与观测到的深度震相走时差进行比对，得到每个扫描深度上所有台站的深度震相匹配成功的总数与相应的走时差绝对值总和。然后选择深度震相匹配成功总数≥90％，而且走时差绝对值总和最小的对应假设起震点作为确定出的起震点。
50.以发生在智利的一个远震事件为例，本发明通过gcmt目录信息自动搜索和下载该事件的地震数据，然后利用直达波自动扫描得到远震深度震相，再通过这些深度震相的走时信息自动完成震源破裂起始点深度扫描，最后获得的起震点深度解非常接近于isc
‑
ehb的定位结果，验证了本专利技术的可行性和实用性。
51.本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。