一种基于“震源-台站”速度模型的微地震事件定位方法与流程

本发明属于微地震监测领域，具体地讲，是涉及一种基于“震源-台站”速度模型的微地震事件定位方法。

背景技术：

震源定位是微地震监测技术相比于其他技术的优势所在，针对微地震震源的定位方法，国内外的学者开展了很多的研究，并且取得了丰富的成果。目前的微地震定位方法主要可以分为基于走时差的定位方法和基于三分量检波器的定位方法两大类(ge,2003a；ge,2003b)，基于三分量传感器的定位方法仅需要一只三分量传感器即可实现对震源的定位，但是该方法需要准确的分离p波和s波并拾取二者的到时，而在微地震监测中，数据普遍噪声大信噪比低，很难实现准确的波场分离，因此该方法在微地震监测中的应用较少，所见报道也不多，如brink等人采用三分量检波器的方法对深部矿井中的岩爆进行监测(brinkando’connor,1983),albright等人将其用于油气储备网络稳定性分析中(albrightandpearson,1982)。

走时差算法是目前微地震震源定位中应用最为广泛的算法，geiger算法是最为经典的走时差算法之一(geiger,1912)，geiger法将非线性问题线性化，以事件的计算到时差和观测到时差的差值作为目标函数，设置初始位置后通过迭代获取最终的定位结果，但该方法对初值的选取有一定的要求，容易发散或陷入局部极值(张唤兰,2014；李楠,2014；田宵,2018)；针对geiger算法的问题，在此基础上也发展了诸如hypo71、hypoinverse、hypocente等多种改进算法(leeetal,1975；lienertetal,1986；nelsonetal,1990)；romney使用距离残差代替时间残差作为目标函数，对震源平面位置、深度和起震时间分别求解，避免了几个参数间相互的影响，但是该方法对震源深度和起震时间的计算并不准确(romney,1957)。2000年，waldhauser和ellsworth提出了双差定位算法(waldhauserandellsworth,2000)，双差定位算法以两个事件之间的到时差作为研究对象，一般来说两个事件中有一个是比较可靠的，以此来降低速度模型对定位结果的影响，提高算法的稳定性和抗噪性，双差定位算法在地震事件重定位，微地震事件定位中取得了广泛的应用(张唤兰,2014；张唤兰,2015；zhang,2003b；徐小明,2008)；丰光亮等人提出了以p波和s波到时差作为目标函数的tdp-s和tdp&s算法，并比较了两种算法性能及其在隧道内应用的效果(fengetal,2017)，不过这两种算法都需要同时拾取p波和s波的到时，在对正演数据进行定位时，可以表现出了较好的可靠性，但是在实际微地震监测中，p波和s波的到时难以被准确拾取，因此会对算法的结果产生一定的影响。

在实际的微地震定位中，应用的最多最广泛的依然是均匀的速度模型，均匀速度模型假设监测区域的介质是各向同性的，各个方向都有相同的传播速度，该速度模型简单易懂，计算方便，但是由于其假设条件的限制，因此在定位中的精度有限，如何实现精准的对微地震震源进行定位是本领域技术人员亟需解决的问题。

技术实现要素：

本发明的目的在于提供一种基于“震源-台站”速度模型的微地震事件定位方法，主要解决现有技术中存在的微地震震源定位精度不高的问题。

为了实现上述目的，本发明采用的技术方案如下：

一种基于“震源-台站”速度模型的微地震事件定位方法，包括如下步骤：

(s1)根据任意数量的主动震源和台站构建任意主动震源到任意台站关于“震源-台站”的速度模型组；

(s2)基于到时排序和旅行时差对速度模型组确定最佳主动震源；

(s3)利用最佳主动震源对应的“震源-台站”速度模型，对待定的微地震事件进行定位。

进一步地，所述步骤(s1)包括如下步骤：

(s11)根据任意震源和任意台站数量求解任意震源到任意台站的速度；

(s12)根据获得的速度求解关于任意主动震源中最佳主动震源到任意台站的速度模型组。

进一步地，所述步骤(s11)中震源数量为m个，台站数量为n个，对任意震源m的关于任意台站n的速度表示为公式(1)：

其中，(xn、yn、zn)表示台站n的空间坐标，(xm、ym、zm)表示震源m的空间坐标，tm表示震源的起震时间，tn表示台站接收到震源信号的时间。

进一步地，所述步骤(s12)中对于任意主动震源si到任意台站tk的速度模型组表示为公式(2)：

其中，表示震源-台站速度模型矩阵，表示震源m到台站n的速度模型。

进一步地，所述步骤(s2)包括如下步骤：

(s21)建立每个主动震源的激发顺序表；

(s22)对待定位的微地震时间激发的台站进行排序；

