首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 358 毫秒
1.
We propose the approach to 3D inversion of airborne electromagnetic data, which is intended for discovering subvertical bodies overlapped by essentially inhomogeneous conductive layers. The approach is based on the geometric inversion in which a geoelectrical medium is parameterized with the use of block structures. During the inversion, the coordinates of the borders between the blocks and the rows of the blocks as well as resistivities inside them are determined. In order to solve the forward problem of the airborne electromagnetic survey, we use the non-conforming optimized mesh with the hexahedral cells, which enables us to reduce the number of degrees of freedom and smoothly approximate the curved borders of a geological medium. For a more reliable discovery of subvertical objects, we propose to carry out 3D inversions at several rotations of block structures relative to the flight lines. The workability of this approach is demonstrated using the data which are synthesized for complex geoelectrical models with topography, inhomogeneous overlapping layers and target subvertical bodies oriented differently relative to the flight lines. The results of this investigation show that, in some way or other, the elongated subvertical object is discovered and its orientation (the direction of its long side) is defined at different rotations of block structures used in 3D inversions. However, the most accurate recovery of the subvertical object length is achieved when the direction of its long side almost coincides with the direction of one of the block structures axes. Thus, the block structures rotations allow not only more reliably discovering a target object in complex geoelectrical conditions, but also more exactly defining its orientation and length.  相似文献   

2.
Many synthetic model studies suggested that the best way to obtain good 3D interpretation results is to distribute the MT sites at a 2D grid array with regular site spacing over the target area. However, MT 3D inversion was very difficult about 10 years ago. A lot of MT data were collected along one profile and then interpreted with 2D inversion. How to apply the state-of-the-art 3D inversion technique to interpret the accumulated mass MT profiles data is an important topic. Some studies on 3D inversion of measured MT profile data suggested that 2D inversions usually had higher resolution for the subsurface than 3D inversions. Meanwhile, they often made their interpretation based on 2D inversion results, and 3D inversion results were only used to evaluate whether the overall resistivity structures were correct. Some researchers thought that 3D inversions could not resolute the local structure well, while 2D inversion results could agree with the surface geologic features much well and interpret the geologic structures easily. But in the present paper, we find that the result of 3D inversion is better than that of 2D inversion in identifying the location of the two local faults, the Shade Fault(SDF)and the Yunongxi Fault(YNXF), and the deep structures. In this paper, we first studied the electrical structure of SDF and YNXF based on a measured magnetotelluric(MT) profile data. Besides, from the point of identifying active faults, we compared the capacity of identifying deep existing faults between 2D inversion models and 3D models with different inversion parameters. The results show that both 2D and 3D inversion of the single-profile data could obtain reasonable and reliable electrical structures on a regional scale. Combining 2D and 3D models, and according to our present data, we find that both SDF and YNXF probably have cut completely the high resistivity layer in the upper crust and extended to the high conductivity layer in the middle crust. In terms of the deep geometry of the faults, at the profile's location, the SDF dips nearly vertically or dips southeast with high dip angle, and the YNXF dips southeast at depth. In addition, according to the results from our measured MT profile, we find that the 3D inversion of single-profile MT data has the capacity of identifying the location and deep geometry of local faults under present computing ability. Finally, this research suggests that appropriate cell size and reasonable smoothing parameters are important factors for the 3D inversion of single-profile MT data, more specifically, too coarse meshes or too large smoothing parameters on horizontal direction of 3D inversion may result in low resolution of 3D inversions that cannot identify the structure of faults. While, for vertical mesh size and data error thresholds, they have limited effect on identifying shallow tectonics as long as their changes are within a reasonable range. 3D inversion results also indicate that, to some extent, adding tippers to the 3D inversion of a MT profile can improve the model's constraint on the deep geometry of the outcropped faults.  相似文献   

3.
Using a subset of the SEG Advanced Modeling Program Phase I controlled‐source electromagnetic data, we apply our standard controlled‐source electromagnetic interpretation workflows to delineate a simulated hydrocarbon reservoir. Experience learned from characterizing such a complicated model offers us an opportunity to refine our workflows to achieve better interpretation quality. The exercise proceeded in a blind test style, where the interpreting geophysicists did not know the true resistivity model until the end of the project. Rather, the interpreters were provided a traditional controlled‐source electromagnetic data package, including electric field measurements, interpreted seismic horizons, and well log data. Based on petrophysical analysis, a background resistivity model was established first. Then, the interpreters started with feasibility studies to establish the recoverability of the prospect and carefully stepped through 1D, 2.5D, and 3D inversions with seismic and well log data integrated at each stage. A high‐resistivity zone is identified with 1D analysis and further characterized with 2.5D inversions. Its lateral distribution is confirmed with a 3D anisotropic inversion. The importance of integrating all available geophysical and petrophysical data to derive more accurate interpretation is demonstrated.  相似文献   

