Fracture Growth Research Articles

Deformation, fracture, and fatigue crack growth behavior in TiB-reinforced near-α Ti matrix composites fabricated using powder metallurgy technique

Titanium matrix composites with TA15 alloy as the matrix and TiB (ranging from 0 to 4 vol%) as reinforcement were produced via powder metallurgy-based hot press sintering, employing AlB2 precursor. The addition of TiB to TA15 led to a substantial reduction in the colony and prior β grain sizes, and an increase in the lath thickness in the composites. Plastic deformation due to planar slip, which traverses across the α/β-ligament phase boundaries within each colony, may cause TiB particles located on colony boundaries to break, potentially more so than those within the colony. The composite containing 1 vol% TiB exhibited simultaneous improvements in strength, ductility, and strain hardening rate, attributed to colony refinement and possibly a sizable amount of intra-boundary TiB. Further increase in the TiB content led to significant drops in ductility, which was further exacerbated by the precipitation of α2-Ti3Al phase. A substantial decrease in the mode I fracture toughness, attributed to reduced colony size and diminished crack deflection, was noted in the composite with 1 vol% TiB. The colony size refinement also contributed to a reduction in the threshold stress intensity factor range, ΔK0, for the fatigue crack growth, owing to the decreased tortuosity of the fatigue crack front. This study elucidates the intricate relationship between TiB location, microstructure, and mechanical properties, providing insights for the design of TiB-reinforced titanium alloys with both α and β phases.

The goodness-of-fit of models for fracture trace length distribution: Does the power law provide a good fit?

Fracture trace length distributions are often assumed to follow a power law, which implies that the distribution is scale-independent. The present study tests this assumption by evaluating the goodness-of-fit of three statistical models—the power law, piecewise power law, and lognormal distribution—upon a dataset of 57 trace maps that cover a range of fracture modes, host rock types, network scales, and topologies. The goodness-of-fit was assessed through the unbiased Kolmogorov-Smirnov (KS) test, which accounts for the fitting procedure and the degrees of freedom of each model. The results show that the power law provides a poor fit to trace length distributions, being rejected in 24 trace maps at a significance level of 0.05. In contrast, the piecewise power law and lognormal distribution demonstrated better fits across the fracture networks, with the piecewise power law performing the best overall. The poor fit of the power law can be attributed to mechanical and chemical controls on fracture growth, mainly fracture abutment, as well as stress shadowing and cementation, which affect growth rate at different length scales and result in scale-dependent trace length distributions. The consistent poor fit of the power law across various fracture networks suggests that these controls are prevalent in natural systems. While the power law remains a simple and effective model for trace length distribution, it should be recognized that it overlooks such controls that can influence the mechanical and hydraulic properties of fracture networks. Meanwhile, the fit of the piecewise power-law suggests the existence of a characteristic length where a transition in fracture growth occurs.

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the "sprain," as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ([Formula: see text]) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Fracture toughness and fatigue crack growth in DMLS Co-Cr-Mo alloy: Unraveling the role of scanning strategies

Co-Cr-Mo alloy is crucial for biomedical implants and aerospace components. These parts often exhibit a high level of geometric intricacy. Direct metal laser sintering (DMLS) is ideal for these complex parts. In DMLS, choosing the right scanning strategies is vital, as it significantly affects the fatigue fracture behavior of the printed components. Thus, the present study investigates the effect of different scanning strategies (stripe, meander, and chessboard) on the fracture toughness and fatigue crack growth behavior of DMLS printed Co-Cr-Mo alloy. For each scanning strategy, fatigue crack growth tests have been performed to evaluate the threshold stress intensity factor and Paris law constants. To corroborate the obtained experimental results, microstructure analyses have been performed using electron backscattered diffraction. Further, failure mechanisms have been identified from fractographs obtained using field emission scanning electron microscopy. It is evident from the obtained test results that scanning strategies caused significant variation in fracture toughness and fatigue crack growth behavior. The stripe scanning strategy has exhibited higher resistance to fracture and fatigue crack growth. However, delayed crack initiation has been observed in the case of the chessboard scanning strategy. The present study provide the background for better selection of scanning strategies to mitigate fatigue fracture in DMLS-printed Co-Cr-Mo alloy designed for specific applications.

