CN111948704A - 一种近区地运动测试信号在地下爆炸监测中的应用方法 - Google Patents

一种近区地运动测试信号在地下爆炸监测中的应用方法 Download PDF

Info

Publication number
CN111948704A
CN111948704A CN202010815655.XA CN202010815655A CN111948704A CN 111948704 A CN111948704 A CN 111948704A CN 202010815655 A CN202010815655 A CN 202010815655A CN 111948704 A CN111948704 A CN 111948704A
Authority
CN
China
Prior art keywords
signal
underground explosion
wavelet packet
ground motion
test signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010815655.XA
Other languages
English (en)
Inventor
吴文溢
王万鹏
熊益波
钟方平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
63653 Troops of PLA
Original Assignee
63653 Troops of PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 63653 Troops of PLA filed Critical 63653 Troops of PLA
Priority to CN202010815655.XA priority Critical patent/CN111948704A/zh
Publication of CN111948704A publication Critical patent/CN111948704A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/02Generating seismic energy
    • G01V1/104Generating seismic energy using explosive charges
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/307Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种近区地运动测试信号在地下爆炸监测中的应用方法。所述近区地运动测试信号在地下爆炸监测中的应用方法包括以下步骤:S1:基于小波包和傅里叶变换建立组合滤波模型;S2:建立基于近区地运动信号的地下爆炸效果的评估算法;S3:通过实验对所述评估算法进行分析和验证本发明提供的近区地运动测试信号在地下爆炸监测中的应用方法,提出基于小波包和傅里叶变换的组合滤波模型,设计了一种利用近区地运动信号来评估爆炸地震效应的算法,其处理结果与地震结果误差精度在5%以内;解决了直接采用地运动信号来评估爆炸地震效应的难题,提高了试验测试数据的利用率,为地下爆炸监测与防护评估提供了一种新的途径,具有重要的工程应用价值。

Description