4.
Electromagnetic methods are routinely applied to image the subsurface from shallow to regional structures. Individual electromagnetic methods differ in their sensitivities towards resistive and conductive structures and in their exploration depths. If a good balance between different electromagnetic data can be be found, joint 3D inversion of multiple electromagnetic datasets can result in significantly better resolution of subsurface structures than the individual inversions. We present a weighting algorithm to combine magnetotelluric, controlled source electromagnetic, and geoelectric data. Magnetotelluric data are generally more sensitive to regional conductive structures, whereas controlled source electromagnetic and geoelectric data are better suited to recover more shallow and resistive structures. Our new scheme is based on weighting individual components of the total data gradient after each model update. Norms of individual data residuals are used to assess how much of the total data gradient must be assigned to each method to achieve a balanced contribution of all datasets for the joint inverse model. Synthetic inversion tests demonstrate advantages of joint inversion in general and also the influence of the weighting. In our tests, the controlled source electromagnetic data gradients are larger than those of the magnetotelluric and geoelectric datasets. Consequently, direct joint inversion of controlled source electromagnetic, magnetotelluric, and geoelectric data results in models that are mostly dominated by structures required by the controlled source electromagnetic data. Applying the new adaptive weighting scheme results in an inversion model that fits the data better and resembles more the original model. We used the modular system electromagnetic as a framework to implement the new joint inversion and briefly describe the new modules for forward modelling and their interfaces to the modular system electromagnetic package.  相似文献   

5.
地面可控源频率测深三维非线性共轭梯度反演   总被引:8,自引:8,他引:0       下载免费PDF全文
讨论了地面可控源电磁勘探三维非线性共轭梯度反演的可行性以及反演过程中考虑场源的必要性.反演采用非线性共轭梯度反演方法.反演过程中,模型响应利用交错网格有限差分技术计算.反演数据采用与发射源平行的电场x分量Ex.利用层状导电模型作为背景,设计了两个理论模型进行数值试验:第一个模型中包含两个电阻率异常,以检验反演的有效性;第二个模型中,在测区外设置了一个低阻异常,以考察源的信息在反演中的作用.两个模型的反演分别从层状背景模型开始,迭代120次后终止.数值试验结果表明,(1)非线性共轭梯度反演所获得的电阻率分布和理论模型吻合较好;(2)非线性共轭梯度算法收敛速度较慢,需要较多的迭代次数完成反演;(3)对于可控源频率电磁勘探,必须考虑源位置信息.因此,本文采用考虑场源信息的地面可控源非线性共轭梯度反演方法能完成真正意义上的可控源频率电磁测深数据的反演.  相似文献   

6.
海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性.  相似文献   

7.
The existence of aquifers extending from land beneath the sea floor up to a distance of several kilometres has been observed and examined all over the world. The coastal aquifer of Israel is a heavily used groundwater reservoir which has to be constantly monitored to ensure the drinking water supply. Former land-based electromagnetic measurements show that it is, in several places, blocked to seawater intrusion and is consequently a candidate for submarine extension. Multicomponent long-offset transient electromagnetic measurements were carried out offshore on the coast of Israel. We deployed a 400-m-long grounded dipole as transmitter and several electric and magnetic receivers on the sea floor up to a distance of 4.8 km from the coast. Altogether, we deployed 8 transmitter positions and received data sets at 14 receiver stations onshore and offshore, with offsets of mostly 400 and 800 m. In this paper, we present the survey and 1D Occam and Marquardt inversions of the offshore horizontal electric components in the broadside and inline configuration. In addition, the vertical magnetic component in the broadside position is also considered. Only single inversions, both single offset and single component, were used to detect the aquifer under sea sediments. We prove the submarine existence of the Israeli coastal aquifer up to a distance to the coast of approximately 3.2 to 3.6 km using all measured long-offset transient electromagnetic components. In addition, we present modelling studies with synthetic data derived from a subsurface model adjusted to our survey area with very shallow water from 10 to 50 m. As part of the planning before the survey, a parameter study of the expected subsurface, the examination of the airwave phenomenon and the justification for our 1D inversion strategy are shown. More detailed eigen parameter analyses are added to explain the measured data.  相似文献   

