Design of Laboratory Models for Study of Miscible Displacement

1963 ◽  
Vol 3 (01) ◽  
pp. 28-40 ◽  
Author(s):  
Anthony L. Pozzi ◽  
Robert J. Blackwell

Abstract Scaled laboratory-model studies provide a powerful method for evaluation of a proposed oil-recovery process. In recent years, models have been used extensively to evaluate processes in which solvents displace oil, both for general cases and for specific reservoir conditions. Since the performance of a miscible flood in a horizontal reservoir can be significantly affected by transverse mixing between solvent and oil, this displacement mechanism must be accurately simulated in the scaled model studies. Unfortunately, precise scaling of transverse dispersion coupled with the requirement of geometric similarity requires impractically large laboratory models and long times for experiments.If scaling requirements for miscible displacements could be relaxed while accurate simulation of essential displacement mechanisms is maintained, the utility of model studies would be greatly enhanced. The purpose of the work reported herein was to evaluate the relative importance of various mechanisms affecting miscible displacement and to ascertain whether the essential features of the displacement process can be simulated even though some scaling groups are not satisfied. These studies were performed with completely miscible systems in linear, horizontal models packed with unconsolidated media.From the experimental results, a set of relaxed scaling criteria was formulated which allows the requirements of geometric similarity and equality of the ratio of viscous to gravity forces to be omitted for specified conditions. The relaxed criteria are valid whether transverse mixing is by molecular diffusion or by convective dispersion.Correlations which permit prediction of vertical sweep efficiencies in linear, horizontal reservoirs were developed from the experimental data when transverse mixing is by molecular diffusion, These same correlations may be used when transverse mixing is by convective dispersion if an empirically defined, effective, transverse dispersion coefficient is used in the description of the mixing process. The effective transverse dispersion coefficient correlation essentially duplicates the dispersion coefficient correlation for equal-viscosity, equal-density fluid systems. Experimental values for the effective transverse dispersion coefficient can be measured readily. Introduction One of the most effective methods for evaluation of miscible-displacement oil-recovery processes is that of displacements in laboratory models scaled to simulate reservoir conditions. For these laboratory studies to be meaningful, however, the essential displacement mechanisms affecting reservoir performance must be accurately simulated.Since the performance of a miscible flood in horizontal reservoirs, or in dipping reservoirs at high rates, can be significantly affected by transverse mixing of solvent and oil, this mechanism must be considered in the design of laboratory experiments. Unfortunately, precise scaling of transverse dispersion coupled with the requirement of geometric similarity requires impracticality large laboratory models and long experiment times. This difficulty seriously limits the utility of laboratory model studies.Craig, et al, demonstrated that geometric similarity is not required when mixing is unimportant. Their experimental data indicate, for the cases studied, that the displacement is sufficiently characterized by scaling the ratio of viscous-to-gravitational forces. This work suggests that relaxation of the requirement of geometric similarity and, possibly, other criteria might also be permissible when mixing is important, provided suitable groups describing the mixing process are scaled.The purpose of the work reported here was to evaluate the relative importance of various mechanisms affecting miscible displacement and to ascertain whether the essential features of the displacement process can be simulated even though some scaling groups are not satisfied. SPEJ P. 28^

1966 ◽  
Vol 6 (01) ◽  
pp. 73-80 ◽  
Author(s):  
J.L. Mahaffey ◽  
W.M. Rutherford ◽  
C.S. Matthews

