CN104749630A - 构建微地震监测速度模型的方法 - Google Patents

Abstract

提供一种构建微地震监测速度模型的方法，包括：选取工区内的多个检波器，并获取对应的初至时间；计算地震波的传播距离；估计出地震波的初始平均速度；基于地震波的初始平均速度、空间径向夹角、射孔点的坐标和当前的约束参数矩阵，得到地震波的平均速度；基于地震波的平均速度和地震波的传播距离，得到地震波的正演走时；计算地震波正演走时的时差的残差；计算偏导数矩阵；基于偏导数矩阵和地震波正演走时的时差的残差，得到约束参数的差值矩阵；在约束参数的差值矩阵满足约束条件时，输出速度模型。上述方法可通过空间建模的方式考虑到地层的各向异性对于地震波传播速度的影响，从而建立了较为复杂的速度模型以提高微地震定位精度。

Description

构建微地震监测速度模型的方法

技术领域

本发明总体来说涉及非常规油气开发的微地震监测技术，更具体地讲，涉及一种基于地表记录中提取的相关信息来构建速度模型的方法。

背景技术

压裂改造是提高非常规油气采收率的关键技术。目前，通过研究压裂改造过程中岩石破裂产生的微地震是监测压裂改造效果的主要技术。利用微地震还可以指导、优化压裂改造，降低对自然环境的污染以及降低诱发地震灾害。井下观测是微地震监测的主要手段，但是，井下观测成本太高，不利于大范围普遍应用。而地表监测的成本较低，使得地表微地震监测成为近几年新兴的微地震监测方式。

速度模型是地震定位的重要先决条件，目前常用的微地震速度模型构建主要有两种方式：一种是根据射孔位置来计算地震波传播的平均速度。这种构建速度模型的方式将地下假设为是一个均匀的介质(即，地下是由同一介质构成)，其计算出的地震波传播的平均速度仅有一个值(即，得到匀速速度模型)。另一种是利用测井资料获得地震波传播的层状速度模型。这种构建速度模型的方式将地下划分为多个地层，假设每个地层具有相同的介质，地震波传播在每个地层中具有相同的速度，因此基于这种方式获得的速度模型是一个关于速度的一维数组(即，水平层状速度模型)。

以上两种方式将地下设想为各向同性介质或横向同性介质，而工区内地震波传播的速度主要受到地层的各向异性影响，因此现有的微地震速度模型构建方式与实际地质结构存在很大的误差，从而导致微地震定位精度不准确。

发明内容

本发明的目的是提出一种构建微地震监测速度模型的方法，以解决现有的速度模型无法真实地反映实际地质结构而导致微地震定位精度差的技术问题。

本发明的一方面提供一种构建微地震监测速度模型的方法，所述方法包括：(a)选取工区内的多个检波器，并获取与所述多个检波器中的每个检波器对应的初至时间；(b)基于射孔点的坐标和所述多个检波器中的每个检波器的坐标，计算所述多个检波器中的每个检波器与射孔点之间地震波的传播距离；(c)对工区的井内速度进行处理，估计出地震波的初始平均速度；(d)基于地震波的初始平均速度、所述多个检波器中的每个检波器相对于射孔点的空间径向夹角、射孔点的坐标和当前的约束参数矩阵，得到所述多个检波器中的每个检波器接收的地震波的平均速度；(e)基于所述地震波的平均速度和所述地震波的传播距离，得到所述多个检波器中的每个检波器接收的地震波的正演走时；(f)从所述多个检波器中选取两组检波器，并计算所述两组检波器之间的初至时间的时差与所述两组检波器之间接收的地震波的正演走时的时差的差值，以得到与所述两组检波器对应的地震波正演走时的时差的残差；(g)基于所述两组检波器接收的地震波的正演走时和当前的约束参数矩阵，得到与所述两组检波器对应的偏导数矩阵；(h)基于所述偏导数矩阵和所述地震波正演走时的时差的残差，得到约束参数的差值矩阵；(i)在约束参数的差值矩阵满足约束条件时，执行步骤(j)：输出速度模型。

可选地，步骤(i)可包括：(i1)计算约束参数的差值矩阵的二范数的平方，并将所述二范数的平方与设定值做比较；(i2)当所述二范数的平方大于设定值时，将当前的约束参数矩阵与约束参数的差值矩阵相加，以更新当前的约束参数矩阵，并返回执行步骤(d)；(i3)当所述二范数的平方不大于设定值时，执行步骤(j)。

