The computational domain is modeled as a fully saturated poroelastic continuum. The continuum is discretized with finite elements with pore pressure and displacement degrees of freedom using coupled pore-pressure-thermal mechanical solid (CPTnnn) elements.

Figure 46.2: Typical Mesh of Reservoir in a Subsurface

The interface between the reservoir and overburden formation is modeled using surface-based contact with both displacement and pore-pressure degrees of freedom. The contact is modeled as bonded and impermeable with respect to displacement and pore-pressure degrees of freedom, respectively.

The following input modifies the real constants to model impermeable contact:

! REAL CONSTANT MODIFICATION (IMPERMEABLE CONTACT)
RMODIF,6,37,0
RMODIF,7,37,0
RMODIF,8,37,0
RMODIF,9,37,0
RMODIF,10,37,0
RMODIF,11,37,0