Coupled Effect of Imbibition Capillary Pressure and Matrix-Fracture Transfer on Oil Recovery from Dual-Permeability Reservoirs

2021 ◽  
Author(s):  
Aymen AlRamadhan ◽  
Yildiray Cinar ◽  
Arshad Hussain ◽  
Nader BuKhamseen

Abstract This paper presents a numerical study to examine how the interplay between the matrix imbibition capillary pressure (Pci) and matrix-fracture transfer affects oil recovery from naturally-fractured reservoirs under waterflooding. We use a dual-porosity, dual-permeability (DPDP) finite difference simulator to investigate the impact of uncertainties in Pci on the waterflood recovery behavior and matrix-fracture transfer. A comprehensive assessment of the factors that control the matrix-fracture transfer, namely Pci, gravity forces, shape factor and fracture-matrix permeabilities is presented. We examine how the use of Pci curves in reservoir simulation can affect the recovery assessment. We present two conceptual scenarios to demonstrate the impact of spontaneous and forced imbibition on the flood-front movement, waterflood recovery processes, and ultimate recovery in the DPDP reservoir systems of varying reservoir quality. The results demonstrate that the inclusion of Pci in reservoir simulation delays the breakthrough time due to a higher displacement efficiency. The study reveals that the matrix-fracture transfer is mainly controlled by the fracture surface area, fracture permeability, shape factor, and the uncertainty in Pci. We underline a discrepancy among various shape factors proposed in the literature due to three main factors: (1) the variations in matrix-block geometries considered, (2) how the physics of imbibition forces that control the multiphase fluid transfer is captured, and (3) how the assumption of pseudo steady-state flow is addressed.

SPE Journal ◽  
2019 ◽  
Vol 24 (06) ◽  
pp. 2653-2670 ◽  
Author(s):  
Didier–Yu Ding

Summary Unconventional shale–gas and tight oil reservoirs are commonly naturally fractured, and developing these kinds of reservoirs requires stimulation by means of hydraulic fracturing to create conductive fluid–flow paths through open–fracture networks for practical exploitation. The presence of the multiscale–fracture network, including hydraulic fractures, stimulated and nonstimulated natural fractures, and microfractures, increases the complexity of the reservoir simulation. The matrix–block sizes are not uniform and can vary in a very wide range, from several tens of centimeters to meters. In such a reservoir, the matrix provides most of the pore volume for storage but makes only a small contribution to the global flow; the fracture supplies the flow, but with negligible contributions to reservoir porosity. The hydrocarbon is mainly produced from matrix/fracture interaction. So, it is essential to accurately model the matrix/fracture transfers with a reservoir simulator. For the fluid–flow simulation in shale–gas and tight oil reservoirs, dual–porosity models are widely used. In a commonly used dual–porosity–reservoir simulator, fractures are homogenized from a discrete–fracture network, and a shape factor based on the homogenized–matrix–block size is applied to model the matrix/fracture transfer. Even for the embedded discrete–fracture model (EDFM), the matrix/fracture interaction is also commonly modeled using the dual–porosity concept with a constant shape factor (or matrix/fracture transmissibility). However, in real cases, the discrete–fracture networks are very complex and nonuniformly distributed. It is difficult to determine an equivalent shape factor to compute matrix/fracture transfer in a multiple–block system. So, a dual–porosity approach might not be accurate for the simulation of shale-gas and tight oil reservoirs because of the presence of complex multiscale–fracture networks. In this paper, we study the multiple–interacting–continua (MINC) method for flow modeling in fractured reservoirs. MINC is commonly considered as an improvement of the dual–porosity model. However, a standard MINC approach, using transmissibilities derived from the MINC proximity function, cannot always correctly handle the matrix/fracture transfers when the matrix–block sizes are not uniformly distributed. To overcome this insufficiency, some new approaches for the MINC subdivision and the transmissibility computations are presented in this paper. Several examples are presented to show that using the new approaches significantly improves the dual–porosity model and the standard MINC method for nonuniform–block–size distributions.