(s23)确定与微地震时间排序相同的主动震源个数，当个数为0，初始化n＝1，执行下一步；当个数大于等于2则执行步骤(s25)；当个数为1则执行步骤(s26)；

(s24)查找排序相差为n的主动震源，如果没有，以n+1继续执行该步骤；如果个数大于等于2则执行步骤(s25)，如果个数为1则执行步骤(s26)；

(s25)计算微地震时间的到时时间，比较每个主动震源和待定位微地震事件到各个台站的旅行时间，从而选取和微地震事件旅行时间最为接近的主动震源；

(s26)记录该主动震源和其速度模型为最佳震源和最佳速度模型。

进一步地，所述步骤(s3)利用最佳主动震源对应的“震源-台站”速度模型，对待定的微地震事件进行定位，其计算表达式为公式(3)：

其中，i表示事件对应的编号，ti表示事件到台站i的到时，t0表示事件的发震时间，vi表示选定的最佳速度模型中，震源到台站i的速度，ri表示事件到台站i的距离。

具体地，所述步骤(s3)利用最佳主动震源对应的“震源-台站”速度模型，对待定的微地震事件进行定位是通过一个震源新型，结合双差定位思想，形成定位约束公式(4)：

其中，i和j表示事件对应的编号，k为台站编号，表示事件i到台站k的到时，表示事件j到台站k的到时，vi表示选定的最佳速度模型中，震源到台站i的速度，表示事件i到台站k的距离，表示事件j到台站k的距离，obs表示括号中的表达式为观测时差。

与现有技术相比，本发明具有以下有益效果：

本发明提供了一种利用震源确定速度模型的方法，并基于此速度模型对微地震的震源进行定位，从而实现了高精度定位微地震震源的目的。并且本发明改变了原有算法中单一速度模型，本发明中的新速度模型更符合地层信息结构，因此可获得更高的微地震事件定位精度

附图说明

图1为本发明的系统流程框图。

图2为本发明震源及速度模型选取方法的流程图。

图3为本发明以为情况下台站激发顺序与震源位置关系图。

图4为本发明二维平面下区域内有一个主动震源的激发顺序与震源位置关系图。

图5为本发明二维平面下区域内有多个主动震源的激发顺序与震源位置关系图。

图6为本发明二维平面下区域内没有主动震源的激发顺序与震源位置关系图。

图7为本发明隧洞内坐标系、传感器与震源的空间分布图。

图8为本发明设备记录的震源s1的剥线及其频谱图。

图9为本发明震源s1激发后传感器记录的数据波形图。

具体实施方式

下面结合附图和实施例对本发明作进一步说明，本发明的实施方式包括但不限于下列实施例。

实施例

如图1至图9所示，一种基于“震源-台站”速度模型的微地震事件定位方法，包括如下步骤：

(s1)根据任意数量的主动震源和台站构建任意主动震源到任意台站关于“震源-台站”的速度模型组；

(s2)对比待定位事件与主动震源到各台站间的到时排序和旅行时差值，确定最佳主动震源；

(s3)利用最佳主动震源对应的“震源-台站”速度模型，对待定的微地震事件进行定位。

具体过程如下：

假设共有m个已知震源，n个台站，对于任意震源m，到任意台站n的速度可用公式(1)进行表达：

那么对于所有的震源和台站，可建立公式(2)所示的矩阵来表示任意主动震源si到任意台站tk的速度模型组，

其中，矩阵中的每一列代表一个主动震源产生的“震源-台站”速度模型。

绝对震源定位方法的原理基本与单点到时差定位算法类似，但是由于系统中加入了若干数量的主动震源，因此需要从多个震源建立的速度模型组中选取合适的速度模型，带入绝对震源定位方法中计算即可。双差单震源定位方法则需要选取合适的主动震源作为事件对中的已知事件，利用其确定的“震源-台站”速度模型的平均速度作为速度带入双差单震源定位方法中计算即可确定震源位置。

陈法兵等(2016)提出了一种以p波到时排序的方法确定微震震源的定位误差，本文中借鉴其思想用于为震源选取最佳的速度模型。当震源作用或者岩爆发生时，地震波会向四周传播，由近及远的被台站所接收，假设监测系统中共有n个台站(t1,t2,…,tn)，介质为各项同性的均匀介质，地震波波速为v，震源si到台站的距离为ri，因此在均匀速度模型下，触发台站的顺序可由震源到台站的距离密切相关。

在一维的均匀速度模型下，如图3，直线上有两个台站t1和t2，假设震源s激发台站的顺序为t2-t1，则震源一定在直线上t1和t2组成线段中点的右侧，如图3中粗线所示。

假设在二维均匀速度模型下，有三个台站t1，t2，t3和t4，震源h的台站激发顺序t2-t1-t4-t3，图中虚线为各台站连线的中间分界线，阴影部分为按照此激发顺序事件的范围，如图4所示。

