一种基于相关性分析的矿山微震信号初至波时刻提取方法技术

一种基于相关性分析的矿山微震信号初至波时刻提取方法，属于矿山微震信号的分析处理方法。其功能为智能识别两信号时间差、实现自动移位对齐和判断矿山微震信号相关性，该方法主要分为五个步骤：一是读取矿山微震信号；二是求解两矿山微震信号的相关函数；三是根据求解出的相关函数，求出最大相关时的时间差；四是根据计算出的时间差对其中一个信号进行移位；五是对对齐后的两矿山微震信号求解其相关系数，判断其相关性。本发明专利技术较好的解决了智能识别矿山微震信号到达时间差及判断矿山微震信号是否来自同一震源的问题。

【技术实现步骤摘要】
一种基于相关性分析的矿山微震信号初至波时刻提取方法
本专利技术涉及矿山微震信号的分析处理方法，特别涉及于一种基于相关性分析的矿山微震信号初至波时刻提取方法。
技术介绍
在自然地震、工程爆破、爆炸等震动事件中，其中有一部分能量必然转化为震动波的形式，它会以震源为中心向周围传播。矿山微震信号初至波拾取是震动研究中非常关键和重要的问题。在震源定位中，准确快捷的拾取到初至波的时刻是进行准确震源定位或信号分析的基础。对采集到的矿山微震信号进行分析和处理，并以此来确定矿山微震信号初至波到达时刻的技术称为初至波拾取技术。初至波拾取技术在军事、民用和工业工程领域都具有非常广泛的应用研究。而矿震是采矿活动诱发的矿井岩体突然失稳破坏的动力现象，它严重威胁着矿井生产与矿工的生命安全，对其进行实时监测和预警具有重要的理论和现实意义。矿震信号相对于地震信号而言，具有震级小、震源浅和影响范围有限等特点，所以可以称之为矿山微震信号。由于可以通过分析矿山微震信号的分布特点及发生频率大致标注出危险地带，实现有效的监测与预警，且因矿山微震信号本身能量释放较小，震波传输距离有限，从而要求矿山微震信号的定位精度相比普通地震信号更高。而定位精度的提高又需要初至波时刻更为快速准确的提取。目前一般采用的矿山微震信号初至波自动拾取技术并无专门针对矿山微震信号的特性和要求的。初至波自动拾取目的是要确定其信号中纯噪声信号和有效信号之间的分界时刻，通常都是根据矿山微震信号的振幅、频率和相位的变化来确定这一时刻。传统的初至拾取方法主要分为两大类：一类是基于地震记录瞬时特征的方法，如极值法(峰值检测)、差分法；这类方法对噪声比较敏感，当地震记录的噪声较严重时，难以准确拾取初至。另一类方法是基于地震记录整体特征的方法，如相关法；这类方法虽然对噪声有较好的抑制作用，但受到地震道之间相似性等因素的影响，对于复杂地震记录，初至拾取的精度也会受到影响。到目前为止，已经提出了许多初至拾取的方法，如人工拾取法、相关法、能量比法、最大振幅法、分形维法及神经网络法等。人工拾取法简单易行，但受人为因素和主观因素影响较大，容易引入人为误差，会直接导致结果误差增大。Gelchinsky和Shtivelman提出了一种用相邻道进行互相关的方法，它假设各道的脉冲形状不发生变化。相关法受初至的续至影响较大，并且与子波的选择关系十分密切，同时要求选择合适的时窗范围，这些在实际工程中都有一定的难度。Hatherrly提出了线性最小平方预测技术与拐点校正相结合的方法，他提出首先识别第一个峰值和拐点，然后估算二者的统计差值。黄成之等采用统计方法将地震初至波记录分成信号和噪音两个部分，并使这两部分统计特征之间的差别为最大。统计特征法收到地震波形相似性影响较大，会对精度造成一定影响。能量比值法使用周期内的信号能量与总时窗能量的比值，对初至比较敏感，续至波衰减比较快，所以将比值的最大值点作为初至的近似值并作适当的时移，即为初至时刻。Coppens提出了在不同大小的时窗内进行能量比较的方法。江玉乐等提出同极性能量比值法，即改进的能量比法。由于能量比值法的抗干扰能力还不够好，所以对于初至波形发生明显变化的地区拾取的初至时刻不够准确。时间域分形维方法拾取初至的过程必须插值，且结果强烈的依靠插值的准确性。FabioBoschetti等提出了一种基于分形维的初至检测算法，该方法是基于地震道随着信号出现其分形维值发生变化的特征来确定地震道初至。但其对时窗和步长的选取十分敏感，稍有不慎就会严重的影响其结果。神经网络法利用多参数特征进行模式识别，充分利用地震记录的瞬时特征和整体特征。由于神经网络不仅具有并行处理、自组织自学习能力，而且具有高度鲁棒性、容错性和高度的映射、计算和分类能力。庄东海等采用将地震记录初至拾取看作一个模式识别过程，充分利用地震记录的瞬时特征和整体特征，用人工神经网络方法进行地震记录初至拾取，其能获得较好的实际效果。但其主要缺点是算法复杂性太高，搜索需要一定时间。地震信号通过小波多分辨分解，可以有效地分离、消除噪声，有利于分形维与神经网络法提高拾取初至的精度。罗光提出了改进型基于小波变换的初至波拾取方法，杨俊峰提出了基于小波变换的三分向震相识别法、能量因子法。但有时其仍需要人工选取数据段来进行拾取以减少拾取时间。除此之外，有一些方法还依赖该道与其近道之间的对比，虽然这一类方法对噪声有一定的压制作用，但受到地震道之间相似性等因素的影响，对于复杂地震记录，初至拾取的精度也会受到影响。综上所述，地震波信号是一个一维的时间序列，由于该序列仅仅表明了时间和幅度之间的关系，且有噪声信号的干扰，对直接进行地震波初至的拾取造成了相当大的困难。解决此类问题可采用如上文所述的智能算法，但由于智能算法普遍具有较高的算法复杂度，不适合用于需要快速、并有较高定位精度的矿山微震环境，而本专利技术较好的解决了对初至波时间的拾取这一问题，并具有较低的算法复杂度。
技术实现思路
本专利技术的目的是针对已有技术中存在的问题，提供一种基于相关性分析的矿山微震信号初至波时刻提取方法，解决矿山微震信号处理中初至波时间智能识别、矿山微震信号对齐和判断矿山微震信号是否来自同一震源的问题。实现本专利技术目的的技术方案：本专利技术的矿山微震信号初至波时刻提取方法包括五个步骤：一是读取矿山微震信号；二是求解两矿山微震信号的相关函数；三是根据求解出的相关函数，求出最大相关时的时间差；四是根据计算出的时间差对其中一个信号进行移位；五是对对齐后的两矿山微震信号求解其相关系数，判断其相关性；具体方法步骤如下：(一)读取矿山微震信号：将两矿山微震信号读取进系统，定义xn为一矿山微震信号序列，yn为另一矿山微震信号序列，要求两序列必须等长，设定其序列长度为L，且两信号采样频率均为fs；(二)求解两矿山微震信号的相关函数：根据下式求解矿山微震信号相关函数，在相关函数取得最大值时，根据采样点及采样率计算出矿山微震信号间的时间差；根据离散信号序列的相关函数公式计算相关函数Rxy为：其中变量m取值范围为0至L；(三)根据求解出的相关函数，求出最大相关时的时间差：根据计算出的时间差对矿山微震信号进行对齐；根据相关函数、采样点与采样率计算最大相关时的时间差，取得相关函数最大值Rmax为：Rmax＝max[|Rxy(m)|](四)根据计算出的时间差对其中一个信号进行移位：对对齐后的矿山微震信号进行相关性分析，并设定阈值以判断矿山微震信号是否来自于同一震源；取得初始时刻至相关函数绝对值最大时刻的采样点总数，记作N，则时间偏差Toffset：Toffset＝N/fs(五)对对齐后的两矿山微震信号求解其相关系数，判断其相关性：两矿山微震信号的相关系数矩阵Rcorr为：其中ρij＝E((Xi-E(Xi))·(Yj-E(Yj)))，其中变量i、变量j取值范围均为1至2，且X为xn序列的总体，Y为yn序列的总体，E为数学期望；求出的相关系数矩阵主对角线表示自相关性，其始终为1；副对角线是两信号本文档来自技高网...

