共查询到20条相似文献，搜索用时 46 毫秒
1.
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. 相似文献
2.
Numerical modeling of hydraulic fracture propagation,closure and reopening using XFEM with application to in‐situ stress estimation 
下载免费PDF全文
下载免费PDF全文 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. 相似文献
3.
新疆油田某地区油藏的储隔层岩性组合复杂,呈现突出的薄互层产状特征,研究合、分压判断条件有利于提高压裂效率,增强储层动用程度与压后改造效果。水力裂缝在薄互层中的穿层与裂缝扩展行为受薄互层地质特征与压裂施工参数的影响。基于此,开展了薄互层物理模型压裂试验,研究界面胶结、岩层分布、岩层厚度、压裂液黏度与注液排量对薄互层中水力裂缝垂向扩展的影响分析。试验结果表明：薄互层的地层特征界面胶结与岩层分布是水力裂缝垂向扩展的主要控制因素,界面胶结强度对裂缝垂向扩展行为的影响强于岩层分布;由于弱胶结界面的存在,水力裂缝垂向扩展穿层时可发生方向偏转,抑制裂缝垂向扩展;提高压裂液黏度与注液排量有利于薄互层中水力裂缝的穿层垂向扩展。 相似文献
4.
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. 相似文献
5.
通过室内水力压裂物理模拟试验系统,对大尺寸原煤进行了水力压裂模拟试验,根据水力裂缝的空间展布形态分析了煤岩储层水力裂缝的延伸规律,揭示了网状裂缝的形成机制。结果表明：水力裂缝易在弱层理处分叉和转向,发育的层理和裂缝系统等结构面为压裂形成裂缝网络提供了前提条件。泵压曲线呈现出的频繁波动是煤岩内产生网状裂缝的一个显著特征。水力裂缝的起裂与延伸有4种基本模式,裂缝网络的形成多为这4种基本模式的组合。地应力差异系数和泵注排量对煤层水力裂缝形态有较大影响。较小的地应力差异系数更利于网状裂缝的形成;较高的压裂液排量易形成相对简单的裂缝形态,导致压裂改造效果较差。该试验方法和试验结果可为现场水力压裂参数设计和优化提供参考和依据。 相似文献
6.
砂砾岩水力压裂裂缝扩展规律的数值模拟分析 总被引：1，自引：0，他引：1
砂砾岩储层一般具有岩性和渗透性变化大、孔隙度低、连通性差、孔隙结构复杂和非均匀性严重等特点,因此,在水力压裂过程中,裂缝扩展形态难以控制,大规模改造难度大。针对国内某典型砂砾岩油藏特征,采用数值计算方法对砂砾岩压裂裂缝的扩展规律进行了研究,包括地应力场、砾石含量和粒径等对裂缝扩展形态及压裂压力的影响。研究表明,砾石的存在增加了压裂裂纹扩展的复杂性,裂纹主要有止裂、偏转、穿透和吸附4种表现模式,但主应力差严格控制着裂纹的走向,随着主应力差的增大,裂纹由总体绕砾扩展转变为总体穿砾扩展,失稳压力随着主应力差的增大而明显减小;砾石含量的多少体现了砂砾岩试样宏观的非均匀性,含量越高均匀性越差,随着砾石含量的提高,裂纹与砾石的相互作用占据主导地位,失稳压力随砾石含量的增加而增大;当砾石体积含量一定时,砾石粒径对压裂压力的影响主要取决于砾石排列的随机性,失稳压力随砾石粒径的增大而略有增大。 相似文献
7.
超临界CO2是一种介于气体和液体之间的特殊状态的CO2流体,具有低黏、高扩散性和零表面张力等独特的性质。利用超临界CO2作为压裂液,有助于裂缝的起裂和扩展,同时可避免储层伤害。通过研究超临界CO2射流破岩和压裂特性,分析得到了超临界CO2岩石致裂机制。研究结果表明,超临界CO2低黏等特性使其更容易进入岩石微孔和微缝之中,在岩石内部建立大小不一的流体压力系统,使岩石发生拉伸和剪切破坏;常规流体压裂起裂压力较高,裂缝一般为单条或多条平直裂缝,大多沿着同一方向贯穿强度较高的胶结颗粒,且裂缝断面光滑、平整;超临界CO2压裂起裂压力相比于常规流体压裂低,在岩石中形成的裂缝网络较为复杂,裂缝互相连通,一般沿着强度较低的胶结物开裂,较少贯穿胶结颗粒,裂缝断面较为粗糙。该研究结果可为超临界CO2压裂技术的实施提供理论支撑。 相似文献
8.
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. 相似文献
9.
Research investigations on three-dimensional (3-D) rectangular hydraulic fracture configurations with varying degrees of fluid lag are reported. This paper demonstrates that a 3-D fracture model coupled with fluid lag (a small region of reduced pressure) at the fracture tip can predict very large excess pressure measurements for hydraulic fracture processes. Predictions of fracture propagation based on critical stress intensity factors are extremely sensitive to the pressure profile at the tip of a propagating fracture. This strong sensitivity to the pressure profile at the tip of a hydraulic fracture is more strongly pronounced in 3-D models versus 2-D models because 3-D fractures are clamped at the top and bottom, and pressures in the 3-D fractures that are far removed from the fracture tip have little effect on the stress intensity factor at the fracture tip. This rationale for the excess pressure mechanism is in marked contrast to the crack tip process damage zone assumptions and attendant high rock fracture toughness value hypotheses advanced in the literature. A comparison with field data is presented to illustrate the proposed fracture fluid pressure sensitivity phenomenon. This paper does not attempt to calculate the length of the fluid lag region in a propagating fracture but instead attempts to show that the pressure profile at the tip of the propagating fracture plays a major role in fracture propagation, and this role is magnified in 3-D models. © 1997 by John Wiley & Sons, Ltd. Int. J. Numer. Anal. Meth. Geomech., vol. 21, 229–240 (1997). 相似文献
10.
Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock 下载免费PDF全文
下载免费PDF全文 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. 相似文献
11.
为了研究煤岩水力压裂的起裂压力和水力压裂裂缝扩展规律,采用型煤试样,利用自主研发的水力压裂实验系统,参照现场压裂施工制定了“施加三向应力-顶部注水”的煤岩水力压裂物理模拟实验方案并开展了水力压裂实验,分析了不同条件下泵注压力和水力压裂裂缝。实验结果表明:压裂液泵注排量越大,起裂压力越大。三向应力满足最大水平主应力σH > 垂向应力σv > 最小水平主应力σh,水力压裂裂缝沿着垂直于σh的方向扩展。σv和σh一定,随着σH的增大,煤岩起裂压力先增大后减小,水力压裂裂缝扩展路径越平直。当σH远大于σv和σh时,水力压裂裂缝扩展路径越复杂,分叉缝角度越大。研究结果可为煤岩水力压裂理论的完善提供一定的参考和借鉴。 相似文献
12.
13.
We present a stabilized extended finite element formulation to simulate the hydraulic fracturing process in an elasto‐plastic medium. The fracture propagation process is governed by a cohesive fracture model, where a trilinear traction‐separation law is used to describe normal contact, cohesion and strength softening on the fracture face. Fluid flow inside the fracture channel is governed by the lubrication equation, and the flow rate is related to the fluid pressure gradient by the ‘cubic’ law. Fluid leak off happens only in the normal direction and is assumed to be governed by the Carter's leak‐off model. We propose a ‘local’ U‐P (displacement‐pressure) formulation to discretize the fluid‐solid coupled system, where volume shape functions are used to interpolate the fluid pressure field on the fracture face. The ‘local’ U‐P approach is compatible with the extended finite element framework, and a separate mesh is not required to describe the fluid flow. The coupled system of equations is solved iteratively by the standard Newton‐Raphson method. We identify instability issues associated with the fluid flow inside the fracture channel, and use the polynomial pressure projection method to reduce the pressure oscillations resulting from the instability. Numerical examples demonstrate that the proposed framework is effective in modeling 3D hydraulic fracture propagation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
14.
Hydraulic fracturing (HF) treatment often involves particle migration and is applied for propping or plugging fractures. Particle migration behaviors, e.g., bridging, packing, and plugging, significantly affect the HF process. Hence, it is crucial to effectively simulate particle migration. In this study, a new numerical approach is developed based on a coupled element partition method (EPM). The EPM is used to model natural and hydraulic fractures, in which a fracture is allowed to propagate across an element, thereby avoiding remeshing in fracture simulations. To characterize the water flow process in a fracture, a fully hydromechanical coupled equation is adopted in the EPM. To model particle transportation in fractures with water flow, each particle is treated as a discrete element. The particles move in the fracture as a result of being dragged by fluid. Their movement, contact, and packing behaviors are simulated using the discrete element method. To reflect the plugging effect, an equivalent aperture approach is proposed. Using this method, the particle migration and its effect on water flow are well simulated. The simulation results show that this method can effectively reproduce particle bridging, plugging, and unblocking in a hydraulic fracture. Furthermore, it is demonstrated that particle plugging significantly affects water flow in a fracture and hence the propagation of hydraulic fracture. This method provides a simple and feasible approach for the simulation of particle migration in a hydraulic fracture. 相似文献
15.
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. 相似文献
16.
酸压改造后的酸蚀裂缝是储层流体的主要流动通道,研究流体在酸蚀裂缝中的流动规律是十分必要的。本文介绍了目前裂隙流动研究中的4种常用控制方程及各自的适用条件,确定了Navier-Stokes方程为酸蚀裂缝流场数值模拟的控制方程。为了分析酸蚀裂缝流场分布规律以及非达西效应对裂缝导流能力的影响,利用逆向工程技术对两种不同刻蚀形态的酸蚀岩样壁面进行了实体重构,并采用有限元数值模拟方法进行了不同流量下的酸蚀裂缝流动试验。结果表明：酸刻蚀形态和裂缝接触关系对流场分布的影响较大。隙宽分布平滑的均匀刻蚀裂缝,其流态稳定、曲折度低,但由于裂缝开度较窄且比表面较高,边界层效应明显,流体流动时产生的黏滞阻力较高;沟槽裂缝曲面粗糙,隙宽分布复杂,流体流动时流态不稳定,曲折度高,在高流量下会产生明显的涡流流动,增大惯性阻力。裂缝通道存在缩颈现象会引起流体绕流和多次加减速,产生额外的压力损耗。随着模拟流量的增大,压降与流量会逐渐偏离线性关系而呈现出非达西流动现象。酸蚀裂缝壁面越粗糙,产生非达西效应时的临界流量和临界雷诺数就越小,在相同流量下非达西效应就越强烈,导流能力下降速率更快。 相似文献
17.
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. 相似文献
18.
Influence of Heterogeneity of Mechanical Properties on Hydraulic Fracturing in Permeable Rocks 总被引：11，自引：4，他引：7
T. H. Yang L. G. Tham C. A. Tang Z. Z. Liang Y. Tsui 《Rock Mechanics and Rock Engineering》2004,37(4):251-275
Summary Numerical simulations of circular holes under internal hydraulic pressure are carried out to investigate the hydraulic fracture initiation, propagation and breakdown behavior in rocks. The hydraulic pressure increases at a constant rate. The heterogeneity of the rocks is taken into account in the study by varying the homogeneity index. In addition, the permeability is varied with the states of stress and fracture. The simulations are conducted by using a finite element code, F-RFPA2D, which couples the flow, stress and damage analyses. The simulation results suggest that the fracture initiation and propagation, the roughness of the fracture path and the breakdown pressure are influenced considerably by the heterogeneity of rocks. The hole diameter elongation and the stress field evolution around the fracture tip during the fracture propagation can also provide useful information for the interpretation of the hydraulic fracturing behaviour. 相似文献
19.
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. 相似文献
20.
为了探究应力阴影效应对交替压裂中压裂间距选取的影响,基于优化后的颗粒流离散元流固耦合计算模型,模拟并分析了双初始水力裂缝下因应力阴影效应产生的诱导应力的分布情况,并与理论解析解进行对比,证明了该数值方法的合理性。在此基础上,分析了应力阴影效应在不同各向异性地应力场及初始压裂间距条件下对新水力裂缝的起裂压力及扩展形态的影响,研究结果表明：初始各向异性应力场不改变裂缝周边的应力场,不影响新水力裂缝的起裂压力;随着初始压裂间距的减小,应力阴影效应增强,新水力裂缝的起裂压力逐渐增加。初始水力裂缝间距与初始各向异性应力场共同影响新水力裂缝的扩展形态,随着初始水力裂缝间距或初始水平地应力场差异系数的增大,应力阴影对新水力裂缝的扩展方向的影响逐步减弱;初始水力裂缝对新水力裂缝的扩展有一定的限制作用,在一定程度上不利于形成复杂的裂缝网络。根据以上分析结果,对交替压裂中压裂间距的优化进行了定性的探讨。 相似文献