Improved Local Wavenumber Methods in the Interpretation of Potential Field Data

We present two new potential-inversion methods for estimating the depth and the nature (structural index) of the source, which use various combinations of different forms of local wavenumbers and the information about the horizontal location to estimate individually the depth and the nature of a magnetic source. The improved local wavenumber methods only use the horizontal offset and vertical offset of local wavenumbers to estimate the depth and the structural index of the source, so they yield more stable results compared with the results obtained by current methods that require the derivatives of local wavenumbers. Tests conducted with synthetic noise-free and noise-corrupted magnetic data show that the proposed methods can successfully estimate the depth and the nature of the geologic body. However, our methods are sensitive to high-wavenumber noise present in the data, and we reduced the noise effect by upward continuing the noise-corrupted magnetic data. The practical application of the new methods is tested on a real magnetic anomaly over a dike whose source parameters are known and the inversion results are consistent with the true values.     

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.     

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.     

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.     

Depth and structural index from normalized local wavenumber of 2D magnetic anomalies

Recent improvements in the local wavenumber approach have made it possible to estimate both the depth and model type of buried bodies from magnetic data. However, these improvements require calculation of third‐order derivatives of the magnetic field, which greatly enhances noise. As a result, the improvements are restricted to data of high quality. We present an alternative method to estimate both the depth and model type using the first‐order local wavenumber approach without the need for third‐order derivatives of the field. Our method is based on normalization of the first‐order local wavenumber anomalies and provides a generalized equation to estimate the depth of some 2D magnetic sources regardless of the source structure. Information about the nature of the sources is obtained after the source location has been estimated. The method was tested using synthetic magnetic anomaly data with random noise and using three field examples.     

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.     

一维波动方程小波逐版本反演

通过将波动方程中的密度及弹性模量参数（函数）投影到一系列波数有上限的Ｌ_２（Ｒ）的子空间Ｖｉ上，并在该子空间（而不是Ｌ_２（Ｒ））中寻找尽量满足条件的密度和弹性模量函数，令Ｊ→－∞．并充分利用小波技术的优越性，相当于对反演指标实施了一系列低通滤波，对解决传统反演问题中局部极值问题是很有希望的途径．逐版本反演不仅可用于波动方程，同样可用于解其他类型的微分方程反问题．     

解释位场全张量数据的张量局部波数法及其与常规局部波数法的比较

全张量探测技术以其信息量大、精度高、干扰小等优点在地球物理领域中得到广泛应用.本文提出采用张量局部波数法来进行位场全张量数据的解释,首先给出了张量局部波数的定义,然后推导出利用张量局部波数法进行反演的基本公式.本文方法在进行张量数据反演时无需事先知道场源体的类型(构造指数)即可获得场源体的位置信息,且可根据位置参数对场源体的类型进行估计.通过理论模型证明张量局部波数法可以很好地完成位场全张量数据的反演工作,并将其与常规局部波数法进行对比,证明全张量局部波数法的反演结果更加准确,即使在测点分布不合理的情况下,张量局部波数法仍可以获得准确的结果.最后应用张量局部波数法对美国得克萨斯州实测重力数据进行了反演,其反演结果与已有的研究成果相一致.     

基于改进的接收点插值算法的频率域海洋可控源电磁法2.5维正演

在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.     

A new concept in Euler deconvolution of isolated gravity anomalies

Euler's homogeneity equation has been used to develop a new technique to interpret the gravity anomalies over some simple geometrical sources, namely a finite horizontal line/vertical line, a finite vertical ribbon, a semicircular dome/basin and an isosceles triangle approximating an anticline/syncline. A linear over-determined system of equations has been solved to compute the depth, the horizontal location and the structural index, all treated as free parameters. The concept of a variable structural index provides better depth estimates and helps to identify the source geometry. Nomograms have been prepared to compute an additional model parameter, namely the horizontal/vertical extent of a line, the vertical extent of a ribbon and the radius of a dome/basin. The efficacy of the proposed method has been evaluated using two real field examples.     

Accuracy‐constrained optimisation methods for staggered‐grid elastic wave modelling

