一种基于主被动震源信号的危险区域检测方法、系统、终端及可读存储介质

摘要

本发明公开了一种基于主被动震源信号的危险区域检测方法、系统、终端及可读存储介质，方法包括：步骤1：获取震源信号；布置微震传感器以及布置主动震源在监测区域内，步骤2：辨识每个震源信号为主动震源信号或是被动震源信号；步骤3：对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻；步骤4：反演出主动震源信号的监测区域波速场、被动震源信号的监测区域波速场；步骤5：基于主动震源信号的监测区域波速场与被动震源信号的监测区域波速场进行比对进行采空区判定。本发明该方法结合主动震源和被动震源进行定位，提高了空区判定的精度。

权利要求

1.一种基于主被动震源信号的危险区域检测方法，其特征在于：包括如下步骤：步骤1：获取震源信号；

其中，布置微震传感器以及在监测区域内布置主动震源，其中，微震传感器接收主动、被动震源信号，并记录每个震源信号触发微震传感器的时刻以及记录主动震源位置和每个主动震源信号的释放时刻；

步骤2：辨识每个震源信号为主动震源信号或是被动震源信号；

步骤3：对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻；

步骤4：基于主动震源位置、被动震源位置、主动震源信号的释放时刻、被动震源信号的释放时刻、触发微震传感器的时刻反演出主动震源信号的波速场、被动震源信号的波速场；

步骤5：基于主动震源信号的监测区域波速场和基于被动震源信号的监测区域波速场进行比对进行采空区判定；

其中，若两种波速场存在相斥情况，则移动主动震源重新反演出主动震源信号的波速场，且以主动震源信号的反演结果为准；

步骤5中进行采空区判定的过程如下：若所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场中均无低波速场，则波速反演区域无采空区；

若所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场中存在相斥情况，则移动主动震源信号重新反演出主动震源信号的波速场，再将重新反演出主动震源信号的监测区域波速场与被动震源信号的监测区域波速场进行比较，若反演结果依旧相斥，以重新反演出主动震源信号的监测区域波速场为准；

若所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场均表示存在采空区，则存在采空区。

2.根据权利要求1所述的方法，其特征在于，步骤4中所述主动震源信号的波速场、被动震源信号的波速场的反演过程如下：将监测区域进行网格划分，并采用路径搜索算法追踪主动震源、被动震源与微震传感器之间震源信号的实际传播路径；

再基于实际传播路径以及主动震源信号的释放时刻、被动震源信号的释放时刻以及触发微振传感器的时刻构建监测区间内主动震源信号、被动震源信号在实际传播路径上的走时误差；

基于所述走时误差反演出主动震源信号、被动震源信号分别在待测区域内每个网格的波速得到所述主动震源信号的波速场、被动震源信号的波速场。

3.根据权利要求2所述的方法，其特征在于：构建的监测区间内主动震源信号、被动震源信号在实际传播路径上的走时误差如下所示：式中， 分别表示监测到的所有主动震源信号、所有被动震源信号在实际传播路径上的走时误差，N1、N2分别为监测到的主动震源信号、被动震源信号的个数，K为监测区域内划分的网格总数， 表示微震传感器i和主动震源j或被动震源j的实际传播路径在一个xyz网格内的传播距离、波速，T′ij表示微震传感器i和主动震源j或被动震源j在传播路径上的实际走时，M为微震传感器的总个数。

4.根据权利要求1所述的方法，其特征在于：步骤3中对任意一个被动震源信号进行定位的过程如下：

首先，基于被动震源信号构建每任意两个微震传感器之间的到时差得到到时差方程组；

然后，基于到时差方程组构建到时差与真实值的偏差方程组；

接着，基于偏差方程组求解偏差平方和最小值时对应的被动震源的位置以及地震波在空区与岩体中传播的速度V；

最后，利用被动震源的位置以及地震波在空区与岩体中传播的速度V计算出被动震源信号的释放时刻。

5.根据权利要求1所述的方法，其特征在于：步骤2中利用训练好的神经网络机器学习模型辨识每个震源信号为主动震源信号或是被动震源信号，其中，所述神经网络机器学习模型的构建过程如下：

首先，利用人工辨识的主动震源信号和被动震源信号建立相应的波形数据库，并分别选取M1组主动震源信号波形和M2组被动震源信号波形作为波形图像样本数据库；

然后，对选取的波形图像样本进行预处理得到每个波形图像对应的二维矩阵；

接着，将每个波形图像对应的二维矩阵以及二维矩阵对应的训练标签输入神经网络机器学习模型中进行训练，其中，依据波形图形为主动震源信号或被动震源信号进行编码标记并作为训练标签。

6.根据权利要求1所述的方法，其特征在于：任意两个微震传感器之间的直线距离在

200m以内。

7.一种基于权利要求1‑6任一项所述方法的检测系统，其特征在于：包括：震源信号获取模块：用于获取主动震源信号和被动震源信号以及每个震源信号触发微震传感器的时刻以及主动震源位置和每个主动震源信号的释放时刻；

辨识模块：用于辨识每个震源信号为主动震源信号或是被动震源信号；

被动震源信号定位模块：用于对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻；

波速场反演模块：用于反演出主动震源信号的监测区域波速场、被动震源信号的监测区域波速场；

采空区判定模块：用于基于所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场进行比对进行采空区判定。

8.一种检测终端，其特征在于：包括处理器和存储器，所述存储器存储了计算机程序，所述处理器调用所述计算机程序执行权利要求1‑6任一项所述方法的步骤。

9.一种可读存储介质，其特征在于：存储了计算机程序，所述计算机程序被处理器调用执行权利要求1‑6任一项所述方法的步骤。

说明书

一种基于主被动震源信号的危险区域检测方法、系统、终端及

可读存储介质

技术领域

[0001] 本发明属于微震监测预警技术领域，具体涉及一种基于主被动震源信号的危险区域检测方法。

背景技术

[0002] 随着浅部资源的逐年减少和枯竭，国内外矿产资源的开采正向深部发展，处于高应力状态下的岩体受机械凿岩、爆破作业等工程扰动的影响，导致周围岩体出现应力集中、
应力释放等应力重分布现象，使得原本稳定的岩体进入不稳定状态，极大地威胁着施工人
员和设备的安全。工程扰动会使岩体变形甚至催生岩体内部裂纹的扩展，在这个过程中岩
体内部储存的弹性势能会在岩体非弹性变形的过程中以波的形式释放，该能量的释放最终
会导致岩体工程失稳并在岩体内部产生微震事件。微震监测技术可以实时采集岩体内部裂
纹萌生、发育、贯通时释放的地震波信息，并通过信号辨识、震源定位及其破坏机制反演等
手段获取丰富的震源信息，并对监测区域内可能存在的高风险区域进行支护和人员撤离。
因此，微震监测技术已成为深部资源开采中有效的监测手段之一。

[0003] 工程扰动对岩体工程失稳的影响根据作用机理大致可以分为两类：一类是通过工程扰动周围岩体应力的重分布间接影响岩体原生裂隙、断层处的应力分布状态典型的如断
层滑移型大震级事件，另一种则为工程扰动直接作用于致灾区域例如采空区塌陷等等。因
此，地下空区的存在极大地影响了矿山开采作业的效率与安全，为了有效探清地下空区分
布位置及其大小并合理安排作业进行，传统方法多采用钻孔探测，对疑似空区位置打孔探
测，根据岩芯分析地质结构构造以作出评价，但该种方法不但要耗费大量的人力物力和财
力，而且需要花费很长的时间。

[0004] 目前地震学中为解决上述问题采取的技术手段主要是波速场成像技术，通过对震源传播路径以及波速进行反演，求解监测区域内波速场模型，以便对危险区域进行圈定，保
证矿山生产的安全进行。按照震源类型的不同，波速场成像技术可以分为主动震源波速成
像技术和被动震源波速成像技术，主动震源的使用虽然保证了震源位置的精确性，但是在
释放主动震源的过程中要停止一切生产活动，极大地限制了波速成像技术在现场中的应
用，另一方面爆破作业、机械凿岩等被动震源虽然能够保证信号的连续释放，但在复杂岩体
结构下缺乏对被动震源精确定位的方法，难以保证反演的精度。因此，如何既保证定位精度
又能实现长期连续检测危险区域成为了深地资源开采中监控防灾的一大难题。

发明内容

[0005] 本发明的目的是提供一种结合主动震源和被动震源的危险区域监测方法，保证了震源定位的精度以及空区判定的精度。

[0006] 本发明提供了一种基于主被动震源信号的危险区域检测方法，包括如下步骤：

[0007] 步骤1：获取震源信号；

[0008] 其中，布置微震传感器以及在监测区域内布置主动震源，其中，微震传感器接收主动、被动震源信号，并记录每个震源信号触发微震传感器的时刻以及记录主动震源位置和
每个主动震源信号的释放时刻；

[0009] 步骤2：辨识每个震源信号为主动震源信号或是被动震源信号；

[0010] 步骤3：对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻；

[0011] 步骤4：基于主动震源位置、被动震源位置、主动震源信号的释放时刻、被动震源信号的释放时刻、触发微震传感器的时刻反演出主动震源信号的监测区域波速场、被动震源
信号的监测区域波速场；

[0012] 步骤5：基于所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场进行比对进行采空区判定；

[0013] 其中，若两种波速场存在相斥情况，则移动主动震源重新反演出主动震源信号的监测区域波速场，且以主动震源信号的反演结果为准。

[0014] 进一步优选，步骤4中所述主动震源信号的监测区域波速场、被动震源信号的监测区域波速场的反演过程如下：

[0015] 将监测区域进行网格划分，并采用路径搜索算法追踪主动震源、被动震源与微震传感器之间震源信号的实际传播路径；

[0016] 再基于实际传播路径以及主动震源信号的释放时刻、被动震源信号的释放时刻以及触发微振传感器的时刻构建监测区间内主动震源信号、被动震源信号在实际传播路径上
的走时误差；

[0017] 基于所述走时误差反演出主动震源信号、被动震源信号分别在待测区域内每个网格的波速得到所述主动震源信号的监测区域波速场、被动震源信号的监测区域波速场。

[0018] 本发明反演波速场的过程是推演出实际传播路径，再基于传播路径中穿过每个网格所需时间之和与实际震源走时之间的误差函数，达到反演出所有传播路径通过的网格波
速场的效果，这是基于根据路径搜索算法可以得到震源地震波的传播路径，通过多条传播
路径对穿过的网格中的波速表达出地震波在每个网格中的走时，将每个网格中的走时相加
与实际震源走时之间构成一个误差函数，满足误差函数值最小的网格波速场即为求得的监
测区域波速场。

[0019] 进一步优选，构建的监测区间内主动震源信号、被动震源信号在实际传播路径上的走时误差如下所示：

[0020]

[0021] 式中， 分别表示监测到的所有主动震源信号、所有被动震源信号在实际传播路径上的走时误差，N1、N2分别为监测到的主动震源信号、被动震源信号的个数，K为
监测区域内划分的网格总数， 表示微震传感器i和主动震源j或被动震源j的实
际传播路径在一个xyz网格内的传播距离、波速，T′ij表示微震传感器i和主动震源j或被动
震源j在传播路径上的实际走时，M为微震传感器的总个数。

[0022] 进一步优选，步骤5中进行采空区判定的过程如下：

[0023] 若所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场中均无低波速场，则波速反演区域无采空区；

