Interpretation of magnetic anomalies by horizontal and vertical derivatives of the analytic signal

Magnetic anomalies are often disturbed by the magnetization direction, so we can’t directly use the original magnetic anomaly to estimate the exact location and geometry of the source. The 2D analytic signal is insensitive to magnetization direction. In this paper, we present an automatic method based on the analytic signal horizontal and vertical derivatives to interpret the magnetic anomaly. We derive a linear equation using the analytic signal properties and we obtain the 2D magnetic body location parameters without giving a priori information. Then we compute the source structural index (expressing the geometry) by the estimated location parameters. The proposed method is demonstrated on synthetic magnetic anomalies with noise. For different models, the proposed technique can both successfully estimate the location parameters and the structural index of the sources and is insensitive to noise. Lastly, we apply it to real magnetic anomalies from China and obtain the distribution of unexploited iron ore. The inversion results are consistent with the parameters of known ore bodies.     

Direct analytic signal (DAS) method in the interpretation of magnetic data

We presented a new method for interpreting 2D magnetic data, called direct analytic signal (DAS) method, which directly used the analytic signal of magnetic anomaly to compute the depth and the structural index of the source. The DAS method needs only the computation of the first order derivatives of magnetic anomaly, so that the inversion results are more stable than the results obtained by the other existing analytic signal methods. The DAS method is tested on synthetic magnetic data with and without noise, and the DAS method can successfully obtain the depth and the structural index of the source. We also applied the DAS method to interpret a real magnetic data over a shallow geological source whose source parameters are known from closely drilling information, and the inversion results are in accord with the true values.     

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

通过分析解析信号概念,指出目前重磁场解析信号事实上是重磁场梯度解析信号. 在借助解析信号分量满足二维希尔伯特变换关系的基础上,提出重磁场直接解析信号的概念,并阐述重力异常及化极磁异常希尔伯特变换——直接解析信号的含义,并给出基于直接解析信号对位场增强的四种处理方法:直接解析信号模、水平分量模、改进Tilt angle和改进Theta map. 理论模型和实际航磁资料处理表明上述方法一定程度上可以解决噪声干扰等问题,对分辨和刻画场源有较好的解释效果,能够用于航空重力和航磁资料处理.     

应用加强解析信号倾斜角进行位场数据的边界检测

边界检测在地球物理位场数据解释中占有重要位置.现有的传统边界识别方法有的不能同时显示不同振幅的异常边界,有的虽然能均衡不同振幅的异常,但识别出来的边界信息中含有一些额外的错误的边界信息,尤其是当测量的异常中同时含有正异常和负异常时.目前已有的去除额外错误边界信息的方法存在着一定的人为主观性.为了解决这些问题,本文定义了加强解析信号倾斜角来进行地质体边界识别.通过模型试验证明了该方法不仅能同时清晰地识别深部和浅部地质体的边界,而且能有效地避免引入一些错误边界信息.最后将该方法应用到四川盆地的重力异常数据中,并取得了良好边界结果.     

Depth and structural index estimation of 2D magnetic source using correlation coefficient of analytic signal

We presented using the correlation coefficient of the analytic signal of real data and the analytic signal of synthetic data generated by the assumed source to estimate the structural index and the depth of the source. First, we assumed that the causative sources are located at different locations in the underground and the structural index of the assumed source is changed from 0 to 3, and then we separately compute the correlation coefficients of the analytic signal of the measured data and the analytic signal of the anomaly generated by each assumed source, the correlation coefficient can get the maximum value when the location and structural index of the assumed source are consistent with the real source. We tested the correlation coefficient method on synthetic noise-free and noise-corrupted magnetic anomalies, and the inversion results indicate that the new method can successfully finish the inversion of magnetic data. We also applied it to measured magnetic data, and we obtain the structural index and the location of the source.     

磁张量梯度数据方向解析信号比值的边界识别和空间位置反演方法

