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.展开更多
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.展开更多
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.展开更多
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.展开更多
Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based ...Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.展开更多
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.展开更多
Horizontal well-stimulation is the key to unconventional resource exploration and development.The development mode of the well plant helps increase the stimulated reservoir volume.Nevertheless,fracture interference be...Horizontal well-stimulation is the key to unconventional resource exploration and development.The development mode of the well plant helps increase the stimulated reservoir volume.Nevertheless,fracture interference between wells reduces the fracturing effect.Here,a 2D hydro-mechanical coupling model describing hydraulic fracture(HF)propagation is established with the extended finite element method,and the effects of several factors on HF propagation during multiple wells fracturing are analyzed.The results show that with an increase in elastic modulus,horizontal principal stress difference and injection fluid displacement,the total fracture area and the reservoir stimulation efficiency are both improved in all three fracturing technologies.After a comparison of the three technologies,the method of improved zipper fracturing is proposed,which avoids mutual interference between HFs,and the reservoir stimulation effect is improved significantly.The study provides guidance for optimizing the fracturing technology of multiple horizontal wells.展开更多
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.展开更多
A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the dr...A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.展开更多
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.展开更多
Shale gas is an important component of unconventional oil and gas resources.Studying the imbibition behavior is helpful to optimize flowback parameters and enhance gas recovery.Recent imbibition studies have focused o...Shale gas is an important component of unconventional oil and gas resources.Studying the imbibition behavior is helpful to optimize flowback parameters and enhance gas recovery.Recent imbibition studies have focused on shale matrix,and the pressure conditions discussed were mostly atmospheric.The initial imbibition behavior begins from propped fractures to matrix,but there are few studies working on explaining the imbibition behavior in propped fractures or the phenomenon of many shale wells exhibit higher productivity after a“soaking”period.Therefore,propped fracture samples were designed for imbibition and migration experiments.In order to accurately study the mechanism and main influencing factors of fracturing fluid imbibition and migration in propped and unpropped shale fractures under high temperature and high pressure,a series of experiments based on nuclear magnetic resonance(NMR)were carried out.Results showed that NMR T_(2) spectra of all samples exhibited a bimodal distribution.The final imbibition volume of fracturing fluid was positively related to pressure and fracture width.The imbibition effect of fracturing fluid was more evident in matrix pores under high pressure.In the migration during soaking stage,the fracturing fluid gradually migrated from large pores to small pores and gradually displaced the shale gas from the matrix,thus allowing the water blocking in propped fractures to self-unlock to some extent.Gas permeability decreased in the imbibition stage,while it recovered in the migration stage to some extent.展开更多
Based on the actual data collected from the tight sandstone development zone, correlation analysis using theSpearman method was conducted to determine the main factors influencing the gas production rate of tightsands...Based on the actual data collected from the tight sandstone development zone, correlation analysis using theSpearman method was conducted to determine the main factors influencing the gas production rate of tightsandstone fracturing. An integrated model combining geological engineering and numerical simulation of fracturepropagation and production was completed. Based on data analysis, the hydraulic fracture parameters wereoptimized to develop a differentiated fracturing treatment adjustment plan. The results indicate that the influenceof geological and engineering factors in the X1 and X2 development zones in the study area differs significantly.Therefore, it is challenging to adopt a uniform development strategy to achieve rapid production increase. Thedata analysis reveals that the variation in gas production rate is primarily affected by the reservoir thickness andpermeability parameters as geological factors. On the other hand, the amount of treatment fluid and proppantaddition significantly impact the gas production rate as engineering factors. Among these factors, the influence ofgeological factors is more pronounced in block X1. Therefore, the main focus should be on further optimizing thefracturing interval and adjusting the geological development well location. Given the existing well location, thereis limited potential for further optimizing fracture parameters to increase production. For block X2, the fracturingparameters should be optimized. Data screening was conducted to identify outliers in the entire dataset, and adata-driven fracturing parameter optimization method was employed to determine the basic adjustment directionfor reservoir stimulation in the target block. This approach provides insights into the influence of geological,stimulation, and completion parameters on gas production rate. Consequently, the subsequent fracturing parameteroptimization design can significantly reduce the modeling and simulation workload and guide field operations toimprove and optimize hydraulic fracturing efficiency.展开更多
For the analysis of the formation damage caused by the compound function of drilling fluid and fracturing fluid,the prediction method for dynamic invasion depth of drilling fluid is developed considering the fracture ...For the analysis of the formation damage caused by the compound function of drilling fluid and fracturing fluid,the prediction method for dynamic invasion depth of drilling fluid is developed considering the fracture extension due to shale minerals erosion by oil-based drilling fluid.With the evaluation for the damage of natural and hydraulic fractures caused by mechanical properties weakening of shale fracture surface,fracture closure and rock powder blocking,the formation damage pattern is proposed with consideration of the compound effect of drilling fluid and fracturing fluid.The formation damage mechanism during drilling and completion process in shale reservoir is revealed,and the protection measures are raised.The drilling fluid can deeply invade into the shale formation through natural and induced fractures,erode shale minerals and weaken the mechanical properties of shale during the drilling process.In the process of hydraulic fracturing,the compound effect of drilling fluid and fracturing fluid further weakens the mechanical properties of shale,results in fracture closure and rock powder shedding,and thus induces stress-sensitive damage and solid blocking damage of natural/hydraulic fractures.The damage can yield significant conductivity decrease of fractures,and restrict the high and stable production of shale oil and gas wells.The measures of anti-collapse and anti-blocking to accelerate the drilling of reservoir section,forming chemical membrane to prevent the weakening of the mechanical properties of shale fracture surface,strengthening the plugging of shale fracture and reducing the invasion range of drilling fluid,optimizing fracturing fluid system to protect fracture conductivity are put forward for reservoir protection.展开更多
Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operation...Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.展开更多
Ultra-low permeability reservoirs are characterized by small pore throats and poor physical properties, which areat the root of well-known problems related to injection and production. In this study, a gas injection f...Ultra-low permeability reservoirs are characterized by small pore throats and poor physical properties, which areat the root of well-known problems related to injection and production. In this study, a gas injection floodingapproach is analyzed in the framework of numerical simulations. In particular, the sequence and timing of fracturechanneling and the related impact on production are considered for horizontal wells with different fracturemorphologies. Useful data and information are provided about the regulation of gas channeling and possible strategiesto delay gas channeling and optimize the gas injection volume and fracture parameters. It is shown that inorder to mitigate gas channeling and ensure high production, fracture length on the sides can be controlled andlonger fractures can be created in the middle by which full gas flooding is obtained at the fracture location in themiddle of the horizontal well. A Differential Evolution (DE) algorithm is provided by which the gas injectionvolume and the fracture parameters of gas injection flooding can be optimized. It is shown that an improvedoil recovery factor as high as 6% can be obtained.展开更多
Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown a...Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown are influenced by hydraulic fractures,which can reflect the geometric features of hydraulic fracture.The shutdown pressure can be used to interpret the hydraulic fracture parameters in a real-time and cost-effective manner.In this paper,a mathematical model for shutdown pressure evolution is developed considering the effects of wellbore friction,perforation friction and fluid loss in fractures.An efficient numerical simulation method is established by using the method of characteristics.Based on this method,the impacts of fracture half-length,fracture height,opened cluster and perforation number,and filtration coefficient on the evolution of shutdown pressure are analyzed.The results indicate that a larger fracture half-length may hasten the decay of shutdown pressure,while a larger fracture height can slow down the decay of shutdown pressure.A smaller number of opened clusters and perforations can significantly increase the perforation friction and decrease the overall level of shutdown pressure.A larger filtration coefficient may accelerate the fluid filtration in the fracture and hasten the drop of the shutdown pressure.The simulation method of shutdown pressure,as well as the analysis results,has important implications for the interpretation of hydraulic fracture parameters.展开更多
A numerical model of hydraulic fracture propagation is introduced for a representative reservoir(Yuanba continental tight sandstone gas reservoir in Northeast Sichuan).Different parameters are considered,i.e.,the inte...A numerical model of hydraulic fracture propagation is introduced for a representative reservoir(Yuanba continental tight sandstone gas reservoir in Northeast Sichuan).Different parameters are considered,i.e.,the interlayer stress difference,the fracturing discharge rate and the fracturing fluid viscosity.The results show that these factors affect the gas and water production by influencing the fracture size.The interlayer stress difference can effectively control the fracture height.The greater the stress difference,the smaller the dimensionless reconstruction volume of the reservoir,while the flowback rate and gas production are lower.A large displacement fracturing construction increases the fracture-forming efficiency and expands the fracture size.The larger the displacement of fracturing construction,the larger the dimensionless reconstruction volume of the reservoir,and the higher the fracture-forming efficiency of fracturing fluid,the flowback rate,and the gas production.Low viscosity fracturing fluid is suitable for long fractures,while high viscosity fracturing fluid is suitable for wide fractures.With an increase in the fracturing fluid viscosity,the dimensionless reconstruction volume and flowback rate of the reservoir display a non-monotonic behavior,however,their changes are relatively small.展开更多
A deep understanding of the geometric impacts of fracture on fracturing fluid flowback efficiency is essential for unconventional oil development. Using nuclear magnetic resonance and 2.5-dimensional matrix-fracture v...A deep understanding of the geometric impacts of fracture on fracturing fluid flowback efficiency is essential for unconventional oil development. Using nuclear magnetic resonance and 2.5-dimensional matrix-fracture visualization microfluidic models, qualitative and quantitative descriptions of the influences of connectivity between primary fracture and secondary fracture on flowback were given from core scale to pore network scale. The flow patterns of oil-gel breaking fluid two-phase flow during flowback under different fracture connectivity were analyzed. We found some counterintuitive results that non-connected secondary fracture (NCSF, not connect with artificial primary fracture and embedded in the matrix) is detrimental to flowbackefficiency. The NCSF accelerates the formation of oil channeling during flowback, resulting in a large amount of fracturing fluid trapped in the matrix, which is not beneficial for flowback. Whereas the connected secondary fracture (CSF, connected with the artificial primary fracture) is conducive to flowback. The walls of CSF become part of primary fracture, which expands the drainage area with low resistance, and delays the formation of the oil flow channel. Thus, CSF increases the high-speed flowback stage duration, thereby enhancing the flowback efficiency. The fracturing fluid flowback efficiency investigated here follows the sequence of the connected secondary fracture model (72%) > the matrix model (66%) > the non-connected secondary fracture model (38%). Our results contribute to hydraulic fracturing design and the prediction of flowback efficiency.展开更多
This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation...This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.展开更多
We aim at the development of a general modelling workflow for design and optimization of the well flowback and startup operation on hydraulically fractured wells.Fracture flowback model developed earlier by the author...We aim at the development of a general modelling workflow for design and optimization of the well flowback and startup operation on hydraulically fractured wells.Fracture flowback model developed earlier by the authors is extended to take into account several new fluid mechanics factors accompanying flowback,namely,viscoplastic rheology of unbroken cross-linked gel and coupled“fracture-reservoir”numerical submodel for influx from rock formation.We also developed models and implemented new geomechanical factors,namely,(i)fracture closure in gaps between proppant pillars and in proppant-free cavity in the vicinity of the well taking into account formation creep;(ii)propagation of plastic deformations due to tensile rock failure from the fracture face into the fluid-saturated reservoir.We carried out parametric calculations to study the dynamics of fracture conductivity during flowback and its effect on well production for the set of parameters typical of oil wells in Achimov formation of Western Siberia,Russia.The first set of calculations is carried out using the flowback model in the reservoir linear flow regime.It is obtained that the typical length of hydraulic fracture zone,in which tensile rock failure at the fracture walls occurs,is insignificant.In the range of rock permeability in between 0.01 mD and 1 D,we studied the effect of non-dimensional governing parameters as well as bottomhole pressure drop dynamics on oil production.We obtained a map of pressure drop regimes(fast,moderate or slow)leading to maximum cumulative oil production.The second set of parametric calculations is carried out using integrated well production modelling workflow,in which the flowback model acts as a missing link in between hydraulic fracturing and reservoir commercial simulators.We evaluated quantitatively effects of initial fracture aperture,proppant diameter,yield stress of fracturing fluid,pressure drop rate and proppant material type(ceramic and sand)on long-term well production beyond formation linear regime.The third set of parametric calculations is carried out using the flowback model history-matched to field data related to production of four multistage hydraulically fractured oil wells in Achimov formation of Western Siberia,Russia.On the basis of the matched model we evaluated geomechanics effects on fracture conductivity degradation.We also performed sensitivity analysis in the framework of the history-matched model to study the impact of geomechanics and fluid rheology parameters on flowback efficiency.展开更多
文摘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.
基金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 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.
基金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.
基金the National Natural Science Foundation of China(No.51974042)National Key Research and Development Program of China(No.2023YFC3009005).
文摘Ground hydraulic fracturing plays a crucial role in controlling the far-field hard roof,making it imperative to identify the most suitable target stratum for effective control.Physical experiments are conducted based on engineering properties to simulate the gradual collapse of the roof during longwall top coal caving(LTCC).A numerical model is established using the material point method(MPM)and the strain-softening damage constitutive model according to the structure of the physical model.Numerical simulations are conducted to analyze the LTCC process under different hard roofs for ground hydraulic fracturing.The results show that ground hydraulic fracturing releases the energy and stress of the target stratum,resulting in a substantial lag in the fracturing of the overburden before collapse occurs in the hydraulic fracturing stratum.Ground hydraulic fracturing of a low hard roof reduces the lag effect of hydraulic fractures,dissipates the energy consumed by the fracture of the hard roof,and reduces the abutment stress.Therefore,it is advisable to prioritize the selection of the lower hard roof as the target stratum.
基金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.
基金funded by Shaanxi Natural Science Basic Research Program Project Study on Liquid Propellant High Energy Gas Fracturing Mechanism in Radial Well Based on Phase Field Method(No.2019JQ-824)NSFC Projects Evolution Mechanism and Effectiveness Evaluation of Fracture Network Produced by Volume Fracturing with Tighter Clusters in Continental Shale Oil Reservoir(No.52274040)+1 种基金Study on Thermal Secondary Pore Evolution and Salt Precipitation Regulation Mechanism in Fire Flooding Reservoirs Based on Multi-field Coupling of Thermal-Flow-Solid-Chemical(No.52274039)Xi’an Shiyou University Youth Scientific Research and Innovation Team Operation Funds in 2018 Flow Mechanism of Complex Reservoirs and High Efficiency Development and Oil Production Technology(No.115080020).
文摘Horizontal well-stimulation is the key to unconventional resource exploration and development.The development mode of the well plant helps increase the stimulated reservoir volume.Nevertheless,fracture interference between wells reduces the fracturing effect.Here,a 2D hydro-mechanical coupling model describing hydraulic fracture(HF)propagation is established with the extended finite element method,and the effects of several factors on HF propagation during multiple wells fracturing are analyzed.The results show that with an increase in elastic modulus,horizontal principal stress difference and injection fluid displacement,the total fracture area and the reservoir stimulation efficiency are both improved in all three fracturing technologies.After a comparison of the three technologies,the method of improved zipper fracturing is proposed,which avoids mutual interference between HFs,and the reservoir stimulation effect is improved significantly.The study provides guidance for optimizing the fracturing technology of multiple horizontal wells.
基金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.
基金Project supported by the National Natural Science Foundation of China(No.42202314)。
文摘A novel phase-field model for the propagation of mixed-mode hydraulic fractures,characterized by the formation of mixed-mode fractures due to the interactions between fluids and solids,is proposed.In this model,the driving force for the phase field consists of both tensile and shear components,with the fluid contribution primarily manifesting in the tension driving force.The displacement and pressure are solved simultaneously by an implicit method.The numerical solution's iterative format is established by the finite element discretization and Newton-Raphson(NR)iterative methods.The correctness of the model is verified through the uniaxial compression physical experiments on fluid-pressurized rocks,and the limitations of the hydraulic fracture expansion phase-field model,which only considers mode I fractures,are revealed.In addition,the influence of matrix mode II fracture toughness value,natural fracture mode II toughness value,and fracturing fluid injection rate on the hydraulic fracture propagation in porous media with natural fractures is studied.
基金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.
基金The authors gratefully acknowledge the support of the National Natural Science Foundation of China(Grant Nos.52174036,51774243,51904257,51874251)the Sichuan Province Science and Technology Program(Grant Nos.2021YJ0345,2022JDJQ0009,2022NSFSC0186).
文摘Shale gas is an important component of unconventional oil and gas resources.Studying the imbibition behavior is helpful to optimize flowback parameters and enhance gas recovery.Recent imbibition studies have focused on shale matrix,and the pressure conditions discussed were mostly atmospheric.The initial imbibition behavior begins from propped fractures to matrix,but there are few studies working on explaining the imbibition behavior in propped fractures or the phenomenon of many shale wells exhibit higher productivity after a“soaking”period.Therefore,propped fracture samples were designed for imbibition and migration experiments.In order to accurately study the mechanism and main influencing factors of fracturing fluid imbibition and migration in propped and unpropped shale fractures under high temperature and high pressure,a series of experiments based on nuclear magnetic resonance(NMR)were carried out.Results showed that NMR T_(2) spectra of all samples exhibited a bimodal distribution.The final imbibition volume of fracturing fluid was positively related to pressure and fracture width.The imbibition effect of fracturing fluid was more evident in matrix pores under high pressure.In the migration during soaking stage,the fracturing fluid gradually migrated from large pores to small pores and gradually displaced the shale gas from the matrix,thus allowing the water blocking in propped fractures to self-unlock to some extent.Gas permeability decreased in the imbibition stage,while it recovered in the migration stage to some extent.
基金Research and Application of Key Technologies for Tight Gas Production Improvement and Rehabilitation of Linxing Shenfu(YXKY-ZL-01-2021)。
文摘Based on the actual data collected from the tight sandstone development zone, correlation analysis using theSpearman method was conducted to determine the main factors influencing the gas production rate of tightsandstone fracturing. An integrated model combining geological engineering and numerical simulation of fracturepropagation and production was completed. Based on data analysis, the hydraulic fracture parameters wereoptimized to develop a differentiated fracturing treatment adjustment plan. The results indicate that the influenceof geological and engineering factors in the X1 and X2 development zones in the study area differs significantly.Therefore, it is challenging to adopt a uniform development strategy to achieve rapid production increase. Thedata analysis reveals that the variation in gas production rate is primarily affected by the reservoir thickness andpermeability parameters as geological factors. On the other hand, the amount of treatment fluid and proppantaddition significantly impact the gas production rate as engineering factors. Among these factors, the influence ofgeological factors is more pronounced in block X1. Therefore, the main focus should be on further optimizing thefracturing interval and adjusting the geological development well location. Given the existing well location, thereis limited potential for further optimizing fracture parameters to increase production. For block X2, the fracturingparameters should be optimized. Data screening was conducted to identify outliers in the entire dataset, and adata-driven fracturing parameter optimization method was employed to determine the basic adjustment directionfor reservoir stimulation in the target block. This approach provides insights into the influence of geological,stimulation, and completion parameters on gas production rate. Consequently, the subsequent fracturing parameteroptimization design can significantly reduce the modeling and simulation workload and guide field operations toimprove and optimize hydraulic fracturing efficiency.
基金Supported by the Key Fund Project of the National Natural Science Foundation of China and Joint Fund of Petrochemical Industry(Class A)(U1762212)National Natural Science Foundation of China(52274009)"14th Five-Year"Forward-looking and Fundamental Major Science and Technology Project of CNPC(2021DJ4402)。
文摘For the analysis of the formation damage caused by the compound function of drilling fluid and fracturing fluid,the prediction method for dynamic invasion depth of drilling fluid is developed considering the fracture extension due to shale minerals erosion by oil-based drilling fluid.With the evaluation for the damage of natural and hydraulic fractures caused by mechanical properties weakening of shale fracture surface,fracture closure and rock powder blocking,the formation damage pattern is proposed with consideration of the compound effect of drilling fluid and fracturing fluid.The formation damage mechanism during drilling and completion process in shale reservoir is revealed,and the protection measures are raised.The drilling fluid can deeply invade into the shale formation through natural and induced fractures,erode shale minerals and weaken the mechanical properties of shale during the drilling process.In the process of hydraulic fracturing,the compound effect of drilling fluid and fracturing fluid further weakens the mechanical properties of shale,results in fracture closure and rock powder shedding,and thus induces stress-sensitive damage and solid blocking damage of natural/hydraulic fractures.The damage can yield significant conductivity decrease of fractures,and restrict the high and stable production of shale oil and gas wells.The measures of anti-collapse and anti-blocking to accelerate the drilling of reservoir section,forming chemical membrane to prevent the weakening of the mechanical properties of shale fracture surface,strengthening the plugging of shale fracture and reducing the invasion range of drilling fluid,optimizing fracturing fluid system to protect fracture conductivity are put forward for reservoir protection.
基金funded by the joint fund of the National Key Research and Development Program of China(No.2021YFC2902101)National Natural Science Foundation of China(Grant No.52374084)+1 种基金Open Foundation of National Energy shale gas R&D(experiment) center(2022-KFKT-12)the 111 Project(B17009)。
文摘Evaluating the physical mechanisms that link hydraulic fracturing(HF) operations to induced earthquakes and the anticipated form of the resulting events is significant in informing subsurface fluid injection operations. Current understanding supports the overriding role of the effective stress magnitude in triggering earthquakes, while the impact of change rate of effective stress has not been systematically addressed. In this work, a modified critical stiffness was brought up to investigate the likelihood, impact,and mitigation of induced seismicity during and after hydraulic fracturing by developing a poroelastic model based on rate-and-state fraction law and linear stability analysis. In the new criterion, the change rate of effective stress was considered a key variable to explore the evolution of this criterion and hence the likelihood of instability slip of fault. A coupled fluid flow-deformation model was used to represent the entire hydraulic fracturing process in COMSOL Multiphysics. The possibility of triggering an earthquake throughout the entire hydraulic fracturing process, from fracturing to cessation, was investigated considering different fault locations, orientations, and positions along the fault. The competition between the effects of the magnitude and change rate of effective stress was notable at each fracturing stage. The effective stress magnitude is a significant controlling factor during fracturing events, with the change rate dominating when fracturing is suddenly started or stopped. Instability dominates when the magnitude of the effective stress increases(constant injection at each fracturing stage) and the change rate of effective stress decreases(the injection process is suddenly stopped). Fracturing with a high injection rate, a fault adjacent to the hydraulic fracturing location and the position of the junction between the reservoir and fault are important to reduce the Coulomb failure stress(CFS) and enhance the critical stiffness as the significant disturbance of stresses at these positions in the coupled process. Therefore,notable attention should be given to the injection rate during fracturing, fault position, and position along faults as important considerations to help reduce the potential for induced seismicity. Our model was verified and confirmed using the case of the Longmaxi Formation in the Sichuan Basin, China, in which the reported microseismic data were correlated with high critical stiffness values. This work supplies new thoughts of the seismic risk associated with HF engineering.
基金supported by the Forward Looking Basic Major Scientific and Technological Projects of CNPC (Grant No.2021DJ2202).
文摘Ultra-low permeability reservoirs are characterized by small pore throats and poor physical properties, which areat the root of well-known problems related to injection and production. In this study, a gas injection floodingapproach is analyzed in the framework of numerical simulations. In particular, the sequence and timing of fracturechanneling and the related impact on production are considered for horizontal wells with different fracturemorphologies. Useful data and information are provided about the regulation of gas channeling and possible strategiesto delay gas channeling and optimize the gas injection volume and fracture parameters. It is shown that inorder to mitigate gas channeling and ensure high production, fracture length on the sides can be controlled andlonger fractures can be created in the middle by which full gas flooding is obtained at the fracture location in themiddle of the horizontal well. A Differential Evolution (DE) algorithm is provided by which the gas injectionvolume and the fracture parameters of gas injection flooding can be optimized. It is shown that an improvedoil recovery factor as high as 6% can be obtained.
基金The work is supported by the Sub-Project of“Research on Key Technologies and Equipment of Reservoir Stimulation”of China National Petroleum Corporation Post–14th Five-Year Plan Forward-Looking Major Science and Technology Project“Research on New Technology of Monitoring and Diagnosis of Horizontal Well Hydraulic Fracture Network Distribution Pattern”(2021DJ4502).
文摘Multistage multi-cluster hydraulic fracturing has enabled the economic exploitation of shale reservoirs,but the interpretation of hydraulic fracture parameters is challenging.The pressure signals after pump shutdown are influenced by hydraulic fractures,which can reflect the geometric features of hydraulic fracture.The shutdown pressure can be used to interpret the hydraulic fracture parameters in a real-time and cost-effective manner.In this paper,a mathematical model for shutdown pressure evolution is developed considering the effects of wellbore friction,perforation friction and fluid loss in fractures.An efficient numerical simulation method is established by using the method of characteristics.Based on this method,the impacts of fracture half-length,fracture height,opened cluster and perforation number,and filtration coefficient on the evolution of shutdown pressure are analyzed.The results indicate that a larger fracture half-length may hasten the decay of shutdown pressure,while a larger fracture height can slow down the decay of shutdown pressure.A smaller number of opened clusters and perforations can significantly increase the perforation friction and decrease the overall level of shutdown pressure.A larger filtration coefficient may accelerate the fluid filtration in the fracture and hasten the drop of the shutdown pressure.The simulation method of shutdown pressure,as well as the analysis results,has important implications for the interpretation of hydraulic fracture parameters.
文摘A numerical model of hydraulic fracture propagation is introduced for a representative reservoir(Yuanba continental tight sandstone gas reservoir in Northeast Sichuan).Different parameters are considered,i.e.,the interlayer stress difference,the fracturing discharge rate and the fracturing fluid viscosity.The results show that these factors affect the gas and water production by influencing the fracture size.The interlayer stress difference can effectively control the fracture height.The greater the stress difference,the smaller the dimensionless reconstruction volume of the reservoir,while the flowback rate and gas production are lower.A large displacement fracturing construction increases the fracture-forming efficiency and expands the fracture size.The larger the displacement of fracturing construction,the larger the dimensionless reconstruction volume of the reservoir,and the higher the fracture-forming efficiency of fracturing fluid,the flowback rate,and the gas production.Low viscosity fracturing fluid is suitable for long fractures,while high viscosity fracturing fluid is suitable for wide fractures.With an increase in the fracturing fluid viscosity,the dimensionless reconstruction volume and flowback rate of the reservoir display a non-monotonic behavior,however,their changes are relatively small.
基金supported by the National Key Research and Development Program of China(Grant No.2019YFA0708700).
文摘A deep understanding of the geometric impacts of fracture on fracturing fluid flowback efficiency is essential for unconventional oil development. Using nuclear magnetic resonance and 2.5-dimensional matrix-fracture visualization microfluidic models, qualitative and quantitative descriptions of the influences of connectivity between primary fracture and secondary fracture on flowback were given from core scale to pore network scale. The flow patterns of oil-gel breaking fluid two-phase flow during flowback under different fracture connectivity were analyzed. We found some counterintuitive results that non-connected secondary fracture (NCSF, not connect with artificial primary fracture and embedded in the matrix) is detrimental to flowbackefficiency. The NCSF accelerates the formation of oil channeling during flowback, resulting in a large amount of fracturing fluid trapped in the matrix, which is not beneficial for flowback. Whereas the connected secondary fracture (CSF, connected with the artificial primary fracture) is conducive to flowback. The walls of CSF become part of primary fracture, which expands the drainage area with low resistance, and delays the formation of the oil flow channel. Thus, CSF increases the high-speed flowback stage duration, thereby enhancing the flowback efficiency. The fracturing fluid flowback efficiency investigated here follows the sequence of the connected secondary fracture model (72%) > the matrix model (66%) > the non-connected secondary fracture model (38%). Our results contribute to hydraulic fracturing design and the prediction of flowback efficiency.
基金supported by the National Natural Science Foundation of China(Nos.52274038,5203401042174143)+1 种基金the Taishan Scholars Project(No.tsqnz20221140)the Open Fund of State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University)of China(No.PLN2020-5)。
文摘This paper presents an integrated study from fracture propagation modeling to gas flow modeling and a correlation analysis to explore the key controlling factors of intensive volume fracturing.The fracture propagation model takes into account the interaction between hydraulic fracture and natural fracture by means of the displacement discontinuity method(DDM)and the Picard iterative method.The shale gas flow considers multiple transport mechanisms,and the flow in the fracture network is handled by the embedded discrete fracture model(EDFM).A series of numerical simulations are conducted to analyze the effects of the cluster number,stage spacing,stress difference coefficient,and natural fracture distribution on the stimulated fracture area,fractal dimension,and cumulative gas production,and their correlation coefficients are obtained.The results show that the most influential factors to the stimulated fracture area are the stress difference ratio,stage spacing,and natural fracture density,while those to the cumulative gas production are the stress difference ratio,natural fracture density,and cluster number.This indicates that the stress condition dominates the gas production,and employing intensive volume fracturing(by properly increasing the cluster number)is beneficial for improving the final cumulative gas production.
文摘We aim at the development of a general modelling workflow for design and optimization of the well flowback and startup operation on hydraulically fractured wells.Fracture flowback model developed earlier by the authors is extended to take into account several new fluid mechanics factors accompanying flowback,namely,viscoplastic rheology of unbroken cross-linked gel and coupled“fracture-reservoir”numerical submodel for influx from rock formation.We also developed models and implemented new geomechanical factors,namely,(i)fracture closure in gaps between proppant pillars and in proppant-free cavity in the vicinity of the well taking into account formation creep;(ii)propagation of plastic deformations due to tensile rock failure from the fracture face into the fluid-saturated reservoir.We carried out parametric calculations to study the dynamics of fracture conductivity during flowback and its effect on well production for the set of parameters typical of oil wells in Achimov formation of Western Siberia,Russia.The first set of calculations is carried out using the flowback model in the reservoir linear flow regime.It is obtained that the typical length of hydraulic fracture zone,in which tensile rock failure at the fracture walls occurs,is insignificant.In the range of rock permeability in between 0.01 mD and 1 D,we studied the effect of non-dimensional governing parameters as well as bottomhole pressure drop dynamics on oil production.We obtained a map of pressure drop regimes(fast,moderate or slow)leading to maximum cumulative oil production.The second set of parametric calculations is carried out using integrated well production modelling workflow,in which the flowback model acts as a missing link in between hydraulic fracturing and reservoir commercial simulators.We evaluated quantitatively effects of initial fracture aperture,proppant diameter,yield stress of fracturing fluid,pressure drop rate and proppant material type(ceramic and sand)on long-term well production beyond formation linear regime.The third set of parametric calculations is carried out using the flowback model history-matched to field data related to production of four multistage hydraulically fractured oil wells in Achimov formation of Western Siberia,Russia.On the basis of the matched model we evaluated geomechanics effects on fracture conductivity degradation.We also performed sensitivity analysis in the framework of the history-matched model to study the impact of geomechanics and fluid rheology parameters on flowback efficiency.