Abstract This paper gives results of an experimental study of the sweep efficiency of a miscible displacement in a five-spot. The study was carried out in a parallel-plate glass model in which effects of diffusion were scaled at or near the molecular diffusion level. The experiments show that very early breakthrough (25 to 35 per cent of pore volume (PV) injected) may be expected in miscible floods because of the unfavorable viscosity ratio. However after 1 PV of displacing fluid is injected, the sweep rises to a reasonable value (50 to 60 per cent). Photographs show that small slugs of less than 10 per cent of PV tend to dissipate before breakthrough. A minimum slug size of 15 per cent of PV would appear to be necessary even in a relatively homogeneous formation. Presence of a slug whose viscosity is intermediate between that of oil and gas increases the sweep efficiency of the oil-gas system. In a typical system the sweep at breakthrough rose from 26 to 37 per cent of PV for a 25 per cent slug. The increase in sweep brought about by use of a large slug could well pay for the extra deferment cost of the additional slug material. Introduction Most miscible displacement processes involve the displacement of oil with fluids of much lower viscosity and density. The displacement process at these adverse viscosity and density ratios is dominated by instability phenomena, i.e., viscous fingers and gravity tongues. These phenomena have highly adverse effects on oil recovery. Although a number of laboratory studies have been made to determine the effect of adverse viscosity ratios on five-spot sweep patters,1,2 the scaling of diffusion effects is uncertain. In the series of scaled model studies reported herein, an attempt was made to scale diffusion. Model studies of miscible displacements in which molecular diffusion predominates are permitted by controlling the parallel plate spacing which reduces convective mixing to arbitrarily small levels. To decide how this scaling relates to any particular field displacement necessitates an estimation of diffusion effects for the natural rock being considered and conditions under which displacement will be conducted. The approach normally taken is to extrapolate data obtained from stable miscible displacements performed in the laboratory, such as those presented by Brigham et al.3 The validity with which such an extrapolation can be applied to an unstable flow system has yet to be established. If this approach is accepted, a family of oil recovery curves can be generated for a single viscosity ratio based on Brigham's observation that the magnitude of the dispersion coefficient is dependent, among other things, upon specific rock properties. Objective of our test was to define the lower limit of this range by presenting the case where dispersion effects were reduced to the molecular diffusion level in both the transverse and longitudinal directions. The scaling of diffusion effects can be handled in two-dimensional systems by the usse of narrow-gap, parallel-plate models. In parallel-plate models the Taylor diffusion coefficient for convective mixing in the direction of flow at low flow rates is given by (following Taylor4):Equation 1 where h is the plate spacing and D is the molecular diffusion coefficient. Clearly, convective mixing can be reduced to arbitrarily small levels by manipulating the gap spacing h. This was the method used in these studies.


1962 ◽  
Vol 2 (04) ◽  
pp. 317-326 ◽  
Author(s):  
C. Van Der Poel

Abstract When oil is displaced from a horizontal formation by another fluid of lower density, the latter tends to override the former in the shape of a tongue owing to gravity segregation. This gravity tongue has an adverse influence on the oil recovery. If the fluids are miscible, diffusion (mixing) takes place at the interfacial boundary of the gravity tongue. This mixing should have a favorable effect on oil recovery. The report describes a laboratory study of the magnitude of the mixing zone under various conditions, so as to assess the effect of diffusion on oil recovery both in laboratory experiments and under actual field conditions. The technique used enables visual observation and measurement of the size of the mixing zone in transparent glass-powder packs. The results show that in experiments in small models and cores the width of the mixing zone may well be of the same order of magnitude as the height of the model. In such cases oil recovery is favorably affected by mixing. It can further be concluded that, under conditions prevailing in the field, mixing of the injected fluid with reservoir oil is equal to that caused by molecular diffusion alone, eddy-mixing not taking place to any appreciable extent. A simple calculation, then, shows that molecular diffusion is too small for a beneficial effect to be expected from the injection of miscible fluids in horizontal or nearly horizontal reservoirs unless pay zones are thin. Introduction This paper gives results of experiments made on the mixing which occurs when a miscible displacement is carried out in a horizontal reservoir. This mixing takes place at the interfacial boundary of the gravity tongue formed when the lighter injected fluid overrides the oil present in the reservoir. The object of the experiments was to simulate the field case where, for example, propane with a viscosity of 0.075 cp under reservoir conditions displaces an oil of viscosity 0.6 cp. In the experiments the lighter fluid had a viscosity of about 1 cp and the heavier one a viscosity of about 8 cp, so as to obtain the same viscosity ratio. In order to enable the results to be compared with those published in the literature, a set of experiments with viscosity ratio equal to one was also performed. EXPERIMENTAL TECHNIQUE The technique developed for the purpose enables the width of the mixing zone to be studied as a function of time and place. A glass-powder pack is saturated with a water-glycerine mixture of suitable viscosity, which represents the reservoir fluid. The pack is then rendered transparent by dissolving sufficient ammonium thiocyanate in the mixture to obtain a solution with a refractive index matching that of the glass powder. Alkaline water to which phenolphthalein has been added is used as the lighter and less-viscous displacing fluid. In those places where mixing or diffusion of the two liquids occurs, the slightly acidic ammonium-thiocyanate solution neutralizes the alkali in the water, and the glass pack shows up white against the red-colored invading water. If a black cloth is hung over the back of the apparatus, the transparent part of the glass pack appears black. In this way the position of the two phases and the transition zone between them is clearly visible, as shown in Fig. 1 (where the red-colored invading water is the grey-shaded zone lying uppermost).In all experiments the amount of alkali in the water was chosen such that the upper contour corresponded to a concentration of 5 per cent of the dense liquid. The lower contour is determined by the size of the glass grains and the thickness of the pack and, consequently, varies from experiment to experiment. SPEJ P. 317^


