Modeling spatial fracture intensity as a control on flow in fractured rock

Spatial fracture intensity (P ₃₂, fracture area by volume) is an important characteristic of a jointed rock mass. Although it can hardly ever be measured, P ₃₂ can be modeled based on available geological information such as spatial data of the fracture network. Flow in a mass composed of low-permeability hard rock is controlled by joints and fractures. In this article, models were developed from a geological data set of fractured andesite in LanYu Island (Taiwan) where a site is investigated for possible disposal of low-level and intermediate-level radionuclide waste. Three different types of conceptual models of spatial fracture intensity distribution were generated, an Enhanced Baecher’s model (EBM), a Levy–Lee Fractal model (LLFM) and a Nearest Neighborhood model (NNM). Modeling was conducted on a 10 × 10 × 10 m synthetic fractured block. Simulated flow was forced by a 1% hydraulic gradient between two vertical x–z faces of the cube (from North to South) with other boundaries set to no-flow conditions. Resulting flow vectors are very sensitive to spatial fracture intensity (P ₃₂). Flow velocity increases with higher fracture intensity (P ₃₂). R-squared values of regression analysis for the variables velocity (V/V _max) and fracture intensity (P ₃₂) are 0.293, 0.353, and 0.408 in linear fit and 0.028, 0.08, and 0.084 in power fit. Higher R ² values are positively linked with structural features but the relation between velocity and fracture intensity is non-linear. Possible flow channels are identified by stream-traces in the Levy–LeeFractal model.     

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.     

基于离散裂隙网络模型的裂隙水渗流计算

离散裂隙网络模型(Discrete Fracture Network（DFN）)是研究裂隙水渗流最为有效的手段之一。文章根据裂隙几何参数和水力参数的统计分布,利用Monte Carlo随机模拟技术生成二维裂隙网络,基于图论无向图的邻接矩阵判断裂隙网络的连通,利用递归算法提取出裂隙网络的主干网或优势流路径。基于立方定律和渗流连续性方程,利用数值解析法建立了二维裂隙网络渗流模型,分析不同边界条件下裂隙网络中的流体流动。结果表明,该方法可以模拟区域宏观水力梯度和边界条件下,裂隙网络水力梯度方向总的流量,以及节点的水位、节点间的流量和流动方向的变化特征,为区域岩溶裂隙水渗流计算提供了一种实用、可行的方法。      

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.     

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.     

An embedded fracture modeling framework for simulation of hydraulic fracturing and shear stimulation

A numerical modeling framework is described that is able to calculate the coupled processes of fluid flow, geomechanics, and rock failure for application to general engineering problems related to reservoir stimulation, including hydraulic fracturing and shear stimulation. The numerical formulation employs the use of an embedded fracture modeling approach, which provides several advantages over more traditional methods in terms of computational complexity and efficiency. Specifically, the embedded fracture modeling strategy avoids the usual requirement that the discretization of the fracture domain conforms to the discretization of the rock volume surrounding the fractures. As fluid is exchanged between the two domains, conservation of mass is guaranteed through a coupling term that appears as a simple source term in the governing mass balance equations. In this manner, as new tensile fractures nucleate and propagate subject to mechanical effects, numerical complexities associated with the introduction of new fracture control volumes are largely negated. In addition, the ability to discretize the fractures and surrounding rock volume independently provides the freedom to choose an acceptable level of discretization for each domain separately. Three numerical examples were performed to demonstrate the utility of the embedded fracture model for application to problems involving fluid flow, mechanical deformation, and rock failure. The results of the numerical examples confirm that the embedded fracture model was able to capture accurately the complex and nonlinear evolution of reservoir permeability as new fractures propagate through the reservoir and as fractures fail in shear.     

Capturing the complete stress–strain behaviour of jointed rock using a numerical approach

