Evaluating the stability and volumetric flowback rate of proppant packs in hydraulic fractures using the lattice Boltzmann-discrete element coupling method

The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.

Pulp health and calcific healing of a complicated crown–root fracture with additional root fracture in a maxillary incisor: A case report

BACKGROUND Complicated crown–root fracture (CRF) involves severe injury to the crown, root,and pulp, and may be accompanied by multiple root fractures. The loss of a toothhas lifelong consequences for children and teenagers, but the maintenance of pulphealth and the calcific healing of multiple root fractures are rarely reported in theliterature.CASE SUMMARY This case reports healing of a permanent tooth with complicated crown–root andadditional root fractures, in which pulp health was maintained. A 10-year-old girlfell and fractured the root of her maxillary left central incisor at the cervical level.After the coronal fragment was repositioned, the tooth was splinted until thetooth was no longer mobile, 2 years later. Eight years after treatment, the toothhas remained asymptomatic with vital pulp and localized gingival overgrowth.Cone-beam computed tomography revealed not only calcified healing of the CRFbut also spontaneous healing in an additional undiagnosed root fracture. Thefracture line on the enamel could not be healed by hard tissue and formed agroove in the cervical crown. It was speculated that the groove was related to thelocalized gingival overgrowth.CONCLUSION This case provides a clinical perspective of the treatment of a tooth with acomplicated CRF and an additional root fracture.

Techniques for improving fracture-controlled stimulated reservoir volume in ultra-deep fractured tight reservoirs: A case study of Kuqa piedmont clastic reservoirs, Tarim Basin, NW China

Based on analysis of the reasons for low efficiency and low production after fracturing of some wells in the ultra-deep fractured tight reservoirs of the Kuqa piedmont zone, Tarim Basin and the matching relationship between the in-situ stress field and natural fractures, technological methods for creating complex fracture networks are proposed. Through theoretical study and large-scale physical simulation experiments, the mechanical conditions for forming complex fracture network in the Kuqa piedmont ultra-deep reservoirs are determined. The effectiveness of temporary plugging and diversion, and multi-stage fracturing to activate natural fractures and consequently realize multi-stage diversion is verified. The coupling effect of hydraulic fractures and natural fractures activating each other and resulting in "fracture swarms" is observed. These insights provide theoretical support for improving fracture-controlled stimulated reservoir volume(FSRV) in ultra-deep tight reservoirs. In addition, following the concept of volume fracturing technology and based on the results of fracture conductivity experiments of different processes, fracturing technologies such as multi-stage fracture-network acid fracturing, "multi-stage temporary plugging + secondary fracturing", fracturing of multiple small layers by vertically softness-and-hardness-oriented subdivision, and weighted-fluid refracturing are proposed to increase the FSRV. New environment-friendly weighted-fluid with low cost and new fracturing fluid system with low viscosity and high proppant-carrying capacity are also developed. These techniques have achieved remarkable results in field application.

Visualization of drained rock volume(DRV) in hydraulically fractured reservoirs with and without natural fractures using complex analysis methods(CAMs)

The drainage areas(and volumes)near hydraulically fractured wells,computed and visualized in our study at high resolution,may be critically affected by the presence of natural fractures.Using a recently developed algorithm based on complex analysis methods(CAMs),the drained rock volume(DRV)is visualized for a range of synthetic constellations of natural fractures near hydraulic fractures.First,flow interference effects near a single hydraulic fracture are systematically investigated for a variety of natural fracture sets.The permeability contrast between the matrix and the natural fractures is increased stepwise in order to better understand the effect on the DRV.Next,a larger-scale model investigates flow interference for a full hydraulically fractured well with a variety of natural fracture sets.The time of flight contours(TOFCs)outlining the DRV are for all cases with natural fractures compared to a base case without any natural fractures.Discrete natural fractures,with different orientations,hydraulic conductivity,and fracture density,may shift the TOFC patterns in the reservoir region drained by the hydraulically fractured well,essentially shifting the location of the well’s drainage area.The CAM-based models provide a computationally efficient method to quantify and visualize the drainage in both naturally and hydraulically fractured reservoirs.

Fracture Network Volume Fracturing Technology in High-temperature Hard Formation of Hot Dry Rock

