Assessment and Inclusion of Capillary Pressure/Relative Permeability Hysteresis Effects in Downhole Water Sink (DWS) Well Technology for Water Coning Control

Author(s):  
Solomon O. Inikori ◽  
Andrew K. Wojtanowicz

Abstract The objective of this study is to assess the effects of capillary pressures and relative permeability hysteresis on the performance of wells using the downhole water sink (DWS) technology for water coning control. In the study a commercial reservoir simulator has been adopted to evaluate well performance under conditions of stabilized oil production/water drainage rates for various combinations of these rates. Operational domain of water-free oil production, Inflow Performance Window (IPW), was used to quantify the effects of capillary pressure transition zone and relative permeability hysteresis on the water coning - control performance of DWS wells. Field data from wells in Canada, West Africa and Louisiana exhibiting severe problems of water coning were used in this study. The simulation results show that the basic concept of the DWS is unchanged by the inclusion of capillary pressure and relative permeability hysteresis. However, these effects may cause considerable reduction in the size of the water-free oil production domain and lead to increase in water production. The results also indicate that, for the same reservoir, converting conventional wells with prior water coning history to DWS application would not be as beneficial as DWS completions on new wells. Thus the effect of drainage-imbibition relative permeability hysteresis should be included in the DWS well design practice.

2021 ◽  
Vol 5 (1) ◽  
pp. 119-131
Author(s):  
Frzan F. Ali ◽  
Maha R. Hamoudi ◽  
Akram H. Abdul Wahab

Water coning is the biggest production problem mechanism in Middle East oil fields, especially in the Kurdistan Region of Iraq. When water production starts to increase, the costs of operations increase. Water production from the coning phenomena results in a reduction in recovery factor from the reservoir. Understanding the key factors impacting this problem can lead to the implementation of efficient methods to prevent and mitigate water coning. The rate of success of any method relies mainly on the ability to identify the mechanism causing the water coning. This is because several reservoir parameters can affect water coning in both homogenous and heterogeneous reservoirs. The objective of this research is to identify the parameters contributing to water coning in both homogenous and heterogeneous reservoirs. A simulation model was created to demonstrate water coning in a single- vertical well in a radial cross-section model in a commercial reservoir simulator. The sensitivity analysis was conducted on a variety of properties separately for both homogenous and heterogeneous reservoirs. The results were categorized by time to water breakthrough, oil production rate and water oil ratio. The results of the simulation work led to a number of conclusions. Firstly, production rate, perforation interval thickness and perforation depth are the most effective parameters on water coning. Secondly, time of water breakthrough is not an adequate indicator on the economic performance of the well, as the water cut is also important. Thirdly, natural fractures have significant contribution on water coning, which leads to less oil production at the end of production time when compared to a conventional reservoir with similar properties.


2002 ◽  
Vol 124 (4) ◽  
pp. 253-261 ◽  
Author(s):  
Andrew K. Wojtanowicz ◽  
Ephim I. Shirman

Dual-completed wells with Downhole Water Sink (DWS) are used for water coning control in oil reservoirs with bottom water drive. In DWS wells, the second (bottom) completion—placed in the water column—is used for draining water. This prevents the water cone invasion and allows free oil inflow in the top completion. The decision on using DWS or a conventional (single-completed) well is based upon deliverability comparison of the two wells. This paper shows how to describe DWS well deliverability in terms of the top and bottom production rates, water cut, and pressure drawdown. Also, the effect of pressure interference between two well completions on deliverability limits has been studied and qualified experimentally. DWS well deliverability depends on two variables, pressure drawdown and water drainage rate, and is described by a three-dimensional Inflow Performance Domain (IPD). Visual-Basic software based on a new analytical model of IPD has been developed to calculate critical (fluid breakthrough) rates for oil and water. The critical rates identify inflow conditions to the well’s completions—single or two-phase inflow. Also calculated are the values of water cut and maximum pressure drawdown at the well. An example demonstrates the procedure and a complete IPD plot. The experimental study, using a Hele-Shaw physical model of DWS well, demonstrates the reduction of well’s deliverability caused by pressure interference from the second (bottom) completion. The experiments have shown, however, that the deliverability decrease is small and over-compensated by the increase of oil rate due to simultaneous reduction of water cut.


Author(s):  
Marcos Faerstein ◽  
Paulo Couto ◽  
José Alves

