Open Access
Issue
Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles
Volume 72, Number 1, January–February 2017
Article Number 2
Number of page(s) 11
DOI https://doi.org/10.2516/ogst/2016023
Published online 26 January 2017

© H. Aguilar-Cisneros et al., published by IFP Energies nouvelles, 2017

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.

Notation

A : Parameter of the equation of state

A1-A3: Coefficients in Equation (5)

a1-a4: Coefficients in Equations (11)-(14)

B : Parameter of the equation of state

b1-b4: Coefficients in Equations (11)-(14)

C k : Carbon number of defined hydrocarbons

c1: Coefficient in Equation (11)

d1: Coefficient in Equation (11)

FG : Number of functional groups types

FObj: Objective Function

H/C: Hydrogen carbon ratio

MW : Molecular Weight

N : Number of data points

P : Pressure

Pc: Critical Pressure

pc: Contribution to Pc by functional groups

PSat: Saturation Pressure

R : Universal constant of gases, 8.31441 J/mol K

RAD%: $ \mathrm{Relative}\enspace \mathrm{Average}\enspace \mathrm{Deviation}=\frac{100}{N}\sum_{i=1}^N\left(\frac{\left|{x}_i^{\mathrm{Exp}}-{x}_i^{\mathrm{Cal}}\right|}{{x}_i^{\mathrm{Exp}}}\right)$

T : Temperature

Tc: Critical temperature

V : Molar volume, m3/mol

Vc: Critical volume, m3/mol

vCi: Contribution to Vc by functional groups

x : Mole fraction

Z : Compressibility factor

Greek letters

α : Coefficient in Equation (9)

β : Coefficient in Equation (9)

χ : Coefficient in Equation (9)

Δv : Group volume increment in Equations (4) and (5)

Φ : Fugacity coefficient

ρ : Mass density, kg/m3

ν : Number of functional groups

ω : Acentric factor

Superscripts

Cal: Calculated data

L: Liquid phase

Exp: Experimental data

Subscripts

c: Critical property

k : Carbon number k

i : Functional group i or pure component i

j : Undefined petroleum fractions j or pure component j

Introduction

Hydrogen solubility predictions are essential part in several processes in the petroleum and chemical industries; for instance, hydrocracking and hydrogenation processes are used to transform heavy oil into more valuable products. Hydrogen solubility data are required for designing and operating such processes and it is also needed in the corresponding kinetic models (Saajanlehto et al., 2014). Several authors have pointed out the importance of taking into account Vapor-Liquid Equilibria (VLE) in reaction modeling and its influence in conversion predictions (Chen et al., 2011; Hook, 1985; Pellegrini et al., 2008). A detailed review about VLE of hydrogen and petroleum fraction systems was presented by Chávez et al. (2014). They presented experimental data and modeling methods, most of which were conducted for systems composed of hydrogen and defined hydrocarbons, with few systems formed by hydrogen and undefined petroleum fractions. This may be attributable to the lack of experimental data, making a hard task to develop accurate prediction methods. Besides, binary systems containing hydrogen all exhibit a type III phase behavior in the classification scheme of Van Konynenburg and Scott (Privat and Jaubert, 2013) which cannot be predicted with cubic equations of state using classical mixing rules because of the quantum nature of hydrogen and the mixture size-asymmetry (Deiters, 2013). Despite these issues, most methods found in literature use cubic Equations Of State (EOS) because of their mathematical simplicity and their capability of modeling multicomponent mixtures over a wide range of temperatures and pressures. To overcome such limitations several authors have proposed to fit EOS parameters to VLE data (Qian et al., 2013); however, in general, their capability of prediction is reduced with such regressions. New approaches have been suggested in order to preserve the predictive features of cubic equations of state, such as that described by Jaubert et al. (2010). In such approach, binary interaction parameters are predicted by using a group contribution method, where classical cubic equations of state and excess free energy models (EOS/gE) are combined, being applied successfully to binary mixtures of hydrogen with defined hydrocarbons (Qian et al., 2013).

A further challenge is the mathematical characterization of the undefined petroleum fraction; that is, the determination of its physical-chemical properties by mathematical means. In such characterization process, it is important to find a convenient number of pseudocomponents in which the undefined petroleum fraction will be divided. According to Lin et al. (1985), better results are obtained as the number of pseudocomponents increases when VLE of systems formed by hydrogen and undefined fractions are modeled. On the other hand, Chávez concluded that 30 pseudocomponents are enough to obtain good agreement between calculated and experimental data. The characterization procedure also involves estimation of the corresponding critical properties entering into cubic equations of state. These properties are frequently calculated from correlations (Firoozabadi, 1988; Lee and Kesler, 1975; Pedersen et al., 2004), but they do not always lead to the same results, because each one was developed for different sets of fluids at limited conditions. Thus, recently, group contribution methods have been recently extended to pseudocomponents, showing a good performance (Carreón-Calderón et al., 2012; Xu et al., 2015).

