大地电磁全张量响应的一维各向异性反演

目前大地电磁(MT)测深资料反演主要基于各向同性介质,但随着MT实际应用的需要,各向异性研究已逐渐引起关注.我们采用广泛应用的广义逆法对一维MT水平层状各向异性介质模型反演进行了探索性研究,并实现了MT全张量响应(即所有的阻抗张量的视电阻率和相位)的一维各向异性反演.理论模型试验表明,无论理论观测值中是否含有噪声,这种方法都能够较好地恢复真实模型,验证了其正确性和有效性.将此方法用于MT实测资料时,能够同时拟合4对视电阻率和阻抗相位曲线,说明本方法可以用于实测资料的处理解释,具有一定的实用价值.     

Magnetotelluric inversion for azimuthally anisotropic resistivities employing artificial neural networks

An extension of an artificial neural network (ANN) approach to solve the magnetotelluric (MT) inverse problem for azimuthally anisotropic resistivities is presented and applied for a real dataset. Three different model classes, containing general 1-D and 2-D azimuthally anisotropic features, have been considered. For each model class, characteristics of three-layer feed forward ANNs trained through an error back propagation algorithm have been adjusted to approximate the inverse modeling function. It appears that, at least for synthetic models, reasonable results would be obtained by applying the amplitudes of the complex impedance tensor elements as inputs. Furthermore, the Levenberg-Marquart algorithm possesses optimal performance as a learning paradigm for this problem. The evaluation of applicability of the trained ANNs for unknown data sets excluded from the learning procedure reveals that the trained ANNs possess acceptable interpolation and extrapolation abilities to estimate model parameters accurately. This method was also successfully used for a field dataset wherein anisotropy had been previously recognized.     

大地电磁资料精细处理和二维反演解释技术研究(一)——阻抗张量分解与构造维性分析

本文对大地电磁观测阻抗实施一种数学变换--"共轭阻抗变换",发现转换后的观测阻抗和区域阻抗之间存在特定的关系,这种特定关系不受电场局部畸变的影响,而且同样不需要关于地下区域结构维性的假设,在区域结构是三维的情况下也是成立的.对转换后的观测阻抗采用已有的Swift旋转方法即可求得区域主轴方位角,然后采用最优化方法求取区域阻抗相位、振幅以及畸变因子.同时根据转换后的观测阻抗重新定义了不受电场局部畸变影响的构造维性参数.采用合成理论数据验证了新算法的正确性,和Swift、Bahr、GB、相位张量、WAL方法进行了对比分析,并将新方法应用于实测资料的解释,发展了一套消除局部畸变和进行构造维性分析的MT精细资料处理技术.进一步的工作是将新算法推广到多点多频,并采用最优化技术分解,发展一种稳定性好、多测点多频点的阻抗张量分解技术.     

大地电磁测深数据处理——4元素和6元素阻抗张量元素的比较研究(英文)

阻抗张量元素的计算是在大地电磁测深数据处理的重要一步。按照常规,阻抗张量被定义为以Zxx,Zxy,Zyx,和Zyy为元素的2×2矩阵。在本次研究中,6个元素的阻抗张量的计算使用了一个含有Zxx,Zxy,Zyx,Zyy,Zxz和Zyz分量的2x3矩阵。对上述两类阻抗张量元素的属性进行了分析。利用由印度古吉拉特邦卡奇沉积盆地采集的5个分量大地电磁数据测试了文中的方法。从视电阻率和相位的计算中我们是观察到在大部分的频带范围内4个元素阻抗和6个元素阻抗Zxy和Zyx两类元素区别不大。然而,较长周期时间的数据,如超过100秒,观察到视电阻率的增加和相位的减少。我们还注意到,倾子幅度在大部分时间几乎是零,但较长周期(超过100秒),逐渐呈增加的趋势。卡奇沉积盆地的地电断面表明在较长的周期内浅层近水平层和深层异常高电导性的不均质层都可能是引起大的Hz分量的原因。这表明,磁场垂直分量Hz对在大的2D/3D结构区域内的电场参数估计发挥的重要作用。     

