An Evaluation Of The Improved Dual Porosity Model For The Simulation Of Gravity Effects In Naturally Fractured Reservoirs

1991 ◽  
Vol 30 (03) ◽  
Author(s):  
L.S.K. Fung ◽  
D.A. Collins
2014 ◽  
Vol 17 (01) ◽  
pp. 82-97 ◽  
Author(s):  
Fikri Kuchuk ◽  
Denis Biryukov

Summary Fractures are common features of many well-known reservoirs. Naturally fractured reservoirs contain fractures in igneous, metamorphic, and sedimentary formations. Faults in many naturally fractured carbonate reservoirs often have high-permeability zones, and are connected to many fractures with varying conductivities. Furthermore, in many naturally fractured reservoirs, faults and fractures can be discrete (i.e., not a connected-network fracture system). New semianalytical solutions are used to understand the pressure behavior of naturally fractured reservoirs containing a network of discrete and/or connected (continuous) finite- and infinite-conductivity fractures. We present an extensive literature review of the pressure-transient behavior of fractured reservoirs. First, we show that the Warren and Root (1963) dual-porosity model is a fictitious homogeneous porous medium because it does not contain any fractures. Second, by use of the new solutions, we show that for most naturally fractured reservoirs, the Warren and Root (1963) dual-porosity model is inappropriate and fundamentally incomplete for the interpretation of pressure-transient well tests because it does not capture the behavior of these reservoirs. We examined many field well tests published in the literature. With few exceptions, none of them shows the behavior of the Warren and Root (1963) dual-porosity model. These examples exhibit very diverse pressure behaviors of discretely and continuously fractured reservoirs. Unlike the single derivative shape of the Warren and Root (1963) model, the derivatives of these examples exhibit many different flow regimes depending on fracture distribution and on their intensity and conductivity. We show these flow regimes with our new model for discretely and continuously fractured reservoirs. Most well tests published in the literature do not exhibit the Warren and Root (1963) dual-porosity reservoir-model behavior. If we interpret them by use of this dual-porosity model, then the estimated permeability, skin factor, interporosity flow coefficient (λ), and storativity ratio (ω) will not represent the actual reservoir parameters.


SPE Journal ◽  
2007 ◽  
Vol 12 (03) ◽  
pp. 367-381 ◽  
Author(s):  
Reza Naimi-Tajdar ◽  
Choongyong Han ◽  
Kamy Sepehrnoori ◽  
Todd James Arbogast ◽  
Mark A. Miller

Summary Naturally fractured reservoirs contain a significant amount of the world oil reserves. A number of these reservoirs contain several billion barrels of oil. Accurate and efficient reservoir simulation of naturally fractured reservoirs is one of the most important, challenging, and computationally intensive problems in reservoir engineering. Parallel reservoir simulators developed for naturally fractured reservoirs can effectively address the computational problem. A new accurate parallel simulator for large-scale naturally fractured reservoirs, capable of modeling fluid flow in both rock matrix and fractures, has been developed. The simulator is a parallel, 3D, fully implicit, equation-of-state compositional model that solves very large, sparse linear systems arising from discretization of the governing partial differential equations. A generalized dual-porosity model, the multiple-interacting-continua (MINC), has been implemented in this simulator. The matrix blocks are discretized into subgrids in both horizontal and vertical directions to offer a more accurate transient flow description in matrix blocks. We believe this implementation has led to a unique and powerful reservoir simulator that can be used by small and large oil producers to help them in the design and prediction of complex gas and waterflooding processes on their desktops or a cluster of computers. Some features of this simulator, such as modeling both gas and water processes and the ability of 2D matrix subgridding are not available in any commercial simulator to the best of our knowledge. The code was developed on a cluster of processors, which has proven to be a very efficient and convenient resource for developing parallel programs. The results were successfully verified against analytical solutions and commercial simulators (ECLIPSE and GEM). Excellent results were achieved for a variety of reservoir case studies. Applications of this model for several IOR processes (including gas and water injection) are demonstrated. Results from using the simulator on a cluster of processors are also presented. Excellent speedup ratios were obtained. Introduction The dual-porosity model is one of the most widely used conceptual models for simulating naturally fractured reservoirs. In the dual-porosity model, two types of porosity are present in a rock volume: fracture and matrix. Matrix blocks are surrounded by fractures and the system is visualized as a set of stacked volumes, representing matrix blocks separated by fractures (Fig. 1). There is no communication between matrix blocks in this model, and the fracture network is continuous. Matrix blocks do communicate with the fractures that surround them. A mass balance for each of the media yields two continuity equations that are connected by matrix-fracture transfer functions which characterize fluid flow between matrix blocks and fractures. The performance of dual-porosity simulators is largely determined by the accuracy of this transfer function. The dual-porosity continuum approach was first proposed by Barenblatt et al. (1960) for a single-phase system. Later, Warren and Root (1963) used this approach to develop a pressure-transient analysis method for naturally fractured reservoirs. Kazemi et al. (1976) extended the Warren and Root method to multiphase flow using a 2D, two-phase, black-oil formulation. The two equations were then linked by means of a matrix-fracture transfer function. Since the publication of Kazemi et al. (1976), the dual-porosity approach has been widely used in the industry to develop field-scale reservoir simulation models for naturally fractured reservoir performance (Thomas et al. 1983; Gilman and Kazemi 1983; Dean and Lo 1988; Beckner et al. 1988; Rossen and Shen 1989). In simulating a fractured reservoir, we are faced with the fact that matrix blocks may contain well over 90% of the total oil reserve. The primary problem of oil recovery from a fractured reservoir is essentially that of extracting oil from these matrix blocks. Therefore it is crucial to understand the mechanisms that take place in matrix blocks and to simulate these processes within their container as accurately as possible. Discretizing the matrix blocks into subgrids or subdomains is a very good solution to accurately take into account transient and spatially nonlinear flow behavior in the matrix blocks. The resulting finite-difference equations are solved along with the fracture equations to calculate matrix-fracture transfer flow. The way that matrix blocks are discretized varies in the proposed models, but the objective is to accurately model pressure and saturation gradients in the matrix blocks (Saidi 1975; Gilman and Kazemi 1983; Gilman 1986; Pruess and Narasimhan 1985; Wu and Pruess 1988; Chen et al. 1987; Douglas et al. 1989; Beckner et al. 1991; Aldejain 1999).


