微地震震源定位方法、装置、介质、设备及定位系统

技术领域

本发明涉及微地震监测技术领域，特别是涉及一种微地震震源定位方法、装置、介质、设备及定位系统。

背景技术

微地震震源定位是压裂监测中的重要技术，通过定位破裂震源位置可对水力压裂过程，尤其是裂缝发育形态进行监测，并通过监测结果进一步指导水力压裂施工。

目前的微地震定位方法主要是初至反演定位方法及波形扫描方法。初至反演定位方法，依赖于地震波初至拾取的准确性，但是近年来，由于压裂储层深度增加，震源能量经过地层传播后衰减较大，地震信号信噪比较低的问题较为凸显，初至反演定位方法较难应对地震信号信噪比差的问题，导致定位准确性低。波形扫描方法主要利用波形信息进行微地震震源的定位，但对微地震震源深度方向的定位精度不高。

发明内容

本发明要解决的技术问题是：现有技术中微地震震源的定位准确性低。

为解决上述技术问题，本发明提供了一种微地震震源定位方法、装置、介质、设备及定位系统。

一种微地震震源定位方法，包括：

获取多个检波器检测目标监测区域的微地震信号所得到的微地震记录、纵波旅行时以及横波旅行时；

通过所述纵波旅行时或所述横波旅行时对所述微地震记录进行叠加，获得波形扫描成像剖面；

基于所述纵波旅行时和所述横波旅行时，分析得到纵波初至时刻及横波初至时刻；

利用所述纵波旅行时、所述横波旅行时以及所述纵波初至时刻与所述横波初至时刻的时差，获得时间匹配扫描成像剖面；

将所述波形扫描成像剖面与所述时间匹配扫描成像剖面相乘得到最终成像剖面；

拾取所述最终成像剖面的最大成像值的位置点，得到所述目标监测区域的微地震震源点的位置。

可选地，所述波形扫描成像剖面包括纵波扫描成像值形成的纵波扫描成像剖面和横波扫描成像值形成的横波扫描成像剖面中的一种；所述通过所述纵波旅行时或所述横波旅行时对所述微地震记录进行叠加，获得波形扫描成像剖面，包括：

或/>

其中，

可选地，所述基于所述纵波旅行时和所述横波旅行时，分析得到纵波初至时刻及横波初至时刻，包括：

基于所述纵波旅行时计算得到叠加道的纵波叠加值；

基于所述横波旅行时计算得到叠加道的横波叠加值；

分别根据所述纵波叠加值和所述横波叠加值，通过能量比法获得纵波初至时刻及横波初至时刻。

可选地，所述基于所述纵波旅行时计算得到叠加道的纵波叠加值，包括：

所述基于所述横波旅行时计算得到叠加道的横波叠加值，包括：

其中，

其中，(xo,yo,zo)为所述目标监测区域的中心点；检波器r为所述目标监测区域的正上方检波器。

可选地，所述分别根据所述纵波叠加值和所述横波叠加值，通过能量比法获得纵波初至时刻及横波初至时刻，包括：

分别根据所述纵波叠加值和所述横波叠加值，采用长短时窗能量比法计算得到多个纵波叠加能量比值和多个横波叠加能量比值；

拾取所述纵波叠加能量比值中的最大值，根据所述纵波叠加能量比值中的最大值所对应的时刻得到所述纵波初至时刻；

拾取所述横波叠加能量比值中的最大值，根据所述横波叠加能量比值中的最大值所对应的时刻得到横波初至时刻。

可选地，所述分别根据所述纵波叠加值和所述横波叠加值，采用长短时窗能量比法计算得到多个纵波叠加能量比值和多个横波叠加能量比值，包括：

其中，

为第n个采样点的横波叠加能量比值，/>

可选地，所述时间匹配扫描成像剖面为时间匹配成像值形成的成像剖面；所述利用所述纵波旅行时、所述横波旅行时以及所述纵波初至时刻与所述横波初至时刻的时差，获得时间匹配扫描成像剖面，包括：

根据所述纵波旅行时、所述横波旅行时、所述纵波初至时刻与所述横波初至时刻的时差，计算所述目标监测区域的中心点垂向上各点的时间匹配成像值；

