[发明专利]一种地震波S波自动拾取的算法

说明书

技术领域

本发明涉及一种在微震事件波形中进行S波自动拾取的算法方法，更具体地说，本发明涉及一种地震波S波自动拾取的算法。

背景技术

矿山生产安全，需要一套能够针对危险区域预警的系统，能够根据预警区域提前做好防护措施；页岩气的压裂开采需要一套能够指导施工的系统，能够提高施工效率、节省施工材料、提高压裂效果等。都可以通过微震监测系统来辅助完成。微震系统一个主要的功能是它可以提供微震事件3D空间的定位，通过微震事件的3D空间位置和能量来预测矿山危险区域和页岩气压裂效果。微震事件的定位现在主要依据人工拾取的初至和检测区域的走时信息来进行3D定位的。

目前，手动拾取消耗较多的人力和时间，本申请提案旨在设计一种自动拾取的算法，能够自动的完成P波和S波初至的拾取，通过自动拾取和人工纠正相结合的方式，来达到节约人力和时间的目的

发明内容

本发明所要解决的问题是提供一种地震波S波自动拾取的算法，。

为了解决上述问题，本发明采取的技术方案为：一种地震波S波自动拾取的算法，包括如下步骤：

(1)首先对波形自动拾取P波初至。

(2)以P波初至为开始选取固定长度的数据，作为构造求解P-SH-SV 映射系数的协方差矩阵数据；

(3)构造协方差矩阵，求解P-SH-SV映射系数； 

(4)波形数据进行P-SH-SV映射；

(5)在映射后的SH或SV波形中拾取S波初至。

有益效果：通过自动拾取、定位，可以快速的给出微震事件的概要位置，粗略估计预警区域或压裂效果，可以节约人力和时间。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

一种地震波S波自动拾取的算法，包括如下步骤：

(1)首先对波形自动拾取P波初至。

(2)以P波初至为开始选取固定长度的数据，作为构造求解P-SH-SV映射系数的协方差矩阵数据；

(3)构造协方差矩阵，求解P-SH-SV映射系数； 

(4)波形数据进行P-SH-SV映射；

(5)在映射后的SH或SV波形中拾取S波初至。

1.构造协方差矩阵

选取从P波初至开始到S波结束波形的一段数据(最好要包含一段S波)来构造协方差矩阵C

设这段时间内的数据有m个点

北分量的数据为：N1，N2……Nm

东分量的数据为：E1，E2……Em

垂直分量的数据为：Z1，Z2……Zm

可得到下列值。

N‾=(Σi=1mNi)/(m-1)]]>

E‾=(Σi=1mEi)/(m-1)]]>