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

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

A regularized Wiener–Hopf filter for inverting models with magnetic susceptibility

For a magnetic target, the spatial magnetic signal can be expressed as a convolutional integral over Green's function of an assumed model with susceptibility as its parameter. A filter can be used to obtain the susceptibility by minimizing the mismatch between observed and the computed magnetic anomalies. In this perspective, we report the development of an advanced digital filter, which is efficient and can be used to map rock susceptibility from the acquired magnetic data. To design the new filter, we modified the space‐domain standard Wiener–Hopf filter by imposing two different constraints: (i) the filter energy constraint; and (ii) normalization of the filter coefficients. These constraints make it capable to characterize source bodies from their produced magnetic anomalies. We assume that the magnetic data are produced by induced magnetization only and interpretation can be as good as the subsurface model. Our technique is less sensitive to the data noise, which makes it efficient in enhancing the interpretation model. The modified filter demonstrates its applicability over the synthetic data with additive white Gaussian noise. In order to check the efficacy and adaptivity of this tool in a more realistic perspective, it is also tested on the real magnetic data acquired over a kimberlitic district adjoining to the western margin of the Cuddapah Basin in India to identify the source bodies from the anomalies. Our result shows that the modified Wiener–Hopf filter with the constraint for the magnetic data is more stable and efficient than the standard Wiener–Hopf filter.     

南海重磁异常特征及火成岩分布(英文)

南海火成岩油气藏具有广阔的勘探前景,综合利用地球物理方法圈划与识别火成岩体、研究火成岩分布是火成岩油气藏研究的基础。针对南海重磁场特征,采用低纬度、变倾角化极技术进行了磁异常化极处理,利用优选延拓方法实现重磁异常分离并提取南海海域浅部火成岩重磁异常信息,利用磁异常三维相关成像给出南海火成岩的三维空间等效分布,在重磁梯度突出局部异常边界信息的基础上,通过梯度加权的重磁相关分析勾画不同类型火成岩的平面展布,火成岩的分布特征显示出受地壳深部结构及断裂构造的控制与影响。     

基于重磁异常相关分析的场源位置及属性识别方法

针对传统的重磁对应分析在场源区域外趋于高相关度的缺点,本文提出了一种基于相关系数和垂向导数的重磁场源位置及属性识别方法,给出了一个新的重磁场源平面位置及属性判别参数(GMCP),该参数能够有效地缩小重磁场源的识别范围.判别参数GMCP值非零分布范围反映了场源的规模.GMCP值的正负反映了场源属性,当GMCP大于零时为正相关,代表了该区域存在高磁高密度或低磁低密度的重磁同源地质体;当GMCP小于零时为负相关,代表了该区域存在高磁低密度或低磁高密度的重磁同源地质体;GMCP趋于零,则表示无重磁同源地质体.两组不同干扰水平的同源复杂组合模型实验、两组不同源组合模型试验以及南海中央海盆实际重磁数据处理结果均验证了该重磁场源位置及属性识别方法的正确性和有效性.     

补偿向下延拓方法研究及应用

位场向下延拓是重、磁处理和解释的常用方法,但其不稳定性限制了其在资料处理及反演中的应用.本文基于补偿圆滑滤波思想以及空间域向下延拓迭代法,通过逐次补偿的办法实现位场的稳定向下延拓.同时,在频率域空间给出了该下延方法的频率域响应因子,并讨论了其低通滤波特性,理论模型和实际位场资料试验表明该方法向下延拓稳定性具有较高的延拓精度.将其应用于重力密度界面反演中,改进反演的稳定性,实际莫霍界面反演表明下延因子具备实用性.     

磁异常ΔT三维相关成像

本文将重力和重力梯度数据三维相关成像方法推广到磁力勘探领域,推导并建立了磁异常ΔT的三维相关成像方法,同时提出了基于熵滤波分离异常的三维相关成像算法来提高成像分辨率.组合模型磁异常ΔT数据和实际磁测资料试验分析表明,本文方法能成像出地下地质体的空间赋存状态和等效磁性分布,具有良好的横向和纵向分辨率.     

Properties of an improved Gabor wavelet transform and its applications to seismic signal processing and interpretation

This paper presents an analytical study of the complete transform of improved Gabor wavelets (IGWs), and discusses its application to the processing and interpretation of seismic signals. The complete Gabor wavelet transform has the following properties. First, unlike the conventional transform, the improved Gabor wavelet transform (IGWT) maps time domain signals to the time-frequency domain instead of the time-scale domain. Second, the IGW’s dominant frequency is fixed, so the transform can perform signal frequency division, where the dominant frequency components of the extracted sub-band signal carry essentially the same information as the corresponding components of the original signal, and the subband signal bandwidth can be regulated effectively by the transform’s resolution factor. Third, a time-frequency filter consisting of an IGWT and its inverse transform can accurately locate target areas in the time-frequency field and perform filtering in a given time-frequency range. The complete IGW transform’s properties are investigated using simulation experiments and test cases, showing positive results for seismic signal processing and interpretation, such as enhancing seismic signal resolution, permitting signal frequency division, and allowing small faults to be identified.     

