冻结立井爆破近区井壁振动信号基线漂移校正和消噪方法

付晓强 杨仁树 刘纪峰 张会芝 张仁巍

引用本文:
Citation:

冻结立井爆破近区井壁振动信号基线漂移校正和消噪方法

    作者简介: 付晓强(1984- ),男,博士,讲师,fuxiaoqiang1984@163.com;
  • 中图分类号: O389; TD235.1

Baseline drift correction and de-noising method of shaft lining vibration signal in near field of freezing vertical shaft blasting

  • CLC number: O389; TD235.1

  • 摘要: 冻结立井爆破过程中,近区监测信号中含有的基线漂零及噪声成分对其局部特征精细化提取影响显著。在对近区井壁振动信号有效采集基础上,通过互补总体经验模态分解(complementary ensemble empirical mode decomposition, CEEMD)方法、稀疏化基线估计消噪(baseline estimation and de-noising with sparsity, BEADS)方法和隐马尔可夫模型消噪(hidden Markov model de-noising, HMMD)方法等,解决了信号中基线漂移和随机噪声消除难题,并采用交叉小波变换对校正和消噪效果进行了相关性评价。实例分析结果表明:信号中缓变的基线成分遍历信号各个模态分量的整个过程,且主要集中于低频分量中,而噪声则集中在高频分量。组合分析方法对低频基线漂零和高频噪声的处理效果好,是一种高效且相对保幅的信号分析方法,可用于批量信号数据的预处理过程。
  • 图 1  基线校正及消噪流程

    Figure 1.  Baseline correction and noise reduction process

    图 2  传感器井壁预埋法

    Figure 2.  Pre-embedding method of vibration instrument

    图 3  炮眼布置(单位:mm)

    Figure 3.  Borehole layout (unit: mm)

    图 4  爆破近区井壁振动信号

    Figure 4.  Vibration signal of shaft lining near blasting area

    图 5  各IMF分量筛分迭代次数关系

    Figure 5.  The relationship of iterations number and shift number

    图 6  信号CEEMD分解各分量及残余项R

    Figure 6.  Component of CEEMD decomposition and residual signal

    图 7  模态分量处理结果

    Figure 7.  Results of baseline correction processing

    图 8  各分量基线成分

    Figure 8.  Baseline of each component

    图 9  罚函数值与迭代次数关系

    Figure 9.  The relation of cost function and the iteration number

    图 10  基线校正后信号

    Figure 10.  Baseline corrected signal

    图 11  CEEMD重构信号

    Figure 11.  reconstructed signals of CEEMD

    图 12  消噪后的校正信号

    Figure 12.  Baseline corrected signal after de-noising

    图 13  信号功率谱密度(wf

    Figure 13.  Power spectral density (wf) of signal

    图 14  信号相关性凝聚谱

    Figure 14.  Correlation condensation spectra of signals

    表 1  模态分量与原信号相关度

    Table 1.  Correlation between components and original signals

    IMF相关性系数IMF相关性系数IMF相关性系数
    10.420250.502390.1083
    20.448460.1419100.3707
    30.375470.1333110.4540
    40.412480.1364120.2940
    下载: 导出CSV

    表 2  不同方法消噪性能指标

    Table 2.  Indexes of noise reduction performance

    去噪方法评价指标
    SNRRMSECCPE
    MD3.4420.5810.8161.758
    SVD20.9830.3790.8961.175
    WED21.3970.3740.9210.492
    HMMD43.3710.0560.9920.244
    下载: 导出CSV
  • [1] 杨仁树, 付晓强, 杨立云, 等. 冻结立井爆破冻结壁成形控制与井壁减振研究 [J]. 煤炭学报, 2016, 41(12): 2975–2985. DOI: 10.13225/j.cnki.jccs.2016.1487.
    YANG R S, FU X Q, YANG L Y, et al. Research on the shaping control of frozen wall and blasting vibrationmitigation of shaft wall effect in mine vertical shaft [J]. Journal of China Coal Society, 2016, 41(12): 2975–2985. DOI: 10.13225/j.cnki.jccs.2016.1487.
    [2] 程桦, 林键, 姚直书, 等. 我国西部地区孔隙型含水基岩段立井单层井壁外荷载研究 [J]. 岩石力学与工程学报, 2019, 38(3): 542–550. DOI: 10.13722/j.cnki.jrme.2018.1155.
    CHENG H, LIN J, YAO Z S, et al. Study on the external load of a single layer wall of the vertical shaft in the pore type water-bearing bedrock section in Western China [J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(3): 542–550. DOI: 10.13722/j.cnki.jrme.2018.1155.
    [3] 马芹永, 袁璞, 张经双, 等. 立井直眼微差爆破模型试验振动测试与分析 [J]. 振动与冲击, 2015, 34(6): 172–176. DOI: 10.13465/j.cnki.jvs.2015.06.033.
    MAQ Y, YUAN P, ZHANG J S, et al. Blasting vibration measurement and analyses of millisecondblasting models for vertical shaft blasting [J]. Journal of Vibration and Shock, 2015, 34(6): 172–176. DOI: 10.13465/j.cnki.jvs.2015.06.033.
    [4] 张军, 潘泽鑫, 郑玉新, 等. 振动信号趋势项提取方法研究 [J]. 电子学报, 2017, 45(1): 22–28. DOI: 10.3969/j.issn.0372-2112.2017.01.004.
    ZHANG J, PAN Z X, ZHENG Y X, et al. Research on vibration signal trend extraction [J]. Acta Electronica Sinica, 2017, 45(1): 22–28. DOI: 10.3969/j.issn.0372-2112.2017.01.004.
    [5] 张景元, 何玉珠. 基于形态学的自动驾驶仪振动信号基线漂移去噪 [J]. 北京航空航天大学学报, 2018, 44(5): 907–913. DOI: 10.13700/j.bh.1001-5965.2017.0371.
    ZHANG J Y, HE Y Z. Removing baseline drift in vibration signal of autopilot based on morphology [J]. Journal ofBeijing University of Aeronautics and Astronautics, 2018, 44(5): 907–913. DOI: 10.13700/j.bh.1001-5965.2017.0371.
    [6] 刘艳丽, 赵为松, 李海坤, 等. 基于形态滤波的脉搏波信号基线漂移消除方法研究 [J]. 合肥工业大学学报(自然科学版), 2011, 34(4): 525–528. DOI: 10.3969/j.issn.1003-5060.2011.04.011.
    LIU Y L, ZHAO W S, LI H K, et al. Research on removing baseline wandering of pulse wave signal based on morphologicalfilter [J]. Journal of Hefei University of Technology, 2011, 34(4): 525–528. DOI: 10.3969/j.issn.1003-5060.2011.04.011.
    [7] 谢芳娟, 朱淑云. 基于空域追踪算法的基线漂移信号噪声修正 [J]. 沈阳工业大学学报, 2016, 38(6): 692–696. DOI: 10.7688/j.issn.1000-1646.2016.06.17.
    XIE F J, ZHU S Y. Baseline drift noise correction based on null space pursuit algorithm [J]. Journal of Shenyang University of Technology, 2016, 38(6): 692–696. DOI: 10.7688/j.issn.1000-1646.2016.06.17.
    [8] 张胜, 凌同华, 曹峰, 等. 模式自适应连续小波去除趋势项方法在爆破振动信号分析中的应用 [J]. 爆炸与冲击, 2017, 37(2): 255–261. DOI: 10.11883/1001-1455(2017)02-0255-07.
    ZHANG S, LING T H, CAO F, et al. Application of removal trend method of pattern adapted continuouswavelet to blast vibration signal analysis [J]. Explosion and Shock Waves, 2017, 37(2): 255–261. DOI: 10.11883/1001-1455(2017)02-0255-07.
    [9] 韩亮, 刘殿书, 辛崇伟, 等. 深孔台阶爆破近区振动信号趋势项去除方法 [J]. 爆炸与冲击, 2018, 38(5): 1006–1012. DOI: 10.11883/bzycj-2016-0194.
    HAN L, LIU D S, XIN C W, et al. Trend removing methods of vibration signals of deep hole bench blasting in near field [J]. Explosion and Shock Waves, 2018, 38(5): 1006–1012. DOI: 10.11883/bzycj-2016-0194.
    [10] 王志亮, 陈贵豪, 黄佑鹏. EEMD修正爆破加速度零漂信号中的最优白噪声系数 [J]. 爆炸与冲击, 2019, 39(8): 084201. DOI: 10.11883/bzycj-2019-0154.
    WANG Z L, CHEN G H, HUANG Y P. Optimal white noise coefficient in EEMD corrected zero drift signal of blastingacceleration [J]. Explosion and Shock Waves, 2019, 39(8): 084201. DOI: 10.11883/bzycj-2019-0154.
    [11] YEH J R, SHIEH J S, HUANG N E. Complementary ensemble empirical mode decomposition: a novel noise enhanced data analysis method [J]. Advances in Adaptive Data Analysis, 2010, 2(2): 135–156. DOI: 10.1142/S1793536910000422.
    [12] 乐友喜, 杨涛, 曾贤德. CEEMD与KSVD字典训练相结合的去噪方法 [J]. 石油地球物理勘探, 2019, 54(4): 729–736, 721. DOI: 10.13810/j.cnki.issn.1000-7210.2019.04.001.
    YUEY X, YANG T, ZENG X D. Seismic denoising with CEEMD and KSVD dictionary combined training [J]. OilGeophysicalProspecting, 2019, 54(4): 729–736, 721. DOI: 10.13810/j.cnki.issn.1000-7210.2019.04.001.
    [13] 杨仁树, 付晓强, 杨国梁, 等. 基于CEEMD与TQWT组合方法的爆破振动信号精细化特征提取 [J]. 振动与冲击, 2017, 36(3): 38–45.
    YANG R S, FU X Q, YANG G L, et al. Precise feature extraction of blasting vibration signals based on combined method ofCEEMD and TQWT [J]. Journal of Vibration and Shock, 2017, 36(3): 38–45.
    [14] NING X R, SELESNICK I W, DUVAL L. Chromatogram baseline estimation and denoising using sparsity (BEADS) [J]. Chemometrics and IntelligentLaboratory Systems, 2014, 139: 156–167. DOI: 10.1016/j.chemolab.2014.09.014.
    [15] 张进, 冯志鹏, 卢文秀, 等. 交叉小波变换在水轮机非平稳信号分析中的应用 [J]. 中国电机工程学报, 2010, 30(23): 84–89. DOI: 10.13334/j.0258-8013.pcsee.2010.23.017.
    ZHANG J, FENG Z P, LU W X, et al. Application of cross-wavelet transform to hydraulic turbine nonstationary signal analysis [J]. Proceeding of the CSEE, 2010, 30(23): 84–89. DOI: 10.13334/j.0258-8013.pcsee.2010.23.017.
    [16] 程知, 何枫, 靖旭, 等. 基于集合经验模态分解和奇异值分解的激光雷达信号去噪 [J]. 光子学报, 2017, 46(12): 1201003. DOI: 10.3788/gzxb20174612.1201003.
    CHENG Z, HE F, JING X, et al. Denoising lidar signal based on ensemble empirical mode decomposition and singular value decomposition [J]. Acta Photonica Sinica, 2017, 46(12): 1201003. DOI: 10.3788/gzxb20174612.1201003.
    [17] ROSSO O A, BLANCO S, YORDANOVA J, et al. Wavelet entropy: a new tool for analysis of short duration brain electricalsignals [J]. Journal of Neuroscience Methods, 2001, 105(1): 65–75. DOI: 10.1016/s0165-0270(00)00356-3.
    [18] 王东, 王新晴, 梁升, 等. 基于数学形态学的流量传感器信号去噪研究 [J]. 传感技术学报, 2012, 25(8): 1097–1101. DOI: 10.3969/j.issn.1004-1699.2012.08.016.
    WANG D, WANG X Q, LIANG S, et al. Study on signal denoising of flow sensor based on mathematical morphology [J]. ChineseJournal of Sensors and Actuators, 2012, 25(8): 1097–1101. DOI: 10.3969/j.issn.1004-1699.2012.08.016.
    [19] CROUSE M S, NOWAK R D, BARANIUK R G. Wavelet-based statistical signal processing using hidden Markov models [J]. IEEE Transactions onSignal Processing, 1998, 46(4): 886–902. DOI: 10.1109/78.668544.
    [20] 张毅刚, 郁惟镛, 黄成军, 等. 基于小波包及隐式马尔科夫模型的局放信号去噪 [J]. 上海交通大学学报, 2004, 38(8): 1269–1272. DOI: 10.3321/j.issn:1006-2467.2004.08.010.
    ZHANG Y G, YU W Y, HUANG C J, et al. Wavelet package-based partial discharge denoising using hidden markov model [J]. Journal of Shanghai Jiaotong University, 2004, 38(8): 1269–1272. DOI: 10.3321/j.issn:1006-2467.2004.08.010.
    [21] GRINSTED A, MOORE J C, JEVREJEVA S. Application of the cross wavelet transform and wavelet coherence to geophysical time series [J]. Nonlinear Processes in Geophysics, 2004, 11(5/6): 561–566. DOI: 10.5194/npg-11-561-2004.
    [22] 邵骏. 基于交叉小波变换的水文多尺度相关分析 [J]. 水力发电学报, 2013, 32(2): 22–26, 42.
    SHAO J. Multi-scale correlation analysis ofhydrological time series based on cross wavelet transform [J]. Journal of Hydroelectric Engineering, 2013, 32(2): 22–26, 42.
  • [1] 中国生房营光谷任国赵奎 . 基于小波分析的建(构)筑物爆破振动安全评估. 爆炸与冲击, doi: 10.11883/1001-1455(2009)01-0035-06
    [2] 费鸿禄曾翔宇杨智广 . 隧道掘进爆破振动对地表影响的小波包分析. 爆炸与冲击, doi: 10.11883/1001-1455(2017)01-0077-07
    [3] 张胜凌同华曹峰黄戡 . 模式自适应连续小波去除趋势项方法在爆破振动信号分析中的应用. 爆炸与冲击, doi: 10.11883/1001-1455(2017)02-0255-07
    [4] 邱贤阳史秀志周健黄丹陈新 . 基于HHT能量谱的高精度雷管短微差爆破降振效果分析. 爆炸与冲击, doi: 10.11883/1001-1455(2017)01-0107-07
    [5] 李维光张继春 . 爆破振动作用下顺层岩质边坡稳定性分析. 爆炸与冲击, doi: 10.11883/1001-1455(2007)05-0426-05
    [6] 杨年华张乐 . 爆破振动波叠加数值预测方法. 爆炸与冲击, doi: 10.11883/1001-1455(2012)01-0084-07
    [7] 路亮龙源谢全民李兴华纪冲赵长啸 . 爆破振动信号的提升小波包分解及能量分布特征. 爆炸与冲击, doi: 10.11883/1001-1455(2013)02-0140-08
    [8] 中国生敖丽萍赵奎 , . 基于小波包能量谱爆炸参量对爆破振动信号能量分布的影响. 爆炸与冲击, doi: 10.11883/1001-1455(2009)03-0300-06
    [9] 李洪涛卢文波舒大强朱传云 . P波作用下衬砌混凝土的爆破安全振动速度研究. 爆炸与冲击, doi: 10.11883/1001-1455(2007)01-0034-06
    [10] 杨招伟卢文波陈明严鹏胡英国刘美山吴新霞冷振东 . 基于实测爆破振动计算岩体介质P波品质因子. 爆炸与冲击, doi: 10.11883/bzycj-2019-0333
    [11] 杨建华卢文波陈明周创兵 . 岩石爆破开挖诱发振动的等效模拟方法. 爆炸与冲击, doi: 10.11883/1001-1455(2012)02-0157-07
    [12] 武旭张云鹏郭奇峰 . 台阶地形爆破振动放大与衰减效应研究. 爆炸与冲击, doi: 10.11883/1001-1455(2017)06-1017-06
    [13] 杨润强严鹏王高辉卢文波陈明 . 地应力水平对深埋隧洞爆破振动频谱结构的影响. 爆炸与冲击, doi: 10.11883/bzycj-2017-0366
    [14] 史秀志林大能陈寿如 . 基于粗糙集模糊神经网络的爆破振动危害预测. 爆炸与冲击, doi: 10.11883/1001-1455(2009)04-0401-07
    [15] 刘爱文俞言祥付长华吕红山陈鲲赵纪生周正华王伟 . 50t科学爆破振动的衰减特征及地形效应. 爆炸与冲击, doi: 10.11883/1001-1455(2010)01-0021-06
    [16] 马冲詹红兵姚文敏余海兵 . 爆破振动作用下含软弱夹层边坡稳定性及安全判据. 爆炸与冲击, doi: 10.11883/bzycj-2016-0275
    [17] 柴修伟 . 不同炸药性能对爆破振动效应的影响. 爆炸与冲击, doi: 10.11883/1001-1455(2011)05-0548-05
    [18] 吴从师徐荣文张庆彬 . 自由面对爆破振动信号能量分布特征的影响. 爆炸与冲击, doi: 10.11883/1001-1455(2017)06-0907-08
    [19] 施建俊李庆亚张琪卫星王辉 . 基于Matlab和BP神经网络的爆破振动预测系统. 爆炸与冲击, doi: 10.11883/1001-1455(2017)06-1087-06
    [20] 楼晓明周文海简文彬郑俊杰 . 微差爆破振动波速度峰值-位移分布特征的延时控制. 爆炸与冲击, doi: 10.11883/1001-1455(2016)06-0839-08
  • 加载中
图(14)表(2)
计量
  • 文章访问数:  179
  • HTML全文浏览量:  179
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-09-23
  • 录用日期:  2020-02-11
  • 网络出版日期:  2020-08-25

冻结立井爆破近区井壁振动信号基线漂移校正和消噪方法

    作者简介:付晓强(1984- ),男,博士,讲师,fuxiaoqiang1984@163.com
  • 1. 三明学院建筑工程学院,福建 三明 365004
  • 2. 三明学院工程材料与结构加固福建省高等学校重点实验室,福建 三明 365004
  • 3. 北京科技大学土木与资源工程学院,北京 100083

摘要: 冻结立井爆破过程中,近区监测信号中含有的基线漂零及噪声成分对其局部特征精细化提取影响显著。在对近区井壁振动信号有效采集基础上,通过互补总体经验模态分解(complementary ensemble empirical mode decomposition, CEEMD)方法、稀疏化基线估计消噪(baseline estimation and de-noising with sparsity, BEADS)方法和隐马尔可夫模型消噪(hidden Markov model de-noising, HMMD)方法等,解决了信号中基线漂移和随机噪声消除难题,并采用交叉小波变换对校正和消噪效果进行了相关性评价。实例分析结果表明:信号中缓变的基线成分遍历信号各个模态分量的整个过程,且主要集中于低频分量中,而噪声则集中在高频分量。组合分析方法对低频基线漂零和高频噪声的处理效果好,是一种高效且相对保幅的信号分析方法,可用于批量信号数据的预处理过程。

English Abstract

  • 冻结井筒采用钻爆法施工时,新浇混凝土强度较低,复杂受力下局部极易产生应力集中,形成不同方向、不同尺度的裂隙和缺陷,这些累积损伤会造成混凝土井壁强度降低、各项物理力学参数与设计初值相比发生较大差异等问题,最终造成井壁承载能力下降,致使井壁位移量过大,易诱发冻结壁破坏、冻结管断裂等安全事故[1-3]。开展井壁振动监测是评估爆破条件下井壁振动响应特征的有效手段之一。振动测试过程中,由于测试环境影响、仪器系统自身原因及监测方案存在缺陷,监测振动波形会局部存在偏离基线中心的现象,即信号基线漂零。尤其在距离爆心位置相对较近的爆破近区,采集信号基线偏离问题更加显著。根据经典的萨道夫斯基公式,对于同一类型爆破工程而言,在理想测试环境下,若振速幅值相当,则最大段药量越大,信号中存在趋势项的概率也相对较大,相应的“近区”的范围也越大。更进一步地,若信号同时存在高频噪声,则两者的不利因素相互叠加,若不采用有效的方法消除上述影响,则整个信号会受到污染,严重影响信号特征的准确提取。因此,在信号预处理时必须予以基线校正和噪声消除,避免对信号相关特征分析的干扰。

    针对非线性振动信号存在干扰成分的问题,相关学者开展了具体的研究,如:张军等[4]针对车载武器振动信号中趋势项消除的难题,根据信号自身特点,研究了基于经验模态分解(empirical mode decomposition, EMD)的信号趋势项提取和判别方法;张景元等[5]通过对导弹自动驾驶仪振动信号进行分析,提出了基于形态学的解决振动信号中基线漂移问题的滤波方法,有效抑制了信号中的基线成分;刘艳丽等[6]基于形态学滤波方法,对脉搏波信号的基线漂移进行了校正,为准确诊断病人的亚健康状态提供了技术手段;谢芳娟等[7]为了消除基线漂移及噪声对信号分析的干扰,提出了基于空域追踪的修正算法,并通过实例对算法的可靠性进行了验证。

    对爆破振动信号干扰性处理的研究,主要集中在趋势项的修正和滤波去噪方面,如:张胜等[8]通过对实测爆破信号进行分析,提出了以时域积分后的振速信号构造模式自适应小波基的方法,并成功地去除了将爆破信号中的趋势项;韩亮等[9]利用集合经验模态分解(ensemble empirical mode decomposition, EEMD)及小波相关方法,提出了以固有模态函数频带分布为指标进行人工判别趋势项的去除方法,以及基于自相关分析识别噪声特征的小波阈值去噪方法,并应用该方法进行露天深孔台阶爆破信号趋势项去除,收到了良好的效果;王志亮等[10]采用集合经验模态分解与高低频处理相结合的方法,对花岗岩爆破振动加速度信号中的漂零信号进行了修正,并对该算法的局限性进行了进一步的讨论。

    爆破振动信号存在基线漂零和噪声的主要因素包括仪器自身原因、测试环境(温度变化、飞石及环境干扰)及测试方法等,其中仪器自身原因、测试环境属于测试中的不可控因素,而测试方法的优化能很大程度上降低甚至避免上述现象对信号波形的干扰。本文针对冻结立井的井壁振动测试,提出传感器井壁预埋振动监测方法,以实现冻结立井爆破近区井壁振动信号的实时准确监测;对爆破近区井壁振动信号中存在的明显基线漂零现象进行校正,并采用最优化去噪方法对校正信号进行消噪分析;通过几种去噪方法评价指标对比和交互小波变换相关性分析,检验基线修正和消噪效果,以期为爆破振动等非线性信号基线漂零校正和噪声消除提供新思路。

    • 在经典的希尔伯特-黄变换基础上,Yeh等[11]提出互补总体经验模态分解(complementary ensemble empirical mode decomposition, CEEMD)方法。该方法是在传统EMD分析的基础上改进的一种噪声自适应的完备算法[12],利用噪声进行辅助分析,可有效解决EMD、EEMD算法中的模态混叠问题。该方法的核心在于将原始信号中添加一对相反的白噪声信号,对生成的信号分别进行EMD分解,将分解得到的固有模态函数(intrinsic mode function, IMF)进行优选重组,以得到最终信号。该算法在保证分解效果的同时,有效抑制了由白噪声引起的重构误差。

      该算法步骤为[13]

      (1) 向原始爆破信号S中加入正、负成对的辅助白噪声,从而生成两个信号S+S-

      $ \left[ \begin{array}{l} {S_ + }\\ {S_ - } \end{array} \right] = \left[ {\begin{array}{*{20}{c}} 1&1\\ 1&{ - 1} \end{array}} \right]\left[ \begin{array}{l} S\\ N \end{array} \right] $

      式中:N为辅助噪声;

      (2) 对S+S-做EMD分解,分别得到两组IMF函数,记为Si+Si-

      $ {S_ + } = \frac{1}{n}\sum\nolimits_{i = 1}^n {{S_{i + }}} ,\quad\quad\quad{S_ - } = \frac{1}{n}\sum\nolimits_{i = 1}^n {{S_{i - }}} $

      (3) 通过将两组分量组合便得到信号S分解形式为:

      $ S=\frac{1}{2}\left({S}_{1}+{S}_{2}\right)=\sum\nolimits_{i=1}^{n}{S}_{i}+R $

      式中:Si为信号S的模态函数分量IMF iR为残余项。

    • 稀疏化基线估计消噪(baseline estimation and de-noising with sparsity, BEADS)算法具体过程如下。

      对于任意给定的稀疏化信号s,若其中含有n点随机基线成分,则该信号可视为特征波形信号和缓变基线漂移信号的叠加[14],即:

      ${{s}} = {{x}} + {{f}}\;\;\;\;\;{{s}} \in {{\bf{R}}^n}$

      式中:x为多峰值且稀疏可导的信号矩阵;f为信号中含有的基线成分,为低通信号。可将观察到的受基线漂移和噪声影响的信号进一步分解为:

      ${{y}} = {{s}} + {{w}} = {{x}} + {{f}} + {{w}}\;\;\;\;\;{{y}} \in {{\bf{R}}^n}$

      式中:w是方差为${\sigma ^2}$的平稳高斯白噪声,信号分析最终目的是剔除信号中含有的基线成分并保留含有峰值点的特征信号。若假定信号不存在峰值,则基线可通过低通滤波器从受噪声污染的观察信号中近似恢复,即:

      ${{f}} \approx {{L}}\left( {{{y}} + {{\hat x}}} \right)$

      式中:L为与信号匹配的低通滤波变换,$\hat x$为最终得到优化后的特征信号估计值。

    • 交叉小波变换(cross wavelet transform, CWT)融合了小波变换和交叉谱分析的功能,可用以揭示两个随机信号在时域和频域中的相关关系。假定$ {{W}}_{n}^{{Y}}\left(t\right) $$ {{W}}_{n}^{{X}}\left(t\right) $和为两个随机信号X={x1x2,···,xn}与Y={y1y2,···,yn}的连续小波变换,则两个信号间的交叉小波变换可表示为[15]

      ${{W}}_n^{{{XY}}}(t) = {{W}}_n^{{X}}(t){{W}}_n^{{{Y}}*}(t)$

      式中:${{W}}_n^{{{Y}}{\rm{*}}}(t)$${{W}}_n^{{Y}}(t)$的复共轭,t为信号时移。交叉小波谱可定义为$ {{W}}_{n}^{{XY}}\left(t\right) $的范数${\rm{||}}{{W}}_n^{{{XY}}}(t){\rm{||}}$,它清晰展示了两个信号在时频空间域上相关程度。在时频域上的值越大,表明两个信号之间的相关程度越高。对于任意给定的两个非平稳随机信号,交叉小波变换的互相关系数为:

      $r({{X}},{{Y}}) = \frac{{\displaystyle\sum\limits_{i = 1}^n {\left[ {{{W}}_i^{{X}}(t) - \overline {{{W}}_i^{{X}}(t)} } \right]\left[ {{{W}}_i^{{Y}}(t) - \overline {{{W}}_i^{{Y}}(t)} } \right]} }}{{\sqrt {\displaystyle\sum\limits_{i = 1}^n {{{\left[ {{{W}}_i^{{X}}(t) - \overline {{{W}}_i^{{X}}(t)} } \right]}^2}} } \sqrt {\displaystyle\sum\limits_{i = 1}^n {{{\left[ {{{W}}_i^{{Y}}(t) - \overline {{{W}}_i^{{Y}}(t)} } \right]}^2}} } }}$

      式中:上标的“—”表示对所有分量求平均值。

      这里,选取Morlet小波作为小波基函数,以反映两个爆破振动信号经小波变换后在时频域上相关性振荡的时频结构特征。

    • 对于给定的两个时间序列信号XY,假定其期望谱为傅里叶红噪声谱$P_k^X$$P_k^Y$,则其交叉小波功率谱分布为[15]

      $D\left( {\frac{{{\rm{|}}{{W}}_n^{{X}}(t){{W}}_n^{{{Y}}*}(t)|}}{{{\sigma _{{X}}}{\sigma _{{Y}}}}} {\text{<}} p} \right) = \frac{{{Z_v}(p)}}{\nu }\sqrt {P_k^{{X}}P_k^{{Y}}} $

      式中:|W|表示矩阵W的行列式;σxσy分别为时间序列的标准差;Zv(p)是与概率p有关的置信度;对于实小波,$v = 1$,对于复小波,$v = 2$。在显著性水平α=0.05条件下,Z1(95%)=2.182,Z2(95%)=3.999,则认为通过显著性水平α=0.05条件下红噪声标准谱的检验,即两者之间显著相关。

    • 根据上述算法理论,整理得到冻结立井爆破近区井壁振动信号基线校正和噪声消除方法,利用matlab软件编制程序,可实现信号的批量化处理,具体如图1所示。

      图  1  基线校正及消噪流程

      Figure 1.  Baseline correction and noise reduction process

    • 与巷道有所不同,冻结立井施工作业空间狭小,井下环境恶劣,尤其是冻结低温的条件和高强混凝土初凝期的水化热现象会对测振仪器的性能和准确性产生负面影响。上述因素对立井爆破振动的现场监测带来很大挑战。针对测试过程中存在的问题,本次现场试验对传统测试方法进行改进:采用传感器预埋安装(绑扎固定在井壁钢筋交叉点处)、大容量专用防爆电源长期供电、专用隔爆保温防护箱等技术措施,集成开发了冻结立井爆破振动监测系统,减少电线敷设、检查和回收等大量工作,避免施工对振动测试的干扰,以实现冻结立井井壁爆破振动信号的长期、稳定和实时监测。具体布置方案如图2所示。

      图  2  传感器井壁预埋法

      Figure 2.  Pre-embedding method of vibration instrument

      采用井壁预埋法在兖矿集团万福煤矿对主立井爆破掘进过程中井壁振动进行了有效监测。该井筒净直径为5.5 m,深度886 m,冻结深度894 m,井筒施工需穿过753 m的表土层,是地质条件极为特殊的井筒,表土层厚度目前是世界第一。测试段为双层井壁结构,钢纤维高强混凝土井壁强度为CF90,单侧厚度1.9 m,测振传感器固定在内层井壁。爆破采用直眼掏槽形式,1~5段毫秒电雷管,掏槽眼深度4.2 m,其余炮眼深度4 m,单循环炮眼总数135个,总装药量339.9 kg,炮眼具体布置如图3所示。

      图  3  炮眼布置(单位:mm)

      Figure 3.  Borehole layout (unit: mm)

      通过对不同测距处的信号对比发现,在爆破近区处采集到的井壁振动信号基线漂零现象显著,如距离立井爆破掌子面3 m处,采集到的井壁垂向振动信号时程曲线如图4所示。该振动信号波峰值为6.35 cm/s,波谷值为−4.78 cm/s,峰峰值为11.13 cm/s。从图中可知:监测信号中包含一定的杂波成分,并且信号波形存在明显的基线漂零(“甩尾”)现象。信号在时间轴上0.1 s内监测波形稳定,随后存在严重的基线偏离问题,信号历时波形在时间轴0.3 s后逐渐回归基线中心并最终趋于稳定。由于爆破网络中周边眼所用的MS5段雷管标定延期时间为(110±15) ms,超出了稳定波形所在时间轴范围。因此,基线漂零会对爆破信号分析产生不利影响。

      图  4  爆破近区井壁振动信号

      Figure 4.  Vibration signal of shaft lining near blasting area

    • 由于信号分析过程中的不同处理阶段的误差会逐步累积,所以分析过程中的各个步骤分析方法均需优选,确定最适合信号特点的分析方法,从而最大程度保留和还原信号的固有特征属性。本次分析设置初始白噪声标准差为0.2,筛分次数为5,最大筛分迭代次数为100。经过5次筛分,各IMF分量迭代次数分布情况如图5所示。

      图  5  各IMF分量筛分迭代次数关系

      Figure 5.  The relationship of iterations number and shift number

      图5表明,信号的IMF分量越复杂,迭代次数也越多,包含的信息越丰富,对原始信号信息的继承度也越高。最终,信号经CEEMD运算被分解为11个IMF分量和1个残余项R,具体如图6所示。

      图  6  信号CEEMD分解各分量及残余项R

      Figure 6.  Component of CEEMD decomposition and residual signal

      本次CEEMD算法中,按照信号的复杂程度和频率高低将其分解为若干个独立分量,分解得到的IMF1~IMF7分量端点振荡效应和模态混叠已基本消除,但具有典型的含噪信号特点,即多峰值、多频带随机特征。而剩余的其余低频分量在波形初始零时刻(或波形尾部截止时刻),均较基线中心位置出现显著的偏离(椭圆线所标记处)。由于信号中存在缓变的基线趋势项,需要建立低通滤波器予以校正和消除。为了客观评价各IMF分量与原信号的相关度,这里将趋势项R也作为一个分量(IMF12),利用互相关函数确定其对应的相关度值,见表1

      IMF相关性系数IMF相关性系数IMF相关性系数
      10.420250.502390.1083
      20.448460.1419100.3707
      30.375470.1333110.4540
      40.412480.1364120.2940

      表 1  模态分量与原信号相关度

      Table 1.  Correlation between components and original signals

      表1中可知:IMF1~IMF5分量的相关性系数远大于IMF6~IMF9分量,同时应注意到IMF10、IMF11分量及残余项R(趋势项)与原信号的相关性也较大,如IMF11分量值甚至超过IMF1~IMF4分量,说明传统的CEEMD分析中直接将残余项及相关分量舍弃的信号重构方法,会导致信号有效成分的严重丢失。从图6中也可看出趋势项幅值较大,因此在分析过程中不能完全舍弃,应采用合理的算法将其中包含的有效信息进行提取,避免导致信号成分的缺失。

    • 对于给定的基线偏移和含噪信号y,可将其建模为m个基本波形(或称为信号原子)的线性组合,从而得到其稀疏化的表达形式:

      ${{y}} = {{\Phi x}} + {{w}} = \sum\limits_{j = 1}^m {{x_j}{{{\varphi}} _j}} + {{w}}$

      式中:x为待提取信号y的稀疏化形式;w为信号中包含的标准差σ>0的高斯白噪声;Φn×m的过完备稀疏表示系数矩阵;xj (j=1, 2, 3,···, m)为y在字典Φ=[φ1,φ2,···,φm]下的稀疏表示系数;信号原子φjΦ的列向量,一般将其归一化到单位${\ell _2}$范数,即:

      ${\rm{||}}{{{\varphi }}_j}|{|_2} = \sum\limits_{i = 1}^n {|\varphi _{ij}^2|} = 1$

      式中:ϕijΦ的分量,i=1, 2, 3,···, n

      根据式(10),字典Φ下稀疏化信号y可以通过信号原子簇φj中的部分原子的叠加而准确描述。

      在分析过程中,为了降低基线校正过程对信号分量幅值的影响,引入决定该过程中信号稀疏化成都的正则化系数$\lambda $和具有非对称补偿罚值功能的对称罚函数$\phi :{\bf{R}} \to {\bf{R}}$,则问题转化为下述优化问题,也称为基追踪去噪(basic pursuit de-noising, BPD)过程[14]

      $ {\arg _x}\min \frac{1}{2}||{{y}} - {{\Phi x}}||_2^2 + \lambda R({{x}})\;\;\;\;R({{x}}) = \sum\limits_m {\phi ({x_m})} $

      由于基线成分主要位于信号低频分量中,分析时设置信号截止频率fc为0.002 Hz,滤波器阶数d取为0~2,罚函数非对称性参数r取为6。正则化参数幅值为0.8,不同阶正则化参数${\lambda _0}{\text{~}}{\lambda _2}$分别为:0.4、3.2和4。

      根据上述理论,对图中分解得到的12个固有模态分量分别进行BEADS算法处理。限于篇幅,这里仅给出IMF3、IMF8~IMF11和残余分量R共6个分量基线成分和校正后的信号曲线,见图7图7表明,BEADS算法可准确识别各模态分量信号中包含的缓变的低通基线杂波成分,且信号不存在明显的畸变和局部细节特征缺失。在相对高频段分量(IMF1~IMF4)中的基线成分主要位于主振时刻内(小于0.1 s)且幅值均较小;相对中频段分量(IMF5~IMF9)中的基线成分主要位于主振时刻内,幅值与高频段相比较大;而相对低频段分量(IMF10~IMF11及R)中含有的基线成分幅值较大且持续时间较长,为信号中趋势项的主要来源。

      图  7  模态分量处理结果

      Figure 7.  Results of baseline correction processing

      图8为各模态分量分离出的其中含有的基线成分。说明基线成分存在并遍历CEEMD分解的整个过程。低频分量中基线偏移幅值比高频分量相对较大,具备典型的低通特性,这与图7的分析结果一致。图9给出了各分量校正过程中罚函数值与迭代次数的历史关系曲线。从图9可以看出,在迭代有限次数(小于20次)罚函数值便趋于收敛,这验证了算法运算速度和特征提取的有效性。

      图  8  各分量基线成分

      Figure 8.  Baseline of each component

      图  9  罚函数值与迭代次数关系

      Figure 9.  The relation of cost function and the iteration number

      采用BEADS算法进行基线校正后的波形时程曲线如图10所示。参考相关文献并结合相关性分析结果,剔除信号趋势项R(IMF12分量)并将相关性系数较大的IMF1~5、IMF10~11分量重新组合便得到了CEEMD方法处理后的重构信号,如图11所示。

      图  10  基线校正后信号

      Figure 10.  Baseline corrected signal

      图  11  CEEMD重构信号

      Figure 11.  reconstructed signals of CEEMD

      对比图10图11可知:剔除残余项及相关低频分量的信号处理方法,使得重组信号幅值降低,波峰值和波谷值分别为3.68 cm/s、4.35 cm/s,峰峰值为8.01 cm/s,与原信号幅值出现了很大偏差;同时,高频分量干扰噪声的残留及低频有效成分的缺失现象均较为严重,在实际分析中需要进行相应的处理,这增加了分析难度和复杂性。

    • BEADS算法对信号中的基线漂零和低频噪声处理效果较好,但无法完美地滤除信号中的高频噪声。通常认为信号中包含的高频噪声满足关于时间的泊松分布,高频噪声的幅值与信号的强度(峰值+基线)的平方根成正比。为了寻求最优信号消噪效果,本文分别采用形态学消噪(morphological de-noising, MD)、奇异值消噪(singular value de-noising, SVD)、小波熵消噪(wavelet entropy de-noising, WED)和隐马尔可夫模型消噪(hidden Markov models de-noising, HMMD)四种方法对图中的基线校正信号进行消噪处理,以上消噪方法的相关理论见文献[16-20]。MD方法充分利用信号波形的形态特征实现其消噪过程,具有运算速度快,自适应强的特点;SVD方法通过构建信号奇异值矩阵,利用信号和噪声的能量可分性,在信号重构过程中将光滑信号所产生的奇异值保留、噪声信号奇异值置零,从而实现消噪信号特征的最优化估计;WED方法中,选取“db8”小波基函数,进行3层分解,计算得到各层的熵值,最后利用wdencmp函数实现对信号噪声的压制;HMMD方法是基于统计数据而提出的新的信号处理方法,该方法不受人为门限设置的影响,已在信号处理相关领域得到成功应用。为了获得最佳分析信号,本文采用信噪比(signal to noise rate, SNR)、均方根误差(root mean square error, RMSE)、相关性系数(coefficient of correlation, CC)和峰值误差(peak error, PE)四个指标来综合评价不同去噪方法的消噪效果。不同方法消噪效果如表2所示。

      去噪方法评价指标
      SNRRMSECCPE
      MD3.4420.5810.8161.758
      SVD20.9830.3790.8961.175
      WED21.3970.3740.9210.492
      HMMD43.3710.0560.9920.244

      表 2  不同方法消噪性能指标

      Table 2.  Indexes of noise reduction performance

      表2中各项指标的对比可知:MD方法避免了小波阈值选取及信号重构,具有自适应特点,但迭代次数选择不合理易导致消噪信号波形失真;SVD方法中奇异值数目的选择对去噪结果影响加大,当奇异值数目较少时,降噪阶次低导致信号信息不完整,反之,降噪阶次高,噪声无法充分消除;WED方法存在小波基选取和分解层数的限制,消噪效果具有一定的盲目性。相比于前三种消噪方法,HMMD方法提高了运算过程中的收敛性,误差小精度高,使得原始信号的信息损失量最小,对局部信号的消噪效果好,消噪后的校正波形见图12

      图  12  消噪后的校正信号

      Figure 12.  Baseline corrected signal after de-noising

      为了验证去噪效果,分别求取原始信号、CEEMD重组信号和HMMD消噪后的校正信号的功率谱密度(wf)曲线,结果如图13所示。图13中原始信号中含有的低频基线成分致使信号在3 Hz处出现峰值突高的奇异点,幅值为0.67 dB/Hz,远大于图13(c)中的主频峰值,且信号含有的基线成分主要位于f≤15 Hz低频段内;图13(b)中,采用CEEMD重组信号时,虽然对相关分量进行了相关性判别和取舍,但对高、频相关分量的处置过于随意,导致在信号的相对高频部分(大于500 Hz)中广泛分布有微幅的随机噪声,同时在信号主振频域范围内(小于500 Hz)噪声特征也较为明显;图13(c)中经过HMMD方法消噪后的校正信号,频谱特征清晰完整,谱峰值与图13(b)幅值相当,主频由低频缓慢过渡,其中第一优势主频为16 Hz,第二优势主频为48 Hz,第三优势主频为103 Hz,具有多段别雷管起爆多次能量输入的典型频谱特征。综上所述,信号中包含的基线漂零现象的存在,会导致对信号主频峰值的误判,同时高频噪声问题对优势主频数量的确定具有显著影响。

      图  13  信号功率谱密度(wf

      Figure 13.  Power spectral density (wf) of signal

    • 小波多尺度相关分析可以展现不同信号总体上的相关程度,清晰描述信号在不同时域和频域上的局部细节相关特征[21]。小波多尺度相关谱综合反映了不同信号之间的相关性在时域和频域上的依赖关系,揭示不同去噪信号与原始信号在不同时间和频率尺度上的相关程度和细部特征。

      为了便于分析,标识原始信号、CEEMD重构信号和基线校正消噪信号为X1X2X3图14分别为X2X3X1的小波相关性凝聚谱,图中的粗实线区域表示通过显著性水平$\alpha $=0.05条件下的红噪声标准谱的检验,箭头表示两个信号之间的位相关系,其中→表示两个信号之间为同位相,说明两者为正相关关系;←表示反位相,说明两者之间为负相关关系[22]。为避免边界效应及小波高频虚假信息,小波影响锥区域(细实线)以内为有效谱值。谱分析计算过程中,对相关性数据做归一化处理。

      图  14  信号相关性凝聚谱

      Figure 14.  Correlation condensation spectra of signals

      相关性凝聚谱说明:X2X3X1在不同的时域和频域尺度上均存在一定相关度。X2X1在高频部分(大于512 Hz)正相关性显著,信号中的高频噪声未得到有效的滤除;同时X2重构过程中直接舍弃相关低频分量,导致信号低频部分信息缺失;相关性在局部出现间断、不连续现象,对信号高频噪声和低频干扰项的保留度均较高。采用文中方法得到X3与原始信号X1的相关性在时频域连续且完整,在低频尺度正相关性最为明显,显著性检验基本贯穿小波影响椎范围内的整个时频域。采用传统的互相关计算X2X3与原始信号X1之间的相关系数分别为:0.996、0.853。验证了组合分析方法的可靠度和有效性。综上分析可知,采用交叉小波变换进行不同信号序列间的相关性分析是可行的,且能够更加清晰描述两信号序列相关关系在时域和频域变化的细部特征和共振位相差异。

    • 文中主要针对深厚表土层冻结立井基岩段爆破近区井壁振动信号开展研究,由于爆破参数的差异和井筒地质条件的复杂性,对于不同爆破条件下信号的基线漂移校正应进行相关参数的调整,以期达到最优化的分析效果。同时,随着信号各分量中基线成分的弱化,其中含有的噪声分量会被一定程度上被放大,因此,滤波消噪方法的选择对分析结果具有显著影响。在后续的研究中应注重从仪器选择、参数设置和测点布置等方面入手,从根源上避免基线漂零现象和环境噪声的产生,深入开展多种类型爆破信号的研究,以便得到更具普适性的研究成果。

    • 通过对冻结立井爆破近区井壁振动信号基线和噪声分析处理,得到以下结论。

      (1) 测试方法的优化可一定程度上降低和控制基线漂零现象对波形特征提取的影响。传感器井壁预埋方法克服了以往井壁振动监测方法的诸多弊端,实现了井壁结构振动响应持续监测。在实际测点布置时,应重点做好仪器相关配件的防水工作,另外,对于有瓦斯涌出的矿井,应确保仪器配套供电系统的防爆性。

      (2) 冻结立井爆破近区井壁振动信号中基线漂零现象严重,基线成分分布在模态分解的各个分量中,且主要位于f≤15 Hz低频段内,基线偏移会导致信号主频的误判,BEADS方法可有效修正信号中包含的低频基线成分。基线分量的弱化会使得信号中的高频噪声在一定程度上得到放大,应采用合理的去噪方法予以消除。HMMD消噪模型自适应性强、收敛速度快,可很大程度上消除信号中的高频噪声。

      (3) 交叉小波变换相关性分析可揭示两信号序列相关关系在时域和频域变化的细部特征和共振位相差异。CEEMD重构信号易导致信号有效成分的缺失,进而导致信号振幅降低,同时高频成分保留较多。组合方法有效削弱了干扰成分对信号特征提取的影响,在信号主振时频域与原始信号关联性高,适合用于该类信号的分析处理。

参考文献 (22)

目录

    /

    返回文章
    返回