基于所述目标监测区域的中心点垂向上各点的时间匹配成像值，监测所述目标监测区域内其余点的时间匹配成像值。

可选地，所述根据所述纵波旅行时、所述横波旅行时、所述纵波初至时刻与所述横波初至时刻的时差，计算所述目标监测区域的中心点垂向上各点的时间匹配成像值，包括：

其中，

其中，(xo,yo,z)为所述目标监测区域的中心点垂向上的点，检波器r为所述目标监测区域的正上方检波器。

可选地，所述基于所述目标监测区域的中心点垂向上各点的时间匹配成像值，监测所述目标监测区域内其余点的时间匹配成像值，包括：

其中，x

一种微地震震源定位装置，包括：

信息获取模块，用于获取多个检波器检测目标监测区域的微地震信号所得到的微地震记录、纵波旅行时以及横波旅行时；

第一成像模块，用于通过所述纵波旅行时或所述横波旅行时对所述微地震记录进行叠加，获得波形扫描成像剖面；

时刻分析模块，用于基于所述纵波旅行时和所述横波旅行时，分析得到纵波初至时刻及横波初至时刻；

第二成像模块，用于利用所述纵波旅行时、所述横波旅行时以及所述纵波初至时刻与所述横波初至时刻的时差，获得时间匹配扫描成像剖面；

成像处理模块，用于将所述波形扫描成像剖面与所述时间匹配扫描成像剖面相乘得到最终成像剖面；

震源定位模块，用于拾取所述最终成像剖面的最大成像值的位置点，得到所述目标监测区域的微地震震源点的位置。

一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如下步骤：

获取多个检波器检测目标监测区域的微地震信号所得到的微地震记录、纵波旅行时以及横波旅行时；

通过所述纵波旅行时或所述横波旅行时对所述微地震记录进行叠加，获得波形扫描成像剖面；

基于所述纵波旅行时和所述横波旅行时，分析得到纵波初至时刻及横波初至时刻；

利用所述纵波旅行时、所述横波旅行时以及所述纵波初至时刻与所述横波初至时刻的时差，获得时间匹配扫描成像剖面；

将所述波形扫描成像剖面与所述时间匹配扫描成像剖面相乘得到最终成像剖面；

拾取所述最终成像剖面的最大成像值的位置点，得到所述目标监测区域的微地震震源点的位置。

一种定位设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现如下步骤：

获取多个检波器检测目标监测区域的微地震信号所得到的微地震记录、纵波旅行时以及横波旅行时；

通过所述纵波旅行时或所述横波旅行时对所述微地震记录进行叠加，获得波形扫描成像剖面；

基于所述纵波旅行时和所述横波旅行时，分析得到纵波初至时刻及横波初至时刻；

利用所述纵波旅行时、所述横波旅行时以及所述纵波初至时刻与所述横波初至时刻的时差，获得时间匹配扫描成像剖面；

将所述波形扫描成像剖面与所述时间匹配扫描成像剖面相乘得到最终成像剖面；

拾取所述最终成像剖面的最大成像值的位置点，得到所述目标监测区域的微地震震源点的位置。

一种定位系统，包括检波器和上述的定位设备，所述检波器连接所述定位设备。

通过基于纵波旅行时或横波旅行时对微地震记录进行叠加获得波形扫描成像剖面，提高微地震信号的成像效果，对微地震信号信噪比的要求低，通过基于纵波初至时刻与横波初至时刻的时差、纵波旅行时、横波旅行时获得时间扫描成像剖面，充分了利用纵横波信息，可以提高微地震震源深度方向上的定位精度，因此，可以整体提高微地震震源的定位准确性。

附图说明

通过结合附图阅读下文示例性实施例的详细描述可更好地理解本公开的范围。其中所包括的附图是：

图1为一个实施例中微地震震源定位方法的流程示意图；

图2为一个实施例中微地震震源定位装置的结构框图；

图3为应用例一中中等强度地面浅井微地震信号图；

图4为应用例一中纵波扫描成像剖面图；

图5为应用例一中横波扫描成像剖面图；

