Open Access
Issue
Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles
Volume 74, 2019
Article Number 64
Number of page(s) 11
DOI https://doi.org/10.2516/ogst/2019035
Published online 30 July 2019

© F.B. Delforouz et al., published by IFP Energies nouvelles, 2019

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Minimum Miscible Pressure (MMP) is the minimum pressure at which first or multi-contact miscible displacement takes place. This parameter plays an important role in selecting miscible flooding method for Enhanced Oil Recovery (EOR) process according to the type and characteristics of oil reservoirs. An accurate estimation of MMP results in appropriate surface facilities design for gas injection, management of costs, and optimized injection pattern.

Various gases are utilized for the injection purpose, including natural gas, flue gas, nitrogen, and supercritical CO2 resulting in various levels of success in operational and economic aspects [1, 2]. Among mentioned gases, high solubility of CO2 in oil reservoir results in an extreme mass transfer between the phases [1, 2], interfacial tension reduction which increase oil sweep by reducing viscosity between the phases [3] and final recovery up to 90% [4, 5]. Injecting CO2 as a greenhouse gas and removing it from atmosphere has also environmental benefits by storing this detrimental gas under the ground [4, 68]. Another advantage of CO2 injection is the considerable reduction of MMP and higher number of potential strategies for CO2 flooding, in comparison to other gases [2]. This is caused because of CO2’s higher molecular weight in comparison to other usually used hydrocarbon gases such as Methane and Ethane and also no corrosion problems in comparison to H2S [3].

To obtain MMP, different experimental and computational methods are available. Experimental methods include slim tube test, Vanishing Interfacial Tension (VIT) technique, multi-contact mixing-cell experiment and rising bubble apparatus [1]. Computation methods contain two main groups: (1) Equation of State (EoS), (2) empirical correlations. For multicomponent injections, semi-analytical and multiple-mixing-cell methods implementing EoS would be appropriate. However, lots of steps for calibrating the EoS with respect to laboratory data must be carried out, which would be complicated. For pure gas injection or injection with small degree of impurities, empirical correlations would perform properly.

One of the most accurate experimental methods for determining MMP is the slim tube test where the oil and gas displacement process in porous medium is simulated. Due to horizontal position of slim tube and low pressure drop across the tube surface because of its small diameter, fingering phenomena and gravitational effects are eliminated. This will result in more accurate MMP measurements. On the other hand, performing slim tube experiment is substantially time and money consuming [1, 9]. Therefore, application of this in miscible injection design where a great number of MMPs should be determined would not be feasible. In fact, the major application of slim tube test is to calibrate EoS for phase equilibria calculations and to develop empirical correlations for MMP prediction in pure and impure gas injections.

In this study, we focused on empirical correlations applied for predicting MMP in pure and impure gas injections. Empirical correlations are mathematical models developed with respect to experimental data. In empirical correlations, MMP is correlated to physical parameters of oil, gas and thermodynamic conditions. Application of these models has eliminated the need for repeating time-consuming and costly experiments for each injection [10, 11].

Cronquist [12] considered reservoir temperature, C5+ molecular weight, and light components mole fraction as the key parameters of MMP prediction. Although, Lee [13], Yellig and Metcalfe [14] and Orr and Jensen [15] considered reservoir temperature (T) as the key parameter of MMP prediction.

To account for impurities in MMP calculations, correction factors were applied to MMPs obtained for pure gases. Various correction factors were introduced by researchers. Sebastian et al. [16] presented the molar average critical temperature of the mixture (Tcm) as the most accurate parameter to correlate impure MMP.

Alston et al. [5] presented a correlation for pure MMP based on reservoir temperature, C5+ molecular weight, and the ratio of light to intermediate components of the reservoir oil. Impure MMP was calculated by multiplying a correction factor based on pseudo-critical temperature. Emera and Sarma [17] presented a correlation using genetic algorithm based on reservoir temperature, C5+ molecular weight, and the ratio of light to intermediate components, similar to Alston et al. [5]. Liao predicted MMP for low permeable reservoirs. In addition, impure MMP was obtained by a parameter called the relative MMP, which was the ratio of impure to pure MMP.

Fathinasab and Ayatollahi [1] introduced a correlation for MMP prediction combining genetic programming with the multivariate search method based on reservoir temperature, C5+ molecular weight, injected gas pseudo-critical temperature, and the ratio of light to intermediate components. Using gene expression programming, Ahmadi et al. [4] predicted MMP based on parameters proposed by Fathinasab and Ayatollahi [1], Liao et al. [18], and Alston et al. [5] (T, Tcm, light to intermediate components ratio, and C5+ molecular weight). However, implementing these implicit methods in operational applications would be complex.

The above mentioned empirical correlations have significant errors for high-temperature reservoirs. Furthermore, applying correction factors to pure MMP in order to predict impure MMP would be very erroneous, since the error of pure MMP calculation can also affect impure MMP results and causes additional error.