一种近区地运动测试信号在地下爆炸监测中的应用方法
技术领域
本发明涉及爆炸监测技术技术领域,尤其涉及一种近区地运动测试信号在地下爆炸监测中的应用方法。
背景技术
地下爆炸能瞬间释放出极大的能量,并在爆室周围的岩土介质中产生一个向四周传播的冲击波,当其强度迅速衰减到不足以使硬岩介质产生塑性变形时,便以地震波的形式向外传播。由于大地对地震波高频成分的吸收,百公里外的地震P波卓越频率约在几十Hz以下,因此,在地震监测中,经常通过对比百公里外同一地震台站的地震波频谱比值来评价地下爆炸的地震效应。地运动通常指地下爆炸引起的地表面运动,试验近区地运动参数测量是爆炸试验测试工作的重要环节,利用多量程多频段传感器获取的大量地运动数据,包含了丰富的爆炸信息。然而,由于近区地运动信号噪声成分的干扰,难以直接适用于爆炸效果的评估,且试验获取的大量地运动数据缺少充分的挖掘,未能给工程决策提供更加丰富、有效的技术支持。
因此,有必要提供一种新的近区地运动测试信号在地下爆炸监测中的应用方法解决上述技术问题。
发明内容
本发明解决的技术问题是提供一种监测精度高的近区地运动测试信号在地下爆炸监测中的应用方法。
为解决上述技术问题,本发明提供的近区地运动测试信号在地下爆炸监测中的应用方法包括以下步骤:
S1:基于小波包和傅里叶变换建立组合滤波模型;
S2:建立基于近区地运动信号的地下爆炸效果的评估算法;
S3:通过实验对所述评估算法进行分析和验证。
优选的,所述步骤S1中基于小波包和傅里叶变换建立组合滤波模型包括以下步骤:
S101:小波包去噪:小波包分析能对信号中的低、高频成分均进行多层次划分,并自适应的匹配信号频谱,具有很好的时-频分辨率;
S102:建立组合滤波模型:首先采用小波包去噪方法滤去地运动信号中的噪声,再利用傅里叶分析提炼出相关有效频谱信息。
优选的,所述步骤S2中的所述评估算法包括以下步骤:
S201:信号转换系数计算;
S202:小波包去噪;
S203:FFT和频谱比值计算;
S204:异常值处理。
优选的,所述步骤S201中信号转换系数计算的具体步骤为:通过地运动传感器(速度、加速度)采集得到的电信号乘以一定的放大系数转化成速度或加速度信号,放大系数的计算与传感器的量程参数有关。
优选的,所述步骤S202中小波包去噪的具体步骤为:选用去噪效果良好的db3小波包基,采取四阶小波包分解,选择Shannon熵标准来计算最优分解树,并选定阈值进行软阈值去噪。
优选的,所述步骤S203中FFT和频谱比值计算的具体步骤为:对于同一测点的地运动信号进行傅里叶分析时,需要保持分析时长、信号起点的一致,并对齐波形的峰值;有效频段的范围与远区地震的频段范围保持一致(30Hz以内),最低有效频率根据采样定理,由采样时长决定。
优选的,所述步骤S204中异常值处理的具体步骤为:对各个测点获取的谱比值数据进行正态分布检验,结果基本服从正态分布,采用莱因达准则进行谱比值结果的异常值处理。
优选的,所述步骤S3中通过实验对所述评估算法进行分析和验证的具体步骤如下:
S301:在距爆心百米范围内布置N个地运动信号测点,分别安装加速度传感器和速度传感器。其中,在测量加速度信号时,各测点放置两个加速度计,用于测量竖直和水平方向的信号;同样,在相同的位置布设速度传感器;
S302:按照所述评估算法步骤,对不同爆炸方式试验获取的近区地运动信号进行处理,试验代号用K表示;
S303:进一步验证算法的效果,对近区地运动谱比值结果和地震频谱比值结果做了比较,并对比了时域峰值比值。
与相关技术相比较,本发明提供的近区地运动测试信号在地下爆炸监测中的应用方法具有如下有益效果:
本发明提供一种近区地运动测试信号在地下爆炸监测中的应用方法,提出基于小波包和傅里叶变换的组合滤波模型,设计了一种利用近区地运动信号来评估爆炸地震效应的算法,其处理结果与地震结果误差精度在5%以内;该方法解决了直接采用地运动信号来评估爆炸地震效应的难题,提高了试验测试数据的利用率,为地下爆炸监测与防护评估提供了一种新的途径,具有重要的工程应用价值。
附图说明
图1为本发明提供的近区地运动测试信号在地下爆炸监测中的应用方法的组合滤波模型原理框图;
图2为部分实验正态分布检验结果图;
图3为对比试验K7/K19的加速度谱比值结果图;
图4为对比试验K7/K19的速度谱比值结果图;
图5为比值结果图。
具体实施方式
下面结合附图和实施方式对本发明作进一步说明。
请结合参阅图1-图5,近区地运动测试信号在地下爆炸监测中的应用方法包括以下步骤:
S1:基于小波包和傅里叶变换建立组合滤波模型;
S2:建立基于近区地运动信号的地下爆炸效果的评估算法;
S3:通过实验对所述评估算法进行分析和验证。
所述步骤S1中基于小波包和傅里叶变换建立组合滤波模型包括以下步骤:
S101:小波包去噪:小波包分析能对信号中的低、高频成分均进行多层次划分,并自适应的匹配信号频谱,具有很好的时-频分辨率;小波包基的选择很大程度上影响了信号的去噪效果;
S102:建立组合滤波模型:傅里叶变换能实现信号的时域和频域的特征关联,具有良好的频域分析能力,但时-频分辨率差,对于傅里叶谱中的某一频率,并不能确定该频率发生的时刻;此外,针对地下爆炸近区地运动信号这类对低频精度要求较高的信号,由于噪声成分的干扰,容易造成有效成分频谱的失真,从而影响对爆炸地震效应的评估;考虑到上述问题,本文利用小波包去噪能同时处理高、低频信号的优点,提出一种组合滤波模型;首先采用小波包去噪方法滤去地运动信号中的噪声,再利用傅里叶分析提炼出相关有效频谱信息,模型原理框图如图1所示。
所述步骤S2中的地下爆炸释放的能量以地震波的形式向外传播出去,由于大地对高频地震波的吸收,百公里至几千公里外的地震信号有效频率集中在几赫兹至几十赫兹;通过监测对比百公里外同一地震台站的地震波频谱比值,可以有效的评估不同措施下地下爆炸试验的爆炸效果;而爆炸近区地运动信号噪声干扰严重,直接进行时域幅值的比较所得结果与地震频谱比值结果关联性较差;基于提出的组合滤波模型,利用近区地运动信号进行地下爆炸效果评估的算法可归纳如下:
S201:信号转换系数计算;
S202:小波包去噪;
S203:FFT和频谱比值计算;
S204:异常值处理。
所述步骤S201中信号转换系数计算的具体步骤为:通过地运动传感器(速度、加速度)采集得到的电信号乘以一定的放大系数转化成速度或加速度信号,放大系数的计算与传感器的量程参数有关。
所述步骤S202中小波包去噪的具体步骤为:选用去噪效果良好的db3小波包基,采取四阶小波包分解,选择Shannon熵标准来计算最优分解树,并选定阈值进行软阈值去噪。
所述步骤S203中FFT和频谱比值计算的具体步骤为:对于同一测点的地运动信号进行傅里叶分析时,需要保持分析时长、信号起点的一致,并对齐波形的峰值;有效频段的范围与远区地震的频段范围保持一致(30Hz以内),最低有效频率根据采样定理,由采样时长决定。
所述步骤S204中异常值处理的具体步骤为:对各个测点获取的谱比值数据进行正态分布检验,结果基本服从正态分布,采用莱因达准则进行谱比值结果的异常值处理,检验结果见图2。
所述步骤S3中通过实验对所述评估算法进行分析和验证的具体步骤如下:
S301:在距爆心百米范围内布置3个地运动信号测点,分别安装加速度传感器和速度传感器;其中,在测量加速度信号时,各测点放置两个加速度计,用于测量竖直和水平方向的信号;同样,在相同的位置布设磁电式速度传感器,频响:0.5~200Hz,测速范围:<24cm/s,灵敏度240mV/cm/s,传感器的量程和技术指标如表1所示;采用EPS-2-M6Q便携式数字地震仪测量更远地区地震信号;
Figure BDA0002632553290000051
表1各测点加速度传感器量程和技术指标
S302:按照所述评估算法步骤,对不同爆炸方式试验获取的近区地运动信号进行处理,试验代号用K表示;图3和图4分别表示一组典型对比试验K7/K19的加速度信号处理结果和速度信号处理结果。
S303:进一步验证算法的效果,对近区地运动谱比值结果和地震频谱比值结果做了比较,并对比了时域峰值比值,结果见表2所示。频段的选择参照算法中的步骤3,范围在30Hz以内;为更加直观的对比,对加速度和速度信号的比值结果做了平均处理,结果如图5所示。
Figure BDA0002632553290000061
表2部分试验加速度、速度及地震波形频谱比值(频段:~30Hz)
综合分析上述结果,可以得出:
1)在时域峰值比较法中,由于噪声的干扰,同一组对比试验的速度和加速度信号会出现比值结果相差较大的现象,并且和地震频谱比值结果关联性较差;
2)从图3和图4的结果来看,部分测点的地运动信号频谱幅值波动较大,但从平均值的结果来看,在低频段(30Hz以下)的波动相对平稳;
3)结合表2和图5中多组对比试验的结果可以看出,采用算法处理近区地运动信号所获得频谱比值结果与地震结果更为接近,误差精度在5%以内。
与相关技术相比较,本发明提供的近区地运动测试信号在地下爆炸监测中的应用方法具有如下有益效果:
本发明提供一种近区地运动测试信号在地下爆炸监测中的应用方法,提出基于小波包和傅里叶变换的组合滤波模型,设计了一种利用近区地运动信号来评估爆炸地震效应的算法,其处理结果与地震结果误差精度在5%以内;该方法解决了直接采用地运动信号来评估爆炸地震效应的难题,提高了试验测试数据的利用率,为地下爆炸监测与防护评估提供了一种新的途径,具有重要的工程应用价值。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其它相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (8)

