Petaquake - Large-Scale Parallel Nonlinear Optimization for High Resolution 3D-Seismic Imaging
People
(Responsible)
Abstract
Large-Scale Parallel Nonlinear Optimization for High Resolution 3D-Seismic Imaging Current methods in global or local-scale seismic tomography rely on approximate descriptions of wave propagation with the result of severely limiting the resolution of tomographic images. However, to truly understand the dynamics of our planet, we need to be able to seismically map its deep structure at resolutions much higher than it is nowadays possible. Major geophysical questions that require high resolution 3D imaging at the planetary scale include a better understanding, e.g., of the nature of mantle plumes and sinking tectonic plates. At the regional scale, reliable seismic images are crucial for more accurate earthquake location and the compilation of seismic hazard maps. Recent advances in algorithms, software development, and high performance computing systems have resulted in PDE-based solvers that scale up to millions of variables, make use of thousands of processors, and accommodate complex multiple-physics. As partial differential equations (PDE) solvers also mature in the Earth Sciences, there is an increasing interest in solving nonlinear seismic inversion problems governed by PDE-based models. Larger computer architectures and new algorithms for optimization and wave propagation now provide the computational ability to address the geophysical issues mentioned above in a more rigorous way: namely, to abandon asymptotic ray-theory approximations in favor of time-dependent PDE-based models, and replace linearized inversions by truly nonlinear optimization. To achieve this goal, it will be necessary to combine recent developments in computational methods for nonlinear optimization and wave propagation, such as high-order finite element discretizations, local time-stepping, iterative methods, and inexact parallel interior-point methods. More specifically, the scientific goals of the project are: to develop parallel numerical methods for forward wave propagation and large-scale nonlinear optimization, to explore the performance of such methods on emerging petascale architectures and to develop a new generation of a seismic inversion code for 3D Earth imaging.
The project was supported by the High Performance and High Productivity Computing initiative (HP2C).
Additional information
Publications
- Basini P., Nissen-Meyer T., Boschi L., Casarotti E., Verbeke J., Schenk O., Giardini D. (2013) The influence of nonuniform ambient noise on crustal tomography in Europe, Geochemistry, Geophysics, Geosystems:1471-1492
- Huber J., Schenk O., Naumann U., Varnik E., Wächter A. (2012) Algorithmic Differentiation and Nonlinear Optimization for an Inverse Medium Problem. Combinatorial Scientific Computing. Chapman and Hall/CRC, 203-231
- Rietmann M., Messmer P., Nissen-Meyer T., Peter D., Basini P., Komatitsch D., Schenk O., Tromp J., Boschi L., Giardini D. (2012) Forward and adjoint simulations of seismic wave propagation on emerging large-scale GPU architectures. International Conference for High Performance Computing, Networking, Storage and Analysis (SC'12). International Conference for High Performance Computing, Networking, Storage and Analysis (SC'12). Salt Lake City, Utah, USA. 10-16 November 2012
- Rietmann M., Messmer P., Nissen-Meyer T., Peter D., Basini P., Komatitsch D., Schenk O., Tromp J., Boschi L., Giardini D. (2012) Forward and adjoint simulations of seismic wave propagation on emerging large-scale GPU architectures. International Conference for High Performance Computing, Networking, Storage and Analysis (SC'12). International Conference for High Performance Computing, Networking, Storage and Analysis (SC'12). Salt Lake City, Utah, USA. 10-16 November 2012
- Christen M., Schenk O., Cui Y. (2012) Patus for convenient high-performance stencils: Evaluation in earthquake simulations. International Conference for High Performance Computing, Networking, Storage and Analysis (SC'12). International Conference for High Performance Computing, Networking, Storage and Analysis. Salt Lake City, Utah, USA. 10-16 November 2012
- Schenk O., Sathe M., Ucar B., Sameh A. (2012) Towards A Scalable Hybrid Linear Solver Based On Combinatorial Algorithms. Combinatorial Scientific Computing. Chapman and Hall/CR, 96-127
- Christen M., Schenk O., Burkhart H. (2011) Automatic code generation and tuning for stencil kernels on modern shared memory architectures, Computer Science - Research and Development:205-210
- Grote M., Huber J., Schenk O. (2011) Inexact Interior-Point Methods for Time-Harmonic Inverse Medium Problems. 10th International Conference on Mathematical and Numerical Aspects of Waves. 10th International Conference on Mathematical and Numerical Aspects of Waves. Vancouver, Canada. July 25-29, 2011
- Christen M., Schenk O., Burkhart H. (2011) PATUS: A Code Generation and Autotuning Framework for Parallel Iterative Stencil Computations on Modern Microarchitectures. 2011 IEEE International Parallel Distributed Processing Symposium. IEEE International Parallel Distributed Processing Symposium. Sheraton Anchorage Hotel and Spa Anchorage, AK, USA. 6 May - 20 May 2011
- Arbenz P., Saad Y., Sameh A., Schenk O. (2011) Special issue on Parallel Matrix Algorithms and Applications (PMAA'10), Parallel Computing:731-732
- Curtis F. E., Schenk O., Wächter A. (2010) An Interior-Point Algorithm for Large-Scale Nonlinear Optimization with Inexact Step Computations, SIAM Journal on Scientific Computing:3447-3475
- Bollhöfer M., Grote M. J., Schenk O. (2009) Algebraic Multilevel Preconditioner for the Helmholtz Equation in Heterogeneous Media, SIAM Journal on Scientific Computing:3781-3805