可选地，步骤(c)可包括：基于测井曲线来估计地震波的初始平均速度，或者依据经验来人为设定地震波的初始平均速度。

可选地，步骤(d)中，可利用下面的公式计算地震波的平均速度， ${\stackrel{\‾}{v}}_r={\stackrel{\‾}{v}}_0(1+m_{1}X+m_2{X}^2)(1+m_{3}Y+m_4{Y}^2)(1+m_{5}Z+m_6{Z}^2)(1+m_7{\mathrm{\θ}}_r+m_8{\mathrm{\θ}}_r^2)$

其中，为n个检波器中的第r个检波器接收的地震波的平均速度，r∈[1,n]，n为所述多个检波器的数量，为地震波的初始平均速度，X,Y,Z为射孔点的坐标，θ_r为n个检波器中的第r个检波器相对于射孔点的空间径向夹角，m_l为由约束参数形成的矩阵，m_l＝[m₁,m₂,m₃,m₄,m₅,m₆,m₇,m₈]^T，该矩阵中的每个元素对应于一个约束参数。

可选地，在步骤(d)中，在初次执行步骤(d)时，可通过以下步骤获取初始约束参数矩阵，基于与所述多个检波器中的每个检波器对应的初至时间和所述多个检波器中的每个检波器与射孔点之间地震波的传播距离，估算出地震波的平均速度的估计值；基于地震波的平均速度的估计值、射孔点的坐标、所述多个检波器中的每个检波器相对于射孔点的空间径向夹角、地震波的初始平均速度，得到初始约束参数矩阵。

可选地，在步骤(e)中，可利用下面的公式计算所述多个检波器中的每个检波器接收的地震波的正演走时，

$T_r=\frac{d_r}{{\stackrel{\‾}{v}}_r}$

其中，T_r为n个检波器中的第r个检波器接收的地震波的正演走时，r∈[1,n]，n为所述多个检波器的数量，d_r为n个检波器中的第r个检波器与射孔点之间地震波的传播距离，为n个检波器中的第r个检波器接收的地震波的平均速度。

可选地，步骤(f)可包括：(f1)从步骤(a)中获取的初至时间中选取与所述两组检波器中的每个检波器对应的初至时间，并计算所述两组检波器之间的初至时间的时差；(f2)从步骤(e)中得到的地震波的正演走时中选取与所述两组检波器中的每个检波器对应的地震波的正演走时，并计算所述两组检波器之间接收的地震波的正演走时的时差；(f3)计算所述初至时间的时差与所述地震波的正演走时的时差的差值，以得到与所述两组检波器对应的地震波正演走时的时差的残差。

可选地，在步骤(h)中，可利用下面的公式计算约束参数的差值矩阵，

ΔT＝J·Δm_l

其中，ΔT为地震波正演走时的时差的残差，J为偏导数矩阵，Δm_l为由约束参数的差值形成的矩阵，Δm_l＝[Δm₁,Δm₂,Δm₃,Δm₄,Δm₅,Δm₆,Δm₇,Δm₈]^T，该矩阵中的每个元素对应于一个约束参数的差值。

可选地，所述多个检波器的数量为大于等于8。

采用本发明所述的构建微地震监测速度模型的方法，通过采用空间建模的方式考虑到地层的各向异性对于地震波传播速度的影响，从而建立了较为复杂的速度模型以提高微地震定位精度。

附图说明

通过下面结合附图进行的详细描述，本发明的上述和其它目的、特点和优点将会变得更加清楚，其中：

图1是示出根据本发明的示例性实施例的构建微地震监测速度模型的方法的流程图；

图2是示出根据本发明的示例性实施例的与空间径向夹角有关的速度模型的示意图。

具体实施方式

下面，将参照附图详细描述本发明的示例性实施例。

提供参照附图的以下描述以帮助对由权利要求及其等同物限定的本发明的示例性实施例的全面理解。包括各种特定细节以帮助理解，但这些细节仅被视为是示例性的。因此，本领域的普通技术人员将认识到在不脱离本发明的范围和精神的情况下，可对描述于此的示例性实施例进行各种改变和修改。此外，为了清楚和简洁，省略对公知的功能和结构的描述。

