Multiphase poroelastic modeling in semi-space for deformable reservoirs

Field observations suggest that while pressure depletion is a local process that occurs dominantly in reservoirs being exploited, it triggers a redistribution of effective stress in a more extensive domain that involves not just the adjacent strata, but the entire overburden and surrounding rocks out to considerable distances. In an infinite or semi-infinite domain, a petroleum reservoir can be simulated as a displacement discontinuity, permitting use of an efficient displacement discontinuity boundary element method to calculate surrounding stress and displacement changes. A 3D poroelastic finite element method is used to account for the local reservoir model. By coupling the displacement discontinuity and finite element models, a 3D poroelastic reservoir in an infinite or semi-infinite domain is simulated. The numerical model has been verified and numerical examples are given. Results show that the relative deformation properties of surrounding strata have significant effects on pressure drawdown and reservoir deformation during production, and thus should be taken into account.