Influence of frequency-dependent anisotropy on seismic amplitude-versus-offset signatures for fractured poroelastic rocks

Frequency-dependent amplitude variation with offset offers an effective method for hydrocarbon detections and analysis of fluid flow during production of oil and natural gas within a fractured reservoir. An appropriate representation for the frequency dependency of seismic amplitude variation with offset signatures should incorporate influences of dispersive and attenuating properties of a reservoir and the layered structure for either isotropic or anisotropic dispersion analysis. In this study, we use an equivalent medium permeated with aligned fractures that simulates frequency-dependent anisotropy, which is sensitive to the filled fluid of fractures. The model, where pores and fractures are filled with two different fluids, considers velocity dispersion and attenuation due to mesoscopic wave-induced fluid flow. We have introduced an improved scheme seamlessly linking rock physics modelling and calculations for frequency-dependent reflection coefficients based on the propagator matrix technique. The modelling scheme is performed in the frequency-slowness domain and can properly incorporate effects of both bedded structure of the reservoir and velocity dispersion quantified with frequency-dependent stiffness. Therefore, for a dispersive and attenuated layered model, seismic signatures represent a combined contribution of impedance contrast, layer thickness, anisotropic dispersion of the fractured media and tuning and interference of thin layers, which has been avoided by current conventional methods. Frequency-dependent amplitude variation with offset responses was studied via considering the influences of fracture fills, layer thicknesses and fracture weaknesses for three classes amplitude variation with offset reservoirs. Modelling results show the applicability of the introduced procedure for interpretations of frequency-dependent seismic anomalies associated with both layered structure and velocity dispersion of an equivalent anisotropic medium. The implications indicate that anisotropic velocity dispersion should be incorporated accurately to obtain enhanced amplitude variation with offset interpretations. The presented frequency-dependent amplitude variation with offset modelling procedure offers a useful tool for fracture fluid detections in an anisotropic dispersive reservoir with layered structures.     

Enhancement of Signal-to-noise Ratio in Natural-source Transient Magnetotelluric Data with Wavelet Transform

—For audio-frequency magnetotelluric surveys where the signals are lightning-stroke transients, the conventional Fourier transform method often fails to produce a high quality impedance tensor. An alternative approach is to use the wavelet transform method which is capable of localizing target information simultaneously in both the temporal and frequency domains. Unlike Fourier analysis that yields an average amplitude and phase, the wavelet transform produces an instantaneous estimate of the amplitude and phase of a signal. In this paper a complex well-localized wavelet, the Morlet wavelet, has been used to transform and analyze audio-frequency magnetotelluric data. With the Morlet wavelet, the magnetotelluric impedance tensor can be computed directly in the wavelet transform domain. The lightning-stroke transients are easily identified on the dilation-translation plane. Choosing those wavelet transform values where the signals are located, a higher signal-to-noise ratio estimation of the impedance tensor can be obtained. ? In a test using real data, the wavelet transform showed a significant improvement in the signal-to-noise ratio over the conventional Fourier transform.     

Improvement of MT data processing using stationary and coherence tests

The great majority of data processing methods for magnetotelluric measurements are based on an average Fourier spectral analysis to estimate either the transfer function or the coherence function. This assumes that the time-series data are stationary. The aim of this paper is to present a preprocessing method able to extract the stationary and coherent part of the original signals. The practitioner can then apply the magnetotelluric method of his choice to these new data. This preliminary data sorting is done in four steps: (i) slow drifts are eliminated using a high-pass filter whose cut-off frequency is determined by an iterative procedure. Based on run tests, this procedure is also able to remove segments with non-independent samples in the time series; (ii) non-stationary segments are eliminated after band-pass filtering; (iii) non-coherent segments are eliminated before spectral analysis; (iv) the impedance tensor value is then retained, at a given frequency, only if the signals are coherent. This preprocessing method was tested on the simplest, but still used, magnetotelluric method which uses only two field components, and it was found that the average resistivity standard deviation decreased significantly from 14.6 Ωm without sorting to 8.6 Ωm after sorting.     

