A Hybrid Approach to Rock Pre-conditioning Using Non-explosive Demolition Agents and Hydraulic Stimulation

This study presents a novel approach to rock pre-conditioning to promote the sustainability of low-grade ore mining applications such as in-situ recovery and cave mining. The proposed method involves a two-stage hybrid approach, utilizing soundless cracking demolition agents (SCDAs) to initiate radial fractures in a predrilled host rock, followed by hydraulic stimulation to extend the fractures. SCDA injection in the first stage creates multiple radial fractures around the injection well. However, the extent of fractures is limited to the near vicinity of the injection well. To overcome this limitation, the second stage involves the application of hydraulic stimulation to extend the initiated fractures, which produces a greater fracture density compared to pure hydraulic stimulation. The concept was assessed using a fully coupled hydro-mechanical discrete element model that simulated the hybrid fracturing method on crystalline rock at the grain scale. The results indicate that the proposed method can create a high density of fractures around the injection well. Additionally, we identify and evaluate the key factors affecting the performance of the proposed method, including rock mass heterogeneity, stress anisotropy, and pre-existing defects, providing valuable insights for further experimental design and execution. Overall, the study offers promising results for a potential solution to enhance the efficiency of low-grade ore mining through the hybrid rock pre-conditioning method. A novel hybrid two-stage pre-conditioning technique for low-grade ore recovery using soundless cracking demolition agents and hydraulic stimulation. A fully coupled hydro-mechanical model in the grain scale using the discrete element method to simulate the hybrid pre-conditioning technique. The hybrid method produces multiple radial fractures around an injection well, regardless of in-situ stress anisotropy. The hybrid method reduces the risk of uncontrolled fracture propagation by lowering the breakdown pressure of hydraulic stimulation. A novel hybrid two-stage pre-conditioning technique for low-grade ore recovery using soundless cracking demolition agents and hydraulic stimulation. A fully coupled hydro-mechanical model in the grain scale using the discrete element method to simulate the hybrid pre-conditioning technique. The hybrid method produces multiple radial fractures around an injection well, regardless of in-situ stress anisotropy. The hybrid method reduces the risk of uncontrolled fracture propagation by lowering the breakdown pressure of hydraulic stimulation.

• A novel hybrid two-stage pre-conditioning technique for low-grade ore recovery using soundless cracking demolition agents and hydraulic stimulation.• A fully coupled hydro-mechanical model in the grain scale using the discrete element method to simulate the hybrid pre-conditioning technique.• The hybrid method produces multiple radial fractures around an injection well, regardless of in-situ stress anisotropy.• The hybrid method reduces the risk of uncontrolled fracture propagation by lowering the breakdown pressure of hydraulic stimulation.

Introduction
The global decline in ore grades (Haque 2022) has led to an increased interest in mining technologies for processing low-grade ore.Cave mining also referred to as block caving (Morrison 2022) and In-Situ Recovery (ISR) (Wang et al. 2022b) have been identified as cost-effective methods to exploit low-grade ore.The widescale industry adaptation of these technologies relies on the general applicability of the techniques in different ore bodies.Fluid-based preconditioning or hydraulic fracturing is being explored as a potential technique to increase the range of applications of low-grade mining methods.

Hydraulic Stimulation in Cave Mining
The use of hydraulic stimulation (HS) for rock mass pre-conditioning in cave mining applications has seen an increase in recent years following its success in the oil and gas industry (He et al. 2016).Hydraulic stress stimulation introduces additional joint sets in a target rock formation which improves the fragmentation of a target ore body.Although cave mining is typically adopted in formations that are sufficiently weak to unravel under gravity loading, hydraulic stimulation methods have the potential to expand the application horizon of cave mining by artificially weakening (preconditioning) a target rock formation.Several hydraulic stimulation projects were performed over the years to improve block caving which includes Northparkes mines copper/gold deposits (Jeffrey et al. 2009) and Newcrest Cadia East Mine (Amorer et al. 2022;Kaiser et al. 2013) in Australia, Tashan coal mine (Huang et al. 2018) and Anshan iron mine (Li et al. 2020) in China, Berezovskaya mine (Lekontsev and Sazhin 2008) in Russia, Salvador copper mine (Chacon et al. 2004) and El Teniente Mine in Morales et al. (2007) in Chile.These case studies indicate that hydraulic stimulation is innately dependent on the in-situ stress condition.The resultant fracture due to the hydraulic loading of an injection well propagates in the major principal stress direction (Chen et al. 2018).In addition, the fracturing fluid viscosity and the injection rate are two other parameters governing the fracture pattern in HS pre-conditioning applications (Bunger and Lecampion 2017).In the presence of natural fractures, HS was found to increase the number of fractures of an existing fracture network, whereas, for cave mining, it is favorable to induce additional fractures that intersect natural fractures leading to a blocky ore body.

Hydraulic Stimulation in ISR Applications
Hydraulic stimulation has been considered a potential preconditioning method to improve host-rock permeability in ISR.ISR is a non-invasive mining method that utilizes the intrinsic permeability of an ore deposit to inject a lixiviant that dissolves and extracts target minerals with minimal ground disturbance.In ISR, a series of injection and extraction wells are drilled into the target formation for the lixiviant circulation.HS was also adopted for artificial permeability enhancement in some ISR mining applications.Field trials of ISR projects included Mount Isa and Gunpowder copper mines in Australia (Seredkin et al. 2016), Kimbley pit (Ahlness and Pojar 1983), Van Dyke Copper mine (Huff et al. 1988) in the USA, and Chuqui Sur Deposit (Pallauta 1985) in Chile.However, these trials lack an in-depth understanding of the mechanics associated with HS pre-conditioning for mining applications.

