Generalized linear inversion of 2½‐dimensional gravity and magnetic anomalies

Author(s):  
S. F. Lai ◽  
W. J. Peeples ◽  
C. L. V. Aiken
2014 ◽  
Vol 65 (2) ◽  
pp. 163-174 ◽  
Author(s):  
Ilya Prutkin ◽  
Peter Vajda ◽  
Miroslav Bielik ◽  
Vladimír Bezák ◽  
Robert Tenzer

Abstract We present a new interpretation of the Kolárovo gravity and magnetic anomalies in the Danube Basin based on an inversion methodology that comprises the following numerical procedures: removal of regional trend, depth-wise separation of signal of sources, approximation of multiple sources by 3D line segments, non-linear inversion based on local corrections resulting in found sources specified as 3D star-convex homogenous bodies and/or 3D contrasting structural contact surfaces. This inversion methodology produces several admissible solutions from the viewpoint of potential field data. These solutions are then studied in terms of their feasibility taking into consideration all available tectono-geological information. By this inversion methodology we interpret here the Kolárovo gravity and magnetic anomalies jointly. Our inversion generates several admissible solutions in terms of the shape, size and location of a basic intrusion into the upper crust, or the shape and depth of the upper/lower crust interface, or an intrusion into the crystalline crust above a rise of the mafic lower crust. Our intrusive bodies lie at depths between 5 and 12 km. Our lower crust elevation rises to 12 km with and 8 km without the accompanying intrusion into the upper crust, respectively. Our solutions are in reasonable agreement with various previous interpretations of the Kolárovo anomaly, but yield a better and more realistic geometrical resolution for the source bodies. These admissible solutions are next discussed in the context of geological and tectonic considerations, mainly in relation to the fault systems.


Geophysics ◽  
1977 ◽  
Vol 42 (3) ◽  
pp. 610-622 ◽  
Author(s):  
Chao C. Ku

A computational method, which combines the Gaussian quadrature formula for numerical integration and a cubic spline for interpolation in evaluating the limits of integration, is employed to compute directly the gravity and magnetic anomalies caused by 2-dimensional and 3-dimensional bodies of arbitrary shape and arbitrary magnetic polarization. The mathematics involved in this method is indeed old and well known. Furthermore, the physical concept of the Gaussian quadrature integration leads us back to the old concept of equivalent point masses or equivalent magnetic point dipoles: namely, the gravity or magnetic anomaly due to a body can be evaluated simply by a number of equivalent points which are distributed in the “Gaussian way” within the body. As an illustration, explicit formulas are given for dikes and prisms using 2 × 2 and 2 × 2 × 2 point Gaussian quadrature formulas. The basic limitation in the equivalent‐point method is that the distance between the point of observation and the equivalent points must be larger than the distance between the equivalent points within the body. By using a reasonable number of equivalent points or dividing the body into a number of smaller subbodies, the method might provide a useful alternative for computing in gravity and magnetic methods. The use of a simplified cubic spline enables us to compute the gravity and magnetic anomalies due to bodies of arbitrary shape and arbitrary magnetic polarization with ease and a certain degree of accuracy. This method also appears to be quite attractive for terrain corrections in gravity and possibly in magnetic surveys.


Sign in / Sign up

Export Citation Format

Share Document