基于三维正演的音频大地电磁阻抗相位不变量校正技术

复杂地形、地质条件的大地电磁数据解释容易出现假象,采用三维正演技术模拟地形和地表不均匀体的背景响应,对实测数据阻抗相位不变量进行校正,实现更准确的定性分析;对三维异常体模型的合成数据进行一维、二维多参数反演试算,以确定地形剧变区选择反演技术的最佳方案.合成数据的试反演结果显示一维反演水平切片假异常较多,二维反演能压制测向假异常,但不能压制走向的假异常,水平切片多出现测向条带.使用本文提出的阻抗相位不变量校正法扣除地形、地表背景响应,结合一维、二维反演,能使实际资料解释成果更加可靠.     

海岸效应对近海地区大地电磁测深数据畸变作用研究

在近海地区采集的大地电磁测深数据通常受到海岸效应的影响,使得大地电磁测深数据发生畸变,因而很难利用大地电磁测深资料较为可靠地获得地下深部的电性结构.本文通过正演模拟方法,分析和总结海水深度变化和海底地形变化对近海地区大地电磁测深数据的畸变影响.当测区与海岸线的距离小于目标频率的大地电磁场趋肤深度时,高导海洋的存在会严重影响测区内电磁场的分布.由于海岸效应的影响,大地电磁测深视电阻率曲线和相位曲线均会发生不同程度的畸变,在低频部分,这种畸变作用尤为明显.大地电磁测深一维Occam反演方法和二维非线性共轭梯度反演方法,对近海地区浅部地层具有较好的反演效果.随着海水深度的增加和海底地形的复杂变化,两种反演方法均会出现不同程度的假异常,为地质解释工作造成了影响.近渤海地区的实测大地电磁测深数据在低频部分可能受到海岸效应的影响而导致视电阻率曲线的严重畸变.     

起伏地形下大地电磁L-BFGS三维反演方法

为推进大地电磁三维反演的实用化,本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法,然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式,解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略,提升了反演的稳定性.利用开发的算法对多个带地形地电模型（山峰地形下的单个异常模型、峰-谷地形下的棋盘模型）的合成数据进行了三维反演,并与已有大地电磁三维反演程序（ModEM）进行对比,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了华南某山区大地电磁实测数据,得到该区三维电性结构,揭示了研究区以高阻介质为基底,中间以低阻不整合面和相对低阻介质连续分布,浅部覆盖高阻介质的电性结构特征,进一步验证了本文算法的实用性.     

江西相山火山盆地地质结构研究:来自大地电磁测深及三维地质建模的证据

江西省相山火山盆地是我国第一大、世界第三大火山岩型铀矿田,其西部牛头山一带铀矿勘探中发现深部有大垂幅的Pb-Zn-Ag矿化.60多年来,以该矿田为对象开展的研究取得了一系列丰硕成果,但对火山机构的认识仍不确定.我们采集了涵盖该火山盆地主要地质体的1386块钻孔岩芯标本和243块地表岩石标本,开展了电阻率、磁化率、密度等物性参数测量,并在火山盆地中施测了19条MT剖面(2条骨干剖面和17条精细剖面),对3000 m以浅主要地层、岩体和断裂带等目标地质体的三维展布特征进行了解译和三维建模.研究结果表明:(1)相山火山盆地具有变质岩-花岗岩双基底.基底变质岩系顶界面表现为南北分带(三隆间两凹)和东西分块(两垒夹一堑)的三维地质格局;南西部有加里东期花岗岩侵人,具有似层状的空间展布特征;盆地基底变质岩系与上覆火山-沉积岩盖层之间呈连续的水平低阻异常带,不整合界面清晰.(2)打鼓顶组火山岩呈似层状产出,主要分布于盆地西部;在河元背一船坑一杏树下一带识别出近东西走向厚层的流纹英安岩凹槽,相山铀矿田西部探明的主要铀矿床分布在该凹槽内或其边缘.鹅湖岭组火山岩总体形态呈蘑菇状,在盆地中部厚度较大.在相山主峰半径约2 km的范围内,发现自下而上贯通式的低阻异常,推测是鹅湖岭组碎斑熔岩喷发的通道相(火山颈相),其火山颈呈陡立管状,深部向南东倾伏,浅部向南东撒开.后期花岗斑岩呈岩墙-岩床组合状,总体构成一个向西开口的环形岩体.打鼓顶期主要岩浆通道位于相山顶一巴山之间,次岩浆通道位于河元背;鹅湖岭期火山活动主岩浆通道也位于相山顶一巴山之间,次岩浆通道位于河元背、阳家山(芙蓉山)、严坑和柏昌.(3)火山盆地中断裂构造发育,MT测量结果显示存在7条北东向、4条北西向和1条南北向格架性断裂构造(其中一条新发现的北东向断裂隐伏于白垩纪红盆之下),盆地北部发育1条弧形火山塌陷构造,表现为大规模延续的低阻异常带.     

