CSAMT中 Ex视电阻率与 Ex/Hy比值视电阻率的对比
万伟1,2, 唐新功1,2
1.油气资源与勘探技术教育部重点实验室(长江大学),湖北 武汉 430100
2.长江大学 地球物理与石油资源学院,湖北 武汉 430100
通讯作者: 唐新功(1968-),男,教授,博士生导师。E-mail:tangxg@yangtzeu.edu.cn

作者简介: 万伟(1989-),男,硕士研究生在读,研究方向为CSAMT正反演。E-mail:beyondwanwei@163.com

摘要

常规CSAMT勘探是利用电场与磁场的比值来计算卡尼亚视电阻率。通过均匀半空间及层状介质模型,对CSAMT各场分量进行了计算,对视电阻率定义方式进行了分析。研究结果表明:采用电场 Ex单分量定义的视电阻率,可以减轻比值视电阻率近区畸变严重的问题。实测资料分析结果进一步验证了 Ex分量视电阻率的这一优点。在野外CSAMT资料采集时不必拘泥比值视电阻率这一种定义方式,可以根据野外实际条件灵活地调整采集方式,提高工作质量和效率。

关键词: CSAMT; 电场 Ex分量; Ex/Hy比值视电阻率; 卡尼亚视电阻率; 水平层状介质模型; 近区畸变
中图分类号:P631 文献标志码:A 文章编号:1000-8918(2015)03-0567-05
A comparative analysis of Ex apparent resistivity and Ex/Hy ratio apparent resistivity in CSAMT
WAN Wei1,2, TANG Xin-Gong1,2
1. Key Laboratory of Exploration Technologies for Oil and Gas Resources, Yangtze University, Wuhan 430100, China
2. Department of Geophysics and Petroleum Resource, Yangtze University, Wuhan 430100, China
Abstract

Currently, the adoption of the ratio of electric field and magnetic field to calculate the Cagniard resistivity is widely used in traditional CSAMT exploration. However, an analysis of the electromagnetic field components and the apparent resistivity definition method in homogeneous half space and layered media indicates that the apparent resistivity defined by electric field Ex component can improve the problem of the near zone distortion of ratio apparent resistivity. Moreover, the advantage of apparent resistivity defined Ex component has been further verified by the measured field data. Therefore it's not necessary to rigidly stick to using the ratio apparent resistivity in the field CSAMT data acquisition, and the acquisition mode flexibility can be adjusted according to actual conditions so as to improve work quality and efficiency.

Keyword: CSAMT; electric field component Ex; Ex/Hy ratio apparent resistivity; comparison; Cagniard resistivity; horizontal layered media model; near zone distortion

在可控源音频大地电磁法(CSAMT)中, 目前仍广泛采用比值电阻率定义, 即通过测量两个相互正交的电场(Ex)和磁场(Hy)来得到视电阻率。比值视电阻率的优点是计算简便, 在波区能较客观地反映地电断面的垂向变化, 然而在非波区会产生严重畸变, 影响测深曲线的解释 13; 其次, 由于通常多个电道共用一个磁道, 一旦磁道的数据丢失, 或测量出现偏差, 或噪声严重, 会影响多个测点卡尼亚视电阻率的使用[4, 5]。虽然工作时可在不影响信号质量的前提下尽可能地增大收发距, 但部分频点仍然不可避免地会进入近区, 导致曲线畸变, 造成资料解释不准确。许多学者基于这个原因, 一方面研究针对近区畸变校正的全区视电阻率法[1, 6, 7], 另一方面开始研究比值视电阻率之外的单分量视电阻率定义方法。林长佑等[8]研究分析了电场和磁场中三种单分量视电阻率— 相位联合反演方法, 并论述了实际应用的效果; 汤井田等[9]探讨了CSAMT电场y方向定义的视电阻率具有计算简单、改善近区畸变的优点; 邱卫忠等[10]运用实测资料分析了CSAMT标量观测中各分量的特征, 认为单分量视电阻率与比值视电阻率在纵向上具有相同的分辨率; 陈明生等[11]认为在电偶源频率电磁测深中的Ex分量定义视电阻率分辨率强于比值视电阻率。这些学者从不同的角度积极探索了不同视电阻率的定义方式。笔者将从场分量的特点出发来讨论不同视电阻率的定义方式的差别。

