[发明专利]一种基于“震源-台站”速度模型的微地震事件定位方法

权利要求书

1.一种基于“震源-台站”速度模型的微地震事件定位方法，其特征在于，包括如下步骤：

(S1)对数量为M的主动震源和数量为N的台站构建“震源-台站”的速度模型组；

(S2)基于到时排序和旅行时差确定最佳主动震源和最佳速度模型；该步骤包括如下流程：

(S21)建立每个主动震源的激发顺序表；

(S22)对待定位的微地震事件激发的台站进行排序；

(S23)确定与微地震事件排序相同的主动震源个数，当个数为0时，初始化p＝1，执行下一步；当个数大于等于2则执行步骤(S25)；当个数为1则执行步骤(S26)；

(S24)查找排序相差为p的主动震源，如果没有，以p+1继续执行该步骤；如果个数大于等于2则执行步骤(S25)，如果个数为1则执行步骤(S26)；

(S25)计算微地震事件的到时时间，比较每个主动震源和待定位微地震事件到各个台站的旅行时间，从而选取和微地震事件旅行时间最为接近的主动震源；

(S26)记录该主动震源和其速度模型为最佳震源和最佳速度模型；

(S3)利用最佳主动震源对应的“震源-台站”速度模型，对待定位的微地震事件进行定位；该步骤中，利用最佳主动震源对应的“震源-台站”速度模型，对待定位的微地震事件进行定位，其计算表达式为公式(3)：

其中，i表示台站编号，共有1、2……N，N个编号；t_i表示事件到编号为i的台站的到时；t₀表示事件的发震时间；v_i表示选定的最佳速度模型中，震源到编号为i的台站的速度；R_i表示事件到编号为i的台站的距离；

或者，利用最佳主动震源对应的“震源-台站”速度模型，对待定位的微地震事件进行定位是通过一个震源新型，结合双差定位思想，形成定位约束公式(4)：

其中，h和j表示事件对应的编号；i表示台站编号，共有1、2……N，N个编号；表示事件h到编号为i的台站的到时；表示事件j到编号为i的台站的到时；v_i表示选定的最佳速度模型中，震源到编号为i的台站的速度；表示事件h到编号为i的台站的距离，表示事件j到编号为i的台站的距离，obs表示括号中的表达式为观测时差。

2.根据权利要求1所述的一种基于“震源-台站”速度模型的微地震事件定位方法，其特征在于，所述步骤(S1)中，任意主动震源m到任意台站n的速度表示为公式(1)：

其中，(x_n、y_n、z_n)表示台站n的空间坐标，(x_m、y_m、z_m)表示主动震源m的空间坐标，t_m表示主动震源m的起震时间，t_n表示台站n接收到震源信号的时间；

按照公式(1)的表示方式对所有的主动震源和台站构建速度模型组，用矩阵表示：