E-mail Alert Rss
 

物探与化探, 2019, 43(1): 143-147 doi: 10.11720/wtyht.2019.1159

方法研究·信息处理·仪器研制

地—井瞬变电磁关键技术问题研究

杜庆丰, 冯晓兰, 黄跃

中国地质科学院 地球物理地球化学勘查研究所,河北 廊坊 065000

Research on key technical problems of surface-borehole TEM

DU Qing-Feng, FENG Xiao-Lan, HUANG Yue

Institute of Geophysical and Geochemical Exploration, CAGS,Langfang 065000,China

责任编辑: 沈效群

收稿日期: 2018-04-20   修回日期: 2018-07-11   网络出版日期: 2019-02-20

基金资助: 国土资源公益性行业科研专项经费项目“深井地—井TEM系统实用化与推广应用”.  201511031

Received: 2018-04-20   Revised: 2018-07-11   Online: 2019-02-20

作者简介 About authors

杜庆丰(1963-),男,教授级高级工程师,主要从事瞬变电磁法研究工作。Email:duqingfeng@igge.cn

摘要

地—井瞬变电磁法在加拿大、澳大利亚等国的勘查实践中已取得了很好的效果,国内由于缺乏从仪器设备到资料处理解释的整套技术,使得该方法未能被广泛应用。为此,研制了地—井瞬变电磁仪器,同时系统地研究了地—井三分量坐标系的确定、探头角度系统偏差的测定、三分量的校正运算、勘查目标体方位及距离的分析与计算等,并在实际的地—井瞬变电磁观测中获得了良好的应用效果。

关键词: 地—井瞬变电磁法; ; 探头角度偏差测定 ; 三分量校正 ; 目标体定位

Abstract

Surface-borehole TEM method has achieved good results in exploration in Canada and Australia. Due to the lack of complete technology from equipment to data processing and interpretation, the method has not been widely used in China. The authors developed the surface-borehole TEM instrument and, at the same time, systematically studied the determination of the three component coordinates of the surface-borehole TEM, the measurement of the deviation of the probe angle, the correction of the three components, and the analysis and calculation of the orientation and distance of the target body, which achieved good application results in actual measurement.

Keywords: surface-borehole TEM ; measuring the deviation of the angle of the probe ; three component correction ; target body location

