Use of a Reservoir Simulator To Interpret Laboratory Waterflood Data

1973 ◽  
Vol 13 (06) ◽  
pp. 343-347 ◽  
Author(s):  
John S. Archer ◽  
S.W. Wong

Abstract Relative permeability curves calculated from laboratory waterflood history by the method of Johnson, Bossler and Naumann (JBN) are often poorly defined or anomalous at low and intermediate poorly defined or anomalous at low and intermediate water saturations. Poor definition can be encountered with strongly water-wet homogeneous cores when the displacement is piston-like. Anomalous curve shapes are associated with laboratory-observed water breakthrough ahead of the main flood front and are common in cores that have contrasting permeability streaks. The JBN technique, although permeability streaks. The JBN technique, although valid for the conditions assumed in its development, is unsatisfactory for the conditions specified above. A reservoir simulator has been used to model laboratory tests and thereby provide an alternative interpretation procedure. The simulation uses core properties and trial-and-error relative permeabilities. properties and trial-and-error relative permeabilities. The shapes of the relative permeability curves are adjusted until calculated oil recovery and relative injectivity curves match those obtained from the laboratory displacement tests. The technique has been used successfully to obtain meaningful relative permeability curves for piston-like displacement, mixed wettability systems, piston-like displacement, mixed wettability systems, and heterogeneous carbonates. The technique has also been used in evaluating empirical equations for calculating relative permeability. Introduction Numerical reservoir simulators are finding increasing application in production history matching and performance predictions. Because of the degree of sophistication reached with these models, it is mandatory that the fluid flow properties be of the highest possible quality. Of all the rock and fluid properties required in predicting performance, it is properties required in predicting performance, it is often the relative permeability characteristics that are the most critically important. These data are usually obtained from laboratory waterflood tests using reservoir core samples. The laboratory waterflood test is an attempt to represent the linear displacement behavior of the oil/water/reservoir-rock system. The wettability properties of the rock system should be preserved properties of the rock system should be preserved in the laboratory core sample if reliable results are to be obtained. Furthermore, the viscosity ratio and surface tension of the oil/water system in the laboratory test should ideally be made the same as those in the reservoir. In interpreting laboratory waterflood tests the unsteady-state equations are usually solved by methods of Buckley-Leverett, Welge and Johnson, Bossler and Neumann (JBN). These interpretations are sometimes inadequate for defining relative permeability curves for heterogeneous reservoir permeability curves for heterogeneous reservoir rock systems or for water displacing a very light oil in a homogeneous sandstone. For example, a number of writers have observed anomalous changes in the relative permeability to water during the flooding of heterogeneous carbonate core samples. The relative permeability to water does not increase smoothly with increasing water saturation, but increases stepwise or even humps. Such behavior appears to reflect small-scale local heterogeneity in the core sample and is likely to be insignificant on a field scale. The heterogeneity is often indicated in the laboratory by an observed water breakthrough at the core-sample production face ahead of the main flood front. The time of water breakthrough is an important measurement used in the calculation of relative permeability by the JBN method. If the breakthrough permeability by the JBN method. If the breakthrough time observed is not that of the main flood front but is a little early, then the relative permeabilities calculated will not represent the properties of the bulk of the core sample. It is under these conditions that anomalous relative permeability curves usually occur. We suggest in this paper that, in many cases, the small changes in pressure and in oil and water production rate that accompany anomalous relative production rate that accompany anomalous relative permeability curves can be smoothed to reflect permeability curves can be smoothed to reflect properties more consistent with the bulk behavior properties more consistent with the bulk behavior of the core sample. In essence, we are saying that together the smoothed oil and water production history and the pressure history of the laboratory core sample represent a unique property of that sample. SPEJ P. 343

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.


1972 ◽  
Vol 12 (05) ◽  
pp. 398-402 ◽  
Author(s):  
Necmettin Mungan

