Research Article  Open Access
Compositional Modeling for Optimum Design of WaterAlternating CO_{2}LPG EOR under Complicated Wettability Conditions
Abstract
The addition of LPG to the CO_{2} stream leads to minimum miscible pressure (MMP) reduction that causes more oil swelling and interfacial tension reduction compared to CO_{2} EOR, resulting in improved oil recovery. Numerical study based on compositional simulation has been performed to examine the injectivity efficiency and transport behavior of wateralternating CO_{2}LPG EOR. Based on oil, CO_{2}, and LPG prices, optimum LPG concentration and composition were designed for different wettability conditions. Results from this study indicate how injected LPG mole fraction and butane content in LPG affect lowering of interfacial tension. Interfacial tension reduction by supplement of LPG components leads to miscible condition causing more enhanced oil recovery. The maximum enhancement of oil recovery for oilwet reservoir is 50% which is greater than 22% for waterwet reservoir. According to the result of net present value (NPV) analysis at designated oil, CO_{2}, propane, and butane prices, the optimal injected LPG mole fraction and composition exist for maximum NPV. At the case of maximum NPV for oilwet reservoir, the LPG fraction is about 25% in which compositions of propane and butane are 37% and 63%, respectively. For waterwet reservoir, the LPG fraction is 20% and compositions of propane and butane are 0% and 100%.
1. Introduction
CO_{2} injection has been found to be an efficient method for oil recovery worldwide through a miscible or an immiscible displacement process. Mechanism of CO_{2} enhanced oil recovery (EOR) is divided into two different processes, miscible flood and immiscible flood. Although miscible gas injection is a widely applied EOR process, it can be only applied when the reservoir pressure is higher than minimum miscible pressure (MMP). The main process of miscible gas injection is displacement efficiency improvement by oil viscosity reduction and swelling effect to reduce residual oil saturation. When reservoir pressure is higher than MMP, the injected CO_{2} and reservoir oil are completely miscible and the displacement efficiency can be enhanced by zero interfacial tension [1]. Immiscible flood is usually applied when reservoir pressure is insufficient to miscible flood or reservoir oil contains many heavy components. The effects of immiscible flood are similar to miscible flood, but one major disadvantage is the limited solubility of CO_{2} in oil, resulting in the restricted swelling effect and viscosity reduction.
Injected CO_{2} and reservoir oil can be miscible by continuous contact. At the foreend of injected fluid, CO_{2} is persistently contacted with fresh oil following flow direction, and they are eventually miscible by the vaporizinggas drive process. In contrast, at the backend of injected CO_{2}, near injection well, reservoir oil is continuously contacted with fresh CO_{2} that causes the miscible state by the condensinggas drive process [2]. CO_{2} miscible flood making high enhanced oil recovery effect has a limit that it can be only applied when the reservoir pressure is higher than MMP. It can be settled by the application of CO_{2}LPG EOR that is able to lower MMP less than that from the application of only CO_{2} EOR. The addition of alkane solvents to the CO_{2} injection generally accelerates swelling oil, reducing oil viscosity and decreasing the interfacial tension that can lead to better performance in enhancing oil recovery [3]. The effects of CO_{2}LPG injection are verified by the experiment [4].
Figure 1 indicates the ternary diagram of phase behavior of reservoir oil and injected solvent. J and I signify injected fluid and reservoir oil. In the inner area of the ternary diagram curve, two phases of the reservoir fluid exist. In case of J_{1}I_{2}, only CO_{2} is injected into reservoir oil I_{2}. J_{1} and I_{2} cannot be miscible at the first contact because the line passes through the twophase area. However, they arrive at miscible condition by multiple contact miscibility process. The J_{1}I_{1} line lies on the twophase territory and both points J_{1} and I_{1} are located in the same side on the basis of limiting tie line. Therefore, J_{1}I_{1} cannot be miscible by first and multiple contact miscibility process. At J_{2}I_{1} and J_{3}I_{1} cases, first and miscible contact miscibility process is available. By the addition of LPG to the CO_{2} stream, the location of injected solvent is moved from J_{1} to J_{2} or J_{3} depending on the amount of injected LPG. It makes miscible condition from J_{1}I_{1} case that was not supposed to be miscible.
To improve sweep efficiency, WAG (wateralternatinggas) process is applied to CO_{2}LPG EOR method in this research. At the same WAG condition, injected LPG amount and composition are the variables considered in the study. Many experimental researches about the effects of LPG and impurities on MMP with oil have been actively developed [6]. Several established researches demonstrate the MMP reduction and oil recovery enhancement by CO_{2}LPG EOR through only experimental ways [7, 8], but numerical approach to analyze the effectiveness of CO_{2}LPG EOR was not included. Shokir [9] developed ACE algorithm model to analyze the effects of impurities on MMP between injected fluid and oil, but it could not explain how lower MMP affects oil recovery. Talbi et al. [3] conducted experimental research on oil swelling, viscosity reduction, interfacial tension reduction, and oil recovery improvement that resulted from injecting solvents into CO_{2}. However, if it is applied to field scale, it should be time consuming, so reservoir simulation model for CO_{2}LPG flood is positively necessary. Recently, Teklu et al. [10] showed MMP reduction effect by CO_{2}LPG flood in various reservoir scenarios using simulation model, but it focused on only the relationship between pore confinement, permeability, and MMP in shale reservoirs. It also does not make connection between MMP reduction and oil recovery in consideration of gas transport in porous media. Recent studies based on the modeling of spontaneous imbibition also indicated that transport properties of oil are affected by wettability condition [11, 12]. For this reason, different wettability conditions are applied for analyzing the performance of CO_{2}LPG EOR.
It has been identified that CO_{2}LPG flood is an effective method for MMP reduction causing oil recovery enhancement through many experimental studies. Compositional model for CO_{2}LPG EOR is necessary to investigate how gas transport affects MMP reduction and oil recovery enhancement. In this research, compositional fluid and multiphase simulation models are developed and injected LPG mole fraction and composition are optimized based on recent oil, CO_{2}, propane, and butane prices for maximum net present value (NPV).
2. Numerical Simulation
2.1. Fluid Modeling
Fluid data of Weyburn reservoir is referred for NPV based solvent injection simulation. Weyburn reservoir, located in southeast Saskatchewan and operated by PanCanadian Petroleum Ltd., has reached its economic limit of production by waterflooding. The reservoir is a target for CO_{2} miscible flooding to enhance oil recovery. The oil composition is shown in Table 1 and comparison between computed fluid model properties and actual fluid data of Weyburn reservoir is given in Table 2 [13]. Oil gravity, formation volume factor, and gasoil ratio are calculated through regression process to match separator experimental data. Saturation pressure is also computed by regression process. Details of the calculation techniques for saturation pressure can be found in [14]. Acceptable match of computed properties from fluid model and Weyburn’s data increases reliability of the fluid model for compositional simulation.


