一种微震震源自动定位方法技术

本发明专利技术涉及一种微震震源自动定位方法。包括：根据获取的条有效通道的信息构建出个四元一次方程组；为有效通道数，；对所述个四元一次方程组中每一个四元一次方程组求解，得到组解；把所述组解分别带入震源定位目标方程中，得到个对应每条有效通道的偏差值；判断所述个偏差值是否均满足误差范围，如果是，输出所述震源定位目标方程中的偏差目标函数值，否则，改变个偏差值中不满足误差范围的偏差值对应的有效通道信息，然后执行根据所述个有效通道的信息构建个四元一次方程组的步骤。本发明专利技术的定位方法通过计算机来实现，省时省力。

【技术实现步骤摘要】
一种微震震源自动定位方法
本专利技术涉及微震检测领域，尤其涉及一种微震震源自动定位方法。
技术介绍
在煤矿开采活动中有裂缝等岩层破坏事件时，岩体的瞬间破裂会激发弹性波。这些弹性波携带着震源的信息，依赖岩体弹性介质向四周传播。矿山微震监测就是利用震动传感器在远处测量这些弹性波信号，并根据所监测的微震信号特征来确定破裂的发生时间、空间位置、尺度、强度及性质，做出有效合理的预警，具有远程、三维和实时的特点，因此在矿山、岩土、水利水电等工程中具有广阔的应用前景，日益受到人们的重视。而矿山微震检测技术的核心是真元定位，能否进行快速、准确的确定为是微震检测系统能否发挥作用的关键。传统的震源定位计算会将震源定位的目标方程转换成对四元一次方程求解的问题，理论上当给出已知的检波器坐标和标尺位置，就能求出方程组的解，实际上，虽然构造的多个方程组能求出多个解，但是未必有满足目标函数的解，这就要求检波器的坐标和我们对标尺位置的标注尽量准确，这样，当标出位置后，发现某一条通道目标函数较大时，就要重新对该条通道进行标注。对于多条通道，每一条通道都要标出准确位置，既费时又费力。
技术实现思路
为解决上述技术问题，本专利技术的目的是提供一种微震震源自动定位方法，通过计算机来实现定位计算，达到省时省力的目的。本专利技术的一种微震震源自动定位方法，包括：根据获取的条有效通道的信息构建出个四元一次方程组；为有效通道数，；对所述个四元一次方程组中每一个四元一次方程组求解，得到组解；把所述组解分别带入震源定位目标方程中，得到个对应每条有效通道的偏差值；判断所述个偏差值是否均满足误差范围,如果是，输出所述震源定位目标方程中的偏差目标函数值，否则，改变所述个偏差值中不满足误差范围的偏差值对应的有效通道信息，然后执行根据所述个有效通道的信息构建个四元一次方程组的步骤。其中，有效通道信息包括：标尺相对位置、标尺相对时间和检波器坐标。进一步地，得到个对应每条有效通道的偏差值之后，输出所述震源定位目标方程中的偏差目标函数值之前，还包括：根据每一条有效通道偏差值计算个对应每条有效通道的均方值；判断所述个均方值中最小均方值是否小于预设阈值，如果是，则执行判断所述个偏差值是否均满足误差范围的步骤，否则，所述震源定位目标方程无解。具体地，根据获取的条有效通道的信息构建出个四元一次方程组，包括：根据获取的条有效通道的信息构建出个定位计算公式；分别对所述个定位计算公式两边平方，得到个平方后的定位计算公式；将所述个平方后的定位计算公式两两相减，得到个四元一次方程；构建出个四元一次方程组。其中定位计算公式为，为震源坐标，为第个测点的坐标，为震源震动时间，为P波到达第个测点的时间，为P波的传播速度。其中，震源定位目标方程为，其中，为偏差目标函数，为有效通道数，为震源坐标，为震源震动时间，为P波到达第个测点的时间，为第个测点的坐标，为P波的传播速度。具体地，改变个偏差值中不满足误差范围的偏差值对应的有效通道信息，包括：判断个偏差值中不满足误差范围的偏差值是否大于0，如果是，则减小不满足误差范围的偏差值对应的有效通道标尺位置，否则，增大不满足误差范围的偏差值对应的有效通道标尺位置。具体地，对个四元一次方程组中每一个四元一次方程组采用全选主元高斯－约当消去法求解。借由上述方案，本专利技术至少具有以下优点：本专利技术实施例只须动标注有效通道信息一次，然后就可以通过计算机代替手工方法寻找最优解，在某条有效通道信息偏差值不满足误差范围时，计算机也自动修改有效通道信息，重新寻找最优解，即偏差目标函数值也就是震源坐标实际和理论的差值，相对于现有技术的手工寻找最优解省时省力。上述说明仅是本专利技术技术方案的概述，为了能够更清楚了解本专利技术的技术手段，并可依照说明书的内容予以实施，以下以本专利技术的较佳实施例并配合附图详细说明如后。附图说明图1是本专利技术实施例提供的一种微震震源自动定位方法示意图。具体实施方式下面结合附图和实施例，对本专利技术的具体实施方式作进一步详细描述。以下实施例用于说明本专利技术，但不用来限制本专利技术的范围。参见图1，一种微震震源自动定位方法，包括如下步骤：1、根据获取的条有效通道信息构建出个四元一次方程组；为有效通道数，；有效通道信息可以是通过自动标注获取的信息，每条有效通道信息包括：标尺相对位置、标尺相对时间和检波器坐标。标尺是指发生震动事件的那一刻所在的时间段（本专利技术实施例中定位5秒钟）中，标尺标注的位置就是这5秒中微震动起震的那一刻。假如事件发生在2014.1.2.15：20：30，该时间段即为15：20：28~15：20：32这5秒数据，标尺就应该标在微震动起震位置15：20：30这一刻。获取上述条有效通道信息后，构建个定位计算公式：具体地，定位基本公式为：（1）其中，为震源坐标，为第个测点的坐标，为震源震动时间，为P波到达第个测点的时间，为P波均匀的传播速度。有条有效通道，则有个定位计算公式（1）。分别对上述个定位计算公式(1)两边平方,得到个平方后的定位计算公式；定位基本公式为公式（1）时，公式（1）两边平方后则变为(2)有个定位计算公式（1），则对应有个公式（2）。将上述个公式（2）两两相减，得到个四元一次方程的标准形式，于是，可以构建出（表示平面定位）个四元一次方程组。其中，C是求组合个数。本专利技术实施例以为例进行说明，根据4条有效通道信息构建4个定位计算公式（1）,则对应得到4个公式（2）。将4个公式（2）两两相减，得到个四元一次方程，因此，构建出个四元一次方程组。2、对上述构建的个方程组中每一个方程组进行求解，得到组解。具体地,可以采用全选主元高斯－约当消去法、复系数方程组的全选主元高斯消去法、复系数方程组的全选主元高斯－约当消去法、对角线方程组的追赶法、大型稀疏方程组的全选主元高斯－约去消去法、求解托伯利兹方程组的列文逊方法等。本专利技术实施例以采用全选主元高斯－约当消去法进行求解。3、把步骤2中得到的组解分别带入震源定位目标方程中，得到个对应每条有效通道的偏差值；4、判断上述个偏差值是否均满足误差范围,如果是，则执行步骤5；否则，执行步骤6；5、输出震源定位目标方程中的偏差目标函数值，偏差目标函数值也就是震源坐标实际和理论的差值。6、改变个偏差值中不满足误差范围的偏差值对应的有效通道信息，然重新执行步骤1。其中，震源定位目标方程为（3）公式（3）的约束条件为。公式（3）中为偏差目标函数，为有效通道数，为震源坐标，为震源震动时间，为P波到达第个测点的时间，为第个测点的坐标，为P波的传播速度。根据公式（3）可得到每一条有效通道的偏差值，即。进一步地，在步骤3之后，步骤4之前，还包括：根据每一条有效通道偏差值计算个对应每条有效通道的均方值；判断上述个均方值中最小均方值是否小于预设阈值，如果是，则执行步骤4，否则，震源定位目标方程无解。若每一条有效通道的偏差值分别为，则每条有效通道的均方值为。通过判断上述个均方值中最小均方值是否小于预设阈值，可以确定公式（3）是否存在最优解，如果是，则公式（3）存在最优解；如果否，则公式（3）不存在满足条件的解，直接结束算法，这样可节省计算时间。进一步地，步骤4中改变个偏差值中不满足误差范围的偏差值对应的有效通道信息包括：判断个偏差值中不满足误差范围的偏差本文档来自技高网...

