A case study of forward calculations of the gravity anomaly by spectral method for a three-dimensional parameterised fault model

Spectral methods provide many advantages for calculating gravity anomalies. In this paper, we derive a kernel function for a three-dimensional (3D) fault model in the wave number domain, and present the full Fortran source code developed for the forward computation of the gravity anomalies and related derivatives obtained from the model. The numerical error and computing speed obtained using the proposed spectral method are compared with those obtained using a 3D rectangular prism model solved in the space domain. The error obtained using the spectral method is shown to be dependent on the sequence length employed in the fast Fourier transform. The spectral method is applied to some examples of 3D fault models, and is demonstrated to be a straightforward and alternative computational approach to enhance computational speed and simplify the procedures for solving many gravitational potential forward problems involving complicated geological models. The proposed method can generate a great number of feasible geophysical interpretations based on a 3D model with only a few variables, and can thereby improve the efficiency of inversion. A kernel function for a parameterized fault model in the wave number domain is derived.Multi-combined fault models as example have been forwarded.The accuracy of spectral method can be controllable by the adequate FFT sequence length.

[1]  William J. Hinze,et al.  The Effect of Finite Data Length in the Spectral Analysis of Ideal Gravity Anomalies , 1976 .

[2]  B. Bhattacharyya,et al.  A FAST FOURIER TRANSFORM METHOD FOR RAPID COMPUTATION OF GRAVITY AND MAGNETIC ANOMALIES DUE TO ARBITRARY BODIES , 1976 .

[3]  William J. Hinze,et al.  Gravity and Magnetic Exploration: Principles, Practices, and Applications , 2013 .

[4]  T. K. Bose,et al.  A GENERAL EXPRESSION FOR THE FOURIER TRANSFORM OF THE GRAVITY ANOMALY DUE TO A FAULT , 1977 .

[5]  Robert L. Parker,et al.  The Rapid Calculation of Potential Anomalies , 1973 .

[6]  D. E. Gill,et al.  Interpretation of gravity anomalies of dike-like bodies by Fourier transformation , 1970 .

[7]  B. Bhattacharyya,et al.  Spectral analysis of gravity and magnetic anomalies due to rectangular prismatic bodies , 1977 .

[8]  William J. Hinze,et al.  THEORETICAL TRANSFORMS OF THE GRAVITY ANOMALIES OF TWO IDEALIZED BODIES , 1978 .

[9]  William J. Hinze,et al.  Fourier transforms of finite length theoretical gravity anomalies. , 1977 .

[10]  Manik Talwani,et al.  Rapid computation of gravitational attraction of three-dimensional bodies of arbitrary shape , 1960 .

[11]  Javier Fullea,et al.  FA2BOUG - A FORTRAN 90 code to compute Bouguer gravity anomalies from gridded free-air anomalies: Application to the Atlantic-Mediterranean transition zone , 2008, Comput. Geosci..

[12]  D. Nagy.  The gravitational attraction of a right rectangular prism , 1966 .

[13]  Fabio Rocca,et al.  Interpretation of Magnetic Anomalies Using Spectral Estimation TECHNIQUES , 1975 .

[14]  S. K. Singh,et al.  Spectral analysis of gravity and magnetic anomalies due to a vertical circular cylinder , 1983 .

[15]  Sonali Das,et al.  Interpretation of the gravitational effect of a two-dimensional dike by Fourier transform , 1977 .

[16]  Houze Xu,et al.  Three-dimensional forward and inverse models for gravity fields based on the Fast Fourier Transform , 2006, Comput. Geosci..

[17]  S. Chacko,et al.  A method for analysing gravity anomalies due to a geologic contact by Fourier transform , 1980 .

[18]  David Gómez-Ortiz,et al.  3DINVER.M: a MATLAB program to invert the gravity anomaly over a 3D horizontal density interface by Parker-Oldenburg's algorithm , 2005, Comput. Geosci..

[19]  D. Oldenburg.  The inversion and interpretation of gravity anomalies , 1974 .

[20]  B. Bhattacharyya,et al.  Some general properties of potential fields in space and frequency domain: A review , 1967 .