SEARCH · Engineering Papers

Results for “Fracturing”

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.

SourceResearch connectionFrom yearThrough year

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

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2019

Femoral Fracture in Primary Soft-Tissue Sarcoma of the Thigh and Groin Treated with Intensity-Modulated Radiation Therapy: Observed versus Expected Risk

Purpose: This study was designed to compare the observed risk of femoral fracture in primary soft-tissue sarcoma (STS) of the thigh/groin treated with intensity-modulated radiation therapy (IMRT) to expected risk calculated using the Princess Margaret Hospital (PMH) nomogram. Methods: Expected femoral fracture risk was calculated by using the PMH nomogram. Cumulative risk of fracture was estimated by using Kaplan–Meier statistics. Prognostic factors were assessed with univariate and multivariate analysis using Cox’s stepwise regression. Results: Between February 2002 and December 2010, 92 consecutive eligible patients were assessed. Median follow-up was 73 months (106 months in surviving patients). IMRT was delivered preoperatively (50 Gy) in 13 (14%) patients and postoperatively in 79 (86%) patients (median dose, 63 Gy; range, 59.4–66.6 Gy). The observed crude risk of fractures was 6.5% compared with 25.6% expected risk from the nomogram; the cumulative risk of fracture using IMRT at 5 years was 6.7% (95% CI 2.8–16.0%). The median time to fracture was 23 months (range, 6.9–88.6). Significant predictors of fracture on univariate analysis were age ≥ 60 years (p = 0.03), tumor location in the anterior thigh (p = 0.008), and periosteal stripping to > 20 cm (p < 0.0001). On multivariate analysis, age ≥ 60 years and periosteal stripping > 20 cm retained significance (p = 0.04 and p = 0.009, respectively). Conclusions: In this study, the cumulative risk of femur fracture in patients treated with IMRT (6.7%) is less than the expected risk using the PMH nomogram (25.6%). Established predictors of femur fracture, such as gender, tumor size, and dose of RT, seem to have less impact on fracture risk when using IMRT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2017

Three-Dimensional Modeling of Fracture Clusters in Geothermal Reservoirs

The objective of this is to develop a 3-D numerical model for simulating mode I, II, and III (tensile, shear, and out-of-plane) propagation of multiple fractures and fracture clusters to accurately predict geothermal reservoir stimulation using the virtual multi-dimensional internal bond (VMIB). Effective development of enhanced geothermal systems can significantly benefit from improved modeling of hydraulic fracturing. In geothermal reservoirs, where the temperature can reach or exceed 350oC, thermal and poro-mechanical processes play an important role in fracture initiation and propagation. In this project hydraulic fracturing of hot subsurface rock mass will be numerically modeled by extending the virtual multiple internal bond theory and implementing it in a finite element code, WARP3D, a three-dimensional finite element code for solid mechanics. The new constitutive model along with the poro-thermoelastic computational algorithms will allow modeling the initiation and propagation of clusters of fractures, and extension of pre-existing fractures. The work will enable the industry to realistically model stimulation of geothermal reservoirs. The project addresses the Geothermal Technologies Office objective of accurately predicting geothermal reservoir stimulation (GTO technology priority item). The project goal will be attained by: (i) development of the VMIB method for application to 3D analysis of fracture clusters; (ii) development of poro- and thermoelastic material sub-routines for use in 3D finite element code WARP3D; (iii) implementation of VMIB and the new material routines in WARP3D to enable simulation of clusters of fractures while accounting for the effects of the pore pressure, thermal stress and inelastic deformation; (iv) simulation of 3D fracture propagation and coalescence and formation of clusters, and comparison with laboratory compression tests; and (v) application of the model to interpretation of injection experiments (planned by our industrial partner) with reference to the impact of the variations in injection rate and temperature, rock properties, and in-situ stress.

15 GEOTHERMAL ENERGY↗

NASA NTRS1980

Fracture modes of high modulus graphite/epoxy angleplied laminates subjected to off-axis tensile loads

