基于PFC2D数值模拟的交替压裂中应力阴影效应研究

为了探究应力阴影效应对交替压裂中压裂间距选取的影响,基于优化后的颗粒流离散元流固耦合计算模型,模拟并分析了双初始水力裂缝下因应力阴影效应产生的诱导应力的分布情况,并与理论解析解进行对比,证明了该数值方法的合理性。在此基础上,分析了应力阴影效应在不同各向异性地应力场及初始压裂间距条件下对新水力裂缝的起裂压力及扩展形态的影响,研究结果表明：初始各向异性应力场不改变裂缝周边的应力场,不影响新水力裂缝的起裂压力;随着初始压裂间距的减小,应力阴影效应增强,新水力裂缝的起裂压力逐渐增加。初始水力裂缝间距与初始各向异性应力场共同影响新水力裂缝的扩展形态,随着初始水力裂缝间距或初始水平地应力场差异系数的增大,应力阴影对新水力裂缝的扩展方向的影响逐步减弱;初始水力裂缝对新水力裂缝的扩展有一定的限制作用,在一定程度上不利于形成复杂的裂缝网络。根据以上分析结果,对交替压裂中压裂间距的优化进行了定性的探讨。     

Numerical models of hydraulic fracturing and the interpretation of syntectonic veins

The hydraulic fracture is modelled as an ellipse in an infinite elastic medium with an internal fluid pressure and loaded under biaxial stresses at infinity. The available stress function for this model has been evaluated numerically, and the magnitudes of the stresses generated around the crack calculated for a variety of loading conditions and crack orientations. Fracture initiation is predicted from the Griffith maximum tensile stress criterion. The location of the maximum tensile stress around the crack is recorded and it is found that for many conditions of applied stresses and crack fluid pressure, the hydraulic shear fracture has a symmetrically developed maximum tensile stress and fracture initiation will occur by growth along the direction of the crack. It is also predicted that fracture initiation will occur when the ratio of fluid pressure to applied least principal stress is considerably less than one. The elastic strain energy fields around elliptical hydraulic flaws have been calculated, and in particular, the change in strain energy upon introduction of a small flaw, and the change in strain energy upon growth of this flaw, have been investigated. The results allow an evaluation of the second part of the Griffith criterion-that fracture growth is accompanied by a decrease in strain energy-for hydraulic fractures. Changes in strain energy with small increases in fluid pressure provide a physical basis for dilatancy hardening and fracture instability. Quasi-static growth from a flaw is modelled by calculating changes in strain energy for unit increases in half length. The distinction between fractures which show an increasing and a decreasing rate of change in strain energy with increasing length, and between fractures which may only extend spontaneously for short distances and those which may show extensive spontaneous growth on the basis of the rate of change of strain energy with length, is made. A gradual drop in crack fluid pressure once the threshold for fracture initiation has been passed may promote the extent to which spontaneous crack growth occurs.The formation of syntectonic veins, particularly in rocks being deformed under low grade metamorphic conditions, is often the most abundant evidence of natural hydraulic fracturing in rocks. Commonly observed geometric features of syntectonic veins-length, simple tapering, symmetric and asymmetric forking, branching, irregular zig-zag traces, en échelon patterns—are discussed primarily with reference to the strain energy model for growth established, and the geometric variation is interpreted in terms of variation in applied stress and fluid pressure conditions and the rate of change of stored strain energy with crack growth. In particular, terminal branching arises when the minimum stress changes from a symmetric to an asymmetric location at the tip of a growing shear fracture, and terminal forking results when there is an increase in the energy release rate during crack growth, and may be symmetric or asymmetric depending on the location of the minimum stress at the crack tip at the time of forking.     

Numerical Modeling of Hydraulic Fractures Interaction in Complex Naturally Fractured Formations

A recently developed unconventional fracture model (UFM) is able to simulate complex fracture network propagation in a formation with pre-existing natural fractures. A method for computing the stress shadow from fracture branches in a complex hydraulic fracture network (HFN) based on an enhanced 2D displacement discontinuity method with correction for finite fracture height is implemented in UFM and is presented in detail including approach validation and examples. The influence of stress shadow effect from the HFN generated at previous treatment stage on the HFN propagation and shape at new stage is also discussed.     

