Study of hydro-mechanical behaviours of rough rock fracture with shear dilatancy and asperities using shear-flow model

The geometric properties of fracture surfaces significantly influence shear-seepage in rock fractures,introducing complexities to fracture modelling.The present study focuses on the hydro-mechanical behaviours of rough rock fractures during shear-seepage processes to reveal how dilatancy and fracture asperities affect these phenomena.To achieve this,an improved shear-flow model(SFM)is proposed with the incorporation of dilatancy effect and asperities.In particular,shear dilatancy is accounted for in both the elastic and plastic stages,in contrast to some existing models that only consider it in the elastic stage.Depending on the computation approaches for the peak dilatancy angle,three different versions of the SFM are derived based on Mohr-Coulomb,joint roughness coefficient-joint compressive strength(JRC-JCS),and Grasselli’s theories.Notably,this is a new attempt that utilizes Grasselli’s model in shearseepage analysis.An advanced parameter optimization method is introduced to accurately determine model parameters,addressing the issue of local optima inherent in some conventional methods.Then,model performance is evaluated against existing experimental results.The findings demonstrate that the SFM effectively reproduces the shear-seepage characteristics of rock fracture across a wide range of stress levels.Further sensitivity analysis reveals how dilatancy and asperity affect hydraulic properties.The relation between hydro-mechanical properties(dilatancy displacement and hydraulic conductivity)and asperity parameters is analysed.Several profound understandings of the shear-seepage process are obtained by exploring the phenomenon under various conditions.

Discontinuity development patterns and the challenges for 3D discrete fracture network modeling on complicated exposed rock surfaces

Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.

Using fracture-based continuum modeling of coupled geomechanical-hydrological processes for numerical simulation of hydraulic fracturing

This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.

Algorithmic approach to discrete fracture network flow modeling in consideration of realistic connections in large-scale fracture networks

Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.

A new elastoplastic model for bolt-grouted fractured rock

Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.

Fourth-order phase-field modeling for brittle fracture in piezoelectric materials

Failure analyses of piezoelectric structures and devices are of engineering and scientific significance.In this paper,a fourth-order phase-field fracture model for piezoelectric solids is developed based on the Hamilton principle.Three typical electric boundary conditions are involved in the present model to characterize the fracture behaviors in various physical situations.A staggered algorithm is used to simulate the crack propagation.The polynomial splines over hierarchical T-meshes(PHT-splines)are adopted as the basis function,which owns the C1continuity.Systematic numerical simulations are performed to study the influence of the electric boundary conditions and the applied electric field on the fracture behaviors of piezoelectric materials.The electric boundary conditions may influence crack paths and fracture loads significantly.The present research may be helpful for the reliability evaluation of the piezoelectric structure in the future applications.

A multi-mechanism numerical simulation model for CO_(2)-EOR and storage in fractured shale oil reservoirs

Under the policy background and advocacy of carbon capture,utilization,and storage(CCUS),CO_(2)-EOR has become a promising direction in the shale oil reservoir industry.The multi-scale pore structure distribution and fracture structure lead to complex multiphase flow,comprehensively considering multiple mechanisms is crucial for development and CO_(2) storage in fractured shale reservoirs.In this paper,a multi-mechanism coupled model is developed by MATLAB.Compared to the traditional Eclipse300 and MATLAB Reservoir Simulation Toolbox(MRST),this model considers the impact of pore structure on fluid phase behavior by the modified Peng—Robinson equation of state(PR-EOS),and the effect simultaneously radiate to Maxwell—Stefan(M—S)diffusion,stress sensitivity,the nano-confinement(NC)effect.Moreover,a modified embedded discrete fracture model(EDFM)is used to model the complex fractures,which optimizes connection types and half-transmissibility calculation approaches between non-neighboring connections(NNCs).The full implicit equation adopts the finite volume method(FVM)and Newton—Raphson iteration for discretization and solution.The model verification with the Eclipse300 and MRST is satisfactory.The results show that the interaction between the mechanisms significantly affects the production performance and storage characteristics.The effect of molecular diffusion may be overestimated in oil-dominated(liquid-dominated)shale reservoirs.The well spacing and injection gas rate are the most crucial factors affecting the production by sensitivity analysis.Moreover,the potential gas invasion risk is mentioned.This model provides a reliable theoretical basis for CO_(2)-EOR and sequestration in shale oil reservoirs.

