TAU-P SEISMIC DATA FOR VISCOELASTIC MEDIA .2. LINEARIZED INVERSION

An approach to extraction of viscoelastic parameters from seismic data is implemented and successfully tested. Viscoelastic inversion is performed using adaptive damping factors to control the sensitivity of the viscoelastic parameters in relation to the tau-p seismic data. A priori information is incorporated through the damping factors as standard deviations of the data and of the viscoelastic model parameters. The stability of the inversion process is controlled by the variation of the damping factors as a function of the residual errors and parameter updates at each iteration. Tests on synthetic and real data show that P- and S-wave quality factors, Q(p) and Q(s), in addition to P- and S-wave velocities and density C(p), C(s) and rho, can be extracted from tau-p seismic information. Singular value decomposition analysis demonstrates that estimated Q(p) and Q(s) values are more affected by the presence of data inaccuracies and noise than are those of C(p) and rho. C(s) and Q(s) are not uniquely recovered due to the limited contribution of P-S converted waves. Knowledge of the viscoelastic parameters is of particular importance in accurately describing petrophysical properties of rocks and pore fluids existing in the subsurface; this is demonstrated with real data from the Gulf of Mexico.