【技术保护点】
一种微震震源自动定位方法，其特征在于：根据获取的条有效通道的信息构建出个四元一次方程组；为有效通道数，；对所述个四元一次方程组中每一个四元一次方程组求解，得到组解；把所述组解分别带入震源定位目标方程中，得到个对应每条有效通道的偏差值；判断所述个偏差值是否均满足误差范围,如果是，输出所述震源定位目标方程中的偏差目标函数值，否则，改变所述个偏差值中不满足误差范围的偏差值对应的有效通道信息，然后执行根据所述个有效通道的信息构建个四元一次方程组的步骤。

【技术特征摘要】
1.一种微震震源自动定位方法，其特征在于：根据获取的n条有效通道的信息构建出k个四元一次方程组；n为有效通道数，n≥4；对所述k个四元一次方程组中每一个四元一次方程组求解，得到l组解，且l满足条件l≤k；把所述l组解分别带入震源定位目标方程中，得到n个对应每条有效通道的偏差值，且l满足条件l≤k；判断所述个偏差值是否均满足误差范围,如果是，输出所述震源定位目标方程中的偏差目标函数F(x0,y0,z0)值，否则，改变所述n个偏差值中不满足误差范围的偏差值对应的有效通道信息，然后执行根据所述n个有效通道的信息构建k个四元一次方程组的步骤；得到n个对应每条有效通道的偏差值之后，输出所述震源定位目标方程中的偏差目标函数值之前，还包括：根据每一条有效通道偏差值计算n个对应每条有效通道的均方值;其中，每一条有效通道的所述偏差值分别为F(1),F(2)…F(n)，则每条有效通道的均方值为[F(1)2＋F(2)2＋...F(n)2]/n；判断所述n个均方值中最小均方值是否小于预设阈值，如果是，则执行判断所述n个偏差值是否均满足误差范围的步骤，否则，所述震源定位目标方程无解。2.根据权利要求1所述的微震震源自动定位方法，其特征在于，所述有效通道信息包括：标尺相对位置、标尺相对时间和检波器坐标。3.根据权利要求1所述的微震震源自动定位方法，其特征在于，根据获取的n条有效通道的信息构建出k个四元一次方程组，包括：根据获取的n条...

【专利技术属性】
技术研发人员：石杰，朱尚嵩，蒋建清，王飞，
申请(专利权)人：江苏三恒科技股份有限公司，
类型：发明
国别省市：江苏;32