The behavior of naturally fractured reservoirs including fluid flow in matrix blocks

1987 ◽  
Vol 2 (2) ◽  
pp. 145-163 ◽  
Author(s):  
Chen Zhongxiang ◽  
You Jun
1965 ◽  
Vol 5 (01) ◽  
pp. 60-66 ◽  
Author(s):  
A.S. Odeh

Abstract A simplified model was employed to develop mathematically equations that describe the unsteady-state behavior of naturally fractured reservoirs. The analysis resulted in an equation of flow of radial symmetry whose solution, for the infinite case, is identical in form and function to that describing the unsteady-state behavior of homogeneous reservoirs. Accepting the assumed model, for all practical purposes one cannot distinguish between fractured and homogeneous reservoirs from pressure build-up and/or drawdown plots. Introduction The bulk of reservoir engineering research and techniques has been directed toward homogeneous reservoirs, whose physical characteristics, such as porosity and permeability, are considered, on the average, to be constant. However, many prolific reservoirs, especially in the Middle East, are naturally fractured. These reservoirs consist of two distinct elements, namely fractures and matrix, each of which contains its characteristic porosity and permeability. Because of this, the extension of conventional methods of reservoir engineering analysis to fractured reservoirs without mathematical justification could lead to results of uncertain value. The early reported work on artificially and naturally fractured reservoirs consists mainly of papers by Pollard, Freeman and Natanson, and Samara. The most familiar method is that of Pollard. A more recent paper by Warren and Root showed how the Pollard method could lead to erroneous results. Warren and Root analyzed a plausible two-dimensional model of fractured reservoirs. They concluded that a Horner-type pressure build-up plot of a well producing from a factured reservoir may be characterized by two parallel linear segments. These segments form the early and the late portions of the build-up plot and are connected by a transitional curve. In our analysis of pressure build-up and drawdown data obtained on several wells from various fractured reservoirs, two parallel straight lines were not observed. In fact, the build-up and drawdown plots were similar in shape to those obtained on homogeneous reservoirs. Fractured reservoirs, due to their complexity, could be represented by various mathematical models, none of which may be completely descriptive and satisfactory for all systems. This is so because the fractures and matrix blocks can be diverse in pattern, size, and geometry not only between one reservoir and another but also within a single reservoir. Therefore, one mathematical model may lead to a satisfactory solution in one case and fail in another. To understand the behavior of the pressure build-up and drawdown data that were studied, and to explain the shape of the resulting plots, a fractured reservoir model was employed and analyzed mathematically. The model is based on the following assumptions:1. The matrix blocks act like sources which feed the fractures with fluid;2. The net fluid movement toward the wellbore obtains only in the fractures; and3. The fractures' flow capacity and the degree of fracturing of the reservoir are uniform. By the degree of fracturing is meant the fractures' bulk volume per unit reservoir bulk volume. Assumption 3 does not stipulate that either the fractures or the matrix blocks should possess certain size, uniformity, geometric pattern, spacing, or direction. Moreover, this assumption of uniform flow capacity and degree of fracturing should be taken in the same general sense as one accepts uniform permeability and porosity assumptions in a homogeneous reservoir when deriving the unsteady-state fluid flow equation. Thus, the assumption may not be unreasonable, especially if one considers the evidence obtained from examining samples of fractured outcrops and reservoirs. Such samples show that the matrix usually consists of numerous blocks, all of which are small compared to the reservoir dimensions and well spacings. Therefore, the model could be described to represent a "homogeneously" fractured reservoir. SPEJ P. 60ˆ


2020 ◽  
Vol 134 (2) ◽  
pp. 399-434
Author(s):  
Daniel Lorng Yon Wong ◽  
Florian Doster ◽  
Sebastian Geiger ◽  
Eddy Francot ◽  
François Gouth

1985 ◽  
Vol 25 (03) ◽  
pp. 445-450 ◽  
Author(s):  
I. Ershaghi ◽  
R. Aflaki

