基于时序贡献图分解的微震信号自适应解算方法技术

技术编号:39655988 阅读:9 留言:0更新日期:2023-12-09 11:24
本申请提出基于时序贡献图分解的微震信号自适应解算方法,所述方法包括:当监测区域待解算时刻的微震状态为波动状态时,利用时序贡献图算法进行自适应解算得到待解算时刻各通道对应的波动状态总贡献值;根据所述待解算时刻各通道对应的波动状态总贡献值确定待解算时刻各通道对应的波动状态贡献率;根据所述待解算时刻各通道对应的波动状态贡献率确定监测区域待解算时刻的微震根源信号通道

【技术实现步骤摘要】
基于时序贡献图分解的微震信号自适应解算方法


[0001]本申请涉及微震监测
,尤其涉及基于时序贡献图分解的微震信号自适应解算方法


技术介绍

[0002]随着煤矿采深逐年增加,各类监测系统已成为保障煤炭安全生产的重要组成部分,其中微震监测系统作为一种区域

实时

连续的监测系统,可对各类震动信号进行实时采集分析,广泛应用在煤矿安全监测领域

但是,现有技术中对微震信号的监测研究主要集中在两部分:(1)从煤岩体的力学特性入手,借助统计参数和综合处理参数,形成统计参量来分析微震能量变化和趋势

众多学者依据微震特性进行微震事件识别的尝试,形成了如震源参数

震动波波速层析成像技术等理论

(2)通过拟合预测和识别预测的方法对微震数据进行判识,利用多元线性回归

相关系数

机器学习等相关统计学方法建立微震数据预测

判识模型

这两方面的研究主要集中在对煤岩力学特征和对微震单通道信号的研究上,这些方法利用微震监测数据进行微震状态识别做了许多创新性研究,但煤矿系统作为一种复杂的动态系统,其微震信号包含大量噪声且煤矿环境复杂

多变,现有的微震信号识别技术精度较低,进而不能保证煤矿工作人员的安全性


技术实现思路

[0003]本申请提供基于时序贡献图分解的微震信号自适应解算方法,以至少解决现有的微震信号识别技术精度较低的技术问题

[0004]本申请第一方面实施例提出一种基于时序贡献图分解的微震信号自适应解算方法,所述方法包括:当监测区域待解算时刻的微震状态为波动状态时,利用时序贡献图算法进行自适应解算得到待解算时刻各通道对应的波动状态总贡献值;根据所述待解算时刻各通道对应的波动状态总贡献值确定待解算时刻各通道对应的波动状态贡献率;根据所述待解算时刻各通道对应的波动状态贡献率确定监测区域待解算时刻的微震根源信号通道;其中,当所述待解算时刻微震信号数据对应的监测统计量值大于预先确定的监测区域的监测阈值时,判定所述监测区域待解算时刻的微震状态为波动状态

[0005]优选的,所述监测区域的监测阈值的确定过程包括:获取监测区域历史时段内各时刻各通道的微震信号数据;根据所述历史时段内各时刻各通道的微震信号数据和奇异值分解法确定所述监测区域历史时段对应的第一映射矩阵

第二映射矩阵;根据所述监测区域历史时段对应的第一映射矩阵

第二映射矩阵确定所述监测区域历史时段的子状态空间;
根据所述历史时段的子状态空间和
Hotelling

T2检验定理确定所述监测区域的监测阈值;其中,所述监测区域历史时段的子状态空间由第一规范矩阵和第二规范矩阵构成

[0006]进一步的,所述第一规范矩阵的计算式如下:
[0007]式中,
C
为第一规范矩阵,
J
为第五映射矩阵,
P
为第一映射矩阵;所述第二规范矩阵的计算式如下:
[0008]式中,
D
为第二规范矩阵,
L
为第六映射矩阵,为第二映射矩阵

[0009]进一步的,所述待解算时刻微震信号数据对应的监测统计量值的确定过程包括:获取监测区域待解算时刻各通道的微震信号数据;根据所述待解算时刻各通道的微震信号数据和奇异值分解法确定所述待解算时刻微震信号数据对应的第三映射矩阵

第四映射矩阵;根据所述第三映射矩阵

所述第四映射矩阵确定所述待解算时刻微震信号数据对应的子状态空间;根据所述待解算时刻微震信号数据对应的子状态空间确定所述待解算时刻微震信号数据对应的监测统计量值;其中,所述待解算时刻微震信号数据对应的子状态空间由第三规范矩阵和第四规范矩阵构成

[0010]进一步的,所述待解算时刻微震信号数据对应的监测统计量值的计算式如下:
[0011]式中,为在待解算时刻
t
下微震信号数据对应的监测统计量值,为第三规范矩阵
,
为待解算时刻
t
至时刻
t+r
的微震信号数据构成的矩阵,为第五映射矩阵的前
k
列构成的矩阵