Angleplied laminates of high modulus graphite fiber/epoxy were examined in several ply configurations at various tensile loading angles to the zero ply direction to determine the effects of ply orientations on tensile properties, fracture modes, and fracture surface characteristics of the various plies. Experimental results consist of stress-strain data, selected plots, fracture stresses and strains, and scanning electron microscope (SEM) photographs of fracture surfaces. It was found that the stress-strain curves were linear to fracture, and that although fracture surface characteristics for a given fracture mode are similar to those for the same fracture mode in uniaxial specimens, no simple load angle range can be associated with a given fracture mode. It was also concluded that SEM results must be supplemented with ply stress calculations in order to identify ranges of fracture modes occurring as a function of ply orientation with respect to the load direction.

Sinclair, J. H.↗

NASA NTRS1979

Mechanics of fracture - Fundamentals and some recent developments

An overview is presented of the fundamental aspects of and recent developments in fracture mechanics. Reference is made to linear elastic fracture mechanics including the state of stresses and displacements in the vicinity of cracks, effects of crack geometry and orientation on stress intensity factors, energy balance of Griffith, Irwin's stress intensity concept, and linear elastic fracture mechanics testing for fracture toughness. Other aspects of this paper include the non-linear behavior of materials and their influence on fracture mechanics parameters, consideration of viscoelasticity and plasticity, non-linear fracture toughness parameters as C.O.D., R-curve and J-integral, and a non-linear energy method, proposed by Liebowitz. Finite element methods applied to fracture mechanics problems are indicated. Also, consideration has been given to slow crack growth, dynamic effects on K(IC), Sih's criterion for fracture, Lee and Liebowitz's criterion relating crack growth with plastic energy, and applications of fracture mechanics to aircraft design. Suggestions are offered for future research efforts to be undertaken in fracture mechanics.

Liebowitz, H.↗

NASA NTRS1990

Deformation and fracture of thin sheet aluminum-lithium alloys: The effect of cryogenic temperatures

The objective is to characterize the fracture behavior and to define the fracture mechanisms for new Al-Li-Cu alloys, with emphasis on the role of indium additions and cryogenic temperatures. Three alloys were investigated in rolled product form: 2090 baseline and 2090 + indium produced by Reynolds Metals, and commercial AA 2090-T81 produced by Alcoa. The experimental 2090 + In alloy exhibited increases in hardness and ultimate strength, but no change in tensile yield strength, compared to the baseline 2090 composition in the unstretched T6 condition. The reason for this behavior is not understood. Based on hardness and preliminary Kahn Tear fracture experiments, a nominally peak-aged condition was employed for detailed fracture studies. Crack initiation and growth fracture toughness were examined as a function of stress state and microstructure using J(delta a) methods applied to precracked compact tension specimens in the LT orientation. To date, J(delta a) experiments have been limited to 23 C. Alcoa 2090-T81 exhibited the highest toughness regardless of stress state. Fracture was accompanied by extensive delamination associated with high angle grain boundaries normal to the fatigue precrack surface and progressed microscopically by a transgranular shear mechanism. In contrast the two peak-aged Reynolds alloys had lower toughness and fracture was intersubgranular without substantial delamination. The influences of cryogenic temperature, microstructure, boundary precipitate structure, and deformation mode in governing the competing fracture mechanisms will be determined in future experiments. Results contribute to the development of predictive micromechanical models for fracture modes in Al-Li alloys, and to fracture resistant materials.

Wagner, John A.↗

NASA NTRS1992

Fracture mechanics; Proceedings of the 22nd National Symposium, Atlanta, GA, June 26-28, 1990. Vols. 1 & 2

