• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

基于改进MP-WVD算法的核电厂建设爆破振动信号处理方法

梅比 汪旭光 杨仁树

引用本文:
Citation:

基于改进MP-WVD算法的核电厂建设爆破振动信号处理方法

    作者简介: 梅 比(1981- ),男,博士,高级工程师,blastmei2017@163.com;
  • 中图分类号: O389

Blasting vibration signal analysis technology of construction of nuclear power plant based on improved MP-WVD algorithm

  • CLC number: O389

  • 摘要: 针对爆破施工产生的振动问题,结合现有信号分析方法优缺点,引入MP-WVD组合算法解决信号分析过程中交叉项干扰问题,可提高时频分辨率,满足精确提取核电爆破振动响应信号时频分布特征的需求;将HHT算法引入MP算法数据的预处理,成功降低了MP算法运算复杂度,可为大数据量分析奠定基础。运用该算法对漳州核电一期工程场地平整土石方爆破工程振动监测信号进行分析,精确提取了爆破振动信号的时频特征。
  • 图 1  分解过程示意图

    Figure 1.  Schematic diagram of decomposition process

    图 3  WVD时频分布图

    Figure 3.  WVD time-frequency distribution

    图 2  示例信号x(t)

    Figure 2.  Example signal x(t)

    图 4  MP算法分解结果

    Figure 4.  Decomposition results of MP algorithm

    图 5  各个子波的WVD分布

    Figure 5.  WVD distribution of each wavelet

    图 6  子波分解的WVD时频分布图

    Figure 6.  WVD time-frequency distribution of wavelet decomposition

    图 7  程序框图

    Figure 7.  Program diagram

    图 8  原始信号

    Figure 8.  Original signal

    图 9  子波集合

    Figure 9.  Collection of wavelets

    图 10  滤波效果对比

    Figure 10.  Comparison of filtering effects

    图 11  MP-WVD算法结果

    Figure 11.  Result of MP-WVD algorithm

    图 12  小波变换结果

    Figure 12.  Result of wavelet transform

    图 13  水平径向分量(X轴)

    Figure 13.  Horizontal radial component

    图 14  水平切向分量(Y轴)

    Figure 14.  Horizontal tangential component

    图 15  瞬时能量谱

    Figure 15.  Instantaneous energy spectrum

    图 16  频率边际谱

    Figure 16.  Marginal spectrum

    图 17  南山村测点数据分析结果

    Figure 17.  Data analysis result of Nanshan village

    图 18  人家村测点数据分析结果

    Figure 18.  Data analysis result of Renjia village

    表 1  两种方法所用机时对比

    Table 1.  Comparison of machine time between two methods

    合成信号类型信号波形传统 MP 所用机时/s改进 MP 所用机时/s
    单正弦96.111.90
    双正弦77.131.90
    三正弦82.961.71
    典型爆破振动信号82.081.58
     注:信号长度,N=512,残差截断阈值,σ=0.05,CPU 型号,Intel Core i7-4710MQ,主频 2.5~3.1 GHz。
    下载: 导出CSV

    表 2  爆破参数表

    Table 2.  Blasting parameters

    孔径/mm孔数孔深/m孔距/m排距/m抵抗线/m孔药量/kg段药量/kg总药量/kg
    1154713.05.53.22.285.0170.03 984
    下载: 导出CSV
  • [1] 方群才, 郑国纲. 漳州核电厂一期工程厂址区域核应急方案编制介绍 [J]. 科技传播, 2016, 8(5): 168–169. DOI: 10.16607/j.cnki.1674-6708.2016.05.097
    FANG Quncai, ZHENG Guogang. Introduction of nuclear emergency plan for site area of phase I project of Zhangzhou nuclear power plant [J]. Public Communication of Science and Technology, 2016, 8(5): 168–169. DOI: 10.16607/j.cnki.1674-6708.2016.05.097
    [2] 郭涛, 方向, 谢全民, 等. 频率切片小波变换在爆破振动信号时频特征精确提取中应用 [J]. 振动与冲击, 2013, 32(22): 73–78. DOI: 10.13465/j.cnki.jvs.2013.22.011
    GUO Tao, FANG Xiang, XIE Quanmin, et al. Application of FSWT in accurate extraction of time-frequency features for blasting vibration signals [J]. Journal of Vibration and Shock, 2013, 32(22): 73–78. DOI: 10.13465/j.cnki.jvs.2013.22.011
    [3] 马华原, 龙源, 郭涛. 田湾核电站基础开挖爆破振动检测与数据分析 [J]. 爆破, 2015, 32(3): 85–90 doi: 10.3963/j.issn.1001-487X.2015.03.015
    MA Huayuan, LONG Yuan, GUO Tao. Monitoring and analysis of blasting vibration of foundation excavation of tianwan nuclear power Plant [J]. Blasting, 2015, 32(3): 85–90 doi: 10.3963/j.issn.1001-487X.2015.03.015
    [4] 李夕兵, 张义平, 左宇军, 等. 岩石爆破振动信号的EMD滤波与消噪 [J]. 中南大学学报(自然科学版), 2006(1): 150–154 doi: 10.3969/j.issn.1672-7207.2006.01.029
    LI Xibing, ZHANG Yiping, ZUO Yujun, et al. EMD filtering and denoising of rock blasting vibration signal [J]. Journal of Central South University (Science and Technology), 2006(1): 150–154 doi: 10.3969/j.issn.1672-7207.2006.01.029
    [5] 魏新江, 谢超, 丁玉琴. 基于平均频率和HHT变换的隧道爆破震动信号研究 [J]. 矿业研究与开发, 2017, 37(7): 13–18
    WEI Xinjiang, XIE Chao, DING Yuqin. Study on the signals of tunnel blasting vibration based on average frequency and HHT method [J]. Mining Research and Development, 2017, 37(7): 13–18
    [6] 周辉, 龙源, 钟明寿, 等. 基于双参数MP算法的不同孔深爆炸地震波特性研究 [J]. 振动与冲击, 2016, 35(18): 76–81. DOI: 10.13465/j.cnki.jvs.2016.18.013
    ZHOU Hui, LONG Yuan, ZHONG Mingshou, et al. Characteristics analysis of explosion seismic waves with different hole depth based on the method of double parameters matching pursuit [J]. Journal of Vibration and Shock, 2016, 35(18): 76–81. DOI: 10.13465/j.cnki.jvs.2016.18.013
    [7] 贾亚飞, 朱永利, 王刘旺. 基于VMD和Wigner-Ville分布的局放信号时频分析 [J]. 系统仿真学报, 2018, 30(2): 569–578
    JIA Yafei, ZHU Yongli, WANG Liuwang. Time-frequency analysis of partial discharge signal based on VMD and Wigner-Ville distribution [J]. Journal of System Simulation, 2018, 30(2): 569–578
    [8] 井爱雯, 刘云, 马轶丽. 基于MP算法的语音信号稀疏分解 [J]. 计算机工程与应用, 2009, 45(5): 144–146 doi: 10.3778/j.issn.1002-8331.2009.05.042
    JING Aiwen, LIU Yun, MA Yili. Speech signal sparse decomposition based on matching pursuit algorithm [J]. Computer Engineering and Applications, 2009, 45(5): 144–146 doi: 10.3778/j.issn.1002-8331.2009.05.042
    [9] 张汛汛, 张繁昌, 刘汉卿. 基于快速匹配追踪算法的地震道集剩余时差校正 [J]. 石油物探, 2015, 54(4): 420–426 doi: 10.3969/j.issn.1000-1441.2015.04.008
    ZHANG Xunxun, ZHANG Fanchang, LIU Hanqing. Seismic gathers residual moveout correction based on fast matching pursuit algorithm [J]. Geophysical Prospecting for Petroleum, 2015, 54(4): 420–426 doi: 10.3969/j.issn.1000-1441.2015.04.008
    [10] 蒋雷龙, 曾小平, 汪胜金, 等. 台山核电厂安全控制爆破技术研究 [J]. 铁道科学与工程学报, 2010, 7(5): 122–125 doi: 10.3969/j.issn.1672-7029.2010.05.024
    JIANG Leilong, ZENG Xiaoping, WANG Shengjin, et al. Research on safety control blasting technology of Taishan nuclear power plant [J]. Journal of Railway Science and engineering, 2010, 7(5): 122–125 doi: 10.3969/j.issn.1672-7029.2010.05.024
    [11] 范兴利, 成谷. 基于Morlet小波尺度参数寻优的匹配追踪时频分析 [J]. 中山大学学报(自然科学版), 2014, 53(6): 85–92
    FAN Xingli, CHENG Gu. Matching pursuit time-frequency analysis based on morlet wavelet scale parameter optimization [J]. Acta Scientiarum Naturalium Universitatis Sunyatseni, 2014, 53(6): 85–92
    [12] 钟明寿, 周辉, 刘影, 等. 基于改进匹配追踪算法的化爆地震波信号时频特征提取 [J]. 爆炸与冲击, 2017, 37(6): 931–938
    ZHONG Mingshou, ZHOU Hui, LIU Ying, et al. Time-frequency analysis of explosion seismic signal based on improved matching pursuit [J]. Explosion and Shock Waves, 2017, 37(6): 931–938
    [13] 魏建国, 张忠, 张心斌. 核电站筏基大体积混凝土温控监测及仿真分析 [J]. 工业建筑, 2017, 37(6): 931–938
    WEI Jianguo, ZHANG Zhong, ZHANG Xinbin. Temperature-controlling monitoring and emulation analysis of mass-concrete of raft foundation of nuclear power plant [J]. Industrial Construction, 2017, 37(6): 931–938
    [14] 方群才, 郑国纲. 漳州核电厂一期工程厂址区域核应急方案编制介绍 [J]. 科技传播, 2016, 8(5): 168–169
    FANG Quncai, ZHENG Guogang. Introduction of regional nuclear emergency plan for Zhangzhou nuclear power plant phase I project site [J]. Public Communication of Science & Technology, 2016, 8(5): 168–169
    [15] 杨真真, 杨震, 孙林慧. 信号压缩重构的正交匹配追踪类算法综述 [J]. 信号处理, 2013, 29(4): 486–496 doi: 10.3969/j.issn.1003-0530.2013.04.011
    YANG Zhenzhen, YANG Zhen, SUN Linhui. A survey on orthogonal matching pursuit type algorithms for signal compression and reconstruction [J]. Journal of Signal Processing, 2013, 29(4): 486–496 doi: 10.3969/j.issn.1003-0530.2013.04.011
    [16] 孙强, 王梦晓, 徐玉山, 等. 岩巷爆破振动信号的HHT分析与应用[J]. 工程爆破, 2016, 22(1): 1-7.
    SUN Qiang, WANG Mengxiao, XU Yushan, et al. HHT analysis and application of blasting vibration in rock roadway excavation [J]. Engineering Blasting, 2016, 22 (1): 1-7.
  • [1] 禹海涛柳献袁勇王庆国程国荣吴晓峰 . 地下墙施工爆破对既有隧道的振动测试与分析. 爆炸与冲击, 2009, 29(6): 566-572. doi: 10.11883/1001-1455(2009)06-0566-07
    [2] 钟明寿周辉刘影龙源郭涛 . 基于改进匹配追踪算法的化爆地震波信号时频特征提取. 爆炸与冲击, 2017, 37(6): 931-938. doi: 10.11883/1001-1455(2017)06-0931-08
    [3] 李夕兵张义平刘志祥左宇军王卫华 . 爆破震动信号的小波分析与HHT变换. 爆炸与冲击, 2005, 25(6): 528-535. doi: 10.11883/1001-1455(2005)06-0528-08
    [4] 刘爱文俞言祥付长华吕红山陈鲲赵纪生周正华王伟 . 50t科学爆破振动的衰减特征及地形效应. 爆炸与冲击, 2010, 30(1): 21-26. doi: 10.11883/1001-1455(2010)01-0021-06
    [5] 李云鹏艾传志韩常领霍明 . 小间距隧道爆破开挖动力效应数值模拟研究. 爆炸与冲击, 2007, 27(1): 75-81. doi: 10.11883/1001-1455(2007)01-0075-07
    [6] 管永红胡八一黄超 . 基于小波包的爆炸容器振动分析. 爆炸与冲击, 2010, 30(5): 551-555. doi: 10.11883/1001-1455(2010)05-00551-05
    [7] 关晓磊颜景龙 . 爆破振动信号的HHT时频能量谱分析. 爆炸与冲击, 2012, 32(5): 535-541. doi: 10.11883/1001-1455(2012)05-0535-07
    [8] 林大超施惠基白春华张奇 . 爆炸地震效应的时频分析. 爆炸与冲击, 2003, 23(1): 31-36.
    [9] 马华原龙源李兴华钟明寿周辉吴建源 . 全尺寸干线输气管道爆炸地振动的空间分布和时频分布特性. 爆炸与冲击, 2019, 39(4): 042201-1-042201-9. doi: 10.11883/bzycj-2018-0123
    [10] 程擂韩焱王鉴杜娟 . 基于改进HHT的水中爆炸冲击波信号时频特性分析方法. 爆炸与冲击, 2011, 31(3): 326-331. doi: 10.11883/1001-1455(2011)03-0326-06
    [11] 晏俊伟龙源方向周春华 . 基于小波变换的爆破振动信号能量分布特征分析. 爆炸与冲击, 2007, 27(5): 405-410. doi: 10.11883/1001-1455(2007)05-0405-06
    [12] 凌同华李夕兵 . 地下工程爆破振动信号能量分布特征的小波包分析. 爆炸与冲击, 2004, 24(1): 63-68.
    [13] 周浩然 . 城市建设中的控制爆破实践. 爆炸与冲击, 1986, 6(2): 176-182.
    [14] 申永康邵建华陈建锋 . 建筑结构爆破地震反应弹塑性精细时程分析. 爆炸与冲击, 2008, 28(1): 92-96. doi: 10.11883/1001-1455(2008)01-0092-05
    [15] 栾龙发王文龙 . 台阶爆破破碎度的分布规律. 爆炸与冲击, 1988, 8(3): 273-280.
    [16] 夏祥李海波张大岩王新远李俊如 . 红沿河核电站基岩爆破的控制标准. 爆炸与冲击, 2010, 30(1): 27-32. doi: 10.11883/1001-1455(2010)01-0027-06
    [17] 杨国梁杨仁树姜琳 . 轴向间隔装药爆破沿炮孔的压力分布. 爆炸与冲击, 2012, 32(6): 653-657. doi: 10.11883/1001-1455(2012)06-0653-05
    [18] 李万张志华周峰张涛 . 水下目标在水下爆炸作用下冲击响应的时频特征. 爆炸与冲击, 2012, 32(3): 309-315. doi: 10.11883/1001-1455(2012)03-0309-07
    [19] 刘慧冯叔瑜 . 炸药单耗对爆破块度分布影响的理论探讨. 爆炸与冲击, 1997, 17(4): 359-362.
    [20] 张小军汪旭光王尹军于亚伦吴春平杨德强 . 基于正态分布的爆破振动评价与安全药量计算. 爆炸与冲击, 2018, 38(5): 1115-1120. doi: 10.11883/bzycj-2017-0091
  • 加载中