PDF (0KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

杜庆丰, 冯晓兰, 黄跃. 地—井瞬变电磁关键技术问题研究. 物探与化探[J], 2019, 43(1): 143-147 doi:10.11720/wtyht.2019.1159

DU Qing-Feng, FENG Xiao-Lan, HUANG Yue. Research on key technical problems of surface-borehole TEM. Geophysical and Geochemical Exploration[J], 2019, 43(1): 143-147 doi:10.11720/wtyht.2019.1159

0 引言

地—井瞬变电磁法是在钻孔上方或钻孔附近的地面上发射,用接收探头在钻孔中逐点测量发射电流关断后地下地质体中感应的瞬变电磁二次场,它不是测井方法,而是一种地—井物探方法。用该方法测量时,接收探头与勘查对象距离较近,受导电覆盖层和电磁干扰小,故较地面TEM更易发现深部的隐伏矿体;该方法可在钻孔周围一定范围内,最大限度地发现井旁及井底矿体,还可对钻孔已见的良导目标体进行定位并追踪其延伸方向。

在国外,地—井TEM方法的研究始于20世纪70年代,80年代获得快速发展,90年代初期三分量仪器便得以商品化,目前该方法在加拿大、澳大利亚等国已被广泛应用。近二十多年来,国内在地—井TEM方法上亦开展了较多的研究工作,研究内容主要集中在方法理论方面,国内地—井TEM工作所使用的仪器设备均需从国外进口,但进口设备的高价格以及下井设备的高风险,制约了该项技术的开展和应用。为此,依托国土资源部公益性行业科研专项经费项目,我们从仪器设备、测量方法技术、处理解释软件等方面进行了系统的研究,希望能推动该方法技术在国内走向实用化。

1 IGGETEM4.0三分量地—井TEM系统

研制的地—井TEM系统包括:适用于井深不超过2 000 m的三分量TEM井下探头,三分量TEM接收机,大功率TEM发射机,5 kW整流器(图1),1 000 m和2 000 m塑胶电缆专用绞车(图2)以及配套的数据处理软件。

图1

图1   三分量TEM整流器、发射机、接收机及探头

Fig.1   Rectifier, transmitter, receiver and probes of the thre-component TEM system


图2

图2   1 000 m、2 000 m塑胶电缆专用绞车

Fig.2   Special winches for 1 000 m and 2 000 m plastic cables


2 三分量坐标系的确定

地—井TEM坐标系目前尚不统一,本文定义的坐标系如图3所示。如此定义坐标系的理由有三:x为北、y为东,这样可使其与大地坐标的北(x)、东(y)一致,亦即与地形图的坐标一致;z为下,与瞬变电磁教科书的多数示意图中的磁场方向一致;xyz三个分量的正方向适用于右手定则(右手笛卡尔坐标系)。

图3

图3   地—井TEM坐标系

Fig.3   Coordinate system of surface-borehole TEM


3 探头角度系统偏差的测定

TEM探头在井中是可以自由转动的,即探头中的xy线圈的方向是不可控的,探头在井中各测点的转角、倾角及方位角等参数由安装在同一探头中的定向器给出。在安装探头中的TEM线圈和定向器时,难以使得TEM线圈的北(x)、东(y)与定向器相应的北、东轴完全一致,为此,每个探头成为产品之前均须通过实际测量,测定出TEM线圈的xy线圈的法线方向与定向器北、东轴的角度偏差。

为测得这些角度数据,制作了专门装置。装置上包含一个与该装置固定在一起的可发射磁场的线圈,当该线圈的角度被确定时,通过旋转探头使TEM线圈接收的值达到最大或最小的方式确定出TEM线圈的角度,而此时同时读取测向器的重、磁参数并计算出相应角度,TEM线圈的角度与测向器相应角度之差便是我们需要的角度偏差值。

TEM线圈的北、东、下各轴与测向器相应轴的夹角在下文中计作:φxφyφz。角度偏差测定的具体装置情况及具体测量过程可参见参考文献[3]

4 三分量的校正运算

4.1 井中定向器角度的计算

各厂家的定向器其坐标系定义不尽相同,由定向器测量结果计算的滚角、倾角、方位角的表达式便不同。文中用图3中的地—井TEM坐标系来定义定向器的xyz坐标,在此坐标系下,整理的各角度的计算公式如下。

1)由重力参数计算重力滚角θ:

cosθ=gy(gx2+gy2)1/2,sinθ=gx(gx2+gy2)1/2,tanθ=gxgy

2)由磁性参数计算磁滚角θm:

cosθm=-Hy(Hx2+Hy2)1/2,sinθm=-Hx(Hx2+Hy2)1/2,tanθm=HxHy

3)由重力参数计算倾角ε:

g=(gx2+gy2+gz2)1/2,cosε=gzg,sinε=(gx2+gy2)1/2g,tanε=(gx2+gy2)1/2gz

4)由重力、磁性参数计算的方位角ϕ:

Hx'= Hx·gy+Hy·gx(gx2+gy2)1/2,

Hy'= Hx·gx+Hy·gy+Hz·gzg,

Hz'= Hz·(gx2+gy2)-Hx·gx·gz-Hy·gz·gyg·(gx2+gy2)1/2,

cosϕ= Hz'(Hz'2+Hx'2)1/2,

sinϕ= -Hx'(Hz'2+Hx'2)1/2,

tanϕ= -Hx'Hz'

上述各式中,HxHyHzgxgygz分别为测向器测得的磁性及重力参数;gHx'Hy'Hz'Vx'Vy'Vz'Vzy'等为中间变量。定向器各角度的定义见图4

图4

图4   定向器各角度示意

Fig.4   Schematic diagrams of the orientation angles