2002 ◽  
Vol 5 (02) ◽  
pp. 154-162 ◽  
Author(s):  
S. Sarda ◽  
L. Jeannin ◽  
R. Basquet ◽  
B. Bourbiaux

Summary Advanced characterization methodology and software are now able to provide realistic pictures of fracture networks. However, these pictures must be validated against dynamic data like flowmeter, well-test, interference-test, or production data and calibrated in terms of hydraulic properties. This calibration and validation step is based on the simulation of those dynamic tests. What has to be overcome is the challenge of both accurately representing large and complex fracture networks and simulating matrix/ fracture exchanges with a minimum number of gridblocks. This paper presents an efficient, patented solution to tackle this problem. First, a method derived from the well-known dual-porosity concept is presented. The approach consists of developing an optimized, explicit representation of the fractured medium and specific treatments of matrix/fracture exchanges and matrix/matrix flows. In this approach, matrix blocks of different volumes and shapes are associated with each fracture cell depending on the local geometry of the surrounding fractures. The matrix-block geometry is determined with a rapid image-processing algorithm. The great advantage of this approach is that it can simulate local matrix/fracture exchanges on large fractured media in a much faster and more appropriate way. Indeed, the simulation can be carried out with a much smaller number of cells compared to a fully explicit discretization of both matrix and fracture media. The proposed approach presents other advantages owing to its great flexibility. Indeed, it accurately handles the cases in which flows are not controlled by fractures alone; either the fracture network may be not hydraulically connected from one well to another, or the matrix may have a high permeability in some places. Finally, well-test cases demonstrate the reliability of the method and its range of application. Introduction In recent years, numerous research programs have been focusing on the topic of fractured reservoirs. Major advances were made, and oil companies now benefit from efficient methodologies, tools, and software for fractured reservoir studies. Nowadays, a study of a fractured reservoir, from fracture detection to full-field simulation, includes the following main steps: geological fracture characterization, hydraulic characterization of fractures, upscaling of fracture properties, and fractured reservoir simulation. Research on fractured reservoir simulation has a long history. In the early 1960s, Barenblatt and Zheltov1 first introduced the dual-porosity concept, followed by Warren and Root,2 who proposed a simplified representation of fracture networks to be used in dual-porosity simulators. Based on this concept, reservoir simulators3 are now able to correctly reproduce the main driving mechanisms occurring in fractured reservoirs, such as water imbibition, gas/oil and water/oil gravity drainage, molecular diffusion, and convection in fractures. Even single-medium simulators can perform fractured reservoir simulation when adequate pseudocapillary pressure curves and pseudorelative permeability curves can be input. Indeed, except for particular cases such as thermal recovery processes, full-field simulation of fractured reservoirs is no longer a problem. Geological characterization of fractures progressed considerably in the 1990s. The challenge was to analyze and integrate all the available fracture data to provide a reliable description of the fracture network both at field scale and at local reservoir cell scale. Tools have been developed for merging seismic, borehole imaging, lithological, and outcrop data together with the help of geological and geomechanical rules.3 These tools benefited from the progress of seismic acquisition and borehole imaging. Indeed, accurate seismic data lead to reliable models of large-scale fracture networks, and borehole imaging gives the actual fracture description along the wells, which enables a reliable statistical determination of fracture attributes. Finally, these tools provide realistic pictures of fracture networks. They are applied successfully in numerous fractured-reservoir studies. The upscaling of fracture properties is the problem of translating the geological description of fracture networks into reservoir simulation parameters. Two approaches are possible. In the first one, the fractured reservoir is considered as a very heterogeneous matrix reservoir; therefore, one applies the classical techniques available for heterogeneous single-medium upscaling. The second approach is based on the dual-porosity concept and consists of upscaling the matrix and the fracture separately. Based on this second approach, methodologies and software were developed in the 1990s to calculate equivalent fracture parameters with respect to the dual-porosity concept (i.e., a fracture-permeability tensor with main flow directions and anisotropy and a shape factor that controls the matrix/fracture exchange kinetics3–5). For a given reservoir grid cell, the upscaling procedures consist of generating the corresponding 3D discrete fracture network and computing the equivalent parameters from this network. In particular, the permeability tensor is computed from the results of steady-state flow simulations in the discrete fracture network alone (without the matrix).