Abstract Water-oil relative permeabilities measured using reservoir fluids and fresh, preserved cores are shown to differ considerably from those obtained routinely using relined fluids and extracted cores. Upon saturating the extracted cores with the reservoir fluids and allowing them to come to equilibrium, in this case for 6 days, the original relative permeability curves were reestablished. Introduction Ordinarily, relative permeabilities are measured using refined fluids and restored-state cores. On occasion, fresh core samples that had been flushed with a refined oil, without any other cleaning or extraction, have been used. To date, no relative permeability measurements have been reported permeability measurements have been reported using reservoir fluids and fresh preserved cores, although a number of published studies have dealt with the influence of core handling and other laboratory experimental factors on results of displacement tests. This is not surprising. Special procedures have to be followed for obtaining, procedures have to be followed for obtaining, preserving, and handling native-state cores; and preserving, and handling native-state cores; and companies are generally reluctant to permit any procedures that may increase the expense or the procedures that may increase the expense or the time required to complete a well. Consequently, suitable cores often are simply not available. Furthermore, only a very few of the petroleum production research laboratories have the facilities production research laboratories have the facilities for measuring interfacial tension, contact angle and relative permeabilities at the elevated temperatures and pressures normally encountered in reservoirs. The following study provides, for the first time, relative permeability data obtained with fresh, preserved cores and reservoir fluids at reservoir preserved cores and reservoir fluids at reservoir pressure and temperature. Measurements were also pressure and temperature. Measurements were also made routinely with refined fluids and extracted cores to afford comparisons. Advancing and receding contact angles were measured as a function of time with the actual reservoir fluids on a solid surface representative of the reservoir rock to characterize the wetting equilibrium during the tests. Finally, a procedure was devised which, for the extracted procedure was devised which, for the extracted cores, yielded the original set of relative permeability curves. permeability curves. EXPERIMENTAL THE CORE The core was cut using lease crude oil in a Pennsylvanian sandstone reservoir. Crude oh was Pennsylvanian sandstone reservoir. Crude oh was considered as the best coring fluid to preserve both reservoir wettability and the interstitial water. The core was obtained from a pumping well that was to be deepened 50 ft for improving the productivity index. The production string was 4 1/2-in. casing and had been cemented and perforated with the total depth at 5,200 ft. Initially, a three-cone, hard-rock rotary bit was run on 2-3/8-in. drill pipe to drill out the cement plug and the casing flow shoe. Drilling was continued for an additional 10 ft to cut new formation and clean out the hole. For the coring operation, the rotary bit was replaced by a diamond core barrel having 3 1/2-in. OD. De-gassed crude oil containing no additives was used as the drilling fluid during both drilling and coring. This crude had an API gravity of 35.5 degrees (i.e., a specific gravity of 0.8473 gm/cc) and provided sufficient pressure against the formation. The intake pipe for pressure against the formation. The intake pipe for the mud pump was positioned several feet above the bottom of a production-storage tank. This was done to avoid picking up any sediments or asphalts that may have accumulated at the tank bottom. The mud return line from the well was diverted into a pit, and the return crude was not recirculated to avoid introduction of any oxygen into the well. The total quantity of the crude oil used was approximately 400 bbl during the entire operation. When brought to the surface, the core was quickly canned in crude oil to minimize exposure to air. In the laboratory, the core container was placed inside an airtight lucite box, which was fitted with sealed gloves for sample manipulation from the outside. A nitrogen atmosphere was maintained inside the box. An approximately 1-ft-long core piece was chosen which, by visual examination, appeared to be uniform. From each end, a 1-in. piece was cut for the x-ray diffraction and mercury injection studies. SPEJ P. 398


1965 ◽  
Vol 5 (04) ◽  
pp. 329-332 ◽  
Author(s):  
Larman J. Heath