总体说来，根据本发明的示例性实施例的构建微地震监测速度模型的方法需首先确定一个工区，并将工区内的压裂井划分为多段，所述多段中的每段压裂井包含多个射孔点，然后从所述多个射孔点中选取初至数据的波形好且信噪比高的射孔点，针对选取的射孔点采用下述方法来构建速度模型。

图1是示出根据本发明的示例性实施例的构建微地震监测速度模型的方法的流程图。

在步骤S10中，选取工区内的多个检波器，并获取与所述多个检波器中的每个检波器对应的初至时间。可选地，可选取n个检波器，n为所述多个检波器的数量，且n∈[8,∞)。

具体说来，可通过下述方式来选取检波器：拾取工区的初至数据，由人工从拾取的初至数据中挑选出具有稳定合理初至时间点的初至数据，然后选取与挑选的初至数据相对应的检波器。这里，一般认为可用于构建速度模型的初至数据越多，则构建出的速度模型相应地更精确，因此，在选取检波器时，可不对选取检波器的数量上限进行限制，将所有满足挑选条件的检波器都选取出来。这里，可利用现有的各种装置和方法来拾取工区的初至数据，例如，可由第三方初至拾取软件进行拾取获得。

在步骤S20中，基于射孔点的坐标和所述多个检波器中的每个检波器的坐标，计算所述多个检波器中的每个检波器与射孔点之间地震波的传播距离。

具体说来，由于在工区***孔点的坐标和n个检波器中的每个检波器的坐标已知，因此可基于射孔点的坐标和n个检波器中的每个检波器的坐标，计算出传播距离。这里，计算得到的传播距离为射孔点与每个检波器之间的直线距离。

在步骤S30中，对工区的井内速度进行处理，估计出地震波的初始平均速度。这里，估计得到的是地震波的整体的初始平均速度，与地表设置的检波器无关，可利用现有的各种方法来估计地震波的初始平均速度。作为示例，可基于测井曲线来估计地震波的初始平均速度，或者可依据经验来人为设定地震波的初始平均速度。

在步骤S40中，基于地震波的初始平均速度、所述多个检波器中的每个检波器相对于射孔点的空间径向夹角、射孔点的坐标和当前的约束参数矩阵，得到所述多个检波器中的每个检波器接收的地震波的平均速度。

这里，地震波的速度对于空间径向夹角比较敏感，且由于射孔点的坐标和n个检波器中的每个检波器的坐标已知，因此可根据射孔点的坐标和n个检波器中的每个检波器的坐标来计算出n个检波器中的每个检波器相对于射孔点的空间径向夹角，进而得到与不同的空间径向夹角对应的地震波的平均速度。

可选地，可利用下面的公式计算地震波的平均速度， ${\stackrel{\‾}{v}}_r={\stackrel{\‾}{v}}_0(1+m_{1}X+m_2{X}^2)(1+m_{3}Y+m_4{Y}^2)(1+m_{5}Z+m_6{Z}^2)(1+m_7{\mathrm{\θ}}_r+m_8{\mathrm{\θ}}_r^2)---\left(1\right)$

公式(1)中，为n个检波器中的第r个检波器接收的地震波的平均速度，r∈[1,n]，为地震波的初始平均速度，X,Y,Z为射孔点的坐标，θ_r为n个检波器中的第r个检波器相对于射孔点的空间径向夹角，m_l为由约束参数形成的矩阵，m_l＝[m₁,m₂,m₃,m₄,m₅,m₆,m₇,m₈]^T，该矩阵中的每个元素对应于一个约束参数。

可选地，在初次执行步骤S40时，可通过以下步骤获取初始约束参数矩阵：基于与所述多个检波器中的每个检波器对应的初至时间和所述多个检波器中的每个检波器与射孔点之间地震波的传播距离，估算出地震波的平均速度的估计值；基于地震波的平均速度的估计值、射孔点的坐标、所述多个检波器中的每个检波器相对于射孔点的空间径向夹角、地震波的初始平均速度，得到初始约束参数矩阵。

具体说来，由于在步骤S10中已经得到与所述多个检波器中的每个检波器对应的初至时间，且在步骤S20中也已经得到所述多个检波器中的每个检波器与射孔点之间地震波的传播距离，因此可利用下面的公式估算出所述多个检波器中的每个检波器接收的地震波的平均速度的估计值，

${\stackrel{\‾}{v}}_r^{\′}=\frac{d_r}{t_r}---\left(2\right)$

