A review of reservoir damage during hydraulic fracturing of deep and ultra-deep reservoirs

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.

Evaluating the stability and volumetric flowback rate of proppant packs in hydraulic fractures using the lattice Boltzmann-discrete element coupling method

The stability and mobility of proppant packs in hydraulic fractures during hydrocarbon production are numerically investigated by the lattice Boltzmann-discrete element coupling method(LB-DEM).This study starts with a preliminary proppant settling test,from which a solid volume fraction of 0.575 is calibrated for the proppant pack in the fracture.In the established workflow to investigate proppant flowback,a displacement is applied to the fracture surfaces to compact the generated proppant pack as well as further mimicking proppant embedment under closure stress.When a pressure gradient is applied to drive the fluid-particle flow,a critical aperture-to-diameter ratio of 4 is observed,above which the proppant pack would collapse.The results also show that the volumetric proppant flowback rate increases quadratically with the fracture aperture,while a linear variation between the particle flux and the pressure gradient is exhibited for a fixed fracture aperture.The research outcome contributes towards an improved understanding of proppant flowback in hydraulic fractures,which also supports an optimised proppant size selection for hydraulic fracturing operations.

Extreme massive hydraulic fracturing in deep coalbed methane horizontal wells:A case study of the Linxing Block,eastern Ordos Basin,NW China

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.

Hydraulic properties and drought response of a tropical bamboo (Cephalostachyum pergracile)

Bamboo plants are an essential component of tropical ecosystems,yet their vulnerability to climate extremes,such as drought,is poorly understood due to limited knowledge of their hydraulic properties.Cephalostachyum pergracile,a commonly used tropical bamboo species,exhibited a substantially higher mortality rate than other co-occurring bamboos during a severe drought event in 2019,but the underlying mechanisms remain unclear.This study investigated the leaf and stem hydraulic traits related to drought responses,including leaf-stem embolism resistance(P50leaf;P50stem) estimated using optical and X-ray microtomography methods,leaf pressure-volume and water-releasing curves.Additionally,we investigated the seasonal water potentials,native embolism level(PLC) and xylem water source using stable isotope.We found that C.pergracile exhibited strong resistance to embolism,showing low P50leaf,P50stem,and turgor loss point,despite its rapid leaf water loss.Interestingly,its leaves displayed greater resistance to embolism than its stem,suggesting a lack of effective hydraulic vulnerability segmentation(HVS) to protect the stem from excessive xylem tension.During the dry season,approximately 49% of the water was absorbed from the upper 20-cm-deep soil layer.Consequently,significant diurnal variation in leaf water potentials and an increase in midday PLC from 5.87±2.33% in the wet season to 12.87±4.09%in the dry season were observed.In summary,this study demonstrated that the rapid leaf water loss,high reliance on surface water,and a lack of effective HVS in C.pergracile accelerated water depletion and increased xylem embolism even in the typical dry season,which may explain its high mortality rate during extreme drought events in 2019.

Stress tensor determination by modified hydraulic tests on pre-existing fractures:Method and stress constraints

The hydraulic testing of pre-existing fractures(HTPF)is one of the most promising in situ stress measurement methods,particularly for three-dimensional stress tensor determination.However,the stress tensor determination based on the HTPF method requires at least six tests or a minimum of 14-15 tests(under different conditions)for reliable results.In this study,we modified the HTPF method by considering the shear stress on each pre-existing fracture,which increased the number of equations for the stress tensor determination and decreased the number of tests required.Different shear stresses were attributed to different fractures by random sampling;therefore,the stress tensors were obtained by searching for the optimal solution using the least squares criterion based on the Monte Carlo method.Thereafter,we constrained the stress tensor based on the tensile strength criterion,compressive strength criterion,and vertical stress constraints.The inverted stress tensors were presented and analyzed based on the tensorial nature of the stress using the Euclidean mean stress tensor.Two stress-measurement campaigns in Weifang(Shandong Province,China)and Mercantour road tunnel(France)were implemented to highlight the validity and efficiency of the modified HTPF(M-HTPF)method.The results showed that the M-HTPF method can be applied for stress tensor inversion using only three to four tests on pre-existing fractures,neglecting the stress gradient.The inversion results were confined to relatively small distribution dispersions and were significantly reliable and stable due to the shear stresses on the fractures and the stress constraints employed.The M-HTPF method is highly feasible and efficient for complete stress tensor determination in a single borehole.

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

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

Influence of Blade Number on the Performance of Hydraulic Turbines in the Transition Stage

