核磁共振与瞬变电磁三维联合解释方法

传统核磁共振地下含水量解释多采用基于均匀半空间或层状导电模型的一维反演,分层给出地下含水信息.然而,这些方法忽略了地下复杂电阻率分布信息对结果的影响,也不能很好地反映局部三维含水构造.本文从三维电介质中核磁共振响应的正演理论出发,提出首先利用瞬变电磁数据进行基于等效导电平面法的快速电阻率成像,然后将成像结果作为核磁共振三维反演的电性模型,进行联合解释.激发磁场的分布采用有限元法直接求解,通过引入伪δ源实现电流源的加载,并强加散度条件排除了三维磁场模拟中"弱解"的影响.针对核磁共振灵敏度矩阵的病态性和数据中存在的干扰信号,提出考虑罚项的非线性拟合目标函数,利用线性化方法进行核磁共振反演.模型数据表明该方法能较准确反映地下三维含水构造,实测算例进一步证明了方法的有效性.本研究将促使核磁共振方法在岩溶、裂隙水、孤立水体等复杂水文地质条件及隧道、矿井灾害水源探测等方面得到有效应用.     

二度半长方体组合模型的重力模拟退火反演

针对重力异常反演的现状,采用矩形单元体网格剖分的二度半体组合模型,模拟任意的多层密度界面或形体,实现了重力正演.该建模方法适应性强,适合于实际工作中常见的密度横向变化的复杂模型.使用改进的全局寻优的快速模拟退火算法,对重力异常进行反演,结合这种灵活的密度建模方法,反演过程中只需要反演各矩形单元体的密度参数,即可同时得到地质体的界面或形态以及密度值分布.通过复杂界面和形体的模型试验,证明了方法的效果.在江苏下扬子地质构造复杂区的实际应用说明,在地震等先验信息约束下,该重力反演方法提高了反演精度并减少了多解性,可有效解决古生界泥盆系-志留系目的层分布等地质问题.     

基于非结构网格的电阻率三维带地形反演

地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.     

同时考虑退磁和剩磁的有限体积法正演模拟

为分析同时考虑退磁和剩磁对磁测数据解释的影响,探讨了利用有限体积法求静磁场数值解的方法.从静磁场中的麦克斯韦方程出发,导出了有限体积法控制方程的离散表达式,对边界条件近似处理后求解方程组得到磁异常.通过与退磁改正计算结果对比,验证了方法的正确性,并分析得到忽略剩磁的相对误差与科尼斯布格比（Q）相关;利用有限体积法计算长方体模型在无地磁场情况下的磁异常和内部磁化强度,从数值模拟上说明剩磁也需要进行退磁改正,并表明退磁作用对剩磁的影响不仅与磁化率相关,而且与剩磁的方向和磁性体的形态相关;组合模型的计算结果对比表明,退磁作用对剩磁的影响还会因为临近强磁性体的作用发生改变.在青海灶火河西工区的应用说明,开展同时考虑退磁和剩磁的解释方法对准确识别强磁性岩体具有实用价值.     

时间域航空电磁三维并行反演研究（英文）

为提高复杂条件下时间域航空电磁数据解释精度,本文开展了时间域航空电磁三维并行反演算法研究。该算法中的三维正演是基于有限差分技术,并采用"移动脚印"技术来减小实际计算模型尺寸;三维反演基于Gauss-Newton反演方法,并采用显式灵敏度矩阵计算技术减少反演过程中的正演次数。为提高三维反演的效率,本文基于OpenMP并行库实现了三维反演的并行化。从理论和实测数据的三维并行反演结果可以看出本文的并行化策略明显地提高了三维反演的速度,能够胜任大数据量时间域航空电磁实测资料三维反演解释任务。     

隧道磁共振全弛豫信号三维反演研究

磁共振探测反演近些年得到较多的关注.传统的磁共振反演多为基于地面大线圈采集,基于初始振幅的反演,这种方法未考虑地下复杂的多孔地质环境,且限制了磁共振方法的应用范围.针对上述问题,本文尝试建立了基于隧道磁共振全弛豫信号的三维反演方法,在传统的磁共振三维正演核函数基础上,加入指数衰减项,在充分考虑地下电阻率不均匀体的情况下,采用有限元方法模拟了磁共振全弛豫信号响应.针对三维磁共振灵敏度矩阵的病态性和隧道中的电磁干扰,建立了考虑罚项的最小二乘目标函数,并结合实际对函数的解加以约束,利用基于不等式约束的最优化算法实现了三维反演.模型计算表明:本文的三维反演效果良好,不仅能实现掌子面前方含水异常体定位,而且能给出不同弛豫时间下的含水量信息,从而间接获知灾害水赋存状态,为隧道磁共振数据处理提供了一个可供借鉴的方法.     