Modified MK model combined with ductile fracture criterion and its application in warm hydroforming

A modified MK model combined with ductile fracture criterion（DFC-MK model） is proposed to compute the forming limit diagrams（FLDs） of 5A06-O aluminum alloy sheet at different temperatures.The material constant（C） of ductile fracture criterion and initial thickness imperfection parameter（f0） at various temperatures are determined by using a new computing method based on wide sheet bending test.The FLDs at 20 and 200 °C are calculated through the DFC-MK model.The DFC-MK model,which includes the influence of through-thickness normal stress,is written into the subroutine VUMAT embedded in Abaqus/ Explicit.The cylindrical cup hydroforming tests are carried out to verify the model.The results show that compared with experimental observations,the predicted FLDs based on DFC-MK model are more accurate than the conventional MK model;the errors between the simulations and experiments in warm hydroforming are 8.23% at 20 °C and 9.24% at 200 °C,which verify the effectiveness of the proposed model.

Seismic wavefield modeling in media with fluid-filled fractures and surface topography

We present a finite difference （FD） method for the simulation of seismic wave fields in fractured medium with an irregular （non-fiat） free surface which is beneficial for interpreting exploration data acquired in mountainous regions. Fractures are introduced through the Coates-Schoenberg approach into the FD scheme which leads to local anisotropic properties of the media where fractures are embedded. To implement surface topography, we take advantage of the boundary-conforming grid and map a rectangular grid onto a curved one. We use a stable and explicit second-order accurate finite difference scheme to discretize the elastic wave equations （in a curvilinear coordinate system） in a 2D heterogeneous transversely isotropic medium with a horizontal axis of symmetry （HTI）. Efficiency tests performed by different numerical experiments clearly illustrate the influence of an irregular free surface on seismic wave propagation in fractured media which may be significant to mountain seismic exploration. The tests also illustrate that the scattered waves induced by the tips of the fracture are re-scattered by the features of the free surface topography. The scattered waves provoked by the topography are re-scattered by the fractures, especially Rayleigh wave scattering whose amplitudes are much larger than others and making it very difficult to identify effective information from the fractures.

3D DEM simulation of hard rock fracture in deep tunnel excavation induced by changes in principal stress magnitude and orientation

To achieve the loading of the stress path of hard rock,the spherical discrete element model(DEM)and the new flexible membrane technology were utilized to realize the transient loading of three principal stresses with arbitrary magnitudes and orientations.Furthermore,based on the deep tunnel of China Jinping Underground Laboratory II(CJPL-II),the deformation and fracture evolution characteristics of deep hard rock induced by excavation stress path were analyzed,and the mechanisms of transient loading-unloading and stress rotation-induced fractures were revealed from a mesoscopic perspective.The results indicated that the stressestrain curve exhibits different trends and degrees of sudden changes when subjected to transient changes in principal stress,accompanied by sudden changes in strain rate.Stress rotation induces spatially directional deformation,resulting in fractures of different degrees and orientations,and increasing the degree of deformation anisotropy.The correlation between the degree of induced fracture and the unloading magnitude of minimum principal stress,as well as its initial level is significant and positive.The process of mechanical response during transient unloading exhibits clear nonlinearity and directivity.After transient unloading,both the minimum principal stress and minimum principal strain rate decrease sharply and then tend to stabilize.This occurs from the edge to the interior and from the direction of the minimum principal stress to the direction of the maximum principal stress on theε1-ε3 plane.Transient unloading will induce a tensile stress wave.The ability to induce fractures due to changes in principal stress magnitude,orientation and rotation paths gradually increases.The analysis indicates a positive correlation between the abrupt change amplitude of strain rate and the maximum unloading magnitude,which is determined by the magnitude and rotation of principal stress.A high tensile strain rate is more likely to induce fractures under low minimum principal stress.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

The tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy: A case study in the Wangyao Oilfield of Ordos Basin, China

