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

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.

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.

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.

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.

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.

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.

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.

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.

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.

A modified 3D mean strain energy density criterion for predicting shale mixed-mode Ⅰ/Ⅲ fracture toughness

The fracture toughness of rocks is a critical fracturing parameter in geo-energy exploitation playing a significant role in fracture mechanics and hydraulic fracturing.The edge-notched disk bending(ENDB)specimens are employed to measure the entire range of mixed-modeⅠ/Ⅲfracture toughness of Longmaxi shale.To theoretically interpret the fracture mechanisms,this research first introduces the detailed derivations of three established fracture criteria.By distinguishing the volumetric and distortional strain energy densities,an improved three-dimensional mean strain energy density(MSED)criterion is proposed.As the critical volumetric to distortional MSED ratio decreases,the transition from tensiondominated fracture to shear-dominated fracture is observed.Our results indicate that both peak load and applied energy increase significantly with the transition from pure mode I(i.e.,tension)to pure modeⅢ(i.e.,torsion or tearing)since mode-Ⅲcracking happens in a twisted manner and mode-Ⅰcracking occurs in a coplanar manner.The macroscopic fracture signatures are consistent with those of triaxial hydraulic fracturing.The average ratio of pure mode-Ⅲfracture toughness to pure mode-Ⅰfracture toughness is 0.68,indicating that the obtained mode-Ⅲfracture resistance for a tensionbased loading system is apparent rather than true.Compared to the three mainstream fracture criteria,the present fracture criterion exhibits greater competitiveness and can successfully evaluate and predict mixed-modeⅠ/Ⅲfracture toughness of distinct materials and loading methods.

Influences of clean fracturing fluid viscosity and horizontal in-situ stress difference on hydraulic fracture propagation and morphology in coal seam

The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.

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.

Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

Borehole instability in naturally fractured rocks poses significant challenges to drilling.Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure(P w)and pore pressure(P p)during drilling,which may cause wellbore instability.However,the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses,which may yield significant errors and misleading predictions.In addition,only limited factors were analyzed,and the fracture distribution was oversimplified.In this paper,the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM(distinct element method)and DFN(discrete fracture network)method.It provides estimates of the effect of fracture strength weakening,wellbore pressure,in situ stresses,and sealing efficiency on borehole stability.The results show that mud intrusion and weakening of fracture strength can damage the borehole.This is demonstrated by the large displacement around the borehole,shear displacement on natural fractures,and the generation of fracture at shear limit.Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure,which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations.A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore.Moreover,the effect of sealing natural fractures on maintaining borehole stability is verified in this study,and the increase in sealing efficiency reduces the radial invasion distance of drilling mud.This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations,which can consider the discrete fracture network,mud intrusion,and associated weakening of fracture strength.The information provided by the numerical approach(e.g.displacement around the borehole,shear displacement on fracture,and fracture at shear limit)is helpful for managing wellbore stability and designing wellbore-strengthening operations.

Identifying the real fracture hidden in rock microcrack zone by acoustic emission energy

Identifying the real fracture of rock hidden in acoustic emission(AE)source clusters(AE-depicted microcrack zone)remains challenging and crucial.Here we revealed the AE energy(representing dissipated energy)distribution rule in the rock microcrack zone and proposed an AE-energy-based method for identifying the real fracture.(1)A set of fracture experiments were performed on granite using wedgeloading,and the fracture process was detected and recorded by AE.The microcrack zone associated with the energy dissipation was characterized by AE sources and energy distribution,utilizing our selfdeveloped AE analysis program(RockAE).(2)The accumulated AE energy,an index representing energy dissipation,across the AE-depicted microcrack zone followed the normal distribution model(the mean and variance relate to the real fracture path and the microcrack zone width).This result implies that the nucleation and coalescence of massive cracks(i.e.,real fracture generation process)are supposed to follow a normal distribution.(3)Then,we obtained the real fracture extension path by joining the peak positions of the AE energy normal distribution curve at different cross-sections of the microcrack zone.Consequently,we distinguished between the microcrack zone and the concealed real fracture within it.The deviation was validated as slight as 1–3 mm.

Visualization and characterization of experimental hydraulic fractures interacting with karst fracture-cavity distributions

Karst fracture-cavity carbonate reservoirs,in which natural cavities are connected by natural fractures to form cavity clusters in many circumstances,have become significant fields of oil and gas exploration and exploitation.Proppant fracturing is considered as the best method for exploiting carbonate reservoirs;however,previous studies primarily focused on the effects of individual types of geological formations,such as natural fractures or cavities,on fracture propagation.In this study,true-triaxial physical simulation experiments were systematically performed under four types of stress difference conditions after the accurate prefabrication of four types of different fracture-cavity distributions in artificial samples.Subsequently,the interaction mechanism between the hydraulic fractures and fracture-cavity structures was systematically analyzed in combination with the stress distribution,cross-sectional morphology of the main propagation path,and three-dimensional visualization of the overall fracture network.It was found that the propagation of hydraulic fractures near the cavity was inhibited by the stress concentration surrounding the cavity.In contrast,a natural fracture with a smaller approach angle(0°and 30°)around the cavity can alleviate the stress concentration and significantly facilitate the connection with the cavity.In addition,the hydraulic fracture crossed the natural fracture at the 45°approach angle and bypassed the cavity under higher stress difference conditions.A new stimulation effectiveness evaluation index was established based on the stimulated reservoir area(SRA),tortuosity of the hydraulic fractures(T),and connectivity index(CI)of the cavities.These findings provide new insights into the fracturing design of carbonate reservoirs.

Effect mechanism of seepage force on the hydraulic fracture propagation

The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.

Fracture geometry and breakdown pressure of radial borehole fracturing in multiple layers

Radial borehole fracturing that combines radial boreholes with hydraulic fracturing is anticipated to improve the output of tight oil and gas reservoirs.This paper aims to investigate fracture propagation and pressure characteristics of radial borehole fracturing in multiple layers.A series of laboratory experiments with artificial rock samples(395 mm×395 mm×395 mm)was conducted using a true triaxial fracturing device.Three crucial factors corresponding to the vertical distance of adjacent radial borehole layers(vertical distance),the azimuth and diameter of the radial borehole are examined.Experimental results show that radial borehole fracturing in multiple layers generates diverse fracture geometries.Four types of fractures are identified based on the connectivity between hydraulic fractures and radial boreholes.The vertical distance significantly influences fracture propagation perpendicular to the radial borehole axis.An increase in the vertical distance impedes fracture connection across multiple radial borehole layers and reduces the fracture propagation distance along the radial borehole axis.The azimuth also influences fracture propagation along the radial borehole axis.Increasing the azimuth reduces the guiding ability of radial boreholes,which makes the fracture quickly curve to the maximum horizontal stress direction.The breakdown pressure correlates with diverse fracture geometries observed.When the fractures connect multi-layer radial boreholes,increasing the vertical distance decreases the breakdown pressure.Decreasing the azimuth and increasing the diameter also decrease the breakdown pressure.The extrusion force exists between the adjacent fractures generated in radial boreholes in multiple rows,which plays a crucial role in enhancing the guiding ability of radial boreholes and results in higher breakdown pressure.The research provides valuable theoretical insights for the field application of radial borehole fracturing technology in tight oil and gas reservoirs.

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.