图(18)表(2)
计量
  • 文章访问数:  291
  • HTML全文浏览量:  748
  • PDF下载量:  32
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-02-09
  • 录用日期:  2018-03-19
  • 网络出版日期:  2019-03-25
  • 刊出日期:  2019-04-01

基于改进MP-WVD算法的核电厂建设爆破振动信号处理方法

    作者简介:梅 比(1981- ),男,博士,高级工程师,blastmei2017@163.com
  • 1. 中国矿业大学(北京)力学与建筑工程学院,北京 100083
  • 2. 北京矿冶研究总院,北京 100160

摘要: 针对爆破施工产生的振动问题,结合现有信号分析方法优缺点,引入MP-WVD组合算法解决信号分析过程中交叉项干扰问题,可提高时频分辨率,满足精确提取核电爆破振动响应信号时频分布特征的需求;将HHT算法引入MP算法数据的预处理,成功降低了MP算法运算复杂度,可为大数据量分析奠定基础。运用该算法对漳州核电一期工程场地平整土石方爆破工程振动监测信号进行分析,精确提取了爆破振动信号的时频特征。

English Abstract

  • 我国核电建设进入加速期,由于核电站规模较大,无法多机组同时建设,因此在后期机组施工时就要考虑到已建机组的安全稳定,尤其是爆破施工产生的振动影响[1]。研究表明,建(构)筑物在地震波作用下产生的动力响应特征与波的时频能量特性密切相关。目前时-频分析算法有各自优缺点:STFT算法(short-time Fourier transform)最初为傅里叶变换应用于时频分析提供了可能,但是郭涛等[2]通过对传统的STFT算法与小波算法进行对比分析得到了传统的STFT以三角函数作为基函数,无法准确与爆破振动这种非稳态随机过程相匹配的结论。马华原等[3]将小波变化方法运用于核电施工爆破振动分析取得了良好的效果,但也发现小波变换虽然引入了更加适配的小波基,但时频分辨率受测不准原理限制。李夕兵等[4]、魏新江等[5]改进了HHT算法(Hilbert-Huang transform)并将其运用于爆炸振动分析中,获得了具有比小波变换分辨率更高的时频分布图谱。但随着研究的深入,HHT变换提出的自适应基底导致算法边界效应较大,这一点目前仍是个难题。另一方面,周辉等[6]对传统的匹配追踪算法进行了改进,提供了消除经典WVD算法(Wigner-Ville distribution)交叉干扰项的思路,在此基础上,本文将改进型MP算法(matching pursuit algorithm)与传统WVD算法相结合,成功解决了交叉项干扰的问题,同时很好地控制了算法的计算复杂度。结合核电爆破信号实例,取得了较好的分析效果。

    • WVD算法最初是由Ville将其应用到信号的时-频分析领域的[7],针对一个时间序列$W(t,\omega )$根据特征函数方法推导出WVD表达式:

      $W(t,\omega ) = \frac{1}{{2\text{π} }}\int {x\left(t - \frac{1}{2}\tau \right)} \cdot x\left(t + \frac{1}{2}\tau \right){{\rm{e}}^{ - j\tau \omega }}{\rm{d}}\tau $

      式中,t为时间,ω为角频率。由式(1)可知,在信号WVD表达中,不存在任何形式的窗函数。因此,WVD分相对STFT、小波变换而言有着更高的时间和频率分辨率灵活性。

    • MP算法是由Mallat于1993年提出的一种信号分解算法[8]。其核心原理是将信号以字典原子为基,进行分解。过程如图1所示:选取与信号Xn匹配程度最高的原子Ψn,并求出投影值an和差值信号Xn+1。此时得到残差信号Xn+1,重复进行原子匹配过程,将其投影到与其最相近的原子Ψn+1上,得到差值信号Xn+2。以此类推,直到残余信号的能量小于设定的阈值为止。

      图  1  分解过程示意图

      Figure 1.  Schematic diagram of decomposition process

      MP算法是一种贪婪算法,该方法与统计学中使用的投影追踪算法和波形增益矢量量化有密切关系。信号的分解需要在超完备子波库中进行,其中,超完备意为信号的分解目标在信号所组成的空间中足够密集,这也就是最终无法以一组正交基进行描述。对于任意一个有限维 Hilbert 空间HD为此空间内的一个超完备词典,设信号为fH,长度为ND中的元素满足:

      ${{D}} = \{ {{ g}_{\rm{\gamma }} }:\gamma \in \Gamma \} \left\| {{{ g}_{\rm{\gamma }} }} \right\| = 1$

      式中:gγ为子波分解算子,Γ为伽马函数。

      MP算法就是将信号f垂直投影到子波库D的子波上。设${{g}}_{\rm{\gamma }} ^{(0)} \in {{D}}$,则f可以表示为:

      ${{f}} = \left. {\left\langle {{{f}},{{g}}_{\rm{\gamma }}^{(0)}} \right.} \right\rangle {{g}}_{\rm{\gamma }}^{(0)} + {{{R}}_{\rm{f}}}$

      式中:Rf的意义为原始信号f通过子波匹配${{g}}_{\rm{\gamma }} ^{(0)}$进行一次分解的残余[9]。为了达到最好的分解效果,就必须使残余信号尽可能小,因此就必然要使内积项$\left\langle {{{f}},{{g}}_{\rm{\gamma }} ^{(0)}} \right\rangle $取最大。很显然,${{g}}_{\rm{\gamma }} ^{(0)}$Rf是正交的,因此:

      ${\left\| {{f}} \right\|^2} = {\left| {\left\langle {\left. {{{f}},{{g}}_{\rm{\gamma }} ^{(0)}} \right\rangle } \right.} \right|^2}{{g}}_{\rm{\gamma }} ^{(0)} + {\left\| {{{{R}}_{\rm{f}}}} \right\|^2}$

      ${{R}}_{\rm{f}}^{(n)} = {{f}}$,经过了n次迭代(n>>0)得到残余信号${{R}}_{\rm{f}}^{(n)}$,此时再选择一个匹配子波${{g}}_{\rm{\gamma }} ^{(n)} \in {{D}}$,使其匹配${{R}}_{\rm{f}}^{(n)}$,即:

      ${{R}}_{\rm{f}}^{(n)} = \left\langle {\left. {{{R}}_{\rm{f}}^{(n)},{{g}}_{\rm{\gamma }} ^{(n)}} \right\rangle } \right.{{g}}_{\rm{\gamma }} ^{(n)} + {{R}}_{\rm{f}}^{(n + 1)}$

      ${{R}}_{\rm{f}}^{(n + 1)}$就是进行了n+1次迭代得到的差值。因此MP算法即是利用式(5)描述的一个重复迭代过程,若迭代m次,则f可以表示为:

      ${{f}} = \sum\limits_{n = 0}^{m - 1} {\left\langle {\left. {{{R}}_{\rm{f}}^{(n)},{{g}}_{\rm{\gamma }} ^{(n)}} \right\rangle } \right.{{g}}_{\rm{\gamma }} ^{(n)} + {{R}}_{\rm{f}}^{(m)}} $

    • WVD分布为非线性时频分布,也就是说其不具备线性分布的可加性,即两信号和的 WVD并不等于每一个信号的 WVD之和[10]

      x(t)=x1(t)+x2(t),则:

      $\begin{split} W(t,\omega ) =& \frac{1}{{2\text{π} }}\int {\left[{x_1}\left(t + \frac{1}{2}\tau \right) + {x_2}\left(t + \frac{1}{2}\tau \right)\right ]} \cdot \left[{x_1}\left(t - \frac{1}{2}\tau \right) + {x_2}\left(t - \frac{1}{2}\tau \right)\right]{{\rm{e}}^{ - j\tau \omega }}{\rm{d}}\tau =\\& {W_{{x_1}}}(t,\omega ) + {W_{{x_2}}}(t,\omega ) + 2{\mathop{\rm Re}\nolimits} [{W_{{x_1} + {x_2}}}(t,\omega )] \end{split}$

      式中:$2{\mathop{\rm Re}\nolimits} [{W_{{x_1} + {x_2}}}(t,\omega )]$x1(t)和x2(t)的交叉项。t为时间,ω为角频率。

      由式(7)可以看出,交叉项的存在对信号的时-频分布产生了很大的干扰。因此,将MP算法与WVD结合起来,将原始信号细分成基本原子,然后再对每个信号做WVD,将得到的结果叠加,在消除干扰项的基础上得到更清晰的时频分布图谱。具体示例如下。

      示例信号x(t),如图2所示。对示例信号x直接进行WVD运算,得到的时频分布情况如图3所示。可以看出图中不仅有4个真实分量的时频能量分布,还在任意2个分量的时频中心连线的中点处出现了交叉项,严重干扰了对信号真实时频分布的判读。

      图  3  WVD时频分布图

      Figure 3.  WVD time-frequency distribution

      图  2  示例信号x(t)

      Figure 2.  Example signal x(t)

      现将示例信号x(t)代入MP算法,分解为x1x2x3x4共4个子项,如图4所示。

      图  4  MP算法分解结果

      Figure 4.  Decomposition results of MP algorithm

      再分别对4个子波进行WVD运算,得到其各自的WVD时频分布结果,如图5所示;将各子波的时频分布合并,即得到了原示例信号的时频分布情况,如图6所示。

      图  5  各个子波的WVD分布

      Figure 5.  WVD distribution of each wavelet

      图  6  子波分解的WVD时频分布图

      Figure 6.  WVD time-frequency distribution of wavelet decomposition

      可见,WVD变换在任意两个有效信号之间产生一个交叉项。图3中信号有4个分量共产生6个交叉项(其中交叉项5和6重叠),对N个分量则会产生N×(N-1)/2个交叉项 。而图456中采用MP子波分解方法,成功剔除了交叉项干扰。

      另一方面,MP算法也存在着缺陷。其计算量过大,需耗费很长机时来完成一次信号分解[11]。而通常时频分析数据量都较大,因此研究者希望通过改进得到一种更高效的算法。

      算法选用Gabor子项的控制参量有:振幅、频率、中心时间和相位[12]。对过完备子波库的扫描过程利用穷举法对这上述参量进行优化选择。引入HHT算法,通过EMD(empirical mode decomposition)分解得到信号的固有模态IMF(intrinsic mode function)分量。对IMF分量进行Hilbert变换,得到信号的瞬时优势频率以及相位,并将其代入MP算法中,即可将4参数扫描运算降维至2参数扫描,大大减少了程序循环步数。程序流程如图7所示。

      图  7  程序框图

      Figure 7.  Program diagram

      与传统算法进行效率对比,分别用常规MP算法与改进后得MP算法对合成信号进行处理:

      表1中对比分析可以看出,经过改进后的算法效率有了明显提升,对同一信号分解所用机时有着50倍左右的差距。

      合成信号类型信号波形传统 MP 所用机时/s改进 MP 所用机时/s
      单正弦96.111.90
      双正弦77.131.90
      三正弦82.961.71
      典型爆破振动信号82.081.58
       注:信号长度,N=512,残差截断阈值,σ=0.05,CPU 型号,Intel Core i7-4710MQ,主频 2.5~3.1 GHz。

      表 1  两种方法所用机时对比

      Table 1.  Comparison of machine time between two methods

    • 结合福建漳州核电厂一期工程场地平整土石方工程施工实例中监测到的爆破振动信号进行分析。本工程施工区域以5#、6#机组及其厂区西侧为界,功能设施为5#机组与6#机组主生产区和厂前区[13-14]。开挖后场平标高,13.5 m;土方量,603 200m3;石方量,5 878 300m3;回填方量,3 573 300 m3;挖沟槽土方,648 m3;石方,1 513 m3;边坡预裂面积,39 600 m2

      爆破振动监测点设置3个测点:人家村村委会东侧民房设置测点1,距离爆破施工位置859 m;南山村地标设置测点2,距离爆破施工位置678 m;水坝闸门设置测点3,距离爆破施工位置291 m。本文中选取三标段爆破实例,深孔爆破,采用乳化90炸药,其爆破技术参数如表2所示。

      孔径/mm孔数孔深/m孔距/m排距/m抵抗线/m孔药量/kg段药量/kg总药量/kg
      1154713.05.53.22.285.0170.03 984

      表 2  爆破参数表

      Table 2.  Blasting parameters

    • 以水坝闸门测点的振动信号的铅垂分量(Z轴)数据为例,利用改进MP-WVD算法进行分析。本次爆破振动测试采用了TC-4850型爆破测振仪,采样频率为10 000 Hz,采集到的典型爆破振动信号如图8所示。第一步做MP分解,所得到的子波集合如图9所示。分解得到的子波均具有良好紧支性[15],可良好地呈现出信号的细节。

      图  8  原始信号

      Figure 8.  Original signal

      图  9  子波集合

      Figure 9.  Collection of wavelets

      数据处理过程中发现,MP算法用于爆破振动信号高频去燥取得了十分理想的效果,与其他滤波算法进行对比研究,结果如图10所示。

      图  10  滤波效果对比

      Figure 10.  Comparison of filtering effects

      图10可看出,MP算法得到的重构信号保留了几乎所有的原始信号振动细节,同时剔除了高频噪声。小波变换得到的重构信号也能保持与原始信号的高度吻合,但有较大噪声残留,信噪比不如MP算法。EMD算法由于采用了自适应基底[16],其算法速度非常快,但自适应基底导致滤波效果不稳定,容易将有用信息一同剔除;其边界效应也会引入多余信号,重构信号和原始信号差值较大。

      将所得原子矩阵代入WVD算法,分别计算每个原子的时频分布,再逐一叠加,得到信号的总体时频分布如图11所示,基于小波变换的时频分析结果如图12所示。

      图  11  MP-WVD算法结果

      Figure 11.  Result of MP-WVD algorithm

      图  12  小波变换结果

      Figure 12.  Result of wavelet transform

      图12可知,小波变换结果不仅频率分辨率较低,且时间分辨率也不如MP-WVD算法得到结果,且存在一定的边界溢出问题。MP-WVD算法不仅有着较高的频率分辨率,同时较好地反映了地震波的频率分布随时间的变化而发生改变。从图11中可以清晰地分辨出此振动信号的特征:振动峰值在0.13 s左右到达,频率中心为13 Hz。此时刻地震波频率分布最宽,在10~40 Hz的范围内都有分布。最后高频成分迅速衰减,只剩余13 Hz分量持续了较长时间。

      将该测点的水平径向分量(X轴)以及水平切向分量(Y轴)信号分别代入MP-WVD算法,得到其时频分布情况,分别如图1314所示。

      图  13  水平径向分量(X轴)

      Figure 13.  Horizontal radial component

      图  14  水平切向分量(Y轴)

      Figure 14.  Horizontal tangential component

      图1314中可以看出,水平切向振动信号的频率中心为21 Hz,水平径向振动信号在21 Hz和13 Hz处都有能量聚集。3个方向的振动信号均具有初始频带宽、高频衰减快的特点。不同点在于:水平切向的信号没有出现持续时间较长的13 Hz分量,且水平切向的能量分布明显比其他2个方向的更集中。

      对信号的时频分布函数分别进行时间积分和频率积分,得到信号的瞬时能量谱如图15所示,频率边际谱如图16所示。

      图  15  瞬时能量谱

      Figure 15.  Instantaneous energy spectrum

      图  16  频率边际谱

      Figure 16.  Marginal spectrum

      频率边际谱中纵轴E·f 1表示能量在频率轴的分布密度, 量纲为J/Hz。

      不同方向的传感器收集的地震波类型不同,其中瑞利波的质点运动方向对应X轴和Z轴,勒夫波的质点运动方向对应Y轴和X轴。边际谱纵坐标E·f −1表示能量在频率轴上的分布密度。由瞬时能量谱可以看出,水平切向(Y轴)的振动峰值率先到达,铅垂方向(Z轴)的振动峰值有较大的时延,而水平径向(X轴)的振动有多个峰值。4个方向的瞬时能量曲线反映出了在此次爆破引起的振动中,勒夫波率先到达测点,瑞利波相较勒夫波有约0.05 s的时延。

      频率边际谱反映出信号的各个频率成分在时间全局上的累加,与傅氏幅频谱意义不同的是,傅氏谱只能反映信号的频率成分存在的概率而不能描述每个分量在整个振动过程中的能量份额。由图16可以看出,水平切向和铅垂方向振动频率较为集中,其频率中心分别为13 Hz和21 Hz。水平径向的振动由于混杂了勒夫波和瑞利波两个成分,所以具有多个频率中心。尤其需注意的是在36 Hz左右,X轴信号的边际谱有较强峰值出现,但在时频谱上并未见较强的时频能量集中,表明能量在此频率处集中但在时间上并不集中,边际谱的峰值是由于全局时间累积产生的,这种峰值小但加载时间长的特殊能量加载形式也需高度重视,防止产生损伤累积效应。

      将南山村测点以及人家村测点所采集到的数据代入算法进行处理,得到的结果如图1718所示。

      图  17  南山村测点数据分析结果

      Figure 17.  Data analysis result of Nanshan village

      图  18  人家村测点数据分析结果

      Figure 18.  Data analysis result of Renjia village

      图17、18中可以看出,由于测试点距离起爆点较远,瑞利波在传递过程中的色散效应产生了到达时间差,导致振动能量在时域上的分布更分散。

      在频域上,能量分布更加向低频集中,绝大部分能量分布在10~20 Hz频带,20 Hz以上频带几乎没有分布。可以看出在爆破地震波中,高频成分衰减较快,而低频成分传播较远,爆破远区的抗震设计当以应对10~20 Hz频段的低频振动为主。

    • (1)将MP算法与WVD分布结合,有效地消除了交叉项的干扰,进而发挥出WVD分布对瞬态信号敏感以及高分辨率的特点。引入HHT算法对信号的瞬时频率以及瞬时相位先行确定,可有效降低计算复杂度,机时平均可缩短至之前的约2%。(2)改进型MP算法用于处理地震波数据有着很好的去燥效果,能在尽量保留信号局部特征的基础上剔除噪声信号,相比小波滤波以及EMD滤波性能较为突出。(3)MP-WVD算法对于核电爆破施工振动信号的分析结果较之小波变换方法具有更高的分辨率以及细节刻画能力。(4)由WVD时频谱衍生出的瞬时能量谱能够清晰地反映出信号的能量峰值到达时间以及加载次数。频率边际谱能够描述任意频率成分在时间全局上的累积,这些特性在研究结构动态响应过程中均可纳入参考。

参考文献 (16)

目录

    /

    返回文章
    返回