Structural mapping in basin-and-range-like geology by electromagnetic methods: a powerful aid to seismic

A case history is presented where electromagnetic (EM) methods were applied as a complement to seismic, for structural mapping in basin-and-range-like geology: 366 five-component magnetotelluric (MT) soundings were carried out together with 331 transient soundings (TDEM) along seismic lines. Due to high structural complexity, seismic shows a number of limitations. For the same reasons, MT is highly perturbed and three specific interpretation techniques were comprehensively applied: 1. a classical correction of static effect using TDEM sounding, to determine the high-frequency nondistorted apparent resistivities and thus the corrected tensor; 2. a so-called regional correction based upon the same concept as the static effect, to transform distorted resistivity curves due to the horst/graben situation into plausible 1D curves, through the use of nomograms built for 2D H-polarization situations; 3. a stripping technique which made it possible to map areas where a deep conductive Mesozoic shale was present below carbonates, at a depth of 3 km. After the best MT interpretation was obtained along each line, it was integrated with seismic and with the results from two boreholes. A crude empirical law relating resistivity and acoustic velocity was established and the MT horizons were plotted on the two-way traveltime seismic sections. The final integrated cross-sections obtained are undoubtedly of greater use to the explorationist than the initial seismic sections alone and two wells were accurately predicted.     

松嫩地块东缘和佳木斯地块西缘电性结构

横过松嫩地块东缘和佳木斯地块西缘的大地电磁测深剖面揭示了两块体结合带附近的深部电性结构.本文对剖面测点做了标准化数据处理,并对二维偏离度、构造走向进行了计算和分析,采用非线性共轭梯度(NLCG)二维反演方法对TM模式的数据进行了反演,获得了该剖面的地壳、上地幔电性结构模型,划分出三个典型构造单元:松嫩地块东缘、碰撞拼合带和佳木斯地块西缘.研究结果表明,研究区上地壳基本呈高阻特征,可能为岩浆岩,代表其经历了多期次岩浆作用,而松嫩地块东缘和佳木斯地块西缘的中下地壳的高导体可能与地幔物质的上涌有关;拼合带下方存在西倾的高导体和高阻体,可能是佳木斯地块向西俯冲到松嫩地块下方的构造遗迹;研究区可能发生了拆沉作用,与之伴随的地幔物质上涌可能是后期伸展作用的一个动力.     

Mapping deep targets based on integrated 3D MT-gravity interpretation: a case study

Mapping deep geological hydrocarbon targets is of significant importance in basin exploration. In areas lacking reliable seismic data, magnetotelluric (MT) and gravity explorations are helpful to delineate the distribution of potential deep geological hydrocarbon targets. Here we investigate the effectiveness of the integrated 3D MT and gravity explorations for mapping the potential deep hydrocarbon source rocks. The result based on the data from the W Basin (part of the Ordes Basin) of China demonstrates that the method is efficient and economical for basin exploration. The method is particularly useful in target areas which are of great interest for oil and gas exploration but lack high quality seismic data. In our method, we first use the high-precision 3D small-bin MT data acquisition to improve the data accuracy. Then we perform datum static correction method and apply 3D inversion to obtain the3D resistivity distribution. We also develop a layered resistivity model based on resistivity logging to assist the interpretation of the inverted 3D resistivity data so as to derive an initial 3D geological model. Starting from the initial model, we use 2D gravity data to update the model via 2D inversion line by line, and then pass the updated model for the next round of the 3D MT inversion. The integrated inversion is implemented iteratively so the model converges to satisfy the need of final geological analysis. The application to the W Basin shows that we could successfully delineate the geological distribution of the potential deep hydrocarbon source rocks within the basin and map the thickness of the upper Paleozoic.     

