Theoretical Investigation of the Transition From Spontaneous to Forced Imbibition

SPE Journal ◽  
2018 ◽  
Vol 24 (01) ◽  
pp. 215-229 ◽  
Author(s):  
Lichi Deng ◽  
Michael J. King

Summary Spontaneous and forced imbibition are recognized as important recovery mechanisms in naturally fractured reservoirs because the capillary force controls the movement of the fluid between the matrix and the fracture. For unconventional reservoirs, imbibition is also important because the capillary pressure is more dominant in these tighter formations, and a theoretical understanding of the flow mechanism for the imbibition process will benefit the understanding of important multiphase-flow phenomena such as waterblocking. In this paper, a new semianalytic method is presented to examine the interaction between spontaneous and forced imbibition and to quantitatively represent the transient imbibition process. The methodology solves the partial-differential equation (PDE) of unsteady-state immiscible, incompressible flow with arbitrary saturation-dependent functions using the normalized water flux concept, which is identical to the fractional-flow terminology used in the traditional Buckley-Leverett analysis. The result gives a universal inherent relationship between time, normalized water flux, saturation profile, and the ratio between cocurrent and total flux. The current analysis also develops a novel stability envelope outside of which the flow becomes unstable caused by strong capillary forces, and the characteristic dimensionless parameter shown in the envelope is derived from the intrinsic properties of the rock and fluid system, and it can describe the relative magnitude of capillary and viscous forces at the continuum scale. This dimensionless parameter is consistently applicable in both capillary-dominated and viscous-dominated flow conditions.

SPE Journal ◽  
2012 ◽  
Vol 17 (02) ◽  
pp. 340-351 ◽  
Author(s):  
E.. Ashoori ◽  
W.R.. R. Rossen

Summary Foam is a promising means of increasing sweep in miscible- and immiscible-gas enhanced oil recovery (EOR). Surfactant alternating gas (SAG) is a preferred method of injection. Numerous studies verify that the water relative permeability function krw(Sw) is unaffected by foam. Studies of foam have used a variety of krw functions. This paper shows a connection between the krw(Sw) function and SAG foam effectiveness that is independent of the details of how foam reduces gas mobility. For simplicity, we analyze SAG processes in the absence of mobile oil; success without oil is a precondition to success with oil, and our analysis also applies to a miscible-gas process with oil in 1D in the absence of dispersion. Fractional-flow methods have proved useful and accurate for modeling foam EOR processes. The success of SAG depends on total mobility at a point of tangency to the fractional-flow curve, which defines the shock front at the leading edge of the foam bank. One can determine total mobility directly from the coordinates of this point (Sw, fw) if the function krw(Sw) is known. Geometric constraints limit the region in the fractional-flow diagram in which this point of tangency can occur. For a given krw(Sw) function, this limits the mobility reduction achievable for any possible SAG process. We examine the implications of this limitation for different krw functions. These implications include the following. Increasing nonlinearity of the krw function is advantageous for SAG processes, regardless of how foam reduces gas mobility. SAG is inappropriate for naturally fractured reservoirs if straight-line relative permeabilities apply, even if extremely strong foam can be stabilized in fractures. It is important to measure krw(Sw) separately for any formation for which a SAG process is envisioned.


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ˆ


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


SPE Journal ◽  
2011 ◽  
Vol 16 (04) ◽  
pp. 795-811 ◽  
Author(s):  
A.. Jamili ◽  
G.P.. P. Willhite ◽  
D.W.. W. Green

