Tie-Simplex-Based Phase-Behavior Modeling in an IMPEC Reservoir Simulator

SPE Journal ◽  
2013 ◽  
Vol 19 (02) ◽  
pp. 327-339 ◽  
Author(s):  
M.. Rezaveisi ◽  
K.. Sepehrnoori ◽  
R.T.. T. Johns

Summary Recently, tie-simplex-based phase-behavior modeling in reservoir simulators has been applied and investigated as a potential method for improving the computational speed of equation-of-state (EOS) -based reservoir simulators. We implemented compositional-space adaptive tabulation (CSAT), the most promising tie-simplex-based method, in UTCOMP, the University of Texas' in-house IMPEC compositional reservoir simulator, to investigate its computational efficiency compared with the phase-behavior algorithm in UTCOMP. The results show that applying CSAT only to skip stability analysis does improve computational time, but only when a significant portion of the gridblocks are in the single-phase region and no other technique for avoiding stability analysis is used. However, in most cases, there is little or no computational advantage to use of CSAT when the simple option in UTCOMP is used where stability analysis is skipped for blocks surrounded by single-phase regions. We explore in detail the performance of CSAT, which depends significantly on the specific gas flood modeled, and the number of tie-lines generated during adaptive tabulation. The results shed light on applicability of CSAT in the IMPEC-type compositional reservoir simulators and show that the advantages of CSAT in this type of simulator are not as great as are reported in the literature for fully implicit or adaptive implicit formulations.

1978 ◽  
Vol 18 (01) ◽  
pp. 20-26 ◽  
Author(s):  
Ralph Simon ◽  
A. Rosman ◽  
Erdinc Zana

February 1978 Original manuscript received in Society of Petroleum Engineers office Jan. 14, 1977. Paper accepted for publication Aug. 15, 1977. Revised manuscript received Sept. 21, 1977. Paper (SPE 6387) was presented at the SPE-AIME Permian Basin Oil and Gas Recovery Conference, held in Midland, Tex., March 10-11, 1977. Abstract This paper presents experimental phase behavior data on two CO2-reservoir oil systems at reservoir pressures and temperatures. pressures and temperatures. The data includepressure-composition diagrams with bubble points, dew points, and critical points;vapor-liquid equilibrium compositions and related K values;vapor and liquid densities compared with values calculated by the Redlich-Kwong equation of state;vapor and liquid viscosities compared with predictions by the Lobrenz-Bray-Clark correlation; andinterfacial tensions for six vapor-liquid mixtures compared with values calculated by the Weinaug-Katz parachor equation. These and other published data contribute to development of the generalized correlations needed by reservoir and production engineers for evaluating, designing, and efficiently operating CO2-injection projects. projects Introduction This paper presents experimental phase behavior data for two CO2-reservoir oil systems. These data are used in predicting the performance of CO2 floods with a compositional simulator. The simulator calculates vapor and liquid compositions, densities, viscosities, and interfacial tensions to describe the phase behavior as the injected CO2 advances through phase behavior as the injected CO2 advances through the reservoir. The simulator predictions are used to evaluate proposed projects and to design and efficiently operate approved ones. The data in this paper consist of pressure-composition diagrams with bubble points, pressure-composition diagrams with bubble points, dew points, and critical points; and compositions, densities, viscosities, and interfacial tensions of vapors and liquids in equilibrium in the two-phase region. These data were obtained by the experimental procedure shown in Fig. 1. procedure shown in Fig. 1. We have compared our measured data with values calculated by existing methods: Redlich-Kwong equation for densities, Lohrenz-Bray-Clark correlation for viscosities, and the Weinaug-Katz parachor equation for interfacial tension. We found parachor equation for interfacial tension. We found that these published methods give acceptable agreement in some areas, but in general, they are not satisfactory for engineering purposes. Therefore, we conclude that improved calculation methods are needed for CO2 systems. For the special case of compositional simulator applications, we devised a technique for obtaining satisfactory calculated density, viscosity, and interfacial tension values. This technique is discussed in the section on "Measurements vs Calculations." We believe that our data, along with previously published information and information yet to come, published information and information yet to come, will advance the development of satisfactory correlations, thus reducing the need for extensive laboratory studies of individual systems. PRESSURE-COMPOSITION DIAGRAMS PRESSURE-COMPOSITION DIAGRAMS OIL A Ten mixtures of CO2 and Reservoir Oil A were prepared. These mixtures contained CO2 concentrations prepared. These mixtures contained CO2 concentrations of 0, 20, 40, 55, 60, 65. 70, 75, 80, and 90 mol percent. At 130 degrees F, pressure traverses were made with each mixture. These traverses started in the single-phase region at a pressure above the bubble (or dew) points and lowered the pressure in discrete steps, passing from the single-phase into the two-phase region. At each step, the vapor and liquid volumes were measured. The results are described in Fig. 2A. At 130 degrees F, the critical point of the CO2-Reservoir Oil A system (where intensive properties of the gas and liquid phases were equal) properties of the gas and liquid phases were equal) is 2,570 psia and 60-mol percent CO2. OIL B Eight mixtures of CO2 and Reservoir Oil B also were prepared and studied in the visual cell at 255 degrees F. CO2 concentrations for these mixtures were 0, 20, 40, 55, 65, 75, 80, and 85 mol percent. The pressure was varied from 800 to 6,100 psia, and the pressure was varied from 800 to 6,100 psia, and the relative vapor and liquid volumes measured. The results are given in Fig. 2B. The critical point of the CO2-Reservoir Oil B system at 255 degrees F is 4,890 psia and 74-mol percent CO2. psia and 74-mol percent CO2. SPEJ P. 20


