Pressure-Transient Behavior of Double-Porosity Reservoirs with Transient Interporosity Transfer with Fractal Matrix Blocks

SPE Journal ◽  
2020 ◽  
pp. 1-24
Author(s):  
Alex Valdes-Perez ◽  
Thomas A. Blasingame

Summary Double-porosity/naturally fractured reservoir models have traditionally been used to represent the flow and pressure behavior for highly fractured carbonate reservoirs. Given that unconventional reservoirs such as shale-oil/gas reservoirs might not be considered to be multiporosity media, the use of the traditional/classical “double-porosity” models might not be adequate (or appropriate). The recent development of anomalous diffusion models has opened the possibility of adapting double-porosity models to estimate reservoir (and related) parameters for unconventional reservoirs. The primary objective of this work is to develop and demonstrate analytical reservoir models that provide (possible) physical explanations for the anomalous diffusion phenomenon. The models considering anomalous diffusion in reservoirs with Euclidean shape are developed using a convolved (i.e., time-dependent) version of Darcy's law. The use of these models can yield a power-law (straight-line) behavior for the pressure and/or rate performance, similar to the fractal reservoir models. The main advantage of using anomalous diffusion models compared with models considering fractal geometry is the reduction from two parameters (i.e., the fractal dimension and the conductivity index) to only one parameter (i.e., the anomalous diffusion exponent). However, the anomalous diffusion exponent does not provide information regarding the geometry or spatial distribution of the reservoir properties. To provide an alternative explanation for the anomalous diffusion phenomenon in petroleum reservoirs, we have developed double-porosity models considering matrix blocks with fractal geometry and fracture networks with either radial or fractal fracture networks. The flows inside the matrix blocks and the fractal fracture network assume that Darcy’s law is valid in its space-dependent (fractal) form, whereas the classical version of Darcy’s law is assumed for the radial-fracture-network case. The transient interporosity transfer is modeled using the classical convolution schemes given in the literature. We have defined the matrix blocks to be “infinite-acting” to represent the nano/micropermeability of shale reservoir. For the system defined by a fractal fracture network and infinite-acting fractal matrix blocks, we have investigated the influence of the fractal parameters (both matrix and fracture network) in the pressure- and rate-transient performance behaviors. We have defined the flow periods that can be observed in these sorts of systems and we have developed analytical solutions for pressure-transient analysis. We demonstrate that the use of the convolved version of Darcy’s law results in a model very similar to the diffusivity equation for double-porosity systems (which incorporates transient interporosity flow). In performing this work, we establish the following observations/conclusions derived from our new solutions: We find that the assumption of a well producing at variable rate (time-dependent inner-boundary condition) has a more-significant effect on the pressure (and derivative) functions and obscures the effects of the properties of the reservoir. We demonstrate that the anomalous-diffusion-phenomena model proposed for unconventional reservoirs can be directly related to the multiporosity concept model. Pressure and pressure-derivative responses can be used in the diagnosis of flow periods and in the evaluation/estimation of reservoir parameters in unconventional reservoirs.

1985 ◽  
Vol 25 (05) ◽  
pp. 743-756 ◽  
Author(s):  
Dimitrie Bossie-Codreanu ◽  
Paul R. Bia ◽  
Jean-Claude Sabathier

