I12-JEEP: Joint Engineering, Environmental and Processing
|
Diamond Proposal Number(s):
[30413]
Abstract: The migration and deposition of fine particles in porous materials is critical in industries such as energy, pharmaceuticals, and environmental engineering. Using 3D time-lapse synchrotron X-ray imaging, we observe fine particles invading porous media, analyzing the effects of pore size and heterogeneity at both pore and macro scales. Glass beads model homogeneous and heterogeneous conditions, revealing a sequence of deposition processes: surface attachment, throat bridging, blocking, pore filling, compaction, and migration. A critical throat-to-particle size ratio of 1.7 governs deposition behavior. At the macro-scale, heterogeneities like beddings and flow pathways influence fines migration and deposition. Based on dynamic 3D imaging, we propose a mechanism for fines behavior in heterogeneous porous media. These findings enhance understanding of fines migration, offering a predictive framework for managing formation damage and optimizing filter cake design in drilling and clean energy applications.
|
Mar 2025
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[30455]
Open Access
Abstract: Permeability is a key rock property important for scientific applications that require simulation of fluid flow. Although permeability is determined using core flooding experiments, recent advancements in micro-CT imaging and pore scale fluid flow simulations have made it possible to constrain permeability honoring pore scale rock structure. Previous studies have reported that complex association of pores and solid grains often results in preferential flow paths which influence the resulting velocity field and, hence, the upscaled permeability value. Additionally, the pore structure may change due to geochemical processes such as microbial growth, mineral precipitation and dissolution. This could result in a flow field which dynamically evolves spatially and temporally. It would require numerous experiments or full physics simulations to determine the resultant upscaled Darcy permeability for such dynamically changing systems. This study presents a graph theory-based approach to upscale permeability from pore-to-Darcy scale for changing pore structure. The method involves transforming a given micro-CT rock image to a graph network map followed by the identification of the least resistance path between the inlet and the outlet faces using Dijkstra's algorithm where resistance is quantified as a function of pore sizes. The least resistance path is equivalent to the path of lowest resistance within the domain. The method was tested on micro-CT images of the samples of Sherwood Sandstone, Ketton Limestone and Estaillades Limestone. The three micro-CT images were used to generate 30 synthetic scenarios for geochemically induced pore structure changes covering a range of pore and solid grain growth. The least resistance value obtained from Dijkstra's algorithm was observed to correlate with upscaled permeability value determined from full physics simulations, while improving computational efficiency by a factor of 250. This provides confidence in using graph theory method as a proxy for full physics simulations for determining effective permeability for samples with changing pore structure.
|
Sep 2024
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[28557]
Open Access
Abstract: Swelling of shale in response to interaction with water is an important consideration within subsurface energy systems. In the case of waste disposal, swelling can provide important barriers around the waste and enhance the sealing ability of rocks. For shale gas exploration purpose, however, swelling may cause wellbore instability. Therefore, a careful study of shale swelling is critical for subsurface energy related applications. Here, the swelling effects of shale were imaged at nanoscale using an advanced synchrotron Transmission X-ray Microscopy (TXM) imaging technique for the first time, with a spatial resolution down to 40.9 nm. Organic matter and clays within the analysed sample were observed to display large swelling effects which resulted in a 50% reduction in porosity. Strain maps generated using Digital Volume Correlation (DVC) show deformation and significant strain were mostly localized to between the contact boundaries of sharp brittle minerals and softer organic matter and clays. This is the first study, to our knowledge, to directly image the swelling deformation of shale at the tens of nanometer scale and provide local information on the strain evolution.
|
Jul 2023
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[15506]
Open Access
Abstract: Determination of mineral texture and diagenetic features in mudstones is crucial to reveal the history of their pore systems and provides key information to predict their future sealing ability, reactivity and storage capacity for sequestered CO2, hydrogen storage or nuclear waste disposal. To understand the spatial transport and storage of fluids, it is necessary to map the distribution of minerals and fractures in three dimensions (3D). This study proposes a novel, multi-scale three-dimensional (3D) imaging method, i.e., a combination of synchrotron- sourced micro- x-ray tomography and lab- sourced nano-tomography, to investigate the sedimentology and diagenetic features of the Bowland Shale, one of the most volumetrically important mudstone-dominated systems in the UK. Diagenetic minerals have been identified and characterised, including pyrite, calcite, kaolinite, illite, chlorite, dolomite, ankerite and authigenic quartz (micro-sized quartz and quartz overgrowths). Multi-scale 3D images provide detailed information about dolomite-ankerite zonation and carbonate dissolution pores. These features cannot be observed or quantified by conventional 2D methods, and they have not been reported in this subject area before. Using these results, potential reactions during carbon storage and other subsurface storage applications are predicted.
|
Apr 2023
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[19904]
Open Access
Abstract: The development of pore and fracture networks at the nano-scale as a response to heating can reveal coupled physical relationships relevant to several energy applications. A combination of time-lapse 3D imaging and finite-element modelling (FEM) was performed on two typical thermally immature shale samples, Kimmeridge Clay and Akrabou shale, to investigate thermal response at the nm-scale for the first time. Samples were imaged using Transmission X-ray Microscopy (TXM) with a voxel resolution of 34 nm at the I13–2 beamline at Diamond Light source, UK. Images were taken after heating to temperatures of 20 °C, 300 °C, 350 °C and 400 °C. The initiation of nano-pores within individual minerals and organic matter particles were observed and quantified alongside the evolution from nano-pores to micro-fractures. The major expansion of pore-volume occurred between 300 and 350 °C in both samples, with the pores elongating rapidly along the organic-rich bedding. The internal pressures induced by organic matter transformation influenced the development of microfractures. Mechanical properties and strain distributions within these two samples were modelled under a range of axial stresses using FEM. The results show that the overall stiffness of the shale reduced during heating, despite organic matter becoming stiffer. The varied roles of ductile (e.g., clay minerals, organic matter) and brittle materials (e.g., calcite, pyrite) within the rock matrix are also modelled and discussed. The configurations of organic matter, mineral components, porosity and connectivity impact elastic deformation during shale pyrolysis. This work extends our understanding of dynamic coupled processes of microstructure and elastic deformation in shales to the nm-scale, which also has applications to other subsurface energy systems such as carbon sequestration, geothermal and nuclear waste disposal.
|
Jun 2021
|
|
I12-JEEP: Joint Engineering, Environmental and Processing
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[13824, 9866, 17606]
Open Access
Abstract: Injection of CO2 into shale reservoirs to enhance gas recovery and simultaneously sequester greenhouse gases is a potential contributor towards the carbon-neutral target. It offers a low-carbon, low-cost, low-waste and large-scale solution during the energy transition period. A precondition to efficient gas storage and flow is a sound understanding of how the shale's micro-scale impacts on these phenomena. However, the heterogeneous and complex nature of shales limits the understanding of microstructure and pore systems, making feasibility analysis challenging. This study qualitatively and quantitatively investigates the Bowland shale microstructure in 3D at five length scales: artificial fractures at 10–100 μm scale, matrix fabric at 1–10 μm-scale, individual mineral grains and organic matter particles at 100 nm–1 μm scale, macropores and micro-cracks at 10–100 nm scale and organic matter and mineral pores at 1–10 nm-scale. For each feature, the volume fraction variations along the bedding normal orientation, the fractal dimensions and the degrees of anisotropy were analysed at all corresponding scales for a multi-scale heterogeneity analysis. The results are combined with other bulk laboratory measurements, including supercritical CO2 and CH4 adsorption at reservoir conditions, pressure-dependent permeability and nitrogen adsorption pore size distribution, to perform a comprehensive analysis on the storage space and flow pathways. A cross-scale pore size distribution, ranging from 2 nm to 3 μm, was calculated with quantified microstructure. The cumulative porosity is calculated to be 8%. The cumulative surface area is 17.6 m2 g−1. A model of CH4 and CO2 flow pathways and storage with quantified microstructure is presented and discussed. The feasibility of simultaneously enhanced gas recovery and subsurface CO2 storage in Bowland shale, the largest shale gas potential formation in the UK, was assessed based using multi-scale microstructure analysis. The potential is estimated to store 19.0–21.2 Gt CO2 as free molecules, together with 18.3–28.5 Gt CO2 adsorbed onto pore surfaces, implying a theoretical maximum of 47.5–49.5 Gt carbon storage in the current estimate of 38 trillion cubic metres (∼1300 trillion cubic feet) of Bowland shale. Simple estimates suggest 6.0–15.8 Gt CO2 may be stored in practice.
|
Jun 2021
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[15506]
Abstract: The feasibility and advantages of synchrotron imaging have been demonstrated to effectively characterise fracture initiation and propagation in shales during indentation tests. These include 1) fast (minute-scale) and high-resolution (μm-scale) imaging of fracture initiation, 2) concurrent spatial and temporal information (4D) about fracture development, 3) quantification and modelling of shale deformation prior to fracture. Imaging experiments were performed on four shale samples with different laminations and compositions in different orientations, representative of three key variables in shale microstructure. Fracture initiation and propagation were successfully captured in 3D over time, and strain maps were generated using digital volume correlation (DVC). Subsequently, post-experimental fracture geometries were characterized at nano-scale using complementary SEM imaging. Characterisation results highlight the influence of microstructural and anisotropy variations on the mechanical properties of shales. The fractures tend to kink at the interface of two different textures at both macroscale and microscale due to deformation incompatibility. The average composition appears to provide the major control on hardness and fracture initiation load; while the material texture and the orientation of the indentation to bedding combine to control the fracture propagation direction and geometry. This improved understanding of fracture development in shales is potentially significant in the clean energy applications.
|
Oct 2020
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[17314]
Open Access
Abstract: Barite precipitation in fractures and nanopores within a shale sample is analysed in situ, in 3D, and over time. Diffusion of barium and sulphate from opposite sides of the sample creates a supersaturated zone where barium sulphate crystals precipitate. Time-lapse synchrotron-based computed tomography was used to track the growth of precipitates over time, even within the shale’s matrix where the nanopores are much smaller than the resolution of the technique. We observed that the kinetics of precipitation is limited by the type and size of the confinement where crystals are growing, i.e., nanopores and fractures. This has a major impact on the ion transport at the growth front, which determines the extent of precipitation within wider fractures (fast and localised precipitation), thinner fractures (non-localised and slowing precipitation) and nanopores (precipitation spread as a front moving at an approximately constant velocity of 10 ± 3 µm/h). A general sequence of events during precipitation in rocks containing pores and fractures of different sizes is proposed and its possible implications to earth sciences and subsurface engineering, e.g., fracking and mineral sequestration, are discussed.
|
Aug 2019
|
|
I13-2-Diamond Manchester Imaging
|
Diamond Proposal Number(s):
[4022]
Abstract: Microstructures and pore systems in shales are key to understanding the role of shale in many energy applications. This study proposes a novel multi-stage upscaling procedure to comprehensively investigate the heterogeneous and complex microstructures and pore systems in a laminated and microfractured shale, utilizing 3D multi-scale imaging data. Five imaging techniques were used for characterisation from sub-nanoscale to macroscale (core-scale), spanning four orders of magnitude. Image data collected using X-ray tomography, Focused Ion Beam, and Electron Tomography techniques range in voxel size from 0.6 nm to 13 μm.
Prior to upscaling, a novel two-step analysis was performed to ensure sub-samples were representative. Following this, a three-step procedure, based on homogenising descriptors and computed volume coefficients, was used to upscale the quantified microstructure and pore system. At the highest resolution (nanoscale), four distinct pore types were identified. At the sub-micron scale equations were derived for three pore-associated phases. At the microscale, the volume coefficients were recalculated to upscale the pore system to the macroscale (millimetre). The accuracy of the upscaling methodology was verified, predicting the total porosity within 7.2% discrepancy. The results provide a unique perspective to understand heterogeneous rock types, breaking though prior scale limitations in the pore system.
|
Jun 2019
|
|
I12-JEEP: Joint Engineering, Environmental and Processing
|
Michael R.
Chandler
,
Anne-Laure
Fauchille
,
Ho Kyeom
Kim
,
Lin
Ma
,
Julian
Mecklenburgh
,
Roberto
Rizzo
,
Mahmoud
Mostafavi
,
Sebastian
Marussi
,
Robert
Atwood
,
Steven
May
,
Mohammed
Azeem
,
Ernie
Rutter
,
Kevin
Taylor
,
Peter
Lee
Diamond Proposal Number(s):
[13824]
Open Access
Abstract: Mode‐I Fracture Toughness, KIc, was measured in six shale materials using the double torsion technique. During loading, crack propagation was imaged both using twin optical cameras, and with fast X‐Ray radiograph acquisition. Samples of Bowland, Haynesville, Kimmeridge, Mancos, Middlecliff and Whitby shales were tested in a range of orientations. The measured fracture toughness values were found to be in good agreement with existing literature values. The two imaging techniques improve our understanding of local conditions around the fracture tip, through in‐situ correlation of mechanical data, inelastic zone size and fracture‐tip velocity. The optical Digital Image Correlation (DIC) technique proved useful as a means of determining the validity of individual experiments, by identifying experiments during which strains had developed in the two "rigid" specimen halves. Strain maps determined through DIC of the optical images suggest that the scale of the inelastic zone is an order of magnitude smaller than the classically used approximation suggests. This smaller damage region suggests a narrower region of enhanced permeability around artificially generated fractures in shales. The resolvable crack‐tip was tracked using radiograph data and found to travel at a velocity around 470μm.s−1 during failure, with little variation in speed between materials and orientations. Fracture pathways in the bedding parallel orientations were observed to deviate from linearity, commonly following layer boundaries. This suggests that while a fracture travelling parallel to bedding may travel at a similar speed to a bedding perpendicular fracture, it may have a more tortuous pathway, and therefore access a larger surface area.
|
Dec 2018
|
|