Carreón-Calderón et al. (2012) presented a new characterization procedure based on group contribution methods. They used selected functional groups to assign a hypothetical chemical structure to an undefined petroleum fraction by a minimization process of its corresponding Gibbs free energy. This methodology was applied to simulate VLE of heavy (Carreón-Calderón et al., 2014) and light hydrocarbons (Uribe-Vargas et al., 2016). Good results were achieved in both cases without making use of correlations to calculate critical properties or adjusted binary interaction parameters. In particular, they modeled vapor-liquid equilibria of heavy petroleum fluids, where no significant differences were shown regardless of the number of pseudocomponents used in modeling. These results show the feasibility of applying this characterization procedure to VLE of hydrogen with heavy undefined petroleum fraction. In this work, we add hydrogen functional group to the set of functional groups original proposed and suggest modifications to the characterization procedure in order to determine the hydrogen solubility properly. The entire methodology was tested with experimental data found in literature: Venezuelan Heavy Coking Gas Oil (HGCO) (Ji et al., 2013), Karamay Atmospheric Residue (KRAR) (Ji et al., 2013), Liaohe Atmospheric Residue (LHAR) (Ji et al., 2013), Venezuelan Atmospheric Residue (VNAR) (Ji et al., 2013), Canadian Light Virgin Gas Oil (CLVGO) (Cai et al., 2001), Canadian Heavy Virgin Gas Oil (CHVGO) (Cai et al., 2001) and Athabasca Bitumen (AB) (Lal et al., 1999).

1 Characterization Procedure

The approach presented by Carreón-Calderón et al. (2012) is used to determine the properties of the undefined petroleum fractions. In that methodology, a hypothetical structure is calculated with a set of selected functional groups, which depend on the group contribution method used. Carreón-Calderón et al. (2012, 2014) have suggested two group contribution methods; Joback-Reid (JR) (Joback and Reid, 1987) and Marrero-Gani (MG) (Marrero and Gani, 2001). JR method was suggested to simulate VLE of light petroleum fluids (Uribe-Vargas et al., 2016), while MG method showed better functionality in heavy petroleum fluids heavier than 14 API (Carreón-Calderón et al., 2014). In the latter, the C7+ fraction, which has unknown composition, was split in 3, 6 and 9 pseudocomponents, founding no significant difference between calculations and experiments regardless of the fraction splitting. Accordingly, in this work, the undefined petroleum fractions are treated as a single pseudocomponent in order to have a simpler characterization, where their corresponding properties are calculated through MG method.

Table 1 shows the functional groups of the MG method and the equivalent groups of the GCVOL (Elbro et al., 1991; Ihmels and Gmehling, 2003) method for determining liquid densities, as originally proposed by Carreón-Calderón et al. (2012). A quick review of these functional groups was made and it was found that the last functional groups are not equivalent each other (aC ≠ AC − C). Therefore, in this work the last functional group (aC) was changed to aC − C functional group from the same MG method.

Table 1

Original functional groups.

An optimization process is used to estimate a hypothetical chemical structure of the undefined petroleum fraction, where its Gibbs free energy is the objective function to be minimized. Since the undefined fraction is treated like a pure component, minimize the liquid fugacity coefficient $ \left({\phi }_i^L\right)$ is equivalent to minimize the Gibbs energy. Thus, the objective function is: min υ 1 , υ 2 , , υ FG     0   ϕ i L [ T C i ( υ 1 , υ 2 , , υ FG ) , P C i ( υ 1 , υ 2 , , υ FG ) , ω i ( υ 1 , υ 2 , , υ FG ) ] $$ \underset{{\upsilon }_1,{\upsilon }_2,\cdots,{\upsilon }_{{FG}}\enspace \ge \enspace 0}{\mathrm{min}}\enspace {\phi }_i^L\left[{T}_{\mathrm{C}i}\left({\upsilon }_1,{\upsilon }_2,\cdots,{\upsilon }_{{FG}}\right),{P}_{\mathrm{C}i}\left({\upsilon }_1,{\upsilon }_2,\cdots,{\upsilon }_{{FG}}\right),{\omega }_i\left({\upsilon }_1,{\upsilon }_2,\cdots,{\upsilon }_{{FG}}\right)\right] $$(1)where TC, PC and ω are the critical temperature, critical pressure and acentric factor, respectively, with FG being the number of functional groups types. The $ {\phi }_i^L$ expression depends on the cubic EOS selected; in this work, the Peng-Robinson cubic Equation Of State PR-EOS (Peng and Robinson, 1976) is used to determine the fugacity coefficient as follows: ln   ( ϕ i L ) = ( Z L - 1 )   ln   ( Z L - B ) - A 2 2   B ln ( Z L + ( 2 + 1 ) B Z L - ( 2 - 1 ) B ) $$ \mathrm{ln}\enspace \left({\phi }_i^L\right)=\left({Z}^L-1\right)\enspace \mathrm{ln}\enspace \left({Z}^L-B\right)-\frac{A}{2\sqrt{2}\enspace B}\mathrm{ln}\left(\frac{{Z}^L+\left(\sqrt{2}+1\right)B}{{Z}^L-\left(\sqrt{2}-1\right)B}\right) $$(2)where Z L  = PV L ⁄(RT) is the compressibility factor of the liquid phase, A and B are PR-EOS parameters. The properties of undefined petroleum fractions are usually given at standard conditions; hence, the pressure P and temperature T are set equal to 0.101325 MPa and 288.15 K, respectively; v L is the molar liquid volume and R is the gas constant.

This minimization problem is constrained by the molecular weight and density of the corresponding undefined petroleum fraction: j = 1 FG υ j   MW j - MW i = 0 $$ \sum_{j=1}^{{FG}}{\upsilon }_j\enspace {{MW}}_j-{{MW}}_i=0 $$(3) j = 1 FG υ j   Δ v j - MW i ρ i L = 0 $$ \sum_{j=1}^{{FG}}{\upsilon }_j\enspace \Delta {v}_j-\frac{{{MW}}_i}{{\rho }_i^L}=0 $$(4)

In the above expression, subscripts j and i indicate the functional group j of the undefined fraction i, MW is the molecular weight, υj, the number (coefficient) of a functional group, Δvj the volume increment, $ {\rho }_i^L\enspace $the liquid density. The Δvj is expressed as a function of temperature by means of Δ v j = A 1 + A 2 T + A 3 T 2 $$ \Delta {v}_j={A}_1+{A}_2T+{A}_3{T}^2 $$(5)where A1, A2 and A3 are the parameters of the GCVOL approach (Ihmels and Gmehling, 2003).

Finally, to solve the minimization problem, Equation (2) requires critical properties (TC, PC and ω). These properties are calculated using Equations from (6) to (8), where, the temperature is in K, pressures in MPa and volume in m3/mol (Marrero and Gani, 2001). T C = 231.239   ln ( i = 1 FG υ i   t C i ) $$ {T}_{\mathrm{C}}=231.239\enspace \mathrm{ln}\left(\sum_{i=1}^{{FG}}{\upsilon }_i\enspace {t}_{\mathrm{C}i}\right) $$(6) P C = 0.1 ( i = 1 FG υ i   p C i + 0.108998 ) - 2 + 0.59827 $$ {P}_{\mathrm{C}}=0.1{\left(\sum_{i=1}^{{FG}}{\upsilon }_i\enspace {p}_{\mathrm{C}i}+0.108998\right)}^{-2}+0.59827 $$(7) V C = ( 7.95 + i = 1 FG υ i   v C i ) × 10 - 6 $$ {V}_{\mathrm{C}}=\left(7.95+\sum_{i=1}^{{FG}}{\upsilon }_i\enspace {v}_{\mathrm{C}i}\right)\times {10}^{-6} $$(8)

In previous equations tCi, pCi and vCi represent contributions to the critical temperature, the critical pressure, and the critical molar volume because of the presence of a defined functional group i in a molecule.

One of the major contributions of the characterization procedure described above is the assignment of a hypothetical chemical structure to an undefined petroleum fraction. This allows us to handle the undefined fraction as hypothetical pure component and, in principle, to make use of any thermodynamic model exclusively developed for known components and their mixtures. Thus, in addition to PR-EOS, Huron-Vidal (Michelsen, 1990) mixing rule together with the UNIversal quasi-chemical Functional group Activity Coefficients (UNIFAC) approach, as activity coefficient model, are used also in this work for calculating VLE of hydrogen and heavy petroleum fluids, where the group-interaction parameters Ψnm in the UNIFAC approach are given by Ψ nm = e xp - α nm + β nm T + χ nm T 2 T $$ {\Psi }_{{nm}}=\mathrm{e}\mathrm{xp}-\frac{{\alpha }_{{nm}}+{\beta }_{{nm}}T+{\chi }_{{nm}}{T}^2}{T} $$(9)

The parameters αnm, βnm and χnm are specific parameters for a couple of functional groups n and m, whose values are obtained from literature (Horstmann et al., 2005). Notice that αnm ≠ αmn, βnm ≠ βmn, χnm ≠ χmn. It is worth pointing out that hydrogen is handled as a single functional group.

