AN AUTOMATIC LEAST‐SQUARES MULTIMODEL METHOD FOR MAGNETIC INTERPRETATION

Geophysics ◽  
1973 ◽  
Vol 38 (2) ◽  
pp. 349-358 ◽  
Author(s):  
Peter H. McGrath ◽  
Peter J. Hood

The magnetic anomalies caused by such diverse model shapes as the finite strike length thick dike, the vertical prism, thes loping step, the parallelepiped body, etc., may be obtained through an appropriate numerical integration of the expression for the magnetic effect produced by a finite thin plate. Using models generated in this manner, an automatic computer method has been developed at the Geological Survey of Canada for the interpretation of magnetic data. Because the magnetic anomalies produced by the various model shapes are nonlinear in parameters of shape and position, it is necessary to use an iterative procedure to obtain the values for the various model parameters which yield a least‐squares best‐fit anomaly curve to a set of discrete observed data. The interpretation method described in this paper uses the Powell algorithm for this purpose. The procedure can sometimes be made more efficient using a Marquardt modification to the Powell algorithm. Examples of the use of the method are presented for an elongated anomaly in the Moose River basin of the Hudson Bay lowlands in northern Ontario, and for an areally large elliptical anomaly in the Sverdrup basin of the Canadian Arctic Islands.

1977 ◽  
Vol 17 (1) ◽  
pp. 85
Author(s):  
Robert J. Whiteley ◽  
Barry F. Long ◽  
David A. Pratt

The magnetic method is used at many stages of a modern petroleum exploration program. Effective interpretation techniques are required to extract maximum geological information from magnetic data. Those techniques which provide the greatest flexibility and make full use of the talents of experienced interpreters are generally of a semi-automated and interactive nature.There are several practical methods for semi-automated quantitative magnetic interpretation in sedimentary basins. Initial interpretation can be achieved by automatic calculation of characteristic anomaly parameters continuously along original or processed magnetic data profiles. Detailed interpretation of more subtle magnetic features can then follow by theoretical anomaly comparison with field anomalies using interactive portfolio modelling or by direct computation.Examples of the use of these semi-automated techniques in the interpretation of basement and intra-sedimentary magnetic anomalies show that combined magnetic and seismic interpretations can provide considerable insight into the structural processes which have operated in a sedimentary basin.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. L35-L42 ◽  
Author(s):  
Mark Pilkington

Gravity and magnetic data are inverted jointly in terms of a model consisting of an interface separating two layers having a constant density and magnetization contrast. A damped least-squares inversion is used to determine the topography of the interface. The inversion requires knowledge of the physical property contrasts across the interface and its average depth. Since the relationship between model parameters and data is weakly nonlinear, a constant damped least-squares inverse is used during the iterative solution search. The effect of this inverse is closely related to a downward continuation of the field to the average interface depth. The method is used to map the base of the Sept-Iles mafic intrusion, Quebec, Canada, and the shape of the central uplift at the Chicxulub impact crater, Yucatan, Mexico. At Sept-Iles, the intrusion reaches a thickness of [Formula: see text], coincident with the maximum gravity anomaly, south of the intrusion center. At Chicxulub, the top of the central uplift is modeled to be [Formula: see text] deep and has a single peak form.


Geophysics ◽  
1978 ◽  
Vol 43 (3) ◽  
pp. 634-636 ◽  
Author(s):  
M. V. Ramanaiah Chowdary

A great deal of interest has been shown in the frequency analysis of gravity and magnetic data originally suggested by Dean (1958). The application of this method for potential field problems has met with considerable success. The purpose of this note is to show that the interpretation of total magnetic anomalies due to a sloping step model, which represents a contact between zones having different magnetic properties in terms of model parameters, is less complicated in the frequency domain than in the spatial domain.


2004 ◽  
Vol 92 (6) ◽  
pp. 1036-1053 ◽  
Author(s):  
PAUL H. GLASER ◽  
BARBARA C. S. HANSEN ◽  
DONALD I. SIEGEL ◽  
ANDREW S. REEVE ◽  
PAUL J. MORIN

Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. L7-L15 ◽  
Author(s):  
Mark Pilkington