Fully coupled simulation of a hydraulic fracture interacting with natural fractures with a hybrid discrete‐continuum method

A hybrid discrete‐continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity‐dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T‐shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs. Copyright © 2017 John Wiley & Sons, Ltd.     

水平井中多条裂缝同步扩展时裂缝竞争机制

水平井中多条水力裂缝间的应力干扰行为,造成了压裂液排量的非均匀分配,影响了水力裂缝的几何形态。采用边界元法研究岩体在压裂液作用下的变形程度,以幂律流体泊肃叶平板流动方程来研究水力裂缝内部的压裂液流场,考虑了多条裂缝间应力干扰和压裂液流量分配,建立了流-固耦合的水平井多条水力裂缝同步扩展模型。模型可模拟水平井多条水力裂缝几何形态、应力干扰情况和压裂液排量的分配情况,可解释水力裂缝之间的竞争机制。多条裂缝同步扩展时,压裂液排量并非均等地分配到各个裂缝之中,进入到内部裂缝的压裂液流量最小,内部裂缝宽度最小;内部的水力裂缝增长一定长度后,停止增长,并且在应力干扰下逐渐闭合。     

青海共和盆地干热岩地热储层水力压裂物理模拟和裂缝起裂与扩展形态研究

水力压裂是青海共和盆地干热岩地热资源开发的难点技术问题之一。本文基于升级改造的大尺寸真三轴水力压裂物理模拟实验系统模拟干热岩储层高温高压环境,利用青海共和盆地露头岩心进行水力压裂物理模拟实验,揭示干热岩储层水力裂缝的起裂和扩展规律。通过物理模拟实验发现：干热岩储层裂缝起裂可以通过文中提出的起裂模型判断起裂方式和预测起裂压力;水力裂缝在岩石基质中的扩展形态简单,仅沿最大主应力方向延伸;但是水力裂缝会受到岩石中弱面的影响,发生转向沿弱面延伸,形成较复杂的裂缝形态。因此,建议在干热岩储层实际施工中,在天然裂缝发育较丰富的层段开展水力压裂,以实现复杂裂缝网络提取地热能。     

Analysis of a deformable fracture in permeable material

The response of deformable fractures to changes in fluid pressure controls phenomena ranging from the flow of fluids near wells to the propagation of hydraulic fractures. We developed an analysis designed to simulate fluid flows in the vicinity of asperity‐supported fractures at rest, or fully open fractures that might be propagating. Transitions between at‐rest and propagating fractures can also be simulated. This is accomplished by defining contact aperture as the aperture when asperities on a closing fracture first make contact. Locations on a fracture where the aperture is less than the contact aperture are loaded by both fluid pressure and effective stress, whereas locations where the aperture exceeds the contact aperture are loaded only by fluid pressure. Fluid pressure and effective stress on the fracture are determined as functions of time by solving equations of continuity in the fracture and matrix, and by matching the global displacements of the fracture walls to the local deformation of asperities. The resulting analysis is implemented in a numerical code that can simulate well tests or hydraulic fracturing operations. Aperture changes during hydraulic well tests can be measured in the field, and the results predicted using this analysis are similar to field observations. The hydraulic fracturing process can be simulated from the inflation of a pre‐existing crack, to the propagation of a fracture, and the closure of the fracture to rest on asperities or proppant. Two‐dimensional, multi‐phase fluid flow in the matrix is included to provide details that are obscured by simplifications of the leakoff process (Carter‐type assumptions) used in many hydraulic fracture models. Execution times are relatively short, so it is practical to implement this code with parameter estimation algorithms to facilitate interpretation of field data. Copyright © 2006 John Wiley & Sons, Ltd.     

川南深层页岩水力压裂缝网扩展规律数值模拟研究