2 Solubility of Hydrogen in Defined Hydrocarbon Components

Before proceeding with the calculation of VLE of systems composed of hydrogen and heavy petroleum fluids, the hydrogen solubility in defined hydrocarbons was estimated in order to evaluate the thermodynamic models and their corresponding binary parameters described above.

Table 2 shows the 13 binary systems involved in this study, which are formed by nine paraffinic and four aromatic compounds. The range of experimental temperatures and pressures are also included and the comparison of calculated with experimental values are given in Table 3.

Table 2

List of 13 binary mixtures used in this study.

Table 3

RAD of hydrogen solubility in defined hydrocarbons.

As observed, hydrogen solubility estimation is smaller than the experimental values in all cases. The difference is higher as the carbon number of the component increases; the Relative Average Deviation (RAD) increases from 13% for pentane to 75% for hexatetracontane, whereas, in the case of aromatic compounds increases from 32% to 52% for benzene and pyrene, respectively. Because of these results, group-interaction parameters Ψnm were calculated again by fitting calculations to the available VLE experimental data in Table 2. This fitting process requires to identify the functional groups forming the defined components. In case of paraffinic compounds, they can be represented by -CH3 and -CH2 functional groups, and aromatic compounds by ACH and AC functional groups. Because -CH3 and AC belongs to the main group 1 (-CH2) and 3 (ACH) respectively, only interaction parameters between -CH2 and H2, ACH and H2, and ACH and CH2 are required.

Table 4 shows the original αnm parameters between main functional groups.

Table 4

Original group-interaction parameters (αnm, αmn) between functional groups used in this study (Horstmann et al., 2005).

In the fitting process, for simplicity, parameter αnm is adjusted while βnm and χnm parameters keep their original values. The fitted group-interaction parameters are obtained by a minimization process, where the objective function (FObj) is given by F Obj = i = 1 N [ P Sat i Exp - P Sat i Cal ( α nm , α mn ) ] 2 N $$ {F}_{\mathrm{Obj}}=\frac{\sqrt{\sum_{i=1}^N{\left[{P}_{{\mathrm{Sat}}_i^{\mathrm{Exp}}}-{P}_{{\mathrm{Sat}}_i^{\mathrm{Cal}}}({\alpha }_{{nm}},{\alpha }_{{mn}})\right]}^2}}{N} $$(10)

Here $ {P}_{{\mathrm{Sat}}_i^{\mathrm{Exp}}}$ and $ {P}_{{\mathrm{Sat}}_i^{\mathrm{Cal}}}$ are the experimental and calculated saturation pressure on the bubble point respectively and N represent the number of data points. The result of this fitting process is a group of new parameters (αnm, αmn) for each binary mixtures (H2-hydrocarbon). Figures 1a and 1b show the trend of αnm and αmn parameters for paraffinic and aromatic compounds as function of the carbon number. The correlated parameters take the form of: α H 2 - C H 2 = a 1 ln ( b 1 C k + c 1 ) + d 1 $$ {\mathrm{\alpha }}_{{\mathrm{H}}_2-\mathrm{C}{\mathrm{H}}_2}={a}_1\mathrm{ln}\left({b}_1{\mathrm{C}}_k+{c}_1\right)+{d}_1 $$(11) α C H 2 - H 2 = a 2 ln ( C k ) + b 2 $$ {\mathrm{\alpha }}_{\mathrm{C}{\mathrm{H}}_2-{\mathrm{H}}_2}={a}_2\mathrm{ln}\left({\mathrm{C}}_k\right)+{b}_2 $$(12) α H 2 - ACH = a 3 ln ( C k ) + b 3 $$ {\mathrm{\alpha }}_{{\mathrm{H}}_2-\mathrm{ACH}}={a}_3\mathrm{ln}\left({\mathrm{C}}_k\right)+{b}_3 $$(13) α ACH - H 2 = a 4 ln ( C k ) + b 4 $$ {\mathrm{\alpha }}_{\mathrm{ACH}-{\mathrm{H}}_2}={a}_4\mathrm{ln}\left({\mathrm{C}}_k\right)+{b}_4 $$(14)where, Ck is the carbon number of the defined compounds. $ {\alpha }_{{\mathrm{H}}_2-\mathrm{C}{\mathrm{H}}_2}$ and $ {\alpha }_{{\mathrm{CH}}_2-{\mathrm{H}}_2}$are the group-interaction parameters between -CH2 and H2 functional groups, $ {\alpha }_{{\mathrm{H}}_2-\mathrm{ACH}}$ and $ {\alpha }_{\mathrm{ACH}-{\mathrm{H}}_2}$ are group-interaction parameters between ACH and H2 functional groups. Finally, coefficients an, bn, cn, dn are shown in Table 5.

