Linking up pressure, chemical potential and thermal gradients

Petroleum reservoirs are remarkable illustrations of the impact of a thermal gradient on fluid pressure and composition. This topic has been extensively studied during the last decades to build tools that are required by reservoir engineers to populate their models. However, one can get only a very limited number of representative samples from a given reservoir and assessing connectivity between all sampling points is often a key issue. In some extreme cases, the whole reservoir fluid properties must be derived from a single point to define the field development plan. To do so, available models are usually not satisfactory as they need too many parameters and so cannot be considered as predictive tools. We propose in this work a comprehensive approach based on the irreversible thermodynamics principles to derive the relationships between pressure, chemical potentials and thermal gradients in porous media. It appears that there is no need for additional assumptions, it is just a matter of a making the right choices along theoretical developments. One of the most important steps is to express the full pressure gradient. As a final result, we obtain the chemical potential gradients for all components of the mixtures that can be easily translated in term of compositions through Equation of State modelling. The most important features of the final expressions are: i) the species relative separation in a thermal field is sensitive to the relative diffusion coefficients at stationary state. In porous media, the separation is sensitive to the permeability when the overall mobility is similar to diffusive mobility; ii) the magnitude of the separation depends on the residual entropy of the species; iii) the separation is not simply balanced by the average residual entropy. The balance is modified by the relative diffusion mobility of the components; iv) in low permeability porous media, the thermal gradient induces a pressure gradient proportional to the fluid residual entropy. As a validation, the proposed approach has been applied on a reservoir fluid subjected to a geothermal gradient and compared with non-equilibrium molecular dynamics simulation results at the stationary state.