Forward and Adjoint Modeling of Three Viscoacoustic Equations Based on Rheological Models Using DEVITO

Abstract

Disregarding the effect of seismic signal attenuation in gas and hydrocarbon zones leads to a final image with low resolution. However, a range of equations called viscoacoustics overcome such limitations. We use three secondorder equations based on Maxwell, Kelvin-Voigt (KV), and Standard Linear Solid (SLS) models. We analyze the dissipation and dispersion effects on each of them through seismograms. We also perform Reverse Time Migration (RTM) using the exact adjoint operators (Q-RTAM). All numerical experiments were implemented using Devito - a domain-specific language (DSL) and code generation framework to design highly optimized finite difference kernels for use in inversion methods.