Thermo-mechanical modelling of spalling around the deposition boreholes in an underground nuclear waste repository during its thermal phase

This paper presents a three-dimensional numerical analysis of multiple fracture growth leading to the development of excavation disturbed zones and spalling around deposition boreholes in a geological disposal facility. The development of fracture patterns is simulated with the Imperial College Geomechanics Toolkit, a finite-element based simulator that can model the simultaneous nucleation, growth, and coalescence of multiple fractures in quasi-brittle rock. In these simulations, fractures develop due to the stress concentrations around the borehole wall, caused by the local in situ stresses, and due to the thermal stresses caused by the radioactive decay of the waste. Fracture patterns, and the extent of the spalled zone, are computed after the borehole drilling, heating, and cooling stages, at the Forsmark repository site in Sweden. The effect of temperature on the nucleation and growth of spalling fractures, as well as on the reactivation of pre-existing fractures, is assessed qualitatively, by comparing fracture patterns, and quantitatively, in terms of the maximum spalling depth, width, and increase in the total fractured surface area. Overall, the simulations presented herein indicate that thermal spalling will increase the depths (away from the borehole) and angular widths of the spalled zone, but is not likely to lead to major increases in fracture aperture, and concomitant increases in hydraulic transmissivity and permeability of the spalled zone, above that which has already been caused by mechanical spalling.

Effect of Low Temperature on the Fatigue Crack Propagation Behavior of Underwater Manned Vehicle Rudder Materials in Arctic Environments.

During polar navigation, adverse environmental conditions like cold temperatures, fatigue, and corrosion can affect surface and underwater manned vehicles (UMVs). Understanding the fatigue fracture growth behavior of polar ship steel is crucial for ensuring safety. This study investigates the mechanical properties and fatigue fracture propagation of steel used in underwater vehicle rudders under various low-temperature conditions through experimental research. It compares and analyzes the static mechanical characteristics, fatigue crack growth rate, and fracture morphology of underwater manned vehicle rudder steels at different low temperatures. Findings show enhancements in yield strength, tensile strength, elastic modulus, and fatigue crack propagation life of steel 925A, steel 20#, and their welded parts under low-temperature conditions. The tensile strength of 925A steel, 20# steel, and their welded parts increases by 6.87%, 14.61%, and 12.55%, respectively, as the temperature decreases from 20 to -60 °C. The yield strength also increases by 14.17%, 29.09%, and 15.76%, respectively. Fatigue crack propagation rate experiments were conducted under different constant low-temperature conditions. This study offers direction for future modeling and experimental testing.

Mechanisms of near-wellbore fracture growth considering the presence of cement sheath microcracks and their implications on wellbore stability

Placement of intended perforations is usually extensive enough into the reservoir to minimize any near wellbore influence on the expected fracture growth. However, some unplanned fractures do not necessarily originate from the intended shot holes. A sufficient volume of published work exists on the processes of crack growth from the near wellbore regions such as from the cement sheath or from fractures observed on the formation wall surrounding the wellbore. A robust quantification of wellbore stability issues should not detach such fracture development, which can potentially result into issues such as cement sheath damage and unnecessary formation damage in the near-wellbore regions. Such fractures may become fluid injection points during stimulation related wellbore pressurization, leading to pronounced development. In such scenarios, fracture development may become a nuisance if it jeopardises the intended structural support and zonal isolation integrity of the cement sheath around the wellbore, leading to wellbore failure, and a potential abandonment. The pressure transmission across this structural seal affects the zonal isolation integrity and its intended structural support leading to wellbore instability. Understanding crack growth behaviour in such cement-rock structural framework is still illusional yet it poses the wellbore to a risky state of instability. Our study provides insight into mechanical and elastic behavior of a tightly bonded cement-rock framework with presence of initial cracks. It also provides insight into how the acting far-field stresses, applied wellbore pressures, crack length, and orientation angle influence the crack growth morphologies including planar fractures, curved and cyclic fractures. The study is important not only for a proper selection of the cement slurry but also for the definition of allowable pressures during the wellbore lifespan. Inappropriately selected parameters might lead to the loss of wellbore integrity to a point as early as the start of a well’s life.

