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

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

基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演

各向异性介质对大地电磁观测数据的影响往往不可忽略,因此需要提高大地电磁各向异性三维反演的可靠性和有效性.为了满足大地电磁各向异性三维反演的需求,本文研究了一种基于交叉梯度结构约束的大地电磁主轴各向异性并行三维反演算法.根据大地电磁平面波理论假设,正演方程采用背景场与二次场分离的计算方式,二次场利用交错网格有限差分法求解.由于各向异性反演的多解性,本文将各向异性介质简化为主轴各向异性,并在此基础上进一步采用有限内存拟牛顿LBFGS法实现三维各向异性反演.为了提高各向异性反演的分辨率,反演目标函数中引入交叉梯度项,利用先验的结构信息,对三个方向的电阻率参数进行结构约束,最终的反演进一步利用MPI(Message Passing Interface,消息传递接口)技术实现分频并行计算,测试结果显示并行接近线性加速比.     

Are the crustal and mantle conductive zones isotropic or anisotropic?

One of the significant problems of modern deep magnetotellurics is the recognition of anisotropy in the crustal and mantle conductive zones. In the paper we perform numerical experiment comparing several 2D models of crustal and mantle isotropic and anisotropic prismatic conductors. Anisotropy is modeled by alternating horizontal or vertical thin layers of different resistivities (the vertical layers are parallel to the prism strike). Using these models, we examine conditions under which the magnetotelluric and magnetovariational response functions distinguish between isotropy and anisotropy. The resolution of MT and MV studies depends on the sediments conductance, lithosphere resistance and deep conductor width. Calculations show that the most favorable conditions for anisotropy studies are observed in the active regions characterized by small sediments conductance (10–20 S) and moderate lithosphere resistance (10⁸ Ohm·m²). However, in the stable regions, where sediments conductance exceeds 50–100 S and the lithosphere resistance comes up to 10⁹ Ohm·m², the crustal and mantle anisotropic and isotropic conductors manifest themselves in the equivalent magnetotelluric and magnetovariational functions, which cannot distinguish between anisotropy and isotropy and admit both the interpretations.     

IMPEDANCE OF A PLANE ELECTROMAGNETIC WAVE AT THE SURFACE OF A LAYERED CONDUCTING EARTH WITH DIPPING ANISOTROPY*

The problem of a plane wave incident on a non-isotropic dipping layer lying over an isotropic conducting substratum has been studied and some numerical results are presented to show the effects of

• 1) degree of anisotropy m,
• 2) conductivity contrast between the substratum to the upper layer b,
• 3) angle of inclination of the axis of anisotropy α,

on the variation of amplitude and phase difference of the wave impedance with the normalized skin-depth (d/δ₁).     

基于有限体积法的二维大地电磁各向异性数值模拟

为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法.     

基于中速-高速摩擦实验研究含碳断层带的电导率特征

野外地质调查结果显示,断层带常富集碳质.断层带中碳的分布结构是影响断层带电导率特征的一种重要参数.本文在室温、室内湿度和2MPa正应力条件下,对不同石墨含量(3,5,6和7wt%)的石英-石墨混合断层泥模拟样品开展了滑动速率介于500μm·s-1～1m·s-1的摩擦实验及相应的电导率测量,以期研究断层运动对碳分布结构的影响以及断层带电性特征对碳含量及分布的响应情况.结果显示,摩擦滑动能够显著地改变样品的电性特征(电导率大小及其各向异性).在平行滑动面方向(径向),样品电导率随着滑动位移的增加快速增加,在滑动约数十厘米之后,其电导率基本达到稳定状态;在垂直滑动面方向(轴向),样品电导率基本不随摩擦滑动速率和滑动距离而变化.SEM显微结构观测显示,摩擦滑动所引起的电导率各向异性直接反映了石墨分布结构的变化.该研究结果深化了对地震断裂带浅部电性特征的认识,为野外断层带大地电磁测深资料的解释提供了约束,同时对于了解含碳断层的力学性质和弱矿物相在剪切变形中的分布特征及其演化过程等方面也具有重要意义.     

Annual variations in the electromagnetic field of the earth and electrical conductivity of the geological environment

Synchronous annual variations in the geoelectric and geomagnetic field are studied on the basis of long-term electromagnetic monitoring. It is shown that the annual geoelectric variations have intraterrestrial origin and are not related to the annual geomagnetic variations. Temporal variations in the magnetotelluric impedance and magnetic tipper, which characterize the electrical conductivity of the geological environment, are analyzed. It is established that annual variations in the magnetotelluric impedance mainly describe the variations in the electrical conductivity of surface crustal layers and are less sensitive to the deep electrical conductivity of the Earth. The annual variations in the imaginary magnetic tipper at the periods of 1000–3000 s probably reflect the changes in conductivity of a deep transversal low-resistive zone (the fault). It is suggested that annual variations in the geoelectrical and geomagnetic fields, as well as in the electrical conductivity of the geological environment, arise as a response to the changes in the geodynamical processes caused by the revolution of the Earth around the Sun.

     

