The interaction between multiple fractures is important in the analysis of rock fracture propagation,fracture network evolution and stability and integrity of rocks under hydro-mechanical(HM)coupling conditions.At pre...The interaction between multiple fractures is important in the analysis of rock fracture propagation,fracture network evolution and stability and integrity of rocks under hydro-mechanical(HM)coupling conditions.At present,modeling the mechanical behavior of multiple fractures is still challenging.Under the condition of multiple fractures,the opening,closing,sliding,propagation and penetration of fractures become more complicated.In order to simulate the HM coupling behavior of multi-fracture system,the paper presents a novel numerical scheme,including mesh reconstruction and topology generation algorithm,to efficiently and accurately represent fractures and their propagation process,and a potential function-based algorithm to address contact problem.The fracture contact algorithm does not need to set contact pairs and thus is suitable for complex contact situations from small to large deformations induced by HM loading.The topology of fracture interfaces is constructed by the dynamic adding algorithm,which makes the mesh reconstruction more rapid in the modeling of fracturing process,especially in the case of multiple fractures intersections.The numerical scheme is implemented in CASRock,a self-developed numerical code,to simulate the propagation process of rock fractures and the interaction of multiple fractures under the condition of HM coupling.To verify the suitability of the code,a series of tests were performed.The code was then applied to simulate hydraulic fracture propagation and fracture interactions caused by fluid injection.The ability of this method to study fracture propagation,multi-fracture interaction and fracture network evolution under hydro-mechanical coupling conditions is demonstrated.展开更多
To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3 D multi-fracture propagation model was proposed. The model considers flu...To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3 D multi-fracture propagation model was proposed. The model considers fluid flow in the wellbore-perforation-fracture system and fluid leak-off into the rock matrix, and uses a 3 D boundary integral equation to describe the solid deformation. The solid-fluid coupling equation is solved by an explicit integration algorithm, and the fracture front is determined by the uniform tip asymptotic solutions and shortest path algorithm. The accuracy of the algorithm is verified by the analytical solution of radial fracture, results of the implicit level set algorithm, and results of organic glass fracturing experiment. Compared with the implicit level set algorithm(ILSA), the new algorithm is much higher in computation speed. The numerical case study is conducted based on a horizontal well in shale gas formation of Zhejiang oilfield. The impact of stress heterogeneity among multiple clusters and perforation number distribution on multi-fracture growth and fluid distribution among multiple fractures are analyzed by numerical simulation. The results show that reducing perforation number in each cluster can counteract the effect of stress contrast among perforation clusters. Adjusting perforation number in each cluster can promote uniform flux among clusters, and the perforation number difference should better be 1-2 among clusters. Increasing perforation number in the cluster with high in situ stress is conducive to uniform fluid partitioning. However, uniform fluid partitioning is not equivalent to uniform fracture geometry. The fracture geometry is controlled by the stress interference and horizontal principal stress profile jointly.展开更多
Due to the reservoir heterogeneity and the stress shadow effect, multiple hydraulic fractures within one fracturing segment cannot be initiated simultaneously and propagate evenly, which will cause a low effectiveness...Due to the reservoir heterogeneity and the stress shadow effect, multiple hydraulic fractures within one fracturing segment cannot be initiated simultaneously and propagate evenly, which will cause a low effectiveness of reservoir stimulation. Temporary plugging and diverting fracturing(TPDF) is considered to be a potential uniform-stimulation method for creating multiple fractures simultaneously in the oilfield. However, the multi-fracture propagation morphology during TPDF is not clear now. The purpose of this study is to quantitatively investigate the multi-fracture propagation morphology during TPDF through true tri-axial fracturing experiments and CT scanning. Critical parameters such as fracture spacing, number of perforation clusters, the viscosity of fracturing fluid, and the in-situ stress have been investigated. The fracture geometry before and after diversion have been quantitively analyzed based on the two-dimensional CT slices and three-dimensional reconstruction method. The main conclusions are as follows:(1) When injecting the high viscosity fluid or perforating at the location with low in-situ stress, multiple hydraulic fractures would simultaneously propagate. Otherwise, only one hydraulic fracture was created during the initial fracturing stage(IFS) for most tests.(2) The perforation cluster effectiveness(PCE) has increased from 26.62% during the IFS to 88.86% after using diverters.(3) The diverted fracture volume has no apparent correlation with the pressure peak and peak frequency during the diversion fracturing stage(DFS) but is positively correlated with water-work.(4) Four types of plugging behavior in shale could be controlled by adjusting the diverter recipe and diverter injection time, and the plugging behavior includes plugging the natural fracture in the wellbore, plugging the previous hydraulic fractures, plugging the fracture tip and plugging the bedding.展开更多
A forward model for optical fiber strain was established based on a planar 3D multi-fracture model. Then the forward method calculating distributed fiber strain induced by multi-fracture growth was proposed. Based on ...A forward model for optical fiber strain was established based on a planar 3D multi-fracture model. Then the forward method calculating distributed fiber strain induced by multi-fracture growth was proposed. Based on this method, fiber strain evolution during fracturing of the horizontal well was numerically simulated. Fiber strain evolution induced by fracture growth can be divided into three stages: strain increasing, shrinkage convergence, and straight-line convergence, whereas the evolution of fiber strain rate has four stages: strain rate increasing, shrinkage convergence, straight-line convergence, and strain rate reversal after pumping stops. Fiber strain does not flip after pumping stop, while the strain rate flips after pumping stop so that strain rate can reflect injection dynamics. The time when the fracture extends to the fiber and inter-well pressure channeling can be identified by the straight-line convergence band of distributed fiber strain or strain rate, and the non-uniform growth of multiple fractures can be evaluated by using the instants of fractures reaching the fiber monitoring well.When the horizontal section of the fiber monitoring well is within the height range of a hydraulic fracture, the instant of the fracture reaching the fiber can be identified;otherwise, the converging band is not apparent. In multi-stage fracturing, under the influence of stress shadow from previous fracturing stages, the tensile region of fiber strain may not appear, but the fiber strain rate can effectively show the fracture growth behavior in each stage. The evolution law of fiber strain rate in single-stage fracturing can be applied to multi-stage fracturing.展开更多
Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present u...Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present unique challenges due to their deep burial depth(4500-8882 m),low matrix permeability,complex crustal stress conditions,high temperature and pressure(HTHP,150-200℃,105-155 MPa),coupled with high salinity of formation water.Consequently,the costs associated with their exploitation and development are exceptionally high.In deep and ultra-deep reservoirs,hydraulic fracturing is commonly used to achieve high and stable production.During hydraulic fracturing,a substantial volume of fluid is injected into the reservoir.However,statistical analysis reveals that the flowback rate is typically less than 30%,leaving the majority of the fluid trapped within the reservoir.Therefore,hydraulic fracturing in deep reservoirs not only enhances the reservoir permeability by creating artificial fractures but also damages reservoirs due to the fracturing fluids involved.The challenging“three-high”environment of a deep reservoir,characterized by high temperature,high pressure,and high salinity,exacerbates conventional forms of damage,including water sensitivity,retention of fracturing fluids,rock creep,and proppant breakage.In addition,specific damage mechanisms come into play,such as fracturing fluid decomposition at elevated temperatures and proppant diagenetic reactions at HTHP conditions.Presently,the foremost concern in deep oil and gas development lies in effectively assessing the damage inflicted on these reservoirs by hydraulic fracturing,comprehending the underlying mechanisms,and selecting appropriate solutions.It's noteworthy that the majority of existing studies on reservoir damage primarily focus on conventional reservoirs,with limited attention given to deep reservoirs and a lack of systematic summaries.In light of this,our approach entails initially summarizing the current knowledge pertaining to the types of fracturing fluids employed in deep and ultra-deep reservoirs.Subsequently,we delve into a systematic examination of the damage processes and mechanisms caused by fracturing fluids within the context of hydraulic fracturing in deep reservoirs,taking into account the unique reservoir characteristics of high temperature,high pressure,and high in-situ stress.In addition,we provide an overview of research progress related to high-temperature deep reservoir fracturing fluid and the damage of aqueous fracturing fluids to rock matrix,both artificial and natural fractures,and sand-packed fractures.We conclude by offering a summary of current research advancements and future directions,which hold significant potential for facilitating the efficient development of deep oil and gas reservoirs while effectively mitigating reservoir damage.展开更多
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 C...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.展开更多
Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the...Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.展开更多
Polypropylene(PP)fiber-reinforced cement-based tailings backfill(FRCTB)is a green compound material with superior crack resistance and has good prospects for application in underground mining.However,FRCTB exhibits su...Polypropylene(PP)fiber-reinforced cement-based tailings backfill(FRCTB)is a green compound material with superior crack resistance and has good prospects for application in underground mining.However,FRCTB exhibits susceptibility to dynamic events,such as impact ground pressure and blast vibrations.This paper investigates the energy and crack distribution behavior of FRCTB under dynamic impact,considering the height/diameter(H/D)effect.Split Hopkinson pressure bar,industrial computed tomography scan,and scanning electron microscopy(SEM)experiments were carried out on six types of FRCTB.Laboratory outcomes confirmed fiber aggregation at the bottom of specimens.When H/D was less than 0.8,the proportion of PP fibers distributed along theθangle direction of80°-90°increased.For the total energy,all samples presented similar energy absorption,reflectance,and transmittance.However,a rise in H/D may cause a rise in the energy absorption rate of FRCTB during the peak phase.A positive correlation existed between the average strain rate and absorbed energy per unit volume.The increase in H/D resulted in a decreased crack volume fraction of FRCTB.When the H/D was greater than or equal to 0.7,the maximum crack volume fraction of FRCTB was observed close to the incidence plane.Radial cracks were present only in the FRCTB with an H/D ratio of 0.5.Samples with H/D ratios of 0.5 and 0.6 showed similar distributions of weakly and heavily damaged areas.PP fibers can limit the emergence and expansion of cracks by influencing their path.SEM observations revealed considerable differences in the bonding strengths between fibers and the FRCTB.Fibers that adhered particularly well to the substrate were attracted together with the hydration products adhering to surfaces.These results show that FRCTB is promising as a sustainable and green backfill for determining the design properties of mining with backfill.展开更多
In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The prim...In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The primary mechanism for gas breakthrough is the opening of micro-fractures due to high gas pressure.This occurs at gas pressures lower than the combined strength of the rock and its minimum principal stress under external loading conditions.To investigate the mechanism of microscale mode-I ruptures,it is essential to incorporate a multiscale approach that includes subcritical microcracks in the modeling framework.In this contribution,we derive the model from microstructures that contain periodically distributed microcracks within a porous material.The damage evolution law is coupled with the macroscopic poroelastic system by employing the asymptotic homogenization method and considering the inherent hydro-mechanical(HM)anisotropy at the microscale.The resulting permeability change induced by fracture opening is implicitly integrated into the gas flow equation.Verification examples are presented to validate the developed model step by step.An analysis of local macroscopic response is undertaken to underscore the influence of factors such as strain rate,initial damage,and applied stress,on the gas migration process.Numerical examples of direct tension tests are used to demonstrate the model’s efficacy in describing localized failure characteristics.Finally,the simulation results for preferential gas flow reveal the robustness of the two-scale model in explicitly depicting gas-induced fracturing in anisotropic clayey rocks.The model successfully captures the common behaviors observed in laboratory experiments,such as a sudden drop in gas injection pressure,rapid build-up of downstream gas pressure,and steady-state gas flow following gas breakthrough.展开更多
The Linxing area within the Ordos Basin exhibits pronounced reservoir heterogeneity and intricate micro-pore structures,rendering it susceptible to water-blocking damage during imbibition extraction.This study delved ...The Linxing area within the Ordos Basin exhibits pronounced reservoir heterogeneity and intricate micro-pore structures,rendering it susceptible to water-blocking damage during imbibition extraction.This study delved into the traits of tight sandstone reservoirs in the 8th member of the Shihezi Formation(also referred to as the He 8 Member)in the study area,as well as their effects on fracturing fluid imbibition.Utilizing experimental techniques such as nuclear magnetic resonance(NMR),high-pressure mercury intrusion(HPMI),and gas adsorption,this study elucidated the reservoir characteristics and examined the factors affecting the imbibition through imbibition experiments.The findings reveal that:①The reservoir,with average porosity of 8.40%and average permeability of 0.642×10^(-3)μm^(2),consists principally of quartz,feldspar,and lithic fragments,with feldspathic litharenite serving as the primary rock type and illite as the chief clay mineral;②Nano-scale micro-pores and throats dominate the reservoir,with dissolution pores and intercrystalline pores serving as predominant pore types,exhibiting relatively high pore connectivity;③Imbibition efficiency is influenced by petrophysical properties,clay mineral content,and microscopic pore structure.Due to the heterogeneity of the tight sandstone reservoir,microscopic factors have a more significant impact on the imbibition efficiency of fracturing fluids;④A comparative analysis shows that average pore size correlates most strongly with imbibition efficiency,followed by petrophysical properties and clay mineral content.In contrast,the pore type has minimal impact.Micropores are vital in the imbibition process,while meso-pores and macro-pores offer primary spaces for imbibition.This study offers theoretical insights and guidance for enhancing the post-fracturing production of tight sandstone reservoirs by examining the effects of these factors on the imbibition efficiency of fracturing fluids in tight sandstones.展开更多
Tight oil reservoirs are complex geological materials composed of solid matrix,pore structure,and mixed multiple phases of fluids,particularly for oil reservoirs suffering from high content of in situ pressurized wate...Tight oil reservoirs are complex geological materials composed of solid matrix,pore structure,and mixed multiple phases of fluids,particularly for oil reservoirs suffering from high content of in situ pressurized water found in China.In this regard,a coupled model considering two-phase flow of oil and water,as well as deformation and damage evolution of porous media,is proposed and validated using associated results,including the oil depletion process,analytical solution of stress shadow effect,and physical experiments on multi-fracture interactions and fracture propagation in unsaturated seepage fields.Then,the proposed model is used to study the behavior of multi-fracture interactions in an unsaturated reservoir in presence of water and oil.The results show that conspicuous interactions exist among multiple induced fractures.Interaction behavior varies from extracted geological profiles of the reservoir due to in situ stress anisotropy.The differential pressures of water and that of oil in different regions of reservoir affect interactions and trajectories of multi-fractures to a considerable degree.The absolute value of reservoir average pressure is a dominant factor affecting fracture interactions and in favor of enhancing fracture network complexity.In addition,difference of reservoir average pressures in different regions of reservoir would promote the fracturing effectiveness.Factors affecting fracture interactions and reservoir treatment effectiveness are quantitatively estimated through stimulated reservoir area.This study confirms the significance of incorporating the two-phase flow process in analyses of multifracture interactions and fracture trajectory predictions during tight sandstone oil reservoir developments.展开更多
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 fr...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.展开更多
Production decline analysis has been considered as an important method to obtain the flow parameters, reservoir properties and original gas in place. Although advanced Blasingame production decline analysis methods fo...Production decline analysis has been considered as an important method to obtain the flow parameters, reservoir properties and original gas in place. Although advanced Blasingame production decline analysis methods for vertical wells, fractured wells and horizontal wells are widely used, limited study has conducted on Blasingame production decline type curves for multi-fractured horizontal well(MFHW). Based on the perpendicular bisection(PEBI) grids, a numerical model was developed and the solution was obtained using control volume finite element method and the fully implicit method. Blasingame production decline-type curves of the infinitely conductive MFHW were plotted through computer programming. A field case was presented to analyse and verify the model developed. Five flow regimes, including early formation linear flow, early radial flow, compound linear flow, transient flow and pseudo-radial flow, are recognized. Fracture spacing is the main factor that affects early radial flow, compound linear flow and transient flow, the distance from the well to the circular boundary affects the pseudo-radial flow, and the type curves are also significantly affected by the formation permeability, fracture number and fracture half-length. The validation of field case suggests that the Blasingame production decline type curves proposed in this work can be applied to the production decline analysis for MFHW in tight gas reservoirs.展开更多
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 ...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.展开更多
Fracture propagation in shale under in situ conditions is a critical but poorly understood mechanical process in hydraulic fracturing for deep shale gas reservoirs. To address this, hydraulic fracturing experiments we...Fracture propagation in shale under in situ conditions is a critical but poorly understood mechanical process in hydraulic fracturing for deep shale gas reservoirs. To address this, hydraulic fracturing experiments were conducted on hollow double-wing crack specimens of the Longmaxi shale under conditions simulating the ground surface(confining pressure σ_(cp)=0, room temperature(Tr)) and at depths of 1600 m(σ_(cp)=40 MPa, Ti=70 ℃) and 3300 m(σ_(cp)=80 MPa, high temperature Ti=110 ℃) in the study area.High in situ stress was found to significantly increase fracture toughness through constrained microcracking and particle frictional bridging mechanisms. Increasing the temperature enhances rather than weakens the fracture resistance because it increases the grain debonding length, which dissipates more plastic energy and enlarges grains to close microdefects and generate compressive stress to inhibit microcracking. Interestingly, the fracture toughness anisotropy in the shale was found to be nearly constant across burial depths, despite reported variations with increasing confining pressure. Heated water was not found to be as important as the in situ environment in influencing shale fracture. These findings emphasize the need to test the fracture toughness of deep shales under coupled in situ stress and temperature conditions rather than focusing on either in situ stress or temperature alone.展开更多
Pulsating hydraulic fracturing(PHF)is a promising fracturing method and can generate a dynamic periodic pressure.The periodic pressure can induce fatigue failure of rocks and decrease initiation pressure of fracture.I...Pulsating hydraulic fracturing(PHF)is a promising fracturing method and can generate a dynamic periodic pressure.The periodic pressure can induce fatigue failure of rocks and decrease initiation pressure of fracture.If the frequency of periodic pressure exceeds 10 Hz,the distribution of pressure along the main fracture will be heterogeneous,which is much different from the one induced by the common fracturing method.In this study,the impact of this special spatial feature of pressure on hydraulic fracture is mainly investigated.A coupled numerical simulation model is first proposed and verified through experimental and theoretical solutions.The mechanism of secondary fracture initiation around the main fracture is then discovered.In addition,sensitivity studies are conducted to find out the application potential of this new method.The results show that(1)this coupled numerical simulation model is accurate.Through comparison with experimental and theoretical data,the average error of this coupled model is less than 1.01%.(2)Even if a reservoir has no natural fracture,this heterogeneous distribution pressure can also cause many secondary fractures around the main fracture.(3)The mechanism of secondary fracture initiation is that this heterogeneous distribution pressure causes tensile stress at many locations along the main fracture.(4)Through adjusting the stimulation parameters,the stimulation efficiency can be improved.The average and amplitude of pressure can increase possibility of secondary fracture initiation.The frequency of this periodic pressure can increase number of secondary fractures.Even 6 secondary fractures along a 100 m-length main fracture can be generated.(5)The influence magnitudes of stimulation parameters are larger than ones of geomechanical properties,therefore,this new fracturing method has a wide application potential.展开更多
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...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.展开更多
Determining the timing of fracturing is crucial for understanding reservoir evolution and hydrocarbon accumulation in foreland basins.Using fracturing data from cores,borehole images,and outcrops,combined with the clu...Determining the timing of fracturing is crucial for understanding reservoir evolution and hydrocarbon accumulation in foreland basins.Using fracturing data from cores,borehole images,and outcrops,combined with the clumped isotope(D47)and fluid inclusion analyses of carbonate minerals filled in pores and fractures,this study ascertained the fracturing timing of the Jurassic reservoirs in the Dibei-Tuziluoke Gas Field,Kuqa Foreland Basin.Data from outcrops and borehole images show two dominant fracture sets in the study area:W-E and NE-SW striking fractures.Some W-E striking fractures are carbonate-filled,while NE-SW striking fractures lack mineral fillings.Bitumen veins,not easy to be identified in borehole images,are prevalent in cores.The petrographic analysis reveals that these bitumen veins formed before the calcite cementation in pores and display high viscosity and low maturity.Homogenization temperatures(T_(h))from primary fluid inclusion assemblages in two representative calcite vein samples were notably lower than T_(△47) values from corresponding samples.This suggests the △_(47) signature underwent alteration due to partial reordering during burial.Thus,△_(47)-derived temperatures(apparent temperatures)may not faithfully represent the mineral precipitation temperatures.When plotting these apparent temperatures vs.the burial history,only the possible latest ages of fracturing emerged.These ages were further refined by considering petroleum charging,tectonic evolution,and stress orientation.Bitumen-filled fractures likely resulted from the Late Cretaceous uplift,marking the migration of low-maturity hydrocarbons in the study area.Carbonate-filled E-W striking fractures emerged during the late Miocene(~13-6.5 Ma)alongside fold development.NE-striking fractures that crosscut W-E ones possibly formed recently due to stress reorientation.展开更多
An optimization method of fracturing fluid volume strength was introduced taking well X-1 in Biyang Sag of Nanxiang Basin as an example.The characteristic curves of capillary pressure and relative permeability were ob...An optimization method of fracturing fluid volume strength was introduced taking well X-1 in Biyang Sag of Nanxiang Basin as an example.The characteristic curves of capillary pressure and relative permeability were obtained from history matching between forced imbibition experimental data and core-scale reservoir simulation results and taken into a large scale reservoir model to mimic the forced imbibition behavior during the well shut-in period after fracturing.The optimization of the stimulated reservoir volume(SRV)fracturing fluid volume strength should meet the requirements of estimated ultimate recovery(EUR),increased oil recovery by forced imbibition and enhancement of formation pressure and the fluid volume strength of fracturing fluid should be controlled around a critical value to avoid either insufficiency of imbibition displacement caused by insufficient fluid amount or increase of costs and potential formation damage caused by excessive fluid amount.Reservoir simulation results showed that SRV fracturing fluid volume strength positively correlated with single-well EUR and an optimal fluid volume strength existed,above which the single-well EUR increase rate kept decreasing.An optimized increase of SRV fracturing fluid volume and shut-in time would effectively increase the formation pressure and enhance well production.Field test results of well X-1 proved the practicality of established optimization method of SRV fracturing fluid volume strength on significant enhancement of shale oil well production.展开更多
This study introduces a novel method integrating CO_(2)flooding with radial borehole fracturing for enhanced oil recovery and CO_(2)underground storage,a solution to the limited vertical stimulation reservoir volume i...This study introduces a novel method integrating CO_(2)flooding with radial borehole fracturing for enhanced oil recovery and CO_(2)underground storage,a solution to the limited vertical stimulation reservoir volume in horizontal well fracturing.A numerical model is established to investigate the production rate,reservoir pressure field,and CO_(2)saturation distribution corresponding to changing time of CO_(2)flooding with radial borehole fracturing.A sensitivity analysis on the influence of CO_(2)injection location,layer spacing,pressure difference,borehole number,and hydraulic fractures on oil production and CO_(2)storage is conducted.The CO_(2)flooding process is divided into four stages.Reductions in layer spacing will significantly improve oil production rate and gas storage capacity.However,serious gas channeling can occur when the spacing is lower than 20 m.Increasing the pressure difference between the producer and injector,the borehole number,the hydraulic fracture height,and the fracture width can also increase the oil production rate and gas storage rate.Sensitivity analysis shows that layer spacing and fracture height greatly influence gas storage and oil production.Research outcomes are expected to provide a theoretical basis for the efficient development of shale oil reservoirs in the vertical direction.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos.52125903).
文摘The interaction between multiple fractures is important in the analysis of rock fracture propagation,fracture network evolution and stability and integrity of rocks under hydro-mechanical(HM)coupling conditions.At present,modeling the mechanical behavior of multiple fractures is still challenging.Under the condition of multiple fractures,the opening,closing,sliding,propagation and penetration of fractures become more complicated.In order to simulate the HM coupling behavior of multi-fracture system,the paper presents a novel numerical scheme,including mesh reconstruction and topology generation algorithm,to efficiently and accurately represent fractures and their propagation process,and a potential function-based algorithm to address contact problem.The fracture contact algorithm does not need to set contact pairs and thus is suitable for complex contact situations from small to large deformations induced by HM loading.The topology of fracture interfaces is constructed by the dynamic adding algorithm,which makes the mesh reconstruction more rapid in the modeling of fracturing process,especially in the case of multiple fractures intersections.The numerical scheme is implemented in CASRock,a self-developed numerical code,to simulate the propagation process of rock fractures and the interaction of multiple fractures under the condition of HM coupling.To verify the suitability of the code,a series of tests were performed.The code was then applied to simulate hydraulic fracture propagation and fracture interactions caused by fluid injection.The ability of this method to study fracture propagation,multi-fracture interaction and fracture network evolution under hydro-mechanical coupling conditions is demonstrated.
基金China National Science and Technology Major Project(2016ZX05023).
文摘To resolve the issue of design for multi-stage and multi-cluster fracturing in multi-zone reservoirs, a new efficient algorithm for the planar 3 D multi-fracture propagation model was proposed. The model considers fluid flow in the wellbore-perforation-fracture system and fluid leak-off into the rock matrix, and uses a 3 D boundary integral equation to describe the solid deformation. The solid-fluid coupling equation is solved by an explicit integration algorithm, and the fracture front is determined by the uniform tip asymptotic solutions and shortest path algorithm. The accuracy of the algorithm is verified by the analytical solution of radial fracture, results of the implicit level set algorithm, and results of organic glass fracturing experiment. Compared with the implicit level set algorithm(ILSA), the new algorithm is much higher in computation speed. The numerical case study is conducted based on a horizontal well in shale gas formation of Zhejiang oilfield. The impact of stress heterogeneity among multiple clusters and perforation number distribution on multi-fracture growth and fluid distribution among multiple fractures are analyzed by numerical simulation. The results show that reducing perforation number in each cluster can counteract the effect of stress contrast among perforation clusters. Adjusting perforation number in each cluster can promote uniform flux among clusters, and the perforation number difference should better be 1-2 among clusters. Increasing perforation number in the cluster with high in situ stress is conducive to uniform fluid partitioning. However, uniform fluid partitioning is not equivalent to uniform fracture geometry. The fracture geometry is controlled by the stress interference and horizontal principal stress profile jointly.
基金the National Natural Science Foundation of China fund (Project number: 52174045 and No. 52104011)Research Foundation of China University of Petroleum-Beijing at Karamay (No. XQZX20210001)PetroChina Innovation Foundation (2020D50070207)。
文摘Due to the reservoir heterogeneity and the stress shadow effect, multiple hydraulic fractures within one fracturing segment cannot be initiated simultaneously and propagate evenly, which will cause a low effectiveness of reservoir stimulation. Temporary plugging and diverting fracturing(TPDF) is considered to be a potential uniform-stimulation method for creating multiple fractures simultaneously in the oilfield. However, the multi-fracture propagation morphology during TPDF is not clear now. The purpose of this study is to quantitatively investigate the multi-fracture propagation morphology during TPDF through true tri-axial fracturing experiments and CT scanning. Critical parameters such as fracture spacing, number of perforation clusters, the viscosity of fracturing fluid, and the in-situ stress have been investigated. The fracture geometry before and after diversion have been quantitively analyzed based on the two-dimensional CT slices and three-dimensional reconstruction method. The main conclusions are as follows:(1) When injecting the high viscosity fluid or perforating at the location with low in-situ stress, multiple hydraulic fractures would simultaneously propagate. Otherwise, only one hydraulic fracture was created during the initial fracturing stage(IFS) for most tests.(2) The perforation cluster effectiveness(PCE) has increased from 26.62% during the IFS to 88.86% after using diverters.(3) The diverted fracture volume has no apparent correlation with the pressure peak and peak frequency during the diversion fracturing stage(DFS) but is positively correlated with water-work.(4) Four types of plugging behavior in shale could be controlled by adjusting the diverter recipe and diverter injection time, and the plugging behavior includes plugging the natural fracture in the wellbore, plugging the previous hydraulic fractures, plugging the fracture tip and plugging the bedding.
基金Supported by the National Natural Science Foundation of China(52104060)Project of Natural Science Outstanding Youth Foundation of Shandong(ZR2019QEE005)。
文摘A forward model for optical fiber strain was established based on a planar 3D multi-fracture model. Then the forward method calculating distributed fiber strain induced by multi-fracture growth was proposed. Based on this method, fiber strain evolution during fracturing of the horizontal well was numerically simulated. Fiber strain evolution induced by fracture growth can be divided into three stages: strain increasing, shrinkage convergence, and straight-line convergence, whereas the evolution of fiber strain rate has four stages: strain rate increasing, shrinkage convergence, straight-line convergence, and strain rate reversal after pumping stops. Fiber strain does not flip after pumping stop, while the strain rate flips after pumping stop so that strain rate can reflect injection dynamics. The time when the fracture extends to the fiber and inter-well pressure channeling can be identified by the straight-line convergence band of distributed fiber strain or strain rate, and the non-uniform growth of multiple fractures can be evaluated by using the instants of fractures reaching the fiber monitoring well.When the horizontal section of the fiber monitoring well is within the height range of a hydraulic fracture, the instant of the fracture reaching the fiber can be identified;otherwise, the converging band is not apparent. In multi-stage fracturing, under the influence of stress shadow from previous fracturing stages, the tensile region of fiber strain may not appear, but the fiber strain rate can effectively show the fracture growth behavior in each stage. The evolution law of fiber strain rate in single-stage fracturing can be applied to multi-stage fracturing.
基金Dao-Bing Wang was supported by the Beijing Natural Science Foundation Project(No.3222030)the National Natural Science Foundation of China(No.52274002)+1 种基金the PetroChina Science and Technology Innovation Foundation Project(No.2021DQ02-0201)Fu-Jian Zhou was supported by the National Natural Science Foundation of China(No.52174045).
文摘Deep and ultra-deep reservoirs have gradually become the primary focus of hydrocarbon exploration as a result of a series of significant discoveries in deep hydrocarbon exploration worldwide.These reservoirs present unique challenges due to their deep burial depth(4500-8882 m),low matrix permeability,complex crustal stress conditions,high temperature and pressure(HTHP,150-200℃,105-155 MPa),coupled with high salinity of formation water.Consequently,the costs associated with their exploitation and development are exceptionally high.In deep and ultra-deep reservoirs,hydraulic fracturing is commonly used to achieve high and stable production.During hydraulic fracturing,a substantial volume of fluid is injected into the reservoir.However,statistical analysis reveals that the flowback rate is typically less than 30%,leaving the majority of the fluid trapped within the reservoir.Therefore,hydraulic fracturing in deep reservoirs not only enhances the reservoir permeability by creating artificial fractures but also damages reservoirs due to the fracturing fluids involved.The challenging“three-high”environment of a deep reservoir,characterized by high temperature,high pressure,and high salinity,exacerbates conventional forms of damage,including water sensitivity,retention of fracturing fluids,rock creep,and proppant breakage.In addition,specific damage mechanisms come into play,such as fracturing fluid decomposition at elevated temperatures and proppant diagenetic reactions at HTHP conditions.Presently,the foremost concern in deep oil and gas development lies in effectively assessing the damage inflicted on these reservoirs by hydraulic fracturing,comprehending the underlying mechanisms,and selecting appropriate solutions.It's noteworthy that the majority of existing studies on reservoir damage primarily focus on conventional reservoirs,with limited attention given to deep reservoirs and a lack of systematic summaries.In light of this,our approach entails initially summarizing the current knowledge pertaining to the types of fracturing fluids employed in deep and ultra-deep reservoirs.Subsequently,we delve into a systematic examination of the damage processes and mechanisms caused by fracturing fluids within the context of hydraulic fracturing in deep reservoirs,taking into account the unique reservoir characteristics of high temperature,high pressure,and high in-situ stress.In addition,we provide an overview of research progress related to high-temperature deep reservoir fracturing fluid and the damage of aqueous fracturing fluids to rock matrix,both artificial and natural fractures,and sand-packed fractures.We conclude by offering a summary of current research advancements and future directions,which hold significant potential for facilitating the efficient development of deep oil and gas reservoirs while effectively mitigating reservoir damage.
基金the funding support from the National Natural Science Foundation of China(Nos.52274093 and 52004236).
文摘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.
基金Supported by the National Natural Science Foundation of China Project(52274014)Comprehensive Scientific Research Project of China National Offshore Oil Corporation(KJZH-2023-2303)。
文摘Deep coal seams show low permeability,low elastic modulus,high Poisson’s ratio,strong plasticity,high fracture initiation pressure,difficulty in fracture extension,and difficulty in proppants addition.We proposed the concept of large-scale stimulation by fracture network,balanced propagation and effective support of fracture network in fracturing design and developed the extreme massive hydraulic fracturing technique for deep coalbed methane(CBM)horizontal wells.This technique involves massive injection with high pumping rate+high-intensity proppant injection+perforation with equal apertures and limited flow+temporary plugging and diverting fractures+slick water with integrated variable viscosity+graded proppants with multiple sizes.The technique was applied in the pioneering test of a multi-stage fracturing horizontal well in deep CBM of Linxing Block,eastern margin of the Ordos Basin.The injection flow rate is 18 m^(3)/min,proppant intensity is 2.1 m^(3)/m,and fracturing fluid intensity is 16.5 m^(3)/m.After fracturing,a complex fracture network was formed,with an average fracture length of 205 m.The stimulated reservoir volume was 1987×10^(4)m^(3),and the peak gas production rate reached 6.0×10^(4)m^(3)/d,which achieved efficient development of deep CBM.
基金financially supported by the National Key Research and Development Program of China(No.2022YFC2905004)the China Postdoctoral Science Foundation(No.2023M742134)。
文摘Polypropylene(PP)fiber-reinforced cement-based tailings backfill(FRCTB)is a green compound material with superior crack resistance and has good prospects for application in underground mining.However,FRCTB exhibits susceptibility to dynamic events,such as impact ground pressure and blast vibrations.This paper investigates the energy and crack distribution behavior of FRCTB under dynamic impact,considering the height/diameter(H/D)effect.Split Hopkinson pressure bar,industrial computed tomography scan,and scanning electron microscopy(SEM)experiments were carried out on six types of FRCTB.Laboratory outcomes confirmed fiber aggregation at the bottom of specimens.When H/D was less than 0.8,the proportion of PP fibers distributed along theθangle direction of80°-90°increased.For the total energy,all samples presented similar energy absorption,reflectance,and transmittance.However,a rise in H/D may cause a rise in the energy absorption rate of FRCTB during the peak phase.A positive correlation existed between the average strain rate and absorbed energy per unit volume.The increase in H/D resulted in a decreased crack volume fraction of FRCTB.When the H/D was greater than or equal to 0.7,the maximum crack volume fraction of FRCTB was observed close to the incidence plane.Radial cracks were present only in the FRCTB with an H/D ratio of 0.5.Samples with H/D ratios of 0.5 and 0.6 showed similar distributions of weakly and heavily damaged areas.PP fibers can limit the emergence and expansion of cracks by influencing their path.SEM observations revealed considerable differences in the bonding strengths between fibers and the FRCTB.Fibers that adhered particularly well to the substrate were attracted together with the hydration products adhering to surfaces.These results show that FRCTB is promising as a sustainable and green backfill for determining the design properties of mining with backfill.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.12302503 and U20A20266)Scientific and Technological Research Projects in Sichuan Province,China(Grant No.2023ZYD0154).
文摘In the context of repositories for nuclear waste,understanding the behavior of gas migration through clayey rocks with inherent anisotropy is crucial for assessing the safety of geological disposal facilities.The primary mechanism for gas breakthrough is the opening of micro-fractures due to high gas pressure.This occurs at gas pressures lower than the combined strength of the rock and its minimum principal stress under external loading conditions.To investigate the mechanism of microscale mode-I ruptures,it is essential to incorporate a multiscale approach that includes subcritical microcracks in the modeling framework.In this contribution,we derive the model from microstructures that contain periodically distributed microcracks within a porous material.The damage evolution law is coupled with the macroscopic poroelastic system by employing the asymptotic homogenization method and considering the inherent hydro-mechanical(HM)anisotropy at the microscale.The resulting permeability change induced by fracture opening is implicitly integrated into the gas flow equation.Verification examples are presented to validate the developed model step by step.An analysis of local macroscopic response is undertaken to underscore the influence of factors such as strain rate,initial damage,and applied stress,on the gas migration process.Numerical examples of direct tension tests are used to demonstrate the model’s efficacy in describing localized failure characteristics.Finally,the simulation results for preferential gas flow reveal the robustness of the two-scale model in explicitly depicting gas-induced fracturing in anisotropic clayey rocks.The model successfully captures the common behaviors observed in laboratory experiments,such as a sudden drop in gas injection pressure,rapid build-up of downstream gas pressure,and steady-state gas flow following gas breakthrough.
基金funded by the National key R&D Program of China(No.2023YFE0120700)the National Natural Science Foundation of China(No.51934005)+2 种基金the Shaanxi Province 2023 Innovation Capability Support Plan(No.2023KJXX-122)the Technology Innovation Leading Program of Shaanxi(No.2022 PT-08)the Project of Youth Innovation Team of Shaanxi Universities(No.22JP063).
文摘The Linxing area within the Ordos Basin exhibits pronounced reservoir heterogeneity and intricate micro-pore structures,rendering it susceptible to water-blocking damage during imbibition extraction.This study delved into the traits of tight sandstone reservoirs in the 8th member of the Shihezi Formation(also referred to as the He 8 Member)in the study area,as well as their effects on fracturing fluid imbibition.Utilizing experimental techniques such as nuclear magnetic resonance(NMR),high-pressure mercury intrusion(HPMI),and gas adsorption,this study elucidated the reservoir characteristics and examined the factors affecting the imbibition through imbibition experiments.The findings reveal that:①The reservoir,with average porosity of 8.40%and average permeability of 0.642×10^(-3)μm^(2),consists principally of quartz,feldspar,and lithic fragments,with feldspathic litharenite serving as the primary rock type and illite as the chief clay mineral;②Nano-scale micro-pores and throats dominate the reservoir,with dissolution pores and intercrystalline pores serving as predominant pore types,exhibiting relatively high pore connectivity;③Imbibition efficiency is influenced by petrophysical properties,clay mineral content,and microscopic pore structure.Due to the heterogeneity of the tight sandstone reservoir,microscopic factors have a more significant impact on the imbibition efficiency of fracturing fluids;④A comparative analysis shows that average pore size correlates most strongly with imbibition efficiency,followed by petrophysical properties and clay mineral content.In contrast,the pore type has minimal impact.Micropores are vital in the imbibition process,while meso-pores and macro-pores offer primary spaces for imbibition.This study offers theoretical insights and guidance for enhancing the post-fracturing production of tight sandstone reservoirs by examining the effects of these factors on the imbibition efficiency of fracturing fluids in tight sandstones.
基金funded by National Natural Science Foundation of China(Grant Nos.51761135102 and 51525402)the Fundamental Research Funds for the Central Universities(Grant No.N180105029)。
文摘Tight oil reservoirs are complex geological materials composed of solid matrix,pore structure,and mixed multiple phases of fluids,particularly for oil reservoirs suffering from high content of in situ pressurized water found in China.In this regard,a coupled model considering two-phase flow of oil and water,as well as deformation and damage evolution of porous media,is proposed and validated using associated results,including the oil depletion process,analytical solution of stress shadow effect,and physical experiments on multi-fracture interactions and fracture propagation in unsaturated seepage fields.Then,the proposed model is used to study the behavior of multi-fracture interactions in an unsaturated reservoir in presence of water and oil.The results show that conspicuous interactions exist among multiple induced fractures.Interaction behavior varies from extracted geological profiles of the reservoir due to in situ stress anisotropy.The differential pressures of water and that of oil in different regions of reservoir affect interactions and trajectories of multi-fractures to a considerable degree.The absolute value of reservoir average pressure is a dominant factor affecting fracture interactions and in favor of enhancing fracture network complexity.In addition,difference of reservoir average pressures in different regions of reservoir would promote the fracturing effectiveness.Factors affecting fracture interactions and reservoir treatment effectiveness are quantitatively estimated through stimulated reservoir area.This study confirms the significance of incorporating the two-phase flow process in analyses of multifracture interactions and fracture trajectory predictions during tight sandstone oil reservoir developments.
文摘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.
基金Project(2013CB228005)supported by the National Basic Research Program of China
文摘Production decline analysis has been considered as an important method to obtain the flow parameters, reservoir properties and original gas in place. Although advanced Blasingame production decline analysis methods for vertical wells, fractured wells and horizontal wells are widely used, limited study has conducted on Blasingame production decline type curves for multi-fractured horizontal well(MFHW). Based on the perpendicular bisection(PEBI) grids, a numerical model was developed and the solution was obtained using control volume finite element method and the fully implicit method. Blasingame production decline-type curves of the infinitely conductive MFHW were plotted through computer programming. A field case was presented to analyse and verify the model developed. Five flow regimes, including early formation linear flow, early radial flow, compound linear flow, transient flow and pseudo-radial flow, are recognized. Fracture spacing is the main factor that affects early radial flow, compound linear flow and transient flow, the distance from the well to the circular boundary affects the pseudo-radial flow, and the type curves are also significantly affected by the formation permeability, fracture number and fracture half-length. The validation of field case suggests that the Blasingame production decline type curves proposed in this work can be applied to the production decline analysis for MFHW in tight gas reservoirs.
基金National Natural Science Foundation of China (51974176, 52174194, 51934004)Shandong Provincial Colleges and Universities Youth Innovation and Technology Support Program (2019KJH006)+1 种基金Taishan Scholars Project (TS20190935)Shandong outstanding youth fund (ZR2020JQ22).
文摘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.
基金supported by the National Natural Science Foundation of China(No.12172240).
文摘Fracture propagation in shale under in situ conditions is a critical but poorly understood mechanical process in hydraulic fracturing for deep shale gas reservoirs. To address this, hydraulic fracturing experiments were conducted on hollow double-wing crack specimens of the Longmaxi shale under conditions simulating the ground surface(confining pressure σ_(cp)=0, room temperature(Tr)) and at depths of 1600 m(σ_(cp)=40 MPa, Ti=70 ℃) and 3300 m(σ_(cp)=80 MPa, high temperature Ti=110 ℃) in the study area.High in situ stress was found to significantly increase fracture toughness through constrained microcracking and particle frictional bridging mechanisms. Increasing the temperature enhances rather than weakens the fracture resistance because it increases the grain debonding length, which dissipates more plastic energy and enlarges grains to close microdefects and generate compressive stress to inhibit microcracking. Interestingly, the fracture toughness anisotropy in the shale was found to be nearly constant across burial depths, despite reported variations with increasing confining pressure. Heated water was not found to be as important as the in situ environment in influencing shale fracture. These findings emphasize the need to test the fracture toughness of deep shales under coupled in situ stress and temperature conditions rather than focusing on either in situ stress or temperature alone.
基金supported by the National Natural Science Foundation of China (Grant No.52004302)Science Foundation of China University of Petroleum,Beijing (No.2462021YXZZ012)the Strategic Cooperation Technology Projects of CNPC and CUPB (ZLZX2020-01)。
文摘Pulsating hydraulic fracturing(PHF)is a promising fracturing method and can generate a dynamic periodic pressure.The periodic pressure can induce fatigue failure of rocks and decrease initiation pressure of fracture.If the frequency of periodic pressure exceeds 10 Hz,the distribution of pressure along the main fracture will be heterogeneous,which is much different from the one induced by the common fracturing method.In this study,the impact of this special spatial feature of pressure on hydraulic fracture is mainly investigated.A coupled numerical simulation model is first proposed and verified through experimental and theoretical solutions.The mechanism of secondary fracture initiation around the main fracture is then discovered.In addition,sensitivity studies are conducted to find out the application potential of this new method.The results show that(1)this coupled numerical simulation model is accurate.Through comparison with experimental and theoretical data,the average error of this coupled model is less than 1.01%.(2)Even if a reservoir has no natural fracture,this heterogeneous distribution pressure can also cause many secondary fractures around the main fracture.(3)The mechanism of secondary fracture initiation is that this heterogeneous distribution pressure causes tensile stress at many locations along the main fracture.(4)Through adjusting the stimulation parameters,the stimulation efficiency can be improved.The average and amplitude of pressure can increase possibility of secondary fracture initiation.The frequency of this periodic pressure can increase number of secondary fractures.Even 6 secondary fractures along a 100 m-length main fracture can be generated.(5)The influence magnitudes of stimulation parameters are larger than ones of geomechanical properties,therefore,this new fracturing method has a wide application potential.
基金The authors gratefully acknowledge the financial support received from the Strategic Cooperation Technology Projects of CNPC and CUPB(No.ZLZX2020-01)the National Key Scientific Research Instrument Research Project of NSFC(No.51827804).
文摘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.
基金funded by the PetroChina Major Research Program on Deep Petroleum System in the Tarim Basin(No.ZD 2019-183-01-003)the Major Research Project on the Tethys Geodynamic System from the National Natural Science Foundation of China(No.92055204)the National Natural Science Foundation of China(No.42072134).
文摘Determining the timing of fracturing is crucial for understanding reservoir evolution and hydrocarbon accumulation in foreland basins.Using fracturing data from cores,borehole images,and outcrops,combined with the clumped isotope(D47)and fluid inclusion analyses of carbonate minerals filled in pores and fractures,this study ascertained the fracturing timing of the Jurassic reservoirs in the Dibei-Tuziluoke Gas Field,Kuqa Foreland Basin.Data from outcrops and borehole images show two dominant fracture sets in the study area:W-E and NE-SW striking fractures.Some W-E striking fractures are carbonate-filled,while NE-SW striking fractures lack mineral fillings.Bitumen veins,not easy to be identified in borehole images,are prevalent in cores.The petrographic analysis reveals that these bitumen veins formed before the calcite cementation in pores and display high viscosity and low maturity.Homogenization temperatures(T_(h))from primary fluid inclusion assemblages in two representative calcite vein samples were notably lower than T_(△47) values from corresponding samples.This suggests the △_(47) signature underwent alteration due to partial reordering during burial.Thus,△_(47)-derived temperatures(apparent temperatures)may not faithfully represent the mineral precipitation temperatures.When plotting these apparent temperatures vs.the burial history,only the possible latest ages of fracturing emerged.These ages were further refined by considering petroleum charging,tectonic evolution,and stress orientation.Bitumen-filled fractures likely resulted from the Late Cretaceous uplift,marking the migration of low-maturity hydrocarbons in the study area.Carbonate-filled E-W striking fractures emerged during the late Miocene(~13-6.5 Ma)alongside fold development.NE-striking fractures that crosscut W-E ones possibly formed recently due to stress reorientation.
文摘An optimization method of fracturing fluid volume strength was introduced taking well X-1 in Biyang Sag of Nanxiang Basin as an example.The characteristic curves of capillary pressure and relative permeability were obtained from history matching between forced imbibition experimental data and core-scale reservoir simulation results and taken into a large scale reservoir model to mimic the forced imbibition behavior during the well shut-in period after fracturing.The optimization of the stimulated reservoir volume(SRV)fracturing fluid volume strength should meet the requirements of estimated ultimate recovery(EUR),increased oil recovery by forced imbibition and enhancement of formation pressure and the fluid volume strength of fracturing fluid should be controlled around a critical value to avoid either insufficiency of imbibition displacement caused by insufficient fluid amount or increase of costs and potential formation damage caused by excessive fluid amount.Reservoir simulation results showed that SRV fracturing fluid volume strength positively correlated with single-well EUR and an optimal fluid volume strength existed,above which the single-well EUR increase rate kept decreasing.An optimized increase of SRV fracturing fluid volume and shut-in time would effectively increase the formation pressure and enhance well production.Field test results of well X-1 proved the practicality of established optimization method of SRV fracturing fluid volume strength on significant enhancement of shale oil well production.
基金This study has been funded by the National Science Fund for Distinguished Young Scholars(No.52204063)Science Foundation of China University of Petroleum,Beijing(No.2462023BJRC025).Moreover,we would like to express our heartfelt appreciation to the Computational Geosciences group in the Department of Mathematics and Cybernetics at SINTEF Digital for developing and providing the free open-source MATLAB Reservoir Simulation Toolbox(MRST)used in this research.
文摘This study introduces a novel method integrating CO_(2)flooding with radial borehole fracturing for enhanced oil recovery and CO_(2)underground storage,a solution to the limited vertical stimulation reservoir volume in horizontal well fracturing.A numerical model is established to investigate the production rate,reservoir pressure field,and CO_(2)saturation distribution corresponding to changing time of CO_(2)flooding with radial borehole fracturing.A sensitivity analysis on the influence of CO_(2)injection location,layer spacing,pressure difference,borehole number,and hydraulic fractures on oil production and CO_(2)storage is conducted.The CO_(2)flooding process is divided into four stages.Reductions in layer spacing will significantly improve oil production rate and gas storage capacity.However,serious gas channeling can occur when the spacing is lower than 20 m.Increasing the pressure difference between the producer and injector,the borehole number,the hydraulic fracture height,and the fracture width can also increase the oil production rate and gas storage rate.Sensitivity analysis shows that layer spacing and fracture height greatly influence gas storage and oil production.Research outcomes are expected to provide a theoretical basis for the efficient development of shale oil reservoirs in the vertical direction.