A Trial for Introducing 6-Element Tensor Impedance in Magnetotelluric Method and Its Application
-
摘要: 在经典大地电磁(MT)理论中,张量阻抗[ Z ]定义为电场分量和磁场分量之间的线性关系.估算张量阻抗[ Z ]及和它有关的其他参数(例如视电阻率、相位、倾子等)是MT数据处理中的一个重要环节.引入了六元素张量阻抗[ R ]的全新概念,并开发了相应的处理方法.为检验本方法的特征和抗噪性能,对采集自朝鲜的MT野外资料进行了分析.分析结果表明在MT资料处理中新定义的六元素张量阻抗[ R ]比传统的四分量张量阻抗[ Z ]提高测深曲线的相干度至少0.1以上,并且改善了大地电磁资料处理的质量.Abstract: The conventional magnetotelluric method is theoretically based upon definition of impedance tensor[ Z ] and the relationship between the electrical and magnetic components Ex, Ey, Hx, Hy, Hz of MT field, the main purpose of MT data processing is to obtain inpedance tensor, apparent resistivity, phase, tipper, skewness and other parameters. In this paper, we propose the definition of 6-elements tensor impedance and briefly describe its some characteristics and determination techniques in the comparison with the former impedance tensor[ Z ]. Furthermore, we explain the necessity of the proposed method and demonstrate its applicability by some field tests conducted in Democratic People's Republic of Korea.
-
Key words:
- magnetotelluric exploration /
- impedance tensor /
- signal to noise ratio /
- geophysics
-
大地电磁测深法(简称MT法)是一种成本较低、探测深度较大的地球物理勘探方法(刘文才等, 2015).经过几十年的发展, MT法在石油勘探(Tripp, 2005; 曹中林等, 2006)、地热勘探(Muñoz, 2014; 李德威和王焰新,2015)、水文勘探(Auken et al., 2017)、海洋勘查(de Groot-Hedlin and Constable, 2004; Linde and Pedersen, 2004)等方面均获得了广泛应用.相应地,MT数据处理、正演和反演算法也引起了许多地球物理学家的兴趣(García and Jones, 2002; Yin, 2003; 胡祖志和胡祥云, 2005; 张莹和张胜业, 2005).然而由于MT法采用天然场源,采集到的数据很容易受到种类繁多、性质复杂的噪声影响, 因此提高信噪比(S/N)是MT法研究的一个重要课题(王家映和徐义贤, 1998).
实际上, MT法中观测的物理量一般是地面上的3个磁场分量和2个水平电场分量, 利用这些分量获得阻抗张量和倾子矢量的估计值, 并评价资料质量.获得可靠的阻抗和倾子数据是大地电磁数据处理的主要任务.除了自然噪声会引起测量数据中较大的误差棒外,铁路、变压器等工业设备也是不容忽视的电磁噪声源;前者的噪声在500 kHz~1 Hz频率范围内的影响最强, 后者对更低频有较大影响, 其底限依赖于当地的电性特征(Pedersen and Engels, 2005).
为了解决上述问题, 目前已提出了许多方法技术, 可以概括为两个部分:第一方面是改善观测方法及装备.其中效果最好的方法就是远点参考观测方式(Gamble et al., 1979), 该方法已在一些论文当中介绍, 并且诸多实例已证明它可以改善噪声严重的环境中采集到的MT资料质量.通过最近几十年的研究, 该方法已从近点参考改进到远点参考, 从两点同步测量发展到三点同步测量法(Ri et al., 2002).此外, 1~5 kHz频率域内AMT信号的能量可以随日随季节变化, 白天磁场的振幅经常低于磁场探头的噪声水平.针对这一问题, García and Jones(2005)在综合频率域地电-地电观测(T-T)和地电-地磁观测(T-MT)的基础上提出了一种新的观测处理方式, RMT(14~250 kHz)法和可控源MT(CSMT)(1~12 kHz)结合可应用于厚度几十米左右的粘土透镜下面砂岩中的地下水勘探.
第二方面是和改进观测数据的数学处理方法有关, 学者们也已提出了许多不同方法, 例如单点观测数据处理中利用权函数消除噪声(Yee et al., 1988), 时间域阻抗张量的自适应重建方法(Gómez-Treviño and Mondragón, 1995), 利用张量阻抗的实部和虚部计算视电阻率(Egbert and Livelybrooks, 1996), 张量阻抗的Robust估计方法(王书明和王家映, 2004)等.其中, Robust资料处理方法被认为是一种特别有效的资料处理方法.还有学者提出了利用高阶统计量MT信号的高斯性检验(Xu and Wang, 2000)和基于小波变换(Gupta and Choubey, 2015)的时变大地电磁信号频谱估计方法.
另一方面, Lilley(1998)提出了一种新方法, 他把张量阻抗分解为全局部分和局部部分, 此方法的目的并不是为了克服人文噪声, 而是为了减小由于三维局部近地表地质体引起的电流畸变效应.Groom and Bailey(1991)曾提出的阻抗张量分解也是和上述方法类似.
然而即使有上述各种途径,目前MT法的信噪比问题仍有待提高,这意味着还需要开发除改进观测方法和处理方法以外的解决方法.
从此问题出发, 笔者再次考虑MT法数据处理的基本方程组.经典MT法处理基于电场E与磁场H之间的线性方程组,
$$ {E_x} = {Z_{xx}}{H_x} + {Z_{xy}}{H_y}, $$ (1) $$ {E_y} = {Z_{yx}}{H_x} + {Z_{yy}}{H_y}, $$ (2) $$ {H_z} = {T_{zx}}{H_x} + {T_{zy}}{H_y}, $$ (3) 式中Ex、Ey、Hx、Hy、Hz分别表示大地电磁场的x, y方向电场分量和x、y、z方向磁场分量, 并且:
$$ {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{Z_{xx}}}&{{Z_{xy}}}\\ {{Z_{yx}}}&{{Z_{yy}}} \end{array}} \right], $$ (4) $$ {\boldsymbol{T}} = \left[ \begin{array}{l} {T_{xx}}\\ {T_{zy}} \end{array} \right], $$ (5) 分别为张量阻抗及倾子矢量(Xu, 1987; García and Jones, 2002).
需要指出的是, (1)~(5)式并不是表示大地电磁场中5个分量之间互相关系的唯一表达式. Berdichevsky(1968)详细解释了这些方程组的意义,他先把大地电磁场中的Ex、Ey、Hx、Hy、Hz分量表示为相互独立的线极化Hxi、Hyi的总和, 即:
$$ {E_x} = {a_{1x}}{H_y}^i + {a_{2x}}{H_x}^i, $$ (6) $$ {E_y} = {a_{1y}}{H_y}^i + {a_{2y}}{H_x}^i, $$ (7) $$ {E_z} = {a_{1z}}{H_y}^i + {a_{2z}}{H_x}^i, $$ (8) $$ {H_x} = {b_{1x}}{H_y}^i + {b_{2x}}{H_x}^i, $$ (9) $$ {H_y} = {b_{1y}}{H_y}^i + {b_{2y}}{H_x}^i, $$ (10) $$ {H_z} = {b_{1z}}{H_y}^i + {b_{2z}}{H_x}^i, $$ (11) 上述式中, a1x、a2x、…、b1z、b2z是复数系数.
从(9)、(10)式, 求解入射波振幅Hxi、Hyi后, 再分别代入到(6)、(7), 可以获得如下:
$$ \begin{array}{l} {E_x} = {Z_{xx}}{H_x} + {Z_{xy}}{H_y},\\ {E_y} = {Z_{yx}}{H_x} + {Z_{yy}}{H_y}, \end{array} $$ 式中:
$$ \begin{array}{c} {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{Z_{xx}}}&{{Z_{xy}}}\\ {{Z_{yx}}}&{{Z_{yy}}} \end{array}} \right] = \\ \frac{1}{{{D_1}}}\left[ {\begin{array}{*{20}{l}} {{a_{1x}}{b_{2y}} - {a_{2x}}{b_{1y}}}&{{a_{2x}}{b_{1x}} - {a_{1x}}{b_{2x}}}\\ {{a_{1y}}{b_{2y}} - {a_{2y}}{b_{1y}}}&{{a_{2y}}{b_{1x}}{a_{1y}}{b_{2x}}} \end{array}} \right], \end{array} $$ 此时,式中的D1该满足D1=b1xb2y-b2xb1y≠0.由此可得到形如1、2、4式的关系式.
从另一方面看, 如果我们假定的D2=b1yb2x-b2xb1y≠0条件, 那经过对(10)、(11)式进行类似的推导可获得如下方程组
$$ {E_x} = {Z_{xy}}{H_y} + {Z_{xz}}{H_z}, $$ (1') $$ {E_y} = {Z_{yy}}{H_y} + {Z_{yz}}{H_z}, $$ (2') $$ {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{Z_{xy}}}&{{Z_{xz}}}\\ {{Z_{yy}}}&{{Z_{yz}}} \end{array}} \right]. $$ (4') 同样地, 从(9)、(11)式出发,假定D3=b1xb2x-b2xb1z≠0的条件,又可以如下关于张量阻抗的定义:
$$ {E_x} = {Z_{xx}}{H_x} + {Z_{xz}}{H_z}, $$ (1") $$ {E_y} = {Z_{yx}}{H_x} + {Z_{yz}}{H_z}, $$ (2") $$ {\boldsymbol{Z}} = \left[ {\begin{array}{*{20}{c}} {{Z_{xx}}}&{{Z_{xz}}}\\ {{Z_{yx}}}&{{Z_{yz}}} \end{array}} \right], $$ (4") 从对上述(1)、(2)、(4)式, (1′)、(2′)、(4′)式, 和(1″)、(2″)、(4″)式分析可以看出,MT张量阻抗的定义可以有三种2×2阶的矩阵形式.换句话说,用2×2阶矩阵形式来表示张量阻抗并不是唯一的.
张量阻抗就是磁场分量转换到电场分量的一种变换矩阵, 显然更一般的形式是在(4)式中把张量阻抗表示为3×3阶矩阵,而非2×2阶矩阵.
最初将MT阻抗张量定义为4个元素是基于假定:大地电磁场5个分量中Hz分量是由于地下电阻率的水平不均匀性及其他噪声影响而产生的.这一假定在一维模型条件下有其合理性,然而在实际处理中忽视Hz分量对Ex、Ey分量的影响(即按照1、2、3式处理)可能会导致张量阻抗估计计算的误差.
基于上述考虑, 本文发展了一种新的理论——大地电磁六元素张量阻抗的资料处理方法.
1. 六元素张量阻抗的定义及其特征
任何电磁场中电场与磁场之间的矢量关系可以用如下线性方程组式来表述:
$$ {E_x} = {R_{xx}}{H_x} + {R_{xy}}{H_y} + {R_{xz}}{H_z}, $$ (12) $$ {E_y} = {R_{yx}}{H_x} + {R_{yy}}{H_y} + {R_{yz}}{H_z}, $$ (13) $$ {E_z} = {R_{zx}}{H_x} + {R_{zy}}{H_y} + {R_{zz}}{H_z}, $$ (14) 因MT测量中总有Ez=0, 所以不必考虑(14)式, 只考虑前两式, 从而可以定义张量阻抗如下:
$$ {\boldsymbol{R}} = \left[ {\begin{array}{*{20}{c}} {{R_{xx}}}&{{R_{xy}}}&{{R_{xz}}}\\ {{R_{yx}}}&{{R_{yy}}}&{{R_{yz}}} \end{array}} \right]. $$ (15) 将[R]称为六元素张量阻抗.以字符[R]表示该张量阻抗就是为了强调它也具有电阻的单位, 并且与常规张量阻抗[Z]有所区别.
首先讨论张量阻抗[R]与张量阻抗[Z]的相互关系.将式(3)分别代入到式(12)和式(13), 然后将它们与(1)、(2)式对比, 得到以下方程:
$$ \begin{array}{l} {Z_{xx}} = {R_{xx}} + {R_{xz}}{T_{zx}},\\ {Z_{xy}} = {R_{xy}} + {R_{xz}}{T_{zy}},\\ {Z_{yx}} = {R_{yx}} + {R_{yz}}{T_{zx}},\\ {Z_{yy}} = {R_{yy}} + {R_{yz}}{T_{zy}}, \end{array} $$ (16) 上式可以再改写为下面矩阵形式:
$$ \left[ {\begin{array}{*{20}{c}} {{Z_{xx}}}&{{Z_{xy}}}\\ {{Z_{yx}}}&{{Z_{yy}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{R_{xx}}}&{{R_{xy}}}\\ {{R_{yx}}}&{{R_{yy}}} \end{array}} \right] + {R_{xz}}\;\;{R_{yz}}\left[ {\begin{array}{*{20}{c}} {{T_{zx}}}&{{T_{zy}}} \end{array}} \right]. $$ (16') 上述两组方程表示常规张量阻抗[Z]、重新定义的张量阻抗[R]和倾子[T]之间的互相关系.由此可以讨论在一维或二维的情况下的张量阻抗[R]的特点.
在一维水平均匀介质的情况下, 电场与磁场矢量总相互垂直, 即:
$$ {\boldsymbol{E}} \cdot {\boldsymbol{H}} = {E_x}{H_x} + {E_y}{H_y} + {E_z}{H_z} = 0, $$ 于是, 可以得到如下方程:
$$ \frac{{{E_x}}}{{{H_y}}} = - \frac{{{E_y}}}{{{H_x}}} = R, $$ (17) 即张量阻抗[R]变成常规张量阻抗[Z]的一个模量.
在二维TE模式下(Ey=Ez=Hz=0), 从(12)式可以得到以下方程组
$$ \begin{array}{l} {E_x} = {R_{xy}}{H_y} + {R_{xz}}{H_z},\\ {Z_{xy}} = {R_{xy}} + {R_{xz}}{T_{zy}}, \end{array} $$ (18) 同样, 在二维介质中的TM模式中, 考虑Ex=Ey=Hz=0, 由(13)式可得
$$ {E_y} = {R_{yx}}{H_x},{Z_{yx}} = {R_{yx}}. $$ (19) 从上面的推导可以获得如下结论:张量阻抗[R]在一般情况下与张量阻抗[Z]不一致, 但如果已知T与[R], 就较容易地可以取得[Z].
另外还要考虑六元素张量阻抗[R]的坐标转换问题.为此, 假设[R]为定义在XOY坐标系中的六元素张量阻抗, [R]’为将[R]向[Z]轴周围旋转θ而得到的X′OY′坐标系内的张量阻抗.那经过简单推导之后可以得到如下方程
$$ \begin{array}{l} {R^\prime }_{xx} = {R_1} - {R_2}\cos2\theta - {R_3}\sin2\theta ,\\ {R^\prime }_{xy} = {R_4} + {R_3}\cos2\theta - {R_2}\sin2\theta ,\\ {R^\prime }_{yx} = - {R_4} + {R_3}\cos2\theta - {R_2}\sin2\theta ,\\ {R^\prime }_{yy} = {R_1} + {R_2}\cos2\theta + {R_3}\sin2\theta ,\\ {R^\prime }_{xz} = {R_{xz}}\cos\theta - {R_{yz}}\sin\theta ,\\ {R^\prime }_{yz} = {R_{xz}}\sin\theta + {R_{yz}}\cos\theta , \end{array} $$ (20) 式中,
$$ \begin{array}{l} 2{R_1} = {R_{xx}} + {R_{yy}},2{R_2} = {R_{yy}} - {R_{xx}},\\ 2{R_3} = {R_{xy}} + {R_{yx}},{R_4} = {R_{xy}} - {R_{yx}}. \end{array} $$ (21) 在此可以证明, 与张量阻抗[Z]类似地, (21)式中的R1、R4和R0=RxxRyy-RxyRyx也是张量不变量.
2. 六元素张量阻抗的估计
假设任意频率MT场5个分量的n组复数振幅测量数据为(右上角参数表示测量次数, T表示转置矩阵)
$$ \begin{array}{l} {\boldsymbol{E}} = \left[ {\begin{array}{*{20}{c}} {{E_x}}&{{E_y}} \end{array}} \right],{\boldsymbol{H }}= \left[ {\begin{array}{*{20}{c}} {{H_x}}&{{H_y}}&{{H_z}} \end{array}} \right],\\ {E_\tau } = [{E_\tau }^1, \ldots ,{E_\tau ^n}]{^T},\tau = x,y,\\ {H_\tau } = [{H_\tau }^1, \ldots ,{H_\tau^n }]{^T},\tau = x,y,z. \end{array} $$ (22) 若将式(12)和(13)作为回归方程, 可以建立关于回归系数Rτx, Rτy, Rτz(τ=x, y)的最小二乘问题.
$$ \begin{array}{l} {Q_\tau } = \sum\limits_{i = 1}^n {{{\left| {{E_\tau }^i - ({R_{\tau x}}{H_x}^i + {R_{\tau y}}{H_y}^i + {R_{\tau z}}{H_z}^i)} \right|}^2}} = \\ \min,\tau = x,y, \end{array} $$ (23) 从式(23)分离实、虚部后, 取
$$ \begin{array}{l} \frac{{\partial {Q_\tau }}}{{\partial {\rm{Re}}{R_{\tau x}}}} = \frac{{\partial {Q_\tau }}}{{\partial {\rm{Im}}{R_{\tau x}}}} = \frac{{\partial {Q_\tau }}}{{\partial {\rm{Re}}{R_{\tau y}}}} = \frac{{\partial {Q_\tau }}}{{\partial {\rm{Im}}{R_{\tau y}}}} = \\ \frac{{\partial {Q_\tau }}}{{\partial {\rm{Re}}{R_{\tau z}}}} = \frac{{\partial {Q_\tau }}}{{\partial {\rm{Im}}{R_{\tau z}}}} = 0,\tau = x,y, \end{array} $$ 即可获得对于6个未知数ReRτx, ImRτx, ReRτy, ImRτy, ReRτz, ImRτz的线性方程.
在此, 任何两个复变量的平均功率谱表示为
$$ \overline {AB} = \frac{1}{n}\sum\limits_{i = 1}^n {\overline {{A^i}} {B^i}} . $$ (25) 结合上述六个表达式后, 可得下面的矩阵方程(上横线表示复数共轭)
$$ \begin{array}{l} {\left[ {\begin{array}{*{20}{c}} {\overline {{H_x}{H_x}} }&{\overline {{H_x}{H_y}} }&{\overline {{H_x}{H_z}} }\\ {\overline {{H_y}{H_x}} }&{\overline {{H_y}{H_y}} }&{\overline {{H_y}{H_z}} }\\ {\overline {{H_z}{H_x}} }&{\overline {{H_z}{H_y}} }&{\overline {{H_z}{H_z}} } \end{array}} \right]^{ - 1}} \cdot \\ \left[ {\begin{array}{*{20}{c}} {{R_{xx}}}&{{R_{yx}}}\\ {{R_{xy}}}&{{R_{yx}}}\\ {{R_{xz}}}&{{R_{yx}}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {\overline {{H_x}{E_x}} }&{\overline {{H_x}{E_y}} }\\ {\overline {{H_y}{E_x}} }&{\overline {{H_y}{E_y}} }\\ {\overline {{H_z}{E_x}} }&{\overline {{H_z}{E_y}} } \end{array}} \right], \end{array} $$ (26) 这样六元素张量阻抗的估计公式为
$$ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {{R_{xx}}}&{{R_{yx}}}\\ {{R_{xy}}}&{{R_{yx}}}\\ {{R_{xz}}}&{{R_{yx}}} \end{array}} \right] = \\ {\left[ {\begin{array}{*{20}{c}} {\overline {{H_x}{H_x}} }&{\overline {{H_x}{H_y}} }&{\overline {{H_x}{H_z}} }\\ {\overline {{H_y}{H_x}} }&{\overline {{H_y}{H_y}} }&{\overline {{H_y}{H_z}} }\\ {\overline {{H_z}{H_x}} }&{\overline {{H_z}{H_y}} }&{\overline {{H_z}{H_z}} } \end{array}} \right]^{ - 1}}\left[ {\begin{array}{*{20}{c}} {\overline {{H_x}{E_x}} }&{\overline {{H_x}{E_y}} }\\ {\overline {{H_y}{E_x}} }&{\overline {{H_y}{E_y}} }\\ {\overline {{H_z}{E_x}} }&{\overline {{H_z}{E_y}} } \end{array}} \right], \end{array} $$ (27) 展开上式, 可得每个Rij=(i=x, y; j=x, y, z)的表达式.
由此可以知道为了估计[Z]需要最少两次测量资料, 但为了估计[R]至少需要3次测量.
方程(27)一般可以写为
$$ \begin{array}{l} \left[ {\begin{array}{*{20}{c}} {{R_{xx}}}&{{R_{yx}}}\\ {{R_{xy}}}&{{R_{yx}}}\\ {{R_{xz}}}&{{R_{yx}}} \end{array}} \right] = \\ {\left[ {\begin{array}{*{20}{c}} {\overline {A{H_x}} }&{\overline {A{H_y}} }&{\overline {A{H_z}} }\\ {\overline {B{H_x}} }&{\overline {B{H_y}} }&{\overline {B{H_z}} }\\ {\overline {C{H_x}} }&{\overline {C{H_y}} }&{\overline {C{H_z}} } \end{array}} \right]^{ - 1}}\left[ {\begin{array}{*{20}{c}} {\overline {A{E_x}} }&{\overline {A{E_y}} }\\ {\overline {B{E_x}} }&{\overline {B{E_y}} }\\ {\overline {C{E_x}} }&{\overline {C{E_y}} } \end{array}} \right]. \end{array} $$ (28) 可以将Ex, Ey, Hx, Hy, Hz中的任意三个量选择为上式中的(A, B, C), 于是,张量阻抗[R]一共有$C_5^2 = \frac{{5!}}{{3! \cdot 2!}} = 10$种估计方法, 但常规张量阻抗[Z]的估计方法只有6种.
考虑(22)和(25)式, (26)和(27)式可以简化为(*表示Hermitian矩阵):
$$ \left[ {H*H} \right]\left[ R \right] = \left[ {H*E} \right], $$ (26') $$ \left[ R \right] = {\left[ {H*H} \right]^{ - 1}}\left[ {H*E} \right], $$ (27') 远点参考MT测量中计算六元素张量阻抗仍然可行.如果将在测点1, 远点2同时测量的电场, 磁场分别用E1,H1, E2, H2表示, 两个测点的六元素张量阻分别用R1, R2表示, R1, R2的估计计算式为如下:
$$ \begin{array}{l} \left[ {{R_1}} \right] = {\left( {{H_2}*{H_1}} \right)^{ - 1}}\left[ {{H_2}*{E_1}} \right],\\ \left[ {{R_2}} \right] = {\left[ {{H_1}*{H_2}} \right]^{ - 1}}\left[ {{H_1}*{E_2}} \right]. \end{array} $$ (29) 3. 六元素张量阻抗的应用实例
本文中引用的MT数据均由朝鲜开发的五通道KTU-MT-1系统(384-0.001Hz频段)采集.
首先对比分析MT资料中磁场垂直分量和水平分量的平均振幅, 此时, 时间域中各个分量的平均绝对振幅可以按照如下公式计算
$$ \overline{H_i} = \frac{1}{{{N_s}}}\sum\limits_{s = 1}^{Ns} {\left( {\frac{1}{N}\sum\limits_{n = 1}^{Nt} {} \left| {{H_i}\left( {s,n} \right)} \right|} \right)} ,i = x,y,z $$ 式中Ns为MT测量点的总数, Nt为在测点上观测时间列的总采样数, Hi(s, n)为在第S个MT测点第i次磁场分量的第n个采样.
表 1为磁场三个分量的振幅分析结果, 此试验中使用的大地电磁测点总数为70个以上.
表 1 MT时间序列中磁场分量的振幅比较Table Supplementary Table Comparison of amplitudes of magnetic components from MT records频率段 采样率t(s) 磁场分量平均振幅(mA/m) Hx Hy Hz 1 6.510 4×10-4 0.562 1.893 0.560 2 1.016 7×10-2 0.078 0.041 0.063 3 0.166 667 0.183 0.149 0.350 4 2.666 67 6.279 10.748 17.410 平均 1.775 3.208 4.596 从表 1可知, 在4个频率范围内垂直磁场分量的振幅几乎等于水平磁场分量, 因此不能忽视垂直磁场分量对水平电场分量的影响.图 1为在一个人文噪声严重的MT测点随频率变化的电场与磁场分量之间的相干度曲线, 相干度计算公式为:
$$ \begin{array}{l} coh({E_i}{H_j}) = \frac{{\overline {\left| {{E_i}{H_j}} \right|} }}{{{{[{E_i}{E_i}\cdot{H_j}{H_j}]}^{1/2}}}},i = x,y,j = \\ x,y,z. \end{array} $$ (30) 表 2为上述测点中算出的Ei(i=x, y)分量和Hj(j=x, y, z)分量之间的平均相干度.
表 2 图 1中的测点Ei(i=x, y)和Hj(j=x, y, z)之间的平均相干度Table Supplementary Table Coherencies between Ei(i=x, y) and Hj(j=x, y, z) componentsCoh(Ex, Hx) Coh(Ex, Hy) Coh(Ex, Hz) Coh(Ey, Hx) Coh(Ey, Hy) Coh(Ey, Hz) 0.488 0.456 0.409 0.473 0.419 0.405 从图 1和表 2可以看出, 在大多数频点中Ei和Hz之间相干度coh(EiHz)(i=x, y)不小于Ei和Hx(或者Hy)之间相干度coh(EiHx)(或者coh(EiHy))(i=x, y), 这意味着不能忽视Hz分量对Ex和Ey分量的影响,使用传统的四元张量阻抗将带来显著的误差.
表 3 图 1中的MT测点中用[Z]和[R]时Ei, Eip之间平均相干度Table Supplementary Table Comparison of coherencies between Ei, Eip components when using [Z] and [R]相干度 Coh(Ex) Coh(Ey) [Z] 0.659 0.641 [R] 0.692 0.692 通常MT法中为了评估观测数据的质量, 可以利用实测电场与预测电场之间的相干度,即:
$$ coh({E_i}) = \frac{{\overline {\left| {{E_i}{E_i}^p} \right|} }}{{{{[{E_i}{E_i}\cdot{E_i}^p{E_i}^p]}^{1/2}}}},i = x,y. $$ 图 2为相同的MT测点中用上式算出来的coh[Ei(R)]和coh[Ei(Z)]曲线的对比例子, 此时利用张量阻抗[Z]和[R]预测到的电场分别为:
$$ \begin{array}{l} {E_i}^p\left( Z \right) = {Z_{ix}}{H_x} + {Z_{iy}}{H_y},\\ {E_i}^p\left( R \right) = {R_{ix}}{H_x} + {R_{iy}}{H_y} + {R_{iz}}{H_z}. \end{array} $$ 表 3为用张量阻抗[Z]和[R]得到的Ei和Eip之间的平均相干度.
通过实际数据处理可以明显地看到对于大多频率点coh[Eip(R)]≥coh[Eip(Z)], 即利用六元素张量阻抗[R]比利用常规四元素张量阻抗[Z]可以提高实测电场和预测电场之间的相干度,从而改善大地电磁资料的信噪比.
图 3为在和图 1、2相同的MT测点上利用常规阻抗张量[Z]和新的阻抗张量[R]得到的视电阻率ρxy、ρyx曲线.实验表明, 虽然两条曲线基本上符合, 但由[Z]得到的曲线比由[R]得到的曲线更分散,常规方法的抗噪性能不如六元素张量阻抗方法.
图 4表示由图 3的曲线计算出的等效视电阻率曲线和由它的反演获得的深度-电阻率曲线.此时, 图 4a的1, 2分别表示用张量阻抗[Z]和[R]计算的等效视电阻率曲线, 它们反演出来的结果分别为图 4b中的深度-电阻率曲线(层状电阻率模型)3和4;图 4a中的实线表示对应于图 4b中电阻率模型3, 4的正演拟合视电阻率曲线.
为了实现MT视电阻率曲线的反演,利用与直流电阻率法中的Zohdy反演(1989)类似的方法, 即先利用Bostick公式获得深度及电阻率的初始值, 按照两条视电阻率曲线(实测与正演计算曲线)的符合程度, 改正深度和电阻率, 该方法速度较快, 收敛稳定.由此可见, 虽然在高频、中频段里两条曲线的一致性较差, 但基本上得到了较满意的反演结果.
4. 结论
在MT资料处理中新定义的六元素张量阻抗[R]比传统的四分量张量阻抗[Z]提高测深曲线的相干度至少为0.1以上, 并且改善了大地电磁资料处理的质量.但是, 本文所提出的方法需要较多的观测数据,其数据量至少为传统方法的1.5倍, 这是该方法的一个不足之处.
目前所有的MT法均以四元素张量阻抗的定义作为理论基础, 本文提出的六元张量阻抗方法还在试验阶段, 需要进一步的研究.为了利用已有成熟的四元素张量阻抗处理和反演方法,我们可以将六元素张量阻抗[R]作为计算常规张量阻抗[Z]的中间辅助变量,同样可以提高数据质量.
-
表 1 MT时间序列中磁场分量的振幅比较
Table 1. Comparison of amplitudes of magnetic components from MT records
频率段 采样率t(s) 磁场分量平均振幅(mA/m) Hx Hy Hz 1 6.510 4×10-4 0.562 1.893 0.560 2 1.016 7×10-2 0.078 0.041 0.063 3 0.166 667 0.183 0.149 0.350 4 2.666 67 6.279 10.748 17.410 平均 1.775 3.208 4.596 表 2 图 1中的测点Ei(i=x, y)和Hj(j=x, y, z)之间的平均相干度
Table 2. Coherencies between Ei(i=x, y) and Hj(j=x, y, z) components
Coh(Ex, Hx) Coh(Ex, Hy) Coh(Ex, Hz) Coh(Ey, Hx) Coh(Ey, Hy) Coh(Ey, Hz) 0.488 0.456 0.409 0.473 0.419 0.405 表 3 图 1中的MT测点中用[Z]和[R]时Ei, Eip之间平均相干度
Table 3. Comparison of coherencies between Ei, Eip components when using [Z] and [R]
相干度 Coh(Ex) Coh(Ey) [Z] 0.659 0.641 [R] 0.692 0.692 -
[1] Auken, E., Boesen, T., Christiansen, A.V., 2017.A Review of Airborne Electromagnetic Methods with Focus on Geotechnical and Hydrological Applications from 2007 to 2017.Advances in Geophysics. [2] Berdichevsky, M.N., 1968.Electrical Prospecting with Magnetotelluric Profiling.Nedra, Moscow(in Russia). [3] Cao, Z.L., He, Z.X., Chang, Y.J., 2006.A Simulation Study of Induced Polarization Effect of Magnetotelluric and Its Application in Oil and Gas Detection.Progress in Geophysics, 21(4):1252-1257(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DQWJ200604031.htm [4] de Groot-Hedlin, C., Constable, S., 2004.Inversion of Magnetotelluric Data for 2D Structure with Sharp Resistivity Contrasts.Geophysics, 69(1):78-86. https://doi.org/10.1190/1.1649377 [5] Egbert, G.D., Livelybrooks, D.W., 1996.Single Station Magnetotelluric Impedance Estimation:Coherence Weighting and the Regression M-Estimate.Geophysics, 61(4):964-970. https://doi.org/10.1190/1.1444045 [6] Gamble, T.D., Goubau, W.M., Clarke, J., 1979.Magnetotellurics with a Remote Magnetic Reference.Geophysics, 44(1):53-68. https://doi.org/10.1190/1.1440923 [7] García, X., Jones, A.G., 2002.Atmospheric Sources for Audio-Magnetotelluric (AMT) Sounding.Geophysics, 67(2):448-458. https://doi.org/10.1190/1.1468604 [8] García, X., Jones, A.G., 2005.A New Methodology for the Acquisition and Processing of Audio-Magnetotelluric (AMT) Data in the AMT Dead Band.Geophysics, 70(5):G119-G126. https://doi.org/10.1190/1.2073889 [9] Gómez-Treviño, E., Mondragón, M., 1995.Uneven Effect of Random Noise in Magnetotelluric Apparent Resistivity Definitions.Geophysics, 60(4):1238-1242. https://doi.org/10.1190/1.1443853 [10] Groom, R.W., Bailey, R.C., 1991.Analytic Investigations of the Effects of Near-Surface Three-Dimensional Galvanic Scatterers on MT Tensor Decompositions.Geophysics, 56(4):496-518. https://doi.org/10.1190/1.1443066 [11] Gupta, D., Choubey, S., 2015.Discrete Wavelet Transform for Image Processing.International Journal of Emerging Technology and Advanced Engineering, 4(3):598-602. http://d.old.wanfangdata.com.cn/OAPaper/oai_doaj-articles_36d77fb15fc197b8636242fbc09e7525 [12] Hu, Z.Z., Hu, X.Y., 2005.Review of Three Dimensional Magnetotelluric Inversion Methods.Progress in Geophysics, 20(1):214-220(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DQWJ2005010BA.htm [13] Li, D.W., Wang, Y.X., 2015.Major Issues of Research and Development of Hot Rock Geothermal Energy.Earth Science, 40(11):1858-1869 (in Chinese with English abstract). https://doi.org/10.3799/dqkx.2015.166 [14] Lilley, F.E.M., 1998.Magnetotelluric Tensor Decomposition:Part I, Theory for a Basic Procedure.Geophysics, 63(6):1885-1897. https://doi.org/10.1190/1.1444481 [15] Linde, N., Pedersen, L.B., 2004.Characterization of a Fractured Granite Using Radio Magnetotelluric (RMT) Data.Geophysics, 69(5):1155-1165. https://doi.org/10.1190/1.1801933 [16] Liu, W.C., Zhang, S.Y., Yang, L.B., et al., 2015.Three-Dimensional Electrical and Deep Structure Features of Akebasitao Area in Western Junggar by AMT Data.Earth Science, 40(3):441-447 (in Chinese with English abstract). https://doi.org/10.3799/dqkx.2015.036 [17] Muñoz, G., 2014.Exploring for Geothermal Resources with Electromagnetic Methods.Surveys in Geophysics, 35(1):101-122. https://doi.org/10.1007/s10712-013-9236-0 [18] Pedersen, L.B., Engels, M., 2005.Routine 2D Inversion of Magnetotelluric Data Using the Determinant of the Impedance Tensor.Geophysics, 70(2):G33-G41. https://doi.org/10.1190/1.1897032 [19] Ri, G., Kim, K., Ri, Y., 2002.The Magnetotelluric Exploration.Publisher of Industry, PyongYang. [20] Tripp, A.C., 2005.Acheron's Rainbow:Free Associations on 75 Years of Exploration Geo-Electromagnetics.Geophysics, 70(6):25-31. https://doi.org/10.1190/1.2127107 [21] Wang, J.Y., Xu, Y.X., 1998.Methods and Advances for Estimation of Magnetotelluric Response Function Abroad.Earth Science Frontiers, 5(2):217-222(in Chinese with English abstract). http://en.cnki.com.cn/article_en/cjfdtotal-dxqy802.006.htm [22] Wang, S.M., Wang, J.Y., 2004.Application of Higher-Order Statistics in Magnetotelluric Data Processing.Chinese Journal of Geophysics, 47(5):928-934(in Chinese with English abstract). doi: 10.1002/cjg2.584/full [23] Xu, Y.X., Wang, J.Y., 2000.A Spectrum Estimation Method for Magnetotelluric Signal Based on Continuous Wavelet Transforms.Chinese Journal of Geophysics, 43(5):717-723. https://doi.org/10.1002/cjg2.86 [24] Xu, Z., 1987.The Processing of Five-Component Magnetotelluric Sounding Data.Oil Geophysical Prospecting, 22(4):435-444. http://en.cnki.com.cn/Article_en/CJFDTotal-SYDQ198704010.htm [25] Yee, E., Kosteniuk, P.R., Paulson, K.V., 1988.The Reconstruction of the Magnetotelluric Impedance Tensor:An Adaptive Parametric Time-Domain Approach.Geophysics, 53(8):1080-1087. https://doi.org/10.1190/1.1442544 [26] Yin, C.C., 2003.Inherent Nonuniqueness in Magnetotelluric Inversion for 1D Anisotropic Models.Geophysics, 68(1):138-146. https://doi.org/10.1190/1.1543201 [27] Zhang, Y., Zhang, S.Y., 2005.A study of the EH-4 Processing and Interpertation System.Chinese Journal of Engineering Geophysics, 2(4):311-315(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-GCDQ200504010.htm [28] Zohdy, A.A.R., 1989.A New Method for the Automatic Interpretation of Schlumberger and Wenner Sounding Curves.Geophysics, 54(2):245-253. https://doi.org/10.1190/1.1442648 [29] 曹中林, 何展翔, 昌彦君, 2006.MT激电效应的模拟研究及在油气检测中的应用.地球物理学进展, 21(4):1252-1257. doi: 10.3969/j.issn.1004-2903.2006.04.032 [30] 胡祖志, 胡祥云, 2005.大地电磁三维反演方法综述.地球物理学进展, 20(1):214-220. doi: 10.3969/j.issn.1004-2903.2005.01.037 [31] 李德威, 王焰新, 2015.干热岩地热能研究与开发的若干重大问题.地球科学, 40(11):1858-1869. https://doi.org/10.3799/dqkx.2015.166 [32] 刘文才, 张胜业, 杨龙彬, 等, 2015.西准噶尔阿克巴斯陶地区三维电性结构和深部地质特征.地球科学, 40(3):441-447. https://doi.org/10.3799/dqkx.2015.036 [33] 王家映, 徐义贤, 1998.国外大地电磁响应函数估计方法.地学前缘, 5(2):217-222. doi: 10.3321/j.issn:1005-2321.1998.02.006 [34] 王书明, 王家映, 2004.高阶统计量在大地电磁测深数据处理中的应用研究.地球物理学报, 47(5):928-934. doi: 10.3321/j.issn:0001-5733.2004.05.027 [35] 张莹, 张胜业, 2005.EH-4资料处理解释系统的研究.工程地球物理学报, 2(4):311-315. doi: 10.3969/j.issn.1672-7940.2005.04.011 期刊类型引用(2)
1. 王维强,宋洪伟,尚铭森,母海东,闫丰录. 铁路隧道水热害不良地质体大地电磁场特征研究——以川藏铁路茶洛隧道为例. 钻探工程. 2021(05): 9-17 . 百度学术
2. 陈斌,汪耀,胡祥云,李倩婕. 大湾区珠江口海上高密度电法探测. 地球科学. 2020(12): 4550-4562 . 本站查看
其他类型引用(3)
-