Summary Gas injection in naturally fractured reservoirs maintains the reservoir pressure and increases oil recovery primarily by gravity drainage and to a lesser extent by mass transfer between the flowing gas in the fracture and the porous matrix. Although gravity drainage has been studied extensively, there has been limited research on mass-transfer mechanisms between the gas flowing in the fracture and fluids in the porous matrix. This paper presents a mathematical model that describes the mass transfer between a gas flowing in a fracture and a matrix block. The model accounts for diffusion and convection mechanisms in both gas and liquid phases in the porous matrix. The injected gas diffuses into the porous matrix through gas and liquid phases, causing the vaporization of oil in the porous matrix, which is transported by convection and diffusion to the gas flowing in the fracture. Compositions of equilibrium phases are computed using the Peng-Robinson EOS. The mathematical model was validated by comparing calculations to two sets of experimental data reported in the literature (Morel et. al. 1990; Le Romancer et. al. 1994), one involving nitrogen (N2) flow in the fracture and the second with carbon dioxide (CO2) flow. The matrix was a chalk. The resident fluid in the porous matrix was a mixture of methane and pentane. In the N2-diffusion experiment, liquid and vapor phases were initially present, while in the CO2 experiment, the matrix was saturated with liquid-hydrocarbon and water phases. Calculated results were compared with the experimental data, including recovery of each component, saturation profiles, and pressure gradient between matrix and fracture. Agreement was generally good. The simulation revealed the presence of countercurrent flow inside the block. Diffusion was the main mass-transfer mechanism between matrix and fracture during N2 injection. In the CO2 experiment, diffusion and convection were both important.


2009 ◽  
Vol 12 (02) ◽  
pp. 200-210 ◽  
Author(s):  
Benjamin Ramirez ◽  
Hossein Kazemi ◽  
Mohammed Al-kobaisi ◽  
Erdal Ozkan ◽  
Safian Atan

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.


Author(s):  
Samir Prasun ◽  
Andrew K. Wojtanowicz

Abstract Reliable predictions of well recovery are crucial for designing reservoir development. In the bottom-water naturally-fractured reservoirs (NFRs), comprising a network of distributed fracture “corridors,” spacing (and apertures) of the corridors varies throughout the reservoir. This makes oil well’s recovery a probabilistic variable as it depends upon uncertain well’s location in the network. The uncertainty is two-fold; it concerns well’s location within corridor network and well’s possible intersection with the nearest corridor. In any network’s location (with closely- or sparsely–spaced corridors), wells may intercept fracture corridors (fracture well) or go in-between two corridors in a matrix block (matrix-well). A simplified way of estimating well recovery is to ignore well’s location within corridor network and consider only probability and performance of fracture well and matrix well in a statistically-equivalent reservoir with uniform spacing and aperture equal to their expected values derived from their known statistics. Another (fully probabilistic) method considers the combined probabilities of the well’s location in the network and being a fracture well or matrix well. The study evaluates discrepancy between the two methods, explains its statistical nature, and demonstrates their implementation in a corridor-type NFR described in the literature. In the study, recovery process is simulated by coupling the inner (near-well) zone’s discrete single-porosity flow model with the outer zone Dual Porosity Dual Permeability (DPDP) simulator. The matrix well’s inner zone extends from the well to the nearest corridor and for the fracture well inner zone covers the corridor and adjacent matrix blocks. In the simulations, matrix and fracture-wells are operated at maximum rate constrained by minimum downhole flowing pressure and the surface handling limit. The study is performed using statistical data from a corridor-type NFR with power-law-distributed spacing size from 19 ft to 153 ft and corridor apertures varying from 8ft to 31ft correlated with the spacing. The simplified method gives recovery values ranging from 28% to 37%, and the single value of total recovery 33% — normalized by the matrix and corridor size fractions of the total reservoir area. Alternatively, the probabilistic method gives two separate distributions of the fracture and matrix wells’ recoveries that are weighted by their probability and converted to a single distribution of total recovery using statistical concept of weighted average. The probabilistic estimation gives higher values of recovery — from 32% to 38% with the expected value of 36.6%. Moreover, there is a considerable 30% probability of having recovery greater than 36.6%. A mathematical proof provides explanation why the probabilistic method gives recovery estimate greater than that from the simplified method. Another advantage of the method is the cumulative probability plot of well recovery that, in practical applications, would let operators make reservoir development decisions based upon the risk-benefit consideration.


Author(s):  
Anuj Gupta

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.


Author(s):  
J. Salvador Flores Mondragón ◽  
César Andres Bernal Huicochea ◽  
Luis Silvestre Zamudio Rivera ◽  
Eduardo Buenrostro González ◽  
Luis Manuel Perera Pérez ◽  
...  

