Simulation of Steamflooding With Distillation and Solution Gas

1976 ◽  
Vol 16 (05) ◽  
pp. 235-247 ◽  
Author(s):  
K.H. Coats

Abstract This paper describes a three-dimensional numerical model for simulating steam-injection processes. The model accounts for solution gas and steam distillation of oil. The relative-permeability treatment presented includes a flexible but simple representation of temperature dependence and a history-dependent hysteresis in gas relative permeability. Since computational stability is a major difficulty in steamflood simulation, an implicit treatment of transmissibilities and capillary pressure is presented in detail. Model applications include comparisons with laboratory data, sensitivity experiments, and a field steam-injection test. Introduction Shutler and Abdalla and Coats described two-dimensional, three-phase flow numerical models for simulating steam-injection processes. Weinstein et al. described a one-dimensional model that accounted for steam distillation of oil. Coats et al. described a three-dimensional steamflood model that neglected steam distillation of oil, release of solution gas at elevated temperatures, and temperature dependence of relative permeability. This paper describes an extended formulation that includes these three phenomena and uses a more implicit treatment of capillary pressures and transmissibilities in the fluid-saturation calculations. The extended formulation represents a step toward a fully compositional thermal model without incurring the computational expense of the latter. The relative-permeability treatment described includes a rather flexible but simple representation of temperature dependence and incorporates a hysteresis in gas-phase relative permeability that varies with the historical maximum grid-block gas saturation. The phase-behavior representation is the weakest element of this work. We have found insufficient data relative to PVT behavior of a heavy-oil/steam system to justify sophisticated schemes of the type used in isothermal hydrocarbon systems. The PVT treatment presented is the simplest we could construct subject to the objectives of "directional correctness," reasonable quantitative accuracy, and ability to obtain required parameters from laboratory data either normally parameters from laboratory data either normally available or readily determinable. Model results presented include a comparison with laboratory data for a steamflood of a distillable oil; sensitivity results indicating effects and relative importance of various types of input data; and a comparison between calculated and observed injection rates for a Cold Lake (Alta.) steam-injection test. The latter is of interest in regard to reservations we have had regarding a model's ability to predict steam-injection rates into virtually immobile oil (100,000 cp). The field-test data showed initial and sustained steam-injection rates of 1,400 STB/D (cold-water equivalent). We discuss several reservoir-fluid parameters that had little effect and one independently measured parameter that had a pronounced effect on the calculated injection rate. pronounced effect on the calculated injection rate. MODEL DESCRIPTION The model consists and sewn equations expressing conservation of energy, conservation of mass, and constraints on sums of liquid and gas phase mol fractions. The mass-conservation equations apply to water and to each of three hydrocarbon components. In finite-difference form these equations are the following. SPEJ P. 235

1981 ◽  
Vol 21 (02) ◽  
pp. 162-178 ◽  
Author(s):  
Y.C. Yortsos ◽  
G.R. Gavalas

Abstract This paper deals with the derivation of upper bounds for the growth of the steam zone in steam injection processes for one- or multidimensional reservoirs at constant or variable injection rates. The bounds are derived from the integral balances describing a reservoir of arbitrary geometry by introducing lower bounds for the heat losses to the surrounding area and the hot liquid zone. In this way, the effect of preheating in the hot liquid zone is estimated to determine the recovery efficiency of a steam drive. The growth rate of a one-dimensional steam zone at variable injection rates is subject to two upper bounds resulting from the total thermal energy and the latent heat balances, respectively. Each of the bounds controls the rate of growth of the steam zone in a certain time interval, depending on the dominant mode of heat transfer in the hot liquid zone. At constant injection rates, the steam zone growth at large times is controlled by the bound based on the latent heat balance. This balance depends on a dimensionless parameter, F, defined as the ratio of the latent heat to the total heat injected. Based on the relative magnitude of F with respect to the critical value F= 2/pi, the region of validity of the Marx-Langenheim solution is delineated on a Ts vs. fs diagram. The Marx-Langenheim solution is satisfactory at large times when F greater than 2/pi and becomes less satisfactory as F assumes smaller values. Similar upper bounds are obtained for a two-dimensional steam drive (thin reservoirs). In three-dimensional reservoirs, on the other hand, bounds are derived only for a special form of displacement (separable front). These bounds depend on the models for the steam front shape, K can be determined in terms of the physical variables of the process. Introduction Injection of steam (steamflood or steam drive) is an important thermal recovery method that is applied on a commercial scale in many parts of the world. The main elements of continuous steam injection, as a displacement process, were analyzed thoroughly by experimental studies under both laboratory and field conditions. Along with laboratory and field tests, mathematical models are sought to aid in understanding and designing the process. The engineering evaluation of a steam drive often is based on a simplified mathematical description of reservoir heating by hot fluid injection presented by Marx and Langenheim and subsequently modified by Mandl and Volek. This theory was combined further with simple fluid flow considerations to determine the oil recovery rates in one-dimensional reservoirs. To account for the important effects of gravity override in three-dimensional geometries, Neuman and van Lookeren (following different approaches) derived simple analytical formulas for the calculation of the performance of a steamflood in three-dimensional reservoirs. An increasing number of investigators also have concentrated on the development of reliable numerical models. Three-phase numerical simulators were derived by Shutter for one- and two-dimensional flow; by Abdalla and Coats for two-dimensional flow; and by Coats et al., Coats, and Weinstein et al. for three-dimensional flow. The last two models also account for steam distillation of oil. SPEJ P. 162^


2013 ◽  
Vol 860-863 ◽  
pp. 1357-1361
Author(s):  
Zhou Rong Zhang ◽  
Qian Wang ◽  
Zhi Xia He ◽  
Li Ming Dai

