Discrete fracture model for simulating waterflooding processes under fracturing conditions

In our study, we develop a model for simulating fracturing processes in a poroelastic medium. The proposed approach combines the discrete fracture model enriched with contact plane mechanics. The model captures mechanical interactions of fractures and a deformable medium, fluid, and heat transfer in fractures and in a porous medium. Both effects of poroelasticity and thermoelasticity are accounted in our model. Mass and heat conservation equations are approximated by the finite volume method, and mechanical equilibrium equations are discretized by means of the Galerkin finite element approach. Two‐dimensional grid facets between 3‐dimensional finite elements are considered as possible fracture surfaces. Most of these facets are inactive from the beginning and are activated throughout the simulation. A fracture propagation criterion, based on Irwin's approach, is verified on each nonlinear iteration. When the criterion is satisfied, additional contact elements are added into finite element and discrete fracture model formulations respectively. The proposed approach allows modeling of existing natural and artificially created fractures within one framework. The model is tested on single‐ and multiple‐phase fluid flow examples for both isothermal and thermal conditions and verified against existing semianalytical solutions. The applicability of the approach is demonstrated on an example of practical interests where a sector model of an oil reservoir is simulated with different injection and production regimes.     

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.     

Numerical Modeling of Natural Fracture Distributions in Shale

The production efficiency of shale gas is affected by the interaction between hydraulic and natural fractures. This study presents a simulation of natural fractures in shale reservoirs, based on a discrete fracture network (DFN) method for hydraulic fracturing engineering. Fracture properties of the model are calculated from core fracture data, according to statistical mathematical analysis. The calculation results make full use of the quantitative information of core fracture orientation, density, opening and length, which constitute the direct and extensive data of mining engineering. The reliability and applicability of the model are analyzed with regard to model size and density, a calculation method for dominant size and density being proposed. Then, finite element analysis is applied to a hydraulic fracturing numerical simulation of a shale fractured reservoir in southeastern Chongqing. The hydraulic pressure distribution, fracture propagation, acoustic emission information and in situ stress changes during fracturing are analyzed. The results show the application of fracture statistics in fracture modeling and the influence of fracture distribution on hydraulic fracturing engineering. The present analysis may provide a reference for shale gas exploitation.     

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

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

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.     

基于ABAQUS平台的水力裂缝扩展有限元模拟研究

随着扩展有限元理论的深入研究,利用扩展有限元方法模拟水力压裂具有了一定的可操作性。相比于常规有限元方法,XFEM方法具有计算结果精度高和计算量小的优点。但是,如何模拟射孔孔眼、如何模拟流体与岩石相互作用以及分析水力裂缝的扩展规律仍然是难题。以研究水力压裂裂缝扩展规律为目的,建立了岩石多孔介质应力平衡方程、流体渗流连续性方程和边界条件。通过有限元离散化方法对耦合方程矩阵进行处理。通过富集函数定义初始裂缝（射孔孔眼）,选择最大主应力及损伤变量D分别作为裂缝起裂和扩展判定准则,利用水平集方法模拟水力裂缝扩展过程。数值模拟结果显示：增加射孔方位角、压裂液排量和减小水平地应力差,起裂压力上升;黏度对起裂压力无明显影响。增加射孔方位角、压裂液排量、黏度和减小水平地应力差值有助于裂缝宽度的增加。增加水平地应力差值、压裂液排量和减小射孔方位角以及压裂液黏度有助于裂缝长度增加,反之亦然。基于ABAQUS的水力裂缝扩展有限元法可对不同井型和诸多储层物性参数及压裂施工参数进行分析,且裂缝形态逼真,裂缝面凹凸程度清晰,结果准确。此研究可作为一种简便有效研究水力压裂裂缝扩展规律的方法为油田水力压裂设计与施工提供参考与依据。     

Anisotropic porothermoelastic solution and hydro‐thermal effects on fracture width in hydraulic fracturing

