IMMISCIBLE DISPLACEMENT IN VERTICALLY FRACTURED RESERVOIRS

A dual-porosity model is defined for saturated, two-phase, compressible, immiscible flow in a vertically fractured reservoir or aquifer. This model allows detailed simulation of the matrix-fracture interaction as well as the matrix flow itself. This is accomplished by directly coupling the matrix and fracture systems along the vertical faces of the matrix blocks, incorporating gravitational effects directly, and simulating flow inside the block. Thus fluid segregation due to gravitational effects and heterogeneities can be simulated. We show that our model can be derived via homogenization techniques. The model (in incompressible form for simplicity of exposition) is then approximated by a computationally efficient finite difference scheme. Calculations are presented to show the convergence of the scheme and to indicate the behavior of the model as a function of several parameters.