通过半空间中电偶源激发的各个电磁场分量在近区和远区的近似表达式 1214不难发现:在远区, 相对于EzHxHy分量, ExEyHz分量与地层电性关系更为密切, 但是Hz分量衰减相对更快, 且在采集过程中易受干扰; 在近区, 除ExEy外, 其余分量均不反映地下电性信息。所以, 在实际测量时, 采集ExEy两个分量较为有利。考虑到CSAMT野外经常采用的是赤道偶极装置进行标量观测, 此次仅对层状介质模型中电偶源频率测深中的Ex分量定义的视电阻率与比值视电阻率定义进行对比研究。

1 CSAMT场分量及视电阻率定义

视电阻率是空间中介质真电阻率的复杂加权平均, 一般通过均匀半空间场等效复杂地电模型场的方式来定义。在准静态极限下, 水平谐变电偶极源在均匀半空间表面产生的电磁场分量为[13]

Ex= ρPE2πr2[eikr(1-ikr)+(3cos2θ -2)](1)

Hy= ρPE2πr2I1K1+ iPE8πriμ0wρ{(cos2θ -1)[I0K1-I1K0]}(2)

式中:r为收-发距, θ 为观测点与偶极源连线方向夹角(如图1 所示), PE为电偶极矩, k为波数, ρ 为半空间中的电阻率, I0、I1、K0、K1为以-ikr/2为虚宗量的第一类、第二类贝塞尔函数。

图1 电偶源及场分量示意

CSAMT观测区域一般限制在如图1所示的梯形或扇形测区内。电偶极源激发的电磁场由两部分组成:一部分在地面传播, 并折射穿到地下, 当|kr|≪1时, 这部分波影响较小, 可以忽略; 另一部分直接从地中传播过去并按指数衰减, 当|kr|≫1时, 进入远区, 可将电磁波视为平面波。当介质由均匀半空间变为层状介质时, 对应的式(1)、(2)改写为

Ex= PE2π(iω μ 0) 0mm+m1/R* J0(mr)dm+

1r0( m1k12R- 1m+m1/R* )J1(mr)dm]sin2θ +

PE2π(iω μ 0)[ 0mm1k12RJ0(mr)dm-

1r0( m1k12R- 1m+m1/R* )J1(mr)dm]cos2θ , (3)

Hy= PE2π[ 0m1R* mm+m1/R* J0(mr)dm-

1r0mm+m1/R* J1(mr)dm]sin2θ +

PE2π[ 1r0mm+m1/R* J1(mr)dm]cos2θ ; (4)

R* =coth m1H1+arcothm1m2·

coth m2H2++arcothmN-1mN, RN* =1, (5)

R=coth m1H1+arcothm1ρ1m2ρ2·

coth m2H2++arcothmN-1ρN-1mNρN, RN=1。 (6)

式中:mi= m2+ki2, m的物理意义是空间频率, 具有距离倒数的量纲, ki2=-iω μ /ρ i; RR* 为电性参数。

|kr|≫1时, 有如下近似关系

Exp=ρPE2πr3(3cos2θ-2), (7)Hyp=(1+i)PE4πr32ρωμ0(3cos2θ-2)(8)

ExpHyp可视为均匀半空间中的平面波场。

计算层状介质的平面波场, 只需将半空间中平面波场公式乘上一个与地层有关的因子。当场强与ρ 成正比时, 乘上1/ G02(G0为地层因子), 当场强与 ρ成正比时, 乘上1/ G014。对于两层介质, 有

