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

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

考虑电流型畸变的大地电磁三维反演

大地电磁测深(MT)的观测数据易受到由近地表小尺度非均匀体或地形起伏引起的电流型畸变干扰,消除或压制这种干扰对获取可靠的深部电性结构至关重要.当区域结构为二维时,电流型畸变可采用张量分解等方法予以消除或压制.当区域结构为三维时,畸变问题更加复杂和严重,传统张量分解方法往往效果不佳或无效,严重地制约了MT三维反演技术的实用性.对此,本文提出一种考虑电流型畸变的MT三维反演算法,将完整的电流型畸变参数引入到目标函数,并采用非线性共轭梯度法与电阻率参数同时反演,从而达到压制畸变的目的.该算法有两个关键点:一是通过分析实测数据所遭受畸变的分布特征,在目标函数中对其进行有效约束;二是在迭代过程中,通过自适应地调整双正则化因子保障算法的稳定和效率.理论模型测试结果显示,常规三维反演算法不能合理解释数据中的畸变成分,而只能通过引入虚假异常体强制地拟合受畸变数据,从而造成电阻率模型严重失真.与之相比,本文算法能够在反演中自动求解各测点所受到的畸变,获得更接近真实的电阻率模型.     

Numerical modeling of 3-D terrain effect on MT field

Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special net division. The result of modeling test for 2-D terrain by this modeling technique is basically coincident with that by 2-D modeling technique, but there is a great difference between the results of 3-D and 2-D modeling for 3-D terrain. Project supported by the National Natural Science Foundation of China.     

Gravity Anomalies of 2.5-D Multiple Prismatic Structures with Variable Density: A Marquardt Inversion

We present an inversion technique based on the Marquardt algorithm to estimate the depth of a 2.5-D sedimentary basin in addition to the regional gravity anomaly that is associated with the residual gravity anomaly, wherein the density contrast varies parabolically with depth. Forward modeling is carried out through a derived analytical gravity expression of a 2.5-D vertical prism. Inversion of a theoretical gravity anomaly with and without a regional gravity anomaly illustrates the procedure that it is found to be insensitive to the regional gravity effect. Furthermore, the algorithm is exemplified with the gravity anomalies of the derived density-depth model of the Godavari subbasin, India with a parabolic density profile resulting in a more consistent geological model rather than a constant density profile. The main advantage of this method is that it works well even when the profile of interpretation does not bisect the strike length of the sedimentary basin.     

庐枞矿集区矾山-将军庙地区AMT三维反演及地质结构解释

庐江-枞阳矿集区是长江中下游成矿带内的一个重要的多矿种成矿区.为揭示矾山-将军庙地区的地质结构以指导找矿,首先,在庐枞矿集区北部的矾山-将军庙地区开展了三维音频大地电磁数据采集工作,采用Rhoplus方法进行AMT数据死频带校正保证数据质量,利用相位张量分析地下介质的维性特征;其次,利用模块化三维反演并行代码ModEM开展了实测数据的三维反演,获得了一个典型的具有双层结构特征的火山岩盆地三维地下电性模型;最后通过电性模型及其他地质资料的综合解释,依据反演模型中近地表的低阻层和深部的高阻隆起区,分别勾画出了火山岩地层的厚度和深部侵入岩的分布范围,并在小岭地区发现地表浅部下方存在一个含矿次级火山机构的高导体.     

THE IMPORTANCE OF TOPOGRAPHIC CORRECTIONS ON MAGNETOTELLURIC RESPONSE DATA FROM RUGGED REGIONS OF ANATOLIA

Topographic irregularities cause some distortions of magnetotelluric (MT) fields. In the vicinity of a topographic feature, the TM-mode distortion increases with the height and inclination of the slope. It is well-known that TM-mode ( E _⊥) topographic effects are much greater than TE-mode ( E _ı) distortions.   We have made a study of MT anomalies in TM-mode due to two-dimensional topography. In order to reduce these effects, the distortion tensor stripping technique was used. After corrections, the resulting data can be interpreted as if they were obtained over a flat surface and depend only on the subsurface structure. However, this technique sometimes causes some geometrical distortions of the real subsurface structure.   One of our aims is to overcome this failure. We have modified the correction coefficients by considering the actual one-dimensional geology. Model studies showed that our approach is especially useful in removing the terrain effects on complex 2D subsurface structures.   The other purpose of this study is to emphasize the importance of a proper terrain correction for data from sites having mountainous topography over complex geology, e.g. strike-slip faults, suture zones and rift valleys. Some examples of MT data sets collected from the North Anatolian Fault Zone and from the thrust regions of the Western Taurides will be presented.     

大地电磁二维层状模型参数化反演