SPE Journal ◽  
2011 ◽  
Vol 16 (04) ◽  
pp. 751-767 ◽  
Author(s):  
R.. Okuno ◽  
R.T.. T. Johns ◽  
K.. Sepehrnoori

Summary CO2 floods at temperatures typically below 120°F can involve complex phase behavior, where a third CO2-rich liquid (L2) phase coexists with the oleic (L1) and gaseous (V) phases. Results of slimtube measurements in the literature show that an oil displacement by CO2 can achieve high displacement efficiency of more than 90% when three hydrocarbon phases coexist during the displacement. However, the mechanism for the high-displacement efficiency is uncertain because the complex interaction of phase behavior with flow during the displacement is not fully understood. In this paper, we present the first detailed study of three-phase behavior predictions and displacement efficiency for low-temperature CO2 floods. Four-component EOS models are initially used to investigate systematically the effects of pressure, temperature, and oil properties on development of three-phase regions and displacement efficiency. Multicomponent oil displacements by CO2 are then considered. We use a compositional reservoir simulator capable of robust three-phase equilibrium calculations. Results show that high displacement efficiency of low-temperature CO2 floods is a consequence of both condensing and vaporizing behavior. The L2 phase serves as a buffer between the immiscible V and L1 phases within the three-phase region. Components in the L1 phase first transfer efficiently to the L2 phase near a lower critical endpoint (LCEP). These oil components then transfer to the V phase near an upper critical endpoint (UCEP) at the trailing edge of the three-phase region. The CEPs are defined where two of the three coexisting phases merge in the presence of the other immiscible phase. Unlike two-phase displacements, condensation and vaporization of intermediate components occur simultaneously within the three-phase region. The simultaneous condensing/vaporizing behavior involving the CEPs is also confirmed for simulations of several west Texas oil displacements. Quaternary fluid models can predict qualitatively the complex displacements because four is the minimum number of components to develop CEP behavior in composition space at a fixed temperature and pressure.


1982 ◽  
Vol 22 (01) ◽  
pp. 28-36 ◽  
Author(s):  
M. Bourrel ◽  
C. Chambu ◽  
R.S. Schechter ◽  
W.H. Wade

Abstract Surfactant/oil/water phase diagrams have become the most important screening tool used to select microemulsion systems for enhanced oil recovery. The number of phases coexisting at a given salinity, the extent of the single-phase region, and the position of the phase boundaries all have relevance with respect to oil displacement efficiency. It is shown that the phase diagrams can be made to take on different configurations depending on the alcohol cosurfactant, the salinity, the impurities present in the surfactant, and the dispersity of the surfactant mixture. Besides the importance of the phase boundary shape, this study provides further insight into factors determining the height of the binodal surface on the pseudoternary phase diagram. Results show the effect of salinity as well as the surfactant, alcohol, and hydrocarbon types on the height of the binodal surface. It is shown that salinity is the main factor; other parameters have little or no influence once a surfactant has been selected. Finally the microemulsion viscosity is shown to be related to the proximity of the formulation to phase boundaries. Extensive data for one system are presented. Introduction It is now recognized that formulating surfactant/oil/brine systems that exhibit desirable phase behavior is an important step in optimizing performance of microemulsion systems for enhanced oil recovery. Oil is displaced by a combination of mechanisms-miscible displacement, swelling of the oil phase, and low tension displacement all of which are related to the topology of the phase boundaries in composition space. To predict the outcome of a particular project, a representation of the phase boundaries and their evolution when diluted with oil or brines having various proportions of divalent ions is required. For example, successful application of the salinity gradient concept demands phase relationships specially structured to accommodate the variations in salinity experienced by the surfactant slug during the course of the flood. Recent publications have dealt with the optimal salinity as a function of total amphiphile concentration (surfactant plus cosurfactant), and reported trends that are quite different from those found if the cosurfactant (alcohol) concentration is held constant. One purpose of this paper is to demonstrate that contorted phase boundaries found by Glover et al are caused by the variation of alcohol concentration when the concentration of total amphiphile is varied and because the direction that the phase boundaries twist or rotate is controlled by the nature of the alcohol. Another important factor is the extent of the single-phase region. More precisely, the height of the demixing curve in the pseudoternary representation should be minimized. This would permit, in principle, the amount of surfactant and cosurfactant in the micellar slug to be minimized. A correlation permitting the determination of the oil, salinity, alcohol, and surfactant at which the height of the demixing curve is minimized has been reported, but few data giving the value of the minimum height have been presented. This height is an important feature of the phase boundary topology and extensive measurements are reported here. The microemulsion viscosity must be high enough to help maintain mobility control. It is sometimes difficult to achieve the required levels of viscosity. Studies of microemulsion viscosity have been reported. We provide further data here and have related the microemulsion viscosities to phase behavior. Materials and Experimental Techniques The phase diagrams have been established by two techniques: a titration procedure and a grid-point technique. SPEJ P. 28^


