CN113156536A - 用于地下转露天开采隐伏采空区的探测方法 - Google Patents

用于地下转露天开采隐伏采空区的探测方法 Download PDF

Info

Publication number
CN113156536A
CN113156536A CN202011619678.XA CN202011619678A CN113156536A CN 113156536 A CN113156536 A CN 113156536A CN 202011619678 A CN202011619678 A CN 202011619678A CN 113156536 A CN113156536 A CN 113156536A
Authority
CN
China
Prior art keywords
goaf
seismic source
elastic wave
underground
concealed
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
CN202011619678.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.)
Changsha Institute of Mining Research Co Ltd
Original Assignee
Changsha Institute of Mining Research Co Ltd
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 Changsha Institute of Mining Research Co Ltd filed Critical Changsha Institute of Mining Research Co Ltd
Priority to CN202011619678.XA priority Critical patent/CN113156536A/zh
Publication of CN113156536A publication Critical patent/CN113156536A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种用于地下转露天开采隐伏采空区的探测方法，包括如下步骤：S1、在隐伏采空区建立弹性波CT系统；S2、采用弹性波CT系统对隐伏采空区进行初步探测，确定是否存在采空区及圈定疑似采空区的大致分布区域；S3、采用三维激光扫描系统对圈定的疑似采空区进行精确探测，确定采空区的准确位置及分布形态。本发明结合工程探测法和地球物理探测法各自的优势，可探测出隐伏采空区的精确位置和形态，具备安全性好、探测精度高、成本低的优点。

Description