One-dimensional magnetotelluric inversion using an adaptation of Zohdy's resistivity method

An iterative refinement method for determining a layered resistivity model from a Schlumberger or Wenner sounding curve is adapted to determine a layered resistivity model by using apparent resistivity and phase derived from the magnetotelluric impedance. Magnetotelluric observations presented as a function of period are first converted to an approximate resistivity–depth profile using Schmucker's transformation and this is used to construct an initial guess (starting) model. A two-stage procedure is then invoked. Keeping resistivities constant, layer boundaries are first adjusted to give a minimum misfit between measured data and responses and this is followed by resistivity adjustments with fixed layer boundaries to reduce the misfit further. The method is illustrated by application to some synthetic data both exact and with added noise, to a real field data set and to some magnetotelluric profile data obtained in a survey over the Carnmenellis granites in south Cornwall. The method is validated by recovering conductivity models from the exact and noisy 1D synthetic data. For complicated three-dimensional data at a single site and along a profile of stations, the method is shown to produce acceptable solutions which may be used as starting models in further two- or three-dimensional studies.     

Three-dimensional electrical structure of the northwestern margin of the Karamay region,China, revealed by magnetotelluric data

In order to investigate the three-dimensional structures of intrusive granite and the deep structure of the Darbut fault in the northwestern margin of the Karamay region, western Junggar Basin, China, new magnetotelluric data were collected along six profiles across the Darbut fault. The magnetotelluric data were processed using a robust estimation technique to obtain the magnetotelluric impedance. Then the off-diagonal components of the impedance tensor were inverted using a three-dimensional nonlinear conjugate gradient inversion technique, which was performed using open-source three-dimensional electromagnetic inversion software. The final three-dimensional model includes two major low-resistivity anomalies and two major high-resistivity anomalies. The first low-resistivity anomaly corresponds to the location of the Darbut fault, which indicates that metallic elements are abnormally enriched there. The second one is located beneath the Darbut fault, and it is most likely a magma channel in the middle crust. The two major high-resistivity anomalies are distributed on either side of the Darbut fault. We interpret them as Karamay rocks and Akebasitao rocks and suggest that they were formed in an extensional setting. The cross sections of three-dimensional magnetotelluric inversion result reveal that the Darbut fault has been reformed by the later magmatism, leading to the variation of its downward extent along its strike. Moreover, our inversion result also indicates that a magma channel exists in the middle crust of the region.     

大地电磁资料的灵敏度研究

本文研究了大地电磁资料对地球电性结构模型参数的灵敏度,它在反演解释中具有重要意义.通过研究均匀半空间中视电导率资料的灵敏度函数与计算分析实际观测资料的Fréchet导数和斜率表明,视电导率资料的灵敏度在地表最大,从地表向下呈指数衰减,频率愈高,电导率愈大,衰减也愈快.视电阻率资料和阻抗相位资料对高频、浅层和低阻层具有更高的灵敏度,几何参数的平均灵敏度高于电性参数的平均灵敏度,通常阻抗相位资料比视电阻率资料更灵敏.同时,根据阻抗实部定义的视电阻率与Kunetz定义的时间域大地电磁响应,反映地层电性差异的灵敏度最高. 文章还指出,在实际的反演解释中,如何正确有效地应用阻抗相位资料,提高分辨能力,仍是急待解决的问题.     

TOWARDS UNIFICATION OF GEOPHYSICAL PROBLEMS FOR HORIZONTALLY STRATIFIED MEDIA*