4.2 TEM三分量的校正计算

TEM探头测得的分别是xyz各线圈法线方向上的感应电动势的值,为将全井的观测结果校正至图3的地—井TEM坐标系,推导出校正计算公式。

1)倾角ε<3.0度时的校正计算:

Vxy0=( Vx02+ Vy02)1/2 ,

Vz=Vz0·cos(ε-φz)+Vxy0·sin(ε-φz) ,

Vx=Vx0·cos(θm-φx)-Vy0·sin(θm-φy) ,

Vy=Vx0·sin(θm-φx)+Vy0·cos(θm-φy) 。

2)倾角ε>=3.0度时的校正计算:

Vx'=Vx0·sin(θ-φx)+Vy0·cos(θ-φy) ,

Vy'=Vx0·cos(θ-φx)-Vy0·sin(θ-φy) ,

Vzy'=Vz0·sin(ε-φz)- Vy'·cos(ε-φz) ,

Vx= Vx'·cosϕ+ Vzy'·sinϕ ,

Vy= Vx'·sinϕ- Vzy'·cosϕ ,

Vz=Vz0·cos(ε-φz)+ Vy'·sin(ε-φz) 。

上述各式中:Vx0Vy0Vz0为实测的感应电动势数据;φxφyφz分别为TEM线圈的北、东、下轴至测向器相应轴的夹角(逆时针为正,顺时针为负); Vx'Vy'Vzy'Vxy0为中间变量;VxVyVz为经三分量校正计算后的感应电动势结果数据。三分量校正计算的具体计算过程参见文献[4]

5 数据处理

地—井瞬变电磁数据处理内容包括:数据编辑、滤波,三分量校正运算,异常识别及分析处理。

6 应用效果

运用自行研制的IGGETEM4.0地—井瞬变电磁仪器系统,在新疆哈密白石泉506钻孔中开展了测量工作,下述三组图分别为xyz三个分量的感应电动势剖面。

图5图6图7的三组图所示的感应电动势数据质量良好,除进行了必要的三分量校正运算外未作圆滑、滤波等处理,即这三组图是由原始数据直接绘制而成的。图中主要异常段位于320~390 m,通过对图件的定性分析及使用我们研制的处理软件计算,对该井的异常分析结果如下。

图5

图5   白石泉506钻孔x分量早、中、晚期感应电动势剖面

Fig.5   X-component profiles of induced electromotive force in early,middle and late stages of Baishiquan 506 borehole


图6

图6   白石泉506钻孔y分量早、中、晚期感应电动势剖面

Fig.6   Y-component profiles of induced electromotive force in early, middle and late stages of Baishiquan 506 borehole


图7

图7   白石泉506钻孔z分量早、中、晚期感应电动势剖面

Fig.7   Z-component profiles of induced electromotive force in early, middle and late stages of Baishiquan 506 borehole


z方向:在360 m深度处,早期为弱的正异常,钻孔在该处穿过低阻体(边缘);晚期为负异常,表明钻孔旁侧有低阻体,由异常衬度可见,该低阻体规模较大。

xy方向:由早期xy感应电动势剖面曲线可见,钻孔穿过低阻体,该低阻体向南倾斜,其主体位于钻孔西侧偏北;由晚期xy感应电动势剖面曲线可见,钻孔旁侧有较大规模的低阻体,该低阻体向北倾斜,其主体位于钻孔东侧偏南。

xyz的中期道,由于钻孔所穿过低阻体异常(正异常)与钻孔东侧较远处的低阻体的异常(负异常)反向,两种异常在中期道相互抵消,故中期道感应电动势相对早、晚期而言异常不够明显。

7 结论与建议

经实际勘查工作检验,这套地—井瞬变电磁法系统性能稳定,测量结果是正确的;经多次测试,开发的三分量校正计算软件,计算结果正确,研制的数据处理程序可基本满足勘查工作的需要。

这套自主研发的地—井瞬变电磁法系统已基本接近实用,并初步形成为一套较完整的应用技术,建议对该系统做进一步的完善并加大推广使用的力度,以使我国可更快、更好、更广泛地开展地—井瞬变电磁法的勘查工作。

