Simultaneous estimation of total magnetization direction and 3-D spatial orientation

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1365-1377 ◽  
Author(s):  
Walter E. Medeiros ◽  
João B.C. Silva

Magnetic interpretations are usually carried out either by assuming induced magnetization and estimating the model geometry, or by presuming a known source spatial orientation to estimate the total magnetization. We present a 3-D magnetic interpretation method that estimates simultaneously the total magnetization direction and the spatial orientation of the source. It is based on the approximation of the anomaly by the series derived from expanding the magnetic potential into multipoles and retaining source moments up to second order. The moments and linear combinations of moments appearing in the series are then inverted from the magnetic anomaly. The total magnetization is assumed constant in direction but not in magnitude. It is also presumed implicitly that the anomalous distribution of magnetization intensity has nonzero values in a finite‐volume region, is far from the observation points, and presents three othogonal planes of symmetry intersecting at the center of the dipole moment. The method is essentially linear and requires no a priori explicit assumption of a fixed geometry for the sources. The method is particularly suited to interpret compact, isolated or disjoint, but spatially correlated sources. If the source satisfies all assumptions presumed by the method, it is possible to obtain accurate, stable estimates of the total dipole moment vector, the position of the center of dipole moment, and the directions of all three principal axes of symmetry. If the source is not far from the observation plane and/or if the total magnetization direction is not constant, it is still possible to obtain accurate and stable estimates of the direction of the mean total magnetization and the projection, on the observation plane, of the center of dipole moment. The method is applied to magnetic data from the Gulf of Guinea Seamount. The estimated magnetic palaeopole is at 50°48′S and 74°54′E which is in good agreement with estimates published by other authors.

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.


Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1342-1353 ◽  
Author(s):  
Walter E. Medeiros ◽  
João B.C. Silva

We present a 3-D gravity interpretation method based on the inversion of source moments employing the series derived from the gravity anomaly expansion in multipoles and retaining moments up to second order only. It presents the advantages of being essentially linear and allowing a straightforward computer implementation. The method requires neither an explicit assumption about the source geometry nor a priori knowledge about the density contrast distribution, which may even be nonuniform. The method assumes implicitly: (1) that the source fits in the interior of an imaginary sphere whose center coincides with the source center of mass and whose radius is smaller than the depth to the source center of mass relative to the measuring plane, and (2) that the spatial density distribution presents three ortogonal planes of symmetry intersecting at the source center of mass. The first assumption can be met by upward continuing the observed anomaly. When both assumptions are met, the method produces reliable and stable estimates of the total anomalous mass, the coordinates of the center of mass, the three principal axes directions of the anomalous body, and the relative importance among the axes. The method is particularly suited for interpreting compact, isolated or disjoint, but spatially correlated sources. The method is sensitive to an incomplete removal of the regional field (presumably superimposed on the residual anomaly of interest) because the unremoved part of the regional field may strongly affect the estimates of the residual source moments.


2014 ◽  
Vol 644-650 ◽  
pp. 3459-3462 ◽  
Author(s):  
Lei Shi ◽  
Liang Hui Guo ◽  
Feng Yi Guo

Processing and interpretation of magnetic data usually require information of total magnetization direction. However, under the effects of remanent magnetization, total magnetization direction is different from induced magnetization direction, which makes data processing and interpretation complexity. In this paper, we present a new method by cross-correlation of magnetic dipole source for determination of magnetization direction from relatively isolated and approximate equiaxial-shape magnetic total field anomaly. This method calculates cross-correlation coefficient between observed magnetic total field anomaly and theoretical magnetic total field anomaly caused by a magnetic dipole source, by using a set of varying parameters of positions and total magnetization direction of dipole source for trial and error. The corresponding magnetization direction of maximum correlation coefficient is regarded as estimated total magnetization direction. Test on synthetic data indicates that this method reliably and effectively estimates the magnetization direction from relatively isolated and approximate equiaxial-shape magnetic total field anomaly.


2014 ◽  
Vol 644-650 ◽  
pp. 3793-3796
Author(s):  
Liang Hui Guo ◽  
Rui Gao ◽  
Guo Li Zhang

Under the effects of remanent magnetization, total magnetization direction is different from geomagnetic field direction, which makes magnetic data processing and interpretation complexity. In this paper, we present a new approach for estimating the total magnetization direction of sources via cross-correlation between the reduced-to-pole anomaly and the normalized source strength (who is less sensitive to remanent magnetization). The geomagnetic field direction is used to calculated the normalized source strength, while various assumed total magnetization directions are used to calculated the RTP anomalies. The maximum correlation between the RTP anomalies and the normalized corresponds to the estimated total magnetization direction. Test on synthetic data showed that the new approach is simple and effective.


Geophysics ◽  
2014 ◽  
Vol 79 (2) ◽  
pp. J11-J19 ◽  
Author(s):  
Shu-Ling Li ◽  
Yaoguo Li

We study the inversion of magnetic data acquired over a rugged observation surface and where the buried source bodies have strong remanent magnetization that leads to unknown total magnetization directions. These factors pose significant challenges for processing and inversion of such data. To tackle the challenges from both a rugged observation surface and an unknown magnetization direction, we propose a strategy through the joint use of the equivalent source technique and 3D amplitude inversion to obtain 3D magnetization strength. We use equivalent source processing to calculate the amplitude data in the space domain because the use of the wavenumber-domain method is invalid due to large variations in the data elevation. We then carried out an amplitude inversion to generate a 3D subsurface distribution of the magnitude of the total magnetization vector. The results from a synthetic example and aeromagnetic data in Daye Mine in China showed that this approach is effective and images the magnetic units whose contact zones with the limestone country rock host the mineralization. The method is general and can be applied to a variety of cases with similar challenges.


Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. L69-L73 ◽  
Author(s):  
Neal Dannemiller ◽  
Yaoguo Li

