A two-scale approach for fluid flow in fractured porous media

A two-scale numerical model is developed for fluid flow in fractured, deforming porous media. At the microscale the flow in the cavity of a fracture is modelled as a viscous fluid. From the micromechanics of the flow in the cavity, coupling equations are derived for the momentum and the mass couplings to the equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. The finite element equations are derived for this two-scale approach and integrated over time. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures is independent from the underlying discretization. The resulting discrete equations are non-linear due to the non-linearity of the coupling terms. A consistent linearization is given for use within a Newton–Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach, and show that faults in a deforming porous medium can have a significant effect on the local as well as on the overall flow and deformation patterns. Copyright q 2006 John Wiley & Sons, Ltd