CN104459763A - 紧支集小波探测地下空洞位置的方法及系统 - Google Patents
紧支集小波探测地下空洞位置的方法及系统 Download PDFInfo
- Publication number
- CN104459763A CN104459763A CN201410737899.5A CN201410737899A CN104459763A CN 104459763 A CN104459763 A CN 104459763A CN 201410737899 A CN201410737899 A CN 201410737899A CN 104459763 A CN104459763 A CN 104459763A
- Authority
- CN
- China
- Prior art keywords
- phase velocity
- frequency division
- rayleigh waves
- underground
- underground cavity
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种紧支集小波探测地下空洞位置的方法,涉及地下空洞探测技术领域,该方法包括确定监测区域,在监测区域中确定监测点,在监测点的具体位置安装数字检波器,根据紧支集小波分析法对采集的瑞雷波地震数据进行分析得到时-频分布图,根据时-频分布图分别计算出瑞雷波相速度和深度,以相速度为横坐标,深度为纵坐标,绘制H-V关系图,H-V关系图中相速度发生突变的地方即为地下空洞的位置。本发明的有益效果:因为瑞雷波是时变信号,选用的紧支集小波能兼顾小波函数的紧支集和平滑性,能高精度的模仿和分析瑞雷波地震信号,通过紧支集小波分析法探测地下空洞位置能准确判断地下空洞的具体的位置,提高判断地下空洞位置的准确性。
Description
技术领域
本发明涉及地下空洞探测技术领域,特别涉及一种紧支集小波探测地下空洞位置的方法及系统。
背景技术
天然形成或人工挖凿形成的地下空洞,对城市基础建设具有潜在巨大隐患。弄清地下空洞、空穴的分布情况是进行工程施工前需要解决的重要问题之一。地震勘探以其高分辨率特点在空洞探测中相较于其他地球物理方法具有显著优势,其中瑞雷面波勘探技术已经成为该方面的研究热点,比如MASW、ReMi、AARW和Spiking Fileter Analysis在空洞、地下管道探测中的应用,究其原因是因为瑞利面波在近地表的能量要远大于其他弹性波。近地表存在不均匀体时,沿近地表传播的瑞雷面波会发生散射现象,散射的强度取决于不均匀体与周围介质的物性差异。现有的技术探测地下空洞位置根据瑞雷波的波长不同,其穿透深度相应不同的特性,频散曲线会出现异常跳动,可以判断地下空洞的位置。如果不是空洞,是其他物体时,也会存在异常跳动,因此,现有技术存在容易误判的问题。
发明内容
有鉴于此,本发明要解决的技术问题在于提供一种紧支集小波探测地下空洞位置的方法及系统,能提高探测地下空洞位置准确度。
本发明通过以下技术手段解决上述技术问题:
本发明提供了一种紧支集小波探测地下空洞位置的方法,包括以下步骤:
S1:确定监测区域,在监测区域中确定监测点,在监测点的具体位置安装数字检波器,所述数字检波器检测原始瑞雷波地震数据;
S2:对所述数字检波器的检测数据利用紧支集小波分析法获取时-频分布图,然后根据所述时-频分布图得到时间差,得到瑞雷波的相速度;
S3:根据所述瑞雷波的相速度计算得到瑞雷波从地面传到地下的深度,以所述相速度为横坐标,所述深度为纵坐标,绘制深度-相速度关系图;
S4:判断所述关系图中是否存在相速度发生突变,如果存在,相速度发生突变处为空洞位置。
进一步,所述S2根据所述探测数据利用紧支集小波分析法获取时-频分布图,再根据所述时-频分布图得到时间差,得到瑞雷波的相速度的具体步骤包括,
S201:根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小波法绘制时-频分布图;
S202:在所述时-频分布图上确定一个固定频率,在所述时-频分布图上找到所述频率能量极值点所对应的时间,计算出所述频率在相邻两个时-频分布图上的时间差t;
S203:根据V=x/t计算出所述频率的瑞雷波相速度。
进一步,所述S3中根据瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的具体步骤:将所述V代入H=V/(2f)计算出H,其中,H代表深度即瑞雷波从地面传到地下的距离,V表示相速度,以V为横坐标,H为纵坐标,绘制出H-V关系图。
本发明还包括一种紧支集小波探测地下空洞位置的系统,包括传感器组,用于采集监测点的地震波数据;
地震波数据处理模块,用于对所述传感器组采集的数据利用紧支集小波分析法获取时-频分布图,然后根据所述时-频分布图得到时间差,计算出瑞雷波的相速度,再计算出瑞雷波从地面传到地下的深度,以相速度为横坐标,深度为纵坐标绘制深度-相速度关系图,判断关系图中是否存在相速度发生突变的地方,如果存在,所述相速度发生突变处为地下空洞位置;
输出装置,用于输出地震波数据处理模块判断出地下空洞所在的位置。
进一步,所述传感器组采用多个数字检波器。
进一步,所述地震波数据处理模块计算出瑞雷波的相速度的具体步骤包括:
S201:根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小波法绘制时-频分布图;
S202:在所述时-频分布图上确定一个固定频率,在所述时-频分布图上找到所述频率能量极值点所对应的时间,计算出所述频率在相邻两个时-频分布图上的时间差t;
S203:根据V=x/t计算出所述频率的瑞雷波相速度。
进一步,所述地震波数据处理模块根据瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的具体步骤:将所述V代入H=V/(2f)计算出H,其中,H代表深度即瑞雷波从地面传到地下的距离,V表示相速度,以V为横坐标,H为纵坐标,绘制出H-V关系图。
本发明的有益效果:
本发明的紧支集小波探测地下空洞位置的方法,确定监测区域,在监测区域中确定监测点,在监测点的具体位置安装数字检波器,根据紧支集小波分析法对采集的瑞雷波地震数据进行分析得到时-频分布图,根据时-频分布图分别计算出瑞雷波相速度和深度,以相速度为横坐标,深度为纵坐标,绘制H-V关系图,H-V关系图中相速度发生突变的地方即为地下空洞的位置。因为瑞雷波是时变信号,选用的紧支集小波能兼顾小波函数的紧支集和平滑性,能高精度的模仿和分析瑞雷波地震信号,通过紧支集小波分析法探测地下空洞位置能准确判断地下空洞的具体的位置,提高判断地下空洞位置的准确性。
紧支集小波探测地下空洞位置的系统采用的硬件简单,只需在待检测点安装多个数字检波器,即可采集瑞雷波原始实测数据,根据实测数据得出波形图,再根据紧支集小波分析法得到时-频分布图,再计算出相速度和深度,绘制H-V关系图,判断H-V关系图中是否存在相速度发生突变,如果相速度发生突变,则可以确定相速度发生突变的位置为地下空洞的位置。
附图说明
下面结合附图和实施例对本发明作进一步描述。
图1为本发明实施例1的流程图;
图2为本发明的一个数字检波器根据实测的数据绘制的波形图;
图3为本发明的另一个数字检波器根据实测的数据绘制的波形图;
图4为本发明的H-V关系图;
图5为本发明实施例2的结构示意图。
具体实施方式
以下将结合附图对本发明进行详细说明,
实施例1
本发明包括一种紧支集小波探测地下空洞位置的方法,具体包括以下步骤:
S1:确定监测区域,在监测区域中确定监测点,在监测点的具体位置安装数字检波器,所述数字检波器检测原始瑞雷波地震数据;
S2:对所述数字检波器的检测数据利用紧支集小波分析法获取时-频分布图,然后根据所述时-频分布图得到时间差,得到瑞雷波的相速度;
S3:根据所述瑞雷波的相速度计算得到瑞雷波从地面传到地下的深度,以所述相速度为横坐标,所述深度为纵坐标,绘制深度-相速度关系图;
S4:判断所述关系图中是否存在相速度发生突变,如果存在,相速度发生突变处为空洞位置。
确定监测区域,在监测区域中确定监测点,在监测点的具体位置安装数字检波器,如图1所示,在激震源附近设置多个数字检波器,数字检波器之间的间距为x。根据间距为x的2个数字检波器实测瑞雷波地震数据绘制的波形图,如图2和图3所示。根据数字检波器采集的地震数据,选取间距为x的两道数字检波器的数据利用紧支集小波函数进行分析,分析后得到时-频分布图。
利用紧支集小波分析求瑞雷波相速度的过程:
设f∈L2(R),紧支集小波变换的定义为 式(1)中,b平移参数,a为尺度参数。式(1)我们可以看出,正是窗函数(2),式(2)中的尺度参数a和位置参数b的变化使紧支集小波变换具有了局域特性,它既是时间的函数又是频率的函数。式(1)即为信号f(t)的时-频分布函数,利用这个时-频分布函数,可以求得瑞雷波的相速度:
(1)对相距为x的N个检波器接收到的信号做多尺度紧支集小波分析根据不同尺度所对应的频率段,选取需要的尺度下的信号。对选定的信号再做多尺度的紧支集小波分析得到一系列时-频分布图;
(2)对于每一个固定频率f,分别在每个时-频分布图上寻找该频率能量极值点所对应的时间t1、t2、…tn,进而可以计算出该频率在相邻两个时-频分布图上的时间差t=tn-tn-1(0<n<N),这个时差就是该频率以其相速度在相邻两个检波器之间传播形成的时间差;
(3)相速度的计算利用公式
可计算出每一固定频率f在相邻两道传播的相速度V。
(4)频散曲线的绘制。根据第(3)步的计算结果,可以绘出瑞雷波相速度随频率的变化曲线。
(5)根据半波长解释法其中,H代表深度即波从地面传到地下的距离,V表示相速度。根据由式(3)、式(4)计算出的H,V,以横坐标为V,纵坐标为H,绘制出H与V的关系图。如图4所示,在地底下0.05~0.15m处出现相速度发生突变,说明该处存在空洞。因此,通过H-V即深度与相速度之间的关系图就可以确定地下空洞所在的位置。
由于瑞雷波是时变信号,选用的紧支集小波能兼顾小波函数的紧支集和平滑性,能高精度的模仿和分析信号,通过紧支集小波分析法探测地下空洞位置能准确判断地下空洞的具体的位置,提高判断地下空洞位置的准确性。
实施例2
如图5所示,本发明还包括一种紧支集小波变换探测地下空洞位置的系统,包括传感器组,用于采集监测点的地震波数据;地震波数据处理模块,用于对所述传感器组采集的数据利用紧支集小波分析法获取时-频分布图,然后根据所述时-频分布图得到时间差,计算出瑞雷波的相速度,再计算出瑞雷波从地面传到地下的深度,以相速度为横坐标,深度为纵坐标,绘制H-V关系图,判断H-V关系图中所示相速度是否发生突变,如果相速度发生突变,H-V关系图中示出相速度发生突变处为地下空洞位置;输出装置,用于输出地震波数据处理模块判断出地下空洞的位置。传感器组与地震波数据处理模块连接,地震波数据处理模块与输出装置连接。所述传感器组采用多个数字检波器,数字检波器之间的间距为x。
地震波数据处理模块计算出瑞雷波的相速度的具体步骤包括:
S201:根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小波法绘制时-频分布图;
S202:对于每一个固定频率f,分别在每个时-频分布图上寻找该频率能量极值点所对应的时间t1、t2、…tn,进而可以计算出该频率在相邻两个时-频分布图上的时间差t=tn-tn-1(0<n<N),这个时差就是该频率以其相速度在相邻两个数字检波器之间传播形成的时间差;
S203:根据V=x/t计算出所述频率的相速度可计算出每一固定频率f在相邻两道传播的相速度V。
地震波数据处理模块判断空洞位置的具体方法:将所述V代入H=V/(2f)计算出H,其中,H代表深度即波从地面传到地下的距离,V表示相速度,以V为横坐标,H为纵坐标,绘制出H-V关系图,判断H-V关系图中是否存在相速度发生突变的地方如果存在,则关系图中相速度发生突变的地方为地下空洞位置,地震波数据处理模块判断出地下空洞的具体位置。如图4所示,在地底下0.05~0.15m处出现相速度发生突变,说明该处存在空洞。
输出装置输出地震波数据处理模块处理后判断出空洞的位置,供人们参考。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.紧支集小波探测地下空洞位置的方法,其特征在于:包括以下步骤:
S1:确定监测区域,在监测区域中确定监测点,在监测点的具体位置安装数字检波器,所述数字检波器检测原始瑞雷波地震数据;
S2:对所述数字检波器的检测数据利用紧支集小波分析法获取时-频分布图,然后根据所述时-频分布图得到时间差,得到瑞雷波的相速度;
S3:根据所述瑞雷波的相速度计算得到瑞雷波从地面传到地下的深度,以所述相速度为横坐标,所述深度为纵坐标,绘制深度-相速度关系图;
S4:判断所述关系图中是否存在相速度发生突变,如果存在,相速度发生突变处为空洞位置。
2.如权利要求1所述的紧支集小波探测地下空洞位置的方法,其特征在于:所述S2根据所述探测数据利用紧支集小波分析法获取时-频分布图,再根据所述时-频分布图得到时间差,得到瑞雷波的相速度的具体步骤包括,
S201:根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小波法绘制时-频分布图;
S202:在所述时-频分布图上确定一个固定频率,在所述时-频分布图上找到所述频率能量极值点所对应的时间,计算出所述频率在相邻两个时-频分布图上的时间差t;
S203:根据V=x/t计算出所述频率的相速度。
3.如权利要求1所述的紧支集小波探测地下空洞位置的方法,其特征在于:所述S3中根据瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的具体步骤:将所述V代入H=V/(2f)计算出H,其中,H代表深度即瑞雷波从地面传到地下的距离,V表示相速度,以V为横坐标,H为纵坐标,绘制出H-V关系图。
4.紧支集小波探测地下空洞位置的系统,其特征在于:包括传感器组,用于采集监测点的地震波数据;
地震波数据处理模块,用于对所述传感器组采集的数据利用紧支集小波分析法获取时-频分布图,然后根据所述时-频分布图得到时间差,计算出瑞雷波的相速度,再计算出瑞雷波从地面传到地下的深度,以相速度为横坐标,深度为纵坐标绘制图形,判断所述绘制图形中是否存在相速度发生突变的地方,所述相速度发生突变处为地下空洞位置;
输出装置,用于输出地震波数据处理模块判断出地下空洞所在的位置。
5.如权利要求4所述的紧支集小波探测地下空洞位置的系统,其特征在于:所述传感器组采用多个数字检波器。
6.如权利要求4所述的紧支集小波探测地下空洞位置的系统,其特征在于:所述地震波数据处理模块计算出瑞雷波的相速度的具体步骤包括:
S201:根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小波分析法绘制时-频分布图;
S202:在所述时-频分布图上确定一个固定频率,在所述时-频分布图上找到所述频率能量极值点所对应的时间,计算出所述频率在相邻两个时-频分布图上的时间差t;
S203:根据V=x/t计算出所述频率的瑞雷波相速度。
7.如权利要求4所述的紧支集小波探测地下空洞位置的系统,其特征在于:所述根据瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的具体步骤:将所述V代入H=V/(2f)计算出H,其中,H代表深度即瑞雷波从地面传到地下的距离,V表示相速度,以V为横坐标,H为纵坐标,绘制出H-V关系图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410737899.5A CN104459763A (zh) | 2014-12-05 | 2014-12-05 | 紧支集小波探测地下空洞位置的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410737899.5A CN104459763A (zh) | 2014-12-05 | 2014-12-05 | 紧支集小波探测地下空洞位置的方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104459763A true CN104459763A (zh) | 2015-03-25 |
Family
ID=52906108
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410737899.5A Pending CN104459763A (zh) | 2014-12-05 | 2014-12-05 | 紧支集小波探测地下空洞位置的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104459763A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106131210A (zh) * | 2015-08-21 | 2016-11-16 | 朱荣华 | 一种采集地下环境检测数据的系统及方法 |
CN106646586A (zh) * | 2016-12-26 | 2017-05-10 | 重庆市第中学校 | 一种利用瑞利波检测结构体内部裂缝位置的方法 |
CN107092030A (zh) * | 2017-05-12 | 2017-08-25 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
CN108037187A (zh) * | 2017-12-05 | 2018-05-15 | 南京航空航天大学 | 基于瑞雷波的机场跑道脱空监测方法 |
CN109188504A (zh) * | 2018-08-08 | 2019-01-11 | 武汉市工程科学技术研究院 | 一种地下管廊探测方法 |
CN112379403A (zh) * | 2020-12-14 | 2021-02-19 | 北京华晖探测科技股份有限公司 | 一种地下采空区的探测方法及系统 |
CN113125559A (zh) * | 2021-04-16 | 2021-07-16 | 山东正元地球物理信息技术有限公司 | 一种基于微动信息识别管线位置的方法 |
CN116776645A (zh) * | 2023-08-21 | 2023-09-19 | 成都协致科技有限责任公司 | 基于小波分析的环境空气监测站布点方法及系统 |
-
2014
- 2014-12-05 CN CN201410737899.5A patent/CN104459763A/zh active Pending
Non-Patent Citations (3)
Title |
---|
何世聪: "瑞雷波频散曲线提取方法的研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 * |
崔竹刚 等: "瑞雷面波在地下空洞探测中的应用", 《科学技术与工程》 * |
章哲辉: "瑞雷面波勘探在空洞探测中的应用", 《湖南水利水电》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106131210A (zh) * | 2015-08-21 | 2016-11-16 | 朱荣华 | 一种采集地下环境检测数据的系统及方法 |
CN106646586A (zh) * | 2016-12-26 | 2017-05-10 | 重庆市第中学校 | 一种利用瑞利波检测结构体内部裂缝位置的方法 |
CN107092030B (zh) * | 2017-05-12 | 2019-04-12 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
CN107092030A (zh) * | 2017-05-12 | 2017-08-25 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
CN108037187A (zh) * | 2017-12-05 | 2018-05-15 | 南京航空航天大学 | 基于瑞雷波的机场跑道脱空监测方法 |
CN109188504B (zh) * | 2018-08-08 | 2020-04-17 | 武汉市工程科学技术研究院 | 一种地下管廊探测方法 |
CN109188504A (zh) * | 2018-08-08 | 2019-01-11 | 武汉市工程科学技术研究院 | 一种地下管廊探测方法 |
CN112379403A (zh) * | 2020-12-14 | 2021-02-19 | 北京华晖探测科技股份有限公司 | 一种地下采空区的探测方法及系统 |
CN112379403B (zh) * | 2020-12-14 | 2024-01-16 | 北京华晖探测科技股份有限公司 | 一种地下采空区的探测方法及系统 |
CN113125559A (zh) * | 2021-04-16 | 2021-07-16 | 山东正元地球物理信息技术有限公司 | 一种基于微动信息识别管线位置的方法 |
CN113125559B (zh) * | 2021-04-16 | 2022-07-26 | 山东正元地球物理信息技术有限公司 | 一种基于微动信息识别管线位置的方法 |
CN116776645A (zh) * | 2023-08-21 | 2023-09-19 | 成都协致科技有限责任公司 | 基于小波分析的环境空气监测站布点方法及系统 |
CN116776645B (zh) * | 2023-08-21 | 2023-10-27 | 成都协致科技有限责任公司 | 基于小波分析的环境空气监测站布点方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104459763A (zh) | 紧支集小波探测地下空洞位置的方法及系统 | |
CN102879805B (zh) | 一种基于钻孔与地面相结合的地震波空间探测方法 | |
CN102373923B (zh) | 一种储层识别方法 | |
CN102565855B (zh) | 油田压裂地面微地震数据处理方法 | |
KR20200014387A (ko) | 지하 구조물의 검출 | |
CN105510880A (zh) | 一种基于双差法的微地震震源定位方法 | |
CN106501848B (zh) | 一种隧道掘进过程中隐性断层超前物探方法 | |
CA2876743C (en) | Quality control for broadband sweeps | |
CN104267442A (zh) | 一种用于煤矿井下的瞬变电磁拟地震探测方法 | |
JPH07325158A (ja) | 垂直弾性波探層を準備する方法 | |
Bellino et al. | An automatic method for data processing of seismic data in tunneling | |
CN103513273B (zh) | 陆用四分量数字地震检波器 | |
CN105277982A (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN104570082A (zh) | 一种基于格林函数表征的全波形反演梯度算子的提取方法 | |
CN103913768A (zh) | 基于地震波资料对地表中浅层进行建模的方法及装置 | |
CN109188528B (zh) | 井间弹性波层析成像系统和方法 | |
CN102914797B (zh) | 一种获得地层各向异性系数的方法及装置 | |
CN104570090A (zh) | 全波形反演噪音滤波算子的提取及使用其噪音滤波的方法 | |
CN117590466A (zh) | 一种微小断层检测方法及断层检测装置 | |
Zhang et al. | Tunnel disturbance events monitoring and recognition with distributed acoustic sensing (DAS) | |
CN107741459A (zh) | 利用孔间声波层析成像技术探测岩性异常体的方法及装置 | |
CN202204937U (zh) | 一种地震勘探仪的检测传输模块 | |
CN109188504B (zh) | 一种地下管廊探测方法 | |
Porras et al. | A semblance-based microseismic event detector for DAS data | |
US20050162974A1 (en) | Resonance scattering seismic method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150325 |
|
RJ01 | Rejection of invention patent application after publication |