Search

Erosional scars, slumps, slides and debrites (mass-transport deposits) on submarine slopes form relief that influences turbidity current behaviour. However, the interaction of mass-transport deposit emplacement kinematics (i.e. rapid emplacement versus creep), the morphology of the evolving seafloor topography and subsequent flow types is complicated. This study describes two outcrop examples of deep-water, predominantly turbiditic, deposits overlying mass-transport deposits, from the Eocene slope succession of the Aínsa Basin (Spanish Pyrenees). In both examples, the mass-transport deposit substrate continued to creep contemporaneously with turbidity current deposition and bypass. In the first case study, structures in the mass-transport deposits are extensional and oriented parallel to flow. In the second, structures are compressional and oriented perpendicular to flow. Mudstones dominate the slope succession, but deposits overlying mass-transport deposits form sandstone-prone accumulations. Lateral confinement by flow-parallel extensional faults enhanced channel-formation. Channel incision occurred close to the exhumed fault plane on the hangingwall. Incision ceased as the fault gradually locked-up, and channels avulsed to the hangingwall of a newly active normal fault, while the abandoned channel was filled by a thinning and fining-upward succession. Barriers formed where the long axes of compressional anticlines developed in mass-transport deposits formed perpendicular to flow. Here, turbidites filled the bathymetric lows in the axes of synclines. Continued tightening of synclinal depocentres led to repeated stoss side trapping and upslope accretion of coarse-grained sediment accumulation. This generated thickening and coarsening-up, followed by thinning and fining-upward successions, tracking the decreasing, then increasing flow bypass associated with topographic expression on the mass-transport deposit. This study shows how post-emplacement creep, and orientation of topographic features on mass-transport deposits, influence the routing and deposition of contemporaneous turbidity currents, and illustrates examples of facies successions that could be misinterpreted as the product of other autogenic submarine slope processes, such as lateral migration of sinuous channels, or compensational lobe stacking.

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.

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 CO 2, 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.

Mud dominates volumetrically the fraction of sediment delivered and deposited in deep-water environments, and mudstone is a major component of basin-floor successions. However, studies of basin-floor deposits have mainly focused on their proximal sandstone-prone part. A consequent bias therefore remains in the understanding of depositional processes and stratigraphic architecture in mudstone-prone distal settings beyond the sandstone pinchouts of basin-floor fans. This study uses macroscopic and microscopic descriptions of over 500 m of continuous cores from research boreholes from the Permian Skoorsteenberg Formation of the Karoo Basin, South Africa, to document the sedimentology, stratigraphy, and ichnology of a distal mudstone-prone basin-floor succession. Very thin- to thin-bedded mudstones, deposited by low-density turbidity currents, stack to form bedsets bounded by thin packages (< 0.7 m thick) of background mudstones. Genetically related bedsets stack to form bedset packages, which are bounded by thicker (> 0.7 m thick) background mudstones. Stratigraphic correlation between cores suggests that bedsets represent the distal fringes of submarine fan lobe elements and/or lobes, and bedset packages represent the distal fringes of lobe complexes and/or lobe complex sets. The internal stacking pattern of bedsets and bedset packages is highly variable vertically and laterally, which records dominantly autogenic processes (e.g., compensational stacking, avulsion of feeder channels). The background mudstones are characterized by remnant tractional structures and outsize particles, and are interpreted as deposited from low-density turbidity currents and debris flows before intense biogenic reworking. These observations challenge the idea that mud accumulates only from hemipelagic suspension fallout in distal basin-floor environments. Thin background mudstones separating bedsets (< 0.7 m thick) are interpreted to mainly represent autogenically driven lobe abandonment due to up-dip channel avulsion. The thicker background mudstones separating bedset packages (> 0.7 m thick) are interpreted to dominantly mark allogenically driven regional decrease of sand supply to the basin floor. The recognition of sandstone-prone basin-floor fans passing into genetically linked distal fringe mudstones suggests that submarine lobes are at least ∼ 20 km longer than previously estimated. This study provides sedimentological, stratigraphic, and ichnological criteria to differentiate mudstones deposited in different sub-environments in distal deep-water basin-floor settings, with implications for the accurate characterization of basin-floor fan architecture, and their use as archives of paleoenvironmental change.

Purpose The proportion of people living in urbanised areas is predicted to rise to > 65% by 2050, and therefore, more humans than ever will be exposed to urban environmental pollution. Accumulation of organic and inorganic substances on street and road surfaces is a major global challenge requiring scientifically robust methods of establishing risk that inform management strategies. This aim of this contribution is to critically review the global literature on urban road–deposited sediment contamination with a specific focus on variability in sampling and analytical methods. Materials and methods In order to assess the concentration of contaminants in global road-deposited sediment (RDS), a comprehensive search of published RDS studies was completed. We review methodological approaches used in RDS studies to highlight the variability in datasets as a result of sampling technique, grain size fractionation, geochemical and mineralogical characterisation methods and establishing the influence of local geology on contaminant concentrations. We also consider emerging contaminants in RDS, and we provide a workflow diagram which promotes a standardised sampling and analysis regime that we believe can reduce data variability and promote collaboration when it comes to tackling the important issue of RDS contamination. Results and discussion Across the literature, Asia (except China) and Africa are underrepresented in RDS studies despite these continents having the largest and fastest growing populations, respectively. The removal of tetraethyl lead from gasoline produced a noticeable decrease in lead concentrations in global RDS, and platinum group element (PGE) concentrations in RDS were consistent with catalytic converter usage. Research into the impact of electric vehicles on non-exhaust emissions suggests other contaminants such as zinc may become more prominent in the future. Most RDS studies consider grain size fractions larger than > 20 μm due to sampling constraints despite RDS Conclusions This review highlights considerable variability in sampling and analytical approach which makes it difficult to identify broad global patterns in RDS contamination. To remove this variability from future RDS research, this review suggests a workflow plan which attempts to improve the comparability between RDS studies. Such comparability is crucial in identifying more discrete RDS trends and informing future emission policy.

Due to the fine-grained nature of shale, organic matter particles are generally micro- and nano-scale in size. Functional groups differ between different organic matter types and as such provide unique chemical information for organic matter. Micro-FTIR can provide direct measurement to characterize sample features at the micrometer scale. However, optical diffraction limits its application at the nanometer scale. As a non-destructive high-resolution scanning probe technique, atomic force microscopy (AFM) is very powerful in nanoscale research and has been widely used in the fields of polymers, semiconductors, electrochemistry and biology. To provide a better combination of AFM’s unique advantages with nanoscale chemical analysis, the AFM-IR technique has been developed in recent years and also attracted the attention of geologists to explore the application in geological materials.In this research, AFM-IR which is a quite new technique in geological research was used to investigate the in-situ geochemical characteristics of organic matter in shale. Nanoscale molecular composition of individual organic particles was captured nondestructively, and the distribution of typical functional groups was displayed via 2D IR mapping. In our samples, both alginite and inertinite display chemical homogeneity. The former is dominated by oxygenated and aliphatic contents which indicates a higher hydrocarbon generation potential, whereas the latter is dominated by aromatic carbon. In contrast, migrated solid bitumen particles show compositional heterogeneities at the nanometer scale as some are aromatic-rich and others are aliphatic-rich. Finally, linking this advanced nanochemical technique to potential applications in subsurface energy was explored. This research demonstrates that AFM-IR is a powerful tool to examine the in-situ nanoscale geochemical characteristics of different organic matter types, which can also provide implications for energy applications.

The composition and characteristics of organic matter within shales can reveal important information on hydrocarbon generation, migration and accumulation processes during catagenesis. It may also impact upon the chemical properties and mechanical behavior of shales acting as hydrocarbon reservoirs or caprocks/seals for carbon and subsurface energy storage. Here we apply advanced atomic force microscopy-based infrared spectroscopy (AFM-IR) to capture the in-situ chemical structure of nanoscale molecular composition of individual organic particles and map the distribution of typical functional groups nondestructively; something that cannot be achieved via conventional analysis since some methods damage the sample, some only collect bulk measurements and some do not achieve nano-scale resolution. The nanoscale geochemical characteristics of organic matter in thermally mature samples of a Baltic Basin shale, from Lithuania and the Bowland Shale, from the UK, obtained from AFM-IR, are compared in this study. Significant chemical heterogeneity can be observed in migrated solid bitumen. In contrast, graptolite, alginite and inertinite generally display homogeneous molecular composition. The potential role of laminae as barriers to the vertical migration of hydrocarbons is also apparent. This research provides insights into nanoscale geochemical information on different organic matter types in thermally-mature shales to enhance the understanding of molecular structure transformation of organic matter during hydrocarbon expulsion and migration processes, and also sheds light on implications of nanogeochemistry for enhanced gas recovery, carbon sequestration and underground hydrogen storage-related energy applications. We also characterised specific materials which have not been reported in other geological studies using AFM-IR. Some of them display similar characteristics to other organic matter particles in the samples, but the possibilities of unknown organic contaminants introduced to the sample during sample storage or preparation stages cannot be ruled out. A comprehensive analysis on this undefined organic material is provided. Standardized approaches during sample preparation and storage are required for future research to dislodge potential surface-attached contaminants.

Due to the complexity of sedimentary rocks, especially the fine-grained and highly-heterogeneous nature of some argillaceous caprocks in the subsurface, it is challenging to directly characterise the macroscale changes during CO 2 storage. The investigation of microstructural changes in CO 2-fluid-rock reactions helps provide a fundamental understanding of macroscopic properties as the microstructure constrains their hydraulic, mechanical and geochemical behaviour. Because of distinct petrophysical properties, the effect of geochemical reaction on caprocks and reservoir rocks is dramatically different. In order to provide a comprehensive knowledge of microstructural changes induced by geochemical reactions and facilitate the application to large scale carbon storage research, this review summarizes representative laboratory measurements conducted on typical caprocks (e.g. mudstone and evaporite) and reservoir rocks (e.g. sandstone and carbonates). Key findings include: (1) Water can generally facilitate CO 2-rock reactions. The effect of dissolution and precipitation on the pore structure largely depends on the inherent heterogeneity of rock (including primary structure, mineralogy, thermal maturity), and the environmental conditions (temperature, pressure, fluid). (2) Adsorption-induced swelling of clay minerals and organic matter in mudstone may increase or reduce the pore space, which is primarily controlled by the treatment pressure of CO 2 within the same caprock system. (3) Fines migration and the re-precipitation of minerals are the major factors causing pore throat blockage of sandstone during CO 2 flooding. (4) Highly conductive wormholes formed in carbonates during CO 2 percolation can significantly enhance fluid injectivity, whereas a long-term thermodynamic equilibrium in the subsurface tends to counteract the increase in porosity due to the healing mechanism of mineral precipitation. Further, how to link laboratory work with modelling to realise large-scale and temporal-scale prediction are summarised. Finally, promising developments and research directions are proposed.

Mineral authigenesis in mudstones responds differently across the compositional spectrum of fine-grained rocks, resulting in a complex array of possible rock fabrics that affect the mechanical capacity of mudstones. Within carbonate-rich mudstones, geochemical controls on silicate and carbonate mineral reactions are variable and poorly understood. Here we analyse the diagenetic minerals of the carbonate-rich Hodder Mudstone from the Carboniferous Bowland Basin, UK, using SEM petrography combined with X-ray-based inorganic geochemical data. Our findings show that up to 90% of quartz cements in the samples are diagenetic, originating mainly from biogenic silica dissolution and clay mineral reactions. It is also evident that due to the varying diagenetic silica-yielding reactions, authigenic silica overgrowths and isolated authigenic quartz crystals are localized in argillaceous samples, while calcareous samples are characterized by silica replacement textures and quartz/calcite intergrowths. Moreover, euhedral dolomite crystals are concentrated within argillaceous units, whereas calcareous units are characterized by anhedral dolomite precipitation and replacement textures. These finding presents a case for facies-selective cementation, as both early and burial diagenetic alterations were observed to be controlled by primary depositional components of biogenic debris and extrabasinal silicate detritus that resulted in complex and variable precipitation of authigenic minerals.

A novel type of burial dolomitization related to alteration of intermediate volcaniclastic rocks is discussed based on mineralogical, petrological and geochemical analyses. Lacustrine, fine-grained (

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.

Carbonate cementation of sandstones at channel bases is commonly interpreted to be a result of solutes generated from surrounding rocks in a burial-diagenetic regime. However, in rapidly buried and under-compacted conglomerates presented in this study, we interpret carbonate cements formed during early-diagenesis (1000–1600 m) to be precipitated directly from Ca- and CO2-rich saline formation water, derived from hypersaline lake water and driven by gravity and compaction. In the Bayindulan Sag of the Er'lian Basin, which is a continental rift basin, the sandstones and conglomerates of the A'ershan Formation (K1ba) were deposited in fan delta front settings during the syn-rift phase. Within these, tightly non-ferroan calcite-cemented conglomerates are present at the bases of distributary channels. In addition, some Fe-dolomite cements are present in sandstones with high present-day porosity, which are located in the upper parts of distributary channels. A petrographic analysis and fluid-inclusion microthermometry study show that the calcite cement predates the Fe-dolomite. Stable isotopic values indicate that both calcite and Fe-dolomite are characterized by enrichment of 13C and depletion of 18O (δ18OV-PDB: −15‰ to −8‰), and the value of δ13C V-PDB in calcite (+1.6‰ to +3.8‰) is less positive than that in Fe-dolomite (+3.1‰ to +8.4‰). Integrated analyses imply that the studied carbonate cements were probably formed within an open diagenetic system and were influenced by multiple diagenetic fluids. During deposition of the A'ershan Formation (K1ba), the lake water was brackish to hypersaline (Sr/Ba>1 and Sr content > 300 ppm). During deposition of the Tenggeer Formation (K1bt), the K1ba progressed rapidly into a deep burial phase (90 m/Ma), diagenetic fluids were controlled by hypersaline lake water, incursion of meteoric water occurred along faults, and methanogenic fermentation. Calcium was probably derived from the saline formation water and the transformation of smectite to illite in associated mudstones, leading to the precipitation of calcite (54–78 °C). In the Saihantala period, the K1ba experienced a slight uplift, the fault activity ceased, and the diagenetic fluids were derived from the saline formation water and were influenced by methanogenic fermentation, and Fe-dolomite began to precipitate. In deeper burial detrital feldspar grains and some calcite cement were probably dissolved by organic acid discharged from mudstones. In conclusion, the early diagenetic characteristics in the continental rift basins were controlled by the composition of the depositional waters and tectonics that controlled the opening and closing the diagenetic system thus influencing fluid composition.

Deep-water mudstones are often considered as background sediments, deposited by vertical suspension fallout, and the range of transport and depositional processes are poorly understood compared to their shallow-marine counterparts. Here, we present a dataset from a 538.50 m-thick coredsuccession through the Permian muddy lower Ecca Group of the Tanqua depocentre (southwest Karoo Basin, South Africa). This study aims to characterize the range of mudstone facies, transport and depositional processes, and stacking patterns recorded in deep-water environments prior to deposition of the Tanqua Karoo sandy basin-floor fans. A combination of macroscopic and microscopic description techniques and ichnological analysis has defined nine sedimentary facies that stack in a repeated pattern to produce 2–26-m-thick depositional units. The lower part of each unit is characterized by bedded mudstone deposited by dilute, low-density turbidity currents with evidence for hyperpycnal-flowprocesses and sediment remobilization. The upper part of each unit is dominated by more organic-rich bedded mudstone with common mudstone intraclasts, deposited by debris flows and transitional flows, with scarce indicators of suspension fallout. The intensity of bioturbation and burrow size increases upward through each depositional unit, consistent with a decrease in physicochemically stressed conditions, linked to a lower sediment accumulation rate. This vertical facies transition in the single well dataset can be interpreted to represent relative sea level variations, where the hyperpycnal stressed conditions in the lower part of the units were driven by sea level fall and the more bioturbated upper part of the unit represents backstepping, related to sea level rise. Alternatively, this facies transition may represent autogenic compensational stacking. The prevalence of sediment density flow deposits, even in positions distal or lateral to the sediment entry point, challenges the idea that deep-water mudstones are primarily the deposits of passive rainout along continental margins.

This paper quantitatively characterises the microstructure in shales across five scales in 3D, builds a multi-scale model of CH 4 and CO 2 flow pathways and storage, and assesses the potential of enhanced gas recovery and CO 2 storage simultaneously. Injection of CO 2 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 CO 2 and CH 4 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 m 2 g −1 . A model of CH 4 and CO 2 flow pathways and storage with quantified microstructure is presented and discussed. The feasibility of simultaneously enhanced gas recovery and subsurface CO 2 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 CO 2 as free molecules, together with 18.3–28.5 Gt CO 2 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 CO 2 may be stored in practice.

The Cenomanian-Turonian Oceanic Anoxic Event (OAE2) is a significant global event that has been linked to major source rock development in marine environments. The aim of this study is to integrate biostratigraphic, sedimentological and geochemical data (including stable carbon isotopes and trace elements) to document the timing and controls of organic matter enrichment across the Cenomanian-Turonian boundary.Morocco has excellent exposures of Tethyan carbonate platform and basinal facies and this paper details a high-resolution analysis of organic-rich Late Cretaceous intervals exposed in the Errachidia-Goulmima Basin, to investigate their relationship to the OAE2. Integrating planktic foraminiferal and ammonite biostratigraphy with carbon isotope correlations indicate that organic-rich black mudstones identified were deposited from latest Cenomanian to early Turonian, suggesting a late-OAE2 to post-OAE2 age.The sedimentology, redox water conditions and palaeoproductivity of the palaeoenvironments were studied to characterise the lithofacies and geochemical signatures, and to analyse the controlling factors for organic–enrichment. In the study area during the lower OAE2 interval, a shallow water mid-ramp environment prevailed, with predominantly bivalve-rich limestone deposition. Subsequent organic-rich mudstone development can be correlated with the latest Cenomanian (late OAE2 interval) and early Turonian (post-OAE interval) transgression, which allowed development of anoxic to euxinic conditions in the outer-ramp environment, characterized by high TOC values, from 2.7 to 17 wt%. Shallow platform locations continued to experience limestone deposition with no organic enrichment.Trace and major element data from the C/T black mudstones show an extremely high concentration in palaeoproductivity-sensitive (P, Zn, Ni and Cu) and redox-sensitive elements (V, Mo, and U). The redox environments evolved from dysoxia to anoxia and sulfidic euxinia, which is interpreted to have provided optimum conditions for organic carbon preservation in the Errachidia area.The results highlight the important interplay of local paleogeography and eustatic sea-level (the late Cenomanian/early Turonian transgressions) as the main controls on organic mudstone deposition, rather than simply global OAE2 control. The development of organic-rich mudstones reflects both increased productivity at this time and enhanced organic matter preservation due to the change in redox bottom water conditions.

Incorporation of reduced inorganic sulfur (S) into organic matter (OM) in euxinic environments is considered a vital route through which OM is preserved for millions of years. However, precipitation of iron sulfides is believed to compete with OM preservation via sulfurisation. High resolution geochemical, petrographic, and electron microprobe (EMPA) analyses were employed to better understand S controls on OM preservation with particular focus on iron sulfides precipitation in organic S-rich mudstones from the Kimmeridge Clay (KCF), Monterey (MF), and Whitby Mudstone [Grey Shale (GS) and Jet Rock (JR) Members] Formations. Analyses indicated that S in the KCF and MF is mostly bound to OM as shown by the relatively high abundance of organic sulfur compounds (OSCs) in their pyrolysates. The relatively low amounts of iron (Fe) in these mudstones permitted extensive reaction of S and OM to occur efficiently, promoting preservation of high amounts OM through sulfurisation. State-of-the-art EMPA imaging of total S and Fe indicated that the S was mainly bound to organics, supporting the geochemical analyses. In contrast, S in the GS and JR was predominantly sequestered as pyrite (FeS2) rather than organically bound, as indicated by the relatively low quantity/absence of OSCs in their pyrolysates, hence, preservation through sulfurisation was less significant in these settings. EMPA imaging of S in these mudstones shows discrete shapes ranging in geometry and size and in most cases correlated with Fe in the corresponding Fe EMPA images, confirming that the S exists in discrete authigenic FeS2. Overall, in line with previous research, this study indicates that preservation of OM through sulfurisation is an important pathway for sequestration of OM. However, an abundance of reactive Fe could hinder this process by preferentially reacting with reduced S, to form iron sulfides (mainly FeS2), thus impeding reactions between OM and S to occur efficiently.

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 characterised 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.

The Bell Island and Wabana groups together comprise a ca. 150-m-thick succession of interbedded clastic andchemical sedimentary rocks composed of eight distinct lithofacies that accumulated along the northern marginof the Rheic Ocean. Lithofacies stacking patterns indicate that deposition occurred during a marine transgressionwith superimposed small-scale sea level fluctuations producing at least six parasequences. Parasequences containingironstone are 10 to 20-m-thick and composed of hummocky cross-stratified sandstone interbeddedwith organic-rich mudstone and phosphatic Fe-silicate-bearing siltstone,which is overlain by hematitic granularironstone capped by an erosive flooding surface.This lithofacies association is interpreted to record the deposition of upwelling-related ironstone on a stormdominatedshelf. The close relationship between Fe-silicates and phosphorite typical of upwelling systemssuggests that Fe was delivered from deep, anoxic, nutrient-rich seawater that also stimulated high surface productivities.The result was the precipitation of authigenic sedimentary apatite in anoxic organic-rich sedimentsthat accumulated near the upwelling front. The gradual advection of Fe-rich waters away from the upwellingfront, initiated precipitation of Fe-silicate coated grains and cements in suboxic pore-waters. Iron pumpedinto shallower environments through advection and Fe-redox cycling is interpreted to have precipitated Fe-(oxyhydr)oxide grains in sediment of the oxygenated middle shelf. These coated grains were subsequentlyconcentrated by fairweather and storm currents on the shoreface to create granular economic Fe deposits.We challenge conventional models of Paleozoic ironstone deposition that rely on a continental source of Fe byproposing a hydrothermal source that supplied Fe2+ to the shelf through upwelling. It also highlights the potentialconnection between the delivery of anoxic, ferruginous seawater to the margins of the Rheic Ocean and theEarly Ordovician extinctions that punctuated the beginning of the Great Ordovician Biodiversification Event.

Produced waters from hydraulically fractured shale formations give insight into the microbial ecology and biogeochemical conditions down-well. This study explores the potential for sulfide production by persistent microorganisms recovered from produced water samples collected from the Marcellus shale formation. Hydrogen sulfide is highly toxic and corrosive, and can lead to the formation of "sour gas" which is costly to refine. Furthermore, microbial colonization of hydraulically fractured shale could result in formation plugging and a reduction in well productivity. It is vital to assess the potential for sulfide production in persistent microbial taxa, especially when considering the trend of reusing produced waters as input fluids, potentially enriching for problematic microorganisms. Using most probable number (MPN) counts and 16S rRNA gene sequencing, multiple viable strains of bacteria were identified from stored produced waters, mostly belonging to the Genus Halanaerobium, that were capable of growth via fermentation, and produced sulfide when supplied with thiosulfate. No sulfate-reducing bacteria (SRB) were detected through culturing, despite the detection of relatively low numbers of sulfate-reducing lineages by high-throughput 16S rRNA gene sequencing. These results demonstrate that sulfidogenic produced water populations remain viable for years post production and, if left unchecked, have the potential to lead to natural gas souring during shale gas extraction.

An appreciation of the diagenetic processes impacting basinal mudstones is important for understanding the properties of mudstones and the nature of fluids which might have been expelled onto adjacent carbonate platforms. This study utilizes petrographic and mineralogical data from mudstones from the Serpukhovian Morridge Formation in the Widmerpool Gulf, UK to investigate the paragenetic succession and associated mineral precipitation during diagenesis of mudstones. Authigenic minerals in the studied samples contain kaolinite, calcite, non-ferroan and ferroan dolomite, pyrite, ankerite, and quartz. The Si/Al ratios in most samples are above 2.5, indicating that most of the studied samples contain excess silica. The widespread presence of authigenic kaolinite and quartz provides important information on the sources for dissolved silica in pore water. The presence of pyritized radiolarians indicates that the dissolution of siliceous radiolarian skeletons is the most likely principle silica source within the system. The illitization of kaolinite also provided a possible silica source for quartz overgrowth during burial diagenesis. It has been long-proposed that the supply of base metals (i.e. Ba, Zn and Mg) for mineralization on the adjacent Derbyshire Platform came from juxtaposed basins, but evidence for this has been limited. The results of this study show that the diagenetic evolution of mudstones in the late Carboniferous links closely to the mineralization observed on the Derbyshire Platform. In the Widmerpool Gulf, the source of magnesium incorporated into dolomite is likely from basinal brine that also provided a source of barium for the barite mineralization. Based on the burial history model, the studied samples were buried to 1.3 km with a temperature of 82 °C in the late Carboniferous. Hydrocarbon was generated during maturation of organic matter and migrated with basinal brines through factures into surrounding formations during basin inversion. This mixed fluid likely contributed to the mineralization on the Derbyshire Platform.

Atomic force microscopy (AFM) has been widely used in materials science and chemistry research since it was invented in the 1980s. It is a relatively new tool for geosciences and is very powerful for micro- and nano-scale surface characterization with a resolution down to atomic scale. One of the major advantages of this microscopy technique is that it can perform continuous non-destructive measurements on a sample surface. Surface studies provide important structural information for geological materials, especially in energy geosciences. In this review, recent applications of AFM in energy geosciences are critically reviewed, including surface topography, mechanical properties, surface dynamics, and state-of-the-art AFM-based infrared spectroscopy (AFM-IR) and the AFM-Raman technique. AFM imaging enables the mapping and recording of surface information and enhances the understanding of morphology and mechanical properties of geological materials, and is especially useful for those displaying a strong heterogeneity. This review aims to demonstrate the accessibility, versatility, limitations, and promising developments of AFM techniques in energy geoscience. We highlight the major advantages and disadvantages of AFM analysis of geological materials and highlight where factors can cause potential deviations between different experimental results. We further highlight where we believe are the most promising potential uses of AFM. It is foreseen that with the improvement of experimental processes and hardware, AFM will provide indispensable information on surface characterization in a wide range of geological studies.

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 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 Finite-element modelling (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.

Analyzing the development of fracture networks in shale is important to understand both hydrocarbon migration pathways within and from source rocks and the effectiveness of hydraulic stimulation upon shale reservoirs. Here we use time-resolved synchrotron X-ray tomography to quantify in four dimensions (3-D plus time) the development of fractures during the accelerated maturation of an organic-rich mudstone (the UK Kimmeridge Clay), with the aim of determining the nature and timing of crack initiation. Electron microscopy (EM, both scanning backscattered and energy dispersive) was used to correlatively characterize the microstructure of the sample preheating and postheating. The tomographic data were analyzed by using digital volume correlation (DVC) to measure the three-dimensional displacements between subsequent time/heating steps allowing the strain fields surrounding each crack to be calculated, enabling crack opening modes to be determined. Quantification of the strain eigenvectors just before crack propagation suggests that the main mode driving crack initiation is the opening displacement perpendicular to the bedding, mode I. Further, detailed investigation of the DVC measured strain evolution revealed the complex interaction of the laminar clay matrix and the maximum principal strain on incipient crack nucleation. Full field DVC also allowed accurate calculation of the coefficients of thermal expansion (8 × 10−5/°C perpendicular and 6.2 × 10−5/°C parallel to the bedding plane). These results demonstrate how correlative imaging (using synchrotron tomography, DVC, and EM) can be used to elucidate the influence of shale microstructure on its anisotropic mechanical behavior.

Silica diagenesis has the potential to drastically change the physical and fluid flow properties of its host strata and therefore plays a key role in the development of sedimentary basins. The specific processes involved in silica diagenesis are, however, still poorly explained by existing models. This knowledge gap is addressed by investigating the effect of silica diagenesis on the porosity of Cenozoic mudstones of the North Viking Graben, northern North Sea through a multiple linear regression analysis. First, we identify and quantify the mineralogy of these rocks by scanning electron microscopy and X-ray diffraction, respectively. Mineral contents and host rock porosity data inferred from wireline data of two exploration wells are then analyzed by multiple linear regressions. This robust statistical analysis reveals that biogenic opal-A is a significant control and authigenic opal-CT is a minor influence on the porosity of these rocks. These results suggest that the initial porosity of siliceous mudstones increases with biogenic opal-A production during deposition and that the porosity reduction during opal-A/CT transformation results from opal-A dissolution. These findings advance our understanding of compaction, dewatering, and lithification of siliceous sediments and rocks. Moreover, this study provides a recipe for the derivation of the key controls (e.g., composition) on a rock property (e.g., porosity) that can be applied to a variety of problems in rock physics.

Organofacies analysis, a fundamental component within source rock appraisal based on the study of kerogen within a source rock, is typically produced from microscopy (palynological) and geochemical (kerogen kinetic) data, both of which are costly to acquire. One-Run-Fixed-Arrhenius (ORFA) kerogen kinetic analysis based on Rock–Eval pyrolysis offers a substantially cheaper kinetic dataset. Here, ORFA and palynological analyses are compared in organofacies characterization of a prospective Mississippian source rock reservoir (Bowland Shale, UK). Two-end-member organofacies were determined based on the abundance of the 56 kcal/mol activation energy peak derived from ORFA data: absence (< 5%) indicating ‘organofacies A’ containing the highest proportion of algal material (Type I kerogen); and presence (> 15%) indicating ‘organofacies B’ containing the highest proportion of sporomorphs (Type II kerogen). A mud-dominated slope setting for the rock reservoir was also used to test the accuracy of organofacies analysis in determining depositional environment. Organofacies A found within lithofacies deposited from dilute waning density flows and hemipelagic suspension settling occurred between shelf edge, slope and basin. Organofacies B found within lithofacies deposited from dilute waning density flows, and low-strength cohesive debrites occurred only within the lower slope. This study demonstrates that ORFA kerogen kinetic analysis provides comparable net results to palynological analysis, enabling cheaper and faster organic characterization during initial source rock appraisal. However, caution must be exercised in drawing interpretations as to biological source(s), organic matter mixing and preservation state(s) without additional investigation using data from detailed palynological analysis.

Here we describe the first automated fully integrated in-microscope broad ion beam (BIB) system. Ar+-BIB has several advantages over Ga+ focused ion beam (FIB) and Xe+ plasma-FIB (PFIB) methods inducing less beam damage, especially for ion beam sensitive materials. It can mill areas several orders of magnitude larger (up to millimetre scale), and is not confined to the edge of the sample with associated curtaining issues. BIB is shown to have sputter rates up to five times higher than comparable FIB techniques. This new coupled BIB–SEM system (commercial name ‘iPrep™II’) enables in-microscope surface polishing to remove contaminants or damage for two dimensional (2D) imaging, as well as automated serial section tomography (SST) by milling and imaging hundreds of slices, cost and time efficiently. The milled slice thickness can be controlled from a few nanometers up to a micrometre. A novel sample transfer, handling and interlock system allows automated and sequential BIB polishing, scanning electron microscopy (SEM) and analysis by secondary electron (SE) imaging, electron back scatter diffraction (EBSD) and energy dispersive spectroscopy (EDS) for 3D microstructure analysis. Furthermore, insulating surfaces can be sputter coated after milling each slice to reduce charging during SEM analysis. The performance of the instrument is demonstrated through a series of case studies across the materials, earth and life sciences exploiting the imaging, crystallographic and chemical mapping capabilities. These include the study of butterfly defects in bearing steels, meta-stable intermetallic phases in bronze bearings, shale gas rock, aluminium plasma electrolytic oxide (PEO) coatings as well as liver and mouse brain tissues.

Despite recent advances, diagenetic processes in fine-grained sediments are still relatively poorly understood. Key questions still to be resolved include the types of diagenetic minerals present in mudstones and the extent of element mobility in these low permeability systems. This study utilizes data from the Haynesville-Bossier Shale, USA, to analyze lithologies, discriminate authigenic phase, and identify mobile elements during diagenesis. It has implications for understanding how authigenic minerals develop and the sources of those authigenic minerals in fine-grained sediments.On the basis of grain-size and mineralogy five lithologies are designated: (i) silica-rich argillaceous mudstones; (ii) argillaceous siliceous mudstones; (iii) mixed siliceous mudstones; (iv) mixed mudstones; and (v) authigenically-dominated mudstones. The diagenetic development of the Haynesville-Bossier Shale can be divided in to early and late diagenesis. Ferroan and non-ferroan dolomite, framboidal pyrite and bioclast pore-filling kaolinite and calcite grain replacements and cements all formed during early diagenesis. Late diagenetic minerals phases include illite formed by the illitisation of smectite, replacive and displacive chlorite, calcite-replacive albite, quartz-replacive calcite and replacive and/or displacive quartz.The presence of extensive late diagenetic mineral precipitates indicates that there was a degree of element mobility on at least the local scale. Aluminium present in albite is most likely to have resulted from the illitization of smectite. Quartz overgrowths probably resulted from illitization and the pressure dissolution at quartz silt grain boundaries. Externally, hydrothermal fluids resulting from regional-scale igneous activity appear to have played a role in the formation of chlorite and possibly albite. The work indicates that extensive mineral development and element mobility occurred during late diagenesis.

Understanding the speciation and fate of radium during operational discharge from the offshore oil and gas industry into the marine environment is important in assessing its long term environmental impact. In the current work, 226Ra concentrations in marine sediments contaminated by produced water discharge from a site in the UK were analysed using gamma spectroscopy. Radium was present in field samples (0.1–0.3 Bq g−1) within International Atomic Energy Agency activity thresholds and was found to be primarily associated with micron sized radiobarite particles (≤2 μm). Experimental studies of synthetic/field produced water and seawater mixing under laboratory conditions showed that a significant proportion of radium (up to 97%) co-precipitated with barite confirming the radiobarite fate pathway. The results showed that produced water discharge into the marine environment results in the formation of radiobarite particles which incorporate a significant portion of radium and can be deposited in marine sediments.

As the fastest growing energy sector globally, shale and shale reservoirs have attracted the attention of both industry and scholars. However, the strong heterogeneity at different scales and the extremely fine-grained nature of shales makes macroscopic and microscopic characterization highly challenging. Recent advances in imaging techniques have provided many novel characterization opportunities of shale components and microstructures at multiple scales. Correlative imaging, where multiple techniques are combined, is playing an increasingly important role in the imaging and quantification of shale microstructures (e.g. one can combine optical microscopy, scanning electron microscopy/transmission electron microscopy and X-ray radiography in 2D, or X-ray computed tomography and electron microscopy in 3D). Combined utilization of these techniques can characterize the heterogeneity of shale microstructures over a large range of scales, from macroscale to nanoscale (c. 100–10−9 m). Other chemical and physical measurements can be correlated to imaging techniques to provide complementary information for minerals, organic matter and pores. These imaging techniques and subsequent quantification methods are critically reviewed to provide an overview of the correlative imaging workflow. Applications of the above techniques for imaging particular features in different shales are demonstrated, and key limitations and benefits summarized. Current challenges and future perspectives in shale imaging techniques and their applications are discussed.

Biogenic CH4 generation has been observed in many shallow, low temperature shale gas basins and coal seams. The depletion of conventional resources and the increasing demand of natural gas for human consumption have spurred the development of so-called unconventional gas resources such as shale gas (SG) and coal-bed methane (CBM). Such unconventional systems represent the opportunity for the stimulation of biogenic CH4 generation. Biogenic CH4 in shale and coal is produced by anaerobic biodegradation of organic matter (OM): methanogenic Archaea represent only the final step of biogenic CH4 generation. Several communities of microorganisms are involved in the initial breakdown of complex geopolymers and the production of intermediate compounds used by methanogens. There are several key knowledge gaps on biogenic CH4 production in unconventional gas systems, such as the exact fraction of bioavailable OM, the microbial communities involved and how they can be stimulated to enhance microbial methanogenesis. Progress on biodegradation studies, isotopic signatures, as well as DNA analyses and proteomics could help unravel interactions within the syntrophic community involved in the methanogenic biodegradation of OM. Questions also remain regarding the environmental impact of unconventional gas production, such as water quality and the mobility of toxic metals and radionuclides. The answers to these questions might have implications for both recovery practices and a sustainable development of unconventional resources. This review summarises the current knowledge regarding biogenic CH4 in SG and CBM: from the nature of the rocks to the producing microbial community and the indicators of biogenic CH4, illustrating how these two environments show remarkably similar opportunities for the stimulation of biogenic CH4 generation.

Silica cementation exerts a key control on the compaction and geotechnical properties of mudstones, and byextension, the style of hydrocarbon and/or mineral systems present in a given sedimentary basin. Integratedmicroscopic and bulk geochemical observations demonstrate that siliceous mudstones in the Bowland ShaleFormation, a target for UK shale gas extraction, exhibit abundant dispersed, discrete, μm-scale quartz cements,and exhibit silica enrichment (‘excess’) above a local detrital Si/Al threshold of 2.5. Dissolution of siliceousradiolarian tests during early diagenesis is identified as the main source of silica (opal A) required for quartzprecipitation, either via opal CT or directly to quartz, and potentially generated as a product of anoxic marine‘weathering’ (dissolution) of reactive silicates during early diagenesis. Excess silica correlates with free hydrocarbons(S1) normalised to total organic carbon (oil saturation index; OSI); we propose early diagenetic quartzprecipitation suppressed pore collapse (‘buttress effect’), retaining the pore space capacity to host oil. Quartzprecipitation was likely catalysed, for example via low porewater pH, elevated Al and/or Fe oxide content, and/or abundant labile organic matter. Juxtaposition of siliceous mudstones and mudstones lacking quartz cementindicates silica was immobile beyond the bed scale. Thus metre-scale siliceous packages likely represent moreprospective units within the Bowland Shale (in terms of unconventional hydrocarbons), on the basis of earlydiagenetic biogenic-derived quartz cementation leading to improved hydrocarbon storage capacity coupled toenhanced brittleness. These findings are relevant for shale oil and shale gas systems, specifically where oilretained in pores subsequently cracks to generate gas. These findings also suggest the Bowland Shale is a subclassof black shale, defined by the potential to host a relatively large volume of early diagenetic fluids,derived from anoxic bottom waters, which were potentially S- and/or metal-bearing. This is potentially relevantfor understanding the genesis of adjacent and related Pb-Zn mineral deposits.

Summary Despite recent advances, diagenetic processes in fine-grained sediments remain relatively poorly understood. Key questions still to be resolved include how mineral cements develop in low porosity and permeability systems and the extent of element mobility during the diagenesis of fine-grained sediments. This study uses petrographic, mineralogical and geochemical data from the Haynesville-Bossier Shale identify early and late diagenetic mineral development. On the basis of grainsize and mineralogy five facies are designated: (i) silica-rich argillaceous mudstones, (ii) argillaceous siliceous mudstones, (iii) mixed siliceous mudstones, (iv) mixed mudstones and (v) cemented mudstones. The diagenetic development of the Haynesville-Bossier Shale can be divided in to early and late diagenesis. Ferroan and non-ferroan dolomite, framboidal pyrite, and microfossil pore-filling kaolinite cements all formed during early diagenesis. Late diagenetic cements include microcrystalline calcite, euhedral pyrite, replacive and displacive chlorite, calcite-replacive albite, and replacive and/or displacive quartz. Late diagenetic cements would have required a wide range of geochemically mobile elements, (bed scale); illitisation of smectite during late diagenesis is a potential source for these elements. The work impacts porosity and permeability evolution and fracture susceptibility in unconventional shales.

In this research the Haynesville-Bossier Shale play (USA), which is a clay-bearing siliceous, organic-rich, gas-mature shale play, was characterized over pore diameters ranging 2 nm to 3000 nm. Three advanced imaging techniques were utilized correlatively, including the application of Xe plasma focused ion beam scanning electron microscopy (plasma FIB or PFIB), complemented by the Ga FIB method which is now frequently used to characterise porosity and organic/inorganic phases, together with transmission electron microscope tomography of the nano-scale pores. This is the first study to comprehensively develop pore geometric and network models in 3D based on experimental characterisation over these critical length scales

The Lower Permian Collyhurst Sandstone Formation is an important gas reservoir both offshore and onshore in the East Irish Sea Basin of north-western UK, and is considered to be a tight gas reservoir. Diagenetic processes have significantly influenced the reservoir quality of the sandstone. This study employs petrographic, mineralogical, and stable isotopic analyses on two different onshore core successions to investigate the composition and texture, the diagenetic history, and the impact of lithology and diagenesis on reservoir quality. The Collyhurst Sandstone Formation consists of fine- to very coarse-grained, poorly to well-sorted, sub-arkosic to lithic arenite sands. The reservoir properties of the sandstone are relatively poor, with average porosity and permeability of 6.6% and 0.387 mD respectively. On the basis of petrographic analyses, five sandstone lithofacies were identified: fine-grained, conglomeratic, carbonate-cemented, quartz-cemented and matrix-free. Contrasting characteristics between the reservoir units suggest the development of divergent diagenetic pathways. Features observed during early diagenesis include compaction, growth of pore-filling early carbonate cements (calcite I, dolomite and siderite I) and vermicular kaolinite. Late diagenetic features include: further carbonate cementations (ferroan dolomite/ankerite, calcite II and siderite II-III), clay mineral cementation (illite and chlorite) with both grain-coating and pore-filling textures. Also common in late diagenesis are syntaxial quartz overgrowths. Albite, barite and anhydrite are minor late diagenetic cement phases. Compaction processes rather than cementation accounts for the major destruction of the initial sandstone porosity. But cementation still accounts for up to 32% of porosity loss. Carbonate cements mostly fill intergranular pores, thus block pore throat and prevent migration of fluid. Clay minerals bridge and occlude pores, resulting in further destruction of the reservoir quality. The study provides insights into how lithofacies (primary composition, texture and sorting) and diagenesis can influence reservoir qualities of tight gas sandstones.

The microstructure of a highly laminated Lower Bowland Shale sample is characterized at the micron-to millimeter scale, to investigate how such characterization can be utilized for microstructure-based modelling of the shale's geomechanical behavior. A mosaic of scanning electron microscope (SEM) back-scattered electron (BSE) images was studied. Mineral and organic content and their anisotropy vary between laminae, with a high variability in fracturing and multi-micrometer aggregates of feldspars, carbonates, quartz and organics. The different microstructural interface types and heterogeneities were located and quantified, demonstrating the microstructural complexity of the Bowland Shale, and defining possible pathways for fracture propagation. A combination of counting-box, dispersion, covariance and 2D mapping approaches were used to determine that the total surface of each lamina is 3 to 11 times larger than the scale of heterogeneities relative to mineral proportion and size. The dispersion approach seems to be the preferential technique for determining the representative elementary area (REA) of phase area fraction for these highly heterogeneous large samples, supported by 2D quantitative mapping of the same parameter. Representative microstructural models were developed using Voronoï tessellation using these characteristic scales. These models encapsulate the microstructural features required to simulate fluid flow through these porous Bowland Shales at the mesoscale.

Pore characterization in shales is challenging owing to the wide range of pore sizes and types present. Haynesville-Bossier shale (USA) was sampled as a typical clay-bearing siliceous, organic-rich, gas-mature shale and characterized over pore diameters ranging 2 nm to 3000 nm. Three advanced imaging techniques were utilized correlatively, including the application of Xe+ plasma focused ion beam scanning electron microscopy (plasma FIB or PFIB), complemented by the Ga+ FIB method which is now frequently used to characterise porosity and organic/inorganic phases, together with transmission electron microscope tomography of the nano-scale pores (voxel size 0.6 nm; resolution 1–2 nm). The three pore-size scales each contribute differently to the pore network. Those 100 nm (greatest total volume hence determines fluid storativity). Four distinct pore types were found: intra-organic, organic-mineral interface, inter-mineral and intra-mineral pores were recognized, with characteristic geometries. The whole pore network comprises a globally-connected system between phyllosilicate mineral grains (diameter: 6–50 nm), and locally-clustered connected pores within porous organic matter (diameter: 200–800 nm). Integrated predictions of pore geometry, connectivity, and roles in controlling petrophysical properties were verified through experimental permeability measurements.

À l’échelle du laboratoire, les roches argileuses sont des matériaux hétérogènes dont le comportement thermo-hydromécanique est en grande partie contrôlé par la microstructure. Le choix du nombre et de la taille des échantillons à étudier en laboratoire est déterminant pour appréhender la variabilité des propriétés de la roche argileuse à petite échelle. Cet article présente une méthode statistique permettant de préciser la surface (ou le volume) et le nombre d’échantillons à prendre en compte pour qu’une propriétépchoisie caractérisant la microstructure, soit statistiquement représentative. Initialement établie dans un cas général par Kanitet al.(2003. Determination of the size of the representative volume element for random composites: statistical and numerical approach.Int J Solids Struct40(13–14): 3647–3679), cette méthode consiste à partitionner un échantillon de propriété moyenne [see formula in PDF] connue, en sous-échantillons de surfaceD × Dafin de calculer l’écart-type et l’erreur relative de la mesure depen fonction deD. Cette méthode permet ainsi de définir des surfaces élémentaires représentatives depen tenant compte de l’erreur relative par rapport à [see formula in PDF]. La méthode est d’abord présentée dans des cas généraux en 2D et 3D, et un exemple type est ensuite développé en 2D pour caractériser la fraction argileuse d’une lamine sédimentaire de Bowland (Royaume-Uni). La fraction surfacique argileuse est choisie comme propriétép, à partir d’une image grand-champ en microscopie électronique à balayage. La méthode est applicable en 2D et 3D sur les matériaux finement divisés autant sur les roches que sur les sols argileux, tant que l’échantillon considéré contient suffisamment d’éléments figurés (inclusions rigides ou pores dans une matrice par exemple) pour permettre l’utilisation des statistiques. L’apport principal visé pour la communauté des ingénieurs est dans la mesure du possible un meilleur ciblage de la quantité d’échantillons à prélever en forage pour mieux évaluer la variabilité des paramètres macroscopiques des roches argileuses. Les limites de la méthode sont ensuite discutées.

An understanding of the nature and scales of diagenetic variability within organic-rich mudstones is critical to the accurate assessment of shale-gas reservoir properties, as well as for elucidating chemical evolution pathways within mudstones. Here we integrate field observations with thin section descriptions (optical and electron optical techniques) and mineralogical data for the Blackhawk Member time-equivalent Mancos Shale in Book Cliffs, Utah, to determine the impacts of early and burial diagenesis on this mudstone succession.The detrital assemblage in the Mancos Shale comprises quartz-silt, feldspar, clay minerals, dolomite and organic matter (TOC of 1 to 2.5%). Biogenic silica is negligible. Field mapping reveals laterally continuous (km scale), ferroan dolomite cemented units up to 0.3 m thick, are present. These cemented units cap both coarsening-upward units (1 to 3 m thick), and stacked successions of coarsening-upward units (5 to 15 m thick). These upward-coarsening sediment packages, capped by dolomite cemented strata, correlate to bedsets and parasequences in updip settings. Pervasive cementation in these dolomite-cemented units is likely to have occurred prior to compaction as a result of bacterially mediated respiratory processes. Cementation at these levels is particularly evident because cement precipitation occurred during breaks in sediment accumulation below marine flooding surfaces. The abundance of dolomite cements highlights the importance of macroscopic-scale diagenetic carbonate mobility in these mudstones.In addition to carbonate-cements, diagenetic alteration and precipitation of quartz and alumina-silicate minerals are also important in these mudstones. Kaolinite is present both in uncompacted test of organisms and as vein fills in septarian concretions. Kaolinite precipitation is interpreted to have occurred prior to significant compaction and indicates that both silicon and aluminium were mobile during early diagenesis. We interpret the abundance of early diagenetic kaolinite cement to be the result of Al-mobilization by organic acids generated during organic matter oxidation reactions, with the Al sourced from poorly crystalline detrital aluminium oxides and clay minerals. There is also indirect evidence for burial diagenetic kaolinitization of feldspar grains. Quartz cement takes the form of quartz overgrowths and microcrystalline quartz crystals. Textures and CL spectra for the quartz microcrystalline cement suggests that recrystallization of biogenic silica (opal-A) was likely to have been an important source for quartz cements, although smectite-to-illite transformation may have contributed some. These mineral phases highlight that microscopic-scale diagenetic mobility of silica is important, even within mudstones lacking obvious sources of biogenic silica and is likely to be an important processes in a wide range of mudstones.

Diagenesis significantly impacts mudstone lithofacies. Processes operating to control diagenetic pathways in mudstones are poorly known compared to analogous processes occurring in other sedimentary rocks. Selected organic-carbon-rich mudstones, from the Kimmeridge Clay and Monterey Formations, have been investigated to determine how varying starting compositions influence diagenesis. The sampled Kimmeridge Clay Formation mudstones are organized into thin homogenous beds, composed mainly of siliciclastic detritus, with some constituents derived from water-column production (e.g., coccoliths, S-depleted type-II kerogen, as much as 52.6% total organic carbon [TOC]) and others from diagenesis (e.g., pyrite, carbonate, and kaolinite). The sampled Monterey Formation mudstones are organized into thin beds that exhibit pelleted wavy lamination, and are predominantly composed of production-derived components including diatoms, coccoliths, and foraminifera, in addition to type-IIS kerogen (as much as 16.5% TOC), and apatite and silica cements. During early burial of the studied Kimmeridge Clay Formation mudstones, the availability of detrital Fe(III) and reactive clay minerals caused carbonate- and silicate-buffering reactions to operate effectively and the pore waters to be Fe(II) rich. These conditions led to pyrite, iron-poor carbonates, and kaolinite cements precipitating, preserved organic carbon being S-depleted, and sweet hydrocarbons being generated. In contrast, during the diagenesis of the sampled Monterey Formation mudstones, sulfide oxidation, coupled with opal dissolution and the reduced availability of both Fe(III) and reactive siliciclastic detritus, meant that the pore waters were poorly buffered and locally acidic. These conditions resulted in local carbonate dissolution, apatite and silica cements precipitation, natural kerogen sulfurization, and sour hydrocarbons generation. Differences in mud composition at deposition significantly influence subsequent diagenesis. These differences impact their source rock attributes and mechanical properties.

Pore systems of tight-gas reservoirs control important parameters such as flow capacity, producible pore volumes, hydrocarbon flow rates and storage which are complex to characterise and to quantify. Detailed understanding of 3D pore structure is vital in predicting tight-gas reservoir productivity. This study investigates the pore architecture of a sandstone sample from Collyhurst Sandstone Formation, northern England, a tight gas reservoir. X-ray Computed Tomography (XCT) was employed to quantify the pore system (pore size, shape, geometry, and connectivity). Three types of pores were identified: inter-grain pores which are at the boundary of detrital grains and/or authigenic minerals, intra-grains pores which are pores within the detrital mineral grains, and intra-cement pores which are pores within the cement phase (authigenic minerals). The first pore type is primary in origin, whereas the latter two are secondary pores. A centerline tree and a pore network model (PNM) algorithm were applied for pore network analysis and a porosity value was estimated. The image-based porosity value was validated with laboratory experimental measurement. 3D image estimates 3.1% volume porosity, while 4.6% was recorded from the laboratory experiment. Spatial quantification of pore types and analysis of pore connectivity of the studied sandstone improves understanding on gas flow behaviour and gas storage. Furthermore, the results are useful in predicting permeability and multi-flow phase properties in a tight gas reservoir.

Source-rock reservoirs are fine-grained petroleum source rocks (“shales” or “mudstones”) having geomechanical properties that allow those rocks to produce hydrocarbons at economic rates after stimulation by hydraulic fracturing. Many of the assumptions commonly adopted by geophysicists to characterize shales cannot be applied to source-rock reservoirs. For example, the mineralogies of many source-rock reservoirs are not dominated by clay minerals and so mathematical and/or conceptual models developed for clay-dominated mudstones are not appropriate and cannot be applied to them. Instead, mudstones of shale plays are generally dominated by biogenic calcite and/or quartz. We use terminology of sedimentary geology to show that anisotropy is scale-dependent in source-rock reservoirs, and we discuss the depositional and diagenetic processes that control these and other geophysical properties of interest. The mudstones of source-rock reservoirs may or may not be anisotropic at the lamination scale (i.e., millimeters), the scale commonly used to measure anisotropic parameters via core plugs, but they are nearly always anisotropic at the bedset (centimeters to several meters) and member (tens of meters) scales. Because of the anisotropic nature of mudstones, elastic properties are not scalars at the length/thickness scales that can be defined using seismic methods. Properties of interest are likely to be different parallel to bedding compared to perpendicular to bedding. Because of the subseismic scale of much of this variability, thin-bed effects are likely to influence the AVO behavior of source-rock reservoirs.

In sandstone reservoirs, despite grain rearrangement during compaction, significant pore volumes may be retained prior to the onset of late diagenetic quartz cementation. In mudstone reservoirs, grain rearrangement during compaction results in significant pore volume reduction prior to late diagenesis. Where quartz overgrowths have been previously reported in post-compaction mudstones they have been in volumetrically low concentrations and interpreted as anomalous occurrences. Quartz cementation alters rock brittleness resulting in changes to mechanical fracture properties. Quartz overgrowths reduce reservoir porosity and permeability. We present petrographic evidence of two phases of quartz cement in the Haynesville-Bossier Shale: (i) grain replacive and (ii) quartz overgrowths. Carbonate grain replacement is volumetrically low (< 1%). Quartz overgrowths identified from SEM-CL imaging are volumetrically more significant (8–13%). Quartz overgrowths were most commonly observed in the coarse mudstone microfacies, but are present in both medium and fine mudstone microfacies. Petrographic evidence indicates three processes in the development of quartz overgrowths. Mica and pyrite are (i) engulfed and (ii) displaced by quartz overgrowth cement. The absence of a supportive, primary granular framework surrounding engulfed detrital and early authigenic minerals would indicate that quartz overgrowths are also replacive. Pressure dissolution of quartz silt grains and smectite-to-illite transformation are likely sources of silica for quartz cement. This study is the first to document large-scale, replacive, authigenic quartz overgrowth development within a producing mudstone.

The Lambda-Mu-Rho (LMR) approach of characterizing reservoir intervals in seismic data is a useful tool for efficient well placement and completion. The successful exploitation of unconventional hydrocarbons in the UK could play a fundamental role in helping to satisfy the growing energy demand. To achieve this, it has been proved necessary to differentiate the more brittle from ductile rocks to maximize reservoir recovery. Economic gas production in unconventional plays requires hydraulic fracture stimulation. This study utilizes the impedance, density and the Lame parameters (Lambda, λ and Mu, μ) logs derived from the Amplitude Variation with Offset (AVO) inversion to compute the more seismically intuitive LMR logs and the elastic properties Poisson’s Ratio and Young’s Modulus. The optimal Carboniferous tight sandstone reservoir intervals at Kirby Misperton in the Cleveland Basin were successfully characterized using the LMR crossplot method. Tight gas-filled sandstones are characterized by a low λρ and a high μρ value as a result of their low incompressibility, high rigidity and intrinsic brittleness. The analysis confirmed that the Lower Bowland sequence is prospective from an unconventional perspective due to its high quartz content and therefore is likely to respond positively to hydraulic fracturing. Unconventional resources currently account for 43% of US gas production with gas consumption expected to increase between 2010 and 2030 (Rogers, 2011). Tight gas sandstones represent approximately 70% of all global unconventional gas production and reserves that are yet to be developed (Khlaifat et al., 2011). Although commercial exploration and production of unconventional gas in the UK is in its infancy, it offers significant potential for future gas production growth and will become increasingly important as the easily exploitable conventional gas reserves are exhausted (Law and Curtis, 2002). Most tight gas reservoirs are often characterized by a multi-layered reservoir system that are often hydraulically fractured to enhance recovery rates from the reservoir intervals (Meckel and Thomasson, 2008).