[0024] 若所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场中存在相斥情况，则移动主动震源信号重新反演出主动震源信号的监测区域波速场，再将
重新反演出主动震源信号的监测区域波速场与被动震源信号的监测区域波速场进行比较，
若反演结果依旧相斥，以重新反演出主动震源信号的监测区域波速场为准；

[0025] 若所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场均表示存在采空区，则存在采空区，并在施工区域和采空区连接处建立临时支护措施。

[0026] 进一步优选，步骤3中对任意一个被动震源信号进行定位的过程如下：

[0027] 首先，基于被动震源信号构建每任意两个微震传感器之间的到时差得到到时差方程组；

[0028] 然后，基于到时差方程组构建到时差与真实值的偏差方程组；

[0029] 接着，基于偏差方程组求解偏差平方和最小值时对应的被动震源的位置以及地震波在空区与岩体中传播的速度V；

[0030] 最后，利用被动震源的位置以及地震波在空区与岩体中传播的速度V计算出被动震源信号的释放时刻。

[0031] 进一步优选，步骤2中利用训练好的神经网络机器学习模型辨识每个震源信号为主动震源信号或是被动震源信号，其中，所述神经网络机器学习模型的构建过程如下：

[0032] 首先，利用人工辨识的主动震源信号和被动震源信号建立相应的波形数据库，并分别选取M1组主动震源信号波形和M2组被动震源信号波形作为波形图像样本数据库；

[0033] 然后，对选取的波形图像样本进行预处理得到每个波形图像对应的二维矩阵；

[0034] 接着，将每个波形图像对应的二维矩阵以及二维矩阵对应的训练标签输入神经网络机器学习模型中进行训练，其中，依据波形图形为主动震源信号或被动震源信号进行编
码标记并作为训练标签。

[0035] 进一步优选，任意两个微震传感器之间的直线距离在200m以内。

[0036] 第二方面，本发明还提供一种基于上述方法的检测系统，包括：

[0037] 震源信号获取模块：用于获取主动震源信号和被动震源信号以及每个震源信号触发微震传感器的时刻以及主动震源位置和每个主动震源信号的释放时刻；

[0038] 辨识模块：用于辨识每个震源信号为主动震源信号或是被动震源信号；

[0039] 被动震源信号定位模块：用于对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻；

[0040] 波速场反演模块：用于反演出主动震源信号的监测区域波速场、被动震源信号的监测区域波速场；

[0041] 采空区判定模块：用于将所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场进行比对进行采空区判定。

[0042] 第三方面，本发明提供一种检测终端，包括处理器和存储器，所述存储器存储了计算机程序，所述处理器调用所述计算机程序执行上述检测方法的步骤。

[0043] 第四方面，本发明提供一种可读存储介质，存储了计算机程序，所述计算机程序被处理器调用执行上述检测方法的步骤。

[0044] 有益效果

[0045] 本发明提供的一种基于主被动震源信号的危险区域检测方法，其将主动震源与被动震源巧妙地的结合在一体，充分利用主动震源和被动震源监测的优势，同时规避各自的
缺点，一方面，由于主动震源信号可以自动识别出来，因此使用本发明所述方法时，主动震
源释放过程中并不需要停止生产活动，二方面，本发明同时考虑了主动震源的监测区域波
速场和被动震源的监测区域波速场，由于主动震源与被动震源之间信号的不可兼容性，导
致在目前的实际应用中以单类震源作为信号的输入，但主动震源和被动震源各有其自身的
优势，同时两类信号反演出的波速场又可互相验证，为了利用两类信号并进行正确的反演，
本发明克服了对两类信号辨识的技术障碍，同时实现了按照网格对监测区域内的波速场进
行主被动震源信号反演，明显提高监测结果的可靠性。

附图说明

[0046] 图1为一种基于主被动信号的危险区域检测方法实施场景图；

[0047] 图2为一种基于主被动信号的危险区域检测方法原理图。

具体实施方式

