[发明专利]一种基于改进程函方程的微地震定位方法

权利要求书

1.一种基于改进程函方程的微地震定位方法，其特征在于，包括以下步骤：

S1、将圈定的监测区域划分为二维离散网格，并在二维离散网格点上布设检波器；

S2、依次将每个检波器作为新震源，采用球面波近似算法计算新震源附近网格点的走时；

S3、根据步骤S2得到的新震源附近网格点的走时，采用多模板快速行进算法计算其余网格点的走时；

S4、利用时间残差函数计算所有网格点走时的函数值，根据计算的函数值定位地震震源位置。

2.根据权利要求1所述的一种基于改进程函方程的微地震定位方法，其特征在于，步骤S1具体包括：

根据微地震监测任务，圈定待监测的区域，建立二维直角坐标系，以微地震事件平面投影到原点的距离作为横坐标，以微地震事件的深度作为纵坐标，将圈定的监测区域划分为二维离散网格，每个网格的大小为d，并在二维离散网格点上以竖井方式布设检波器。

3.根据权利要求1所述的一种基于改进程函方程的微地震定位方法，其特征在于，步骤S2具体包括：

S21、将检波器(x_k,z_k)作为新震源；

S22、根据新震源到达网格点t₁(x,z₁)和网格点t₂(x,z₂)的走时，在两个网格点之间插入一个网格点t₀(x,z₀)，采用Lagrange插值法计算新震源到达网格点t₀(x,z₀)的走时，并计算新震源经过网格点t₀(x,z₀)到达网格点t(x+Δx,z₂)的走时；

S23、最小化网格点t(x+Δx,z₂)的走时，计算网格点t(x+Δx,z₂)的走时对新震源到达网格点t₀(x,z₀)的深度z₀的偏导逼近于0的点，得到新震源到达网格点t₀(x,z₀)的深度z₀，进而求出当网格点t(x+Δx,z₂)的走时最小时，网格点t₀(x,z₀)和网格点t(x+Δx,z₂)的走时大小。