CN103105624A - 基于数据库技术的纵横波时差定位方法 - Google Patents

Abstract

本发明涉及地球物理勘探微地震监测技术，尤其涉及一种结合数据库技术以在保证计算精度的情况下，极大地提高计算效率的纵横波时差定位方法。本发明的技术方案包括：(1)微地震发生空间区域的建立及网格剖分；(2)利用射线追踪算法将每个网格的P波和S波正演结果写入数据库；(3)结合实际微地震发生空间的方位角，判断Volume_i(x，y，z)是否在微地震事件方位角所确定的范围内(θ_m-θ′，θ_m+θ′)，若其方位角在(θ_m-θ′，θ_m+θ′)范围内，得微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)；(4)根据

在微地震发生空间区域Volume中，所有方位角满足在(θ_m-θ′，θ_m+θ′)内的网格点Volume_i(x，y，z，θ_m±θ′)进行最优震源位置搜索。本发明结合数据库技术进行纵横波时差法的微地震震源定位技术，能满足定位精度并极大地提高速度，该技术在微地震监测定位技术中属于首创技术。

Description

基于数据库技术的纵横波时差定位方法

技术领域

本发明涉及地球物理勘探微地震监测技术，更具体地讲涉及一种结合数据库技术的纵横波时差定位方法。

背景技术

微地震监测就是通过观测、分析生产活动中所产生的微小地震事件来监测生产活动之影响、效果及地下状态的地球物理技术。该方法是在油气藏开采开发的压裂过程中，通过在地面或井中布置的检波器而接收压裂施工中地层因压裂而产生人工裂缝时的微地震波(通称为微地震事件)。通过对接收到的微地震事件(P波和S波)进行一系列数据处理技术，将发生人工裂缝的空间位置反演出来，从而预测人工裂缝的发育趋势及空间展布，为压裂效果的评判及后期的油气田开采开发提供指导。现有的微地震监测定位反演的总体思路主要采取以下步骤：

步骤一：将检波器接收到的微地震信号P波和S波进行波至时间拾取；

步骤二：根据测井、射孔资料建立区域精细速度模型；

步骤三：根据步骤二建立的速度模型，利用非线性反演算法进行一个震源的定位；

步骤四：重复步骤三，最终完成对一个微地震事件定位。

上述方法中，对于一个记录的微地震事件，需要给定一个初始点(即初始空间位置估计)，若初始震源位置与当前记录的微地震事件空间的偏差很大，则在进行上述步骤三中，因为判断标准而导致出现局部极值，使得定位出的震源空间位置与实际微地震事件的空间位置误差较大。

若选用逐点搜索对地下空间位置进行震源扫描，则对每个记录的微地震事件，都会引起数量庞大的正演初至迭代计算，过大的计算量必将导致处理时间的延长。

微地震监测最重要的一点是在压裂施工的现场进行实时监测(在压裂施工的同时，能够在很短的时间内，根据接收到的微地震事件，反演出该微地震事件的空间坐标)，只有这样，才能真正地发挥微地震监测对压裂效果进行实时的评判和指导作用。因此，微地震监测震源位置的反演是计算速度和精度必须兼顾的技术问题。若保证不了反演定位精度，则可能得出错误的评判；若保证不了反演定位的速度，则不能保证微地震监测的实时性。

发明内容

为了克服在现有技术中的上述或其他缺点为此本发明提供一种基于数据库技术的纵横波时差定位方法。

本发明一种基于数据库技术的纵横波时差定位方法，其特征在于，包括如下步骤：(a)将微地震发生空间区域的网格的模拟的P波和S波的初至时间集写入初至时间数据库；(b)对监测的微地震事件进行波至时间拾取，以得到微地震事件的P波和S波的波至时间集；(c)对监测的微地震事件进行方位角求取，结合方位角误差范围及空间网格区域，得到满足方位角范围内的空间网格；(d)用已有的P波和S波的初至时间数据库中，搜索出该方位角的空间网格区域，结合微地震事件的波至事件集，进行微地震事件的空间位置的反演。