A 3D coupled numerical simulation of energised fracturing with CO2: Impact of CO2 phase on fracturing process

Engineered fractures play a critical role in enhancing energy extraction efficiency. In this study, energised fracturing with CO2, as an alternative approach to conventional water-based hydraulic fracturing, is investigated via numerical simulations. We validated the CO2 finite element-based fracturing model against analytical as well as CO2-fracturing laboratory experiments, then utilised the model to investigate the effects of pressure-temperature dependent properties of CO2 on energised fracturing process. To account for the temperatures expected in a real field, four cases with injection temperatures of CO2 varying between 250 K and 350K, under both isothermal and adiabatic conditions have been considered. In the adiabatic conditions, the temperature variation during compression of CO2 is captured using the Joule-Thompson coefficient, assuming no thermal exchange between the CO2 and the surrounding medium. The results highlight the significant influence of CO2 phase on the fracturing process, during the pressurisation stage, as well as post-breakdown, the speed of fracture growth after the breakdown and subsequent depressurisation and associated cooling of CO2. In the designed cases, the phase-change from gas to liquid or supercritical occurs during the pressurisation and prior to breakdown, while the phase remains unchanged post breakdown and during fracture propagation. Liquid CO2 presents a fast-pressurising process while gaseous CO2 undergoes a lengthy compression stage. Supercritical CO2 is the best performing as the pressurisation is not too lengthy, while the instantaneous post breakdown fracturing is significant. Results show that higher temperature of supercritical CO2 is causing larger instantaneous fracture propagation as it has lower viscosity for the given in situ stresses (>10 MPa).

Effects of fluid salinity and composition on subcritical crack growth in marble

Subcritical fracture growth of rocks in fluid-chemical environments played a crucial role in the weathering process of rocks and the preservation of stone building heritage. Fracture mechanics properties of marble under different fluid compositions and salinity conditions were investigated by carrying out double torsion tests. The subcritical crack growth index (SCI) exhibited different salinity dependencies in two salt solutions. In NaCl solution, the SCI decreased by 1.97 % – 40.50 % with increasing salinity due to the increase in solubility of calcite. The SCI decreased by 4.60 % and 14.01 % respectively as the Na2SO4 salinity increases from 0.01 mol/L to 1 mol/L, and then increased with the increase in salinity, which was related to the relationship between calcite solubility and salinity. Compared to Na2SO4 solution, the weakening effect of NaCl solution on Ka was smaller. In contrast, the fracture toughness (KIC) was insensitive to fluid composition and salinity in the short term, with KIC ranging from 0.65 to 0.79 MPa m1/2 in all salt solutions. The three-dimensional morphology results of the fracture surfaces indicated that the effect of salinity on the roughness of both fracture types was negligible. However, the microstructure of the marble was significantly affected by salinity, and apparent dissolution was observed. The results highlighted the dependence of the subcritical fracture behavior of marble on fluid chemistry, which had implications for the weathering assessment of stone masonry buildings and the permanent protection of stone heritage monuments. In these applications, fluid chemistry could either strengthen or weaken the rock, depending on variations in fluid composition and salinity.

Exploring the Nb-Ti-W-O System for Use as Li-Ion Battery Anodes

