Next Article in Journal
Homodyne Spectroscopy with Broadband Terahertz Power Detector Based on 90-nm Silicon CMOS Transistor
Next Article in Special Issue
Effect of Capillary Number on the Residual Saturation of Colloidal Dispersions Stabilized by a Zwitterionic Surfactant
Previous Article in Journal
Stability and Manoeuvrability Simulation of a Semi-Autonomous Submarine Free-Running Model SUBOFF with an Autopilot System
Previous Article in Special Issue
The Effect of Carbonyl and Hydroxyl Compounds Addition on CO2 Injection through Hydrocarbon Extraction Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compositional Modeling of Dimethyl Ether–CO2 Mixed Solvent for Enhanced Oil Recovery

1
Department of Earth Resources and Environmental Engineering, Hanyang University, 222 Wangsimni-ro, Seongdong-gu, Seoul 04763, Korea
2
Center for Climate/Environment Change Prediction Research, Ewha Womans University, 52 Ewhayeodae-gil, Seodaemun-gu, Seoul 03760, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(1), 406; https://doi.org/10.3390/app11010406
Submission received: 16 December 2020 / Revised: 26 December 2020 / Accepted: 28 December 2020 / Published: 4 January 2021
(This article belongs to the Special Issue Advanced Oil Recovery Technologies)

Abstract

:
Dimethyl ether (DME) is a compound first introduced by Shell as a chemical solvent for enhanced oil recovery (EOR). This study aims to investigate the efficiency of EOR using the minimum miscible pressure (MMP) and viscous gravity number when a mixed solvent of CO2 and DME is injected. Adding DME to the CO2 water-alternating-gas process reduces the MMP and viscous gravity number. Reduction in MMP results in miscible conditions at lower pressures, which has a favorable effect on oil swelling and viscosity reduction, leading to improved mobility of the oil. In addition, the viscous gravity number decreases, increasing the sweep efficiency by 26.6%. Numerical studies were conducted through a series of multi-phase, multi-component simulations. At a DME content of 25%, the MMP decreased by 30.1% and the viscous gravity number decreased by 66.4% compared with the injection of CO2 only. As a result, the maximum oil recovery rate increased by 31% with simultaneous injection of DME and CO2 compared with only using CO2.

1. Introduction

As the discovery of new oil fields becomes increasingly difficult, most oil companies are focusing on maintaining economic oil prices through enhanced oil recovery (EOR) technology to maximize the recovery from oil fields [1]. Solvent injection is a mature EOR method [2]. CO2 and hydrocarbon gas are widely used as injection fluids [3], wherein their contact with reservoir fluid causes viscosity reduction and oil swelling. The effect is greater under miscible conditions than immiscible conditions [4]. Research using hydrocarbon gas or CO2 as a solvent has been conducted, from the laboratory to the field, for decades [3,5,6]. Additionally, research is being conducted to improve oil recovery and economic efficiency using a mixed solvent of CO2 and hydrocarbon gas or a surfactant [7,8,9,10]. In particular, according to previous studies that added liquefied petroleum gas (LPG) to CO2 flooding, LPG increases the displacement efficiency by reducing the minimum miscible pressure (MMP), which accelerates the effects of oil viscosity reduction, interfacial tension (IFT) reduction, and oil swelling, resulting in an increase in oil recovery [6,11,12]. Dimethyl ether (DME), an industrial chemical used as a fuel additive and aerosol propellant, is being investigated as a solvent for EOR. Recently, DME was mixed with LPG and used as fuel for transportation as well as for cooking and heating fuel [13,14]. As a solvent, DME has the advantage of being soluble in both water and oil. DME can be mixed nearly ideally in the oleic phase and is partially soluble in water because of its slight polarity [15]. These properties enable a unique and improved oil recovery mechanism. DME remaining near the injection well is extracted by the chasing water, which acts as a DME carrier, improving the oil mobility near the production well [13]. Additionally, DME can be extracted from the produced fluid and reused through distillation and adsorption processes [14]. DME was first introduced as a solvent for EOR by Shell [15,16,17,18], and this method has proven effectiveness in improving oil recovery [14,15,16,19,20,21]. Recent research on DME in the field of EOR has focused on modeling and experiments of the DME–water–oil phase behavior or studies modeling the core scale [15,16,21,22,23,24,25]. Additionally, most of the EOR using DME has been studied for the DME-enhanced water flood method, and to the best of our knowledge, no studies have been reported on a CO2–DME mixed solvent. Though DME has been proven as an effective solvent in improving oil recovery, it is an expensive solvent around $565/m3 [26]. Therefore, the optimum amount of DME into the CO2 flood should be determined through extensive compositional simulations to improve the economics of the project. Therefore, this study evaluated the use of a mixed solvent to reduce DME and increase the oil recovery rate by adding DME to the CO2 flooding. In addition, the amount of DME that was recovered by chased water was examined with compositional analysis of the produced fluid. In this study, a compositional model was developed to investigate oil recovery when injecting a CO2–DME mixed solvent under the WAG scenario. Based on the results from simulations, the sweep efficiency and displacement efficiency by WAG using a CO2–DME mixed solvent were analyzed in terms of oil viscosity, viscous gravity number, and MMP.

