Thermal conductivity model function of porosity: review and fitting using experimental data
IFP Energies nouvelles, 1–4 avenue du Bois-Préau, 92852 Rueil-Malmaison, France
* Corresponding author: Christophe.firstname.lastname@example.org
Accepted: 9 August 2021
Thermal conductivity of porous rocks depends on a large variety of proper to rock parameters as well as external influences. Thus, it can generate difficulties in determining accurate thermal behavior of the rock. The rock parameters which influence the thermal conductivity are principally the porosity, the microstructure  and the mineral composition. However, these parameters, in turn, can be impacted by external influences such as temperature and pressure.
An accurate determining of the thermal conductivity is crucial in oil and gas engineering or in geothermal application. For example, during thermal EOR or geothermal application, the porosity and/or the microstructure of the sedimentary rocks can vary due to the increase of temperature and pressure, and this modification must be quantified to be accounted for the thermal behavior of rocks.
Many efforts have been done to estimate the thermal conductivity of sedimentary rocks in parallel to the experimental methods for its determination. These estimations have always been the subject of intensive studies, and a lot of data [2, 3] are obtained as well as models and methodologies to characterize the thermal conductivity of rocks [4–7]. Moreover, this type of estimation is well-known by other research communities. Indeed, we find the same formal analogy between Fourier, Ohm’s law, Darcy’s laws and thermal conductivity. For example, considering Darcy’s laws, the same problem is well-known and termed “upscaling” [8, 9] and consists in computation of the effective permeability considering a heterogeneous rock. Classically, the upscaling process can be related to percolation theory , which describes connectivity of objects within for example, a porous structure. We can also determine effects of this connectivity on macroscale properties such as thermal conductivity . In particular, the fundamental contributions of Torquato who proposed strategies via rigorous microstructure-property relations [1, 12].
Finally, many technics are based on a porosity dependence and a link between the conceptual thermal conductivity of the non-porous rock, λR, and the thermal conductivity of the fluid saturated the porous rock λf. These technics are simple to implement especially when there is no precise information about the microstructure.
In order to predict accurately the thermal efficiency of the geothermal installation or the oil recovery of a thermal EOR process, such as, for example, Steam Assisted Gravity Drainage (SAGD), very often the engineers invoke numerical simulations. Numerous reservoir simulators [13–16] allow to estimate the thermal conductivity as function of porosity, but these solutions are often based on a mixing laws which are quite simplistic models.
The purpose of this paper is to propose a better methodology to predict a thermal conductivity of reservoir rocks depending on the porosity for a reservoir simulator which satisfy the following conditions
easy to implement in a software, precise and accurate,
it should have finite values to avoid numerical problems such as non-convergence,
have ideally correct limits: and ,
have a limited number of parameters in order to minimize a required number of laboratory experiments.
First for all, we perform a review of many means or methods available in the literature, starting by the classical mixing laws and finishing by empirical and theoretical methods. For each method, the limits for φ → 1 and φ → 0 are examined and the corresponding parameters are detailed.
Then, we compare these methods to experimental data. Two experimental datasets are considered, water and air-saturated rock of the same properties for the both cases (sandstone). The purpose is to match these two datasets using the same model. Indeed, in reservoir simulation (oil and gas or geothermal), various fluids may flow in the same rock (oil, gas, CO2, water, steam). Thus, it is crucial to propose a model which yields correct thermal conductivity independently on the fluids flowing in the pore volume.
The heat transfer in the reservoir is mainly dominated by conduction. Thermal conductivity is the property that quantifies the ability of a geological formation for heat transfer. There are three types of methods to define this conductivity of multicomponent systems, mixing laws, empirical models and theoretical models. The following denotations are used
λR, the rock thermal conductivity,
λf, the fluid thermal conductivity,
φ, the rock porosity.
The mixing laws are very common and are also used for Ohm’s law or Darcy’s laws effective parameters determination [8, 9]. But they don’t consider the structural characteristics of rocks and, therefore, their applications can be constrained.
First, consider the weighted arithmetic and harmonic means, reviewed by Zimmerman  which are also known as parallel, or linear, and series models, respectively. They correspond to the classical Wiener’s bound used to compute effective permittivity.
The weighted arithmetic or so-called parallel model is written as
The weighted harmonic or so-called series model, is given by
It can be noticed that the parallel model yields the highest values of thermal conductivity and the series model yields the lowest ones . It was also shown by Wiener, that these two models represent upper and lower bounds on the effective conductivity .
Then, the mixing law is reviewed by Beck 18] as the weighted geometric mean. It has no physical background but seems to work better than the parallel and series means; it is written as
There is another mixing law, proposed by Hashin and Shtrikman , which yields the bounds (HS bounds) always tighter than the Wiener bounds
In , the mean of both bounds is often used as best approximation of rock bulk thermal conductivity:
However, it can be easily demonstrated that this formulation corresponds to the upper Hashin–Shtrikman bound λeff, max(φ). Therefore, in this paper, only HS bounds model is referred.
To conclude, all presented mixing laws are tend to expected values at the porosity limits, i.e., and . Moreover, no additional parameter is required to set these models. Notice that Hashin and Shtrikman bounds are classical and appear also in the computation of effective permeability of heterogeneous porous media .
This section addresses a review of the existing empirical models to determine the thermal conductivity of the porous rock. These models determine the thermal conductivity through the application of a regression technics. Generally, they are limited to a particular type of rocks. The purpose is to demonstrate the link in-between different models, empirical or mixing laws, as well as verify that the conductivity values at the porosity limits satisfy the expected values, λf for φ → 1 and λR for φ → 0.
First, Asaad  proposes a thermal conductivity function obtained by the mean which is very similar to weighted geometric model
where c is an empirical exponent and its value can be fitted using experimental data. When c = 1, this model becomes identical to the weighted geometric mean model. Moreover, ∀c > 0 only one of the conditions at the porosity limits is satisfied
with n (> 0) is an empirical exponent which needs to be fitted using experimental data. Notice that if n = −1, one obtains the weighted arithmetic mean (parallel model). Moreover ∀n > 0,
Zeb and Maqsood  used exponential decay trial to predict thermal conductivity of consolidated porous media at room temperature and normal pressure,
where z is an empirical exponent which can be fitted using experimental data of thermal conductivity λeff and corresponding values of φ and λR as proposed by Zeb and Maqsood . This means that the thermal conductivity of the non-porous rock is not given and is not a constant. Notice that this conductivity, λR, is very difficult to obtain because the existence of a non-porous rock with the same mineral composition as a porous one is almost impossible. However, for a reservoir simulator, only correct limits and a limited number of undefined parameters are identified as most important.
Here, is directly verified. For φ = 1, the required limit condition is satisfied only if , which yields the weighted geometric mean.
A model proposed by Zeb et al.  for prediction of thermal conductivity of porous consolidated igneous rocks is written as
where m is an empirical coefficient and should be fitted using experimental data of thermal conductivity λeff and corresponding values of φ and λf. It means that the thermal conductivity of the fluids is not given and is not a constant.
For the zero porosity, the model verifies , as required. However, for φ = 1 the condition is fulfilled only if , i.e. when the model degenerates toward the weighted harmonic mean. This could represent a problem for a reservoir simulator since λR is not easy to determine contrary to λf which can be measured.
where and with . The term is to account for a high thermal conductivity ratio , ψ is an empirical coefficient which is fitted using experimental data. According to Veerendra and Chaudhary, there is a corrective term, , which may be either added or subtracted from the precedent equation. For the porosity limits, φ = 0 and φ = 1, the required conditions are verified if = 0.
This section addresses a review of the existing theoretical models to determine the thermal conductivity of the porous rock. These models are based on physical consideration such as pore structures.
The first theoretical model to predict the thermal conductivity of fluid saturated rocks is the Krupiczka model . It is a semi-empirical equation derived from the numerical calculation of heat transfer through a bundle of cylinders for the effective conductivity of a packed bed of spheres, therefore, it is indexed in theoretical models,
were the coefficients , and are obtained for this particular case of packed bed of spheres . It should be noted that the required condition at the zero porosity limit, , cannot be satisfied for this model.
The last considered here model is based on the concept of the porous rocks represented as a composite material with the spheroidal non-connected pores. This formulation represents some mathematical difficulties due to arbitrary aspect ratio of the spheroidal inclusions. The solution has been proposed by Fricke  for the electrical conductivity  and can be successfully applied to the thermal conductivity. Thus, for small porosity, Fricke showed that the effective conductivity can be given by
where and M is a factor that depends on the pore aspect ratio . For the two regimes of oblate (α < 1) and prolate (α > 1) spheroids, this factor is given by
There are three limiting cases which are distinguished by Zimmerman  as thin cracks for α → 0, spherical pores for α = 0 and needle-like pores for α → ∞, with
To consider rock with any value of porosity and, therefore, to consider realistic rocks we must evoque the “effective medium theory” for taking into account the interactions of neighboring pores as explained by Zimmerman . Then, applying the effective medium theory proposed by Maxwell , the actual Fricke model  can be written as
where r and β are defined as previously. This model is tested in the next section.
The limit for the zero porosity verifies and for φ = 1, the model tends to λf as required. Moreover, Fricke’s equation always satisfies the Hashin–Shtrikman bound.
To test all these models, the published experimental data is used from samples of sandstone collected in five wells  in Perth Basin, Western Australia and Soultz-sous-Forets Basin, astern France using  to test thermal conductivity calculation from P-wave velocity. Two sets of data are available
Measured thermal conductivity for dry sample, λf = λair,
Measured thermal conductivity for saturated sample, λf = λwater.
These values are used together with the laboratory data and summarized in Table 1.
Experimental data and weights; and denote the effective conductivities for dry and saturated media, respectively.
The purpose is to match different models to the experimental data. Independently on the fluid (air or water), the value of the thermal conductivity of the rock λR is assumed to be a unique constant within the same model. Thus, the both datasets, air and water, are considered together to search λR. Therefore, the matching methodology is constrained as follows.
First, an Objective Function (OF) based on least squares errors is built
where N is a number of experimental entries. Each term is weighted by ωi which is assumed to be equal to one for all the experimental data. To impose the convergence at the two known points λAir and λwater for φ = 1, their weight is set to ω = 0.1 (Tab. 1).
In order to get optimal model parameters to approach the data, the OF is minimized using the Powell’s dog leg minimization method.
Then, the standard deviation is used to estimate the veracity of the results
In order to bound a maximal thermal conductivity, the study of Orlander et al.  for several sandstone outcrops originate from Fontainebleau (France), Castlegate (USA), Bentheim (Germany), Obernkirchen (Germany) and Berea (USA), is used. A maximal thermal conductivity for sandstone obtained by Orlander et al.  is near 8 Wm−1 K−1 for φ → 0. However, to release the parameter restrictions, the bounds for λR are set to
The mixing laws are given as simple means, which depend only on the thermal rock conductivity λR. As mentioned before, the measure of λR is very difficult, thus, for each model, λR is obtained in order to match as close as possible the experimental data.
In Figure 1, the resulting mixing laws are compared to the experimental data. The precision of each model is characterized by the corresponding OF and the standard deviation errors (SDEs) given in Table 2 together with resulting thermal rock conductivity.
The Objective Function (OF), Standard Deviation Error (SDE) and resulting from the minimization of thermal rock conductivity λR.
The worst fit is obtained with the weighted harmonic mean. Moreover, the optimal value of λR corresponds to the set maximal bound which implies a bad convergence of the minimization.
Then, the Hashin and Shtrikman maximal and minimal bounds and mean models also yield high OF and SDE. The corresponding fitting curve for minimal bound matches with the data for the saturated sample, while for the non-saturated sample the minimal bound is found extensively low compared to the data. Moreover, the obtained optimal value of λR is higher than value proposed by Orlander et al. . The maximal bound match with non-saturated sample data. The Hashin and Shtrikman mean model yields small SDE, however the objective function keeps high.
The weighted arithmetic mean has also a high error compared to the experimental data. This model is a linear correlation and it is not adapted for the used set of experimental results which has a non-linear trend. It can be noted that the SDE is very close to those of the Hashin and Shtrikman mean model.
Finally, the best fit is obtained with the weighted geometric mean model. It yields the smallest objective function and the standard deviation for the thermal conductivity of rock λR = 6.74 . Moreover, this model yields the required bounds for porosity limits, and .
The empirical and theoretical models are more complex as the mixing laws. Additionally to the thermal conductivity of rock, these models have other fitting parameters. Similarly to the mixing laws, these parameters can be determined through minimizing the objective function given before or defined by the pore structure. The minimization parameters should satisfy the following conditions.
The Asaad model  tends to the weighted geometrical mean when c = 1. Moreover, it was demonstrated that the weighted geometrical mean yields the best fit with the experimental data. Therefore, the empirical exponent c of the Asaad model is assumed to be closed to 1.
Zeb and Maqsood model  found the values of the empirical exponent z between 0.1 and 5.
Zeb et al.  found the values of the model parameter m in the range [0.1, 5].
The parameter ψ of the Veerendra and Chaudhary model is to account only for a high thermal conductivity ratio . Moreover, in order to maintain the condition on the thermal conductivity at porosity limits, and , it should be close to zero.
For the Krupiczka model, we can think that the parameters values must be closed to those found by the author, a1 = 0.28, a2 = −0.757 and B = −0.057.
The parameter M of Fricke model is a factor which depends on the aspect ratio of the pore. For spherical pore, M = 0. Consider that for the studied sample of sandstone the pores are nearly spherical. Thus, the parameter M should be found close to zero.
Table 3 summarizes the resulting optimal parameters as well as the corresponding objective function values and the SDEs. For all the models, the minimization parameters are obtained in agreement with the required conditions. The fits obtained for each empirical model are compared to the experimental data in Figure 2.
The Objective Function (OF), Standard Deviation Error (SDE), resulting from the minimization of thermal rock conductivity λR and optimal parameters.
The worst results are obtained with the Zeb and Maqsood model and followed by Zeb, Gurmani, Ali and Maqsood and Veerendra and Chaudhary models. Similarly to the weighted harmonic mean, the Zeb, Gurmani, Ali and Maqsood model yields the thermal rock conductivity equal to the imposed upper bound, λR = 16, which means the poor convergence of the solution. Then, the Sugawara and Yoshizawa and Krupiczka models don’t yield a good approximation as well. Moreover, at zero porosity limit, the Krupiczka model has not a finite solution, therefore it is not of interest for the present study. Indeed, a not finite solution can bring a divergence of the reservoir simulator. For these five models, it seems impossible to match the data using common minimization parameters for saturated (water) and non-saturated (air) cases.
Finally, the Asaad and Fricke models give the best match with the experimental results. Applying the optimal value of the Asaad model parameter c = 0.99992, this model degenerates to the weighted geometric mean. The errors obtained with Fricke model are insignificantly higher than those of Asaad and geometric mean models. However, it is the only model which has a strong physical background compared to the others mostly pure mathematical models.
Consider four empirical models which are revealed unsatisfactory for the data approximation using unique parameters independent on the fluid, the Sugawara and Yoshizawa, Zeb and Maqsood, Zeb, Gurmani, Ali and Maqsood and Veerendra and Chaudhary models. The Krupiczka model is not considered here since it doesn’t yield the finite thermal conductivity values at zero porosity limit. Suppose that the fitting parameters are fluid dependent and λR keeps a unique constant for water and air curves. Thus, for the corresponding objective function there is one minimization parameter more than before. The results are summarized in Table 4.
The Objective Function (OF), Standard Deviation Error (SDE), resulting from the minimization of thermal rock conductivity λR and optimal air and water parameters.
As previously, the Zeb, Gurmani, Ali and Maqsood model is not correctly converged since the resulting value of λR is found on its upper limit, λR = 16. The Veerendra and Chaudhary models are found the same for fluid dependent parameters as for unique ones (Fig. 3).
Empirical models obtained with fluid dependent minimization parameters.
Considering the fluid dependent parameters for Sugawara and Yoshizawa and Zeb and Maqsood, allow to obtain the better approximation of the data. The corresponding objective functions and the SDEs are as small as for the geometric mean (Tab. 4).
The dependency on type of fluid increases the number of the model parameters. Moreover, in a reservoir simulation, the algorithm must choose the good fitting parameter function of saturated fluids in each cell. It becomes especially problematic when the rock is saturated with many fluids (steam and water, for example). The fitting parameter is available only for one fluid. Therefore, when the cell is saturated of several fluids, the correlation is not available anymore. Thus, the choice of the fluid dependent parameters makes sense only if the rock is saturated of one fluid and don’t evaluate with time (change porosity, saturation, etc...).
The analysis of various models is summarized in Table 5.
Summary of the mixing laws and empirical/theoretical models analysis.
Combining the model accuracy and required for the numerical simulator conditions, the Fricke and Asaad models reveal to be the most appropriate. Both of them yield a good precision (low OF and SDE), they have finite bounds for porosity limits (φ = 0, 1) and require only one additional to λR fitting parameter.
For both models, at zero porosity, the condition is always satisfied, and for porosity φ = 1, the required condition is also fulfilled.
It should be noted that the Asaad model is comparable to the geometric mean, but an additional fitting parameter makes the Asaad model more flexible which may be interesting if dealing with another rock type, for example.
Finally, it should be kept in the mind that the Fricke model has a strong physical background compared to other purely mathematical models. In view of the results, the Fricke model seems to be the most suitable: correct limits, only one parameter, very good matching and strong physical background.
- Torquato S., Labo F. (1986) Effective properties of two-phase disordered composite media: II. Evaluation of bounds on the conductivity and bulk modulus of dispersions of impenetrable spheres, Phys. Rev. B 33, 9, 6428. [CrossRef] [Google Scholar]
- El Sayed N.A., El Sayed A.M.A. (2019) Thermal conductivity calculation from P-wave velocity and porosity assessment for sandstone reservoir rock, Geothermics 82, 91–96. [CrossRef] [Google Scholar]
- Esteban L., Pimienta L., Sarout J., Piane C.D., Haffen S., Géraud Y., Timms N.E. (2015) Study cases of thermal conductivity prediction from P-wave velocity and porosity, Geothermics 53, 255–269. [CrossRef] [Google Scholar]
- Zeb A., Maqsood A. (2007) Modeling of the effective thermal conductivity of consolidated porous media with different saturants: a test case of Gabbro rocks, Int. J. Thermophys. 28, 4, 1317–1386. [Google Scholar]
- Zimmerman R.W. (1989) Thermal conductivity of fluid-saturated rocks, J. Pet. Sci. Eng. 3, 219–227. [CrossRef] [Google Scholar]
- Anand J., Somerton W.H., Gomaa E. (1973) Predicting thermal conductivities of formations from other known properties, SPE J 13, 5, 267–273. [Google Scholar]
- Sugawara A., Yoshizawa Y. (1962) An experimental investigation on the thermal conductivity of consolidated porous materials, J. Appl. Phys. 33, 10, 3135–3138. [CrossRef] [Google Scholar]
- Renard Ph., de Marsily G. (1997) Calculating equivalent permeability: a review, Adv. Water Resour. 20, 5–6, 253–278. [CrossRef] [Google Scholar]
- Preux C., Le Ravalec M., Enchéry G. (2016) Selecting an appropriate upscaled reservoir model based on connectivity analysis, Oil Gas Sci. Technol. – Rev. IFP Energies nouvelles 71, 5, 60. [Google Scholar]
- Sahimi M. (2009) Applications of Percolation Theory, Taylor & Francis [Google Scholar]
- Chatterjee A., Verma R., Umashankar H.P., Kasthurirengan S., Shivaprakash N.C., Behera U. (2019) Heat conduction model based on percolation theory for thermalconductivity of composites with high volume fraction offiller in base matrix, Int. J. Therm. Sci. 136, 389–395. [CrossRef] [Google Scholar]
- Torquato S. (2020) Predicting transport characteristics of hyperuniform porous media via rigorous microstructure-property relations, Adv. Water Resour. 140, 103565. [CrossRef] [Google Scholar]
- IFPEN (2020) PumaFlow: Technical Manual. [Google Scholar]
- Schlumberger (2016) Eclipse: reference Manual. [Google Scholar]
- Schlumberger (2015) Intersect Version 2015.1 Beta. [Google Scholar]
- CMG (2015) Stars User Guide – Advanced processes and thermal reservoir simulator. [Google Scholar]
- Böttcher C.J.F. (1952) Theory of Electric Polarization, Elsevier, Amsterdam. [Google Scholar]
- Beck A.E. (1976) An improved method of computing the thermal conductivity of fluid-filled sedimentary rocks, Geophysics 41, 133–144. [Google Scholar]
- Hashin Z., Shtrikman H. (1962) A variational approach to the theory of the effective magnetic permeability of multiphasic materials, J. Appl. Phys. 33, 3125–3131. [CrossRef] [Google Scholar]
- Fuchs S., Schütz F., Förster H.-J., Förster A. (2013) Evaluation of common mixing models for calculating bulk thermal conductivity of sedimentary rocks: Correction charts and new conversion equations, Geothermics 47, 40–52. [Google Scholar]
- Zeb A., Gurmani S.F., Ali Z., Maqsood A. (2006) Simultaneous measurement of thermal conductivity, thermal diffusivity and prediction of effective thermal conductivity of porous consolidated igneous rocks at room temperature, J. Phys. D Appl. Phys. 39, 17, 3876–3881. [Google Scholar]
- Colecchio I., Boschan A., Otero A.D., Noetinger B. (2020) On the multiscale characterization of effective hydraulic conductivity in random heterogeneous media: a historical survey and some new perspectives, Adv. Water Resour. 140, 103594. [Google Scholar]
- Asaad Y. (1955) A study of the thermal conductivity of fluid bearing porous rocks, Ph.D. dissertation, Ed. Univ. California, Berkeley, CA. [Google Scholar]
- Sugawara A., Yoshizawa Y. (1961) An Investigation on the thermal conductivity of porous materials and its application to porous rook, Australian J. Phys. 14, 4, 469–480. [Google Scholar]
- Veerendra K., Chaudhary D.R. (1980) Ind. J. Pure Appl. Phys. 18, 984. [Google Scholar]
- Krupiczka R. (1967) Analysis of thermal conductivity in granular materials, Int. Chem. Eng. 7, 122–144. [Google Scholar]
- Fricke H. (1924) A mathematical treatment of the electric conductivity and capacity of disperse systems, Phys. Rev. 24, 575–587. [Google Scholar]
- Maxwell J.C. (1892) A treatise on electricity and magnetism, Vol. 1, 3rd ed. Clarendon Press, Ed. Oxford. [Google Scholar]
- Clauser C., Huenges E. (1995) Thermal conductivity of rocks and minerals in rock physics and phase relations: a handbook of physical constants, Vol. 3, T.J. Ahrens (ed.), American Geophysical Union, Washington, USA. [Google Scholar]
- Orlander T., Adamopoulou E., Jerver Asmussen J., Marczyski A.A., Milsch H., Pasquinelli L., Fabricius I.L. (2018) Thermal conductivity of sandstones from Biot’s coefficient, Geophysics 83, 5, 173–185. [Google Scholar]
© C. Preux & I. Malinouskaya, published by IFP Energies nouvelles, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.