【技术保护点】
一种基于相关性分析的矿山微震信号初至波时刻提取方法，其特征是：矿山微震信号初至波时刻提取方法包括五个步骤：一是读取矿山微震信号；二是求解两矿山微震信号的相关函数；三是根据求解出的相关函数，求出最大相关时的时间差；四是根据计算出的时间差对其中一个信号进行移位；五是对对齐后的两矿山微震信号求解其相关系数，判断其相关性；具体方法步骤如下：(一)读取矿山微震信号：将两矿山微震信号读取进系统，定义xn为一矿山微震信号序列，yn为另一矿山微震信号序列，要求两序列必须等长，设定其序列长度为L，且两信号采样频率均为fs；(二)求解两矿山微震信号的相关函数：根据下式求解矿山微震信号相关函数，在相关函数取得最大值时，根据采样点及采样率计算出矿山微震信号间的时间差；根据离散信号序列的相关函数公式计算相关函数Rxy为：Rxy(m)=Σn=-∞∞x(n)y(n+m)]]>其中变量m取值范围为0至L；(三)根据求解出的相关函数，求出最大相关时的时间差：根据计算出的时间差对矿山微震信号进行对齐；根据相关函数、采样点与采样率计算最大相关时的时间差，取得相关函数最大值Rmax为：Rmax＝max[|Rxy(m)|](四)根据计算出的时间差对其中一个信号进行移位：对对齐后的矿山微震信号进行相关性分析，并设定阈值以判断矿山微震信号是否来自于同一震源；取得初始时刻至相关函数绝对值最大时刻的采样点总数，记作N，则时间偏差Toffset：Toffset＝N/fs(五)对对齐后的两矿山微震信号求解其相关系数，判断其相关性：两矿山微震信号的相关系数矩阵Rcorr为：Rcorr=ρ11ρ12ρ21ρ22]]>其中ρij＝E((Xi‑E(Xi))·(Yj‑E(Yj)))，其中变量i、变量j取值范围均为1至2，且X为xn序列的总体，Y为yn序列的总体，E为数学期望；求出的相关系数矩阵主对角线表示自相关性，其始终为1；副对角线是两信号的相关性，绝对值越大，相关性越好，其取值范围为[0，1]；分界值可根据实际情况不同设定为0到1之间的任何值；最终判断程序是否需要结束，若需结束则结束方法，若需继续运行则返回读取数据步骤，读取新数据进行移位对齐和相关性比较。...

【技术特征摘要】
1.一种基于相关性分析的矿山微震信号初至波时刻提取方法，其特征是：矿山微震信号初至波时刻提取方法包括五个步骤：一是读取矿山微震信号；二是求解两矿山微震信号的相关函数；三是根据求解出的相关函数，求出最大相关时的时间差；四是根据计算出的时间差对其中一个信号进行移位；五是对对齐后的两矿山微震信号求解其相关系数，判断其相关性；具体方法步骤如下：(一)读取矿山微震信号：将两矿山微震信号读取进系统，定义xn为一矿山微震信号序列，yn为另一矿山微震信号序列，要求两序列必须等长，设定其序列长度为L，且两信号采样频率均为fs；(二)求解两矿山微震信号的相关函数：根据下式求解矿山微震信号相关函数，在相关函数取得最大值时，根据采样点及采样率计算出矿山微震信号间的时间差；根据离散信号序列的相关函数公式计算相关函数Rxy为：其中变量m取值范围为0至L；(三)根据求解出的相关函数，求出最大相关时的时间差：根据计算出的时间差对矿山微震信号进行对齐；根据相关函数、采样点与采样率计...

【专利技术属性】
技术研发人员：张申，张然，程婷婷，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：江苏;32