大地电磁局部畸变异常特征三维数值模拟研究

局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。      

基于矢量有限元的带地形大地电磁三维反演研究

研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.     

大地电磁资料精细处理和二维反演解释技术研究(三)——构建二维反演初始模型的印模法

在主流的线性最优化大地电磁二维反演中,如何合理构建初始模型是一个亟待解决的问题.常用的是采用均匀半空间或一维反演结果构建初始模型,不易获得稳定可靠的反演效果.实践表明,尽管基于不同初始模型的大地电磁二维反演结果差别较大,但均较初始模型更为接近真实模型.基于这样一种认识,经过反复的理论和实践探索,我们提出构建大地电磁二维反演初始模型的印模法.印模法的基本思想是依据已有反演结果和均匀半空间模型之间的加权来确定下一步二维反演的初始模型,它一方面保留了已有反演结果中关于真实模型的宏观轮廓信息,另一方面,保证了深部电性结构的均匀性,从而满足大地电磁二维正演所要求的底边界条件.基于印模法,本文进一步提出了迭代重构的反演思想.通过多个理论模型和实测数据的反演计算,验证了上述方法可在很大程度上压制初始模型对反演结果的影响.     

稀疏测线大地电磁资料三维反演研究:合成算例(英文)

受勘探成本和工区环境等因素的影响,当前大多数大地电磁实际工作采取布置稀疏测线采集数据和使用二维反演方法解释这些稀疏测线数据的方式。然而,二维反演方法在解释三维地电构造数据时存在局限性,有时甚至做出错误的地质解释。本文尝试了使用三维反演方法对大地电磁稀疏测线数据进行反演解释。使用大地电磁全信息资料三维共轭梯度反演程序对理论模型合成稀疏测线数据进行了三维反演。结果表明:这种反演方案是可行与有效的。同时,我们发现在不同数据的三维反演结果中,四个张量阻抗元素和两个倾子数据同时反演的结果相对更为准确,更接近理论模型。     

SOTEM数据一维OCCAM反演及其应用于三维模型的效果

本文基于垂直磁场分量研究了SOTEM数据的一维OCCAM反演方法,并将其应用于理论三维数据及野外实测数据的反演.对于大部分一维模型,OCCAM反演可取得较好的反演效果,且反演结果不依赖于偏移距;噪声对SOTEM数据的OCCAM反演具有较大影响,但当信号含噪水平不超过5%时,反演结果仍具有较好的准确性;若浅层存在较厚的低阻层,OCCAM反演结果对下部地层的分辨能力下降,仅能获得具有平均效应的电阻率.将一维算法应用于SOTEM三维数据的反演,会产生较大的误差,尤其是在异常体边缘地带影响最为严重.该影响程度与异常体和背景电阻率之间的差异有关,对于大多数电性近似呈连续变化的真实大地而言,一维OCCAM反演算法仍可获得较好的效果.最后通过陕西某煤田深部富水性调查的实测SOTEM数据反演验证了本文的研究成果.     

阿尔山火山区地壳上地幔电性结构初探

在阿尔山活火山区一条北北西向测线上进行了7个大地电磁测深点的观测。2-D解释结果表明,研究区内新、老两条火山带可能存在通往深部的岩浆通道。其中,新发现的活火山带地下在10~12km还保持着较高的热状态,很可能富含流体,在30~50km处可能是地幔向上的供热通道;而老火山带30km以上可能存在一个正在冷却的岩浆通道,两条火山条带的深部可能是同源的。     

The 3D profile inversion of the MTS–MVS data by the example of the Shiveluch volcanic zone in the Kamchatka Peninsula

The joint use of the magnetovariational (MV) and magnetotelluric (MT) data and the a priori information about the areal distribution of the deep electrical conductivity within the studied area and the integral conductance of the sedimentary cover enabled us to conduct a three-dimensional (3D) inversion of these data and construct the 3D model of the northeastern part of Kamchatka. The results show that within the Shiveluch volcano, north and south of it, there are northwesterly striking conductive faults penetrating up to the mantle which ascends from a depth of 41 to 35 km. The electrical conductivity of these crustal structures which accommodate high seismic and volcanic activity is probably accounted for by fluid saturation.     

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.     

