No abstract

a b s t r a c tThe internal structure and petrophysical property distribution of fault zones are commonly exceedingly complex compared to the surrounding host rock from which they are derived. This in turn produces highly complex fluid flow patterns which affect petroleum migration and trapping as well as reservoir behavior during production and injection. Detailed rendering and forecasting of fluid flow inside fault zones require high-resolution, explicit models of fault zone structure and properties. A fundamental requirement for achieving this is the ability to create volumetric grids in which modeling of fault zone structures and properties can be performed. Answering this need, a method for generating volumetric fault zone grids which can be seamlessly integrated into existing standard reservoir modeling tools is presented. The algorithm has been tested on a wide range of fault configurations of varying complexity, providing flexible modeling grids which in turn can be populated with fault zone structures and properties.

Fault damage zones in porous sandstones commonly exhibit networks of deformation bands reflecting crushing and reorganization of grains associated with small-scale, localized displacement. Deformation bands introduce anisotropic, orderof-magnitude reduction of effective permeability, which will affect fluid flow in reservoir rocks. We here present a method for incorporating these features in industrial-type reservoir models. The method involves the use of a three-dimensional fault zone grid generation technique that allows property modeling on a discrete high-resolution fault zone grid without refining the entire reservoir model. Deformation band data from 106 outcrop scan lines of fault damage zones were classified into discrete fault facies defined according to deformation band density. The distributional pattern of fault facies in the data exhibits recurrent spatial relationships, which could be reproduced using truncated Gaussian simulation in the modeling process. The frequency distribution of deformation band density for each facies was analyzed, and average density values were assigned to each facies for calculating cell permeability. Permeability anisotropy was handled by approximating the relationship between deformation band densities in different directions based on published high-resolution fault zone maps and cross sections. Fluid-flow simulations were carried out on several damage zones models, and results were benchmarked against models with conventional fault rendering without damage zones. Simulation results show that flow paths, remaining oil distribution, and reservoir responses in models incorporating damage zones deviate from models employing conventional fault representation without damage zones, and these differences increase as deformation band permeability decreases.

Access to 3D descriptions of fault zone architectures and recent development of modeling techniques allowing explicit rendering of these features in reservoir models, provide a new tool for detailed implementation of fault zone properties. Our aim is to assess how explicit rendering of fault zone architecture and properties affects performance of fluid flow simulation models. The test models use a fault with a maximum 100 m displacement and a fault damage zone with petrophysical heterogeneity caused by the presence of deformation bands. The distribution pattern of deformation bands in fault damage zones is well-documented, which allows generation of realistic models. A multiscale modeling workflow is applied to incorporate these features into reservoir models. Model input parameters were modulated to provide a range of property distributions, and the interplay between the modeling parameters and reservoir performance was analyzed. The influence of deformation-band damage zone on reservoir performance in the presence of different fault core transmissibility-multipliers was investigated. Two configurations are considered: one in which the fault terminates inside the model domain, representing a case in which the fluid can flow around the fault, and one in which the fault dissects the entire model domain, representing a case in which the fluid is forced to cross the fault. We observed that the impact of deformation-band fault damage zone on reservoir performance changes when the fault core transmissibility multiplier is changed. Reservoir performance is insensitive to changing damage zone heterogeneity in a configuration in which flow can move around the fault. Where flow cannot bypass the fault, the influence of fault damage zone heterogeneity on reservoir performance is significant even when the fault core transmissibility multiplier is low.