当再次发生另一微地震事件时，重复步骤(b)至步骤(d)，从而完成对该另一微地震事件的定位。

本发明的步骤(a)中的P波和S波初至时间集通过如下步骤而得出：(a1)根据微地震监测的测井资料和射孔数据建立微地震事件的P波速度模型V_P(x，y，z)和S波精细速度模型V_S(x，y，z)，其中，x，y，z分别表示笛卡尔坐标系中的x，y，z方向的坐标；(a2)建立微地震事件发生空间区域，对该微地震事件发生空间区域按照x，y，z方向进行网格剖分，从而得到I个网格，其中，第i个网格表示为Volume_i(x，y，z)，I是自然数；(a3)利用射线追踪算法和布置在地面或井中的检波器，模拟出第i个网格Volume_i(x，y，z)所表示的微地震点的P波和S波初至时间集

i是小于或等于I的自然数。

本发明步骤(c)中求取方位角的过程包括：(c1)对第m个事件的第j个检波器所接收到的P波或S波的两个水平分量取一定长度的时窗；(c2)对该时窗内的P波或S波的x和y分量进行偏振分析，可得第m个微地震事件的第j个检波器所表示的方位角；(c3)针对J个检波器的方位角求平均，则得到第m个微地震事件的方位角，其中m是自然数。

当选用P波的方位角时，则求出的P波方位角是微地震事件的真实方位角。

当选用S波的方位角时，则微地震事件的方位角与该S波的求出的方位角有90°的角度差。

本发明的步骤(d)包括：(d1)根据式1进行微地震发生方位区间的判别，

|θ[Volume_i(x，y，z)]-θ_m|＜θ′......式1，其中，θ[Volume_i(x，y，z)]表示网格Volume_i(x，y，z)的方位角，θ_m表示第m个微地震事件的方位角，θ′为预先给定的方位角误差，若网格Voluum_i(x，y，z)的方位角不在(θ_m-θ′，θ_m+θ′)内，则跳过该网格点；(d2)在微地震事件发生空间区域中，对所有方位角满足(θ_m-θ′，θ_m+θ′)的空间网格Volume_i(x，y，z，θ_m±θ′)，根据式2并应用P波和S波的初至时间数据库和波至时间集进行最优微地震空间位置的求取，

