To enable a stable and efficient simulation,
the coupling terms in the two-phase flow
and the deformation processes are handled
by a mixture approach of monolithic and staggered schemes.
Moreover a pressure-pressure scheme is applied to the two-phase flow model for convenience in stress calculation in porous media [36].
With the present assumption,
the changes in flow and deformation fields of the model are simulated by using the standard Galerkin finite elementmethod,
which is realized in an object oriented scientific tool,
Open-GeoSys (an open source finite element code) developed by the authors.
Finally, we present some simulation results and discussion.