thumbnail Figure 1

Group-interaction parameters as function of carbon number a) paraffinic compounds, b) aromatic compounds.

Table 5

Group-interaction parameters (αnm) between functional groups used in this study.

Similar results have been obtained by Florusse et al. (2003); they used a Statistical Associating Fluid Theory (SAFT) model to predict hydrogen solubility in heavy alkanes, where the corresponding cross-binary interaction parameters were obtained as a function of the carbon number as well.

Figures 2a and 2b show experimental (♦) and calculated hydrogen solubility data for pyrene and hexatetracontane, respectively. The original group-interaction parameters underestimate hydrogen solubility (-----) whereas the predicted values are in better agreement with experimental data (—).

thumbnail Figure 2

Solubility of hydrogen a) pyrene at 433.2 K, b) hexatetracontane at 402.4 K.

Table 3 summarizes the RAD of all defined components with the original group-interaction parameters and with those given by Equations from (11) to (14).

3 Solubility of Hydrogen in Undefined Petroleum Fractions

As first step, critical properties and molecular pseudostructures of the seven undefined petroleum fractions studied in this work were determined using the characterization procedure described in Section 1. Table 6 shows density, molecular weight and hydrogen-carbon ratio of the seven undefined fractions used in this work. As will be discussed later, the last two properties play an important role in modeling such systems with hydrogen. In all cases, the experimental error reported in the literature for hydrogen solubility was ±5% (Cai et al., 2001; Ji et al., 2013), except for AB system where such experimental error was not reported (Lal et al., 1999).

Table 6

Relevant physical properties of the undefined petroleum fractions used in this work.

Once the molecular pseudostructure is formed with the functional groups from Table 1, it is possible to make a distinction between paraffinic, naphthenic and aromatic contributions through to the no-integer coefficients (υj) resulting from the minimization process. Notice that although nine functional groups are used in the characterization procedure (see Tab. 1), there are only two main functional groupS according to the UNIFAC method: first six functional belong to the main -CH2 functional group, whereas the rest of them belong to the main ACH functional group. Accordingly, the calculated coefficients υj can be added in order to obtain the carbon number according to these two main functional groups. Therefore, Equations from (11) to (14) can be used to estimate group-interaction parameters between the H2 and the functional groups calculated by the characterization procedure.

As illustration, Figure 3 shows the experimental and predicted hydrogen solubility in CHVGO (at 603.15 K) which has a molecular weight of 350 g/mol. The black diamonds (♦) represent the experimental values and the crossed bars represent the experimental reported error of ±5%. Calculations using the original group-interaction parameters (— —) underestimates the hydrogen solubility, as it is also observed for the defined hydrocarbons from Table 2; in this case the RAD is 51%. When our proposed correlations, Equations from (11) to (14), are used to estimate group-interaction parameters, the predicted values (⋅⋅⋅⋅) remains below experimental data, but the RAD is reduced up to 25%. To improve the accuracy of hydrogen solubility prediction, a modification of the original minimization process of characterization is suggested. This modification consists of including another constraint, in addition to Equations (3) and (4). We believed that this constraint should somehow account for the chemical nature of the undefined fractional. Hence, the Hydrogen-Carbon ratio (H/C) is proposed as a third constraint. According to experimental observations made by Ji et al. (2013), this parameter may be related to the aromaticity of undefined fraction.

thumbnail Figure 3

Solubility of hydrogen in CHVGO petroleum fraction at 603.15 K as a function of group-interaction parameters.

Figure 4 shows the coefficients (υj) obtained from the original characterization procedure and the modified one with the hydrogen-carbon ratio as constraint. The latter gives smaller values for the aromatic functional groups and larger values for all other functional groups. According to Park et al. (1996) hydrogen solubility decreases as the number of aromatic rings increases. Hence, we believe that the inclusion of the H/C ratio may help to construct a pseudostructure in better agreement with its aromatic nature. Figure 3 shows the predicted hydrogen solubility (—) using the modified characterization procedure. As seen, predicted hydrogen solubilities are in better agreement with experimental data (RAD = 5%).

thumbnail Figure 4

Coefficients calculated by characterization process.

The modified characterization procedure is summarized in Figure 5, where the stages are described in a flowchart. Here (H⁄C)j is the hydrogen-carbon ratio of a functional group j and (H⁄C)Exp is the experimental hydrogen-carbon ratio of a i undefined petroleum fraction.

thumbnail Figure 5

Flow chart of the modified characterization process.