Current research on fracture mechanics is reviewed, focusing on ductile fracture; high-temperature and time-dependent fracture; 3D problems; interface fracture; microstructural aspects of fatigue and fracture; and fracture predictions and applications. Particular attention is given to the determination and comparison of crack resistance curves from wide plates and fracture mechanics specimens; a relationship between R-curves in contained and uncontained yield; the creep crack growth behavior of titanium alloy Ti-6242; a crack growth response in three heat resistant materials at elevated temperature; a crack-surface-contact model for determining effective-stress-intensity factors; interfacial dislocations in anisotropic bimaterials; an effect of intergranular crack branching on fracture toughness evaluation; the fracture toughness behavior of exservice chromium-molybdenum steels; the application of fracture mechanics to assess the significance of proof loading; and a load ratio method for estimating crack extension.

Ernst, Hugo A.↗

NASA NTRS2004

Tensile Creep Fracture of Polycrystalline Near-Stoichiometric NiAl

Tensile creep fracture behavior of polycrystalline near-stoichiometric NiAl has been studied between 700 and 1200 K under initial applied stresses varying between 10 and 200 MPa. The stress exponent for fracture varied between 5.0 and 10.7 while the activation energy for fracture was 250 +/- 22 kJ/mol. The fracture life was inversely proportional to the secondary creep rate in accordance with the Monkman-Grant relation although there was extensive scatter in the data. This observation suggests that the fracture life for near-stoichiometric NiAl was influenced by creep under these stress and temperature conditions. Several different fracture morphologies were observed. Transgranular ductile cleavage fracture occurs at 700 K and at the higher stresses at 800 K. The fracture mode transitions to transgranular creep fracture at 900 and 1000 K and at lower stresses at 800 K, while plastic rupture and grain boundary cavitation occur at 1100 and 1200 K. An experimental fracture mechanism map is constructed for near-stoichiometric NiAl.

Raj, Sai V.↗

DOE OSTI2025

Utah FORGE: Well 16B(78)-32 Drill Core Fracture Analysis Images and Data

This dataset contains drilling core data from well 16B(78)-32, including PDF documents with flattened core images annotated by feature type and core interval, as well as spreadsheets detailing feature morphologies by depth, planar feature measurements, and planar feature orientations rotated to in situ conditions. Core was recovered from three intervals, one per stimulation stage, in the crystalline rocks affected by the stimulation of well 16A(78)-32. Seven core runs were conducted, yielding 135.8 feet of recovered core. Features in the core were categorized into planar fractures, semi-planar fractures, unbroken mineralized fractures, rough fractures, curviplanar fractures, concave-convex surfaces, and planar compositional features such as mylonite or dike-like structures. Planar features were measured while the core was positioned horizontally, with the core axis aligned to a downhole azimuth of 42 degrees. Planar core measurements from stimulations 2 and 3 that could be confidently correlated with FMI data were rotated to in situ orientations. This was done by rotating the planes along vertical and horizontal axes to match the azimuth and inclination data recorded in the directional survey of well 16B(78)-32, as well as applying an axial rotation to resemble the fracture orientations observed in the FMI log at corresponding depths. Coherent sets of planar fracture measurements were made by aligning the core within each 3-foot section of the dissected core barrel, and between adjacent 3-foot sections within a core run by matching rock fabrics, saw cuts and/or tool marks. Where coherent fracture measurements could not be made within a core run, data sets are denoted by a subscript (i.e. 2-Ta and 2-Tb both come from tangent core run number 2).

15 GEOTHERMAL ENERGY↗

DOE OSTI2016

High-stress triaxial direct-shear fracturing of Utica shale and in situ X-ray microtomography with permeability measurement