致谢:在地—井瞬变电磁法系统的实用化应用中,得到了新疆地质矿产勘查开发局第六地质大队、青海省第三地质矿产勘查院、内蒙古国土资源勘查开发院的大力支持,在此表示感谢。

The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。

参考文献

蒋邦远 . 实用近区磁源瞬变电磁法勘探[M]. 北京: 地质出版社, 1998.

Jiang B Y. Applid near zone magnetic source transient electromagnetic exploration[M]. Beijing: Geological Publishing House, 1998.

杜庆丰 .

瞬变电磁 (TEM)数据组合滤波

[J]. 地质与勘探, 2002,38(s):129-131.

DOI:10.3969/j.issn.0495-5331.2002.z1.036     

随着TEM方法理论的不断进步,其正反演计算及解释在理论上已基本成熟,为取得良好的勘测效果,如何从观测数据中去除干扰,获得良好、可靠的原始信息变得至关重要,故就此问题进行探讨.

Du Q F .

TEM data composite filter

[J]. Geology and Prospecting, 2002,38(s):129-131.

黄跃, 沈钟, 冯晓兰 , .

TEM井中磁探头的校准装置及校准方法:发明专利证书号第2231480号

[P].2016-09-07.

[本文引用: 1]

Huang Y, Shen Z, Feng X L , et al. Calibration device and calibration method of surface-borehole TEM magnetic probe: Invention Patent Certificate No.2231480[P].2016 -09-07.

[本文引用: 1]

中国地质科学院地球物理地球化学勘查研究所

.地—井瞬变电磁(TEM)三分量校正软件1.0:软著登字第0663591号[CP].2013-12-26.

[本文引用: 1]

Institute of Geophysical and Geochemical Exploration

Surface-borehole TEM three-component correction software 1.0: Computer software copyright registration No.0663591 [CP].2013 -12-26.

[本文引用: 1]

张杰, 邓晓红, 谭捍东 , .

地—井瞬变电磁资料矢量交会解释方法

[J]. 物探与化探, 2015,39(3):572-579.

DOI:10.11720/wtyht.2015.3.23      URL     Magsci    

<p>以往地&mdash;井瞬变电磁法资料解释研究主要为定性解释。本文在瞬变电磁等效涡流的理论基础上,通过研究等效场矢量的空间特征,发现感应二次场矢量的方向与一次场强度、观测点距离及观测时间均无关,只与观测点的方位有关,经简单的修正,二次场特征矢量就可以交会于异常中心位置。研究的地&mdash;井瞬变电磁矢量交会法是一种新的定量解释方法,该方法通过绘制每个测点的特征矢量,根据异常区矢量交会于一点的规律,确定异常中心位置。分别采用理论模型和实测资料对该方法进行了验证,结果表明矢量交会法能够对井旁盲矿中心位置进行较准确地定位。</p>

Zhang J, Deng X H, Tan H D , et al.

A study of vector intersection for borehole transient electromagnetic method

[J]. Geophysical and Geochemical Exploration, 2015,39(3):572-579.

Magsci    

武军杰, 李貅, 智庆全 , .

电性源地—井瞬变电磁法三分量响应特征分析

[J]. 地球物理学进展, 2017,32(3).

URL    

地-井瞬变电磁法是在地面发射,井(钻孔)中接收的装置形式,能够利用已有钻孔使接收探头深入地下更加接近矿体,获得更加可靠的目标体信息.文中针对复杂地形地区矿区深部找矿面临的实际问题,提出了电性源地-井瞬变电磁法的思路,并以FDTD正演模拟来研究均匀半空间瞬变场的空间分布以及一维、三维模型钻孔中的三分量瞬变响应特征,以期获得理论模型三分量响应特征规律,为实际勘查提供参考.理论模型计算结果表明,电性源地-井TEM方法三分量瞬变响应曲线对于电性界面、异常体具有良好的反应.

