微地震事件有效信号拾取方法

1.本技术涉及微地震数据处理技术领域，例如涉及一种微地震事件有效信号拾取方法。


背景技术：

2.目前,微地震监测技术理论都是建立在对有效信号的处理以及分析基础之上的。由于微地震信号通常具有能量弱、信噪比低及持续时间短等特点，早期的一些拾取方法虽然可以大致识别出有效信号，但对噪声敏感度较高，在低信噪比下，其识别的准确性及完整性仍难以取得令人满意的效果。
3.在识别微地震事件有效信号的方法中，长短时窗能量比法是当前较为主流的一种方法。长短时窗能量比法(sta/lta),通过给定一个滑动的长时间窗,在此窗口内再取一个短时间窗，两窗口终点或起始点重合，用短时窗信号平均值(sta)和长时窗信号平均值(lta)之比来反映信号振幅或能量的变化。sta主要反映微震信号的平均值,lta主要反映背景噪音的平均值。在微震信号到达的时间处,sta要比lta变化快，相应的sta/lta值会有一个明显的增加,当其比值大于某个阈值时,即可判定有微震事件发生，从而达到自动检测和拾取微震初至的目的。
4.该算法对于低信噪比的情况，sta/lta方法识别的准确率较低。同时，对该算法起到主要影响因素sta和lta时窗长度、触发阈值、特征函数等因素，容易受人为因素影响，导致微震信号的识别和拾取效果差，使得最终的结果不准确。


技术实现要素：