2010 ◽  
Vol 13 (01) ◽  
pp. 109-117 ◽  
Author(s):  
V.. Er ◽  
T.. Babadagli

Summary CO2 injection has been applied in naturally fractured reservoirs (NFRs) for the purpose of enhanced oil recovery (i.e., the Wey-burn and Midale fields, Canada; the Wasson and Slaughter fields, USA; and the Bati Raman field, Turkey). The matrix part of these types of reservoirs could potentially be a good storage medium as well. Understanding the matrix/fracture interaction during this process and the dynamics of the flow in this dual-porosity system requires visual analyses. We mimicked fully miscible CO2 injection in NFRs using 2D models with a single fracture and oil (solute)/hydrocarbon solvent pairs. The focus was on the visual pore-scale analysis of miscibility interaction, breakthrough of solvent through fracture, transfer between matrix and fracture, and the dynamics of miscible displacement inside the matrix. First, matrix/fracture interaction was studied intensively using 2D glass-bead models experimentally. The model was prepared using acrylic sheets and glass beads saturated with oil as a porous medium while a narrow gap of 1-mm size containing filter paper served as a fracture. The first contact miscible solvent (pentane) was injected into the fracture, and the flow distribution was monitored using an image-acquisition and -processing system. The produced solvent and solute were continuously analyzed for compositional study. The effects of several parameters, such as flow rate, viscosity ratio (oil/solvent), and gravity, were studied. Next, the process was modeled numerically using a commercial compositional simulator, and the saturation distribution in the matrix was matched to experimental data. The key parameters in the matching process were the effective diffusion coefficients and the longitudinal and the transverse dispersivities. The diffusion coefficients were specified for each fluid, and dispersivities were assigned into gridblocks separately for the fracture and the matrix.


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.


2009 ◽  
Vol 13 (01) ◽  
pp. 44-55 ◽  
Author(s):  
Hamidreza Salimi ◽  
Johannes Bruining

Summary Most simulations of waterflooding in fractured media are based on the Warren and Root (WR) approach, which uses an empirical transfer function between the fracture and matrix block. We use homogenization to obtain an improved flow model in fractured media, leading to an integro-differential equation; also called the boundary-condition (BC) approach. We formulate a well-posed numerical 3D model for the BC approach. This paper derives this numerical model to solve full 3D integro-differential equations in a field reservoir simulation. We compare the results of the upscaled model with ECLIPSETM results. For the interpretation, it is useful to define three dimensionless parameters that characterize the oil production in fractured media. The most important of these parameters is a Peclet number, defined as the ratio between the time required to imbibe water into the matrix block and the travel time of water in the fracture system. The results of the WR approach and the BC approach are in good agreement when the travel time is of the same order of magnitude as the imbibition time. However, if the travel time is shorter or longer than the imbibition time, the approaches give different results. The BC approach allows the use of transfer functions based on fundamental principles (e.g., the use of a rate-dependent capillary pressure function). When implemented, it can be used to improve recovery predictions for waterflooded fractured reservoirs.


2021 ◽  
Vol 137 (3) ◽  
pp. 739-767
Author(s):  
Abdul Salam Abd ◽  
Na Zhang ◽  
Ahmad S. Abushaikha

