A method is described for monitoring a stimulated reservoir volume (SRV) including receiving simulation parameters, performing 3D fully coupled quasi-static poro-elastic finite difference modeling using the simulation parameters, wherein the 3D fully coupled quasi-static poro-elastic finite difference modeling is based on a rescaling of solid rock and fluid flow density parameters and generates simulated temporal quasi-static stresses, and pore pressure. In addition, simulated stresses may be used for performing calculation of the 3D rotation of the simulated stresses to principal directions; performing calculation of the temporal 3D Mohr-Coulomb (MC) failure criteria from the calculated principal stresses and the simulated pore pressure for all or selected time steps; and displaying the computed temporal MC failure criteria results on a graphical display. The method may also be used in time-lapse monitoring of the reservoir for microseismic depletion delineation.