基于多级多尺度网格相似性系数计算的微震震源定位方法
【技术领域】
[0001] 本发明是一种基于多级多尺度网格划分的利用相似性系数的微地震定位反演方 法，涉及地球物理探测技术领域，特别涉及微地震监测定位反演方法。
【背景技术】
[0002] 地下的微小震动作为地下介质大型动力灾害的前兆，对于微震震源的实时定位， 可以有效的监测地下介质的破裂情况以及对于大型动力灾害的预防提供信息。在微地震监 测工作中，如何准确的对微震震源进行定位是监测工作成功的关键。微震监测因其可以实 时监控介质破裂程度的特点，已大量运用在工程领域中，在煤矿系统利用微震监测判断煤 层顶底板破裂情况，在工程物探中利用微震来判断粧基介质的内部情况。
[0003] 目前成熟的微震震源定位算法大多是基于盖格理论的方法，通过检波器空间位 置、地下介质体速度和地震波到达检波器的时间的关系建立方程组，基于最优化理论求得 定位结果。基于盖格理论的定位过程一种是利用非线性反演算法中目标函数的几何性质寻 求迭代方向寻求最优值。但是这种手段往往使最终结果寻求到局部的最优值导致定位结果 不精确。另一种则是利用在搜索空间的随机搜寻来寻找全局的最优值来减低反演算法陷入 局部最优的危险。
[0004] 基于盖格类的定位方法在地震数据信噪比高的情况下可以进行准确的定位，但是 微地震监测中得到的数据信噪比较低，利用盖格方法对数据进行定位处理往往得不到理想 的定位结果。另外，面对大量地震数据，盖格类方法只是利用了地震波的旅行时间，丢弃了 地震数据中的波形数据。
[0005] 基于相似性系数计算的定位方法是一种同时利用地震波旅行时和地震波形数据 计算相似性系数进行定位的方法，在低信噪比或速度模型较复杂时该方法的定位精度较之 盖格类方法有所提高，在天然地震定位中得到了很好的应用。微震震源定位借鉴天然地震 定位的方法，对基于相似性系数计算的定位方法开展了研究并取得了一定的应用成果。然 而微震震源定位在指导资源开采时，对方法的实时性和定位精度都要求较高。为了同时达 到较好的定位精度和实时性，本发明提出一种基于多级多尺度网格划分的相似性系数计算 方法。首先通过对定位空间进行稀疏且均匀的网格划分得到对应的网格点，计算网格点中 相似性系数最高的网格点;然后对相似性系数计算结果最大值的网格点附近的空间进行密 集的网格划分，再一次对网格点进行相似性系数的计算，直至定位精度达到要求，最终得到 的相似性系数值最高的空间点就可认定是最终的微震震源位置。
[0006] 本算法的目的在于多级多尺度的网格划分同时提高了定位方法的计算效率和计 算精度，利用相似性系数定位的方法通过利用地震波形数据进行精确的微震震源定位。此 外，本方法不用拾取地震波从震源到检波器的旅行时间，避免了盖格法在地震信号信噪比 不高的情况下，拾取地震波旅行时不准确造成定位不理想。
【发明内容】
[0007] 本发明是为了解决传统基于相似性系数计算的震源定位算法计算效率和计算精 度无法同时提高，应用在微震震源定位时精度和实时性无法同时满足，导致应用范围受限 的问题，通过多级多尺度网格划分，较之传统均匀网格划分时相似性系数计算的定位方法， 同时提高了定位方法的计算效率和计算精度，满足微震监测实时性的同时，实现震源的精 确定位。
[0008] 基于多级网格划分和相似性系数的微震震源定位方法的基本实现步骤如下：
[0009] 步骤一:将微震监测区域均匀离散化为i行Xj列的网格模型，从原点处逐行进行 编号，将对应的网格分割点坐标位置存储在i行X j列的二维数组Lo中，同时建立i行X j列 的二维数组U，数组U中元素的值为数组Lo中对应网格分割点的平均传播速度，1^的编号顺 序同Lo;
[0010] 步骤二:根据数组LdPU中的值，计算i行Xj列的网格模型中每一个网格位置到每 一个检波器的地震波旅行时，并存储在i行X j列Xm道的三维数组L3中，其中m道表示检波 器道数；
[0011] 步骤三:选取长度为N的时窗，根据数组L3中对应的地震波旅行时对时窗进行移 动，选取视窗内的地震振幅参数；
[0012] 步骤四：根据时窗内选取的振幅信息，将相同行数和列数处对应的所有道的振幅 进行相加，计算相似性系数，并将结果存在存储在i行Xj列的二维数组1^中的对应位置；
[0013] 步骤五:选取步骤四二维数组L2中相似性系数最大的点，以这个点对应的Lo空间坐 标位置作为暂时的定位结果；
[0014] 步骤六:对步骤五中得到的相似性系数最大的网格进行密集地划分，将其均匀离 散化为i行X j列的网格模型，将每个网格分割点对应的空间坐标和空间平均传播速度分别 存入步骤一中提到的二维数组Lo和U中，重复步骤二至对步骤五，直至定位精度满足要求；
[0015] 步骤七:将步骤六中计算得到的相似性系数最大的网格点，对应到Lo数组对应的 真实物理空间位置处，并将该位置是最终的微震震源定位结果。
[0016] 进一步的，在所述步骤一中，根据监测区域的大小将其等效成长方形剖面，在对监 测区域的大小和网格的大小进行设定后，以坐标轴原点处的网格为起始，编号为1，以测线 方向为行，垂直测线方向为列，逐行依次编号，从而将剖面划分为i行和j列，并将对应行数 和列数处网格分割点对应的空间坐标存储在数组L〇中，网格内平均速度的值存入二维数组 U中；若监测空间是三维空间，则可将该区域等效成长方体，二维数组LdPU需要扩充为三 维数组。
[0017] 进一步的，所述步骤二中，每一个网格划分点到每一个检波器的地震波旅行时的 具体表示为：
[0018] 其中/；；Hi表示三维数组U中第i行和第j列所对应的网格点所对应到第m道检波器 的地震波旅行时;Xl，^表示Lo中第i行和第j列对应空间位置的坐标;Xm，ym表示第m道检波器 的坐标。为二维数组1^对应的第i行和第j列网格点的空间物理位置的平均传播速度。
[0019] 进一步的，所述步骤三中，选取长度为N的时窗，根据数组L3中对应的地震波旅行 时对时窗进行移动，选取时窗内的地震振幅参数，其中时窗的长度N至少要大于权利要求1 步骤二中计算得到的数组U中的最大值。
[0020] 进一步的，所述步骤四中，据时窗内选取的振幅信息，将相同行数和列数处对应的 所有道的振幅进行相加，计算相似性系数，并将结果存在存储在i行Xj列的二维数组1^中 的对应位置，其具体形式为：
[0021] 其中是表示二维数组1^中，第i行和第j列对应的网格点的相似性系数，N表示时 窗的长度，Μ表示检波器的个数，表示对第m个地震道的数据，将时窗进行长度 为的移动后，时窗内第η个点的地震数据幅值。
[0022] 进一步的，所述步骤五中，选取步骤四二维数组1^2中相似性系数最大的点，以这个 点对应的Lo空间坐标位置作为暂时的定位结果;若监测区域出现多个微震震源，则需要建 立一个阈值，超过阈值的相似性系数则被认定为暂时的定位结果，对应的空间网格处都需 要再次加密细分。
[0023] 进一步的，所述步骤六中，对步骤五中得到的相似性系数最大的网格进行密集地 划分，将每个网格分割点对应的空间坐标和空间平均传播速度分别存入步骤一中提到的二 维数组Lo和U中，重复步骤二至对步骤五，直至定位精度满足要求。
[0024] 进一步的，所述步骤七中，其特征在于将步骤六中计算得到的相似性系数最大的 网格点，对应到数组Lo中对应的真实物理空间位置处，并将该位置作为最终的微震震源定 位结果。
[0025] 本发明的技术方案，首先通过对监测区域进行稀疏且均匀的网格划分，计算网格 点对应的相似性系数，得到一个距离真实位置比较接近的空间位置，再以此位置为中心，以 一定的范围进行密集的网格划分，再次通过计算相似性系数，选取相似性系数值最大的点 作为最终的定位结果。这样，能够避免监测区域整体的密集网格划分造成的计算时间过长 的缺点，也能够避免由于整体的稀疏网格划分导致的定位精度过低的缺点，最终实现尚效、 高精度的微震震源定位。 【附图说明】 图1是基于多级多尺度网格相似性系数计算的微震震源定位方法的流程框图。 【具体实施方式】 基于多级网格划分和相似性系数的微震震源定位方法的基本实现步骤如下： 步骤一:将微震监测区域均匀离散化为i行Xj列的网格模型，将对应的网格分割点的 坐标位置存储在i行X j列的二维数组Lo中，从原点处逐行进行编号，建立i行X j列的二维 数组U，数组U中元素的值为对