In this study, the effects of the temperature difference between hydraulic fracturing fluid and rock formation on the time‐dependent evolution of fracture width were investigated using a newly derived one‐dimensional anisotropic porothermoelastic analytical solution. The solution is shown to correctly reproduce existing solutions for special cases and corrections for an earlier publication are provided. An analysis of time‐dependent fracture width evolution using Woodford Shale data was also presented. It was found that when the fracturing fluid has the same temperature as the shale formation, the fracture gradually closes back after the initial opening due to the invasion of the fracturing fluid. Practically, in this scenario, proppants should be pumped into the fracture as soon as possible to obtain maximum fracture conductivity. On the other hand, with a fracturing fluid 60 °C colder than the formation, the thermal contraction of the rock dominates the fracture aperture evolution, resulting in a fracture aperture approximately 70% larger than that produced by the hotter fracturing fluid. Consequently, in this case, it is beneficial to delay proppant placement to take advantage of the widening fractures. Finally, it was found that the fracture aperture is directly controlled by the spacing of natural fractures. Therefore, the presence of natural fractures in the shale formation and their spacing influence not only the type of hydraulic fractures created but also what kind and size of proppants should be used to keep them open. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical modeling of stress-permeability coupling in rough fractures

A numerical model is described for coupled flow and mechanical deformation in fractured rock. The mechanical response of rock joints to changes in hydraulic pressure is strongly influenced by the geometric characteristics of the joint surfaces. The concept of this work is to combine straightforward finite element solutions with complex and realistic fracture surface geometry in order to reproduce the non-linear stress-deformation-permeability coupling that is commonly observed in fractures. Building on the numerous studies that have expanded the understanding of the key parameters needed to describe natural rough-walled fractures, new methods have been developed to generate a finite element mesh representing discrete fractures with realistic rough surface geometries embedded in a rock matrix. The finite element code GeoSys/Rockflow was then used to simulate the coupled effects of hydraulic stress, mechanical stress, and surface geometry on the evolving permeability of a single discrete fracture. The modeling concept was experimentally verified against examples from the literature. Modeling results were also compared to a simple interpenetration model.     

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

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

Simulation of non‐planar three‐dimensional hydraulic fracture propagation

Hydraulic fracturing is the method of choice to enhance reservoir permeability and well efficiency for extraction of shale gas. Multi‐stranded non‐planar hydraulic fractures are often observed in stimulation sites. Non‐planar fractures propagating from wellbores inclined from the direction of maximum horizontal stress have also been reported. The pressure required to propagate non‐planar fractures is in general higher than in the case of planar fractures. Current computational methods for the simulation of hydraulic fractures generally assume single, symmetric, and planar crack geometries. In order to better understand hydraulic fracturing in complex‐layered naturally fractured reservoirs, fully 3D models need to be developed. In this paper, we present simulations of 3D non‐planar fracture propagation using an adaptive generalized FEM. This method greatly facilitates the discretization of complex 3D fractures, as finite element faces are not required to fit the crack surfaces. A solution strategy for fully automatic propagation of arbitrary 3D cracks is presented. The fracture surface on which pressure is applied is also automatically updated at each step. An efficient technique to numerically integrate boundary conditions on crack surfaces is also proposed and implemented. Strongly graded localized refinement and analytical asymptotic expansions are used as enrichment functions in the neighborhood of fracture fronts to increase the computational accuracy and efficiency of the method. Stress intensity factors with pressure on crack faces are extracted using the contour integral method. Various non‐planar crack geometries are investigated to demonstrate the robustness and flexibility of the proposed simulation methodology. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

危岩水力劈裂分析的扩展有限元法

危岩是三峡库区典型的地质灾害类型之一,而主控结构面受荷断裂扩展是危岩发育成灾的关键核心。将危岩主控结构面类比为宏观裂纹,利用扩展有限元法在模拟裂纹扩展方面的优势,基于考虑裂纹面水压力作用的虚功原理推导出了采用扩展有限元法分析水力劈裂问题的控制方程,给出了危岩主控结构面水力劈裂问题的扩展有限元实现方法,对重庆万州太白岩危岩主控结构面的水力劈裂进行了数值模拟分析。计算结果表明：暴雨是威胁危岩稳定性的最敏感因素,随着裂隙水压力上升,裂端拉应力会急剧升高,危岩的稳定性降低;I型裂纹扩展是危岩主要的结构面扩展形式,结构面一旦发生开裂,将处于非稳定扩展状态。     

水力劈裂对岩体中自然裂纹的影响研究