This paper presents the results of a series of numerical experiments using the synthetic rock mass (SRM) approach to quantify the behaviour of jointed rock masses. Field data from a massive sulphide rock mass, at the Brunswick mine, were used to develop a discrete fracture network (DFN). The constructed DFN model was subsequently subjected to random sampling whereby 40 cubic samples, of height to width ratio of two, and of varying widths (0.05 to 10 m) were isolated. The discrete fracture samples were linked to 3D bonded particle models to generate representative SRM models for each sample size. This approach simulated the jointed rock mass as an assembly of fractures embedded into the rock matrix. The SRM samples were submitted to uniaxial loading, and the complete stress–strain behaviour of each specimen was recorded. This approach provided a way to determine the complex constitutive behaviour of large‐scale rock mass samples. This is often difficult or not possible to achieve in the laboratory. The numerical experiments suggested that higher post‐peak modulus values were obtained for smaller samples and lower values for larger sample sizes. Furthermore, the observed deviation of the recorded post‐peak modulus values decreased with sample size. The ratio of residual strength of rock mass samples per uniaxial compressive strength intact increases moderately with sample size. Consequently, for the investigated massive sulphide rock mass, the pre‐peak and post‐peak representative elemental volume size was found to be the same (7 × 7 × 14 m). Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

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.

     

非常低延展性裂隙岩体REV存在性研究

裂隙岩体渗流问题一直是水文地质与工程地质的前沿课题之一,等效多孔介质模型是裂隙岩体渗流研究和工程应用的主要方法,此方法的重点在于裂隙岩体渗流模型即典型单元体(REV)的确定。本文采用作者自行开发的离散元软件FractureToKarst,根据国际岩石力学协会对岩石的分类,讨论了非常低延展性裂隙岩体典型单元体(REV)的存在性及大小,并给出确定其存在性的一般方法。通过研究可知,非常低延展性裂隙岩体REV存在的条件就是裂隙平均间距在0.2m以下,即极密间距、很密间距和密间距裂隙岩体存在REV。裂隙平均间距大于0.2 m,即中等间距、宽间距、很宽间距和极宽间距裂隙岩体不存在REV。     

岩体裂隙网络非线性渗流特性研究

基于人工交叉裂隙模型,通过室内透水试验,利用电荷耦合元件（CCD）照相机可视化技术,对流体在裂隙交叉点内的非线性流动特性进行研究。建立两种离散裂隙网络（DFN）模型,考虑两种边界条件,改变模型进口和出口之间的压力,直接求解Navier-Stokes（简称N-S）方程,对DFN的非线性渗流特性进行研究。结果表明,室内试验可以观测到与出口3相连的裂隙单元内发生了明显的非线性流动,且通过模型的流量Q和模型两端的压力P具有非线性关系。数值计算结果也表明,在水力梯度J较大时（比如J > 0.1）,通过DFN的Q和P具有非线性关系,而当J较小时（比如J < 10-4）,Q与P线性相关;根据文中的算例,建议利用局部立方定律求解DFN内每条裂隙的渗流特性的临界条件为J ≤10-4;裂隙表面粗糙会造成通过DFN渗流量的降低,但对相对流量误差的影响可忽略不计。     

A Connectivity Index for Discrete Fracture Networks

Connectivity is an important measure for assessing flow transport in rock, especially through fractures. In this paper, rock fracture systems are modelled by a discrete fracture model simulated by a marked point process. A connectivity index is then introduced to quantify the connectivity between any two points in space. Monte Carlo simulation is used to evaluate the connectivity index for stationary cases and relationships between the connectivity index and the parameters of the discrete fracture model are analysed. The average number of intersections per fracture, X_f, and the fracture intensity, P12 (P32), are calculated and the relationships between these parameters and the connectivity index are investigated, concluding that X_f is the more suitable parameter for the classification of rock mass flow properties. The relationships between the connectivity index and the percolation state of the fractured medium are also discussed. An edge correction is briefly discussed and a practical example is used to demonstrate the method of computing the connectivity index.     

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.     

Numerical analysis of the effect of natural microcracks on the supercritical CO2 fracturing crack network of shale rock based on bonded particle models