Phase behavior of fluid model was determined by PengRobinson EOS [15] with the reservoir oil and injected fluid composition. The PR EOS is given byor in terms of Z factor,and .
The EOS constants for pure components are given bywhere , , and
Robinson and Peng [16] proposed a modified for heavier components () as follows:
Fugacity expressions are given bywhere mixing rules are used for multicomponent fugacity expression as follows:where is binaryinteraction parameters.
Multiple mixing cell method [17] was applied to fluid model to estimate MMP between injected CO_{2} and reservoir oil. Multiple mixing cell method follows the order below.(1)Specify the reservoir temperature and an initial pressure.(2)Calculate the tieline length for each pressure step by using the equation below: where is the number of components and and are liquid and gas equilibrium compositions, respectively.(3)Draw a tieline length graph as a function of pressures.(4)Perform a multipleparameter regression of the minimum tieline lengths to determine the exponent in (powerlaw extrapolation). These parameters are determined when correlation coefficient exceeds 0.999.(5)Determine the MMP when the powerlaw extrapolation gives zero of minimum tieline length.
After generating the fluid model which has approximate MMP to Weyburn fluid, MMPs were computed between oil and LPGs. The composition of LPG is propane 63% and butane 37%, and the calculated MMPs are indicated in Table 3. MMPs of LPG (composition: propane 100% and butane 0%) mole fraction 20% and 25% are 1,747 psi and 1,614 psi.

