Comparison of Chemical-Component Transport in Naturally Fractured Reservoirs Using Dual-Porosity and Multiple-Interacting-Continua Models

SPE Journal ◽  
2020 ◽  
Vol 25 (04) ◽  
pp. 1964-1980
Author(s):  
Ali Al-Rudaini ◽  
Sebastian Geiger ◽  
Eric Mackay ◽  
Christine Maier ◽  
Jackson Pola
Keyword(s):  
Oil Recovery ◽  
Fractured Reservoirs ◽  
Transport Processes ◽  
Chemical Component ◽  
Naturally Fractured Reservoirs ◽  
Spontaneous Imbibition ◽  
Dual Porosity ◽  
Fine Scale ◽  
Naturally Fractured ◽  
The Matrix

Summary We propose a workflow to optimize the configuration of multiple-interacting-continua (MINC) models and overcome the limitations of the classical dual-porosity (DP) model when simulating chemical-component-transport processes during two-phase flow. Our new approach captures the evolution of the saturation and concentration fronts inside the matrix, which is key to design more effective chemical enhanced-oil-recovery (CEOR) projects in naturally fractured reservoirs. Our workflow is intuitive and derived from the simple concept that fine-scale single-porosity (SP) models capture fracture/matrix interaction accurately; it can hence be easily applied in any reservoir simulator with MINC capabilities. Results from the fine-scale SP model are translated into an equivalent MINC model that yields more accurate results compared with a classical DP model for oil recovery by spontaneous imbibition; for example, in a water-wet (WW) case, the root-mean-square error (RMSE) improves from 0.123 to 0.034. In general, improved simulation results can be obtained when selecting five or fewer shells in the MINC model. However, the actual number of shells is case specific. The largest improvement in accuracy is observed for cases where the matrix permeability is low and fracture/matrix transfer remains in a transient state for a prolonged time. The novelty of our approach is the simplicity of defining shells for a MINC model such that the chemical-component-transport process in naturally fractured reservoirs can be predicted more accurately, especially in cases where the matrix has low permeability. Hence, the improved MINC model is particularly suitable to model chemical-component transport, key to many CEOR processes, in (tight) fractured carbonates.

Verification and Proper Use of Water-Oil Transfer Function for Dual-Porosity and Dual-Permeability Reservoirs

2009 ◽  
Vol 12 (02) ◽  
pp. 189-199 ◽  
Author(s):  
Adetayo S. Balogun ◽  
Hossein Kazemi ◽  
Erdal Ozkan ◽  
Mohammed Al-kobaisi ◽  
Benjamin Ramirez
Keyword(s):  
Transfer Function ◽  
Oil Recovery ◽  
Transfer Functions ◽  
Fractured Reservoirs ◽  
Fracture Network ◽  
Naturally Fractured Reservoirs ◽  
Dual Porosity ◽  
Fluid Exchange ◽  
Matrix Block ◽  
Naturally Fractured

Summary Accurate calculation of multiphase fluid transfer between the fracture and matrix in naturally fractured reservoirs is a very crucial issue. In this paper, we will present the viability of the use of a simple transfer function to accurately account for fluid exchange resulting from capillary and gravity forces between fracture and matrix in dual-porosity and dual-permeability numerical models. With this approach, fracture- and matrix-flow calculations can be decoupled and solved sequentially, improving the speed and ease of computation. In fact, the transfer-function equations can be used easily to calculate the expected oil recovery from a matrix block of any dimension without the use of a simulator or oil-recovery correlations. The study was accomplished by conducting a 3-D fine-grid simulation of a typical matrix block and comparing the results with those obtained through the use of a single-node simple transfer function for a water-oil system. This study was similar to a previous study (Alkandari 2002) we had conducted for a 1D gas-oil system. The transfer functions of this paper are specifically for the sugar-cube idealization of a matrix block, which can be extended to simulation of a match-stick idealization in reservoir modeling. The basic data required are: matrix capillary-pressure curves, densities of the flowing fluids, and matrix block dimensions. Introduction Naturally fractured reservoirs contain a significant amount of the known petroleum hydrocarbons worldwide and, hence, are an important source of energy fuels. However, the oil recovery from these reservoirs has been rather low. For example, the Circle Ridge Field in Wind River Reservation, Wyoming, has been producing for 50 years, but the oil recovery is less than 15% (Golder Associates 2004). This low level of oil recovery points to the need for accurate reservoir characterization, realistic geological modeling, and accurate flow simulation of naturally fractured reservoirs to determine the locations of bypassed oil. Reservoir simulation is the most practical method of studying flow problems in porous media when dealing with heterogeneity and the simultaneous flow of different fluids. In modeling fractured systems, a dual-porosity or dual-permeability concept typically is used to idealize the reservoir on the global scale. In the dual-porosity concept, fluids transfer between the matrix and fractures in the grid-cells while flowing through the fracture network to the wellbore. Furthermore, the bulk of the fluids are stored in the matrix. On the other hand, in the dual-permeability concept, fluids flow through the fracture network and between matrix blocks. In both the dual-porosity and dual-permeability formulations, the fractures and matrices are linked by transfer functions. The transfer functions account for fluid exchanges between both media. To understand the details of this fluid exchange, an elaborate method is used in this study to model flow in a single matrix block with fractures as boundaries. Our goal is to develop a technique to produce accurate results for use in large-scale modeling work.

