A derivative-based interpretation approach to estimating source parameters of simple 2D magnetic sources from Euler deconvolution, the analytic-signal method and analytical expressions of the anomalies

The major advantage of using either the analytic‐signal or the Euler‐deconvolution technique is that we can determine magnetic‐source locations and depths independently of the ambient earth magnetic parameters. In this study, we propose adopting a joint analysis of the analytic signal and Euler deconvolution to estimate the parameters of 2D magnetic sources. The results can avoid solution bias from an inappropriate magnetic datum level and can determine the horizontal locations, depths, structural types (indices), magnetization contrasts and/or structural dips. We have demonstrated the feasibility of the proposed method on 2D synthetic models, such as magnetic contacts (faults), thin dikes and cylinders. However, the method fails to solve the parameters of magnetic sources if there is severe interference between the anomalies of two adjacent magnetic sources.     

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.     

Automatic Interpretation of Magnetic Data Using Euler Deconvolution with Nonlinear Background

The voluminous gravity and magnetic data sets demand automatic interpretation techniques like Naudy, Euler and Werner deconvolution. Of these techniques, the Euler deconvolution has become a popular choice because the method assumes no particular geological model. However, the conventional approach to solving Euler equation requires tentative values of the structural index preventing it from being fully automatic and assumes a constant background that can be easily violated if the singular points are close to each other. We propose a possible solution to these problems by simultaneously estimating the source location, depth and structural index assuming nonlinear background. The Euler equation is solved in a nonlinear fashion using the optimization technique like conjugate gradient. This technique is applied to a published synthetic data set where the magnetic anomalies were modeled for a complex assemblage of simple magnetic bodies. The results for close by singular points are superior to those obtained by assuming linear background. We also applied the technique to a magnetic data set collected along the western continental margin of India. The results are in agreement with the regional magnetic interpretation and the bathymetric expressions.     

An Improved Analytic Signal Technique for the Depth and Structural Index from 2D Magnetic Anomaly Data

This paper presents a new inversion method for the interpretation of 2D magnetic anomaly data, which uses the combination of the analytic signal and its total gradient to estimate the depth and the nature (structural index) of an isolated magnetic source. However, our proposed method is sensitive to noise. In order to lower the effect of noise, we apply upward continuation technique to smooth the anomaly. Tests on synthetic noise-free and noise corrupted magnetic data show that the new method can successfully estimate the depth and the nature of the causative source. The practical application of the technique is applied to measured magnetic anomaly data from Jurh area, northeast China, and the inversion results are in agreement with the inversion results from Euler deconvolution of the analytic signal.     

基于垂向一阶导数与解析信号比值的欧拉反演方法

常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。     

Elimination of the water-layer response from multi-component source and receiver marine electromagnetic data

This paper presents the theory to eliminate from the recorded multi‐component source, multi‐component receiver marine electromagnetic measurements the effect of the physical source radiation pattern and the scattering response of the water‐layer. The multi‐component sources are assumed to be orthogonally aligned above the receivers at the seabottom. Other than the position of the sources, no source characteristics are required. The integral equation method, which for short is denoted by Lorentz water‐layer elimination, follows from Lorentz' reciprocity theorem. It requires information only of the electromagnetic parameters at the receiver level to decompose the electromagnetic measurements into upgoing and downgoing constituents. Lorentz water‐layer elimination replaces the water layer with a homogeneous half‐space with properties equal to those of the sea‐bed. The source is redatumed to the receiver depth. When the subsurface is arbitrary anisotropic but horizontally layered, the Lorentz water‐layer elimination scheme greatly simplifies and can be implemented as deterministic multi‐component source, multi‐component receiver multidimensional deconvolution of common source gathers. The Lorentz deconvolved data can be further decomposed into scattering responses that would be recorded from idealized transverse electric and transverse magnetic mode sources and receivers. This combined electromagnetic field decomposition on the source and receiver side gives data equivalent to data from a hypothetical survey with the water‐layer absent, with idealized single component transverse electric and transverse magnetic mode sources and idealized single component transverse electric and transverse magnetic mode receivers. When the subsurface is isotropic or transverse isotropic and horizontally layered, the Lorentz deconvolution decouples into pure transverse electric and transverse magnetic mode data processing problems, where a scalar field formulation of the multidimensional Lorentz deconvolution is sufficient. In this case single‐component source data are sufficient to eliminate the water‐layer effect. We demonstrate the Lorentz deconvolution by using numerically modeled data over a simple isotropic layered model illustrating controlled‐source electromagnetic hydrocarbon exploration. In shallow water there is a decrease in controlled‐source electromagnetic sensitivity to thin resistors at depth. The Lorentz deconvolution scheme is designed to overcome this effect by eliminating the water‐layer scattering, including the field's interaction with air.     