建立了求解自然裂纹和水力裂纹扩展的扩展有限元法,对裂纹附近区域的节点采用广义形函数,并采用线增函数消除混合单元,以提高裂纹附近的精度。引入水力劈裂的非耦合模型,即假设裂纹中的水压力为均布力;用砂浆法（线段-线段接触法）结合增广型拉格朗日乘子法处理受压裂纹段的接触条件。并通过算例分析了以下内容：计算了受压裂纹和裂纹面分布均布水压力的水力裂纹的应力强度因子,并与解析解进行了比较,结果表明,提出的方法具有很高的精度;模拟了水力裂纹对自然裂纹面的影响,并分析了自然裂纹面上的接触力和接触状态。     

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

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

An algorithm for the simulation of curvilinear plane-strain and axisymmetric hydraulic fractures with lag using the universal meshes

We present an algorithm to simulate curvilinear hydraulic fractures in plane strain and axisymmetry. We restrict our attention to sharp fractures propagating in an isotropic, linear elastic medium and driven by the injection of a laminar, Newtonian fluid governed by lubrication theory, and we require the existence of a finite lag region between the fluid front and the crack tip. The key novelty of our approach is in how we discretize the evolving crack and fluid domains: we utilize universal meshes (UMs), a technique to create conforming triangulations of a problem domain by only perturbing nodes of a universal background mesh in the vicinity of the boundary. In this way, we construct meshes, which conform to the crack and to the fluid front. This allows us to build standard piecewise linear finite element spaces and to monolithically solve the quasistatic hydraulic fracture problem for the displacement field in the rock and the pressure in the fluid. We demonstrate the performance of our algorithms through three examples: a convergence study in plane strain, a comparison with experiments in axisymmetry, and a novel case of a fracture in a narrow pay zone.     

不同尺度岩体结构面对页岩气储层水力压裂裂缝扩展的影响

储层岩体中的天然结构面对水力压裂缝网改造具有重要的影响。本文采用真实破裂过程分析软件RFPA2D-Flow,在考虑岩体非均质性和岩体渗流-应力-损伤破裂特性的基础上,对不同尺度天然结构面影响的水力压裂裂缝扩展与演化行为进行了模拟分析和讨论,研究结果表明:（1）当水力裂缝遇天然非闭合裂隙时,在水力裂缝靠近非闭合裂隙区间形成拉张应力区,水力裂缝与区间非闭合裂隙间微元体累进性张拉破坏是导致水力裂缝与非闭合裂隙贯通的主要机制;（2）层理等优势结构对水力压裂裂缝扩展及缝网形态影响十分显著,当最大主应力方向与层理面走向小角度相交时,层理结构面对水力裂隙的扩展起主要作用,当最大主应力方向与层理面走向大角度相交时,最大主压应力与层理面共同对缝网扩展起主导作用,随着优势结构面的增多和差应力的增大,水力压裂形成的缝网范围和复杂性程度随之增大;（3）储层水力压裂是一种局部范围内的短暂动力扰动过程,尽管断层的存在可以极大地影响水力裂缝的扩展模式,增大水力裂隙扩展高度,但相比于储层埋深,水力压裂对断层封闭性的破坏范围和断层活动性的扰动程度十分有限。     

Simulations of hydro-fracking in rock mass at meso-scale using fully coupled DEM/CFD approach

The paper deals with two-dimensional (2D) numerical modelling of hydro-fracking (hydraulic fracturing) in rocks at the meso-scale. A numerical model was developed to characterize the properties of fluid-driven fractures in rocks by combining the discrete element method (DEM) with computational fluid dynamics (CFD). The mechanical behaviour of the rock matrix was simulated with DEM and the behaviour of the fracturing fluid flow in newly developed and pre-existing fractures with CFD. The changes in the void geometry in the rock matrix were taken into account. The initial 2D hydro-fracking simulation tests were carried out for a rock segment under biaxial compression with one injection slot in order to validate the numerical model. The qualitative effect of several parameters on the propagation of a hydraulic fracture was studied: initial porosity of the rock matrix, dynamic viscosity of the fracking fluid, rock strength and pre-existing fracture. The characteristic features of a fractured rock mass due to a high-pressure injection of fluid were realistically modelled by the proposed coupled approach.

     