2013 ◽  
Vol 16 (01) ◽  
pp. 85-96 ◽  
Author(s):  
Olaoluwa O. Adepoju ◽  
Larry W. Lake ◽  
Russell T. Johns

Summary Dispersion (or local mixing) degrades miscibility in miscible-flood displacements by interfering with the transfer of intermediate components that develop miscibility. Dispersion, however, also can improve oil recovery by increasing sweep efficiency. Either way, dispersion is an important factor in understanding miscible-flood performance. This paper investigates longitudinal and transverse local mixing in a finite-difference compositional simulator at different scales (both fine and coarse scale) using a 2D convection-dispersion model. All simulations were of constant-mobility and -density, first-contact miscible flow. The model allows for variations of velocity in both directions. We analyzed local (gridblock) concentration profiles for various miscible-displacement models with different scales of heterogeneity and permeability autocorrelation lengths. To infer dispersivity, we fitted an analytical 2D convection-dispersion model to the local concentration profile to determine local longitudinal and transverse dispersivities simultaneously. Streamlines of simulation models were traced using the algorithm proposed by Pollock (1988). To our knowledge, this is the first systematic attempt to numerically study local transverse dispersivity. The results show that transverse mixing, which is usually neglected in the 1D convection-dispersion model of dispersion, is significant when the flow direction changes locally as a result of heterogeneity. The computed streamlines, which highlight the variation in flow directions, agree with the computed transverse-dispersivity trends. We find that both transverse and longitudinal dispersion can grow with travel distance and that there are several instances in which transverse dispersion is the larger of the two. Often, the variations in the streamlines are suppressed (homogenized) during upscaling. This paper gives a quantitative and systematic procedure to estimate the degree of transverse mixing (dispersivity) in any model. We conclude that local mixing, including transverse mixing, should be considered when upscaling a fine-scale model for miscible displacement to ensure proper preservation of fine-scale sweep and displacement efficiency and ultimate oil recovery for miscible-displacement simulations.


1982 ◽  
Vol 22 (05) ◽  
pp. 647-657 ◽  
Author(s):  
J.P. Batycky ◽  
B.B. Maini ◽  
D.B. Fisher

Abstract Miscible gas displacement data obtained from full-diameter carbonate reservoir cores have been fitted to a modified miscible flow dispersion-capacitance model. Starting with earlier approaches, we have synthesized an algorithm that provides rapid and accurate determination of the three parameters included in the model: the dispersion coefficient, the flowing fraction of displaceable volume, and the rate constant for mass transfer between flowing and stagnant volumes. Quality of fit is verified with a finite-difference simulation. The dependencies of the three parameters have been evaluated as functions of the displacement velocity and of the water saturation within four carbonate cores composed of various amounts of matrix, vug, and fracture porosity. Numerical simulation of a composite core made by stacking three of the individual cores has been compared with the experimental data. For comparison, an analysis of Berea sandstone gas displacement also has been provided. Although the sandstone displays a minor dependence of gas recovery on water saturation, we found that the carbonate cores are strongly affected by water content. Such behavior would not be measurable if small carbonate samples that can reflect only matrix properties were used. This study therefore represents a significant assessment of the dispersion-capacitance model for carbonate cores and its ability to reflect changes in pore interconnectivity that accompany water saturation alteration. Introduction Miscible displacement processes are used widely in various aspects of oil recovery. A solvent slug injected into a reservoir can be used to displace miscibly either oil or gas. The necessary slug size is determined by the rate at which deterioration can occur as the slug is Another commonly used miscible process involves addition of a small slug within the injected fluids or gases to determine the nature and extent of inter well communication. The quantity of tracer material used is dictated by analytical detection capabilities and by an understanding of the miscible displacement properties of the reservoir. We can develop such understanding by performing one-dimensional (1D) step-change miscible displacement experiments within the laboratory with selected reservoir core material. The effluent profiles derived from the experiments then are fitted to a suitable mathematical model to express the behavior of each rock type through the use of a relatively small number of parameters. This paper illustrates the efficient application of the three-parameter, dispersion-capacitance model. Its application previously has been limited to use with small homogeneous plugs normally composed of intergranular and intencrystalline porosity, and its suitability for use with cores displaying macroscopic heterogeneity has been questioned. Consequently, in addition to illustrating its use with a homogeneous sandstone, we fit data derived from previously reported full-diameter carbonate cores. As noted earlier, these cores were heterogeneous, and each of them displayed different dual or multiple types of porosity characteristic of vugular and fractured carbonate rocks. Dispersion-Capacitance Model The displacement efficiency of one fluid by a second immiscible fluid within a porous medium depends on the complexity of rock and fluid properties. SPEJ P. 647^


