瞬变电磁法三分量测量方法研究

席振铢，刘剑，龙霞，侯海涛

(中南大学 信息物理工程学院，湖南 长沙，410083)

摘  要：在不考虑本征电流之间互感的情况下，研究用本征电流代替涡流计算地表瞬变电磁场的模拟计算方法，计算均匀半空间中定源装置下薄板导体在不同倾角时的X，Y和Z 3个分量的响应特征，并与1条已知地质剖面的实测结果进行对比分析。研究结果表明：除板状体直立外，薄板导体X分量的过零点基本上是板状体的中心在地面的投影位置；当板状体水平或直立时，Y分量曲线是对称的，Z分量曲线是反对称的；当板状体倾斜时，Z分量和Y分量曲线都不对称，在倾斜方向曲线平缓，响应范围大，在反方向曲线陡立，响应范围小，因此，用X分量可以确定异常的中心，用Z和Y分量可以定性判断板状体的倾向；实测结果与计算结果一致，说明综合分析X，Y和­Z方向3个分量的方法比单分量分析的方法更有助于对异常体的解释。

关键词：瞬变电磁法；三分量；本征电流；薄板导体

中图分类号：P319.1+2          文献标志码：A         文章编号：1672-7207(2010)01-0272-05

Three-component measurement in transient electromagnetic method

XI Zhen-zhu, LIU Jian, LONG Xia, HOU Hai-tao

(School of Info-physics and Geomatics Engineering, Central South University, Changsha 410083, China)

Abstract: Without considering the mutual inductance between eigencurrent and eddy current, the eigencurrent which was instead of eddy current was used to calculate the fixed-loop X, Y and Z component transient responses over a thin plate conductor in homogeneous half-space with various plate dip angles. The results show that except for a vertical target the cross-over of X-component is approximately the projection position of the plate’s center. For a horizontal plate, both Y and Z-component plots are symmetrical. For a vertical target Y-component is symmetrical and Z-component is anti-symmetrical. While for an inclined plate, neither component is symmetrical, it is broad and gentle in the inclination direction but narrow and steep in the other direction. So target location and orientation can be defined qualitatively by three-component measurement. The observed results are coincident with the calculated ones, which proves that the method of unitarily analyzing three-component is better than that for the interpretation of conductor body.

Key words: transient electromagnetic method; three-component; eigencurrent; thin plate conductor

瞬变电磁法^[1-4](Transient electromagnetic method，简称TEM^[5-6])是属于时间域的电磁勘探方法，它是利用不接地回线或接地线源向地下发送一次脉冲磁  场，在一次脉冲磁场的间歇利用线圈或接地电极观测二次场的方法。人们对瞬变电磁法的研究大多是测量垂直分量并对其响应特征进行分析，但由于受仪器限制以及易受噪声等因素的影响导致对水平分量的研究还不够深入。Gallagher等^[7]计算了均匀半空间中导电薄板在不同倾角时垂直分量的偶极瞬变响应特征；Adhidjaja等^[8]计算了自由空间中大定源组合装置下垂直分量在导电薄板不同走向长度、不同电导率时的响应特征；Zonge等^[9]讨论了水平分量在探测地下管线中的作用；McNill等^[10-13]阐述了测量三分量的意义，但并没有进一步分析异常体产生的三分量响应特征和规律。由于水平分量提供了一些仅仅通过垂直分量无法得到的有关异常体的信息，因此，对水平分量和垂直分量的曲线响应特征及其规律进行综合分析很有必要。尽管均匀半空间是一种理想的地电结构模型，但对野外情况下地电结构进行分析仍具有借鉴意义。鉴于板状体是野外常见的地质模型，本文作者用本征电流代替涡流计算均匀半空间中薄板导体产生的磁场三分量曲线响应特征，并总结不同板状体参数下三分量曲线的响应规律，然后，结合1条已知地质剖面进行验证。

1  三分量响应模拟计算

1.1  计算方法

