Simulated Annealing in Seismic Inversion
Dr. Mrinal K. Sen
Institute for Geophysics
John A.& Katherine G. Jackson School of Geosciences
The University of Texas at Austin

Seismic inversion involves estimating elastic properties or some attributes of the subsurface rocks by iterative fitting of observed seismograms with the theoretical seismograms. For many applications, optimal fitting of data with model may become computationally intractable due to nonlinearity, ill-posedness and expensive forward modeling. I will discuss several efficient approaches to address these issues. The non-linear optimization problem is solved very efficiently using a method called simulated annealing (SA) that does not usually get trapped in a local minimum of the fitness function that measures the similarities or differences between observed and synthetic data. The SA approach can be applied to post-stack inversion for impedances, wavelet estimation, and pre-stack waveform inversion. It can also be used for estimation of uncertainties in the derived results. The efficiency can be increased further by a combination of SA with a local optimization. I will show several field data examples to demonstrate the applicability and usefulness of our algorithm.