可控源音频大地电磁三维共轭梯度反演研究

可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性.     

地面矩形大定源电磁法频率域三维非线性共轭梯度反演

本文实现了地面矩形大定源三维频率域反演.矩形大定源三维模型响应计算采用交错网格有限差分技术.正演的微分方程为异常电场满足的非齐次Helmholtz方程,方程右手边源项中的大定源产生的背景格林函数由虚界面法结合虚框法计算.频率域三维反演采用非线性共轭梯度反演技术.反演的数据类型为垂直磁场的频率域响应Hz的实部和虚部分量.数值结果表明,(1)三维模型正演模拟数值结果与前人一致,为三维反演奠定基础;(2)针对两个三维导电模型,分别进行了三维反演数值试算.反演结果可以清晰恢复出异常体的电阻率和位置信息,表明地面矩形大定源三维频率域非线性共轭梯度反演具有可行性.本文研究的意义在于,在电磁响应时频转换技术的基础上,如果将野外实测的瞬变电磁数据变换为对应的频率响应,则结合本文提出的三维反演技术,可以为矩形大定源瞬变电磁数据的三维解释提供一个新的思路.     

数据类型对三维地面可控源电磁勘探效果的影响

本文基于考虑人工场源辐射的可控源地面电磁勘探三维数值模拟技术,以目前常用的标量可控源音频大地电磁勘探测量方式为依据,讨论地面可控源三维电磁勘探中,为达到最佳的探测效果,如何选择合适的反演数据类型.地面可控源数据的三维反演采用有限内存拟牛顿方法.反演过程中,三维可控源频率域响应数值计算采用交错网格有限差分法,求解基于二次电场的Helmholtz方程.发射装置采用长度为1000m的有限长直导线源,测量频率为10Hz;测点个数200个,分布在10条剖面上.在异常体分布区分别观测(1)电场Ex分量的振幅和相位,(2)Ex振幅,(3)Ex与Ey分量的振幅和相位,(4)Ex与Hy的振幅和相位.反演的数据类型分别为上述4种数据以及导出的阻抗Zxy振幅和相位.反演模型由30×30×20个网格组成,测区内水平方向大小为50m×50m,垂直方向厚度为50m,最后5层厚度倍增.反演都从均匀半空间开始,迭代120次结束.数值模拟结果发现,(1)单个电场分量Ex,其相位信息对异常体信息提取非常重要,若只反演该电场振幅,深部电阻率分辨率低;(2)观测正交电场的效果比单个电场分量效果好,其浅部异常的边缘效应明显减弱,深部异常形态完整;(3)同时观测正交的电场和磁场,反演效果与只观测电场振幅和相位的相当;(4)从阻抗数据反演得到的异常位置和电阻率分布均有较大的改善,反演模型分辨率最好.因此,在理论上,建议在野外数据采集时,最好测量正交电磁场;次之,同时测量正交的电场;最次的,可以只观测电场Ex的振幅和相位,亦可取得较足够的信息.但如果只能获得电场的振幅信息,反演结果的深部将出现较大的不确定性.     

时间域航空电磁2.5维非线性共轭梯度反演

对于时间域航空电磁法二维和三维反演来说,最大的困难在于有效的算法和大的计算量需求.本文利用非线性共轭梯度法实现了时间域航空电磁法2.5维反演方法,着重解决了迭代反演过程中灵敏度矩阵计算、最佳迭代步长计算、初始模型选取等问题.在正演计算中,我们采用有限元法求解拉式傅氏域中的电磁场偏微分方程,再通过逆拉氏和逆傅氏变换高精度数值算法得到时间域电磁响应.在灵敏度矩阵计算中,采用了基于拉式傅氏双变换的伴随方程法,时间消耗只需计算两次正演,从而节约了大量计算时间.对于最佳步长计算,二次插值向后追踪法能够保证反演迭代的稳定性.设计两个理论模型,检验反演算法的有效性,并讨论了选择不同初始模型对反演结果的影响.模型算例表明:非线性共轭梯度方法应用于时间域航空电磁2.5维反演中稳定可靠,反演结果能够有效地反映地下真实电性结构.当选择的初始模型电阻率值与真实背景电阻率值接近时,能得到较好的反演结果,当初始模型电阻率远大于或远小于真实背景电阻率值时反演效果就会变差.     

基于印模法的地面磁电阻率数据三维非线性共轭梯度反演

本文提出了能提高异常体分辨能力,同时得到绝对电导率的地面磁电阻率数据三维反演方法.磁电阻率响应用准直流的低频磁场代替;数值模拟由频率域电场满足的Helmholtz方程出发,采用三维交错网格有限差分法;长直导线源作为发射源,其中源的计算包含在背景场中;结合地面磁电阻率数据各分量的特点,选择y分量进行反演研究;反演采用三维非线性共轭梯度反演技术,为了提高异常体的深度分辨能力,进行迭代重构反演;用印模法对初始模型进行重构,采用的是辅模型在浅部,元模型在深部的组合方式.从合成数据和实际数据的反演结果可以得到以下的认识:(1)由频率域麦克斯韦方程组出发,低频磁场数据反演可以直接得到电导率,而不是相对电导率之比;(2)采用印模法组合初始模型,进行迭代重构反演,可以提高地面磁电阻率数据反演对异常体的分辨能力,确定埋深位置,同时不会丧失对于浅部异常体的分辨能力;(3)在结合印模法的地面磁电阻率数据三维反演中,深部异常体的分辨能力受地表不均匀导电体影响较小;(4)确定印模深度可以采用上一次重构反演结束时的模型变化量,通过相邻两次重构反演结束时的模型变化量之差来确定迭代重构是否终止.因为静磁场与重力场在数学上的相似性,本文的反演方法可以被运用到重力场等位场的地面数据的反演中.     