[0048] 本发明提供的一种基于主被动震源信号的危险区域检测方法是结合主动震源和被动震源进行采空区检测的手段，具体是在矿山采场周围附近布置数个微震传感器使监测
区域被微震传感器网络所覆盖；其中，在监测区域中心位置处放置主动震源，定时对测点产
生冲击释放微震波；其次由于主动震源释放时伴随着回采、掘进等作业，故需要对主动震源
信号和被动(工程扰动)信号进行辨识；然后对被动震源信号震源定位并对主动震源信号和
被动震源信号分别进行震动波速度成像；最后综合分析两类震源信号波速成像结果以评估
监测区域岩体稳定性状态，圈定采空区域指导施工作业。下面将结合实施例对本发明做进
一步的说明，本发明实施例中所述方法包括如下步骤：

[0049] S1：在监测区域附件布置微震传感器以及在监测区域内布置主动震源，并利用微震传感器实时获取震源信号，并记录每个震源信号触发微震传感器的时刻以及记录主动震
源位置和每个主动震源信号的释放时刻。

[0050] 其中，微震传感器网络需要覆盖监测区域，使得微震传感器既可以接收主动震源信号也可以接收被动震源信号，同时能保证24小时连续采集。为保证被动震源定位的精确，
要使尽可能多的微震传感器接收震源信号(微震传感器数量n，其根据实际需要可以设定)，
且所有微震传感器之间的直线距离优选在200m以内。

[0051] 本实施例中，主动震源于监测区域内，并每隔一段时间自动冲击测点释放地震波，本发明记录主动震源位置为(xa,ya,za)与每次的主动释放时间ta。

[0052] S2：辨识主动震源信号和被动震源信号。

[0053] 本发明实施例中选择训练的好的神经网络机器学习模型进行辨识。其中，神经网络机器学习模型的构建过程如下：

[0054] 首先，利用人工辨识的主动震源信号和被动震源信号建立相应的波形数据库，并分别选取M1组主动震源信号波形和M2组被动震源信号波形作为波形图像样本数据库。优选
M1＝M2；

[0055] 然后，对选取的波形图像样本进行预处理得到每个波形图像对应的二维矩阵。本实施例中，预处理的过程为：将两类波形图像尺寸大小调整至相同，并对每幅波形图像中的
每个像素点的灰度值(范围0‑255)进行提取得到一个二维矩阵，将得到的二维矩阵进行归
一化处理。

[0056] 接着，将每个波形图像对应的二维矩阵以及二维矩阵对应的训练标签输入神经网络机器学习模型中进行训练，其中，依据波形图形为主动震源信号或被动震源信号进行编
码标记并作为训练标签。

[0057] 本实施例中，选用的神经网络机器学习模型为现有模型，二维矩阵(特征图)经过神经网络机器学习模型中卷积核进行卷积操作得到新的特征矩阵，并对其进行下采样操
作，经过数个卷积‑池化操作后降低了特征图的维度，最后将二维特征矩阵转化为一维行向
量。其中，模型架构是现有技术，本发明对此不进行具体的阐述。

[0058] 应当理解，其他可可行的实施例中，也不限于使用上述架构的模型，也不限于使用机器学习模型进行辨识。

[0059] S3：对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻。其中，定位原理如下：

[0060] 假定被动震源位置(x0,y0,z0)和地震波在空区与岩体中传播的速度V均未知，那么被动震源信号经过传播到达第i(i＝1,2,…,n)个微震传感器的到时ti和传播路径Li(i＝1,
2,…,n)之间的关系为：

[0061]

[0062] 其中，t0为震源释放地震波的时刻，ti为触发微震传感器的时刻，(xi,yi,zi)分别为第i个微震传感器位置的空间坐标。

[0063] 同理，将所有接收到被动震源信号的微震传感器对应的传播路径、波速、初始时间和到达时间值之间的关系建立方程组：

[0064]