2.2. Interfacial Tension Calculations
The equation for calculating interfacial tension in multicomponent systems is as follows [18]:where is the interfacial tension between liquid and gas phases () and and are molar densities of liquid and gas phases (), respectively. The parachor () is defined as follows:where and CN is the carbon number of the components .
2.3. Reservoir Modeling
The reservoir model was assumed as 2D model which is discretized into 33 × 33 × 1 grid blocks. Each grid block has dimension as 10 ft × 10 ft × 20 ft as shown in Figure 2. The model size is general one injectorone producer scale of 10acre fivespot model [19]. This simulation study utilized homogeneous 2D areal model not considering heterogeneity and gas overriding effect. Without these effects, oil recovery can be governed only by displacement efficiency from LPG addition and can be expected near 100% [20].
Contact angle which is a determinant for wettability is defined by Young’s equation as follows: where , , and are oilsolid, watersolid, and oilwater interfacial tensions. As indicated in the above equation, if is greater than , is smaller than 90°, so the reservoir rock exhibits waterwet solid. The inverse case is oilwet condition. Waterwet and oilwet reservoirs have the constant porosity and isotropic permeability is also assumed. Reservoir initial conditions are shown in Table 4. The porosity, permeability, and relative permeability were gained from the same reservoir, and two different relative permeability curves (Figure 3) are used in this simulation for establishing different residual oil saturation and mobility [21]. The relative permeability curves are predicted by simulations. Simulation methods to predict relative permeability are already verified by previous studies [22, 23] and similar water relative curve can be found. Residual oil saturations of oil and waterwet reservoirs are 18% and 15%, respectively.

After waterflooding for three years, wateralternating CO_{2} EOR and CO_{2}LPG EOR were applied to water and oilwet reservoirs for ten years. WAG cycle of CO_{2} and CO_{2}LPG EOR is 1 : 1, and one cycle period is 6 months. Production pressure is 1,500 psi which is within a limitation of miscible condition by first or multiple contact miscibility process when added LPG concentration is larger than 20% (Figure 4). Injected LPG mole fraction and composition are indicated in Table 5.

2.4. Net Present Value
The NPV of a time series of cash flows is defined as the sum of the present values. NPV considering prices of oil, CO_{2}, propane, and butane and costs of water injection and produced water handling is calculated by the following equation [24]:where is total production period (day), is time, is net profit at time , and is daily discount rate. is estimated by yearly discount rate aswhere yearly discount rate is 10% and is defined by the difference between the profit from oil production and total investment costs at time :where , , , , , and are oil production rate (bbl/day), CO_{2} injection rate (lb/day), propane injection rate (lb/day), butane injection rate (lb/day), water injection rate (bbl/day), and water production rate (bbl/day). Parameters , , , , , and are oil price ($/bbl), CO_{2} price ($/lb), propane price ($/lb), butane price ($/lb), water injection cost ($/bbl), and produced water handling cost ($/bbl). All values of parameters for NPV calculation are shown in Table 6 [25, 26].