Due to inherent limits of data acquisition and geophysical data resolution, there are large uncertainties in the characterization of subsurface fractures. However, outcrop analogies can provide qualitative and quantitative information on a large number of fractures, based on which the accuracy of subsurface fracture characterization can be improved. Here we take the tectonic fracture modeling of an ultra-low permeability sandstone reservoir based on an outcrop analogy, a case study of the Chang6t~ Formation of the Upper Triassic Yanchang Group of the Wangyao Oilfield in the Ordos Basin of China. An outcrop at the edge of the basin is a suitable analog for the reservoir, but the prerequisite is that they must have equivalent previous stress fields, similar final structural characteristics, relative timing and an identical depositional environment and diagenesis. The relationship among fracture density, rock type and bed thickness based on the outcrop is one of the most important fracture distribution models, and can be used to interpret fracture density in individual wells quantitatively. Fracture orientation, dip, geometry and scale, also should be described and measured in the outcrop, and can be used together with structure restoration and single well fracture density interpretation to guide fracture intensity prediction on bed surfaces and to constrain the construction of the 3D fracture geometry model of the subsurface reservoir. The application of the above principles shows the outcrop-based tectonic fracture models of the target ultra-low permeability sandstone reservoir are consistent with fractures inferred from microseismic interpretation and tracer tests. This illustrated that the fracture modeling based on the outcrop analogy is reliable and can reduce the uncertainty in stochastic fracture modeling.

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.

Establishment and validation of standardized animal models of spinal cord injury by normal external force-caused fracture dislocation

The duplication of animal models plays a key role in spinal cord injury research; however, there has been limited study into normal, external force-derived fracture dislocation. This study adopted experimental devices, designed in-house, to construct standardized ventral and dorsal spinal cord injury animal models of 6 g and 17 g falling from a height of 2, 4, and 10 cm, and 15, 30 or 50 g transversal compression on the spinal cord. The results showed that gradual increases in the degree of histopathological injury led to decreased Tarlov and Basso, Beattie and Bresnahan scores for the behavioral test, and increased Ashworth scores for the hind limb. Furthermore, there was a gradual decline in the slope test in the rats with dorsal spinal cord injury that correlated to increases in the falling substance weight or falling height. Similar alterations were observed in the ventral spinal cord injured rats, proportional to the increase in compression weight. Our experimental findings indicate that the standardized experimental rat models of dorsal and ventral spinal cord injury are stable, reliable and reproducible.

Experimental study and modeling of hydromechanical behavior of concrete fracture

In this study, the hydromechanical behavior of a concrete fracture under coupled compressive and shear stresses was investigated. A special experimental device was designed to create a planar fracture in a cylindrical sample and to carry out different kinds of hydromechanical tests on the fracture. Four series of laboratory tests were performed on an ordinary concrete sample. Hydrostatic compression tests were first conducted to characterize the normal compressibility of the fracture. In the second series, direct shear tests were conducted on the fracture under different normal stresses. The maximal shear stress of the fracture was determined as a function of the normal stress. In the third series, fluid flow tests were carried out in view of characterizing the overall hydraulic conductivity of the fracture as a function of its opening and closure. Shear tests with a constant fluid pressure were finally performed to investigate the influence of fluid pressure on the deformation behavior of concrete fractures. Based on the experimental investigation, an elastoplastic model is proposed. This model takes into account the nonlinear elastic behavior of a fracture under normal compression and the plastic deformation and failure due to shear stress. The model was coupled with the classical Darcy's law to describe the fluid flow along the fracture by considering the variation of permeability with fracture aperture. Numerical results agree with experimental data from various laboratory tests.

Fractured reservoir modeling by discrete fracture network and seismic modeling in the Tarim Basin,China