The classical finite‐difference methods for seismic wave modelling are very accurate at low wavenumbers but suffer from inaccuracies at high wavenumbers, particularly at Nyquist wavenumber. In contrast, the optimisation finite‐difference methods reduce inaccuracies at high wavenumbers but suffer from inaccuracies at low wavenumbers, particularly at zero wavenumber when the operator length is not long and the whole range of wavenumbers is considered. Inaccuracy at zero wavenumber means that the optimisation methods only have a zeroth‐order accuracy of truncation and thus are not rigorously convergent. To guarantee the rigorous convergence of the optimisation methods, we have developed accuracy‐constrained optimisation methods. Different‐order accuracy‐constrained optimisation methods are presented. These methods not only guarantee the rigorous convergence but also reduce inaccuracies at low wavenumbers. Accuracy‐constrained optimisation methods are applied to staggered‐grid elastic wave modelling.     

Wavenumber Sampling Issues in 2.5D Frequency Domain Seismic Modelling

There are several important wavenumber sampling issues associated with 2.5D seismic modelling in the frequency domain, which need careful attention if accurate results are to be obtained. At certain critical wavenumbers there exist rapid disruptions in the mainly smooth oscillatory spectra. The amplitudes of these disruptions can be very large, and this affects the accuracy of the inverse Fourier transformed frequency-space domain solution. In anisotropic elastic media there are critical wavenumbers associated with each wave mode—the quasi-P (qP) wave, and the two quasi-shear (qS1 and qS2) waves. A small wavenumber sampling interval is desirable in order to capture the highly oscillatory nature of the wavenumber spectrum, especially at increasing distance from the source. Obviously a small wavenumber sampling interval adds greatly to the computational effort because a 2D problem must be solved for every wavenumber and every frequency. The discretisation should be carried out up to some maximum wavenumber, beyond which the field becomes evanescent (exponentially decaying or diffusive). For receivers close to the source, activity persists beyond the critical wavenumber associated with the minimum shear wave velocity in the model. Fortunately, for receivers well removed from the source, the contribution from the evanescent energy is negligible and so there is no need to sample beyond this critical wavenumber. Sampling at Gauss–Legendre spacings is a satisfactory approach for acoustic media, but it is not practical in elastic media due to the difficulty of partitioning the integration around the different critical wavenumbers. We found to our surprise that in transversely isotropic media, the critical wavenumbers are independent of wave direction, but always occur at those wavenumbers corresponding to the maximum phase velocities of the three wave modes (qP, qS1 and qS2), which depend only on the elastic constants and the density. Additionally, we have observed that intermediate layers between source and receiver can filter out to a large degree, the sharp irregularities around the critical wavenumbers in the ω–k _y spectra. We have found that, using the spectral element method, the singularities (poles) at the critical wavenumbers which exist with analytic solutions, do not arise. However, the troublesome spike-like behaviour still occurs and can be damped out without distorting the spectrum elsewhere, through the introduction of slight attenuation.     

低纬度磁异常化极方法应用效果对比

低纬度磁异常化极的方法很多,包括在空间域化极和在频率域化极,比较起来,频率域方法计算简单快速,目前化极主要在频率域进行.然而大部分低纬度化极方法只压制了噪声并没有针对化极因子做改动,也就没有真正达到压制南北向条带状拉长的目的.本文针对其中4种在频率域对化极因子做处理的方法(双曲正弦法,压制因子法,直接阻尼法,伪倾角法),设计模型进行对比,并用南海某地区的实测数据进行4种方法的化极效果对比,选出一种最符合南海地区实际情况的化极方法.并对原有的直接阻尼法和伪倾角法做了适当的改进.     

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.     

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

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

A stable iterative downward continuation of potential field data

Downward continuation is a useful tool in the processing of potential field data, which can effectively enhance weak anomalies and identify overlap anomalies, but we all know that the computation of downward continuation is unstable, and easily distorts the true feature of potential field data. Because the computation of upward continuation and horizontal derivatives is stable, we proposed using the combination of upward continuation and horizontal derivative to accomplish the downward continuation of potential field data. The proposed method is demonstrated on synthetic potential field data, and the results show that the proposed method can finish the downward continuation of the data stably and precisely, and the precision of the proposed method is higher than the traditional method. We also apply it to real potential field data, and the results show that the proposed method accomplishes the downward continuation of the real data stably.     

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.     

重力异常视深度滤波及应用

为了解观测面以下位场场源体产生的异常,设计了一种波数域的滤波器,该滤波器可在无需岩石物性资料的前提下用于位场数据的反演．对位场资料用解析法和差分法上、下延拓,以消除浅表干扰源的影响,并最后下延至观测面以下,根据下延后场的极值确定场源的质心深度．理论模型的试算和实际资料的处理均获得了令人满意的效果．