G0(ρ 1, ρ 2, h1)= (1/ρ1+1/ρ2)-(1/ρ1-1/ρ2)e2ik1h1(1/ρ1+1/ρ2)+(1/ρ1-1/ρ2)e2ik1h1, (9)

ρ 1ρ 2分别为第一、二层的电阻率, h1为第一层的厚度。对于层介质, 只需要将 1/ρn-1变成 1/ρn-1G0(σ n-1, σ n, hn-1), 再代入式(9)即可递推得到。显然对于均匀半空间有G0=1。

由式(7)、(8)可分别得到各场分量及对应比值视电阻率(卡尼亚视电阻率)的定义式

ρ|Ex|=2πr3PE(1-3sin2θ)|Ex|, (10)ρ|Hy|=4ωμ0πr3PE(3cos2θ-2)2|Hy|2, (11)ρExHy=1ωμ0ExHy2(12)

2 Ex分量及视电阻率计算结果分析

赤道偶极装置(θ =π /2)是CSAMT最常采用的测量方式。均匀半空间作为最简单的模型, 可以直观地反映出各视电阻率定义方法的优劣。图2a分别为均匀半空间(ρ =100 Ω · m)中各视电阻率随频率变化的曲线, 图2b为各视电阻率值相对真电阻率的偏差δ 。为了方便比较, 均使用对数刻度计算, 计算式定义为

δ=lg(Ff)-lg(F)lg(F)×100%, (13)

式中, F表示真电阻率或平面波场幅值, Ff表示各视电阻率或全区场幅值。由图2a、b可以看出, 由磁场分量Hy定义的视电阻率偏差最大, 其次是卡尼亚阻抗视电阻率, 偏差最小的是用电场Ex定义的视电阻率。视电阻率曲线的形态与场分量直接相关, 可从场的幅值随频率的变化来探讨其对视电阻率曲线形态的影响。

图2 固定收发距下均匀半空间各视电阻率曲线及其分析结果

CSAMT视电阻率的定义是基于平面波的假设, 而在收发距r一定时, 电偶极源产生的电磁波在低频段通常并不是平面波而导致计算的比值视电阻率发生畸变。图2c给出了对数坐标系中平面波场与非平面波场的电场与磁场幅值随频率的变化规律, 可以看出磁场幅值变化跨度远大于电场, 使得在计算阻抗视电阻率时磁场的影响占有较大的比重; 在图2d中, 磁场偏离平面波场的幅度也远大于电场, 进一步导致阻抗(Ex/Hy)的偏离。

综合分析图2可见, 比值视电阻率在近区发生严重的畸变与磁场分量密切相关, 虽然电场与磁场都有非平面波效应, 但磁场的影响比电场大得多, 继而影响到比值视电阻率的形态。基于这个原因, 在实际测量时, 可以考虑不使用磁场分量, 这样做一方面可以在一定程度上降低近场效应的影响, 提高数据资料的可利用率, 另一方面也能减少野外工作量, 降低成本, 提高工作效率。

图3 三层介质中的各视电阻率曲线

下面继续以水平层状介质模型来验证Ex分量定义视电阻率的优点。图3为三层介质中Ex分量视电阻率、比值视电阻率及平面波场视电阻率(ρ s)的对比, 分别计算了基底为高阻和基底为低阻情况下的各视电阻率。由图可见, 以ρ s为参考, 在低频段(近区), 在各模型中Ex视电阻率均优于比值视电阻率, 即Ex视电阻率比比值视电阻率在近区发生严重畸变的情况要轻得多。这种情况对于基底为低阻模型时优势更为明显, Ex视电阻率整体形态上更加接近于平面波场视电阻率。所以层状模型的计算结果说明了Ex分量定义的视电阻率相对比值视电阻率有减轻近区畸变的优点。

3 实测资料验证

