CN104391335A - 瞬变电磁纯异常信号提取方法 - Google Patents
瞬变电磁纯异常信号提取方法 Download PDFInfo
- Publication number
- CN104391335A CN104391335A CN201410746340.9A CN201410746340A CN104391335A CN 104391335 A CN104391335 A CN 104391335A CN 201410746340 A CN201410746340 A CN 201410746340A CN 104391335 A CN104391335 A CN 104391335A
- Authority
- CN
- China
- Prior art keywords
- field
- induced voltage
- voltage data
- measuring point
- data
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种瞬变电磁纯异常信号提取方法,包括如下步骤:确定背景场感应电压数据;确定测点感应电压数据;将测线上每一个测点感应电压数据与背景场感应电压数据进行反卷积处理,剔除二次场背景场数据,得到不含二次场背景场的二次场纯异常场数据;依据二次场纯异常场数据确定地下地质目标体性质。本发明瞬变电磁纯异常信号提取方法,通过将地下介质产生的测点感应电压数据中的二次场背景场剔除,即提取出二次场的异常场部分,然后根据二次场的异常场进行资料判定,对地下目标性质解释的精准度大幅度提高。
Description
技术领域
本发明涉及一种地质信息提取方法,具体的说,是涉及一种瞬变电磁纯异常信号提取方法。
背景技术
瞬变电磁场法(Transient Electromagnetic Method,简称TEM)是一种建立在电磁感应原理基础上的时间域人工源电磁探测方法。它是利用阶跃形波电磁脉冲激发,利用不接地回线向地下发射一次场,在一次场关断后,在测点测量由地下介质产生的感应二次场随时间的变化,来达到寻找各种地质目标体的一种地球物理信息提取方法。
瞬变电磁法大多采用空芯线圈和有芯线圈形式,用以测量与磁场随时间变化率dB/dt成正比的二次感应电压归一化参数V(t)/I,并根据以下公式由二次感觉电压计算视电阻率ρτ(t):
其中μ0为真空中的磁导率,M为发送回线的磁矩,q为接收线圈的有效面积。进一步根据装置参数换算出深度数据后,可形成视电阻率一深度剖面。这种解释方法仍处于定性和“看图识字”的阶段。对瞬变电磁法测深资料定量解释还局限于单点一维反演,很多情况下是靠解释人员的工作经验对测深结果做出大致判断,人为性较大。对地下目标体的评价精度低。远不能满足实际工程地质的需要,这一状况也限制了方法的推广应用。
发明内容
针对上述现有瞬变电磁法解释技术中的不足,本发明提供一种能够精确确定地下目标性质的瞬变电磁纯异常信号提取方法。
本发明所采取的技术方案是:
一种瞬变电磁纯异常信号提取方法,包括如下步骤:
确定背景场感应电压数据;
对需要进行探测的地下地质目标体,利用不接地回线向地下发射一次场;在一次场关断后,测量由地下介质产生的测点感应电压数据;
将测线上每一个测点感应电压数据与背景场感应电压数据进行反卷积处理,剔除二次场背景场数据,得到不含二次场背景场的二次场纯异常场数据;
依据二次场纯异常场数据确定地下目标性质。
所述确定背景场感应电压数据的步骤为:选取一条没有异常体存在的剖面,利用不接地回线向地下发射一次场;在一次场关断后,在测点测量由地下介质产生的感应电压数据;将感应电压数据作为二次场背景场感应电压数据。
所述确定背景场感应电压数据的步骤为:任意选取一条剖面,利用不接地回线向地下发射一次场;在一次场关断后,在测点测量由地下介质产生的感应电压数据;选取10-99个测点,将所有测点的感应电压数据的平均值作为二次场背景场感应电压数据,将各测点的观测数据作为测点感应电压数据。
所述测点感应电压数据与背景场感应电压数据进行反卷积处理过程为:将时间域测点感应电压数据与背景场感应电压数据通过快速傅里叶变换转换到频率域。
所述测点感应电压数据与背景场感应电压数据进行反卷积处理过程为:将频率域二次场纯异常场响应转换到时间域;得到二次场纯异常数据。
本发明相对现有技术的有益效果:
本发明瞬变电磁纯异常信号提取方法,通过将地下介质产生的测点感应电压数据中的二次场背景场剔除,即提取出二次场的异常场部分,然后根据二次场的异常场进行资料判定,对地下目标性质解释的精准度大幅度提高。
附图说明
图1是本发明瞬变电磁纯异常信号提取方法的二次场背景场响应图;
图2是本发明瞬变电磁纯异常信号提取方法的二次场异常场响应图;
图3是本发明瞬变电磁纯异常信号提取方法的归一化纯二次异常场响应对数坐标图;
图4是本发明瞬变电磁纯异常信号提取方法的实测数据二次场背景场响应图;
图5是本发明瞬变电磁纯异常信号提取方法的实测数据二次场异常场响应图;
图6是本发明瞬变电磁纯异常信号提取方法的归一化二次场纯异常场图。
具体实施方式
以下参照附图及实施例对本发明进行详细的说明:
图1、2、4和5坐标系说明:横坐标表示时间(单位ms),纵坐标表示接收到的电压响应幅值(单位V);
由于数值模拟响应采用单位偶极源,其响应幅度与实测响应存在区别。
图3和6坐标系说明:横坐标表示时间(单位ms),纵坐标表示归一化二次场纯异常响应(无量纲)。
附图1-6可知,一种瞬变电磁纯异常信号提取方法,包括如下步骤:
确定背景场感应电压数据;
对需要进行探测的地下地质目标体,利用不接地回线向地下发射一次场;在一次场关断后,测量由地下介质产生的测点感应电压数据;
将测线上每一个测点感应电压数据与背景场感应电压数据进行反卷积处理,剔除二次场背景场数据,得到不含二次场背景场的二次场纯异常场数据;
依据二次场纯异常场数据确定地下目标性质。
所述确定背景场感应电压数据的步骤为:选取一条没有异常体存在的剖面,利用不接地回线向地下发射一次场;在一次场关断后,在测点测量由地下介质产生的感应电压数据;将感应电压数据作为二次场背景场感应电压数据。
所述确定背景场感应电压数据的步骤为:任意选取一条剖面,利用不接地回线向地下发射一次场;在一次场关断后,在测点测量由地下介质产生的感应电压数据;选取10-99个测点,将所有测点的感应电压数据的平均值作为二次场背景场感应电压数据
所述测点感应电压数据与背景场感应电压数据进行反卷积处理过程为:将时间域测点感应电压数据与背景场感应电压数据通过快速傅里叶变换转换到频率域。由于测点感应电压数据等于背景场感应电压数据与二次场纯异常场响应的卷积,时间域的卷积等价于频率域的乘积,将提取二次场纯异常场的反卷积过程转换到频率域进行计算后,它等价于频率域测点感应电压数据除以频率域背景场感应电压数据。
所述测点感应电压数据与背景场感应电压数据进行反卷积处理过程为:将频率域二次场纯异常场响应转换到时间域;得到二次场纯异常数据。
本发明瞬变电磁纯异常信号提取方法,通过将地下介质产生的测点感应电压数据中的二次场背景场剔除,即提取出二次场的异常场部分,然后根据二次场的异常场进行资料判定,对地下目标性质解释的精准度大幅度提高。
瞬变电磁观测的是二次场,二次场信号包括二次场背景场和二次场异常场,将二次场背景场剔除出去,即提取出二次场的异常场部分,然后根据二次场的异常场进行资料解释,那么,解释的精准度会大大提高。
将所接收到的大地二次场响应看作是一个时变系统,系统输出的二次场总响应包括系统背景二次场正常响应和目标体的纯二次场异常响应,通过反卷积的办法提取来自目标体的二次场纯异常响应信号。具体过程如下:
根据线性时变系统特性,输出信号可表示为:
ak(xr,yr,t)=s(t)*g(xr,yr,t)+n(xr,yr,t) (2)
式中,(xr,yr)表示测点坐标,t表示观测时间,ak(xr,yr,t)表示所接收到的总二次场响应,
s(t)表示与大地有关的二次场背景正常场响应
g(xr,yr,t)表示目标体的二次场纯异常响应
n(xr,yr,t)表示噪声。
提取二次场的异常场步骤如下:
通过反卷积获得自目标体的纯二次场响应,以得到的目标体纯二次场响应剖面比传统的视电阻率剖面更容易反映出地下目标体。
具体过程如下:
(1)将等式(2)两端进行傅里叶变换。根据卷积定理,时间域的卷积等价于频率域的乘积,从而卷积在频率域的表达式为:
A(xr,yr,ω)=S(xr,yr,ω)G(xr,yr,ω)+N(xr,yr,ω) (3)
(2)将等式两端同时除以背景正常场响应的频谱S(xs,xr,ω),即乘以乘积因子1/S(xs,xr,ω)。由于噪声项的存在,在背景正常场响应的某些频点幅度很低时,噪声项将被放大从而干扰反卷积的结果。为此在除以背景正常场响应的频谱S(xs,xr,ω)时引入常数ε,常数ε的存在避免了频域反卷积结果出现偏差较大的频点。引入常数因子很之后,原乘积因子转换转换为:
其中,ε表示值很小的正常数,*表示复数的共轭。该因子的分母为实数,常数ε的存在使得当|S(xr,yr,ω)|很小时,|S(xr,yr,ω)|作为分母引入的不稳定性得到了消除。
(3)上述反卷积过程中并没有考虑观测系统的有限带宽性质,因此在等式两端同时乘以如下零相位高斯滤波器以适应有限带宽系统:
式中ω0为D(ω)衰减为1/e时的角频率值,可以得到如下等式:
F(ω)D(ω)V(ω)=F(ω)D(ω)S(ω)G(ω)+F(ω)D(ω)N(ω)
≈D(ω)G(ω)+F(ω)D(ω)N(ω)
(6)
零相位滤波器的时间域形式为:
实际上这是对delt脉冲的有限带宽的近似,其在零时刻有极大值,且其在整个时间轴上的积分为1。
(4)将等式(6)两端进行整理可以得到经过整理,可以得到频率域的二次场纯异常响应,然后通过反傅里叶变换即可得到时间域二次场纯异常响应,其表达式为:
g(xs,xr,t)=f(t)*d(t)*a(xs,xr,t)-f(t)*n(xr,t) (8)
二次场纯异常响应的提取模拟:
设背景模型为电阻率为100Ω·m的均匀半空间,其瞬变电磁响应设定为背景场异常。背景场响应曲线如图1所示。
异常模型为在深度范围500~700m赋存有电阻率为1Ω.m后的低阻层,其二次场异常场响应如图2所示。图2异常场曲线与图1背景场曲线形状基本上很相似,如果仅依靠图2的曲线(二次场异常场曲线)是无法明确判断是否有异常体存在。所以需要进行纯二次场异常场数据信息提取。
采用上述反卷积的方法,利用公式(8),得到的纯二次异常场响应曲线图3,从图3可以看出,曲线的早时间和晚时间道的值都比较小,说明其对应二次场的背景场,而在中间时间道的响应值比较大,明显对应二次场纯异常。如果图3曲线比较平缓,说明不存在纯二次场,即地下不存在目标体;如果图3曲线峰值交小,说明纯二次异常场较弱,即地下目标体规模较小或者埋藏深度较大;如果图3曲线峰值较大大,说明纯二次异常场较强,即地下目标体规模较大或者埋藏深度较浅。
本发明野外观测数据处理实例:
利用瞬变电磁法对河北省隆化县某金矿进行探测,二次场背景场响应实测区线为图4。二次场异常场响应(实测数据)为图5。图5异常场曲线与图4背景场曲线形状基本上很相似,如果仅依靠图2的曲线(二次场异常场曲线)是无法明确判断是否有异常体存在。
采用上述反卷积的方法,得到的归一化纯二次异常场响应图6。从图6可以直接判断矿体异常情况,在图6中的晚期时间道上,出现两个负尖脉冲。说明在这两处位置上可能存在矿体异常,这一判断与实际相符。
Claims (5)
1.一种瞬变电磁纯异常信号提取方法,其特征在于,包括如下步骤:
确定背景场感应电压数据;
对需要进行探测的地下地质目标体,利用不接地回线向地下发射一次场;在一次场关断后,测量由地下介质产生的测点感应电压数据;
将测线上每一个测点感应电压数据与背景场感应电压数据进行反卷积处理,剔除二次场背景场数据,得到不含二次场背景场的二次场纯异常场数据;
依据二次场纯异常场数据确定地下目标性质。
2.根据权利要求1所述瞬变电磁纯异常信号提取方法,其特征在于:所述确定背景场感应电压数据的步骤为:选取一条没有异常体存在的剖面,利用不接地回线向地下发射一次场;在一次场关断后,在测点测量由地下介质产生的感应电压数据;将感应电压数据作为二次场背景场感应电压数据。
3.根据权利要求1所述瞬变电磁纯异常信号提取方法,其特征在于:所述确定背景场感应电压数据的步骤为:任意选取一条剖面,利用不接地回线向地下发射一次场;在一次场关断后,在测点测量由地下介质产生的感应电压数据;选取10-99个测点,将所有测点的感应电压数据的平均值作为二次场背景场感应电压数据,将各测点的观测数据作为测点感应电压数据。
4.根据权利要求1所述瞬变电磁纯异常信号提取方法,其特征在于:所述测点感应电压数据与背景场感应电压数据进行反卷积处理过程为:将时间域测点感应电压数据与背景场感应电压数据通过快速傅里叶变换转换到频率域。
5.根据权利要求1所述瞬变电磁纯异常信号提取方法,其特征在于:所述测点感应电压数据与背景场感应电压数据进行反卷积处理过程为:将频率域二次场纯异常场响应转换到时间域;得到二次场纯异常数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410746340.9A CN104391335A (zh) | 2014-12-10 | 2014-12-10 | 瞬变电磁纯异常信号提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410746340.9A CN104391335A (zh) | 2014-12-10 | 2014-12-10 | 瞬变电磁纯异常信号提取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104391335A true CN104391335A (zh) | 2015-03-04 |
Family
ID=52609259
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410746340.9A Pending CN104391335A (zh) | 2014-12-10 | 2014-12-10 | 瞬变电磁纯异常信号提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104391335A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105277993A (zh) * | 2015-11-21 | 2016-01-27 | 山东省第一地质矿产勘查院 | 一种海域金矿找矿方法 |
CN105824052A (zh) * | 2016-02-25 | 2016-08-03 | 西安西北有色物化探总队有限公司 | 一种地下目标体瞬变电磁多分量合成方法和装置 |
CN105842738A (zh) * | 2015-12-25 | 2016-08-10 | 中国科学院地质与地球物理研究所 | 一种地下目标体异常定量判定方法和装置 |
CN106405663A (zh) * | 2016-10-09 | 2017-02-15 | 煤炭科学技术研究院有限公司 | 一种基于背景场反馈的矿井瞬变电磁数据采集装置及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040232917A1 (en) * | 2001-09-07 | 2004-11-25 | Wright David A | Detection of subsurface resistivity contrasts with application to location of fluids |
CN101176017A (zh) * | 2005-03-14 | 2008-05-07 | 睦泰姆有限责任公司 | 真实幅度瞬时电磁系统响应测量 |
US20090309599A1 (en) * | 2008-06-11 | 2009-12-17 | Antoni Marjan Ziolkowski | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
CA2689870A1 (en) * | 2009-01-27 | 2010-07-27 | Mtem Ltd. | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
-
2014
- 2014-12-10 CN CN201410746340.9A patent/CN104391335A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040232917A1 (en) * | 2001-09-07 | 2004-11-25 | Wright David A | Detection of subsurface resistivity contrasts with application to location of fluids |
CN101176017A (zh) * | 2005-03-14 | 2008-05-07 | 睦泰姆有限责任公司 | 真实幅度瞬时电磁系统响应测量 |
US20090309599A1 (en) * | 2008-06-11 | 2009-12-17 | Antoni Marjan Ziolkowski | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
CA2689870A1 (en) * | 2009-01-27 | 2010-07-27 | Mtem Ltd. | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
Non-Patent Citations (6)
Title |
---|
ANTON ZIOLKOWSKI,ET AL.: "Multitransient electromagnetic demonstration survey in France", 《GEOPHYSICS》 * |
ANTON ZIOLKOWSKI,ET AL.: "Multi-transient electromagnetic repeatability experiment over the North Sea Harding field", 《GEOPHYSICAL PROSPECTING》 * |
CHRIS ANDERSON,ET AL.: "Multi-transient EM technology in practice", 《FIRST BREAK》 * |
DAVID A.WRIGHT,ET AL.: "Hydrocarbon detection with a multi-channel transient electromagnetic survey", 《SEG INT’L EXPOSITION AND ANNUAL MEETING》 * |
DAVID WRIGHT,ET AL.: "Hydrocarbon detection and monitoring with a multicomponent transient electromagnetic (MTEM) survey", 《THE LEADING EDGE》 * |
蒋帮远: "《实用近区磁源瞬变电磁法勘探》", 31 December 1998 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105277993A (zh) * | 2015-11-21 | 2016-01-27 | 山东省第一地质矿产勘查院 | 一种海域金矿找矿方法 |
CN105842738A (zh) * | 2015-12-25 | 2016-08-10 | 中国科学院地质与地球物理研究所 | 一种地下目标体异常定量判定方法和装置 |
CN105842738B (zh) * | 2015-12-25 | 2018-05-04 | 中国科学院地质与地球物理研究所 | 一种地下目标体异常定量判定方法和装置 |
CN105824052A (zh) * | 2016-02-25 | 2016-08-03 | 西安西北有色物化探总队有限公司 | 一种地下目标体瞬变电磁多分量合成方法和装置 |
CN105824052B (zh) * | 2016-02-25 | 2019-01-18 | 西安西北有色物化探总队有限公司 | 一种地下目标体的瞬变电磁多分量合成方法和装置 |
CN106405663A (zh) * | 2016-10-09 | 2017-02-15 | 煤炭科学技术研究院有限公司 | 一种基于背景场反馈的矿井瞬变电磁数据采集装置及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Salem et al. | Analytic signal approach and its applicability in environmental magnetic investigations | |
Smith et al. | Using an induction coil sensor to indirectly measure the B-field response in the bandwidth of the transient electromagnetic method | |
CN104360401B (zh) | 一种瞬变电磁b场确定地下目标体地质信息方法 | |
Podgorski et al. | Processing and inversion of commercial helicopter time-domain electromagnetic data for environmental assessments and geologic and hydrologic mapping | |
CN103995301A (zh) | 一种评价页岩气储层中总有机碳含量的方法及装置 | |
CN104391335A (zh) | 瞬变电磁纯异常信号提取方法 | |
CN109541695B (zh) | 人工场源频率域电场梯度远区视电阻率快速成像方法 | |
Prouty et al. | MetalMapper: A multi-sensor TEM system for UXO detection and classification | |
CN108169802B (zh) | 一种粗糙介质模型的时域电磁数据慢扩散成像方法 | |
JP2939334B2 (ja) | 可聴周波以下の低周波磁気測定装置 | |
CN104793268B (zh) | 一种瞬变电磁探测的盲深度测量方法和装置 | |
Kang et al. | Forward modeling and analyzing for frequency domain semi-airborne EM method | |
RU2733095C2 (ru) | Способ поиска трехмерных объектов методами геоэлектрики тм-поляризации | |
He et al. | A fusion approach for suppression of environmental noise in spread spectrum induced polarization data | |
CN107942373A (zh) | 基于裂缝性油气储层断裂系统检测的相干算法 | |
CN106610509A (zh) | 一种瞬变电磁资料的时域处理方法 | |
CN107193048A (zh) | 一种隧道超前地质预报瞬变电磁响应的金属干扰消除方法 | |
Wang et al. | Characterization of excitation effects and data interpretation of combined time-domain multiwaveform transmission currents | |
CN105629314A (zh) | 一种瞬变电磁发射线圈的自感信号的剔除方法和装置 | |
CN104536062B (zh) | 一种基于地下结构反演的自动拾取磁性层埋深的方法 | |
AU2022202145B2 (en) | Comprehensive detection method for coal mine spontaneous combustion fire zone | |
Prospect | Report on a Helicopter-Borne Versatile Time Domain Electromagnetic (VTEM Plus) and Horizontal Magnetic Gradiometer Geophysical Survey | |
Zhou et al. | Divergence of tipper real induction vector in tensor frequency-domain ground-airborne electromagnetic method | |
Casteel et al. | Estimating subsurface permeability with 3d seismic attributes: A neural net approach | |
Su et al. | Analysis of response characteristics of ground-well transient electromagnetic method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information |
Inventor after: Xue Guoqiang Inventor after: Li Hai Inventor after: Di Qingyun Inventor after: Chen Weiying Inventor after: Zhong Huasen Inventor before: Xue Guoqiang Inventor before: Li Hai Inventor before: Chen Weiying Inventor before: Zhong Huasen |
|
COR | Change of bibliographic data | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20150304 |
|
RJ01 | Rejection of invention patent application after publication |