As a study of natural gas engine, three-dimensional numerical simulations of diesel injection rates were conducted by using AVL FIRE code. And particularly the effects of the injection rates in the pilot diesel jets on the engine combustion and emission were analyzed. It is found that injection rate in the pilot diesel jets plays an important role on the emission, and slow ending rate of the injection keeps lower NO and Soot.


2017 ◽  
pp. 63-67
Author(s):  
L. A. Vaganov ◽  
A. Yu. Sencov ◽  
A. A. Ankudinov ◽  
N. S. Polyakova

The article presents a description of the settlement method of necessary injection rates calculation, which is depended on the injected water migration into the surrounding wells and their mutual location. On the basis of the settlement method the targeted program of geological and technical measures for regulating the work of the injection well stock was created and implemented by the example of the BV7 formation of the Uzhno-Vyintoiskoe oil field.


1982 ◽  
Vol 14 (3) ◽  
pp. 33-39
Author(s):  
C Y Kuo

An existing, three-dimensional, Eulerian-Lagrangian finite-difference model was modified and used to examine the far-field transport processes of dumped sewage sludge in the New York Bight. Both in situ and laboratory data were utilized in an attempt to approximate model inputs such as mean current speed, vertical and horizontal diffusion coefficients, particle size distributions, and specific gravities. Concentrations of the sludge near the sea surface predicted from the computer model were compared qualitatively with those remotely sensed.


2021 ◽  
Vol 503 (3) ◽  
pp. 3309-3325
Author(s):  
Sabine Bellstedt ◽  
Aaron S G Robotham ◽  
Simon P Driver ◽  
Jessica E Thorne ◽  
Luke J M Davies ◽  
...  

ABSTRACT We analyse the metallicity histories of ∼4500 galaxies from the GAMA survey at z < 0.06 modelled by the SED-fitting code ProSpect using an evolving metallicity implementation. These metallicity histories, in combination with the associated star formation histories, allow us to analyse the inferred gas-phase mass–metallicity relation. Furthermore, we extract the mass–metallicity relation at a sequence of epochs in cosmic history, to track the evolving mass–metallicity relation with time. Through comparison with observations of gas-phase metallicity over a large range of redshifts, we show that, remarkably, our forensic SED analysis has produced an evolving mass–metallicity relationship that is consistent with observations at all epochs. We additionally analyse the three-dimensional mass–metallicity–SFR space, showing that galaxies occupy a clearly defined plane. This plane is shown to be subtly evolving, displaying an increased tilt with time caused by general enrichment, and also the slowing down of star formation with cosmic time. This evolution is most apparent at lookback times greater than 7 Gyr. The trends in metallicity recovered in this work highlight that the evolving metallicity implementation used within the SED-fitting code ProSpect produces reasonable metallicity results over the history of a galaxy. This is expected to provide a significant improvement to the accuracy of the SED-fitting outputs.


Author(s):  
Tianyu Jin ◽  
Yu Sun ◽  
Chuqiao Wang ◽  
Adams Moro ◽  
Xiwen Wu ◽  
...  

Abstract The stringent emission regulations diesel engines are required to meet has resulted in the usage of multi-hole and ultra-multi-hole injectors, nowadays. In this research study, a double layered 8-hole diesel injection nozzle was investigated both numerically and experimentally. A three-dimensional model of the nozzle which was validated with experimental results was used to analyze the injection characteristics of each hole. The validation was conducted by comparing experiment and simulation injection rate results, acquired simultaneously from all the holes of the injector and the model. The fuel flow rates of the lower layered holes are higher than those of the upper layered holes. Two different needle eccentricity models were established. The first model only included the lateral displacement of the needle during needle lift. The needle reached maximum displacement at full needle lift. The second model considered the needle inelastic deformation into consideration. The needle radially displaces and glides along with the needle seat surface during needle lift. When the eccentricity reached maximum in the radial direction, the needle began to lift upwards vertically. The differences in injection characteristics under the different eccentricity models were apparent. The results indicated that the cycle injection quantity, fuel injection rate and cavitation of each hole were affected during the initial lifting stages of the needle lift. As the eccentricity of the needle increases, the injection rate uniformity from the nozzle hole deteriorates. The result showed that the upper layered holes were affected by the needle eccentricity during needle lift.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Gen-ichiro Yamamoto ◽  
Atsushi Kyono ◽  
Satoru Okada

AbstractMineral trapping through the precipitation of carbonate minerals is a potential approach to reduce CO2 accumulation in the atmosphere. The temperature dependence of amorphous magnesium carbonate (AMC), a precursor of crystalline magnesium carbonate hydrates, was investigated using synchrotron X-ray scattering experiments with atomic pair distribution function (PDF) and X-ray absorption fine structure analysis. PDF analysis revealed that there were no substantial structural differences among the AMC samples synthesized at 20, 60, and 80 °C. In addition, the medium-range order of all three AMC samples was very similar to that of hydromagnesite. Stirring in aqueous solution at room temperature caused the AMC sample to hydrate immediately and form a three-dimensional hydrogen-bonding network. Consequently, it crystallized with the long-range structural order of nesquehonite. The Mg K-edge X-ray absorption near-edge structure spectrum of AMC prepared at 20 °C was very similar to that of nesquehonite, implying that the electronic structure and coordination geometry of Mg atoms in AMC synthesized at 20 °C are highly similar to those in nesquehonite. Therefore, the short-range order (coordination environment) around the Mg atoms was slightly modified with temperature, but the medium-range order of AMC remained unchanged between 20 and 80 °C.


Sign in / Sign up

Export Citation Format

Share Document