Application of Hilbert‐like transforms for enhanced processing of full tensor magnetic gradient data

Commonly, geomagnetic prospection is performed via scalar magnetometers that measure values of the total magnetic intensity. Recent developments of superconducting quantum interference devices have led to their integration in full tensor magnetic gradiometry systems consisting of planar‐type first‐order gradiometers and magnetometers fabricated in thin‐film technology. With these systems measuring directly the magnetic gradient tensor and field vector, a significantly higher magnetic and spatial resolution of the magnetic maps is yield than those produced via conventional magnetometers. In order to preserve the high data quality in this work, we develop a workflow containing all the necessary steps for generating the gradient tensor and field vector quantities from the raw measurement data up to their integration into high­resolution, low­noise, and artefactless two‐dimensional maps of the magnetic field vector. The gradient tensor components are processed by superposition of the balanced gradiometer signals and rotation into an Earth‐centred Earth‐fixed coordinate frame. As the magnetometers have sensitivity lower than that of gradiometers and the total magnetic intensity is not directly recorded, we employ Hilbert‐like transforms, e.g., integration of the gradient tensor components or the conversion of the total magnetic intensity derived by calibrated magnetometer readings to obtain these values. This can lead to a better interpretation of the measured magnetic anomalies of the Earth's magnetic field that is possible from scalar total magnetic intensity measurements. Our conclusions are drawn from the application of these algorithms on a survey acquired in South Africa containing full tensor magnetic gradiometry data.     

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

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

Estimating magnetic dike parameters using a non‐linear constrained inversion technique: an example from the Särna area,west central Sweden

In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated g_zz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of g_zz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.     

Applying the tilt‐depth and contact‐depth methods to the magnetic anomalies of thin dykes

A new method for the calculation of the depth, location, and dip of thin dykes from pole‐reduced magnetic data is introduced. The depth can be obtained by measuring the distance between chosen values of a tilt angle that is based upon the ratio of the magnetic field and its Hilbert transform over the dyke. Alternatively, it can be obtained from the horizontal derivative of the ratio of the Hilbert transform of the field to the field itself, over the dyke. The latter method also allows the dip of the dyke to be estimated from the gradient of the depth estimates.     

Feasibility study for an anisotropic full waveform inversion of cross‐well seismic data

Anisotropy is often observed due to the thin layering or aligned micro‐structures, like small fractures. At the scale of cross‐well tomography, the anisotropic effects cannot be neglected. In this paper, we propose a method of full‐wave inversion for transversely isotropic media and we test its robustness against structured noisy data. Optimization inversion techniques based on a least‐square formalism are used. In this framework, analytical expressions of the misfit function gradient, based on the adjoint technique in the time domain, allow one to solve the inverse problem with a high number of parameters and for a completely heterogeneous medium. The wave propagation equation for transversely isotropic media with vertical symmetry axis is solved using the finite difference method on the cylindrical system of coordinates. This system allows one to model the 3D propagation in a 2D medium with a revolution symmetry. In case of approximately horizontal layering, this approximation is sufficient. The full‐wave inversion method is applied to a crosswell synthetic 2‐component (radial and vertical) dataset generated using a 2D model with three different anisotropic regions. Complex noise has been added to these synthetic observed data. This noise is Gaussian and has the same amplitude f?k spectrum as the data. Part of the noise is localized as a coda of arrivals, the other part is not localized. Five parameter fields are estimated, (vertical) P‐wave velocity, (vertical) S‐wave velocity, volumetric mass and the Thomsen anisotropic parameters epsilon and delta. Horizontal exponential correlations have been used. The results show that the full‐wave inversion of cross‐well data is relatively robust for high‐level noise even for second‐order parameters such as Thomsen epsilon and delta anisotropic parameters.     

High‐accuracy practical spline‐based 3D and 2D integral transformations in potential‐field geophysics