Wu J J, Li X, Zhi Q Q , et al.

Analysis of three component TEM response characteristic of electric source dill hole TEM

[J]. Progress in Geophysics (in Chinese), 2017,32(3).

孟庆鑫, 潘和平, 牛峥 .

大地介质影响下地—井瞬变电磁的正演模拟分析

[J]. 中国矿业大学学报, 2014,43(6).

URL    

对大地介质影响下地-井瞬变电磁响应特征规律进行了研究.基于时域有限差分法和Mur吸收边界条件对不同观测方式('定源-异井位观测'和'方位观测')下典型地电模型(均质半空间和含浅部低阻层的半空间介质)中三维体的地-井瞬变电磁响应进行了正演模拟,认为观测所得响应是围岩背景场与目标体异常场共同作用的结果,提出大地介质会影响观测结果的观点.研究结果表明:围岩背景场越强,总响应的异常特征越弱;随着观测延时的增加,背景场对总响应的影响减弱;受导电介质影响的观测曲线在剖面形态和衰减特征上与高阻围岩的典型解释模型存在差别.研究结果为相关资料解释工作提供了理论依据.

Meng Q X, Pan H P, Niu Z .

Forward simulation of surface-borehole TEM in geological medium effect

[J]. Journal of China University of Mining & Technology, 2014,43(6).

李建慧, 刘树才, 焦险峰 , .

地—井瞬变电磁法三维正演研究

[J]. 石油地球物理勘探, 2015,50(3):556-564.

DOI:10.13810/j.cnki.issn.1000-7210.2015.03.025      URL     Magsci    

地&mdash;井瞬变电磁法是一种使用较多的瞬变电磁法工作装置。鉴于现今国内外对其正演研究还比较少, 在此实现了基于一维解析法和矢量有限单元法的地&mdash;井瞬变电磁法一维和三维正演。通过对比层状介质的一维解析解和水平板状体地电模型的数值解, 验证了矢量有限单元法用于模拟地&mdash;井瞬变电磁法的有效性; 通过对低阻块状体地电模型进行三维正演, 分析总结了矩形回线源激发瞬变电磁场随时间和深度变化的扩散规律。

Li J H, Liu S C, Jiao X F , et al.

Three-dimensional forward modeling for surface-borehole transient electromagnetic method

[J]. Oil Geophysical Prospecting, 2015,50(3):556-564.

Magsci    

Xue G Q, Qin K Z, Li X , et al.

Discovery of a largescale porphyry molybdenum deposit in tibet through a modified TEM Exploration Method

[J]. JEEG, 2012,17(1):19-25.

DOI:10.2113/JEEG17.1.19      URL    

uring the last decade,there have been many exploration achievements in the Tibetan Gangdese metallogenic bell.The Sharang area of the Tibetan region is covered by a lowtemperature mineralized alterable clay that is considered to be a low-grade ore.Although small intervals of rich molybdenum(Mo) mineralization have been discovered,the ore deposit scale is limited and the condition of deep ore is still unknown.To explore these deeper targets,a modified large-loop TEM system was used in the Sharang area.The TEM receiver configuration is redesigned and the late-time resistivity equation of large-loop TEM has also been defined.During data processing,two regions with low resistivity anomalies were discovered.The interpreted results indicate that the main ore deposit is buried 200 m beneath the surface,and extends 600 m vertically.The total anomalous area associated with the ore deposit is estimated at 3.77 km~2.The interpretation results are consistent with drilling data acquired after the geophysical survey.The results show that it is the first ultra-large porphyry molybdenum deposit that has been found in Tibet.

Guo J L, Zeng Y Q, Li X.

Numerical simulation analysis of surface-to-borehole TEM based on the finite difference method

[C]// International Workshop and Gravity,Electrical & Magnetic Methods and their Applications, 2015: 115-118.

Swidinsky A, Nabighian M .

On smoke rings produced by a loop buried in a conductive half-space

[J]. Geophysics, 2015,80(4):E225-E236.

DOI:10.1190/geo2015-0002.1      URL    