To analyze the effect of blade number on the performance of hydraulic turbines during the transient stage in which theflow rate is not constant,six hydraulic turbines with different blade numbers are considered.The instantaneous hydraulic performance of the turbine and the pressure pulsation acting on the impeller are investigated numerically by using the ANSYS CFX software.The ensuing results are compared with the outcomes of experimental tests.It is shown that thefluctuation range of the pressure coefficient increases with time,but the corresponding range for the transient hydraulic efficiency decreases gradually when theflow velocity transits to larger values.During the transition to smallflow velocity,thefluctuation range of the pressure coefficient gradually decreases as time passes,but the correspondingfluctuation range of its transient hydraulic efficiency gradually becomes larger.Thefluctuation range in the Z9 case is small during the transition.The main frequency of transient hydraulic efficiency pulsation is equal to the blade frequency.At the main frequency,Z7 has the largest amplitude of the hydraulic efficiency pulsation,Z10 has the smallest amplitude,and the difference between Z7 and Z9 is limited.As the number of blades grows,the pressure pulsation during the transition process gradually decreases,but the pressure pulsation of Z10 at the volute tongue is larger.In the steady state,Z9 has the highest efficiency and in the transient stage,the pressure coefficientfluctuation range is small.Accordingly,for the hydraulic turbine Z9,the performance is optimal.

Hydraulic fracturing behaviors of shale under coupled stress and temperature conditions simulating different burial depths

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.

Pumping-induced Well Hydraulics and Groundwater Budget in a Leaky Aquifer System with Vertical Heterogeneity in Aquitard Hydraulic Properties

In groundwater hydrology,aquitard heterogeneity is often less considered compared to aquifers,despite its significant impact on groundwater hydraulics and groundwater resources evaluation.A semi-analytical solution is derived for pumping-induced well hydraulics and groundwater budget with consideration of vertical heterogeneity in aquitard hydraulic conductivity(K)and specific storage(S_(s)).The proposed new solution is innovative in its partitioning of the aquitard into multiple homogeneous sub-layers to enable consideration of various forms of vertically heterogeneous K or S_(s).Two scenarios of analytical investigations are explored:one is the presence of aquitard interlayers with distinct K or S_(s) values,a common field-scale occurrence;another is an exponentially depth-decaying aquitard S_(s),a regional-scale phenomenon supported by statistical analysis.Analytical investigations reveal that a low-K interlayer can significantly increase aquifer drawdown and enhance aquifer/aquitard depletion;a high-S_(s) interlayer can noticeably reduce aquifer drawdown and increase aquitard depletion.Locations of low-K or high-S_(s) interlayers also significantly impact well hydraulics and groundwater budget.In the context of an exponentially depth-decaying aquitard S_(s),a larger decay exponent can enhance aquifer drawdown.When using current models with a vertically homogeneous aquitard,half the sum of the geometric and harmonic means of exponentially depth-decaying aquitard S_(s) should be used to calculate aquitard depletion and unconfined aquifer leakage.

Experimental Investigation on Vertical Hydraulic Transport of Ores in Deepsea Mining

Deepsea mining has been proposed since the 1960s to alleviate the lack of resources on land.Vertical hydraulic transport of collected ores from the seabed to the sea surface is considered the most promising method for industrial applications.In the present study,an indoor model test of the vertical hydraulic transport of particles was conducted.A noncontact optical method has been proposed to measure the local characteristics of the particles inside a vertical pipe,including the local concentration and particle velocity.The hydraulic gradient of ore transport was evaluated with various particle size distributions,particle densities,feeding concentrations and mixture flow velocities.During transport,the local concentration is larger than the feeding concentration,whereas the particle velocity is less than the mixture velocity.The qualitative effects of the local concentration and local fluid velocity on the particle velocity and slip velocity were investigated.The local fluid velocity contributes significantly to particle velocity and slip velocity,whereas the effect of the local concentration is marginal.A higher feeding concentration and mixture flow velocity result in an increased hydraulic gradient.The effect of the particle size gradation is slight,whereas the particle density plays a crucial role in the transport.

Precision Motion Control of Hydraulic Actuator Using Adaptive Back-Stepping Sliding Mode Controller

Hydraulic actuators are highly nonlinear when they are subjected to different types of model uncertainties and dynamic disturbances.These unfavorable factors adversely affect the control performance of the hydraulic actuator.Although various control methods have been employed to improve the tracking precision of the dynamic system,optimizing and adjusting control gain to mitigate the hydraulic actuator model uncertainties remains elusive.This study presents an adaptive back-stepping sliding mode controller(ABSMC)to enhance the trajectory tracking precision,where the virtual control law is constructed to replace the position error.The adaptive control theory is introduced in back-stepping controller design to compensate for the model uncertainties and time-varying disturbances.Based on Lyapunov theory,the finite-time convergence of the position tracking errors is proved.Furthermore,the effectiveness of the developed control scheme is conducted via extensive comparative experiments.

