Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fracture control”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Autonomous Fracture Conductivity Using Expandable Proppants in Enhanced Geothermal Systems

Summary Early thermal breakthrough in enhanced geothermal systems (EGS) due to the presence of preferential flow channels is a major challenge that endangers efficient and economic heat extraction in such systems. Previous studies mainly focused on adjusting circulation rates of the working fluid, which still leaves significant amounts of untapped heat behind. Currently, there is a lack of technologies for altering flow distribution within the fracture network to achieve uniform heat sweeping in the reservoir. This work presents a novel concept for making proppants to autonomously control fracture conductivity based on the surrounding temperature. Here, proppants with negative thermal expansion coefficients have demonstrated the capability for appropriate fracture conductivity adjustment as a function of temperature to achieve uniform flow across the fracture network. Particle-particle interactions governing such functions are explicitly modeled, and then the Lattice Boltzmann methods (LBM) is used to determine the potential impact of closure stress and temperature changes on the permeability of the proposed proppant packs. Microscale analyses are further used to determine the required material properties to achieve a certain improvement in the permeability of the proppant pack. Our analyses show an enhancement in permeability and the associated fracture conductivity by half of their initial values. Field-scale analysis further confirms the effectiveness of the proposed concept as 31.4% more heat can be extracted from EGS over 50 years of production when the proposed proppants are used. Such novel proppants may effectively delay thermal breakthrough, sweep heat from larger rock volumes, and elongate the life span of the EGS project.

Engineering↗
DOE OSTI2020

Application of an Embedded Fracture and Borehole Modeling Approach to the Understanding of EGS Collab Experiment 1

EGS Collab is a series of meso-scale experiments and associated numerical simulation activities being funded by the United States Department of Energy, Geothermal Technologies Office (GTO) to investigate enhanced geothermal system processes under in-situ stress and slightly elevated temperature conditions. This project is designed to provide scientists and engineers with immediate access to impermeable rock at scales larger than possible in the laboratory, but generally smaller than those for commercial production. Immediate access to rock is provided via the existing drifts of the former Homestake Gold Mine, now operated as the Sanford Underground Research Facility in Lead, South Dakota. The objectives of the EGS Collab project are to develop well controlled fracture networks between injector and producer boreholes using normal and shear stimulation for permeability enhancement. The first experimental site was located off the West Access drift on the 4850 Level (4850 feet below ground surface) in phyllite of the Precambrian Poorman formation, and involved the creation of a fracture network comprising a combination of hydraulic and natural fractures. A second experimental site is now being considered near the battery alcove on the 4100 Level in amphibolite and rhyolite of the Yates Unit. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. Experiment 1 has comprised a series of successful tests, including a long-term chilled-water circulation test, but the testbed has two atypical EGS elements. Active ventilation in the adjacent drift over a 50-year period cooled the testbed rock mass and hydraulic fracturing intersected monitoring boreholes making them conduits for fluid flow. Numerical simulations executed in support of the design of the EGS Collab Experiment 1, computed a radial temperature distribution and stress gradient orthogonal to the drift axis, resulting in the forecast of an oblong hydraulic fracture geometry, extended in the direction of the drift from the stimulation borehole. This paper describes the numerical simulation of long-term chilled-water circulation test with an embedded fracture and borehole modeling approach. The principal objective of the simulation work is show agreement between experimental observations in terms of production fluid temperatures, monitoring borehole temperatures, injection and production pressures, and tracer recoveries using a single conceptual model for the fracture network, using characterization and monitoring data generally available for EGS. A secondary objective is to improve the numerical simulation result comparisons with additional data available from the broader monitoring equipment within the testbed versus deeper, hotter, and more remote EGS. One complicating factor for comparisons between the numerical simulation results and experimental observations, is the Joule-Thomson heating associated with large pressure drops across the fracture network.

numerical simulation, embedded fractures and boreh↗
DOE OSTI2022

Well-Controlled Hydro-thermo-chemical Stimulation of a Coso Geothermal Well to Enhance Productivity

The goal of this project was to demonstrate controlled fracturing and permeability increase to enhance resource production and power generation at Coso. It was planned to stimulate at least two intervals in the Well 73-18RD with zonal isolation using advanced high-temperature packer technology. The anticipation was that a hydraulic fracture would form in one zone and in the second zone natural fractures would reopen, shear and extend, thereby increasing fluid communication in the field. In both cases, the intent was to carefully design the stimulation to take advantage of the relative roles of hydraulic pressure, thermal stress and chemical (“hydro-thermo-chemical”) dissolution. The project had planned to use the Welltec Annular Barrier (WAB) system developed as part of another DOE-funded project which was on schedule to be ready for field deployment. If successful, the Coso Operating Company would incorporate this approach into their ongoing well maintenance in other wells increasing the total MWs produced. The project had multiple tasks and the objective of Task 1 was to develop a workover plan, to build a geologic model, to design a stimulation program, and to work with the operating and services companies to develop a detailed budget and obtain the necessary permits and sundry notices. In Phase 2 the workover plan would be implemented, the designed stimulation would be carried out, and subsequently the well would logged and tested to measure the impact of stimulation. Per sponsor request, a Pre-Task 1 activity with a Go/No Go decision point was scheduled in Task 1.