While much attention in recent years has been given to high performance Li-ion batteries for electrical vehicles, other applications such as medical devices require a unique set of properties that demand a different set of material considerations such as long lifetimes and extremely high safety. Anodes for such batteries must be carefully selected to properly achieve these goals. The most commonly industrially utilized anode material, graphite, is not suitable for medical applications due to particle fracture and Li dendrite growth at high rates. Currently utilized medical device anode materials obtain this high performance by sacrificing energy density.Herein a new set of Nb based materials for use as anodes are presented that do not need to make such a sacrifice as they display high longevity along with high energy densities. The Nb-Ti-W pseudoternary features many promising single-phase compositions that could be used as anodes for such devices. These compositions derive from crystal structures that are found on both the Nb2O5-WO3 and Nb2O5-TiO2 pseudobinaries1, though they posses a more diverse range of chemical compositions. These compositions feature high performance discharge capacities of up to 300 mAh/g. Some materials also display robust capacity retentions of up to 98% after 10 cycles. The materials were also tested at 37 °C to simulate operation within the human body. At these temperatures, improvements were found in electrochemical performance for nearly all compositions, notably demonstrating a large improvement in capacity retention regardless of crystal structure. References (1) Rehman, S.; Sieffert, J. M.; Lang, C. J.; McCalla, E. NbyW1-yOz and NbxTi1-xOz pseudobinaries as anodes for Li-ion batteries. Electrochimica Acta 2023, 439. DOI: 10.1016/j.electacta.2022.141665.

Effects of Notches on Breakdown Pressures and Fracture Evolution in Hydraulic Fracturing

AbstractSustainable ore extraction in cave mining heavily relies on the effective fragmentation or caveability of the orebody. Since cave mining offers substantial benefits, it has gained popularity after preconditioning was introduced to help improve caveability. Therefore, hydraulic fracturing serves as a vital technique for risk management and cave stimulation. The increased rock competency and high stress levels in the rock mass around the orebody significantly influence fracturing and, thus, cavability processes. In order to improve the çefficiency of preconditioning by hydraulic fracturing to specific parts of the non-caving or poorly caving formations, the use of notches as artificial flaws offers an influence on the directionality of fracture propagation; this approach also has the potential to decrease the necessary breakdown pressures, thereby lifting limitations on the design and mechanical capabilities of fracturing campaigns and reducing required breakdown pressures, which could improve hydraulic fracturing capabilities. In this study, we studied the effects of notches on hydraulic fracturing performance under varying stress conditions. A number of hydraulic fracturing experiments were conducted using different notch quantities and spacings. Notches were created parallel to the axis of confining stress, and specimens were then subjected to constant axial loads of 40 MPa under varying confining pressures ranging from 5 to 40 MPa. A supplementary 3-D discrete element method using 3DEC was performed, and the results were compared with the hydraulic fracturing experiment. The 3DEC models incorporated Darcy's Law to describe fluid flow through fractures, and the Mohr–Coulomb softening yield criterion was used to simulate failures on predefined surfaces, providing a thorough hydro-mechanical coupled solution. We found that introducing notches can effectively reduce the pressures needed for fracture initiation and growth. Moreover, with appropriate spacing, fracture direction can be controlled. This knowledge, combined with the use of numerical modelling, has advanced our understanding of fracture behaviour and the influence of notches on propagation paths under different stress regimes. These findings could potentially revolutionise the field of hydraulic fracturing, making it more efficient and sustainable.

Local quadratic convergence of the SQP method for an optimal control problem governed by a regularized fracture propagation model

We prove local quadratic convergence of the sequential quadratic programming (SQP) method for an optimal control problem of tracking type governed by one time step of the Euler-Lagrange equation of a time discrete regularized fracture or damage energy minimization problem. This lower-level energy minimization problem % describing the fracture process contains a penalization term for violation of the irreversibility condition in the fracture growth process and a viscous regularization term. Conditions on the latter, corresponding to a time step restriction, guarantee strict convexity and thus unique solvability of the Euler Lagrange equations. Nonetheless, these are quasilinear and the control problem is nonconvex. For the convergence proof with $L^\infty$ localization of the SQP-method, we follow the approach from \cite{Troeltzsch:1999}, utilizing strong regularity of generalized equations and arguments from \cite{HoppeNeitzel:2020} for $L^2$-localization. 2020 Mathematics Subject Classification 90C55, 49M41, 49M15.

The Fluvial Battering Ram: Collisional Experiments Reveal the Importance of Particle Impact Energies on Bedrock Erosional Efficiency

AbstractThe battering of bedrock by bedload collisions is the primary mechanism by which bedrock rivers erode and landscapes evolve. The energy imparted via impacts acts to detach bedrock via the growth and intersection of surface fractures. We present impact experiments designed to test the influence of particle impact energy on bedrock erosion rates. We found that erosional efficiency increased with increasing impact energy. Notably, these increases in efficiency are not captured by a widely‐used mechanistic bedrock erosion model. Observed increases in erosional efficiency were linked with enhanced elastic energy dissipation captured by differences in the coefficient of restitution. We suggest that this increase in energy dissipation is indicative of enhanced crack extension for high velocity impacts. Our experiments indicate a clear energy‐dependence for bedrock detachment processes that is not yet captured by bedrock incision models but may be integrated into long‐term erosion rates and landscape evolution.

The Potential for Fracture Growth in Stepped Subglacial Topography as a Quarrying Mechanism

AbstractUnderstanding the rates and mechanisms of erosion by subglacial quarrying is a major unsolved problem in geomorphology. Stress enhancement due to load concentration on bedrock ledges between cavities is hypothesized to drive the growth of fractures. Prior work assumed the formation of vertically oriented tensile fractures at the downstream margins of cavities as the controlling process, but did not account for the evolution of the stress field as fractures lengthen, and in particular the dominance of the shearing mode at fracture tips. We used 2D finite element analysis and J‐integral methods to analyze stress intensity factors and fracture growth potentials at the tips of preexisting fractures in loaded bedrock steps, taking into account normal and shear components and measured rock strengths. By examining different step heights, step riser angles, rock types, prior fracture locations and orientations, and extents of ice‐rock contact zones, we identified some situations favorable for fracture growth, especially in brittle rock types. Typically, however, the growth direction will not be vertically downward but angled up‐glacier away from the step riser, a situation unfavorable for quarrying. Moreover, in many situations, the normal stress across fracture planes will be compressive. Non‐vertical step risers buttress the bedrock and also suppress fracture growth. In contrast, reducing the sizes of ice‐rock contact zones not only increases the loading magnitude, as previously recognized, but also increases intensification of tensile stress at the tips of fractures located just up‐glacier. Thus, larger cavities, and hence, fast sliding and low effective pressures, favor quarrying more strongly than previously recognized.

Evaluation of the combined influence of geological layer property and in-situ stresses on fracture height growth for layered formations

Fracture geometry is important when stimulating low-permeability reservoirs for natural gas or oil production. The geological layer (GL) properties and contrasts in in-situ stress are the two most important parameters for determination of the vertical fracture growth extent and containment in layered rocks. However, the method for assessing the cumulative impact on growth in height remains ambiguous. In this research, a 3D model based on the cohesive zone method is used to simulate the evolution of hydraulic fracture (HF) height in layered reservoirs. The model incorporates fluid flow and elastic deformation, considering the friction between the contacting fracture surfaces and the interaction between fracture components. First, an analytical solution that was readily available was used to validate the model. Afterwards, a quantitative analysis was performed on the combined impacts of the layer interface strength, coefficient of interlayer stress difference, and coefficient of vertical stress difference. The results indicate that the observed fracture height geometries can be categorized into three distinct regions within the parametric space: blunted fracture, crossed fracture, and T-shaped fracture. Furthermore, the results explained the formation mechanism of the low fracture height in the deep shale reservoir of the Sichuan Basin, China, as well as the distinction between fracture network patterns in mid-depth and deep shale reservoirs.

Utilizing well-reservoir pseudo-connections for multi-stage hydraulic fracturing modeling in tight gas saturated formations