Abstract Synthetic rock with predictable porosity and permeability bas been prepared from mixtures of sand, cement and water. Three series of mixes were investigated primarily for the relation between porosity and permeability for certain grain sizes and proportions. Synthetic rock prepared of 65 per cent large grains, 27 per cent small grains and 8 per cent Portland cement, gave measurable results ranging in porosity from 22.5 to 40 per cent and in permeability from 0.1 darcies to 6 darcies. This variation in porosity and permeability was caused by varying the amount of blending water. Drainage- cycle relative permeability characteristics of the synthetic rock were similar to those of natural reservoir rock. Introduction The fundamental behavior characteristics of fluids flowing through porous media have been described in the literature. Practical application of these flow characteristics to field conditions is too complicated except where assumptions are overly simplified. The use of dimensionally scaled models to simulate oil reservoirs has been described in the literature. These and other papers have presented the theoretical and experimental justification for model design. Others have presented elements of model construction and their operation. In most investigations the porous media have consisted of either unconsolidated sand, glass beads, broken glass or plastic-impregnated granular substances-materials in which the flow behavior is not identical to that in natural reservoir rock. The relative permeability curves for unconsolidated sands differ from those for consolidated sandstone. The effect of saturation history on relative permeability measurements A discussed by Geffen, et al. Wygal has shown quite conclusively that a process of artificial cementation can be used to render unconsolidated packs into synthetic sandstones having properties similar to those of natural rock. Many theoretical and experimental studies have been made in attempts to determine the structure and properties of unconsolidated sand, the most notable being by Naar and Wygal. Others have theorized and experimented with the fundamental characteristics of reservoir rocks. This study was conducted to determine if some general relationship could be established between the size of sand grains and the porosity and permeability in consolidated binary packs. This paper presents the results obtained by changing some of the factors which affect the porosity and permeability of synthetically prepared sandstone. In addition, drainage relative permeability curves are presented. EXPERIMENTAL PROCEDURE Mixtures of Portland cement with water and aggregate generally are designed to have certain characteristics, but essentially all are planned to be impervious to water or other liquids. Synthetic sandstone simulating oil reservoir rock, however, must be designed to have a given permeability (sometimes several darcies), a porosity which is primarily the effective porosity but quantitatively similar to natural rock, and other characteristics comparable to reservoir rock, such as wettability, pore geometry, tortuosity, etc. Unconsolidated ternary mixtures of spheres gave both a theoretically computed and an experimentally observed minimum porosity of about 25 per cent. By using a particle-distribution system, one-size particle packs had reproducible porosities in the reproducible range of 35 to 37 per cent. For model reservoir studies of the prototype system, a synthetic rock having a porosity of 25 per cent or less and a permeability of 2 darcies was required. The rock bad to be uniform and competent enough to handle. Synthetic sandstone cores mere prepared utilizing the technique developed by Wygal. Some tight variations in the procedure were incorporated. The sand was sieved through U.S. Standard sieves. SPEJ P. 329ˆ


Author(s):  
Pouriya Esmaeilzadeh ◽  
Mohammad Taghi Sadeghi ◽  
Alireza Bahramian

Many gas condensate reservoirs suffer a loss in productivity owing to accumulation of liquid in near-wellbore region. Wettability alteration of reservoir rock from liquid-wetting to gas-wetting appears to be a promising technique for elimination of the condensate blockage. In this paper, we report use of a superamphiphobic nanofluid containing TiO2 nanoparticles and low surface energy materials as polytetrafluoroethylene and trichloro(1H,1H,2H,2H-perfluorooctyl)silane to change the wettability of the carbonate reservoir rock to ultra gas-wetting. The utilization of nanofluid in the wettability alteration of carbonate rocks to gas-wetting in core scale has not been reported already and is still an ongoing issue. Contact angle measurements was conducted to investigate the wettability of carbonate core plugs in presence of nanofluid. It was found that the novel formulated nanofluid used in this work can remarkably change the wettability of the rock from both strongly water- and oil-wetting to highly gas-wetting condition. The adsorption of nanoparticles on the rock and formation of nano/submicron surface roughness was verified by Scanning Electron Microscope (SEM) and Stylus Profilometer (SP) analyses. Using free imbibition test, we showed that the nanofluid can imbibe interestingly into the core sample, resulting in notable ultimate gas-condensate liquid recovery. Moreover, we studied the effect of nanofluid on relative permeability and recovery performance of gas/water and gas/oil systems for a carbonate core. The result of coreflooding tests demonstrates that the relative permeability of both gas and liquid phase increased significantly as well as the liquid phase recovery enhanced greatly after the wettability alteration to gas-wetting.


1982 ◽  
Vol 22 (01) ◽  
pp. 108-116 ◽  
Author(s):  
John R. Counsil ◽  
Henry J. Ramey