Unraveling the hydraulic properties of loess for landslide prediction:A study on variations in loess landslides in Lanzhou,Dingxi,and Tianshui,China

Loess has distinctive characteristics,leading to frequent landslide disasters and posing serious threats to the lives and properties of local re sidents.The involvement of water repre sents a critical factor in inducing loess landslides.This study focuses on three neighboring cities sequentially situated on the Loess Plateau along the direction of aeolian deposition of loess,namely Lanzhou,Dingxi,and Tianshui,which are densely populated and prone to landslide disasters.The variations in hydraulic properties,including water retention capacity and permeability,are investigated through Soil Water Characteristic Curve(SWCC)test and hydraulic conductivity test.The experimental findings revealed that Tianshui loess exhibited the highest water retention capacity,followed by Dingxi loess,while Lanzhou loess demonstrated the lowest water retention capacity.Contrastingly,the results for the saturated permeability coefficient were found to be the opposite:Tianshui loess showed the lowest permeability,whereas Lanzhou loess displayed the highest permeability.These results are supported and analyzed by scanning electron microscopy(SEM)observation.In addition,the water retention capacity is mathematically expressed using the van Genuchten model and extended to predict unsaturated hydraulic properties of loess.The experimental results exhibit a strong accordance with one another and align with the regional distribution patterns of disasters.

Jet Characteristics and Optimization of a Cavitation Nozzle for Hydraulic Fracturing Applications

Hydraulic jetting is a form of fracturing that involves using a high-pressure jet of water to create fractures in the reservoir rock with a nozzle serving as the central component of the hydraulic sandblasting perforation tool.In this study,the flow behavior of the nozzle is simulated numerically in the framework of a SST k-ωturbulence model.The results show that the nozzle structure can significantly influence the jet performance and related cavitation effect.Through orthogonal experiments,the nozzle geometric parameters are optimized,and the following configuration is found accordingly:contraction angle 20°,contraction segment length 6 mm,cylindrical segment diameter 6 mm,cylindrical segment length 12 mm,spread segment length 10 mm,and spread angle 55°.

A phase-field model for simulating the propagation behavior of mixed-mode cracks during the hydraulic fracturing process in fractured reservoirs

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.

Spatiotemporal variations of sand hydraulic conductivity by microbial application methods

The spatiotemporal distributions of microbes in soil by different methods could affect the efficacy of the microbes to reduce the soil hydraulic conductivity.In this study,the specimens of bio-mediated sands were prepared using three different methods,i.e.injecting,mixing,and pouring a given microbial so-lution onto compacted sand specimens.The hydraulic conductivity was measured by constant-head tests,while any soil microstructural changes due to addition of the microbes were observed by scan-ning electron microscope(SEM)and mercury intrusion porosimetry(MIP)tests.The amount of dextran concentration produced by microbes in each type of specimen was quantified by a refractometer.Results show that dextran production increased exponentially after 5-7 d of microbial settling with the supply of culture medium.The injection and mixing methods resulted in a similar amount and uniform dis-tribution of dextran in the specimens.The pouring method,however,produced a nonuniform distri-bution,with a higher concentration near the specimen surface.As the supply of culture medium discontinued,the dextran content near the surface produced by the pouring method decreased dramatically due to high competition for nutrients with foreign colonies.Average dextran concentration was negatively and correlated with hydraulic conductivity of bio-mediated soils exponentially,due to the clogging of large soil pores by dextran.The hydraulic conductivity of the injection and mixing cases did not change significantly when the supply of culture medium was absent.

Influence of High-Density Bedding Plane Characteristics on Hydraulic Fracture Propagation in Shale Oil Reservoir

The existence of high-density bedding planes is a typical characteristic of shale oil reservoirs.Understanding the behavior of hydraulic fracturing in high-density laminated rocks is significant for promoting shale oil production.In this study,a hydraulic fracturing model considering tensile failure and frictional slip of the bedding planes is established within the framework of the unified pipe-interface element method(UP-IEM).The model developed for simulating the interaction between the hydraulic fracture and the bedding plane is validated by comparison with experimental results.The hydraulic fracturing patterns in sealed and unsealed bedding planes are compared.Additionally,the effects of differential stress,bedding plane permeability,spacing,and the friction coefficient of the bedding plane are investigated.The results showed that a single main fracture crossing the bedding planes is more likely to form in sealed bedding planes under high differential stress.The decrease in bedding plane permeability and the increase in the friction coefficient also promote the fracture propagating perpendicular to the bedding planes.Shale with high-density bedding planes has a poorer fracturing effect than that with low-density bedding planes,as the hydraulic fracture is prone to initiate and propagate along the bedding planes.Moreover,higher injection pressure is needed to maintain fracture propagation along the bedding.An increase in bedding density will lead to a smaller fracturing area.Fracturing fluid seepage into the bedding planes slows shale fracturing.It is recommended that increasing the injection flow rate,selecting alternative fracturing fluids,and employing multi-well/multi-cluster fracturing may be efficient methods to improve energy production in shale oil reservoirs.

