Geology from geophysics; or, converting magnetic data to depth, dip, and susceptibility contrast using the source parameter imaging (SPITM) method

1996 ◽  
Author(s):  
Jeffrey B. Thurston ◽  
Richard S. Smith
Geophysics ◽  
1997 ◽  
Vol 62 (3) ◽  
pp. 807-813 ◽  
Author(s):  
Jeffrey B. Thurston ◽  
Richard S. Smith

The Source Parameter Imaging (SPI™) method computes source parameters from gridded magnetic data. The method assumes either a 2-D sloping contact or a 2-D dipping thin‐sheet model and is based on the complex analytic signal. Solution grids show the edge locations, depths, dips, and susceptibility contrasts. The estimate of the depth is independent of the magnetic inclination, declination, dip, strike and any remanent magnetization; however, the dip and the susceptibility estimates do assume that there is no remanent magnetization. Image processing of the source‐parameter grids enhances detail and provides maps that facilitate interpretation by nonspecialists. The SPI method tests successfully on synthetic profile and gridded data. SPI maps derived from aeromagnetic data acquired over the Peace River Arch area of northwestern Canada correlate well with known basement structure and furthermore show that the Ksituan Magmatic Arc can be divided into several susceptibility subdomains.


Geophysics ◽  
1984 ◽  
Vol 49 (8) ◽  
pp. 1354-1360 ◽  
Author(s):  
A. Guillen ◽  
V. Menichetti

The nonuniqueness of gravity or magnetic data inversion is well known. In order to remove ambiguity, some authors have sought solutions minimizing a functional describing geometrical or physical properties. Last and Kubik (1983), in particular, developed a method explaining the observed anomaly by structures of minimum volume. In this method the domain where anomalous sources are searched is divided into elementary prisms of a constant density or susceptibility contrast. Each elementary contrast is allowed to vary individually. Thus a contrast distribution is computed. The search for this kind of solution leads in general to geologically more appropriate bodies, but exceptions do occur. In this paper, the technique is broadened to include the search for solutions minimizing the moment of inertia with respect to the center of gravity or with respect to a given dip line passing through it. The resulting structures are both deeper and more compact, precisely as is required in specific cases. Theoretical and actual examples illustrate this flexible inversion technique.


2020 ◽  
Vol 10 (2) ◽  
pp. 171
Author(s):  
Sehah Sehah ◽  
Sukmaji Anom Raharjo ◽  
Urip Nurwijayanto Prabowo

Two dimensional modeling to basaltic rocks intrusion in Pekuncen and Karanglewas Villages Jatilawang District, Banyumas Regency, Central Java based on the local magnetic anomalies data has been carried out in March – June 2020. The amount of magnetic data obtained from the acquisition in the field was 239 data stretching in position of 109.107222° – 109.134944°E and 7.561361° – 7.577306°S, with the local magnetic anomalies values ranging of -2,961.11 – 1,516.31 nT. To model anomalous sources in the subsurface in two dimensions, then the local magnetic anomalies data is transformed into pseudogravity anomalies data, so that anomalous value can be obtained as -27.815 – 41.087 mGal. Based on the pseudogravity anomalous map, the basaltic rock intrusion is interpreted to be located in the eastern part of the research area, so modeling of anomalous sources is conducted in this area. The results of 2D-modeling to local magnetic anomalies data indicate the presence of anomalous object interpreted as basaltic rock intrusion with magnetic susceptibility contrast value of 0.0223 cgs, located at depth of 52.61 – 505.97 m and a lateral length of 1777.94 m. This rock intrudes sediment rock from the Halang Formation and is connected to other basaltic rock near the surface with magnetic susceptibility contrast value of 0.0165 cgs, located at depth of 1.94 – 80.90 m and lateral length of 751.83 m. The results of lithological interpretation are in accordance with the geological information of the research area.


1989 ◽  
Vol 20 (2) ◽  
pp. 25 ◽  
Author(s):  
P.M. Smith ◽  
M. Whitehead