Theorems and relations describing a system of horizontal layers are considered as the elements of a meta-system having a layered structure. This implies a layered structure of a solution of geophysical problems for horizontally stratified medium. Wave equation, Laplace's equation, and Maxwell's equations are the special cases of a more general differential equation. By separation of variables the basic equation for layered structures is obtained. A multi-layer medium may result from accumulation of successive layers. The descending or ascending accumulation is possible, depending whether the successive layers are inserted at the top or at the bottom of the medium. The discretization of a basic equation provides two sets of formulas for descending or ascending accumulation, respectively. The relationship between the mathematical expressions of the seismic transfer function and the magnetotelluric input impedance is derived. The deduction leads to the possibility of converting the magnetotelluric observations into an equivalent synthetic seismogram.     

Processing of magnetotelluric data

The processing of magnetotelluric data involves concepts from electromagnetic theory, time series analysis and linear systems theory for reducing natural electric and magnetic field variations recorded at the earth's surface to forms suitable for studying the electrical properties of the earth's interior.The electromagnetic field relations lead to either a scalar transfer impedance which couples an electric component to an orthogonal magnetic component at the surface of a plane-layered earth, or a tensor transfer impedance which couples each electric component to both magnetic components in the vicinity of a lateral inhomogeneity.A number of time series spectral analysis methods can be used for estimating the complex spectral coefficients of the various field quantities. These in turn are used for estimating the nature of the transfer function or tensor impedance. For two dimensional situations, the tensor impedance can be rotated to determine the principal directions of the electrical structure.In general for real data, estimates of the apparent resistivity are more stable when calculated from the tensor elements rather than from simple orthogonal field ratios (Cagniard estimates), even when the fields are measured in the principal coordinates.     

3-D inversion of MT data from the Sabalan geothermal field,Ardabil, Iran

Since the true Earth is 3-D in nature, a three-dimensional (3-D) inversion has clear advantages over lower dimensional inversions. We utilized a 3-D magnetotelluric (MT) inversion code, WSINV3DMT, to obtain a realistic resistivity model using a long period MT data set collected in the Northwest Sabalan geothermal field in Ardabil, Iran. The apparent resistivity and phase curves, the magnetic induction vectors, the impedance polar diagrams and the rotational invariant of impedance tensor, indicate a complex 3-D conductivity structure. After setting up the model parameters and designing the appropriate block discretization, we performed the 3-D inversions for two sets of observed data; one set includes the full MT impedance tensor and another set contains only off-diagonal elements of the MT impedance. The final model was selected according to the relative magnitude of the data misfit and the model norm with respect to various Lagrangian multipliers. The results of this study illustrate the 3-D inversion of the off-diagonal elements of MT impedance tensor is precisely enough to explain the structures related to the geothermal source. The obtained results were compared with the results of available 2-D models and they are then interpreted using all of the geological and drilling data of the area. The main outcome of this study is the precise delineation of the geometry of geothermal source that is located at the center of the study area with a surface coverage of about 7 km².     

大地电磁资料精细处理和二维反演解释技术研究(五)——利用阻抗张量成像识别大地线性构造

大地线性构造往往与块体或物性的边界、断裂的几何延伸相关联,是地学研究的重要目标构造.在大地电磁法中,传统的方法是通过反演得到的电性结构的差异来解译线性构造,是一种间接分析方法.由于线性构造发育处走向和倾向明显,二维性较其他地方一般会更强,故而在阻抗张量成像分析中,相应的二维有效因子就会比较大,最佳主轴统计分布会更加明显.据此,本文利用大地电磁阻抗张量成像(包括最佳主轴统计成像、构造维性成像)来探测识别大地线性构造.本文首先通过三维理论模型的研究结果来验证这一技术在线性结构探测中的有效性,然后将其应用于郯庐断裂带南段实测剖面的具体分析中,展示了这一技术在线性构造探测中的实用性及其数据处理与分析过程.本文研究表明,大地电磁阻抗张量成像技术可独立于反演为地质解释尤其是线性构造的解释提供诸多参考信息.     

利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性