Inversion of potential field data using the structural index as weighting function rate decay

Nonparametric inverse methods provide a general framework for solving potential‐field problems. The use of weighted norms leads to a general regularization problem of Tikhonov form. We present an alternative procedure to estimate the source susceptibility distribution from potential field measurements exploiting inversion methods by means of a flexible depth‐weighting function in the Tikhonov formulation. Our approach improves the formulation proposed by Li and Oldenburg (1996, 1998) , differing significantly in the definition of the depth‐weighting function. In our formalism the depth weighting function is associated not to the field decay of a single block (which can be representative of just a part of the source) but to the field decay of the whole source, thus implying that the data inversion is independent on the cell shape. So, in our procedure, the depth‐weighting function is not given with a fixed exponent but with the structural index N of the source as the exponent. Differently than previous methods, our choice gives a substantial objectivity to the form of the depth‐weighting function and to the consequent solutions. The allowed values for the exponent of the depth‐weighting function depend on the range of N for sources: 0 ≤N≤ 3 (magnetic case). The analysis regarding the cases of simple sources such as dipoles, dipole lines, dykes or contacts, validate our hypothesis. The study of a complex synthetic case also proves that the depth‐weighting decay cannot be necessarily assumed as equal to 3. Moreover it should not be kept constant for multi‐source models but should instead depend on the structural indices of the different sources. In this way we are able to successfully invert the magnetic data of the Vulture area, Southern Italy. An original aspect of the proposed inversion scheme is that it brings an explicit link between two widely used types of interpretation methods, namely those assuming homogeneous fields, such as Euler deconvolution or depth from extreme points transformation and the inversion under the Tikhonov‐form including a depth‐weighting function. The availability of further constraints, from drillings or known geology, will definitely improve the quality of the solution.     

欧拉反褶积与解析信号相结合的位场反演方法

由于解析信号具有不受(二维)或少受磁化方向影响,能够较好反映磁性体边界的特性,因此受到人们的重视.欧拉反褶积法可以确定场源的位置和深度以及形状因子,具有较强的适应性.因此前人提出将二者相结合的方法.针对前人提出的方法中存在受高频干扰严重的问题,本文提出低阶的欧拉反褶积与解析信号相结合的位场反演方法.本方法在反演中只需计...     

基于最小反演拟合差的重磁场源深度计算方法

以等效源及位场物性反演为基础,本文提出一种新的求取重磁场源深度的方法.该方法将一层等效源以一定的间隔从浅部向深部移动,并将等效源作为初始模型进行反演,当反演拟合差最小时,停止反演,此时的等效源底深即为所求场源的中心深度.由于仅需要反演一层等效源,比传统的物性反演计算时间大大减少,并且不需要进行深度加权约束.理论模型数据处理结果表明该方法能够获得较准确的场源深度:以长宽比为7.5的薄板模型为例,深度计算误差约为1个点距(25 m);以长宽比为0.5~1.5的厚板模型为例,深度计算误差小于1个点距(25m).将该方法应用于实测航磁梯度数据,计算的磁源中心深度在200~250m之间,钻井资料显示该异常由埋藏深度在200~300m的闪长岩引起,计算结果与钻井资料较吻合.     

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.     

TSVD analysis of Euler deconvolution to improve estimating magnetic source parameters: An example from the Åsele area,Sweden