公式(2)中，为n个检波器中的第r个检波器接收的地震波的平均速度的估计值，t_r为与n个检波器中的第r个检波器对应的初至时间，d_r为n个检波器中的第r个检波器与射孔点之间地震波的传播距离。

然后，再将地震波的平均速度的估计值射孔点的坐标X,Y,Z、空间径向夹角θ_r、地震波的初始平均速度代入公式(1)中，此时可得到一个关于初始约束参数的方程组，求解该方程组即可得到初始约束参数矩阵。这里，每段压裂井中包含的多个射孔点具有相同的初始约束参数矩阵。

图2是示出根据本发明的示例性实施例的与空间径向夹角有关的速度模型的示意图。

如图2所示，S点表示射孔点，A点、B点、C点为空间上的任意三个点，这里，可理解为是在地表上安放的检波器的位置。作为示例，A点、B点在同一个圆弧上，则位于A点、B点的检波器相对于射孔点S具有相同的空间径向夹角，即，θ_A＝θ_B＝θ₁，位于C点的检波器相对于射孔点S的空间径向夹角为θ₂，由于A点、B点、C点、S点的坐标以及各检波器相对于射孔点S的空间径向夹角均已知，且约束参数矩阵也已知，则可将上述已知的多个参数代入公式(1)来计算出地震波的平均速度。这里，可利用现有的各种方法和装置来获取每个检波器和射孔点的坐标。

在步骤S50中，基于所述地震波的平均速度和所述地震波的传播距离，得到所述多个检波器中的每个检波器接收的地震波的正演走时。

可选地，可利用下面的公式计算所述多个检波器中的每个检波器接收的地震波的正演走时，

$T_r=\frac{d_r}{{\stackrel{\‾}{v}}_r}---\left(3\right)$

公式(3)中，T_r为n个检波器中的第r个检波器接收的地震波的正演走时。这里，地震波的正演走时T_r为与约束参数矩阵有关的矩阵。

在步骤S60中，从所述多个检波器中选取两组检波器，并计算所述两组检波器之间的初至时间的时差与所述两组检波器之间接收的地震波的正演走时的时差的差值，以得到与所述两组检波器对应的地震波正演走时的时差的残差。可选地，在从所述多个检波器中选取检波器时，可从n个检波器中选取任意两个检波器配成一对，共选取多对，每一对的两个检波器分别属于第i组和第j组检波器。优选地，可将彼此之间的距离相对远的两个检波器配成一对。

具体说来，步骤S60可包括：首先从步骤S10中获取的初至时间中选取与所述两组检波器中的每个检波器对应的初至时间，并计算所述两组检波器之间的初至时间的时差，然后从步骤S50中得到的地震波的正演走时中选取与所述两组检波器中的每个检波器对应的地震波的正演走时，并计算所述两组检波器之间接收的地震波的正演走时的时差，再计算所述初至时间的时差与所述地震波的正演走时的时差的差值，以得到与所述两组检波器对应的地震波正演走时的时差的残差。

可选地，可利用下面的公式计算地震波正演走时的时差的残差，

ΔT＝(t_i-t_j)-[T_i-T_j]     (4)

公式(4)中，ΔT为地震波正演走时的时差的残差，t_i为第i组检波器中的所有检波器对应的初至时间形成的矩阵，该矩阵中的每个元素对应于一个检波器对应的初至时间，t_j为与第j组检波器中的所有检波器对应的初至时间形成的矩阵，该矩阵中的每个元素对应于一个检波器对应的初至时间，T_i为与第i组检波器中的所有检波器对应的地震波的正演走时形成的矩阵，该矩阵中的每个元素对应于一个检波器对应的地震波的正演走时，T_j为与第j组检波器中的所有检波器对应的地震波的正演走时形成的矩阵，该矩阵中的每个元素对应于一个检波器对应的地震波的正演走时。

这里，两组检波器中包括的检波器的数量应相同。作为示例，当第i组检波器中包括8个检波器时，矩阵t_i可表示为[t_i1,t_i2,t_i3,t_i4,t_i5,t_i6,t_i7,t_i8]^T，相应地矩阵t_j可表示为[t_j1,t_j2,t_j3,t_j4,t_j5,t_j6,t_j7,t_j8]^T。这里，t_i1与t_j1即为配成一对的两个检波器。

在步骤S70中，基于所述两组检波器接收的地震波的正演走时和当前的约束参数矩阵，得到与所述两组检波器对应的偏导数矩阵。