The Gap
According to Bunger and Lecampion (2017), the four critical issues that require addressing for conventional HS applications in oil and gas reservoirs are: (1) controlling fracture height growth, (2) reducing fracture tortuosity near-wellbore, (3) predicting/modelling localized growth geometry, and (4) stimulation of multiple fractures simultaneously.From, a mining perspective, HS pre-conditioning is much smaller in scale compared to the oil and gas industry.For instance, the injection volume, rate and fracture extent for mining applications vary in the range of 8-20 m 3 , 5-10 l/s and 1-2.5 m respectively while in the oil and gas industry, a typical well requires about 1000 m 3 fluid at a rate of 250 l/s, and the fracture extent is in the order of 100 s of meters (Adams and Rowe 2013;Bunger et al. 2012).At smaller-scale mining applications, issues (3) and (4) take precedence.He et al. (2017) and Zeeb and Konietzky (2015) proposed several options to overcome some aspects such as stress shadowing and producing prescribed fractures in a formation using directional hydraulic stimulation and hydraulic stimulation with perforations.However, technologies for pre-conditioning methods for mining applications require further research.
In this paper, we present a conceptual hybrid pre-conditioning method that addresses the issue of inducing multiple fractures.A two-stage fracturing process is proposed that utilizes a soundless cracking demolition agent (SCDA) (De Silva et al. 2018a) to initiate fractures in an injection well, followed by HS to simultaneously extend the initiated fractures.First, a grain scale numerical model using the 3 Dimensional Discrete Element Modelling software 3DEC 7.0 (Itasca 2022) is presented considering the material heterogeneity of a target rock mass.This is used to simulate SCDA-charged fracturing.Next, a fully coupled hydromechanical model simulating hydraulic fracturing at the grain scale is presented.This simulates hydraulic fracture extension.Finally, the model is used to assess the fracturing performance of the proposed hybrid fracturing method under varying in-situ stress conditions and micro-mechanical discontinuities.The proposed method shows promising results that can be used as a precursor to design and execute laboratory/field trials in unconventional low-grade ore mining methods such as cave mining, and permeability-enhanced in-situ recovery (ISR).

The Hybrid SCDA-HS Method
The proposed method consists of injecting SCDA slurry into a predrilled borehole to initiate fractures followed by HS to extend the initiated fractures (Fig. 1a).SCDA is a non-explosive demolition agent capable of generating expansive pressure inside a confined space.SCDAs have also been modified, specifically for subsurface applications below the water table which extends its usefulness to pre-conditioning for permeability-enhanced ISR (De Silva et al. 2018a).It mainly consists of calcium oxide (CaO), which forms an injectable slurry when mixed with water.During hydration, the SCDA slurry solidifies and expands volumetrically in an exothermic reaction producing Ca(OH) 2 (Eq.1).The expansive pressure induced by SCDA in a pre-drilled injection well produces radial compressive stresses and tangential tensile stresses around the injection well (over 180 MPa according to experimental results of Natanzi et al. ( 2016)).Fractures are initiated around the surrounding rock when the tensile capacity of the rock is exceeded.After the first crack initiation, additional radial tensile cracks are formed due to the dissipation of strain energy by SCDA expansion (Fig. 1b).However, the extent of the fracture growth is primarily dependent on the in-situ stress conditions (De Silva et al. 2023) and is typically limited to the near vicinity of the injection well.Therefore, practical applications of SCDA charging require a secondary stage of fracture extension.
The hybrid method considers the extension of multiple radial fractures initiated by SCDA charging using HS.During the volumetric expansion process, the SCDA slurry transforms into a solid state.Following the first stage of fracturing, the solidified SCDA may be drilled out from the borehole.During the second stage of fracture extension, the SCDA-charged area may be isolated by a packer system and hydraulically restimulated to extend the initiated radial fractures.This process was numerically simulated in this study to assess its viability as a potential pre-conditioning technique.

Numerical Simulation of the Hybrid SCDA-HS Method
Micro-mechanical properties affecting fracture initiation and propagation must be well understood for the intended applications of block caving and ISR.Therefore, a grain scale model was produced to understand the fracture mechanics in the rock mass that falls below the representative elementary volume, where the fracture initiation and evolution are heavily dependent on material heterogeneity and pre-existing defects (Fig. 2).The numerical model in 3DEC was developed in four steps.(1) An intact rock model was produced and calibrated at the grain scale using Voronoi tessellation (2) The permeability of the grain assembly was matched to that of intact granite matrix permeability by making the grains impermeable and allowing fluid flow only between grain contact.(3) SCDA-charged fracturing was simulated on the calibrated numerical assembly.(4) Hydraulic stimulation was mimicked on the intact grain assembly by considering compatible injection flow rates for the calibrated hydraulic apertures of the model.SCDA-charged fracturing was simulated in this study using the method proposed by De Silva et al. (2023), where the physical expansion of SCDA was simulated using deformable blocks in an injection well.The expansion of the SCDA was capped by either a peak expansive pressure or a maximum strain, depending on which parameter was reached first in the simulation.For HS pre-conditioning, the typical extended leak-off test (XLOT) was adopted.The XLOT method is an extension of the standard leak-off test where the injection flow is sustained post-fracturing to ensure fracture growth (Hamidi and Mortazavi 2014).Therefore, in the simulations conducted, the injection flow rate was maintained beyond the initial breakdown pressure until continuous fracture extension was achieved in the model.