Electromagnetic measurements in the vicinity of the KTB drill site. Part II. Magnetotelluric results

Summary In addition to the magnetovariational measurements across an array in Western Bohemia, close to the KTB ultradeep borehole (Germany), discussed in part I of this paper [1], magnetotelluric results from pivot sounding point Ostrvek within the array are presented here. Good quality of long-period magnetotelluric data (period range from 30 s to about 1 hour) allowed structural dimensionality of the medium to be analysed in detail. The geoelectrical structure was identified as a slightly distorted two-dimensional regional substratum, with dominating E - W strike, overlaid by a heterogeneous subsurface layer with extremely strong and anisotropic galvanic distortion effect on the magnetotelluric data. Estimating the total static shift distortion tensor by fitting the local magnetotelluric curve to the curve of the global magnetovariational soundings (for the European continent), the static distortions were identified as of generally multidirectional origin. The resulting telluric ellipse is, however, strongly anisotropic, indicating an approximately SW - NE apparent local strike, which is in the approximate agreement with remote reference magnetovariational results. Finally, the magnetotelluric results from the station Ostrvek are compared with long-period data from the immediate neighbourhood of the KTB borehole on the German territory.     

用有限元法正演京津唐地区二维构造的大地电磁测深资料

本文阐述在大地电磁测深工作中,用有限单元法,模拟计算京津唐地区二维地质构造的大地电磁响应。 文章前部分,论述了用有限单元法正演二维大地电磁测深曲线的原理及其数学公式推导。文章后部分,根据两条测线实测的大地电磁测深结果并结合其它地球物理探测资料,给出了两个二维剖面模型,并计算了两种极化的视电阻率曲线。两条实测剖面分别是横切唐山地震带的马兰峪—柏各庄剖面及横过北京地区北北东构造走向的文化营—王辛庄剖面     

3D geoelectric model of the Mutnov steam hydrothermal deposit

The article discusses the results of the interpretation of magnetotelluric sounding (MTS) data generated in the vicinity of the Mutnov steam hydrothermal deposit during the period from 2004 to 2006. Approaches to an interpretation of magnetotelluric (MT) data in a situation characterized by significant 3D distortions and the coastal effect phenomenon are discussed. The resulting 3D geoelectric model is presented. Based on analysis of the created geoelectric model, the principal geoelectric horizons are characterized, which includes the parameters of the integral conductivity and the morphology of the geoelectric boundaries. Zones of the highest values of the electrical conductivity of the watered tuffaceous sediment horizon have been identified, as well as zones of the maximum lateral electrical conductivity contrasts that are interpreted as discharge zones of hydrotherms. The results of the drilling carried out based on MTS data are presented.     

Inversions for MT data in 2D symmetrical anisotropic media

Introduction Ready and Renkin (1971) were the first to make the research on anisotropy problems in magnetotellurics (MT). The progress in the research is not evident because it is more complex and difficult than isotropic problems. Now, the one-dimensional (1D) anisotropy problems in MT have been well solved, while for the two-dimensional (2D) cases, the numerical solutions have only been obtained for some particular conditions (Ready and Renkin, 1975). As to the three-dimensional (3D) ani…     

INDUCTIVE SOUNDING OF A STRATIFIED EARTH WITH TRANSITION LAYER RESTING ON DIPPING ANISOTROPIC BEDS *

The generalised three layer boundary value problem with a transition layer sand-witched between an isotropic overburden and dipping anisotropic substratum is discussed assuming that plane electro-magnetic waves are incident normally over the air-earth inter-face. The tangential electric (E_y) and magnetic (H_x) fields and the expression for surface-impedance (E_y/H_x) have been evaluated at the earth's surface. Through numerical analysis it is shown that changes in the values of the parameters m (coefficient of anisotropy), h (thickness of the transition layer), α (angle of inclination of the dipping beds), and b (conductivity ratio between substratum and upper layer) modify the amplitude and phase-variation curves (with skindepth) significantly.     

波西米亚块体西部长周期电磁研究(英文)

在波西米亚块体西南部做了一条长周期电磁剖面 ,它靠近德国超深钻的位置。用一个二维合成模型对数据做了大体的拟合 ,模型走向东西 ,近地表盖层高度地不均匀畸变。畸变层显示强烈的各向异性 ,优势电导方向为NW -SE至NNW -SSE。讨论了深部区域结构的几个不同的模型 ,其中一个模型可以解释由MT和GDS数据显示的主方向的偏差 ,还讨论了沿特别剖面得到的结果与区域地电结构的关系     

Comparative Magnetotelluric Modeling of Smooth 2D and 3D Conducting Bodies Using Rayleigh-Fourier Codes