磁异常△T三维相关成像

本文将重力和重力梯度数据三维相关成像方法推广到磁力勘探领域,推导并建立了磁异常△T的三维相关成像方法,同时提出了基于熵滤波分离异常的三维相关成像算法来提高成像分辨率.组合模型磁异常△T数据和实际磁测资料试验分析表明,本文方法能成像出地下地质体的空间赋存状态和等效磁性分布,具有良好的横向和纵向分辨率.     

21世纪重力与磁法勘探的展望

对21世纪重力与磁法勘探的仪器，数据处理技术，解释理论与方法，应用领域等方面的发展方向进行了分析与展望，发展航空标量，矢量，梯度重力测量和航空全梯度磁力测量，三分量磁力测量，提高综合信息采集能力；开展卫星重磁测量，综合卫星，航空，地面重磁测量资料研究地球结构与构造；发展高精度数据处理技术；重磁异常弱信号的提取，不同深度重磁异常的划分，低纬底变倾角化磁极以及位场面延拓；发展复杂条件下三维重磁场多参数综合反演可视化技术以及快速自动反演技术；探索磁性多参数的应用新领域，充分发挥磁法在环境污染调查中的作用并开拓应用新领域。     

地磁匹配用于水下载体导航的初步分析

随着地磁场模型的日趋完善和滤波技术的不断成熟,地磁匹配方法在近几年得到迅速发展,它一般用于近地卫星,还未用于海洋环境.针对这一问题,本文结合地磁场特点和海洋地磁测量的现状,分析了水下地磁匹配的可行性.分别从匹配的地磁模型和测量值中存在干扰磁场两方面,指出了制约该技术发展的几个因素.在此基础上,提出了地磁用于水下载体定位的两条技术途径,一是用小波变换滤去地磁异常,与大尺度地磁模型匹配定位;另外,也提出了利用地磁异常,结合同时定位与构图(SLAM)算法实现自主导航.     

位场异常分离的递减半径迭代法

本文提出了位场区域-剩余异常分离的空间域递减半径迭代法.由给定半径的圆周上八点位场值的算术平均导出一个新的八点圆周平均公式,它是一个由重磁数据计算区域异常的滤波器.该滤波器的传递函数有一主瓣和多个旁瓣,且半径越大,旁瓣越多,滤波器特性越差.由从大到小不同半径对应的传递函数的连乘积构造了递减半径迭代传递函数,它以大半径为参变量.递减半径迭代传递函数类似低通滤波器,其截止波数与大半径成反比.由递减半径迭代传递函数,给出空间域分离区域异常的递减半径线性迭代法,由重磁数据减去区域异常求得剩余异常.进一步,通过构造非线性修正系数,把空间域递减半径线性迭代法中线性迭代公式变成非线性迭代公式,得到空间域递减半径非线性迭代法.重磁理论模型及安徽泥河铁矿重磁资料的试验表明,空间域递减半径非线性迭代法可有效地压制假异常和高频干扰以及减少异常畸变,并能有效地分离出区域异常和剩余异常.     

磁总场异常垂直梯度三维相关成像

本文提出了磁总场异常垂直梯度三维相关成像方法,用于成像地下等效磁源分布.它首先将地下待成像空间划分为三维规则网格,然后直接计算每个网格节点磁偶板子在观测面理论磁总场异常垂直梯度与实测磁总场异常垂直梯度的互相关,其相关系数值表征等效磁偶极子分布(即磁偶极子发生的概率).理论长方体组合模型数据和实际某矿区磁测资料试验结果表明本文方法计算得到的相关系数值能基本反映地下的磁源分布,且分辨率明显高于磁总场异常三维相关成像的分辨率,也高于基于熵滤波分离异常的磁总场异常三维相关成像的分辨率.     

利用加强水平方向总水平导数对位场全张量数据进行边界识别

位场全张量梯度数据以其信息量大、含有更高频的信号成分,能更好地描述小的异常特征等优点在地球物理领域中得到广泛应用.边界检测是位场解释中不可缺少的任务,需要新的边界探测器来处理位场梯度张量数据.为了充分利用位场梯度张量数据的多信息成分,本文定义了方向总水平导数和加强方向总水平导数,并利用其定义新的边界检测器.为了能同时显示不同振幅大小异常的边界,本文对其进行了归一化处理.通过模型试验,证明了归一化方法能更加清晰准确地显示浅部和深部的地质体边界信息.最后将该边界检测方法用于加拿大圣乔治湾实际测得全张量重力梯度数据和中国朱日和地区的磁异常数据中,并得到了较好的边界检测结果.     

磁赤道处化极方法