02 PETROLEUM↗
DOE OSTI2025

Geochemical and Hydromechanical Stimulation for Reaction Acceleration (GeoHydRA) of Serpentinization for In Situ Hydrogen Production

The project team aims to develop an innovative approach to stimulate geological hydrogen generation by enhancing in-situ serpentinization of ultramafic rocks. This method can maximize hydrogen production by promoting reactive surface areas through controlled fracturing and optimizing fluid chemistry and flow rates. Key objectives include establishing optimal conditions for hydrogen generation, investigating subsurface limitations, and sustaining production via hierarchical crack systems.

08 HYDROGEN↗
DOE OSTI2021

Deep reinforcement learning control of hydraulic fracturing

Hydraulic fracturing is a technique to extract oil and gas from shale formations, and obtaining a uniform proppant concentration along the fracture is key to its productivity. Recently, various model predictive control schemes have been proposed to achieve this objective. But such controllers require an accurate and computationally efficient model which is difficult to obtain given the complexity of the process and uncertainties in the rock formation properties. In this article, we design a model-free data-based reinforcement learning controller which learns an optimal control policy through interactions with the process. Deep reinforcement learning (DRL) controller is based on the Deep Deterministic Policy Gradient algorithm that combines Deep-Q-network with actor-critic framework. In addition, we utilize dimensionality reduction and transfer learning to quicken the learning process. We show that the controller learns an optimal policy to obtain uniform proppant concentration despite the complex nature of the process while satisfying various input constraints.

42 ENGINEERING↗
DOE OSTI2026

3D seismic imaging of a fracture damage zone controlling reservoir compartmentalization at the Raft River EGS using multi-azimuth walkaway VSP

Accurate imaging of steeply dipping fracture zones in crystalline enhanced geothermal systems (EGS) is critical for constraining permeability architecture and guiding stimulation design. However, such structures remain poorly resolved by conventional surface seismic methods. We present a fully three-dimensional (3D) elastic-waveform inversion-migration workflow applied to multi-azimuth walkaway vertical seismic profiling (VSP) data acquired in a deviated borehole at the Raft River EGS. The workflow integrates first-arrival traveltime tomography, multi-scale elastic waveform inversion (EWI), and elastic least-squares reverse-time migration (ELSRTM) to recover high-resolution compressional-and shear-wave velocity models and to image structural discontinuities in the crystalline basement. The results reveal a laterally continuous low-velocity anomaly, with shear-wave velocity reductions of 25-30%, consistent with fractureinduced mechanical weakening. Two steeply dipping discontinuities bound a 50-80 m wide fracture damage zone. Independent constraints from microseismic clustering and geochemical compartmentalization corroborate the geometry and structural significance of this feature. Synthetic modeling further confirms that structures of this scale are resolvable (~30 m). These findings indicate that the Narrows structure is a distributed fracture damage zone rather than a discrete fault plane. By resolving fracture-zone geometry at the tens-of-meters scale using a single borehole, this workflow provides a practical and transferable approach for improving structural characterization, reducing uncertainty in permeability architecture, and supporting reservoir modeling and stimulation design in fractured crystalline EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2023

Gravity-driven controls on fluid and carbonate precipitation distributions in fractures

Abstract Many challenges related to carbon-dioxide ( $$\hbox {CO}_2$$ CO 2 ) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90 $$^\circ$$ ∘ (vertical fracture plane) to 30 $$^\circ$$ ∘ . The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0 $$^\circ$$ ∘ ). However, for fracture inclinations greater than 10 $$^\circ$$ ∘ , the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester $$\hbox {CO}_2$$ CO 2 through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Planarizing Spalled GaAs(100) Surfaces by MOVPE Growth