用于地下转露天开采隐伏采空区的探测方法
技术领域
本发明涉及矿山安全技术领域，尤其涉及一种用于地下转露天开采隐伏采 空区的探测方法。
背景技术
受前期无序开采影响，导致我国许多矿山地下采空区具有隐伏性强、空间 分布不明等特征。这些地下隐伏采空区已经成为制约矿山发展的一个重要难题， 尤其对地下转露天开采的矿山影响很大，如河南栾川钼矿、广东大宝山矿、鞍 钢弓长岭矿、太原袁家村铁矿等许多矿山都存在大量的地下隐伏采空区。由于 对地下采空区空间分布不明确，导致对采空区顶板冒落塌陷情况难以掌控，对 露天开采带来了严重的安全隐患。因此，如何对地下转露天开采隐伏采空区的 分布特征进行探测是保证这类矿山安全生产的前提和关键技术难题。
对于隐伏采空区的探测，目前主要采用工程探测法和地球物理探测法两大 类。工程探测法是利用钻机对疑似采空区的区域实施高密度钻孔，一般同时会 结合三维激光扫描技术对探测到的采空区进行分布形态扫描，为采空区处理提 供准确数据。该方法探测结果较为准确，但是成本较高，一般多用于对已经初 步圈定采空区的确认。地球物理探测方法具有安全性好、成本低、探测工作条 件好等特点，但是探测效果受实际地形、地质条件和岩体介质物理力学性质的 限制和影响，对于隐伏采空区分布形态的精确探测效果相对较差，一般多用于 隐伏采空区的初步探测和疑似采空区区域的圈定，目前常用的方法主要有常规 电法、浅层地震法、微重力法、探地雷达法、瞬变电磁法、高密度电法、弹性 波CT技术等。
发明内容
本发明的目的在于针对现有技术的不足之处，提供一种可以实现采空区分 布位置和形态精确探测的用于地下转露天开采隐伏采空区的探测方法。
本发明提供的这种用于地下转露天开采隐伏采空区的探测方法，包括如下 步骤：
S1、在隐伏采空区建立弹性波CT系统；
S2、采用弹性波CT系统对隐伏采空区进行初步探测，确定是否存在采空区 及圈定疑似采空区的大致分布区域；
S3、采用三维激光扫描系统对圈定的疑似采空区进行精确探测，确定采空 区的准确位置及分布形态。
在所述步骤S1中，弹性波CT系统的建立过程如下：
S11、在地面上布置网格状的钻孔；
S12、选择任意两个钻孔作为震源发射钻孔和信号接收钻孔；
S13、在震源发射钻孔和信号接收钻孔内按相同的间距布置震源和接收器， 将震源、接收器和数据采集器连接好后，构建弹性波CT系统。
在步骤S11中，钻孔在地面的露天台阶上按照36m*36m的网度进行分布， 单个钻孔深度36m，相邻钻孔间组成观测线跨度为36m的弹性波CT系统，对角 钻孔间组成观测线跨度为51m的弹性波CT系统。
在步骤S13中，震源发射钻孔内每间隔1m激发一次震源，信号接收钻孔内 每间隔1m布置一个检波器串。
所述步骤S2的具体操作过程如下：
S21、激发震源，通过接收器接收震源信号，采用数据采集器对接收器接收 的信号进行采集，完成钻孔间区域岩体弹性波CT数据的采集；
S22、对采集的弹性波CT数据进行数据校对、初值拾取、计算模型建立和 反演成像，给出数据分析结果；
S23、根据数据分析结果确定是否存在疑似为采空区的低波速异常区域并圈 定低波速异常区域的大致分布区域。
所述步骤S3的具体操作过程如下：对圈定的低波速异常区域再次进行三维 激光扫描，最终确定采空区的准确分布位置和形态。
在步骤S21中，激发震源的具体过程如下：现场将弹性波CT系统连接好后， 开始对电火花震源进行充电，电压达到预定值后激发震源，激发一次后将震源 位置上升或下降1m，再次激发震源，循环操作直至全部采集完毕。
本发明结合工程探测法和地球物理探测法各自的优势，可探测出隐伏采空 区的精确位置和形态，具备安全性好、探测精度高、成本低的优点。
附图说明
图1为本发明中弹性波CT系统钻孔布置示意图。
图2为本发明中弹性波CT系统结构示意图。
图3为本发明中弹性波CT数据分析示意图。
图4为采空区三维激光扫描示意图。
图5采空区三维激光扫描结果示意图。
图6为本发明中弹性波CT波速模型。
图7为本发明中三维激光扫描原理示意图。
具体实施方式
本发明这种用于地下转露天开采隐伏采空区的探测方法，包括如下步骤：
S1、在隐伏采空区建立弹性波CT系统
S11、在地面上钻取n个呈网格状布置的钻孔；
S12、选择任意两个钻孔作为震源发射钻孔和信号接收钻孔；
S13、在震源发射钻孔和信号接收钻孔内按相同的间距布置震源和接收器， 将震源、接收器和数据采集器连接好后，构建弹性波CT系统；
S2、采用弹性波CT系统对隐伏采空区进行初步探测，确定是否存在采空区 及圈定疑似采空区的大致分布区域
S21、激发震源，通过接收器接收震源信号，采用数据采集器对接收器接收 的信号进行采集，完成钻孔间区域岩体弹性波CT数据的采集；
S22、对采集的弹性波CT数据进行数据校对、初值拾取、计算模型建立和 反演成像，给出数据分析结果；
S23、根据数据分析结果确定是否存在疑似为采空区的低波速异常区域并圈 定低波速异常区域的大致分布区域；
S3、采用三维激光扫描系统对圈定的疑似采空区进行精确探测，确定采空 区的准确位置及分布形态
对圈定的低波速异常区域再次进行三维激光扫描，最终确定采空区的准确 分布位置和形态。
从图1可以看出，在本发明的的步骤S11中，钻孔有四个，四个钻孔在地 面的露天台阶上按照36m*36m的网度进行分布，单个钻孔深度36m，相邻钻孔间 组成观测线1跨度为36m的弹性波CT系统，对角钻孔间组成观测线2跨度为51m 的弹性波CT系统。
在本发明的步骤S13中，震源发射钻孔内每间隔1m激发一次震源，信号接 收钻孔内每间隔1m布置一个检波器串。
弹性波CT系统主要包括震源、接收器、数据采集器和数据分析软件。弹性 波CT系统要求震源能量大、重复性好、计时精确，并且发射的震动信号频带宽、 高频成份丰富，这些要求比较难以完全满足，根据不同的任务要求，可选用空 气枪、电火花、炸药、雷管束、超磁致伸缩、专用可控震源等，目前最常用的 是电火花震源。接收器一般采用压电陶瓷传感器，压电陶瓷具有高灵敏度、宽 频带、大动态范围、相位一致性好、谐波失真小、耦合性好稳定性好、寿命长、 抗干扰能力强等特性。为了提高效率，常采用等间距的检波器串同时接收。数 据采集器常用12道或24道地震仪或类似的多道工程探测仪。如图2所示，在 本发明中，弹性波CT系统的震源采用电火花震源，数据采集器采用24道地震 仪。
在本发明的步骤S21中，激发震源的具体过程如下：如图2所示，在现场 将弹性波CT系统连接好后，开始对电火花震源进行充电，电压达到预定值后激 发震源，激发一次后将震源位置上升或下降1m，再次激发震源，循环操作直至 全部采集完毕。
在本发明的步骤S22中，对数据是否完整、重复及错误进行校对；如图3 左半图所示，校对后开始初值拾取，可采用手动、半自动、全自动三种拾取方 式，一般先选择全自动拾取，再进行手动拾取调整；计算模型建立包括网格划 分、速度模型确定、迭代算法及收敛条件确定，完成后进行反演成像计算，如 图3右半图所示，为典型的弹性波CT反演结果，可以看出存在低波速异常区域， 该区域基本可以判断存在采空区。
在本发明中，三维激光扫描系统采用C-ALS采空区三维激光自动扫描系统， 如图4所示，该系统可以将激光的扫描头通过钻孔进入到空区内部进行扫描， 扫描控制仪内置导航模块，可测量空区内扫描头的位置和方向。扫描控制仪可 以360°全方位无盲区扫描，设计的水平和垂直扫描方式可以确保对采空区整个 空间的扫描，如图5所示，通过自带的软件可以得到采空区的三维模型，并可 以和CAD、3Dimine和Surpac等常用数字矿山软件进行对接。
透射型波速CT成像原理如图6所示，假设弹性波在二维空间上的波速分布 为V(x,z)，那么沿射线l的走时为：
Figure BDA0002873847670000051
如图6所示，S为震源，R为接收点，1、2、…16为网格编号，从震源到接 收点发出的第i条射线穿过第j个网格的长度用li,j表示。
若将二维平面划分成网格，每一小网格的速度为Vj(j＝1,2,3,…,n＝M×N)，那么沿着射线l的的走时可以写为：
Figure BDA0002873847670000052
其中xj是网格j的慢度值，cij为第j个像元对走时yi的贡献，亦即第i条射线 在第j个像元上的射线长度，如果有n个像元，m条射线，那么它可以表示为：
CX＝Y (3)
其中C为m×n的矩阵，每个元素为cij，X＝(x1,x2,…,xn)，Y＝(y1,y2,…,yn)。对 方程组进行求解，得到测区内各节点慢度分布，最后得到速度分布图像。
三维激光扫描技术是建立在激光测距原理基础上，扫描过程中，从三维激 光扫描仪的扫描头发出一束激光脉冲，发射的激光到达空区边界后反射回仪器 可记录激光到达时间T，同时激光扫描仪可以任意旋转激光脉冲横向扫描角度值 β和纵向扫描角度值α。通过建立坐标系统，由此可以得到目标点P点的坐标 计算公式。
如图7所示，激光发射点O距离目标点P的距离L为：
Figure BDA0002873847670000053
，式中：c为激光在大气中的传播速度；T为激光发射后到达目标点反射返 回的传播时间；
待测目标点的坐标为：
Figure BDA0002873847670000061
，式中：L为激光发射点距离待测目标点的距离；α为激光纵向扫描角度 值；β为激光横向扫描角度值。