I have developed an inversion approach that determines a 3D susceptibility distribution that produces a given magnetic anomaly. The subsurface model consists of a 3D, equally spaced array of dipoles. The inversion incorporates a model norm that enforces sparseness and depth weighting of the solution. Sparseness is imposed by using the Cauchy norm on model parameters. The inverse problem is posed in the data space, leading to a linear system of equations with dimensions based on the number of data, [Formula: see text]. This contrasts with the standard least-squares solution, derived through operations within the [Formula: see text]-dimensional model space ([Formula: see text] being the number of model parameters). Hence, the data-space method combined with a conjugate gradient algorithm leads to computational efficiency by dealing with an [Formula: see text] system versus an [Formula: see text] one, where [Formula: see text]. Tests on synthetic data show that sparse inversion produces a much more focused solution compared with a standard model-space, least-squares inversion. The inversion of aeromagnetic data collected over a Precambrian Shield area again shows that including the sparseness constraint leads to a simpler and better resolved solution. The degree of improvement in model resolution for the sparse case is quantified using the resolution matrix.


2021 ◽  
pp. 1-14
Author(s):  
Marissa A. Davies ◽  
Jerome Blewett ◽  
B. David A. Naafs ◽  
Sarah A. Finkelstein

Abstract A multiproxy Holocene record from a bog in the Hudson Bay Lowlands, northern Ontario, Canada, was used to evaluate how ecohydrology relates to carbon accumulation. The study site is located at a somewhat higher elevation and on coarser grained deposits than the surrounding peatlands. This promotes better drainage and thus a slower rate of carbon accumulation relative to sites with similar initiation age. The rate of peat vertical accretion was initially low as the site transitioned from a marsh to a rich fen. These lower rates took place during the warmer temperatures of the Holocene thermal maximum, confirming the importance of hydrological controls limiting peat accretion at the local scale. Testate amoebae, pollen, and plant macrofossils indicate a transition to a poor fen and then a bog during the late Holocene, as the carbon accumulation rate and reconstructed water table depth increased. The bacterial membrane lipid biomarker indices used to infer paleotemperature show a summer temperature bias and appear sensitive to changes in peat type. The bacterial membrane lipid biomarker pH proxy indicates a rich to a poor fen and a subsequent fen to bog transition, which are supported by pollen, macrofossil, and testate amoeba records.


2012 ◽  
Vol 78 (2) ◽  
pp. 275-284 ◽  
Author(s):  
Joan Bunbury ◽  
Sarah A. Finkelstein ◽  
Jörg Bollmann

AbstractMultiple proxies from a 319-cm peat core collected from the Hudson Bay Lowlands, northern Ontario, Canada were analyzed to determine how carbon accumulation has varied as a function of paleohydrology and paleoclimate. Testate amoeba assemblages, analysis of peat composition and humification, and a pollen record from a nearby lake suggest that isostatic rebound and climate may have influenced peatland growth and carbon dynamics over the past 6700 cal yr BP. Long-term apparent rates of carbon accumulation ranged between 8.1 and 36.7 g C m− 2 yr− 1 (average = 18.9 g C m− 2 yr− 1). The highest carbon accumulation estimates were recorded prior to 5400 cal yr BP when a fen existed at this site, however following the fen-to-bog transition carbon accumulation stabilized. Carbon accumulation remained relatively constant through the Neoglacial period after 2400 cal yr BP when pollen-based paleoclimate reconstructions from a nearby lake (McAndrews et al., 1982) and reconstructions of the depth to the water table derived from testate amoeba data suggest a wetter climate. More carbon accumulated per unit time between 1000 and 600 cal yr BP, coinciding in part with the Medieval Climate Anomaly.


2019 ◽  
Vol 6 (1.) ◽  
Author(s):  
Nasar Aldian Ambark Shashoa ◽  
Abdurrezag S. Elmezughi ◽  
Ibrahim N. Jleta ◽  
Nasser B. Ekreem

In this paper, a new-type recursive least squares algorithm is proposed for identifying the system model parameters and the noise model parameters of Box–Jenkins Systems. The basic idea is based on replacing the unmeasurable variables in the information vectors with their estimates. The proposed algorithm has high computational efficiency because the dimensions of the involved covariance matrices in each subsystem become small. Validation of the model is evaluated using some statistical methods, Which, best-fit criterion and Histogram. Simulation results are presented to illustrate the effectiveness of the proposed algorithm.


Geophysics ◽  
1993 ◽  
Vol 58 (8) ◽  
pp. 1214-1214 ◽  
Author(s):  
N. L. Mohan

It is quite interesting to learn the 3-D analytic signal interpretation of Roest et al. (1992), using the vector addition. However, I am quite skeptical about their objective of determining the depth under the assumption that the magnetic anomalies are caused by vertical contacts from gridded magnetic data which, it appears to me, is nothing but an oversimplification of interpretation.


Sign in / Sign up

Export Citation Format

Share Document