Fractured reservoirs are an important target for oil and gas exploration in the Tarim Basin and the prediction of this type of reservoir is challenging.Due to the complicated fracture system in the Tarim Basin,the conventional AVO inversion method based on HTI theory to predict fracture development will result in some errors.Thus,an integrated research concept for fractured reservoir prediction is put forward in this paper.Seismic modeling plays a bridging role in this concept,and the establishment of an anisotropic fracture model by Discrete Fracture Network （DFN） is the key part.Because the fracture system in the Tarim Basin shows complex anisotropic characteristics,it is vital to build an effective anisotropic model.Based on geological,well logging and seismic data,an effective anisotropic model of complex fracture systems can be set up with the DFN method.The effective elastic coefficients,and the input data for seismic modeling can be calculated.Then seismic modeling based on this model is performed,and the seismic response characteristics are analyzed.The modeling results can be used in the following AVO inversion for fracture detection.

Evaluating the mechanical properties of anisotropic shale containing bedding and natural fractures with discrete element modeling

Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.

Evaluation of mountain slope stability considering the impact of geological interfaces using discrete fractures model

The stability of rock slope is often controlled by the existing discontinuous surfaces, such as discrete fractures, which are ubiquitously distributing in a geological medium. In contrast with the traditional approaches used in soil slope with a continuous assumption, the simulation methods of jointed rock slope are different from that of in soil slope. This paper presents a study on jointed rock slope stability using the proposed discontinuous approach, which considers the effects of discrete fractures. Comparing with traditional methods to model fractures in an implicit way, the presented approach provides a method to simulate fractures in an explicit way, where grids between rock matrix and fractures are independent. To complete geometric components generation and mesh partition for the model, the corresponding algorithms were devised. To evaluate the stability state of rock slope quantitatively, the strength reduction method was integrated into our analysis framework. A benchmark example was used to verify the validation of the approach. A jointed rock slope, which contains natural fractures, was selected as a case study and was simulated regarding the workflow of our framework. It was set up in the light of the geological condition of the site. Slope stability was evaluated under different loading conditions with various fracture patterns. Numerical results show that fractures have significant contributions to slope stability, and different fracture patterns would lead to different shapes of the slip surface. The devised method has the ability to calculate a non-circular slip surface, which is different from a circular slip surface obtained by classical methods.

A discrete model for prediction of radon flux from fractured rocks

Prediction of radon flux from the fractured zone of a propagating cave mine is basically associated with uncertainty and complexity. For instance, there is restricted access to these zones for field measure- ments, and it is quite difficult to replicate the complex nature of both natural and induced fractures in these zones in laboratory studies. Hence, a technique for predicting radon flux from a fractured rock using a discrete fracture network （DFN） model is developed to address these difficulties. This model quantifies the contribution of fractures to the total radon flux, and estimates the fracture density from a measured radon flux considering the effects of advection, diffusion, as well as radon generation and decay. Radon generation and decay are classified as reaction processes. Therefore, the equation solved is termed as the advection-diffusion-reaction equation （ADRE）. Peclet number （Pe）, a conventional dimensionless parameter that indicates the ratio of mass transport by advection to diffusion, is used to classify the transport regimes. The results show that the proposed model effectively predicts radon flux from a fractured rock. An increase in fracture density for a rock sample with uniformly distributed radon generation rate can elevate radon flux significantly compared with another rock sample with an equivalent increase in radon generation rate. In addition to Pe, two other independent dimensionless parameters （derived for radon transport through fractures） significantly affect radon dimensionless flux. Findings provide insight into radon transport through fractured rocks and can be used to improve radon control measures for proactive mitigation.

Dolomite fracture modeling using the Johnson-Holmquist concrete material model:Parameter determination and validation

In this paper,the Johnson-Holmquist concrete(JHC)constitutive model is adopted for modeling and simulating the fracture of dolomite.A detailed step-by-step procedure for determining all required parameters,based on a series of experiments under quasi-static and dynamic regimes,is proposed.Strain rate coefficients,failure surfaces,equations of state and damage/failure constants are acquired based on the experimental data and finite element analyses.The JHC model with the obtained parameters for dolomite is subsequently validated using quasi-static uniaxial and triaxial compression tests as well as dynamic split Hopkinson pressure bar(SHPB)tests.The influence of mesh size is also analyzed.It shows that the simulated fracture behavior and waveform data are in good agreement with the experimental data for all tests under both quasi-static and dynamic loading conditions.Future studies will implement the validated JHC model in small-and large-scale blasting simulations.