图6为应用例一中时间匹配扫描成像剖面图；

图7为应用例一中纵波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面图；

图8为应用例一中横波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面图；

图9为应用例二中弱地面浅井微地震信号图；

图10为应用例二中纵波扫描成像剖面图；

图11为应用例二中横波扫描成像剖面图；

图12为应用例二中时间匹配扫描成像剖面图；

图13为应用例二中纵波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面图；

图14为应用例二中横波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，以下将结合附图及实施例来详细说明本发明的实施方法，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。

微地震震源定位方法源自天然地震定位方法，最早的是经典的Gieger定位方法(1912)，其对地震信号初至进行拾取，而后通过基于构建目标函数的反演迭代方法实现地震震源的定位。其后，Douglas(1967)提出的联合震中测定方法(JED)通过多事件定位，降低了速度模型误差对震源定位的影响。而后，国内外专家和学者(Dewey，1972；王椿镛等人，1993；Spence，1980；周仕勇等人，1999)在JED方法基础上进行了改进，形成了相对定位算法。上述方法均可划归初至反演定位方法，其依赖于地震波初至拾取的准确性。

近年来，由于压裂储层深度增加，震源能量经过地层传播后衰减较大，地震信号信噪比较低的问题较为凸显。王维波等人(2012)提出地面微地震发射层析成像(SET)定位算法，利用走时信息对地面微地震信号进行偏移及多道相关处理，降低了噪音对微地震定位的影响。Gajewski等人(2005，2007)通过对微地震记录进行逆时反向延拓波场定位微地震震源。Saenger(2010)通过声波方程进行了微地震事件的逆时定位。Sava(2008，2011)通过干涉成像条件降低了逆时延拓波场中的串扰问题。王晨龙等人(2013)和李振春等人(2014)进行了微地震事件的经地联合逆时定位研究。上述方法可划归波形扫描方法，主要利用波形信息进行微地震震源的定位，但在定位过程中难以进行纵横波信息的结合，对微地震震源深度方向的定位精度不高。

综上所述，目前的微地震定位方法主要是初至反演定位方法及波形扫描方法，初至反演定位方法较难应对地震信号信噪比差的问题，波形扫描方法无法充分利用纵横波信息，两种方法都较难应对微地震微弱信号的三维准确定位。

基于此，本发明提供了一种可以提高微地震震源定位准确性的方案。

实施例一

如图1所示，提供了一种微地震震源定位方法，包括如下步骤：

S110：获取多个检波器检测目标监测区域的微地震信号所得到的微地震记录、纵波旅行时以及横波旅行时。

目标监测区域是需要检测微地震震源的区域，可位于地面浅井，可以包括多个监测点；即，目标监测区域可以为地面浅井的区域。检波器布置于地面浅井之上的地面。其中，纵波旅行时是指震源点到检波器的纵波走时，横波旅行时是指震源点到检波器的横波走时。

S120：通过纵波旅行时或横波旅行时对微地震记录进行叠加，获得波形扫描成像剖面。

具体地，可以是通过纵波旅行时对微地震记录进行叠加，获得波形扫描成像剖面，也可以是通过横波旅行时对微地震记录进行叠加，获得波形扫描成像剖面。

S130：基于纵波旅行时和横波旅行时，分析得到纵波初至时刻及横波初至时刻。

S140：利用纵波旅行时、横波旅行时以及纵波初至时刻与横波初至时刻的时差，获得时间匹配扫描成像剖面。

纵波初至时刻与横波初至时刻的时差是纵波初至时刻与横波初至时刻的差值，比如纵波初至时刻减去横波初至时刻的值。

S150：将波形扫描成像剖面与时间匹配扫描成像剖面相乘得到最终成像剖面。

将波形扫描成像剖面与时间匹配扫描成像剖面相乘，可以是将波形扫描成像剖面各个位置点的成像值与时间匹配扫描成像剖面对应位置点的成像值相乘，得到最终成像剖面各位置点的成像值。

S160：拾取最终成像剖面的最大成像值的位置点，得到目标监测区域的微地震震源点的位置。

