[发明专利]一种基于台站对EDT面概率分布函数的微震定位方法

权利要求书

1.一种基于台站对EDT面概率分布函数的微震定位方法，其特征在于，包括：

对所有台站的微震信号进行滤波预处理；

对滤波预处理后的所述微震信号进行拾取，得到每个所述台站的微震信号初至数据；

根据所述初至数据，计算获得所有所述台站的微震信号信噪比，并对所述微震信号信噪比进行归一化计算；

根据所述初至数据，得到所有台站中任意两个不同所述台站的观测到时差数据，并保留两个不同所述台站总所述微震信号信噪比小的信噪比数据；

对微震监测区域进行三维网格划分，并计算每个网格点到每个所述台站的旅行走时，得到计算走时表数据；

将任意两个不同的所述计算走时表数据相减，得到计算走时差数据；

将所述观测到时差数据、所述微震信号信噪比小的信噪比数据和所述计算走时差数据代入到EDT面概率分布函数中，得到每个网格点的概率值；

最大概率值的网格点位置为微震事件定位结果。

2.如权利要求1所述的基于台站对EDT面概率分布函数的微震定位方法，其特征在于，根据所述初至数据，计算获得所有所述台站的微震信号信噪比包括：

信噪比公式为：

式中：为微震事件i在台站k的信噪比；W₁为从时刻往前取长度为W的时窗内微震信号振幅值；W₂为从时刻往后取长度为W的时窗内背景噪音振幅值；和分别表示为W₁和W₂窗口内波形振幅的均方根值。

3.如权利要求2所述的基于台站对EDT面概率分布函数的微震定位方法，其特征在于，对所述微震信号信噪比进行归一化计算包括：

归一化计算公式为：

式中：为微震事件i在台站k的归一化信噪比；maxⁱ(SNR)为微震事件i的信噪比的最大值。

4.如权利要求3所述的基于台站对EDT面概率分布函数的微震定位方法，其特征在于，根据所述初至数据，得到所有台站中任意两个不同所述台站的观测到时差数据包括：

式中：和分别为微震事件i在台站k和台站l的观测初至；为微震事件i在台站k和台站l的观测到时差，即为微震事件的台站对到时差。

5.如权利要求4所述的基于台站对EDT面概率分布函数的微震定位方法，其特征在于，保留两个不同所述台站总所述微震信号信噪比小的信噪比数据包括：

式中：为保留微震事件i在台站对k，l中微震信号信噪比较小的信噪比数据；为取较小值运算。

6.如权利要求5所述的基于台站对EDT面概率分布函数的微震定位方法，其特征在于，将所述观测到时差数据、所述微震信号信噪比小的信噪比数据和所述计算走时差数据代入到EDT面概率分布函数中，得到每个网格点的概率值包括：

计算公式为：

式中：L^g为网格点g的概率值；TERR为容忍度参数；为台站对k，l的计算走时差与观测到时差的差；为EDT面概率分布函数。

7.如权利要求6所述的基于台站对EDT面概率分布函数的微震定位方法，其特征在于，最大概率值的网格点位置为微震事件定位结果包括：

S(G)＝max(L^g)；

式中：S(G)为概率最大值，G为微震事件定位结果。