Gallagher等^[7]用15个本征电流环代替涡流来计算和模拟了薄板导体的电磁响应，这种计算方法的优点是计算方法简单，速度快而且准确。由于不考虑本征电流之间的互感，故薄板导体产生的瞬变响应为所有本征电流各自产生的响应的总和。在早期，所有的本征电流都会产生瞬变响应，即；在晚期，瞬变响应几乎都是由最内层的本征电流引起，即。其中：n为本征电流环的数量；t为测量时间；τ_i为第i个本征电流的时间常数；c_i为与板的几何参数和时间常数等相关的脉冲响应的幅度因子。本文作者引用上述算法，用Geonics公司的MOTEM^[13]程序和部分本文作者编制的Matlab程序计算了均匀半空间中定源装置下薄板导体产生的磁场三分量曲线响应特征。

1.2  模拟条件及计算结果

设均匀半空间电阻率为500 ??m，发送框面积  (长×宽)为300 m×600 m，以框中心为原点建立直角坐标系，发送框中心坐标用(x, y, z)表示。板走向为0?，沿走向长为400 m，向下延伸400 m，纵向电导为10 S，板顶部中心坐标用(x, y, z)表示。倾角沿X轴正方向顺时针方向旋转，Z方向以竖直向上为正，发送电流为15 A，t=277 ms，对薄板导体在不同倾角(0?，45?，90?和135?)下的三分量曲线响应特征进行模拟计算，结果见图1~4。其中：纵坐标dB/dt为磁感应强度的微分即本征电流。由图1~4可见：当板倾角为0?时，产生的3个分量的响应比其他倾角情况下所产生的相应分量的响应要强得多，且Z分量的响应最强，X分量的响应次之，Y分量的响应最小；板倾角为90?时产生的响应最小；当板水平时，X分量响应曲线过零点处对应板的中心在地面的投影位置，且以其过零点处呈反对称状分布，Y分量在板的上方表现为单峰异常，并以X分量曲线的过零点处呈对称状分布；Z分量在板的上方表现为1个明显的单峰异常极大值，并在靠近板边缘的处各出现1个较小的负异常极小值，其曲线形态也以X分量的过零点处呈对称状分布。当板直立时，X和Y分量在板体上方表现为单峰异常，Z分量的过零点位置对应板的顶部中心在地面的投影位置，且以其过零点处呈反对称分布；当板倾斜时，X分量曲线的过零点位置仍大致对应板的中心在地面的投影位置，Y分量曲线分布形态在倾斜侧变化较平缓，反倾斜侧变化较陡，据此可以判定板的倾向。虽然Z分量响应具有在板倾斜侧变化相对平缓，在反倾斜侧变化相对陡的特征，但不如Y分量对板体倾角变化的反映敏感，且曲线形态比X和Y分量的曲线形态复杂。

1—X分量；2—Y分量；3—Z分量

图1  板体倾角为0?时的三分量响应

Fig.1  Fixed-loop response of three-component calculated from eigencurrent over a thin plate conductor in homogeneous half-space when plate dip is 0?

1—X分量；2—Y分量；3—Z分量

图 2  板体倾角为45?时的三分量响应

Fig.2  Fixed-loop response of three-component calculated from eigencurrent over a thin plate conductor in homogeneous half-space when plate dip is 45?

1—X分量；2—Y分量；3—Z分量

图 3  板体倾角为90?时的三分量响应

Fig.3  Fixed-loop response of three-component calculated from eigencurrent over a thin plate conductor in homogeneous half-space when plate dip is 90?

1—X分量；2—Y分量；3—Z分量

图 4  板体倾角为135?的三分量响应

Fig.4  Fixed-loop response of three-component calculated from eigencurrent over a thin plate conductor in homogeneous half-space when plate dip is 135?

2  实例

2.1  地质及地球物理概况