以山西大同帽帽山地区2013年CSAMT的L1线实测资料为例。野外采集方式如图1所示, 通过实测的电场Ex和磁场Hy分量计算了由电场Ex和比值定义的两种视电阻率, 图4给出了部分测点的计算结果。图中显示, 相对阻抗比值视电阻率, 电场Ex定义的视电阻率均有改善近区畸变的效果。

实测资料的计算结果进一步说明, 在改善近区畸变方面, 电场Ex定义视电阻率要优于比值视电阻率。因而, 在野外实际工作时可以只采集一个电场分量, 这样可以大大减少野外工作量, 提高工作效率。

图4 帽帽山地区CSAMT测量L1线部分测点的 ρ|Ex|ρ|Ex/Hy|曲线

4 结论

传统的CSAMT阻抗视电阻率曲线在近区畸变严重, 较大程度上是受到磁场Hy分量的影响, 因此在数据采集时可以考虑放弃测量或放弃使用磁场Hy分量, 即单独采集Ex分量。均匀半空间和层状模型的计算结果表明, 用Ex分量定义的视电阻率比比值视电阻率的近区畸变程度小, 曲线形态更加接近平面波场视电阻率, 能够更准确地反映地下介质的电性变化。野外实测资料的计算也表明了这一点。使用Ex分量定义的视电阻率突出的优点是使得野外数据采集可以仅观测一个电场分量, 从而大大减少野外工作量, 降低成本, 提高工作效率。

The authors have declared that no competing interests exist.

参考文献
[1] 汤井田, 何继善. 水平电偶源频率测深中全区视电阻率定义的新方法[J]. 地球物理学报, 1994, 37(4): 543-552. [本文引用:1]
[2] 汤井田, 何继善. 可控源音频大地电磁法及其应用[M]. 长沙: 中南大学出版社, 2005. [本文引用:1]
[3] 林威. CSAMT法过渡区电磁场的特征[J]. 物探与化探, 2009, 33(2): 148-150 [本文引用:1]
[4] 王坤坤, 廖全涛. 高压线对CSAMT探测结果的影响及数据处理[J]. 物探与化探, 2014, 38(5): 1051-1054. [本文引用:1]
[5] 杨建军, 钟仁. CSAMT中磁探头摆放位置对卡尼亚视电阻率的影响[J]. 物探与化探, 2013, 37(3): 521-523. [本文引用:1]
[6] 毛先近, 鲍光淑. 水平电偶源频率电磁测深全区视电阻率的直接算法[J]. 中南工业大学学报, 1996, 27(3): 253-256. [本文引用:1]
[7] 詹少全, 钱美平, 冯戋戋. CSAMT全区视电阻率电场正演迭代拟合近场校正方法[J]. 物探与化探, 2011, 35(5): 663-665. [本文引用:1]
[8] 林长佑, 武玉霞, 罗东山, . 频率测深中单分量视电阻率—相位联合反演[J]. 煤田地质与勘探, 1993, 21(3): 46-52. [本文引用:1]
[9] 汤井田, 周聪, 张林成. CSAMT电场 y方向视电阻率的定义及研究[J]. 吉林大学学报: 地球科学版, 2011, 41(2): 552-558. [本文引用:1]
[10] 邱卫忠, 闫述, 薛国强, . CSAMT的各分量在山地精细勘探中的作用[J]. 地球物理学进展, 2011, 26(2): 664-668. [本文引用:1]
[11] 陈明生, 阎述, 陶冬琴. 电偶源频率电磁测深中的 Ex分量[J]. 煤田地质与勘探, 1998, 26: 60-66. [本文引用:1]
[12] 朴华荣. 电磁法测深原理[M]. 北京: 地质出版社, 1990. [本文引用:1]
[13] 底青云, 王若. 可控源音频大地电磁数据正反演及方法应用[M]. 北京: 科学出版社, 2008: 9-27. [本文引用:1]
[14] 石昆法. 可控源音频大地电磁法理论与应用[M]. 北京: 科学出版社, 1999: 10-12. [本文引用:1]