A Critical Review for Proper Use of Water/Oil/Gas Transfer Functions in Dual-Porosity Naturally Fractured Reservoirs: Part I

2009 ◽  
Vol 12 (02) ◽  
pp. 200-210 ◽  
Author(s):  
Benjamin Ramirez ◽  
Hossein Kazemi ◽  
Mohammed Al-kobaisi ◽  
Erdal Ozkan ◽  
Safian Atan
Keyword(s):  
Mass Transfer ◽  
Oil Recovery ◽  
Transfer Functions ◽  
Fractured Reservoirs ◽  
Gas Injection ◽  
Naturally Fractured Reservoirs ◽  
Dual Porosity ◽  
Fine Grid ◽  
Matrix Block ◽  
The Matrix

Summary Accurate calculation of multiphase-fluid transfer between the fracture and matrix in naturally fractured reservoirs is a crucial issue. In this paper, we will present the viability of the use of simple transfer functions to account accurately for fluid exchange resulting from capillary, gravity, and diffusion mass transfer for immiscible flow between fracture and matrix in dual-porosity numerical models. The transfer functions are designed for sugar-cube or match-stick idealizations of matrix blocks. The study relies on numerical experiments involving fine-grid simulation of oil recovery from a typical matrix block by water or gas in an adjacent fracture. The fine-grid results for water/oil and gas/oil systems were compared with results obtained with transfer functions. In both water and gas injection, the simulations emphasize the interaction of capillary and gravity forces to produce oil, depending on the wettability of the matrix. In gas injection, the thermodynamic phase equilibrium, aided by gravity/capillary interaction and, to a lesser extent, by molecular diffusion, is a major contributor to interphase mass transfer. For miscible flow, the fracture/matrix mass transfer is less complicated because there are no capillary forces associated with solvent and oil; nevertheless, gravity contrast between solvent in the fracture and oil in the matrix creates convective mass transfer and drainage of oil. Using the transfer functions presented in this paper, fracture- and matrix-flow calculations can be decoupled and solved sequentially--reducing the complexity of the computation. Furthermore, the transfer-function equations can be used independently to calculate oil recovery from a matrix block.

Characterization of Matrix Wettability and Mass Transfer From Matrix to Fractures in Carbonate Reservoirs

2012 ◽  
Author(s):  
Anuj Gupta
Keyword(s):  
Oil Recovery ◽  
History Matching ◽  
Oil And Gas ◽  
Fractured Reservoirs ◽  
Naturally Fractured Reservoirs ◽  
Carbonate Reservoirs ◽  
Shape Factors ◽  
Naturally Fractured ◽  
Fractured Carbonate ◽  
The Matrix

This paper presents results of an experimental investigation, supported by numerical analysis, to characterize oil recovery from fractured carbonate reservoirs. Imbibition recovery of oil is measured as a function of time for samples with varying wettability and shape factors. One of the objectives of this study is to verify the validity of exponential transfer function for matrix-fracture systems with varying wettability and flow-boundary conditions. Another objective is to establish the possibility of quantitatively determining the wettability of a sample based on history-matching of cumulative imbibition recovery and recovery rate data. The productivity of most carbonate oil and gas reservoirs is closely tied to the natural or stimulated fracture system present in the reservoir. Further, the recovery from naturally fractured reservoirs, in presence of aquifer drive or waterflooding is closely tied to the wettability of the matrix. The approach presented in this paper offers means to evaluate how recovery factor in a fractured system can be affected by wettability. A detailed understanding of rock-fluid interactions and wettability alterations at the fracturing face should help design improved strategies for exploiting naturally fractured carbonate reservoirs.

A Fully Implicit, Compositional, Parallel Simulator for IOR Processes in Fractured Reservoirs

SPE Journal ◽  
2007 ◽  
Vol 12 (03) ◽  
pp. 367-381 ◽  
Author(s):  
Reza Naimi-Tajdar ◽  
Choongyong Han ◽  
Kamy Sepehrnoori ◽  
Todd James Arbogast ◽  
Mark A. Miller
Keyword(s):  
Fractured Reservoirs ◽  
Naturally Fractured Reservoirs ◽  
Dual Porosity ◽  
Fractured Reservoir ◽  
Matrix Fracture ◽  
Naturally Fractured ◽  
Matrix Blocks ◽  
Dual Porosity Model ◽  
The Matrix ◽  
Porosity Model

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).

Pressure-Transient Tests and Flow Regimes in Fractured Reservoirs

2015 ◽  
Vol 18 (02) ◽  
pp. 187-204 ◽  
Author(s):  
Fikri Kuchuk ◽  
Denis Biryukov
Keyword(s):  
History Matching ◽  
Fractured Reservoirs ◽  
Flow Regimes ◽  
Naturally Fractured Reservoirs ◽  
Dual Porosity ◽  
Hydraulic Fractures ◽  
Fractured Reservoir ◽  
Pressure Transient ◽  
Matrix Block ◽  
Naturally Fractured

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