共查询到20条相似文献，搜索用时 187 毫秒
1.
局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。 相似文献
2.
大地电磁测深(MT)的观测数据易受到由近地表小尺度非均匀体或地形起伏引起的电流型畸变干扰,消除或压制这种干扰对获取可靠的深部电性结构至关重要.当区域结构为二维时,电流型畸变可采用张量分解等方法予以消除或压制.当区域结构为三维时,畸变问题更加复杂和严重,传统张量分解方法往往效果不佳或无效,严重地制约了MT三维反演技术的实用性.对此,本文提出一种考虑电流型畸变的MT三维反演算法,将完整的电流型畸变参数引入到目标函数,并采用非线性共轭梯度法与电阻率参数同时反演,从而达到压制畸变的目的.该算法有两个关键点:一是通过分析实测数据所遭受畸变的分布特征,在目标函数中对其进行有效约束;二是在迭代过程中,通过自适应地调整双正则化因子保障算法的稳定和效率.理论模型测试结果显示,常规三维反演算法不能合理解释数据中的畸变成分,而只能通过引入虚假异常体强制地拟合受畸变数据,从而造成电阻率模型严重失真.与之相比,本文算法能够在反演中自动求解各测点所受到的畸变,获得更接近真实的电阻率模型. 相似文献
3.
本文对大地电磁观测阻抗实施一种数学变换--"共轭阻抗变换",发现转换后的观测阻抗和区域阻抗之间存在特定的关系,这种特定关系不受电场局部畸变的影响,而且同样不需要关于地下区域结构维性的假设,在区域结构是三维的情况下也是成立的.对转换后的观测阻抗采用已有的Swift旋转方法即可求得区域主轴方位角,然后采用最优化方法求取区域阻抗相位、振幅以及畸变因子.同时根据转换后的观测阻抗重新定义了不受电场局部畸变影响的构造维性参数.采用合成理论数据验证了新算法的正确性,和Swift、Bahr、GB、相位张量、WAL方法进行了对比分析,并将新方法应用于实测资料的解释,发展了一套消除局部畸变和进行构造维性分析的MT精细资料处理技术.进一步的工作是将新算法推广到多点多频,并采用最优化技术分解,发展一种稳定性好、多测点多频点的阻抗张量分解技术. 相似文献
4.
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. 相似文献
5.
Yu. F. Moroz T. A. Moroz V. A. Loginov A. G. Nurmukhamedov D. A. Alekseev 《Izvestiya Physics of the Solid Earth》2016,52(1):29-44
The results of the magnetotelluric (MT) soundings before and after an earthquake are analyzed. The interpretation is based on the longitudinal and transverse MT curves (along and across the strike of the main tectonic elements, respectively). The MT curves are distorted by the ρ- and coast effect. The distortions due to the coast effect are estimated by the testing three-dimensional (3D) model. It is established that the coast effect distortion at the periods up to 1000 s is small and can be disregarded. The divergence of the longitudinal and transverse MT curves, which points to the presence of the deep faults, is thoroughly studied. The inversion of the MTsounding curves is carried out by the REBOCC program of the numerical two-dimensional modeling. This program implements the procedures of elimination of the ρ-effect and the joint inversion of the longitudinal and transverse MT sounding curves. The obtained geoelectrical cross sections provide an insight into the structure of electrical conductivity of the lithosphere before and after the earthquake. The more intense variations in the electric conductivity are observed in the zone of the deep faults. These variations are related to the changes in the porosity and saturation of the rocks by the highly mineralized fluids. 相似文献
6.
大地电磁测深二维反演方法求解复杂电性结构问题的适应性研究 总被引：3，自引：1，他引：2
为探讨二维反演方法在三维电性结构中的适应性问题,本研究中设计了一系列的二维/三维正演模型进行实验计算,分析了三维高阻/低阻异常体对模型响应的畸变作用,并从反演模式选择和数据旋转方向两个方面进行模型二维反演的对比分析,与三维反演的结果进行了比较,最后采用了实测数据进一步进行了二维和三维反演的比较实验.实验和研究结果表明,在剖面选择方面,在剖面方向与垂直主构造方向相差不大的情况下,截取剖面方向,将电性主轴旋转到垂直剖面方向的二维反演结果与垂直主构造方向的反演结果都可以较好地还原正演模型,在大的构造的反映上并无太大差异.在地下为二维或近三维条件时,正演模型的主要结构都可以较好地被二维和三维反演解析出来.二维的反演结果可能甚至会比三维的反演结果的边界更清晰,更精确.然而,对于具有较强的三维结构的模型而言,其二维反演结果与原始模型可能仍然存在较大差异,其中TM+TP或TM模式的二维反演结果相对更接近原始模型,而TE模式的结果往往会有较大误差,需要在解释时特别注意以免得出错误结论. 相似文献
7.
The Magnetotelluric Phase Tensor: A Critical Review 总被引：1，自引：0，他引：1
John R. Booker 《Surveys in Geophysics》2014,35(1):7-40
8.
提出了一种简单的大地电磁阻抗张量畸变分解方法，引入3d/2d模型，使得全畸变的大地电磁阻抗张量被生意人经，通过对比Swift旋转和Bahr分解后的资料，用联系二的畸变矩阵表示浅层三维电性异常体对区域二维构造的局部畸变影响，通过研究畸变张量矩阵，并对其进行分解，可将局部三维电性不均匀体产生的畸变与区域构造分离开来，从而了解区域构造的维数特征和局部畸变类型及其对区域构造的影响，最后，应用3个测点的实际资料，对这种简单的畸变分解方法进行了检验，结果表明，该方法简单且应用方便，解释效果较好。 相似文献
9.
We propose a system for the analysis of magnetotelluric (MT) data, which makes use of the invariant characteristics of the impedance tensor such as the maximum and minimum induction curves and the phase tensor. We examine the coefficients of the appearance and normalization of principal values of the impedance tensor. By the case study for Koryakiya, it is shown that the three-dimensional (3D) mathematical modeling and the Wiese-Parkinson vectors allow one to correct the results of one-dimensional (1D) and two-dimensional (2D) inversion of MT curves. Comparison between model and observed data based on the 1D inversion of MTS curves provides a pictorial view of the distortions of MT curves and their sensitivity to the parameters of a geological cross section. 相似文献
10.
George R. Jiracek 《Surveys in Geophysics》1990,11(2-3):163-203
11.
可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性. 相似文献
12.
当前无人机频率域半航空电磁方法(SAEM)成为地球物理勘探中的新兴技术,该方法通过空中无人飞行器测量地面上单个或多个可控源的垂直磁场.本文为无人机频率域SAEM开发了三维反演程序,正演采用交错网格有限差分.由于无人机采集的数据量巨大,因此使用了有限内存拟牛顿法(LBFGS)实现快速的三维反演,以避免计算和存储巨大的灵敏度矩阵.为了提升无人机频率域SAEM的三维反演效果,本文进一步实现了无人机频率域SAEM与大地电磁(MT)的三维联合反演,并对联合反演实施了地震结构约束,结构约束采用交叉梯度法实现.MT可以增加勘探深度,而结构约束可以提高电磁法反演的分辨率.最后通过建立四个理论模型,证明了本文程序的有效性. 相似文献
13.
基于大地电磁阻抗张量分解技术,本文提出了两种电性主轴方位的统计描述图像：随频率变化的统计分布成像（频率分布云图）和随测点序列变化的统计分布成像（测点分布云图）.这两种图像与传统的统计玫瑰图一起,较全面地描述了最佳主轴的分布特征.在进行构造维性分析过程中,通过定义二维有效因子e2d,来压制一维结构和三维结构、突出纯二维结构的影响.e2d被用于电性主轴的统计加权,有效地起到了滤波的作用;同时,统计成像中还考虑了数据质量的影响.为了得到稳定、高质量的区域阻抗张量数据,提出并实现了共主轴的多测点-多频点阻抗张量分解新算法.最终,完成了以上各项处理手段的可视化实现.本文通过两个理论模型和一个实测算例,以共轭阻抗法（CCZ法）为基础,展示了这一新技术的有效性. 相似文献
14.
A magnetotelluric acquisition and modelling campaign, comprising 400 full tensor soundings, was carried out in the Sub-Andean Foothills in the area located south of Santa Cruz de la Sierra, Bolivia. The objective of the survey was to improve the imaging of the Paleozoic section and provide insight to the overall structure of the Devonian Silurian complex. The acquired data were inverted assuming both two-dimensional and three-dimensional dimensionality of the subsurface structure, adopting a multi-step iterative workflow, which began applying unconstrained inversion and introduced increasing constraints determined at each following step based on the interpretation of the latest inversion results. The geological interpretation of the final inversion results allowed to better image the expected complexity in the structures, particularly highlighting several different levels of detachment, mainly in the area of Kirusillas Formation. Furthermore, an oblique sense was identified of the main tendency of N-S deformation, related to a dextral transgressional movement that generates discrete dislocations, at a length scale of the order of a few kilometres. 相似文献
15.
In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated gzz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of gzz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area. 相似文献
16.
Zhi-Yong Zhang Han-Dong Tan Kun-Peng Wang Chang-Hong Lin Bin Zhang Mao-Bi Xie 《应用地球物理》2016,13(1):13-24
Traditional two-dimensional (2D) complex resistivity forward modeling is based on Poisson’s equation but spectral induced polarization (SIP) data are the coproducts of the induced polarization (IP) and the electromagnetic induction (EMI) effects. This is especially true under high frequencies, where the EMI effect can exceed the IP effect. 2D inversion that only considers the IP effect reduces the reliability of the inversion data. In this paper, we derive differential equations using Maxwell’s equations. With the introduction of the Cole–Cole model, we use the finite-element method to conduct 2D SIP forward modeling that considers the EMI and IP effects simultaneously. The data-space Occam method, in which different constraints to the model smoothness and parametric boundaries are introduced, is then used to simultaneously obtain the four parameters of the Cole—Cole model using multi-array electric field data. This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity. To improve the computational efficiency, message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion. Synthetic datasets were tested using both serial and parallel algorithms, and the tests suggest that the proposed parallel algorithm is robust and efficient. 相似文献
17.
This paper presents a short theoretical summary of the series expansion-based 2.5D combined geoelectric weighted inversion (CGWI) method and highlights the advantageous way with which the number of unknowns can be decreased due to the simultaneous characteristic of this inversion. 2.5D CGWI is an approximate inversion method for the determination of 3D structures, which uses the joint 2D forward modeling of dip and strike direction data. In the inversion procedure, the Steiner’s most frequent value method is applied to the automatic separation of dip and strike direction data and outliers. The workflow of inversion and its practical application are presented in the study. For conventional vertical electrical sounding (VES) measurements, this method can determine the parameters of complex structures more accurately than the single inversion method. Field data show that the 2.5D CGWI which was developed can determine the optimal location for drilling an exploratory thermal water prospecting well. The novelty of this research is that the measured VES data in dip and strike direction are jointly inverted by the 2.5D CGWI method. 相似文献
18.
G. Michael Hoversten 《Geophysical Prospecting》1996,44(6):935-961
Hydrocarbon exploration in the Papuan fold belt is made extremely difficult by mountainous terrain, equatorial jungle and thick karstified Miocene limestones at the surface. The high-velocity karstified limestones at or near the surface often render the seismic technique useless for imaging the subsurface. In such areas magnetotellurics (MT) provides a valuable capability for mapping subsurface structure. The main structural interface which can be mapped with MT, due to the large electrical contrast, is the contact between the resistive Darai limestone and the underlying conductive sediments of the Ieru Formation. In some areas the base of the Darai can be mapped with reasonable accuracy by fitting 1D models to the observed MT data. However, in many cases where 2D and 3D effects are severe, 1D interpretations can yield dramatically incorrect results. Numerical and field data examples are presented which demonstrate the severity of the 1D errors and the improvements in accuracy which can be achieved using a 2D inverse solution. Two MT lines over adjacent anticlines, both with well control and seismic data, are used to demonstrate the application of 1D and 2D inversions for structural models. In both cases the seismic data provide no aid in the interpretations. The example over the Hides anticline illustrates a situation where 1D inversion of either TE or TM mode provides essentially the same depth to base of Darai as 2D inversion of both TE and TM. Both models provide base Darai depth estimates which are within 10% of that measured in the well. The example over the Angore anticline illustrates the inadequacy of 1D inversion in structurally complex geology complicated by electrical statics. The TE mode fits a 1D Darai thickness of 800 metres while the TM mode fits a 1D Darai thickness of 3500 metres, bracketing the thickness of 2450 metres observed in the well. The final 2D inversion model provides a depth estimate of 2250 metres. Four MT lines along the Angore anticline have been interpreted using 2D inversion. A high degree of correlation exists between lineaments observed on an airborne radar image and zones of low resistivity within the high-resistivity material interpreted as Darai limestone. These low-resistivity zones are interpreted as fault zones. Three-dimensional modelling has been used to simulate 3D statics in an otherwise 2D earth. These data were used to test the Groom-Bailey (GB) decomposition for possible benefits in reducing static effects and estimating geoelectric strike in the Papua New Guinea (PNG) field data. It has been found that the GB decomposition can provide improved regional 2D strike estimates in 3D contaminated data. However, in situations such as PNG, where the regional 2D strike is well established and hence can be fixed, the GB decomposition provides apparent resistivities identical to those simply rotated to strike. 相似文献
19.
在前人研究基础上,对Groom-Bailey(GB)张量分解畸变因子和区域阻抗的求解方法进行了改进.首先,通过Swift旋转与GB分解的扭变和剪切矩阵的求逆变换,利用变换后区域阻抗主对角元素为0的条件获得关于扭变因子和剪切因子的超定方程组,采用模拟退火全局优化算法进行求解.其次,由得到的扭变因子和剪切因子,结合Swift旋转确定的走向角和区域阻抗元素的估计,作为非线性最小二乘局部优化算法的初始值,对GB分解定义式的超定方程组进行求解,得到各畸变参数和区域阻抗的解.通过模型试验验证了方法的正确,对方法的稳定性进行了比较与评价,并通过与已有结果的对比和实际资料的应用,表明了方法实际应用的效果. 相似文献
20.
Fernando A. Monteiro Santos Antnio Trota Antnio Soares Rafael Luzio Nuno Loureno Liliana Matos Eugnio Almeida Joo L. Gaspar Jorge M. Miranda 《Journal of Applied Geophysics》2006,59(4):314-323
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. 相似文献