CN104267442A - 一种用于煤矿井下的瞬变电磁拟地震探测方法 - Google Patents
一种用于煤矿井下的瞬变电磁拟地震探测方法 Download PDFInfo
- Publication number
- CN104267442A CN104267442A CN201410490742.7A CN201410490742A CN104267442A CN 104267442 A CN104267442 A CN 104267442A CN 201410490742 A CN201410490742 A CN 201410490742A CN 104267442 A CN104267442 A CN 104267442A
- Authority
- CN
- China
- Prior art keywords
- transient
- coal mine
- data
- field signal
- transient electromagnetic
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种探测方法,主要用于煤矿井下探测,具体涉及一种用于用于煤矿井下的瞬变电磁拟地震探测方法。该方法对每一探测网格增加了测量数据,大幅度减弱巷道内环境因素带来的影响,提高了通过异常体的探测数据密度,并将瞬变场数据转换为了虚拟波场数据,消弱了电磁法的体积效应影响,同时,该方法采用地震的前沿精细处理技术,有利于地质异常的精细解释,能显著提高对矿井水害隐患的探查精度。
Description
技术领域
本发明涉及一种探测方法,主要用于煤矿井下探测,具体涉及一种用于煤矿井下的瞬变电磁拟地震探测方法。
背景技术
瞬变电磁法是利用不接地回线或接地线源向地下发射一次脉冲磁场,在一次脉冲磁场的间隙期间,利用线圈或接地电极观测二次涡流场的方法。
小线框瞬变电磁法是近年来广泛应用于煤矿井下超前探、顶底板水害探测、侧帮及工作面内部水害探测等领域的一种电磁类勘探方法。但是目前,现有技术中的该方法存在以下问题:
1、线框激发能量有限:由于井下施工空间限制,该方法不能采用地面常规使用的大发射线框(一般边长大于100m),往往使用的是边长小于2m的多匝小回线,激发能量有限;
2、设备运行不稳定:单点垂直测线方向激发,同方向接收的探测方式造成穿过异常体的测点数目很少,且受巷道内空间、金属物件、电力器件等的影响,往往出现异常无法判断或者出现大量假异常等情况;
3、隐患探查精度低:资料处理目前主要靠移植地面瞬变电磁一维定性处理技术,解释精度较低,对资料的精细解释不利,影响煤矿井下的水害隐患探查精度,对矿井安全生产带来威胁。
发明内容
本发明主要是解决现有技术所存在的发射线框激发能量有限,设备运行不稳定,隐患探查精度低等问题,提出了一种用于煤矿井下的瞬变电磁拟地震探测方法,该方法对每一探测网格增加了测量数据,大幅度减弱巷道内环境因素带来的影响,提高了通过异常体的探测数据密度,并将瞬变场数据转换为了虚拟波场数据,消弱了电磁法的体积效应影响,同时,该方法采用地震的前沿精细处理技术,有利于地质异常的精细解释,能显著提高对矿井水害隐患的探查精度。
本发明的上述技术问题主要是通过下述技术方案得以解决的:
一种用于煤矿井下的瞬变电磁拟地震探测方法,包括:
步骤1,根据探测区域设计测线方向,在测线上按照一定点距布置测点,并在每一测点上均布置接收装置;
步骤2,依次在每一个测点上布置发射线框并发射瞬变信号,并在其余所有测点同时接收采集瞬变场信号;
步骤3,将步骤2中采集到的瞬变场信号转换为虚拟波动场信号,并对转换后的虚拟波动场信号进行波形压缩;
步骤4,以步骤3中的虚拟波动场信号的传播时间为横坐标,以其波形幅值为纵坐标,通过描点得到测线波形图;
步骤5,根据瞬变电磁定性算法处理步骤2中采集到的瞬变场信号,以得到视电阻率分布,并将所述视电阻率分布作为步骤3中转换后的虚拟波动场信号的初始速度场;
步骤6,在所述测线波形图上选取参考子波,并以步骤5中的初始速度场为初始模型,采用波形反演技术对步骤3中的压缩后的虚拟波动场信号进行二维或三维反演;
步骤7,对反演后得到的虚拟波场信号的速度数据进行网格化,绘制等值线图并成像;
步骤8,依据成像结果分析,将速度值出现明显的低值分布的区域作为地质异常体。
优化的,上述的一种用于煤矿井下的瞬变电磁拟地震探测方法,所述步骤3中基于以下公式将瞬变场所数据转换为波动数据:
式中hz(t)为瞬变场数据,U(τ)为波动场数据,t为采样时间,τ为与t对应的转换后的虚拟波场传播时间,K=tn为波形均衡因子,n为(0,2.5)区间内的任意值。
优化的,上述的一种用于煤矿井下的瞬变电磁拟地震探测方法,所述步骤1中的点距为5米。
优化的,上述的一种用于煤矿井下的瞬变电磁拟地震探测方法,其特征在于,所述步骤1中的接收装置为探头或线圈。
因此,本发明具有如下优点:1、对每一探测网格增加了测量数据,可大幅度减弱巷道内环境因素带来的影响;2、提高了通过异常体的探测数据密度,并将瞬变场数据转换为了虚拟波场数据,消弱了电磁法的体积效应影响;3、采用地震的前沿精细处理技术,有利于地质异常的精细解释,能显著提高对矿井水害隐患的探查精度。
附图说明
图1为传统煤矿井下瞬变电磁探测方法示意图。
图2为煤矿井下瞬变电磁拟地震探测方法示意图。
图3为直接转换得到的波场数据。
图4为经过反褶积压缩后的锐化波形。
图5为瞬变电磁实测衰减曲线图。
图6经过本发明方法得到的虚拟波场数据。
图7是常规瞬变电磁方法处理得到视电阻率剖面。
图8是对图7进行精细反演得到的虚拟波场速度剖面。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明。图中,煤层1,煤板岩层2,巷道空间3,发射线框4,地质异常体5,探测射线6,接收点7。
实施例:
图1是传统煤矿井下瞬变电磁探测方法示意图。传统的瞬变电磁探测方法采用单点垂直测线方向激发,同方向接收的探测方式造成穿过异常体的测点数目很少。
图2所示是改进的用于煤矿井下的瞬变电磁拟地震探测方法。从图中可以看出该方法对每一探测网格增加了测量数据,提高了通过异常体的探测数据密度,并将瞬变场数据转换为了虚拟波场数据,消弱了电磁法的体积效应影响。
本实施例通过以下步骤实现:
(1)、根据探测区域设计测线方向,在测线上按照一定点距(如5m)布置测点;
(2)、在每一测点上均布置接收装置(探头、线圈均可);
(3)、在每一个测点上依次布置发射线框并发射瞬变信号,其余所有测点同时接收采集;
(4)、将数据从仪器导入电脑;
(5)、根据以下公式将采集到的瞬变场信号转换为虚拟波动场信号,并对转换后的虚拟波动场信号进行波形压缩;
式中hz(t)为瞬变场数据,U(τ)为波动场数据,t为采样时间,τ为与t对应的转换后的虚拟波场传播时间,K=tn为波形均衡因子,n为(0,2.5)区间内的任意值。
由于采集到的信号为瞬变场信号,通过公式可转换为虚拟波动场信号,瞬变场为扩散场,该转换相当于将扩散场中的波动特征提取了出来。
公式中的波形均衡因子K可以对瞬变场转换过来的波形做均衡,直接转换会出现早期的波形能看见,晚期的波形幅值太小看不见,通过调整K可以压制前期的波形幅值,抬升晚期的波形幅值,使前后波形差异缩小,有利于深部信息的反演成像。n取0-2.5与瞬变电磁法的基本公式有关,是根据下面公式分母上t的幂指数确定:
式中:μ0—磁导率;
M—发送线圈磁矩;
q—接收线圈等效面积;
ρ—地层电阻率;
t—时间;
V—瞬变电磁二次场信号。
(6)、对转换后的波场数据进行波形压缩;直接转换得到的波场数据中波的特征不明显,不同时间的波峰波谷混在一起,如图3所示,经过反褶积压缩后能提取出锐化波形,如图4所示。这一压缩有利于提高反演的分辨率。
(7)、根据瞬变电磁定性算法给出视电阻率分布作为初始速度场;
(8)、以虚拟波动场信号的传播时间为横坐标,以其波形幅值为纵坐标,通过描点得到测线波形图,在测线波形图上选取参考子波;
(9)、采用波形匹配反演技术对数据进行二维或三维反演;
(10)、对反演后得到的虚拟波场速度数据进行网格化,绘制等值线图成像;
(11)、依据成像结果解释可能存在的地质异常体。
采用本实施例处理后的图形如图5-8所示。图5为瞬变电磁实测衰减曲线图,经过本发明处理后,可以得到如图6所示的虚拟波场数据。图7是常规瞬变电磁方法处理得到视电阻率剖面(探测目标为钻孔,图中圆点所示)。对图7进行精细反演可以得到图8所示的虚拟波场速度剖面(探测目标为钻孔,图中圆点所示)。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
尽管本文较多地使用了煤层1,煤板岩层2,巷道空间3,发射线框4,地质异常体5,探测射线6,接收点7术语,但并不排除使用其它术语的可能性。使用这些术语,是为了更方便地描述和解释本发明的本质;把它们解释成任何一种附加的限制都是与本发明精神相违背的。
Claims (4)
1.一种用于煤矿井下的瞬变电磁拟地震探测方法,其特征在于,包括:
步骤1,根据探测区域设计测线方向,在测线上按照一定点距布置测点,并在每一测点上均布置接收装置;
步骤2,依次在每一个测点上布置发射线框并发射瞬变信号,并在其余所有测点同时接收采集瞬变场信号;
步骤3,将步骤2中采集到的瞬变场信号转换为虚拟波动场信号,并对转换后的虚拟波动场信号进行波形压缩;
步骤4,以步骤3中的虚拟波动场信号的传播时间为横坐标,以其波形幅值为纵坐标,通过描点得到测线波形图;
步骤5,根据瞬变电磁定性算法处理步骤2中采集到的瞬变场信号,以得到视电阻率分布,并将所述视电阻率分布作为步骤3中转换后的虚拟波动场信号的初始速度场;
步骤6,在所述测线波形图上选取参考子波,并以步骤5中的初始速度场为初始模型,采用波形反演技术对步骤3中的压缩后的虚拟波动场信号进行二维或三维反演;
步骤7,对反演后得到的虚拟波场信号的速度数据进行网格化,绘制等值线图并成像;
步骤8,依据成像结果分析,将速度值出现明显的低值分布的区域作为地质异常体。
2.根据权利要求1所述的一种用于煤矿井下的瞬变电磁拟地震探测方法,其特征在于,所述步骤3中基于以下公式将瞬变场所数据转换为波动数据:
式中hz(t)为瞬变场数据,U(τ)为波动场数据,t为采样时间,τ为与t对应的转换后的虚拟波场传播时间,K=tn为波形均衡因子,n为(0,2.5)区间内的任意值。
3.根据权利要求1所述的一种用于煤矿井下的瞬变电磁拟地震探测方法,其特征在于,所述步骤1中的点距为5米。
4.根据权利要求1所述的一种用于煤矿井下的瞬变电磁拟地震探测方法,其特征在于,所述步骤1中的接收装置为探头或线圈。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410490742.7A CN104267442B (zh) | 2014-09-23 | 2014-09-23 | 一种用于煤矿井下的瞬变电磁拟地震探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410490742.7A CN104267442B (zh) | 2014-09-23 | 2014-09-23 | 一种用于煤矿井下的瞬变电磁拟地震探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104267442A true CN104267442A (zh) | 2015-01-07 |
CN104267442B CN104267442B (zh) | 2016-08-24 |
Family
ID=52158981
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410490742.7A Active CN104267442B (zh) | 2014-09-23 | 2014-09-23 | 一种用于煤矿井下的瞬变电磁拟地震探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104267442B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104614779A (zh) * | 2015-01-16 | 2015-05-13 | 中国矿业大学 | 一种多参数电磁法动态监测系统及其方法 |
CN104793251A (zh) * | 2015-04-10 | 2015-07-22 | 河南理工大学 | 中远距离保护层开采保护范围的测试装置及其测试方法 |
CN106054260A (zh) * | 2016-05-18 | 2016-10-26 | 吉林大学 | 一种借助平行导硐的实时隧道超前预报方法 |
CN106199733A (zh) * | 2016-06-30 | 2016-12-07 | 江苏省有色金属华东地质勘查局地球化学勘查与海洋地质调查研究院 | 用于评价地热富水地质异常体的方法及装置 |
CN106199732A (zh) * | 2016-06-27 | 2016-12-07 | 中国科学院地质与地球物理研究所 | 一种瞬变电磁多道覆盖观测装置和方法 |
CN106970424A (zh) * | 2017-03-17 | 2017-07-21 | 中煤科工集团西安研究院有限公司 | 一种煤矿井下孔巷瞬变电磁叠加超前探测装置及方法 |
CN109752762A (zh) * | 2018-12-11 | 2019-05-14 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN110133726A (zh) * | 2019-04-25 | 2019-08-16 | 中铁二院工程集团有限责任公司 | 铁路隧道航空电磁法勘探测线布置方法 |
CN112147707A (zh) * | 2020-08-06 | 2020-12-29 | 中煤科工集团西安研究院有限公司 | 钻孔瞬变电磁一维反演电性信息立体成像方法 |
CN112394418A (zh) * | 2020-11-06 | 2021-02-23 | 天津大学 | 一种近地表瞬变电磁感应高分辨率检测系统 |
CN112799139A (zh) * | 2020-12-29 | 2021-05-14 | 河北煤炭科学研究院有限公司 | 一种赋水异常区域三维探测方法及装置 |
US11567230B1 (en) * | 2021-09-03 | 2023-01-31 | University Of Science And Technology Beijing | Direction-finding and positioning system of electromagnetic emission of coal or rock fracture |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
CN102508308A (zh) * | 2011-11-03 | 2012-06-20 | 山东大学 | 一种隧道瞬变电磁多点阵列式探测方法 |
US20140218017A1 (en) * | 2012-01-27 | 2014-08-07 | Texas A&M University System | Transient electromagnetic geophysics system with cooperative source(s) |
-
2014
- 2014-09-23 CN CN201410490742.7A patent/CN104267442B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102419453A (zh) * | 2011-07-15 | 2012-04-18 | 中国科学院地质与地球物理研究所 | 长导线源瞬变电磁地空探测方法 |
CN102508308A (zh) * | 2011-11-03 | 2012-06-20 | 山东大学 | 一种隧道瞬变电磁多点阵列式探测方法 |
US20140218017A1 (en) * | 2012-01-27 | 2014-08-07 | Texas A&M University System | Transient electromagnetic geophysics system with cooperative source(s) |
Non-Patent Citations (4)
Title |
---|
G.Q.XUE ET AL.: "《3D pseudo-seismic imaging of transient electromagnetic data—a feasibility study》", 《GEOPHYSICAL PROSPECTING》 * |
G.Q.XUE ET AL.: "《Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geolopy exploration》", 《GEOPHYSICAL RESEARCH LETTERS》 * |
朱宏伟等: "《瞬变电磁法三维拟地震成像信息提取技术》", 《地球物理学进展》 * |
范涛等: "《瞬变电磁方法在探测煤矿浅层高阻采空区中的应用》", 《西北地质》 * |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104614779A (zh) * | 2015-01-16 | 2015-05-13 | 中国矿业大学 | 一种多参数电磁法动态监测系统及其方法 |
CN104793251A (zh) * | 2015-04-10 | 2015-07-22 | 河南理工大学 | 中远距离保护层开采保护范围的测试装置及其测试方法 |
CN104793251B (zh) * | 2015-04-10 | 2017-05-17 | 河南理工大学 | 中远距离保护层开采保护范围的测试方法 |
CN106054260A (zh) * | 2016-05-18 | 2016-10-26 | 吉林大学 | 一种借助平行导硐的实时隧道超前预报方法 |
CN106199732A (zh) * | 2016-06-27 | 2016-12-07 | 中国科学院地质与地球物理研究所 | 一种瞬变电磁多道覆盖观测装置和方法 |
CN106199733A (zh) * | 2016-06-30 | 2016-12-07 | 江苏省有色金属华东地质勘查局地球化学勘查与海洋地质调查研究院 | 用于评价地热富水地质异常体的方法及装置 |
CN106970424B (zh) * | 2017-03-17 | 2019-05-28 | 中煤科工集团西安研究院有限公司 | 一种煤矿井下孔巷瞬变电磁叠加超前探测装置及方法 |
CN106970424A (zh) * | 2017-03-17 | 2017-07-21 | 中煤科工集团西安研究院有限公司 | 一种煤矿井下孔巷瞬变电磁叠加超前探测装置及方法 |
CN109752762A (zh) * | 2018-12-11 | 2019-05-14 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN109752762B (zh) * | 2018-12-11 | 2020-10-16 | 中煤科工集团西安研究院有限公司 | 单发多收观测装置瞬变电场数据动校正方法及装置 |
CN110133726A (zh) * | 2019-04-25 | 2019-08-16 | 中铁二院工程集团有限责任公司 | 铁路隧道航空电磁法勘探测线布置方法 |
CN110133726B (zh) * | 2019-04-25 | 2024-03-22 | 中铁二院工程集团有限责任公司 | 铁路隧道航空电磁法勘探测线布置方法 |
CN112147707A (zh) * | 2020-08-06 | 2020-12-29 | 中煤科工集团西安研究院有限公司 | 钻孔瞬变电磁一维反演电性信息立体成像方法 |
CN112147707B (zh) * | 2020-08-06 | 2023-01-24 | 中煤科工集团西安研究院有限公司 | 钻孔瞬变电磁一维反演电性信息立体成像方法 |
CN112394418A (zh) * | 2020-11-06 | 2021-02-23 | 天津大学 | 一种近地表瞬变电磁感应高分辨率检测系统 |
CN112394418B (zh) * | 2020-11-06 | 2023-03-17 | 天津大学 | 一种近地表瞬变电磁感应高分辨率检测系统 |
CN112799139A (zh) * | 2020-12-29 | 2021-05-14 | 河北煤炭科学研究院有限公司 | 一种赋水异常区域三维探测方法及装置 |
US11567230B1 (en) * | 2021-09-03 | 2023-01-31 | University Of Science And Technology Beijing | Direction-finding and positioning system of electromagnetic emission of coal or rock fracture |
Also Published As
Publication number | Publication date |
---|---|
CN104267442B (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104267442A (zh) | 一种用于煤矿井下的瞬变电磁拟地震探测方法 | |
US10613245B2 (en) | Mine TEM three-component detection method | |
CN101603423B (zh) | 一种在煤矿巷道内顺层超前探测含水构造的直流电法方法 | |
CN109613616B (zh) | 一种浅层岩溶隧道地质综合立体探测方法 | |
CN103995296B (zh) | 瞬变电磁法地孔探测方法与装置 | |
CN102788991B (zh) | 基于z分量的透射槽波高密度快速探测方法 | |
CN103995295B (zh) | 直流电法地孔探测方法与装置 | |
CN103777247B (zh) | 一种瞬变电磁雷达探测系统及探测方法 | |
CN103245977B (zh) | 一种矿井回采区灾害源的地质雷达层析探测方法 | |
CN104614774B (zh) | 一种瞬变电磁探测方法、装置和系统 | |
CN102426384A (zh) | 一种探测地下采空区和岩溶分布的方法 | |
CN104035137A (zh) | 地下全空间瞬变电磁探测仪及探测方法 | |
CN106154310A (zh) | 一种便携式城市地质异常地震探测系统及方法 | |
CN105510981A (zh) | 一种磁铁矿采空区地球物理判定方法和装置 | |
CN103630946A (zh) | 单孔电磁波层析成像超前探测装置及方法 | |
CN104765072A (zh) | 一种用环形天线旋转进行磁共振超前探测的方法 | |
Lv et al. | Application of GPR reverse time migration in tunnel lining cavity imaging | |
CN104459763A (zh) | 紧支集小波探测地下空洞位置的方法及系统 | |
CN105277976A (zh) | 基于岩石露头雷达探测的地震正演模拟方法 | |
CN111538093A (zh) | 一种用于浅层地表探测方法及瞬变电磁仪器 | |
CN106371137A (zh) | 用于浅层异常体地震探测装置及其三维观测方法 | |
CN208334667U (zh) | 一种用于软岩隧道长距离三维超前地质预报装置 | |
CN113703058A (zh) | 一种利用视电导率和相对介电常数探测地下障碍物的方法 | |
CN104142522B (zh) | 一种城市隐伏断裂的探测方法 | |
CN113050085A (zh) | 一种超前地质预报的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |