Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture characterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Field Experiments and Model Validation: The EGS Collab Project

Implementing enhanced geothermal systems (EGS) will require improvements in understanding stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The US Department of Energy (DOE) is addressing these and other challenges at multiple scales. The EGS Collab project, addressed here, is performing tests and modeling at the 10 m scale. The FORGE project is performing tests at the full reservoir scale. The EGS Collab team created an underground testbed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently building a second testbed at SURF at a depth of about 1.25 km aimed at investigating shear stimulation (Experiment 2). In Experiment 1, we drilled eight boreholes in a well-characterized phyllite and installed geophysical sensors in six of them to create a well-instrumented testbed to allow careful monitoring of stimulation events and flow tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results both in near-real-time and in detailed process studies. Stimulations performed in this testbed allowed quantification of processes occurring during stimulation and the examination of dynamic flow occurrences. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Our second testbed, targeted at shear stimulation, is currently being built at the SURF Facility at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1.

artificial intelligence↗

Hydraulic fracturing experiments at 1500 m depth in a deep mine: Highlights from the kISMET project

In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.

Oldenburg, C↗

Characteristics and mechanisms of hydrogen-induced quasi-cleavage fracture of lath martensitic steel

Here this study presents an in-depth characterization of the microstructures, crystallographic orientations, and dislocation characteristics beneath hydrogen-induced quasi-cleavage fracture features of an as-quenched, lath martensitic (α') 22MnB5 steel. The fracture surfaces of gaseous hydrogen-embrittled martensitic specimens were analyzed by a combination of multiple analytical tools: scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), transmission electron microscopy (TEM), and transmission Kikuchi diffraction (TKD). The dominant fracture mode in the hydrogen-affected zones was quasi-cleavage fracture, which involved significant plasticity, evidenced by plastic zones near tear ridges and a high density of dislocations beneath the quasi-cleavage facets. The martensite constituent sizes, variant orientations, and boundaries influenced the quasi-cleavage surface morphologies. Fractography revealed the occurrence of {100} α' -type cleavage across martensitic laths, developing relatively “flat” quasi-cleavage surfaces, in addition to {110} α' -type cracking likely along lath and block boundaries, and fracture along non-cleavage planes. The likelihood of the formation of the relatively “flat” quasi-cleavage surfaces increased with increasing martensitic constituent sizes. Substantial dislocation bands formed on intersecting {112} α' slip planes within a martensite block below the hydrogen-induced cleavage fractures on {100} α' planes. River markings on the quasi-cleavage surfaces were found to originate from the complex, hierarchical lath martensitic microstructure. Steps and ridges on the quasi-cleavage facets generally linked with the various martensitic boundaries, suggesting that they were produced as a result of crack deviations at those boundaries. The fracture paths and martensite quasi-cleavage mechanisms are discussed in the context of the role of hydrogen and Cottrell cleavage model.

36 MATERIALS SCIENCE↗

Newly acquired core enhances geologic understanding of the northern Paradox Basin Cane Creek play, southeastern Utah

The Cane Creek interval is an emerging unconventional heterolithic play within the Pennsylvanian Paradox Formation of the northern Paradox Basin, southeastern Utah. Structural controls (basement uplifts, salt movement, faults, fractures) likely control production from the Cane Creek, where total oil production of ~10 MMBO is just a fraction of current estimates of the undiscovered resource (up to 1.2 BBOE). Oil production has been most successful from unstimulated horizontal wells in the Cane Creek Unit in the central play area, whereas areas to the north (Greentown and Gunnison Valley) have to date proven to be largely unproductive from vertical and horizontal wells drilled without the benefit of 3D seismic or the deployment of modern hydraulic fracturing techniques. However, the limited historical drilling has shown significant promise (e.g., significant shows). The lack of core data in the northern part of the play has impeded thorough assessment of reservoir quality, character, and lithological facies heterogeneity. The Cane Creek is informally divided into three distinct zones. The upper A and basal C zones are composed of fabric destructive and wavy bedded anhydrite, dolomitic mudstone, and organic-rich algal laminated mudstones (source rock). The middle B zone is generally wave rippled and burrowed, low-permeability sandstone-siltstone (reservoir) with subordinate displacive anhydrite. Recent potash cores collected from the Salt Valley anticline near Crescent Junction, Utah (but east of the main play area), show thicker reservoir sandstone/siltstone packages (~40 ft) in the B zone, and parts of the C zone, in the northern part of the Paradox Basin in comparison to the central (~35 ft) and southern (~30 ft) areas. Source rock analyses from northern core/cuttings also indicate that it is well positioned within the oil window (VRo ~0.8), with thin organic-rich mudstone beds containing up to 20 wt% TOC. Because the siliciclastic reservoir rocks have low permeabilities (0.009–0.202 mD) and variable porosities (6–17%) due to dolomite-anhydrite cement, naturally occurring and possibly stimulated fractures may be essential for hydrocarbon recovery. To better understand the northern part of the Cane Creek play, a new core was obtained from the research stratigraphic well State 16-2 (January 2021). Details from this new core will help characterize reservoir quality via core analyses, petrography, fracture orientations, geomechanics, and correlation of facies and depositional trends. Additionally, an evaluation of matrix porosity/permeability will be conducted to determine its contribution to productivity.

02 PETROLEUM↗

Low-Rate Cycling Regeneration Sustains NMC811–Graphite Battery Performance through Extreme Ultrafast Charging

Urban air mobility (UAM) systems, such as electric vertical takeoff and landing (eVTOL) aircraft, require batteries capable of extreme ultrafast charging to support high-frequency operations. Here, in this work, we evaluate the electrochemical stability of NMC811–graphite full cells using a complex mixed-modal protocol designed to simulate the aggressive turnaround demands of UAM missions. We demonstrate that incorporating a periodic low-rate (0.3C) “regeneration” cycle enables cells to maintain 96% capacity retention over 900 cycles, despite repeated exposure to ultrafast pulses. To validate the efficacy of this regeneration mechanism, we performed simplified control experiments comparing continuous 10C charging against a 10C–0.3C interleaved protocol. Advanced characterization reveals that while mechanical particle fracture is pervasive at these rates, the regeneration cycles decelerate the electrochemical degradation by resolving spatially heterogeneous Ni oxidation states and relaxing lithium-ion concentration gradients. These findings provide a vital mechanistic framework for battery operation, proving that restorative charging can effectively delay and mitigate failure in demanding UAM applications.

batteries↗

EGS Collab Experiment 1: SIMFIP Notch-164 GRL Paper

Characterizing the stimulation mode of a fracture is critical to assess the hydraulic efficiency and the seismic risk related to deep fluid manipulations. We have monitored the three-dimensional displacements of a fluid-driven fracture during water injections in a borehole at ~1.5 km depth in the crystalline rock of the Sanford Underground Research Facility (USA). The fracture initiates at 61% of the minimum horizontal stress by micro-shearing of the borehole on a foliation plane. As the fluid pressure increases further, borehole axial and radial displacements increase with injection time highlighting the opening and sliding of a new hydrofracture growing ~10 m away from the borehole, in accordance with the ambient normal stress regime and in alignment with the microseismicity. Our study reveals how fluid-driven fracture stimulation can be facilitated by a mixed-mode process controlled by the complex hydromechanical evolution of the growing fracture. The data presented in this submission refer to the SIMFIP measurements and analyses of the stimulation tests conducted on the 164 ft (50 m) notch of the Sanford Underground Research Facility (SURF), during the EGS-Collab test 1. In addition to the datafiles, there is the draft of a manuscript submitted to Geophysical Research Letters (GRL).

15 GEOTHERMAL ENERGY↗

Mechanism of Microseismic Generation During Hydraulic Fracturing - With Evidence From HFTS 2 Observations

The objective of this study is to understand how microseismic events are generated during hydraulic fracturing, as well as the role of geomechanical conditions (i.e., stress and mechanical stratigraphy) in this process. In the industry, microseismic event clouds have been generally used as an “outer-boundary” of the “stimulated reservoir volume” (SRV). However, by comparing with other surveillance data (Low frequency Distributed Acoustic Sensing, or LF-DAS strain) in the Hydraulic Fracturing Test Site (HFTS) 2 experiment, we show that this assumption is fundamentally flawed. The HFTS 2 data has three unique observations that have not been commonly observed in other datasets: 1. Due to influence of offset pad depletion, microseismic data shows that hydraulic fractures from the child well can propagate over 3000 feet into the depleted low stress zone. 2. By comparing microseismic and horizontal fiber LF-DAS strain data, we observe that microseismic event cloud does not necessarily reflect the created hydraulic fracture volume. Particularly, the extent of microseismic event clouds near heel stages are much shorter than what is shown with LF-DAS strain data. 3. Microseismic event magnitudes are larger in the depleted regions. Through geomechanical analysis, we demonstrate that the “bedding-plane-slip” model is likely the mechanism for microseismic generation during hydraulic fracturing. Furthermore, this model successfully explains the above field observations from HFTS 2 experiment. We also provide a quantitative relationship connecting the microseismic event magnitude with fracture width increment and layer mechanical property contrast.

58 GEOSCIENCES↗

Towards an Era of Low-Temperature Integral Critical Experiments: Surrogate Testing of Low-Temperature TEX Configurations

United States and IAEA regulations for the transportation of radioactive materials require that fissile materials must be subcritical under expected ambient temperatures ranging from -40°C to 38°C. While historical experiments have been conducted at room and elevated temperatures, no criticality integral benchmarks exist for temperatures below room temperature. The Low-Temperature Thermal Epithermal eXperiments (LT-TEX) aim to provide benchmark evaluations down to -40°C. The LT-TEX design incorporates stacked HEU plates interspersed with varying thicknesses of polyethylene disks to provide varying levels of moderation, with five configurations at 20°C having already been benchmarked. A novel vacuum/cryogenic chamber was fabricated to contain the LT-TEX configurations in a thermally isolated environment and cool it to -40°C. Prior to testing with fissile materials, testing with surrogate fuel plates was performed to characterize the thermal properties of the system, ensure no condensation forms on the stack, and to test that the uranium plates would not warp, or fracture during the thermal cycling. Surrogate testing has characterized the cryogenic chamber and verified its viability to perform the benchmark experiments, allowing for the next steps in the LT-TEX experimental campaign. The results of the surrogate testing are discussed and compared to heat transfer calculations performed as part of the experiment design process, to establish confidence in the thermal behavior of the system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An adaptive space-time phase field formulation for dynamic fracture of brittle shells based on LR NURBS

Here, we present an adaptive space-time phase field formulation for dynamic fracture of brittle shells. Their deformation is characterized by the Kirchhoff–Love thin shell theory using a curvilinear surface description. All kinematical objects are defined on the shell’s mid-plane. The evolution equation for the phase field is determined by the minimization of an energy functional based on Griffith’s theory of brittle fracture. Membrane and bending contributions to the fracture process are modeled separately and a thickness integration is established for the latter. The coupled system consists of two nonlinear fourth-order PDEs and all quantities are defined on an evolving two-dimensional manifold. Since the weak form requires C 1 -continuity, isogeometric shape functions are used. The mesh is adaptively refined based on the phase field using Locally Refinable (LR) NURBS. Time is discretized based on a generalized-α method using adaptive time-stepping, and the discretized coupled system is solved with a monolithic Newton–Raphson scheme. Finally, the interaction between surface deformation and crack evolution is demonstrated by several numerical examples showing dynamic crack propagation and branching.

42 ENGINEERING↗

Two-Phase Fluid Flow Properties of Rough Fractures With Heterogeneous Wettability: Analysis With Lattice Boltzmann Simulations

Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.

58 GEOSCIENCES↗

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗

Characterizing Pore-Scale Geochemical Alterations in Eagle Ford and Barnett Shale from Exposure to Hydraulic Fracturing Fluid and CO 2 /H 2 O

As demand increases for an affordable energy source that is tied to an environmental obligation to reduce greenhouse gas emissions and water usage, there is a growing consideration in shale production utilizing processes such as 1) enhancing hydrocarbon recovery via carbon dioxide (CO 2 ) flooding, 2) using CO 2 as a fracturing agent to minimize water use, and 3) storing CO 2 in depleted shale formations to mitigate emissions to the atmosphere. Understanding the geochemical reactions and alterations that occur as shale is exposed to fluids and CO 2 is necessary to develop and optimize each of these processes for field applications. While the majority of shale formations are stimulated using traditional fracturing fluid, some may be fractured using CO 2 or other non-traditional means. We examine the effect fracturing fluid has on shale and how it behaves with secondary exposure to dry CO 2 or CO 2 -saturated water using in situ Fourier Transform infrared spectroscopy (FTIR), feature relocation scanning electron microscopy (SEM), and surface area and pore size analysis using volumetric gas sorption. These techniques were performed on Eagle Ford and Barnett shale samples that were exposed to fracturing fluid and unexposed (as received). Shales that have been exposed to traditional fracturing fluid experienced two reaction fronts. The first reaction front was formed during exposure to the fracturing fluid (pH of ~1.4). A secondary reaction front was formed as a result of CO 2 -saturated fluid exposure in the form of carbonic acid (pH ~5.6). These two different reaction mechanisms drove multiple dissolution and precipitation cycles which altered petrophysical properties of the shale and could lead to a significant impact on flow pathways. FTIR showed that equilibration of carbonate dissolution and precipitation cycles could take as long as 35 days. Samples exposed to fracturing fluid showed significantly less carbonate reactivity compared to those exposed to water. Pore size analysis results indicate exposure to fracturing fluid blocked small nanopores (10-15 nm) reducing BET surface area and total pore volume. SEM results show barite precipitated heavily during exposure to fracturing fluid. It appeared that carbonic acid was able to extract sulfur from organic matter to form gypsum evaporites. The mineralogical (barite precipitation and calcite dissolution/precipitation) and pore-scale alterations observed in these samples may lead to enhancement of flow pathways for injected CO 2 or produced hydrocarbons.

04 OIL SHALES AND TAR SANDS↗

Acoustic velocity and permeability of acidized and propped fractures in shale

We report that from geochemical reactions to proppant emplacement, hydraulic fracturing induces various chemomechanical fracture alterations in shale reservoirs. Hydraulic fracturing through the injection of a vast amount and variety of fluids and proppants has substantial impacts on fluid flow and hydrocarbon production. There is a strong need to improve our understanding on how fracture alterations affect flow pathways within the stimulated rock volume and develop monitoring tools. We have conducted time-lapse rock-physics experiments on clay-rich (carbonate-poor) Marcellus shales to characterize the acoustic velocity and permeability responses to fracture acidizing and propping. Acoustic P- and S-wave velocities and fracture permeability were measured before and after laboratory-induced fracture alterations along with microstructural imaging through X-ray computed tomography and scanning electron microscopy. Our experiments indicate that the S-wave velocity is an important geophysical observable, particularly the S-wave polarized perpendicular to fractures because it is sensitive to fracture stiffness. The acidizing and propping of a fracture decrease its elastic stiffness. This effect is stronger for acidizing, so it is possible that proppant monitoring will be masked by chemical alteration except when propping is highly efficient (i.e., most fractures are propped). However, fracture permeability is undermined by the softening of fracture surfaces due to acidizing, while being greatly enhanced by propping. These contrasting effects on fluid flow in combination with similar seismic attributes indicate the importance of experiments to improve existing rock-physics models, which must include changes to the rock frame. Such improvements are necessary for a correct interpretation of seismic velocity monitoring of flow pathways in stimulated shales.

58 GEOSCIENCES↗

Overview of Hydraulic Fracturing Test Site 2 in the Permian Delaware Basin (HFTS-2)

Here, the Hydraulic Fracturing Test Site 2 (HFTS-2) is a large collaborative field-based R&D program in the Permian Delaware Basin, funded by the US Department of Energy through the National Energy Technology Laboratory (NETL) and the E&P industry, with support from academia. The projects' main objective is to improve the understating of the hydraulic fracturing process through utilization of advanced diagnostics and collection of through-fracture cores to provide undisputable evidence and attributes of the created hydraulic fractures. At the HFTS-2, in excess of $30 million was used to perform hydraulic fracturing research focusing on the Wolfcamp formation at a field site hosted and operated by Occidental. In addition to the research data collected by the project, Occidental provided a significant amount of background data for about a dozen existing wells in the test area as well as access to previously collected core. Additional technical and laboratory support was provided by the program members. Building on learnings and unanswered questions from HFTS-1 in the Permian Midland basin, the HFTS-2 used eight new producing wells and two existing (parent) wells to perform hydraulic fracturing research. Multiple science wells were drilled to sample and characterize the subsurface, including the collection of 540 feet of core in a vertical pilot hole and 948 feet of high-angle through-fracture core. The project installed permanent fiber optic cables in 3 wells to monitor near wellbore signals during fracturing and to collect cross-well strain measurements. Additional advanced diagnostics included a significant formation evaluation program on the vertical whole core, multi array moment tensor inversion capable microseismic survey, multi-well time-lapse geochemistry analysis, analysis of proppant distribution in producing child and slant core well, and others. We will provide an overview of the HFTS-2 project, including list of the consortium members, details of the test site, experiments performed, and technologies tested.

58 GEOSCIENCES↗

Evaluating proxies for the drivers of natural gas productivity using machine-learning models

We report the extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

58 GEOSCIENCES↗

Multiscale Characterization of the Caney Shale — An Emerging Play in Oklahoma

From a hydrocarbon perspective, the Caney Shale has historically been evaluated as a sealing unit, which resulted in limited studies characterizing the rock properties of the Caney Shale and its suitability for hydraulic fracturing. The objective of our research is to help bridge the current knowledge gap through the integration of multiscale laboratory techniques and to characterize the macro- and microscale rock properties of the Caney Shale. We employed an integrated approach for the characterization of the Caney using 200 ft (61 m) of Caney core from a target well in southern Oklahoma. Core observation and petrographic analysis of thin sections were combined to characterize the general rock types and associated fabrics and textures. Mineralogical composition, pore system architecture, and rock fabric were analyzed using x-ray diffraction (XRD), scanning electron microscopy/energy dispersive x-ray spectroscopy (SEM/EDS), and focused ion beam (FIB)-SEM. In addition, rebound hardness and indentation testing were carried out to determine rock hardness (brittleness) and elasticity, respectively. With the integrated multiscale characterization, three mixed carbonate-siliciclastic rock types were identified — mudstone, calcareous siltstone, and silty carbonate — likely representing a spectrum of deposition from low to relatively high energy environments in the distal portions of a ramp system. Silty carbonate contains mostly interparticle pores. The calcareous siltstones and silty mudstones contain a combination of organic matter pores and interparticle pores. Each of the rock types shows unique mineralogical compositions based on XRD. The mudstone lithofacies has the highest clay content and the least carbonate content. Calcareous siltstones show moderate carbonate and clay content. Silty carbonate indicates the highest carbonate content with the least clay content. In an order of mudstone, calcareous siltstone, and silty carbonate, rebound hardness and Young’s modulus show an increasing trend. As a result of rock-fluid interactions, there are potential scaling reactions during completion and production that could ultimately affect permeability and production rates. Overall, the proposed multiscale integration approach is critical for the geologic characterization of most rocks. However, in shale reservoirs dominated by microporosity and microstructure where engineered fractures are expected to provide permeability at a reservoir scale, successful integration is essential. An optimized, integrated geological characterization of the Caney Shale that is well aligned with the engineering designs in drilling, completing, and producing wellbores will ultimately lead to optimal production while providing safe and environmentally responsible operations.

Wang, Yulun↗

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