共查询到19条相似文献，搜索用时 953 毫秒
1.
为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释. 相似文献
2.
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征. 相似文献
3.
CSAMT在山区金属矿勘查中,采用各种滤波和相位积分之类的处理方法,校正因地形起伏和局部电性不均匀引起的静态效应,往往难保奏效,开发消除静态效应的新方法是提高CSAMT资料处理与解释水平和方法应用效果的重要研究课题. 本文以如何消除地形影响为重点,对起伏地形下CSAMT二维大地三维源地电模型,采用加权余弦数值积分法,进行波数域电磁场二维有限单元法正演. 为模拟复杂地形地电模型,选取交叉对称网格三角形剖分法,实现了在国内常用赤道电偶极装置的CSAMT二维正演计算;在二维正演的基础上,开发了基于奥克姆反演法的CSAMT二维反演技术,研制出一套起伏地形下CSAMT二维正反演处理与解释方法技术系统. 通过理论模型试算和实测数据处理证实,本系统能有效地削减起伏地形影响. 在找矿应用中,该系统反演的电阻率断面,极大地消除了起伏地形影响和静态效应,突显出清晰的控矿构造和矿体的异常,取得了重要成效. 相似文献
4.
5.
6.
在山区进行MT勘探时,用规则网格有限元方法模拟起伏地形会受到限制.本文采用非结构三角网格可以有效地模拟任意二维地质结构,如起伏地形、倾斜岩层和多尺度构造等.正演引入自适应有限元方法,其在网格剖分过程中能根据单元误差自动细化网格,保证了正演结果的精度.将自适应有限元与Occam算法结合,且引用并行处理技术提高正反演计算速度.通过对比两个理论模型,讨论了地形对MT正演响应的影响;其次进行了不同地电模型带地形反演展示了本文算法的正确性和适用性;最后将该方法应用于实测MT数据处理,证明了自适应非结构有限元方法是复杂地形下处理MT数据的有力工具. 相似文献
7.
本文实现了非平坦地形条件下电阻率三维反演,讨论了几种消除反演中地形影响的方法. 结果表明,只有将地形直接带入反演算法中,进行带地形电阻率三维反演才能有效消除地形影响及其对反演结果的偏差,得到与地下电性结构相符的反演结果. 相似文献
8.
航空瞬变电磁法以其速度快、成本低、通行性好等的优势能够有效的应用于地质地形条件复杂的地区.目前对于航空瞬变电磁法的研究主要基于平坦地形的理想情况,对于地形效应的研究相对较少,然而实际应用中地形不可避免,若忽略地形影响将对资料解释造成较大的误差,从而制约航空电磁方法的进一步发展.本文基于交错网格的时域有限差分方法对三维起伏地形条件下航空瞬变电磁进行正演模拟,在保证算法准确性的前提下给出大量模型算例.以经典地形模型为例,利用所给方法计算三维正演响应,结果显示起伏对于航空瞬变电磁数据有着显著的影响且影响主要集中在早期.而后,以实际地质资料为基础,构建起伏地形条件下包含多个异常体的三维复杂模型,计算了复杂模型的航空瞬变电磁响应,并给出三维全域视电阻率曲线,从而对地形效应的影响有了更加直观的认知.最后,通过大量模型讨论了地形的尺寸参数、电性参数、飞行轨迹与飞行高度等因素变化对于航空瞬变电磁数据的影响情况,并得出有价值的结论. 相似文献
9.
海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性. 相似文献
10.
11.
Myeong-Jong Yi Jung-Ho Kim Yoonho Song Seong-Jun Cho Seung-Hwan Chung Jung-Hee Suh 《Geophysical Prospecting》2001,49(4):483-497
We have developed a three-dimensional inverse scheme for carrying out DC resistivity surveys, incorporating complicated topography as well as arbitrary electrode arrays. The algorithm is based on the finite-element approximation to the forward problem, so that the effect of topographic variation on the resistivity data is effectively evaluated and incorporated in the inversion. Furthermore, we have enhanced the resolving power of the inversion using the active constraint balancing method. Numerical verifications show that a correct earth image can be derived even when complicated topographic variation exists. By inverting the real field data acquired at a site for an underground sewage disposal plant, we obtained a reasonable image of the subsurface structures, which correlates well with the surface geology and drill log data. 相似文献
12.
Ahmad Neyamadpour 《Studia Geophysica et Geodaetica》2011,55(4):755-769
The applicability of three kinds of electrode configurations used to delineate a buried horizontal pipe was studied. A 3D
resistivity imaging survey was carried out along eight parallel lines using pole-pole, pole-dipole, and dipole-dipole arrays
with 1m minimum electrode spacings. Roll-along measurements were carried out to cover a rectangular grid. The 2D and 3D least
squares algorithms based on the robust inversion method were used in the inversion of the apparent resistivity data sets.
The 2D inversion of data sets could not delineate the orientation and dimension of the subsurface anomalies clearly. To obtain
more accurate results, a 3D joint inversion of the pole-pole and pole-dipole data sets was performed, as well as of pole-pole
and dipole-dipole data sets. In this case, both horizontal and vertical dimensions of subsurface structures were resolved.
The resulting model obtained from each array was compared to those of joint inversion method. The result showed that the horizontal
resolution does not improve so much as that in the vertical direction when joint inversion is applied. 相似文献
13.
适用于近地表(2000m以内)勘探的频率域电磁法主要有音频大地电磁法(audio-frequency magnetotellurics,AMT),无线电大地电磁法(radio-magnetotellurics,RMT),可控源音频大地电磁法(controlled source audiofrequency magnetotellurics,CSAMT),广域电磁法(Wide Field Electromagnetic Method,WFEM).本文拟从最新的数据采集技术、数据处理技术、正反演算法、实例等四个方面,论述适用于浅部勘探的AMT,RMT,CSAMT和WFEM方法的国内外最新进展,总结目前AMT,RMT,CSAMT和WFEM方法遇到的困难,对潜在的发展方向提出建议.综述表明:(1)张量测量、多站阵列、多站叠加可提高AMT、RMT和CSAMT数据的质量.利用近区数据WFEM法可获得良好的效果.国产与国外仪器在质量方面的差距正在逐步缩小.(2)数学形态滤波技术、Hilbert-Huang变换等可有效分离出有用的数据,局部畸变仍然是亟待解决的难题,需要更为深入的研究.(3)矢量有限元与非结构网格的出现大幅度提高了有限元处理复杂电磁问题模拟的精度与应用范围,成为目前电磁正演的首选工具.完全非线性反演算法仍然局限于1D、2D问题,共轭梯度法和高斯牛顿算法等为解决3D问题的发展趋势.地质约束的引入和多数据联合反演可以减小反演的非唯一性.各向异性的反演为目前反演研究的热点之一.(4)野外数据解释的正确性严重依赖于对地下结构先期的维性判别,在2D特性不明显、3D特性明显时,需要采用3D进行反演解释. 相似文献
14.
Irfan Akca 《Acta Geophysica》2016,64(2):443-462
ELRIS2D is an open source code written in MATLAB for the two-dimensional inversion of direct current resistivity (DCR) and time domain induced polarization (IP) data. The user interface of the program is designed for functionality and ease of use. All available settings of the program can be reached from the main window. The subsurface is discre-tized using a hybrid mesh generated by the combination of structured and unstructured meshes, which reduces the computational cost of the whole inversion procedure. The inversion routine is based on the smoothness constrained least squares method. In order to verify the program, responses of two test models and field data sets were inverted. The models inverted from the synthetic data sets are consistent with the original test models in both DC resistivity and IP cases. A field data set acquired in an archaeological site is also used for the verification of outcomes of the program in comparison with the excavation results. 相似文献
15.
M. E. Sholpo 《Izvestiya Physics of the Solid Earth》2013,49(1):113-119
The sensitivity of transversal apparent resistivity to the changes in the electrical resistivity of elements of a two-dimensional (2D) geoelectrical structure is studied by numerical modeling. This sensitivity is found to have a series of specific features, due to which the monitoring of transversal apparent resistivity can serve as a helpful addition to the monitoring of longitudinal apparent resistivity to trace the dynamics of subsurface and shallow crustal elements of the medium. It is shown that the method previously suggested for the inversion of relative changes in apparent resistivity into relative changes in electrical resistivities of the elements of the geoelectrical structure is applicable to the transversal electrical resistivity. 相似文献
16.
The travel time and amplitude of ground-penetrating radar (GPR) waves are closely related to medium parameters such as water content, porosity, and dielectric permittivity. However, conventional estimation methods, which are mostly based on wave velocity, are not suitable for real complex media because of limited resolution. Impedance inversion uses the reflection coefficient of radar waves to directly calculate GPR impedance and other parameters of subsurface media. We construct a 3D multiscale stochastic medium model and use the mixed Gaussian and exponential autocorrelation function to describe the distribution of parameters in real subsurface media. We introduce an elliptical Gaussian function to describe local random anomalies. The tapering function is also introduced to reduce calculation errors caused by the numerical simulation of discrete grids. We derive the impedance inversion workflow and test the calculation precision in complex media. Finally, we use impedance inversion to process GPR field data in a polluted site in Mongolia. The inversion results were constrained using borehole data and validated by resistivity data. 相似文献
17.
Radio-magnetotelluric(RMT)是以无线电发射机为信号源的一种地球物理勘探方法,近年来被广泛应用于数米至数十米内的近地表工程和环境地球物理勘探.目前,各类电磁资料的反演均是以寻求满足目标拟合差的地下介质电阻率分布为目的.然而,对于勘探频率为10~300kHz的RMT数据,由介电常数所引起的波动场在总场中的比例可达20%以上,在这种情况下,忽略介电常数,仅通过电阻率参数的反演来进行数据拟合势必降低反演资料解释的准确性.为解决这一问题,本文研究了基于电阻率-介电常数的双参数同步反演算法.构建了一个全新的双参数目标函数,并推导了双参数反演迭代方程组;通过灵敏度分析,研究了电阻率和介电常数对正演响应的影响,并据此提出相对电导率的概念,统一了反演参数的灵敏度;通过理论模型分析了参考频率、双参数正则化因子对反演结果的影响,并给出了一般性的参数优选方案.此外,为了能够灵活处理复杂地形,本文采用非结构的正反演双网格进行模型离散,并通过局部加密技术保证反演的速度和精度.最后,对一带地形的理论模型分别进行了单参数和双参数反演,结果表明单参数反演无法正确反映出地电信息,而双参数反演能够准确得到异常的分布,验证了本文所开发的双参数反演程序的有效性. 相似文献
18.
Vikas C. Baranwal Antje Franke Ralph-Uwe Börner Klaus Spitzer 《Journal of Applied Geophysics》2011,75(2):363-372
We present a 2-D inversion code incorporating a damped least-squares and a minimum-model approach for plane wave electromagnetic (EM) methods using an adaptive unstructured grid finite element forward operator. Unstructured triangular grids permit efficient discretization of arbitrary 2-D model geometries and, hence, allow for modeling arbitrary topography. The inversion model is parameterized on a coarse parameter grid which constitutes a subset of the forward modeling grid. The mapping from parameter to forward modeling grid is obtained by adaptive mesh refinement. Sensitivities are determined by solving a modified sensitivity equation system arising from the derivative of the finite element equations with respect to the model parameters. Firstly, we demonstrate that surface topography may induce significant effects on the EM response and in the inversion result, and that it cannot be ignored when the scale length of topographic variations is in the order of magnitude of the skin depth. Secondly, the dependency of the inversion on the starting model is discussed for VLF and VLF-R data. Thirdly, we demonstrate the inversion of a synthetic data set obtained from a model with topography. Finally, the inversion approach is applied to field data collected in a region with undulating topography. 相似文献
19.
Ahmad Neyamadpour W. A. T. Wan Abdullah Samsudin Taib 《Studia Geophysica et Geodaetica》2010,54(2):299-311
The objective of this paper is to investigate the applicability of four-electrode arrays in 3D electrical resistivity imaging
survey. A 3D resistivity imaging survey was carried out along fourteen parallel lines using dipole-dipole, Wenner-Schlumberger,
and Wenner arrays with 2 m minimum electrode spacings. Roll-along measurements using a line spacing of 1 m were carried out
covering a grid of 20 × 14 electrodes. The 3D least squares algorithm, based on the robust inversion method, was used in the
inversion of the 3D apparent resistivity data sets. The results show that the 3D electrical resistivity imaging survey using
the Wenner-Schlumberger and the dipole-dipole arrays, or the Wenner and the dipole-dipole arrays, in combination with an appropriate
3D inversion method, can be highly useful when the site conditions do not allow using the pole-pole or pole-dipole arrays. 相似文献