Author(s):  
Teresa B. Hoberg ◽  
Kenshiro Muramatsu ◽  
Erica M. Cherry ◽  
John K. Eaton

Open-cell metal foams are of interest for a variety of thermal engineering applications because of their high surface-to-volume ratio and high convective heat transfer coefficients relative to conventional fins. The tortuous flow path through the foam promotes rapid transverse mixing, a fact that is important in heat exchanger applications. Transverse mixing acts to spread heat away from a heated surface, bringing cooler fluid to the foam elements that are in direct contact with the surface. Heat is also spread by conduction in the foam ligaments. The present work addresses fully-coupled thermal dispersion in a metal foam. Dispersion of the thermal wake of a line source was measured. A conjugate heat transfer model was developed which showed good agreement with the data. The validated model was used to examine the complementary effects of the mechanical dispersion, molecular diffusion in the gas, and conduction in the solid.


2019 ◽  
Vol 7 (9) ◽  
pp. 291 ◽  
Author(s):  
Xiang Cui ◽  
Changqi Zhu ◽  
Mingjian Hu ◽  
Xinzhi Wang ◽  
Haifeng Liu

Dispersion characteristics are important factors affecting groundwater solute transport in porous media. In marine environments, solute dispersion leads to the formation of freshwater aquifers under islands. In this study, a series of model tests were designed to explore the relationship between the dispersion characteristics of solute in calcareous sands and the particle size, degree of compactness, and gradation of porous media, with a discussion of the types of dispersion mechanisms in coral sands. It was found that the particle size of coral sands was an important parameter affecting the dispersion coefficient, with the dispersion coefficient increasing with particle size. Gradation was also an important factor affecting the dispersion coefficient of coral sands, with the dispersion coefficient increasing with increasing d10. The dispersion coefficient of coral sands decreased approximately linearly with increasing compactness. The rate of decrease was −0.7244 for single-grained coral sands of particle size 0.25–0.5 mm. When the solute concentrations and particle sizes increased, the limiting concentration gradients at equilibrium decreased. In this study, based on the relative weights of molecular diffusion versus mechanical dispersion under different flow velocity conditions, the dispersion mechanisms were classified into five types, and for each type, a corresponding flow velocity limit was derived.


2019 ◽  
Vol 2019 ◽  
pp. 1-25
Author(s):  
Shaorui Sun ◽  
Feng Zhu ◽  
Jihong Wei ◽  
Wuchao Wang ◽  
Huilin Le

Soil-rock mixture (SRM) is a special geological material that has the unique properties of rock or soil. Studies on the strength characteristics of the SRM have very important theoretical significance and practical value. In this study, the gravel proportions, gravel sizes, gravel shapes, and repetitive results of shear experiments are considered in laboratory experiments and for the identification of strength parameters. To the gravel shapes, from the angle of the composition materials of SRM, the experimental samples include samples with breccia gravels and with subrounded gravels. And, in this study, the laboratory model experiment is used to research the strength characteristics of SRM. In addition, the shear failure mechanism is used to establish the relationship between the microfailure mechanism and the macrostrength parameters identification. Taking gravel proportion, gravel size, gravel shape, and repetitive shear process as the influencing factors of the SRM, the laboratory models have been remolded, and laboratory direct shear experiments have been carried out. The shear deformation laws of the SRM are researched on the basis of the analysis of the curves of shear stress and the horizontal displacement.


1965 ◽  
Vol 5 (02) ◽  
pp. 160-166 ◽  
Author(s):  
A.M. Rowe ◽  
I.H. Silberberg