Hydro-mechanical Simulation of the Voronoi Grain Structure
Strathbogie granite was used as the reference material to calibrate the mechanical and hydraulic properties of the numerical grain assembly.Individual grains were modelled as an assembly of discrete Voronoi blocks with minimum and maximum edge lengths of 3 mm and 8 mm respectively.The grain sizes were selected to accurately represent material heterogeneity while being able to simulate fluid flow with limited computational resources and time.Individual grains were internally meshed as finite difference elements consisting of tetrahedral zones with a zone edge length of 2 mm, allowing them to deform as linear elastic material with Young's modulus of E and Poisson's ratio of ν.Since 3DEC only supports parallel plate geometries for fluid flow calculations in joints, grain contacts are parallel in the Voronoi grain assembly.Contacts between grains were modelled using the Coulomb slip model allowing contacts to fail under both shear and tension (Fig. 3a).The shear strength of the joint is governed by contact cohesion c and the friction angle φ.Tensile failure is simulated by defining a tensile cut-off, σ t between contacts.The shear and normal stiffnesses, k s and k n govern the displacement between joints and the joint apertures during mechanical deformation in coupled hydro-mechanical simulations.
A crack is formed between grains when the applied stresses exceed the tensile or shear capacity of the bond.Joint apertures vary depending on the mechanical loading and the fluid flow through grain contacts evolve with the mechanical loading of the model.As shown in Fig. 3b, grain contact apertures vary linearly depending on the applied stress and fluid pressure.The cubic law, given by Eq. ( 2) is used to calculate the flow q, through contact planes in the numerical model.The fluid flow is restricted to grain contacts only, enabling non-linear pressure dependence of joint apertures.The fluid simulations were run under the small-strain condition, which assumes grain positions and grid points are not changed during the simulation.This condition inhibits significant deformations between grain boundaries and allows the application of the cubic law by considering grain boundaries as parallel plates.Therefore, a fully coupled hydro-mechanical model in 3DEC requires meshing to be done internally as flow planes are superimposed on mechanical contact planes.Due to computational limitations, only intergranular fracturing was simulated in this study to achieve results within a practical timeframe.
In Eq. (2), a is the hydraulic aperture, Δp is the fluid pressure drop along the selected boundary, l is the boundary length, and μ is the fluid viscosity.For flow calculations, the hydraulic aperture is assumed to be equal to the mechanical aperture and the average aperture between contact planes is used for flow calculations for non-parallel contact planes.The fluid pressures, p are updated based on net flow into the domain and volume change due to the deformation of grain boundaries and joints according to Eq. (3).Readers are directed to 3DEC documentation for an in-depth description of fluid flow simulations in joints (Itasca 2022).In Eq. (3), p 0 is the initial fluid pressure, K w is the fluid bulk modulus, Q is the net flow rate into the domain, and V is the volume.V = V 0 + ΔV, where ΔV is the volume change in the volume change in the domain at time Δt, and V 0 is the previous volume in the domain.
Upper and lower limits are imposed on the hydraulic aperture, used for the calculation of the fluid flow (Fig. 3b) so that the cubic law stands valid in the model.The upper limit is imposed to avoid turbulent flow regimes and the lower limit is imposed to account for non-zero joint apertures in real life.The residual aperture, a res permits a minute fluid flow through contacts that can be calibrated to simulate the intrinsic matrix permeability of the intact material (i.e., leak-off).

HS Laboratory Experiment
Past research has shown that the direction of hydraulic fracture propagation originating from a central injection well is directly influenced by the initial principal stress state of (3) the rock mass (Zeng et al. 2018).Material heterogeneity also influences the direction of fracture initiation around an injection well where the fractures orient in the direction of planes with the least fracture toughness (Sesetty and Ghassemi 2018).The existence of fissures has also been shown to increase the possibility of crack branching during preconditioning using HS (Taleghani 2010).The fluid viscosity influences the fracturing mode where low-viscosity fracture fluids were found to initiate intergranular (IG) fracturing (Chen et al. 2015).
A laboratory-scale hydraulic stimulation experiment was conducted on Strathbogie granite (see Table 1 for petrographic and mechanical properties of the specimen) to assess the nature of the induced fractures.A 22.5 mm diameter cylindrical test specimen was cored from an undamaged granite block.The specimen ends were finished using a diamond face grinder to produce a level and smooth surfaces with minimum surface irregularities (< 10 µm).A 3-mm hole was drilled into the center of the specimen for fluid injection.A schematic of the experimental setup and the boundary conditions used for hydraulic stimulation is shown in Fig. 4 and water was injected into the drilled borehole at a constant rate to initiate a hydraulic fracture.The readers are directed to Kumari et al. (2018) for a detailed description of the experimental procedure.

Calibration of Mechanical Properties
The heterogeneity of Strathbogie granite was introduced to the model by considering varying mineralogical strength properties of grains and different grain contacts.Table 1 shows the experimentally observed petrographic and mechanical properties of granite.
For numerical simulations, the mineralogy in Table 1 was simplified to quartz, feldspar and mica representing the dominant mineral phases in the numerical assembly.
The mineralogy was randomly assigned to individual grains until the target percentages were achieved by the volume of the material using the Random Distribution of Properties (RDP) method (Wang et al. 2022a).Table 2 shows the volume percentage, calibrated elastic properties of individual grains and the different contact properties used to simulate interparticle bonds in the model.These micro-mechanical parameters for the grains and grain contacts were assigned by considering a Weibull random distribution given by Eq. ( 4).
In Eq. ( 4), x is the Weibull distributed random variable, λ (λ = 1) is the scale parameter and k (k = 7) is the shape parameter.
The numerical model was calibrated by adjusting the micro-mechanical properties until the macro-mechanical responses for the Uniaxial Compressive Strength (UCS) and the Brazilian Disk tensile Strength (BTS) tests of the numerical model matched that of the granite specimens.The BTS σ t was calculated using Eq. ( 5).
In Eq. ( 5), P is the peak diametrical load, and D and T are the diameter and the thickness of the specimen respectively.
A comparison between the experimental results and the numerical simulations of the calibrated model for the mechanical response of granite is shown in Table 3 and Fig. 5.The calibrated model was verified by calculating the fracture toughness K IC of the material (Eq.6) using a Notched Semi-Circular Bending (NSCB) test (Kuruppu et al. 2014).
(4)  In Eq. ( 6), P is the peak bending load, A is the notch length at the center of the specimen, and Y′ is given by Eq. ( 7).

