Jin, Luchao (University of Oklahoma) | Budhathoki, Mahesh (University of Oklahoma) | Jamili, Ahmad (University of Oklahoma) | Li, Zhitao (The University of Texas at Austin) | Luo, Haishan (The University of Texas at Austin) | Delshad, Mojdeh (The University of Texas at Austin) | Shiau, Ben (University of Oklahoma) | Harwell, Jeffrey H. (University of Oklahoma)
The surfactant screening process to develop an optimum formulation under reservoir conditions is typically time consuming and expensive. Theories and correlations like HLB, R-ratio and packing parameters have been developed. But none of them can quantitatively consider both the effect of oil type, salinity, hardness and temperature, and model microemulsion phase behavior.
This paper uses the physics based Hydrophilic Lipophilic Difference (HLD) Net Average Curvature (NAC) model, and comprehensively demonstrated its capabilities in predicting the optimum formulation and microemulsion phase behavior based on the ambient conditions and surfactant structures. By using HLD equation and quantitatively characterized parameters, four optimum surfactant formulations are designed for target reservoir with high accuracy compared to experimental results. The microemulsion phase behavior is further predicted, and well matched the measured equilibrium interfacial tension. Its predictability is then reinforced by comparing to the empirical Hand's rule phase behavior model. Surfactant flooding sandpack laboratory tests are also interpreted by UTCHEM chemical flooding simulator coupled with the HLD-NAC phase behavior model.
The results indicate the significance of HLD-NAC equation of state in not only shorten the surfactant screening processes for formulators, but also predicting microemulsion phase behavior based on surfactant structure. A compositional reservoir simulator with such an equation of state will increase its predictability and hence help with the design of surfactant formulation.
Ghasemi, M. (Petrostreamz AS) | Astutik, W. (Petrostreamz AS) | Alavian, S. A. (PERA AS) | Whitson, C. H. (PERA AS/NTNU) | Sigalas, L. (Geological Survey of Denmark and Greenland) | Olsen, D. (Geological Survey of Denmark and Greenland) | Suicmez, V. S. (Maersk Oil & Gas A/S)
This paper presents a novel technique to determine multi-component diffusion coefficients for CO2 injection in a North Sea Chalk Field (NSCF) at reservoir conditions. Constant volume diffusion (CVD) method is used, consisting of an oil-saturated chalk core in contact with an overlying free-space, which is filled with the CO2. The experimental data are matched with an EOS-based compositional model.
Transport by diffusion controls the dynamics of the constant-volume system, together with phase equilibria, allowing a consistent estimation of diffusion coefficients needed to describe the observed changes in system pressure.
We conduct two experiments at reservoir condition: one utilizes a core plug saturated with live-oil, and the other with stock tank oil (STO). Once the experiments are completed, EOS-based compositional simulation is performed to match the experimental data using the oil and gas diffusion coefficients as history matching parameters. The modeling work is conducted with a commercial reservoir simulator using a two dimensional radial grid model to describe the experimental setup.
The experiment utilizes a vertically-oriented core holder with a height of 92 mm and 37.6 mm in diameter. An outcrop chalk core with a sealing sleeve is mounted in the core holder, which has the same diameter and a height of 64.6 mm, thus resulting in an overlying void space. The system is initially saturated with oil at reservoir condition. CO2 is then injected from the top, forming an overlying CO2 chamber, and displacing oil towards the bottom of the core holder. Once CO2 fills the overlying bulk space, the system is isolated with no further injection or production.
The CO2 and oil reach and remain in equilibrium locally at the gas-oil interface throughout the test, initiating and maintaining the diffusion mechanism. Diffusion of CO2 into the oil results in a decreasing pressure, which is the main history matching parameter.
The multi-component diffusion coefficients are found to match the model pressure-time prediction to the experimental data. This suggests the modelling workflow incorporates a representative EOS model and the main transport dynamics controlled by diffusion are being treated properly.
The two main challenges in the modeling are (1) the limitation on setting an appropriately-high permeability for the CO2 chamber, and (2) the reservoir simulator neglects compositional dependency of diffusion coefficients.
Proper simulation of CO2 injection in fractured chalk reservoirs requires the ability to model multi-component diffusion accurately. The proposed CVD-method provides such modeling capabilities. Our modeling and experimental work indicate the novelty of the CVD method to determine the diffusion coefficients of a system where diffusion is the dominant displacement mechanism. The fact that the oil is contained within a low-permeability chalk sample reduces density-driven convection that could result due to non-monotonic oil density changes as CO2 dissolves into the oil.
The use of isenthalpic flash has become of interest for the simulation of some heavy oil recovery processes where large temperature changes are experienced. For these thermal simulations energy can be used as a primary variable. This leads to thousands or millions of individual multiphase isenthalpic flash calculations. Robust and efficient algorithms for multiple-phase isenthalpic flash are required to improve the efficiency of compositional simulations for thermal recovery.
The general framework on state function based flash specifications proposed by
Narrow boiling mixtures can be dealt with in the majority of cases without any significant difficulty. This is true of the direct substitution algorithm and the proposed solution procedure. The vast majority of examples can be solved without using Q function maximisation. The challenges associated with multiphase calculations in the Newton steps are investigated. In particular, inadequate initial estimate of the equilibrium type may lead to non-convergent iteration. This can usually be solved by introduction of a new phase and/or elimination of an existing phase. The speed of the method is analysed for a large number of specifications and is found to be only slightly more expensive than isothermal flash in the majority of cases.
Compositional reservoir simulation plays a vital role in the development of conventional and unconventional reservoirs. Two major building blocks of compositional simulation are phase behavior and fluid transport computations. The oil and gas reserves and flow of reservoir fluids are strongly dependent on phase behavior. In conventional reservoirs, capillary pressure is relatively small and is typically ignored in phase behavior calculations. However, large capillary pressure values are encountered in tight formations such as shales; and therefore, its effects should not be ignored in phase equilibria calculations. Neglecting the effects of capillary pressure on phase behavior can lead to an inaccurate estimation of original oil and gas in place as well as recovery performance. In spite of this, the effect of capillary pressure on phase behavior in tight reservoirs has not been well studied using compositional simulation, especially for hydraulically-fractured reservoirs.
In this paper, we develop a new compositional reservoir simulator capable of modeling discrete fractures and incorporating the effect of capillary pressure on phase behavior. Large-scale natural and hydraulic fractures in tight rocks and shales are modeled with a technique called Embedded Discrete Fracture Model (EDFM) where fractures are modeled explicitly without using local grid refinement or an unstructured grid. Flow of hydrocarbons occurs simultaneously within similar and different porosity types. Capillary pressure is considered in both flow and flash calculations, where simulations also include variable pore size as a function of gas saturation in each grid block. We examine the impact of capillary pressure on the original oil in place and cumulative oil production for different initial reservoir pressures (above and below the bubble-point pressure) on Bakken and Eagle Ford fluids. The importance of capillary pressure on both flow and flash calculations from hydraulically fractured horizontal wells during primary depletion in fractured tight reservoirs using Bakken fluid composition is demonstrated.
Phase behavior calculations show that bubble-point pressure is suppressed allowing the production to remain in the single-phase region for a longer period of time and altering phase compositions and fluid properties such as density and viscosity of equilibrium liquid and vapor. The results show that bubble-point suppression is larger in the Eagle Ford shale than for Bakken. When capillary pressure is considered, we found an increase in original oil in place up to 4.1% for Bakken and 46.33% for the Eagle Ford crude. Depending on the initial reservoir pressure, cumulative primary production after one year increases owing to capillary pressure by approximately 9.0 – 38.2% for Bakken oil and 7.2 – 154% for Eagle Ford oil. The recovery increase caused by capillary pressure becomes more significant when reservoir pressure is far below bubble-point pressure. The simulation results with hydraulically fractured wells give similar recovery differences; cumulative oil production after 1 year is 3.5 – 5.2% greater when capillary pressure is considered in phase behavior calculations for Bakken.