假设区域内有且只有一个主动震源，那么就以该主动震源的速度模型作为计算微震事件的速度模型，直接计算即可。

如果区域内有两个或两个以上主动震源存在，如图5所示，假设区域内有k个主动震源，首先利用k个主动震源的平均速度计算微地震事件h的起震时间t0进行第一次求解，则微地震事件h到各个台站的旅行时为公式(5)：

其中，表示事件h到台站tj的旅行时，为台站tj的观测到时，t0为事件的起震时间。

微地震事件h与同区域内主动震源si的到相同台站tj的旅行时差为公式(6)：

其中，表示事件h到主动震源si的旅行时差。

从公式(6)中可以看出，微地震事件h与不同主动震源si的差值dt的大小关系不受t0的影响，可将差值之和最小的主动震源si为选为最佳震源，即公式(7)对应最小值的si：

如果区域内没有主动震源存在，则需要根据前几个传感器的触发顺序，判断触发顺序接近的主动震源。如图6中主动震源s1的激发顺序为t1-t2-t4-t3，主动震源s2的激发顺序为t2-t4-t1-t3，两组震源均只有一个排序与震源h不符，此时，可带入两个速度模型分别计算事件h的位置，与多个震源在同一区域的处理方法相同，选取旅行时差值最小的主动震源为最佳震源。

速度模型选取和震源的选取流程可概括如下，如图2所示：

(s21)建立每个主动震源的激发顺序表；

(s22)对待定位的微地震事件激发的台站进行排序；

(s23)确定与微地震事件排序相同的主动震源个数，个数为0，初始化n＝1，转入(s24)；个数大于等于2转入(s25)；个数为1转入(s26)；

(s24)寻找排序相差n的主动震源，如没有，n+1，继续执行该步骤；如果个数大于等于2，转入(s25)；个数为1转入(s26)；

(s25)计算微地震事件的到时时间，比较每个主动震源和微地震事件到各个台站的旅行时，选取和微地震事件旅行时最为接近的主动震源；

(s26)记此震源和其速度模型为最佳震源和最佳速度模型。

从算法的流程中的旅行时差比较的过程，即利用公式(7)选取最佳震源的过程，选取的主动震源点为与微地震事件激发顺序相同或类似的主动震源中，距离微地震事件最近的主动震源点的速度模型。

在如图7在隧洞实验室宽约2.6米，长约8米，高约3米的隧道内共布设6个震源(s1～s6)和6个传感器(t1～t6)，以正门左侧地面位置为原点，建立如图7所示的坐标系，传感器和震源的坐标分别如表1和表2所示。

表1传感器坐标(m)

表2震源坐标(m)

计算波形模型：试验中采用了三组已知震源，因此共有三个速度模型，针对每个已知震源点，多次敲击计算震源点到各个传感器的平均速度模型，由于监测区域较小，监测中仅直达波波至参与计算，因此速度模型中也仅包含直达波的波速模型。

试验时，利用震源记录仪记录主动震源的起震时间，采用swm-aic算法拾取到时后，相减即可得到震源到传感器的旅行时，用于计算震源到台站的速度模型；每组震源至少保证有10次有效激发，震源-台站速度模型计算结果如表3所示，设备记录的震源波形及频谱如图8所示，传感器记录的数据如图9所示。

表3“震源-台站”速度模型组(m/s)

震源定位：进行震源定位时，由于有六个震源事件，可依据震源激发台站的顺序对其两两分组，分组完成后，将两个事件中的任意一个当做已知，通过其确定的“震源-台站”速度模型采用绝对定位算法或双差单震源定位算法对另外一个事件进行定位。震源激发台站的顺序如表4所示，从表中可知，震源s2和震源s5具有相同的台站激发顺序，因此可将二者作为一组事件来进行定位；另外四组中，s3和s4，s1和s6均是只相差一组便可构成相同的激发震序，因此s3和s4作为一组，s1和s6作为一组。

表4震源触发台站次序表

方法一：

绝对定位算法采用公式(3)为目标函数，其中的n取值为2，即，

公式的速度v采用合适震源处的“震源-台站”速度模型，定位结果和误差如表所示。

表5-5震源定位结果及误差

方法二：

双差单震源定位算法采用公式(4)，取其平方作为目标函数，即

公式的速度v采用合适震源处的均匀“震源-台站”速度模型，事件i的位置选用合适已知震源位置。定位结果和误差如表6所示。

表6震源定位结果及误差

通过使用上述两种方式中的一种方法，即可较为精确的确定微地震震源的位置。

上述实施例仅为本发明的优选实施例，并非对本发明保护范围的限制，但凡采用本发明的设计原理，以及在此基础上进行非创造性劳动而做出的变化，均应属于本发明的保护范围之内。