Mesoscopic Wave-Induced Fluid Flow Effect Extraction by Using Frequency-Dependent Prestack Waveform Inversion

Patchy saturation of gas and brine within the porous rock can induce significant attenuation and velocity dispersion effects, which in turn have a profound impact on seismic data. Thus, there should be a close relationship between dispersion and reservoir parameters, such as saturation, porosity, and permeability. Hence, using seismic data to determine dispersion information has been the subject of intensive research in recent years. It can help quantitatively predicting gas saturation and monitoring CO<sub>2</sub> sequestration, a key strategy for mitigation of global warming. Here, a frequency-dependent prestack waveform inversion workflow was proposed to extract the wave-induced fluid flow (WIFF) effect from seismic data (WIFF prestack waveform inversion strategy). The proposed approach consists of three essential parts, comprising spectral decomposition, <inline-formula> <tex-math notation="LaTeX">$Q$ </tex-math></inline-formula>-compensated prestack waveform inversion, and frequency-dependent prestack waveform inversion. First, <inline-formula> <tex-math notation="LaTeX">$\ell _{1}$ </tex-math></inline-formula> norm constrained inverse spectral decomposition is exploited to provide time-frequency amplitude and phase spectra with high resolution and reliable accuracy. Then, constant-<inline-formula> <tex-math notation="LaTeX">$Q$ </tex-math></inline-formula> compensated prestack waveform inversion is used to generate relatively precise <inline-formula> <tex-math notation="LaTeX">$P$ </tex-math></inline-formula>- and <inline-formula> <tex-math notation="LaTeX">$S$ </tex-math></inline-formula>-wave velocities, density, and Fréchet derivative for the subsequent dispersion inversion. Finally, frequency-divided seismic data and the estimated parameters and Fréchet derivative are used to invert the <inline-formula> <tex-math notation="LaTeX">$P$ </tex-math></inline-formula>-wave velocity dispersion. A comparison between the inversion results and band-limited rock physics analysis results shows that the proposed method can provide a quantitative inversion of the velocity dispersion to some extent. The proposed strategy provides a possibility for quantifying dispersion and further estimating gas saturation.