[发明专利]一种微地震干涉成像的方法

权利要求书

1.一种微地震干涉成像的方法，其特征在于，所述方法包括：

采用一定大小的网格，将地下模型进行离散，并在所述地下模型中设置n个检波器，n为自然数；

建立地下速度模型，根据所述速度模型，分别计算每个网格点到所述n个检波器的每个检波器中的初至波的走时，获取走时表；

分别获取微地震的主事件在所述n个检波器中的第一记录信息，所述第一记录信息包括所述微地震的主事件的初至波的走时以及波形；

分别获取所述微地震的目标事件在所述n个检波器中的第二记录信息，所述第二记录信息包括所述目标事件的波形；

根据所述第一记录信息以及所述第二记录信息，对所述微地震的主事件和目标事件进行互相关计算，获取互相关值；计算与所述互相关值相对应的包络值，获取互相关包络道集，所述互相关包络道集包括所述主事件和目标事件之间的初至波的走时差值和所述目标事件的激发时刻信息；

在每一个所述网格点，按照特定的时间步长，设定不同的激发时刻值，对所述互相关包络道集进行加时窗干涉成像处理；

对所述加时窗干涉成像处理后的图像进行剖面处理，获取不同剖面图像中的最大值；

对所述不同剖面图像的最大值进行比对，获取目标最大值。

2.根据权利要求1所述的方法，其特征在于，所述目标最大值对应的时间和位置为所述微地震的目标事件的激发时刻和位置。

3.根据权利要求1所述的方法，其特征在于，所述分别获取微地震的主事件在所述n个检波器中的第一记录信息具体包括：获取所述微地震的主事件在所述n个检波器的每个检波器中的初至波的走时以及波形，所述初至波的走时包括直达纵波的走时和直达横波的走时。

4.根据权利要求1所述的方法，其特征在于，根据所述第一记录信息以及所述第二记录信息，对所述微地震的主事件和目标事件进行互相关计算，获取互相关值具体包括：

计算所述互相关值如下式所示：

C(tc)=Σtx(t)y(t-tc)]]>

其中，C(tc)为互相关系数，x(t)为所述第一记录信息中所述微地震的主事件波形的时间离散序列，y(t)为所述第二记录信息中所述微地震目标事件波形的时间离散序列，tc为所述微地震目标事件波形的时间离散序列的时间轴移位值，y(t-tc)为所述微地震目标事件波形的时间离散序列y(t)经过移位长度tc后的离散序列，t∈[1，N]，tc∈[-N，N]，N为每一个所述检波器中的地震波波形的时间长度。

5.根据权利要求1所述的方法，其特征在于，所述计算与所述互相关值相对应的包络值具体包括：

计算所述包络值如下式所示：

E(t)=C2(t)+H2{C(t)}]]>

其中，E(t)为所述微地震的主事件和目标事件进行互相关计算后，获取的与互相关值C(t)相对应的包络值，H²(C(t))为所述互相关值C(t)的希尔伯特变换。

6.根据权利要求1所述的方法，其特征在于，所述对所述互相关包络道集进行加时窗干涉成像处理具体包括：

进行加时窗干涉成像处理具体如下式所示：

其中，M(x，t₀)为不同激发时刻t₀对应的干涉偏移成像值，x为所述网格点的位置向量，i为所述目标事件中的第i个检波器，为所述主事件与目标事件的互相关包络道集，t_w为[-W，W]大小的时窗，i^m为所述主事件中的第i个检波器，m代表所述主事件的位置向量，τ_im和τ_ix分别为所述主事件的第i个检波器中所检测到的初至波的走时和所述目标事件到第i个检波器的理论初至波的走时，t₀为设定的不同的激发时刻值。