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.

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.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Stress tensors and quantification of fracture patterns to analyze connectivity and potential fluid flow in a mesa landform of the Northern Andes

This work applies stress tensors inversions and quantification of fracture patterns along the Mesa de Los Santos, in the Eastern Cordillera of Colombia, to better understand the potential fluid flow. It thus contributes to the conceptual hydrogeological model. The area was subdivided into three blocks, separated by the NW Potreros and the Los Santos faults, having minor inner faults of different orientations. This separation facilitates the analysis of the fractures measured in the field, which in general show high dip angles and a conjugate geometry in the northern block, tension fractures(Mode I) in the central block, and a random distribution in the southwestern block. WinTensor treatment of slickensides yielded a maximum horizontal stress(SHmax) of 111o, which coincides with the WNW-ESE tensor observed from the conjugate and tension joints. We then used Frac Pa Q to generate interpolation maps of fracture intensity and density. The maps show the highest values in the central block and the lowest in the northern block,where the precipitation is higher, causing intensive rock weathering and homogenization of the fracture planes. Although the highest values of connectivity by line are found to the south of the mesa, we suggest the possibility of greater flow from the recharge zone(NE) along bedding planes and open NW-SE fractures.

Damage Characteristics Analysis and Fractal Study of Shale With Prefabricated Fractures under Thermal-mechanical Coupling

To study the damage and failure of shale with different fracture inclination angles under uniaxial compression loading,in this work,RFPA2D-Thermal,a two-dimensional real failure process analysis software,was used for numerical simulation.Numerical simulation results show that quartz in shale mainly affects the tensile and compressive strength of shale by increasing rock brittleness.The coupling of temperature and pressure will cause lateral and volume destruction of shale,which enables the shale body to be more easily broken.Fracture inclination is the key factor affecting shale damage patterns.The failure mode of shale with low-and high-angle fractures is mainly shear failure,and the compressive strength does not vary with crack inclination.The damage mode of obliquely intersecting fractured shale is slip damage along the fracture face,the compressive strength decreases and then increases with the fracture inclination,and a minimum value exists.The acoustic emission simulation results of the damage process effectively reflect the accumulated internal damage and macroscopic crack appearance until fracture instability when the prefabricated fractured shale is subjected to uniaxial compressive loading.The crack inclinations of 0°and 120℃ corresponds to the most complex"N"shape damage mode.The crack inclinations of 30°and 60°,and the damage mode is an inverted"λ"shape.

Mechanical properties and fracture surface roughness of thermally damaged granite under dynamic splitting

In order to understand the mechanical properties and the fracture surface roughness characteristics of thermally damaged granite under dynamic splitting,dynamic Brazilian splitting tests were conducted on granite samples after thermal treatment at 25,200,400,and 600℃.Results show that the dynamic peak splitting strength of thermally damaged granite samples increases with increasing strain rate,showing obvious strain‐rate sensitivity.With increasing temperature,thermally induced cracks in granite transform from intergranular cracks to intragranular cracks,and to a transgranular crack network.Thermally induced damages reduce the dynamic peak splitting strength and the maximum absorbed energy while increasing the peak radial strain.The fracture mode of the thermally damaged granite under dynamic loads is mode Ⅱ splitting failure.By using the axial roughness index Z2 a,the distribution ranges of the wedge‐shaped failure zones and the tensile failure zones in the fracture surfaces under dynamic Brazilian splitting can be effectively identified.The radial roughness index Z_(2)^(r)is sensitive to the strain rate and temperature.It shows a linear correlation with the peak splitting strength and the maximum absorbed energy and a linear negative correlation with the peak radial strain.Z_(2)^(r)can be used to quantitatively estimate the dynamic parameters based on the models proposed.

Factors influencing further vertebral height loss following percutaneous vertebroplasty in osteoporotic vertebral compression fractures:A 1-year follow-up study