This paper discusses the impacts that rock wettability may have upon the production and recovery of oil with waterflooding in carbonate reservoirs and how it should be modeled. A broad review of the state of the art has been conducted surveying existing disagreements and knowledge gaps, basic definitions, as well as the correct understanding of the physical phenomena and identification of the characteristics of the various wettability scenarios. Case studies conducted with a black oil reservoir simulator evaluated the impact of different wettability scenarios on oil production and recovery. A comprehensive approach considering all the parameters involved in the wettability modeling was applied to the case studies, showing how the behavior of the reservoir varies as a function of their wettability. This paper shows how relative permeability and capillary pressure should be varied to correctly represent different wettability scenarios and consequently assess its impacts on oil production and recovery. The case studies show that the evaluation of the volume of oil in the reservoir is impacted by wettability through the irreducible water saturation and primary drainage capillary pressure and must be considered in the analyses. In long term analyses, mixed-wet scenarios have a higher oil production and recovery. In medium and short term, the water-wet scenarios have the higher recovery, but in relation to oil production, these scenarios are negatively influenced by the smaller volume of oil in place. The main contribution of this paper is the simultaneous analyses of all the parameters involved in the modeling of wettability showing how they impact the behavior of a reservoir. It shows how the parameters must be varied in a heterogeneous reservoir and how heterogeneity impacts the relevance of wettability in the studies.


SPE Journal ◽  
2010 ◽  
Vol 15 (04) ◽  
pp. 1003-1019 ◽  
Author(s):  
Odd Steve Hustad ◽  
David John Browning

Summary A coupled formulation for three-phase capillary pressure and relative permeability for implicit compositional reservoir simulation is presented. The formulation incorporates primary, secondary, and tertiary saturation functions. Hysteresis and miscibility are applied simultaneously to both capillary pressure and relative permeability. Two alternative three-phase capillary pressure formulations are presented: the first as described by Hustad (2002) and the second that incorporates six representative two-phase capillary pressures in a saturation-weighting scheme. Consistency is ensured for all three two-phase boundary conditions through the application of two-phase data and normalized saturations. Simulation examples of water-alternating-gas (WAG) injection are presented for water-wet and mixed-wet saturation functions. 1D homogeneous and 2D and 3D heterogeneous examples are employed to demonstrate some model features and performance.


1977 ◽  
Vol 17 (04) ◽  
pp. 251-262 ◽  
Author(s):  
E.G. Woods ◽  
A.K. Khurana