It is more difficult for a hot dry rock to form a fracture network system than shale due to its special lithology, physical and mechanical properties under high temperature. The essential characteristics, rock mechanics and in-situ stress characteristics of a hot rock mass have been systematically studied by means of laboratory tests and true tri-axial physical simulation. The fracture initiation and propagation characteristics under different geological and engineering conditions are physically simulated, and the main controlling factors for the formation of a complex fracture network are revealed. The technology of low displacement for enhancing thermal cracking, gel fluid for expanding fracture and variable displacement cyclic injection for increasing a fracture network has been applied in the field, and good results have been achieved. Microseismic monitoring results demonstrate that complex fractures were formed in the field test, and the stimulation volume for heat exchanging reaches more than 3 million cubic meters. The research results play an important role in the stimulation technology of an enhanced geothermal system(EGS) and realize a breakthrough for power generation.

Analysis of fracture propagation and shale gas production by intensive volume fracturing

This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.

Drained rock volume around hydraulic fractures in porous media:planar fractures versus fractal networks

This study applies the Lindenmayer system based on fractal theory to generate synthetic fracture networks in hydraulically fractured wells.The applied flow model is based on complex analysis methods,which can quantify the flow near the fractures,and being gridless,is computationally faster than traditional discrete volume simulations.The representation of hydraulic fractures as fractals is a more realistic representation than planar bi-wing fractures used in most reservoir models.Fluid withdrawal from the reservoir with evenly spaced hydraulic fractures may leave dead zones between planar fractures.Complex fractal networks will drain the reservoir matrix more effectively,due to the mitigation of stagnation flow zones.The flow velocities,pressure response,and drained rock volume(DRV)are visualized for a variety of fractal fracture networks in a single-fracture treatment stage.The major advancement of this study is the improved representation of hydraulic fractures as complex fractals rather than restricting to planar fracture geometries.Our models indicate that when the complexity of hydraulic fracture networks increases,this will suppress the occurrence of dead flow zones.In order to increase the DRV and improve ultimate recovery,our flow models suggest that fracture treatment programs must find ways to create more complex fracture networks.

Revolutionizing Tight Reservoir Production: A Novel Dual-Medium Unsteady Seepage Model for Optimizing Volumetrically Fractured Horizontal Wells

This study presents an avant-garde approach for predicting and optimizing production in tight reservoirs,employing a dual-medium unsteady seepage model specifically fashioned for volumetrically fractured horizontal wells.Traditional models often fail to fully capture the complex dynamics associated with these unconventional reservoirs.In a significant departure from these models,our approach incorporates an initiation pressure gradient and a discrete fracture seepage network,providing a more realistic representation of the seepage process.The model also integrates an enhanced fluid-solid interaction,which allows for a more comprehensive understanding of the fluid-structure interactions in the reservoir.This is achieved through the incorporation of improved permeability and stress coupling,leading to more precise predictions of reservoir behavior.The numerical solutions derived from the model are obtained through the sophisticated finite element method,ensuring high accuracy and computational efficiency.To ensure the model’s reliability and accuracy,the outcomes were tested against a real-world case,with results demonstrating strong alignment.A key revelation from the study is the significant difference between uncoupled and fully coupled volumetrically fractured horizontal wells,challenging conventional wisdom in the field.Additionally,the study delves into the effects of stress,fracture length,and fracture number on reservoir production,contributing valuable insights for the design and optimization of tight reservoirs.The findings from this study have the potential to revolutionize the field of tight reservoir prediction and management,offering significant advancements in petroleum engineering.The proposed approach brings forth a more nuanced understanding of tight reservoir systems and opens up new avenues for optimizing reservoir management and production.

FRACTURE序列在中轴型脊柱关节炎骶髂关节结构性病变中的诊断价值

目的:评估限制回波间隔的快速梯度回波类CT成像(FRACTURE)序列在检测中轴型脊柱关节炎(SpA)骶髂关节结构病变中的诊断价值。方法:回顾性分析83例确诊中轴型SpA患者在2021年12月-2022年8月进行骶髂关节MRI和CT检查的影像资料。评估图像包括FRACTURE、T_(1)WI序列和CT平扫图像。对骶髂关节的影像图像进行结构性病变评分:关节间隙变化(0~5分)、骨质侵蚀(0~3分)和骨质硬化(0~2分)。根据修改后的纽约标准评分系统对骶髂关节炎进行综合评分。结果:总共有166个骶髂关节图像(83名受试者)可供分析。以CT图像作为参考标准,FRACTURE序列在关节间隙改变、骨质侵蚀、骨质硬化和骶髂关节炎综合评分与CT结果表现出较高的一致性,AUC分别为0.908、0.943、0.918和0.944。与T_(1)WI相比,FRACTURE在关节间隙变化(91.8%vs.70.5%)、骨质侵蚀(93.3%vs.56.4%)、骨质硬化(94%vs.77.6%)和骶髂关节炎综合评分(98.9%vs.75.3%)方面具有更高的诊断准确性。结论:FRACTURE成像可以评估中轴型SpA患者骶髂关节的结构性病变,并显示出良好的诊断性能。

A method to interpret fracture aperture of rock slope using adaptive shape and unmanned aerial vehicle multi-angle nap-of-the-object photogrammetry

The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.

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.

Complication rates after direct anterior vs posterior approach for hip hemiarthroplasty in elderly individuals with femoral neck fractures

BACKGROUND Dislocation rates after hemiarthroplasty reportedly vary from 1%to 17%.This serious complication is associated with increased morbidity and mortality rates.Approaches to this surgery are still debated,with no consensus regarding the superiority of any single approach.AIM To compare early postoperative complications after implementing the direct anterior and posterior approaches(PL)for hip hemiarthroplasty after femoral neck fractures.METHODS This is a comparative,retrospective,single-center cohort study conducted at a university hospital.Between March 2008 and December 2018,273 patients(a total of 280 hips)underwent bipolar hemiarthroplasties(n=280)for displaced femoral neck fractures using either the PL(n=171)or the minimally invasive direct anterior approach(DAA)(n=109).The choice of approach was related to the surgeons’practices;the implant types were similar and unrelated to the approach.Dislocation rates and other complications were reviewed after a minimum followup of 6 mo.RESULTS Both treatment groups had similarly aged patients(mean age:82 years),sex ratios,patient body mass indexes,and patient comorbidities.Surgical data(surgery delay time,operative time,and blood loss volume)did not differ significantly between the groups.The 30 d mortality rate was higher in the PL group(9.9%)than in the DAA group(3.7%),but the difference was not statistically significant(P=0.052).Among the one-month survivors,a significantly higher rate of dislocation was observed in the PL group(14/154;9.1%)than in the DAA group(0/105;0%)(P=0.002).Of the 14 patients with dislocation,8 underwent revision surgery for recurrent instability(posterior group),and one of them had 2 additional procedures due to a deep infection.The rate of other complications(e.g.,perioperative and early postoperative periprosthetic fractures and infection-related complications)did not differ significantly between the groups.CONCLUSION These findings suggest that the DAA to bipolar hemiarthroplasty for patients with femoral neck fractures is associated with a lower dislocation rate(<1%)than the PL.

Production induced fracture closure of deep shale gas well under thermo-hydro-mechanical conditions

Deep shale gas reservoirs have geological characteristics of high temperature,high pressure,high stress,and inferior ability to pass through fluids.The multi-stage fractured horizontal well is the key to exploiting the deep shale gas reservoir.However,during the production process,the effectiveness of the hydraulic fracture network decreases with the closure of fractures,which accelerates the decline of shale gas production.In this paper,we addressed the problems of unclear fracture closure mechanisms and low accuracy of shale gas production prediction during deep shale gas production.Then we established the fluid—solid—heat coupled model coupling the deformation and fluid flow among the fracture surface,proppant and the shale matrix.When the fluid—solid—heat coupled model was applied to the fracture network,it was well solved by our numerical method named discontinuous discrete fracture method.Compared with the conventional discrete fracture method,the discontinuous discrete fracture method can describe the three-dimensional morphology of the fracture while considering the effect of the change of fracture surface permeation coefficient on the coupled fracture—matrix flow and describing the displacement discontinuity across the fracture.Numerical simulations revealed that the degree of fracture closure increases as the production time proceeds,and the degree of closure of the secondary fractures is higher than that of the primary fractures.Shale creep and proppant embedment both increase the degree of fracture closure.The reduction in fracture surface permeability due to proppant embedment reduces the rate of fluid transfer between matrix and fracture,which has often been overlooked in the past.However,it significantly impacts shale gas production,with calculations showing a 24.7%cumulative three-year yield reduction.This study is helpful to understand the mechanism of hydraulic fracture closure.Therefore,it provides the theoretical guidance for maintaining the long-term effectiveness of hydraulic fractures.

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.

