Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Dec 4;115(49):12365-12370.
doi: 10.1073/pnas.1808402115. Epub 2018 Nov 15.

Mesoscale structure, mechanics, and transport properties of source rocks' organic pore networks

Affiliations

Mesoscale structure, mechanics, and transport properties of source rocks' organic pore networks

Jeremie Berthonneau et al. Proc Natl Acad Sci U S A. .

Abstract

Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores' topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.

Keywords: electron tomography; fluid transport; mechanics; mesoscale; porous media.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

Fig. 1.
Fig. 1.
Morphological characterization of the organic mesoporous networks imaged in electron tomography. Aperture maps of the organic pore networks of (A) LEF (φmeso = 9.1%), (B) HAY (φmeso = 19.3%), and (C) MAR (φmeso = 16.9%). (D) Pore-size distributions from aperture maps calculation. (E) In-pore CLDs. (F) Specific surface areas (As) calculation from Eq. 1 as a function of the cutoff distance (rc) of the in-pore C.L.D.
Fig. 2.
Fig. 2.
Elastic simulations of the two porous mesostructures originating from LEF and MAR under two types of agent: (a) imposed strain boundary, and (b) adsorbed fluid with free strain boundary. The averaged strain of the simulation box is ϵii = 0.01 for all cases. The size is 700 × 700 pixels with a resolution of 0.42 nm for LEF and 0.35 nm for MAR. Mesopore phase is in white. (Top) Elastic energy density fields (106 J/m3). (Bottom) Distribution of the elastic energy density (106 J/m3) for each simulation.
Fig. 3.
Fig. 3.
Homogenized fluid transport properties of the mesoscale lattice models. (A) Homogenized diffusion coefficient of the MAR reconstruction scaled by Dnano as function of the diffusion contrast rD=Dmeso/Dnano for different concentration ratios rρ=ρmeso/ρnano indicated by colors. The corresponding curves stand for the effective medium theory (Eq. 2) where the value of γ(rρ)2 is D¯/Dnano when Dmeso=Dnano (i.e., rD=1). (B) Concentration-dependent obstruction factor γ(rρ)2 as function of the probability p(rρ) for a random walker to cross the interface from the nano- to the mesoporous phase (red: LEF, green: HAY, blue: MAR). The corresponding curves display the analytical model (Eq. 4) where γ02 is simply taken as γ02=γ(p=0)2. (C) Homogenized diffusion coefficient of methane for LEF (red) and MAR (blue) as function of the bulk fluid pressure. The dashed curves correspond to the diffusion in the nanoporous phase only, while the solid lines represent the homogenized model with mesoporosity. Standard geological conditions (T = 400 K and lithostatic pressure of 25 MPa) have been considered.

References

    1. Kerr RA. Energy. Natural gas from shale bursts onto the scene. Science. 2010;328:1624–1626. - PubMed
    1. Cueto-Felgueroso L, Juanes R. Forecasting long-term gas production from shale. Proc Natl Acad Sci USA. 2013;110:19660–19661. - PMC - PubMed
    1. Valenza JJ, Drenzek N, Marques F, Pagels M, Mastarlerz M. Geochemical controls on shale microstructure. Geology. 2013;41:611–614.
    1. Monteiro PJM, Rycroft CH, Barenblatt GI. A mathematical model of fluid and gas flow in nanoporous media. Proc Natl Acad Sci USA. 2012;109:20309–20313. - PMC - PubMed
    1. Silin D, Kneafsey T. Shale gas: Nanometer-scale observations and well modelling. J Can Pet Technol. 2012;51:464–475.

Publication types

LinkOut - more resources