AbstractCapillary dominated flow or imbibition—whether spontaneous or forced—is an important physical phenomena in understanding the behavior of naturally fractured water-driven reservoirs (NFR’s). When the water flows through the fractures, it imbibes into the matrix and pushes the oil out of the pores due to the difference in the capillary pressure. In this paper, we focus on modeling and quantifying the oil recovered from NFR’s through the imbibition processes using a novel fully implicit mimetic finite difference (MFD) approach coupled with discrete fracture/discrete matrix (DFDM) technique. The investigation is carried out in the light of different wetting states of the porous media (i.e., varying capillary pressure curves) and a full tensor representation of the permeability. The produced results proved the MFD to be robust in preserving the physics of the problem, and accurately mapping the flow path in the investigated domains. The wetting state of the rock affects greatly the oil recovery factors along with the orientation of the fractures and the principal direction of the permeability tensor. We can conclude that our novel MFD method can handle the fluid flow problems in discrete-fractured reservoirs. Future works will be focused on the extension of MFD method to more complex multi-physics simulations.


Author(s):  
Mahdi Abbasi ◽  
Alireza Kazemi ◽  
Mohammad Sharifi

Fractured reservoirs contain most of the oil in the world’s reserves. The existence of two systems of matrix and fracture with completely different characteristics has caused the modeling of the mechanisms of fractured reservoirs to be more complex than conventional ones. Modeling of this type of reservoirs is possible using two methods of single and dual porosity model. Modeling via single porosity scheme is very time-consuming as it takes into account huge matrix blocks (low permeability and high porosity) and small fractures (high permeability and low porosity) alongside each other explicitly. The dual porosity model, however, attempts to solve this problem using the concept of shape factor, which is defined as the amount of fluid transferred from the matrix to the fracture. The shape factor coefficients expressed so far have been derived via simplifying assumptions which keep them away from real conditions prevailing in fractured reservoirs. In this paper, shape factor is calculated more realistically with consideration of the quadratic pressure gradient in the diffusivity equation, the heterogeneity of the matrix block and the change of the rock properties by pressure change. For these three cases, the analytical modeling of the flow of fluid from the matrix to the fracture system has been discussed and its results with previous models have been compared. In addition, the dependence of shape factor on the stated parameters was evaluated and in order to validate the results of the proposed analytical model, its results were compared with the results of a commercial simulator. Investigating the shape factor with the assumptions about the physics of the fractured reservoirs will improve our understanding of the fluid transfer between the matrix and the fracture, and this capability will allow numerical and commercial simulators to predict the behavior of fractured reservoirs more accurately.


2008 ◽  
Vol 11 (03) ◽  
pp. 633-640 ◽  
Author(s):  
Martin Stoll ◽  
Jan Hofman ◽  
Dick J. Ligthelm ◽  
Marinus J. Faber ◽  
Paul van den Hoek