Claims (7)

1.一种用于地下转露天开采隐伏采空区的探测方法，其特征在于，包括如下步骤：
S1、在隐伏采空区建立弹性波CT系统；
S2、采用弹性波CT系统对隐伏采空区进行初步探测，确定是否存在采空区及圈定疑似采空区的大致分布区域；
S3、采用三维激光扫描系统对圈定的疑似采空区进行精确探测，确定采空区的准确位置及分布形态。
2.根据权利要求1所述的用于地下转露天开采隐伏采空区的探测方法，其特征在于，在所述步骤S1中，弹性波CT系统的建立过程如下：
S11、在地面上布置网格状的钻孔；
S12、选择任意两个钻孔作为震源发射钻孔和信号接收钻孔；
S13、在震源发射钻孔和信号接收钻孔内按相同的间距布置震源和接收器，将震源、接收器和数据采集器连接好后，构建弹性波CT系统。
3.根据权利要求2所述的用于地下转露天开采隐伏采空区的探测方法，其特征在于：在步骤S11中，钻孔在地面的露天台阶上按照36m*36m的网度进行分布，单个钻孔深度36m，相邻钻孔间组成观测线跨度为36m的弹性波CT系统，对角钻孔间组成观测线跨度为51m的弹性波CT系统。
4.根据权利要求2所述的用于地下转露天开采隐伏采空区的探测方法，其特征在于：在步骤S13中，震源发射钻孔内每间隔1m激发一次震源，信号接收钻孔内每间隔1m布置一个检波器串。
5.根据权利要求2所述的用于地下转露天开采隐伏采空区的探测方法，其特征在于，所述步骤S2的具体操作过程如下：
S21、激发震源，通过接收器接收震源信号，采用数据采集器对接收器接收的信号进行采集，完成钻孔间区域岩体弹性波CT数据的采集；
S22、对采集的弹性波CT数据进行数据校对、初值拾取、计算模型建立和反演成像，给出数据分析结果；
S23、根据数据分析结果确定是否存在疑似为采空区的低波速异常区域并圈定低波速异常区域的大致分布区域。
6.根据权利要求5所述的用于地下转露天开采隐伏采空区的探测方法，其特征在于，所述步骤S3的具体操作过程如下：对圈定的低波速异常区域再次进行三维激光扫描，最终确定采空区的准确分布位置和形态。
7.根据权利要求5所述的用于地下转露天开采隐伏采空区的探测方法，其特征在于，在步骤S21中，激发震源的具体过程如下：现场将弹性波CT系统连接好后，开始对电火花震源进行充电，电压达到预定值后激发震源，激发一次后将震源位置上升或下降1m，再次激发震源，循环操作直至全部采集完毕。
CN202011619678.XA 2020-12-30 2020-12-30 用于地下转露天开采隐伏采空区的探测方法 Pending CN113156536A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011619678.XA CN113156536A (zh) 2020-12-30 2020-12-30 用于地下转露天开采隐伏采空区的探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011619678.XA CN113156536A (zh) 2020-12-30 2020-12-30 用于地下转露天开采隐伏采空区的探测方法