The problem of predicting the geometric structure of induced fractures is highly complex and significant in the fracturing stimulation of rock reservoirs. In the traditional continuous fracturing models, the mechanical properties of reservoir rock are input as macroscopic quantities. These models neglect the microcracks and discontinuous characteristics of rock, which are important factors influencing the geometric structure of the induced fractures. In this paper, we simulate supercritical CO₂ fracturing based on the bonded particle model to investigate the effect of original natural microcracks on the induced‐fracture network distribution. The microcracks are simulated explicitly as broken bonds that form and coalesce into macroscopic fractures in the supercritical CO₂ fracturing process. A calculation method for the distribution uniformity index (DUI) is proposed. The influence of the total number and DUI of initial microcracks on the mechanical properties of the rock sample is studied. The DUI of the induced fractures of supercritical CO₂ fracturing and hydraulic fracturing for different DUIs of initial microcracks are compared, holding other conditions constant. The sensitivity of the DUI of the induced fractures to that of initial natural microcracks under different horizontal stress ratios is also probed. The numerical results indicate that the distribution of induced fractures of supercritical CO₂ fracturing is more uniform than that of common hydraulic fracturing when the horizontal stress ratio is small.     

Damage characterization during laboratory strength testing: A 3D-finite-discrete element approach

A combined finite-discrete element approach is used to simulate the complete 3D fracture process during conventional laboratory testing, including Brazilian indirect tension and uniaxial and biaxial compression. A typical granite rock type (based on the Lac du Bonnet granite) was simulated to investigate the fracture pattern and mechanical strength of brittle rock in the laboratory. Damage intensity parameters (D₂₁ and D₃₂) are introduced and utilized to characterize the induced damage in the models. These parameters provide an improved representation of the cumulative associated damage and facilitate a quantitative characterization of crack intensity during testing. The numerical simulations included both 3D and 2D models, and show that there is a good agreement between the strength response derived from simulations both in 3D and 2D and the considered rock material. A good correlation also exists between the fracture pattern in 3D and the equivalent 2D models. The influence of confinement on the biaxial strength and the associated damage in compression is investigated. While axial splitting is the dominant failure mode at low confinement, finite-discrete element simulations show that a shear failure mode tends to dominate as the confinement increases. The dependency of dilation upon the confining pressure is also demonstrated, the dilation angle decreasing with increased confinement.     

评价页岩压裂形成缝网能力的新方法

页岩储层的“体积压裂”,使美国页岩气产业取得巨大成功,有效评价压裂裂缝网络形成的难易程度,是压裂开采的首要目标,目前国内外尚未发现有效的评价方法,为此开发了一种新的测试方法。针对10种岩芯,测试岩石力学参数,并对比分析常用的3种页岩脆性评价方法。采用压后裂隙结构面迹长分布的分维值和面密度对裂缝进行定量表征,并对压后崩落碎块进行对比分析。通过实验认为,杨氏模量和泊松比判别法与塑性系数判别法用于评价岩石脆性,精确度更高;脆性岩石通常表现为高杨氏模量或（和）低泊松比的特征,与单轴抗压强度、抗张强度和压入硬度没有对应关系;压裂裂缝的分布具有统计意义上的分形特征,分维可用于定量评价压后裂缝网络复杂度;硬度越高,压后裂缝密度越小;脆性越强,压后裂缝密度越大。新方法是岩石脆性、硬度和天然裂缝系统（和沉积层理）特征的综合体现,用于评价页岩压后形成缝网的能力,不仅直观可靠,而且简单有效,有利于现场推广应用,对于今后页岩气或致密砂岩气开发的理论研究和现场应用具有一定的指导意义。     

Numerical investigation of hydraulic fracture network propagation in naturally fractured shale formations