This study aims to provide an accurate, explicit and simple empirical correlation with less computational errors compared to prior correlations for MMP prediction. The necessity to present correction factors for impure MMP prediction have been also eliminated and two separate correlations were developed to predict pure and impure CO2 MMP using multi-variable optimization algorithms. Furthermore, as a unique feature of this study, two new parameters considering the ratio of reservoir temperature to pseudo critical temperature as one parameter and acentric factor (ѡ) for the other one, have been implemented in impure MMP correlation. These two parameters can effectively determine the impact of injected gas impurities on MMP predictions.

A large data bank of oil reservoirs collected from authenticated articles [5, 7, 1624] as well as two Iranian oil fields (Darkhovin and Yadavaran) were collected and applied to develop the new correlations.

2 Data analysis/Experimental

2.1 Data analysis

Data are divided in two categories. The first category includes 126 data points which will be used to develop the correlation for pure CO2 injection. The second category includes 126 data points which will be used to develop the correlation for impure CO2 injection.

Most part of data that are used in this study have been collected from the literature [5, 7, 1624]. It is worth mentioning that previous correlations had also been developed based on these set of data [5, 7, 1624]. Therefore, a unique and new dataset is not considered to develop the correlation, except for some data points, which are added from Iranian oil fields in both pure CO2 injection (two points) and impure CO2 injection (four points). The reason to add Iranian oil fields dataset was to increase the range of temperature and because Iranian oil fields are involved with relatively high-temperature deep reservoirs.

2.2 Experimental

Iranian datasets are obtained through slim tube experimental result. Tables 1 and 2 represent the oil properties of these reservoirs.

Table 1

Properties of field oil Darkhovin.

Table 2

Properties of field oil Yadavaran.

The characteristics of the slim tube used in this study are presented in Table 3.

Table 3

Slim-tube features.

Slim tube is a stainless steel, packed with glass beds that fairly simulates one-dimensional flow through pore geometry. Before starting the test, toluene is injected in the slim tube in order to clean it and after that N2 is injected to remove the remaining amount of toluene. Moreover, a vacuum pump evacuates the porous media for several hours.

At the beginning of the fluid displacement tests, the slim tube system is saturated by the oil (reservoir fluid) at the reservoir temperature and a pressure above the bubble pressure. Then the gas with a constant flow rate is injected into the tube (1.2 pore volume) by an injection pump for miscibility process to occur. This process is repeated in several pressures and a sight glass is imbedded for the flow/process observation. There are also an accumulator and measuring systems at the end of the tube to measure gas breakthrough through, checking the producing gas–oil ratio and composition as functions of the injected volume. The schematic diagram of the slim-tube test used in this study is shown in Figure 1.

thumbnail Fig. 1

Schematic of slim-tube apparatus.

The common experimental procedure for determining MMP of CO2/Crude oil system, once the injection CO2 becomes miscible with crude oil, an inflection point is observed in the curve of recovery factor with respect of displacement pressure and the recovery will not improve as much above with a step change in pressure (Fig. 2).

thumbnail Fig. 2

Schematic graph of oil recovery versus injection pressure obtained from slim-tube test.

Table 4 represents the experimental MMP data obtained in this study. As can be seen in this table, three experimental points were obtained for each fluid sample using the slim-tube testing.

Table 4

Experimental obtained data from slim-tube testing.

3 Theory

The Group Method of Data Handling (GMDH) was developed by a Russian cybernet specialist, Prof. Alexey Ivakhnenko, in 1966. In standard regression models, the only criterion is the least squared error, and thus it cannot be determined whether the final model is simple or complex. However, using in the Ivakhnenko polynomials, one can obtain a polynomial with optimal complexity [25, 26].

The GMDH algorithm is robust and gives unique answers and produce linear and explicit correlations. This method is very suitable for solving complex and multi-dimensional problems with limited data [27] as with the case we have encountered in this article.

The GMDH algorithm creates a process for developing higher order polynomials as in equation (1): y = a + i = 1 m a i u i + i = 1 m j = 1 m a ij u i u j + i = 1 m j = 1 m k = 1 m a ijk u i u j u k + i = 1 m j = 1 m k = 1 m l = 1 m a ijkl u i u j u k u l + , $$ y=a+\sum_{i=1}^m{a}_i{u}_i+\sum_{i=1}^m\sum_{j=1}^m{a}_{{ij}}{u}_i{u}_j+\sum_{i=1}^m\sum_{j=1}^m\sum_{k=1}^m{a}_{{ijk}}{u}_i{u}_j{u}_k+\sum_{i=1}^m\sum_{j=1}^m\sum_{k=1}^m\sum_{l=1}^m{a}_{{ijkl}}{u}_i{u}_j{u}_k{u}_l+\dots, $$(1)which relates m input parameters u 1, u 2, ..., u m to a single target parameter called y. In the GMDH algorithm, it is not necessary to use all formats of the summations (double, triple, quadruple, etc.) in equation (1). In fact, it depends on the difficulties involved in modeling of a system. In systems with more variables involved, a higher order summation might be essential for accurate modeling. In this study we intended to make the correlations as simple (consisting of less constants) as possible while keeping a significant accuracy. In this regard, based on the optimization process, some variables from the regular and double summation terms are utilized for the modeling purpose as can be seen for pure CO2 MMP system (Sect. 4.1) and for impure CO2 MMP system (Sect. 4.2). Hence, implementing the triple and quadruple summations makes the correlation more complex while not improving its accuracy noticeably here.

We have divided the data into training and testing datasets. Training datasets will be implemented in developing the correlation and testing datasets will be applied for validation. To this end, 70% of data are randomly used for training and 30% of data are assigned to the test subset.

In the first step, all possible first-order polynomials are created for all existing input parameters. In this case, the constants are determined in such a way that the resulting polynomials have the least sum of squared error compared to the training data. For example, assuming only linear relationship, the following polynomials are constructed for a three-parameter function as can be seen in equation (2): y 1 = a 0 ( 1 ) + a 1 ( 1 ) u 1 y 2 = a 0 ( 2 ) + a 2 ( 2 ) u 2 y 3 = a 0 ( 3 ) + a 1 ( 3 ) u 1 + a 2 ( 3 ) u 2 y 4 = a 0 ( 4 ) + a 3 ( 4 ) u 3 y 5 = a 0 ( 5 ) + a 1 ( 5 ) u 1 + a 2 ( 5 ) u 2 y 6 = a 0 ( 6 ) + a 2 ( 6 ) u 2 + a 3 ( 6 ) u 3 y 7 = a 0 ( 7 ) + a 1 ( 7 ) u 1 + a 2 ( 7 ) u 2 + a 3 ( 7 ) u 3 . $$ \begin{array}{c}{y}_1={a}_0^{(1)}+{a}_1^{(1)}{u}_1\\ {y}_2={a}_0^{(2)}+{a}_2^{(2)}{u}_2\\ \begin{array}{c}{y}_3={a}_0^{(3)}+{a}_1^{(3)}{u}_1+{a}_2^{(3)}{u}_2\\ {y}_4={a}_0^{(4)}+{a}_3^{(4)}{u}_3\\ \begin{array}{c}{y}_5={a}_0^{(5)}+{a}_1^{(5)}{u}_1+{a}_2^{(5)}{u}_2\\ {y}_6={a}_0^{(6)}+{a}_2^{(6)}{u}_2+{a}_3^{(6)}{u}_3\\ {y}_7={a}_0^{(7)}+{a}_1^{(7)}{u}_1+{a}_2^{(7)}{u}_2+{a}_3^{(7)}{u}_3.\end{array}\end{array}\end{array} $$(2)where y k is the estimated output of kth unit, k = 1, 2, ...,7; and a i (k), i = 0, 1, 2, 3 are their weight constants. In general, M = (2 m  − 1) polynomials are constructed for m input parameters.

In the second step, for each of these polynomials, the sum of least squared error is calculated for the test data as mentioned in equation (3): d j 2 = i = 1 Number   of   data ( y ij - z i ) 2   j = 1 ,   2 ,   ,   2 m - 1 , $$ {d}_j^2=\sum_{i=1}^{\mathrm{Number}\enspace \mathrm{of}\enspace \mathrm{data}}{({y}_{{ij}}-{z}_i)}^2\enspace {j}=1,\enspace 2,\enspace \dots,\enspace {2}^m-1, $$(3)where z i represents the ith test data and y ij represents the ith corresponding prediction of the jth equation. d j is the sum of least squared error of the jth polynomial. In the third step, the polynomial with least squared error is selected as the solution. If the solution is undesired, other modes such as the division or multiplication of the current parameters can be selected as new parameters to be added to the previous parameters. Thus, the number of parameters varies and therefore the number of equations. Then the process starts again from the first step.

4 Result and discussion

4.1 Empirical correlation for pure CO2 MMP

Initial input parameters for pure MMP determination are reservoir temperature (T), MwC5+ and volatile to intermediate oil fraction (x) which is defined in equation (4) (mol stands for mole fraction of each component). Table 5 shows the range of changes in reservoir temperature (T) in °C, molar ratio of volatile (C1, N2) to intermediate components (C2, C3, C4, CO2, H2S), average molecular weight of components heavier than pentane (MwC5+), and pure MMP in MPa: x = mol C 1 + mol N 2 mol C 2 + mol C 3 + mol C 4 + mol CO 2 + mol H 2 S . $$ x=\frac{{\mathrm{mol}}_{{\mathrm{C}}_1}+{\mathrm{mol}}_{{\mathrm{N}}_2}}{{\mathrm{mol}}_{{\mathrm{C}}_2}+{\mathrm{mol}}_{{\mathrm{C}}_3}+{\mathrm{mol}}_{{\mathrm{C}}_4}+{\mathrm{mol}}_{{\mathrm{CO}}_2}+{\mathrm{mol}}_{{\mathrm{H}}_2\mathrm{S}}}. $$(4)

Table 5

Range of oil properties for pure CO2 injection.