Calibration of Intact Rock Permeability
Following the calibration of the macro-mechanical properties of the rock, the permeability of the intact matrix was calibrated considering different residual joint apertures in the model as previously discussed.The residual joint aperture ( 6) 2A D 2 a res (Fig. 3b) was changed in the model and the matrix permeability was calculated using constant head permeability tests.The permeability, k of the matrix was calculated using Eq. ( 8).
where Q is the cumulative steady-state flow rate observed at the injection point normal to the specimen cross-sectional area A, and ΔP∕H is the pressure gradient.Model dimen- sions and boundary parameters used in this simulation are identical to those used by Kumari et al. (2018) to measure the intrinsic permeability of Strathbogie granite.The model parameters used and the average flowrate q observed in the flow planes at the injection plane for the selected parameters are shown in Table 4.The fluid pressure and the flow rate distributions of the specimen at a steady state are shown in Fig. 6a and b, respectively.The steady state was achieved by running the model until the pore pressures stabilized  within the specimen (Fig. 6c).The cumulative flow rate Q, in this simulation corresponds to the summation of flowrates through flow planes between grain boundaries at the injection face as indicated by red line segments in Fig. 6b ( Q in = ∑ q i l i = qL ).q i , the flowrate per unit width is cal- culated using Eq.(2).A similar methodology was adopted by Tan and Konietzky (2019) to simulate the intact permeability of Granite in 2D. Figure 6d shows the variation of the specimen permeability with grain contact plane aperture a 0 which indicates a non-linear variation.Considering the experimental permeability of intact Strathbogie granite of 1.00 × 10 -19 m 2 at 10 MPa confining pressure (Kumari et al. 2018), a 0 of 2 × 10 -7 m yielded a reasonable result.

SCDA Charging and Hydraulic Fracture Numerical Simulation
SCDA is typically introduced into a rock mass as an injectable slurry through pre-drilled boreholes.Once, injected, the SCDA transforms into a solid state during the hydration process and results in the volumetric expansion explained in Sect.1.4.The elastic properties of SCDA evolve with the hydration of SCDA (De Silva et al. 2017;Xu et al. 2021) and the stiffness of the SCDA influences the expansive pressure exerted on the injection well boundary (Tang et al. 2022).
However, the numerical model in this study simulates SCDA expansion from the point of transition to a solid state, and constant material properties were assumed for simulating the solid SCDA in the numerical model (Table 5).
The expansion of SCDA was modelled by generating a uniform radial distribution of tetrahedral expansive elements within the borehole.Experimental observations of SCDA expansion yield an exponential decay of the strain rate applied by SCDA on the surrounding material (De Silva et al. 2023).This can be translated into a normal velocity boundary in SCDA elements for implementation in 3DEC that follows the decay function given by Eq. ( 9).
where, v is the velocity magnitude, f and b are shape factors governing the rate of expansion in SCDA and t is the total simulation time in 3DEC.SCDA-charged fracture propagation was terminated in the model in one of two ways; either by v asymptotically reaching zero as the simulation progresses (corresponding to the finite strain energy in SCDA expansion) or when the expansive pressure in the model reaches a predefined peak value based on experimental data.In this simulation, for the calibrated rock mass, the peak expansive pressure was limited to 45 MPa, f = 1 and b = 20 in Eq. ( 9).For further details of the simulation method, readers are directed to De Silva et al. (2023).
In the coupled numerical model presented, the grain boundaries can open and slip as a function of effective stress and their bond strength.A fracture is registered between two contacts (grain boundaries) when the tensile or shear strength of the bond is exceeded as a result of block  (Kumari et al. 2017) and diametrical splitting for the BTS specimen (Kumari et al. 2019) deformation and/or pore-pressure build-up in the contact during fluid injection.Hydraulic fracturing is simulated by applying a constant flowrate injection boundary to a predefined injection well in the model (See Sect.4.3).During the fracture simulations, the bonds between grains are allowed to progressively fail representing hydraulic fracture propagation.As the fractures propagate, the pore pressure dissipation within the fluid planes is simulated with the coupled hydro-mechanical simulations (as explained in Eq. 3).Since the grain assembly represents a heterogeneous material, in addition to the initial stress state of the model, the fracture evolution is also dictated by the varying bond strength between different grain contact boundaries.