Summary Densely-fractured oil-wet carbonate fields pose a true challenge for oil recovery that traditional primary and secondary processes fail to meet. The difficulty arises from the combination of two unfavorable characteristics: First, the dense fracturing frustrates an efficient waterflood; second, because of the oil-wetness, the water pressure exceeds the oil pressure inside the matrix blocks, thus inhibiting spontaneous imbibition of water. In the past decade, using a new class of surfactants, enhanced oil recovery (EOR) researchers have studied the options to chemically revert the wettability of carbonate rock without drastically decreasing the oil-water interfacial tension. These chemicals, termed "wettability modifiers" (WMs), effectively reverse the sign of capillary pressure at the prevalent saturation. With the oil pressure exceeding the water pressure, the capillary pressure becomes the driving force for oil expulsion from the matrix and into the fracture system. Previous publications on chemical wettability modification focused on the performance of different chemical wettability modifiers for a chosen rock/oil/brine system. In some cases, they demonstrated an almost full oil recovery from core plugs. Little attention, however, has been given to the mechanism underlying the transport of the chemical into the matrix block and to the proper scaling of laboratory results to reservoir size. The present study aims to demonstrate that imbibition after wettability modification is diffusion-limited. To this end, the recovery profiles for spontaneous capillary imbibition, as well as for imbibition after wettability modification, are calculated. The results are then used to compare with the data of Amott cell imbibition experiments. It is confirmed that in both cases, the cumulative recovery is initially proportional to the square root of time. Imbibition after wettability modification, however, takes approximately 1,000 times longer than spontaneous capillary imbibition into a water-wet medium. The slow recovery observed in the case of imbibition after wettability modification is in excellent agreement with the assumption that, in the absence of significant spontaneous imbibition, the WM, to unfold its action, must first diffuse into the porous medium. In any diffusion process, the time scale is linked to the square of the length scale of the medium. Therefore, it would take up to 1,000 times longer (an equivalent of 200 years) before the same recovery is obtained from a meter-scale matrix block as is obtained from a centimeter-scale plug in a laboratory in 100 days. Consequently, unless a significantly faster transport mechanism for the wettability modifier is identified, or unless viscous forces or buoyancy enable forced imbibition, the chemical wettability modification of fractured oil-wet carbonate rock does not provide an economically interesting opportunity. Introduction Rock fractures provide comparatively highly permeable flow paths through oil reservoirs. In a densely fractured reservoir, the permeability contrast between the fracture network and the oil-bearing matrix can be significant. In that case, the viscous pressure differential across individual matrix blocks can be too small to release oil from the blocks under waterflood, thus leading to a poor recovery. Depending on the wetting state of the matrix and its initial water saturation, Swi, capillary action can cause imbibition of water up to a "spontaneous" equilibrium saturation, commonly denoted as Sspw. At this saturation, however, the capillary pressure inside the matrix block coincides with that in the fracture, and the recovery ceases. Experience has shown that carbonate fields often range from intermediate-wet to preferentially oil-wet (Treiber et al. 1972; Chilingar and Yen 1983), which is synonymous with Sspw being close or equal to Swi ; thus, they exhibit very limited recovery during primary and secondary production. Recently, a new EOR technique, designed specifically to tackle the challenges outlined previously, has been suggested by Austad and coworkers (Austad and Milter 1997; Standnes and Austad 2000a, b). In their pioneering work, these authors show that certain chemicals, when dissolved in the surrounding brine, can initiate water imbibition into oil-saturated core plugs and, hence, lead to the recovery of oil. One possible mechanism that explains these observations is the solubilization of adsorbed hydrocarbon components from the pore surface—as demonstrated by an atomic force microscopy study by Kumar et al. (2005), this exposes the intrinsically hydrophilic matrix. Another possibility is the formation of an additional chemical layer covering the adsorbed hydrophobic material. In either case, the pore surface becomes more hydrophilic, and the wettability of the matrix is thus modified. In a capillary rise experiment into parallel plates, Kumar et al. also observed different time scales for different types of wettability-modifying chemicals (2005). Using the cationic wettability modifier dodecyl trimethyl ammonium bromide (DTAB, also known as C12TAB), Standnes and Austad deduced that wettability modification was achieved through the comparatively slow process of partitioning the chemical into the oil phase, followed by desorption and solubilization of anionic hydrocarbon components (2000a, b). Shen et al. (2006) and Rao et al. (2006) measured the effect of surfactants on the relative water/oil permeabilities at different interfacial tensions. Wu et al. (2006) studied the properties and ranked the efficiency of chemical model compounds, based on their chemical structure, to modify the wettability and enhance recoveries. Several groups have taken initiative to model wettability modification in numerical simulators (Adibhatla et al. 2005; Delshad et al. 2006). So far, no significant attention has been given to time dependence and to the subsequent upscaling of the laboratory results to matrix block scale. This subject will be addressed in the present work. The structure of the article is as follows: In the Theory section, the basic results for countercurrent capillary imbibition will be briefly reviewed and compared to Fick's law of molecular diffusion. The oil recovery as a function of time for both capillary imbibition and imbibition after wettability modification will be predicted. The experimental approach to imbibition at different wetting situations will be described in the section Materials and Preparation. The recovery results will then be analyzed using the previously derived equations. Finally, tentative conclusions for the upscaling will be drawn.