Potential, potential field and potential‐field gradient data are supplemental to each other for resolving sources of interest in both exploration and solid Earth studies. We propose flexible high‐accuracy practical techniques to perform 3D and 2D integral transformations from potential field components to potential and from potential‐field gradient components to potential field components in the space domain using cubic B‐splines. The spline techniques are applicable to either uniform or non‐uniform rectangular grids for the 3D case, and applicable to either regular or irregular grids for the 2D case. The spline‐based indefinite integrations can be computed at any point in the computational domain. In our synthetic 3D gravity and magnetic transformation examples, we show that the spline techniques are substantially more accurate than the Fourier transform techniques, and demonstrate that harmonicity is confirmed substantially better for the spline method than the Fourier transform method and that spline‐based integration and differentiation are invertible. The cost of the increase in accuracy is an increase in computing time. Our real data examples of 3D transformations show that the spline‐based results agree substantially better or better with the observed data than do the Fourier‐based results. The spline techniques would therefore be very useful for data quality control through comparisons of the computed and observed components. If certain desired components of the potential field or gradient data are not measured, they can be obtained using the spline‐based transformations as alternatives to the Fourier transform techniques.     

Multiridge Euler deconvolution

Potential field interpretation can be carried out using multiscale methods. This class of methods analyses a multiscale data set, which is built by upward continuation of the original data to a number of altitudes conveniently chosen. Euler deconvolution can be cast into this multiscale environment by analysing data along ridges of potential fields, e.g., at those points along lines across scales where the field or its horizontal or vertical derivative respectively is zero. Previous work has shown that Euler equations are notably simplified along any of these ridges. Since a given anomaly may generate one or more ridges we describe in this paper how Euler deconvolution may be used to jointly invert data along all of them, so performing a multiridge Euler deconvolution. The method enjoys the stable and high‐resolution properties of multiscale methods, due to the composite upward continuation/vertical differentiation filter used. Such a physically‐based field transformation can have a positive effect on reducing both high‐wavenumber noise and interference or regional field effects. Multiridge Euler deconvolution can also be applied to the modulus of an analytic signal, gravity/magnetic gradient tensor components or Hilbert transform components. The advantages of using multiridge Euler deconvolution compared to single ridge Euler deconvolution include improved solution clustering, increased number of solutions, improvement of accuracy of the results obtainable from some types of ridges and greater ease in the selection of ridges to invert. The multiscale approach is particularly well suited to deal with non‐ideal sources. In these cases, our strategy is to find the optimal combination of upward continuation altitude range and data differentiation order, such that the field could be sensed as approximately homogeneous and then characterized by a structural index close to an integer value. This allows us to estimate depths related to the top or the centre of the structure.     

Seismic channel edge detection using 3D shearlets—a study on synthetic and real channelised 3D seismic data

Automatic feature detection from seismic data is a demanding task in today's interpretation workstations. Channels are among important stratigraphic features in seismic data both due to their reservoir capability or drilling hazard potential. Shearlet transform as a multi‐scale and multi‐directional transformation is capable of detecting anisotropic singularities in two and higher dimensional data. Channels occur as edges in seismic data, which can be detected based on maximizing the shearlet coefficients through all sub‐volumes at the finest scale of decomposition. The detected edges may require further refinement through the application of a thinning methodology. In this study, a three‐dimensional, pyramid‐adapted, compactly supported shearlet transform was applied to synthetic and real channelised, three‐dimensional post‐stack seismic data in order to decompose the data into different scales and directions for the purpose of channel boundary detection. In order to be able to compare the edge detection results based on three‐dimensional shearlet transform with some famous gradient‐based edge detectors, such as Sobel and Canny, a thresholding scheme is necessary. In both synthetic and real data examples, the three‐dimensional shearlet edge detection algorithm outperformed Sobel and Canny operators even in the presence of Gaussian random noise.     

重磁场二维希尔伯特变换——直接解析信号解释方法