3500 m以深页岩气资源量占整个川南地区总资源量的比例高达86.5%,该区深层页岩气藏构造复杂,压裂形成复杂缝网的难度大,有必要通过数值模拟研究深层页岩气复杂缝网主控因素,对实现川南地区深层页岩气的效益开发具有重要意义。在对川南地区页岩气气田某井的岩芯进行细观尺度下的观察并构建二维裂缝模型的基础上,利用位移间断边界元法(DDM)模拟深层页岩水力压裂过程中水力裂缝与天然裂缝相互作用的物理力学过程,研究主应力、应力差和压裂液排量对裂缝扩展的影响。结果表明:在高应力差条件下缝网的复杂程度和总长度急剧降低,缝网的平均宽度增大,且平均宽度随排量增加而增大的能力变得有限。在高应力差条件下提升压裂液排量,缝网长度的增加以产生新生裂缝为主,同时提升排量对于激活天然裂缝有一定的提升作用,但是效果有限。相比于拉张裂缝,剪切裂缝的形成受主应力和压裂液排量的影响更显著,在高应力差条件下缝网中剪切裂缝的长度急剧降低。随着压裂液的注入,在较低应力差和相同压裂液注入量的情况下,低排量工况下的裂缝长度逐渐大于高排量工况下的裂缝长度。在应力差较高的情况下裂缝扩展的速率较低,同时会使提升排量而形成更多新生裂缝的能力变得有限。基于显微镜观察构建的裂缝模型计算出的结果能够较好地符合场地实际,为深层水力压裂设计和施工提供参考。     

Flow–stress coupled permeability tensor for fractured rock masses

In this paper a new analytical model is proposed to determine the permeability tensor for fractured rock masses based on the superposition principle of liquid dissipation energy. This model relies on the geometrical characteristics of rock fractures and the corresponding fracture network, and demonstrates the coupling effect between fluid flow and stress/deformation. This model empirically considers the effect of pre‐peak shear dilation and shear contraction on the hydraulic behavior of rock fractures and can be used to determine the applicability of the continuum approach to hydro‐mechanical coupling analysis. Results of numerical analysis presented in this paper show that the new model can effectively describe the permeability of fractured rock masses, and can be applied to the coupling analysis of seepage and stress fields. Copyright © 2008 John Wiley & Sons, Ltd.     

岩石裂隙的非饱和渗透特性及其演化规律研究

岩石裂隙的非饱和渗透特性是岩土、能源和环境等领域科学研究中的热点问题。采用三维激光扫描获取花岗岩裂隙的表面形貌特征,分析裂隙微观形貌特征对非饱和渗透特性的影响。研究在张拉、压缩、剪切等复杂荷载作用下裂隙开度分布的演化规律,建立复杂荷载作用下岩石裂隙非饱和毛细压力曲线演化模型。基于裂隙的微观形貌特征推导了岩石裂隙非饱和相对渗透系数模型,通过与试验数据对比,验证了模型的准确性和有效性,并在此基础上建立了复杂荷载作用下岩石裂隙非饱和相对渗透系数演化模型。研究成果对非饱和条件下裂隙岩体的水-力耦合机制研究具有一定指导意义。     

基于扩展有限元法的碳酸盐岩地热储层岩体裂缝扩展行为

水力压裂作为一种主要的地热能开采手段,其压裂效果除与岩体基本物理力学性质有关外,还与裂隙分布、地应力状态、压裂工程参数等密切相关.为了探究以上因素对水力压裂过程中裂缝扩展行为的影响,以冀中坳陷碳酸盐岩储层岩体为研究对象,基于扩展有限元法,建立裂缝扩展流固耦合模型,分析了水平应力差、射孔方位角、注入液排量和压裂液黏度等参数对裂缝扩展行为的影响.结果表明:单裂缝扩展时,射孔方位角越小、注入量越大、越有利于裂缝扩展;双裂缝扩展时,水平应力差增大,裂缝偏转程度变小;水力裂缝与天然裂缝相交时,较小水平应力差有利于天然裂缝开启.      

Analysis of refracturing in horizontal wells: Insights from the poroelastic displacement discontinuity method

In this paper, a fully coupled 2‐dimensional poroelastic displacement discontinuity method is used to investigate the refracturing process in horizontal wells. One of the objectives of refracturing is to access new reserves by adding new hydraulic fractures in zones that were bypassed in the initial fracturing attempt. Pore pressure depletion in the vicinity of old fractures directly affects the state of stress and eventually the propagation of newly created hydraulic fractures. Thus, a poroelastic analysis is required to identify guidelines for the refracturing process, in particular to understand the extension of the pore pressure depletion, and eventually, the orientation of new as well as old fractures. We propose a fully coupled approach to model the whole process of child fracture propagation in a depleted area between 2 parent fractures in the same wellbore. This approach omits the need of using multistep workflow that is regularly used to model the process. The maximum tensile stress criterion (σ criterion) is used for hydraulic fracture propagation. The proposed method is verified using available analytical solutions for total stress and pore pressure loading modes on a line fracture in drained and undrained conditions. Then, test cases of multifractured horizontal wells are studied to calculate the time evolution of the stress and pore pressure fields around old fractures and to understand the effect of these fields on the propagation path of newly created fractures. Finally, the effect of the pore pressure depletion on the propagation path of the newly created fractures in the bypassed area of the same wellbore is studied. The results show that the depleted areas around old fractures are highly affected by the extent and severity of the stress redistribution and pore pressure depletion. It is observed that a successful creation of new fractures may only happen in certain time frames. The results of this study provide new insights on the behavior of newly created fractures in depleted zones. They also clarify the relationship between stress change and pore pressure depletion in horizontal wells.     

Development of a full 3D numerical model to investigate the hydraulic fracture propagation under the impact of orthogonal natural fractures

Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.

     

Overview of numerical models for interactions between hydraulic fractures and natural fractures: Challenges and limitations

The intersection of natural fractures with hydraulic fractures results in formation of complex fracture networks, including non-planar fractures or multi-stranded fractures. On one hand, opening of these natural fractures improves productivity of the formation; on the other hand, coalescence of these fractures into a hydraulic fracture makes pressure analysis and prediction of fracture growth very complicated. Overall, interactions between natural fractures and hydraulic fractures pose more challenges in the fracturing design and its execution. Investigation and understanding of their interaction are crucial in achieving successful fracture treatments in formations with pre-existing natural fracture network. In this paper, we will review the numerical works that have been done in the last decade to model opening of natural fractures during hydraulic fracturing, focusing especially on mechanical models that address propagation of hydraulic fractures in naturally fractures reservoirs. Linear elastic fracture mechanics, cohesive element methods and continuum damage mechanics techniques utilized to understand interaction of hydraulic fractures with natural fractures are discussed here based on their capability to reproduce experimental results and field observations.     

金湖凹陷天96断块阜二段储层裂缝多参数定量预测

基于应变能、表面能理论, 采用格里菲斯岩石破裂准则, 建立了金湖凹陷天96断块造缝期古地应力、现今地应力与裂缝参数之间的关系; 基于构造分析确定了造缝时期古地应力的类型和方向, 用声速法和水力压裂法确定了现今地应力的方向和大小, 将古今地应力场的数值模拟结果导入裂缝参数计算模型中, 模拟计算了天96断块阜二段裂缝的密度、开度、孔渗参数, 并结合古地应力场分析以及岩心裂缝统计数据对裂缝预测结果进行验证; 最后, 综合分析了影响天96断块裂缝发育的因素。      

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

水力压裂扩展特性的数值模拟研究

采用ABAQUS建立了水力压裂计算模型,模拟了地应力、岩石力学特性、压裂液流体特性等各种复杂因素对水力压裂扩展的影响。通过计算分析得到一些有益结论：（1）在注入压力一定的情况下起裂压力与最小水平地应力、临界应力、初始孔隙压力成正比,而与压裂液黏度、最大水平地应力、弹性模量无关;（2）裂缝扩展长度和最大缝宽与最小水平地应力、初始孔隙压力、弹性模量成反比,而与最大水平地应力无关;（3）水力压裂作业中,缝长的扩展过程可分为无扩展阶段、快速扩展阶段、稳定扩展阶段以及缓慢扩展阶段等4个阶段。研究结论对于水力压裂作业优化具有参考价值。     

页岩储层体积压裂裂缝扩展机制研究

页岩储层天然裂缝、水平层理发育,水力压裂过程中可能形成复杂的体积裂缝。针对页岩储层体积裂缝扩展问题,基于流-固耦合基本方程和损伤力学原理,建立了页岩储层水力压裂体积裂缝扩展的三维有限元模型。将数值模型的模拟结果与页岩储层裂缝扩展室内试验结果进行对比,二者吻合较好,从而证明了数值模型的可靠性。通过一系列数值模拟发现：（1）水力压裂过程中水平层理可能张开,形成水平缝,水平与垂直缝相互交错,形成复杂的体积裂缝网络;（2）水平主应力差增大,体积裂缝的分布长度（水平最大主应力方向压裂裂缝的展布距离）增加、分布宽度（水平最小主应力方向压裂裂缝的展布距离）减小,体积裂缝的长宽比增加;（3）压裂施工排量增大,体积裂缝的分布长度减小、宽度增加,压裂裂缝的长宽比降低;（4）天然裂缝的残余抗张强度增大,体积裂缝分布宽度减小、分布长度增加,体积裂缝的长宽比增加。研究成果可以为国内的页岩气的压裂设计和施工提供一定的参考和借鉴。     

Numerical modeling of hydraulic fracture propagation,closure and reopening using XFEM with application to in‐situ stress estimation

In this paper, a fully coupled model is developed for numerical modeling of hydraulic fracturing in partially saturated weak porous formations using the extended finite element method, which provides an effective means to simulate the coupled hydro‐mechanical processes occurring during hydraulic fracturing. The developed model is for short fractures where plane strain assumptions are valid. The propagation of the hydraulic fracture is governed by the cohesive crack model, which accounts for crack closure and reopening. The developed model allows for fluid flow within the open part of the crack and crack face contact resulting from fracture closure. To prevent the unphysical crack face interpenetration during the closing mode, the crack face contact or self‐contact condition is enforced using the penalty method. Along the open part of the crack, the leakage flux through the crack faces is obtained directly as a part of the solution without introducing any simplifying assumption. If the crack undergoes the closing mode, zero leakage flux condition is imposed along the contact zone. An application of the developed model is shown in numerical modeling of pump‐in/shut‐in test. It is illustrated that the developed model is able to capture the salient features bottomhole pressure/time records exhibit and can extract the confining stress perpendicular to the direction of the hydraulic fracture propagation from the fracture closure pressure. Copyright © 2016 John Wiley & Sons, Ltd.     

An explicitly coupled hydro‐geomechanical model for simulating hydraulic fracturing in arbitrary discrete fracture networks

Modeling hydraulic fracturing in the presence of a natural fracture network is a challenging task, owing to the complex interactions between fluid, rock matrix, and rock interfaces, as well as the interactions between propagating fractures and existing natural interfaces. Understanding these complex interactions through numerical modeling is critical to the design of optimum stimulation strategies. In this paper, we present an explicitly integrated, fully coupled discrete‐finite element approach for the simulation of hydraulic fracturing in arbitrary fracture networks. The individual physical processes involved in hydraulic fracturing are identified and addressed as separate modules: a finite element approach for geomechanics in the rock matrix, a finite volume approach for resolving hydrodynamics, a geomechanical joint model for interfacial resolution, and an adaptive remeshing module. The model is verified against the Khristianovich–Geertsma–DeKlerk closed‐form solution for the propagation of a single hydraulic fracture and validated against laboratory testing results on the interaction between a propagating hydraulic fracture and an existing fracture. Preliminary results of simulating hydraulic fracturing in a natural fracture system consisting of multiple fractures are also presented. Copyright © 2012 John Wiley & Sons, Ltd.