化向地磁极(化极)是最基本的磁测资料处理方法之一,化极能消除或减少斜磁化影响,提高对磁测资料的认识程度和解释水平,对研究地壳产生的磁异常具有重要意义.但低纬度地区特别是磁赤道处,化极处理很不稳定甚至奇异,一直是位场研究的难点.针对地磁纬度较低特别是磁赤道地区磁异常化极的困难,利用从磁北极处垂直磁化向低纬度地区水平磁化方向转换稳定的特点,提出"狭义化赤"概念,并将其与低纬度磁异常"倒相"解释方法结合,提出专门用于磁赤道处化极的方法.该方法扩展了现有的化极理论,实现了磁赤道处的稳定化极.区别于目前任何方法,专门用于(近)水平磁化条件下的化极计算,具有原理简单,实现方便,收敛速度快等特点.对理论模型和实际资料计算表明这种针对磁赤道地区磁异常的化极处理方法是稳定、可靠的.     

强磁性体ΔT异常计算的误差分析研究

本文主要针对当前磁法勘探中高精度处理解释的需求,对强磁性体ΔT异常计算存在的误差进行分析研究.我们首先通过理论模型计算试验,证明常规计算采用的投影关系的ΔT与实际测量的模量差ΔT之间的误差E在磁异常幅值大时是明显存在的,其影响不容忽视.其次,当磁性强且剩磁存在时,投影ΔT曲线及其误差曲线在磁化方向与地磁场方向改变时具有一定的对称性;地磁场T0、磁性体形态(如二度水平圆柱体模型的半径r、柱体埋深R)和磁性参数(如磁化率κ)等参数确定的情况下,最大误差值出现在磁性体正上方,且其大小与磁性参数(κ)和模型体规模(如r/R)之间皆是指数关系;另外,研究还发现ΔT的计算误差曲线的一些其他规律特点,如在各纬度带上,ΔT计算误差的最大值Emax曲线的极值主要分布在中纬度地区;磁异常矢量Ta与地磁场T0的夹角θ逐渐变化时,随θ变化Emax曲线的极值分布在θ=90°~120°范围内;当磁异常幅值小于10000nT时,最大误差近似为磁异常矢量垂直于地磁场方向的测点附近的误差值;另外,磁性体(圆柱体为例)的半径(即尺度)与埋深的比值r/R超过0.5,且磁化率超过0.1SI时误差已达到3.9nT,磁化率增大与对应的Emax的值呈指数增长特点.因此,我们的研究表明,在强磁性体、磁异常幅值大的数据处理、反演及解释时,现有方法会产生较大的误差,应该基于严格的模量差ΔT,完善相应的处理以及反演方法.     

Edge Detection in Potential-Field Data by Enhanced Mathematical Morphology Filter

Enhancement on the edges of the causative source is an indispensable tool in the interpretation of potential-field data. There are a number of methods for recognizing the edges, most of which involve high-pass filters based on derivatives of potential-field data. A new edge-detection method is presented, called the enhanced mathematical morphology (EMM) filter. The EMM filter uses the ratio of the erosion of the total horizontal derivative to the dilation of total horizontal derivative to recognize the edges of the sources, and can display the edges of the shallow and deep bodies simultaneously. The EMM filter does not require the computation of vertical derivatives, which makes this method computationally stable. The EMM filter is tested on synthetic and real potential field data in China. Compared to other edge-detection filters, the new method is able to recognize the source edges more clearly, and the outputs are more insensitive to noise.     

位场向下延拓组合滤波器的设计和应用

一 重、磁异常向下延拓能将水平方向迭加的异常分开,并能突出异常细节,从而丰富异常的解释内容,提高异常解释的可靠性.本文从对引起向下延拓不稳定的各种因素的讨论入手,结合对一些向下延拓方法在频率域的分析,提出了一种与已知资料频谱特征相“匹配”的组合滤波器。经理论模型和实例考核,说明该滤波器既能保证向下延拓有较高的精度,又有较强的抗干扰能力。     

Application of balanced edge detection filters to estimate the location parameters of the causative sources using potential field data

Balanced edge detection filters can recognize the edges of the shallow and deep bodies simultaneously, and are commonly used in the edge detection of potential field data. In this paper, we present using the balanced edge detection filters to estimate source locations, and derive two linear equations based on the balanced edge detection filters that can estimate the locations of the source without any priori information about the nature (structural index) of the source. The proposed methods are demonstrated on synthetic gravity anomalies, and the inversion results show that the proposed methods can successfully estimate location parameters of the sources. I also apply the proposed methods to real magnetic data, and the inversion results estimated by the proposed methods are consistent with the results estimated by the other similar method.     

LINEAR FILTERING OF VLF DIP-ANGLE MEASUREMENTS*

The suggested interpretation technique is based on discrete linear filtering of VLF data. The output of the described filtering results is expressed in terms of an equivalent current density at a specified depth that would cause the measured magnetic field. The most practical six-point filter gives an accuracy of 8%. The filter is an extension of the commonly used Fraser filter to process VLF dip-angle data. Filtering the same data set for various depths gives an idea about the change of current density with depth. Areas with high current-density correspond to good conductors. The conductor dip can also be determined. The use of the technique is illustrated on theoretical and field examples. In all cases a good correlation with original models and other types of geophysical measurements was obtained. As shown in the last example, the filtering technique is also applicable in interpretation of other electromagnetic methods.     

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.