This paper presents the results of numerical simulation of dry, forward combustion tube experiments. The kinetic aspects of in-situ combustion processes also are discussed. The goals of the study are to investigate processes also are discussed. The goals of the study are to investigate the fuel deposition mechanism and to identify the key parameters affecting the performance of in-situ combustion processes. The thermal simulator developed at Gulf R and D Co. was used in the study. It was modified to include the capillary outlet effects for a more realistic description of the oil and water productions. The following experimental data were matched: cumulative water and oil productions, position of the combustion front as a function of time, fuel consumption, position of the combustion front as a function of time, fuel consumption, temperature as a function of time and position, and the pressure drop across the tube. History matches were performed for two crude oils with distinctly different physical properties (gravities of 26.5 and 13 API [0. 896 and 0. 979 g/cm3]). The agreements between experimental data and simulation results were excellent. Results indicate that the component equilibrium K-values and the kinetics of cracking reactions are the most important parameters affecting the fuel deposition, and that the fuel deposition mechanism, the fuel composition, and the locations and sizes of the transient zones depend on the crude oil and reservoir rock properties. Simulation results are always sensitive to the K-values of the light oil component but insensitive to the K-values of the heavy oil component. Results are sensitive to the kinetics of cracking reaction only if the cracking reaction is catalytic or the peak temperature and the fuel consumption are sufficiently high. Furthermore, the fuel available may or may not be solely in the form of coke. Our study suggests that further investigations of the catalytic effect of reservoir rocks and reaction kinetics of the cracking reaction are needed. Also, more than two crude oil components may be required to simulate the evaporation effect of crude oil accurately. Introduction In in-situ combustion processes, many physical changes as well as chemical reactions take place simultaneously or sequentially in the vicinity of the combustion front. It is generally believed that the combustion zone is preceded by a cracking or superheated steam zone, where coke is formed and preceded by a cracking or superheated steam zone, where coke is formed and deposited on the sand grains, and some lighter crude oil components evaporate and move forward with the flowing gas phase. The kinetics of combustion and cracking reactions in the combustion zone and the cracking zone has been discussed widely in the literature. The mechanisms of the physical changes and chemical reactions occurring around the combustion zone can be studied effectively through numerical simulation by using a thermal simulator. Although a number of numerical simulations of combustion tube experiments have been performed with different thermal simulators, no conclusions regarding the mechanism of fuel deposition can be drawn from these studies. The mentioned simulations either neglect the formation of coke from cracking reaction or use a high cracking rate so that no residual oil will be present in the combustion zone. The mechanism of fuel deposition is controlled by two important processes: the evaporation of crude oil components and the kinetics of the processes: the evaporation of crude oil components and the kinetics of the cracking reaction. These two processes determine how much fuel eventually will be burned and how much fuel will be in the form of coke. It has been reported, that low-temperature oxidation can have a significant effect on the fuel deposition and fuel characteristics. However, this reaction is important only when oxygen is available downstream of the combustion front. If oxygen is used completely in a combustion tube experiment, low-temperature oxidation will not play an important role in the fuel deposition mechanism. For a system with a high cracking reaction rate, it is likely that all of the crude oil in the cracking zone will be either evaporated or coked so that coke is the sole source of fuel. However, if the cracking rate is so low that only a portion of crude oil in the cracking zone is evaporated or coked, then some residual crude oil also will be burned in the combustion zone. This is supported strongly by the experimental data of Hildebrand who conducted a number of combustion tube experiments using clean, crushed Berea sandpacks with a variety of crude oils. SPEJ p. 657
A new fully implicit formulation for compositional simulators is presented. Rather than solving for pressure, saturations and phase compositions, the new formulation solves for pressure, overall concentrations and K-values. This change of primary unknowns to be solved improves numerical stability by yielding a more diagonally dominant Jacobian. The use of K-values as primary unknowns also allows the composition constraints to be solved separately from the other equations. The solution of the constraint equations is very efficient because they can be reformulated as monotonic functions*. The primary unknowns are further classified into reduced-unknowns and pivotal-unknowns. The latters are eliminated from the mass conservation equations using the fugacity equalities and constraints. The pivotal-unknowns are selected according to the sensitivity of the equations to the unknowns. This selection conforms to the partial pivoting strategy of Gaussian elimination and enhances numerical stability. Finally, a partial solution method is introduced to eliminate unnecessary calculations of those equations which have met the convergence criteria. For a simple system, it is shown that the new formulation results in a more diagonally dominant Jacobian matrix. Numerical experiments involving one dimensional multicontact-miscibile (MCM) and immiscible problems, a two dimensional MCM problem and a three dimensional MCM problem are conducted to illustrate the capability and efficiency of the new formulation. Results indicate that the number of iterations per time step required by the new formulation is only about 40% – 60% of an existing fully implicit scheme. It is almost as low as that required by an IMPES scheme.
Experimental phase equilibrium data are presented for three reservoir oils at conditions approximating those encountered in in-situ thermal recovery processes. The fluid systems involved consist of three major groups of components: flue gas, water, and crude oil. Data were measured at temperatures from 204.4 to 371.1°C (400 to 700°F) and pressures from 6996.0 to 20785.6 kPa (1,000 to 3,000 psia). Experimental phase equilibrium data were used to develop a correlation of binary interaction coefficients of crude-oil fractions required for the Peng-Robinson equation of state. Phase equilibrium data predicted using the Peng-Robinson equation of state, using our interaction coefficients, are compared with experimental data. Generally, the Peng-Robinson equation of state predictions were in close agreement with the experimental data. Effect of feed gas/oil ratio and water/oil ratio on the equilibrium coefficients was examined through the Peng-Robinson equation of state. A study on the feasibility of representing the crude oil by only two fractions was made also. This study includes a procedure for lumping the crude-oil fractions and examples showing the importance of mixing rules in determining the pseudo critical properties of lumped fractions. Introduction The steady growth of commercial thermal recovery processes1 has created a need for basic data on phase equilibria that involve water and hydrocarbons ranging from methane to high boiling-point fractions. The in-situ thermal recovery processes often are operated at pressures above 6800 kPa (1,000 psia) and temperatures above 200°C (400°F). Experimental data and theoretical correlations on phase equilibria approximating these systems are virtually nonexistent. Early work by White and Brown2 dealt with high boiling-point hydrocarbon phase equilibria. However, the highest pressure studied was 6894.8 kPa (1,000 psia) and the lightest component was pentane. Poettmann and Mayland,3 on the basis of an empirical correlation,4 constructed charts of equilibrium coefficients, or K values, as functions of pressure and temperature for various boiling-point fractions. But the maximum pressure studied was 6894.8 kPa (1,000 psia). Later, Hoffmann et al.5 studied phase behavior of a gas-condensate system with the highest pressure reaching 20 684.3 kPa (3,000 psia) but the highest temperature investigated was only 94.2°C (201°F). In 1963, Grayson and Streed6 reported experimental vapor/liquid equilibrium data for high-temperature and high-pressure hydrocarbon systems. They also extended the Chao-Seader correlation to cover the higher temperature ranges. However, the. major light component in Grayson and Streed's system was hydrogen. Recently, because of the increasing activity in carbon dioxide flooding processes, the phase equilibria of systems involving carbon dioxide and crude oil has received attention. Simon et al.7 studied phase behavior and other properties of carbon-dioxide/reservoir-oil systems. Shelton and Yarborough8 examined phase behavior in porous media during carbon dioxide or rich-gas flooding. No extensive data on equilibrium coefficients were reported in those papers, and the temperature ranges (out of physical reality) were below 93.5°C (200°F). None of these papers surveyed included water as a component.
A new mathematical model is developed for analyzing single-well pulse test data from a reservoir consisting of two permeable zones separated by a less permeable zone. The model considers well-bore storage and skin effects and allows the three zones to have different rock properties. The line-source assumption used in developing the analytic solution is validated via comparison of results with another model also developed in this work. This latter model assumes identical rock properties for the two permeable zones but allows the wellbore radius to be finite. The numerical Laplace inversion methods available in the literature are evaluated with respect to some simple functions and the analytic Laplace space solution of the problem of interest. It is found that the Fourier-series and Zakian's methods can serve as alternatives to the well-known Stehfest algorithm. The effects of wellbore storage and skin on single-well pulse testing of layered reservoirs are investigated. Consistent with the findings in the literature for multiwell pulse testing, it is shown that the influence of wellbore storage and skin is to reduce the response amplitude and to increase the time lag of the pulse response. Among the wellbore storage and skin factors involved in a single well pulse testing, the skin associated with the observation interval is least influential on the pressure response. Using field data, it is shown that the mathematical model presented can be used, in conjunction with an optimization subroutine, to estimated the parameters of interests. It is further shown numerically that neglecting the effects of wellbore storage and skin could result in an underestimation of the formation permeabilities of interest. The problem of non-uniqueness of solution for parameter estimation using optimization techniques is also discussed for the subject problem.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.