As observed in Table 6, three undefined fractions have molecular weights below 360 g/mol, whereas the rest of them have values above 500 g/mol. Although good results were achieved in the former, hydrogen solubility was overestimated in the latter. When the calculated pseudostructures are reviewed, it is found that they are formed mainly by paraffinic and cyclic functional groups and by a small portion of aromatic ones. This seems to be the reason why the hydrogen solubility is overestimated. Hence, a sensitivity analysis was made to know the effect of the H/C ratio in the characterization process from Figure 5.

In Figure 6a, three different pseudostructures of AB depending on the H/C ratio are shown. When H/C parameter is set equal to the experimental value, the pseudostructure is mainly formed by paraffinic and naphthenic functional groups and hydrogen solubility is overestimated (RAD > 50%). In contrast, hydrogen solubility is underestimated (RAD > 20%) when H/C parameter is set equal to 30% below the experimental value. In latter case, aromatic functional groups predominate over paraffinic and naphthenic ones. Figure 6b shows the trend of the RAD as function of the used H/C ratio. It can be seen that there is minimum value when H/C parameter is set around 80% of the experimental value; its corresponding pseudostructure is shown in Figure 6a. At this point, it is important to emphasize that the characterization procedure do not intend to make a characterization at molecular level, but look for an average structure that allow to reproduce bulk properties (Carreón-Calderón et al., 2012, 2014; Uribe-Vargas et al., 2016).

thumbnail Figure 6

Sensitivity analysis of H/C parameter used in the characterization process for AB. a) Pseudostructure arrangements of AB depending on percent of the H/C ratio, b) RAD vs. percent of H/C ratio used in the minimization process.

As with the A/B undefined petroleum fraction, the same sensitivity analysis was made for the VNAR, KRAR, and LAHR fractions. Figure 7 illustrates that minimum values of the RAD are achieved when H/C parameter lies between 80 and 85% of the experimental value. Therefore, we suggested to set the H/C ratio in this interval in order to obtained better results for heavy undefined petroleum fraction with molecular weights above 500 g/mol.

thumbnail Figure 7

Sensitivity analysis of RAD as a function of H/C parameter in the characterization process for different cuts.

The results of hydrogen solubility prediction in Athabasca Bitumen (AB) at different temperatures and pressures are illustrated in Figures 8a and 8b. Here, the experimental value of the H/C ratio is 1.51, but, according to that previously discussed, it was set equal to 1.208, obtaining an overall RAD of 6.8%. In Table 7 all results are summarized, the overall RAD being equal to 14.12%.

thumbnail Figure 8

Prediction of hydrogen solubility in Athabasca Bitumen a) at 373.15 K, 473.15 K and 573.15 K, b) at 323.15 K, 423.5 K and 523.15 K.

Table 7

RAD on prediction of H2 solubility in undefined hydrocarbon fractions.

Conclusions

The characterization procedure presented by Carreón-Calderón et al. (2012, 2014) was extended to predict hydrogen solubility in heavy undefined petroleum fractions including H2 as another functional group, and the hydrogen-carbon ratio as third constraint in the characterization process. Group-interaction parameters between hydrogen and the main functional groups of the UNIFAC approach were calculated through a set of correlations based on experimental vapor-liquid equilibria data of binary mixtures of hydrogen and defined hydrocarbons. These modifications were enough for prediction purposes with molecular weights below 500 g/mol. For undefined petroleum fraction with molecular weights above 500 g/mol, the hydrogen-carbon ratio was modified to 80% of its experimental value. The inclusion of the hydrogen-carbon ratio was revealed as important property in order to construct the hypothetical structures in better accordance with its molecular nature.

Acknowledgments

The authors express their gratitude to the Mexican Institute of Petroleum and SENER-CONACyT for financial support through the Y.61006 research project.

