Forward modelling of spectral depths using 3D Fourier convolution
Roger CliftonNT Geological Survey, GPO Box 4550, Darwin, NT 0801, Australia, and Centre for Exploration Targeting, The University of Western Australia, 35 Stirling Highway, Crawley, WA 6009, Australia. Email: Roger.Clifton@NT.gov.au
Exploration Geophysics 48(2) 166-176 https://doi.org/10.1071/EG15092
Submitted: 3 September 2015 Accepted: 5 December 2015 Published: 29 January 2016
Abstract
A procedure for creating depth transects across the Northern Territory (Australia) has been established using magnetic spectral depths. Although these transects showed depths to layers, they failed to show depths in the upper few hundred metres. Because the magnetic depths had been derived on the basis of heterogeneity of the bodies, a method of forward modelling based on heterogeneity is needed to explain this limitation and other issues. Spatial convolutions based on heterogeneity suffice for primitive models but are too slow for detailed work.
This paper demonstrates fast forward modelling using Fourier convolution, that is convolution of three-dimensional (3D) arrays via the frequency domain, to obtain total magnetic intensity grids and magnetic depth profiles for hypothetical structures. Randomly located dipoles are used to simulate the heterogeneity of the material of modelled bodies.
The loss of shallow depth signal in the magnetic transects is shown to arise mainly from the limitation of the line spacing of the underlying surveys. Depths of bodies at less than half the line spacing of the survey are not resolved at all and depths less than the line spacing itself appear deeper than the actual source depth.
Fourier convolution works equally well with non-layered, non-prismatic bodies. Modelling of an inclined, elliptical body is demonstrated by way of example. The associated depth profile shows a clear equivalent layer at a depth representative for such a body. The result allows interpretation of a characteristic pattern in magnetic depth transects as indicating the depth to a relatively compact non-layered body.
Fourier convolution showed a considerable speed advantage over spatial convolution at all array sizes used in the study. Convolutions of model arrays of 1000 × 1000 × 500 were calculated within a few minutes.
Key words: Fourier modelling, heterogeneity, magnetic depths, random dipoles.
References
Blakely, R. J., 1995, Potential theory in gravity and magnetic applications: Cambridge University Press.Caratori Tontini, F., 2012, Rapid interactive modelling of 3-D magnetic anomalies: Computers & Geosciences, 48, 308–315
| Rapid interactive modelling of 3-D magnetic anomalies:Crossref | GoogleScholarGoogle Scholar |
Caratori Tontini, F., Cocchi, L., and Carmisciano, C., 2009, Rapid 3-D forward model of potential fields with application to the Palinuro Seamount magnetic anomaly (Southern Tyrrhenian Sea, Italy): Journal of Geophysical Research, 114, B02103
| Rapid 3-D forward model of potential fields with application to the Palinuro Seamount magnetic anomaly (Southern Tyrrhenian Sea, Italy):Crossref | GoogleScholarGoogle Scholar |
Clifton, R., 2011, NT wide geophysical stitch, magnetics: Northern Territory Geological Survey.
Clifton, R., 2013, Magnetic depth profiles across the Northern Territory: NTGS Technical Note 2013–001, Northern Territory Geological Survey.
Clifton, R., 2015, Magnetic depths to basalts: extension of spectral depths method: Exploration Geophysics, 46, 284–296
| Magnetic depths to basalts: extension of spectral depths method:Crossref | GoogleScholarGoogle Scholar |
Krahenbuhl, R. A., and Li, Y., 2007, Influence of self-demagnetization effect on data interpretation in strongly magnetic environments: ASEG Extended Abstracts, 1–4.
Kruse, P. D., and Maier, R. C., 2010, Frew River, Northern Territory (second edition revised). 1: 250 000 geological map series, SF 53–03: Northern Territory Geological Survey, Darwin.
Phillips, J. D., 2014, Using vertical Fourier transforms to invert potential-field data to magnetization or density models in the presence of topography: SEG Technical Program Expanded Abstracts, 1339–1343.
Reid, A. B., 1980, Aeromagnetic survey design: Geophysics, 45, 973–976
| Aeromagnetic survey design:Crossref | GoogleScholarGoogle Scholar |
Sheriff, R. E., and Geldart, L. P., 1983, Exploration seismology, vol. 2: Cambridge University Press.
Spector, A., and Grant, F. S., 1970, Statistical models for interpreting aeromagnetic data: Geophysics, 35, 293–302
| Statistical models for interpreting aeromagnetic data:Crossref | GoogleScholarGoogle Scholar |
University of British Columbia, 2005, MAG3D software. Available at: http://www.eos.ubc.ca/research/ubcgif/iag/sftwrdocs/mag3d/index.htm
Vacquier, V., 1972, Geomagnetism in marine geology: Elsevier.