III-V photovoltaic devices have demonstrated exceptional performance across various applications, with controlled crystal fracturing, known as controlled spalling, emerging as a promising method to reduce costs by enabling substrate reuse. Spalling GaAs(100) substrates, a commonly used substrate in III-V photovoltaics, results in faceted ridges that must be planarized to grow high-quality photovoltaic devices. Here, in this study, we demonstrate that a GaAs(100) wafer offcut toward [$0\bar{11}$] and spalled toward [$011$] can be efficiently planarized by growing C:GaAs by metal-organic vapor phase epitaxy (MOVPE) on the surface, with up to 95% of the nominally deposited material used to fill the valleys between ridges. We find that reducing the offcut to 2° enhances the planarizing capability of C:GaAs. A surface morphology model indicates that the density of surface dangling bonds significantly influences the growth evolution of undoped GaAs surfaces. In contrast, the model suggests that the effectiveness of C:GaAs as a smoothing layer stems from modifying the atomic surface structure and, consequently, the associated sticking coefficients of the facets, which can alter the evolution of surface morphology. Our findings provide guidelines for the epitaxial planarization of semiconductor surfaces and improve the understanding of MOVPE growth on nonplanar surfaces.

14 SOLAR ENERGY↗
DOE OSTI2021

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↗
DOE OSTI2020

Experimental Investigation of Hydraulic Fracturing and Stress Sensitivity of Fracture Permeability Under Changing Polyaxial Stress Conditions

Understanding and predicting fracture propagation and subsequent fluid flow characteristics is critical to geoenergy technologies that engineer and/or utilize favorable geological conditions to store or extract fluids from the subsurface. Fracture permeability decreases nonlinearly with increasing normal stress, but the relationship between shear displacement and fracture permeability is less well understood. We utilize the new Geo–Reservoir Experimental Analogue Technology (GREAT cell), which can apply polyaxial stress states and realistic reservoir temperatures and pressures to cylindrical samples and has the unique capability to alter both the magnitude and orientation of the radial stress field by increments of 11.25° during an experiment. We load synthetic analogue materials and real rock samples to stress conditions representative of 500–1,000 m depth, investigate the hydraulic stimulation process, and then conduct flow experiments while changing the fluid pressure and the orientation of the intermediate and minimum principal stresses. High–resolution circumferential strain measurements combined with fluid pressure data indicate fracture propagation can be both stable (no fluid pressure drop) and unstable (fluid pressure drop). The induced fractures exhibit both opening and shear displacements during their creation and/or during fluid flow with changing radial stress states. Flow tests during radial stress field rotation reveal that fracture normal effective stress has first–order control on fracture permeability but increasing fracture offset can lead to elevated permeabilities at maximum shear stress. Here, the results have implications for our conceptual understanding of fracture propagation as well as fluid flow and deformation around fractures.

58 GEOSCIENCES↗
DOE OSTI2022

Controlling Sustainability of Hydraulic Fracture Permeability in Ductile Shales

Hydraulic fracturing technology, along with horizontal drilling and associated completion technologies, dramatically increased US oil and gas production, by increasing the permeability of tight shales and allowing for cost-effective extraction of trapped hydrocarbons. For efficient and economical production, the targeted shale has to be fractured easily (i.e., good ‘frackability’), but also the permeability of induced fractures have to be sustained during production (good ‘sustainability’). For these reasons, currently, ductile shale with high clay content (>~40%) is difficult to exploit as a hydrocarbon resource, although hydrocarbons can still be found in it. Good frackability, high-TOC, low-clay-content resource shales—the ‘low-hanging fruit’—will be depleted, and we need to develop tools that allow tapping into currently underutilized, clay-rich ductile shales. This research project aims (1) to understand the behavior of fractures in clay-rich, ductile (and sometimes swelling) shales and (2) to begin to develop technologies for efficient and economical production from such shales. With these objectives, we examine time-dependent, coupled mechanical-hydrological behavior of open and proppant-filled fractures within different types of shales. In the preceding project, the focus of the research was to understand the behavior of the fractures and to predict their permeability changes due to fracture closure from shale matrix deformation and proppant embedment. In the current project, we investigate the possibility of chemically manipulating the sustainability of hydraulic fractures in ductile shales—specifically by altering the proppant-embedding behavior—via chemical means, at different stages of oil and gas production.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2025

Towards machine-learning a fully-coupled constitutive model for thermal-hydraulic fracture in geothermal systems: phase I (Final Report)

This project, entitled “Towards machine-learning a fully-coupled constitutive model for thermal-hydraulic fracture in geothermal systems: phase I,” addresses challenges in understanding and controlling subsurface fracture networks, which are crucial for applications like deep geothermal heat mining and deep-crustal minerals/metals/hydrogen extraction. The research focuses on advancing the understanding of coupled thermal-hydro-mechanical-chemical (THMC) processes in geologic materials, particularly under the high temperature and pressure conditions found in the deep crust. This seed grant focused specifically on thermal cracking and the development of new constitutive models. Significant progress was made in both experimental and theoretical domains. To study micro-scale fracture formation, the project demonstrated the ability to create thermal cracking under stress in granite samples using a Paterson Gas-medium Deformation Apparatus.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Numerical study of proppant transport and settling processes in fractures