Purpose. Research is aimed at integrating multi-stage hydraulic fracturing in horizontal wells with hydrodynamic simulation as a mandatory part of planning the mining of any shale oil or gas reservoir. Methods. Geological and hydrodynamic reservoir modeling is part of the research. The properties and geometries of the hydraulic fracture network and its representation in the dynamic reservoir model were assessed. The comparative characterization was carried out based on the two methods of fracture modeling: cell dimension reduction for explicit fracture modeling (LGR – local grid refinement) and implicit fracture modeling method, presented in this paper, with additional pseudo-connections between well and reservoir. Findings. A hydrodynamic model for low-permeable reservoir, produced by horizontal well, hydraulically fractured with 5 stages, has been generated. This model is calibrated to the production history and flowing bottom hole pressure by applying two methods of fracture modeling. Modeling results show that it is possible to replicate historical well production by using both methods. However, the proposed method with pseudo connections has several advantages compared to the generally accepted, local grid refinement (LGR) method. Originality. For the first time, a system of pseudo connections between well and reservoir was constructed to model a multi-stage hydraulic fracturing for a hydrodynamic model of tight reservoir. Hydrodynamic simulation results were refined and calibrated to the history of hydrocarbon production and flowing bottom hole pressure data using the pseudo-connections and LGR methods. The similarity of the results by applying LGR and pseudo-connections methods was revealed. Practical implications. The use of pseudo connections for hydraulic fracturing modeling can reduce simulation run time for cases where multi-stage hydraulic fracturing has already been carried out or is planned in the future. Additionally, the use of this method allows testing a larger number of realizations and scenarios, including hydraulic fracturing design (number of stages, size and conductivity of resulted fracture systems, fracture orientation, etc.), well placement and fracture growth relative to well trajectory. Also, there is no need to rebuild a model every time for each realization, as is the case with the LGR method.

The interplay between dolomitizing fluids, tectonically-controlled saddle dolomite and calcite cements in lower Cambrian to Furongian strata in the Tazhong Uplift

The interplay between dolomitizing fluids, tectonically-controlled saddle dolomite and calcite cements in Lower-Cambrian to Furongian strata in the Tazhong Uplift, Tarim Basin were interpreted from petrography, geochemistry and fluid inclusions microthermometry. Stochiometric replacement dolomites (RD1, RD2 and RD3) and saddle dolomite (DC) cement were identified. Equant-drusy calcite (ECAL) and blocky-to columnar calcite (LCAL) cements were identified. Weak relationships exist between δ13C, δ18O, 87Sr/86Sr, Mg/Ca versus Mn/Sr ratios, which are tools used to evaluate the impact of diagenetic environment on C–O–Sr isotope composition of the carbonates. The replacement dolomites and DC with high Mg/Ca ratios (0.97–1.21) have low depleted δ18O, δ13C values, compared to the ECAL and LCAL with more depleted δ18O, δ13C values, indicating that the C–O–Sr isotope composition suffered diagenetic alteration by higher-temperature and equilibrium isotopic fractionation. The carbonates fluid system is characterized by negative shift in C–O-isotopic composition with hypersaline signatures compared to Early Paleozoic evaporative seawater. Moreover, the increasing trend in (Cl, Br, I)/(CaO + MgO) ratios with Mg/Ca ratios, couple with positive relationship between Cl and Br, suggests fluids contamination by halogen-rich inclusions. The RD1, RD2 and RD3 were precipitated by relatively higher-temperature saline fluids at depth during shallow-deep burial dolomitization. The DC have average Th-values (141.5–199.2 °C) mostly higher than the estimated ambient temperatures (165 °C) of the studied strata, meaning their growth in fractures resulted from higher-temperature Mg-enriched saline basinal fluids conveyed by hydrothermal convection and short-lived geothermal squeegee-type flow. The LCAL were precipitated by Ca/Mg-enriched basinal brine fluids diluted by influx of meteoric water, cementing fractures in host rocks, as dolomite was replaced by calcite (i.e., calcitization). The conceptual model reveals that dolomitizing fluids and calcitization fluids were conveyed through permeable horizons resulting from strong extensional tectonics and weak compressional tectonics, forming irregular networks of saddle dolomite and calcite cement in host rocks.

Role of Fracture Roughness in Fracture Surface Evolution During Acidized Corefloods of Calcareous Shales