Effect of fracture fluid flowback on shale microfractures using CT scanning

The field data of shale fracturing demonstrate that the flowback performance of fracturing fluid is different from that of conventional reservoirs,where the flowback rate of shale fracturing fluid is lower than that of conventional reservoirs.At the early stage of flowback,there is no single-phase flow of the liquid phase in shale,but rather a gas-water two-phase flow,such that the single-phase flow model for tight oil and gas reservoirs is not applicable.In this study,pores and microfractures are extracted based on the experimental results of computed tomography(CT)scanning,and a spatial model of microfractures is established.Then,the influence of rough microfracture surfaces on the flow is corrected using the modified cubic law,which was modified by introducing the average deviation of the microfracture height as a roughness factor to consider the influence of microfracture surface roughness.The flow in the fracture network is simulated using the modified cubic law and the lattice Boltzmann method(LBM).The results obtained demonstrate that most of the fracturing fluid is retained in the shale microfractures,which explains the low fracturing fluid flowback rate in shale hydraulic fracturing.

Systemic modulation of skeletal mineralization by magnesium implant promoting fracture healing: Radiological exploration enhanced with PCA-based machine learning in a rat femoral model

The clinical application of magnesium(Mg)and its alloys for bone fractures has been well supported by in vitro and in vivo trials.However,there were studies indicating negative effects of high dose Mg intake and sustained local release of Mg ions on bone metabolism or repair,which should not be ignored when developing Mg-based implants.Thus,it remains necessary to assess the biological effects of Mg implants in animal models relevant to clinical treatment modalities.The primary purpose of this study was to validate the beneficial effects of intramedullary Mg implants on the healing outcome of femoral fractures in a modified rat model.In addition,the mineralization parameters at multiple anatomical sites were evaluated,to investigate their association with healing outcome and potential clinical applications.Compared to the control group without Mg implantation,postoperative imaging at week 12 demonstrated better healing outcomes in the Mg group,with more stable unions in 3D analysis and high-mineralized bridging in 2D evaluation.The bone tissue mineral density(TMD)was higher in the Mg group at the non-operated femur and lumbar vertebra,while no differences between groups were identified regarding the bone tissue volume(TV),TMD and bone mineral content(BMC)in humerus.In the surgical femur,the Mg group presented higher TMD,but lower TV and BMC in the distal metaphyseal region,as well as reduced BMC at the osteotomy site.Principal component analysis(PCA)-based machine learning revealed that by selecting clinically relevant parameters,radiological markers could be constructed for differentiation of healing outcomes,with better performance than 2D scoring.The study provides insights and preclinical evidence for the rational investigation of bioactive materials,the identification of potential adverse effects,and the promotion of diagnostic capabilities for fracture healing.

Experimental investigation on coal pore-fracture variation and fractal characteristics synergistically affected by solvents for improving clean gas extraction

Chemical solvents instead of pure water being as hydraulic fracturing fluid could effectively increase permeability and improve clean methane extraction efficiency.However,pore-fracture variation features of lean coal synergistically affected by solvents have not been fully understood.Ultrasonic testing,nuclear magnetic resonance analysis,liquid phase mass spectrometry was adopted to comprehensively analyze pore-fracture change characteristics of lean coal treated by combined solvent(NMP and CS_(2)).Meanwhile,quantitative characterization of above changing properties was conducted using geometric fractal theory.Relationship model between permeability,fractal dimension and porosity were established.Results indicate that the end face fractures of coal are well developed after CS2and combined solvent treatments,of which,end face box-counting fractal dimensions range from 1.1227 to 1.4767.Maximum decreases in ultrasonic longitudinal wave velocity of coal affected by NMP,CS_(2)and combined solvent are 2.700%,20.521%,22.454%,respectively.Solvent treatments could lead to increasing amount of both mesopores and macropores.Decrease ratio of fractal dimension Dsis 0.259%–2.159%,while permeability increases ratio of NMR ranges from 0.1904 to 6.4486.Meanwhile,combined solvent could dissolve coal polar and non-polar small molecules and expand flow space.Results could provide reference for solvent selection and parameter optimization of permeability-enhancement technology.