基于非结构网格剖分的海洋可控源电磁三维正则化反演研究

为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释.     

2.5维海洋可控源电磁反演算法及影响参数研究

海洋可控源电磁法作为一种较新的地球物理勘探方法,在国外已被成功应用于海洋油气资源与天然气水合物的探测中,而我国在这方面的研究开始较晚.本文编写了基于非线性共轭梯度算法的海洋可控源电磁2.5维反演程序,结合理论数据对反演代码的正确性进行验证并讨论了反演参数对反演运算速度及效果的影响.研究结果表明,1)本文编写的2.5维反演程序正确可靠,计算结果与理论模型一致,2)在非线性共轭梯度反演计算中,正则化因子、预条件矩阵、线性搜索及共轭梯度更新因子等参数对实际的反演速度及精度都存在一定影响,可根据不同的勘探需求,调整反演参数以达到较好的反演结果.     

基于VTI各向异性介质的频率域海洋可控源电磁三维约束反演

近年来,海洋可控源电磁法(MCSEM)被引入油气勘探领域以降低勘探风险.在海洋环境中,受沉积因素所造成的电阻率各向异性的影响,地电模型往往会非常复杂.为更好地反映地下电性结构,本文实现了基于VTI各向异性介质的频率域海洋可控源电磁三维反演.其中,正演采用基于Yee氏交错网格的三维有限差分算法,所形成的离散线性系统通过大规模并行矩阵直接求解器(MUMPS)进行求解.反演采用基于不等式约束的有限内存BFGS(L-BFGS)算法.最后,利用VTI各向异性介质合成数据,分别进行了电阻率各向异性覆盖层和电阻率各向异性高阻层的三维反演,结果表明:(1)基于并行直接法的MCSEM非常适用于海洋电磁所特有的多场源问题;(2)针对各向异性覆盖层模型进行三维各向异性约束反演,提高了解的可靠性;(3)针对电阻率各向异性高阻层,Inline和broadside数据覆盖的反演结果对异常体位置有很好的反映.     

三维频率域可控源电磁反演研究

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

基于全区视电阻率的线源FCSEM二维正演模拟

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.     

基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟

从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.     

3D LBFGS inversion of controlled source extremely low frequency electromagnetic data

The controlled source extremely low frequency (CSELF) electromagnetic method is characterized by extremely long and powerful sources and a huge measurement range. Its electromagnetic field can therefore be affected by the ionosphere and displacement current. Research on 3D forward modeling and inversion of CSELF electromagnetic data is currently in its infancy. This paper makes exploratory attempts to firstly calculate the 1D extremely low frequency electromagnetic field under ionosphere-air-earth coupling circumstances, and secondly analyze the propagation characteristics of the background electromagnetic field. The 3D staggered-grid finite difference scheme for solving for the secondary electric field is adopted and incorporated with the 1D modeling algorithm to complete 3D forward modeling. Considering that surveys can be carried out in the near field and transition zone for lower frequencies, the 3D Limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) inversion of CSELF electromagnetic data is presented (in which the sources, or primary fields, are included), with the aim of directly inverting the impedance data, regardless of where it is acquired. Derivation of the objective functional gradient is the core component in the inversion. Synthetic tests indicate that the well-chosen approximation to the Hessian can significantly speed up the inversion. The model responses corresponding to the coexistence of conductive and resistive blocks show that the off-diagonal components of tensor impedance are much more sensitive to the resistivity variation than the diagonal components. In comparison with conventional scalar inversion, tensor inversion is superior in the recoveries of electric anomalies and background resistivity.     

频率域2维电磁测深资料的RRI反演

本文将大地电磁反演中的RRI方法应用于线源频率测深模型资料的反演中.当背景电导率变化很小时,可用变化前的电场来近似代替变化后的电场,从而在反演方程的推导过程中将有源电磁场中的源项消掉,得到和大地电磁场相同的反演方程,使有源电磁波的反演也可以应用RRI方法.反演过程中所需要的模型资料通过有限元方法得到,该资料不需做近场校正,直接用适合于有源电磁场的RRI方法反演,避免了近场校正带来的误差.数值模型结果证明该方法是可行的.最后,用RRI方法讨论了当源和目标区间存在低阻异常体时只对目标区反演的可行性,对实际工作的解释有一定的指导意义.     

Field test of sub‐basalt hydrocarbon exploration with marine controlled source electromagnetic and magnetotelluric data