Reservoir stimulation by creating hydraulically conductive fractures is the key step for enabling enhanced geothermal systems (EGS). The effectiveness of stimulation is significantly influenced by the deposition of proppant inside induced fractures. The transportation and settling of proppant in a propagating fracture is controlled by a multitude of operational and physical parameters, including the fracturing fluid rheology, injection rate, proppant concentration, fracture length/aperture evolution, proppant size/density/shape, etc. A numerical tool that robustly and efficiently accounts for all important attributes can facilitate the design and optimization of reservoir stimulation. This study presents the novel computational tool ELK (ELectrical fracKing) developed for the numerical simulation of proppant-fluid mixture circulation in a fractured geothermal reservoir. We enriched the MOOSE-based PorousFlow module with a suite of equations to consider the fluid-proppant mixture with particle-particle/fluid interactions, which include gravitational settling, particle convection, particle hampering, and strong density and viscosity contrasts. The computational tool is validated by comparing the predicted proppant bed evolution against two different laboratory scale experiments of proppant transport in a fixed aperture channel. Further parameter studies were performed, and the modeling results show that the proppant deposition is determined by the mixing characteristics and settling of the particles from the slurry. Concentration-dependent density and viscosity lead to an inhomogeneous distribution of the proppant, particle collision, and enhanced settling at the bottom of the fractures. Preliminary coupling with dynamic fracture propagation shows promising results and will be further developed to simulate hydraulic stimulation at high fidelity.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Fracture Network Localization Preceding Catastrophic Failure in Triaxial Compression Experiments on Rocks

We quantify the spatial distribution of fracture networks throughout six in situ X-ray tomography triaxial compression experiments on crystalline rocks at confining stresses of 5–35 MPa in order to quantify how fracture development controls the final macroscopic failure of the rock, a process analogous to those that control geohazards such as earthquakes and landslides. Tracking the proportion of the cumulative volume of fractures with volumes >90th percentile to the total fracture volume, ∑ v 90 / v t o t indicates that the fracture networks tend to increase in localization toward these largest fractures for up to 80% of the applied differential stress. The evolution of this metric also matches the evolution of the Gini coefficient, which measures the deviation of a population from uniformity. These results are consistent with observations of localizing low magnitude seismicity before large earthquakes in southern California. In both this analysis and the present work, phases of delocalization interrupt the general increase in localization preceding catastrophic failure, indicating that delocalization does not necessarily indicate a reduction of seismic hazard. However, the proportion of the maximum fracture volume to the total fracture volume does not increase monotonically. Experiments with higher confining stress tend to experience greater localization. To further quantify localization, we compare the geometry of the largest fractures, with volumes >90th percentile, to the best fit plane through these fractures immediately preceding failure. The r 2 scores and the mean distance of the fractures to the plane indicate greater localization in monzonite than in granite. The smaller mean mineral diameter and lower confining stress in the granite experiments may contribute to this result. Tracking these various metrics of localization reveals a close association between macroscopic yielding and the acceleration of fracture network localization. Near yielding, ∑ v 90 / v t o t and the Gini coefficient increase while the mean distance to the final failure plane decreases. Macroscopic yielding thus occurs when the rate of fracture network localization increases.

58 GEOSCIENCES↗
DOE OSTI2026

Aperture and roughness govern iron oxide passivation in olivine fractures during carbon mineralization

Efficient carbon dioxide mineralization in silicate-rich formations depends on reaction rates within fracture networks, yet the role of fracture microstructure remains poorly constrained. We investigated how fracture aperture and surface roughness influence mineralization by reacting synthetic forsteritic olivine fractures with controlled apertures (0.5–1.6 millimeters) and roughness (6–16 micrometers) with supercritical carbon dioxide at 13.7 megapascals and 185 degrees Celsius. Optical microscopy and Raman spectroscopy showed that rougher fracture surfaces promoted iron oxide precipitation, which passivated olivine surfaces, limiting magnesium availability for carbonate formation while removing iron as a reactant. These effects were more pronounced in smaller-aperture fractures. A geochemical model incorporating surface topography and passivation reproduced the observed trends, confirming that increased roughness can reduce reactivity despite greater surface area. Our results demonstrate that fracture microstructure strongly influences iron oxide passivation and carbon dioxide mineralization, with important implications for predicting long-term performance of engineered carbon dioxide storage systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

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↗