1887

Abstract

Summary

It is present an efficient parallel algorithm for the inversion of 3-D gravity data, which goal is to estimate the depth of a sedimentary basin in which the density contrast varies parabolically with depth. The efficiency of the gravity inversion methods applied to the interpretation of sedimentary basins depends on the number of data and model parameters to be estimated, making it very poor when the number of parameters is very large. We present the simulation results with a synthetic model of a sedimentary basin inspired in a real situation, taking advantage of a parallel Levenberg-Marquardt algorithm implemented using both MPI and OpenMP. Lanczos bidiagonalization method has been used to obtain the solution for the linearized subproblem at each iteration. The idea of obtaining the solution of a large system of equations using the bidiagonalization procedure is quite useful in practical problems, and allows to implement selection methods for the optimal regularization parameter in an easy way, like the weighted generalized cross validation method, adopted in this work. The hybrid parallel implementation combined with Lanczos bidiagonalization allows us to achieve a significant reduction of the computational cost, which is otherwise very high due to the scale of the problem.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201412757
2015-06-01
2024-03-29
Loading full text...

Full text loading...

References

  1. Abedi, M., Golhami, A. and Fathianpour, N.
    [2013] Fast inversion of magnetic data using lanczos bidiagonalization method. Journal of Applied Geophysics, 90, 126–137.
    [Google Scholar]
  2. Chakravarthi, V., Raghuram, H.M. and Singh, S.
    [2002] 3D forward gravity modeling of density interfaces above which the density contrast varies continuously with depth. Computers and Geosciencies, 28, 53–57, doi:10.1016/S0098‑3004(01)00080‑2.
    https://doi.org/10.1016/S0098-3004(01)00080-2 [Google Scholar]
  3. Chakravarthi, V. and Sundararajan, N.
    [2007] 3D gravity inversion of basement relief - a depth-dependent density approach. Geophysics, 72, I23–I32, doi:10.1190/1.2431634.
    https://doi.org/10.1190/1.2431634 [Google Scholar]
  4. Chung, J., Nagy, J.G. and O’Leary, D.P.
    [2008] A weighted gcv method for lanczos hybrid regularization. Electronic Transactions on Numerical Analysis, 28, 149–167.
    [Google Scholar]
  5. Marquardt, D.
    [1963] An algorithm for least squares estimation of non-linear parameters. Journal of the society of Industrial and Applied Mathematics, 11, 431–441.
    [Google Scholar]
  6. Martins, C.M., Barbosa, V.C.F. and Silva, J.B.C.
    [2010] Simultaneous 3D depth-tobasement and density-contrast estimates using gravity data and depth control at few points. Geophysics, 75, I21–I28, doi:10.1190/1.3380225.
    https://doi.org/10.1190/1.3380225 [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201412757
Loading
/content/papers/10.3997/2214-4609.201412757
Loading

Data & Media loading...

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