The proposed correlation for pure CO2-MMP is shown in equation (5): MMP pure = A 1 + A 2 x Mw c 5 + + A 3 Mw c 5 + T + A 4 Mw c 5 + e x + A 5 ln ( T ) + A 6 e x ln ( T ) . $$ {\mathrm{MMP}}_{\mathrm{pure}}={A}_1+{A}_2\frac{x}{{\mathrm{Mw}}_{\mathrm{c}5+}}+{A}_3\frac{{\mathrm{Mw}}_{\mathrm{c}5+}}{T}+{A}_4{\mathrm{Mw}}_{\mathrm{c}5+}{e}^x+{A}_5\mathrm{ln}(T)+{A}_6\frac{{e}^x}{\mathrm{ln}(T)}. $$(5)

A1–A6 are constants shown in Table 6.

Table 6

Correlating parameters for pure CO2-MMP prediction.

In this study, Mean Squared Error (MSE) is used for estimating absolute deviation. Average Absolute Percentage Relative Error (AAPRE) is applied for estimating error precisely. These deviation and error measurement methods are defined as follows: E MSE = i = 1 n ( MMP i observed - MMP i calculated ) 2 n , $$ {E}_{\mathrm{MSE}}=\frac{\sum_{i=1}^n{\left({\mathrm{MMP}}_i^{\mathrm{observed}}-{\mathrm{MMP}}_i^{\mathrm{calculated}}\right)}^2}{n}, $$(6) E AAPRE = i = 1 n ( | MMP i observed - MMP i calculated MMP i calculated | ) n × 100 , $$ {E}_{\mathrm{AAPRE}}=\frac{\sum_{i=1}^n\left(\left|\frac{{\mathrm{MMP}}_i^{\mathrm{observed}}-{\mathrm{MMP}}_i^{\mathrm{calculated}}}{{\mathrm{MMP}}_i^{\mathrm{calculated}}}\right|\right)}{n}\times 100, $$(7)where MMP i observed $ {\mathrm{MMP}}_i^{\mathrm{observed}}$ is the ith observed (experimental) MMP value. MMP i calculated $ {\mathrm{MMP}}_i^{\mathrm{calculated}}$ is the ith calculated MMP value and n is the number of data points.

For pure CO2, 126 data points were available. 70% of these data (including 88 data points) were used to develop the corresponding correlation. The remaining 30% of dataset (including 38 data points) were used to evaluate and test the obtained correlation. Performance of the proposed correlation is evaluated based on each dataset (train, test and total data) as presented in Table 7 which provides results of simulation based on AAPRE and MSE.

Table 7

Performance evaluation of proposed pure MMP correlation.

On the other hand, there is no certainty if the previous authors have selected a specific part of data for test or if they do, which part of data have been used for test. Therefore, it would be reasonable to compare the performance of developed correlation with other correlations based on total data points available, not just the test data.

Thus, MSE and AAPRE calculated based on total pure data points for well-known correlations in predicting pure MMPs are given in Table 8. This table is sorted based on descending values of AAPRE.

Table 8

MSE and AAPRE error for pure MMP.

As seen in, Liao correlation [18] shows the highest error. In contrast, the correlation presented in this study has the lowest error, reducing error by at least 3.3% compared to other existing correlations.

In Figure 3 the experimental MMP graphs for three correlations with higher accuracy and minimum Error Fathinasab and Ayatollahi [1], Ahmadi et al. [4] and Emera and Sarma [17], are compared with the proposed correlation in this paper. In these graphs, the vertical and horizontal axes show the experimental MMP and the corresponding calculated MMP, respectively. Accumulation of data around the diagonal line in these graphs indicates the accuracy of each correlation in prediction of MMP. As can be seen, the graph plotted for the present study shows the best accumulation around the diagonal line y = x.

thumbnail Fig. 3

Comparison of experimental versus calculated MMP graphs for three correlations with minimum error.

4.2 Empirical correlation for impure CO2 MMP

Input parameters for impure MMP determination include C5+ molecular weight (MwC5+) of the oil, volatile to intermediate ratio (Eq. (4)), pseudo-critical pressure of the injected gas (P pc), relative pseudo-reduced temperature (reservoir temperature to pseudo critical temperature ratio) (T pr) and average acentric factor of the injected gas (ω). Where the parameters T pc, P pc are defined in equations (8) and (9) in accordance with the Kay’s rule [28]: T pc = i n y i T ci , $$ {T}_{\mathrm{pc}}=\sum_i^n{y}_i{T}_{\mathrm{ci}}, $$(8) P pc = i n y i P ci . $$ {P}_{\mathrm{pc}}=\sum_i^n{y}_i{P}_{\mathrm{ci}}. $$(9)

T pr is also defined in equation (10): T pr = T T pc . $$ {T}_{\mathrm{pr}}=\frac{T}{{T}_{\mathrm{pc}}}. $$(10)

ω is molar averaged based on equation (11) [29]: ω = i n y i ω i i n y i , $$ \omega =\frac{\sum_i^n{y}_i{\omega }_i}{\sum_i^n{y}_i}, $$(11)where y i stand for the ith component mole fraction in the gas in equations (8)(11). T ci and P ci represent critical temperature and pressure of the ith component, respectively. As mentioned earlier, 126 data points were used to develop the impure CO2 MMP empirical correlation. Range of data used for developing this correlation is given in Table 9.