可选地，可利用下面的公式计算偏导数矩阵J，

$J=[\frac{\∂T_i}{\∂m_l}-\frac{\∂T_j}{\∂m_l}]---\left(5\right)$

在步骤S80中，基于所述偏导数矩阵和所述地震波正演走时的时差的残差，得到约束参数的差值矩阵。

可选地，可利用下面的公式计算约束参数的差值矩阵，

ΔT＝J·Δm_l   (6)

公式(6)中，ΔT为地震波正演走时的时差的残差，Δm_l为由约束参数的差值形成的矩阵，Δm_l＝[Δm₁,Δm₂,Δm₃,Δm₄,Δm₅,Δm₆,Δm₇,Δm₈]^T，该矩阵中的每个元素对应于一个约束参数的差值。

在步骤S90中，计算约束参数的差值矩阵的二范数的平方，并将约束参数的差值矩阵的二范数的平方与设定值做比较。

具体说来，可利用下面的公式计算约束参数的差值矩阵的二范数的平方，

${\left|\right|\mathrm{\Δ}{m}_l\left|\right|}_2^2=\mathrm{\Δ}{m}_1^2+\mathrm{\Δ}{m}_2^2+\mathrm{\Δ}{m}_3^2+\mathrm{\Δ}{m}_4^2+\mathrm{\Δ}{m}_5^2+\mathrm{\Δ}{m}_6^2+\mathrm{\Δ}{m}_7^2+\mathrm{\Δ}{m}_8^2---\left(7\right)$

公式(7)中，为约束参数的差值矩阵的二范数的平方，Δm₁,Δm₂,Δm₃,Δm₄,Δm₅,Δm₆,Δm₇,Δm₈为约束参数的差值矩阵Δm_l中的元素。

当约束参数的差值矩阵的二范数的平方大于设定值时，执行步骤S100：将当前的约束参数矩阵与约束参数的差值矩阵相加，以更新当前的约束参数矩阵，并返回执行步骤S40。

当约束参数的差值矩阵的二范数的平方不大于(即，小于等于)设定值时，执行步骤S110：输出速度模型。

当对选取的射孔点执行上述方法输出速度模型之后，将该输出的速度模型作为该选取的射孔点所在的这段压裂井(即，当前段压裂井)的速度模型。在此基础上，为了确保最终获得的每段压裂井的速度模型的准确性，在针对当前段压裂井输出速度模型后，在对下一段压裂井的数据进行处理时，可使用当前段压裂井的速度模型对下一段压裂井的射孔位置(即，从下一段压裂井包含的多个射孔点中选取的射孔点)进行定位，并将该定位结果与该射孔点的实际坐标进行对比，若误差小于等于设定值，则无需调整速度模型，即，将当前段压裂井输出的速度模型作为下一段压裂井的速度模型，如果误差大于设定值，则需要重复上述步骤，以确定出下一段压裂井的速度模型，从而保证每段压裂井的速度模型的准确性。

根据本发明示例性实施例的构建微地震监测速度模型的方法使用反演的方法进行速度模型的构建，并考虑到各向异性对于地震波的传播速度的影响，建立了相对复杂的速度模型，从而在很大程度上提高了微地震定位精度。

尽管已经参照其示例性实施例具体显示和描述了本发明，但是本领域的技术人员应该理解，在不脱离权利要求所限定的本发明的精神和范围的情况下，可以对其进行形式和细节上的各种改变。

Claims (9)

1.一种构建微地震监测速度模型的方法，所述方法包括：

(a)选取工区内的多个检波器，并获取与所述多个检波器中的每个检波器对应的初至时间；

(b)基于射孔点的坐标和所述多个检波器中的每个检波器的坐标，计算所述多个检波器中的每个检波器与射孔点之间地震波的传播距离；

(c)对工区的井内速度进行处理，估计出地震波的初始平均速度；

(d)基于地震波的初始平均速度、所述多个检波器中的每个检波器相对于射孔点的空间径向夹角、射孔点的坐标和当前的约束参数矩阵，得到所述多个检波器中的每个检波器接收的地震波的平均速度；

(e)基于所述地震波的平均速度和所述地震波的传播距离，得到所述多个检波器中的每个检波器接收的地震波的正演走时；

