Q-interface imaging using accumulative attenuation estimation
A high-resolution Q model is beneficial for more accurate attenuation compensation and preferable for gas-related interpretation. Given an accurate velocity model, viscoacoustic/viscoelastic full-waveform inversion ( Q-FWI) could reconstruct a high-resolution Q model, but it requires significant computational cost due to the iterative process of solving viscoacoustic/viscoelastic wave equations. We have proposed an efficient high-resolution Q-interface imaging method through the following steps. First, we estimate the attenuated traveltime via inversion of the dynamic match filter between synthetic acoustic and observed viscoacoustic prestack records. Second, we derive virtual Q reflectivities via piecewise linear regression on the attenuated traveltime estimation. Finally, by convolving a source wavelet on the virtual Q reflectivities, we generate the virtual Q reflection gathers and migrate them through reverse time migration (RTM) to image the Q interfaces. The Q-interface information is essentially derived by comparing the accumulative attenuation effects estimated from near-offset primary reflections arriving at the same receiver successively in time, and the high resolution is assured by the piecewise linear regression based on prior knowledge of the Q-interface number along the depth. The key insight of our method is to use accumulative attenuation effects to derive immediate effects of Q interfaces (virtual Q reflections) in the prestack data domain, which are readily applicable for Q-interface imaging through simple acoustic RTM. Numerical examples demonstrate that our method produces unprecedented high-resolution images of Q interfaces along the vertical direction with satisfying positioning and interpretable polarity.