通过分析解析信号概念,指出目前重磁场解析信号事实上是重磁场梯度解析信号.在借助解析信号分量满足二维希尔伯特变换关系的基础上,提出重磁场直接解析信号的概念,并阐述重力异常及化极磁异常希尔伯特变换——直接解析信号的含义,并给出基于直接解析信号对位场增强的四种处理方法:直接解析信号模、水平分量模、改进Tilt angle和改...     

利用Hilbert变换计算重力归一化总梯度

针对提高重力勘探正反演解释的分辨率问题,提出利用Hilbert变换计算和研究重力归一化总梯度.文中从理论上详细地证明了方法的可行性,给出了适合计算机实现的计算方法.为探讨该方法对油气藏的分辨能力,在模型计算中,分别利用Fourier级数法、Fourier变换法和本文提出的Hilbert变换法计算模型的GH场值,发现当三度体储油球冠模型（模拟似三度贮油气藏背斜模型）油气藏部分（低密度体）厚度减小到低于球冠厚度的十分之一时,Fourier级数及变换法不能分辨出低密度体所产生的异常,而Hilbert变换法仍然可以清晰地识别,这说明用Hilbert变换法计算的GH场对异常识别的分辨率优于其他两种方法.     

Surface nuclear magnetic resonance signals recovery by integration of a non‐linear decomposition method with statistical analysis

Presence of noise in the acquisition of surface nuclear magnetic resonance data is inevitable. There are various types of noise, including Gaussian noise, spiky events, and harmonic noise that affect the signal quality of surface nuclear magnetic resonance measurements. In this paper, we describe an application of a two‐step noise suppression approach based on a non‐linear adaptive decomposition technique called complete ensemble empirical mode decomposition in conjunction with a statistical optimization process for enhancing the signal‐to‐noise ratio of the surface nuclear magnetic resonance signal. The filtering procedure starts with applying the complete ensemble empirical mode decomposition method to decompose the noisy surface nuclear magnetic resonance signal into a finite number of intrinsic mode functions. Afterwards, a threshold region based on de‐trended fluctuation analysis is defined to identify the noisy intrinsic mode functions, and then the no‐noise intrinsic mode functions are used to recover the partially de‐noised signal. In the second stage, we applied a statistical method based on the variance criterion to the signal obtained from the initial phase to mitigate the remaining noise. To demonstrate the functionality of the proposed strategy, the method was evaluated on an added‐noise synthetic surface nuclear magnetic resonance signal and on field data. The results show that the proposed procedure allows us to improve the signal‐to‐noise ratio significantly and, consequently, extract the signal parameters (i.e., and V₀) from noisy surface nuclear magnetic resonance data efficiently.     

The theory of the continuous wavelet transform in the interpretation of potential fields: a review

We consider the use of the continuous wavelet transform in the interpretation of potential field data. We report its development since the publication of the first paper by Moreau et al . in 1997. Basically, it consists in the interpretation in the upward continued domain since dilation of the wavelet transform is the upward continuation altitude. Thus within a range of altitudes, the wavelet transform of the noise is decreased faster than the wavelet transform of the potential field caused by underground sources; this means that the signal-to-noise ratio is much better than those involved in other enhancing methods (e.g., Euler deconvolution, gradient analysis, or the analytic signals). Similarly to the Euler deconvolution, its first target parameters were the source positions and shape. The method has then been developed to estimate size and directions of extended sources (e.g., faults and dikes of finite dimensions) and also the magnetization direction in the case of magnetic data. Latest developments show that when combined with a Radon transform, the continuous wavelet transform can help in the automatic detection of elongated structures in 3D, simultaneously to the estimation of their strike direction, shape and depth. Several applications to real case studies have been shown before; however for clarity's sake in the present paper, only synthetic cases have been reproduced to clearly sum up the development of the methodology.     

Interpretation of magnetic and gravity gradient tensor data using normalized source strength – A case study from McFaulds Lake,Northern Ontario,Canada