2. Materials and Methods

2.1. Fluid Modeling

Fluid modeling was performed using the properties of Weyburn W3 fluid. The components of the Weyburn W3 fluid and the parameters of the equation of state (EOS) for fluid modeling are listed in Table 1 [27]. The acentric factor was proposed by Pitzer [28] as a comprehensive expression of the boiling point, molecular weight and polarity characteristics of pure ingredients. Parachor is a material-specific value that relates surface tension to chemical structure.
Regression matching is promoted using the experimental data of the Weyburn W3 fluid. Results from the fluid model adequately matched the experimental data and confirmed the reliability for composition simulation. The data values compared with the modeling fluid are shown in Table 2.
The Winprop software, developed by the Computer Modeling Group (CMG), was used for fluid modeling. The Peng–Robinson EOS (PR EOS) [29,30] was used to calculate the fugacity of a component in the oil and gas phases. The concept of fugacity is usually used for multi-component equilibrium involving solid, liquid and gas phases. It is used to determine the amount of ingredients in each stage of a multi-component mixture at various temperatures and pressures [31]. Therefore, the Peng–Robinson equation was used to predict the phase behavior of the mixed fluid and to tune the experimental values. The PR EOS is as follows:
p = R T v b a v ( v + b ) + b ( v b )
or, in terms of the Z factor:
Z 3 ( 1 B ) Z 2 + ( A 3 B 2 2 B ) Z ( A B B 2 B 3 ) = 0
where a is the attraction parameter, b is the repulsion parameter, p is the pressure (Pa), T is the temperature (K), v is the volume ( m 3 mol ), and R is the universal gas constant ( m 3 · Pa K · mol ). The EOS constants for pure components are given by:
A = a p ( R T ) 2
B = b p R T
a = Ω a 0 R 2 T c 2 p c α
b = Ω b 0 R T c p c
α = [ 1 + m ( 1 T r ) ] 2
where α is a parameter for temperature, p c is the critical pressure, T c is the critical temperature, and T r is the reduced temperature. Additionally, values of the constants are given as follows:
α = [ 1 + m ( 1 T r ) ] 2   α = [ 1 + m ( 1 T r ) ] 2
The parameter m is expressed as a function of the acentric factor ω :
m = 0.37464 + 1.54226 ω 0.26992 ω 2
Robinson and Peng [30] proposed a modified expression for m for heavier components ( ω > 0.49 ):
m = 0.3796 + 1.485 ω 0.1644 ω 2 + 0.01667 ω 3
In PR EOS, the fugacity of component i of the hydrocarbon mixture, Φ i , is expressed through the following equation:
ln Φ i = B i B ( Z 1 ) ln ( Z B ) + A 2 2 B ( B i B 2 A j = 1 N y j A i j ) ln [ Z + ( 1 + 2 ) B Z ( 1 2 ) B ]
where y j is the gas mole fraction of component i. The mixing rules for expressing multi-component fugacity are as follows:
A = i = 1 N j = 1 N y i y j A i j
B = i = 1 N y i B i
A i j = ( 1 K i j ) A i A j
where K i j is the binary interaction parameter.
The parameters of DME used in PR EOS were obtained from literature (Table 3) [13,32].
DME is also soluble in the aqueous phase, and Henry’s law was used to determine the fugacity in the aqueous phase. The phase of the gas that dissolves in the water is obtained through Henry’s law as [33]:
f i w = y i w H i ,   i = 1 , . , n c
where H i is Henry’s constant of component i and y i w is the mole fraction in the aqueous phase of component i. Henry’s constant for the pressure change is calculated by:
ln H i = ln H i * + v i ¯ ( p p * ) R T
where H i * is Henry’s constant of component i at a reference pressure (Pa) and temperature (K) p * ,   T and v i ¯ is the molar volume of component i. The values of the solubility properties of DME in the aqueous phase were obtained from literature [34]. The phase behavior and solubility in the aqueous solution were determined through fugacity obtained by PR EOS and Henry’s law.

2.2. Reservoir Modeling