上述微地震震源定位方法，通过基于纵波旅行时或横波旅行时对微地震记录进行叠加获得波形扫描成像剖面，提高微地震信号的成像效果，对微地震信号信噪比的要求低，通过基于纵波初至时刻与横波初至时刻的时差、纵波旅行时、横波旅行时获得时间扫描成像剖面，充分了利用纵横波信息，可以提高微地震震源深度方向上的定位精度，因此，可以整体提高微地震震源的定位准确性。

应用上述微地震震源定位方法，在微弱微地震信号条件下，仍能实现微地震震源的三维准确定位，与初至反演定位方法相比，对微地震信号信噪比的要求更低；与波形扫描方法相比，充分利用纵横波信息，提高了微地震震源深度方向上的定位精度。

可选地，波形扫描成像剖面包括纵波扫描成像值形成的纵波扫描成像剖面和横波扫描成像值形成的横波扫描成像剖面中的一种。其中，纵波扫描成像剖面为通过纵波旅行时对微地震记录进行叠加获得的波形扫描成像剖面；横波扫描成像剖面为通过横波旅行时对微地震记录进行叠加获得的波形扫描成像剖面。

步骤S120包括：

或

其中，

利用纵波旅行时对微地震记录进行叠加，获得各个点的纵波扫描成像值，扫描监测区域所有点形成纵波扫描成像剖面；利用横波旅行时对微地震记录进行叠加，获得各个点的横波扫描成像值，扫描监测区域所有点形成横波扫描成像剖面。

可选地，步骤S130包括步骤(a1)至步骤(a3)。

步骤(a1)：基于纵波旅行时计算得到叠加道的纵波叠加值。

具体地，利用纵波旅行时校正微地震记录形成叠加道。

步骤(a2)：基于横波旅行时计算得到叠加道的横波叠加值。

具体地，利用横波旅行时校正微地震记录形成叠加道。

步骤(a3)：分别根据纵波叠加值和横波叠加值，通过能量比法获得纵波初至时刻及横波初至时刻。

具体地，根据纵波叠加值，通过能量比法获得纵波初至时刻；根据横波叠加值，通过能量比法获得横波初至时刻。

可选地，步骤(a1)包括：

步骤(a2)包括：

其中，

可选地，步骤(a3)包括步骤(a31)至步骤(a33)。

步骤(a31)：分别根据纵波叠加值和横波叠加值，采用长短时窗能量比法计算得到多个纵波叠加能量比值和多个横波叠加能量比值。

具体地，根据纵波叠加值，采用长短时窗能量比法计算得到多个纵波叠加能量比值；根据横波叠加值，采用长短时窗能量比法计算得到多个横波叠加能量比值。

步骤(a32)：拾取纵波叠加能量比值中的最大值，根据纵波叠加能量比值中的最大值所对应的时刻得到纵波初至时刻。

例如，纵波初至时刻

步骤(a33)：拾取横波叠加能量比值中的最大值，根据横波叠加能量比值中的最大值所对应的时刻得到横波初至时刻。

例如，横波初至时刻＝记录初始时刻+横波叠加能量比值中的最大值*采样间隔。

实际应用中，可以通过长短时窗能量比法获得相应的长短时窗能量比曲线，拾取相应的长短时窗能量比曲线中的最大值时刻分别获得纵波初至时刻及横波初至时刻。

可选地，步骤(a31)包括：

其中，

可选地，时间匹配扫描成像剖面为时间匹配成像值形成的成像剖面。步骤S140包括步骤(b1)至步骤(b2)。

步骤(b1)：根据纵波旅行时、横波旅行时、纵波初至时刻与横波初至时刻的时差，计算目标监测区域的中心点垂向上各点的时间匹配成像值。

步骤(b2)：基于目标监测区域的中心点垂向上各点的时间匹配成像值，监测目标监测区域内其余点的时间匹配成像值。

进一步地，步骤(b1)包括：

其中，

进一步地，步骤(b2)包括：

其中，x

可选地，以波形扫描成像剖面为纵波扫描成像剖面为例，步骤S150中最终成像剖面的成像公式如下：

其中，S