3. Results and Discussion
3.1. Oil Production
The aim of this study is to confirm the effectiveness of wateralternating CO_{2}LPG EOR process in oil recovery for different reservoirs. The performance of CO_{2}LPG injection process has been compared with that of CO_{2} WAG process. LPG is composed of 63% propane and 37% butane. Results of oil recovery with various LPG concentrations are indicated in Figure 5. Increased oil recoveries for oil and waterwet reservoirs by CO_{2}LPG flood are 46% and 22%. For both wettability conditions, the higher LPG mole fraction is injected, the more oil is produced. However, significant differences are not found if LPG mole fraction is greater than 25%. The tendency is also identified by experimental results in the literature [7]. To detect the influence of LPG composition in the injected fluid, oil recoveries with different ratio of propane and butane (LPG 15%) are shown in Figure 6. Figure 6 shows that the higher fraction of butane causes more enhanced oil recovery. Increments of oil recovery for oil and waterwet reservoirs are 25% and 15% as compared with CO_{2} EOR. This phenomenon was already revealed by the experimental study and it was explained that the result is because of higher mole weight [27].
(a) Oilwet
(b) Waterwet
(a) Oilwet
(b) Waterwet
When reservoir oil and injected gas are miscible, gas saturation decreases further than immiscible condition (Figure 7). Injected gas reached production well at around 2004, so the gas saturation of WAG CO_{2} case increased abruptly. However, in case of WAG CO_{2} + LPG 30%, the gas saturation did not increase even though injected gas already reached production well. It indicates that miscible condition reduces gas saturation. The reduction of gas saturation causes a decrease in gas relative permeability (Figure 8). As both gas saturation and relative permeability decline, liquid saturation and relative permeability increase, which leads to the enhancement of oil recovery.
The addition of LPG to CO_{2} stream is more effective to lower interfacial tension between oil and gas phases. In particular, if the reservoir is in miscibility condition, interfacial tension reaches zero [1]. As shown in Figure 9(a), the swept zone is left in nonzero interfacial area, so reservoir is not in miscible condition by only CO_{2} injection. In contrast, the addition of LPG to CO_{2} stream as 20% brings the swept area into zero interfacial zone indicating miscible condition (Figure 9(b)). Zero interfacial tension indicates that oil and gas become singlephase, so it flows easier than twophase fluid.
(a) CO2 100%
(b) LPG 20%
If more butane content is injected than propane, more oil recovery is expected because of its higher molecular weight. It was proved that butane is much more effective in MMP reduction [27]. The addition of alkane solvents to the CO_{2} stream accelerates the process of reducing oil viscosity; thus, it leads to higher oil recovery [4]. To compare the aspect of oil recovery by injected LPG composition, oil saturation in reservoir is shown in Figure 10. Figure 10 indicates the oil saturation when LPG mole fraction is 25% for oilwet reservoir after 6 months from the end of waterflooding. In case of 100% propane, oil saturation near injection well is zero because of immaculate expulsion and it is 0.5 at the foreend of injected fluid (Figure 10(a)). In case of 100% butane, zero zone of oil saturation is more widespread with near wellbore as a center. Furthermore, oil saturation at the foreend is 0.7 which is higher than that in the case of 100% propane because more oil is displaced from wider area (Figure 10(b)).
(a) 100% propane
(b) 100% butane
Tables 7 and 8 indicate the amount of increase in maximum NPV after waterflooding for different wettability conditions. NPVs are calculated according to LPG concentration and composition. The maximum NPV increments by CO_{2} WAG are 12% and 13% for oil and waterwet reservoirs. As shown in Tables 7 and 8, the maximum value is 24.1% (LPG 25%: propane 63%, butane 37%) for oilwet reservoir and 17.0% (LPG 20%: propane 0%, butane 100%) for waterwet reservoir. When LPG mole fraction is less than 15% and 20% (propane 100%), maximum increase in NPV is less than CO_{2} WAG. These cases are in immiscible condition, so oil recovery is not high compared to the economic feasibility of LPG. Injected fluid and reservoir oil are in miscible condition, and maximum NPV improvements are higher than those of CO_{2} WAG cases. Maximum NPV increment by CO_{2}LPG EOR occurred for two different wettability conditions, but the effect in oilwet reservoir is better than in waterwet reservoir because of higher residual oil saturation after waterflooding.


