[发明专利]一种震源定位方法、计算机设备及计算机可读存储介质

权利要求书

1.一种震源定位方法，其特征在于，包括如下步骤：

S1、获取各个震波采集器的经纬度坐标信息，转换为X、Y坐标；

S2、根据震源属性选择滤波器参数，对各个震波采集器采集的波形信号进行滤波处理，并提取目标信号到达时刻；

S3、对每个震波采集器，截取包含目标信号到达时刻的、设定时间长度的波形信号，并计算峰度；

S4、根据各个震波采集器对应的峰度筛选可用的震波采集器，构成可用的震波采集器集合，且初始化循环次数；

S5、令循环次数加一，设组合参数为m，在n个可用的震波采集器中提取所有包含M个震波采集器的不同采集器组合；其中n-m≤M≤n，n-m≥4，m取2,3或4；

S6、计算步骤S5中提取的每组采集器组合的系统误差，筛选出系统误差最小的一组采集器组合，并将该组采集器组合计算得到的结果作为当前的最优解保存；

S7、对步骤S6中筛选出的采集器组合，计算其中每个震波采集器的误差，将误差最大的震波采集器从可用的震波采集器集合中删去，令可用的震波采集器个数减一；

S8、判断是否满足结束条件，是则结束循环执行步骤S9，否则返回执行步骤S5；

S9、将当前的最优解作为输出结果，完成震源定位；

所述步骤S3中，计算峰度时，表达式为：

其中，g为峰度，s_i为截取的波形信号在第i个采样点上的信号值，为截取的波形信号的信号值均值，n为截取的波形信号的采样点总数。

2.根据权利要求1所述的震源定位方法，其特征在于：

所述步骤S4中，根据各个震波采集器对应的峰度筛选可用的震波采集器时，将所有峰度大于峰度阈值的震波采集器加入可用的震波采集器集合。

3.根据权利要求1所述的震源定位方法，其特征在于：

所述步骤S6包括如下步骤：

S6-1、对步骤S5中提取的每一组采集器组合，通过采集器组合中各个震波采集器的X、Y坐标及提取到的目标信号到达时刻，求解采集器组合对应的震源位置X、Y坐标，震波传播速度和震动发生时刻，表达式为：

其中，k为采集器组合中的震波采集器个数，e_i表示采集器组合中第i个震波采集器引入的误差，x_i、y_i分别为采集器组合中第i个震波采集器的X、Y坐标，t_i为采集器组合中第i个震波采集器对应的目标信号到达时刻，x、y分别为采集器组合计算得到的震源位置X、Y坐标，v为采集器组合计算得到的震波传播速度，t为采集器组合计算得到的震动发生时刻；

S6-2、比较所有采集器组合，筛选出系统误差最小的一组采集器组合，将该采集器组合对应的震源位置X、Y坐标，震波传播速度和震动发生时刻作为当前计算得到的最优解(x_b，y_b，v_b，t_b)。

4.根据权利要求3所述的震源定位方法，其特征在于：

所述步骤S7中，对步骤S6中筛选出的采集器组合，计算其中每个震波采集器引入的误差时，表达式为：

其中，e_i′为采集器组合中第i个震波采集器引入的误差。

5.根据权利要求1所述的震源定位方法，其特征在于，

所述步骤S8中的结束条件为：

可用的震波采集器个数小于预设个数阈值，或循环次数达到预设次数阈值。

6.根据权利要求5所述的震源定位方法，其特征在于：

所述预设个数阈值N的取值范围为N4。

7.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6中任一项所述震源定位方法的步骤。

8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至6中任一项所述震源定位方法的步骤。