Traditional AVO forward modeling only considers the impact of reflection coefficients at the interface on seismic wave field amplitude and ignores various propagation effects. Introducing wave propagation effects including geometric spreading, transmission loss, attenuation into seismic wave propagation, multi-wave amplitude-preserved AVO forward modeling for horizontally layered media based on ray theory is proposed in this paper. We derived the multi-wave geometric spreading correction formulas for horizontally layered media in order to describe the geometric spreading effect of multi-wave propagation. Introducing the complex traveltime directly, we built the relationship between complex traveltime and quality factor without the help of complex velocity to describe the attenuation of viscoelastic media. Multi-wave transmission coefficients, obtained by solving the Zoeppritz equations directly, is used to describe the transmission loss. Numerical results show that the effects of geometric spreading, attenuation, and transmission loss on multi-wave amplitude varies with offset and multi-wave amplitude-preserved AVO forward modeling should consider the reconstructive effect of wave propagation on reflection amplitude.
Traditionally, fluid substitutions are often conducted on log data for calculating reservoir elastic properties with different pore fluids. Their corresponding seismic responses are computed by seismic forward modeling for direct gas reservoir identification. The workflow provides us with the information about reservoir and seismic but just at the well. For real reservoirs, the reservoir parameters such as porosity, clay content, and thickness vary with location. So the information from traditional fluid substitution just at the well is limited. By assuming a rock physics model linking the elastic properties to porosity and mineralogy, we conducted seismic forward modeling and AVO attributes computation on a three-layer earth model with varying porosity, clay content, and formation thickness. Then we analyzed the relations between AVO attributes at wet reservoirs and those at the same but gas reservoirs. We arrived at their linear relations within the assumption framework used in the forward modeling. Their linear relations make it possible to directly conduct fluid substitution on seismic AVO attributes. Finally, we applied these linear relations for fluid substitution on seismic data and identified gas reservoirs by the cross-plot between the AVO attributes from seismic data and those from seismic data after direct fluid substitution.