In this paper, we present a case study on the use of the normalized source strength (NSS) for interpretation of magnetic and gravity gradient tensors data. This application arises in exploration of nickel, copper and platinum group element (Ni‐Cu‐PGE) deposits in the McFaulds Lake area, Northern Ontario, Canada. In this study, we have used the normalized source strength function derived from recent high resolution aeromagnetic and gravity gradiometry data for locating geological bodies. In our algorithm, we use maxima of the normalized source strength for estimating the horizontal location of the causative body. Then we estimate depth to the source and structural index at that point using the ratio between the normalized source strength and its vertical derivative calculated at two levels; the measurement level and a height h above the measurement level. To discriminate more reliable solutions from spurious ones, we reject solutions with unreasonable estimated structural indices. This method uses an upward continuation filter which reduces the effect of high frequency noise. In the magnetic case, the advantage is that, in general, the normalized magnetic source strength is relatively insensitive to magnetization direction, thus it provides more reliable information than standard techniques when geologic bodies carry remanent magnetization. For dipping gravity sources, the calculated normalized source strength yields a reliable estimate of the source location by peaking right above the top surface. Application of the method on aeromagnetic and gravity gradient tensor data sets from McFaulds Lake area indicates that most of the gravity and magnetic sources are located just beneath a 20 m thick (on average) overburden and delineated magnetic and gravity sources which can be probably approximated by geological contacts and thin dikes, come up to the overburden.     

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

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

The feasibility of electromagnetic gradiometer measurements

The quantities measured in transient electromagnetic (TEM) surveys are usually either magnetic field components or their time derivatives. Alternatively it might be advantageous to measure the spatial derivatives of these quantities. Such gradiometer measurements are expected to have lower noise levels due to the negative interference of ambient noise recorded by the two receiver coils. Error propagation models are used to compare quantitatively the noise sensitivities of conventional and gradiometer TEM data. To achieve this, eigenvalue decomposition is applied on synthetic data to derive the parameter uncertainties of layered-earth models. The results indicate that near-surface gradient measurements give a superior definition of the shallow conductivity structure, provided noise levels are 20–40 times smaller than those recorded by conventional EM instruments. For a fixed-wing towed-bird gradiometer system to be feasible, a noise reduction factor of at least 50–100 is required. One field test showed that noise reduction factors in excess of 60 are achievable with gradiometer measurements. However, other collected data indicate that the effectiveness of noise reduction can be hampered by the spatial variability of noise such as that encountered in built-up areas. Synthetic data calculated for a vertical plate model confirm the limited depth of detection of vertical gradient data but also indicate some spatial derivatives which offer better lateral resolution than conventional EM data. This high sensitivity to the near-surface conductivity structure suggests the application of EM gradiometers in areas such as environmental and archaeological mapping.     

Continuous wavelet transform,theoretical aspects and application to aeromagnetic data at the Huanghua Depression,Dagang Oilfield,China

We use the continuous wavelet transform based on complex Morlet wavelets, which has been developed to estimate the source distribution of potential fields. For magnetic anomalies of adjacent sources, they always superimpose upon each other in space and wavenumber, making the identification of magnetic sources problematic. Therefore, a scale normalization factor, a^?n, is introduced on the wavelet coefficients to improve resolution in the scalogram. By theoretical modelling, we set up an approximate linear relationship between the pseudo‐wavenumber and source depth. The influences of background field, random noise and magnetization inclination on the continuous wavelet transform of magnetic anomalies are also discussed and compared with the short‐time Fourier transform results. Synthetic examples indicate that the regional trend has little effect on our method, while the influence of random noise is mainly imposed on shallower sources with higher wavenumbers. The source horizontal position will be affected by the change of magnetization direction, whereas the source depth remains unchanged. After discussing the performance of our method by showing the results of various synthetic tests, we use this method on the aeromagnetic data of the Huanghua depression in central China to define the distribution of volcanic rocks. The spectrum slices in different scales are used to determine horizontal positions of volcanic rocks and their source depths are estimated from the modulus maxima of complex coefficients, which is in good accordance with drilling results.     