Finite element simulations of 3D planar hydraulic fracture propagation using a coupled hydro-mechanical interface element

Two-dimensional hydraulic fracturing simulations using the cohesive zone model (CZM) can be readily found in the literature; however, to our knowledge, verified 3D cohesive zone modeling is not available. We present the development of a 3D fully coupled hydro-mechanical finite element method (FEM) model (with parallel computation framework) and its application to hydraulic fracturing. A special zero-thickness interface element based on the CZM is developed for modeling fracture propagation and fluid flow. A local traction-separation law with strain softening is used to capture tensile cracking. The model is verified by considering penny-shaped hydraulic fracture and plain strain Kristianovich‑Geertsma‑de Klerk hydraulic fracture (in 3D) in the viscosity- and toughness-dominated regimes. Good agreement between numerical results and analytical solutions has been achieved. The model is used to investigate the influence of rock and fluid properties on hydraulic fracturing. Lower stiffness tip cohesive elements tend to yield a larger elastic deformation around the fracture tips before the tensile strength is reached, generating a larger fracture length and lower fracture pressure compared with higher stiffness elements. It is found that the energy release rate has almost no influence on hydraulic fracturing in the viscosity-dominated regime because the energy spent in creating new fractures is too small when compared with the total input energy. For the toughness-dominated regime, the released energy during fracturing should be accurately captured; relatively large tensile strength should be used in order to match numerical results to the asymptotic analytical solutions. It requires smaller elements when compared with those used in the viscosity-dominated regime.     

不同水力加载条件下非均质储层缝网扩展规律研究

水力压裂可显著提高页岩气等致密储层岩体的渗透性以增加油气产量,然而受多因素影响,水力压裂形成缝网结构的机理和压裂优化设计一直是研究的焦点和难点。本研究基于渗流-应力-破坏耦合计算模拟方法,对不同水力加载条件下的非均质储层水力压裂过程进行了模拟和对比研究。研究结果表明:水力压裂过程中起始注水压力和增量大小对水力压裂缝网扩展和改造区域形态有着显著的影响。当起始注水压力小于等于模型材料体抗拉强度,并缓慢增压致裂时,压裂过程可近似视为稳态应力-破坏-渗流耦合作用过程的不同阶段,这种情况下仅在压裂井孔周围形成两组对称式的伞状水力裂缝带。当对模型体施加高于模型材料体破裂压力的注水压力时,相当于对压裂孔快速施加高动水压力,水力裂缝沿压裂孔全方位迅速萌生并快速扩展,当注水压力值高于破裂压力一定幅值时,压裂改造可形成围绕压裂井全方位的放射状裂缝网络,使压裂储层得以最大范围改造。在拟静力和拟动力两种加载条件下,不同水岩相互作用机理是造成不同水力加载条件出现不同缝网结构的力学机制,而对于实际的页岩气储层改造,压裂产生围绕压裂井全方位放射状的缝网结构则是一种最优的体积压裂改造。     

Experiments and Analysis of the Hydraulic Fracture Propagation Behaviors of the Granite with Structural Planes in the Gonghe Basin

Hydraulic fracture (HF) propagation behavior is signi?cant when building enhanced geothermal systems (EGS). HF geometry is closely related to the structural planes (SPs) in hot dry rock (HDR), such as natural fractures (NFs), quartz veins (QVs) and lithologic interfaces (LIs). However, the HF behaviors in HDR have not been well understood, especially the influence of multiple SPs on the HF geometry. To clarify this mechanism, several groups of physical simulation experiments of hydraulic fracturing were conducted to investigate the intersection relationship between the HFs and the SPs. Results show that the HF geometry shows great differences when intersecting with different SPs. In summary, the HF geometry displays four basic patterns, namely, propagation along the SPs, branching, capture, penetration/non-dilation. The fluctuation degree of the pressure-time curve and the HF complexity show a positive correlation. The cementing strength of the SP and their different mechanical properties from rock matrix influence the HF behaviors significantly. Therefore, the HF shows diverse geometries when intersecting with the NFs and LIs, while propagating along the QV when intersecting with it. For the complex networks, it is favorable for the HF to penetrate through and dilate several SPs, rather than simply cross or propagate along the SP.