Geophysical modeling in oilfield: fractures wellbores and other fine-scaled problems.
Abstract not provided.
Weiss, Chester Joseph↗
SEARCH · Engineering Papers
Results for “fracture modeling”
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.
Deterministic Discrete Fracture Network (DFN) Model for the EGS Collab Project on the 4850 Level of the Sanford Underground Research Facility (SURF).
Abstract not provided.
Schwering, Paul Conrad↗
Analysis and Modeling of Explosive Fracturing Process in a Transparent Surrogate of Jointed Rock
Explore the source record for details and available documents.
Geosciences↗
Generating Multiphase Fluid Configurations in Fractures using Diffusion Models
Explore the source record for details and available documents.
Computer Science↗
The Role of Chemistry in Fracture Pattern Development and Opportunities to Advance Interpretations of Geological Materials
Abstract Fracture pattern development has been a challenging area of research in the Earth sciences for more than 100 years. Much has been learned about the spatial and temporal complexity inherent to these systems, but severe challenges remain. Future advances will require new approaches. Chemical processes play a larger role in opening‐mode fracture pattern development than has hitherto been appreciated. This review examines relationships between mechanical and geochemical processes that influence the fracture patterns recorded in natural settings. For fractures formed in diagenetic settings (~50 to 200 °C), we review evidence of chemical reactions in fractures and show how a chemical perspective helps solve problems in fracture analysis. We also outline impediments to subsurface pattern measurement and interpretation, assess implications of discoveries in fracture history reconstruction for process‐based models, review models of fracture cementation and chemically assisted fracture growth, and discuss promising paths for future work. To accurately predict the mechanical and fluid flow properties of fracture systems, a processes‐based approach is needed. Progress is possible using observational, experimental, and modeling approaches that view fracture patterns and properties as the result of coupled mechanical and chemical processes. A critical area is reconstructing patterns through time. Such data sets are essential for developing and testing predictive models. Other topics that need work include models of crystal growth and dissolution rates under geological conditions, cement mechanical effects, and subcritical crack propagation. Advances in machine learning and 3‐D imaging present opportunities for a mechanistic understanding of fracture formation and development, enabling prediction of spatial and temporal complexity over geologic timescales. Geophysical research with a chemical perspective is needed to correctly identify and interpret fractures from geophysical measurements during site characterization and monitoring of subsurface engineering activities.
58 GEOSCIENCES↗
Stochastic reconstruction of fracture network pattern using spatial point processes
Fracture spatial patterns can strongly affect fluid flow in the subsurface. Proximity and distribution of fractures control reservoir flow behavior over various length scales. In many studies, however, simplified geometrical patterns are generated for fractures which may lead to unrealistic subsurface models. Here we introduce a new method for characterization and modeling of fracture spatial patterns based on outcrop observations. We use Ripley's K-function to characterize the arrangement of fracture barycenters and intersection points over various length scales. In addition, we use semivariograms to quantify spatial correlation in fracture intensity maps. Using this information, we develop a stochastic algorithm that generates two-dimensional fracture network realizations with spatial properties similar to those of a real fracture network measured in the field. Numerical simulation models indicate that the generated fracture realizations exhibit similar flow behaviors as that of the original fracture network. Furthermore, such modeling tools expand and improve our capability in building representative fracture models and in quantification of uncertainty in naturally and hydraulically fractured reservoirs.
58 GEOSCIENCES↗
Modeling the corrosion of steel casing and the damage of well cement in a borehole system
This work presents the effect of chloride-induced corrosion in an underground borehole system, specifically considering the steel casing and further damage to the surrounding cement. In the past four decades, oil well cement has been successfully used for underground construction. However, high concentrations of chloride ions gradually degrade the functionality and durability of well cement. Gasses and liquids can leak from deterioration caused by chloride ions, and thus motivate this study. The chemo-mechanical coupling diffusion model and a classical fracture mechanics model are applied. The model is validated with OPC concrete corrosion data from literature, and then applied to well cement using experimentally obtained material property data. The time when crack initiation, peak pressure, and complete fracture is reached in the model of the cement sheath provides a quantification of the expected service life of the borehole system.
36 MATERIALS SCIENCE↗
Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation
This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.
42 ENGINEERING↗
Modeling Metallic Fuel using Peridynamics
Based on available modeling and simulation capabilities of peridynamics module in MOOSE framework for oxide fuel, the overall goal of this project is to further develop the peridynamics capabilities for modeling metallic fuel. It includes two major tasks: 1) develop validated scheme to handle the shape tensor singularity due to insufficient active neighbors of a material particle in the peridynamic correspondence model for fracture problems, and 2) develop failure modeling scheme including failure criterion for metallic fuels. Before the peridynamics can be applied to model metallic fuel, the formulation instability of the peridynamic correspondence model should be addressed. The PI first worked on developing new stabilization method to improve the performance of the peridynamic correspondence model and reduce the possibility of getting a singular shape tensor while applying the model for fracture problems. The new stabilization scheme uses bond-associated weight function rather than bond-associated horizon. Compared to bondassociated horizon stabilized method, this new stabilization scheme has better performance with improved prediction accuracy and reduced free surface effect. Using this newly developed stabilization, materials models from BISON can be directly used in peridynamics for metallic fuels, such as fission rate and burnup dependent creep and swell models. Publication of this work is under preparation.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
A multiscale anisotropic polymer network model coupled with phase field fracture
The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.
42 ENGINEERING↗
Mechanisms for the Generation of Complex Fracture Networks: Observations From Slant Core, Analog Models, and Outcrop
We use observations of hydraulic fractures in core, outcrop attributes of natural hydraulic fractures, and analogue models, to address how hydraulic fracture networks evolve. A slant core from the Wolfcamp Formation—an unconventional shale hydrocarbon reservoir in the Permian Basin of West Texas—collected within 18 m and 30 m of two hydraulically stimulated horizontal wells, provided an opportunity to examine hydraulic fractures directly. In approximately 183 m of core, 309 calcite-sealed natural opening-mode fractures and 375 hydraulic fractures were identified. Many hydraulic fractures in the core show complex morphology, including twist-hackle segmentation, diversion, and bifurcation; these structures most commonly develop at lithological bed boundaries and mechanical heterogeneities such as natural fractures and concretions. An outcrop of bed-parallel pavements in the Cretaceous Boquillas Formation in West Texas contains opening-mode fractures that likely formed by natural hydraulic fracturing. Fracture traces provide evidence of twist-hackle segmentation, and are typically associated with bed boundaries and preexisting bed-parallel stylolites. A laboratory study of hydraulic fracturing of 33 synthetic blocks of gypsum and hydrostone revealed fracture steps, diversions, twist hackles, and multiple overlapping fractures together with information on fracture growth directions. These complexities in the fracture network were dominantly nucleated at inclusions used to simulate pre-existing fractures, and as a result of mechanical heterogeneity introduced by the wellbore and perforations. Collectively, our results show that complex fracture networks are produced in hydraulic fracturing of self-sourced reservoir strata. Mechanical stratigraphic boundaries and other heterogeneities are likely to enhance fracture network complexity through the processes of segmentation, diversion, and bifurcation. These processes create multiple fracture strands, resulting in an increased number of hydraulic fractures over those initiated, thereby increasing total fracture surface area. Our study provides insight into hydraulic fracture network propagation, and has applications for evaluation, completion, production, and fracture modeling of unconventional reservoirs.
54 ENVIRONMENTAL SCIENCES↗
A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport
Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.
58 GEOSCIENCES↗
Investigation of Proppant Distributions in Rock Fractures with Applications to Hydraulic Fracturing
Numerical modeling of fluid and proppant transport. Shows impact of transport rate and proppant size on the extent of fracture coverage by proppant.
Mofakham, Amir A.↗
A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation
Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.
36 MATERIALS SCIENCE↗
Phase-field modeling of rate-dependent fluid-driven fracture initiation and propagation
The rate-dependent behavior associated with deformation and fracturing of materials, such as natural rocks, poses significant challenges for modeling. In addition to the complications of the viscoelastic response, the speed of fracture propagation reflects micromechanical mechanisms in the fracture process zone(FPZ). In order to represent these complicated behaviors, a thermodynamically consistent, rate-dependent fracture model is required. Based on rigorous thermodynamic principles, we derive a rate-dependent phase-field mechanical model coupled with single-phase fluid flow in both the matrix and the fracture. The model is guaranteed to satisfy energy conservation during fracture propagation. Here, the system of equations is solved using the introduced solution procedure and a novel preconditioner that accounts for the complex fluid-structure interaction. The proposed phase-field model is tested against several benchmark problems on solid-fluid coupling, fluid-driven fracture propagation and rate-dependent viscoelastic deformation. The model serves as a strong basis for investigating rate-dependent fracturing experiments and for making predictions of material behaviors under new conditions.
42 ENGINEERING↗
Colloid-facilitated transport of 238 Pu, 233 U and 137 Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal
The influence of montmorillonite colloids on the mobility of 238 Pu, 233 U and 137 Cs through a chalk fracture was investigated to assess the transport potential for radioactive waste. Radioisotopes of each element, along with the conservative tracer tritium, were injected in the presence and absence of montmorillonite colloids into a naturally fractured chalk core. In parallel, batch experiments were conducted to obtain experimental sorption coefficients (K d , mL/g) for both montmorillonite colloids and the chalk fracture material. Breakthrough curves were modelled to determine diffusivity and sorption of each radionuclide to the chalk and the colloids under advective conditions. Uranium sorbed sparingly to chalk (log K d = 0.7 ± 0.2) in batch sorption experiments. 233 U(VI) breakthrough was controlled primarily by the matrix diffusion and sorption to chalk (15 and 25% recovery with and without colloids, respectively). Cesium, in contrast, sorbed strongly to both the montmorillonite colloids and chalk (batch log K d = 3.2 ± 0.01 and 3.9 ± 0.01, respectively). The high affinity to chalk and low colloid concentrations overwhelmed any colloidal Cs transport, resulting in very low 137 Cs breakthrough (1.1–5.5% mass recovery). Batch and fracture transport results, and the associated modelling revealed that Pu migrates both as Pu (IV) sorbed to montmorillonite colloids and as dissolved Pu(V) (7% recovery). Transport experiments revealed differences in Pu(IV) and Pu(V) transport behavior that could not be quantified in simple batch experiments but are critical to effectively predict transport behavior of redox-sensitive radionuclides. Finally, a brackish groundwater solution was injected after completion of the fracture flow experiments and resulted in remobilization and recovery of 2.2% of the total sorbed radionuclides which remained in the core from previous experiments. In general, our study demonstrates consistency in sorption behavior between batch and advective fracture transport. Finally, the results suggest that colloid-facilitated radionuclide transport will enhance radionuclide migration in fractured chalk for those radionuclides with exceedingly high affinity for colloids.
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states
Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.
42 ENGINEERING↗