1887
Volume 48, Issue 3
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

[

We describe a new fast gravity inversion method to recover a 3D density model from gravity data. This gravity inversion method introduces a stabiliser model norm with a depth weighting function to produce smooth models and uses a new symmetric successive over-relaxation (SSOR) iterative conjugate gradient (CG) algorithm.

,

The subsurface three-dimensional (3D) model of density distribution is obtained by solving an under-determined linear equation that is established by gravity data. Here, we describe a new fast gravity inversion method to recover a 3D density model from gravity data. The subsurface will be divided into a large number of rectangular blocks, each with an unknown constant density. The gravity inversion method introduces a stabiliser model norm with a depth weighting function to produce smooth models. The depth weighting function is combined with the model norm to counteract the skin effect of the gravity potential field. As the numbers of density model parameters is NZ (the number of layers in the vertical subsurface domain) times greater than the observed gravity data parameters, the inverse density parameter is larger than the observed gravity data parameters. Solving the full set of gravity inversion equations is very time-consuming, and applying a new algorithm to estimate gravity inversion can significantly reduce the number of iterations and the computational time. In this paper, a new symmetric successive over-relaxation (SSOR) iterative conjugate gradient (CG) method is shown to be an appropriate algorithm to solve this Tikhonov cost function (gravity inversion equation). The new, faster method is applied on Gaussian noise-contaminated synthetic data to demonstrate its suitability for 3D gravity inversion. To demonstrate the performance of the new algorithm on actual gravity data, we provide a case study that includes ground-based measurement of residual Bouguer gravity anomalies over the Humble salt dome near Houston, Gulf Coast Basin, off the shore of Louisiana. A 3D distribution of salt rock concentration is used to evaluate the inversion results recovered by the new SSOR iterative method. In the test model, the density values in the constructed model coincide with the known location and depth of the salt dome.

]
ASEG
Loading

Article metrics loading...

/content/journals/10.1071/EG15041
2017-09-01
2026-01-18

  • From This Site

    /content/journals/10.1071/EG15041
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Abdelrahman E. M. Bayoumi A. I. El-Araby H. M. 1991 A least-squares minimization approach to invert gravity data: Geophysics 56 115 118 10.1190/1.1442946
    https://doi.org/10.1190/1.1442946 [Google Scholar]
  2. Abdelrahman E. S. M. El-Araby T. M. El-Araby H. M. Abo-Ezz E. R. 2001 a A new method for shape and depth determinations from gravity data: Geophysics 66 1774 1780 10.1190/1.1487119
    https://doi.org/10.1190/1.1487119 [Google Scholar]
  3. Abdelrahman E. M. El-Araby H. M. El-Araby T. M. Abo-Ezz E. R. 2001 b Three least-squares minimization approaches to depth, shape, and amplitude coefficient determination from gravity data: Geophysics 66 1105 1109 10.1190/1.1487058
    https://doi.org/10.1190/1.1487058 [Google Scholar]
  4. Abedi M. Gholami A. Norouzi G. H. Fathianpour N. 2013 Fast inversion of magnetic data using Lanczos bidiagonalization method: Journal of Applied Geophysics 90 126 137 10.1016/j.jappgeo.2013.01.008
    https://doi.org/10.1016/j.jappgeo.2013.01.008 [Google Scholar]
  5. Bear G. W. Al-Shukri H. J. Rudman A. J. 1995 Linear inversion of gravity data for 3-D density distributions: Geophysics 60 1354 1364 10.1190/1.1443871
    https://doi.org/10.1190/1.1443871 [Google Scholar]
  6. Bosch M McGaughey J 2001 Joint inversion of gravity and magnetic data under lithologic constraints: The Leading Edge 20 877 881
    [Google Scholar]
  7. Botros Y. Y. Volakis J. L. 1999 Preconditioned generalized minimal residual iterative scheme for perfectly matched layer terminated applications: IEEE Microwave and Guided Wave Letters 9 45 47 10.1109/75.755039
    https://doi.org/10.1109/75.755039 [Google Scholar]
  8. Braile L. W. Keller G. R. Peeples W. J. 1974 Inversion of gravity data for two‐dimensional density distributions: Journal of Geophysical Research 79 2017 2021 10.1029/JB079i014p02017
    https://doi.org/10.1029/JB079i014p02017 [Google Scholar]
  9. Canning F. X. Scholl J. F. 1996 Diagonal preconditioners for the EFIE using a wavelet basis: IEEE Transactions on Antennas and Propagation 44 1239 1246 10.1109/8.535382
    https://doi.org/10.1109/8.535382 [Google Scholar]
  10. Caratori Tontini F Cocchi L Carmisciano C 2006 Depth-to-the-bottom optimization for magnetic data inversion: magnetic structure of the Latium volcanic region, Italy: Journal of Geophysical Research: Solid Earth 111 B11104
    [Google Scholar]
  11. Chasseriau P Chouteau M 2003 3D gravity inversion using a model of parameter covariance: Journal of Applied Geophysics 52 59 74
    [Google Scholar]
  12. Chen R. S. Yung E. K. Chan C. H. Fang D. G. 2000 Application of preconditioned CG–FFT technique to method of lines for analysis of the infinite-plane metallic grating: Microwave and Optical Technology Letters 24 170 175 10.1002/(SICI)1098‑2760(20000205)24:3<170::AID‑MOP8>3.0.CO;2‑S
    https://doi.org/10.1002/(SICI)1098-2760(20000205)24:3<170::AID-MOP8>3.0.CO;2-S [Google Scholar]
  13. Chen R. S. Yung E. K. N. Chan C. H. Wang D. X. Fang D. G. 2002 Application of the SSOR preconditioned CG algorithm to the vector FEM for 3D full-wave analysis of electromagnetic-field boundary-value problems: IEEE Transactions on Microwave Theory and Techniques 50 1165 1172 10.1109/22.993420
    https://doi.org/10.1109/22.993420 [Google Scholar]
  14. Chen, X., 2005, Preconditioners for iterative solutions of large-scale linear systems arising from Biot’s consolidation equations: Ph.D. thesis, National University of Singapore.
  15. Čuma M. Zhdanov M. S. 2014 Massively parallel regularized 3D inversion of potential fields on CPUs and GPUs: Computers & Geosciences 62 80 87 10.1016/j.cageo.2013.10.004
    https://doi.org/10.1016/j.cageo.2013.10.004 [Google Scholar]
  16. Ellis R. G. Oldenburg D. W. 1994 The pole-pole 3-D Dc-resistivity inverse problem: a conjugate gradient approach: Geophysical Journal International 119 187 194 10.1111/j.1365‑246X.1994.tb00921.x
    https://doi.org/10.1111/j.1365-246X.1994.tb00921.x [Google Scholar]
  17. Essa K. S. 2007 Gravity data interpretation using the s-curves method: Journal of Geophysics and Engineering 4 204 213 10.1088/1742‑2132/4/2/009
    https://doi.org/10.1088/1742-2132/4/2/009 [Google Scholar]
  18. Fisher N. J. Howard L. E. 1980 Gravity interpretation with the aid of quadratic programming: Geophysics 45 403 419 10.1190/1.1441090
    https://doi.org/10.1190/1.1441090 [Google Scholar]
  19. Golub, G. H., and Van Loan, C. F., 1996, Matrix computations (3rd edition): Johns Hopkins University Press.
  20. Haáz I. B. 1953 Relationship between the potential of the attraction of the mass contained in a finite rectangular prism and its first and second derivatives: Geophysical Transactions II 57 66
    [Google Scholar]
  21. Hinze W. J. 1990 The role of gravity and magnetic methods in engineering and environmental studies: Geotechnical and Environmental Geophysics 1 75 126
    [Google Scholar]
  22. Jackson D. D. 1979 The use of a priori data to resolve non-uniqueness in linear inversion: Geophysical Journal International 57 137 157 10.1111/j.1365‑246X.1979.tb03777.x
    https://doi.org/10.1111/j.1365-246X.1979.tb03777.x [Google Scholar]
  23. Li X. Chouteau M. 1998 Three-dimensional gravity modeling in all space: Surveys in Geophysics 19 339 368 10.1023/A:1006554408567
    https://doi.org/10.1023/A:1006554408567 [Google Scholar]
  24. Li Y. Oldenburg D. W. 1996 3-D inversion of magnetic data: Geophysics 61 394 408 10.1190/1.1443968
    https://doi.org/10.1190/1.1443968 [Google Scholar]
  25. Li Y. Oldenburg D. W. 1998 3-D inversion of gravity data: Geophysics 63 109 119 10.1190/1.1444302
    https://doi.org/10.1190/1.1444302 [Google Scholar]
  26. Mackie R. L. Madden T. R. 1993 Conjugate direction relaxation solutions for 3-D magnetotelluric modeling: Geophysics 58 1052 1057 10.1190/1.1443481
    https://doi.org/10.1190/1.1443481 [Google Scholar]
  27. Malehmir A. Thunehed H. Tryggvason A. 2009 Case history: the Paleoproterozoic Kristineberg mining area, northern Sweden: results from integrated 3D geophysical and geologic modeling, and implications for targeting ore deposits: Geophysics 74 B9 B22 10.1190/1.3008053
    https://doi.org/10.1190/1.3008053 [Google Scholar]
  28. Mareschal J. C. 1985 Inversion of potential field data in Fourier transform domain: Geophysics 50 685 691 10.1190/1.1441943
    https://doi.org/10.1190/1.1441943 [Google Scholar]
  29. Mohan N. L. Anandababu L. Rao S. S. 1986 Gravity interpretation using the Mellin transform: Geophysics 51 114 122 10.1190/1.1442024
    https://doi.org/10.1190/1.1442024 [Google Scholar]
  30. Najafi H. S. Edalatpanah S. A. 2014 A new modified SSOR iteration method for solving augmented linear systems: International Journal of Computer Mathematics 91 539 552 10.1080/00207160.2013.792923
    https://doi.org/10.1080/00207160.2013.792923 [Google Scholar]
  31. Nakatsuka T. 1995 Minimum norm inversion of magnetic anomalies with application to aeromagnetic data in the Tanna area, central Japan: Journal of Geomagnetism and Geoelectricity 47 295 311 10.5636/jgg.47.295
    https://doi.org/10.5636/jgg.47.295 [Google Scholar]
  32. Namaki L. Gholami A. Hafizi M. A. 2011 Edge-preserved 2-D inversion of magnetic data: an application to the Makran arc-trench complex: Geophysical Journal International 184 1058 1068 10.1111/j.1365‑246X.2010.04877.x
    https://doi.org/10.1111/j.1365-246X.2010.04877.x [Google Scholar]
  33. Nettleton L. L. 1962 Gravity and magnetics for geologists and seismologists: AAPG Bulletin 46 1815 1838
    [Google Scholar]
  34. Nettleton, L. L., 1976, Gravity and magnetics in oil prospecting: McGraw-Hill Companies.
  35. Nocedal, J., and Wright, S., 2006, Numerical optimization: Springer Science and Business Media.
  36. Nolet G. 1985 Solving or resolving inadequate and noisy tomographic systems: Journal of Computational Physics 61 463 482 10.1016/0021‑9991(85)90075‑0
    https://doi.org/10.1016/0021-9991(85)90075-0 [Google Scholar]
  37. Nolet, G., 1993, Solving large linearized tomographic problems, in H. M. Iyer, and K. Hirahara, eds., Seismic tomography: theory and practice: Chapman and Hall, 227–247.
  38. Oldenburg D. W. Li Y. 1994 Inversion of induced polarization data: Geophysics 59 1327 1341 10.1190/1.1443692
    https://doi.org/10.1190/1.1443692 [Google Scholar]
  39. Oruç B. 2010 Depth estimation of simple causative sources from gravity gradient tensor invariants and vertical component: Pure and Applied Geophysics 167 1259 1272 10.1007/s00024‑009‑0021‑4
    https://doi.org/10.1007/s00024-009-0021-4 [Google Scholar]
  40. Paterson N. R. Reeves C. V. 1985 Applications of gravity and magnetic surveys: the state-of-the-art in 1985: Geophysics 50 2558 2594 10.1190/1.1441884
    https://doi.org/10.1190/1.1441884 [Google Scholar]
  41. Pignatelli A. Nicolosi I. Chiappini M. 2006 An alternative 3D source inversion method for magnetic anomalies with depth resolution: Annals of Geophysics 49 1021 1027
    [Google Scholar]
  42. Pilkington M. 1997 3-D magnetic imaging using conjugate gradients: Geophysics 62 1132 1142 10.1190/1.1444214
    https://doi.org/10.1190/1.1444214 [Google Scholar]
  43. Portniaguine O. Zhdanov M. S. 1999 Focusing geophysical inversion images: Geophysics 64 874 887 10.1190/1.1444596
    https://doi.org/10.1190/1.1444596 [Google Scholar]
  44. Portniaguine O. Zhdanov M. S. 2002 3-D magnetic inversion with data compression and image focusing: Geophysics 67 1532 1541 10.1190/1.1512749
    https://doi.org/10.1190/1.1512749 [Google Scholar]
  45. Saad, Y., 2003, Iterative methods for sparse linear systems: SIAM.
  46. Safon C Vasseur G Cuer M 1977 Some applications of linear programming to the inverse gravity problem: Geophysics 42 1215 1229
    [Google Scholar]
  47. Salem A. Ravat D. Mushayandebvu M. F. Ushijima K. 2004 Linearized least-squares method for interpretation of potential-field data from sources of simple geometry: Geophysics 69 783 788 10.1190/1.1759464
    https://doi.org/10.1190/1.1759464 [Google Scholar]
  48. Sarkar T Arvas E 1985 On a class of finite step iterative methods (conjugate directions) for the solution of an operator equation arising in electromagnetics: IEEE Transactions on Antennas and Propagation 33 1058 1066
    [Google Scholar]
  49. Scales J. A. 1987 Tomographic inversion via the conjugate gradient method: Geophysics 52 179 185
    [Google Scholar]
  50. Shamsipour P Marcotte D Chouteau M Keating P 2010 3D stochastic inversion of gravity data using cokriging and cosimulation: Geophysics 75 I1 I10
    [Google Scholar]
  51. Shamsipour P. Marcotte D. Chouteau M. Rivest M. Bouchedda A. 2013 3D stochastic gravity inversion using nonstationary covariances: Geophysics 78 G15 G24 10.1190/geo2012‑0122.1
    https://doi.org/10.1190/geo2012-0122.1 [Google Scholar]
  52. Shaw R. K. Agarwal B. N. P. 1997 A generalized concept of resultant gradient to interpret potential field maps: Geophysical Prospecting 45 1003 1011 10.1046/j.1365‑2478.1997.740310.x
    https://doi.org/10.1046/j.1365-2478.1997.740310.x [Google Scholar]
  53. Shi, Y. F., Chen, R. S., and Xia, M. Y., 2008, SSOR preconditioner accelerated time domain finite element boundary integral method, in Asia-Pacific Microwave Conference 2008: IEEE, 1–4.
  54. Smith, G. D., 1985, Numerical solution of partial differential equations: finite difference methods: Oxford University Press.
  55. Tikhonov, A. N., and Arsenin, V. I., 1977, Solutions of ill-posed problems: Winston.
  56. VanDecar J. C. Snieder R. 1994 Obtaining smooth solutions to large, linear, inverse problems: Geophysics 59 818 829 10.1190/1.1443640
    https://doi.org/10.1190/1.1443640 [Google Scholar]
  57. Ward, S. H., 1990, Geotechnical and environmental geophysics: SEG.
  58. Zhang J. Mackie R. L. Madden T. R. 1995 3-D resistivity forward modeling and inversion using conjugate gradients: Geophysics 60 1313 1325 10.1190/1.1443868
    https://doi.org/10.1190/1.1443868 [Google Scholar]
  59. Zhdanov, M. S., 1988, Integral transforms in geophysics: Springer Verlag.
  60. Zhdanov, M. S., 2002, Geophysical inverse theory and regularization problems: Elsevier.
  61. Zhdanov, M. S., 2009, Geophysical electromagnetic theory and methods: Elsevier.
  62. Zhdanov, M. S., and Wannamaker, P. E., 2002, Three-dimensional electromagnetics: Elsevier.
/content/journals/10.1071/EG15041
Loading
Fast inversion of gravity data using the symmetric successive over-relaxation (SSOR) preconditioned conjugate gradient algorithm
Exploration Geophysics 48, 294 (2017); https://doi.org/10.1071/EG15041
/content/journals/10.1071/EG15041
/content/journals/10.1071/EG15041
Loading

Data & Media loading...

  • Article Type: Research Article
Keyword(s): conjugate gradient, depth of the weighting function, gravity inversion, symmetric successive over-relaxation.

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error