Table 9

Range of oil properties for impure CO2 injection.

It is noteworthy saying that while the temperature of the published data mostly ranged from 32.2 to 118.3 °C, utilizing Iran’s data in this paper increased the temperature range according to the higher depth and temperature of these reservoirs (up to 143 °C). Moreover, the range of molecular weight of components heavier than pentane was in former studies Liao et al. [18] and Alston et al. [5], was extended between 154 and 350.3 g/gmole range of C5+ molecular weight is extended considerably using the Iran’s data. Amount of injected gas impurities for impure MMP dataset is given in Table 10.

Table 10

Range of impurities along with injected CO2.

It is worth noting that in previous studies [1, 4], T and T pc were considered as correlating parameters. In this study, a ratio of these parameters called relative pseudo-reduced temperature is chosen as a correlating parameter. This choice was due to higher observed coefficient of determination between T pr and MMP relative to MMP and T pc or T.

A higher linear coefficient of determination can be measured as a representative of the correlation between the target parameter and the input parameters. It can be observed that the data correlation for T pr and MMP is the highest, compared to T C and T (Fig. 4). Therefore, it is more suitable to develop the correlation based on T pr (Eq. (10)).

thumbnail Fig. 4

Comparison of data correlations for T, T pc, T pr.

The Pitzer’s acentric factor [30] was introduced in 1955 with the aim of developing the corresponding states theorem; increasing its reliability and accuracy in fluid properties modeling and prediction. This coefficient is defined as follows: ω = - log 10 ( P r sat ) - 1   at   T r = 0.7 , $$ \omega =-{\mathrm{log}}_{10}\left({P}_r^{\mathrm{sat}}\right)-1\enspace \mathrm{at}\enspace {T}_r=0.7, $$(12)where T r = T T c $ {T}_r=\frac{T}{{T}_c}$ is the reduced temperature and P r sat = P sat P c $ {P}_r^{\mathrm{sat}}=\frac{{P}^{\mathrm{sat}}}{{P}_c}$ is the reduced saturation vapor pressure.

Employment of acentric factor has significantly ameliorated the prediction of fluid phase behavior and calculation of reservoir fluid properties [3134].

In evolutionary process of equations of state, Soave [31] proposed a correction factor as a function of acentric factor on the attractive term of Redlich and Kwong [35] EoS which was previously introduced merely as function of temperature. The proposed format of attractive term temperature dependency is subsequently incorporated in development of equations of state. Moreover, Pitzer’s acentric factor has also played an important role in developing three-parameter equations of state, which led to significant improvements in prediction of fluid volumetric data. This impact can be clearly observed in equations of state such as Schmidt and Wenzel [33], Esmaeilzadeh and Roshanfekr [32] and Patel and Teja [34].

Success of foregone applications of acentric factor in phase behavior predictions, gave an idea to use this parameter as a correlating parameter to predict MMP alongside other previously alluded parameters.

The GMDH-based MMP function is presented as follows (Eq. (13)): MMP impure P pc = A 1 + A 2 T pr + A 3 T pr ln ( T pr ) + A 4 ω Mw c 5 + + A 5 ω x + A 6 x ω + A 7 Mw c 5 + x + A 8 Mw c 5 + ln ( T pr ) + A 9 x Mw c 5 + . $$ \frac{{\mathrm{MMP}}_{\mathrm{impure}}}{{P}_{\mathrm{pc}}}={A}_1+{A}_2{T}_{\mathrm{pr}}+\frac{{A}_3}{{T}_{\mathrm{pr}}\mathrm{ln}({T}_{\mathrm{pr}})}+{A}_4\frac{\omega }{{\mathrm{Mw}}_{\mathrm{c}5+}}+{A}_5{\omega x}+{A}_6\frac{x}{\omega }+{A}_7{\mathrm{Mw}}_{\mathrm{c}5+}x+{A}_8\frac{{\mathrm{Mw}}_{\mathrm{c}5+}}{\mathrm{ln}({T}_{\mathrm{pr}})}+{A}_9\frac{x}{{\mathrm{Mw}}_{\mathrm{c}5+}}. $$(13)

Constants A1–A8 are also listed in Table 11.

Table 11

Correlation parameters for impure MMP.

It is worth noting that, the reason for presenting two separate correlations for pure CO2 and impure CO2 is to achieve the ultimate goal to adhere simplicity and accuracy, and to avoid using correction factors for impure MMP predictions which would lead up to additional errors. As with the importance and prevalence of impure CO2 injection scenario in oil and gas industry (since pure CO2 can be hardly accessible), it is preferred to develop a separate correlation for impure CO2 injection to properly handling the simultaneous presence of accuracy and simplicity.

Moreover, using a dimensionless relative temperature (T pr), a dimensionless pressure ratio MMP / P pc $ \mathrm{MMP}/{P}_{\mathrm{pc}}$, relative molar ratio (x) and molecular weight of C5+ (MwC5+) which could be considered dimensionless despite its unit, all result in a semi-dimensionless correlation, and also reduce the number of correlation parameters.