References

  • Cai H.Y., Shaw J.M., Chung K.H. (2001) Hydrogen solubility measurements in heavy oil and bitumen cuts, Fuel 80, 8, 1055–1063. [CrossRef] [Google Scholar]
  • Carreón-Calderón B., Uribe-Vargas V., Ramírez-Jaramillo E., Ramírez-de-Santiago M. (2012) Thermodynamic characterization of undefined petroleum fractions using group contribution methods, Ind. Eng. Chem. Res. 51, 43, 14188–14198. [CrossRef] [Google Scholar]
  • Carreón-Calderón B., Uribe-Vargas V., Ramírez-de-Santiago M., Ramírez-Jaramillo E. (2014) Thermodynamic characterization of heavy petroleum fluids using group contribution methods, Ind. Eng. Chem. Res. 53, 13, 5598–5607. [CrossRef] [Google Scholar]
  • Chávez L.M., Alonso F., Ancheyta J. (2014) Vapor-liquid equilibrium of hydrogen-hydrocarbon systems and its effects on hydroprocessing reactors, Fuel 138, 156–175. [CrossRef] [Google Scholar]
  • Chen J., Mulgundmath V., Wang N. (2011) Accounting for vapor-liquid equilibrium in the modeling and simulation of a commercial hydrotreating reactor, Ind. Eng. Chem. Res. 50, 1571–1579. [CrossRef] [Google Scholar]
  • Connolly J.F., Kandalic G.A. (1986) Gas solubilities, vapor-liquid equilibria, and partial molal volumes in some hydrogen-hydrocarbon systems, J. Chem. Eng. Data 31, 4, 396–406. [CrossRef] [Google Scholar]
  • Deiters U.K. (2013) Comments on the modeling of hydrogen and hydrogen-containing mixtures with cubic equations of state, Fluid Phase Equilib. 352, 93–96. [CrossRef] [Google Scholar]
  • Elbro H.S., Fredenslund A., Rasmussen P. (1991) Group contribution method for the prediction of liquid densities as a function of temperature for solvents, oligomers, and polymers, Ind. Eng. Chem. Res. 30, 12, 2576–2582. [CrossRef] [Google Scholar]
  • Firoozabadi A. (1988) Reservoir-Fluid Phase Behavior and Volumetric Prediction With Equations of State (includes associated papers 18400 and 18579), J. Petrol. Technol. 40, 397–406. [CrossRef] [Google Scholar]
  • Florusse L.J., Peters C.J., Pàmies J.C., Vega L.F., Meijer H. (2003) Solubility of hydrogen in heavy n-alkanes: experiments and SAFT modeling, AIChE J. 49, 12, 3260–3269. [CrossRef] [Google Scholar]
  • Hook B.D. (1985) Effect of feed volatility on conversion in trickle bed reactors, Ind. Eng. Chem. Fundam. 24, 3, 398–401, Available at: http://www.scopus.com/inward/record.url?eid=2-s2.0-0022110777&partnerID=40&md5=a9e63ff4de18848d030986e035bd2866. [CrossRef] [Google Scholar]
  • Horstmann S., Jabłoniec A., Krafczyk J., Fischer K., Gmehling J. (2005) PSRK group contribution equation of state: comprehensive revision and extension IV, including critical constants and α-function parameters for 1000 components, Fluid Phase Equilib. 227, 2, 157–164. [CrossRef] [Google Scholar]
  • Ihmels E.C., Gmehling J. (2003) Extension and revision of the group contribution method GCVOL for the prediction of pure compound liquid densities, Ind. Eng. Chem. Res. 42, 2, 408–412. [CrossRef] [Google Scholar]
  • Jaubert J.N., Privat R., Mutelet F. (2010) Predicting the phase equilibria of synthetic petroleum fluids with the PPR78 approach, AIChE J. 56, 12, 3225–3235. [CrossRef] [Google Scholar]
  • Ji S., Wang Z., Guo A., Zhou Y., Chen K. (2013) Determination of hydrogen solubility in heavy fractions of crude oils by a modified direct method, J. Chem. Eng. Data 58, 12, 3453–3457. [CrossRef] [Google Scholar]
  • Joback K.G., Reid R.C. (1987) Estimation of pure-component properties from group-contributions, Chem. Eng. Commun. 57, 1-6, 233–243. [CrossRef] [Google Scholar]
  • Lal D., Otto F.D., Mather A.E. (1999) Solubility of hydrogen in Athabasca bitumen, Fuel 78, 12, 1437–1441. [CrossRef] [Google Scholar]
  • Lee B.I., Kesler M.G. (1975) A Generalized Thermodynamic Correlation Based on Three-Parameter Corresponding, AIChE J. 21, 3, 510–527. [Google Scholar]
  • Lin H.M., Kim H., Guo T., Chao K.C. (1985) Equilibrium vaporization of a coal liquid from a Kentucky No. 9 coal, Ind. Eng. Chem. Process Des. Dev. 24, 1049–1055. [CrossRef] [Google Scholar]
  • Marrero J., Gani R. (2001) Group-contribution based estimation of pure component properties, Fluid Phase Equilib. 183-184, 183–208. [CrossRef] [Google Scholar]
  • Michelsen M.L. (1990) A modified Huron-Vidal mixing rule for cubic equations of state, Fluid Phase Equilib. 60, 1-2, 213–219. [CrossRef] [Google Scholar]
  • Nichols W.B., Reamer H.H., Sage B.H. (1957) Volumetric and phase behavior in the hydrogen-n-hexane system, AIChE J. 3, 2, 262–267. [CrossRef] [Google Scholar]
  • Park J., Robinson R.L., Gasem K.A.M. (1996) Solubilities of Hydrogen in Aromatic Hydrocarbons from 323 to 433 K and Pressures to 21.7 MPa, J. Chem. Eng. Data 41, 1, 70–73. [CrossRef] [Google Scholar]
  • Pedersen K.S., Milter J., Sørensen H. (2004) Cubic equations of state applied to HT/HP and highly aromatic fluids, SPE J. 9, 186–192. [CrossRef] [Google Scholar]
  • Pellegrini L.A., Gamba S., Calemma V., Bonomi S. (2008) Modelling of hydrocracking with vapour - liquid equilibrium, Chem. Eng. Sci. 63, 4285–4291. [CrossRef] [Google Scholar]
  • Peng D.-Y., Robinson D.B. (1976) A new two-constant equation of state, Ind. Eng. Chem. Fundam. 15, 1, 59–64. [CrossRef] [Google Scholar]
  • Privat R., Jaubert J.N. (2013) Classification of global fluid-phase equilibrium behaviors in binary systems, Chem. Eng. Res. Des. 91, 10, 1807–1839. [CrossRef] [Google Scholar]
  • Qian J.W., Jaubert J.N., Privat R. (2013) Phase equilibria in hydrogen-containing binary systems modeled with the Peng-Robinson equation of state and temperature-dependent binary interaction parameters calculated through a group-contribution method, J. Supercrit. Fluids 75, 58–71. [CrossRef] [Google Scholar]
  • Saajanlehto M., Uusi-Kyyny P., Alopaeus V. (2014) Hydrogen solubility in heavy oil systems: experiments and modeling, Fuel 137, 393–404. [CrossRef] [Google Scholar]
  • Uribe-Vargas V., Carreón-Calderón B., Ramírez-Jaramillo E., Ramírez-de-Santiago M. (2016) Thermodynamic characterization of undefined petroleum fractions of gas condensate using group contribution, Oil Gas Sci. Technol. - Rev. IFP 71, 1, 5. [CrossRef] [EDP Sciences] [Google Scholar]
  • Xu X., Jaubert J.-N., Privat R., Duchet-Suchaux P., Braña-Mulero F. (2015) Predicting binary-interaction parameters of cubic equations of state for petroleum fluids containing pseudo-components, Ind. Eng. Chem. Res. 54, 10, 2816–2824. [CrossRef] [Google Scholar]
  • Young C.L. (1981) Hydrogen and Deuterium (Solubility Data Series), Pergamon Press, Australia, pp. 314–321. [Google Scholar]