Near‐field seismic effects in a homogeneous medium and their removal in vertical seismic profile attenuation estimates

To better understand (and correct for) the factors affecting the estimation of attenuation (Q), we simulate subsurface wave propagation with the Weyl/Sommerfeld integral. The complete spherical wavefield emanating from a P‐wave point source surrounded by a homogeneous, isotropic and attenuative medium is thus computed. In a resulting synthetic vertical seismic profile, we observe near‐field and far‐field responses and a 90° phase rotation between them. Depth dependence of the magnitude spectra in these two depth regions is distinctly different. The logarithm of the magnitude spectra shows a linear dependence on frequency in the far‐field but not in those depth regions where the near‐field becomes significant. Near‐field effects are one possible explanation for large positive and even negative Q‐factors in the shallow section that may be estimated from real vertical seismic profile data when applying the spectral ratio method. We outline a near‐field compensation technique that can reduce errors in the resultant Q estimates.     

航磁垂直梯度调整ΔT水平方法研究

使用中国自主研发、具有完全自主知识产权的航磁全轴梯度勘查系统首次试飞获得的资料,利用航磁垂直梯度数据进行ΔT水平调整研究,讨论了航磁及航磁梯度资料水平调整的原则,提出梯度资料长波水平的调平方法.对经长波水平调整后的垂直梯度数据进行位场转换,获得不受地磁场随时间变化影响的磁异常,并通过恢复磁异常的长波信息,实现了ΔT的水平调整和增强.实际资料处理表明调平过程具有物理意义、效果明显.     

Fast imaging of subsurface conductors using very low‐frequency electromagnetic data

The study presents a fast imaging technique for the very low‐frequency data interpretation. First, an analytical expression was derived to compute the vertical component of the magnetic field at any point on the Earth's surface for a given current density distribution in a rectangular block on the subsurface. Current density is considered as exponentially decreasing with depth, according to the skin depth rule in a particular block. Subsequently, the vertical component of the magnetic field due to the entire subsurface was computed as the sum of the vertical component of the magnetic field due to an individual block. Since the vertical component of the magnetic field is proportional to the real part of very low‐frequency anomaly, an inversion program was developed for imaging of the subsurface conductors using the real very low‐frequency anomaly in terms of apparent current density distribution in the subsurface. Imaging results from the presented formulation were compared with other imaging techniques in terms of apparent current density and resistivity distribution using a standard numerical forward modelling and inversion technique. Efficacy of the developed approach was demonstrated for the interpretation of synthetic and field very low‐frequency data. The presented imaging technique shows improvement with respect to the filtering approaches in depicting subsurface conductors. Further, results obtained using the presented approach are closer to the results of rigorous resistivity inversion. Since the presented approach uses only the real anomaly, which is not sensitive to very small isolated near‐surface conducting features, it depicts prominent conducting features in the subsurface.     

重力场向下延拓Milne法

重力场向下延拓能够突出局部和浅部的异常信息,分离叠加的异常特征.但是向下延拓通常具有过程不稳定、下延深度小、结果不准确等问题.针对向下延拓所存在的不足,本文利用重力场及其垂向一阶导数,基于辛普森(Simpson)求积公式,推导出重力场向下延拓米尔尼(Milne)公式.将本文向下延拓方法应用于模型数据,向下延拓模型结果及误差曲线表明,相对于向下延拓快速傅里叶变换(FFT)法和积分迭代法,向下延拓Milne法的深度更大,相对误差更小;相对模型值,向下延拓Milne法能够获得稳定且准确的结果.对加拿大乃查科(Nechako)盆地地区实测航空重力数据进行本文方法向下延拓验证,处理结果表明,相对于实测异常,本文方法向下延拓结果能够很好还原实测数据,并且在进一步向下延拓中反映原始异常的趋势,增强局部和细小异常信息.