CN107193048A - 一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 - Google Patents
一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 Download PDFInfo
- Publication number
- CN107193048A CN107193048A CN201710328475.7A CN201710328475A CN107193048A CN 107193048 A CN107193048 A CN 107193048A CN 201710328475 A CN201710328475 A CN 201710328475A CN 107193048 A CN107193048 A CN 107193048A
- Authority
- CN
- China
- Prior art keywords
- mtd
- msub
- mrow
- mtr
- resistance
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法，包括如下步骤：(1)通过现场试验获取隧道背景场以及金属体与探测线圈不同耦合情况下的瞬变电磁响应值；(2)分别计算步骤(1)中两种瞬变电磁响应的电阻极限值，并进一步求取金属体与探测线圈不同耦合情况下的电阻极限函数；(3)求取各测点低阻地质异常体瞬变电磁响应的电阻极限值；(4)以低阻地质异常体电阻极限值为特征值建立反演目标函数，求解计算低阻地质异常体的埋深及大小。本发明能够有效剔除金属体对瞬变电磁法探测结果的干扰，减少误判、提高隧道地质预报精度。
Description
技术领域
本发明涉及隧道空间瞬变电磁法探测领域，尤其涉及一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法。
背景技术
隧道开挖时掌子面前方经常会遇见富水断层、充水溶洞等低阻不良地质体，它们是隧道施工过程中重要安全隐患。成功预测掌子面前方低阻不良地质体成为保障隧道安全施工的重要措施，为此人们将瞬变电磁法引入到隧道超前地质预报中。然而隧道空间探测环境复杂，往往存在大量金属干扰体，已有研究表明这些金属体对探测数据及解译结果有较大的影响，有时甚至会导致误判。
目前，已有部分学者针对矿井瞬变电磁法金属干扰体进行了研究，但已有方法现场应用时仍存在局限性，如：对于给定的金属体和探测装置，金属体的瞬变电磁响应时间谱还取决于其与线圈的耦合情况；当有多个低阻体时，各个低阻体之间存在互感作用，瞬变电磁响应值并非单个低阻体的简单加减；屏蔽探头只适合中心回线，对地质异常响应较强的重叠回线不适合，且由于磁场线为闭合曲线，接收磁探头对后方信号的屏蔽效果还有待进一步验证。故需要一套有效的金属干扰消除方法，分离出低阻异常体的响应信号，准确反演低阻地质异常体位置及大小，对提高瞬变电磁法隧道超前地质预报精度具有十分重要的意义。
发明内容
本发明针对上述现有技术的不足，利用不同导体电阻极限值可简单加减的特点，提供了一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法，从而有效剔除金属体对瞬变电磁法探测结果的干扰，减少误判、提高隧道地质预报精度。
为实现上述目的，本发明的技术方案是：一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法，包括以下步骤：
(1)通过现场试验获取隧道背景场以及金属体与探测线圈不同耦合情况下的瞬变电磁响应值；
(2)分别计算步骤(1)中两种瞬变电磁响应的电阻极限值，并进一步求取金属体与探测线圈不同耦合情况下的电阻极限函数；
瞬变电磁响应的电阻极限值计算方法，包括：
首先按式(1)对瞬变电磁响应值进行拟合，计算各展开指数项的振幅系数Aj与时间常数τj，然后按式(2)计算各测点电阻极限值。
式中ti为不同采样时刻，Vi为不同时刻所采集的感应电压值；
式中，RL表示电阻极限值，j为所展开的项数；
所述展开指数项项数取10项，即可满足拟合精度要求；
所述电阻极限函数采用指数函数形式拟合求取；
(3)求取各测点低阻地质异常体瞬变电磁响应的电阻极限值；
(4)以低阻地质异常体电阻极限值为特征值建立反演目标函数，求解计算低阻地质异常体的埋深及大小，
所述反演目标函数是根据低阻地质异常体瞬变电磁响应的电阻极限解析式建立，其形式如下：
式中fk(r)为实际各测点低阻地质异常体的电阻极限值，mk(r,qi)为根据电阻极限解析式计算所得的电阻极限值，r为测点编号，N为测点数；qi为反演参数，包括低阻地质异常体埋深、大小，
所述的反演目标函数求解采用通用全局优化算法UGO求解。
所述步骤(3)中各测点低阻地质异常体瞬变电磁响应的电阻极限值是通过将各测点总电阻极限值减去金属干扰体及隧道地质背景场电阻极限值获得。
本发明的突出优点在于：
利用电阻极限的以下特点：①是在频率趋于零时定义；②导体之间不存在互感；③不同导体的电阻极限值可简单加减，实现金属体干扰信号分离。并将电阻极限值作为特征值进一步实现了低阻地质异常体埋深及大小的快速反演计算。对提高瞬变电磁法隧道超前地质预报精度具有十分重要的意义。
附图说明
图1为本发明所述的隧道超前地质预报瞬变电磁响应的金属干扰消除方法的总流程图。
图2为本发明所述的应用实施例中感应电动势多测道剖面图。
图3为本发明所述的应用实施例中低阻地质异常体电阻极限值曲线。
图4为本发明所述的应用实施例中低阻地质异常体埋深及大小的反演结果。
具体实施方式
下面结合附图和应用实施例对本发明做进一步说明。
如图1所示，本发明所述的隧道超前地质预报瞬变电磁响应的金属干扰消除方法，具体操作步骤如下：
(1)通过现场试验获取隧道背景场以及金属体与探测线圈不同耦合情况下的瞬变电磁响应值；其中隧道背景场是指无金属体及低阻地质异常体时的瞬变电磁响应值，不同耦合情况具体指金属体与探测线圈平面不同距离、不同角度。
(2)分别计算步骤(1)中两种瞬变电磁响应的电阻极限值，并进一步求取金属体与探测线圈不同耦合情况下的电阻极限函数；
瞬变电磁响应的电阻极限值计算方法，包括：
首先按式(1)对瞬变电磁响应值进行拟合，计算各展开指数项的振幅系数Aj与时间常数τj，然后按式(2)计算各测点电阻极限值。
式中ti为不同采样时刻，Vi为不同时刻所采集的感应电压值；
式中，RL表示电阻极限值，j为所展开的项数；
所述展开指数项项数取10项，即可满足拟合精度要求；
所述电阻极限函数采用指数函数形式拟合求取；
(3)求取各测点低阻地质异常体瞬变电磁响应的电阻极限值；
(4)以低阻地质异常体电阻极限值为特征值建立反演目标函数，求解计算低阻地质异常体的埋深及大小，
所述反演目标函数是根据低阻地质异常体瞬变电磁响应的电阻极限解析式建立，其形式如下：
式中fk(r)为实际各测点低阻地质异常体的电阻极限值，mk(r,qi)为根据电阻极限解析式计算所得的电阻极限值，r为测点编号，N为测点数；qi为反演参数，包括低阻地质异常体埋深、大小，
所述的反演目标函数求解采用通用全局优化算法UGO求解。
所述的电阻极限定义为：
在瞬变电磁响应晚期阶段，二次感应电流完全扩散到了导体内部，涡流分布趋于稳定，磁场按指数规律衰减。此时，任意局部导体的瞬变电磁响应可展开为多个按指数规律衰减的响应项总和，其频率域表达式分别为：
式中，j为所展开的项数；τj为时间常数，取决于导体自身的形状、大小与电导率；Aj为振幅系数，与导体的埋深、倾角、观测点位置有关。
从而引入参数电阻极限概念，在频率域，电阻极限定义为频率趋于零时瞬变电磁响应对频率的导数。其表达式为：
所述的实测数据电阻极限值计算过程为：
瞬变电磁仪测量数据，即瞬变电磁响应，常为感应电动势与时间的关系曲线，根据已知的一系列V和t值，按(6)式对瞬变电磁数据进行拟合，以计算各展开指数项的振幅系数Aj与时间常数τj，从而进一步按(5)计算各测点瞬变电磁响应的电阻极限值。
式中ti为不同采样时刻，Vi为不同时刻所采集的感应电压值。显然，展开项数越多，拟合精度就越高，实际拟合时一般展开10项，即可满足精度要求。为减少一次场及晚期电磁场波动对拟合结果的影响，拟合计算前，需根据曲线衰减特征选择合适的时间段。
所述的各测点低阻地质异常体瞬变电磁响应的电阻极限值是通过将各测点总电阻极限值减去金属干扰体及隧道地质背景场电阻极限值获得。
图2为本发明所述的应用实施例中感应电动势多测道剖面图，数据来源于广西岑溪市岭脚隧道左洞掌子面FK18+218处。测点间距为0.5m，线圈距隧道边墙大于1m，金属台车距掌子面20m。
在本实例中采用高阻围岩中非磁性导电球体模型建立反演目标函数，其电阻极限解析表达式为：
其中A(l)为无量纲的振幅，表达式为：
F(x0,l)为几何坐标函数，表达式为：
其中：
式中，X0为球心在线圈平面投影到回线中心的距离，x0＝X0/h；l＝L/h；r＝R/h；L为回线半边长，R为球体半径，h为球体埋深。
图3为采用本发明所述的方法，剔除金属干扰与背景场后，计算所得各测点低阻地质异常体电阻极限值。
图4为低阻地质异常体埋深及大小的反演结果，反演计算所得低阻异常体半径为5.7米、埋深11.7米，低阻体中心靠近隧道右边墙。实际开挖至掌子面FK18+212.5处时遭遇泥浆充填型空腔，空腔为椭圆状，长半径约6.1米，采用管棚支护得以顺利通过，反演结果与实际开挖情况较为吻合。
Claims (2)
1.一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法，其特征在于，该方法包括以下步骤：
(1)通过现场试验获取隧道背景场以及金属体与探测线圈不同耦合情况下的瞬变电磁响应值；
(2)分别计算步骤(1)中两种瞬变电磁响应的电阻极限值，并进一步求取金属体与探测线圈不同耦合情况下的电阻极限函数；
瞬变电磁响应的电阻极限值计算方法，包括：
首先按式(1)对瞬变电磁响应值进行拟合，计算各展开指数项的振幅系数Aj与时间常数τj，然后按式(2)计算各测点电阻极限值，
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mi>m</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>t</mi>
<mi>n</mi>
</msub>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mi>m</mi>
</msub>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>m</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中ti为不同采样时刻，Vi为不同时刻所采集的感应电压值；
<mrow>
<mi>R</mi>
<mi>L</mi>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<msub>
<mi>A</mi>
<mi>j</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，RL表示电阻极限值，j为所展开的项数；
所述展开指数项项数取10项，即可满足拟合精度要求；
所述电阻极限函数采用指数函数形式拟合求取；
(3)求取各测点低阻地质异常体瞬变电磁响应的电阻极限值；
(4)以低阻地质异常体电阻极限值为特征值建立反演目标函数，求解计算低阻地质异常体的埋深及大小，
所述反演目标函数是根据低阻地质异常体瞬变电磁响应的电阻极限解析式建立，其形式如下：
<mrow>
<mi>&Phi;</mi>
<mrow>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>/</mo>
<mi>N</mi>
<mo>&CenterDot;</mo>
<munderover>
<mo>&Sigma;</mo>
<mn>0</mn>
<mi>N</mi>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>f</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>m</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中fk(r)为实际各测点低阻地质异常体的电阻极限值，mk(r,qi)为根据电阻极限解析式计算所得的电阻极限值，r为测点编号，N为测点数；qi为反演参数，包括低阻地质异常体埋深、大小，
所述的反演目标函数求解采用通用全局优化算法UGO求解。
2.根据权利要求1所述的一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法，其特征在于，所述步骤(3)中各测点低阻地质异常体瞬变电磁响应的电阻极限值是通过将各测点总电阻极限值减去金属干扰体及隧道地质背景场电阻极限值获得。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710328475.7A CN107193048B (zh) | 2017-05-11 | 2017-05-11 | 一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710328475.7A CN107193048B (zh) | 2017-05-11 | 2017-05-11 | 一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107193048A true CN107193048A (zh) | 2017-09-22 |
| CN107193048B CN107193048B (zh) | 2018-11-30 |
Family
ID=59872982
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710328475.7A Active CN107193048B (zh) | 2017-05-11 | 2017-05-11 | 一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107193048B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649881A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 一种各向异性介质瞬变电磁场的数值模拟方法及系统 |
| CN113987877A (zh) * | 2021-10-27 | 2022-01-28 | 南通大学 | 一种钢筋网干扰源瞬变电磁响应的有限元仿真简化方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6335624B1 (en) * | 1999-06-02 | 2002-01-01 | Commissariat A L'energie Atomique | Process for electromagnetic detection of conducting objects |
| CN103760614A (zh) * | 2014-02-24 | 2014-04-30 | 中国科学院电子学研究所 | 一种适用于不规则发射波形的瞬变电磁正演方法 |
-
2017
- 2017-05-11 CN CN201710328475.7A patent/CN107193048B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6335624B1 (en) * | 1999-06-02 | 2002-01-01 | Commissariat A L'energie Atomique | Process for electromagnetic detection of conducting objects |
| CN103760614A (zh) * | 2014-02-24 | 2014-04-30 | 中国科学院电子学研究所 | 一种适用于不规则发射波形的瞬变电磁正演方法 |
Non-Patent Citations (2)
| Title |
|---|
| 于景邨 等: "巷道内金属体瞬变电磁响应特征及处理技术", 《煤炭学报》 * |
| 刘宗辉 等: "隧道空间内金属体对瞬变电磁法探测结果影响研究", 《现代隧道技术》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112649881A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 一种各向异性介质瞬变电磁场的数值模拟方法及系统 |
| CN113987877A (zh) * | 2021-10-27 | 2022-01-28 | 南通大学 | 一种钢筋网干扰源瞬变电磁响应的有限元仿真简化方法 |
| CN113987877B (zh) * | 2021-10-27 | 2024-08-09 | 南通大学 | 一种钢筋网干扰源瞬变电磁响应的有限元仿真简化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107193048B (zh) | 2018-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Clark | Methods for determining remanent and total magnetisations of magnetic sources–a review | |
| Wang et al. | A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data | |
| CN104267442A (zh) | 一种用于煤矿井下的瞬变电磁拟地震探测方法 | |
| CN104614774B (zh) | 一种瞬变电磁探测方法、装置和系统 | |
| CN102495430B (zh) | 深埋隧道岩体完整性的探测方法 | |
| CN105510981A (zh) | 一种磁铁矿采空区地球物理判定方法和装置 | |
| CN114236621A (zh) | 一种基于多技术协同的堤防隐患无损探测方法 | |
| CN203930090U (zh) | 三分量tem井中磁探头 | |
| CN114265124B (zh) | 基于时间域瞬变电磁概率反演的不良地质体定位方法 | |
| SONG et al. | Research on transient electromagnetic response of magnetic source in borehole | |
| CN103852795A (zh) | 一种水下小目标的磁异常信号提取方法 | |
| CN111538093A (zh) | 一种用于浅层地表探测方法及瞬变电磁仪器 | |
| Ramazi et al. | Contribution of geophysical inversion theory and geostatistical simulation to determine geoelectrical anomalies | |
| CN104615843A (zh) | 一种识别隐伏火山岩型铀矿化信息的物化探综合方法 | |
| CN107193048B (zh) | 一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 | |
| Geng et al. | Identifying void defects behind Tunnel composite lining based on transient electromagnetic radar method | |
| CN105842738B (zh) | 一种地下目标体异常定量判定方法和装置 | |
| CN103196991A (zh) | 连续诊断管体金属腐蚀与缺陷的全覆盖瞬变电磁检测方法 | |
| CN104391335A (zh) | 瞬变电磁纯异常信号提取方法 | |
| CN104793268B (zh) | 一种瞬变电磁探测的盲深度测量方法和装置 | |
| CN106054257A (zh) | 一种含磁层状介质电磁响应的确定系统 | |
| Xiong et al. | Automatic defect detection in operational high-speed railway tunnels guided by train-mounted ground penetrating radar data | |
| Li et al. | The analysis of the early electromagnetic response of the receiving coil and its application at close-range TEM detection | |
| Thabit et al. | Detection of subsurface cavities by using pole-dipole array (Bristow's method)/hit area-western Iraq | |
| CN106610509A (zh) | 一种瞬变电磁资料的时域处理方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |