Activated desorption at heterogeneous interfaces and long-time kinetics of hydrocarbon recovery from nanoporous media

Hydrocarbon recovery from unconventional reservoirs (shale gas) is debated due to its environmental impact and uncertainties on its predictability. But a lack of scientific knowledge impedes the proposal of reliable alternatives. The requirement of hydrofracking, fast recovery decay and ultra-low permeability—inherent to their nanoporosity—are specificities of these reservoirs, which challenge existing frameworks. Here we use molecular simulation and statistical models to show that recovery is hampered by interfacial effects at the wet kerogen surface. Recovery is shown to be thermally activated with an energy barrier modelled from the interface wetting properties. We build a statistical model of the recovery kinetics with a two-regime decline that is consistent with published data: a short time decay, consistent with Darcy description, followed by a fast algebraic decay resulting from increasingly unreachable energy barriers. Replacing water by CO₂ or propane eliminates the barriers, therefore raising hopes for clean/efficient recovery

Bottom-up model of adsorption and transport in multiscale porous media

We develop a model of transport in multiscale porous media which accounts for adsorption in the different porosity scales. This model employs statistical mechanics to upscale molecular simulation and describe adsorption and transport at larger time and length scales. Using atom-scale simulations, which capture the changes in adsorption and transport with temperature, pressure, pore size, etc., this approach does not assume any adsorption or flow type. Moreover, by relating the local chemical potential μ(r) and density ρ(r), the present model accounts for adsorption effects and possible changes in the confined fluid state upon transport. This model constitutes a bottom-up framework of adsorption and transport in multiscale materials as it (1) describes the adsorption-transport interplay, (2) accounts for the hydrodynamics breakdown at the nm scale, and (3) is multiscale.France. Investissements d'avenir (ICoME2/ANR-11-LABX-0053)France. Investissements d'avenir (A*NUDEX/ANR-11-IDEX-0001-02)Schlumberger FoundationShell Oil Compan

Subcontinuum mass transport of condensed hydrocarbons in nanoporous media

Although hydrocarbon production from unconventional reservoirs, the so-called shale gas, has exploded recently, reliable predictions of resource availability and extraction are missing because conventional tools fail to account for their ultra-low permeability and complexity. Here, we use molecular simulation and statistical mechanics to show that continuum description—Darcy’s law—fails to predict transport in shales nanoporous matrix (kerogen). The non-Darcy behaviour arises from strong adsorption in kerogen and the breakdown of hydrodynamics at the nanoscale, which contradict the assumption of viscous flow. Despite this complexity, all permeances collapse on a master curve with an unexpected dependence on alkane length. We rationalize this non-hydrodynamic behaviour using a molecular description capturing the scaling of permeance with alkane length and density. These results, which stress the need for a change of paradigm from classical descriptions to nanofluidic transport, have implications for shale gas but more generally for transport in nanoporous media.France. Investissements d'avenir (ICoME2 Labex ANR-11-LABX-0053)France. Investissements d'avenir (A*MIDEX ANR-11-IDEX-0001-02)Royal Dutch-Shell GroupSchlumberger Foundatio

Adsorption-based characterization of hierarchical metal–organic frameworks

Nitrogen adsorption at 77 K on metal–organic framework (MOF) is investigated by means of molecular simulations. We consider both regular Cu–BTC crystal and a MOF-based hierarchical porous solid consisting of a mesopore carved out of a Cu–BTC crystal. The t-plot method is applied to these solids by using a non-porous Cu–BTC surface as the reference sample. The values of the mesoporous and external surface areas are determined from the t-plot, and the validity of the method for this type of hierarchical solid is discussed.Institut Carnot (Chimie, Environnement et De´veloppement Durable

Adsorption of carbon dioxide-methane mixtures in porous carbons: effect of surface chemistry

A combined experimental and molecular simulation study of the coadsorption of CO2 and CH4 in porous carbons is reported. We address the effect of surface chemistry by considering a numerical model of disordered porous carbons which has been modified to include heterochemistry (with a chemical composition consistent with that of the experimental sample). We discuss how realistic the numerical sample is by comparing its pore size distribution (PSD), specific surface area, porous volume, and porosity with those for the experimental sample. We also discuss the different criteria used to estimate the latter properties from a geometrical analysis. We demonstrate the ability of the MP method to estimate PSD of porous carbons from nitrogen adsorption isotherms. Both the experimental and simulated coadsorption isotherms resemble those obtained for pure gases (type I in the IUPAC classification). On the other hand, only the porous carbon including the heterogroups allows simulating quantitatively the selectivity of the experimental adsorbent for different carbon dioxide/methane mixtures. This result shows that taking into account the heterochemistry present in porous carbons is crucial to represent correctly adsorption selectivities in such hydrophobic samples. We also show that the adsorbed solution theory describes quantitatively the simulated and experimental coadsorption isotherms without any parameter adjustment

Surface of glassy GeS[subscript 2]: A model based on a first-principles approach

First-principles calculations within the framework of the density functional theory are used to construct realistic models for the surface of glassy GeS[subscript 2](g−GeS[subscript 2]). Both calculations at T = 0 K and at finite temperature (T = 300 K) are considered. This allows for a comparison between the structural and electronic properties of surface and bulk g−GeS[subscript 2]. Although the g−GeS[subscript 2] surface recovers the main tetrahedral structural motif of bulk g−GeS[subscript 2], the number of fourfold coordinated Ge atoms and twofold coordinated S atoms is smaller than in the bulk. On the contrary, the surface system features a larger content of overcoordinated S atoms and threefold coordinated Ge atoms. This effect is more important for the g−GeS[subscript 2] surface relaxed at 0 K. Maximally localized Wannier functions (WF) are used to inspect the nature of the chemical bonds of the structural units present at the g−GeS[subscript 2] surface. We compare the ability of several charge derivation methods to capture the atomic charge variations induced by a coordination change. Our estimate for the charges allows exploiting the first-principles results as a data base to construct a reliable interatomic force field