可以理解，若波形扫描成像剖面为横波扫描成像剖面，则最终成像剖面的成像公式如下：

可选地，步骤S160包括：选择目标监测区域内某点的最终成像剖面的成像值作为起始点，记录所选择的点的位置；依次将目标监测区域内其他点的成像值与起始点的成像值比较大小，若起始点的成像值小，则将大于起始点的成像值的点作为新的起始点，更新起始点的成像值与位置；若起始点的成像值大，则比对下一个点。如此，直到比对完所有的点，挑选出最终成像剖面的成像值中的最大值，最大值的位置点为微地震震源点的位置。

应该理解的是，虽然图1的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。

实施例二

如图2所示，提供了一种微地震震源定位装置，包括信息获取模块210、第一成像模块220、时刻分析模块230、第二成像模块240、成像处理模块250和震源定位模块260。

信息获取模块210用于获取多个检波器检测目标监测区域的微地震信号所得到的微地震记录、纵波旅行时以及横波旅行时。

第一成像模块220用于通过所述纵波旅行时或所述横波旅行时对所述微地震记录进行叠加，获得波形扫描成像剖面。

时刻分析模块230用于基于所述纵波旅行时和所述横波旅行时，分析得到纵波初至时刻及横波初至时刻。

第二成像模块240用于利用所述纵波旅行时、所述横波旅行时以及所述纵波初至时刻与所述横波初至时刻的时差，获得时间匹配扫描成像剖面。

成像处理模块250用于将所述波形扫描成像剖面与所述时间匹配扫描成像剖面相乘得到最终成像剖面。

震源定位模块260用于拾取所述最终成像剖面的最大成像值的位置点，得到所述目标监测区域的微地震震源点的位置。

上述微地震震源定位装置，通过基于纵波旅行时或横波旅行时对微地震记录进行叠加获得波形扫描成像剖面，提高微地震信号的成像效果，对微地震信号信噪比的要求低，通过基于纵波初至时刻与横波初至时刻的时差、纵波旅行时、横波旅行时获得时间扫描成像剖面，充分了利用纵横波信息，可以提高微地震震源深度方向上的定位精度，因此，可以整体提高微地震震源的定位准确性。

可选地，波形扫描成像剖面包括纵波扫描成像值形成的纵波扫描成像剖面和横波扫描成像值形成的横波扫描成像剖面中的一种。其中，纵波扫描成像剖面为通过纵波旅行时对微地震记录进行叠加获得的波形扫描成像剖面；横波扫描成像剖面为通过横波旅行时对微地震记录进行叠加获得的波形扫描成像剖面。

第一成像模块220根据如下公式获得纵波扫描成像值：

或根据如下公式获得横波扫描成像值：

其中，

可选地，时刻分析模块230用于基于纵波旅行时计算得到叠加道的纵波叠加值；基于横波旅行时计算得到叠加道的横波叠加值；分别根据纵波叠加值和横波叠加值，通过能量比法获得纵波初至时刻及横波初至时刻。

可选地，时刻分析模块230根据如下公式计算得到叠加道的纵波叠加值：

时刻分析模块230根据如下公式计算得到叠加道的横波叠加值：

其中，

其中，(xo,yo,zo)为目标监测区域的中心点；检波器r为目标监测区域的正上方检波器。

可选地，时刻分析模块230分别根据纵波叠加值和横波叠加值，通过能量比法获得纵波初至时刻及横波初至时刻，包括：分别根据纵波叠加值和横波叠加值，采用长短时窗能量比法计算得到多个纵波叠加能量比值和多个横波叠加能量比值；拾取纵波叠加能量比值中的最大值，根据纵波叠加能量比值中的最大值所对应的时刻得到纵波初至时刻；拾取横波叠加能量比值中的最大值，根据横波叠加能量比值中的最大值所对应的时刻得到横波初至时刻。

可选地，时刻分析模块230根据如下公式计算得到纵波叠加能量比值和横波叠加能量比值：

其中，