Author(s):  
Qichen Zhang ◽  
Xiaodong Kang ◽  
Huiqing Liu ◽  
Xiaohu Dong ◽  
Jian Wang

AbstractCurrently, the reservoir heterogeneity is a serious challenge for developing oil sands with SAGD method. Nexen’s Long Lake SAGD project reported that breccia interlayer was widely distributed in lower and middle part of reservoir, impeding the steam chamber expansion and heated oil drainage. In this paper, two physical experiments were conducted to study the impact of breccia interlayer on development of steam chamber and production performance. Then, a laboratory scale numerical simulation model was established and a history match was conducted based on the 3D experimental results. Finally, the sensitivity analysis of thickness and permeability of breccia layer was performed. The influence mechanism of breccia layer on SAGD performance was analyzed by comparing the temperature profile of steam chamber and production dynamics. The experimental results indicate that the existence of breccia interlayer causes a thinner steam chamber profile and longer time to reach the peak oil rate. And, the ultimate oil recovery reduced 15.8% due to much oil stuck in breccia interlayer areas. The numerical simulation results show that a lower permeability in breccia layer area has a serious adverse impact on oil recovery if the thickness of breccia layer is larger, whereas the effect of permeability on SAGD performance is limited when the breccia layer is thinner. Besides, a thicker breccia layer can increase the time required to reach the peak oil rate, but has a little impact on the ultimate oil recovery.


SPE Journal ◽  
2018 ◽  
Vol 23 (06) ◽  
pp. 2243-2259 ◽  
Author(s):  
Pengfei Dong ◽  
Maura Puerto ◽  
Guoqing Jian ◽  
Kun Ma ◽  
Khalid Mateen ◽  
...  

Summary Oil recovery in heterogeneous carbonate reservoirs is typically inefficient because of the presence of high-permeability fracture networks and unfavorable capillary forces within the oil-wet matrix. Foam, as a mobility-control agent, has been proposed to mitigate the effect of reservoir heterogeneity by diverting injected fluids from the high-permeability fractured zones into the low-permeability unswept rock matrix, hence improving the sweep efficiency. This paper describes the use of a low-interfacial-tension (low-IFT) foaming formulation to improve oil recovery in highly heterogeneous/fractured oil-wet carbonate reservoirs. This formulation provides both mobility control and oil/water IFT reduction to overcome the unfavorable capillary forces preventing invading fluids from entering an oil-filled matrix. Thus, as expected, the combination of mobility control and low-IFT significantly improves oil recovery compared with either foam or surfactant flooding. A three-component surfactant formulation was tailored using phase-behavior tests with seawater and crude oil from a targeted reservoir. The optimized formulation simultaneously can generate IFT of 10−2 mN/m and strong foam in porous media when oil is present. Foam flooding was investigated in a representative fractured core system, in which a well-defined fracture was created by splitting the core lengthwise and precisely controlling the fracture aperture by applying a specific confining pressure. The foam-flooding experiments reveal that, in an oil-wet fractured Edward Brown dolomite, our low-IFT foaming formulation recovers approximately 72% original oil in place (OOIP), whereas waterflooding recovers only less than 2% OOIP; moreover, the residual oil saturation in the matrix was lowered by more than 20% compared with a foaming formulation lacking a low-IFT property. Coreflood results also indicate that the low-IFT foam diverts primarily the aqueous surfactant solution into the matrix because of (1) mobility reduction caused by foam in the fracture, (2) significantly lower capillary entry pressure for surfactant solution compared with gas, and (3) increasing the water relative permeability in the matrix by decreasing the residual oil. The selective diversion effect of this low-IFT foaming system effectively recovers the trapped oil, which cannot be recovered with single surfactant or high-IFT foaming formulations applied to highly heterogeneous or fractured reservoirs.


Sign in / Sign up

Export Citation Format

Share Document