磁张量梯度测量具有高分辨率、多参量的优点,能更准确地描述磁源体的分布特征,在矿产资源勘探中具有广阔的用途.磁异常解析信号具有受倾斜磁化干扰小的特点,且为了增强深部地质体的分辨能力,本文提出磁张量梯度数据的解析信号比值的均衡边界识别及空间位置反演技术.磁张量梯度数据的均衡边界识别方法为不同方向解析信号比值的反正切函数,在降低倾斜磁化干扰的同时能有效地均衡不同深度地质体的响应,提高了对较深地质体的分辨率;空间位置反演技术是建立解析信号比值与地质体位置参数的对应方程,利用解析信号比值与地质体的对应关系作为约束条件来反演获得地质体的水平位置和深度信息,具有无需已知任何先验信息的优势.通过磁性体张量异常试验表明解析信号比值的边界识别方法能清晰和准确地获得不同深度地质体的边界,所建立的反演方程能准确地计算出地质体的范围和深度,具有较高的水平分辨率和精度.将本文方法应用于实测磁张量梯度数据的解释,获得了地下铁矿的分布特征,为区域矿产资源潜力评价提供了翔实的基础资料.

     

Interpretation of marine magnetic anomalies. Part I. A survey of existing methods and analysis of the analytic signal method

In this work, the problem of interpreting linear marine magnetic anomalies in the context of the spreading model is considered. Analysis and a comparison of the methods proposed earlier for solving this problem are carried out and the analytic signal method is considered in detail. In terms of a simplest example, it is shown that the AS method has significant restrictions associated primarily with the size and/or occurrence depth of the bodies.     

Euler deconvolution of the analytic signal and its application to magnetic interpretation

Euler deconvolution and the analytic signal are both used for semi‐automatic interpretation of magnetic data. They are used mostly to delineate contacts and obtain rapid source depth estimates. For Euler deconvolution, the quality of the depth estimation depends mainly on the choice of the proper structural index, which is a function of the geometry of the causative bodies. Euler deconvolution applies only to functions that are homogeneous. This is the case for the magnetic field due to contacts, thin dikes and poles. Fortunately, many complex geological structures can be approximated by these simple geometries. In practice, the Euler equation is also solved for a background regional field. For the analytic signal, the model used is generally a contact, although other models, such as a thin dike, can be considered. It can be shown that if a function is homogeneous, its analytic signal is also homogeneous. Deconvolution of the analytic signal is then equivalent to Euler deconvolution of the magnetic field with a background field. However, computation of the analytic signal effectively removes the background field from the data. Consequently, it is possible to solve for both the source location and structural index. Once these parameters are determined, the local dip and the susceptibility contrast can be determined from relationships between the analytic signal and the orthogonal gradients of the magnetic field. The major advantage of this technique is that it allows the automatic identification of the type of source. Implementation of this approach is demonstrated for recent high‐resolution survey data from an Archean granite‐greenstone terrane in northern Ontario, Canada.     

An analytic signal approach to the interpretation of total field magnetic anomalies1

The analytic signal (AS) is defined as the square root of the sum of the squares of the vertical and the two horizontal derivatives of the total magnetic field ΔT. This paper verifies theoretically that peaks of the AS correlate directly with their magnetic causative bodies and are positioned symmetrically over them, i.e. the main feature of the AS is that it is independent of the inclination of the magnetic field. This avoids the difficulties that are often faced in the conventional process of reduction to pole for ΔT, when the direction of magnetization of the causative bodies is not known. In addition, the AS has characteristics similar to the derivative features of the magnetic field, so that it is very sensitive to edge effects of the causative magnetic bodies. The theoretical derivations are tested by comparison with calculations on models, and, in a field example from Hunan Province, China, the AS is applied successfully to the interpretation of ΔT, whereas the conventional process of reduction to pole fails, due to the reverse magnetization of the causative body.     

Interpretation of magnetic anomalies using a genetic algorithm

A genetic algorithm (GA) is an artificial intelligence method used for optimization. We applied a GA to the inversion of magnetic anomalies over a thick dike. Inversion of nonlinear geophysical problems using a GA has advantages because it does not require model gradients or well-defined initial model parameters. The evolution process consists of selection, crossover, and mutation genetic operators that look for the best fit to the observed data and a solution consisting of plausible compact sources. The efficiency of a GA on both synthetic and real magnetic anomalies of dikes by estimating model parameters, such as depth to the top of the dike (H), the half-width of the dike (B), the distance from the origin to the reference point (D), the dip of the thick dike (δ), and the susceptibility contrast (k), has been shown. For the synthetic anomaly case, it has been considered for both noise-free and noisy magnetic data. In the real case, the vertical magnetic anomaly from the Pima copper mine in Arizona, USA, and the vertical magnetic anomaly in the Bayburt–Sar?han skarn zone in northeastern Turkey have been inverted and interpreted. We compared the estimated parameters with the results of conventional inversion methods used in previous studies. We can conclude that the GA method used in this study is a useful tool for evaluating magnetic anomalies for dike models.     