The characterization and interpretation of magnetic anomalies rely upon knowledge of the total magnetization direction. Magnetization is usually assumed to consist solely, or primarily, of induced magnetization. The presence of strong remanent magnetization can alter the direction significantly and consequently adversely affect the interpretation, leading to erroneous sizes or shapes of causative bodies. Therefore, it is imperative to have some understanding of the total magnetization direction. We propose a method based upon the correlation between two quantities in magnetic data interpretation: the vertical gradient and the total gradient of the reduced-to-pole (RTP) field. This method is tested on both synthetic and field data sets. The results show that the method is effective in a variety of situations, including those with two-dimensional and three-dimensional dipping bodies and a field example that has a large deviation between the inducing field direction and the total magnetization direction.


Geophysics ◽  
2009 ◽  
Vol 74 (3) ◽  
pp. L21-L30 ◽  
Author(s):  
Peter G. Lelièvre ◽  
Douglas W. Oldenburg

Inversion of magnetic data is complicated by the presence of remanent magnetization. To deal with this problem, we invert magnetic data for a three-component subsurface magnetization vector, as opposed to magnetic susceptibility (a scalar). The magnetization vector can be cast in a Cartesian or spherical framework. In the Cartesian formulation, the total magnetization is split into one component parallel and two components perpendicular to the earth’s field. In the spherical formulation, we invert for magnetization amplitude and the dip and azimuth of the magnetization direction. Our inversion schemes contain flexibility to obtain different types of magnetization models and allow for inclusion of geologic information regarding remanence. Allowing a vector magnetization increases the nonuniqueness of the magnetic inverse problem greatly, but additional information (e.g., knowledge of physical properties or geology) incorporated as constraints can improve the results dramatically. Commonly available information results in complicated nonlinear constraints in the Cartesian formulation. However, moving to a spherical formulation results in simple bound constraints at the expense of a now nonlinear objective function. We test our methods using synthetic and real data from scenarios involving complicated remanence (i.e., many magnetized bodies with many magnetization directions). All tests provide favorable results and our methods compare well against those of other authors.


Geophysics ◽  
2013 ◽  
Vol 78 (6) ◽  
pp. D429-D444 ◽  
Author(s):  
Shuang Liu ◽  
Xiangyun Hu ◽  
Tianyou Liu ◽  
Jie Feng ◽  
Wenli Gao ◽  
...  

Remanent magnetization and self-demagnetization change the magnitude and direction of the magnetization vector, which complicates the interpretation of magnetic data. To deal with this problem, we evaluated a method for inverting the distributions of 2D magnetization vector or effective susceptibility using 3C borehole magnetic data. The basis for this method is the fact that 2D magnitude magnetic anomalies are not sensitive to the magnetization direction. We calculated magnitude anomalies from the measured borehole magnetic data in a spatial domain. The vector distributions of magnetization were inverted methodically in two steps. The distributions of magnetization magnitude were initially solved based on magnitude magnetic anomalies using the preconditioned conjugate gradient method. The preconditioner determined by the distances between the cells and the borehole observation points greatly improved the quality of the magnetization magnitude imaging. With the calculated magnetization magnitude, the distributions of magnetization direction were computed by fitting the component anomalies secondly using the conjugate gradient method. The two-step approach made full use of the amplitude and phase anomalies of the borehole magnetic data. We studied the influence of remanence and demagnetization based on the recovered magnetization intensity and direction distributions. Finally, we tested our method using synthetic and real data from scenarios that involved high susceptibility and complicated remanence, and all tests returned favorable results.


Geophysics ◽  
2021 ◽  
pp. 1-47
Author(s):  
Rukuan Xie ◽  
Shengqing Xiong ◽  
Shuling Duan ◽  
Jinlong Wang ◽  
Ping Wang ◽  
...  

The total-field magnetic anomaly [Formula: see text] is an approximation of the projection [Formula: see text] of the magnetic anomaly vector [Formula: see text] onto the normal geomagnetic field [Formula: see text]. However, for highly magnetic sources, the approximation error of [Formula: see text] cannot be ignored. To reduce the error, we have developed a method for calculating [Formula: see text] by using airborne vector magnetic data based on the vector relationship of geomagnetic field [Formula: see text]. The calculation uses the magnitude of the vectors [Formula: see text], [Formula: see text], and [Formula: see text] through a simple approach. To ensure that each magnitude has the same level, we normalize the magnitude of [Formula: see text] using the total-field magnetic data measured by the scalar magnetic sensor. The method is applied to the measured airborne vector magnetic data at the Qixin area of the East Tianshan Mountains in China. The results indicate that the calculated [Formula: see text] has high precision and can distinguish the approximation error less than 3.5 nT. We also analyze the characteristics of the approximation error that are caused by the effects of different total magnetization inclinations. These error characteristics are used to predict the total magnetization inclination of a 2D magnetic source based on the measured airborne vector magnetic data.


Sign in / Sign up

Export Citation Format

Share Document