A new potential field shape descriptor using continuous wavelet transforms

Geophysics ◽  
2020 ◽  
Vol 85 (5) ◽  
pp. G81-G92
Author(s):  
P. Cavalier ◽  
D. W. O’Hagan

Potential field characterization aims at determining source depths, inclination, and type, preferably without a priori information. For ideal sources, the type is often defined from the field’s degree of homogeneity, derived from its expression in the space domain. We have developed a new shape descriptor for potential field source functions, stemming from spectral-domain parameters, which manifest clearly when using continuous wavelet transforms (CWTs). We generalize the use of the maximum wavelet coefficient points in the CWT diagram for the analysis of all types of potential fields (gravity, magnetic, and self-potential). We interpret the CWT diagram as a similarity diagram between the wavelet and the analyzed signal, which has fewer limitations than its interpretation as a weighted and upward-continued field projection. We develop new formulas for magnetic source depth prediction, as well as for effective inclination estimation, using various kinds of wavelets. We found that the potential field source functions exhibit precise behaviors in the CWT analysis that can be predicted using a single parameter [Formula: see text], which is related to their Fourier transforms. This parameter being scale and rotation-invariant can be used as a source-body shape descriptor similar to the commonly used structural index (SI). An advantage of the new descriptor is an increased level of discrimination between sources because it takes different values to describe the horizontal or the vertical cylinder structures. Our approach is illustrated on synthetic examples and real data. The method can be applied directly with the native form of the CWT without scaling factor modification, negative plane diagram extension, or downward plotting. This framework offers an alternative to existing wavelet-like projection methods or other classic deconvolution techniques relying on SI for determining the source depth, dip, and type without a priori information, with an increased level of differentiation between source structures thanks to the new shape descriptor.

Geophysics ◽  
1997 ◽  
Vol 62 (3) ◽  
pp. 814-830 ◽  
Author(s):  
Maurizio Fedi

The depth to the top, or bottom, and the density of a 3-D homogeneous source can be estimated from its gravity or magnetic anomalies by using a priori information on the maximum and minimum source depths. For the magnetic case, the magnetization direction is assumed to be constant and known. The source is assumed to be within a layer of known depth to the top h and thickness t. A depth model, satisfying both the data and the a priori information is found, together with its associated density/magnetization contrast. The methodology first derives, from the measured data, a set of apparent densities [Formula: see text] (or magnetizations), which do not depend on the layer parameters h and t, but only on source thickness. A nonlinear system of equations based on [Formula: see text], with source thicknesses as unknowns, is constructed. To simplify the solution, a more practical system of equations is formed. Each equation depends on only one value of thickness. Solving for the thicknesses, taking into account the above a priori information, the source depth to the top (or to the bottom) is determined uniquely. Finally, the depth solutions allow a unit‐density gravity model to be computed, which is compared to the observed gravity to determine the density contrast. A similar procedure can be used for magnetic data. Tests on synthetic anomalies and on real data demonstrate the good performance of this method.


Sign in / Sign up

Export Citation Format

Share Document