An enriched–FEM technique for numerical simulation of interacting discontinuities in naturally fractured porous media

In this paper, an extended finite element method is presented for simulation of interaction between hydraulic fracturing and natural fractures in saturated porous media. The well-known u−p formulation is employed in order to obtain the fully coupled set of governing equations. Natural faults are modeled for both opening and closure modes where the fluid inflow and contact conditions are considered at the interface, respectively. The Darcy law is employed in conjunction with an aperture dependent permeability for the fracture channel to describe the interfacial inflow. The contact constraints of both the solid and fluid phases are imposed using the Penalty method. The Heaviside and modified ramp enrichment functions are introduced to model the required strong and weak discontinuities in the displacement and pressure fields, respectively. Furthermore, new enrichment functions are introduced to model the junction and/or intersection of the discontinuities. The system of nonlinear equations is solved using the Newton–Raphson iterative algorithm with an implicit time stepping scheme. Finally, several numerical examples are chosen to illustrate the robustness of the proposed formulation as well as to study the interaction of hydraulic fracturing with natural faults. •An enriched FEM for simulation of hydraulic fracturing process in fractured porous media.•Modeling HF using inflow continuity equation and momentum balance equation of porous domain.•A cohesive crack model for nonlinear fracturing process occurring ahead of hydro-fracture tip.•Applying the contact constraints along the natural fault faces for both solid and fluid phases.•Modeling the junction between discontinuities by enrichment functions for solid/fluid phases.