SPE Journal ◽  
2013 ◽  
Vol 18 (04) ◽  
pp. 670-684 ◽  
Author(s):  
S.. Geiger ◽  
M.. Dentz ◽  
I.. Neuweiler

Summary A major part of the world's remaining oil reserves is in fractured carbonate reservoirs, which are dual-porosity (fracture-matrix) or multiporosity (fracture/vug/matrix) in nature. Fractured reservoirs suffer from poor recovery, high water cut, and generally low performance. They are modeled commonly by use of a dual-porosity approach, which assumes that the high-permeability fractures are mobile and low-permeability matrix is immobile. A single transfer function models the rate at which hydrocarbons migrate from the matrix into the fractures. As shown in many numerical, laboratory, and field experiments, a wide range of transfer rates occurs between the immobile matrix and mobile fractures. These arise, for example, from the different sizes of matrix blocks (yielding a distribution of shape factors), different porosity types, or the inhomogeneous distribution of saturations in the matrix blocks. Thus, accurate models are needed that capture all the transfer rates between immobile matrix and mobile fracture domains, particularly to predict late-time recovery more reliably when the water cut is already high. In this work, we propose a novel multi-rate mass-transfer (MRMT) model for two-phase flow, which accounts for viscous-dominated flow in the fracture domain and capillary flow in the matrix domain. It extends the classical (i.e., single-rate) dual-porosity model to allow us to simulate the wide range of transfer rates occurring in naturally fractured multiporosity rocks. We demonstrate, by use of numerical simulations of waterflooding in naturally fractured rock masses at the gridblock scale, that our MRMT model matches the observed recovery curves more accurately compared with the classical dual-porosity model. We further discuss how our multi-rate dual-porosity model can be parameterized in a predictive manner and how the model could be used to complement traditional commercial reservoir-simulation workflows.


2015 ◽  
Vol 18 (02) ◽  
pp. 187-204 ◽  
Author(s):  
Fikri Kuchuk ◽  
Denis Biryukov

Summary Fractures are common features in many well-known reservoirs. Naturally fractured reservoirs include fractured igneous, metamorphic, and sedimentary rocks (matrix). Faults in many naturally fractured carbonate reservoirs often have high-permeability zones, and are connected to numerous fractures that have varying conductivities. Furthermore, in many naturally fractured reservoirs, faults and fractures can be discrete (rather than connected-network dual-porosity systems). In this paper, we investigate the pressure-transient behavior of continuously and discretely naturally fractured reservoirs with semianalytical solutions. These fractured reservoirs can contain periodically or arbitrarily distributed finite- and/or infinite-conductivity fractures with different lengths and orientations. Unlike the single-derivative shape of the Warren and Root (1963) model, fractured reservoirs exhibit diverse pressure behaviors as well as more than 10 flow regimes. There are seven important factors that dominate the pressure-transient test as well as flow-regime behaviors of fractured reservoirs: (1) fractures intersect the wellbore parallel to its axis, with a dipping angle of 90° (vertical fractures), including hydraulic fractures; (2) fractures intersect the wellbore with dipping angles from 0° to less than 90°; (3) fractures are in the vicinity of the wellbore; (4) fractures have extremely high or low fracture and fault conductivities; (5) fractures have various sizes and distributions; (6) fractures have high and low matrix block permeabilities; and (7) fractures are damaged (skin zone) as a result of drilling and completion operations and fluids. All flow regimes associated with these factors are shown for a number of continuously and discretely fractured reservoirs with different well and fracture configurations. For a few cases, these flow regimes were compared with those from the field data. We performed history matching of the pressure-transient data generated from our discretely and continuously fractured reservoir models with the Warren and Root (1963) dual-porosity-type models, and it is shown that they yield incorrect reservoir parameters.


Sign in / Sign up

Export Citation Format

Share Document