(f)从所述多个检波器中选取两组检波器，并计算所述两组检波器之间的初至时间的时差与所述两组检波器之间接收的地震波的正演走时的时差的差值，以得到与所述两组检波器对应的地震波正演走时的时差的残差；

(g)基于所述两组检波器接收的地震波的正演走时和当前的约束参数矩阵，得到与所述两组检波器对应的偏导数矩阵；

(h)基于所述偏导数矩阵和所述地震波正演走时的时差的残差，得到约束参数的差值矩阵；

(i)在约束参数的差值矩阵满足约束条件时，执行步骤(j)：输出速度模型。

2.根据权利要求1所述的方法，其中，步骤(i)包括：

(i1)计算约束参数的差值矩阵的二范数的平方，并将所述二范数的平方与设定值做比较；

(i2)当所述二范数的平方大于设定值时，将当前的约束参数矩阵与约束参数的差值矩阵相加，以更新当前的约束参数矩阵，并返回执行步骤(d)；

(i3)当所述二范数的平方不大于设定值时，执行步骤(j)。

3.根据权利要求1所述的方法，其中，步骤(c)包括：基于测井曲线来估计地震波的初始平均速度，或者依据经验来人为设定地震波的初始平均速度。

4.根据权利要求1所述的方法，其中，步骤(d)中，利用下面的公式计算地震波的平均速度，

${\stackrel{\‾}{v}}_r={\stackrel{\‾}{v}}_0(1+m_{1}X+m_2{X}^2)(1+m_{3}Y+m_4{Y}^2)(1+m_{5}Z+m_6{Z}^2)(1+m_7{\mathrm{\θ}}_r+m_8{\mathrm{\θ}}_r^2)$

其中，为n个检波器中的第r个检波器接收的地震波的平均速度，r∈[1,n]，n为所述多个检波器的数量，为地震波的初始平均速度，X,Y,Z为射孔点的坐标，θ_r为n个检波器中的第r个检波器相对于射孔点的空间径向夹角，m_l为由约束参数形成的矩阵，m_l＝[m₁,m₂,m₃,m₄,m₅,m₆,m₇,m₈]^T，该矩阵中的每个元素对应于一个约束参数。

5.根据权利要求4所述的方法，其中，在步骤(d)中，在初次执行步骤(d)时，通过以下步骤获取初始约束参数矩阵，

基于与所述多个检波器中的每个检波器对应的初至时间和所述多个检波器中的每个检波器与射孔点之间地震波的传播距离，估算出地震波的平均速度的估计值；

基于地震波的平均速度的估计值、射孔点的坐标、所述多个检波器中的每个检波器相对于射孔点的空间径向夹角、地震波的初始平均速度，得到初始约束参数矩阵。

6.根据权利要求1所述的方法，其中，在步骤(e)中，利用下面的公式计算所述多个检波器中的每个检波器接收的地震波的正演走时，

$T_r=\frac{d_r}{{\stackrel{\‾}{v}}_r}$

其中，T_r为n个检波器中的第r个检波器接收的地震波的正演走时，r∈[1,n]，n为所述多个检波器的数量，d_r为n个检波器中的第r个检波器与射孔点之间地震波的传播距离，为n个检波器中的第r个检波器接收的地震波的平均速度。

7.根据权利要求1所述的方法，其中，步骤(f)包括：

(f1)从步骤(a)中获取的初至时间中选取与所述两组检波器中的每个检波器对应的初至时间，并计算所述两组检波器之间的初至时间的时差；

(f2)从步骤(e)中得到的地震波的正演走时中选取与所述两组检波器中的每个检波器对应的地震波的正演走时，并计算所述两组检波器之间接收的地震波的正演走时的时差；

(f3)计算所述初至时间的时差与所述地震波的正演走时的时差的差值，以得到与所述两组检波器对应的地震波正演走时的时差的残差。

8.根据权利要求1所述的方法，其中，在步骤(h)中，利用下面的公式计算约束参数的差值矩阵，

ΔT＝J·Δm_l

其中，ΔT为地震波正演走时的时差的残差，J为偏导数矩阵，Δm_l为由约束参数的差值形成的矩阵，Δm_l＝[Δm₁,Δm₂,Δm₃,Δm₄,Δm₅,Δm₆,Δm₇,Δm₈]^T，该矩阵中的每个元素对应于一个约束参数的差值。

9.根据权利要求1所述的方法，其中，所述多个检波器的数量为大于等于8。