Bayesian joint residual moveout and amplitude-variation-with-offset inversion

Residual moveout (RMO) can have a severe impact on seismic amplitude-variation-with-offset (AVO) analysis. It is, therefore, common practice to quality control and correct the processed seismic data for RMO before AVO analysis. However, a complicating factor is that AVO effects and tuning may result in up- or down-dipping events that are easily mistaken for events with RMO, e.g., AVO Class 2p response. Flattening these events will lead to an incorrect AVO estimation. We present a new Bayesian joint RMO and AVO inversion to estimate the RMO time shifts and the AVO intercept and gradient. The joint inversion corrects the seismic data based on the RMO and AVO prior models rather than explicitly assuming that the data should be flattened. The prior models will typically guide toward flat gathers but will also allow for up- and down-dipping events when these are possible within the prior model. The method is illustrated on synthetic and real seismic data. For cases wherein flat events are correct, the results are similar to sequential methods with RMO correction before AVO analysis, but in situations wherein dipping seismic reflectors may be misinterpreted as events with RMO, the joint inversion provides better results by evaluating RMO and AVO simultaneously. The inversion results include the posterior covariance matrix, which represents uncertainties for the AVO intercept and gradient, the RMO time shifts, and the correlations between these at all samples. The uncertainty of the RMO time shift varies within seismic gathers and depends on the seismic quality and how clear and well determined the seismic events are. The uncertainty of the RMO is lower for clear and continuous events but increases in the zones with weaker and noisy events. The uncertainty of the RMO time shifts has a low impact on the uncertainty of the AVO intercept but increases the uncertainty of the AVO gradient significantly.