BACKGROUND Osteoporotic vertebral compression fractures(OVCFs)contribute to back pain and functional limitations in older individuals,with percutaneous vertebroplasty(PVP)emerging as a minimally invasive treatment.However,further height loss post-PVP prompts investigation into contributing factors.AIM To investigate the factors associated with further height loss following PVP with cement augmentation in OVCF patients.METHODS A total of 200 OVCF patients who underwent successful PVP between January 2021 and December 2022 were included in this study.“Further height loss”during 1 year of follow-up in OVCF patients with bone edema was defined as a vertical height loss of≥4 mm.The study population was divided into two groups for analysis:The“No Further Height Loss group(n=179)”and the“Further Height Loss group(n=21).”RESULTS In comparing two distinct groups of patients,significant differences existed in bone mineral density(BMD),vertebral compression degree,prevalence of intravertebral cleft(IVF),type of bone cement used,and cement distribution patterns.Results from binary univariate regression analysis revealed that lower BMD,the presence of IVF,cleft distribution of bone cement,and higher vertebral compression degree were all significantly associated with further height loss.Notably,the use of mineralized collagen modified-poly(methyl methacrylate)bone cement was associated with a significant reduction in the risk of further height loss.In multivariate regression analysis,lower BMD and the presence of IVF remained significantly associated with further height loss.CONCLUSION Further height loss following PVP in OVCF patients is influenced by a complex interplay of factors,especially lower BMD and the presence of IVF.These findings underscore the importance of assessing and managing these factors when addressing height loss following PVP in OVCF patients.

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患者骶髂关节的结构性病变,并显示出良好的诊断性能。

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.

Shared genetic architecture highlights the bidirectional association between major depressive disorder and fracture risk

Background There is limited evidence suggesting that osteoporosis might exacerbate depressive symptoms,while more studies demonstrate that depression negatively affects bone density and increases fracture risk.Aims To explore the relationship between major depressive disorder(MDD)and fracture risk.Methods We conducted a nested case-control analysis(32670 patients with fracture and 397017 individuals without fracture)and a matched cohort analysis(16496 patients with MDD and 435492 individuals without MDD)in the same prospective UK Biobank data set.Further,we investigated the shared genetic architecture between MDD and fracture with linkage disequilibrium score regression and the MiXeR statistical tools.We used the conditional/conjunctional false discovery rate approach to identify the specific shared loci.We calculated the weighted genetic risk score for individuals in the UK Biobank and logistic regression was used to confirm the association observed in the prospective study.Results We found that MDD was associated with a 14%increase in fracture risk(hazard ratio(HR)1.14,95%CI 1.14 to 1.15,p<0.001)in the nested case-control analysis,while fracture was associated with a 72%increase in MDD risk(HR 1.72,95%CI 1.64 to 1.79,p<0.001)in the matched cohort analysis,suggesting a longitudinal and bidirectional relationship.Further,genetic summary data suggested a genetic overlap between MDD and fracture.Specifically,we identified four shared genomic loci,with the top signal(rs7554101)near SGIP1.The protein encoded by SGIP1 is involved in cannabinoid receptor type 1 signalling.We found that genetically predicted MDD was associated with a higher risk of fracture and vice versa.In addition,we found that the higher expression level of SGIP1 in the spinal cord and muscle was associated with an increased risk of fracture and MDD.Conclusions The genetic pleiotropy between MDD and fracture highlights the bidirectional association observed in the epidemiological analysis.The shared genetic components(such as SGIP1)between the diseases suggest that modulating the endocannabinoid system could be a potential therapeutic strategy for both MDD and bone loss.

Management and classification of the fracture of lateral process of talus:An overview and literature update

Fracture of the lateral process of the talus(FLPT)is uncommon in clinical practice and can be easily missed or misdiagnosed.In recent years,as researchers from all over the world have further deepened their research on FLPT,there has been a breakthrough in the classification,and the methods and principles of clinical management have changed accordingly;however,there is still no standardized guideline for the diagnosis and management of FLPT,and there have been few relevant literature review articles related to this kind of fracture in the past at least 5 years.In this article,we review the clinical classification,classification-based therapeutic recommendations,and prognosis of FLPT,with the aim of providing a reference for the clinical diagnosis and management of this infrequent fracture.

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.

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.

Creep constitutive model considering nonlinear creep degradation of fractured rock

Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.

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.

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.

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.

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.

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.

Effects of thawing-induced softening on fracture behaviors of frozen rock

Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.