[发明专利]一种基于槽波频散特征的微地震震源定位方法

权利要求书

1.一种基于槽波频散特征的微地震震源定位方法，其特征在于，包括：

计算微地震数据每两道之间的互相关记录，将所述互相关记录转换至时频域以得到时频谱；

将监测区域划分为若干个网格，基于各网格与各接收点的距离将各网格点与各接收点的时频谱变换为互相关频散谱；

计算各网格点对应的所有互相关频散谱之间的相似系数，依据所述相似系数确定震源位置。

2.根据权利要求1所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，包括：

对输入微地震数据进行预处理，所述预处理包括对数据进行去野值、去工频干扰、道均衡处理。

3.根据权利要求1所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，将所述互相关记录转换至时频域以得到时频谱包括：采用多次滤波计算每一道互相关记录X_ij(τ)的对应的二维时频谱其中f_k取值为槽波的有效频带范围[f_min,f_max]上以Δf等间隔划分的第k个离散值，i和j表示计算X_ij(τ)的接收点号。

4.根据权利要求3所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，所述二维时频谱基于下式计算得到：

将一道互相关记录X_ij(τ)利用FFT变换到频率域，得到频率域表示的互相关记录

对每个需要计算的离散频率值f_k，利用下式计算出滤波器函数H(f,f_k):

式中，D为滤波器的相对宽度,取值范围为(0，1)，α为高斯函数峰值的锐度参数，取值范围(1-10)；

基于下式对互相关记录进行频域滤波，得到滤波后记录

对执行FFT反变换，得到滤波后的记录X′_ij(f_k,τ)；

计算所有离散频率值f_k即可得到二维时频谱

5.根据权利要求1所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，将监测区域划分为若干个网格，各网格的边长为λ/4，其中，λ为槽波埃里相的波长。

6.根据权利要求1所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，基于各网格与各接收点的距离将各网格点与各接收点的时频谱变换为互相关频散谱包括：

将每个二维时频谱通过变换得到对应的二维频散谱：

其中| |表示模；m表示速度离散序号，v_m＝v_min+(m-1)×Δv，槽波的速度分布范围[v_min:v_max]，速度离散值Δv，离散速度值总数为Nv。

7.根据权利要求6所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，计算各网格点对应的所有互相关频散谱之间的相似系数包括：对所有的二维频散谱Y_ij(f_k,v_m)，基于下式计算它们的零延时互相关系数val：

其中，∑为求和运算符，∏为求积运算符。

8.根据权利要求6所述的一种基于槽波频散特征的微地震震源定位方法，其特征在于，将零延时互相关系数val最大值对应的网格坐标作为微震事件对应的震源位置。