Cite this article as: H. Aguilar-Cisneros, V. Uribe-Vargas, B. Carreón-Calderón, J.M. Domínguez-Esquivel and M. Ramirez-de-Santiago (2017). Hydrogen Solubility in Heavy Undefined Petroleum Fractions Using Group Contributions Methods. Oil Gas Sci. Technol 72, 2.

All Tables

Table 1

Original functional groups.

Table 2

List of 13 binary mixtures used in this study.

Table 3

RAD of hydrogen solubility in defined hydrocarbons.

Table 4

Original group-interaction parameters (αnm, αmn) between functional groups used in this study (Horstmann et al., 2005).

Table 5

Group-interaction parameters (αnm) between functional groups used in this study.

Table 6

Relevant physical properties of the undefined petroleum fractions used in this work.

Table 7

RAD on prediction of H2 solubility in undefined hydrocarbon fractions.

All Figures

thumbnail Figure 1

Group-interaction parameters as function of carbon number a) paraffinic compounds, b) aromatic compounds.

In the text
thumbnail Figure 2

Solubility of hydrogen a) pyrene at 433.2 K, b) hexatetracontane at 402.4 K.

In the text
thumbnail Figure 3

Solubility of hydrogen in CHVGO petroleum fraction at 603.15 K as a function of group-interaction parameters.

In the text
thumbnail Figure 4

Coefficients calculated by characterization process.

In the text
thumbnail Figure 5

Flow chart of the modified characterization process.

In the text
thumbnail Figure 6

Sensitivity analysis of H/C parameter used in the characterization process for AB. a) Pseudostructure arrangements of AB depending on percent of the H/C ratio, b) RAD vs. percent of H/C ratio used in the minimization process.

In the text
thumbnail Figure 7

Sensitivity analysis of RAD as a function of H/C parameter in the characterization process for different cuts.

In the text
thumbnail Figure 8

Prediction of hydrogen solubility in Athabasca Bitumen a) at 373.15 K, 473.15 K and 573.15 K, b) at 323.15 K, 423.5 K and 523.15 K.

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.