Abstract This paper describes an approach to simulating the flow of water, oil, and gas in fully or partially fractured reservoirs with conventional black-oil models. This approach is based on the dual porosity concept and uses a conventional tridimensional, triphasic, black-oil model with minor modifications. The basic feature is an elementary volume of the fractured reservoir that is simulated by several model cells; the matrix is concentrated into one matrix cell and tee fractures into the adjacent fracture cells. Fracture cells offer a continuous path for fluid flows, while matrix cello are discontinuous ("checker board" display). The matrix-fracture flows are calculated directly by the model. Limitations and applications of this approximate approach are discussed and examples given. Introduction Fractured reservoir models were developed to simulate fluid flows in a system of continuous fractures of high permeability and low porosity that surround discontinuous, porous, oil-saturated matrix blocks of much lower permeability but higher porosity. The use of conventional models that permeability but higher porosity. The use of conventional models that actually simulate the fractures and matrix blocks is restricted to small systems composed of a limited number of matrix blocks. The common approach to simulating a full-field fractured reservoir is to consider a general flow within the fracture network and a local flow (exchange of fluids) between matrix blocks and fractures. This local flow is accounted for by the introduction of source or sink terms (transfer functions). In this formulation, the model is not directly predictive because the source term (transfer function) is, in fact, entered data and is derived from outside the model by one of the following approaches:analytical computation,empirical determination (laboratory experiments), ornumerical simulation of one or several matrix blocks on a conventional model. To derive these transfer functions, imposing some boundary conditions is necessary. Unfortunately, it is generally impossible to foresee all the conditions that will arise in a, matrix block and its surrounding fractures during its field life. It would be helpful, therefore, to have a model that is able to compute directly the local flows according to changing conditions. However, to have low computing times, it is necessary to use an approximate formulation and, thus, to adjust some parameters to match results that are externally (and more accurately) derived in a few basis, well-defined conditions. By current investigative techniques, only a very general description of the matrix blocks and fissures can be obtained, so our knowledge of local flows is very approximate. This paper presents a modeling procedure that is an approximate but helpful approach to the simulation of fractured reservoirs and requires a few, simple modifications of conventional black-oil mathematical models. Review of the Literature Numerous papers related to single- and multiphase flow in fractured porous media have been published over the last three decades. On the basis of data from fractured limestone and sand-stone reservoirs, fractured reservoirs are pictured as stacks of matrix blocks separated by fractures (Figs. 1 and 2). The fractured reservoirs with oil-saturated matrices usually are referred to as "double porosity" systems. Primary porosity is associated with matrix blocks, while secondary porosity is associated with fractures. The porosity of the matrices is generally much greater than that of the fractures, but permeability within fractures may be 100 and even over 10,000 times higher permeability within fractures may be 100 and even over 10,000 times higher than within the matrices. The main difference between flow in a fractured medium and flow in a conventional porous system is that, in a fractured medium, the interconnected fracture network provides the main path for fluid flow through the reservoir, while local flows (exchanges of fluids) occur between the discontinuous matrix blocks and the surrounding fractures. Matrix oil flows into the fractures, and the fractures carry the oil to the wellbore. For single-phase flow, Barenblatt et al constructed a formula based on the dual porosity approach. They consider the reservoir as two overlying continua, the matrices and the fractures. SPEJ p. 743


Lithosphere ◽  
2021 ◽  
Vol 2021 (Special 1) ◽  
Author(s):  
Jiazheng Liu ◽  
Xiaotong Liu ◽  
Hongzhang Zhu ◽  
Xiaofei Ma ◽  
Yuxue Zhang ◽  
...  

Abstract The gridless analytical and semianalytical methodologies can provide credible solutions for describing the well performance of the fracture networks in a homogeneous reservoir. Reservoir heterogeneity, however, is common in unconventional reservoirs, and the productivity can vary significantly along the horizontal wells drilled for producing such reservoirs. It is oversimplified to treat the entire reservoir matrix as homogeneous if there are regions with extremely nonuniform properties in the reservoir. However, the existing analytical and semianalytical methods can only model simple cases involving matrix heterogeneity, such as composite, layered, or compartmentalized reservoirs. A semianalytical methodology, which can model fracture networks in heterogeneous reservoirs, is still absent; in this study, we propose a decomposed fracture network model to fill this gap. We discretize a fractured reservoir into matrix blocks that are bounded by the fractures and/or the reservoir boundary and upscale the local properties to these blocks; therefore, a heterogeneous reservoir can be represented with these blocks that have nonuniform properties. To obtain a general flow equation to characterize the transient flow in the blocks that may exhibit different geometries, we approximate the contours of pressure with the contours of the depth of investigation (DOI) in each block. Additionally, the borders of each matrix block represent the fractures in the reservoir; thus, we can characterize the configurations of complex fracture networks by assembling all the borders of the matrix blocks. This proposed model is validated against a commercial software (Eclipse) on a multistage hydraulic fracture model and a fracture network model; both a homogeneous case and a heterogeneous case are examined in each of these two models. For the heterogeneous case, we assign different permeabilities to the matrix blocks in an attempt to characterize the reservoir heterogeneity. The calculation results demonstrate that our new model can accurately simulate the well performance even when there is a high degree of permeability heterogeneity in the reservoir. Besides, if there are high-permeability regions existing in the fractured reservoir, a BDF may be observed in the early production period, and formation linear flow may be indistinguishable in the early production period because of the influence of reservoir heterogeneity.


1985 ◽  
Vol 25 (01) ◽  
pp. 14-26 ◽  
Author(s):  
Karsten Pruess

