紧支集小波探测地下空洞位置的方法及系统的制作方法
【专利摘要】本发明公开了一种紧支集小波探测地下空洞位置的方法，涉及地下空洞探测【技术领域】，该方法包括确定监测区域，在监测区域中确定监测点，在监测点的具体位置安装数字检波器，根据紧支集小波分析法对采集的瑞雷波地震数据进行分析得到时-频分布图，根据时-频分布图分别计算出瑞雷波相速度和深度，以相速度为横坐标，深度为纵坐标，绘制H-V关系图，H-V关系图中相速度发生突变的地方即为地下空洞的位置。本发明的有益效果：因为瑞雷波是时变信号，选用的紧支集小波能兼顾小波函数的紧支集和平滑性，能高精度的模仿和分析瑞雷波地震信号，通过紧支集小波分析法探测地下空洞位置能准确判断地下空洞的具体的位置，提高判断地下空洞位置的准确性。
【专利说明】紧支集小波探测地下空洞位置的方法及系统
【技术领域】
[0001] 本发明涉及地下空洞探测【技术领域】，特别涉及一种紧支集小波探测地下空洞位置 的方法及系统。
【背景技术】
[0002] 天然形成或人工挖凿形成的地下空洞，对城市基础建设具有潜在巨大隐患。弄清 地下空洞、空穴的分布情况是进行工程施工前需要解决的重要问题之一。地震勘探以其高 分辨率特点在空洞探测中相较于其他地球物理方法具有显著优势，其中瑞雷面波勘探技术 已经成为该方面的研究热点，比如MASW、ReMi、AARW和Spiking Fileter Analysis在空洞、 地下管道探测中的应用，究其原因是因为瑞利面波在近地表的能量要远大于其他弹性波。 近地表存在不均匀体时，沿近地表传播的瑞雷面波会发生散射现象，散射的强度取决于不 均匀体与周围介质的物性差异。现有的技术探测地下空洞位置根据瑞雷波的波长不同，其 穿透深度相应不同的特性，频散曲线会出现异常跳动，可以判断地下空洞的位置。如果不是 空洞，是其他物体时，也会存在异常跳动，因此，现有技术存在容易误判的问题。
【发明内容】
[0003] 有鉴于此，本发明要解决的技术问题在于提供一种紧支集小波探测地下空洞位置 的方法及系统，能提高探测地下空洞位置准确度。
[0004] 本发明通过以下技术手段解决上述技术问题：
[0005] 本发明提供了一种紧支集小波探测地下空洞位置的方法，包括以下步骤：
[0006] Sl :确定监测区域，在监测区域中确定监测点，在监测点的具体位置安装数字检波 器，所述数字检波器检测原始瑞雷波地震数据；
[0007] S2 :对所述数字检波器的检测数据利用紧支集小波分析法获取时-频分布图，然 后根据所述时-频分布图得到时间差，得到瑞雷波的相速度；
[0008] S3 :根据所述瑞雷波的相速度计算得到瑞雷波从地面传到地下的深度，以所述相 速度为横坐标，所述深度为纵坐标，绘制深度-相速度关系图；
[0009] S4 :判断所述关系图中是否存在相速度发生突变，如果存在，相速度发生突变处为 空洞位置。
[0010] 进一步，所述S2根据所述探测数据利用紧支集小波分析法获取时-频分布图，再 根据所述时-频分布图得到时间差，得到瑞雷波的相速度的具体步骤包括，
[0011] S201 :根据所述相距为X的相邻数字检波器的两道地震波数据记录分别利用紧支 集小波法绘制时-频分布图；
[0012] S202 :在所述时-频分布图上确定一个固定频率，在所述时-频分布图上找到所述 频率能量极值点所对应的时间，计算出所述频率在相邻两个时-频分布图上的时间差t ;
[0013] S203 :根据V = x/t计算出所述频率的瑞雷波相速度。
[0014] 进一步，所述S3中根据瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的 具体步骤：将所述V代入H = VA2f)计算出H，其中，H代表深度即瑞雷波从地面传到地下 的距离，V表示相速度，以V为横坐标，H为纵坐标，绘制出H-V关系图。
[0015] 本发明还包括一种紧支集小波探测地下空洞位置的系统，包括传感器组，用于采 集监测点的地震波数据；
[0016] 地震波数据处理模块，用于对所述传感器组采集的数据利用紧支集小波分析法获 取时-频分布图，然后根据所述时-频分布图得到时间差，计算出瑞雷波的相速度，再计算 出瑞雷波从地面传到地下的深度，以相速度为横坐标，深度为纵坐标绘制深度-相速度关 系图，判断关系图中是否存在相速度发生突变的地方，如果存在，所述相速度发生突变处为 地下空洞位置；
[0017] 输出装置，用于输出地震波数据处理模块判断出地下空洞所在的位置。
[0018] 进一步，所述传感器组采用多个数字检波器。
[0019] 进一步，所述地震波数据处理模块计算出瑞雷波的相速度的具体步骤包括：
[0020] S201 :根据所述相距为X的相邻数字检波器的两道地震波数据记录分别利用紧支 集小波法绘制时-频分布图；
[0021] S202 :在所述时-频分布图上确定一个固定频率，在所述时-频分布图上找到所述 频率能量极值点所对应的时间，计算出所述频率在相邻两个时-频分布图上的时间差t ;
[0022] S203 :根据V = x/t计算出所述频率的瑞雷波相速度。
[0023] 进一步，所述地震波数据处理模块根据瑞雷波相速度计算得到瑞雷波从地面传到 地下的深度的具体步骤：将所述V代入H = VA2f)计算出H，其中，H代表深度即瑞雷波从 地面传到地下的距离，V表示相速度，以V为横坐标，H为纵坐标，绘制出H-V关系图。
[0024] 本发明的有益效果：
[0025] 本发明的紧支集小波探测地下空洞位置的方法，确定监测区域，在监测区域中确 定监测点，在监测点的具体位置安装数字检波器，根据紧支集小波分析法对采集的瑞雷波 地震数据进行分析得到时-频分布图，根据时-频分布图分别计算出瑞雷波相速度和深度， 以相速度为横坐标，深度为纵坐标，绘制H-V关系图，H-V关系图中相速度发生突变的地方 即为地下空洞的位置。因为瑞雷波是时变信号，选用的紧支集小波能兼顾小波函数的紧支 集和平滑性，能高精度的模仿和分析瑞雷波地震信号，通过紧支集小波分析法探测地下空 洞位置能准确判断地下空洞的具体的位置，提高判断地下空洞位置的准确性。
[0026] 紧支集小波探测地下空洞位置的系统采用的硬件简单，只需在待检测点安装多个 数字检波器，即可采集瑞雷波原始实测数据，根据实测数据得出波形图，再根据紧支集小波 分析法得到时-频分布图，再计算出相速度和深度，绘制H-V关系图，判断H-V关系图中是 否存在相速度发生突变，如果相速度发生突变，则可以确定相速度发生突变的位置为地下 空洞的位置。
【专利附图】
【附图说明】
[0027] 下面结合附图和实施例对本发明作进一步描述。
[0028] 图1为本发明实施例1的流程图；
[0029] 图2为本发明的一个数字检波器根据实测的数据绘制的波形图；
[0030] 图3为本发明的另一个数字检波器根据实测的数据绘制的波形图；
[0031] 图4为本发明的H-V关系图；
[0032] 图5为本发明实施例2的结构示意图。
【具体实施方式】
[0033] 以下将结合附图对本发明进行详细说明，
[0034] 实施例1
[0035] 本发明包括一种紧支集小波探测地下空洞位置的方法，具体包括以下步骤：
[0036] Sl :确定监测区域，在监测区域中确定监测点，在监测点的具体位置安装数字检波 器，所述数字检波器检测原始瑞雷波地震数据；
[0037] S2 :对所述数字检波器的检测数据利用紧支集小波分析法获取时-频分布图，然 后根据所述时-频分布图得到时间差，得到瑞雷波的相速度；
[0038] S3 :根据所述瑞雷波的相速度计算得到瑞雷波从地面传到地下的深度，以所述相 速度为横坐标，所述深度为纵坐标，绘制深度-相速度关系图；
[0039] S4 :判断所述关系图中是否存在相速度发生突变，如果存在，相速度发生突变处为 空洞位置。
[0040] 确定监测区域，在监测区域中确定监测点，在监测点的具体位置安装数字检波器， 如图1所示，在激震源附近设置多个数字检波器，数字检波器之间的间距为X。根据间距为 X的2个数字检波器实测瑞雷波地震数据绘制的波形图，如图2和图3所示。根据数字检 波器采集的地震数据，选取间距为X的两道数字检波器的数据利用紧支集小波函数进行分 析，分析后得到时-频分布图。
[0041] 利用紧支集小波分析求瑞雷波相速度的过程：
[0042] 设f e L2 (R)，紧支集小波变换的定义为
【权利要求】
1. 紧支集小波探测地下空洞位置的方法，其特征在于：包括以下步骤： 51 :确定监测区域，在监测区域中确定监测点，在监测点的具体位置安装数字检波器， 所述数字检波器检测原始瑞雷波地震数据； 52 :对所述数字检波器的检测数据利用紧支集小波分析法获取时-频分布图，然后根 据所述时-频分布图得到时间差，得到瑞雷波的相速度； 53 :根据所述瑞雷波的相速度计算得到瑞雷波从地面传到地下的深度，以所述相速度 为横坐标，所述深度为纵坐标，绘制深度-相速度关系图； 54 :判断所述关系图中是否存在相速度发生突变，如果存在，相速度发生突变处为空洞 位置。
2. 如权利要求1所述的紧支集小波探测地下空洞位置的方法，其特征在于：所述S2根 据所述探测数据利用紧支集小波分析法获取时-频分布图，再根据所述时-频分布图得到 时间差，得到瑞雷波的相速度的具体步骤包括， 5201 :根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小 波法绘制时-频分布图； 5202 :在所述时-频分布图上确定一个固定频率，在所述时-频分布图上找到所述频率 能量极值点所对应的时间，计算出所述频率在相邻两个时-频分布图上的时间差t ; 5203 :根据V = x/t计算出所述频率的相速度。
3. 如权利要求1所述的紧支集小波探测地下空洞位置的方法，其特征在于：所述S3中 根据瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的具体步骤：将所述V代入H = VA2f)计算出H，其中，H代表深度即瑞雷波从地面传到地下的距离，V表示相速度，以V为 横坐标，H为纵坐标，绘制出H-V关系图。
4. 紧支集小波探测地下空洞位置的系统，其特征在于：包括传感器组，用于采集监测 点的地震波数据； 地震波数据处理模块，用于对所述传感器组采集的数据利用紧支集小波分析法获取 时-频分布图，然后根据所述时-频分布图得到时间差，计算出瑞雷波的相速度，再计算出 瑞雷波从地面传到地下的深度，以相速度为横坐标，深度为纵坐标绘制图形，判断所述绘制 图形中是否存在相速度发生突变的地方，所述相速度发生突变处为地下空洞位置； 输出装置，用于输出地震波数据处理模块判断出地下空洞所在的位置。
5. 如权利要求4所述的紧支集小波探测地下空洞位置的系统，其特征在于：所述传感 器组采用多个数字检波器。
6. 如权利要求4所述的紧支集小波探测地下空洞位置的系统，其特征在于：所述地震 波数据处理模块计算出瑞雷波的相速度的具体步骤包括： 5201 :根据所述相距为x的相邻数字检波器的两道地震波数据记录分别利用紧支集小 波分析法绘制时-频分布图； 5202 :在所述时-频分布图上确定一个固定频率，在所述时-频分布图上找到所述频率 能量极值点所对应的时间，计算出所述频率在相邻两个时-频分布图上的时间差t ; 5203 :根据V = x/t计算出所述频率的瑞雷波相速度。
7. 如权利要求4所述的紧支集小波探测地下空洞位置的系统，其特征在于：所述根据 瑞雷波相速度计算得到瑞雷波从地面传到地下的深度的具体步骤：将所述V代入H = V/ (2f)计算出H，其中，H代表深度即瑞雷波从地面传到地下的距离，V表示相速度，以V为横 坐标，H为纵坐标,绘制出H-V关系图。
【文档编号】G01V1/28GK104459763SQ201410737899
【公开日】2015年3月25日 申请日期:2014年12月5日 优先权日:2014年12月5日
【发明者】柴贺军, 甘平, 孟云伟, 黄河, 李东链, 申晓青, 陈卓, 张宇, 李真真, 余欢, 曾彩兰 申请人:招商局重庆交通科研设计院有限公司, 重庆大学