The GEM compositional simulator developed by CMG was used for numerical simulation. In this paper, the 2D homogenous cross-sectional reservoir model was designed to analyze the sweep efficiency change according to the DME mole fraction. The properties of the reservoir were obtained from the literature (Table 4) [31]. In oil-wet system, water flooding is not effective for enhancing oil recovery [35]. Therefore, the oil-wet system was chosen to show the improved performance of DME-CO2 mixed solvent. The end point values and crossover saturation of relative permeability are measures of wettability [36]. The relative permeability obtained from the literature is shown in Figure 1 [37]. The relative permeability of water is quite high at residual oil saturation and crossover saturation is less than 0.5, which are good qualitative indications of oil wet conditions.
The size of the reservoir was designed as 52 × 1 × 20 grids and each grid block was 3 m × 3 m × 1.5 m. In this paper, WAG was applied to enhance the oil recovery. The total injection process was conducted for 12 years: 3 years for the secondary water flooding, 6 years for WAG (1:1.5 WAG ratio), and 3 years for the water flooding to produce the residual DME (Figure 2).
Four WAG cases depending on the mole fraction of DME (0% to 25%) were analyzed (Table 4).
Other properties, including production pressure, were maintained to examine the effect of the injection fluid. Reservoir properties and the operation conditions are shown in Table 5.

2.3. Viscous Gravity Number (Ngv)

Since the density of the gas used in the WAG was less than that of the reservoir fluid, the difference between the density and gravity causes an override in which the injected gas moves to the upper layer, which reduces the vertical sweep efficiency. The vertical sweep efficiency was affected by the relationship between the viscosity and gravity. The change in the vertical sweep efficiency with the addition of DME was quantitatively investigated using the viscous gravitational number (Ngv), which can be expressed by [38]:
N g v = time   for   horizontal   flow time   for   vertical   flow = k v k r s Δ ρ g cos α A L q s μ s h
where k v is vertical permeability, k r s is relative permeability, Δ ρ is density difference between injected and reservoir fluids, g is gravitational acceleration, A is cross sectional area, L is length of the reservoir, q s is flow rate of the solvent, μ s is solvent viscosity and h is height of the reservoir.
If Ngv = 10, horizontal flow takes ten times longer than vertical flow. The guidelines for determining which gravity or viscous force dominates are as follows:
N g v < 0.1 :   Viscous   forces   dominate ; 0.1   <   N g v   < 10 :   Intermediate ; N g v > 10 :   Gravity   dominate .

2.4. Minimum Miscible Pressure

The MMP is the minimum pressure at which the injected gas and oil can be mixed under a specific reservoir temperature condition. If the reservoir pressure is less than the MMP, the immiscible flooding process is applied. Conversely, if the reservoir pressure is greater than the MMP, the miscible flooding process is applied. Though the main mechanism of the effect of oil expansion and oil viscosity reduction in immiscible flooding and miscible flooding is similar, the effect is greater in miscible flooding [39]. In addition, in the case of miscible flooding, the displacement efficiency increases because of the decrease in IFT between gas and oil [40]. In this study, the MMP between oil and gas according to DME addition was investigated. The slim tube test is recognized to be the most accurate experimental method for determining MMP. However, the slim tube test is expensive and takes a long time to measure MMP [41]. To overcome the limitations, slim tube test simulation has been widely used. Since slim tube test simulation is a cost-effective method providing an accurate prediction [42]. Therefore, we used slim tube simulation in this study.

3. Results

3.1. Effects of DME Injection on Oil Recovery

3.1.1. Vertical Sweep Efficiency

The vertical sweep efficiency according to the mole fraction of DME was compared using 2D reservoir modeling. The front of the sweep area was confirmed through the change in the viscosity of the oil in contact with the injected gas. The steeper the slope of the front contact in the sweep area, the more the injected gas contacts the fluid under the reservoir. Figure 3 shows the oil viscosity when 0.5 PV was injected. When comparing the slope of the front of sweep area, Cases 1, 2, 3 and 4 were 5, 3.75, 5, and 8, respectively. The sweep efficiency increased with the addition of DME, as confirmed by the slope comparison. The contact area at breakthrough and the viscous gravity values were compared to confirm the increase in sweep efficiency. When comparing the area of the swept grids when the breakthrough occurred, Cases 1, 2, 3, and 4 were 6372 m3, 6867.6 m3, 7235.8 m3 and 8071.2 m3, respectively. Case 4 was 26.6% greater than Case 1. Table 6 shows the viscous gravity number of all cases. In Case 1, the viscous gravity number decreased by as much as 66.4% when DME was added. Decreasing the viscous gravity number increases the viscous forces. The addition of DME is more advantageous to sweep the lower reservoir fluid. As a result, the vertical sweep efficiency was significantly improved by the addition of DME.

3.1.2. Displacement Efficiency