—Recently, a method for 3D magnetotelluric modeling was developed, which is based on the application of the Rayleigh scattering theory. This method, RF-3D, is especially capable of modeling multilayered structures with smooth irregular boundaries. The formulation allows inclusion of layers with vertically anisotropic electrical conductivity.¶Using RF-3D, the response of smooth structures of practical interest is calculated and the importance of 3D effects is evaluated. Two models consisting of a 3D conductive body in the presence of a 2D shallow distortion are analyzed. In the first model, the direction of maximum elongation of the body is perpendicular to the strike direction of the 2D upper structure, and in the second one both directions coincide. In addition, the case of a small 3D shallow conductor over a regional 2D structure is also considered.¶3D effects are compared to those generated by 2D models with identical cross sections. In all the cases, the 3D responses differ from those of the 2D, especially directly over the bodies. A good agreement between the 2D transverse magnetic response and the corresponding components of the 3D response, along centrally located transverse profiles, is expected for elongate, prismatic conductors. Then, the differences obtained for the models considered in this study, particularly for the second and third models, are a consequence of the smooth geometry. They can be explained in terms of galvanic effects produced by boundary charges, which are greater near the vertical sides of a prism than on the sides of a body with smooth contours.¶Equivalent 2D models of the first and second structures are also obtained. In these models, the thickness of the conductor is underestimated, respectively, by about 30% and 24%.¶For the third model, when vertical anisotropy is analyzed, it is found that only the anisotropy of the first layer can be detected. This is because the effect of vertical anisotropy decreases strongly with depth and appears to be important only near the 3D anomaly.     

冻土区天然气水合物各向异性储层的AMT响应特征

陆域天然气水合物通常发育于冻土层下方破碎带和岩层裂隙处,其储层会表现明显电性各向异性特征.音频大地电磁法（AMT）能有效探测陆域水合物储藏范围,且不易受高阻冻土层的压制和干扰,可用于陆域天然气水合物探测研究.本文采用AMT对角各向异性二维正演方法,对多种陆域天然气水合物各向异性储层模型及相关参数进行了模拟试算,分析其视电阻率和相位响应特征.结果表明AMT能清晰显示各向异性储层空间位置和分布情况,并对储层的水合物饱和度变化以及储层数量等特征都有所反映.天然气水合物电性各向异性模拟研究为这一新型潜力能源的勘查工作提供了新的思路.     

Estimation of the correlation between magnetotelluric and geothermal data in the Bishkek geodynamic research area

The paper analyzes the correlation between the electrical conductivity and temperature in the upper crust of the Bishkek geodynamic research area (the Northern Tien Shan). Electrical conductivity profiles constructed from magnetotelluric data and thermograms from the boreholes near magnetotelluric sounding (MTS) points are used for estimations. The correlation analysis of conductivity and temperature profiles to depths of 3–4 km showed that, first, the correlation coefficients do not depend on the distance between the borehole and the nearest MTS point; second, the good correlation between the conductivity and temperature observed for the majority of borehole-MTS point pairs is accounted for by the fact that the study parameters vary with depth in a manner normal for laminated sedimentary rocks; and, third, a low correlation is due to specific features of the geological structure between the borehole and MTS point under consideration.     

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

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

用磁流源并矢Green函数的递推矩阵方法计算层状各向异性地层中多分量感应测井响应

采用递推矩阵方法计算各向异性介质的磁流源并矢Green函数,并利用上述并矢Green函数对层状各向异性倾斜地层中多分量感应测井的响应进行数值模拟,分析了线圈距、层厚、倾角和围岩对多分量感应测井响应的影响.计算发现,共面视电导率比共轴视电导率变化规律复杂,不能反映地层电导率的真实情况.各向异性地层共轴视电导率随井眼相对倾角的增加而减小,而共面视电导率随井眼相对倾角的增加而增加.仪器的垂向分辨率、围岩各向异性对目的层响应的影响程度取决于线圈距,围岩各向异性对高电导率目的层中测井响应的影响大于对高电阻率目的层中测井响应的影响.     

MAGNETOTELLURIC EFFECT OF DIPPING ANISOTROPIES*

Reddy , I. K., D. Rankin , 1971, Magnetotelluric Effect of Dipping Anisotropies, Geophysical Prospecting 19, 84–97. A solution is obtained for the magnetotelluric effect in the case of multi-layered earth, each layer containing dipping anisotropy, with the strike of anisotropy oriented at an arbitrary angle to the measuring directions. A technique to compute tensor impedances and tensor apparent resistivities is presented and the importance of anisotropies in magnetotelluric interpretation is discussed.     

Asthenospheric parameters of dipping plate regions by magnetotelluric investigations

The multifrequency resolving capability of the electromagnetic surface impedance parameter, employed in magnetotelluric investigations, has been examined for an upper mantle model incorporating (1) a deep-seated lithosphere-asthenosphere boundary and (2) preferred orientations of olivine crystals in the lithosphere, usually associated with the subducting plate boundary zones. Numerical results display the quantitative errors in the resolution of various anisotropic and dipping plate conditions. It is found that (1) conductivity anisotropy, particularly with higher dips, causes a significant dispersion in the impedance values and higher resolving power, and (2) variations in the geometrical/physical parameters are reflected predominantly in the phase component of the surface impedance. The study has relevance to the understanding of the electrical nature of materials, its relationship with physical properties and associated geothermal and seismic activities in the dipping plate regions.