Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hyd...Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.展开更多
Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System(EGS)design.We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data,where the overa...Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System(EGS)design.We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data,where the overall direction of hydraulic fractures does not show a delineated shape parallel to the maximum principal stress orientation.A field-scale numerical model based on the distinct element method is set up to carry out a fully coupled hydromechanical simulation,with the explicit representation of natural fractures via the discrete fracture network(DFN)approach.The effects of injection parameters and in situ stress on hydraulic fracture patterns are then quantitatively assessed.The study reveals that shear-induced deformation primarily governs the fracturing morphology in the GR1 well,driven by smaller injection rates and viscosities that promote massive activation of natural fractures,ultimately dominating the direction of hydraulic fracturing.Furthermore,the increase of in situ differential stress may promote shear damage of natural fracture surfaces,with the exact influence pattern depending on the combination of specific discontinuity properties and in situ stress state.Finally,we provide recommendations for EGS fracturing based on the influence characteristics of multiple parameters.This study can serve as an effective basis and reference for the design and optimization of EGS in the Gonghe basin and other sites.展开更多
In the concurrent extraction of coal and gas,the quantitative assessment of evolving characteristics in mining-induced fracture networks and mining-enhanced permeability within coal seams serves as the cornerstone for...In the concurrent extraction of coal and gas,the quantitative assessment of evolving characteristics in mining-induced fracture networks and mining-enhanced permeability within coal seams serves as the cornerstone for effective gas extraction.However,representing mining-induced fracture networks from a three-dimensional(3D)sight and developing a comprehensive model to evaluate the anisotropic mining-enhanced permeability characteristics still pose significant challenges.In this investigation,a field experiment was undertaken to systematically monitor the evolution of borehole fractures in the coal mass ahead of the mining face at the Pingdingshan Coal Mining Group in China.Using the testing data of borehole fracture,the mining-induced fracture network at varying distances from the mining face was reconstructed through a statistical reconstruction method.Additionally,utilizing fractal theory,a model for the permeability enhancement rate(PER)induced by mining was established.This model was employed to quantitatively depict the anisotropic evolution patterns of PER as the mining face advanced.The research conclusions are as follows:(1)The progression of the mining-induced fracture network can be classified into the stage of rapid growth,the stage of stable growth,and the stage of weak impact;(2)The PER of mining-induced fracture network exhibited a typical progression that can be characterized with slow growth,rapid growth and significant decline;(3)The anisotropic mining-enhanced permeability of the reconstructed mining-induced fracture networks were significant.The peak PER in the vertical direction of the coal seam is 6.86 times and 4446.38 times greater than the direction perpendicular to the vertical thickness and the direction parallel to the advancement of the mining face,respectively.This investigatione provides a viable approach and methodology for quantitatively assessing the anisotropic PER of fracture networks induced during mining,in the concurrent exploitation of coal and gas.展开更多
Fractures are critical to subsurface activities such as oil and gas extraction,geothermal energy production,and carbon storage.Hydraulic fracturing,a technique that enhances fluid production,creates complex fracture n...Fractures are critical to subsurface activities such as oil and gas extraction,geothermal energy production,and carbon storage.Hydraulic fracturing,a technique that enhances fluid production,creates complex fracture networks within rock formations containing natural discontinuities.Accurately distinguishing between hydraulically induced fractures and pre-existing discontinuities is essential for understanding hydraulic fracture mechanisms.However,this remains challenging due to the interconnected nature of fractures in three-dimensional(3D)space.Manual segmentation,while adaptive,is both labor-intensive and subjective,making it impractical for large-scale 3D datasets.This study introduces a deep learning-based progressive cross-sectional segmentation method to automate the classification of 3D fracture volumes.The proposed method was applied to a 3D hydraulic fracture network in a Montney cube sample,successfully segmenting natural fractures,parted bedding planes,and hydraulic fractures with minimal user intervention.The automated approach achieves a 99.6%reduction in manual image processing workload while maintaining high segmentation accuracy,with test accuracy exceeding 98%and F1-score over 84%.This approach generalizes well to Brazilian disc samples with different fracture patterns,achieving consistently high accuracy in distinguishing between bedding and non-bedding fractures.This automated fracture segmentation method offers an effective tool for enhanced quantitative characterization of fracture networks,which would contribute to a deeper understanding of hydraulic fracturing processes.展开更多
In order to identify the development characteristics of fracture network in tight conglomerate reservoir of Mahu after hydraulic fracturing,a hydraulic fracturing test site was set up in the second and third members o...In order to identify the development characteristics of fracture network in tight conglomerate reservoir of Mahu after hydraulic fracturing,a hydraulic fracturing test site was set up in the second and third members of Triassic Baikouquan Formation(T1b2 and T1b3)in Ma-131 well area,which learned from the successful experience of hydraulic fracturing test sites in North America(HFTS-1).Twelve horizontal wells and a high-angle coring well MaJ02 were drilled.The orientation,connection,propagation law and major controlling factors of hydraulic fractures were analyzed by comparing results of CT scans,imaging logs,direct observation of cores from Well MaJ02,and combined with tracer monitoring data.Results indicate that:(1)Two types of fractures have developed by hydraulic fracturing,i.e.tensile fractures and shear fractures.Tensile fractures are approximately parallel to the direction of the maximum horizontal principal stress,and propagate less than 50 m from perforation clusters.Shear fractures are distributed among tensile fractures and mainly in the strike-slip mode due to the induced stress field among tensile fractures,and some of them are in conjugated pairs.Overall,tensile fractures alternate with shear fractures,with shear fractures dominated and activated after tensile ones.(2)Tracer monitoring results indicate that communication between wells was prevalent in the early stage of production,and the static pressure in the fracture gradually decreased and the connectivity between wells reduced as production progressed.(3)Density of hydraulic fractures is mainly affected by the lithology and fracturing parameters,which is smaller in the mudstone than the conglomerate.Larger fracturing scale and smaller cluster spacing lead to a higher fracture density,which are important directions to improve the well productivity.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
The goal of this research is to develop mine-scale discrete fracture network(DFN)models in which the influence of the spatial heterogeneity of fracture distributions may be investigated on the rock wedge stability of ...The goal of this research is to develop mine-scale discrete fracture network(DFN)models in which the influence of the spatial heterogeneity of fracture distributions may be investigated on the rock wedge stability of an open pit slope.For this purpose,spatially conditioned DFN models were developed for the pit walls at Tasiast mine using comprehensive structural data from the mine.Using Sequential Gaussian Simulation(SGS),volumetric fracture intensities(P32)were modeled across the entire mine site in the form of 3D block models.The simulated P32 block models were used as the input constraints for conditional DFN fracture generation,where the DFN grid dimension is the same as the SGS 3D blocks.The spatially constrained DFN models were further calibrated using aerial fracture intensities(P21)data from the pit walls,obtained by a survey of the pit walls using an unmanned aerial vehicle(UAV)and measured traces of joints from 3D point cloud data.The final DFN model is expected to honor the fracture intensities gathered through different means with optimal model accuracy.Finally,bench-scale and interramp scale rock wedge slope stability analyses were conducted using the calibrated conditional DFN models.This work proves the significance of conditioned DFN models in rock wedge stability analysis.Such models provide detailed information regarding rock wedge stability so that site monitoring and prevention plans can be conducted with higher efficiency.展开更多
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 directi...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.展开更多
The complex network of fractures formed by randomly distributed natural fractures in hot-dry rocks(HDRs)complicates the heat transfer regularity of injected fluid.On the basis of the fracture network,exploring the cha...The complex network of fractures formed by randomly distributed natural fractures in hot-dry rocks(HDRs)complicates the heat transfer regularity of injected fluid.On the basis of the fracture network,exploring the characteristics of the fluid flow and heat transfer as influenced by different parameters helps enable efficient resource extraction and effectively promotes the construction of diversified energy utilization structures.Accordingly,accounting for the effect of the thermal shock on the evolution of the permeability of the rock matrix,a thermo-hydromechanical(THM)coupling model is developed to analyze the influences of fracture network characteristics on the heat extraction performance of HDRs.In addition,a large-scale injection and production physical simulation experiment is performed using a newly developed,in-house,large-scale true triaxial experimental system.The corresponding numerical model is established and validated.The good agreement between the numerical and experimental results verifies the reliability and accuracy of the proposed THM model.Subsequently,a two-dimensional model is established under complex fracture network conditions,taking,as a research object,the natural fracture characteristics of HDR in the Qinghai Gonghe Basin in combination with the regional geological information.The effects of different parameters,including the production well location,rock matrix permeability,injection rate,initial fracture width,and number of fractures,on the production temperature and heat extraction performance are systematically analyzed.The results indicate that an increase in the number of fractures,the distance between the injection well and the production well,or the width of the initial fractures leads to an improved heat extraction performance.The number of fractures increased from 11 horizontal fractures and 22 high-angle fractures to 35 horizontal fractures and 70 high-angle fractures,with a 20%increase in heat extraction rate.While the influence of the rock matrix permeability is not highly significant,it cannot be ignored.It is crucial to select an injection rate that is neither too low nor too high,taking into consideration economic factors.展开更多
The topological attributes of fracture networks in limestone,subject to intense hydrodynamics and intricate geological discontinuities,substantially influence the mechanical and hydraulic characteristics of the rock m...The topological attributes of fracture networks in limestone,subject to intense hydrodynamics and intricate geological discontinuities,substantially influence the mechanical and hydraulic characteristics of the rock mass.The dynamical evolution of fracture networks under stress is crucial for unveiling the interaction patterns among fractures.However,existing models are undirected graphs focused on stationary topology,which need optimization to depict fractures'dynamic development and rupture process.To compensate for the time and destruction terms,we propose the damage network model,which defines the physical interpretation of fractures through the ternary motif.We focus primarily on the evolution of node types,topological attributes,and motifs of the fracture network in limestone under uniaxial stress.Observations expose the varying behavior of the nodes'self-dynamics and neighbors'adjacent dynamics in the fracture network.This approach elucidates the impact of micro-crack behaviors on large brittle shear fractures from a topological perspective and further subdivides the progressive failure stage into four distinct phases(isolated crack growth phase,crack splay phase,damage coalescence phase,and mechanical failure phase)based on the significance profile of the motif.Regression analysis reveals a positive linear and negative power correlation between fracture network density and branch number to the rock damage resistance,respectively.The damage network model introduces a novel methodology for depicting the interaction of two-dimensional(2D)projected fractures,considering the dynamic spatiotemporal development characteristics and fracture geometric variation.It helps dynamically characterize properties such as connectivity,permeability,and damage factors while comprehensively assessing damage in rock mass fracture networks.展开更多
Class III tight oil reservoirs have low porosity and permeability,which are often responsible for low production rates and limited recovery.Extensive repeated fracturing is a well-known technique to fix some of these ...Class III tight oil reservoirs have low porosity and permeability,which are often responsible for low production rates and limited recovery.Extensive repeated fracturing is a well-known technique to fix some of these issues.With such methods,existing fractures are refractured,and/or new fractures are created to facilitate communication with natural fractures.This study explored how different refracturing methods affect horizontal well fracture networks,with a special focus on morphology and related fluid flow changes.In particular,the study relied on the unconventional fracture model(UFM).The evolution of fracture morphology and flow field after the initial fracturing were analyzed accordingly.The simulation results indicated that increased formation energy and reduced reservoir stress differences can promote fracture expansion.It was shown that the length of the fracture network,the width of the fracture network,and the complexity of the fracture can be improved,the oil drainage area can be increased,the distance of oil and gas seepage can be reduced,and the production of a single well can be significantly increased.展开更多
Proppant transport within fractures is one of the most critical tasks in oil,gas and geothermal reservoir stimulation,as it largely determines the ultimate performance of the operating well.Proppant transport in rough...Proppant transport within fractures is one of the most critical tasks in oil,gas and geothermal reservoir stimulation,as it largely determines the ultimate performance of the operating well.Proppant transport in rough fracture networks is still a relatively new area of research and the associated transport mechanisms are still unclear.In this study,representative parameters of rough fracture surfaces formed by supercritical CO_(2) fracturing were used to generate a rough fracture network model based on a spectral synthesis method.Computational fluid dynamics(CFD)coupled with the discrete element method(DEM)was used to study proppant transport in this rough fracture network.To reveal the turning transport mechanism of proppants into branching fractures at the intersections of rough fracture networks,a comparison was made with the behavior within smooth fracture networks,and the effect of key pumping parameters on the proppant placement in a secondary fracture was analyzed.The results show that the transport behavior of proppant in rough fracture networks is very different from that of the one in the smooth fracture networks.The turning transport mechanisms of proppant into secondary fractures in rough fracture networks are gravity-driven sliding,high velocity fluid suspension,and fracture structure induction.Under the same injection conditions,supercritical CO_(2)with high flow Reynolds number still has a weaker ability to transport proppant into secondary fractures than water.Thickening of the supercritical CO_(2)needs to be increased beyond a certain value to have a significant effect on proppant carrying,and under the temperature and pressure conditions of this paper,it needs to be increased more than 20 times(about 0.94 m Pa s).Increasing the injection velocity and decreasing the proppant concentration facilitates the entry of proppant into the branching fractures,which in turn results in a larger stimulated reservoir volume.The results help to understand the proppant transport and placement process in rough fracture networks formed by reservoir stimulation,and provide a theoretical reference for the optimization of proppant pumping parameters in hydraulic fracturing.展开更多
The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective meth...The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective method to improve oil recovery factor from unconventional oil reservoirs. Hydrocarbon gas huff-n-puff becomes preferable when the CO_(2) source is limited. However, the impact of complex fracture networks and well interference on the EOR performance of multiple MFHWs is still unclear. The optimal gas huff-n-puff parameters are significant for enhancing oil recovery. This work aims to optimize the hydrocarbon gas injection and production parameters for multiple MFHWs with complex fracture networks in unconventional oil reservoirs. Firstly, the numerical model based on unstructured grids is developed to characterize the complex fracture networks and capture the dynamic fracture features.Secondly, the PVT phase behavior simulation was carried out to provide the fluid model for numerical simulation. Thirdly, the optimal parameters for hydrocarbon gas huff-n-puff were obtained. Finally, the dominant factors of hydrocarbon gas huff-n-puff under complex fracture networks are obtained by fuzzy mathematical method. Results reveal that the current pressure of hydrocarbon gas injection can achieve miscible displacement. The optimal injection and production parameters are obtained by single-factor analysis to analyze the effect of individual parameter. Gas injection time is the dominant factor of hydrocarbon gas huff-n-puff in unconventional oil reservoirs with complex fracture networks. This work can offer engineers guidance for hydrocarbon gas huff-n-puff of multiple MFHWs considering the complex fracture networks.展开更多
Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual conne...Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.展开更多
Pressure buildup testing can be used to analyze fracture network characteristics and conduct quantitative interpretation of relevant parameters for shale gas wells,thus providing bases for assessing the well productiv...Pressure buildup testing can be used to analyze fracture network characteristics and conduct quantitative interpretation of relevant parameters for shale gas wells,thus providing bases for assessing the well productivity and formulating proper development strategies.This study establishes a new well test interpretation model for fractured horizontal wells based on seepage mechanisms of shale reservoirs and proposes a method for identifying fracturing patterns based on the characteristic slopes of pressure buildup curves and curve combination patterns.The pressure buildup curve patterns are identified to represent three types of shale reservoirs in the Sichuan Basin,namely the moderately deep shale reservoirs with high pressure,deep shale reservoirs with ultra-high pressure,and moderately deep shale reservoirs with normal pressure.Based on this,the relationship between the typical pressure buildup curve patterns and the fracture network types are put forward.Fracturing effects of three types of shale gas reservoir are compared and analyzed.The results show that typical flow patterns of shale reservoirs include bilinear flow in primary and secondary fractures,linear flow in secondary fractures,bilinear flow in secondary fractures and matrix,and linear flow in matrix.The fracture network characteristics can be determined using the characteristic slopes of pressure buildup curves and curve combinations.The linear flow in early secondary fractures is increasingly distinct with an increase in primary fracture conductivity.Moreover,the bilinear flow in secondary fractures and matrix and the subsequent linear flow in the matrix occur as the propping and density of secondary fractures increase.The increase in the burial depth,in-situ stress,and stress difference corresponds to a decrease in the propping of primary fractures that expand along different directions in the shale gas wells in the Sichuan Basin.Four pressure buildup curve patterns exist in the Sichuan Basin and its periphery.The pattern of pressure buildup curves of shale reservoirs in the Yongchuan area can be described as 1/2/→1/4,indicating limited stimulated reservoir volume,poorly propped secondary fractures,and the forming of primary fractures that extend only to certain directions.The pressure buildup curves of shale reservoirs in the main block of the Fuling area show a pattern of 1/4/→1/2 or 1/2,indicating greater stimulated reservoir volume,well propped secondary fractures,and the forming of complex fracture networks.The pattern of pressure buildup curves of shale reservoirs in the Pingqiao area is 1/2/→1/4→/1/2,indicating a fracturing effect somewhere between that of the Fuling and Yongchuan areas.For reservoirs with normal pressure,it is difficult to determine fracture network characteristics from pressure buildup curves due to insufficient formation energy and limited liquid drainage.展开更多
Prepulse combined hydraulic fracturing facilitates the development of fracture networks by integrating prepulse hydraulic loading with conventional hydraulic fracturing.The formation mechanisms of fracture networks be...Prepulse combined hydraulic fracturing facilitates the development of fracture networks by integrating prepulse hydraulic loading with conventional hydraulic fracturing.The formation mechanisms of fracture networks between hydraulic and pre-existing fractures under different prepulse loading parameters remain unclear.This research investigates the impact of prepulse loading parameters,including the prepulse loading number ratio(C),prepulse loading stress ratio(S),and prepulse loading frequency(f),on the formation of fracture networks between hydraulic and pre-existing fractures,using both experimental and numerical methods.The results suggest that low prepulse loading stress ratios and high prepulse loading number ratios are advantageous loading modes.Multiple hydraulic fractures are generated in the specimen under the advantageous loading modes,facilitating the development of a complex fracture network.Fatigue damage occurs in the specimen at the prepulse loading stage.The high water pressure at the secondary conventional hydraulic fracturing promotes the growth of hydraulic fractures along the damage zones.This allows the hydraulic fractures to propagate deeply and interact with pre-existing fractures.Under advantageous loading conditions,multiple hydraulic fractures can extend to pre-existing fractures,and these hydraulic fractures penetrate or propagate along pre-existing fractures.Especially when the approach angle is large,the damage range in the specimen during the prepulse loading stage increases,resulting in the formation of more hydraulic fractures.展开更多
As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness dis...As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness distribution functions of the Bakhtiary dam site and Oskarshamn/Forsmark mountain were fitted using statistical methods.The COMSOL Multiphysics finite element software was utilized to analyze the effects of fracture roughness distribution types and empirical formulas for fracture hydraulic aperture on the seepage field and temperature field of rock masses.The results show that:(1)The fracture roughness at the Bakhtiary dam site and Oskarshamn/Forsmark mountain follows lognormal and normal distributions,respectively;(2)For rock masses with the same expected value and standard deviation of fracture roughness,the outflow from rock masses with lognormal distribution of fracture roughness is significantly larger than that of rock masses with normal distribution of fracture roughness;(3)The fracture hydraulic aperture,outflow,and cold front distance of the Li and Jiang model are significantly larger than those of the Barton model;(4)The outflow,hydraulic pressure distribution,and temperature distribution of the Barton model are more sensitive to the fracture roughness distribution type than those of the Li and Jiang model.展开更多
Curtain grouting projects are characterized by their large scale and complexity,presenting significant challenges for real-time prediction of grout penetration using traditional methods.This study introduces an intell...Curtain grouting projects are characterized by their large scale and complexity,presenting significant challenges for real-time prediction of grout penetration using traditional methods.This study introduces an intelligent prediction method for grouting in fractured rock masses based on three core principles:integration of multi-source input features,fracture voxel modeling,and shortest path in sequential grouting.Three categories of data(geological structure data,grouting environmental data,and grouting operation data in the concept of a grouting geological model)are integrated and served as multi-source structured data in the intelligent prediction of grouting.A voxelization model quantifies the spatial characteristics of fractures,with voxel size optimized for capturing grouting paths.A shortest path algorithm based on a hierarchical solution is then developed to calculate grout penetration distances in the process of sequential grouting.A complete analysis framework is established,from the voxelization of the fracture network model to precise voxel classification,ultimately achieving an accurate prediction of grout penetration.The method demonstrates excellent performance on the test set,with validation against numerical methods in single-fracture and sequential grouting scenarios confirming its accuracy and prediction efficiency as hundreds of times faster than numerical methods.Application to the Dongzhuang hydraulic project’s grouting test area further validates its effectiveness in multi-hole grouting scenarios.展开更多
This study investigated the micro-sliding frictional behavior of shale in fracturing fluids under varying operational conditions using Chang 7 shale oil reservoir core samples.Through systematic micro-sliding friction...This study investigated the micro-sliding frictional behavior of shale in fracturing fluids under varying operational conditions using Chang 7 shale oil reservoir core samples.Through systematic micro-sliding friction experiments,the characteristics and governing mechanisms of shale friction were elucidated.Complementary analyses were conducted to characterize the mineral composition,petrophysical properties,and micromorphology of the shale samples,providing insights into the relationship between microscopic structure and frictional response.In this paper,the characteristics and variation law of shale micro-sliding friction under different types of graphite materials as additives in LGF-80(Low-damage Guar Fluid)oil flooding recoverable fracturing fluid system were mainly studied.In addition,the finite element numerical simulation experiment of hydraulic fracturing was adopted to study the influence of the friction coefficient of natural fracture surfaces on fracture propagation and formation of the fracture network.The geometric complexity of fracture networks was systematically quantified under varying frictional coefficients of natural fracture surfaces through multi-parametric characterization and morphometric analysis.The research results show that graphite micro-particles reduce friction and drag.Based on this,this paper proposes a new idea of graphite micro-particles as an additive in the LGF-80 oil flooding recoverable fracturing fluid system to reduce friction on the fracture surface.展开更多
Based on the finite element-discrete element numerical method,a numerical model of fracture propagation in deflagration fracturing was established by considering the impact of stress wave,quasi-static pressure of expl...Based on the finite element-discrete element numerical method,a numerical model of fracture propagation in deflagration fracturing was established by considering the impact of stress wave,quasi-static pressure of explosive gas,and reflection of stress wave.The model was validated against the results of physical experiments.Taking the shale reservoirs of Silurian Longmaxi Formation in Luzhou area of the Sichuan Basin as an example,the effects of in-situ stress difference,natural fracture parameters,branch wellbore spacing,delay detonation time,and angle between branch wellbore and main wellbore on fracture propagation were identified.The results show that the fracture propagation morphology in deflagration fracturing is less affected by the in-situ stress difference when it is 5-15 MPa,and the tendency of fracture intersection between branch wellbores is significantly weakened when the in-situ stress difference reaches 20 MPa.The increase of natural fracture length promotes the fracture propagation along the natural fracture direction,while the increase of volumetric natural fracture density and angle limits the fracture propagation area and reduces the probability of fracture intersection between branch wells.The larger the branch wellbore spacing,the less probability of the fracture intersection between branch wells,allowing for the fracture propagation in multiple directions.Increasing the delay detonation time decreases the fracture spacing between branch wellbores.When the angle between the branch wellbore and the main wellbore is 45°and 90°,there is a tendency of fracture intersection between branch wellbores.展开更多
基金Australian Research Council Linkage Program(LP200301404)for sponsoring this researchthe financial support provided by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology,SKLGP2021K002)National Natural Science Foundation of China(52374101,32111530138).
文摘Discrete fracture network(DFN)commonly existing in natural rock masses plays an important role in geological complexity which can influence rock fracturing behaviour during fluid injection.This paper simulated the hydraulic fracturing process in lab-scale coal samples with DFNs and the induced seismic activities by the discrete element method(DEM).The effects of DFNs on hydraulic fracturing,induced seismicity and elastic property changes have been concluded.Denser DFNs can comprehensively decrease the peak injection pressure and injection duration.The proportion of strong seismic events increases first and then decreases with increasing DFN density.In addition,the relative modulus of the rock mass is derived innovatively from breakdown pressure,breakdown fracture length and the related initiation time.Increasing DFN densities among large(35–60 degrees)and small(0–30 degrees)fracture dip angles show opposite evolution trends in relative modulus.The transitional point(dip angle)for the opposite trends is also proportionally affected by the friction angle of the rock mass.The modelling results have much practical meaning to infer the density and geometry of pre-existing fractures and the elastic property of rock mass in the field,simply based on the hydraulic fracturing and induced seismicity monitoring data.
基金support from the National Natural Science Foundation of China(Grant Nos.42320104003,42177175,and 42077247)the Fundamental Research Funds for the Central Universities.
文摘Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System(EGS)design.We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data,where the overall direction of hydraulic fractures does not show a delineated shape parallel to the maximum principal stress orientation.A field-scale numerical model based on the distinct element method is set up to carry out a fully coupled hydromechanical simulation,with the explicit representation of natural fractures via the discrete fracture network(DFN)approach.The effects of injection parameters and in situ stress on hydraulic fracture patterns are then quantitatively assessed.The study reveals that shear-induced deformation primarily governs the fracturing morphology in the GR1 well,driven by smaller injection rates and viscosities that promote massive activation of natural fractures,ultimately dominating the direction of hydraulic fracturing.Furthermore,the increase of in situ differential stress may promote shear damage of natural fracture surfaces,with the exact influence pattern depending on the combination of specific discontinuity properties and in situ stress state.Finally,we provide recommendations for EGS fracturing based on the influence characteristics of multiple parameters.This study can serve as an effective basis and reference for the design and optimization of EGS in the Gonghe basin and other sites.
基金supported by the National Natural Science Foundation of China (Grant No.42377143)Sichuan Natural Science Foundation (Grant No.2024NSFSC0097)the Open Fund of State Key Laboratory of Coal Mining and Clean Utilization,China (Grant No.2021-CMCU-KFZD001).
文摘In the concurrent extraction of coal and gas,the quantitative assessment of evolving characteristics in mining-induced fracture networks and mining-enhanced permeability within coal seams serves as the cornerstone for effective gas extraction.However,representing mining-induced fracture networks from a three-dimensional(3D)sight and developing a comprehensive model to evaluate the anisotropic mining-enhanced permeability characteristics still pose significant challenges.In this investigation,a field experiment was undertaken to systematically monitor the evolution of borehole fractures in the coal mass ahead of the mining face at the Pingdingshan Coal Mining Group in China.Using the testing data of borehole fracture,the mining-induced fracture network at varying distances from the mining face was reconstructed through a statistical reconstruction method.Additionally,utilizing fractal theory,a model for the permeability enhancement rate(PER)induced by mining was established.This model was employed to quantitatively depict the anisotropic evolution patterns of PER as the mining face advanced.The research conclusions are as follows:(1)The progression of the mining-induced fracture network can be classified into the stage of rapid growth,the stage of stable growth,and the stage of weak impact;(2)The PER of mining-induced fracture network exhibited a typical progression that can be characterized with slow growth,rapid growth and significant decline;(3)The anisotropic mining-enhanced permeability of the reconstructed mining-induced fracture networks were significant.The peak PER in the vertical direction of the coal seam is 6.86 times and 4446.38 times greater than the direction perpendicular to the vertical thickness and the direction parallel to the advancement of the mining face,respectively.This investigatione provides a viable approach and methodology for quantitatively assessing the anisotropic PER of fracture networks induced during mining,in the concurrent exploitation of coal and gas.
基金supported through the Natural Sciences and Engineering Research Council of Canada(NSERC)Discovery Grants 341275,CRDPJ 543894-19NSERC/Energi Simulation Industrial Research Chair program.
文摘Fractures are critical to subsurface activities such as oil and gas extraction,geothermal energy production,and carbon storage.Hydraulic fracturing,a technique that enhances fluid production,creates complex fracture networks within rock formations containing natural discontinuities.Accurately distinguishing between hydraulically induced fractures and pre-existing discontinuities is essential for understanding hydraulic fracture mechanisms.However,this remains challenging due to the interconnected nature of fractures in three-dimensional(3D)space.Manual segmentation,while adaptive,is both labor-intensive and subjective,making it impractical for large-scale 3D datasets.This study introduces a deep learning-based progressive cross-sectional segmentation method to automate the classification of 3D fracture volumes.The proposed method was applied to a 3D hydraulic fracture network in a Montney cube sample,successfully segmenting natural fractures,parted bedding planes,and hydraulic fractures with minimal user intervention.The automated approach achieves a 99.6%reduction in manual image processing workload while maintaining high segmentation accuracy,with test accuracy exceeding 98%and F1-score over 84%.This approach generalizes well to Brazilian disc samples with different fracture patterns,achieving consistently high accuracy in distinguishing between bedding and non-bedding fractures.This automated fracture segmentation method offers an effective tool for enhanced quantitative characterization of fracture networks,which would contribute to a deeper understanding of hydraulic fracturing processes.
基金Supported by the National Natural Science Foundation of China(52274051)CNPC-China University of Petroleum(Beijing)Strategic Cooperative Project(ZLZX2020-01).
文摘In order to identify the development characteristics of fracture network in tight conglomerate reservoir of Mahu after hydraulic fracturing,a hydraulic fracturing test site was set up in the second and third members of Triassic Baikouquan Formation(T1b2 and T1b3)in Ma-131 well area,which learned from the successful experience of hydraulic fracturing test sites in North America(HFTS-1).Twelve horizontal wells and a high-angle coring well MaJ02 were drilled.The orientation,connection,propagation law and major controlling factors of hydraulic fractures were analyzed by comparing results of CT scans,imaging logs,direct observation of cores from Well MaJ02,and combined with tracer monitoring data.Results indicate that:(1)Two types of fractures have developed by hydraulic fracturing,i.e.tensile fractures and shear fractures.Tensile fractures are approximately parallel to the direction of the maximum horizontal principal stress,and propagate less than 50 m from perforation clusters.Shear fractures are distributed among tensile fractures and mainly in the strike-slip mode due to the induced stress field among tensile fractures,and some of them are in conjugated pairs.Overall,tensile fractures alternate with shear fractures,with shear fractures dominated and activated after tensile ones.(2)Tracer monitoring results indicate that communication between wells was prevalent in the early stage of production,and the static pressure in the fracture gradually decreased and the connectivity between wells reduced as production progressed.(3)Density of hydraulic fractures is mainly affected by the lithology and fracturing parameters,which is smaller in the mudstone than the conglomerate.Larger fracturing scale and smaller cluster spacing lead to a higher fracture density,which are important directions to improve the well productivity.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金Kinross Gold and MITACS for their financial support(Grant No.FR42880).
文摘The goal of this research is to develop mine-scale discrete fracture network(DFN)models in which the influence of the spatial heterogeneity of fracture distributions may be investigated on the rock wedge stability of an open pit slope.For this purpose,spatially conditioned DFN models were developed for the pit walls at Tasiast mine using comprehensive structural data from the mine.Using Sequential Gaussian Simulation(SGS),volumetric fracture intensities(P32)were modeled across the entire mine site in the form of 3D block models.The simulated P32 block models were used as the input constraints for conditional DFN fracture generation,where the DFN grid dimension is the same as the SGS 3D blocks.The spatially constrained DFN models were further calibrated using aerial fracture intensities(P21)data from the pit walls,obtained by a survey of the pit walls using an unmanned aerial vehicle(UAV)and measured traces of joints from 3D point cloud data.The final DFN model is expected to honor the fracture intensities gathered through different means with optimal model accuracy.Finally,bench-scale and interramp scale rock wedge slope stability analyses were conducted using the calibrated conditional DFN models.This work proves the significance of conditioned DFN models in rock wedge stability analysis.Such models provide detailed information regarding rock wedge stability so that site monitoring and prevention plans can be conducted with higher efficiency.
基金supported by the National Natural Science Foundation of China(Nos.42077243,52209148,and 52079062).
文摘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.
基金supported by the Major Program of National Natural Science Foundation of China(No.52192622)the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project(No.SKLGP2022Z018).
文摘The complex network of fractures formed by randomly distributed natural fractures in hot-dry rocks(HDRs)complicates the heat transfer regularity of injected fluid.On the basis of the fracture network,exploring the characteristics of the fluid flow and heat transfer as influenced by different parameters helps enable efficient resource extraction and effectively promotes the construction of diversified energy utilization structures.Accordingly,accounting for the effect of the thermal shock on the evolution of the permeability of the rock matrix,a thermo-hydromechanical(THM)coupling model is developed to analyze the influences of fracture network characteristics on the heat extraction performance of HDRs.In addition,a large-scale injection and production physical simulation experiment is performed using a newly developed,in-house,large-scale true triaxial experimental system.The corresponding numerical model is established and validated.The good agreement between the numerical and experimental results verifies the reliability and accuracy of the proposed THM model.Subsequently,a two-dimensional model is established under complex fracture network conditions,taking,as a research object,the natural fracture characteristics of HDR in the Qinghai Gonghe Basin in combination with the regional geological information.The effects of different parameters,including the production well location,rock matrix permeability,injection rate,initial fracture width,and number of fractures,on the production temperature and heat extraction performance are systematically analyzed.The results indicate that an increase in the number of fractures,the distance between the injection well and the production well,or the width of the initial fractures leads to an improved heat extraction performance.The number of fractures increased from 11 horizontal fractures and 22 high-angle fractures to 35 horizontal fractures and 70 high-angle fractures,with a 20%increase in heat extraction rate.While the influence of the rock matrix permeability is not highly significant,it cannot be ignored.It is crucial to select an injection rate that is neither too low nor too high,taking into consideration economic factors.
基金supported by the National Natural Science Foundation of China(Grant No.52090081)the State Key Laboratory of Hydroscience and Engineering(Grant No.2022-KY-02).
文摘The topological attributes of fracture networks in limestone,subject to intense hydrodynamics and intricate geological discontinuities,substantially influence the mechanical and hydraulic characteristics of the rock mass.The dynamical evolution of fracture networks under stress is crucial for unveiling the interaction patterns among fractures.However,existing models are undirected graphs focused on stationary topology,which need optimization to depict fractures'dynamic development and rupture process.To compensate for the time and destruction terms,we propose the damage network model,which defines the physical interpretation of fractures through the ternary motif.We focus primarily on the evolution of node types,topological attributes,and motifs of the fracture network in limestone under uniaxial stress.Observations expose the varying behavior of the nodes'self-dynamics and neighbors'adjacent dynamics in the fracture network.This approach elucidates the impact of micro-crack behaviors on large brittle shear fractures from a topological perspective and further subdivides the progressive failure stage into four distinct phases(isolated crack growth phase,crack splay phase,damage coalescence phase,and mechanical failure phase)based on the significance profile of the motif.Regression analysis reveals a positive linear and negative power correlation between fracture network density and branch number to the rock damage resistance,respectively.The damage network model introduces a novel methodology for depicting the interaction of two-dimensional(2D)projected fractures,considering the dynamic spatiotemporal development characteristics and fracture geometric variation.It helps dynamically characterize properties such as connectivity,permeability,and damage factors while comprehensively assessing damage in rock mass fracture networks.
基金the China Research and Pilot Test on Key Technology of Efficient Production of Changqing Tight Oil(Grant No.2021DJ2202).
文摘Class III tight oil reservoirs have low porosity and permeability,which are often responsible for low production rates and limited recovery.Extensive repeated fracturing is a well-known technique to fix some of these issues.With such methods,existing fractures are refractured,and/or new fractures are created to facilitate communication with natural fractures.This study explored how different refracturing methods affect horizontal well fracture networks,with a special focus on morphology and related fluid flow changes.In particular,the study relied on the unconventional fracture model(UFM).The evolution of fracture morphology and flow field after the initial fracturing were analyzed accordingly.The simulation results indicated that increased formation energy and reduced reservoir stress differences can promote fracture expansion.It was shown that the length of the fracture network,the width of the fracture network,and the complexity of the fracture can be improved,the oil drainage area can be increased,the distance of oil and gas seepage can be reduced,and the production of a single well can be significantly increased.
基金the support from the National Key Research and Development Program of China(Grant No.2022YFE0137200)the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2024JC-YBQN-0381,2023JC-QN-0403)+2 种基金the Natural Science Basic Research Program of Shaanxi Province,China(Program No.2022JC-37)the Innovation Capability Support Program of Shaanxi(Program No.2023-CX-TD31)the Funded by Open Foundation of Shaanxi Key Laboratory of Carbon Dioxide Sequestration and Enhanced Oil Recovery,and the Youth Innovation Team of Shaanxi Universities。
文摘Proppant transport within fractures is one of the most critical tasks in oil,gas and geothermal reservoir stimulation,as it largely determines the ultimate performance of the operating well.Proppant transport in rough fracture networks is still a relatively new area of research and the associated transport mechanisms are still unclear.In this study,representative parameters of rough fracture surfaces formed by supercritical CO_(2) fracturing were used to generate a rough fracture network model based on a spectral synthesis method.Computational fluid dynamics(CFD)coupled with the discrete element method(DEM)was used to study proppant transport in this rough fracture network.To reveal the turning transport mechanism of proppants into branching fractures at the intersections of rough fracture networks,a comparison was made with the behavior within smooth fracture networks,and the effect of key pumping parameters on the proppant placement in a secondary fracture was analyzed.The results show that the transport behavior of proppant in rough fracture networks is very different from that of the one in the smooth fracture networks.The turning transport mechanisms of proppant into secondary fractures in rough fracture networks are gravity-driven sliding,high velocity fluid suspension,and fracture structure induction.Under the same injection conditions,supercritical CO_(2)with high flow Reynolds number still has a weaker ability to transport proppant into secondary fractures than water.Thickening of the supercritical CO_(2)needs to be increased beyond a certain value to have a significant effect on proppant carrying,and under the temperature and pressure conditions of this paper,it needs to be increased more than 20 times(about 0.94 m Pa s).Increasing the injection velocity and decreasing the proppant concentration facilitates the entry of proppant into the branching fractures,which in turn results in a larger stimulated reservoir volume.The results help to understand the proppant transport and placement process in rough fracture networks formed by reservoir stimulation,and provide a theoretical reference for the optimization of proppant pumping parameters in hydraulic fracturing.
基金funded by the National Natural Science Foundation of China(No.51974268)Open Fund of Key Laboratory of Ministry of Education for Improving Oil and Gas Recovery(NEPUEOR-2022-03)Research and Innovation Fund for Graduate Students of Southwest Petroleum University(No.2022KYCX005)。
文摘The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective method to improve oil recovery factor from unconventional oil reservoirs. Hydrocarbon gas huff-n-puff becomes preferable when the CO_(2) source is limited. However, the impact of complex fracture networks and well interference on the EOR performance of multiple MFHWs is still unclear. The optimal gas huff-n-puff parameters are significant for enhancing oil recovery. This work aims to optimize the hydrocarbon gas injection and production parameters for multiple MFHWs with complex fracture networks in unconventional oil reservoirs. Firstly, the numerical model based on unstructured grids is developed to characterize the complex fracture networks and capture the dynamic fracture features.Secondly, the PVT phase behavior simulation was carried out to provide the fluid model for numerical simulation. Thirdly, the optimal parameters for hydrocarbon gas huff-n-puff were obtained. Finally, the dominant factors of hydrocarbon gas huff-n-puff under complex fracture networks are obtained by fuzzy mathematical method. Results reveal that the current pressure of hydrocarbon gas injection can achieve miscible displacement. The optimal injection and production parameters are obtained by single-factor analysis to analyze the effect of individual parameter. Gas injection time is the dominant factor of hydrocarbon gas huff-n-puff in unconventional oil reservoirs with complex fracture networks. This work can offer engineers guidance for hydrocarbon gas huff-n-puff of multiple MFHWs considering the complex fracture networks.
基金sponsored by the General Program of the National Natural Science Foundation of China(Grant Nos.52079129 and 52209148)the Hubei Provincial General Fund,China(Grant No.2023AFB567)。
文摘Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.
基金SINOPEC's Scientific and Technological Research Project:Research on effective production strategies of Jurassic continental shale oil and gas(No.P21078-5).
文摘Pressure buildup testing can be used to analyze fracture network characteristics and conduct quantitative interpretation of relevant parameters for shale gas wells,thus providing bases for assessing the well productivity and formulating proper development strategies.This study establishes a new well test interpretation model for fractured horizontal wells based on seepage mechanisms of shale reservoirs and proposes a method for identifying fracturing patterns based on the characteristic slopes of pressure buildup curves and curve combination patterns.The pressure buildup curve patterns are identified to represent three types of shale reservoirs in the Sichuan Basin,namely the moderately deep shale reservoirs with high pressure,deep shale reservoirs with ultra-high pressure,and moderately deep shale reservoirs with normal pressure.Based on this,the relationship between the typical pressure buildup curve patterns and the fracture network types are put forward.Fracturing effects of three types of shale gas reservoir are compared and analyzed.The results show that typical flow patterns of shale reservoirs include bilinear flow in primary and secondary fractures,linear flow in secondary fractures,bilinear flow in secondary fractures and matrix,and linear flow in matrix.The fracture network characteristics can be determined using the characteristic slopes of pressure buildup curves and curve combinations.The linear flow in early secondary fractures is increasingly distinct with an increase in primary fracture conductivity.Moreover,the bilinear flow in secondary fractures and matrix and the subsequent linear flow in the matrix occur as the propping and density of secondary fractures increase.The increase in the burial depth,in-situ stress,and stress difference corresponds to a decrease in the propping of primary fractures that expand along different directions in the shale gas wells in the Sichuan Basin.Four pressure buildup curve patterns exist in the Sichuan Basin and its periphery.The pattern of pressure buildup curves of shale reservoirs in the Yongchuan area can be described as 1/2/→1/4,indicating limited stimulated reservoir volume,poorly propped secondary fractures,and the forming of primary fractures that extend only to certain directions.The pressure buildup curves of shale reservoirs in the main block of the Fuling area show a pattern of 1/4/→1/2 or 1/2,indicating greater stimulated reservoir volume,well propped secondary fractures,and the forming of complex fracture networks.The pattern of pressure buildup curves of shale reservoirs in the Pingqiao area is 1/2/→1/4→/1/2,indicating a fracturing effect somewhere between that of the Fuling and Yongchuan areas.For reservoirs with normal pressure,it is difficult to determine fracture network characteristics from pressure buildup curves due to insufficient formation energy and limited liquid drainage.
基金financially supported by,the Fundamental Research Funds for the Central Universities(Grant No.2023QN1064)the China Postdoctoral Science Foundation(Grant No.2023M733772)Jiangsu Funding Program for Excellent Postdoctoral Talent(Grant No.2023ZB847)。
文摘Prepulse combined hydraulic fracturing facilitates the development of fracture networks by integrating prepulse hydraulic loading with conventional hydraulic fracturing.The formation mechanisms of fracture networks between hydraulic and pre-existing fractures under different prepulse loading parameters remain unclear.This research investigates the impact of prepulse loading parameters,including the prepulse loading number ratio(C),prepulse loading stress ratio(S),and prepulse loading frequency(f),on the formation of fracture networks between hydraulic and pre-existing fractures,using both experimental and numerical methods.The results suggest that low prepulse loading stress ratios and high prepulse loading number ratios are advantageous loading modes.Multiple hydraulic fractures are generated in the specimen under the advantageous loading modes,facilitating the development of a complex fracture network.Fatigue damage occurs in the specimen at the prepulse loading stage.The high water pressure at the secondary conventional hydraulic fracturing promotes the growth of hydraulic fractures along the damage zones.This allows the hydraulic fractures to propagate deeply and interact with pre-existing fractures.Under advantageous loading conditions,multiple hydraulic fractures can extend to pre-existing fractures,and these hydraulic fractures penetrate or propagate along pre-existing fractures.Especially when the approach angle is large,the damage range in the specimen during the prepulse loading stage increases,resulting in the formation of more hydraulic fractures.
基金College Students Innovation and Entrepreneurship Project of Guangzhou Railway Polytechnic(2025CXCY015)。
文摘As the dominant seepage channel in rock masses,it is of great significance to study the influence of fracture roughness distribution on seepage and heat transfer in rock masses.In this paper,the fracture roughness distribution functions of the Bakhtiary dam site and Oskarshamn/Forsmark mountain were fitted using statistical methods.The COMSOL Multiphysics finite element software was utilized to analyze the effects of fracture roughness distribution types and empirical formulas for fracture hydraulic aperture on the seepage field and temperature field of rock masses.The results show that:(1)The fracture roughness at the Bakhtiary dam site and Oskarshamn/Forsmark mountain follows lognormal and normal distributions,respectively;(2)For rock masses with the same expected value and standard deviation of fracture roughness,the outflow from rock masses with lognormal distribution of fracture roughness is significantly larger than that of rock masses with normal distribution of fracture roughness;(3)The fracture hydraulic aperture,outflow,and cold front distance of the Li and Jiang model are significantly larger than those of the Barton model;(4)The outflow,hydraulic pressure distribution,and temperature distribution of the Barton model are more sensitive to the fracture roughness distribution type than those of the Li and Jiang model.
基金supported by the National Natural Science Foundation of China(Grant No.U23A6018)Science and Technology Program of Hebei(Grant No.E2022202041,2022HBQZYCXY004,242Q9920Z)the project of“Key technologies of seepage control system for large-scale hydraulic projects”was also gratefully appreciated.
文摘Curtain grouting projects are characterized by their large scale and complexity,presenting significant challenges for real-time prediction of grout penetration using traditional methods.This study introduces an intelligent prediction method for grouting in fractured rock masses based on three core principles:integration of multi-source input features,fracture voxel modeling,and shortest path in sequential grouting.Three categories of data(geological structure data,grouting environmental data,and grouting operation data in the concept of a grouting geological model)are integrated and served as multi-source structured data in the intelligent prediction of grouting.A voxelization model quantifies the spatial characteristics of fractures,with voxel size optimized for capturing grouting paths.A shortest path algorithm based on a hierarchical solution is then developed to calculate grout penetration distances in the process of sequential grouting.A complete analysis framework is established,from the voxelization of the fracture network model to precise voxel classification,ultimately achieving an accurate prediction of grout penetration.The method demonstrates excellent performance on the test set,with validation against numerical methods in single-fracture and sequential grouting scenarios confirming its accuracy and prediction efficiency as hundreds of times faster than numerical methods.Application to the Dongzhuang hydraulic project’s grouting test area further validates its effectiveness in multi-hole grouting scenarios.
文摘This study investigated the micro-sliding frictional behavior of shale in fracturing fluids under varying operational conditions using Chang 7 shale oil reservoir core samples.Through systematic micro-sliding friction experiments,the characteristics and governing mechanisms of shale friction were elucidated.Complementary analyses were conducted to characterize the mineral composition,petrophysical properties,and micromorphology of the shale samples,providing insights into the relationship between microscopic structure and frictional response.In this paper,the characteristics and variation law of shale micro-sliding friction under different types of graphite materials as additives in LGF-80(Low-damage Guar Fluid)oil flooding recoverable fracturing fluid system were mainly studied.In addition,the finite element numerical simulation experiment of hydraulic fracturing was adopted to study the influence of the friction coefficient of natural fracture surfaces on fracture propagation and formation of the fracture network.The geometric complexity of fracture networks was systematically quantified under varying frictional coefficients of natural fracture surfaces through multi-parametric characterization and morphometric analysis.The research results show that graphite micro-particles reduce friction and drag.Based on this,this paper proposes a new idea of graphite micro-particles as an additive in the LGF-80 oil flooding recoverable fracturing fluid system to reduce friction on the fracture surface.
基金Supported by the National Natural Science Foundation of China(52374004)National Key R&D Program of China(2023YFF0614102,2023YFE0110900).
文摘Based on the finite element-discrete element numerical method,a numerical model of fracture propagation in deflagration fracturing was established by considering the impact of stress wave,quasi-static pressure of explosive gas,and reflection of stress wave.The model was validated against the results of physical experiments.Taking the shale reservoirs of Silurian Longmaxi Formation in Luzhou area of the Sichuan Basin as an example,the effects of in-situ stress difference,natural fracture parameters,branch wellbore spacing,delay detonation time,and angle between branch wellbore and main wellbore on fracture propagation were identified.The results show that the fracture propagation morphology in deflagration fracturing is less affected by the in-situ stress difference when it is 5-15 MPa,and the tendency of fracture intersection between branch wellbores is significantly weakened when the in-situ stress difference reaches 20 MPa.The increase of natural fracture length promotes the fracture propagation along the natural fracture direction,while the increase of volumetric natural fracture density and angle limits the fracture propagation area and reduces the probability of fracture intersection between branch wells.The larger the branch wellbore spacing,the less probability of the fracture intersection between branch wells,allowing for the fracture propagation in multiple directions.Increasing the delay detonation time decreases the fracture spacing between branch wellbores.When the angle between the branch wellbore and the main wellbore is 45°and 90°,there is a tendency of fracture intersection between branch wellbores.