HS Model Verification in 3DEC
An initial model was run to verify the validity of planar hydraulic fracture propagation in 3DEC.The closed-form theoretical approximation for the fracture aperture and pressure variation along the fracture plane of a radially growing penny-shaped hydraulic fracture (Fig. 7) in an isotropic elastic medium was considered for the validation (Dontsov 2016;Savitski and Detournay 2002).The solution is only valid under the following assumptions.
1. Incompressible Newtonian fluid injection from a point source at the center of the fracture 2. The lag between the fluid front and the fracture tip during fracture propagation is negligible  Under zero fluid leak-off, fracture propagation is categorized as either in toughness-dominated (K) or viscosity-dominated (M) regime (Dontsov 2016).The limiting regime (for C = 0 case) is determined by the dominance of the dissipative mechanism associated with either fluid viscosity or rock toughness during fracture propagation.The M vertex limit corresponds to a viscosity dominated and zero toughness (K IC = 0) condition, indicating zero (tensile and shear) strength in the hydraulic fracture.Therefore, hydraulic fracture propagation in an isotropic elastic medium as shown in Fig. 7 can be described by: where Q o is the injection flow rate, w is the fracture aperture, and p is the pressure gradient along the fracture radius R, at a given time t (See Fig. 7).F is defined as an integral containing a piecewise function, defined as where, ρ = r/R, s = 1−ρ and the piecewise function M(ρ,s) is defined as In Eq. ( 15), K(•) and E(•) are complete elliptic integrals of the first and second kind respectively.For the M vertex solution, ~ 0.487, and ~ 2/3 (Dontsov 2016).The asymptotic solution for Eqs.(11 and 12) is numerically  A 60 mm cube with a horizontal penny-shaped crack with a 10 mm diameter and zero tensile strength and cohesion was modelled with a single injection point with a flow rate Q o = 10 -5 m 3 /s in 3DEC to replicate the analytical solution.Young's modulus, E of 40 MPa, Poisson's ratio, ν of 0.22 and fluid viscosity, μ of 0.001 Pa s was considered for the model.An initial (and minimum) hydraulic aperture of 3 × 10 -6 m and a maximum hydraulic aperture of 10 -4 m was defined for fluid flow through fractures.To simulate zero leak-off, fluid flow was permitted only in failed contacts in this instance.The analytical solution was compared ( 16) against the numerical results at t = 10 -3 s after starting the fluid injection.
Figure 8 shows the initial model and the observed pressure and fracture aperture distribution comparison with the analytical solution at the end of 10 −3 s.The model simulates the outward propagation of the fracture with reasonable accuracy.Although there is a deviation at the crack tip, this can be attributed to the non-zero initial fracture aperture and the finite bulk modulus (compressible fluid) defined in 3DEC.These were assumed to be zero aperture at the crack tip and incompressible fluid in the analytical solution.

The Hybrid SCDA-HS Method in the Granular Assembly
In this study, first SCDA charging was simulated to initiate multiple radial fractures around a central injection well followed by a constant rate fluid injection applied to the boundary nodes of the injection well.The boundary conditions and the model geometry used to simulate the hybrid SCDA-HS pre-conditioning method are shown in Fig. 9.

SCDA Pre-conditioning
Previous research has shown that the in-situ stress conditions (De Silva et al. 2018b;Guo et al. 2015;Zhai et al. 2018) and material properties strongly influence the fracture initiation pattern of SCDA charging (Lei et al. 2021).In this study, the effects of stress anisotropy and increasing confining pressure on the overall fracture initiation, growth and termination in heterogeneous rock were investigated first.Figure 10 shows the SCDA-charged fracture pattern observed for four different stress regimes shown in Table 6.The stress regimes in Table 6 were used for all subsequent simulations.A comparison of Case 1 and Case 3 exhibits the influence of isotropic confining pressure in the XY plane, while Case 2 and Case 3 exhibit the effect of stress anisotropy on SCDA-charged fracturing.Finally, a comparison of Case 2 and Case 4 reveals how the increased intensity of stress anisotropy influences fracture propagation.
The tangential stresses developed around the injection well during SCDA charging are shown in Fig. 10 (compression is considered positive).Figures 10a-d show the variation of the tangential tensile stress developed around the injection well at fracture growth termination for the four stress conditions considered.The results indicate fracture termination within the tangential tensile stress field.It is evident, that at lower confining pressures (Case 1) the tensile stress field reaches further outward from the injection well, contributing to extended radial fractures, whereas increasing the confinement (Case 3) appears to restrain the extension of SCDA-charged fracturing.This phenomenon can be better explained by observing the evolution of the average tangential stress around the injection well (Fig. 10e).The initial tangential stress around the injection well is compressive (positive) due to the in-situ stress state.Next, the tangential tensile stress created by SCDA expansion in the borehole results in the reorientation of stresses around the injection well.In other words, the stress reorientation transforms the near-borehole tangential stress from a compressive state to a tensile state (state 1 in Fig. 10e).Fracture initiation is inhibited until this tensile stress exceeds the tensile strength between mineral grains (state 2, in Fig. 10e) and the initiated fractures continue to grow within the tensile stress field (state 3, in Fig. 10e).
As the in-situ stress increases, the energy spent by SCDA expansion for stress reorientation in state 1 increases.Therefore, a reduction in the energy spent on fracture extension is expected leading to shorter fracture lengths.Furthermore, the in-situ stress anisotropy skews the otherwise circular tensile stress state around the injection well as shown in Fig. 10b and d.However, the overall geometry of the SCDA-charged fracturing remains relatively unaffected by the in-situ stress distribution at lower stress ratios (Case 2).Comparing Fig. 10b and d shows that the increased stress anisotropy in the XY plane further skews the tensile stress field, contributing to fracture extension in the Y direction (higher stress) and crack growth suppression in the X direction.

HS Pre-conditioning
The hydraulic fracturing test was performed at unconfined conditions to achieve the best visibility of the generated fractures.However, a constant axial load of 2 kN was Given the similar densities of feldspars and quartz, distinguishing these minerals from the CT images is difficult.However, a clear demarcation was observed between mica (lighter spots) and the rest of the rock matrix.Comparing mica boundaries in the fracture path, it is evident that the fracture path mainly follows grain boundaries.Therefore, the IG cracking simulated in the numerical assembly presented in this study would yield a realistic result in terms of the crack path.
The limitations of the small-strain conditions and the maximum apertures of the numerical model suggests, it is more suited to simulate confined rock and may be unrealistic to simulate large deformations observed during HS in unconfined conditions.However, a realistic prediction of the breakdown pressure can be obtained through the simulation process presented (Fig. 11b).By simulating identical conditions to the HS experiment, a comparable fracture pattern and a breakdown pressure of 8.39 MPa (8.35% error) were  observed.This error can be attributed to several factors such as the grain size used in the model, fluid leak-off during the experiment, and the compressibility of the fluid.