大地电磁野外实测数据月前大多为二维剖面数据.如何反演这些二维剖面数据获得较为接近实际地电情况的结果,是多数大地电磁工作者关心的问题.我们通过对理论模型的三维响应进行分析和对合成数据及实测资料的反演结果进行对比研究,讨论了利用三维反演的方法来获得大地电磁二维剖面附近三维电阻率结构的可行性.结果表明:可用三维反演的方法来解...     

Stability of the time‐domain analysis method including a frequency‐dependent soil–foundation system

A number of methods have been proposed that utilize the time‐domain transformations of frequency‐dependent dynamic impedance functions to perform a time‐history analysis. Though these methods have been available in literature for a number of years, the methods exhibit stability issues depending on how the model parameters are calibrated. In this study, a novel method is proposed with which the stability of a numerical integration scheme combined with time‐domain representation of a frequency‐dependent dynamic impedance function can be evaluated. The method is verified with three independent recursive parameter models. The proposed method is expected to be a useful tool in evaluating the potential stability issue of a time‐domain analysis before running a full‐fledged nonlinear time‐domain analysis of a soil–structure system in which the dynamic impedance of a soil–foundation system is represented with a recursive parameter model. Copyright © 2015 John Wiley & Sons, Ltd.     

SOME EFFECTS OF MULTIPLE LATERAL IN HOMOGENEITIES IN MAGNETOTELLURICS1

The analytical solution, for the H-polarization magnetotelluric impedance, of a series of multiple, vertical, conducting slabs (dikes) embedded in a host medium is extended to an infinite array in order to model an anisotropic layer. The solution is used to study the effects of such strongly anisotropic media on the surface impedance. At low frequencies such vertically laminated structures behave as a bulk material. It is shown that the effective bulk parameters are those expected from d.c. theory. However, practical electrode separations may not be long enough or adequately positioned to correctly obtain these bulk parameters from the measured impedance. Thus, such structures can masquerade as quite different one-dimensional structures. A haphazard use of long electrode spacings will not necessarily produce correct results.     

大地电磁全信息资料三维共轭梯度反演研究(英文)

在对张量阻抗数据、倾子数据和共轭梯度算法深入分析的基础上,我们实现了大地电磁全信息资料三维共轭梯度反演算法。基于全信息资料的三维共轭梯度反演研究,探讨了同时利用五个电磁场分量整理得到的大地电磁资料进行三维反演定量解释的方法以及全信息数据在三维反演中的作用。理论模型合成数据的反演结果表明,在三维反演中使用张量阻抗和倾子数据结合的全信息数据的反演结果优于只使用张量阻抗数据(或只使用倾子数据)的反演结果,提高了反演结果的分辨率和可信度。合成数据的反演算例也验证了所实现的大地电磁全信息资料三维共轭梯度反演算法的正确性和稳定性。     

大地电磁各向异性二维模拟及实例分析

经过半个多世纪的发展,国内外利用大地电磁法研究地球内部电性结构取得了令人瞩目的成就,这些研究成果多数是基于电性各向同性理论.然而地球内部普遍存在电性各向异性现象,地壳和上地幔中存在的电性各向异性是地电模型、地下结构和构造模型间一个重要的联系因素.本文首先由麦克斯韦方程出发,引入张量电导率,根据二维电性各向异性结构的特点,得到一组关于平行走向的电场分量E_x和磁场分量H_x的偏微分方程.使用有限差分法求解偏微分方程,求出E_x和H_x的近似解,并以此求得其它场分量;随后,通过对普通及特殊的二维电性各向异性结构做正演模拟,研究其对观测大地电磁场的影响,从而认识在普通及某种特定地质条件下的电磁传播特性,为其后对大地电磁实测资料的处理解释奠定理论基础;最后,以本文的研究成果为基础,将电性各向异性理论引入对实测大地电磁资料的处理解释中,通过对新疆某地的大地电磁资料做二维正演拟合解释,说明了电性各向异性现象的普遍存在,也验证了理论的正确性及算法的实用性,为今后分析解释大地电磁资料中的电性各向异性现象提供理论依据和技术指导,并开拓了对大地电磁实测资料处理的思路和方法.