The transient electromagnetic (TEM) response of a wire loop at the surface of the earth has been studied in detail for several decades. When current in the loop is rapidly switched off, an identical system of current, commonly referred to as a moke ring, is induced in the conductive subsurface and diffuses downward and outward into the earth. What has not been studied as thoroughly is the response of a loop buried within a conductive medium, as might be the case in the marine environment or when TEM prospecting systems are placed in mine tunnels or horizontal boreholes. We have examined the TEM response of a horizontal wire loop buried within a conductive half-space. Expressions were derived for the azimuthal electric field and corresponding vertical and radial magnetic fields. Our results showed that when the loop is far from the earth-air interface, a single smoke ring system diffuses radially outward from the transmitter, while the electric field and corresponding current density decays away in the vertical direction. As the loop approaches the interface, the smoke ring system diffuses radially at early times, but gradually, the complex image of the loop in the air produces a system of secondary azimuthal electric fields, which, when combined with the primary field, adds a vertical component to the field diffusion. At late times, the field behavior reduces to the well-known surface case and the maximum current system diffuses downward at a constant angle of 26 with respect to the plane of the loop. We concluded that it was the effect of the interface that produces the downward migration of a smoke ring system, whereas the outward migration is mainly generated by the primary field.

Annetts D.

Inversion of surface and downhole electromagnetic data for a 3D earth

[C]// 22nd ASEG Geophysical Conference, Extended Abstracts, 2012: 1-3.

Buselli G, Lee S K .

Modelling of drill-hole TEM response from multiple targets covered by a conductive overburden

[J]. Exploration Geophysics, 1996,27(3):141-153.

DOI:10.1071/eg996141      URL    

A number of transient electromagnetic (TEM) analogue modelling projects have been carried out at CSIRO to compile suites of model responses for surface, drill-hole and airborne surveys. The drill-hole modelling described in this paper has been carried out to determine the response of targets in free space and covered by overburden. The objective of this modelling was specifically to investigate the effect of inductive interaction between a conductive overburden and electrically-isolated targets beneath it, and to examine the effect of current channelled from the overburden into a target in electrical contact with it. For the experiments discussed here, a model consisting of two parallel plate targets with a dip angle of 45 and covered by a conductive overburden was employed. Measurements were made both with and without the overburden present. The modelling results show that for cases of simple conductors without overburden, the characteristics of the drill-hole TEM responses can be understood in terms of the geometrical relationship between the transmitter loop, the drill-hole receiver probe, and the conductive targets. However, the modelling experiments demonstrate that the presence of a conductive overburden can significantly alter the drill hole response. Inductive interaction between the overburden and targets can be large enough to cause polarity reversals of the drill-hole response of the targets. Additionally, when the upper target is in electrical contact with the overburden, current channelling produces a galvanic current that may suppress or enhance peaks in the drill-hole profiles. When interpreting field data with the use of superposition of responses from various components of a given model, unexpected polarity reversals and enhanced or suppressed peak responses in the field data could be explained by these two effects.

Gao J, Zhao A B, Peng F , et al.

Inversion of array induction logs and its application

[J]. Petroleum Science, 2007,4(3):31-45.

DOI:10.1007/s12182-007-0005-x      URL    

With the help of the modified geometrical factor theory, the Marquardt method was used to calculate the true electrical parameters of the formation from array induction logs. The inversion results derived from the assumed model and some practical cases show that the rebuilt formation profile determined by 2-ft resolution array induction logs is reasonable when the formation thickness is greater than 1 m, which thus indicates that the inversion method is reliable and can provide quantitative information for the discrimination of oil/gas or water zone.

Hughes N A, Ravenhurst W R .

Three componet DHEM surveying at balcooma

[J]. Exploration Geophysics, 1996,27:77-89.

DOI:10.1071/EG996077      URL    

http://library.seg.org/doi/pdf/10.1071/EG996077

/

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址:北京市学院路29号航遥中心 邮编:100083
电话:010-62060192;62060193 E-mail:whtbjb@sina.com