一种基于环境适应算法的时差微震震源定位方法技术

本发明专利技术提供一种基于环境适应算法的时差微震震源定位方法，(1)微震信号被接收；(2)去噪处理；(3)初至到时拾取；(4)构造目标函数，计算拾取到时与计算到时的时差，求解微震震源位置；(5)初始化IEAM

【技术实现步骤摘要】
一种基于环境适应算法的时差微震震源定位方法


[0001]本专利技术提供一种基于环境适应算法的时差微震震源定位方法，属于微震定位


技术介绍

[0002]微震定位指根据检波器接收到的微地震信号，推演出微地震发生位置的一种地球物理方法，在岩土工程中被广泛应用。微震定位的精度和效率与煤矿资源开采的高效性与安全性息息相关。
[0003]微震定位算法分为类偏移定位算法和时差定位算法，前者计算效率低所以国内外往往采用时差定位方法。时差定位方法起源于Geiger于1912年提出的经典定位算法，其基本思想是将非线性方程组线性化并给出最小二乘解，但该算法受初始值影响较大。随着多学科的交流渗透，遗传算法、差分进化算法、微粒群优化算法等智能计算方法得到了迅速发展，基于此，国内外学者做了广泛研究以提高微震定位的精度和稳定性。例如，粒子群差分进化定位算法(DEPSO)、震源扫描算法(SSA)、网格搜索
‑
牛顿迭代法、MOPSO
‑
SA混合算法、基于斯奈尔定律的布谷鸟算法等，这些较为新颖的方法进一步提高了微震定位的精度与稳定性，但同时也有着各自的缺陷。采用DEPSO与单纯用粒子群算法进行微震定位相比，定位精度有着一定提升，但在采集信号信噪比极低时其定位精度仍然会受影响；SSA和网格搜索法有着更长的定位耗时；MOPSO
‑
SA有效避免了寻优过程陷入局部极值的问题，但由于受微震信号检测通道个数影响，定位性能依然需要进一步优化。
[0004]另外还有改进定位目标函数的差分进化微震定位方法(DE)，基于传统微震P波走时定位方法，在定位目标函数中加入了误差加权系数，通过降低初至拾取精度低的微震信号道对定位误差的影响，减弱了初至拾取精度对微震定位的影响；利用差分进化算法，通过种群初始化、变异、交叉和选择等操作，求解改进的微震定位目标函数，实现微震定位。但是该方法在种群规模小的情况下，收敛速度慢且稳定性差。

技术实现思路