Abstract A computer program was written to predict the phase behavior generated by the enriched-gas-drive process. This program is based, in part, on a new concept of convergence pressure, which is then used to select vapor-liquid equilibrium ratios (K-factors) for performing a series of flash calculations. The results of these calculations are the equilibrium vapor and liquid phase compositions which define the phase envelopes. The program was used to predict phase envelopes for 11 different hydrocarbon systems on which published experimental data were available. The predicted and experimental results compare favorably. Introduction The reservoir engineer is frequently faced with the problem of predicting what will happen if gas is injected into a reservoir. One aspect of this general problem is predicting the phase changes that will occur when a non-equilibrium gas displaces a reservoir fluid. In particular, a "dry" gas, upon displacing a volatile oil will pick up intermediate components from the oil. On the other hand, a "wet" gas, containing a high concentration of intermediates, will lose some of these components to a relatively low-gravity, non-equilibrium crude. It is this latter process, occurring in the enriched-gas displacement, which is treated in this paper. In the past, these phase changes have been determined experimentally and the results incorporated into various modifications of the Buckley-Leverett analysis. Such experimental work is time consuming, and the results are sensitive to numerous experimental errors. With the wide availability of high-speed digital computing equipment and numerous correlations pertaining to the vapor-liquid equilibria of hydrocarbon systems, it is now practical to calculate such phase behavior. This paper describes a computer program for performing these calculations. THE ENRICHED GAS DISPLACEMENT PROCESS Experimental results have shown that oil recovery can be significantly increased by enriching the displacing gas with intermediate hydrocarbon components. The essential features of the phase behavior generated by this enriched-gas-drive process are commonly illustrated with ternary diagrams such as Fig. 1. In this figure, Gas D, which contains a high concentration of intermediate hydrocarbons with respect to the undersaturated Crude A, is injected into the reservoir. When D contacts A, gas goes into solution until the oil becomes saturated (Point. B). Further contacting of Gas D and saturated Oil B results in a Mixture C which separates into Vapor Y(c) and Liquid X(c). Liquid X(c) is contacted by additional Gas D, resulting in Mixture E which separates into Vapor Y(e) and Liquid X(e). Repeated contacts of the liquid by the injected gas will eventually result in Liquid X(d) of maximum enrichment existing in equilibrium with Gas Y(d). The equilibrium tie-line X(d) Y(d), when extended, passes through the Point D representing the enriched injection gas. For systems of more than three components, the predicted equilibrium states are dependent upon not only reservoir temperature and pressure, but also the compositions of the crude oil and injected gas. If the gas is sufficiently enriched, a miscible displacement is generated. Line is tangent to the phase envelope at the critical point (Point Z) and represents the limiting slope of the tie-lines as the critical state is approached. Point I therefore represents the minimum enrichment of injection gas required to generate a miscible displacement. Point G represents the minimum enrichment required for initial miscibility of the injection gas with Crude A.Attra has presented a method to be used for prediction of oil recovery by the enriched gas displacement process. To develop the phase behavior data needed, he designed the experimental procedure described in the following quotation from his paper SPEJ P. 160ˆ


SPE Journal ◽  
2019 ◽  
Vol 24 (03) ◽  
pp. 973-987 ◽  
Author(s):  
Neha Anand ◽  
Brandon Tang ◽  
Bradley (Duong) Nguyen ◽  
Chao-yu Sie ◽  
Marco Verlaan ◽  
...  

Summary Application of thermal and solvent enhanced-oil-recovery (EOR) technologies for viscous heavy-oil recovery in naturally fractured reservoirs is generally challenging because of low permeability, unfavorable wettability and mobility, and considerable heat losses. Vapor/oil gravity drainage (VOGD) is a modified solvent-only injection technology, targeted at improving viscous oil recovery in fractured reservoirs. It uses high fluid conductivity in vertical fractures to rapidly establish a large solvent/oil contact area and eliminates the need for massive energy and water inputs, compared with thermal processes, by operating at significantly lower temperatures with no water requirement. An investigation of the effects of solvent-injection rate, temperature, and solvent type [n-butane and dichloromethane (DCM)] on the recovery profile was performed on a single-fracture core model. This work combines the knowledge obtained from experimental investigation and analytical modeling using the Butler correlation (Das and Butler 1999) with validated fluid-property models to understand the relative importance of various recovery mechanisms behind VOGD—namely, molecular diffusion, asphaltene precipitation and deposition, capillarity, and viscosity reduction. Experimental and analytical model studies indicated that molecular diffusion, convective dispersion, viscosity reduction by means of solvent dissolution, and gravity drainage are dominant phenomena in the recovery process. Material-balance analysis indicated limited asphaltene precipitation. High fluid transmissibility in the fracture along with gravity drainage led to early solvent breakthroughs and oil recoveries as high as 75% of original oil in place (OOIP). Injecting butane at a higher rate and operating temperature enhanced the solvent-vapor rate inside the core, leading to the highest ultimate recovery. Increasing the operating temperature alone did not improve ultimate recovery because of decreased solvent solubility in the oil. Although with DCM, lower asphaltene precipitation should maximize the oil-recovery rate, its higher solvent (vapor)/oil interfacial tension (IFT) resulted in lower ultimate recovery than butane. Ideal density and nonideal double-log viscosity-mixing rules along with molecular diffusivity as a power function of oil viscosity were used to obtain an accurate physical description of the fluids for modeling solvent/oil behavior. With critical phenomena such as capillarity and asphaltene precipitation missing, the Butler analytical model underpredicts recovery rates by as much as 70%.


Sign in / Sign up

Export Citation Format

Share Document