The challenge of characterizing subsurface fluid flow has motivated extensive laboratory studies, yet fluid flow through rock specimens in which fractures are created and maintained at high-stress conditions remains underinvestigated at this time. The studies of this type that do exist do not include in situ fracture geometry measurements acquired at stressed conditions, which would be beneficial for interpreting the flow behavior. Therefore, this study investigates the apparent permeability induced by direct-shear fracture stimulation through Utica shale (a shale gas resource and potential caprock material) at high triaxial stress confinement and for the first time relates these values to simultaneously acquired in situ X-ray radiography and microtomography images. Change in fracture geometry and apparent permeability was also investigated at additional reduced stress states. Finite element and combined finite-discrete element modeling were used to evaluate the in situ observed fracturing process. Results from this study indicate that the increase in apparent permeability through fractures created at high-stress (22.2 MPa) was minimal relative to the intact rock (less than 1 order of magnitude increase), while fractures created at low stress (3.4 MPa) were significantly more permeable (2 to 4 orders of magnitude increase). Here, this study demonstrates the benefit of in situ X-ray observation coupled with apparent permeability measurement to analyze fracture creation in the subsurface. Our results show that the permeability induced by fractures through shale at high stress can be minor and therefore favorable in application to CO 2 sequestration caprock integrity but unfavorable for hydrocarbon recovery from unconventional reservoirs.

58 GEOSCIENCES↗

DOE OSTI2022

Wellbore cement fracture permeability as a function of confining stress and pore pressure

Cemented annulus fractures are a major leakage path in a wellbore system, and their permeability plays an important role in the behavior of fluid flow through a leaky wellbore. The permeability of these fractures is affected by changing conditions including the external stresses acting on the fracture and the fluid pressure within the fracture. Laboratory gas flow experiments were conducted in a triaxial cell to evaluate the permeability of a wellbore cement fracture under a wide range of confining stress and pore pressure conditions. For the first time, an effective stress law that considers the simultaneous effect of confining stress and pore pressure was defined for the wellbore cement fracture permeability. Here the results showed that the effective stress coefficient (λ) for permeability increased linearly with the Terzaghi effective stress ( -p) with an average of λ = 1 in the range of applied pressures. The relationship between the effective stress and fracture permeability was examined using two physical-based models widely used for rock fractures. The results from the experimental work were incorporated into numerical simulations to estimate the impact of effective stress on the interpreted hydraulic aperture and leakage behavior through a fractured annular cement. Accounting for effective stress-dependent permeability through the wellbore length significantly increased the leakage rate at the wellhead compared with the assumption of a constant cemented annulus permeability.

58 GEOSCIENCES↗

DOE OSTI2021

Hydraulic fracture mapping using wavelet-based fusion of wave transmission and emission measurements

Fracture characterization is essential to hydrocarbon and geothermal exploration and production as well as carbon geo-sequestration. Active and passive measurements have been extensively applied to map fractures across various scales. Active and passive methods have inherent advantages and limitations that complement each other. Here, in this study, to improve the mapping of embedded fractures and the surrounding geomechanically altered regions, we integrate active shear-wave transmission measurement with passive acoustic-emission measurement collected during a lab-scale hydraulic-fracturing experiment. The proposed data fusion leverages the advantages of the two modalities of measurements, while minimizing their individual limitations. Polarized shear-wave transmission (active) measurement was collected before and after the hydraulic fracturing, and the acoustic emission (passive) measurement was collected during hydraulic fracturing. Two sets of two-dimensional maps of fracture and fracture-induced damage in axial, median and frontal planes were obtained by separately processing active and passive measurements. These 2D maps are then fed to wavelet-based image fusion technique to integrate the two sources of information for the reliable mapping/imaging of the embedded fractures and the surrounding geomechanically altered regions.

58 GEOSCIENCES↗

DOE OSTI2017

Effect of Unsaturated Flow Modes on Partitioning Dynamics of Gravity-Driven Flow at a Simple Fracture Intersection: Laboratory Study and Three-Dimensional Smoothed Particle Hydrodynamics Simulations