For impure CO2, 126 data points were available. It should be emphasized that these 126 data points are completely distinct from 126 data points which were used for pure CO2. 70% of these data (including 88 data points) were used to develop the corresponding correlation. The remaining 30% of dataset (including 38 data points) were used to evaluate and test the obtained correlation. Performance of the proposed correlation is evaluated based on each dataset (train, test and total data) and presented in Table 12 which provides results of simulation based on AAPRE and MSE.

Table 12

Performance evaluation of proposed for impure correlation.

Again, there is no certainty if the previous authors have selected a specific part of data for test or if they do, which part of data have been used for test. Therefore, it would be reasonable to compare the performance of developed correlation with other correlations based on total data points available, not just the test data.

Thus, AAPRE (Eq. (7)) and MSE (Eq. (6)) errors for the above equation (13) based on total impure data points, is shown and compared with former studies in Table 13.

Table 13

MSE and AAPRE error for impure MMP.

Table 13 is sorted based on descending values of the AAPRE. It should be noted that in this study, correlations proposed by Emera and Sarma [17], Yellig and Metcalfe [14], Orr and Jensen [15], and Lee [13] were applied using the Sebastian correction factor [16] to predict the impure MMP.

Experimental MMP graphs in terms of calculated MMP for three correlations with the least error from Table 13 along with the correlation of this study were compared in Figure 5. In the present study data points are well-accumulated around the diagonal line. This indicates the higher accuracy of the correlation provided for MMP calculation.

thumbnail Fig. 5

Comparison of experimental versus calculated MMP graphs for three correlations with minimum error.

5 Sensitivity analysis

In this study, the sensitivity analysis was performed on parameters affecting both pure and impure correlations using the relevancy factor: r ( input k , MMP ) =   i = 1 n ( input k.i - input ave . k ) ( MMP i - MMP ave ) i = 1 n ( input k.i - input ave .k ) 2 i = 1 n ( MMP i - MMP ave ) 2 . $$ {r}_{\left({\mathrm{input}}_k,\mathrm{MMP}\right)}=\enspace \frac{{\sum }_{i=1}^n({\mathrm{input}}_{{k.i}}-{\mathrm{input}}_{\mathrm{ave}.k})({\mathrm{MMP}}_i-{\mathrm{MMP}}_{\mathrm{ave}})}{\sqrt{{\sum }_{i=1}^n{({\mathrm{input}}_{{k.i}}-{\mathrm{input}}_{\mathrm{ave}{.k}})}^2{\sum }_{i=1}^n{({\mathrm{MMP}}_i-{\mathrm{MMP}}_{\mathrm{ave}})}^2}}. $$(14)

In equation (14), input k.i and inputave.k are the ith value and the average value of the kth input, respectively. The index k refers to each enumerating parameter, e.g. Temperature, volatile to intermediate ratio, etc; MMP i stands for the ith value of predicted MMP and MMPave is the arithmetic average of predicted MMP values. r shows the effect of each parameter (each input k ) on the correlation output (MMP in this study). If r > 0, then the associated parameter has a positive effect; in contrast, if r is negative (r < 0), the associated parameter has a negative effect. The parameter r ranges from −1 to 1, indicating the highest negative or positive effect.

The results of sensitivity analysis on pure MMP correlation are depicted in Figure 6. As can be seen, all parameters (including x, MwC5+, T) have a positive effect on predicted MMP. Temperature has the largest effect and C5+ molecular weight of components has the smallest effect. The impact of these parameters in this study is consistent with the results of foregone studies [7, 18, 36].

thumbnail Fig. 6

Relevancy factor for pure MMP.

The results of sensitivity analysis are rendered for the correlation of impure MMP parameters (ω, P pc, MwC5+, x and T pr) in Figure 7 (T pr consists of two parameters, T and T pc.). Previous researchers reported direct relation of T [7, 18, 36] and inverse relation of T pc [7]. Consequently, the effect of their ratio (T/T pc) reported as T pr is ultimately deduced to have a direct relation with MMP. The effect of the parameter x and C5+ molecular weight is similar to pure MMP. Sensitivity analysis on correlating parameters in this study confirms these results.

thumbnail Fig. 7

Relevancy factor for impure MMP considering acentric factor.

A gas with high molecular weight has a lower P pc. As a result, one can assume that increasing the molecular weight of gas has the same effect as decreasing P pc. By increasing the molecular weight of gas, P pc and hence the MMP decreases. The results of the sensitivity analysis also confirm this conclusion.

In general, with increasing the gas molecular weight, the acentric factor (ω) increases and accordingly, MMP decreases. Therefore, inverse relation of the acentric factor on MMP can also be explained and interpreted. Among the parameters mentioned above, T pr, P pc, and ѡ showed the most while x and MWC5+ showed the least impact.

6 Conclusion