Abstract Liquid vaporization can influence the results of unsteady, external gas-drive relative permeability experiments. At elevated temperatures, liquid vaporization may affectdisplacing gas mixture volume,displacing gas mixture viscosity, andvolumetric liquid saturation calculated from a material balance. Approximate methods are presented to correct laboratory displacement data for the effect of liquid vaporization on displacing gas mixture volume and viscosity. An approximate method also is presented to evaluate the magnitude of liquid saturation reduction caused by liquid vaporization. By use of a modified Jones and Roszelle calculation procedure, equations are developed to describe the dynamic displacement of liquid water by nitrogen gas at elevated temperatures. A conventional analysis of three displacement experiments demonstrated the apparent temperature dependence of gas relative permeability. Use of the proposed method indicated that corrected gas and water relative permeability curves are not strongly temperature dependent for the artificially consolidated sandstone cores used in this study. Introduction Relative permeability curves are required for numerical modelling of multiphase fluid flow through porous media. Although natural reservoir heterogeneity often reduces the utility of laboratory-derived relative permeabilities, laboratory studies are still required to understand basic fluid flow processes. Welge first modified the Buckley-Leverett theory and presented the equations for calculating (relative) permeability ratios from linear displacement data. Johnson et al. later extended this theory, to allow the calculation of individual relative permeabilities. The base permeability was the predrive, displacing-fluid effective permeability at the initial wetting, phase saturation. Jones and Roszelle then presented a simplified graphical technique that yielded individual relative permeabilities with the absolute (brine)permeability as a base. Osoba et al., Geffen et al., Welge, Rapoport and Leas, Stewart et al., Owens et al., Corey and Rathjens. Estes and Fulton. Richardson and Perkins, Craig et al., and others demonstrated the importance of end effects, flow rate, pressure gradient, drainage imbibition hysteresis, viscosity ratio, interfacial tension, contact angle, critical scaling factor, core heterogeneity, gas slippage, and other factors. In addition, temperature-dependent permeability effects were observed by Davidson, Poston et al., Weinbrandt et al., Casse and Ramey, and others. The reasons for the temperature effects were never fully understood or explained. This paper presents a method of eliminating some of the "apparent" gas relative-permeability temperature dependence by correcting approximately for the temperature- and pressure-dependent vapor/liquid phase behavior. Experimental Process and Apparatus The calculation procedure developed in this study models an isothermal, unsteady, linear gas drive. SPEJ P. 108^


Geophysics ◽  
2009 ◽  
Vol 74 (6) ◽  
pp. E251-E262 ◽  
Author(s):  
Marc H. Schneider ◽  
Patrick Tabeling ◽  
Fadhel Rezgui ◽  
Martin G. Lüling ◽  
Aurelien Daynes

Core analysis from reservoir rock plays an important role in oil and gas exploration as it can provide a large number of rock properties. Some of these rock properties can be extracted by image analysis of microscopic rock images in the visible light range. Such properties include the size, shape, and distribution of pores and grains, or more generally the texture, mineral distribution, and so on. A novel laboratory instrument and method allows for easy and reliable core imaging. This method is applicable even when the core sample is in poor shape. The capabilities of this technique can be verified by core images, image interpretation, and dynamic measurements of rock samples during flooding. A microscopic imager instrument is operated in video acquisition mode and can measure additional properties, such as fluid mobility, by detecting the emergence of injected fluids across the core sample.


SPE Journal ◽  
2009 ◽  
Vol 15 (02) ◽  
pp. 509-525 ◽  
Author(s):  
Yudou Wang ◽  
Gaoming Li ◽  
Albert C. Reynolds

Summary With the ensemble Kalman filter (EnKF) or smoother (EnKS), it is easy to adjust a wide variety of model parameters by assimilation of dynamic data. We focus first on the case where realizations and estimates of the depths of the initial fluid contacts, as well as grid- block rock-property fields, are generated by matching production data with the EnKS. Then we add the parameters defining power law relative permeability curves to the set of parameters estimated by assimilating production data with EnKS. The efficiency of EnKF and EnKS arises because data are assimilated sequentially in time and so "history matching data" requires only one forward run of the reservoir simulator for each ensemble member. For EnKS and EnKF to yield reliable characterizations of the uncertainty in model parameters and future performance predictions, the updated reservoir-simulation variables (e.g., saturations and pressures) must be statistically consistent with the realizations of these variables that would be obtained by rerunning the simulator from time zero using the updated model parameters. This statistical consistency can be established only under assumptions of Gaussi- anity and linearity that do not normally hold. Here, we use iterative EnKS methods that are statistically consistent, and show that, for the problems considered here, iteration significantly improves the performance of EnKS.


