共查询到20条相似文献，搜索用时 147 毫秒
1.
Numerical Simulation of 3D Hydraulic Fracturing Based on an Improved Flow-Stress-Damage Model and a Parallel FEM Technique 总被引：4，自引：4，他引：0
L. C. Li C. A. Tang G. Li S. Y. Wang Z. Z. Liang Y. B. Zhang 《Rock Mechanics and Rock Engineering》2012,45(5):801-818
The failure mechanism of hydraulic fractures in heterogeneous geological materials is an important topic in mining and petroleum engineering. A three-dimensional (3D) finite element model that considers the coupled effects of seepage, damage, and the stress field is introduced. This model is based on a previously developed two-dimensional (2D) version of the model (RFPA2D-Rock Failure Process Analysis). The RFPA3D-Parallel model is developed using a parallel finite element method with a message-passing interface library. The constitutive law of this model considers strength and stiffness degradation, stress-dependent permeability for the pre-peak stage, and deformation-dependent permeability for the post-peak stage. Using this model, 3D modelling of progressive failure and associated fluid flow in rock are conducted and used to investigate the hydro-mechanical response of rock samples at laboratory scale. The responses investigated are the axial stress–axial strain together with permeability evolution and fracture patterns at various stages of loading. Then, the hydraulic fracturing process inside a rock specimen is numerically simulated. Three coupled processes are considered: (1) mechanical deformation of the solid medium induced by the fluid pressure acting on the fracture surfaces and the rock skeleton, (2) fluid flow within the fracture, and (3) propagation of the fracture. The numerically simulated results show that the fractures from a vertical wellbore propagate in the maximum principal stress direction without branching, turning, and twisting in the case of a large difference in the magnitude of the far-field stresses. Otherwise, the fracture initiates in a non-preferred direction and plane then turns and twists during propagation to become aligned with the preferred direction and plane. This pattern of fracturing is common when the rock formation contains multiple layers with different material properties. In addition, local heterogeneity of the rock matrix and macro-scale stress fluctuations due to the variability of material properties can cause the branching, turning, and twisting of fractures. 相似文献
2.
S. Y. Wang S. W. Sloan S. G. Fityus D. V. Griffiths C. A. Tang 《Rock Mechanics and Rock Engineering》2013,46(5):1165-1182
Rock is a heterogeneous geological material. When rock is subjected to internal hydraulic pressure and external mechanical loading, the fluid flow properties will be altered by closing, opening, or other interaction of pre-existing weaknesses or by induced new fractures. Meanwhile, the pore pressure can influence the fracture behavior on both a local and global scale. A finite element model that can consider the coupled effects of seepage, damage and stress field in heterogeneous rock is described. First, two series of numerical tests in relatively homogeneous and heterogeneous rocks were performed to investigate the influence of pore pressure magnitude and gradient on initiation and propagation of tensile fractures. Second, to examine the initiation of hydraulic fractures and their subsequent propagation, a series of numerical simulations of the behavior of two injection holes inside a saturated rock mass are carried out. The rock is subjected to different initial in situ stress ratios and to an internal injection (pore) pressure at the two injection holes. Numerically, simulated results indicate that tensile fracture is strongly influenced by both pore pressure magnitude and pore pressure gradient. In addition, the heterogeneity of rock, the initial in situ stress ratio (K), the distance between two injection holes, and the difference of the pore pressure in the two injection holes all play important roles in the initiation and propagation of hydraulic fractures. At relatively close spacing and when the two principal stresses are of similar magnitude, the proximity of adjacent injection holes can cause fracturing to occur in a direction perpendicular to the maximum principal stress. 相似文献
3.
水平井中多条水力裂缝间的应力干扰行为,造成了压裂液排量的非均匀分配,影响了水力裂缝的几何形态。采用边界元法研究岩体在压裂液作用下的变形程度,以幂律流体泊肃叶平板流动方程来研究水力裂缝内部的压裂液流场,考虑了多条裂缝间应力干扰和压裂液流量分配,建立了流-固耦合的水平井多条水力裂缝同步扩展模型。模型可模拟水平井多条水力裂缝几何形态、应力干扰情况和压裂液排量的分配情况,可解释水力裂缝之间的竞争机制。多条裂缝同步扩展时,压裂液排量并非均等地分配到各个裂缝之中,进入到内部裂缝的压裂液流量最小,内部裂缝宽度最小;内部的水力裂缝增长一定长度后,停止增长,并且在应力干扰下逐渐闭合。 相似文献
4.
Magnus Wangen 《Computational Geosciences》2013,17(4):647-659
A procedure based on the finite element method is suggested for modeling of 3D hydraulic fracturing in the subsurface. The proposed formulation partitions the stress field into the initial stress state and an additional stress state caused by pressure buildup. The additional stress is obtained as a solution of the Biot equations for coupled fluid flow and deformations in the rock. The fluid flow in the fracture is represented on a regular finite element grid by means of “fracture” porosity, which is the volume fraction of the fracture. The use of the fracture porosity allows for a uniform finite element formulation for the fracture and the rock, both with respect to fluid pressure and displacement. It is demonstrated how the fracture aperture is obtained from the displacement field. The model has a fracture criterion by means of a strain limit in each element. It is shown how this criterion scales with the element size. Fracturing becomes an intermittent process, and each event is followed by a pressure drop. A procedure is suggested for the computation of the pressure drop. Two examples of hydraulic fracturing are given, when the pressure buildup is from fluid injection by a well. One case is of a homogeneous rock, and the other case is an inhomogeneous rock. The fracture geometry, well pressure, new fracture area, and elastic energy released in each event are computed. The fracture geometry is three orthogonal fracture planes in the homogeneous case, and it is a branched fracture in the inhomogeneous case. 相似文献
5.
A two‐dimensional coupled hydromechanical discontinuum model for simulating rock hydraulic fracturing 
下载免费PDF全文
下载免费PDF全文 Yu‐Yong Jiao Huan‐Qiang Zhang Xiu‐Li Zhang Hai‐Bo Li Qing‐Hui Jiang 《国际地质力学数值与分析法杂志》2015,39(5):457-481
Within the framework of our discontinuous deformation analysis for rock failure algorithm, this paper presents a two‐dimensional coupled hydromechanical discontinuum model for simulating the rock hydraulic fracturing process. In the proposed approach, based on the generated joint network, the calculation of fluid mechanics is performed first to obtain the seepage pressure near the tips of existing cracks, and then the fluid pressure is treated as linearly distributed loads on corresponding block boundaries. The contribution of the hydraulic pressure to the initiation/propagation of the cracks is considered by adding the components of these blocks into the force matrix of the global equilibrium equation. Finally, failure criteria are applied at the crack tips to determine the occurrence of cracking events. Several verification examples are simulated, and the results show that this newly proposed numerical model can simulate the hydraulic fracturing process correctly and effectively. Although the numerical and experimental verifications focus on one unique preexisting crack, because of the capability of discontinuous deformation analysis in simulating block‐like structures, the proposed approach is capable of modeling rock hydraulic fracturing processes. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
6.
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. 相似文献
7.
Yao Yao 《Rock Mechanics and Rock Engineering》2012,45(3):375-387
Hydraulic fracturing technology is being widely used within the oil and gas industry for both waste injection and unconventional
gas production wells. It is essential to predict the behavior of hydraulic fractures accurately based on understanding the
fundamental mechanism(s). The prevailing approach for hydraulic fracture modeling continues to rely on computational methods
based on Linear Elastic Fracture Mechanics (LEFM). Generally, these methods give reasonable predictions for hard rock hydraulic
fracture processes, but still have inherent limitations, especially when fluid injection is performed in soft rock/sand or
other non-conventional formations. These methods typically give very conservative predictions on fracture geometry and inaccurate
estimation of required fracture pressure. One of the reasons the LEFM-based methods fail to give accurate predictions for
these materials is that the fracture process zone ahead of the crack tip and softening effect should not be neglected in ductile
rock fracture analysis. A 3D pore pressure cohesive zone model has been developed and applied to predict hydraulic fracturing
under fluid injection. The cohesive zone method is a numerical tool developed to model crack initiation and growth in quasi-brittle
materials considering the material softening effect. The pore pressure cohesive zone model has been applied to investigate
the hydraulic fracture with different rock properties. The hydraulic fracture predictions of a three-layer water injection
case have been compared using the pore pressure cohesive zone model with revised parameters, LEFM-based pseudo 3D model, a
Perkins-Kern–Nordgren (PKN) model, and an analytical solution. Based on the size of the fracture process zone and its effect
on crack extension in ductile rock, the fundamental mechanical difference of LEFM and cohesive fracture mechanics-based methods
is discussed. An effective fracture toughness method has been proposed to consider the fracture process zone effect on the
ductile rock fracture. 相似文献
8.
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.
相似文献9.
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. 相似文献
10.
A three-phase hydro-mechanical model for hydraulic fracturing is proposed. Three phases include: porous solid, fracturing fluid and host fluid. Discontinuity is handled using extended finite element method (XFEM) while cohesive crack model is used as fracturing criterion. Flow through fracture is defined as one-dimensional laminar flow, and flow through porous medium (host reservoir) is defined as two-dimensional Darcy flow. Coupling between two fluids in each space, fracture and pore, is captured through capillary pressure–saturation relationship, while the identical fluids in fracture and pore are coupled through a so-called leak-off mass transfer term. Coupling between fluids and deformation is captured through compatibility of volumetric strain of fluids within fracture and pore, and volumetric strain of the matrix. Spatial and temporal discretisation is achieved using the standard Galerkin method and the finite difference technique, respectively. The model is verified against analytical solutions available from literature. The leaking of fracturing fluid into the medium and suction of porous fluid into the fracture around the tip, are investigated. Sensitivity analyses are carried out for cases with slow and fast injection rates. It is shown that the results by single-phase flow may underestimate the leak-off. 相似文献
11.
The hydraulic fracturing technique has been widely applied in many fields, such as the enhanced geothermal systems (EGS), the improvement of injection rates for geologic sequestration of CO2, and for the stimulations of oil and gas reservoirs. The key points for the success of hydraulic fracturing operations in unconventional resources are to accurately estimate the redistribution of pore pressure and stresses around the induced fracture and predict the reactivations of preexisting natural fractures. The pore pressure and stress regime around hydraulic fracture are affected by poroelastic and thermoelastic phenomena as well as by fracture opening compression. In this work, a comprehensive semi-analytical model is used to estimate the stress and pore pressure distribution around an injection-induced fracture from a single well in an infinite reservoir. The model allows the leak-off distribution in the formation to be three-dimensional with the pressure transient moving ellipsoidically outward into the reservoir from the fracture surface. The pore pressure and the stress changes in three dimensions at any point around the fracture caused by poroelasticity, thermoelasticity, and fracture compression are investigated. With Mohr-Coulomb failure criterion, we calculate the natural fracture reactivations in the reservoir. Then, two case studies of constant water injection into a hydraulic fracture are presented. This work is of interest in the interpretation of microseismicity in hydraulic fracturing and in the estimation of the fracture spacing for hydraulic fracturing operations. In addition, the results from this study can be very helpful for the selection of stimulated wells and further design of the refracturing operations. 相似文献
12.
Water outbursts from the floor during underground mining, and those from the surrounding rock mass of tunnels, involve the basic principle of hydraulic fracturing. Based on the hydraulic-fracturing mechanism, considered to be dependent on the coupling between seepage and damage, it is deemed that the variation of the pore-fluid pressure coefficient must be taken into account during this coupled process, in order to correctly establish the crack propagation mechanism during hydraulic fracturing. The coupled seepage-damage model is validated using numerical simulations of hydraulic fracturing around one hole and three holes; the model may also enable scientific and reasonable explanation of the dominance of hydraulic gradient on the crack propagation path in permeable rock. Finally, the water outburst from the floor at a coal mining site in Hebei Province, China, is numerically simulated, and the coupled seepage and damage mechanism during the mining-induced rock failure is clarified. The numerical simulation implies that the seepage-damage is the main mechanism for controlling the water outburst. Therefore this mechanism should be considered in the numerical simulation to understand the essence of water outburst induced in mines. 相似文献
13.
This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids. 相似文献
14.
Zhang Liaoyuan Qu Zhanqing Zhai Mingyang Li Aishan Huang Bo Zhang Zilin 《Geotechnical and Geological Engineering》2021,39(2):839-854
Tight heterogeneous glutenite reservoir is typically not easy to form complex hydraulic fracture (HF) due to its poor physical properties, poor matrix seepage capacity, and small limit discharge radius and undeveloped natural fracture system. To improve the HF complexity and the stimulated reservoir volume (SRV), a novel stimulation technology called CO2 miscible fracturing has been introduced and its fracturing mechanism has been studied. The CO2 miscible fracturing modifies the in situ stress field by injecting low viscosity fluid to increase the HF complexity and SRV. Therefore, a series of numerical simulations based on a hydro-mechanical-damage model were carried out to study the effects of low viscosity fluid pre-injection on pore pressure, stress field, and fracturing effect in tight heterogeneous glutenite reservoirs. The results indicate that the low viscosity fluid injection can effectively increase the pore pressure around the wellbore and reduce the effective stress of the glutenite. The FCI and SRV increase with the increase of the pre-injection amount of the low viscosity fluid. The HF complexity and SRV can be improved by pre-injecting low viscosity fluid to transform the in situ stress field. The field application of this technology in a well of Shengli Oilfield showed that low-viscosity fluid pre-injection can effectively increase the width of the fractured zone, improve the SRV, and optimize the fracturing effect.
相似文献15.
超深储层地层起裂压力较高,水力压裂受现场泵注设备的限制严重,文中重点研究了大尺度水力压裂物模实验水泥样品尺寸（762 mm×762 mm×914 mm）在循环和常规两种泵注条件下的起裂扩展和声发射规律。实验结果显示,（1）相对于普通泵注,采用循环泵注方式进行水力压裂可以有效降低起裂压力,类似于单轴和三轴循环加载岩石力学行为,都是由于循环加载引起疲劳损伤;（2）对于螺旋射孔完井方式,水力压裂裂缝只从最薄弱射孔处起裂,一旦起裂后其他射孔孔眼很难再开启,水力压裂现场应合理选择分段距离和簇间距,实现储层改造效率最优化;（3）循环泵注水力压裂存在Kaiser效应（当加载应力到前次加载最高应力值时出现的声发射信息）。因孔隙流体扩散到岩石并导致孔隙压力的局部上升,破坏模式仍然可以由摩尔圆表示。研究成果对循环泵注条件的裂缝扩展规律研究以及发展新型压裂改造技术具有重要意义。 相似文献
16.
基于完全流-固耦合的弹塑性理论给出了水力压裂数值计算的弥散裂缝模型,其中材料的弹性部分采用线弹性本构关系,塑性部分采用摩尔-库仑破坏准则及强化准则。依据当前的有效应力状态修正渗透系数来模拟压裂液在裂缝中的流动。渗透系数的修改使用双曲正切函数,并采用平均有效应力作为水力裂缝的起裂判据。在ABAQUS软件中通过用户自定义程序添加了该模型。根据岩石的切面照片建立了含有硬包裹体分布的非均质岩石的有限元计算模型,模拟了中心点注水条件下的水力压裂传播过程,讨论了在常应力状态下非均质岩石中开裂区域、典型位置的应力路径变化和裂缝传播范围随时间变化的特点。进行了多种条件下含有硬包裹体分布的岩石材料的数值试验,得出了基岩材料的弹性模量、凝聚力和渗透系数以及注水速率对峰值注水压力、平均注水压力和裂缝开度的影响规律。 相似文献
17.
水力压裂是低渗油气藏的主要开发手段,传统数值模型所得到的基质-裂缝窜流量以及断裂参数精度不足.为此以流固耦合理论与断裂力学相结合的压裂模型为基础,模拟了水力裂缝扩展过程.在模型中分别引入离散裂缝模型和广义J积分计算基质-裂缝流量交换和断裂参数,并采用动态网格技术对裂缝尖端进行局部加密,以提高模拟的效率和精度.模型计算结果显示,影响水力压裂过程的主要参数中:基质渗透率和压裂液粘度主要影响水力裂缝的最终形态;岩石弹性模量影响裂缝宽度.对压裂车而言,最高工作压力一般都能够满足压裂增产需求,其最大输出功率和最大输出流量是限制压裂能力的主要因素. 相似文献
18.
Jiao Kaituo Han Dongxu Wang Daobing Chen Yujie Li Jingfa Gong Liang Bai Bofeng Yu Bo 《Computational Geosciences》2022,26(5):1167-1187
Thermal-hydro-mechanical (THM) coupled fracture propagation is common in underground engineering. Rock damage, as an inherent property of rock, significantly affects fracture propagation, but how it influences the THM coupled fracturing remains stubbornly unclear. A pore-scale THM coupling model is developed to study this problem, which combines the lattice Boltzmann method (LBM), the discrete element method (DEM), and rock damage development theory together for the first time. This model can more accurately calculate the exchanged THM information at the fluid-solid boundary and fluid conductivity dependent on fracture and rock damage. Based on the developed model, the synergistic effect of injected temperature difference (fluid temperature below rock temperature) and rock damage (characterized by the parameter “critical fracture energy”, abbreviated as “CFE”) on fracture propagation of shale are investigated particularly. It is found that: (1) the generation of branched cracks is closely related to the temperature response frontier, and the fracture process zone of single bond failure increases in higher CFE. (2) through the analysis of micro failure events, hydraulic fracturing is more pronounced in the low CFE, while thermal fracturing displays the opposite trend. The fluid conductivity of fractured rock increases with a higher injected temperature difference due to the more penetrated cracks and wider fracture aperture. However, this enhancement weakens when rock damage is significant. (3) in the multiple-layered rock with various CFEs, branched cracks propagating to adjacent layers are more difficult to form when the injection hole stays in the layer with significant rock damage than without rock damage.
相似文献19.
Flowback analysis recently has been considered as a potential tool for assessing induced fractures through corresponding pressure analyses and rate transient analysis. In this paper, we study fracture closure mechanisms during the flowback of fracturing fluid after hydraulic fracturing treatments. Although it is known that flowback can be significantly affected by the geometry of the fractures and closure stress, there has not been any effort to understand the problem from the geomechanical point of view; rather, available studies assume that a fracture closes uniformly with constant fracture compressibility. The coupled geomechanics and fluid flow model presented in this paper help to elucidate how geostresses may affect fracturing fluid recovery under different conditions. We perform a scaling analysis to formulate the occurrence of different fracture closure modes and then use numerical analyses to verify our scaling parameters. The factors governing the flowback process include the mechanical and petrophysical properties of the rock as well as preexisting discontinuities such as natural fractures. Three different closure modes for fracture networks are described and numerically verified. The occurrence of each mode can dramatically affect fracturing fluid recovery. The role of fluid leakoff into the formation, fractures conductivity, and drawdown strategy are examined for each fracture closure scenario. 相似文献
20.
随着扩展有限元理论的深入研究,利用扩展有限元方法模拟水力压裂具有了一定的可操作性。相比于常规有限元方法,XFEM方法具有计算结果精度高和计算量小的优点。但是,如何模拟射孔孔眼、如何模拟流体与岩石相互作用以及分析水力裂缝的扩展规律仍然是难题。以研究水力压裂裂缝扩展规律为目的,建立了岩石多孔介质应力平衡方程、流体渗流连续性方程和边界条件。通过有限元离散化方法对耦合方程矩阵进行处理。通过富集函数定义初始裂缝（射孔孔眼）,选择最大主应力及损伤变量D分别作为裂缝起裂和扩展判定准则,利用水平集方法模拟水力裂缝扩展过程。数值模拟结果显示：增加射孔方位角、压裂液排量和减小水平地应力差,起裂压力上升;黏度对起裂压力无明显影响。增加射孔方位角、压裂液排量、黏度和减小水平地应力差值有助于裂缝宽度的增加。增加水平地应力差值、压裂液排量和减小射孔方位角以及压裂液黏度有助于裂缝长度增加,反之亦然。基于ABAQUS的水力裂缝扩展有限元法可对不同井型和诸多储层物性参数及压裂施工参数进行分析,且裂缝形态逼真,裂缝面凹凸程度清晰,结果准确。此研究可作为一种简便有效研究水力压裂裂缝扩展规律的方法为油田水力压裂设计与施工提供参考与依据。 相似文献