$P_m(x,y,z)=\mathrm{Min}\left(\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{\left[\right|(t_{i,j}^P-t_{i,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2\right)$ ......式2

其中，P_m(x，y，z)表示反演的震源位置，

分别表示第j个检波器接收到的第i个震源的P波和S波初至时间；分别表示第j个检波器接收到的第m个微地震事件的P波和S波波至时间，Volume_i(x，y，z，θ_m±θ′)是Volume_i(x，y，z)网格的方位角满足(θ_m-θ′，θ_m+θ′)的限制条件。

反演震源位置P_m(x，y，z)的详细过程如下：

设 $\mathrm{temp}=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{\left[\right|(t_{i,j}^P-t_{i,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$

若 $\mathrm{temp}\≥\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$ 则 $\mathrm{temp}=\underset{j=1}{\overset{J}{\mathrm{\Σ}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2;$

若 $\mathrm{temp}<\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$ 则temp保持不变，

如此循环，直到完成对上述所有满足Volume_i(x，y，z，θ_m±θ′)网格的搜索，最终temp最小时所对应的网格空间Volume_i(x，y，z)即为第m个微地震事件的震源位置。

本发明所应用的射线追踪算法是使用试射法或迭代法，或者是结合试射法和迭代法的算法。

为了实现上述目的，根据本发明的基于数据库技术的纵横波时差定位方法的优点如下：本发明在正式压裂施工前或压裂施工的间隙，利用计算机资源将地下可能的震源区域进行网格划分，用建立好的速度模型，根据射线追踪算法正演地下可能微地震震源的P波和S波初至时间，并以数据库的形式记录P波和S波记录初至时间，在对任意一个记录的微地震事件进行震源位置反演时，只需要读取数据库中的P波和S波记录，而避免了在压裂施工过程中大量重复的利用射线追踪方法进行初至计算，从而保证了定位的速度。

本发明在精度上，进行了两重保证。充分考虑θ_m的误差因素，将θ_m误差范围内进行网格

值最小(即最优)震源空间位置的搜索。避免了θ_m误差导致的P_m(x，y，z)空间误差以及仅使用

准则而导致的最优震源位置的不确定性，最大限度地保证了定位的精度。

附图说明

图1根据本发明的基于初至时间数据库而反演微地震事件的空间位置的流程图；

图2是制作本发明的初至时间数据库的流程图。

具体实施方式

首先，为了便于本领域人员理解，将解释几个名词。并且本发明中使用的i、j、m、M、I均是自然数。

初至时间：地震波波前到达某个观测点时，此点介质的质点开始发生振动的时刻称为波的初至时间，简称初至。

波至时间：有别于初至时间，一个微地震事件的地震波到达观测点的时间。

以下，参照图1和图2详细说明根据本发明的基于数据库技术的纵横波时差定位方法。

首先，建立初至时间数据库，即，将微地震发生空间区域的I个网格的模拟的P波和S波的初至时间集写入初至时间数据库(S101)。

图2示出了根据本发明的建立初至时间数据库的详细过程。为了制作初至时间数据库，按照一定深度给定一个速度值的模式设立地震波的初始速度模型V_k(x，y，z)，再根据射孔数据对初始速度模型V_k(x，y，z)进行校正，从而建立纵波(P波)速度模型V_P(x，y，z)和横波(S波)精细速度模型V_S(x，y，z)(S201)。其中，x，y，z分别表示笛卡尔坐标系中的x、y、z三个坐标。

在压裂施工设计中，会有一个给定的压裂发生区域。将该区域结合地震层位、断层特征进行范围拓展，得到一个更大范围的微地震震源区域。根据前面设定的微地震震源区域建立一个长方体空间，该长方体空间Volume必须将步骤S201建立的微地震空间发生区域全部包含在内，计算出该长方体空间区域Volume的x方向长度为L_x，y方向长度为L_y，z方向长度为L_z。并且沿着x，y，z方向进行网格剖分。若x方向的网格单元为dx，y方向的网格单元为dy，z方向的网格单元为dz，则长度为L_x的x方向可以得到网格数为

则长度为L_y的y方向可以得到网格数为

则长度为L_z的z方向可以得到网格数为则空间区域Volume可以划分的总体网格为I＝N_x·N_y·N_z，则第i个网格表示为Volume_i(x，y，z)(S202)。

$\mathrm{Volume}=\underset{i=1}{\overset{l}{\mathrm{\&Sigma;}}}{\mathrm{Volume}}_i(x,y,z)$ ..................式1

其中，i、I是自然数，Volume表示地层中可能发生震源区域，Volume_i(x，y，z)表示对Volume进行网格剖分后的第i个网格。

假定微地震发生空间区域中第i个网格Volume_i(x，y，z)是一个微地震震源点，结合S201建立的精细P波速度模型V_P(x，y，z)和精细S波速度模型V_S(x，y，z)，利用公知的射线追踪算法RAY(Valume_i(x，y，z)，V_P(x，y，z)，V_S(x，y，z)R(x，y，z))，模拟震源点为Volume_i(x，y，z)的在J个检波器接收到的P波和S波初至时间集

(S203)。

$\mathrm{RAY}({\mathrm{Volume}}_i(x,y,z),V_P(x,y,z),V_S(x,y,z),R(x,y,z))=T_i^{\mathrm{PS}}$ ............式2

$T_i^{\mathrm{PS}}=\{t_{i,1}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^P,t_{i,1}^S,\&CenterDot;\&CenterDot;\&CenterDot;t_{i,j}^S\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^S\}$ ..................式3

其中，Volume_i(x，y，z)是根据步骤S202计算出的第i个震源网格；V_P(x，y，z)、V_S(x，y，z)分别是P波和S波速度模型；R_j(x，y，z)表示第j个检波器的空间坐标；

表示利用射线追踪算法计算的第i个震源网格的P波和S波初至时间集；

分别表示第j个检波器接收到的第i个震源网格的P波和S波初至时间。

最后将I＝N_x·N_y·N_z个网格的P波和S波的初至时间集记录到数据库得到初至时间数据库(S204)。

$\mathrm{Database}\{T_1^{\mathrm{PS}},\&CenterDot;\&CenterDot;\&CenterDot;,T_i^{\mathrm{PS}},\&CenterDot;\&CenterDot;\&CenterDot;T_I^{\mathrm{PS}}\}$ ........................式4

其中 $T_1^{\mathrm{PS}}=\{t_{1,1}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{1,J}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{1,J}^P,t_{1,1}^S,\&CenterDot;\&CenterDot;\&CenterDot;t_{1,j}^S\&CenterDot;\&CenterDot;\&CenterDot;t_{1,J}^S\},$

$T_i^{\mathrm{PS}}=\{t_{i,1}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^P,t_{i,1}^S,\&CenterDot;\&CenterDot;\&CenterDot;t_{i,j}^S\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^S\},$

$T_I^{\mathrm{PS}}=\{t_{I,1}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{I,J}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{I,J}^P,t_{I,1}^S,\&CenterDot;\&CenterDot;\&CenterDot;t_{I,j}^S\&CenterDot;\&CenterDot;\&CenterDot;t_{I,J}^S\}.$

表示I＝N_x·N_y·N_z个微地震源所模拟出的P波和S波初至时间数据库。

接下来，参照回到图1，当发生微地震事件时，通过布置在地面或井中的检波器对监测的微地震事件进行波至时间拾取，可分别得到微地震事件的P波和S波波至时间集，其中，第m个微地震事件的P波和S波波至时间集可被表示为(S102)。其中，

分别表示第j个检波器所接收到的第m个微地震事件的P波和S波波至时间，m是自然数。

当发生微地震事件时，同时也求取该微地震事件的方位角(S103)。以第m个微地震事件为例，对第j个检波器所接收到的P波或S波的两个水平分量取一定长度的时窗(该时窗必须将整个P波或S波包含在内，但窗不能太长，以至干扰信号影响分析的结果)，对该时窗内的P波或S波的x和y分量进行偏振分析，可得第m个微地震事件的第j个检波器所表示的方位角θ_m，j，将J个检波器的方位角θ_m，j求平均，若选用P波的方位角，则求出的P波方位角是微地震事件的真实方位角，若选用S波的方位角，则S波的求出的方位角加上90°，则得到第m个微地震事件的方位角θ_m。利用正演的微地震P波和S波初至时间数据库和波至时间集进行微地震事件空间位置的反演

(S104)。

|θ[Volume_i(x，y，z)]-θ_m|＜θ′......式5

其中，因θ_m在求取时存在一定的误差，因此，将误差范围设定为θ′，可以得到θ_m的真解区域为(θ_m-θ′，θ_m+θ′)，方位角误差θ′通常取小于10°的值。

满足式5时，则表示网格Volume_i(x，y，z)在给定方位角度(θ_m-θ′，θ_m+θ′)范围内，并得微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)，在微地震发生空间区域Volume中只针对其方位角满足(θ_m-θ′，θ_m+θ′)的网格Volume_i(x，y，z)进行最优微地震空间位置的求取，若网格Volume_i(x，y，z)的方位角不在(θ_m-θ′，θ_m+θ′)内，则跳过该网格点。

根据式|θ[Volume_i(x，y，z)]-θ_m|＜θ′，可得到其方位角在给定方位角度(θ_m-θ′，θ_m+θ′)范围内的网格Volume_i(x，y，z)，并得到微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)，并在初至时间数据库中利用所述微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)对应的初至时间进行最优微地震空间位置的求取。在根据式5得到的微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)，根据公式

来反演第m个微地震的地震位置，先求取第i个网格的第j道的P波和S波波至时差

的绝对值和第m个微地震事件的第j道的P波和S波波至时差

的绝对值，设

$\mathrm{temp}=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i,j}^P-t_{i,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$

若 $\mathrm{temp}\&GreaterEqual;\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$ 则 $\mathrm{temp}=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2;$

若 $\mathrm{temp}<\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$ 则temp保持不变。

如此循环，直到对微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)中的所有网格空间进行搜索，最终temp最小时所对应的网格空间Volume_i(x，y，z)即为第m个微地震事件的震源位置。

当又发生一个微地震事件时，利用已经制作好的初至时间数据库，并重复步骤S102至步骤S104，即可完成对该微地震事件的定位。

根据本发明的基于数据库技术的纵横波时差定位方法，本发明在正式压裂施工前或压裂施工的间隙，利用计算机资源将地下可能的震源区域进行网格划分，用建立好的速度模型，根据射线追踪算法正演地下可能微地震震源的P波和S波初至时间，并以数据库的形式记录P波和S波初至时间，在对任意一个记录的微地震事件进行震源位置反演时，只需要读取数据库中的P波和S波记录，而避免了大量重复的利用射线追踪方法进行初至计算，从而保证了定位的速度。对一定的网格数量I＝N_x·N_y·N_z，假定一次压裂监测有m个微地震事件需要进行反演，若在压裂施工过程中，对每个记录的微地震事件进行网格初至的计算，则需要正演初至的次数为I·m＝N_x·N_y·N_z·m次。而用基于数据库技术的纵横波时差定位方法，则仅需要在处理前利用空余的计算机资源对每个网格进行一次正演，从而使得在实际震源定位过程中，完全不需要再进行正演计算，因此，计算效率可提高约N_x·N_y·N_z·m倍。

本发明在实时处理中，根据拾取的微地震事件的P波和S波波至时间，以及该微地震事件的方位角θ_m。考虑到方位角θ_m的求取存在误差，且震源位置离射孔位置距离越远，则方位角θ_m误差导致的震源空间位置越大。首先根据θ_m建立一个微地震发生的角度区域(θ_m-θ′，θ_m+θ′)。然后判断Volume_i(x，y，z)是否在Volume_i(x，y，z)内，可得微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)。并对微地震可能发生空间点Volume_i(x，y，z，θ_m±θ′)中的网格，通过公式

求取取值最小的得网格，该网格即为此次微地震事件的空间位置。本发明在精度上还进行了双重保证。充分考虑方位角θ_m的误差因素，将方位角θ_m误差范围内进行网格

值最小(即最优)震源空间位置的搜索。避免了方位角θ_m误差导致的P_m(x，y，z)空间误差，以及仅使用准则而导致的最优震源位置的不确定性，最大限度地保证了定位的精度。

Claims (9)

1.一种基于数据库技术的纵横波时差定位方法，其特征在于，包括如下步骤：

(a)将通过正演得出的P波和S波的初至时间集写入初至时间数据库，初至时间数据库表示为

其中， $T_i^{\mathrm{PS}}=\{t_{i,1}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^P\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^P,t_{i,1}^S,\&CenterDot;\&CenterDot;\&CenterDot;t_{i,j}^S\&CenterDot;\&CenterDot;\&CenterDot;t_{i,J}^S\},$ $t_{i,j}^P,t_{i,j}^S$ 分别表示第j个检波器接收到的第i个震源的P波和S波初至时间，i、j是自然数；

(b)对监测的微地震事件进行波至时间拾取，建立微地震事件的P波和S波的波至时间集；

(c)求取所述微地震事件的方位角；

(d)利用所述微地震事件的P波和S波的初至时间数据库和波至时间集，进行微地震事件的空间位置的反演。

2.根据权利要求1所述的纵横波时差定位方法，其特征在于，当再次发生另一微地震事件时，重复步骤(b)至(d)，从而完成对该另一微地震事件的定位。

3.根据权利要求1所述的纵横波时差定位方法，其特征在于，步骤(a)中的P波和S波初至时间集通过如下步骤而得出：

(a1)建立微地震事件的P波速度模型和S波精细速度模型；

(a2)建立微地震发生空间区域，针对所述微地震事件发生空间区域，按照笛卡尔坐标系中的x、y、z轴方向进行网格剖分；

(a3)利用射线追踪算法和布置在地面或井中的检波器，模拟出第i个网格所表示的微地震点的P波和S波初至时间集。

4.根据权利要求1所述的纵横波时差定位方法，其特征在于，步骤(c)中求取方位角的过程包括：

(c1)对第m个微地震事件的第j个检波器所接收到的P波或S波的两个水平分量取预定长度的时窗，其中m是自然数；

(c2)对所述时窗内的P波或S波的x和y分量进行偏振分析，取得第m个微地震事件的第j个检波器所表示的方位角；

(c3)求出所有检波器求取的方位角的平均值，得到第m个微地震事件的方位角。

5.根据权利要求4所述的纵横波时差定位方法，其特征在于，当选择利用P波求出的方位角时，将所述方位角确定为所述微地震事件的方位角。

6.根据权利要求4所述的纵横波时差定位方法，其特征在于，当选择利用S波求出的方位角时，将所述方位角加上90°的取值确定为所述微地震事件的方位角。

7.根据权利要求1所述的纵横波时差定位方法，其特征在于，步骤(d)包括：

(d1)根据式1进行微地震发生方位区间的判别，

|θ[Volume_i(x，y，z)]-θ_m|＜θ′......式1，将满足所述式1的网格确定为微地震可能发生空间点，

其中，θ[Volume_i(x，y，z)]表示第i个网格的方位角，θ_m表示第m个微地震事件的方位角，θ′为预先给定的方位角误差；

(d2)当第i个网格的方位角不在(θ_m-θ′，θ_m+θ′)范围内时，则跳过该网格点；

(d3)针对所有微地震可能发生空间点，根据式2应用所述微地震事件的P波和S波的初至时间数据库和波至时间集进行最优微地震空间位置的求取，

$P_m(x,y,z)=\mathrm{Min}\left(\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i,j}^P-t_{i,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2\right)$ ......式2

其中，P_m(x，y，z)表示反演的震源位置，

分别表示第j个检波器接收到的第m个微地震事件的P波和S波波至时间。

8.根据权利要求7所述的纵横波时差定位方法，其特征在于，

反演震源位置P_m(x，y，z)的过程如下：

设 $\mathrm{temp}=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i,j}^P-t_{i,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$

若 $\mathrm{temp}\&GreaterEqual;\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$ 则 $\mathrm{temp}=\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2;$

若 $\mathrm{temp}<\underset{j=1}{\overset{J}{\mathrm{\&Sigma;}}}{\left[\right|(t_{i+1,j}^P-t_{i+1,j}^S)|-|(t_{m,j}^P-t_{m,j}^S)\left|\right]}^2,$ 则temp保持不变，

如此循环，直到针对所有微地震可能发生空间点进行搜索，最终的temp所对应的网格Volume_i(x，y，z)的坐标是第m个微地震事件的震源位置。

9.根据权利要求3所述的纵横波时差定位方法，其特征在于，射线追踪算法是使用试射法或迭代法，或者是结合试射法和迭代法的算法。

Images