8.
为探讨二维反演方法在三维电性结构中的适应性问题,本研究中设计了一系列的二维/三维正演模型进行实验计算,分析了三维高阻/低阻异常体对模型响应的畸变作用,并从反演模式选择和数据旋转方向两个方面进行模型二维反演的对比分析,与三维反演的结果进行了比较,最后采用了实测数据进一步进行了二维和三维反演的比较实验.实验和研究结果表明,在剖面选择方面,在剖面方向与垂直主构造方向相差不大的情况下,截取剖面方向,将电性主轴旋转到垂直剖面方向的二维反演结果与垂直主构造方向的反演结果都可以较好地还原正演模型,在大的构造的反映上并无太大差异.在地下为二维或近三维条件时,正演模型的主要结构都可以较好地被二维和三维反演解析出来.二维的反演结果可能甚至会比三维的反演结果的边界更清晰,更精确.然而,对于具有较强的三维结构的模型而言,其二维反演结果与原始模型可能仍然存在较大差异,其中TM+TP或TM模式的二维反演结果相对更接近原始模型,而TE模式的结果往往会有较大误差,需要在解释时特别注意以免得出错误结论.  相似文献   

9.
Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.  相似文献   

10.
On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained.  相似文献   

11.
The detecting capabilities of some electrical arrays for the estimation of position, size and depth of small-scale targets were examined in view of the results obtained from 2D inversions of apparent-resistivity data. The two-sided three-electrode apparent-resistivity data are obtained by the application of left- and right-hand pole–dipole arrays that also permit the computation of four-electrode and dipole–dipole apparent-resistivity values without actually measuring them. Synthetic apparent-resistivity data sets of the dipole–dipole, four-electrode and two-sided three-electrode arrays are calculated for models that simulate buried tombs. The results of two-dimensional inversions are compared with regard to the resolution in detecting the exact location, size and depth of the target, showing some advantage for the two-sided three-electrode array. A field application was carried out in the archaeological site known as Alaca Hoyuk, a religious temple area of the Hittite period. The two-dimensional inversion of the two-sided three-electrode apparent-resistivity data has led to locating a part of the city wall and a buried small room. The validity of the interpretation has been checked against the results of subsequent archaeological excavations.  相似文献   

12.
On the basis of the dispersion relation of magnetotelluric response functions (MTRF), a filter coefficient algorithm has been made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivity data. The tests of theoretical models and observed magnetotelluric (MT) data show that this algorithm is effective. Comparing the impedance phase estimated using dispersion relation with the observed phase, it can be checked whether the dispersion relation between the observed apparent resistivities and phase data was satisfied. The use of phase data corrected using the dispersion relation in the joint inversion for MT impedance is advantageous to obtain more reliable inversion results. The problems on the one-dimensional joint inversion for the (MT) apparent resistivity and the apparent resistivity of the frequency electromagnetic sounding (FEMS) with horizontal electric dipole, whose observed frequency bands are linked up each other, are studied. The observed data of two kinds of electromagnetic (EM) methods at two sites are used to inverse, the comparison with the drilling data show the results are more reliable. To supply the phase data of FEMS using the dispersion relation, for the apparent resistivity-phase data and impedance real part-imaginary part apparent resistivities of two kinds of EM methods the imitated MT joint inversions are made, and more similar results also are obtained. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 91–96, 1993. The projects sponsored by the Chinese Joint Seismological Science Foundation.  相似文献   

13.
One of the ways to improve the information content of a set of field data is that of combining the interpretation of disparate data sets. Electromagnetic and direct current resistivity methods suffer from inherent equivalence problem. Application of joint inversion for these measurements can overcome the problem of equivalence very well. In the present work, synthetic data from vertical electrical sounding (VES) and horizontal coplanar low-frequency induction sounding (EMHD) are inverted individually and jointly over different types of 1D earth structures. Global optimization with Monte Carlo Multistart algorithm was used in the calculations. The results obtained from the inversions of synthetic data indicate that the joint inversion significantly improves the solution reducing the ambiguity of the models.  相似文献   

14.
井中垂直双极源体积分方程法三维模拟研究   总被引:2,自引:0,他引:2       下载免费PDF全文
类比水平方向电偶极源和垂直方向磁偶极源电磁场求解的原理,又利用偶极沿电流方向积分是双极原理,推导出地下均匀半空间中垂直双极源电场和磁场解析计算公式.采用体积分方程法对井中垂直双极源电磁响应进行了三维模拟.它对提高井中垂直双极源圈定储层范围的精度有着非常重要的意义,同时也为井中垂直双极源三维反演问题奠定了基础.  相似文献   