The presence of a large anomalous structure in the northern part of Permit AC/P2 in the Timor Sea has been recognised ever since seismic data were first acquired in the area. Historically, however, sparse seismic coverage has always prevented a detailed and unambiguous interpretation of the complicated structure. In order to overcome this problem, some 2000 km of 3D seismic data were acquired over the feature. In conjunction with this seismic survey, detailed gravity and magnetic data sets were also recorded over the structure.Interpretation of the new seismic data indicated the presence of a piercement structure which is associated with a small negative Bouguer gravity anomaly and a magnetic intensity anomaly resulting from a positive susceptibility contrast. Modelling of the magnetic data indicated that an acidic or intermediate intrusive body was the most likely cause of the piercement structure. The presence of an acidic intrusive body was consistent with the gravity data which indicated that no large density contrast existed between the material of the piercement structure and the surrounding sediments.The combined interpretation of these three data sets was tested by a well, Paqualin-1, drilled on the flank of the piercement structure. The well encountered a thick evaporite sequence with associated thin bands of magnetitie and intermediate igneous rocks. It was logged with a three component downhole magnetic probe and forward magentic modelling work incorporating the results of the magnetic log gave good agreement with the observed aeromagnetic profiles.


Geophysics ◽  
2000 ◽  
Vol 65 (3) ◽  
pp. 849-860 ◽  
Author(s):  
Jörg Herwanger ◽  
Hansruedi Maurer ◽  
Alan G. Green ◽  
Jürg Leckebusch

A vertical‐gradient magnetic system based on optically pumped Cesium sensors has been used to map subtle magnetic anomalies across infilled pit houses and ditches at a medieval archeological site in northern Switzerland. For estimating the locations and dimensions of these features from the recorded data, we have designed and implemented an appropriate inversion scheme. Tests of this scheme on realistic synthetic data sets suggested that suitable minimum magnetic susceptibility contrasts and smoothing parameters for the inversion may be directly extracted from the data. Inversions with minimum magnetic susceptibility contrasts generated causative bodies with maximum plausible sizes. By using higher magnetic susceptibility contrasts, a complete suite of models that matched the data equally well was produced. To constrain better the magnetic susceptibility constrast within a selected area of the archeological site, shallow samples of topsoil and sediment were analyzed in the laboratory. An inversion based on the measured magnetic susceptibility contrast yielded reliable estimates of the locations, 3-D geometries, and sizes of two small pit houses. The depth extent of one pit house was subsequently verified by shallow drilling. We concluded that inversions of vertical‐gradient magnetic data constrained by magnetic susceptibility or shallow borehole information are rapid and inexpensive means of providing key knowledge on the depth distribution of inductively magnetized bodies.


2020 ◽  
Vol 1 (3) ◽  
Author(s):  
Maysam Abedi

The presented work examines application of an Augmented Iteratively Re-weighted and Refined Least Squares method (AIRRLS) to construct a 3D magnetic susceptibility property from potential field magnetic anomalies. This algorithm replaces an lp minimization problem by a sequence of weighted linear systems in which the retrieved magnetic susceptibility model is successively converged to an optimum solution, while the regularization parameter is the stopping iteration numbers. To avoid the natural tendency of causative magnetic sources to concentrate at shallow depth, a prior depth weighting function is incorporated in the original formulation of the objective function. The speed of lp minimization problem is increased by inserting a pre-conditioner conjugate gradient method (PCCG) to solve the central system of equation in cases of large scale magnetic field data. It is assumed that there is no remanent magnetization since this study focuses on inversion of a geological structure with low magnetic susceptibility property. The method is applied on a multi-source noise-corrupted synthetic magnetic field data to demonstrate its suitability for 3D inversion, and then is applied to a real data pertaining to a geologically plausible porphyry copper unit.  The real case study located in  Semnan province of  Iran  consists  of  an arc-shaped  porphyry  andesite  covered  by  sedimentary  units  which  may  have  potential  of  mineral  occurrences, especially  porphyry copper. It is demonstrated that such structure extends down at depth, and consequently exploratory drilling is highly recommended for acquiring more pieces of information about its potential for ore-bearing mineralization.


Sign in / Sign up

Export Citation Format

Share Document