Publications (1)

Publication Number Publication Date
CN113156536A true CN113156536A (zh) 2021-07-23

Family

ID=76878171

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011619678.XA Pending CN113156536A (zh) 2020-12-30 2020-12-30 用于地下转露天开采隐伏采空区的探测方法

Country Status (1)

Country Link
CN (1) CN113156536A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113805250A (zh) * 2021-08-23 2021-12-17 鞍钢集团矿业有限公司 一种露天矿采空区高密度地震与电法联合探测及识别方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107741459A (zh) * 2017-11-14 2018-02-27 中国电建集团贵阳勘测设计研究院有限公司 利用孔间声波层析成像技术探测岩性异常体的方法及装置
CN110531406A (zh) * 2019-08-20 2019-12-03 北京科技大学 一种地下隐伏病害精准探查方法及装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107741459A (zh) * 2017-11-14 2018-02-27 中国电建集团贵阳勘测设计研究院有限公司 利用孔间声波层析成像技术探测岩性异常体的方法及装置
CN110531406A (zh) * 2019-08-20 2019-12-03 北京科技大学 一种地下隐伏病害精准探查方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
彭府华: "地下转露天开采隐伏采空区综合探测技术应用研究", 《中国钨业》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113805250A (zh) * 2021-08-23 2021-12-17 鞍钢集团矿业有限公司 一种露天矿采空区高密度地震与电法联合探测及识别方法
CN113805250B (zh) * 2021-08-23 2024-01-12 鞍钢集团矿业有限公司 一种露天矿采空区高密度地震与电法联合探测及识别方法

Similar Documents

Publication Publication Date Title
Li et al. An overview of ahead geological prospecting in tunneling
CN109613616B (zh) 一种浅层岩溶隧道地质综合立体探测方法
CN102866417A (zh) 一种地下溶洞地震跨孔ct探测及层析成像装置及方法
CN110988981B (zh) 一种适用于钻爆法隧道的相控阵声波超前预报系统及方法
CN109407144B (zh) 一种基于多波的单孔孤石立体探测方法
CN202837558U (zh) 一种地下溶洞地震跨孔ct探测及层析成像装置
CN108957548B (zh) 一种多波多分量联合观测地震页岩气富集区预测方法
CN113419294A (zh) 一种多维度岩溶特殊地质综合探测方法
Takahashi et al. ISRM suggested methods for borehole geophysics in rock engineering
CN102877828A (zh) 一种三维多井联合井地ct成像方法
CN110632667B (zh) 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法
Luo et al. Tomographic imaging of rock conditions ahead of mining using the shearer as a seismic source—a feasibility study
CN113156536A (zh) 用于地下转露天开采隐伏采空区的探测方法
CN208334667U (zh) 一种用于软岩隧道长距离三维超前地质预报装置
CN113050085A (zh) 一种超前地质预报的方法
CN110850472B (zh) 一种基于冲击波激发震源的可变偏移距超前探测断层方法
CN112130206A (zh) 一种井下巷道全空间地震探测系统及探测方法
CN111158050B (zh) 数据采集系统、方法及隧道地震波超前预报方法
CN112965139A (zh) 一种复杂地质条件隧道超前地质综合预报方法
Sanfirov et al. Frozen wall construction control in mine shafts using land and borehole seismology techniques
CN115182736A (zh) 隧道的施工方法
CN109375251B (zh) 利用城市既有地下空间与地表的探测方法及系统
CN110045414B (zh) 一种矿区深部金属矿的探测方法
CN112414290A (zh) 一种用地质雷达检测铝矿直接顶板厚度的方法
Giroux et al. Geotechnical application of borehole GPR-a case history

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

Application publication date: 20210723

RJ01 Rejection of invention patent application after publication