Seismic numerical modeling technique provide most suitable way for the simulation of wave propagation through subsurface rocks and play an essential role in seismic interpretation, seismic inversion and in evaluating and designing a seismic survey. In order to delineate subsurface more precisely, different approaches are proposed for numerical modeling in which the solutions of Biot’s poroelastic equation are given by using reflectivity, finite element and finite difference algorithms. In current study, a finite difference forward modeling method is proposed to quantify the variation in wave properties during its propagation in a poroelastic media exhibiting mesoscopic heterogeneity due to the existence of three phase fluids. The Biot’s poroelastic wave equations are resolved and a generalized finite difference technique in time domain (FDTD) is used to examine seismic responses in a poroelastic media saturated with three phase fluids. So that the ability to detect geo-fluids through forward modeling technique can be improved. The accuracy of the simulation will lay a good foundation for the follow up work and will enhance our understandings about wave propagation through multi-phase poroelastic media and will also assist in detection and discrimination between fluids of different nature. Also, the proposed approach is of great significance for imaging and prediction of subsurface structures and distribution of fluids.