[0065] 其中，(x1,y1,z1)、(x2,y2,z2)、(xn,yn,zn)分别为第一个、第二个和第n个触发的微震传感器的坐标，t1、t2、tn分别上述触发微震传感器相对应的到时。

[0066] 通过上述建立的方程组可以求解出被动震源的位置坐标，但由于波速在传播过程中的不均匀性以及监测机器在时间同步上的误差导致方程组解出的值为理论值，需要进一
步对方程进行求解以使理论值尽可能逼近真实值。

[0067] 任意两个触发微震传感器之间的到时差可以建立方程组，如下：

[0068]

[0069] 式中，Δt1,2表示微震传感器1和微震传感器2之间的到时差。

[0070] 由于上述公式计算得出的到时差为回归值，其与真实值 之间的偏差为：

[0071]

[0072] 式中，D1,2表示微震传感器1和微震传感器2的到时差与真实值的偏差。

[0073] 回归之与真实值之间偏差越小拟合关系越小，故为了得到接近真实的震源坐标实际是求解使得所有触发传感器之间到时差偏差平方和的最小值，定义到时差偏差平方和F
(x0,y0,z0,V)为：

[0074]

[0075] 对于上述超静定方程，存在一组(x0,y0,z0,V)使得F(x0,y0,z0,V)达到最小，当F(x0,y0,z0,V)达到最小时，满足：

[0076]

[0077]

[0078]

[0079]

[0080] 通过求解上式，可以定位被动震源的位置(x0,y0,z0)以及波速V，再利用上述公式可以计算出被动震源释放地震波的释放时刻t0。上述参数为被动震源反演波速场提供数据
支持。

[0081] 应当理解，利用微震传感器接收震源信号的时间可以鉴别出对应哪个震源的信号。尤其是主动震源可以控制其触发时间。

[0082] S4：反演出主动震源信号的监测区域波速场、被动震源信号的监测区域波速场(震动波速度场成像)。

[0083] 首先，对监测区域按照一定比例划分网格(网格数量为X×Y×Z＝K，网格为立体网格，XYZ分别表示X坐标轴上的坐标、Y坐标轴上的坐标、Z坐标轴上的坐标，K为网格总数)，并
采用A*算法(路径搜索算法)分别追踪主动震源信号、被动震源信号与微震传感器之间的实
际传播路径，记微震传感器i和震源信号j之间的最短传播路径为Lij,那么在xyz的网格内该
路径传播的距离为 波速为 那么第i个微震传感器接收到的第j个震源释放震动
波的走时Tij为：

[0084]

[0085] 假设共有n个微震传感器、N1个主动震源信号和N2个被动震源信号(N1+N2＝N)，第i个微震传感器接收到的第j个震源释放震动波走时的实际值与计算值之间的误差τij为：

[0086] τij＝Tij‑T′ij

[0087] 其中，T′ij为走时的实际值，第N1主动震源中该值为 为第N1主动震源信号的释放时刻， 为第N1主动震源信号触发第i个微震传感器的时刻。同理，被动震
源的走时的实际值也是可以计算出的。监测范围内的N1M+N2M传播路径的走时误差可以表示
为：

[0088]

[0089] 式中， 分别表示监测到的所有主动震源信号、所有被动震源信号在实际传播路径上的走时误差。其中，走时误差 和 越小说明求解得到的每个网格的
波速与现实的出入越小，因此最小的 和 对应的每个网格的波速分别为主动震源
信号和被动震源信号求解得到的波速场。

[0090] S5：基于所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场进行比对实现采空区判定。

[0091] 地震波的传播路径偏向于高速区域介质，故而可通过对监测区域波速场的成像判别采空区等危险区域，依据上述主动震源信号与被动震源信号的成像结果，本实施例中，对
采空区的判定与灾害防控措施如下：

[0092] 施工区域附近无低波速场出现时，可安全施工；