In mature fields, low oil production, increased gas production and water fractional flow of low pressure reservoir combined with the mobility ratio between the gas and oil feed contacts the w/o induced by the oil extraction process are accentuated in naturally fractured reservoirs (NFR) -. It is common N2 injection for pressure maintenance and decline of oil production; however N2 causes channeling towards producing wells. In various fields - NFR, closure of oil production by this mechanism loss value make unprofitable oil production and surface facilities for handling demanding high volume gas and / or water or gas production out of specification. In volumes estimated residual oil trapped in areas invaded by the gas cap and in areas of lower conductivity can be recovered if it has clearly identified enhanced recovery processes. Previous efforts to this work results showed potential benefits in terms of increased oil production and the significant reduction of Gas-Oil Ratio (GOR) and the technical and economic feasibility of using this type of process with surfactants developed specifically for the conditions of this area The technology tested and evaluated under methodological process is based on new supramolecular complexes wettability modifier, corrosion inhibitor and able to generate stable foams under high pressure, temperature and salinity, which penetrate and invade the channels of high conductivity formation cause decreased flow of gas by reducing gas mobility. The product in the liquid phase diffuses into the channels of lower conductivity which cannot penetrate the foam, and by spontaneous imbibition mechanism resulting from the change of wettability of the rock surface and reducing interfacial tension and favors an increase the oil recovery factor in naturally fractured reservoirs. The application of a methodological process allowed the parameters measurement and evaluation of test results, visualizing future opportunities for the new chemicals. This project was approved after evaluation from a process of allocation of federal funds. With the purpose of defining the further steps in the search for the chemicals stability and risk mitigation stages of industrial upgrading for the complexity the NFR, the following discussion is presented. In order to accelerate the knowledge of new technologies and its deployment on the field, PEMEX has diversified the efforts, to achieve the principal goals regarding new technologies. This will provide greater ability to assess best practices and technologies. To evaluate the efforts of companies a performance assessment model was designed and apply from 2008, which takes into account the integral complexity of each technology to attend the specific challenges from an Asset and to be fair in comparing the results obtained for the particular design of the test. The aim of this paper is to describe the results and the methodology used for developing the performance evaluation and identifying the new opportunities in the state of the art of these tests.


SPE Journal ◽  
2020 ◽  
Vol 25 (04) ◽  
pp. 1964-1980
Author(s):  
Ali Al-Rudaini ◽  
Sebastian Geiger ◽  
Eric Mackay ◽  
Christine Maier ◽  
Jackson Pola

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.


SPE Journal ◽  
2012 ◽  
Vol 17 (02) ◽  
pp. 540-554 ◽  
Author(s):  
F.. Civan ◽  
M.L.. L. Rasmussen

Summary Methodology is presented and proved for determination of the best-estimate parameter values affecting the matrix/fracture-interface fluid transfer in naturally fractured reservoirs. Fracture/surface-hindered interface transfer of immiscible fluids is considered between matrix blocks and surrounding natural fractures. Improved matrix/fracture-transfer models are applied on the basis of presumed matrix-block shapes. Analytical solutions and the limiting isotropic-matrix long-time shape factors developed for special boundary conditions are used for interpretation of typical laboratory tests conducted using rectangular- and cylindrical-shaped rock samples. Workable equations and straight-line data-plotting schemes are developed for effective analysis and interpretation of laboratory data obtained from various-shaped oil-saturated reservoir-rock samples immersed into brine. Applications concerning the water/air and water/decane systems in laboratory core tests are also presented. The present approach allows rapid determination of the characteristic parameters of the matrix/fracture-transfer models for various-shaped matrix blocks, which are essential for prediction of petroleum recovery from naturally fractured reservoirs. The methodology is verified using various experimental data, and the values of the characteristic parameters (e.g., the average diffusion-coefficient and the interface-skin-mass-transfer coefficient) are determined. The Arrhenius (1889) equation is shown to represent the temperature dependency of these parameters effectively.


Sign in / Sign up

Export Citation Format

Share Document