Uncertainty analysis for FWI

Uncertainty analysis is important for seismic interpretation. Based on the framework of Bayesian, we can analyse different statistic parameters of our FWI result. However, directly sampling the posterior probability density function (pdf) is computationally intractable. In order to make this problem computationally tractable, in this work, we use Gaussian distribution approximation and low rank approximation to generate the posterior pdf. Simultaneous shots are also used to reduce the computational costs.