1.一种近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,包括以下步骤:
S1:基于小波包和傅里叶变换建立组合滤波模型;
S2:建立基于近区地运动信号的地下爆炸效果的评估算法;
S3:通过实验对所述评估算法进行分析和验证。
2.根据权利要求1所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S1中基于小波包和傅里叶变换建立组合滤波模型包括以下步骤:
S101:小波包去噪:小波包分析能对信号中的低、高频成分均进行多层次划分,并自适应的匹配信号频谱,具有很好的时-频分辨率;
S102:建立组合滤波模型:首先采用小波包去噪方法滤去地运动信号中的噪声,再利用傅里叶分析提炼出相关有效频谱信息。
3.根据权利要求1所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S2中的所述评估算法包括以下步骤:
S201:信号转换系数计算;
S202:小波包去噪;
S203:FFT和频谱比值计算;
S204:异常值处理。
4.根据权利要求3所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S201中信号转换系数计算的具体步骤为:通过地运动传感器(速度、加速度)采集得到的电信号乘以一定的放大系数转化成速度或加速度信号,放大系数的计算与传感器的量程参数有关。
5.根据权利要求3所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S202中小波包去噪的具体步骤为:选用去噪效果良好的db3小波包基,采取四阶小波包分解,选择Shannon熵标准来计算最优分解树,并选定阈值进行软阈值去噪。
6.根据权利要求3所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S203中FFT和频谱比值计算的具体步骤为:对于同一测点的地运动信号进行傅里叶分析时,需要保持分析时长、信号起点的一致,并对齐波形的峰值;有效频段的范围与远区地震的频段范围保持一致(30Hz以内),最低有效频率根据采样定理,由采样时长决定。
7.根据权利要求3所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S204中异常值处理的具体步骤为:对各个测点获取的谱比值数据进行正态分布检验,结果基本服从正态分布,采用莱因达准则进行谱比值结果的异常值处理。
8.根据权利要求1所述的近区地运动测试信号在地下爆炸监测中的应用方法,其特征在于,所述步骤S3中通过实验对所述评估算法进行分析和验证的具体步骤如下:
S301:在距爆心百米范围内布置N个地运动信号测点,分别安装加速度传感器和速度传感器。其中,在测量加速度信号时,各测点放置两个加速度计,用于测量竖直和水平方向的信号;同样,在相同的位置布设速度传感器;
S302:按照所述评估算法步骤,对不同爆炸方式试验获取的近区地运动信号进行处理,试验代号用K表示;
S303:进一步验证算法的效果,对近区地运动谱比值结果和地震频谱比值结果做了比较,并对比了时域峰值比值。
CN202010815655.XA 2020-08-14 2020-08-14 一种近区地运动测试信号在地下爆炸监测中的应用方法 Pending CN111948704A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010815655.XA CN111948704A (zh) 2020-08-14 2020-08-14 一种近区地运动测试信号在地下爆炸监测中的应用方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010815655.XA CN111948704A (zh) 2020-08-14 2020-08-14 一种近区地运动测试信号在地下爆炸监测中的应用方法