Reviewing experimental data in foregone studies and data from two Iranian reservoirs, a database with a higher temperature range was collected. Then, to predict MMP in pure and impure CO2 injection operations, correlations were provided using the GMDH algorithm.

  1. The correlation proposed to predict pure MMP is an explicit correlation based on MWC5+, x, and T parameters. This correlation ameliorates the results compared to other previous correlations and reduces the computational error by at least 2.5%. Remarkable decrease of computational error corroborates robustness of the GMDH algorithm approach.

  2. The correlation proposed to predict impure MMP is presented explicitly without using the correction factor. The effective parameters in this correlation include T pr, x, MwC5+ and ω. Using the dimensionless temperature and presenting the correlation as MMP/P pc, this correlation is developed in a semi-dimensionless form. It should be noted that such non-dimensionalization not only shortens the correlation but also provides its applicability for many reservoirs. Employing the GMDH algorithm approach along with implementing the gas acentric factor, eventuated in at least 2% decrease in computational error compared to previous studies.

  3. Since the data from Iranian reservoirs were used to develop new correlations in this study, they can be used as means for predicting CO2-MMP in reservoirs with high depth and temperatures such as some Middle East reservoirs.

References

  • Fathinasab M., Ayatollahi S. (2016) On the determination of CO2–crude oil minimum miscibility pressure using genetic programming combined with constrained multivariable search methods, Fuel 173, 180–188. [CrossRef] [Google Scholar]
  • Shaw J., Bachu S. (2002) Screening, evaluation, and ranking of oil reservoirs suitable for CO2-flood EOR and carbon dioxide sequestration, J. Can. Pet. Technol. 41, 51–61. [CrossRef] [Google Scholar]
  • Rahimi V., Bidarigh M., Bahrami P. (2017) Experimental study and performance investigation of miscible water-alternating-CO2 flooding for enhancing oil recovery in the Sarvak formation, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles 72, 35. [CrossRef] [Google Scholar]
  • Ahmadi M.A., Zendehboudi S., James L.A. (2017) A reliable strategy to calculate minimum miscibility pressure of CO2-oil system in miscible gas flooding processes, Fuel 208, 117–126. [CrossRef] [Google Scholar]
  • Alston R., Kokolis G., James C. (1985) CO2 minimum miscibility pressure: A correlation for impure CO2 streams and live oil systems, Soc. Pet. Eng. 25, 268–274. [CrossRef] [Google Scholar]
  • Huang Y., Huang G., Dong M., Feng G. (2003) Development of an artificial neural network model for predicting minimum miscibility pressure in CO2 flooding, J. Pet. Sci. Eng. 37, 83–95. [Google Scholar]
  • Kamari A., Arabloo M., Shokrollahi A., Gharagheizi F., Mohammadi A.H. (2015) Rapid method to estimate the minimum miscibility pressure (MMP) in live reservoir oil systems during CO2 flooding, Fuel 153, 310–319. [CrossRef] [Google Scholar]
  • Sun L., Hao X., Dou H., Adenutsi C.D., Li Z., Zhang Y. (2018) Co-optimization of oil recovery and CO2 storage for cyclic CO2 flooding in ultralow permeability reservoirs, Oil Gas Sci. Technol. - Rev. IFP Energies nouvelles 73, 42. [CrossRef] [Google Scholar]
  • Glaso O. (1985) Generalized minimum miscibility pressure correlation (includes associated papers 15845 and 16287), Soc. Pet. Eng. 25, 927–934. [CrossRef] [Google Scholar]
  • Danesh A. (1998) PVT and phase behaviour of petroleum reservoir fluids, Elsevier, The Netherlands. [Google Scholar]
  • Kharrat R., Asadollahi M., Namani M. (2008) Enhanced oil recovery, Bahr Danesh, Tehran, Iran. [Google Scholar]
  • Cronquist C. (1978) Carbon dioxide dynamic miscibility with light reservoir oils), in: Proc. Fourth Annual US DOE Symposium, Tulsa. [Google Scholar]
  • Lee J.I. (1979) Effectiveness of carbon dioxide displacement under miscible and immiscible conditions, Report RR-40, Pet. Rec. Inst. Calg., Canada. [Google Scholar]
  • Yellig W., Metcalfe R. (1980) Determination and prediction of CO2 minimum miscibility pressures (includes associated paper 8876), J. Pet. Technol. 32, 160–168. [CrossRef] [Google Scholar]
  • Orr F. Jr, Jensen C. (1984) Interpretation of pressure-composition phase diagrams for CO2/crude-oil systems, Soc. Pet. Eng. 24, 485–497. [CrossRef] [Google Scholar]
  • Sebastian H., Wenger R., Renner T. (1985) Correlation of minimum miscibility pressure for impure CO2 streams, J. Pet. Technol. 37. 2,076–072,082. [Google Scholar]
  • Emera M.K., Sarma H.K. (2005) Use of genetic algorithm to estimate CO2–oil minimum miscibility pressure – A key parameter in design of CO2 miscible flood, J. Pet. Sci. Eng. 46, 37–52. [Google Scholar]
  • Liao C., Liao X., Chen J., Ye H., Chen X., Wang H. (2014) Correlations of minimum miscibility pressure for pure and impure CO2 in low permeability oil reservoir, J. Energy Inst. 87, 208–214. [CrossRef] [Google Scholar]
  • Yuan H., Johns R., Egwuenu A., Dindoruk B. (2004) Improved MMP correlations for CO2 floods using analytical gas flooding theory, in: SPE/DOE Symposium on Improved Oil Recovery, Society of Petroleum Engineers. [Google Scholar]
  • Shokir E.M.E.-M. (2007) CO2–oil minimum miscibility pressure model for impure and pure CO2 streams, J. Pet. Sci. Eng. 58, 173–185. [Google Scholar]
  • Bon J., Emera M.K., Sarma H.K. (2006) An experimental study and genetic algorithm (GA) correlation to explore the effect of nC5 on impure CO2 minimum miscibility pressure (MMP), in: SPE Asia Pacific Oil & Gas Conference and Exhibition, Society of Petroleum Engineers. [Google Scholar]
  • Eakin B., Mitch F. (1988) Measurement and correlation of miscibility pressures of reservoir oils, in: SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers. [Google Scholar]
  • Moosazadeh M., Keshavarzi B., Ghotbi C. (2017) Investigation of the minimum miscibility pressure for injection of two different gases into two Iranian oil reservoirs: Experimental and theory, Cad. J Chem. Eng. 95, 1014–1020. [CrossRef] [Google Scholar]
  • Zuo Y.-X., Chu J.-Z., Ke S.-L., Guo T.-M. (1993) A study on the minimum miscibility pressure for miscible flooding systems, J. Pet. Sci. Eng. 8, 315–328. [Google Scholar]
  • Farlow S.J. (1981) The GMDH algorithm of Ivakhnenko, Am. Stat. 35, 210–215. [Google Scholar]
  • Madala H.R., Ivakhnenko A.G. (1994) Inductive learning algorithms for complex systems modeling, CRC Press, Boca Raton, FL. [Google Scholar]
  • Ivakhnenko A. (1970) Heuristic self-organization in problems of engineering cybernetics, Automatica 6, 207–219. [Google Scholar]
  • Kay W. (1936) Density of hydrocarbon gases and vapors at high temperature and pressure, Ind. Eng. Chem. 28, 1014–1019. [Google Scholar]
  • Joergensen M., Stenby E. (1995) Optimization of pseudo-component selection for compositional studies of reservoir fluids, Society of Petroleum Engineers. [Google Scholar]
  • Curl R., Pitzer K. (1958) Volumetric and thermodynamic properties of fluids – enthalpy, free energy, and entropy, Ind. Eng. Chem. 50, 265–274. [Google Scholar]
  • Soave G. (1972) Equilibrium constants from a modified Redlich-Kwong equation of state, Chem. Eng. Sci. 27, 1197–1203. [Google Scholar]
  • Esmaeilzadeh F., Roshanfekr M. (2006) A new cubic equation of state for reservoir fluids, Fluid Phase Equil. 239, 83–90. [CrossRef] [Google Scholar]
  • Schmidt G., Wenzel H. (1980) A modified van der Waals type equation of state, Chem. Eng. Sci. 35, 1503–1512. [Google Scholar]
  • Patel N.C., Teja A.S. (1982) A new cubic equation of state for fluids and fluid mixtures, Chem. Eng. Sci. 37, 463–473. [Google Scholar]
  • Redlich O., Kwong J.N. (1949) On the thermodynamics of solutions. V. An equation of state. Fugacities of gaseous solutions, Chem. Rev. 44, 233–244. [CrossRef] [PubMed] [Google Scholar]
  • Chen G., Fu K., Liang Z., Sema T., Li C., Tontiwachwuthikul P., Idem R. (2014) The genetic algorithm based back propagation neural network for MMP prediction in CO2-EOR process, Fuel 126, 202–212. [CrossRef] [Google Scholar]