Interpretation of magnetic anomalies of dikes using correlation factors

The magnetic anomaly due to a buried dike consists of the sum of two easily separated elementary functions. These functions, which have simple symmetry, are called even and odd functions. The correlation factors (r _0,1 for the even andr _0,2 for the odd function) between least-squares residual anomalies from even and odd functions are computed. Correlation values are used to determine the depth to the top and the half-width of the dike. The method also includes the determination of the index parameter and the amplitude coefficient. The validity of the method is tested against a theoretical and a field example where the parameters of the latter were determined by other investigators in comparing the results.     

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

由于解析信号具有不受(二维) 或少受磁化方向影响,能够较好反映磁性体边界的特性, 因此受到人们的重视. 欧拉反褶积法可以确定场源的位置和深度以及形状因子,具有较强的适应性. 因此前人提出将二者相结合的方法.针对前人提出的方法中存在受高频干扰严重的问题,本文提出低阶的欧拉反褶积与解析信号相结合的位场反演方法.本方法在反演中只需计算磁异常的一、二阶导数,这样能够将高频噪声的干扰减少到人们可以控制的水平.文中详细推导了反演计算公式.模型计算证明了方法的正确性,同时探讨了噪声和邻近磁性体干扰影响的问题.对安徽庐江—枞阳火山岩盆地磁异常反演计算的结果证实了方法的实用性.本方法也可以应用于重力异常的反演计算中.     

瞬变电磁场资料的联合时-频分析解释

本文充分利用瞬变电磁（TEM）信号包含的信息,采用联合时-频分析方法解释了TEM资料。文章选取了Wigner-Ville分布和Gabor展开的联合时-频分布（或表示）,将一维时间域信号拓展到时-频二维平面上,不仅将不同地电断面的响应成功地分开来,而且所显示的信号特征与地层电性结构之间有着明确的物理意义,从时间和频率两个方面同时描述了瞬变涡流场在地层中的激发和衰减的过程。利用时间窗数据和频率窗数据所进行的联合时-频反演,缩小了拟合差,提高了定量解释精度。     

Interpretation of tensor gravity data using an adaptive tilt angle method

Full Tensor Gravity Gradiometry (FTG) data are routinely used in exploration programmes to evaluate and explore geological complexities hosting hydrocarbon and mineral resources. FTG data are typically used to map a host structure and locate target responses of interest using a myriad of imaging techniques. Identified anomalies of interest are then examined using 2D and 3D forward and inverse modelling methods for depth estimation. However, such methods tend to be time consuming and reliant on an independent constraint for clarification. This paper presents a semi‐automatic method to interpret FTG data using an adaptive tilt angle approach. The present method uses only the three vertical tensor components of the FTG data (T_zx, T_zy and T_zz) with a scale value that is related to the nature of the source (point anomaly or linear anomaly). With this adaptation, it is possible to estimate the location and depth of simple buried gravity sources such as point masses, line masses and vertical and horizontal thin sheets, provided that these sources exist in isolation and that the FTG data have been sufficiently filtered to minimize the influence of noise. Computation times are fast, producing plausible results of single solution depth estimates t hat relate directly to anomalies. For thick sheets, the method can resolve the thickness of these layers assuming the depth to the top is known from drilling or other independent geophysical data. We demonstrate the practical utility of the method using examples of FTG data acquired over the Vinton Salt Dome, Louisiana, USA and basalt flows in the Faeroe‐Shetland Basin, UK. A major benefit of the method is the ability to quickly construct depth maps. Such results are used to produce best estimate initial depth to source maps that can act as initial models for any detailed quantitative modelling exercises using 2D/3D forward/inverse modelling techniques.     