Distortions of the electromagnetic field by shallow basins and by resistive outcrops

The interpretation of magnetotelluric (MT) measurements carried out on shallow (several hundred metres deep) basins and on the surrounding highly resistive rock outcrops can be difficult due to different forms of distortion, mainly 3D effects. A good example for study of this problem is the Mór Graben, which is the transition zone between the Bakony and Vértes Mountains (W-Hungary), where more than 20 MT soundings have been made with stations spaced about 2 km apart to determine the structure of a conductivity anomaly at a depth of about 3–4 km.

The statistical treatment of the different distortions due to varying sediment cover (the S-effect or static shift) and to the steep wall of the resistive basement outcrops (edge effect), etc., enabled the estimation of the actual parameters of the conductive formation. The interpretation of the field data is supported by numerical modelling.     

APPLICATION OF 3D HIGH PRECISION GRAVITY METHOD TO THE EXPLORATION OF DEEP BURIED HILL

Due to the influence of volcanic rock or weathering crust coverage, the quality of deep seismic data is often not reliable and it is difficult to identify small deep buried hill structure beneath the cover. However, the gravity exploration technology can use the remarkable density differences between the object and the overlying strata to identify this special geological reservoir. Although recently several rounds of seismic exploration have been carried out in Raoyang depression, Huabei Oilfield, North China to determine the existence of Hu 8 north deep buried hill, whether the buried hill really exists or not as well as its scale is still in controversial. In this paper, based on the information of seismic data and formation density, deep processing of 3D high precision ground gravity data has been carried out for 3D forward and inversion computation. The dimensional gravity forward calculation results show that the density of the deep anomaly body forms a relatively low gravity anomaly in the earth's surface. By using the potential processing method of vertical second derivative, and sliding filtering, the residual anomaly is separated from regional or background field, which illustrates the existence of a local high gravity anomaly at Hu 8 north area. According to the amplitude of residual gravity anomaly and formation density difference modeling, through a number of 3D forward calculations and 3D inversion of gravity data, the existence of Hu 8 north buried hill and its possible scale are basically determined. The results prove that the 3D high precision gravity method is effective in determining the deep buried hill structure in case that seismic data is not reliable.     

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

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

Internal structure of Ustica volcanic complex (Italy) based on a 3D inversion of magnetic data

A ground magnetic study of Ustica Island was performed to provide new insights into subsurface tectonic and volcanic structures. The total-intensity anomaly field, obtained after a data-reduction procedure, shows the presence of a W–E-striking magnetic anomaly in the middle of the island and another two intense anomalies, which seem to continue offshore, in the southwestern and the northeastern sides, respectively. The detected anomalies were analyzed by a quadratic programming (QP) algorithm to obtain a 3D subsurface magnetization distribution. The volcano magnetization model reveals the presence of intensely magnetized volumes, interpreted as the feeding systems of the main eruptive centers of the island, which roughly follow the trend of the main regional structural lineaments. These findings highlight how regional tectonics has strongly affected the structural and magmatic evolution of the Ustica volcanic complex producing preferential ways for magma ascent.     

Deep structure of the northeastern margin of the Parnaiba Basin, Brazil, from magnetotelluric imaging

The magnetotelluric (MT) method has been applied to the determination of the deep resistivity structure of the northeastern margin of the Parnaiba Basin. Transient electromagnetic (TEM) and MT data were collected in early 1999 along a 95 km long N–S line, extending from the coast across the projected subcrop position of a discontinuous fault found to the west of the study area that is believed to be a possible basin‐bounding fault. The MT data were processed to yield the TE‐ and TM‐mode responses and then corrected for static shift using central‐loop and single‐loop TEM data, respectively. Regularized 2D MT inversion was subsequently undertaken using a structured initial model with the near‐surface constrained by TEM inversion results. As a consistency check, we performed another set of 2D inversions using different smooth initial models. The various optimal 2D inversion models show clearly the presence of a major basement trough, over 2 km deep, located about 70 km from the coast. We interpret it as possibly marking the main basin margin and suggest that it may have implications for groundwater resource development in the area.