4. Conclusions
In this study, wateralternating CO_{2}LPG EOR simulation model was developed. To examine the efficiency of CO_{2}LPG EOR considering oil, CO_{2}, and LPG prices, extensive simulations have been performed for different wettability conditions and the following conclusions have been drawn.(1)When LPG concentration is 30% and composition of butane is 100%, oil recovery increased by 46% and 25% for oilwet reservoir. When LPG concentration is 30% and butane composition is 100%, the maximum increasing amounts are 22% and 15% in case of waterwet reservoir. As injected LPG concentration and butane composition increased, significantly enhanced oil recovery was observed from the reduction of MMP and interfacial tension. Oil recovery for different wettability by CO_{2}LPG EOR has become close to 100%.(2)When LPG concentration is 25% and butane composition is 37%, maximum NPV improvement is 24.1% for oilwet reservoir. When LPG concentration is 20% and butane composition is 100%, maximum NPV improvement is 17.0% for waterwet reservoir. For both oil and waterwet reservoirs, when LPG concentrations are 10%, 15%, 20% (propane 100%), and 25% (propane 100%), the reservoir condition is immiscible and maximum NPV increment is lower than CO_{2} WAG process. When LPG concentration is higher than 20% (miscible condition), maximum NPV improved and optimum LPG concentration and composition exist for maximum NPV improvement.(3)CO_{2}LPG EOR can be applicable in low pressure reservoirs that CO_{2} is not miscible. LPG addition to CO_{2} stream can appreciably improve oil recovery by zero interfacial tension bringing miscible condition. Moreover, the optimization of LPG concentration and composition is absolutely necessary for economic feasibility. The necessity of optimization is required more in oilwet reservoir due to better performance of displacement efficiency.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by the Energy Efficiency & Resources Core Technology Program of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) granted financial resource from the Ministry of Trade, Industry & Energy, Republic of Korea (no. 20122010200060).
References
 D. Makimura, M. Kunieda, Y. Liang, T. Matsuoka, S. Takahashi, and H. Okabe, “Application of molecular simulations to CO_{2}enhanced oil recovery: phase equilibria and interfacial phenomena,” Society of Petroleum Engineers, vol. 18, no. 2, pp. 319–330, 2013. View at: Publisher Site  Google Scholar
 F. I. Stalkup Jr., Miscible Displacement, Society of Petroleum Engineers of AIME, Dallas, Tex, USA, 2nd edition, 1983.
 K. Talbi, T. M. V. Kaiser, and B. B. Maini, “Experimental investigation of CO_{2}based VAPEX for recovery of heavy oils and bitumen,” Journal of Canadian Petroleum Technology, vol. 47, no. 4, pp. 29–36, 2008. View at: Google Scholar
 H. Li, S. Zheng, and D. Yang, “Enhanced swelling effect and viscosity reduction of solvent(s)/CO_{2}/heavyoil systems,” SPE Journal, vol. 18, no. 4, pp. 695–707, 2013. View at: Publisher Site  Google Scholar
 L. W. Lake, Enhanced Oil Recovery, Society of Petroleum Engineers, Richardson, Tex, USA, 2010.
 R. B. Alston, G. P. Kokolis, and C. F. James, “CO_{2} minimum miscibility pressure: a correlation for impure CO_{2} streams and live oil systems,” Society of Petroleum Engineers Journal, vol. 25, no. 2, pp. 268–274, 1985. View at: Publisher Site  Google Scholar
 N. Kumar and W. D. V. Gonten, “An investigation of oil recovery by injecting CO_{2} and LPG mixtures,” in Proceedings of the 48th Annual Fall Meeting of the Society of Petroleum Engineers of AIME, SPE 4581, Las Vegas, Nev, USA, September 1973. View at: Publisher Site  Google Scholar
 P. Y. Zhang, S. Huang, S. Sayegh, and X. L. Zhou, “Effect of CO_{2} impurities on gasinjection EOR processes,” in Proceedings of the SPE/DOE Symposium on Improved Oil Recovery, SPE 89477, Tulsa, Okla, USA, April 2004. View at: Google Scholar
 E. M. E. Shokir, “CO_{2}oil minimum miscibility pressure model for impure and pure CO_{2} streams,” in Proceedings of the Offshore Mediterranean Conference and Exhibition (OMC '07), Ravenna, Italy, March 2007. View at: Google Scholar
 T. W. Teklu, N. Alharthy, H. Kazemi, X. Yin, and R. M. Graves, “Hydrocarbon and nonhydrocarbon gas miscibility with light oil in shale reservoirs,” in Proceedings of the SPE Improved Oil Recovery Symposium, SPE 169123, Tulsa, Okla, USA, April 2014. View at: Google Scholar
 J. Cai, X. Hu, D. C. Standnes, and L. You, “An analytical model for spontaneous imbibition in fractal porous media including gravity,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 414, pp. 228–233, 2012. View at: Publisher Site  Google Scholar
 J. Cai, E. Perfect, C.L. Cheng, and X. Hu, “Generalized modeling of spontaneous imbibition based on hagenpoiseuille flow in tortuous capillaries with variably shaped apertures,” Langmuir, vol. 30, no. 18, pp. 5142–5151, 2014. View at: Publisher Site  Google Scholar
 R. K. Srivastava, S. S. Huang, and M. Dong, “Laboratory investigation of Weyburn CO_{2} miscible flooding,” Journal of Canadian Petroleum Technology, vol. 39, no. 2, pp. 41–51, 2000. View at: Google Scholar
 L. X. Nghiem, Y.K. Li, and R. A. Heidemann, “Application of the tangent plane criterion to saturation pressure and temperature computations,” Fluid Phase Equilibria, vol. 21, no. 12, pp. 39–60, 1985. View at: Publisher Site  Google Scholar
 D.Y. Peng and D. B. Robinson, “A new twoconstant equation of state,” Industrial and Engineering Chemistry Fundamentals, vol. 15, no. 1, pp. 59–64, 1976. View at: Publisher Site  Google Scholar
 D. B. Robinson and D. Y. Peng, “The characterization of the heptanes and heavier fractions,” Research Report 28, Gas Processors Association, Tulsa, Okla, USA, 1978. View at: Google Scholar
 K. Ahmadi and R. T. Johns, “Multiplemixingcell method for MMP calculations,” SPE Journal, vol. 16, no. 4, pp. 733–742, 2011. View at: Publisher Site  Google Scholar
 R. C. Reid, J. M. Prausnitz, and T. K. Sherwood, The Properties of Gases and Liquids, McGrawHill, New York, NY, USA, 1977.
 G. P. Willhite, Waterflooding, Society of Petroleum Engineers, Richardson, Tex, USA, 1986.
 G. F. Teletzke, P. D. Patel, and A. L. Chen, “Methodology for miscible gas injection EOR screening,” in Proceedings of the SPE International Improved Oil Recovery Conference in Asia Pacific (IIORC '05), SPE 97650, pp. 315–325, Kuala Lumpur, Malaysia, December 2005. View at: Google Scholar
 M. Delshad, N. F. Najafabadi, G. A. Anderson, G. A. Pope, and K. Sepehrnoori, “Modeling wettability alteration in naturally fractured reservoirs,” in Proceedings of the 15th SPE/DOE Improved Oil Recovery Symposium, vol. 2 of SPE 100081, Tulsa, Okla, USA, April 2006. View at: Google Scholar
 B. Xiao, J. Fan, and F. Ding, “Prediction of relative permeability of unsaturated porous media based on fractal theory and Monte Carlo simulation,” Energy and Fuels, vol. 26, no. 11, pp. 6971–6978, 2012. View at: Publisher Site  Google Scholar
 B. Xiao, J. Fan, and F. Ding, “A fractal analytical model for the permeabilities of fibrous gas diffusion layer in proton exchange membrane fuel cells,” Electrochimica Acta, vol. 134, pp. 222–231, 2014. View at: Publisher Site  Google Scholar
 S. Salem and T. Moawad, “Economic study of miscible CO_{2} flooding in a mature waterflooded oil reservoir,” in Proceedings of the SPE Saudi Arabia Section Annual Technical Symposium and Exhibition, SPE 168064, AlKhobar, Saudi Arabia, May 2013. View at: Publisher Site  Google Scholar
 C. L. Liao, X. W. Liao, X. L. Zhao et al., “Study on enhanced oil recovery technology in low permeability heterogeneous reservoir by wateralternategas of CO_{2} flooding,” in Proceedings of the SPE Asia Pacific Oil and Gas Conference and Exhibition, SPE 165907, Jakarta, Indonesia, October 2013. View at: Google Scholar
 OPIS Europe LPG Report, 2014, http://www.opisnet.com.
 R. S. Metcalfe, “effects of impurities on minimum miscibility pressures and minimum enrichment levels for CO_{2} and richgas displa cements,” SPE Journal, vol. 22, no. 2, pp. 219–225, 1982. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Jinhyung Cho et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.