本文针对目前二维反演存在求解参数过多和求其偏导数困难的两个问题，在近年来关于连续介质二维大地电磁测深资料快速反演的理论基础上，对二维层状模型反演的目标函数，进行线性近似，采用加权法来近似计算二维模型参数的偏导数。从而实现比较快速的二维层状模型的参数化反演。理论模型的反演试验表明，反演迭代可以较快地以足够的精度拟合视电阻率和相位资料，并收敛接近理论模型。     

宁洱地震区深部电性结构及发震构造初析

在2007年6月3日宁洱6.4级地震区进行了宽频带(1000～1/4096Hz)大地电磁观测,采用了远参考道技术和Robust技术对观测数据进行了处理,分析了视电阻率和相位曲线、二维偏离度及区域电性结构方向,采用非线性共扼梯度(NLCG)二维反演技术对实测资料进行了反演,得到了与实测曲线拟合程度较高的反演模型。从二维解释结果可以看出:横向上,以12号测点为界,明显地划分为2个不同的电性分区,推测存在具一定延伸的断裂构造。纵向上,30km深度以上大致可分为4个电性层:1)中、新生代盆地沉积层(Mz—Cz),厚约1km。2)古生代沉积层(Pz),厚约5km。3)壳内高导层(HCL)。4)中、下地壳(MLC)。其中壳内高导层埋深在12号点左右存在较大的差异,推测在12—13号点下高导层顶部存在壳内的滑脱构造。宁洱地震震源深度为5km,文中所推测的滑脱构造位于这一层位上。认为该滑脱构造为本次地震的发震构造     

Electromagnetic excitation of asthenospheric conducting zones

The relations between galvanic and inductive electromagnetic excitation of 3-D asthenospheric conducting zones are considered. Mathematical modeling with the application of integral equations shows that the rather intense magnetotelluric anomalies reflecting the topography of the conducting asthenosphere arise via the galvanic channel, which has two branches: (1) the lithospheric branch (an electric current flows from the sedimentary cover into the lithosphere and concentrates in uplifts of the conducting asthenosphere), which is inactive in the case of a very high resistivity of the lithosphere (the screening effect), and (2) the mantle branch (the deep mantle current flows from below into uplifts of the conducting asthenosphere), which is free from the screening effect of a high resistivity lithosphere. In previous works, the mantle branch of the galvanic channel was erroneously associated with local induction in uplifts of the asthenosphere. According to modern estimates, local induction in the asthenosphere generates weak magnetotelluric anomalies that cannot be reliably measured.     

基于非结构化网格的任意复杂2D RMT有限元模拟

射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性.     

用边界单元法模拟二维地形对大地电磁场的影响

本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定.     

On Two-Dimensional Modeling Of Magnetotelluric Field Data

In this paper, some recent topics on the modeling of magnetotelluric data are introduced. The focus is on the handling of real field data for two-dimensional resistivity modeling. First, the removal of the effects of near surface heterogeneity is reviewed. It covers telluric distortions (phase mixing and static shift) and magnetic distortions using conventional Groom-Bailey type 3D/2D model (three-dimensional local anomaly underlain by regional two-dimensional structure). The extension of a 3D/2D distortion model for multi-site, multi-frequency is a new development. Magnetic distortion seems to be less significant for land observations, but significant for sea floor data, where the regional magnetic field is weak due to seawater. In special cases involving for example, distortion due to topography and bathymetry, explicit removal is possible. There are some schemes proposed for a 3D/3D model (three-dimensional local anomaly underlain by regional three-dimensional structure). Along with the removal of the distortion, it is important to recognize the dimensionality of the dataset prior to modeling. A property using strike estimates for each site is an indicator for dimensionality which is unaffected by local distortion. Mapping the local strike or a rose diagram is an effective visualization of the dimensionality.Two-dimensional inversion is becoming routine. For the fast calculation of derivatives, approximate calculation, reciprocity or conjugate gradient methods are used. In order to incorporate a priori information and to overcome the intrinsic ill-posed nature of the inversion problem, imposing constraints on the model structure is important. A proper tradeoff between the data fit and constraints should be optimized to obtain minimum structures that are required by the field data. However, the choice of constraints is rather subjective and depends on the geological situations. For field data, two-dimensional inversion has limits on modes, area, and period range. Special care must be taken for the structure outside the profile. Two-dimensional inversion incorporating anisotropy is interesting and becoming popular, but the structure may not be unique. Future development in three-dimensional inversion for real datasets should take the above points into consideration.     

用可视化技术分析虚拟现实环境下三维地震数据场的研究