Oil saturation, oil viscosity, and MMP were compared to confirm the displacement efficiency. Figure 4 shows the oil saturation at the end of the WAG cycle. After three WAG cycles, the oil saturation in Case 1 remains relatively greater than Case 4. As the content of DME increases, the MMP between the injection gas and oil decreases, which has a greater effect on reducing the oil viscosity and IFT, both of which are advantageous for oil flow. Figure 5 shows the oil viscosity after 0.6 PV injection. When the initial oil viscosity was 2 mPa s, it decreased to 0.56 mPa s in Case 1 and to 0.36 mPa s in Case 4. For Case 4, the lowest oil viscosity value was 35% less than that for Case 1. In addition, in Case 4, the deep blue area, in which oil viscosity was less than 0.5 mPa s, was more widespread. The deep blue area in Case 4 was 2902.8 m3, and in Case 1, the deep blue area did not exist. The values at the middle of the reservoir were compared to evaluate the fluid properties and residual oil component after displacement. As shown in Figure 6, after some point, the oil viscosity of Case 4 was greater than that of Case 1 because the intermediate component was displaced (Figure 7a). The mole fraction of the heavy component (C28+) increases (Figure 7b) because the heavy component remains behind the flood front. MMP was measured using a slim tube (Table 7). Figure 8 shows the IFT of the swept area in the middle of the reservoir for the cases, where the IFT value of Case 4 is less than that of Case 1. The lowest IFT value of Case 4 was 0.012 mM/m. The values for Cases 1, 2, and 3 were 0.054 mM/m, 0.115 mM/m, 0.485 mM/m, which were 350%, 858%, and 3941% greater than Case 4.

3.2. Recovery

The oil recovery factor (RF) can be expressed as follows:
RF = E D E v E A
where, ED is displacement efficiency, Ev is vertical sweep efficiency and EA is areal sweep efficiency. Since is not considered in 2D cross-sectional system, oil recovery can be increased by improving ED and Ev. In Section 3.1.1 and Section 3.1.2, it was confirmed that both ED and Ev increased with the addition of DME. Figure 9 shows the oil recovery for each case according to the DME content. For Case 1, the maximum oil recovery was 55%. Cases 2, 3, and 4 showed recoveries of 60%, 64%, and 72%, respectively. As the mole fraction of DME increased, the oil recovery increased. In Case 4, the recovery increased by 31% compared with Case 1 because the displacement efficiency and sweep efficiency increased as the mole fraction of DME increased. In addition, as DME is easily dissolved in water, it is recovered with the production fluid. To improve the economics of the process, the produced DME can be recovered through distillation and adsorption processes. Figure 10 compares the cumulative amount of injected DME and produced DME in Case 4. It can be confirmed that about 70% of the injected amount has been recovered and can be reused.

4. Conclusions

A model using a mixed solvent of CO2–DME was developed to explain the change in oil recovery when CO2 WAG and DME were added. To evaluate the reliability of the compositional simulation, results from fluid modeling were matched with experimental data from real Weyburn W3 fluids. A cross-sectional 2D model was designed to examine the sweep efficiency considering the gravity override phenomenon during WAG. Since DME is an expensive solvent, a mixed solvent of DME and CO2 was used and the effect was analyzed. Four cases with different DME addition ratios were modeled and the oil recovery was analyzed through changes in the sweep efficiency and displacement efficiency. DME in CO2 decreases the IFT and MMP, which increases the displacement efficiency. Adding DME to the CO2 flood also increases the vertical sweep efficiency. The developed compositional model illustrated that DME increases the sweep efficiency and displacement efficiency. In addition, the vertical sweep efficiency and displacement efficiency was presented quantitatively with viscous gravitational numbers, MMP, and IFT. When the DME content was 25%, the oil recovery was 72%, a 31% increase over CO2 only. This study confirmed the favorable effect of DME addition into the CO2 WAG process over simple CO2 flooding. Since DME is an expensive solvent, the optimum design with an integrative techno-economical model is required for successful field implementation of the process. In addition, a full 3D model is required to analyze the effect of vertical and areal reservoir heterogeneities on the oil recovery from actual fields.

Author Contributions