Summary During the hydraulic fracturing process, an acidic hydraulic fracturing fluid (HFF) is injected at high flow rates to break the rock and enhance its flow potential. This rock-fluid interaction induces both physical and chemical alterations on the fracture surface, resulting in the formation of a “reaction-altered zone.” Recent research has revealed that the depth of reaction penetration is minimal, and most changes occur on the fracture surface. To gain a deeper understanding of how fracture roughness affects fracture aperture change, in this work we adopt an experimental approach. Two similar samples of carbonate-rich Wolfcamp shale with calcite-filled fractures are selected. One sample is cut through the center creating a smooth fracture (SF), while the other is fractured by parting along the calcite-filled fracture, generating a rough fracture (RF). The fracture surface topography, mineral distribution, fracture aperture, and rock hardness are characterized before a reactive coreflood using an equilibrated acidic brine is conducted. The pressure drop across the core is measured, and the effluent is periodically collected and analyzed using mass spectrometry. The temporal changes in the fracture surface are observed by conducting physicochemical surface characterization after the coreflood. The results indicate that calcite dissolution is the primary chemical reaction occurring on the fracture surface, weakening it. Furthermore, this dissolution decreases the fracture roughness, which results in fracture closure and ultimately a decrease in the fracture conductivity. The most significant change in the fracture aperture is observed near the inlet. These results highlight the potential impact of fracture roughness on the mechanism of fracture evolution during acidized corefloods. Higher fracture roughness is associated with increased fines migration and a more significant overall change in fracture aperture during injection. This research provides valuable insights into the intricate processes at play during hydraulic fracturing and aids in understanding the dynamics of fracture growth in such conditions.

Effect of swab and surge pressure on the time-dependent wellbore natural fracture development

Drilling engineers continue to grapple with the persistent challenge of maintaining wellbore stability. Throughout the drilling process, the wellbore pressure experiences fluctuations induced by various factors such as swab and surge pressure, leading to potential instability. This study employs a numerical approach, utilizing Abaqus software, to investigate the impact of swab and surge pressure on the natural fracture growth within the wellbore and the evolving pore pressure. A computational tool developed with MATLAB is then utilized to ascertain a safe operational mud window by assessing the time-dependent collapse and fracture pressures. The findings illustrated a notable increase in fracture width as a function of time in response to swab and surge pressure, with the most significant percentage increase reaching 69.92%. Notably, a 69.16% augmentation in fracture width is observed in the immediate vicinity of the wellbore following the application of swab and surge pressure. However, parameters such as fracture length, loss circulation rate, and pore fluid pressure exhibited marginal changes post-integration of swab and surge effects. The examination of the time-dependent wellbore stability after integration of swab and surge pressure indicated a narrowing of the initial mud window as a function of time, with a 14.33% increase in collapse pressure and a 13.80% decrease in fracture pressure. The numerical model verification against analytical solutions demonstrated good agreement, highlighting its potential utility in optimizing particle size for wellbore reinforcement and mitigating lost circulation through natural fracture sealing.

Numerical Modelling of the Influence of Tidal Stresses on Fracture Patterns on the Surface of Europa

AbstractJupiter’s satellite Europa is covered by an ice shell exhibiting many surface features, including linear structures called lineae, which in this work are treated as fractures. A three-dimensional finite-element simulator is used to numerically model fracture nucleation, growth, and interaction, assuming the ice is an isotropic, linear elastic medium. Tidal stresses are exerted upon the ice through the Jovian orbital relationship. These stresses are calculated using a closed-form model derived from first principles. The fractures grow in response to stress concentration around their tips, and a damage criterion models the weakening of the ice matrix. Three-dimensional non-planar multiple fracture growth is modelled as a function of geometric multi-modal stress intensity factors computed at the fracture tips. Fracturing is evaluated over multi-scale periods, from days to millions of years, thus capturing multiple tidal effects. Fracture behaviour is modelled across the Europan surface in one domain. The patterns are dense clusters of lineae about the stress maxima with diffuse fracturing in outlying regions. Fractures are also modelled in the vicinity of subsurface meltwater lenses, where fractures form parallel to the surface in contrast to the usual perpendicular orientation. The resultant fracture patterns are qualitatively compared against images from NASA’s Galileo mission. This work contributes to the understanding of Europan lineae by illustrating how they behave in a fracture mechanics framework, and suggests interesting results regarding lineae interaction with meltwater lenses. This work is also a proof of concept for this modelling approach, and will serve as the framework for future work.