Publications (1)

Publication Number Publication Date
CN111948704A true CN111948704A (zh) 2020-11-17

Family

ID=73343108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010815655.XA Pending CN111948704A (zh) 2020-08-14 2020-08-14 一种近区地运动测试信号在地下爆炸监测中的应用方法

Country Status (1)

Country Link
CN (1) CN111948704A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114002734A (zh) * 2021-11-02 2022-02-01 中国人民解放军63653部队 一种地运动数据处理方法、装置、存储介质和电子设备

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130046475A1 (en) * 2011-08-19 2013-02-21 Pei-Yang Lin System and method for on-site instant seismic analysis
US20140324356A1 (en) * 2013-04-30 2014-10-30 Republic Of Korea (National Disaster Management Institute) Apparatus for evaluating safety of building using earthquake acceleration measurement
CN105137498A (zh) * 2015-09-17 2015-12-09 鲁东大学 一种基于特征融合的地下目标探测识别系统及方法
CN107144879A (zh) * 2017-05-12 2017-09-08 重庆大学 一种基于自适应滤波与小波变换结合的地震波降噪方法
CN107515424A (zh) * 2017-07-26 2017-12-26 山东科技大学 一种基于vmd与小波包的微震信号降噪滤波方法
CN107515421A (zh) * 2017-08-15 2017-12-26 中国石油化工股份有限公司江汉油田分公司物探研究院 基于小波包变换的频谱成像方法
CN107991706A (zh) * 2017-11-24 2018-05-04 华北科技学院 基于小波包多重阈值和改进经验模态分解的煤层水力压裂微震信号联合降噪方法
CN109828318A (zh) * 2019-01-25 2019-05-31 吉林大学 一种基于变分模态分解的磁共振测深信号噪声滤除方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130046475A1 (en) * 2011-08-19 2013-02-21 Pei-Yang Lin System and method for on-site instant seismic analysis
US20140324356A1 (en) * 2013-04-30 2014-10-30 Republic Of Korea (National Disaster Management Institute) Apparatus for evaluating safety of building using earthquake acceleration measurement
CN105137498A (zh) * 2015-09-17 2015-12-09 鲁东大学 一种基于特征融合的地下目标探测识别系统及方法
CN107144879A (zh) * 2017-05-12 2017-09-08 重庆大学 一种基于自适应滤波与小波变换结合的地震波降噪方法
CN107515424A (zh) * 2017-07-26 2017-12-26 山东科技大学 一种基于vmd与小波包的微震信号降噪滤波方法
CN107515421A (zh) * 2017-08-15 2017-12-26 中国石油化工股份有限公司江汉油田分公司物探研究院 基于小波包变换的频谱成像方法
CN107991706A (zh) * 2017-11-24 2018-05-04 华北科技学院 基于小波包多重阈值和改进经验模态分解的煤层水力压裂微震信号联合降噪方法
CN109828318A (zh) * 2019-01-25 2019-05-31 吉林大学 一种基于变分模态分解的磁共振测深信号噪声滤除方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
LIU SHU CONG: "Seismic signal wavelet-packet denoising and fast spectrum analysis", 《APPLIED MECHANICS AND MATERIALS》 *
宋诗谦: "爆破中实测地运动数据分析及应用", 《测控技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114002734A (zh) * 2021-11-02 2022-02-01 中国人民解放军63653部队 一种地运动数据处理方法、装置、存储介质和电子设备