[0093] 当施工区域附近存在两种波形场反演结果相斥的情况，需将主动震源移动至监测区域附近，使尽地震波的传播方向尽可能通过待测区域，反演结果若与被动震源反演结果
相同，可认定圈定范围内存在采空区，如在其附近进行开采作业则需进行支护措施，若反演
结果依旧与之相悖，则以主动震源的反演结果为准，条件允许的情况下联系井下工作人员
对该区域进行人工核定，根据人工探测结果决定是否继续开采；

[0094] 若两种反演结果皆显示该区域附近确实有采空区的存在，则需在施工区域和危险区域连接处建立临时支护措施，根据后续施工与成像结果考虑是否转为永久支护。

[0095] 基于上述方法，本发明提供一种检测系统，包括相互连接通信的震源信号获取模块，辨识模块，被动震源信号定位模块，波速场反演模块以及采空区判定模块。

[0096] 震源信号获取模块：用于获取主动震源信号和被动震源信号以及每个震源信号触发微震传感器的时刻以及主动震源位置和每个主动震源信号的释放时刻；

[0097] 辨识模块：用于辨识每个震源信号为主动震源信号或是被动震源信号；

[0098] 被动震源信号定位模块：用于对被动震源信号进行定位得到被动震源位置、被动震源信号的释放时刻；

[0099] 波速场反演模块：用于反演出主动震源信号的监测区域波速场、被动震源信号的监测区域波速场；

[0100] 采空区判定模块：用于基于所述主动震源信号的监测区域波速场与所述被动震源信号的监测区域波速场进行比对实现采空区判定。

[0101] 具体各个模块的实现过程请参照方法部分的描述，本实施例对此不进行具体的说明。

[0102] 应当理解，上述单元模块的具体实现过程参照方法内容，本发明在此不进行具体的赘述，且上述功能模块单元的划分仅仅是一种逻辑功能的划分，实际实现时可以有另外
的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以
忽略，或不执行。同时，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能
单元的形式实现。

[0103] 在其他实施例中，本发明还提供一种检测终端，包括处理器和存储器，所述存储器存储了计算机程序，所述处理器调用所述计算机程序执行上述检测方法的步骤。

[0104] 在其他实施例中，本发明还提供一种可读存储介质，存储了计算机程序，所述计算机程序被处理器调用执行上述检测方法的步骤。

[0105] 应当理解，在本发明实施例中，所称处理器可以是中央处理单元(Central Processing Unit，CPU)，该处理器还可以是其他通用处理器、数字信号处理器(Digital 
Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，
ASIC)、现成可编程门阵列(Field‑Programmable GateArray，FPGA)或者其他可编程逻辑器
件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理
器也可以是任何常规的处理器等。存储器可以包括只读存储器和随机存取存储器，并向处
理器提供指令和数据。存储器的一部分还可以包括非易失性随机存取存储器。例如，存储器
还可以存储设备类型的信息。

[0106] 所述可读存储介质为计算机可读存储介质，其可以是前述任一实施例所述的控制器的内部存储单元，例如控制器的硬盘或内存。所述可读存储介质也可以是所述控制器的
外部存储设备，例如所述控制器上配备的插接式硬盘，智能存储卡(Smart Media Card,
SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)等。进一步地，所述可读存储
介质还可以既包括所述控制器的内部存储单元也包括外部存储设备。所述可读存储介质用
于存储所述计算机程序以及所述控制器所需的其他程序和数据。所述可读存储介质还可以
用于暂时地存储已经输出或者将要输出的数据。

[0107] 基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分，或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在
一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或
者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的可读存储介
质包括：U盘、移动硬盘、只读存储器(ROM，Read‑OnlyMemory)、随机存取存储器(RAM，
RandomAccess Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

[0108] 需要强调的是，本发明所述的实例是说明性的，而不是限定性的，因此本发明不限于具体实施方式中所述的实例，凡是由本领域技术人员根据本发明的技术方案得出的其他
实施方式，不脱离本发明宗旨和范围的，不论是修改还是替换，同样属于本发明的保护范
围。