基于非结构有限元的时间域海洋电磁三维反演

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

考虑关断时间的地面瞬变电磁三维带地形反演

起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少每次反演迭代的计算量.针对发射线圈随地表起伏变化的特点,利用基于偶极子离散的场源处理技术模拟发射源的实际形状,采用瞬时电流脉冲技术实现考虑关断时间的地面瞬变电磁三维正演模拟.我们首先将本文开发的三维反演算法应用于理论模型的反演计算中,检验本文算法的可靠性,并分析地形和关断时间对反演结果的影响特征.在此基础上,进一步将本文算法应用于实测数据反演,验证本文算法的实用性.     

强剩磁强退磁条件下的二维井中磁测反演

强剩磁、强退磁改变了总磁化强度的大小和方向,给磁测资料解释带来困难.为此,本文利用二维井中磁测数据反演磁化强度矢量的二维分布.首先利用井中磁测的磁异常模量反演磁化强度大小的分布.然后,在已知磁化强度大小分布的前提下,拟合磁场分量,反演磁化强度方向的分布.其中,磁化强度大小和方向均用共轭梯度法求解,并通过预优矩阵改善磁化强度大小的反演效果.理论模拟说明,该方法能准确获得磁化强度矢量分布.磁化强度矢量反演结果包括感磁、剩磁及退磁的影响,这为研究强剩磁、高磁化率矿床提供了一种有效方法.     

可控源音频大地电磁数据的反演方法

从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一.     

Interactive 3D forward modeling of total field surface and three-component borehole magnetic data for the Daye iron-ore deposit (Central China)

The ground magnetic response of deep ore bodies in the Daye iron-ore deposit is relatively weak, and sometimes concealed by the strong magnetic background of shallower sources. Apart from the low-quality ground magnetic data, another critical problem for reconstructing the deep skarn-type ore bodies is developing a versatile inversion scheme that can simultaneously resolve 3D sources with arbitrary shapes. In this case, we resort to interactive 3D forward modeling solution with the joint use of two data sets-total field surface and three-component borehole magnetic data. Joint inversion of the two data sets is expected to help resolve the ambiguity associated with either data set and greatly reduces the nonuniqueness of the magnetic inversion. Such nonuniqueness is especially severe when a 3-D distribution of magnetic susceptibility, instead of a simple body, is sought from the inversion.In this paper, we calculate the magnetic field on the surface and in the borehole caused by 3D arbitrarily-shaped bodies with the triple integral method. The complex 3D magnetic sources having arbitrary shapes are constructed with cross-sections, termination points and facets in our visualization technology. We specify, interactively and in a user-friendly environment, the outline of the sources in terms of geometric elements and their magnetic parameters. The method automatically fits the observations within a prescribed precision. If dissatisfied, the user can redefine the model parameters and proceed to a new inversion. The method's ability to interpret a complicated 3D geologic environment is demonstrated on synthetic models and real data profiles in the Daye iron-ore deposit in central China. The interactive forward modeling results in all tests demonstrate a good correlation of estimated magnetic sources with corresponding known geologic features.     

Palaeomagnetic distortion modelling and possible recovery by inversion

A robust finite-element technique is presented for computation of both the internal demagnetization effects and magnetic terrain effects in bodies with arbitrary shape and arbitrary susceptibility distribution. This method facilitates a flexible analysis of the palaeomagnetic deflection problem. Tests on geologically realistic settings of highly magnetic rocks demonstrate that deflections of several degrees may occur even for relatively simple two-dimensional models. Similarly, the magnetic intensity may well be biased by 5-15% by demagnetization effects. The present paper focuses on deflections and intensity variations inside the magnetized body, where we find a systematic shallowing of inclination for bodies with a horizontal elongation. Because the bodies sampled at a typical palaeomagnetic site will have a dominant direction of elongation, the magnetic deflection effect will tend to impose a systematic bias which doesn’t average out. An inversion-based procedure for elimination of the deflection effect is presented. It requires that the magnetic body is quite homogeneous and that its surface geometry is known, as may be the case for historical lava flows. Tests demonstrate that in order to recover both ambient palaeofield direction and the effective susceptibility at blocking temperature it is necessary to sample near strong topographic elements in the magnetic body. Since the surface geometry rarely is known it is proposed as an alternative to inversion that an effective susceptibility is assessed and a horizontal slab correction is applied for samples taken far from topographical features. When shape geometry is unknown and no correction applied, palaeomagnetic conclusions must take into account the possible bias from internal demagnetization and magnetic terrain effects.     

强磁性三度体内部磁化强度的分布

