Abstract Design and modeling of solvent-based heavy oil recovery processes require the knowledge of solvent diffusivity in heavy oil. The solvent diffusivity in heavy oil at constant temperature can be measured by using the so-called pressure decay method. In this method, the solvent diffusivity is determined by monitoring the decaying pressure of solvent phase while it is in contact with heavy oil in a closed high-pressure cell. This paper examines the effects of three different boundary conditions at the solvent–heavy oil interface on the determined diffusion coefficient. The so-called equilibrium, quasi-equilibrium and non-equilibrium boundary conditions are applied to represent three rather different interface mass transfer models. An analytical solution to the bulk diffusion equation subjected to each boundary condition is presented to determine the concentration distribution of solvent in the heavy oil. With the determined solvent concentration distribution in the heavy oil, pressure in the solvent phase is then predicted by utilizing the solvent mass conservation and the equation of state for a real gas. Therefore, an average solvent diffusivity, i.e., independent of solvent concentration, is determined by matching the numerically predicted pressures with the experimentally measured data. Experimental data for two different systems, CO2-heavy oil and CH4-heavy oil, are analyzed and the diffusion coefficients are compared with the previously published values. Based on the minimum discrepancy between the predicted and measured pressure data, it is found that the non-equilibrium BC is more applicable to CO2-heavy oil system and that the equilibrium BC is most suitable for CH4-heavy oil system. Thus it is concluded that the obtained solvent diffusivity is sensitive to the boundary condition applied at the solvent–heavy oil interface. Hence, a proper boundary condition for each solvent–heavy oil system should be used to determine the solvent diffusivity in heavy oil.
[1]
Byung-Ik Lee,et al.
A generalized thermodynamic correlation based on three‐parameter corresponding states
,
1975
.
[2]
Mohammad R. Riazi,et al.
A new method for experimental measurement of diffusion coefficients in reservoir fluids
,
1996
.
[3]
Simant R. Upreti,et al.
Experimental measurement of gas diffusivity in bitumen: Results for carbon dioxide
,
2000
.
[4]
J. C. Jaeger,et al.
Conduction of Heat in Solids
,
1952
.
[5]
Stanley M. Walas,et al.
Modeling with differential equations in chemical engineering
,
1991
.
[6]
John Crank,et al.
The Mathematics Of Diffusion
,
1956
.
[7]
C. Whitson,et al.
Estimating Diffusion Coefficients of Dense Fluids
,
1993
.
[8]
M. Rasmussen,et al.
Analysis and Interpretation of Gas Diffusion in Quiescent Reservoir, Drilling, and Completion Fluids: Equilibrium vs. Non-equilibrium Models
,
2003
.
[9]
Simant R. Upreti,et al.
Diffusivity of CO2, CH4, C2H6 and N2 in athabasca bitumen
,
2002
.
[10]
W. V. Pinczewski,et al.
The Role of Molecular Diffusion Processes in Tertiary CO2 Flooding
,
1987
.
[11]
L. Hepler,et al.
Comments on "Solubility of carbon dioxide in tar sand bitumens: experimental determination and modeling"
,
1992
.
[12]
Donald L. Katz,et al.
Density of Natural Gases
,
1942
.
[13]
T. A. Renner,et al.
Measurement and Correlation of Diffusion Coefficients for CO2 and Rich-Gas Applications
,
1988
.
[14]
Phillip M. Sigmund,et al.
Prediction of Molecular Diffusion At Reservoir Conditions. Part 1- Measurement And Prediction of Binary Dense Gas Diffusion Coefficients
,
1976
.
[15]
Y. P. Zhang,et al.
Measurement of Gas Diffusivity In Heavy Oils
,
1998
.
[16]
L. Dake.
Fundamentals of Reservoir Engineering
,
1983
.