15.
为解决地球物理反演中多解性的问题,综合多种地球物理信息的联合反演受到了广泛的关注.本文依据不同地球物理响应可能由相同异常体引起,而不同地球物理分布参数之间存在相关性等特点,提出了一种基于局部Pearson相关系数约束的联合反演方法.该方法假设每个局部区域模型参数的分布具有线性相关特性,在拟合不同类型观测数据时,对局部模型参数施加相关性约束,进行联合反演以减少多解性.本文采用交替迭代联合反演流程,改善了同一目标函数下联合反演收敛性和速度问题.基于新的联合反演方法和流程,我们测试了三维大地电磁和重力仿真数据的联合反演.结果表明,本文提出的基于局部相关性约束的联合反演方法,能充分利用大地电磁和重力观测数据信息,有效改善单一地球物理反演收敛性和多解性的问题,反演效果得到明显提升.  相似文献   

16.
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.  相似文献   

17.
The time-domain controlled source electromagnetic method is a geophysical prospecting tool applied to image the subsurface resistivity distribution on land and in the marine environment. In its most general set-up, a square-wave current is fed into a grounded horizontal electric dipole, and several electric and magnetic field receivers at defined offsets to the imposed current measure the electromagnetic response of the Earth. In the marine environment, the application often uses only inline electric field receivers that, for a 50% duty-cycle current waveform, include both step-on and step-off signals. Here, forward and inverse 1D modelling is used to demonstrate limited sensitivity towards shallow resistive layers in the step-off electric field when transmitter and receivers are surrounded by conductive seawater. This observation is explained by a masking effect of the direct current signal that flows through the seawater and primarily affects step-off data. During a step-off measurement, this direct current is orders of magnitude larger than the inductive response at early and intermediate times, limiting the step-off sensitivity towards shallow resistive layers in the seafloor. Step-on data measure the resistive layer at times preceding the arrival of the direct current signal leading to higher sensitivity compared to step-off data. Such dichotomous behaviour between step-on and step-off data is less obvious in onshore experiments due to the lack of a strong overlying conductive zone and corresponding masking effect from direct current flow. Supported by synthetic 1D inversion studies, we conclude that time-domain controlled source electromagnetic measurements on land should apply both step-on and step-off data in a combined inversion approach to maximize signal-to-noise ratios and utilize the sensitivity characteristics of each signal. In an isotropic marine environment, step-off electric fields have inferior sensitivity towards shallow resistive layers compared to step-on data, resulting in an increase of non-uniqueness when interpreting step-off data in a single or combined inversion.  相似文献   

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

19.
Large‐scale inversion methods have been recently developed and permitted now to considerably reduce the computation time and memory needed for inversions of models with a large amount of parameters and data. In this work, we have applied a deterministic geostatistical inversion algorithm to a hydraulic tomography investigation conducted in an experimental field site situated within an alluvial aquifer in Southern France. This application aims to achieve a 2‐D large‐scale modeling of the spatial transmissivity distribution of the site. The inversion algorithm uses a quasi‐Newton iterative process based on a Bayesian approach. We compared the results obtained by using three different methodologies for sensitivity analysis: an adjoint‐state method, a finite‐difference method, and a principal component geostatistical approach (PCGA). The PCGA is a large‐scale adapted method which was developed for inversions with a large number of parameters by using an approximation of the covariance matrix, and by avoiding the calculation of the full Jacobian sensitivity matrix. We reconstructed high‐resolution transmissivity fields (composed of up to 25,600 cells) which generated good correlations between the measured and computed hydraulic heads. In particular, we show that, by combining the PCGA inversion method and the hydraulic tomography method, we are able to substantially reduce the computation time of the inversions, while still producing high‐quality inversion results as those obtained from the other sensitivity analysis methodologies.  相似文献   

20.
三维频率域可控源电磁反演研究   总被引:1,自引:1,他引:0       下载免费PDF全文
对于三维可控源电磁,反演计算效率、张量测量、旁侧效应以及阴影效应是目前研究的热点.本文正演采用基于库仑规范条件的耦合势有限体积算法,反演采用有限内存BFGS(L-BFGS)算法.合成数据反演结果表明:(1)有限内存BFGS法比非线性共轭梯度法,在反演计算效率上具有一定的优势,更适合求解大规模三维可控源电磁反演问题.(2)张量可控源电磁法相对于标量可控源电磁法,前者在模型分辨率上优于后者.(3)在某个区域无法布置测网的情况下,我们可利用旁侧效应在异常体周围布置测网进行三维反演,从而获得真实异常体的信息.同时,为避免阴影效应,我们应在测网外增加可控源电磁控制点,使得三维反演的数据更加完备.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号