铜山铜矿位于铜陵—安庆成矿带中部，其构造位置位于姥山背斜的南东翼，地层为泥盆系—二叠系地层。姥山背斜的总体走向为北东50?，在矿床产出部位轴向弯转，呈向南东突出的弧形。区内断裂构造发育，以平行褶皱轴向的逆冲断裂、层间伸展滑动引起的构造破碎带和北东向的压扭性断裂为主。花岗闪长斑岩沿上述断裂的交接处侵入，并沿断裂面和层间滑动面贯入，引起广泛的接触交代变质、热液蚀变及矿化作用。铜山铜矿现已发现的矿体主要赋存于铜山花岗闪长斑岩体的接触带和层间破碎带内，主矿体多呈似层状，小矿体多呈透镜状、扁豆状、囊状。沿走向及倾向常出现分支复合、膨胀收缩及尖灭再现现象。主矿体一般连续性较好^[16]。该区闪长岩电阻率为1×10³~1×10⁴ ??m，大理岩电阻率大于1×10⁴ ??m，角页岩电阻率小于1×10³ ??m，金属硫化矿电阻率低于1 k??m。由此可见，测区矿石的电性表现为低阻，围岩的电性恰好相反，表现为高阻。良好的电性条件为瞬变电磁法在该区进行找矿提供了有利的地球物理 条件。

2.2  工作装置

实验采用的仪器是由加拿大Geonics公司生产的PROTEM-67型瞬变电磁仪。采用大定源回线装置对已知剖面即13线进行框内观测。13线地质剖面见图5。测线方位为NW346?，发送电流为20 A，基频为6.25 Hz，采用30个时窗进行采样，高精度石英钟晶振同步，发送框面积(长×宽)为200 m×600 m，关断时间为245 ?s，接收线圈有效面积为200 m²，同时接收X，Y和Z 3个方向的二次场本征电流(即dB/dt)，所得结果见图6~8。

图5  13线地质剖面

Fig.5  Geology profile for line 13

1—采样时间道为5；2—采样时间道为10；3—采样时间道为15；4—采样时间道为2；5—采样时间道为25

图 6  13线X分量的响应

Fig.6  Response of X component for line 13

1—采样时间道为5；2—采样时间道为10；3—采样时间道为15；4—采样时间道为2；5—采样时间道为25

图 7  13线Y分量的响应特征

Fig.7  Response of Y component for line 13

1—采样时间道为5；2—采样时间道为10；3—采样时间道为15；4—采样时间道为2；5—采样时间道为25

图 8  13线Z分量的响应特征

Fig.8  Response of Z component for line 13

2.3  实测结果分析

由图6~8可见：在0~160 m之间由于地表有厚层灰岩覆盖，尽管地下约600 m以下有矿体分布，但由于规模不大且埋深较大，3个分量响应均较弱，但在160~340 m之间X，Y和Z 3个分量的响应曲线均表现出较强的异常。X分量的响应特征表明：在点位260 m处附近出现的过零点大致对应矿体的中心位置，且由于过零点左侧矿体厚度较小且埋深较大，导致过零点左侧的瞬变响应小于右侧响应。由于矿体形态可以近似地看作板状体，其倾向也与图7中Y分量的响应规律一致，特别是在晚期表现得尤为明显，即在倾斜侧相对平缓，反倾斜侧陡立；Z分量的响应虽然较X分量和Y分量都要强，分辨异常的能力也较强，但是，图8所示曲线仅能说明在160~340 m段存在异常，并不能较准确地给出异常体中心和倾向的信息。因此，如果仅测量Z分量并只对其进行分析就可能会丢失一些异常体的重要信息，而这些信息可以通过对X分量和Y分量进行分析获得。

3  结论

(1) 利用 X 分量的过零点位置可以大致确定板状体中心在地面的投影位置，利用 Y 分量可以判断板状体的倾向，而利用 Z 分量便于发现异常，因此，综合 X，Y和Z 3个分量的响应特征进行分析可以准确地对异常进行评价。

(2) 在测量垂直分量的同时，并不需要付出额外的工作量就可以同时获得其他2个水平分量的信息，从而为更深入地对这3个分量的响应特征及其规律进行研究提供了良好的条件。

(3) 在已知板状体上方定源回线三分量瞬变响应的实际观测结果验证了本文对三分量的响应特征及其规律的分析方法是正确的。因此，应注重将综合三分量瞬变响应特征进行分析的方法应用于找矿工作中。

致谢  中南大学牛之琏老师对本文提出了许多宝贵意见，在此深表谢意！

参考文献：

[1] 牛之琏. 时间域电磁法原理[M]. 长沙: 中南大学出版社, 2007.
NIU Zhi-lian. The theory of time domain electromagnetic method[M]. Changsha: Central South University Press, 2007.

