COLLOQUIUM: Seismic Imaging and Inversion Based on Spectral-Element and Adjoint Methods
Harnessing high-performance computers and accurate numerical methods to better constrain physical properties of Earth's interior is rapidly becoming one of the most important research topics in exploration and global seismology. We use spectral-element and adjoint methods to iteratively improve 3D subsurface images ranging from exploration to global scales. The spectral-element method, a high-order finite-element method with the advantage of a diagonal mass matrix, is used to accurately calculate three-component synthetic seismograms in complex 3D Earth models.
An adjoint method is used to numerically compute Frechet derivatives of a misfit function based on the interaction between the wavefield for a reference Earth model and a wavefield obtained by using time-reversed differences between data and synthetics at all receivers as simultaneous sources. In combination with gradient-based optimization methods, such as the L-BFGS method, we are able to iteratively improve 3D images of Earth's interior and gradually minimize discrepancies between observed and simulated seismograms. Various misfit functions may be chosen to quantify these discrepancies, such as cross-correlation traveltime differences, frequency-dependent phase and amplitude anomalies as well as full-waveform differences. Various physical properties of the Earth are constrained based on this method, such as elastic wavespeeds, radial anisotropy, shear attenuation and impedance contrasts. We apply this method to study seismic inverse problems at various scales, from exploration-scale full-waveform inversion to global-scale seismic tomography.
The Princeton Plasma Physics Laboratory 2012-2013 Colloquium Committee is comprised of the following people. Please feel free to contact them by e-mail regarding any possible speakers or topics for future colloquia.