210 likes | 457 Views
Spectral Element Method and GPU Computing for Seismic Imaging. Chaiwoot Boonyasiriwat May 1, 2010. Outline. Introduction to Seismic Imaging Spectral Element Method GPU Computing and Computational Steering Summary. 1. Seismic Imaging. Given a seismic data set.
E N D
Spectral Element Method andGPU Computing for Seismic Imaging Chaiwoot Boonyasiriwat May 1, 2010
Outline • Introduction to Seismic Imaging • Spectral Element Method • GPU Computing and Computational Steering • Summary 1
Seismic Imaging • Given a seismic data set • Find a model that best describes the given data or minimizes the misfit function where CSG Common shot gather m? 2
Seismic Imaging • Use iterative gradient-based methods m? 3
Seismic Imaging in 3D • Improve efficiency • Forward modeling: use better numerical methods • Improve convergence rate: preconditioning • Use many sources simultaneously 4
Forward Modeling • Important component of seismic imaging methods • Numerically solving a wave equation, e.g., acoustic wave equation • Numerical methods: • Finite difference • Spectral/pseudo-spectral • Finite element/spectral element 5
Outline • Introduction to Seismic Imaging • Spectral Element Method • GPU Computing and Computational Steering • Summary 6
Finite Element Formulation Acoustic wave equation Weighted residual method Weak form Linear combination of basis functions Matrix form where 9
Geophysical Applications of FEM • Global-scale • Exploration-scale 10
Outline • Introduction to Seismic Imaging • Spectral Element Method • GPU Computing and Computational Steering • Summary 11
Seismic Imaging on GPUs • High performance and large band-width • General-purpose computing on GPUs using CUDA • Advanced hardware (GPUs) + algorithm (multisource imaging) • Computational steering: computation + visualization 12
Performance of GPU vs CPU 1000 750 Peak GFLOP/s 500 250 0 13 Courtesy of NVIDIA
Memory Bandwidth of GPU vs CPU 120 100 80 Bandwidth GB/s 60 40 20 0 14 Courtesy of NVIDIA
Conventional vs Multisource Method Conventional RTM Image using 200 gathers 15
Conventional vs Multisource Method Multisource RTM Image using 20 supergathers 10x speedup 16
Computational Steering User Interface • Initial model • Line search method • Regularization • Iteration number • Perturb model manually (steering) Misfit value Iteration number 17
Tomogram 1 True model Tomogram 2 18
Summary • Seismic imaging using spectral element methods • Global-scale waveform inversion using SEM • 3D seismic imaging on GPU cluster (CUDA+MPI) • Computational steering on GPU cluster 19