Effect mechanism of seepage force on the hydraulic fracture propagation

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

Estimation of the anisotropy of hydraulic conductivity through 3D fracture networks using the directional geological entropy

With an extension of the geological entropy concept in porous media,the approach called directional entrogram is applied to link hydraulic behavior to the anisotropy of the 3D fracture networks.A metric called directional entropic scale is used to measure the anisotropy of spatial order in different directions.Compared with the traditional connectivity indexes based on the statistics of fracture geometry,the directional entropic scale is capable to quantify the anisotropy of connectivity and hydraulic conductivity in heterogeneous 3D fracture networks.According to the numerical analysis of directional entrogram and fluid flow in a number of the 3D fracture networks,the hydraulic conductivities and entropic scales in different directions both increase with spatial order(i.e.,trace length decreasing and spacing increasing)and are independent of the dip angle.As a result,the nonlinear correlation between the hydraulic conductivities and entropic scales from different directions can be unified as quadratic polynomial function,which can shed light on the anisotropic effect of spatial order and global entropy on the heterogeneous hydraulic behaviors.

Numerical Analysis of Perforation during Hydraulic Fracture Initiation Based on Continuous-Discontinuous Element Method

Perforation is a pivotal technique employed to establish main flow channels within the reservoir formation at the outset of hydraulic fracturing operations.Optimizing perforation designs is critical for augmenting the efficacy of hydraulic fracturing and boosting oil or gas production.In this study,we employ a hybrid finite-discrete element method,known as the continuous–discontinuous element method(CDEM),to simulate the initiation of post-perforation hydraulic fractures and to derive enhanced design parameters.The model incorporates the four most prevalent perforation geometries,as delineated in an engineering technical report.Real-world perforations deviate from the ideal cylindrical shape,exhibiting variable cross-sectional profiles that typically manifest as an initial constriction followed by an expansion,a feature consistent across all four perforation types.Our simulations take into account variations in perforation hole geometries,cross-sectional diameters,and perforation lengths.The findings show that perforations generated by the 39g DP3 HMX perforating bullet yield the lowest breakdown pressure,which inversely correlates with increases in sectional diameter and perforation length.Moreover,this study reveals the relationship between breakdown pressure and fracture degree,providing valuable insights for engineers and designers to refine perforation strategies.

Using fracture-based continuum modeling of coupled geomechanical-hydrological processes for numerical simulation of hydraulic fracturing

This paper describes numerical simulation of hydraulic fracturing using fracture-based continuum modeling(FBCM)of coupled geomechanical-hydrological processes to evaluate a technique for high-density fracturing and fracture caging.The simulations are innovative because of modeling discrete fractures explicitly in continuum analysis.A key advantage of FBCM is that fracture initiation and propagation are modeled explicitly without changing the domain grid(i.e.no re-meshing).Further,multiple realizations of a preexisting fracture distribution can be analyzed using the same domain grid.The simulated hydraulic fracturing technique consists of pressurizing multiple wells simultaneously:initially without permeating fluids into the rock,to seed fractures uniformly and at high density in the wall rock of the wells;followed by fluid injection to propagate the seeded fracture density hydraulically.FBCM combines the ease of continuum modeling with the potential accuracy of modeling discrete fractures and fracturing explicitly.Fractures are modeled as piecewise planar based on intersections with domain elements;fracture geometry stored as continuum properties is used to calculate parameters needed to model individual fractures;and rock behavior is modeled through tensorial aggregation of the behavior of discrete fractures and unfractured rock.Simulations are presented for previously unfractured rock and for rock with preexisting fractures of horizontal,shallow-dipping,steeply dipping,or vertical orientation.Simulations of a single-well model are used to determine the pattern and spacing for a multiple-well design.The results illustrate high-density fracturing and fracture caging through simultaneous fluid injection in multiple wells:for previously unfractured rock or rock with preexisting shallow-dipping or horizontal fractures,and in situ vertical compressive stress greater than horizontal.If preexisting fractures are steeply dipping or vertical,and considering the same in situ stress condition,well pressurization without fluid permeation appears to be the only practical way to induce new fractures and contain fracturing within the target domain.