Effect of Stress Anisotropy and Heterogeneity on HS Fracture Propagation
Figure 12 shows the fracture patterns observed in the granular assemblies for the in-situ stress conditions given in Table 6 at a fluid injection rate of 18 ml/min.Compared to SCDA charging, HS produced a single fracture that reached the specimen boundary, which was driven by the fluid pressure at the fracture tip irrespective of the in-situ stress field.At a stress ratio of 1 (Case 1: σ x = σ y = 1.0 MPa and Case 3: σ x = σ y = 5.0 MPa) in the XY plane, the fracture grew diagonally, which was dictated by the grain arrangement.
Increasing the in-situ stress initiated additional minor cracking around the injection well (from 10,439 cracks in Case 1 to 11,377 cracks in Case 2: σ x = 2.5 MPa, σ y = 5.0 MPa) and increased the tortuosity of the main fracture (Fig. 12a and  c).The Tortuosity was calculated on the horizontal plane by considering the fracture length on the top face of the specimen.The geometric two-dimensional tortuosity was defined as the ratio of fracture length (geodesic length) and the straight line distance between the two ends of the fracture (Euclidean length) (Shanti et al. 2014).Additionally, HS shows a higher sensitivity to the in-situ stress anisotropy compared to SCDA charging.At an XY plane stress ratio of 2 (Case 2), the HS fracture grows in the Y direction with some minor crack branching in the X direction (Fig. 12b), particularly around the injection well.This crack branching is suppressed when the stress ratio was further increased to 5 (Case 4).Furthermore, increasing the stress ratio in the XY plane from 1 to 5 (Fig. 12a, b, and d) exhibits an increase in the shearing effect of the fracture as indicated by the reduction of tensile to shear crack (T/S) ratios from 8.45 to 6.49.However, these additional shear The fluid pressure inside the injection well during the fracturing process for each case is shown in Fig. 12e.As expected, the highest breakdown pressure was recorded for the highest in-situ stress condition (Case 4) and the lowest breakdown pressure was recorded for the lowest in-situ stress state (Case 1).However, the slightly higher breakdown pressures for Case 2 and Case 3 correspond to additional cracking around the injection and the greater overall compressive stress in the X-Y plane.This phenomenon can be alleviated by inducing crack initiation points around the injection well.However, comparing Case 2 and Case 4 indicates that although the breakdown pressures are relatively similar a significantly lower T/S ratio was observed for Case 4, indicating a strong influence of the stress anisotropy on the failure mode during HSs.

The Effect of Material Heterogeneity, Grain Size and Mineral Distribution on HS
Four different Voronoi grain distributions were considered to assess the influence of grain arrangement on fracture initiation and propagation.The grain size distributions  2) and only the spatial distribution of the minerals was randomized.The minimum edge length for meshing the microstructure was kept constant at 1.4 mm for all four assemblies, which eliminated the influence of mesh size on fracture propagation.The four assemblies were hydraulically fractured with an injection rate of 18 ml/min at the central injection well under Case 4 (σ x = 1.0 MPa, σ y = 5.0 MPa) in-situ stress condition.The observed fracture patterns, fluid leak-off into the matrix, and the different grain contacts that failed during the simulation are shown in Fig. 14a-d.Although the overall crack growth is in the Y direction (due to the stress anisotropy), the grain arrangement seems to influence both the fracture path and the extent of fluid leak-off into the rock matrix, despite maintaining identical injection rates and in-situ stress states.This is because the fluid flow during hydraulic stimulation follows the path of least resistance and the energy required to fail contacts depends on the adjoining mineral grain type (Table 2).
Figure 14e shows a comparison between the different contact types as a fraction of total (failed and intact) contacts in the entire assembly (denoted by T) and the contact types as a fraction of total failed contacts (denoted by C) for the four grain assemblies simulated.In other words, T includes both failed and intact contacts and C includes only failed contacts.The two most abundant contact types F-Q and Q-Q in the model show an indication of crack growth in the path of least resistance.As seen in Fig. 14e, an overall increase in the fraction of F-Q contacts and a decrease in the fraction of Q-Q contacts were observed in the failed contacts.This is because F-Q contacts require lower energy to fail than Q-Q contacts, leading to fracture propagation with the least resistance.This relationship was less apparent in the other contact types due to the lower percentage of mica grains in the assemblies.
The fluid pressure development in the injection wells for the four assemblies considered is shown in Fig. 14f.Although the crack path showed dependence on the mineral distribution and the grain arrangement, the breakdown pressure showed little to no variation.This is because the macro-mechanical response of the four grain assemblies remained relatively unchanged since identical mineralogy and comparable grain-size distributions were considered for each of the four grain assemblies.The pressure distribution within the rock matrix due to fluid leak-off is also shown in Fig. 14a-d.This fluid leak-off pressure distribution reflected dependency on grain arrangement and it is fair to assume that pre-existing defects, particularly fissures in granitic rock can also influence the fracture pattern produced by hydraulic stimulation.Therefore, hydraulic stimulation in the presence of grain scale fissures (micro-cracks) was investigated next.Zhou et al. (2020) have shown two types of defects that can develop in granite micro-structure.The first is thermal cracks arising due to thermal expansion and different thermal coefficients between individual grains that form on grain boundaries (IG cracks).The second type is cracking, which forms both at grain boundaries and through individual grains (TG cracks) due to induced mechanical stresses.However, as previously discussed, only pre-existing IG cracks were modelled in this study due to computational limitations.Hydraulic simulations were performed on two defect densities for comparison.The defects were modelled by random selection of grain contact boundaries and setting the tensile strength and cohesion of the contact to zero.The two models produced had a percentage of unbonded contacts (fissures) at 5% and 15%.HS was simulated under stress Case 3 (σ x = 5 MPa, σ y = 5 MPa). Figure 15a and b show the even distribution of the defects in the models.Forming defects around the injection well were deliberately left out as it would substantially alter the fracture initiation and would not form a basis for direct comparison.Therefore, the fluid pressure build-up of the injection well and the breakdown pressure that corresponds to the fracture initiation were similar for the two models considered.However, the fracture propagation and the fluid leak-off following the breakdown point showed dependence on the pre-defined fissures in the model.Increased fissure density (Fig. 15c and d) contributes to fracture extension.At the same time, a greater distribution of fluid pressure was observed due to fluid leak-off in the fissures.Comparing the strain energy dissipation of the two instances two important aspects are revealed.First, the tensile strain energy dissipation in the contacts is relatively unaffected despite the substantial difference in the defect density.Second, the introduction of fissures tends to have a significant impact on the shearing-induced fractures as the shear strain energy and the subsequent joint friction energy (frictional slip after the shear failure of contacts) in failed contacts increased by 18% and 148%, respectively when the fissure density was increased from 5 to 15%.