In this paper, I introduce a new approach based on truncated singular value decomposition (TSVD) analysis for improving implementation of grid-based Euler deconvolution with constraints of quasi 2D magnetic sources. I will show that by using TSVD analysis of the gradient matrix of magnetic field anomaly (reduced to pole) for data points located within a square window centered at the maximum of the analytic signal amplitude, we are able to estimate the strike direction and dip angle of 2D structures from the acquired eigenvectors. It is also shown that implementation of the standard grid-based Euler deconvolution can be considerably improved by solving the Euler's homogeneity equation for source location and structural index, simultaneously, using the TSVD method. The dimensionality of the magnetic anomalies can be indicated from the ratio between the smallest and intermediate eigenvalues acquired from the TSVD analysis of the gradient matrix. For 2D magnetic sources, the uncertainty of the estimated source location and structural index is significantly reduced by truncating the smallest eigenvalue.Application of the method is demonstrated on an aeromagnetic data set from the Åsele area in Sweden. The geology of this area is dominated by several dike swarms. For these dolerite dikes, the introduced method has provided useful information of strike directions and dip angles in addition to the estimated source location and structural index.     

An enhanced method for source parameter imaging of magnetic data collected for mineral exploration

We have developed a method for imaging magnetic data collected for mineral exploration to yield the following structural information: depth, model type (structural index) and susceptibility. The active nature of mineral exploration data requires we derive the structural information from a robust quantity: we propose that the first‐ or second‐order analytic‐signal amplitude is suitably stable. The procedure is to normalize the analytic‐signal amplitude by the peak value and then use non‐linear inversion to estimate the depth and the structural index for each anomaly. In our field example, different results are obtained depending on whether we inverted for the first‐ or second‐order analytic‐signal amplitude. This is probably because the two‐dimensional contact, thin sheet or horizontal cylinder models we have assumed are not appropriate. In cases such as these, when our model assumptions are not correct, the results should not be interpreted quantitatively, but they might be useful for giving a qualitative indication of how the structure might vary. With a priori information, it is possible to assume a model type (i.e. set the structural index) and generate estimates of the depth and susceptibility. These data can then be gridded and imaged. If a contact is assumed, the susceptibility contrast is estimated; for the dike model, the susceptibility‐thickness is estimated; for the horizontal cylinder, the susceptibility‐area is estimated. To emphasize that the results are dependent on our assumed model, we advocate prefixing any derived quantity by the term ‘apparent’.     

The use of potential field data in revealing the basement structure in sub-basaltic settings: an example from the Møre margin, offshore Norway

This study investigates the utility of the potential fields (gravity and magnetics) in volcanic settings as observed on the Møre margin. Synthetic models are used to investigate the effect of volcanics on the gravity and magnetic fields. The focus is on detecting sub-basaltic basement structures. The methods applied to the models are Euler deconvolution on magnetic data, gravity gradients and integrated 3D gravity and magnetic forward modelling. The same methods are used on the Møre margin and the results compared to the synthetic models. The Euler deconvolution on the magnetic signal does provide limited depth solutions in the volcanic environment and the use of different observation levels does not enhance the results. Forward gravity and magnetic models provide a valuable tool to estimate both the basalt and sub-basaltic sedimentary thickness but are limited by the ambiguity inherent in potential field methods. The use of gravity gradients significantly decreases the available model solutions and provides boundary detection even in sub-basaltic settings.     

AUTOMATIC INTERPRETATION OF GRAVITY GRADIOMETRIC DATA IN TWO DIMENSIONS: VERTICAL GRADIENT1

The magnetic and gravity field produced by a given homogeneous source are related through Poisson's equation. Starting from this consideration, it is shown that some 2D interpretation tools, widely applied in the analysis of aeromagnetic data, can also be used for the interpretation of gravity gradiometric data (vertical gradient). This paper deals specifically with the Werner deconvolution, analytic signal and Euler's equation methods. After a short outline of the mathematical development, synthesized examples have been used to discuss the efficiency and limits of these interpretation methods. These tools could be applied directly to airborne gravity gradiometric data as well as ground gravity surveys after transformation of the Bouguer anomalies into vertical gradient anomalies. An example is given of the application of the Werner deconvolution and Euler's equation methods to a microgravity survey.     

The Automatic Determination of the Location and Depth of Contacts and Dykes from Aeromagnetic Data

This paper demonstrates how to automatically obtain the location and depth of magnetic field sources of known structural index, such as contacts and dykes, for profile and map datasets. The method uses the analytic signal amplitude, and does not require that the magnetic data be pole-reduced or have a known magnetisation vector. The method is applied both to synthetic data and to aeromagnetic data from South Africa.     

Can Euler deconvolution outline three-dimensional magnetic sources?