可选地，时间匹配扫描成像剖面为时间匹配成像值形成的成像剖面。第二成像模块240用于：根据纵波旅行时、横波旅行时、纵波初至时刻与横波初至时刻的时差，计算目标监测区域的中心点垂向上各点的时间匹配成像值；基于目标监测区域的中心点垂向上各点的时间匹配成像值，监测目标监测区域内其余点的时间匹配成像值。

进一步地，第二成像模块240根据如下公式计算目标监测区域的中心点垂向上各点的时间匹配成像值：

其中，

进一步地，第二成像模块240根据如下公式获得目标监测区域内其余点的时间匹配成像值：

其中，x

可选地，以波形扫描成像剖面为纵波扫描成像剖面为例，最终成像剖面的成像公式如下：

其中，S

可以理解，若波形扫描成像剖面为横波扫描成像剖面，则最终成像剖面的成像公式如下：

可选地，震源定位模块260选择目标监测区域内某点的最终成像剖面的成像值作为起始点，记录所选择的点的位置；依次将目标监测区域内其他点的成像值与起始点的成像值比较大小，若起始点的成像值小，则将大于起始点的成像值的点作为新的起始点，更新起始点的成像值与位置；若起始点的成像值大，则比对下一个点。如此，直到比对完所有的点，挑选出最终成像剖面的成像值中的最大值，最大值的位置点为微地震震源点的位置。

关于微地震震源定位装置的具体限定可以参见上文中对于微地震震源定位方法的限定，在此不再赘述。上述微地震震源定位装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于定位设备中的处理器中，也可以以软件形式存储于定位设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。需要说明的是，本申请实施例中对模块的划分是示意性的，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式。

实施例三

提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述各实施例中方法的步骤。

上述计算机可读存储介质，由于存储有可以实现上述各实施例方法的计算机程序，同理，可以提高微地震震源定位的准确性。

实施例四

提供了一种定位设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现上述各实施例中方法的步骤。

上述定位设备，由于实现上述各实施例方法的步骤，同理，微地震震源定位的准确性高。

实施例五

提供了一种定位系统，包括检波器和上述实施例中的定位设备，检波器连接定位设备。

上述定位系统，由于采用了前述定位设备，同理，微地震震源定位的准确性高。

为更好地说明本发明的效果，以下结合实际应用例进行说明：

应用例一

应用本发明微地震震源定位方法，进行中等强度地面浅井微地震信号下的扫描定位，如图3所示为中等强度地面浅井微地震信号。图4所示为利用纵波旅行时进行波形扫描叠加获得的纵波扫描成像剖面。图5所示为利用横波旅行时进行波形扫描叠加获得的横波扫描成像剖面。

利用纵波初至时刻与横波初至时刻的时差及纵波旅行时、横波旅行时，得到时间匹配扫描成像剖面，如图6所示。图7所示为通过纵波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面。图8所示为通过横波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面。

可以看出，纵横波时差约束的最终成像剖面在深度方向的定位精度明显优于单纯的纵波或横波扫描成像剖面。

应用例二

应用本发明微地震震源定位方法，进行弱地面浅井微地震信号下的扫描定位，如图9所示为弱地面浅井微地震信号。图10所示为利用纵波旅行时进行波形扫描叠加获得的纵波扫描成像剖面。图11所示为利用横波旅行时进行波形扫描叠加获得的横波扫描成像剖面。

利用纵波初至时刻与横波初至时刻的时差及纵波旅行时、横波旅行时，得到时间匹配扫描成像剖面，如图12所示。图13所示为通过纵波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面。图14所示为通过横波扫描成像剖面及时间匹配扫描成像剖面相乘得到纵横波约束下的最终成像剖面。

同样，纵横波时差约束的最终成像剖面在深度方向的定位精度明显优于单纯的纵波或横波扫描成像剖面。

综上所述，本发明在微弱微地震信号条件下，仍能实现微地震震源的三维准确定位，与初至反演方法相比，对微地震信号信噪比的要求更低；与波形扫描方法相比，充分利用纵横波信息，提高了微地震震源深度方向上的定位精度。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-Only Memory，ROM)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory，RAM)或外部高速缓冲存储器。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static Random Access Memory，SRAM)或动态随机存取存储器(Dynamic Random Access Memory，DRAM)等。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。