[0005]针对上述技术问题，本专利技术提供一种基于环境适应算法的时差微震震源定位方法(IEAM
‑
R)，达到以下两个目的：
[0006](1)基于IEAM
‑
R提出一种所需种群规模更小的方法，减小反演过程的计算量，使整个流程更高效；
[0007](2)加速反演算法的收敛过程，并进一步提高定位结果的准确性与稳定性。
[0008]本专利技术的方案为：
[0009]一种基于环境适应算法的时差微震震源定位方法，包括以下步骤：
[0010](1)在微震震源待测地点的不同位置设置多个检波器，当微震事件发生时，产生的微震信号被检波器接收并传送给计算机；
[0011](2)对监测得到的微震信号进行去噪处理，通过波形相关或其他检测算法提取有
效微震事件，剔除无效微震事件；
[0012](3)利用时窗比法与AIC结合方法进行初至到时拾取；
[0013](4)构造下述目标函数，计算拾取到时与计算到时的时差，利用最优化思想求解微震震源位置：
[0014][0015][0016]d
i
＝|E
i
‑
E
med
|
ꢀꢀ
(3)
[0017][0018]式中，(x0，y0，z0)为震源坐标，t0为发震时刻，N为检波器个数，r
i
为各道误差加权系数，t
ical
表示第i个检波器的模拟波至时间，t
iobs
表示第i个检波器的观测波至时间，E
i
表示第i个检波器的计算误差，E
med
表示第i个检波器的误差中值；
[0019](5)初始化IEAM
‑
R的参数，包括种群规模POPSIZE、最大迭代次数MIT，搜索空间范围[lb，ub]3；
[0020](6)产生初始种群oldPop，在搜索空间中随机生成每个个体位置坐标，并通过公式(1)计算每个个体的适应度值；种群中，每个个体的位置坐标构成一个震源位置的解(x0’
，y0’
，z0’
)，其适应度值用来评估与实际震源位置的接近程度，适应度值越小，则与实际震源位置越接近；
[0021](7)计算当前种群的平均适应度值F
avg
；所有个体按适应度值大小升序排序后，设置序号index来将种群分为适应度优于F
avg
的部分和劣于F
avg
的部分；
[0022](8)从适应度优于F
avg
的个体中随机选择个体I1，从适应度劣于F
avg
的个体中随机选择个体I2；
[0023](9)用I1与I2的差值定义带宽bw；
[0024](10)生成临时种群tempPop，其中每个个体的位置坐标通过以下公式产生；从而使I1中的个体开发更优的潜在解，I2中的个体探索新的搜索空间；
[0025]tempPop(i)＝oldPop(i)+bw*randi＝1，2，...，index
ꢀꢀ
(5)
[0026]tempPop(j)＝(F
j
/F
index
)*oldPop(j)+rand j＝index，index+1，...，POPSIZE
ꢀꢀ
(6)
[0027]其中，tempPop(i)表示临时种群中的第i个体，oldPop(i)表示当前种群中第i个体，rand是(0，1)内的随机数；
[0028](11)对tempPop中超出搜索空间的个体进行越界处理；
[0029](12)通过式(1)计算tempPop中所有个体的适应度值；
[0030](13)将tempPop与oldPop合并，保留适应度更优的POPSIZE个个体作为新的oldPop；
[0031](14)判断是否满足终止条件；满足则输出当前oldPop中最优个体的位置，即反演得到的震源坐标；否则重复步骤(7)
‑
(14)。
[0032]本专利技术首次将环境适应算法IEAM
‑
R应用于微震震源定位问题，取得了比已有的差分进化算法更好的效果。
附图说明
[0033]图1为本专利技术所述方法流程图；
[0034]图2为实施例的IEAM
‑
R与DE的进化曲线；
[0035]图3为实施例的IEAM
‑
R与DE定位结果。
具体实施方式
[0036]由于井下监测的检波器基本在目标层深度附近，且在储层附近地层分层比较明显，因此层状速度模型比较适用于井下监测。下面结合附图和具体实施例进一步说明本专利技术的技术方案。
[0037]如图1所示的流程，包括以下步骤：
[0038]首先，采用地面稀疏观测建立五层的水平层状速度模型，在地面放置埋深1m的检波本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于环境适应算法的时差微震震源定位方法，其特征在于，包括以下步骤：(1)在微震震源待测地点的不同位置设置多个检波器，当微震事件发生时，产生的微震信号被检波器接收并传送给计算机；(2)对监测得到的微震信号进行去噪处理；(3)利用时窗比法与AIC结合方法进行初至到时拾取；(4)构造目标函数，计算拾取到时与计算到时的时差，利用最优化思想求解微震震源位置；(5)初始化IEAM
‑
R的参数；(6)产生初始种群oldPop，在搜索空间中随机生成每个个体位置坐标，并通过步骤(4)构造的目标函数计算每个个体的适应度值；(7)计算当前种群的平均适应度值F
avg
；所有个体按适应度值大小升序排序后，设置序号index来将种群分为适应度优于F
avg
的部分和劣于F
avg
的部分；(8)从适应度优于F
avg
的个体中随机选择个体I1，从适应度劣于F
avg
的个体中随机选择个体I2；(9)用I1与I2的差值定义带宽bw；(10)生成临时种群tempPop，产生每个个体的位置坐标；从而使I1中的个体开发更优的潜在解，I2中的个体探索新的搜索空间；(11)对tempPop中超出搜索空间的个体进行越界处理；(12)通过步骤(4)构造的目标函数计算tempPop中所有个体的适应度值；(13)将tempPop与oldPop合并，保留适应度更优的POPSIZE个个体作为新的oldPop；(14)判断是否满足终止条件；满足则输出当前oldPop中最优个体的位置，即反演得到的震源坐标；否则重复步骤(7)
‑
(14)。2.根据权利要求1所述的一种基于环境适应算法的时差微震震源定位方法，其特征在于，步骤(2)中通过波形相关检测算法提取有效微震事件，提取有效微震事件，剔除无效微震事件。3.根据权利要求1所述的一种基于环境适应算法的时差微震震源定位方法，其特...

【专利技术属性】
技术研发人员：李泽文，张启升，吴祯，乔帅卿，郭丰，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：