CN102426384B - 一种探测地下采空区和岩溶分布的方法 - Google Patents

一种探测地下采空区和岩溶分布的方法 Download PDF

Info

Publication number
CN102426384B
CN102426384B CN201110262208.7A CN201110262208A CN102426384B CN 102426384 B CN102426384 B CN 102426384B CN 201110262208 A CN201110262208 A CN 201110262208A CN 102426384 B CN102426384 B CN 102426384B
Authority
CN
China
Prior art keywords
wave
karst
goaf
velocity
underground
Prior art date
Application number
CN201110262208.7A
Other languages
English (en)
Other versions
CN102426384A (zh
Inventor
赵永贵
蒋辉
赵晓鹏
Original Assignee
赵永贵
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 赵永贵 filed Critical 赵永贵
Priority to CN201110262208.7A priority Critical patent/CN102426384B/zh
Publication of CN102426384A publication Critical patent/CN102426384A/zh
Application granted granted Critical
Publication of CN102426384B publication Critical patent/CN102426384B/zh

Links

Abstract

本发明提供一种探测地下采空区和岩溶分布的方法,步骤如下:步骤1,在待探测区域布设检波器和震源;步骤2,震源激发产生地震波,检波器采集地下震动信息;步骤3,进行空间方向滤波,滤除杂波,得到地下采空区和岩溶地质体的散射波和反射波,计算其能量分布谱;步骤4,对散射波和反射波进行速度扫描,确定岩体波速分布;步骤5,制定地下采空区和岩溶的形态与分布特征图。本发明提供了一种用于探测地下采空区及岩溶分布的方法,克服现有技术中不能有效探查地形和地质结构复杂背景和无水高阻采空区及岩溶的缺陷。本发明综合利用散射波和反射波,应用合成孔径偏移成像技术实现对采空区和岩溶分布的有效探测,为工程施工的安全进行提供有效保证。

Description

一种探测地下采空区和岩溶分布的方法
技术领域
[0001] 本发明涉及工程地质勘查领域,尤其涉及一种探测地下采空区和岩溶地质灾害诊断技术。
背景技术
[0002] 自20世纪末以来,我国矿业开采秩序较为混乱,非法无序的乱采滥挖在一些矿山及其周边留下了大量的采空区,随着时间的推移,地下水和瓦斯气体充填采空区,一旦挖通就会造成透水、瓦斯爆炸等矿难。由于地下采空区具有隐伏性强、空间分布特征规律性差、采空区顶板冒落塌陷情况难以预测等特点,因煤矿采空区引起的煤矿灾害屡屡发生,著名的王家岭矿“3.28”特别重大透水事故就是由于挖通采空区导致巷道被淹。因此,如何勘查地下采空区的分布范围、空间形态特征,一直是工程技术领域的难题。
[0003] 岩溶形态复杂、界面极不规则,对地震波难以形成有效反射,主要以散射形式为主,信号较弱、干扰大、难以探测。在铁路、公路建设中对施工安全和施工质量有极大地威胁。
[0004] 当前,地下采空区和岩溶的常用探测方法主要有三维地震法、高密度电法和瞬变电磁法。
[0005] 其中,三维地震法是建立在反射地震理论基础之上,主要研究水平连续地层,不适合复杂、地质结构横向变化大的山区。三维地震法在山区探测的效果不佳,难于获得清楚的地质图像,更难于判断采空区的位置。高密度电法和瞬变电磁对于有填充的低阻采空区反应灵敏,但是对于未填充的高阻采空区极不敏感,这种采空区往往被漏报。这两种方法都不能有效地解决中小型采空区和岩溶的探测问题,特别是对于未填充的采空区和岩溶的探测效果不佳。
[0006] 本发明建立在散射理论基础上,不但适合大的连续界面,更适合于横向剧烈变化的场合;散射波可以有效的发现孤立的溶洞和采空区。,不管是否有填充。
发明内容
[0007] 本发明所要解决的技术问题是提供一种能准确探测地下采空区和岩溶分布的方法。
[0008] 本发明探测地下采空区和岩溶分布的方法,具体步骤如下:
[0009] 步骤1,在待探测区域内依次间隔地布设多个地震检波器和多个机械震源;且,地震检波器之间的距离不大于所述机械震源产生的地震波波长的1/4 ;布设的多个地震检波器和多个机械震源形成的总排列长度不小于所述地震波的波长;
[0010] 或者,在待探测区域内依次布设多个地震检波器后再布设一个机械震源;或者,在待探测区域内布设一个机械震源后再依次布设多个地震检波器;
[0011] 其中,地震检波器之间的距离不大于所述机械震源产生的地震波波长的1/4 ;且,布设的多个地震检波器形成的总排列长度不小于所述地震波的波长;或地下散射波与反射波的偏移叠加能量按1,当叠加能量取极值时对应的波速就是最
I:
3得到的地下散射波和反射波的能量、以及觉孔径偏移成像,制定所述地下采空区和岩
艮空区及岩溶分布的方法,以克服现有技术7能有效探查采空区及岩溶的缺陷。本发明,更适合于横向剧烈变化的场合;散射波可[0036] 图7为速度扫描结果的示意图;
[0037] 图8为合成孔径偏移成像的结果示意图。
具体实施方式
[0038] 以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
[0039] 为了提高地下采空区和岩溶分布的探测准确度,本发明提供一种综合考虑地质体散射波和反射波进行地层与地质病害(这里以采空区和岩溶为典型代表)探测的方法。如图1所示,具体步骤如下:
[0040] 步骤1,在待探测区域地表上布设地震检波器10和机械震源20,形成检测系统。其中,地震检波器10还要与主控机相通信,将采集的信息发送至主控机进行数字化和数据分析处理。
[0041] 步骤2,机械震源20激发产生地震波,地震检波器10采集包含地下反射与散射信号的地下震动信息,并发送至主控机。
[0042] 步骤3,主控机接收地震检波器10发送来的地下震动信息,对地下震动信息进行空间方向滤波,滤除其中的杂波干扰,得到来自地下采空区和岩溶的散射波和反射波,并进一步通过散射波和反射波的幅值变化计算获得散射波和反射波的能量分布谱。
[0043] 步骤4,对散射波和反射波进行速度扫描,确定地下岩体波速分布。
[0044] 步骤5,使用滤波后的散射波和反射波、扫描得到的速度分布进行合成孔径偏移成像,来确定采空区和岩溶的形态与分布特征图。
[0045] 如图2所示,为了提高主控机处理地下震动信息的速度及精准度,在进行步骤3方向滤波之前,先对地震检波器10采集的地下震动信息进行数字化和预处理。主要包括将地震检波器10发送来的电信号转变为数字信号,然后再进行去除地震检波器10采集的噪声干扰信号、为提高多个地震检波器10采集信息的同步性而进行的延时校正、为更好的进行信号识别而进行的数据增益放大等预处理操作中的一种或多种。
[0046] 下面对本发明的探测地下采空区和岩溶分布的具体过程进行简单介绍。
[0047] 首先,在待探测区内布设检测系统,主要包括用于产生地震波的机械震源20,用于采集地下震动信息的地震检波器10,用于进行信号数字化和数据处理的主控机。其中主控机与地震检波器10之间可以通过线缆相连进行有线通信,二者也可以通过无线方式连接进行无线通信。本发明中针对主控机带载能力的不同,可以按照以下两种方式在待探测区内布设机械震源20和地震检波器10:
[0048] 1.若主控机通道数量小(24道以内),则在待探测区域内依次间隔地布设多个地震检波器10和多个机械震源20,也就是说每间隔一个地震检波器10就布置一个机械震源20,如图3所示。若机械震源20激发的地震波的波长为λ,则地震检波器10之间的距离d应满足1/4 λ ;而多个地震检波器10排列形成的长度1,即布设多个地震检波器10和多个机械震源20形成的总排列长度I应满足I > λ。进一步地,为了达到最佳的探测结果,多个地震检波器10和多个机械震源20应布设在同一水平线上。在进行现场探测时,布设好检测系统后,可保持地震检波器10和机械震源20的相对位置不变,沿方向a同步向前移动进行探测。[0049] 2.若主控机通道数量大(32道以上),则在待探测区域内依次布设多个地震检波器10后再布设一个机械震源20 ;或者,在待探测区域内布设一个机械震源20后再依次布设多个地震检波器10。也就是说间隔布设多个地震检波器10后,在首尾地震检波器10位置上再布设一个机械震源20,形成检测系统,如图4所示。若机械震源20激发的地震波的波长为λ,则地震检波器10之间的距离cT应满足cT ^ 1/4 λ ;而布设的多个地震检波器ίο形成的总排列长度I y应满足I y ^ λ。同样地,为了达到最佳的探测结果,多个地震检波器10和一个机械震源20布设在同一水平线上。在进行现场探测时,布设好检测系统后,可保持地震检波器10和机械震源20的相对位置不变,沿方向a同步向前移动进行探测。
[0050] 其次,机械震源20激发产生地震波,地震检波器10为一组传感器,可以将地下的震动转变为电信号发送至主控机。主控机接收到地震检波器10发送的地下震动信息,先将地下震动信息进行数字化,再对其进行预处理,去除采集到的检测现场的干扰噪声信号,对信息进行同步和增益放大后,对初步处理之后生成的地震记录信息进行保存。然后把记录投影在F-K域里,如图5所示,不同方向来的回波会分布在不同的F-K域区间,通过选择合适的滤波因子,保留从地下来的反射波和散射波,滤除表面波、直达波、侧向回波等其他方向的回波等干扰波。经过方向滤波后,就得到了用于表征地下形态特征的散射波和反射波能量分布谱,如图6所示。至此就得到了用于探测地下采空区和岩溶形态特征图的第一方面数据。
[0051] 接着,主控机对散射波和反射波进行速度扫描,确定不同地段不同深度的围岩波速分布情况。
[0052] 对待探测区域地下散射波与反射波的偏移叠加能量按照以下公式求极值。递进改变波速计算叠加能量,当叠加能量取极值时对应的波速就是最佳波速。按照以下公式计算不同部位的最佳波速:
Figure CN102426384BD00061
[0054] 其中,V为介质波速,为待求量,它反应介质的力学性状。h和X分别代表垂直方向和水平方向的距离,其中,h为距离地面的深度;x为水平方向的坐标。
[0055] 当波速值取真实值时积分曲线与记录到的走时曲线一致,积分取最大值。因而真实的波速可由像点值的极值求取。
[0056] 上式的计算可以通过分区速度扫描完成,计算不同部位、不同速度对应的像点值,选取各部位的最大像点值对应的速度,就可得到各个深度岩石的最佳波速分布谱。至此就得到了用于制定地下采空区和岩溶形态特征图的第二方面数据。
[0057] 最后,对上述获得的两方面数据(散射波和反射波的能量分布谱、不同地段不同深度的最佳波速分布谱)进行合成孔径偏移成像,根据偏移叠加后的图像来识别地下采空区和岩溶的位置、形态和力学性质。其中,若速度低,则代表岩石强度低、岩石破碎程度高,该区域为采空区或者岩溶区。若速度高,则代表岩石强度高、岩石坚硬、完整度好。
[0058] 下面以贵州省毕威高速水塘隧道煤矿采空区地段为例,说明本发明的探测方法应用于实际工程试验时的效果。探测深度为120米,结果发现隧道仰拱下方15米处存在采空区,进行钻孔验证,检测结果与实际结果仅相差I米。误差仅为O. 83%,在工程领域是相当小的误差。具体探测过程如下。
[0059] (1)测线长60米。在地面上按照一定间距布置地震检波器10,每隔一个检波器点由机械震源20激发产生震动波,同时测量地震检波器10记录下来的散射波和反射波,共得到21个记录。
[0060] (2)对散射波和反射波进行方向滤波,滤除面波和直达波,得到地震检波器10记录的来自地下采空区和岩溶地质体的散射波和反射波。
[0061] (3)对散射波和反射波进行速度扫描,确定不同地段不同深度的围岩波速分布,如图7所示。
[0062] (4)对散射波和反射波进行合成孔径偏移成像,得到偏移图像,如图8所示。
[0063] (5)在偏移图像上识别,得到采空区等位置,如图8所示。
[0064] (6)根据上述结果,在现场钻孔验证。钻孔结果显示,在(5, 1930)处地下存在一采空区。
[0065] 以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种探测地下采空区和岩溶分布的方法,其特征在于, 步骤1,在待探测区域内依次间隔地布设多个地震检波器和多个机械震源;且,地震检波器之间的距离不大于所述机械震源产生的地震波波长的1/4 ;布设的多个地震检波器和多个机械震源形成的总排列长度不小于所述地震波的波长; 或者,在待探测区域内依次布设多个地震检波器后再布设一个机械震源;或者,在待探测区域内布设一个机械震源后再依次布设多个地震检波器; 其中,地震检波器之间的距离不大于所述机械震源产生的地震波波长的1/4 ;且,布设的多个地震检波器形成的总排列长度不小于所述地震波的波长; 步骤2,机械震源激发产生地震波,地震检波器采集地下震动信息; 步骤3,滤除所述地下震动信息中的杂波,得到来自地下采空区和岩溶的散射波和反射波,计算获得散射波和反射波的能量分布谱; 步骤4,对所述散射波和反射波进行速度扫描,确定地下岩体波速分布; 步骤5,制定地下采空区和岩溶的形态与分布特征图。
2.按照权利要求1所述的探测地下采空区和岩溶分布的方法,其特征在于, 在所述步骤3滤除杂波之前,先对所述地震检波器采集的地下震动信息进行预处理; 所述预处理包括去除噪声、延时校正、数据增益处理中的一种或多种。
3.按照权利要求1或2所述的探测地下采空区和岩溶分布的方法,其特征在于, 所述步骤3中,对所述地下震`动信息进行空间方向滤波,以滤除所述地下震动信息中的杂波。
4.按照权利要求1或2所述的探测地下采空区和岩溶分布的方法,其特征在于,所述杂波为面波、直达波、侧向回波中的一种或多种。
5.按照权利要求4所述的探测地下采空区和岩溶分布的方法,其特征在于,所述多个地震检波器和多个机械震源布设在同一水平线上。
6.按照权利要求4所述的探测地下采空区和岩溶分布的方法,其特征在于,所述多个地震检波器和一个机械震源布设在同一水平线上。
7.按照权利要求1、2、5或6所述的探测地下采空区和岩溶分布的方法,其特征在于, 所述步骤4中,对待探测区域地下散射波与反射波的偏移叠加能量按照以下公式求极值,递进改变波速计算叠加能量,当叠加能量取极值时对应的波速即为最佳波速; 按照以下公式计算不同部位的最佳波速: ——:-=0 01, 其中,V为介质波速; 11为距离地面的深度; X为水平方向的坐标。
8.按照权利要求7所述的探测地下采空区和岩溶分布的方法,其特征在于, 所述步骤5中,通过对步骤3得到的地下散射波和反射波的能量、以及步骤4得到的各个地层对应的最佳波速进行合成孔径偏移成像,制定所述地下采空区和岩溶的形态与分布特征图。
CN201110262208.7A 2011-09-06 2011-09-06 一种探测地下采空区和岩溶分布的方法 CN102426384B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110262208.7A CN102426384B (zh) 2011-09-06 2011-09-06 一种探测地下采空区和岩溶分布的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110262208.7A CN102426384B (zh) 2011-09-06 2011-09-06 一种探测地下采空区和岩溶分布的方法

Publications (2)

Publication Number Publication Date
CN102426384A CN102426384A (zh) 2012-04-25
CN102426384B true CN102426384B (zh) 2014-06-04

Family

ID=45960386

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110262208.7A CN102426384B (zh) 2011-09-06 2011-09-06 一种探测地下采空区和岩溶分布的方法

Country Status (1)

Country Link
CN (1) CN102426384B (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102862775B (zh) * 2012-04-28 2014-11-05 中国神华能源股份有限公司 一种矿井地下水的分布式存储方法
US9134444B2 (en) * 2012-11-27 2015-09-15 Chevron U.S.A., Inc. System and method for deducing cavern properties
CN103064111B (zh) * 2012-12-12 2015-11-18 中国石油天然气集团公司 一种基于形态滤波的微地震事件识别方法
CN103149274B (zh) * 2013-01-28 2015-10-07 中国科学院声学研究所 一种混凝土缺陷检测方法
CN104614770A (zh) * 2015-01-26 2015-05-13 北京市市政工程研究院 一种资源采空区的探查装置及方法
CN104749637B (zh) * 2015-04-21 2017-10-03 中煤科工集团重庆研究院有限公司 随钻地震钻头震源侧帮地质构造探测方法
CN105155502B (zh) * 2015-09-25 2017-01-25 青岛理工大学 岩溶溶洞型地基塌陷风险的测定方法
CN105589094A (zh) * 2015-12-16 2016-05-18 辽宁工程技术大学 一种不明采空区探测装置与方法
CN107643542B (zh) * 2017-08-28 2019-07-02 中煤科工集团西安研究院有限公司 基于井地联合的煤矿采空区三维弹性波层析探测方法
CN108845030A (zh) * 2018-05-03 2018-11-20 张家港申乾工程技术服务有限公司 一种在大口径管道内检测管道周边介质状态的方法
CN109343107B (zh) * 2018-09-30 2020-05-08 中国葛洲坝集团易普力股份有限公司 一种采空区精确且安全的钻孔探测方法
CN109441541B (zh) * 2018-11-06 2020-01-03 中国矿业大学 一种煤矿采空区充填体承载压缩率监测系统及其监测方法
CN109581486A (zh) * 2018-12-10 2019-04-05 上海宝冶集团有限公司 一种隧道溶洞的检测方法、装置和计算机可读存储介质
CN110531406B (zh) * 2019-08-20 2021-01-15 北京科技大学 一种地下隐伏病害精准探查方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5753818A (en) * 1995-05-15 1998-05-19 North American Geotechnical Company Method and apparatus for measuring scour around bridge foundations
CN101295027B (zh) * 2008-06-18 2011-03-30 上海大学 大堤隐患综合无损检测方法
CN101975089B (zh) * 2010-09-28 2012-11-21 中国水电顾问集团华东勘测设计研究院 一种用探地雷达首波相位法预报地下水的方法

Also Published As

Publication number Publication date
CN102426384A (zh) 2012-04-25

Similar Documents

Publication Publication Date Title
Foti et al. Guidelines for the good practice of surface wave analysis: A product of the InterPACIFIC project
Li et al. An overview of ahead geological prospecting in tunneling
CN102322294B (zh) 岩溶隧道施工用综合性地质预测预报方法
Park et al. Combined use of active and passive surface waves
NL1001100C2 (nl) Inrichting en werkwijze voor het lokaliseren van een bron van akoestische golven.
Xia et al. Comparing shear-wave velocity profiles from MASW with borehole measurements in unconsolidated sediments, Fraser River Delta, BC, Canada
CA2557526C (en) Prediction of shallow drilling hazards using seismic refraction data
Olsson et al. Borehole radar applied to the characterization of hydraulically conductive fracture zones in crystalline rock 1
Park et al. Multichannel analysis of surface waves (MASW)—active and passive methods
CN100487488C (zh) 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法
Krawczyk et al. Sinkholes in the city of Hamburg—New urban shear-wave reflection seismic system enables high-resolution imaging of subrosion structures
Martínez-Moreno et al. Integrated geophysical methods for studying the karst system of Gruta de las Maravillas (Aracena, Southwest Spain)
Cardarelli et al. Evaluation of tunnel stability using integrated geophysical methods
CN101403797B (zh) 地下工程施工超前地质预报系统及其预报方法
Juhlin et al. Reflection seismic imaging of the end-glacial Pärvie Fault system, northern Sweden
CN101907725B (zh) 裂缝预测方法和装置
Jeannin et al. Multiconfiguration GPR measurements for geometric fracture characterization in limestone cliffs (Alps)
CN102455439B (zh) 基于克希霍夫积分法的绕射波场分离方法
CN102353996B (zh) 钻孔内的定向瞬变电磁装置
Vanneste et al. Pseudo-3D imaging of a low-slip-rate, active normal fault using shallow geophysical methods: the Geleen fault in the Belgian Maas River valley
CN107589471B (zh) 一种高速铁路隧道综合超前地质预报方法
Anbazhagan et al. Mapping of average shear wave velocity for Bangalore region: a case study
CN102788991B (zh) 基于z分量的透射槽波高密度快速探测方法
CN102681004A (zh) 以掘进机为震源的巷道随掘地震超前探测装置及方法
Willenberg et al. Internal structure and deformation of an unstable crystalline rock mass above Randa (Switzerland): Part I—Internal structure from integrated geological and geophysical investigations

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
CP02 Change in the address of a patent holder

Address after: 541002 the Guangxi Zhuang Autonomous Region in Xiufeng District of Guilin City Road No. 31

Patentee after: Zhao Yonggui

Address before: 100101, room 1706, building 5, No. 108 East Fourth Ring Road, Chaoyang District, Beijing

Patentee before: Zhao Yonggui

C56 Change in the name or address of the patentee
CP02 Change in the address of a patent holder

Address after: 100101, room 1706, building 5, No. 108 East Fourth Ring Road, Chaoyang District, Beijing

Patentee after: Zhao Yonggui

Address before: 541002 the Guangxi Zhuang Autonomous Region in Xiufeng District of Guilin City Road No. 31

Patentee before: Zhao Yonggui

C56 Change in the name or address of the patentee
TR01 Transfer of patent right

Effective date of registration: 20160919

Address after: 100191 room 2, building 2, Huayuan Road 1112, Haidian District, Beijing, China

Patentee after: BEIJING TONGDU ENGINEERING GEOPHYSICS LTD. CORP.

Address before: 100101, room 1706, building 5, No. 108 East Fourth Ring Road, Chaoyang District, Beijing

Patentee before: Zhao Yonggui

C41 Transfer of patent application or patent right or utility model