The Effect of Injection Flow Rate on HS
Three injection rates, 18 ml/min, 42 ml/min and 66 ml/min were considered for the simulations to assess the effect of injection flow rate on HS fracture propagation at the grain scale.The fluid viscosity was kept constant and the in-situ stress Case 2 (σ x = 2.5 MPa, σ y = 5.0 MPa) was used for the simulations.The observed results are shown in Fig. 16.The injection flow rate corresponds to the rate of energy transfer to the system.The increase in fluid injection rate translates to increases in the pressure development in the injection well, the breakdown pressure and the number of fractures around the injection well.Although there is experimental evidence (Zhuang et al. 2022), the reasons for this phenomenon have not conclusively been proven.However, the increased breakdown pressure can be explained as follows.
As the rate of net energy input to the system (viscosity × injection rate) increases under higher injection rates, a greater amount of strain energy is absorbed by the rock.At lower injection rates, fluid pressure acting on weaker grain boundaries leads to fracturing through grain boundaries and the eventual dissipation of the pore pressure around the injection well.However, at sufficiently high injection flow rates, the fluid pressure is not dissipated and initiates additional cracking around the injection well.Similar observations have been made in other studies where an increased rate of energy input contributes to higher failure stresses and increasingly abrupt drops in the injection pressure (Duan et al. 2018;Fallahzadeh et al. 2017;Kim and Changani 2016;Komurlu 2018;Zhuang et al. 2019).

The Hybrid SCDA-HS Pre-conditioning Method
Given the limitations associated with each pre-conditioning method, namely limited fracture growth in SCDA charging and single fracture propagation in HS, the hybrid fracturing method was assessed to improve the fractures generated by pre-conditioning.As explained in Sect.1.4, the injection well was first charged by SCDA to initiate fractures.Next, the injection well was pressurized by applying a fluid pressure boundary condition to the injection well (see Fig. 9 in Sect.4.3) where a fluid injection rate of 90 ml/min was maintained.The hybrid fracturing method was simulated under two conditions.First, under anisotropic stress conditions (Case 2: σ x = 2.5 MPa, σ y = 5.0 MPa) with no defects and second, under stress Case 2 in the presence of 15% fissures in the model, which represent pre-existing defects at the grain scale.
Figure 17a shows the initial fracture pattern produced by SCDA charging for Case 2 and the expansive pressure development is shown in Fig. 17d.It should be noted that the time taken for the SCDA expansion is not representative of the field or experimental conditions (which would be in the order of hours).But rather, the total strain energy released by SCDA was modelled over a shorter period for ease of simulation.As previously shown, SCDA charging produces multiple radial fractures despite the stress anisotropy around the injection well.The fracture extension produced by subsequent fluid injection is shown in Fig. 17b, where fracture extension was observed in five fractures.In comparison, the hydraulic loading of the injection well without prior SCDA charging only produced 3 cracks (Fig. 17c) but also registered a much higher breakdown pressure of 31.2Mpa.
The expansive pressure development in the injection well during hybrid SCDA charging, and a comparison of fluid pressures in the injection well for the HS and hybrid fracturing are also shown in Fig. 17d.Since the fractures around the injection well are already initiated before fluid injection in the hybrid method, the sudden drop in the injection pressure observed in pure HS pre-conditioning is not observed during the hybrid fracturing process.Instead, multiple peaks corresponding to fracture extension and a gradual reduction in the injection pressure were observed.This reduction in the breakdown pressure and the absence of an abrupt drop in injection pressure could further contribute to the mitigation of uncontrolled fracture propagation and possibly induced seismicity that is often seen in pure HS.
The performance of the proposed hybrid fracturing method in the presence of 15% micro-fractures is illustrated in Fig. 18. Figure 18a shows the SCDA-charged fracture pattern observed for the first stage of hybrid fracturing.A comparison of the expansive pressures and the injection pressures developed during hybrid fracturing for 0% and 15% pre-existing fissures is shown in Fig. 18d.The expansive pressure development during SCDA charging was nearly identical for 0% and 15% fissures (Fig. 18d).Yet, the SCDA-charged fracture exhibited an elongation in the major principal stress direction, which may have been caused by the coalescence of fissures in the rock matrix (See Figs.17a and 18a).However, for the fluid injection stage of hybrid fracturing, inferior fracture propagation results were observed in the presence of fissures (Fig. 18b) compared to the case of no fissures (Fig. 17b).This lack of fracture extension during the fluid injection can be attributed to increased fluid leak-off into the matrix with the introduction of fissures.This dispersed fluid pressure distribution observed within the rock matrix during the fluid injection is shown in Fig. 18c.Furthermore, the comparison shown in Fig. 18d (for 0% and 15% fissures) indicates a lag in the injection pressure build-up and a reduction in the peak injection pressure for the specimen with 15% fissures which is indicative of the increased fluid leak-off into the matrix.
Simulating the hybrid SCDA-HS pre-conditioning method in a numerical rock block with equal dimensions, although computationally intensive on commercially available computer hardware, provides insight into the 3-dimensional morphology of the induced fractures along the axis of the injection well.Material heterogeneity at the grain scale appears to further contribute to the tortuosity of the fractures produced.Figure 19a  Comparing the cross sections with the fracture patterns of the rock 'slab' simulations (Fig. 17b), Fig. 19c exhibits a more complex fracture pattern with multiple crack bifurcations in addition to the main radial cracks.The results of the HS stage of the simulation are consistent with laboratory HS experimental observations, where crack bifurcation was also observed near mica grain boundaries (Fig. 19d).This is due to the 3-dimensional effect of crack propagation in the direction of least resistance considering the strength variability between grain contacts as shown by the crack bifurcations at weak grain contacts in the numerical model.This added