Abstract This paper presents a critical analysis of some recently published papers on naturally fractured reservoirs. These published papers on naturally fractured reservoirs. These publications have pointed out that for a publications have pointed out that for a matrix-to-fracture-gradient flow regime, the transition portion of pressure test data on the semilog plot develops a portion of pressure test data on the semilog plot develops a slope one half that of the late-time data. We show that systems under pseudosteady state also may develop a 1:2 slope ratio. Examples from published case studies are included to show the significant errors associated with the characterization of a naturally fractured system by using the 1:2 slope concept for semicomplete well tests. Introduction Idealistic models of the dual-porosity type often have been recommended for interpretation of a well test in naturally fractured reservoirs. The evolutionary aspects of these models have been reviewed by several authors. Gradual availability of actual field tests and recent developments in analytical and numerical solution techniques have helped to create a better understanding of application and limitation of various proposed models. Two important observations should be made here. First, just as it is now recognized that classical work published by Warren and Root in 1963 was not the end of the line for interpretation of the behavior of naturally fractured systems, the present state of knowledge later may be considered the beginning of the technology. Second, parallel with the ongoing work by various investigators who progressively include more realistic assumptions in their progressively include more realistic assumptions in their analytical modeling, one needs to ponder the implication of these findings and point out the inappropriate impressions that such publications may precipitate in the mind of practicing engineers. practicing engineers. This paper is intended to scrutinize statements published in recent years about certain aspects of the anticipated transition period developed on the semilog plot of pressure-drawdown or pressure-buildup test data. pressure-drawdown or pressure-buildup test data. The Transition Period In the dual-porosity models published to date, a naturally fractured reservoir is assumed to follow the behavior of low-permeability and high-storage matrix blocks in communication with a network of high-permeability and low-storage fractures. The difference among the models has been the assumed geometry of the matrix blocks or the nature of flow between the matrix and the fracture. However, in all cases, it is agreed that a transition period develops that is strictly a function of the matrix period develops that is strictly a function of the matrix properties and matrix-fracture relationship. Fig. 1 shows properties and matrix-fracture relationship. Fig. 1 shows a typical semilog plot depicting the transition period and the parallel lines. Estimation of Warren and Root's proposed and to characterize a naturally fractured proposed and to characterize a naturally fractured system requires the development of the transition period. The Warren and Root model assumes a set of uniformly distributed matrix blocks. Furthermore, the flow from matrix to fracture is assumed to follow a pseudosteady-state regime. Under such conditions, in theory, this period is an S-shaped curve with a point of inflection. Uldrich and Ershaghi developed a technique to use the coordinates of this point of inflection for estimating and under conditions where either the early- or the late-time straight lines were not available. Kazemi and de Swann presented alternative approaches to represent naturally fractured reservoirs. They assumed a geometrical configuration consisting of layered matrix blocks separated by horizontal fractures. Their observation was that for such a system the transition period develops as a straight line with no inflection point. Bourdet and Gringarten identified a semilog straight line during the transition period for unsteady-state matrix-fracture flow. Recent work by Streltsova and Serra et al emphasized the transient nature of flow from matrix to fracture and pointed out the development of a unique slope ratio. These authors, later joined by Cinco-L. and Samaniego-V., stated that under a transient flow condition, the straight-line shape of the transition period develops a slope that is numerically one-half the slope of the parallel straight lines corresponding to the early- or late-time data. It was further pointed out that the transient flow model is a more realistic method of describing the matrix-fracture flow. As such, they implied that in the absence of wellbore-storage-free early-time data, or late-time data in the case of limited-duration tests, one may use the slope of the transition straight line and proceed with the estimation of the reservoir properties. Statement of the Problem The major questions that need to be addressed at this time are as follows. SPEJ P. 445


1998 ◽  
Author(s):  
J. Kamath ◽  
S.H. Lee ◽  
C.L. Jensen ◽  
W. Narr ◽  
H. Wu

SPE Journal ◽  
2013 ◽  
Vol 19 (02) ◽  
pp. 289-303 ◽  
Author(s):  
Ali Moinfar ◽  
Abdoljalil Varavei ◽  
Kamy Sepehrnoori ◽  
Russell T. Johns

Summary Many naturally fractured reservoirs around the world have depleted significantly, and improved-oil-recovery (IOR) processes are necessary for further development. Hence, the modeling of fractured reservoirs has received increased attention recently. Accurate modeling and simulation of naturally fractured reservoirs (NFRs) is still challenging because of permeability anisotropies and contrasts. Nonphysical abstractions inherent in conventional dual-porosity and dual-permeability models make them inadequate for solving different fluid-flow problems in fractured reservoirs. Also, recent technologies for discrete fracture modeling may suffer from large simulation run times, and the industry has not used such approaches widely, even though they give more-accurate representations of fractured reservoirs than dual-continuum models. We developed an embedded discrete fracture model (DFM) for an in-house compositional reservoir simulator that borrows the dual-medium concept from conventional dual-continuum models and also incorporates the effect of each fracture explicitly. The model is compatible with existing finite-difference reservoir simulators. In contrast to dual-continuum models, fractures have arbitrary orientations and can be oblique or vertical, honoring the complexity of a typical NFR. The accuracy of the embedded DFM is confirmed by comparing the results with the fine-grid, explicit-fracture simulations for a case study including orthogonal fractures and a case with a nonaligned fracture. We also perform a grid-sensitivity study to show the convergence of the method as the grid is refined. Our simulations indicate that to achieve accurate results, the embedded discrete fracture model may only require moderate mesh refinement around the fractures and hence offers a computationally efficient approach. Furthermore, examples of waterflooding, gas injection, and primary depletion are presented to demonstrate the performance and applicability of the developed method for simulating fluid flow in NFRs.


Sign in / Sign up

Export Citation Format

Share Document