Fast inversion of magnetic data using Lanczos bidiagonalization method

This paper describes application of a fast inversion method to recover a 3D susceptibility model from magnetic anomalies. For this purpose, the survey area is divided into a large number of rectangular prisms in a mesh with unknown susceptibilities. Solving the full set of equations is substantially time consuming, and applying an algorithm to solve it approximately can reduce the time significantly. It is shown that the Lanczos bidiagonalization method can be an appropriate algorithm to solve a Tikhonov cost function for this purpose. Running time of the inverse modeling significantly decreases by replacing the forward operator matrix with a matrix of lower dimension. A weighted generalized cross validation method is implemented to choose an optimum value of a regularization parameter. To avoid the natural tendency of magnetic structures to concentrate at shallow depth, a depth weighting is applied. This study assumes that there is no remanent magnetization. The method is applied on a noise-corrupted synthetic data to demonstrate its suitability for 3D inversion. A case study including ground based measurement of magnetic anomalies over a porphyry-Cu deposit located in Kerman providence of Iran, Now Chun deposit, is provided to show the performance of the new algorithm on real data. 3D distribution of Cu concentration is used to evaluate the obtained results. The intermediate susceptibility values in the constructed model coincide with the known location of copper mineralization.     

基于航磁矢量数据的磁源总磁化方向估算

在地质体具有较强剩磁或自退磁效应的情况下,获得地质体总磁化方向对磁法勘探的数据处理、反演及解释具有重要意义.与传统航磁总场测量相比,航磁矢量(三分量)测量能够获得地磁场的矢量信息.基于三分量磁异常一阶矩与磁源磁矩的积分关系,我们实现了基于航磁矢量数据的磁源总磁化方向估算.针对该方法易受相邻异常影响,从而使多异常的估算难...     

CSAMT单分量数据解释方法

可控源音频大地电磁法（CSAMT）一直沿用大地电磁法（MT）的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性.

     

Interpretation of gravity and magnetic anomalies near Naples,Italy, using computer techniques

Some aspects of the analytical procedure for an automatic fitting of gravity and magnetic anomalies are discussed. An example of application relative to the Campania volcanic district near Naples, Italy, is reported. Gravity anomalies in that area mainly reflect the depth to the carbonate basement rocks. A small gravity high over Vesuvius may indicate buried volcanic layers. Magnetic anomalies near L. Patria and Castel Volturno are interpreted to be caused by buried intrusive rocks of higher susceptibility. Computer models of subsurface rock units with appropriate densities and susceptibilities have been generated to match the observed anomalies.     

Interpretation of gravity and magnetic data in the Phlegraean fields geothermal area, Naples, Italy

Available gravity and magnetic data of the Phlegraean Fields geothermal area, Naples, Italy, have been interpreted and the obtained structural models discussed in the light of the other available geological, volcanological and geophysical data.On the basis of the results of a previous seismic reflection survey in the Gulf of Naples and in the Pozzuoli Bay, which delineated a basement characterized by a seismic velocity of 4–6 km/s, it has been possible to evaluate the gravity anomaly connected with the morphology of this horizon ( = 2.7 g/cm³).The residual anomaly map, obtained after subtraction of the regional long-wavelength components relative to mantle and deep crustal structures and the computed components relative to the above-mentioned seismic basement, shows up as a circular low with an amplitude of 10 mgal centred in the Pozzuoli Bay. This gravity low has been interpreted as due to the occurrence, in the centre of Pozzuoli Bay, of light (Δ = −0.2 g/cm³) material with a maximum thickness of about 2 km. However, a contribution to the anomaly due to a narrow magmatic body intruded in the basement, as suggested by volcanological and ground deformation data, cannot be excluded.The aeromagnetic map of the Phlegraean Fields is characterized by three main anomalies which have been fitted by superficial tridimensional parallelepipedic bodies, schematically representing lava flows and domes. Their anomalies have been subsequently subtracted from the observed field, obtaining as a residual a large anomaly centred in the southwestern area of the Pozzuoli Bay. It has been interpreted as being due to a lowmagnetized body which, taking into account the thermal state of the area, should represent that part of the pyroclastic sequence which has lost part of its magnetization by thermo-chemical alteration.     

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

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