The recent use of marine electromagnetic technology for exploration geophysics has primarily focused on applying the controlled source electromagnetic method for hydrocarbon mapping. However, this technology also has potential for structural mapping applications, particularly when the relative higher frequency controlled source electromagnetic data are combined with the lower frequencies of naturally occurring magnetotelluric data. This paper reports on an extensive test using data from 84 marine controlled source electromagnetic and magnetotelluric stations for imaging volcanic sections and underlying sediments on a 128‐km‐long profile. The profile extends across the trough between the Faroe and Shetland Islands in the North Sea. Here, we focus on how 2.5D inversion can best recover the volcanic and sedimentary sections. A synthetic test carried out with 3D anisotropic model responses shows that vertically transverse isotropy 2.5D inversion using controlled source electromagnetic and magnetotelluric data provides the most accurate prediction of the resistivity in both volcanic and sedimentary sections. We find the 2.5D inversion works well despite moderate 3D structure in the synthetic model. Triaxial inversion using the combination of controlled source electromagnetic and magnetotelluric data provided a constant resistivity contour that most closely matched the true base of the volcanic flows. For the field survey data, triaxial inversion of controlled source electromagnetic and magnetotelluric data provides the best overall tie to well logs with vertically transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data a close second. Vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data provided the best interpreted base of the volcanic horizon when compared with our best seismic interpretation. The structural boundaries estimated by the 20‐Ω·m contour of the vertical resistivity obtained by vertical transverse isotropy inversion of controlled source electromagnetic and magnetotelluric data gives a maximum geometric location error of 11% with a mean error of 1.2% compared with the interpreted base of the volcanic horizon. Both the model study and field data interpretation indicate that marine electromagnetic technology has the potential to discriminate between low‐resistivity prospective siliciclastic sediments and higher resistivity non‐prospective volcaniclastic sediments beneath the volcanic section.     

郯庐断裂带中段(沂沭断裂带)电性结构研究与孕震环境

穿过郯庐断裂带中段（沂沭断裂带,36°N）所做的大地电磁测深（MT）剖面长约150km.使用Robust技术和远参考道大地电磁方法处理观测数据.通过分析视电阻率、阻抗相位、Swift二维偏离度和区域走向,定性确定测区的电性结构.二维反演解释中选择非线性共轭梯度（NLCG）方法,使用TE、TM两种模式资料联合反演,沿剖面的二维电性结构显示：自西向东,鲁西隆起、郯庐断裂带、胶莱坳陷及鲁东隆起4个电性区块分别对应,鲁东和鲁西隆起区为高阻,郯庐断裂带电性结构复杂,高、低阻相间,胶莱坳陷为低阻（高导）区.沿MT剖面附近曾发生3个地震,其震源区处在电性变化剧烈部位,并在震源区附近存在高导体.     

Multi‐dimensional analyses of the SEAM controlled source electromagnetic data—the story of a blind test of interpretation workflows

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.     

3D modelling of near-surface, environmental effects on AEM data

This study considers the three-dimensional (3D) modelling of compact, at-surface conductive bodies on frequency domain airborne electromagnetic (AEM) survey data. The context is the use of AEM data for environmental and land quality applications. The 3D structures encountered are typically conductive, of limited thickness (<20 m) and form ‘point’ source locations carrying potential environmental risk. The scale of such bodies may generate single-profile, ‘bulls-eye’ anomalies. In attempts to recover geological information, such anomalies may be considered to represent noise. In environmental AEM, the correct interpretation of such features is important. The study uses a combination of theoretical models and trial-fixed-wing survey data obtained in populated areas of the UK. Scale issues are discussed in terms of the volumetric footprints of the induced electric field generated by systems flown at both low and high elevation. One of the primary uses of AEM survey data lies in the assessment of conductivity maps. These are typically obtained using one-dimensional (1D) conductivity models at individual measurement points. In order to investigate the limitations of this approach, 3D modelling of conductive structures with dimensions less than 350×350 m and thicknesses extending to 20 m has been carried out. A 1D half space inversion of the data obtained at each frequency is then used to assess the behaviour of the spatial information. The results demonstrate that half space conductivity values obtained over compact 3D targets generally provide only apparent conductivity results. For thin, at-surface bodies, conductivity values are biased to lower values than the true conductivity except at high frequency. The spatial perturbation to both coupling ratios and 1D conductivity models can be laterally extensive. The results from 3D modelling indicate that the use of horizontal derivatives applied to the conductivity models offers enhanced edge detection. The practical application of such derivatives to both regional- and local-scale survey data is presented.. The special case of a near-surface, metallic pipeline has been modelled. The problem constitutes an inductive limit (current gathering) response in which the perturbation is largely confined to the in-phase coupling ratios. The main perturbations, in data and conductivity models, are within about 40 m of each side of the pipeline. The maximum perturbation to the conductivity model is only a factor of 1.5 above background. Detailed survey data across a former compact landfill (about 100×100 m) are used to compare the model behaviour predicted by the 3D modelling with survey results. The survey, conducted at two separate altitudes, provides a demonstration of 3D effects on 1D survey models as a function of frequency and elevation. Although the nature of the landfill materials and their location are not known precisely, the mapping information appears realistic.