Similar Documents

Publication Publication Date Title
Chen et al. Main frequency band of blast vibration signal based on wavelet packet transform
CN107991706B (zh) 基于小波包多重阈值和改进经验模态分解的煤层水力压裂微震信号联合降噪方法
CN101201407B (zh) 相对无高频泄漏等效n点平滑谱模拟反褶积方法
CN108267784A (zh) 一种地震信号随机噪声压制处理方法
CN110487136B (zh) 一种基于频谱叠加的台阶爆破毫秒延时压制降振方法
CN105242309A (zh) 一种地震数据中规则干扰的压制方法及装置
Zuo et al. Noise suppression method of microseismic signal based on complementary ensemble empirical mode decomposition and wavelet packet threshold
CN104133248B (zh) 一种高保真声波干扰压制方法
CN111007569A (zh) 一种集成经验模态分解的样本熵阈值微地震信号降噪方法
CN110737023A (zh) 一种矿用微震监测信号的处理方法
CN111948704A (zh) 一种近区地运动测试信号在地下爆炸监测中的应用方法
CN114814939A (zh) 一种煤矿微震台网监测效果评估方法
CN103792574A (zh) 一种储层频变含气性检测方法
CN103217709B (zh) 一种提高地震数据信噪比和分辨率的面波衰减方法
US20160245938A1 (en) Method and device for increasing frequency of seismic digital signal
CN103743969A (zh) 一种近场电磁辐射干扰测量信号的测量噪声抑制方法
CN112183407B (zh) 一种基于时频域谱减法的隧道地震波数据去噪方法及系统
CN103675917B (zh) 一种识别微地震监测信号谐振干扰并进行有效压制的方法
Barnaba et al. Site effect of the strong-motion site at Tolmezzo-Ambiesta Dam in northeastern Italy
CN104459772B (zh) 一种地震数字信号的提频方法及装置
Jiang et al. An improved method of local mean decomposition with adaptive noise and its application to microseismic signal processing in rock engineering
CN101285892A (zh) 宽频地震勘探方法
CN102890288B (zh) 一种地震波层速度反演方法
Parolai Determination of dispersive phase velocities by complex seismic trace analysis of surface waves (CASW)
CN109270575A (zh) 一种基于建筑物地震响应等效的爆破地震波模型构造方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20201117