目的 用可视化方法分析三维地震数据，直接观察数据中隐藏的地质现象。方法：将第每个反射子波用一个粒子模型表示，来自同一反射面的子波粒子相互叠合成一个同相面。为加快显示速度，将一个三维地震数据体转换为三维二值图像，用有序扩展算法追踪同相面所在的三维连通域，细化后得到反射面的坐标数据。再在虚拟现实系统下的对反射面进行立体显示，直观地分析地下地层的结构。结果 实际数据试验得到了反射面的三维形态及断层的分布规律，立体显示更深刻地揭示了复杂地质现象及其规律。结论 只有在虚拟现实系统中进行三维地震数据的可视化，才能直观地看到地下地层的真实场景，真正实现三维数据的全三维解释，所提出的方法是有效。     

An audio-magnetotelluric investigation in Terceira Island (Azores)

Ten audio-magnetotelluric soundings have been carried out along a profile crossing the Serra do Cume caldera in the eastern part of the Terceira Island (Azores). The main objectives of this investigation were to detect geoelectrical features related with tectonic structures and to characterize regional hydrological and hydrothermal aspects mainly those related to geothermal fluid dynamics.Three-dimensional numerical investigation showed that the data acquired at periods shorter than 1 s are not significantly affected by ocean effect. The data was analysed using the Smith's decomposition method in order to investigate possible distortions caused by superficial structures and to estimate a global regional strike. The results suggest that in general the soundings were not distorted. A regional N55°W strike was chosen for the two-dimensional data inversion.The low-resistivity zones (10–30 ohm-m) displayed in the central part of the 2-D geoelectrical model have been interpreted as caused by hydrothermal circulation. The low-resistivity anomalies at the ends of the profile might be attributed to alteration zones with interaction of seawater intrusion. High-resistivity (> 300 ohm-m) values have been related with less permeable zones in the SW of Cinco Picos and Guilherme Moniz caldera walls.     

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.     

昌马地震应力场及发震机制研究

本文利用前人大量的地质、地球物理资料，建立了昌马盆地的立体地质模型，利用三维有限元法研究了１９３２年昌马７．６级地震前后应力场的变化及其昌马地震成因机制，结果表明：（１）震源断层走向与利用主压应力推测的断裂走向相差６５°，可能是由于地壳浅部旋转的主压应力与深部北东方向的主压应力所产生的力矩使块体发生旋转所致，其枢纽点即为震源所在；同时，该旋转也导致了地震断层和震源断层在力学性质和几何性质上的不同。（２）昌马地震是在闭锁、贯通和块体旋转联合作用下形成的。（３）区域应力场与局部应力场不一致、应力松弛单元的出现可能与地震的孕育、发生有关。（４）断裂活动的不均一性与平面最大剪应力分布不均匀成正相关，且各断裂均以相对左旋走滑兼挤压为特征，但是其走滑量的分布是不均匀的。     

三维精细构造解释的方法流程和关键技术

详细论述了三维精细构造解释的方法和流程.对层位标定、层位解释、断裂、速度解释及成图每个步骤的方法给以具体论述,重点分析了三维精细构造解释使用的关键技术如可视化技术、相干体技术、变速解释技术和真三维解释技术.本文可为从事三维精细构造解释的人员提供借鉴.     

三维地震勘探数据处理的问题及其解决方法

三维地震勘探在十年内取得了显著的地质效果和经济效益。三维地震数据处理技术也得到了迅速的发展。但是,与二维相比,三维地震数据处理技术有许多不完善之处。本文对其中几个主要问题进行分析并提出解决方法。这些问题是三维叠加速度分析、三维静校正和叠前三维偏移等处理方法。     

Electromagnetic profiling interpretation across vertical faults and dikes

When a vertical fault or dike scatters a normally incident TE-mode plane wave, the horizontal components of the electric and the magnetic fields vary along the direction perpendicular to the strike. This scattering is also responsible for the formation of a vertical component of the magnetic field, which also varies along the direction perpendicular to the strike. Analysis of the lateral variations of the field components permits the identification of the type of geological structure, either fault or dike, as well as an estimate of the position of the geological contacts and the electrical conductivity of each medium. Previous exact analytical solutions have shown that two Fourier cosine integrals can express each field component. A series of functions specified for each model represent the kernel of each integral. From the field components obtained from the first non-zero term of each series, we calculated the following functions: the dip angle and the ellipticity of the vertical polarization ellipse in a plane perpendicular to the strike of the structure, and the azimuth and the ellipticity of the horizontal ellipse. We established master-curves of these functions for the interpretation of vertical faults and dikes for the polarization ellipsoid, for instance VLF or audio-frequency methods. This representation has as variable the induction number     , and as parameters σ ₂/ σ ₁ and     , where a denotes the half-thickness of the dike. The interpretation procedure using curve fitting is possible because the induction number is represented on a logarithmic scale. This procedure represents an important step before automatic interpretation is carried out. Two case histories using field data illustrate the effectiveness of the procedure and the type of quantitative interpretation obtained from it.