In this work, we study unsaturated gravity-driven free-surface flows on a synthetic surface intersected by a horizontal fracture and the importance of flow modes (droplet; rivulet) on the partitioning behavior at the fracture intersection. We present laboratory experiments, three-dimensional smoothed particle hydrodynamics (SPH) simulations using a heavily parallelized code, and an analytical solution. The flow-rate-dependent mode switching from droplets to rivulets is reproduced by the SPH model, and the transition ranges agree with the laboratory experiments. We show that flow modes heavily influence the bypass behavior. Flows favoring the formation of droplets exhibit a much stronger bypass capacity compared to rivulet flows, where nearly the whole fluid mass is initially stored within the horizontal fracture. This behavior is demonstrated for a multi-inlet laboratory setup where the inlet- specific flow rate is chosen so that either a droplet or rivulet flow persists. The effect of fluid buffering within the horizontal fracture is presented in terms of dimensionless fracture in- flow so that characteristic scaling regimes can be recovered. For both cases (rivulets and droplets), the flow within the horizontal fracture transitions into a Washburn regime until a critical threshold is reached and the bypass efficiency increases. For rivulet flows, the initial filling of the horizontal fracture is described by classical plug flow. Meanwhile, for droplet flows, a size-dependent partitioning behavior is observed, and the filling of the fracture takes longer. For the case of rivulet flow, we provide an analytical solution that demonstrates the existence of classical Washburn flow within the horizontal fracture.

unsaturated flow, fractured media, smoothed partic↗

DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

DOE OSTI2019

Post-forming Room Temperature Brittle Fracture in a High-Strength Low-Alloy Steel Sheet After Various Forming Modes

The effects of various sheet forming modes on the post-forming room temperature fracture behavior of a 550-MPa yield strength high-strength low-alloy (HSLA) steel are investigated. In-plane biaxial stretching, plane strain (cold rolling), uniaxial tension, cylindrical cup drawing, and in-plane compression (IPC) are examined up to a von Mises effective prestrain ε{sub eff} = 0.7. Most of the results pertain to sub-size Charpy-type impact specimens prepared from the prestrained material. After sufficiently large prestrains, the fracture behavior is highly anisotropic. Low-energy fracture modes, namely cleavage and intergranular fracture, occur after all modes of prestraining along plane directions that are perpendicular to the principal compressive prestrain. These planes are parallel to the sheet surface after biaxial stretching and cold rolling and correspond to brittle splits extending into the primary, ductile fracture surface. After cup drawing, the brittle planes are oriented perpendicular to the circumferential compressive prestrain. This results in very low energy fractures propagating in the length direction of a fully drawn cup. After in-plane compression, brittle fracture occurs along planes that are perpendicular to the compression direction. The fracture toughness for a crack propagating along such a plane after IPC to ε{sub eff} = 0.38 was very low (12.7 MPa √m), much lower than in the undeformed condition (247 MPa√m). Changes in grain shape and in crystallographic texture caused by the various prestraining modes, as well as the microstructural damage at non-metallic inclusions and carbides, are examined to try to understand the fracture behavior.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Evolution of planar fractures in limestone: The role of flow rate, mineral heterogeneity and local transport processes

In this study, the processes that affect the evolution of fracture geometry and reactive transport in limestone fractures were evaluated using an integrated experimental and modeling approach. Flow-through experiments were performed to investigate the alteration of three artificially made planar limestone fractures, which were exposed to HCl-acidified water flow (pH 3.8 ± 0.1) at three different rates. Reactive infiltration instability developed in the fractures because of small fracture aperture variations and naturally occurring heterogeneity in the rock, i.e., the presence of clay spots of higher porosity. The dissolution regime changed from conical to dominant and ramified wormhole with increasing flow rate, i.e., Péclet number. Maps of fracture aperture obtained from X-ray micro-tomography imaging, together with measurements of permeability and analyses of fluid chemistry, provided experimental constraints for modeling. Numerical simulations were performed using a 2.5D continuum reactive transport model to examine the development of flow patterns and the reactivity and permeability evolution over time, and to interrogate the impacts of fluid reactivity, mineral heterogeneity and diffusive boundary layer. The modeling results showed that fracture alteration is very sensitive to initial fluid reactivity and mineralogical heterogeneity, in addition to flow rate. Finally, the changes of local transport processes during fracture opening also affected the reactivity of the fractured samples.

58 GEOSCIENCES↗

DOE OSTI2024

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗

DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