Conclusions
In this study, we propose a two-stage hybrid rock pre-conditioning method using soundless cracking demolition agents (SCDA) and hydraulic stimulation (HS).W first compared SCDA charging and HS in isolation to evaluate the limitations of each pre-conditioning method.We found that SCDA charging creates multiple radial fractures around an injection well with limited fracture extent, while HS typically produces a single unidirectional fracture.In the proposed hybrid method, SCDA charging is used to initiate multiple radial fractures, followed by HS to extend the fractures.To simulate the SCDA-HS hybrid method on crystalline rock at the grain scale, we employed a fully-coupled hydro-mechanical discrete element model.Our results indicate that the hybrid SCDA-HS method generates additional radial fractures and results in a significant decrease in the breakdown pressure (~ 30%) which could also help mitigate issues such as uncontrolled fracture propagation and induced seismicity.We also examined other key factors that affect the performance of the hybrid method such as rock mass heterogeneity, in-situ stress anisotropy, and pre-existing defects.The findings from these simulations provide a basis for further experimentation to evaluate the viability of the proposed hybrid pre-conditioning method for applications such as rock mass pre-conditioning for low-grade ore mining.otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http:// creat iveco mmons.org/ licen ses/ by/4.0/.

Fig. 1
Fig. 1 The hybrid SCDA-HS injection method.Stage 1 SCDA injection and stage 2 hydraulic stimulation a section along A-B, and b plan view

Fig. 2
Fig. 2 The dependency of rock mass properties on the scale and the selection of the representative elementary volume (modified after Bear (1988) and Farahmand et al. (2018))

Fig. 3
Fig. 3 Joint contact model a Coulomb slip model between contact planes of two grains, and b hydraulic aperture variation in contact (modified after Zhang and Dontsov 2018)

Fig. 4
Fig. 4 Schematic diagram of the experimental setup for hydraulic stimulation

Fig. 5
Fig. 5 Stress-strain response comparison of the numerical model a UCS, and b BTS.The numerical result is shifted in the x-axis to align with the experimental results as indicated by the arrow.Point strain

8Fig. 6
Fig. 6 Steady state constant head permeability test for intact rock a numerical assembly and the pressure distribution with pressure measurement points PT1-PT3, b flowrate distribution within grain contacts during steady state (inflow and outflow pressure boundaries denoted by red and blue points respectively), c pressure stabiliza-

Fig. 8
Fig. 8 Comparison of analytical and numerical results for a penny-shaped hydraulic fracture propagation at t = 10 -3 s. a Initial numerical model, b hydraulic aperture contour c fluid pressure contour, and the variation of d crack width (aperture) and e pressure along the fracture radius

Fig. 9
Fig. 9 Model boundary conditions, a heterogeneous grain assembly and SCDA expansion elements, b stress boundary conditions, and c fluid boundary conditions with permeable boundaries at edges and flow boundary at the injection well

Fig. 11
Fig. 11 Hydraulic stimulation.a CT scanned Granite specimen after hydraulic stimulation under unconfined conditions, and b numerical simulation of the hydraulic fracture.Left, 3-dimensional reconstruc-

Fig. 13 Fig. 14
Fig. 13 Grain assemblies.The corresponding grain assemblies for minimum and maximum seed lengths for the four Voronoi distributions a 3-6 mm, b 3-7 mm, c 3-8 mm, d 3-9 mm, and e the equivalent particle size distributions

Fig. 15
Fig.15HS fracture propagation in a fissured granular assembly.Spatial distribution, pole positions in a stereonet, distribution of the defect contact areas (in m 2 ) for models with (a 5% and b 15% of contacts defined as defects.HS fracture pattern, block displacement field and b illustrate the fracture patterns observed from the first (SCDA) and second (HS) stages of pre-conditioning of a numerical specimen with 100 × 100 × 100 mm 3 dimensions fractured under stress Case 3 (σ x = σ y = 5.0 MPa).Cross sections of the observed fracture profile at planes (1)-(4) in Fig. 19b are shown in Fig. 19c.

Fig. 16
Fig. 16 The influence of energy transmission rate (flow rate) on the HS process.The HS fracture patterns observed for different injection flow rates a 18 ml/min, b 42 ml/min, c 66 ml/min, and d the corresponding injection pressure plots Fig. 18 SCDA-HS hybrid fracture in the presence of pre-existing fissures and anisotropic stress conditions (Case 2). a Fracture pattern produced by SCDA charging, b fracture extension of SCDA-charged specimen by hydraulic loading, (the red points indicate the position of fissures, and the background color field indicates block displacement)

Fig. 19
Fig. 19 Effect of grain heterogeneity on the 3-dimensional fracture propagation during hybrid SCDA-HS fracturing under stress Case 3, a SCDA-charged fracture, b HS fracture extension, c cross sections of

Table 1
Petrographic and mechanical properties of granite a Elastic properties were calculated using cylindrical specimens with 22.5 mm diameter and 45 mm height

Table 2
Calibrated micro-mechanical properties of the numerical assembly

Table 4
Model parameters used for the calibrated intact rock matrix flow

Table 6
In-situ stress states