5.有鉴于此，本发明实施例提供了一种微地震事件有效信号拾取方法，以解决在低信噪比情况下微地震信号识别效果差的问题。
6.为了达到上述目的，本发明实施例提供如下技术方案：
7.一种微地震事件有效信号拾取方法，包括获取微地震事件中的多个微地震信号，并选取其中具有微地震特征的信号作为模板信号，将其余信号设为待识别微地震信号；
8.对所述模板信号进行信号识别，设定微地震有效信号的识别阈值；
9.所述信号识别过程包括：
10.s1对待处理信号进行第一形态变换，去除所述待处理信号中的波峰、波谷尖锐特征，获得平滑处理后的平滑处理信号；
11.s2对所述平滑处理信号进行多尺度放大的第二形态变换，获得所述平滑处理信号的多尺度剖面；
12.s3拾取所述多尺度剖面中显示微地震信号特征的剖面作为识别剖面，并对所述识别剖面进行数据采样，将采样数据与模拟数据比对，获取识别误差；
13.s4根据所述识别误差设定微地震信号的识别阈值；
14.对所述待识别微地震信号，分别执行步骤s1-s2，获取与所述识别剖面同一尺度的
待识别剖面；
15.对所述待识别剖面进行数据采样，并与所述模拟数据比对，获得所述待识别微地震信号的待识别误差；
16.当所述待识别误差在所述识别阈值内时，识别所述待识别微地震信号为微地震事件有效信号。
17.优选地，在上述微地震事件有效信号拾取方法中，所述模板信号为由信噪比大于第一预设值的高信噪比微地震信号中截取；
18.所述模板信号的截取范围为所述高信噪比微地震信号中，具有明显波峰的信号。
19.优选地，在上述微地震事件有效信号拾取方法中，所述第一形态变换为，对所述待处理信号进行多尺度形态学开启运算和多尺度形态学闭合运算。
20.优选地，在上述微地震事件有效信号拾取方法中，所述去除所述待处理信号中的波峰、波谷尖锐特征，包括用所述多尺度形态学开启运算结果的波峰代替所述待处理信号的波峰部分，以及用所述多尺度形态学闭合运算结果的波谷代替所述待处理信号的波谷部分。
21.优选地，在上述微地震事件有效信号拾取方法中，步骤s2包括在所述平滑模板信号上选择一个识别时窗；
22.对所述识别时窗内的信号进行白帽变换处理和黑帽变换处理；
23.计算所述白帽变换处理结果和所述黑帽变换处理结果的平均结果，得到所述平滑处理信号的多尺度剖面。
24.优选地，在上述微地震事件有效信号拾取方法中，步骤s3中所述识别剖面的选取，包括从所述多尺度剖面的多个尺度中，拾取尺度波形趋于稳定的临界尺度剖面，且所述临界尺度剖面上具有高波形的波峰。
25.优选地，在上述微地震事件有效信号拾取方法中，步骤s3中所述获取识别误差为，根据误差判断准则计算所述采样数据与所述模拟数据的识别误差。
26.优选地，在上述微地震事件有效信号拾取方法中，所述模拟数据由弹性波动方程正演模拟而来。
27.优选地，在上述微地震事件有效信号拾取方法中，所述误差判断准则的公式为
28.其中，e(i)为第i道数据的识别误差，t(i)为第i道数据的模拟结果，t(i)为第i道数据的识别，和分别为t(i)和t(i)的平均。
29.本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述微地震事件有效信号拾取方法。
30.在本发明实施例中，从多个微地震信号中选取具有微地震特征的信号作为模板信号，将模板信号经过多尺度形态学运算和多尺度高帽变换处理后，得到模板信号的多尺度剖面并从中选取出识别剖面，根据识别剖面上的识别数据与模拟数据作比对以设定误差阈值；再将其余微地震信号进行多尺度形态学运算和多尺度高帽变换处理后，以模板信号选择识别剖面时的尺度参数和误差阈值去识别其余微地震信号上的有效信号。
31.通过上述识别方法，先对模板信号进行预处理，获得尺度参数和误差阈值，再以此为标准对其余微地震信号进行处理和识别，能显著提升低信噪比情况下，微地震事件有效信号的识别准确度。
附图说明
32.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
33.图1是本发明实施例中微地震事件有效信号拾取方法流程图；
34.图2是本发明实施例中信号识别流程图；
35.图3是本发明实施例中形态学四种基本运算示意图；
36.图4是本发明实施例中多尺度高帽变换示意图；
37.图5是本发明实施例中一种高信噪比合成信号的多尺度剖面图；
38.图6是本发明实施例中一种高信噪比合成信号的第7尺度剖面拾取误差的直方分布图。
具体实施方式
39.本发明实施例提供了一种微地震事件有效信号拾取方法，以解决在低信噪比情况下微地震信号识别效果差的问题。
40.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
41.图1为本发明实施例中微地震事件有效信号拾取方法流程图，该方法包括:
42.步骤101，获取微地震事件中的多个微地震信号，并选取其中具有微地震特征的信号作为模板信号，将其余信号设为待识别微地震信号；
43.步骤102，对所述模板信号进行信号识别，设定微地震有效信号的识别阈值；
44.步骤103，对所述待识别微地震信号，分别执行步骤s1-s2，获取与所述识别剖面同一尺度的待识别剖面；
45.步骤104，对所述待识别剖面进行数据采样，并与所述模拟数据比对，获得所述待识别微地震信号的待识别误差；
46.步骤105，当所述待识别误差在所述识别阈值内时，识别所述待识别微地震信号为微地震事件有效信号。
47.在步骤101中，同一微地震事件信号往往有高信噪比和低信噪比的多组信号，选取一个高信噪比信号作为模板信号，即选取出最典型的微地震信号。对于具有微地震特征的高信噪比微地震事件信号，信号的波形振幅较大，在信号到达时会有较明显的波峰，能肉眼可见地区分于噪声。
48.在步骤102中，信号识别过程如图2所示，为本发明实施例中信号识别流程图，该过
程包括：
49.s1对待处理信号进行第一形态变换，去除所述待处理信号中的波峰、波谷尖锐特征，获得平滑处理后的平滑处理信号；
50.s2对所述平滑处理信号进行多尺度放大的第二形态变换，获得所述平滑处理信号的多尺度剖面；
51.s3拾取所述多尺度剖面中显示微地震信号特征的剖面作为识别剖面，并对所述识别剖面进行数据采样，将采样数据与模拟数据比对，获取识别误差；
52.s4根据所述识别误差设定微地震信号的识别阈值。
53.在步骤s1中，第一形态变换如图3所示，为本发明实施例中形态学四种基本运算示意图，该变换包括：
54.步骤s101，对待处理信号进行多尺度形态学开启运算，采用的计算公式为
[0055][0056]
其中bn为第n个结构元素，f为原始信号，表示膨胀运算，θ表示腐蚀运算。膨胀运算是将结构元素沿信号上方滑动时结构元素顶点所经过的位置的轨迹；同理，腐蚀运算是将结构元素颠倒沿着信号下方滑动时结构元素顶点所经过的位置的轨迹；开启运算是结构元素沿信号上方滑动时结构元素圆心所经过的位置的轨迹；
[0057]
步骤s102，对模板信号进行多尺度形态学闭合运算，采用的公式为
[0058][0059]
同理，闭合运算是将结构元素颠倒沿着信号下方滑动时结构元素圆心所经过的位置的轨迹；
[0060]
步骤s103，在图中红色半椭圆为结构元素，黄色实线为原始信号，绿色实线为膨胀运算结果，紫色实线为腐蚀运算结果，蓝色点线为开启运算结果，黑色点线为闭合运算结果。可以发现，膨胀运算的结果始终位于原始信号上方，腐蚀运算的结果始终位于原始信号的下方，开启和闭合运算大部分结果与原始信号吻合，只是开启运算的结果在波峰附近位于原始信号下方，闭合运算的结果在波谷位置位于信号上方。这一现象表明开启运算可以用来消除原始信号中尖锐的波峰，闭合运算可以消除原始信号中尖锐的波谷。那么采用开启运算结果的波峰部分代替原始信号的波峰，采用闭合运算结果的波谷部分代替原始信号的波谷，就能得到能综合反映模板信号特征的平滑信号。
[0061]
在步骤s2中，第二形态变换如图4所示，为本发明实施例中多尺度高帽变换示意图，该变换包括：
[0062]
步骤s201，在平滑处理信号上选择一个滑动的时窗；
[0063]
步骤s202，对选择的时窗内的信号进行白帽变换处理，采用的公式为
[0064]
步骤s203，对选择的时窗内的信号进行黑帽变换处理，采用的公式为b_hatn(f)＝(f
·bn
)-f；
[0065]
步骤s204，利用白帽变换和黑帽变换之和得到该信号对应的多尺度剖面，即
[0066]
[0067]
在步骤s3中，识别剖面的选取为在多尺度波形图中，从低尺度到高尺度，在前几个尺度中，尺度波形会逐尺度变化，直到某个尺度及其之后，尺度波形趋于稳定，变化不大，且在该尺度下，尺度波形中会出现一个明显的波峰，即为初至波的波峰，区别于其余部分，选择该尺度作为识别剖面；
[0068]
例如，基于一组采用井中接收方式的模型数据，该模型数据模拟了水平放置在巷道内的三分量检波器接收到的煤矿微地震事件，检波器共计两排，每排6个；微地震事件的时间采样间隔为0.5ms，时间采样点为500个，微地震信号子波为主频是50hz的雷克子波，合成数据中各道的有效信号的初至波峰位置见表1.1；
[0069]
表1.1合成数据中各道的有效信号初至的波峰位置
[0070]
trace123456789101112波峰位置(采样点)200215230245260270235250265280295310trace131415161718192021222324波峰位置(采样点)210225240255270285245260275290305320trace252627282930313233343536波峰位置(采样点)220235250265280295240255270285300315
[0071]
向上述合成数据添加高斯噪音，合成信噪比为1db的高信噪比数据，表1.2记录了在高信噪比合成数据中识别到的有效信号波峰的时间采样点位置；
[0072]
表1.2在高信噪比合成数据中识别到的有效信号波峰位置
[0073]
trace123456789101112波峰位置(采样点)199217227244261276235248264276296314trace131415161718192021222324波峰位置(采样点)210226239256272285244259276287304318trace252627282930313233343536波峰位置(采样点)221234253263280297241255273287297313
[0074]
将识别到的高信噪比数据进行步骤s1和步骤s2的形态变换后，生成如图5所示的，为本发明实施例中一种高信噪比合成信号的多尺度剖面图；
[0075]
图5所示的多尺度剖面表明，在前3个尺度剖面中有效信号的特征并不明显；从第4个尺度开始，尺度剖面上的有效信号的特征开始明显；第6个尺度开始，尺度剖面的波形变化不大；于是选择在第7个尺度剖面上进行有效信号识别工作。
[0076]
在步骤s3中，将采样数据与模拟数据比对，获取识别误差为根据误差判断准则计算识别剖面的识别误差，误差判断准则采用的公式为
[0077][0078]
其中e(i)为第i道数据的识别误差，由于是对公式右侧取绝对值，所以都为正，单位为个采样点；t(i)为第i道数据的模拟结果，模拟数据可由弹性波动方程正演模拟而来；t(i)为第i道数据的识别；t和t分别为t(i)和t(i)的平均值。
[0079]
在步骤s4中，根据所述识别误差设定微地震信号的识别阈值；
[0080]
例如，根据上述高信噪比合成数据的识别误差，拾取误差的直方分布图如图6所示，为本发明实施例中一种高信噪比合成信号的第7尺度剖面拾取误差的直方分布图；
[0081]
若将拾取误差阈值设置为5个采样点，则对应相差的时间采样间隔为2.5ms，此差异在图中通常较难察觉，因此若产生这样较小的误差，可以认为该方法的结果是准确的；图6所示最大误差也仅有2个采样点，均小于所设置的阈值，故将误差阈值确定为5个采样点。
[0082]
在步骤103中，对其余所有微地震事件信号，即待识别微地震信号，分别执行步骤s1-s2，获得与其对应的多尺度剖面，再按照模板信号拾取识别剖面时选择的尺度参数，拾取其余待识别微地震信号的待识别剖面；
[0083]
在步骤104中，对待识别剖面进行数据采样，并与之前得到的模拟数据比对，获得所述待识别微地震信号的待识别误差。
[0084]
在步骤105中，将待识别误差与误差阈值进行比较，若待识别误差在所述误差阈值内，则识别其为一个有效信号；若所述误差差值大于所述预设误差阈值，则判定其为非有效信号。
[0085]
本技术的实施例还提供一种计算机可读存储介质，能够实现上述实施例中的微地震事件有效信号拾取方法中全部步骤，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的微地震事件有效信号拾取方法的全部步骤。
[0086]
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。