[2] 蒋邦远. 实用近区磁源瞬变电磁法勘探[M]. 北京: 地质出版社, 1998.
JIANG Bang-yuan. Applied near zone magnetic source exploration of transient electromagnetic method[M]. Beijing: Geology Press, 1998.

[3] 考夫曼 A A, 凯勒 G V. 频率域与时间域电磁测深[M]. 王建谋, 译. 北京: 地质出版社, 1987.
Kaufman A A, Keller G V. Frequency and transient soundings[M]. WANG Jian-mou, trans. Beijing: Geology Press, 1987.

[4] 米萨克N.纳比吉安. 勘查地球物理电磁法[M]. 赵经祥, 译. 北京: 地质出版社, 1992.
Nabighian M N. Electromagnetic methods in applied geophysics[M]. ZHAO Jing-xiang, trans. Beijing: Geology Press, 1992.

[5] 李貅. 瞬变电磁测深的理论与应用[M]. 西安: 陕西科学技术出版社, 2002.
LI Xiu. The theory and application of transient electromagnetic[M]. Xi’an: Shanxi Science and Technology Press, 2002.

[6] 阎述, 陈明生, 傅君眉. 瞬变电磁场的时域数值分析[J]. 地球物理学报, 2002, 45(2): 275-283.
YAN Shu, CHEN Ming-sheng, FU Jun-mei. Direct time-domain numerical analysis of transient electromagnetic field[J]. Chinese J Geophys, 2002, 45(2): 275-283.

[7] Gallagher P R, Ward S H, Hohmann G W. A model study of a thin plate in free space for the EM37 transient electromagnetic system[J]. Geophysics, 1985, 50(6): 1002-1005.

[8] Adhidjaja J I, Hohmann G V, Oristaglio M L. Two-dimension time domain electromagnetic response[C]//Technical Program Abstracts and Biographies. 53rd Annual International Soc of Expl Geophys Meeting. Tuls, 1983: 623-626.

[9] Zonge K L, Carlson N. The utility of horizontal component measurements in random-walk TEM surveys[R]. Tucoson: Zonge Engineering & Research Organization Inc, 2002.

[10] McNill J D. Why should I measure three components in a time-domain electromagnetic survey[R]. Mississauga: Geonics Limited, 1992.

[11] 刘金涛, 顾汉明, 胡祥云. 瞬变电磁法三分量解释剖析[J]. 人民长江, 2008, 39: 114-116.
LIU Jin-tao, GU Han-ming, HU Xiang-yun. Analysis of three-component transient electromagnetic[J]. Yangtze River, 2008, 39: 114-116.

[12] McNill J D, Edwards R N, Levy G M. Approximate calculations of the transient electromagnetic response from buried conductors in a conductive half-space[J]. Geophysics, 1984, 49(7): 918-924.

[13] Chen C S, Chiu W H, Lin C R. Three component time-domain electromagnetic surveying: Modeling and data. Analysis[J]. Piers Online, 2008, 4(4): 475-480.

[14] Wang T L, Tripp A C, Hohmannt G W. Studying the TEM response of a 3-D conductor at a geological contact using the FDTD method[J]. Geophysics, 1995, 60(4): 1265-1269.

[15] Barnett C T. Simple inversion of time-domain electromagnetic data[J]. Geophysics, 1984, 49(7): 925-933.

[16] 席振铢. 安徽省铜陵市铜山铜矿接替资源勘查项目物探设计[R]. 长沙: 中南大学信息物理工程学院, 2006.
XI Zhen-zhu. Tongling City, Anhui Province, the geophysical survey design to replace resource of Tongshan copper project[R]. Changsha: School of Info-physics and Geomatics Engineering, Central South University, 2006.

收稿日期：2008-12-11；修回日期：2009-03-16

基金项目：国家“十一五”科技支撑计划项目(2006BAB01B07)

通信作者：席振铢(1966-)，男，内蒙古商都人，博士研究生，副教授，从事电磁法勘探理论与应用研究；电话：13873150690；E-mail: xizhenzhu@163.com

(编辑 陈灿华)