All Tables

Table 1

Properties of field oil Darkhovin.

Table 2

Properties of field oil Yadavaran.

Table 3

Slim-tube features.

Table 4

Experimental obtained data from slim-tube testing.

Table 5

Range of oil properties for pure CO2 injection.

Table 6

Correlating parameters for pure CO2-MMP prediction.

Table 7

Performance evaluation of proposed pure MMP correlation.

Table 8

MSE and AAPRE error for pure MMP.

Table 9

Range of oil properties for impure CO2 injection.

Table 10

Range of impurities along with injected CO2.

Table 11

Correlation parameters for impure MMP.

Table 12

Performance evaluation of proposed for impure correlation.

Table 13

MSE and AAPRE error for impure MMP.

All Figures

thumbnail Fig. 1

Schematic of slim-tube apparatus.

In the text
thumbnail Fig. 2

Schematic graph of oil recovery versus injection pressure obtained from slim-tube test.

In the text
thumbnail Fig. 3

Comparison of experimental versus calculated MMP graphs for three correlations with minimum error.

In the text
thumbnail Fig. 4

Comparison of data correlations for T, T pc, T pr.

In the text
thumbnail Fig. 5

Comparison of experimental versus calculated MMP graphs for three correlations with minimum error.

In the text
thumbnail Fig. 6

Relevancy factor for pure MMP.

In the text
thumbnail Fig. 7

Relevancy factor for impure MMP considering acentric factor.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.