Numerical modeling of fracture propagation of supercritical CO_(2)compound fracturing

The exploitation of shale gas is promising due to depletion of the conventional energy and intensification of the greenhouse effect.In this paper,we proposed a heat-fluid-solid coupling damage model of supercritical CO_(2)(SC-CO_(2))compound fracturing which is expected to be an efficient and environmentally friendly way to develop shale gas.The coupling model is solved by the finite element method,and the results are in good agreement with the analytical solutions and fracturing experiments.Based on this model,the fracture propagation characteristics at the two stages of compound fracturing are studied and the influence of pressurization rate,in situ stress,bedding angle,and other factors are considered.The results show that at the SC-CO_(2)fracturing stage,a lower pressurization rate is conducive to formation of the branches around main fractures,while a higher pressurization rate inhibits formation of the branches around main fractures and promotes formation of the main fractures.Both bedding and in situ stress play a dominant role in the fracture propagation.When the in situ stress ratio(δ_(x)/δ_(y))is 1,the presence of bedding can reduce the initiation pressure and failure pressure.Nevertheless,it will cause the fracture to propagate along the bedding direction,reducing the fracture complexity.In rocks without bedding,hydraulic fracturing has the lengthening and widening effects for SC-CO_(2)induced fracture.In shale,fractures induced at the hydraulic fracturing stage are more likely to be dominated by in situ stresses and have a shorter reorientation radius.Therefore,fracture branches propagating along the maximum principal stress direction may be generated around the main fractures induced by SC-CO_(2)at the hydraulic fracturing stage.When the branches converge with the main fractures,fracture zones are easily formed,and thus the fracture complexity and damage area can be significantly increased.The results are instructive for the design and application of SC-CO_(2)compound fracturing.

Estimation of the anisotropy of hydraulic conductivity through 3D fracture networks using the directional geological entropy

With an extension of the geological entropy concept in porous media,the approach called directional entrogram is applied to link hydraulic behavior to the anisotropy of the 3D fracture networks.A metric called directional entropic scale is used to measure the anisotropy of spatial order in different directions.Compared with the traditional connectivity indexes based on the statistics of fracture geometry,the directional entropic scale is capable to quantify the anisotropy of connectivity and hydraulic conductivity in heterogeneous 3D fracture networks.According to the numerical analysis of directional entrogram and fluid flow in a number of the 3D fracture networks,the hydraulic conductivities and entropic scales in different directions both increase with spatial order(i.e.,trace length decreasing and spacing increasing)and are independent of the dip angle.As a result,the nonlinear correlation between the hydraulic conductivities and entropic scales from different directions can be unified as quadratic polynomial function,which can shed light on the anisotropic effect of spatial order and global entropy on the heterogeneous hydraulic behaviors.

A coupled thermo-mechanical peridynamic model for fracture behavior of granite subjected to heating and water-cooling processes

Thermal damage and thermal fracture of rocks are two important indicators in geothermal mining projects.This paper investigates the effects of heating and water-cooling on granite specimens at various temperatures.The laboratory uniaxial compression experiments were also conducted.Then,a coupled thermo-mechanical ordinary state-based peridynamic(OSB-PD)model and corresponding numerical scheme were developed to simulate the damage of rocks after the heating and cooling processes,and the change of crack evolution process was predicted.The results demonstrate that elevated heating temperatures exacerbate the thermal damage to the specimens,resulting in a decrease in peak strength and an increase in ductility of granite.The escalating occurrence of thermal-induced cracks significantly affects the crack evolution process during the loading phase.The numerical results accurately reproduce the damage and fracture characteristics of the granite under different final heating temperatures(FHTs),which are consistent with the test results in terms of strength,crack evolution process,and failure mode.