Conceptualization, Y.W.L. and K.S.L.; methodology, Y.W.L., J.C. and K.S.L.; software, Y.W.L. and J.C.; validation, Y.W.L. and K.S.L.; formal analysis, Y.W.L.; investigation, Y.W.L., M.S.J. and H.S.L.; resources, Y.W.L. and K.S.L.; data curation, Y.W.L.; writing—original draft preparation, Y.W.L.; writing—review and editing, Y.W.L. and K.S.L.; visualization, Y.W.L., M.S.J. and H.S.L.; supervision, K.S.L.; project administration, K.S.L.; funding acquisition, K.S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (No. 2020R1F1A1070406). Cho was supported by the NRF grant (No. 2018R1A6A1A08025520 and No. 2020R1I1A1A01067015).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors are grateful to the Computer Modelling Group Ltd. (CMG) for technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Muggeridge, A.H.; Cockin, A.; Webb, K.; Frampton, H.; Collins, I.; Moulds, T.; Salino, P. Recovery rates, enhanced oil recovery and technological limits. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2014, 372, 20120320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Lake, L.W.; Johns, R.; Rossen, W.R.; Pope, G.A. Fundamentals of Enhanced Oil Recovery; Society of Petroleum Engineers: Richardson, TX, USA, 2014. [Google Scholar]
  3. Hu, W.; Wang, Z.; Ding, J.; Wang, Z.; Ma, Q.; Gao, Y. A new integrative evaluation method for candidate reservoirs of hydrocarbon gas drive. Geosyst. Eng. 2015, 18, 38–44. [Google Scholar] [CrossRef]
  4. Makimura, D.; Kunieda, M.; Liang, Y.; Matsuoka, T.; Takahashi, S.; Okabe, H. Application of Molecular Simulations to CO2-Enhanced Oil Recovery: Phase Equilibria and Interfacial Phenomena. SPE J. 2013, 18, 319–330. [Google Scholar] [CrossRef]
  5. Fassihi, M.R.; Gillham, T.H. The use of air injection to improve the double displacement processes. In Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, TX, USA, 3–6 October 1993. [Google Scholar]
  6. Talbi, K.; Kaiser, T.; Maini, B. Experimental Investigation of CO-Based VAPEX for Recovery of Heavy Oils and Bitumen. J. Can. Pet. Technol. 2008, 47, 1–8. [Google Scholar] [CrossRef]
  7. Cho, J.; Jeong, M.S.; Lee, Y.W.; Lee, H.S.; Lee, K.S. Techno-economic analysis of intermediate hydrocarbon injection on coupled CO2 storage and enhanced oil recovery. Energy Explor. Exploit. 2020, 1–19. [Google Scholar] [CrossRef]
  8. Cho, J.; Lee, K.S. Effects of hydrocarbon solvents on simultaneous improvement in displacement and sweep efficiencies during CO2 -enhanced oil recovery. Pet. Sci. Technol. 2016, 34, 359–365. [Google Scholar] [CrossRef]
  9. Cho, J.; Park, G.; Kwon, S.; Lee, K.S.; Lee, H.S.; Min, B. Compositional Modeling to Analyze the Effect of CH4 on Coupled Carbon Storage and Enhanced Oil Recovery Process. Appl. Sci. 2020, 10, 4272. [Google Scholar] [CrossRef]
  10. Kopcak, U.; Mohamed, R.S. Caffeine solubility in supercritical carbon dioxide/co-solvent mixtures. J. Supercrit. Fluids 2005, 34, 209–214. [Google Scholar] [CrossRef]
  11. Alston, R.; Kokolis, G.; James, C. CO2 Minimum Miscibility Pressure: A Correlation for Impure CO2 Streams and Live Oil Systems. Soc. Pet. Eng. J. 1985, 25, 268–274. [Google Scholar] [CrossRef]
  12. Kumar, N.; Von Gonten, W.D. An investigation of oil recovery by injecting CO2 and LPG mixtures. In Proceedings of the Fall Meeting of the Society of Petroleum Engineers of AIME, Las Vegas, NV, USA, 30 September–3 October 1973. [Google Scholar]
  13. Cho, J.; Kim, T.H.; Lee, K.S. Compositional modeling and simulation of dimethyl ether (DME)-enhanced water flood to investigate oil mobility improvement. Pet. Sci. 2018, 15, 297–304. [Google Scholar] [CrossRef] [Green Version]
  14. Riele, P.T.; Parsons, C.; Boerrigter, P.; Plantenberg, J.; Suijkerbuijk, B.; Burggraaf, J.; Chernetsky, A.; Boersma, D.; Broos, R. Implementing a Water Soluble Solvent Based Enhanced Oil Recovery Technology—Aspects of Field Development Planning. In Proceedings of the SPE EOR Conference at Oil and Gas West Asia, Muscat, Oman, 21–23 March 2016. [Google Scholar]
  15. Chahardowli, M.; Farajzadeh, R.; Bruining, H. Experimental investigation of the use of the dimethyl ether/polymer hybrid as a novel enhanced oil recovery method. J. Ind. Eng. Chem. 2016, 38, 50–60. [Google Scholar] [CrossRef]
  16. Chahardowli, M.; Farajzadeh, R.; Masalmeh, S.K.; Mahani, H.; Bruining, H. A Novel Enhanced Oil Recovery Technology Using Dimethyl Ether/Brine: Spontaneous Imbibition in Sandstone and Carbonate Rocks. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dubai, UAE, 26–28 September 2016. [Google Scholar]
  17. Chernetsky, A.; Masalmeh, S.; Eikmans, D.; Boerrigter, P.M.; Fadili, A.; Parsons, C.A.; Parker, A.; Boersma, D.M.; Cui, J.; Dindoruk, B. A novel enhanced oil recovery technique: Experimental results and modelling workflow of the DME enhanced water flood technology. In Proceedings of the Abu Dhabi International Petroleum Exhibition and Conference, Abu Dhabi, UAE, 9–12 November 2015. [Google Scholar]
  18. Parsons, C.; Chernetsky, A.; Eikmans, D.; Te Riele, P.; Boersma, D.; Sersic, I.; Broos, R. Introducing a novel enhanced oil recovery technology. In Proceedings of the IOR 2017-19th European Symposium on Improved Oil Recovery, Tulsa, OK, USA, 11–13 April 2016. [Google Scholar]
  19. Chahardowli, M.; Zholdybayeva, A.; Farajzadeh, R.; Bruining, H. Solvent-enhanced spontaneous imbibition in fractured reservoirs. In Proceedings of the EAGE Annual Conference & Exhibition incorporating SPE Europec, London, UK, 10–13 June 2013. [Google Scholar]
  20. Javanmard, H.; Seyyedi, M.; Nielsen, S.M. On Oil Recovery Mechanisms and Potential of DME–Brine Injection in the North Sea Chalk Oil Reservoirs. Ind. Eng. Chem. Res. 2018, 57, 15898–15908. [Google Scholar] [CrossRef]
  21. Groot, J.A.W.M.; Chernetsky, A.; Riele, P.M.T.; Dindoruk, B.; Cui, J.; Wilson, L.C.; Ratnakar, R.R. Representation of Phase Behavior and PVT Workflow for DME Enhanced Water-Flooding. In Proceedings of the SPE EOR Conference at Oil and Gas West Asia, Muscat, Oman, 21–23 March 2016. [Google Scholar]
  22. Ratnakar, R.R.; Dindoruk, B.; Wilson, L. Use of DME as an EOR Agent: Experimental and Modeling Study to Capture Interactions of DME, Brine and Crudes at Reservoir Conditions. In Proceedings of the SPE Annual Technical Conference and Exhibition, Dubai, UAE, 26–28 September 2016. [Google Scholar]
  23. Ratnakar, R.R.; Dindoruk, B.; Wilson, L. Experimental investigation of DME–water–crude oil phase behavior and PVT modeling for the application of DME-enhanced water flooding. Fuel 2016, 182, 188–197. [Google Scholar] [CrossRef]
  24. Ratnakar, R.R.; Dindoruk, B.; Wilson, L.C. Development of empirical correlation for DME-partitioning between brine and crudes for enhanced water flooding applications. J. Pet. Sci. Eng. 2017, 157, 264–272. [Google Scholar] [CrossRef]
  25. Ratnakar, R.R.; Dindoruk, B.; Wilson, L.C. Phase behavior experiments and PVT modeling of DME-brine-crude oil mixtures based on Huron-Vidal mixing rules for EOR applications. Fluid Phase Equilibria 2017, 434, 49–62. [Google Scholar] [CrossRef]
  26. Karagoz, S. Process Design, Simulation and Integration of Dimethyl Ether (DME) Production from Shale Gas by Direct and Indirect Methods. Master’s Thesis, Texas A & M University, College Station, TX, USA, 2014. [Google Scholar]
  27. Srivastava, R.; Huang, S.; Dong, M. Laboratory Investigation of Weyburn CO2 Miscible Flooding. J. Can. Pet. Technol. 2000, 39, 41–51. [Google Scholar] [CrossRef]
  28. Pitzer, K.S.; Lippmann, D.Z.; Curl, R.F.; Huggins, C.M.; Petersen, D.E. The Volumetric and Thermodynamic Properties of Fluids. II. Compressibility Factor, Vapor Pressure and Entropy of Vaporization1. J. Am. Chem. Soc. 1955, 77, 3433–3440. [Google Scholar] [CrossRef]
  29. Peng, D.-Y.; Robinson, D.B. A New Two-Constant Equation of State. Ind. Eng. Chem. Fundam. 1976, 15, 59–64. [Google Scholar] [CrossRef]
  30. Robinson, D.B.; Peng, D.Y. The Characterization of the Heptanes and Heavier Fractions for the GPA Peng-Robinson Programs; Research Report 28; Gas Processors Association: Tulsa, OK, USA, 1978. [Google Scholar]
  31. Lee, K.S.; Cho, J.; Lee, J.H. CO2 Storage Coupled with Enhanced Oil Recovery; Springer Science and Business Media LLC: Cham, Switzerland, 2020. [Google Scholar]
  32. Ganjdanesh, R.; Rezaveisi, M.; Pope, G.A.; Sepehrnoori, K. Treatment of Condensate and Water Blocks in Hydraulic-Fractured Shale-Gas/Condensate Reservoirs. SPE J. 2016, 21, 665–674. [Google Scholar] [CrossRef]
  33. Li, Y.K.; Nghiem, L.X. Phase equilibria of oil, gas and water/brine mixtures from a cubic equation of state and Henry’s law. Can. J. Chem. Eng. 1986, 64, 486–496. [Google Scholar] [CrossRef]
  34. Sander, R. Compilation of Henry’s law constants (version 4.0) for water as solvent. Atmos. Chem. Phys. Discuss. 2015, 15, 4399–4981. [Google Scholar] [CrossRef] [Green Version]
  35. Kułynycz, V. The influence of wettability on oil recovery. Drilling–Oil–Gas AGH 2015, 32, 493–502. [Google Scholar]
  36. Wheaton, R. Basic Rock and Fluid Properties. In Fundamentals of Applied Reservoir Engineering; Elsevier BV: Amsterdam, The Netherlands, 2016; pp. 5–57. [Google Scholar]
  37. Delshad, M.; Najafabadi, N.F.; Anderson, G.A.; Pope, G.A.; Sepehrnoori, K. Modeling wettability alteration in naturally fractured reservoirs. In Proceedings of the SPE/DOE Symposium on Improved Oil Recovery, Tulsa, OK, USA, 22–26 April 2006. [Google Scholar]
  38. Chugh, S.; Baker, R.; Cooper, L.; Spence, S. Simulation of Horizontal Wells to Mitigate Miscible Solvent Gravity Override in the Virginia Hills Margin. J. Can. Pet. Technol. 2000, 39, 28–34. [Google Scholar] [CrossRef]
  39. Cho, J.; Park, S.S.; Jeong, M.S.; Lee, K.S. Compositional Modeling for Optimum Design of Water-Alternating CO2 -LPG EOR under Complicated Wettability Conditions. J. Chem. 2015, 2015, 1–9. [Google Scholar] [CrossRef] [Green Version]
  40. Pourhadi, S.; Fath, A.H. Performance of the injection of different gases for enhanced oil recovery in a compositionally grading oil reservoir. J. Pet. Explor. Prod. Technol. 2019, 10, 641–661. [Google Scholar] [CrossRef] [Green Version]
  41. Adel, I.A.; Tovar, F.D.; Schechter, D.S. Fast-Slim Tube: A Reliable and Rapid Technique for the Laboratory Determination of MMP in CO2—Light Crude Oil Systems. In Proceedings of the SPE Improved Oil Recovery Conference, Tulsa, OK, USA, 11–13 April 2016. [Google Scholar]
  42. Ameri, M.; Fard, M.A.; Akbari, M.; Zamanzadeh, S.; Nasiri, E. A Comparison of Accuracy and Computational Time for Common and Artificial Methods in Predicting Minimum Miscibility Pressure. Energy Explor. Exploit. 2013, 31, 221–236. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Relative permeability curve for the (a) water–oil and (b) liquid–gas systems.