Hydraulic fracture network (HFN) propagation in naturally fractured shale formations is investigated numerically using a 3D complex fracturing model based on the discrete element method. To account for the plastic deformation behavior of shales, the Drucker–Prager plasticity model is incorporated into the fracturing model. Parametric studies are then conducted for different Young's moduli, horizontal differential stresses, natural fracture (NF) properties, injection rates, and number and spacing of perforation clusters. Numerical results show that horizontal differential stress primarily determines the generation of a complex HFN. The plastic deformation of shale can reduce the stimulated reservoir volume; this is more obvious with Young's modulus of less than 20 GPa. In addition, a higher injection rate could largely increase the fracture complexity index (FCI). Moreover, increasing perforation cluster numbers per fracturing stage is beneficial for increasing the FCI, but it also increases the potential merging of neighboring fractures, which may lead to non-uniform development of HFN in far-wellbore regions. To achieve uniform development of HFN within a fracturing stage, the distribution of NFs should be fully considered. The results presented here may provide improved understanding of HFN generation and are favorable for optimizing fracturing treatment designs for shale formations.     

XFEM-based cohesive zone approach for modeling near-wellbore hydraulic fracture complexity

Near-wellbore fracture tortuosity has important impacts on the productivity of fractured oil and gas wells and the injectivity of CO₂ or solids disposal injectors. Previous models for simulating near-wellbore fracture tortuosity usually assume fracture growth in linear-elastic media, without considering the effects of porous features of the rock. In this paper, a 2D fully coupled model is developed to simulate near-wellbore fracturing using the XFEM-based cohesive segment method. The model takes into account a variety of crucial physical aspects, including fracture extension and turning, fluid flow in the fracture, fluid leak-off through wellbore wall and fracture surfaces, pore fluid flow, and rock deformation. The proposed model was verified against two sets of published experimental results. Numerical examples were carried out to investigate the effects of various parameters on near-wellbore fracture trajectory, injection pressure, and fracture width. Results show that near-wellbore fracture behaviors are not only dependent on rock elastic properties and field stresses, but also greatly influenced by porous properties of the rock, such as permeability and leak-off coefficient. Some field implications were provided based on the simulation results. By overcoming some limitations of the previous models, the proposed model predicts more realistic fracture evolution in the near-wellbore region and provides an attractive tool for design and evaluation of many field operations, for which near-wellbore fracture behaviors play an important role on their successes.

     

Flow and transport predictions during multi-borehole tests in fractured chalk using discrete fracture network models

Multi-borehole pumping and tracer tests on the 10 to 100-m scale were conducted in a fractured chalk aquitard in the Negev Desert, Israel. Outcrop and core fracture surveys, as well as slug tests in packed-off intervals, were carried out at this site to obtain the parameters needed for construction of a stochastic discrete fracture network (DFN). Calibration of stochastic DFNs directly to the multiple borehole test data was inadequate. Instead, two equivalent deterministic DFN flow models were used: the vertical-fractures (VF) model, consisting of only vertical fractures, and the fractures’ intersections (INT) model, consisting of vertical and horizontal fractures with enhanced transmissivity at their intersections. Both models were calibrated against the multi-borehole response of one pumping test and their predictions were tested against three other independent pumping tests. The average accuracies of all transient drawdown predictions of the VF and INT models were 65 and 66%, respectively. In contrast to this equality in average drawdown predictions of both models, the INT model predicted better important breakthrough curve features (e.g., first and peak arrival times), than the VF model. This result is in line with previously assumed channeled flow, derived from analytical analysis of these pumping and tracer tests. Ronit Nativ, deceased, may her memory be blessed.     

洞周岩体内裂纹水力劈裂临界水压分析计算

水力劈裂是深埋隧洞施工涌水重要因素之一,对其破坏机制研究是岩土工程界的热点课题。根据裂纹面的应力状态,从断裂力学角度将岩体的裂纹扩展分为拉剪复合扩展和压剪复合扩展。应用地下岩体复合失稳扩展判据分别推导出两种破坏模式的临界水压计算公式,并对其随裂纹在隧洞洞周围的位置、方位变化规律进行分析。分析结果表明,在发生拉剪复合扩展情况下,裂纹方向与最大地应力方向平行时,最易失稳扩展;在压剪复合扩展情况下,裂纹与最大地应力夹角约介于30°～75°之间时,最易失稳扩展。