Abstract A multiple interacting continua (MINC) method is presented, which is applicable for numerical simulation presented, which is applicable for numerical simulation of heat and multiphase fluid flow in multidimensional, fractured porous media. This method is a generalization of the double-porosity concept. The partitioning of the flow domain into computational volume elements is based on the criterion of approximate thermodynamic equilibrium at all times within each element. The thermodynamic conditions in the rock matrix are assumed to be controlled primarily by the distance from the fractures, which leads to the use of nested gridblocks. The MINC concept is implemented through the integral finite difference (IFD) method. No analytical approximations are made for coupling between the fracture and matrix continua. Instead, the transient flow of fluid and heat between matrix and fractures is treated by a numerical method. The geometric parameters needed in simulation are preprocessed from a specification of fracture spacings and apertures and geometry of the matrix blocks. The numerical implementation of the MINC method is verified by comparison with the analytical solution of Warren and Root. Illustrative applications are given for several geothermal reservoir engineering problems. Introduction In this paper, we present a numerical method for simulating transient nonisothermal, two-phase flow of water in fractured porous medium. The method is base on a generalization of a concept originally proposed by Barenblatt et al. and introduced into the petroleum literature by Warren and Root, Odeh, and others in the form of what has been termed the "double-porosity" model. The essence of this approach is that in a fractured porous medium, fractures are characterized by much porous medium, fractures are characterized by much larger diffusivities (and hence, much smaller response times) than the rock matrix. Therefore, the early system response is influenced by the matrix. In seeking to analytically solve such a system, all fractures were grouped into one continuum and all the matrix blocks into another, resulting in two interacting continua coupled through a mass transfer function determined by the size and shape of the blocks, as well as the local difference in potentials between the two continua. Later, Kazemi and Duguid and Lee incorporated the double-porosity concept into a numerical model. For a more detailed description of the concept and its application, see Refs. 6 through 8. Very little work has been done in investigating nonisothermal, two-phase fluid flow in fractured porous media. Moench and coworkers used the discrete fracture approach to study the behavior of fissured, vapor-dominated geothermal reservoirs. The purpose of our work is first to generalize the double-porosity concept into one of many interacting continua. We then incorporate the MINC model into a simulator for nonisothermal transport of a homogeneous two-phase fluid (water and steam) in multidimensional systems. Our approach is considerably broader in scope and more general than any previous models discussed in the literature. The MINC previous models discussed in the literature. The MINC method permits treatment of multiphase fluids with large and variable compressibility and allows for phase transitions with latent heat effects, as well as for coupling between fluid and heat flow. The transient interaction between matrix and fractures is treated in a realistic way. Although the model can permit alternative formulations for the equation of motion, we shall assume that, macroscopically, each continuum obeys Darcy's law; in particular, we shall use the "cubic law" for the flow of particular, we shall use the "cubic law" for the flow of fluids in fracture. While the methodology presented in this paper is generally applicable to multiphase compositional thermal systems, our illustrative calculations were restricted to geothermal reservoir problems. The numerical method chosen to implement the MINC concept is the IFD method. In this method, all thermophysical and thermodynamic properties are represented by averages over explicitly defined finite subdomains, while fluxes of mass or energy across surface segments are evaluated through finite difference approximations. An important aspect of this method is that the geometric quantities required to evaluate the conductance between two communicating volume elements are provided directly as input data rather than having them generated from data on nodal arrangements and nodal coordinates. Thus, a remarkable flexibility is attained by which one can allow a volume element in any one continuum to communicate with another element in its own or any other continuum. Inasmuch as the interaction between volume elements of different continua is handled as a geometric feature, the IFD methodology does not distinguish between the MINC method and the conventional porous-medium type approaches to modeling. porous-medium type approaches to modeling. SPEJ p. 14


SPE Journal ◽  
2013 ◽  
Vol 18 (05) ◽  
pp. 969-981 ◽  
Author(s):  
Mehmet A. Torcuk ◽  
Basak Kurtoglu ◽  
Najeeb Alharthy ◽  
Hossein Kazemi

Summary In this paper, we present a new method to model heterogeneity and flow channeling in petroleum reservoirs—especially reservoirs containing interconnected microfractures. The method is applicable to both conventional and unconventional reservoirs where the interconnected microfractures form the major flow path. The flow equations, which could include flow contributions from matrix blocks of various size, permeability, and porosities, are solved by the Laplace-transform analytical solutions and finite-difference numerical solutions. The accuracy of flow from and into nanodarcy matrix blocks is of great interest to those dealing with unconventional reservoirs; thus, matrix flow equations are solved by use of both pseudosteady-state (PSS) and unsteady state (USS) formulations and the results are compared. The matrix blocks can be of different size and properties within the representative elementary volume (REV) in the analytical solutions, and within each control volume (CV) in the numerical solutions. Although the analytical solutions were developed for slightly compressible rock/fluid linear systems, the numerical solutions are general and can be used for nonlinear, multiphase, multicomponent flow problems. The mathematical solutions were used to analyze the longterm and short-term performances of two separate wells in an unconventional reservoir. It is concluded that matrix contribution to flow is very slow in a typical low-permeability unconventional reservoir and much of the enhanced production is from the fluids contained in the microfractures rather than in the matrix. In addition to field applications, the mathematical formulations and solution methods are presented in a transparent fashion to allow easy usage of the techniques for reservoir and engineering applications.


2019 ◽  
Vol 129 ◽  
pp. 70-79 ◽  
Author(s):  
Yuhang Wang ◽  
Saman A. Aryana ◽  
Myron B. Allen

Sign in / Sign up

Export Citation Format

Share Document