2021 ◽  
Author(s):  
Vincenzo Tarantini ◽  
Cristian Albertini ◽  
Hana Tfaili ◽  
Andrea Pirondelli ◽  
Francesco Bigoni

Abstract Karst systems heterogeneity may become a nightmare for reservoir modelers in predicting presence, spatial distribution, impact on formation petrophysical characteristics, and particularly in dynamic behaviour prediction. Moreover, the very high resolution required to describe in detail the phenomena does not reconcile with the geo-cellular model resolution typically used for reservoir simulation. The scope of the work is to present an effective approach to predict karst presence and model it dynamically. Karst presence recognition started from the analysis of anomalous well behaviour and potential sources of precursors (logs, drilling evidence, etc.) to derive concepts for karst reservoir model. This first demanding step implies then characterizing each cell classified as karstified in terms of petrophysical parameters. In a two-phase flow, karst brings to fast travelling of water which leaves the matrix almost unswept. This feature was characterized through dedicated fine simulations, leading to an upscaling of relative permeability curves for a single porosity formulation. The workflow was applied to a carbonate giant field with a long production history under waterflood development. Firstly, a machine learning algorithm was trained to recognize karst features based on log response, seismic attributes, and well dynamic evidence, then a karst probability volume was generated and utilized to predict the karst presence in the field. Karst characterization just in terms of porosity and permeability is sufficient to model the reservoir when still in single phase, however it fails to reproduce observed water production. Karst provides a high permeability path for water transport: classical history match approaches, such as the introduction of permeability multipliers, proved to be ineffective in reproducing the water breakthrough timing and growth rate. In fact, the reservoir consists of two systems, matrix, and karst: however, the karst is less known and laboratory analysis shows relative permeability only for the matrix medium. The introduction of equivalent or pseudo-relative permeability curves, accounting for both the media, was crucial for correct modelling of the reservoir underlying dynamics, allowing a proper reproduction of water breakthrough timing and water cut (WCT) trends. The implementation of a dedicated pseudo relative permeability curve dedicated to karstified cells allowed to replicate early water arrival, thus bringing to a correct prediction of oil and water rates, also highlighting the presence of bypassed oil associated with water circuiting, particularly in presence of highly karstified cells.


2021 ◽  
Author(s):  
Subodh Gupta

Abstract The objective of this paper is to present a fundamentals-based, consistent with observation, three-phase flow model that avoids the pitfalls of conventional models such as Stone-II or Baker's three-phase permeability models. While investigating the myth of residual oil saturation in SAGD with comparing model generated results against field data, Gupta et al. (2020) highlighted the difficulty in matching observed residual oil saturation in steamed reservoir with Stone-II and Baker's linear models. Though the use of Stone-II model is very popular for three-phase flow across the industry, one issue in the context of gravity drainage is how it appears to counter-intuitively limit the flow of oil when water is present near its irreducible saturation. The current work begins with describing the problem with existing combinatorial methods such as Stone-II, which in turn combine the water-oil, and gas-oil relative permeability curves to yield the oil relative permeability curve in presence of water and gas. Then starting with the fundamentals of laminar flow in capillaries and with successive analogical formulations, it develops expressions that directly yield the relative permeabilities for all three phases. In this it assumes a pore size distribution approximated by functions used earlier in the literature for deriving two-phase relative permeability curves. The outlined approach by-passes the need for having combinatorial functions such as prescribed by Stone or Baker. The model so developed is simple to use, and it avoids the unnatural phenomenon or discrepancy due to a mathematical artefact described in the context of Stone-II above. The model also explains why in the past some researchers have found relative permeability to be a function of temperature. The new model is also amenable to be determined experimentally, instead of being based on an assumed pore-size distribution. In that context it serves as a set of skeletal functions of known dependencies on various saturations, leaving constants to be determined experimentally. The novelty of the work is in development of a three-phase relative permeability model that is based on fundamentals of flow in fine channels and which explains the observed results in the context of flow in porous media better. The significance of the work includes, aside from predicting results more in line with expectations and an explanation of temperature dependent relative permeabilities of oil, a more reliable time dependent residual oleic-phase saturation in the context of gravity-based oil recovery methods.


Sign in / Sign up

Export Citation Format

Share Document