[0012]进一步的,所述利用时序贡献图算法进行自适应解算得到待解算时刻各通道对应的波动状态总贡献值,包括:根据所述监测区域待解算时刻各通道的微震信号数据和所述第三映射矩阵分别确定所述监测区域内各通道在待解算时刻至待解算时刻加时各个时刻的波动状态贡献值,其中,为预设的时刻长度;
[0013]根据所述波动状态贡献值确定待解算时刻各通道对应的波动状态总贡献值

[0014]进一步的,所述待解算时刻各通道对应的波动状态总贡献值的计算式如下:
[0015]式中,为在待解算时刻
t
下第
i
个通道对应的波动状态总贡献值,为待解算时刻
j

i
通道对应的波动状态总贡献值,其中,
r
为对待解算时刻
t
解算有影响的时刻总数

[0016]进一步的,所述待解算时刻各通道对应的波动状态贡献率计算式如下:
[0017]式中,为在待解算时刻
t
下第
i
个通道对应的波动状态贡献率,为在待解算时刻
t
下微震信号数据对应的监测统计量值

[0018]进一步的,所述根据所述待解算时刻各通道对应的波动状态贡献率确定监测区域待解算时刻的微震根源信号通道,包括:将所述待解算时刻各通道对应的波动状态贡献率按照从大到小的顺序进行排列,构成贡献率序列;将所述贡献率序列中最大值对应的通道作为所述监测区域待解算时刻的微震根源信号通道

[0019]进一步的,所述方法还包括:将待解算时刻各通道对应的波动状态贡献率进行可视化展示

[0020]本申请的实施例提供的技术方案至少带来以下有益效果:本申请提出了一种基于时序贡献图分解的微震信号自适应解算方法,所述方法包括:当监测区域待解算时刻的微震状态为波动状态时,利用时序贡献图算法进行自适应解算得到待解算时刻各通道对应的波动状态总贡献值;根据所述待解算时刻各通道对应的波动状态总贡献值确定待解算时刻各通道对应的波动状态贡献率;根据所述待解算时刻各通道对应的波动状态贡献率确定监测区域待解算时刻的微震根源信号通道;其中,当所述待解算时刻微震信号数据对应的监测统计量值大于预先确定的监测区域的监测阈值时,判定所述监测区域待解算时刻的微震状态为波动状态

本申请提出的技术方案,可以监测各通道的波动程度,丰富了监测信本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.
一种基于时序贡献图分解的微震信号自适应解算方法,其特征在于,所述方法包括:当监测区域待解算时刻的微震状态为波动状态时,利用时序贡献图算法进行自适应解算得到待解算时刻各通道对应的波动状态总贡献值;根据所述待解算时刻各通道对应的波动状态总贡献值确定待解算时刻各通道对应的波动状态贡献率;根据所述待解算时刻各通道对应的波动状态贡献率确定监测区域待解算时刻的微震根源信号通道;其中,当所述待解算时刻微震信号数据对应的监测统计量值大于预先确定的监测区域的监测阈值时,判定所述监测区域待解算时刻的微震状态为波动状态
。2.
如权利要求1所述的方法,其特征在于,所述监测区域的监测阈值的确定过程包括:获取监测区域历史时段内各时刻各通道的微震信号数据;根据所述历史时段内各时刻各通道的微震信号数据和奇异值分解法确定所述监测区域历史时段对应的第一映射矩阵

第二映射矩阵;根据所述监测区域历史时段对应的第一映射矩阵

第二映射矩阵确定所述监测区域历史时段的子状态空间;根据所述历史时段的子状态空间和
Hotelling

T2检验定理确定所述监测区域的监测阈值;其中,所述监测区域历史时段的子状态空间由第一规范矩阵和第二规范矩阵构成
。3.
如权利要求2所述的方法,其特征在于,所述第一规范矩阵的计算式如下:式中,
C
为第一规范矩阵,
J
为第五映射矩阵,
P
为第一映射矩阵;所述第二规范矩阵的计算式如下:式中,
D
为第二规范矩阵,
L
为第六映射矩阵,为第二映射矩阵
。4.
如权利要求2所述的方法,其特征在于,所述待解算时刻微震信号数据对应的监测统计量值的确定过程包括:获取监测区域待解算时刻各通道的微震信号数据;根据所述待解算时刻各通道的微震信号数据和奇异值分解法确定所述待解算时刻微震信号数据对应的第三映射矩阵

第四映射矩阵;根据所述第三映射矩阵

所述第四映射矩阵确定所述待解算时刻微震信号数据对应的子状态空间;根据所述待解算时刻微震信号数据对应的子状态空间确定所述待解算时刻微震信号数据对应的监测统计量值;其中,所述待解算时刻微震信号数据对应的子状态空间由...

【专利技术属性】
技术研发人员:程健石林松骆意杨凌凯周天白张晓雨
申请(专利权)人:天地科技股份有限公司北京技术研究分公司
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1