Abstract Three-dimensional numerical models of bottom-water-drive reservoirs show delayed water breakthrough into individual wells when compared with observed well performance and individual-well coning models. This reservoir-model behavior results from masking of the well coning effect by volume-averaging pressure and saturation profiles around a well over a grid block with a large volume. The reservoir-simulator prediction of well performance can be improved by mathematically performance can be improved by mathematically transforming the production performance of a detailed well-coning model into a set of time-independent pseudorelative-permeability and capillary-pressure curves that then can be used in the reservoir model. A procedure for obtaining the required pseudofunctions is described and the results of their application in simple models and in a large reservoir-simulator model are shown. Introduction The prohibitive cost of numerical reservoir simulation with fine-grid definition models of large reservoirs has led to development of techniques whereby vertical saturation distribution and/or localized flow conditions in the vicinity of individual wells can be approximately accounted for in relatively coarse-grid models at an acceptable incremental cost. In particular, vertical cross-section models under capillary and gravity equilibrium have been used to derive pseudorelative permeabilities and capillary pressures for use in two-dimensional, areal models to simulate the average vertical distribution of flow without having to pay the computing price of a full three-dimensional model. Coats et al. described the use of the vertical equilibrium concept for developing pseudorelative-permeability and capillary-pressure pseudorelative-permeability and capillary-pressure functions for simulating the vertical dimension in a two-dimensional, areal simulator model This method assumes gravity-capillary equilibrium in the vertical direction. Also, Coats et al. developed a dimensionless parameter for estimating when these conditions are valid. Martin formed a mathematical basis for pseudofunctions by reducing the equations for pseudofunctions by reducing the equations for three-phase, three-dimensional, compressible flow to two-dimensional relations by partial integration of the equations of flow. Hearn extended the pseudorelative-permeability concept by adapting it pseudorelative-permeability concept by adapting it to stratified reservoirs where viscous rather than gravity and capillary forces dominate the vertical sweep efficiency. Hawthorne studied the effects of capillary pressure on pseudorelative permeability derived from the Hearn stratified model. Jacks et al. further enlarged thepseudorelative-perrneability concept by developing dynamic pseudorelative permeabilities. (Dynamic pseudos, denoting pseudos permeabilities. (Dynamic pseudos, denoting pseudos determined under flowing rather than static conditions, allow one to account for the interaction between viscous and gravity forces resulting from rate variation in the vertical plane.) Kyte and Berry generalized the work of Jacks et al. by introducing the concept of pseudocapillary pressures and improving dynamic pseudofunction calculations to include varying flow potential gradients. Emanual and Cook expanded the concept of vertical cross-section, pseudorelative permeabilities to include the vertical performance of individual wells. Their procedure combines the effect of coning and well pseudorelative permeabilities for use in a two-dimensional, areal model. Chappelear and Hirasaki used a different approach to including of coning effects in a two-dimensional, areal simulator by developing a functional relationship among water cut, average oil-column thickness, and total rate based on an analytical incompressible, steady-state model. The most sophisticated approach to representing well-coning effects in a reservoir simulator has been taken by Mrosovsky and Ridings and Akbar et al. They incorporated detailed numerical well models into the reservoir-model grid. SPEJ P. 251


1963 ◽  
Vol 3 (02) ◽  
pp. 164-176 ◽  
Author(s):  
Russell L. Nielsen ◽  
M.R. Tek

The scaling laws as formulated by Rapport relate dynamically similar flow systems in porous media each involving two immiscible, incompressible fluids. A two-dimensional numerical technique for solving the differential equations describing systems of this type has been employed to assess the practical value of the scaling laws in light of the virtually unscalable nature of relative permeability and capillary pressure curves and boundary conditions.Two hypothetical systems - a gas reservoir subject to water drive and the laboratory scaled model of that reservoir - were investigated with emphasis placed on water coning near a production well. Comparison of the computed behavior of these particular systems shows that water coning in the reservoir would be more severe than one would expect from an experimental study of a laboratory model scaled within practical limits to the reservoir system.This paper also presents modifications of the scaling laws which are available for systems that can be described adequately in two-dimensional Cartesian coordinates. Introduction Present day digital computing equipment and methods of numerical analysis allow realistic and quantitative studies to be carried out for many two-phase flow systems in porous media. Before these tools became available the anticipated behavior of systems of this type could be inferred only from analytical solutions of simplified mathematical models or from experimental studies performed on laboratory models.To reproduce the behavior of a reservoir system on the laboratory scale, certain relationships must be satisfied between physical and geometric properties of the reservoir and laboratory systems. Where the reservoir fluids may be considered as two immiscible and incompressible phases, the necessary relationships have been formulated by Rapoport and others. Rapoport's scaling laws follow from inspectional analysis of the differential equation describing phase saturation distribution in such systems.It will be recalled that these scaling laws presuppose three conditions:the relative permeability curves must be identical for the model and prototype;the capillary pressure curve (function of phase saturation) for the model must be linearly related to that of the prototype; andboundary conditions imposed on the model must duplicate those existing at the boundaries of the prototype. These three requirements seldom if ever can be satisfied in scaling an actual reservoir to the laboratory system because:The laboratory medium normally will be unconsolidated (glass beads or sand) while the reservoir usually is consolidated. Relative permeability and capillary pressure curves are usually quite different for consolidated and unconsolidated porous media.The reservoir usually will be surrounded by a large aquifer which could be simulated in the laboratory only to a limited extent.Wells present in the reservoir would scale to microscopic dimensions in the laboratory if geometric similarity is to be maintained. In view of these considerations, rigorous scaling of even a totally defined reservoir probably would never be possible.The purpose of this paper is to assess the practical value of the scaling laws in the light of the unscalable variables. This has been done by carrying out numerical solutions in two dimensions to the differential equations describing the flow of two immiscible, incompressible fluids in porous media for a field scale reservoir and a laboratory model of that reservoir. While both the reservoir and the laboratory model were purely fictional, each has been made as realistic and representative as possible.The field problem selected as the basis for the investigation was an inhomogeneous, layered gas reservoir initially at capillary gravitational equilibrium and subsequently produced in the presence of water drive. The laboratory model of this reservoir was designed to utilize oil and water in a glass bead pack. SPEJ P. 164^


SPE Journal ◽  
2014 ◽  
Vol 20 (02) ◽  
pp. 267-276 ◽  
Author(s):  
Xianhui Kong ◽  
Mojdeh Delshad ◽  
Mary F. Wheeler

Summary Numerical modeling and simulation are essential tools for developing a better understanding of the geologic characteristics of aquifers and providing technical support for future carbon dioxide (CO2) storage projects. Modeling CO2 sequestration in underground aquifers requires the implementation of models of multiphase flow and CO2 and brine phase behavior. Capillary pressure and relative permeability need to be consistent with permeability/porosity variations of the rock. It is, therefore, crucial to gain confidence in the numerical models by validating the models and results by use of laboratory and field pilot results. A published CO2/brine laboratory coreflood was selected for our simulation study. The experimental results include subcore porosity and CO2-saturation distributions by means of a computed tomography (CT) scanner along with a CO2-saturation histogram. Data used in this paper are all based on those provided by Krause et al. (2011), with the exception of the CT porosity data. We generated a heterogeneous distribution for the porosity but honoring the mean value provided by Krause et al. (2011). We also generated the permeability distribution with the mean value for the whole core given by Krause et al. (2011). All the other data, such as the core dimensions, injection rate, outlet pressure, temperature, relative permeability, and capillary pressure, are the same as those in Krause et al. (2011). High-resolution coreflood simulations of brine displacement with supercritical CO2 are presented with the compositional reservoir simulator IPARS (Wheeler and Wheeler 1990). A 3D synthetic core model was constructed with permeability and porosity distributions generated by use of the geostatistical software FFTSIM [Jennings et al. (2000)], with cell sizes of 1.27×1.27×6.35 mm. The core was initially saturated with brine. Fluid properties were calibrated with the equation-of-state (EOS) compositional model to match the measured data provided by Krause et al. (2011). We used their measured capillary pressure and relative permeability curves. However, we scaled capillary pressure on the basis of the Leverett J-function (Leverett 1941) for permeability, porosity, and interfacial tension (IFT) in every simulation grid cell. Saturation images provide insight into the role of heterogeneity of CO2 distribution in which a slight variation in porosity gives rise to large variations in CO2-saturation distribution in the core. High-resolution numerical results indicated that accurate representation of capillary pressure at small scales was critical. Residual brine saturation and the subsequent shift in the relative permeability curves showed a significant impact on final CO2 distribution in the core.


Author(s):  
Huiying Li ◽  
Sergio A. Vasquez ◽  
Mohammed Azhar

The present work concerns the development of a comprehensive model capability in ANSYS CFD software FLUENT towards modeling of multiphase flows in porous media with targeted applications in reservoir/well analyses. The modeling approach is based on the Eulerian multifluid model. Porous media are modeled by both Superficial and Physical Velocity formulation with embedded sub-models to account for the resistance sink, relative-permeability and capillary pressure effect. An advanced numerical algorithm has been developed to achieve time-step and mesh independent solutions as well as to satisfy the physical constraints/limits. In particular, the resistance sinks are rearranged and linearized to ensure numerical stability and to handle mathematically infinite resistance caused by possible zero relative permeability. The capillary pressure and body forces are implicitly treated to enhance solver robustness. The multiphase porous medium model is compatible with all the numerical schemes and solvers (iterative and non-iterative) available in FLUENT. The present model has been applied to simulate 1D, 2D and 3D transient oil-water two-phase flows mimicking the conditions in reservoirs and wells. The solutions are time-step and grid independent, and successfully reproduce the flow characteristics and physical limits. The solvers are fast and robust, allowing the time step to be as large as 2 hours for a reservoir setting with the flow physical time in 2–20 years. The model capability shows great promises for reservoir and well performance analysis.


SPE Journal ◽  
2012 ◽  
Vol 17 (04) ◽  
pp. 1047-1055 ◽  
Author(s):  
Romain L. Chassagne ◽  
Paul S. Hammond

Summary We used a commercial reservoir simulator to study, first, the dissipation of aqueous drilling fluid filtrate invasion around a cased observation well in an oil-saturated formation under the action of capillary pressure and, second, the interaction of a waterflood front with the cased well and remaining invaded zone. Hysteretic behavior of the capillary pressure and relative permeabilities is critically important to these processes and is taken into account by the use of the Carlson model, with the various bounding drainage and imbibition curves computed from a pore network model. Filtrate invasion into a hydrocarbon formation influences the readings of well-logging tools. Although this phenomenon has been known, and corrected for, for many years, uncertainty remains with regard to the long-time behavior of invasion around observation wells where no flow in or out of the formation occurs after completion, and with regard to the influence of formation wettability. We find that after sufficient time, the invaded zone dissipates completely in a water-wet formation, but some invasion always remains in the oil/mixed-wet case. Nonwetting-phase trapping, manifested through relative permeability hysteresis, is the cause. Because trapping affects the values and the endpoints of the relative permeability curves, a waterflood front passing across an observation well is more distorted in the oil/mixed-wet case. The simulation results allow us to understand how logging-tool measurements made in cased observation wells are influenced by drilling-fluid invasion and will therefore lead to improved interpretation. This study shows strong links between the wettability of the formation and the persistence of invaded zone saturation and between invaded zone saturation and the distortion of subsequent flood fronts.


Sign in / Sign up

Export Citation Format

Share Document