Author(s):  
Md. Hamidul Kabir ◽  
Ravshan Makhkamov ◽  
Shaila Kabir

The solution properties and phase behavior of ammonium hexylene octyl succinate (HOS) was investigated in water and water-oil system. The critical micelle concentration (CMC) of HOS is lower than that of anionic surfactants having same carbon number in the lipophilic part. The phase diagrams of a water/ HOS system and water/ HOS/ C10EO8/ dodecane system were also constructed. Above critical micelle concentration, the surfactant forms a normal micellar solution (Wm) at a low surfactant concentration whereas a lamellar liquid crystalline phase (La) dominates over a wide region through the formation of a two-phase region (La+W) in the binary system. The lamellar phase is arranged in the form of a biocompatible vesicle which is very significant for the drug delivery system. The surfactant tends to be hydrophilic when it is mixed with C10EO8 and a middle-phase microemulsion (D) is appeared in the water-surfactant-dodecane system where both the water and oil soluble drug ingredient can be incorporated in the form of a dispersion. Hence, mixing can tune the hydrophile-lipophile properties of the surfactant. Key words: Ammonium hexylene octyl succinate, mixed surfactant, lamellar liquid crystal, middle-phase microemulsion. Dhaka Univ. J. Pharm. Sci. Vol.3(1-2) 2004 The full text is of this article is available at the Dhaka Univ. J. Pharm. Sci. website


2007 ◽  
Vol 539-543 ◽  
pp. 2928-2933 ◽  
Author(s):  
V.S. Yusupov ◽  
A.I. Milyaev ◽  
Galia F. Korznikova ◽  
Alexander V. Korznikov ◽  
J.K. Kovneristii

Results of experimental research into evolution of the structure and microhardness of the hard magnetic Fe-30Cr-8Co-0,7Ti-0,5V-0,7Si alloy during complex two-level loading (compression + torsion) in isothermal conditions at various temperatures in single-phase region are reported. It was revealed that the deformation leads to a strong refinement of initial coarse-grained structure in the whole volume of the sample, however the generated structure is non-uniform through the body of the sample. In an active zone of deformation, near to mobile head, there is a microcrystalline layer with a grain size of about 5 microns which thickness poorly depends on the formation. With removal from the active zone of deformation the grain size increases, and microhardness decreases.


2021 ◽  
Author(s):  
Daniele Pampaloni ◽  
Antonio Andreini ◽  
Alessandro Marini ◽  
Giovanni Riccio ◽  
Gianni Ceccherini

Abstract Thermoacoustic characterization of gas turbine combustion systems is of primary importance for successful development of gas turbine technology, to meet the stringent targets on pollutant emissions. In this context, it becomes more and more necessary to develop reliable tools to be used in the industrial design process. The dynamics of a lean-premixed full-annular combustor for heavy-duty applications has been numerically studied in this work. The well-established CFD-SI method has been used to investigate the flame response varying operational parameters such as the flame temperature (global equivalence ratio) and the fuel split between premixed and pilot fuel injections: such a wide range experimental characterization represents an opportunity to validate the employed numerical methods and to give a deeper insight into the flame dynamics. URANS simulations have been performed, due to their affordable computational costs from the industrial perspective, after validating their accuracy through the comparison against LES results. Furthermore, an approach where the pilot and the premixed flame responses are analyzed separately is proposed, exploiting the independence of their evolution. The calculated FTFs have been implemented in a 3D FEM model of the chamber, in order to perform linear stability analysis and to validate the numerical approach. A boundary condition for rotational periodicity based on Bloch-Wave theory has been implemented into the Helmholtz solver and validated against full-annular chamber simulations, allowing a significant reduction in computational time. The reliability of the numerical procedure has been assessed through the comparison against full-annular experimental results.


Sign in / Sign up

Export Citation Format

Share Document