Severe limitations of the standard Euler deconvolution to outline source shapes have been pointed out. However, Euler deconvolution has been widely employed on field data to outline interfaces, as faults and thrust zones. We investigate the limitations of the 3D Euler deconvolution–derived estimates of source dip and volume with the use of reduced-to-the-pole synthetic and field anomalies. The synthetic anomalies are generated by two types of source bodies: (1) uniformly magnetized prisms, presenting either smooth or rough interfaces, and (2) bodies presenting smooth delimiting interfaces but strong internal variation of magnetization intensity. The dip of the first type of body might be estimated from the Euler deconvolution solution cluster if the ratio between the depth to the top and vertical extent is relatively high (>1/4). For the second type of body, besides dip, the source volume can be approximately delimited from the solution cluster envelope, regardless of the referred ratio. We apply Euler deconvolution to two field anomalies which are caused by a curved-shape thrust zone and by a banded iron formation. These anomalies are chosen because they share characteristics with the two types of synthetic bodies. For the thrust zone, the obtained Euler deconvolution solutions show spatial distribution allowing to estimate a source dip that is consistent with the surface geology data, even if the above-mentioned ratio is much less than 1/4. Thus, there are other factors, such as a heterogeneous magnetization, which might be controlling the vertical spreading of the Euler deconvolution solutions in the thrust zone. On the other hand, for the iron-ore formation, the solution cluster spreads out occupying a volume, in accordance with the results obtained with the synthetic sources having internal variation of magnetization intensity. As conclusion, although Euler deconvolution–derived solutions cannot offer accurate estimates of source shapes, they might provide a sufficient degree of reliability in the initial estimates of the source dip and volume, which may be useful in a later phase of more accurate modelling.     

用Euler反褶积方法反演台湾海峡磁异常

用Ｅｕｌｅｒ反褶积方法反演台湾海峡磁异常，选取的窗口不能大，而决定取舍反褶积解的误差限又不能小，这势必对解的质量有较大影响．为此，不仅用了不同延深的单体模型，而且进行了多体模型实验．结果表明，复杂分布的磁性体，用Ｅｕｌｅｒ反褶积方法确定磁性体的轮廓可能是困难的，但可确定磁性体的水平位置和深度，从而降低了对资料精度和计算参数选取的要求．在磁异常比较复杂的地区．     

Imaging magnetic sources using Euler's equation

The conventional Euler deconvolution method has the advantage of being independent of magnetization parameters in locating magnetic sources and estimating their corresponding depths. However, this method has the disadvantage that a suitable structural index must be chosen, which may cause spatial diffusion of the Euler solutions and bias in the estimation of depths to the magnetic sources. This problem becomes more serious when interfering anomalies exist. The interpretation of the Euler depth solutions is effectively related to the model adopted, and different models may have different structural indices. Therefore, I suggest a combined inversion for the structural index and the source location from the Euler deconvolution, by using only the derivatives of the magnetic anomalies. This approach considerably reduces the diffusion problem of the location and depth solutions. Consequently, by averaging the clustered solutions satisfying a given criterion for the solutions, we can image the depths and attributes (or types) of the causative magnetic sources. Magnetic anomalies acquired offshore northern Taiwan are used to test the applicability of the proposed method.     

Avoidable Euler Errors – the use and abuse of Euler deconvolution applied to potential fields

Window‐based Euler deconvolution is commonly applied to magnetic and sometimes to gravity interpretation problems. For the deconvolution to be geologically meaningful, care must be taken to choose parameters properly. The following proposed process design rules are based partly on mathematical analysis and partly on experience.

相似文献   

不同高度数据联合欧拉反褶积法研究

基于不同测量高度重力场及其梯度数据可同时对应同一场源并用于反演场源位置的分析原理,拓展不同高度场数据在欧拉反褶积法中的应用范围.首先,立足于对欧拉反褶积方法的理论研究基础,提出不同高度数据融合联合欧拉反演公式.其次,在理论模型上对多种高度数据联合反演做了测试分析计算,验证了不同高度场数据融合联合欧拉反褶积法能够改善位场解释中单一观测面数据计算带来的解的发散问题,收敛过程由此改善.最后,将本文方法应用于龙门山地区实际重力数据的解释,获得了研究区断裂分布特征.