Figure 1. Relative permeability curve for the (a) water–oil and (b) liquid–gas systems.
Applsci 11 00406 g001
Figure 2. Schematic diagram of the water alternating gas (WAG) injection design.
Figure 2. Schematic diagram of the water alternating gas (WAG) injection design.
Applsci 11 00406 g002
Figure 3. Oil viscosity when 0.5 PV solvent was injected: (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Figure 3. Oil viscosity when 0.5 PV solvent was injected: (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Applsci 11 00406 g003
Figure 4. Residual oil saturation after the last WAG cycle for (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Figure 4. Residual oil saturation after the last WAG cycle for (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Applsci 11 00406 g004
Figure 5. Oil viscosity when 0.6 PV solvent was injected for (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Figure 5. Oil viscosity when 0.6 PV solvent was injected for (a) Case 1, (b) Case 2, (c) Case 3, and (d) Case 4.
Applsci 11 00406 g005
Figure 6. Comparison of the oil viscosity for Cases 1 and 4 at the midpoint of the reservoir.
Figure 6. Comparison of the oil viscosity for Cases 1 and 4 at the midpoint of the reservoir.
Applsci 11 00406 g006
Figure 7. Residual oil content after sweep for Cases 1 and 4: (a) C6–9 mole fraction, (b) C28+ mole fraction at the midpoint of the reservoir.
Figure 7. Residual oil content after sweep for Cases 1 and 4: (a) C6–9 mole fraction, (b) C28+ mole fraction at the midpoint of the reservoir.
Applsci 11 00406 g007
Figure 8. Comparison of interfacial tension (IFT) for all cases at the midpoint of the reservoir.
Figure 8. Comparison of interfacial tension (IFT) for all cases at the midpoint of the reservoir.
Applsci 11 00406 g008
Figure 9. Oil recovery according to DME and the CO2 ratio.
Figure 9. Oil recovery according to DME and the CO2 ratio.
Applsci 11 00406 g009
Figure 10. DME injection and DME production in Case 4.
Figure 10. DME injection and DME production in Case 4.
Applsci 11 00406 g010
Table 1. Oil components and properties of each component for equation of state (EOS) calculation.
Table 1. Oil components and properties of each component for equation of state (EOS) calculation.
ComponentCompositionCritical Pressure
(kPa)
Critical Temperature
(K)
ParachorAcentric FactorMolecular Weight
(g/gmol)
N20.02073394.4126.2410.0428.01
CO20.00747376.5304.2780.22544.01
H2S0.00128936.9373.2800.134.08
C10.07494600.2190.6770.00816.04
C20.04224883.9305.41080.09830.07
C30.07854245.5369.81500.15244.09
DME0.00015442.2400.31320.246.07
C40.06553722.7416.51850.184558.12
C50.04593379.4464.92280.23972.15
C6-90.21553019.6556.42960.331102.5
C10-170.22022017.5692.35050.584184
C18-270.10271327.0808.47680.893306.2
C28+0.12521155.1915.510011.1585.61
Table 2. Comparison of the fluid model and experimental data for Weyburn W3 fluid.
Table 2. Comparison of the fluid model and experimental data for Weyburn W3 fluid.
ParametersFluid ModelExperimental Data
Saturation Pressure (kPa)4784.954920
Oil density at psat (kg/m3)805.73806.4
Oil viscosity at psat (mPa·s)1.751.76
Formation volume factor1.111.12
API (°)34.831
Table 3. Binary interaction coefficients between DME and each component.
Table 3. Binary interaction coefficients between DME and each component.
N2CO2H2SC1C2C3DMEC4C5C6-9C10–17C18–27C28+
N20.00------------
CO20.000.00-----------
H2S0.130.140.00----------
C10.250.110.070.00---------
C20.010.130.090.000.00--------
C30.090.130.080.010.000.00-------
DME0.100.000.000.290.250.250.00------
C40.100.120.080.020.010.000.250.00-----
C50.100.120.070.020.010.000.250.000.00----
C6-90.110.120.050.030.020.010.200.000.000.00---
C10-170.110.120.050.060.040.020.080.020.010.010.00--
C18-270.110.120.050.090.060.050.080.030.030.020.000.00-
C28+0.110.120.050.120.090.070.080.050.050.030.010.000.00
Table 4. Classification of cases according to DME content.
Table 4. Classification of cases according to DME content.
CaseInjection Gas Contents
1100% CO2 + 0% DME
290% CO2 + 10% DME
385% CO2 + 15% DME
475% CO2 + 25% DME
Table 5. Reservoir initial conditions and operating conditions.
Table 5. Reservoir initial conditions and operating conditions.
PropertiesValues
Initial Pressure (kPa)13,789
Initial oil saturation0.88
Porosity (%)30
Permeability (m2)4.9 × 10−14
Producing pressure (kPa)11,031
Total injection (PV)1.5
Table 6. Comparison of the viscous gravity number.
Table 6. Comparison of the viscous gravity number.
CaseViscous Gravity Number (Ngv)
15.06
23.94
33.26
41.70
Table 7. Minimum miscible pressure (MMP) for each case calculated using a slim tube.
Table 7. Minimum miscible pressure (MMP) for each case calculated using a slim tube.
CaseMMP (kPa)
114,127
213,065
312,079
49866
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lee, Y.W.; Lee, H.S.; Jeong, M.S.; Cho, J.; Lee, K.S. Compositional Modeling of Dimethyl Ether–CO2 Mixed Solvent for Enhanced Oil Recovery. Appl. Sci. 2021, 11, 406. https://doi.org/10.3390/app11010406

AMA Style

Lee YW, Lee HS, Jeong MS, Cho J, Lee KS. Compositional Modeling of Dimethyl Ether–CO2 Mixed Solvent for Enhanced Oil Recovery. Applied Sciences. 2021; 11(1):406. https://doi.org/10.3390/app11010406

Chicago/Turabian Style

Lee, Young Woo, Hye Seung Lee, Moon Sik Jeong, Jinhyung Cho, and Kun Sang Lee. 2021. "Compositional Modeling of Dimethyl Ether–CO2 Mixed Solvent for Enhanced Oil Recovery" Applied Sciences 11, no. 1: 406. https://doi.org/10.3390/app11010406

APA Style

Lee, Y. W., Lee, H. S., Jeong, M. S., Cho, J., & Lee, K. S. (2021). Compositional Modeling of Dimethyl Ether–CO2 Mixed Solvent for Enhanced Oil Recovery. Applied Sciences, 11(1), 406. https://doi.org/10.3390/app11010406

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop