Gravity inversion using wavelet-based compression on parallel hybrid cpu/gpu systems: Application to southwest ghana

Geophysical Journal International

Volumn 195, Issue 3, 2013, Pages 1594-1619

  Martin, Roland ^a   Monteiller, Vadim ^b,c   Komatitsch, Dimitri ^b   Perrouty, Stéphane ^a   Jessell, Mark ^a   Bonvalot, Sylvain ^a   Lindsay, Mark ^a,d  

^a CNRS   (France)

^b AIX MARSEILLE UNIVERSITY   (France)

^c UNIVERSITÉ DE PAU ET DES PAYS DE L'ADOUR   (France)

^d MONASH UNIVERSITY   (Australia)

Author keywords

Gravity anomalies and Earth structure; Inverse theory; Numerical approximations and analysis; Satellite gravity; Wavelet transform

Indexed keywords

DENSITY DISTRIBUTIONS; GRAPHICS PROCESSING UNIT; GRAVITY ANOMALIES AND EARTH STRUCTURES; INVERSE THEORY; NUMERICAL APPROXIMATIONS AND ANALYSIS; SATELLITE GRAVITY; THREE ORDERS OF MAGNITUDE; WAVELET-BASED COMPRESSIONS;

COMPUTER GRAPHICS; DISCRETE WAVELET TRANSFORMS; INVERSE PROBLEMS; ITERATIVE METHODS; LINEAR SYSTEMS; MATRIX ALGEBRA; PROGRAM PROCESSORS; WAVELET TRANSFORMS;

PARALLEL PROCESSING SYSTEMS;

GRAVITY ANOMALY; GRAVITY FIELD; INVERSE PROBLEM; LEAST SQUARES METHOD; MATHEMATICAL ANALYSIS; PARALLEL COMPUTING; SOFTWARE; TRANSFORM;

GHANA;

EID: 84887551043     PISSN: 0956540X     EISSN: 1365246X     Source Type: Journal    
DOI: 10.1093/gji/ggt334     Document Type: Article

References (91)

1. Adadey, K., Clarke, B., Th́eveniaut, H., Urien, P., Delor, C., Roig, J.Y. & Feybesse, J.L., 2009. Geological map explanation: map sheet 0503 B (1:100 000), CGS/BRGM/Geoman, Tech. Rep, Geological Survey Department of Ghana (GSD). MSSP/2005/GSD/5a.
2. Agarwal, B.N.P.&Srivastava, S., 2010.AFORTRAN program to implement the method of finite elements to compute regional and residual anomalies from gravity data, Comput. Geosci., 36, 848-852.
3. Allibone, A., McCuaig, T., Harris, D., Etheridge, M., Munroe, S. & Byrne, D., 2002. Structural controls on gold mineralization at the Ashanti Gold Deposit, Obuasi, Ghana, Soc. Econ. Geologist., 9, 65-93. [Special publication]
4. Andersen, O., 2010. The DTU10 global gravity field and mean sea surface: improvements in the Arctic, in Proceedings of the 2nd IGFS Meeting, Fairbanks, Alaska, USA.
5. Argüello, F., Heras, D.B., Bóo, M. & Lamas-Rodríguez, J., 2012. The splitand-merge method in general purpose computation on GPUs, Parallel Comput., 38, 277-288.
6. s and density structure beneath the Colorado Plateau constrained by gravity anomalies and joint inversions of receiver function and phase velocity data, J. geophys. Res., 117, B02313, doi:10.1029/2011JB008522.
7. Balmino, G., Vales, N., Bonvalot, S. & Briais, A., 2012. Spherical harmonic modelling to ultrahigh degree of Bouguer and isostatic anomalies, J. Geod., 86(7), 499-520.
8. Barritt, S. & Kuma, J., 1998. Constrained gravity models and structural evolution of the Ashanti Belt, southwest Ghana, J. Afr. Earth Sci., 26(4), 539-550.
9. Battacharyya, B.K., 1980. A generalized multibody model for inversion of magnetic anomalies, Geophysics, 45, 255-270.
10. Beck, A. & Teboulle, M., 2009. A fast iterative shrinkage-thresholding algorithm for linear inverse problems, SIAM J. Imag. Sci., 2(1), 183- 202.
11. Blakely, R.J., 1995. Potential Theory in Gravity and Magnetic Applications, Cambridge Univ. Press.
12. Bonhomme,M., 1962. Contribution à l'étude géochronologique de la plateforme de l'Ouest Africain, Ann. Fac. Sci., Université de Clermont-Ferrand, France, Géologie, Minéralogie, 5, 1-62.
13. Bonvalot, S. et al. Eds., 2012.World Gravity Map, BGI/CGMW/CNES/IRD.
14. Boschetti, F., Hornby, P. & Horowitz, F.G., 2001. Wavelet based inversion of gravity data, Explor. Geophys., 32, 48-55.
15. Calcagno, P., Chilès, J., Courrioux, G. & Guillen, A., 2008. Geological modelling from field data and geological knowledge. Part I. Modelling method coupling 3D potentiel-field interpolation and geological rules, Phys. Earth planet. Inter., 171, 147-157.
16. Chapman, C. & Pratt, R.G., 1992. Traveltime tomography in anisotropic media-I: theory, Geophys. J. Int., 109, 1-19.
17. Chevrot, S., Martin, R. & Komatitsch, D., 2012. Optimized discrete wavelet transforms in the cubed-sphere with the lifting scheme: implications for global finite-frequency tomography, Geophys. J. Int., 191(3), 1391-1402.
18. Cohen, A., Daubechies, I. & Fauveau, J., 1992. Biorthogonal bases of compactly supported wavelets, Commun. Pure appl. Math., 45(5), 485-560.
19. Commer, M., 2011. Three-dimensional gravity modelling and focusing inversion using rectangular meshes, Geophys. Prospect., 59, 966-979.
20. Cuma, M., Wilson, G. & Zhdanov, M., 2012. Large-scale 3D inversion of potential data, Geophys. Prospect., 60, 1186-1199.
21. Daubechies, I., 1988. Orthonormal bases of compactly supported wavelets, Commun. Pure appl. Math., 41(7), 909-996.
22. Daubechies, I., 1992. Ten lectures on wavelets, in Proceedings of the CBMSNSF Regional Conference Series in Applied Mathematics, Vol. 61, SIAM.
23. Daubechies, I. & Sweldens, W., 1998. Factoring wavelet transforms into lifting steps, J. Fourier Anal. Appl., 4(3), 247-269.
24. DeKock, G., Armstrong, R., Siegfried,H.&Thomas, E., 2011. Geochronology of the Birim Supergroup of the West African craton in the Wa-bolé region of west-central Ghana: implications for the stratigraphic framework, J. Afr. Earth Sci., 59(1), 1-40.
25. Dongarra, J.J., Du Croz, J., Hammarling, S. & Hanson, R.J., 1988. An extended set of FORTRAN basic linear algebra subprograms, ACMTrans. Math. Softw., 14(1), 1-17.
26. Durrheim, R. & Nyblade, A., 2009. Africa array special volume: introduction, S. Afr. J. Geol., 112, 209-212.
27. Farquharson, C. & Oldenburg, D., 2004. A comparison of automatic techniques for estimating the regularization parameter in non-linear inverse problems, Geophys. J. Int., 156, 411-425.
28. Fatahalian, K.&Houston, M., 2008. A closer look at GPUs, Commun. ACM, 51(10), 50-57.
29. Fullagar, P., Pears, G., Hutton, D. & Thmpson, A., 2004. 3D gravity and aeromagnetic inversion for MVT lead-zinc exploration at pillara,Western Australia, Explor. Geophys., 35(2), 142-146.
30. Fullagar, P., Pears, G. & McMonnies, B., 2008. Constrained inversion of geologic surfaces-pushing the boundaries, Leading Edge, 27(1), 98-105.
31. Gallardo, L., 2007. Multiple cross-gradient joint inversion for geospectral imaging, Geophys. Res. Lett., 34, L19301, doi:10.1029/2007GL030409.
32. Gallardo, L. & Meju, M., 2011. Structure-coupled multiphysics imaging in geophysical sciences, Rev. Geophys., 49, 1-20.
33. Garcia-Abdeslem, J., 2000. Nonlinear 3-D inversion of gravity data over a sulfide ore body, Geofisica Internacional, 39(2), 179-188.
34. Garcia-Abdeslem, J., 2008. 3D forward and inverse modelling of total-field magnetic anomalies caused by a uniformly magnetized layer defined by a linear combination of 2D Gaussian functions, Geophysics, 73(1), L11- L18.
35. Garcia-Abdeslem, J., Espinosa-Cardea, J., Mungua-Orozco, L., Wong-Ortega, V. & Ramrez-Hernndez, J., 2001. Crustal structure from 2-D gravity and magnetic data modeling, magnetic power spectrum inversion, and seismotectonics in the Laguna SaladaBasin, northernBaja California, Mexico, Geofis. Int., 40(2), 67-87.
36. Gropp, W., Lusk, E. & Skjellum, A., 1994. Using MPI, Portable Parallel Programming with the Message-Passing Interface, MIT Press.
37. Guillen, A., Calcagno, P., Courrioux, G., Joly, A. & Ledru, P., 2008. Geological modelling from field data and geological knowledge, part II: modelling validation using gravity and magnetic data inversion, Phys. Earth planet. Inter., 171, 158-169.
38. Hansen, P., 2001. The L-curve and its use in the numerical treatment of inverse problems, in Advances in Computational Bioengineering, Computational Inverse Problems in Electrocardiology, Vol. 5, pp. 119-142, ed. Johnston, P., WIT Press.
39. Hastings, D., 1978. The southwest Ghana gravity survey, withWest African regional correlations and implications. Tech. Rep, Unpublished Report, Ghana Geological Survey.
40. Hastings, D., 1982. On the tectonics and metallogenesis of West Africa: a model incorporating new geophysical data, Geoexploration, 20, 295-327.
41. Hildebrand, A., Pilkington, M., Ortiz-Aleman, C., Chavez, R., Urrutia-Fucugauchi, J., Connors, M., Graniel-Castro, E. & Niehaus, D., 1998. Mapping Chicxulub crater structure with gravity and seismic reflection data, meteorites: flux with time and impact effects, Geol. Soc. Lond. Spec. Publ., 140, 155-176.
42. Hornby, P., Boschetti, F. & Horowitz, F.G., 1998. Analysis of potential field data in the wavelet domain, Geophys. J. Int., 137, 175-196.
43. Kirk, D.B. & Hwu, W.-M.W., 2010. Programming Massively Parallel Processors: A Hands-on Approach, Morgan Kaufmann.
44. Komatitsch, D., 2011. Fluid-solid coupling on a cluster of GPU graphics cards for seismic wave propagation, C. R. Acad. Sci., Ser. IIb Mec., 339, 125-135.
45. Komatitsch, D., G̈oddeke, D., Erlebacher, G. & Mich́ea, D., 2010. Modeling the propagation of elastic waves using spectral elements on a cluster of 192 GPUs, Comp. Sci. Res. Dev., 25(1-2), 75-82.
46. Koren, Z.,Mosegaard, K., Landa, E., Thore, P.&Tarantola, A., 1991. Monte Carlo estimation and resolution analysis of seismic background velocities, J. geophys. Res., 96(B12), 20 289-20 299.
47. Lajaunie, C., Courrioux, G. & Manuel, L., 1997. foliation fields and 3D cartography in geology: principles of a method based on potential interpolation, Math. Geol., 29(4), 571-584.
48. Li, Y. & Oldenburg, D.W., 2003. Fast inversion of large-scale magnetic data using wavelet transforms and a logarithmic barrier method, Geophys. J. Int., 152, 251-265.
49. Liu, G., Meng, X. & Chen, Z., 2012. 3D magnetic inversion based on probability tomography and its GPU implement, Comput. Geosci., 48, 86-92.
50. Loris, I., Nolet, G., Daubechies, I. & Dahlen, F., 2007. Tomographic inversion using l1-norm regularization of wavelet coefficients, Geophys. J. Int., 170, 359-370.
51. Loris, I., Douma, H., Nolet, G., Daubechies, I. & Regone, C., 2010. Nonlinear regularization techniques for seismic tomography, J. Comput. Phys., 229, 890-905.
52. Mallet, J.L., 1992. Discrete smooth interpolation in geometric modeling, Comput.-Aided Des., 24, 178-191.
53. Martin, R. & Komatitsch, D., 2009. An unsplit convolutional perfectly matched layer technique improved at grazing incidence for the viscoelastic wave equation, Geophys. J. Int., 179(1), 333-344.
54. Martin, R., Komatitsch, D. & Gedney, S.D., 2008. A variational formulation of a stabilized unsplit convolutional perfectly matched layer for the isotropic or anisotropic seismic wave equation, Comput.Model. Eng. Sci., 37(3), 274-304.
55. Martin, R., Komatitsch, D., Gedney, S.D. & Bruthiaux, E., 2010. A highorder time and space formulation of the unsplit perfectly matched layer for the seismic wave equation using Auxiliary Differential Equations (ADE-PML), Comput. Model. Eng. Sci., 56(1), 17-42.
56. Mauri, G.,Williams-Jones, G.&Saracco, G., 2010. Depth determinations of shallow hydrothermal systems by self-potential and multi-scale wavelet tomography, J. Volc. Geotherm. Res., 191(3-4), 233-244.
57. Mauri, G., Williams-Jones, G. & Saracco, G., 2011. MWTmat application of multiscale wavelet tomography on potential fields, Comput. Geosci., 37(11), 1825-1835.
58. Mich́ea,D.&Komatitsch,D., 2010. Accelerating a 3D finite-differencewave propagation code using GPU graphics cards, Geophys. J. Int., 182(1), 389-402.
59. Moorkamp, M., Jegen, M., Roberts, A. & Hobbs, R., 2010. Massively parallel forward modeling of scalar and tensor gravimetry data, Comput. Geosci., 36(36), 680-686.
60. Moorkamp, M., Heincke, B., Jegen, M., Roberts, A. & Hobbs, R., 2011. A framework for 3-D joint inversion of MT, gravity and seismic refraction data, Geophys. J. Int., 184(1), 477-493.
61. Morgan, J.V.,Warner,M.R., Collins, G.S., Grieve, R.A.F., Christeson, G.L., Gulick, S.P.S. & Barton, P.J., 2011. Full waveform tomographic imagesof the peak ring at the Chicxulub impact crater, J. geophys. Res., 116, B06303, doi:10.1029/2010JB008015.
62. Mosegaard, K. & Tarantola, A., 1995. Monte Carlo sampling of solutions to inverse problems, J. geophys. Res., 10(B7), 12 431-12 447.
63. Mosegaard, K. & Tarantola, A., 2002. Probabilistic Approach to Inverse Problems, Academic Press.
64. NVIDIA Corporation, 2011. NVIDIA CUDA toolkit documentation version 4. NVIDIA Corporation, Santa Clara, California, USA.
65. Ortiz-Aleman, C. & Martin, R., 2005. Two-phase oil-gas pipe flow imaging by simulated annealing, J. geophys. Eng., 2(1), 32-37.
66. Ortiz-Aleman, C. & Urrutia-Fucugauchi, J., 2010. Aeromagnetic anomaly modeling of central zone structure and magnetic sources in the Chicxulub crater, Phys. Earth planet. Inter., 179(3-4), 127-138.
67. Owens, J.D., Houston, M., Luebke, D.P., Green, S., Stone, J.E. & Phillips, J.C., 2008. GPU computing, Proc. IEEE, 96(5), 879-899.
68. Pacheco, P.S., 1997. Parallel Programming with MPI, Morgan Kaufmann.
69. Parker, R., 1994. Geophysical Inverse Theory, Princeton Univ. Press.
70. Pavlis, N.K., Holmes, S.A., Kenyon, S.C. & Factor, J.K., 2008. An earth gravitational model to degree 2160: EGM2008, Tech. Rep, General Assembly of the European Geosciences Union, Vienna, Austria, 2008 April 13-18.
71. Perrouty, S., 2012. Évolution Structurale de la Ceinture Minéralisée d'Ashanti, Sud-Ouest Ghana, PhD thesis, Universit́e de Toulouse III -Paul Sabatier, France.
72. Perrouty, S., Aillères, L., Jessell, M., Baratoux, L.&Bourassa, Y., 2012. Revised Eburnean geodynamic evolution of the gold-rich southern Ashanti Belt, Ghana, with new field and geophysical evidence of pre-Tarkwaian deformations, Precambrian Res., 204-205, 12-39.
73. Peter, D. et al., 2011. Forward and adjoint simulations of seismic wave propagation on fully unstructured hexahedral meshes, Geophys. J. Int., 186(2), 721-739.
74. Pilkington, M., 2009. 3D magnetic data-space inversion with sparseness constraints, Geophysics, 74(1), L7-L15.
75. Pilkington, M. & Hildebrand, A., 2000. Three-dimensional magnetic imaging of the Chicxulub crater, J. geophys. Res., 105, 23 479-23 491.
76. Pilkington, M., Hildebrand, A. & Ortiz-Aleman, C., 1994. Gravity and magnetic field modeling and structure of the Chicxulub crater, Mexico, J. geophys. Res., 99(E6), 13 147-13 162.
77. Pontoise, B., Bonvalot, S., Mascle, J. & Basile, C., 1990. Structure crustale de la marge transformante de Ĉ ote d'Ivoire-Ghana, C.R. Académie des Sciences de Paris, 310(2), 527-534.
78. Pratt, R.G., 1999. Seismic waveforminversion in the frequency domain, part 1: theory and verification in a physical scale model, Geophysics, 64(3), 888-901.
79. Pratt, R.G. & Chapman, C., 1992. Traveltime tomography in anisotropic media -II: application, Geophys. J. Int., 109, 20-37.
80. Pratt, R.G., Shin, C. & Hicks, G.J., 1998. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion, Geophysics, 133(2), 341-362.
81. Shahbahrami, A., 2012. Algorithms and architectures for2Ddiscretewavelet transform, J. Supercomput., 62, 1045-1064.
82. Sirgue, L. & Pratt, R.G., 2004. Efficient waveform inversion and imaging: a strategy for selecting temporal frequencies, Geophysics, 69(1), 231-248.
83. Stollnitz, E., deRose, T.&Salesin,D., 1995.Wavelets for computer graphics: a primer, part 1, IEEE Comput. Graph. Appl., 15(3), 76-86.
84. Talwani, M., Worzel, J.L. & Landisman, M., 1959. Rapid computation for two-dimensional bodies with application to the Mendocino submarine fracture zone, J. geophys. Res., 64, 49-59.
85. The International Association of Geodesy, 2012. The International Gravimetric Bureau, J. Geod., 86(10), 946-949.
86. Uytterhoeven, G., Roose, D.&Bultheel, A., 1997.Wavelet transforms using the lifting scheme. Tech. Rep, Report ITA-Wavelets- WP1.1 (Revised Version).
87. Vermeesch, P.M., Morgan, J.V., Christeson, G.L., Barton, P.J. & Surendra, A., 2009. Three-dimensional joint inversion of traveltime and gravity data across the Chicxulub impact crater, J. geophys. Res., 114, B02105, doi:10.1029/2008JB005776.
88. Wilson, G., Cuma,M.&Zhdanov, M., 2011. Massively parallel 3D inversion of gravity and gradiometry data, Preview, 152, 29-34.
89. Zhang, J., Wang, C., Shiz, Y., Cai, Y., Chi, W., Dreger, D., Cheng, W. & Yuan, Y., 2004. Three-dimensional crustal structure in central Taiwan from gravity inversion with a parallel genetic algorithm, Geophysics, 69(4), 917-924.
90. Zhdanov, M., 2009. New advances in regularized inversion of gravity and electromagnetic data, Geophys. Prospect., 57, 463-478.
91. Zhdanov, M., Liu, X.,Wilson, G. &Wan, L., 2011. Potential field migration for rapid imaging of gravity gradiometry data, Geophys. Prospect., 59(6), 1052-1071