共查询到19条相似文献，搜索用时 328 毫秒
1.
基于FDEM-Flow研究地应力对水力压裂的影响 总被引：2，自引：0，他引：2
《岩土力学》2016,(1):237-246
用提出的FDEM-Flow(考虑流固耦合的离散元-有限元耦合方法)方法作为工具,研究了地应力对水力压裂的影响。通过一个注水圆孔的算例,研究不同地应力状态对压裂裂隙的走向和形态的影响。研究结果表明,起裂压力的大小和压裂裂隙的方位均与地应力有密切关系。在竖向地应力v?保持不变的情况下,且侧压力系数?1时,随着?的增大,起裂压力和失稳压力均减小;?1时且取值较小时,主要产生水平向的裂隙,并有斜向裂隙产生;?1时且取值较大时,裂隙严格按照最大主压应力的方向扩展,不再出现斜向裂隙;?1时,主要产生竖向的裂隙;?=1时,水平地应力和竖向地应力相等,裂隙的扩展不存在优势方向。不同侧压力系数条件下,裂隙的扩展方向与最大主应力的方向一致,水力压裂裂隙的起裂和扩展主要由最大主拉应力控制,裂隙在拉应力集中的区域起裂。这些结果与已有的试验及理论认识是相符的,进一步验证了FDEM-Flow方法用于模拟水力压裂问题的有效性。 相似文献
2.
3.
将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。 相似文献
4.
将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。 相似文献
5.
6.
水力压裂技术是油气藏尤其是页岩气开发中的核心技术,利用数值模拟方法进行压裂优化和产能预测又是水力压裂成功的关键。本文首先介绍了水力压裂技术的发展历程。然后从计算模型(二维模型、拟三维模型和全三维模型)和数值模拟方法(基于连续介质和基于非连续介质)两方面对油气藏开发领域的水力压裂计算模拟技术进行较全面的总结。最后,从以下3个方面指出现今研究的不足并提出了进一步的研究建议:(1)全三维模型的完善-全三维模型应当与真实的工程参数和监测数据结合,用于校正模型本身,而校正后的全三维模型又可预测和优化新的现场水力压裂作业; (2)数值模拟方法的选用-已有的水力压裂数值模拟方法种类繁多,需要针对各种方法的适用范围、计算效率和模拟效果等,进行全面的比较和优化; (3)页岩储层中天然裂缝网络的数值模拟-天然裂隙网络加剧了页岩储层力学性质的各向异性,同时水力裂缝沟通天然裂缝活化扩展是有利于储层的增渗增产,对压裂缝网的形态、尺寸和连通率等起着至关重要的作用。因此,数值计算过程中综合考虑页岩储层中天然裂缝与水力裂缝的相互作用,将是未来水力压裂模拟的热点。 相似文献
7.
水力压裂作为煤层强化增透技术的一种,其应力演化特征及裂隙形态与扩展范围的判断尤为重要。采用离散元数值方法,以导向压裂为背景,建立水力压裂流固耦合模型;通过应力路径、裂纹热点图等手段,探究水力压裂过程中压裂排量、泊松比、天然裂隙密度对应力演化和裂隙演化的影响及其细观规律。结果表明:不同压裂排量下的应力演化方向及最终应力路径曲线形状有着明显的不同,低排量下裂隙附近的应力比值逐渐增大,而在高排量下先增大后减小;煤层泊松比越大,平均压裂半径越低,但对起裂时间及裂隙的扩展形态影响不明显;天然裂隙的发育情况对水力裂隙的扩展起着关键性作用,高裂隙发育煤层水力裂隙扩展的方向性无法预测,应力演化方向会出现反转现象;压裂过程中不同区域的应力演化特征能够反映出裂隙的扩展状态,现场可通过监测压裂区域附近应力变化,判断水力压裂缝网的扩展范围。 相似文献
8.
陶可友 《地质灾害与环境保护》2018,(1)
运用岩石破裂过程分析软件RFPA2D,通过设置不同尺寸、剪切速率,探究裂隙形成、声发射能量以及位移之间的关系。数值模拟结果表明非均质岩石剪切破坏过程裂纹的扩展贯穿模式及声发射能量的对应关系。结果表明:(1)岩石剪切破坏的尺寸效应显著,剪切速率影响较小。(2)大裂隙的生成伴随较强的声发射能量,引起Y方向位移波动明显。(3)岩石剪切破坏剪裂纹是由大量细小张裂隙发育演变而来,对破坏模式起主导作用。 相似文献
9.
对基坑开挖影响范围内土体的应力路径进行平面应变试验离散元数值模拟,以研究结构性与卸荷形式对坑周土体宏微观力学特性的影响。首先,将一个描述土颗粒间胶结效应的简单三维胶结接触模型植入三维离散元软件PFC3D;其次,对初始K0固结状态的重塑土、结构性土试样分别进行常规三轴以及平面应变条件下4种不同卸荷应力路径的离散元模拟;最后,对经历不同卸荷形式的坑底土体单元进行再加荷模拟。模拟结果表明,在卸荷过程中,被动区土体峰值强度以及破坏时竖向应变随卸荷比增大而增大,且其强度小于主动区土体强度;在卸荷、再加荷过程中,被动区土体峰值强度随卸荷比增大而增大,但均小于不卸荷而直接加荷条件下的峰值强度;由于结构性的存在,土体由应变硬化向应变软化过渡,且强度增长;结构性与卸荷形式显著影响土体体积改变。在微观尺度,增大卸荷比或结构性均会增大垂直大主应力方向的平面上的法向接触力,进而提高其强度。 相似文献
10.
地壳中岩石的变形模式受构造应力、流体压力和上覆岩层重力共同作用的影响。岩石组成和构造应力的大小、方向决定着岩石的变形过程,同时岩石的破裂还受先存断裂构造的影响。流体压力增大,岩石可以发生水力破裂,而引起水力引张破裂的条件是σ1-σ3<4T和Pf=σ3+T。随着深度的增加,受地温梯度的影响,岩石由脆性变形向韧性变形转变。在无流体超压影响的情况下,脆韧性转换的温度在300~450℃之间,大约在地壳15km处。当流体压力和应变速率增大时,韧性条件下的岩石变形行为由韧性向脆性变化,脆韧性转变的深度随之增大。从构造角度探讨热液成矿作用,热液矿床形成的深度与流体压力、应变速率、裂隙的发育、介质的渗透率、温度变化等相关。岩石断裂的类型和方向影响岩石的渗透率,提供流体运移的通道和聚集场所,控制矿床形成的深度、位置和矿体产状。 相似文献
11.
页岩储层天然裂缝、水平层理发育,水力压裂过程中可能形成复杂的体积裂缝。针对页岩储层体积裂缝扩展问题,基于流-固耦合基本方程和损伤力学原理,建立了页岩储层水力压裂体积裂缝扩展的三维有限元模型。将数值模型的模拟结果与页岩储层裂缝扩展室内试验结果进行对比,二者吻合较好,从而证明了数值模型的可靠性。通过一系列数值模拟发现：（1）水力压裂过程中水平层理可能张开,形成水平缝,水平与垂直缝相互交错,形成复杂的体积裂缝网络;（2）水平主应力差增大,体积裂缝的分布长度（水平最大主应力方向压裂裂缝的展布距离）增加、分布宽度（水平最小主应力方向压裂裂缝的展布距离）减小,体积裂缝的长宽比增加;（3）压裂施工排量增大,体积裂缝的分布长度减小、宽度增加,压裂裂缝的长宽比降低;（4）天然裂缝的残余抗张强度增大,体积裂缝分布宽度减小、分布长度增加,体积裂缝的长宽比增加。研究成果可以为国内的页岩气的压裂设计和施工提供一定的参考和借鉴。 相似文献
12.
水力压裂可显著提高页岩气等致密储层岩体的渗透性以增加油气产量,然而受多因素影响,水力压裂形成缝网结构的机理和压裂优化设计一直是研究的焦点和难点。本研究基于渗流-应力-破坏耦合计算模拟方法,对不同水力加载条件下的非均质储层水力压裂过程进行了模拟和对比研究。研究结果表明:水力压裂过程中起始注水压力和增量大小对水力压裂缝网扩展和改造区域形态有着显著的影响。当起始注水压力小于等于模型材料体抗拉强度,并缓慢增压致裂时,压裂过程可近似视为稳态应力-破坏-渗流耦合作用过程的不同阶段,这种情况下仅在压裂井孔周围形成两组对称式的伞状水力裂缝带。当对模型体施加高于模型材料体破裂压力的注水压力时,相当于对压裂孔快速施加高动水压力,水力裂缝沿压裂孔全方位迅速萌生并快速扩展,当注水压力值高于破裂压力一定幅值时,压裂改造可形成围绕压裂井全方位的放射状裂缝网络,使压裂储层得以最大范围改造。在拟静力和拟动力两种加载条件下,不同水岩相互作用机理是造成不同水力加载条件出现不同缝网结构的力学机制,而对于实际的页岩气储层改造,压裂产生围绕压裂井全方位放射状的缝网结构则是一种最优的体积压裂改造。 相似文献
13.
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. 相似文献
14.
金湖凹陷天96断块阜二段储层裂缝多参数定量预测 总被引：1，自引：1，他引：0
基于应变能、表面能理论, 采用格里菲斯岩石破裂准则, 建立了金湖凹陷天96断块造缝期古地应力、现今地应力与裂缝参数之间的关系; 基于构造分析确定了造缝时期古地应力的类型和方向, 用声速法和水力压裂法确定了现今地应力的方向和大小, 将古今地应力场的数值模拟结果导入裂缝参数计算模型中, 模拟计算了天96断块阜二段裂缝的密度、开度、孔渗参数, 并结合古地应力场分析以及岩心裂缝统计数据对裂缝预测结果进行验证; 最后, 综合分析了影响天96断块裂缝发育的因素。 相似文献
15.
页岩储层孔隙度和渗透率极低,天然裂缝和水平层理发育,常规压裂增产措施无法满足页岩气的开发要求,水平井多段分簇压裂是页岩气开发的关键技术之一,该技术能够大幅度提升压裂改造的体积、产气量和最终采收率。为确定页岩储层水平井多段分簇射孔压裂的起裂点和起裂压力,采用有限元方法建立了水平井套管完井（考虑水泥环和套管的存在）多段分簇射孔的全三维起裂模型。数值模型的起裂压力与室内试验结果吻合较好,证明了数值模型的准确性和可靠性。利用数值模型研究了页岩水平井多段分簇射孔压裂的起裂点和起裂压力的影响因素,研究发现：射孔孔眼附近无天然裂缝或水平层理影响,起裂点发生在射孔簇孔眼的根部;射孔簇间距越小,中间射孔簇的干扰越大,可能造成中间的射孔簇无法起裂;射孔密度和孔眼长度增大,起裂压力降低;天然裂缝的存在,在某些情况能够降低起裂压力且改变起裂位置,主要与天然裂缝的分布方位及水平主应力差有关;水平层理可能会降低起裂压力,但与垂向主应力与水平最小主应力的差值有关。获得的起裂压力变化规律,可作为进一步研究水平井多段分簇射孔条件下的裂缝扩展规律的基础,可以为压裂设计和施工的射孔参数确定及优化给出具体建议。 相似文献
16.
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.
相似文献17.
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. 相似文献
18.
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. 相似文献
19.
Fully coupled simulation of a hydraulic fracture interacting with natural fractures with a hybrid discrete‐continuum method 
下载免费PDF全文
下载免费PDF全文 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. 相似文献