本文讨论稳定磁化场中强磁性三度体内部磁化强度的数值解法。展示了由于退磁作用而造成的不均匀性,并以长方体内一些典型截面为例,较详细地讨论了磁化强度的特征,作了分布图;推导了较均匀磁化体磁化强度分布的近似公式,并计算了一组长方体的视退磁系数。     

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

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

3D inversion of gravity data by separation of sources and the method of local corrections: Kolarovo gravity high case study

We present a novel methodology for 3D gravity/magnetic data inversion. It combines two algorithms for preliminary separation of sources and an original approach to 3D inverse problem solution. The first algorithm is designed to separate sources in depth and to remove the shallow ones. It is based on subsequent upward and downward data continuation. For separation in the lateral sense, we approximate the given observed data by the field of several 3D line segments. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. It allows retrieving unknown 3D geometry of anomalous objects in terms of restricted bodies of arbitrary shape and contact surfaces. For restricted objects, we apply new integral equations of gravity and magnetic inverse problems. All steps of our methodology are demonstrated on the Kolarovo gravity anomaly in the Danube Basin of Slovakia.     

复杂介质有限元法2.5维可控源音频大地电磁法数值模拟

利用有限元2.5维可控源音频大地电磁法（简称CSAMT）数值模拟方法,对100Ωm均匀半空间介质中有限长度的电偶极源产生的电场、磁场及视电阻率、相位特征进行了数值模拟,研究了场的空间变化规律. 在一个象限中,场的特征存在双叶现象,当收发距大于4个趋肤深度时,电阻率较接近介质真实的电阻率,这些结果为观测系统和收发距的选择提供了依据. 波数域场的特征表明,低波数对源的贡献占较大的比例. 有限元法2.5维CSAMT数值模拟的优势在于能较准确地获得复杂介质结构的波场特征. 本文结合直立异常体、倾斜异常体及断陷模型对CSAMT电阻率、相位剖面特征及频率曲线特征的可靠性进行了研究. 数值模拟结果直观地给出了异常体的剖面异常形态. 通过对比研究异常体的剖面异常形态和半空间场的特征进一步说明本文方法和软件在模拟复杂介质结构场特征时是可靠的. 这为认识观测数据,指导反演解释提供了较好的依据.     

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

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

Magnetic interface forward and inversion method based on Padé approximation

The magnetic interface forward and inversion method is realized using the Taylor series expansion to linearize the Fourier transform of the exponential function. With a large expansion step and unbounded neighborhood, the Taylor series is not convergent, and therefore, this paper presents the magnetic interface forward and inversion method based on Padé approximation instead of the Taylor series expansion. Compared with the Taylor series, Padé’s expansion’s convergence is more stable and its approximation more accurate. Model tests show the validity of the magnetic forward modeling and inversion of Padé approximation proposed in the paper, and when this inversion method is applied to the measured data of the Matagami area in Canada, a stable and reasonable distribution of underground interface is obtained.     

Calculation of magnetic anomalies caused by 2D bodies of arbitrary shape with consideration of demagnetization

Forward calculations of magnetic anomalies caused by two-dimensional bodies of any shape and magnetic properties may be performed either without considering demagnetization as in the equivalent source technique or taking demagnetization into account as in the volume integral equation (VIE) approach, in which, for this purpose, magnetized bodies are divided into a set of rectangular prismatic cells. Ignoring demagnetization may result in distortion of the shape and the amplitude of an anomaly, whereas rectangular cells may not be an optimal representation of the source. Moreover, an inaccurate form approximation in the VIE technique may lead to inconsistent results in the near-body region. In this paper, a method is proposed, based on the VIE approach but differing by applying triangular elementary cells. The method largely overcomes the above-mentioned limitations of the VIE technique. It allows us to delineate large and complex structures exactly and only requires the source to be divided into a few elementary cells to take demagnetization into account satisfactorily. These improvements have been attained through analytical calculation of the Green's function in the complex plane, using the theory of the Cauchy-type integral. Comparing numerical solutions with analytical solutions for homogeneous elliptic cylinders without remanence, the method is found to be consistent with the theory in the range of relative magnetic permeability of 2–20, not only far from but also at subcell distances from the body. The method is appropriate for modelling highly and inhomogeneously magnetized 2D bodies of any shape. It may be of value in interpreting underground measurements or topographic effects, as well as in modelling regional geomagnetic profiles, and it is also a convenient tool for testing questionable geological hypotheses. In the framework of the method, the gravitational anomaly for the same causative bodies can be easily calculated. However, at higher and geologically uncommon values of relative magnetic permeability, the algorithm may become unstable but may be stabilized with SVD regularization. The fact that discrepancies were found with the method employed is a basis for further research.