含层理页岩气藏水力压裂裂纹扩展规律解析分析

页岩气蕴藏在页岩层中,页岩层的层理性构造使其水力压裂裂纹扩展与常规均质储层不同.为研究页岩储层水力压裂的裂纹扩展规律,基于复变函数保角变换,得出裂纹尖端应力集中解,考虑页岩非均质、强度各向异性特点,通过比较裂纹沿各方向扩展所需的裂缝尖端水压力,推导出水力压裂裂纹垂直于最小地应力方向稳定扩展过程中在斜交层理后的扩展判据.分别定义了水力压裂裂纹在层理处起裂和沿层理扩展的弱层和岩石基体临界强度比,根据两个临界强度比确定水力压裂裂纹遇层理时在层理处起裂和沿层理扩展的层理弱面强度范围,以此表示水力压裂裂纹转向层理扩展的难易程度.通过对裂纹扩展判据的分析得出:层理起裂弱层和岩石基体临界强度比随层理走向线与第一主应力夹角和层理倾角的减小以及第三主应力和岩石基体强度的增大而增大;层理走向角小于35.26°时,层理起裂弱层和岩石基体临界强度比随第一主应力的减小以及第二主应力的增大而增大;反之,层理起裂弱层和岩石基体临界强度比随第一主应力的减小以及第二主应力的增大而减小;层理扩展弱层和岩石基体临界强度比随层理走向线与第一主应力夹角、层理倾角和地应力差的减小以及岩石基体抗拉强度的增大而增大.层理起裂条件与层理扩展条件同时满足时,水力压裂裂纹转向层理方向扩展.      

Finite Element study of fracture initiation in flaws subject to internal fluid pressure and vertical stress

Hydraulic fracturing is a method used routinely in oil and gas exploitation and in engineered geothermal systems. While used frequently, there are many aspects of hydraulic fracturing, such as the direction of propagation of the newly-created fractures, which are not very well understood. Even though it is known that the local stress field plays a fundamental role in the orientation of the new fractures, there may be other factors, such as the geometry of the existing fractures and the magnitude of the hydraulic pressure applied, that may play a major role in the path that a new fracture follows when pressurized.The main goal of this study is to numerically analyze the effect of the ratio between a vertical load, or stress, and the hydraulic pressure applied in existing flaws on the stress field in the vicinity of the flaw tips. For that purpose, a double flaw geometry 2a-30-30 was modeled in the Finite Element code ABAQUS, and different vertical loads and internal flaw pressures were applied to the model. The variation of the maximum principal stresses and maximum shear stresses around the flaw tips were analyzed and related to fracture initiation.The study showed that the ratio between the water pressure applied in the flaws and the vertical load/stress (WP/VL) plays a crucial role in the magnitude and shape of the stress field around a flaw tip, and therefore in the location of tensile and shear fracture initiation. As WP/VL increases, the location of initiation of new tensile fractures shifts from the upper face of the studied flaw towards the region right ahead of the flaw tip; simultaneously, the location of initiation of new shear fractures shifts from the region ahead of the flaw tip to the upper face of the analyzed tip.     

Model of stepwise propagation of the tip of a hydraulic fracture in the absence of filtration

Quasi-static stepwise propagation of a hydraulic fracture in rock with a regular structure in the absence of filtration is considered. It is proposed to use a brittle fracture diagram taking into account the hydraulic fracturing fluid pressure and the confining pressure. Fracture curves describing the brittle rock fracture where the hydraulic fracturing fluid partially fills the fracture are constructed and used to predicted the possibility of stepwise propagation of hydraulic fracturing in the case where the fluid gradually flows into the fracturing crack. The regularity of the structure of the brittle rocks fracture is estimated from the results of two full-scale experiments: the critical stress intensity factor and the tensile strength limit of the rock. Experiments on pulsed loading of polymethylmethacrylate samples with stepwise crack propagation along concentric circular arcs were performed. The results of the experiments are consistent with theoretical predictions.     

页岩气藏压裂缝网扩展数值模拟

为探究页岩气藏水力压裂复杂裂缝网络的形成机理,开展了缝网扩展的数值模拟研究.考虑应力阴影和天然裂缝作用,建立了井筒和裂缝中流体流动模型,利用位移不连续方法求解应力与位移不连续量,然后构建了压力与裂缝宽度的迭代方程,并采用牛顿迭代法求解.通过比较数值解经典模型解析解,验证了模型和迭代解法的正确性.多簇裂缝同步扩展时裂缝间距越小,压裂液分配到各条裂缝越不均匀,靠近井筒跟部的裂缝的分流量越大,从而裂缝宽度越大;考虑天然裂缝作用时,逼近角越小或者应力各向异性越弱,水力裂缝越容易发生转向扩展,裂缝网络越复杂.      

梁在纯弯曲力矩作用下的动态塑性断裂

本文采用刚塑性分析方法,研究了带裂纹有限长韧性材料梁在纯弯曲力矩作用下的动态塑性断裂过程,本文考虑了由于断裂引起的附加轴向力对断裂过程的影响,并在分析中考虑了应变率对断裂过程的影响。文章给出了断裂过程中断裂截面上弯曲力矩和轴向力随时间的变化规律,以及裂纹长度,裂纹扩展速度和加速度随时间的变化规律。     

页岩储层体积压裂的地应力变化研究

页岩储层属于致密超低渗透储层,需改造形成复杂缝网才有经济产能.体积压裂是页岩储层增产改造的主要措施,而地应力场特别是水平主应力差值是体积压裂的关键控制因素. 理论研究表明:(1)当初始两向水平主应力差较小时,容易形成缝网,反之不易产生缝网;(2)人工裂缝的形成能够改变地层初始应力场. 因此应在前人研究的基础上优化设计压裂方式,以克服和翻转初始水平主应力差值,产生体积缝网.基于此,建立了页岩气藏水平井体积压裂数值模型,模型中采用多孔介质流固耦合单元模拟页岩基质的行为,采用带有孔压的"cohesive"单元描述水力裂缝的性质,模型对"Texas Two-Step" 压裂方法进行了数值模拟,模拟结果得到了压裂过程中地层应力场的分布及其变化,模拟结果和解析公式计算结果吻合良好.模拟结果表明:(1)裂缝的产生减弱了地层应力场的各向异性;(2 对于低水平应力差页岩储层,采用"Texas Two-Step"压裂方法可以产生缝网. 对于采用"Texas Two-Step"压裂方法无法产生缝网的高应力差页岩储层,提出了三次应力"共振" 和四次应力"共振" 压裂方法并进行了数值模拟,模拟结果得到了压裂过程中页岩储层应力场的分布及其变化,得到了缝网形成的区域,模拟结果表明:(1)对于高应力差页岩储层,采用"Texas Two-Step" 压裂方法无法产生缝网;(2)对于高应力差页岩储层,三次应力"共振" 和四次应力"共振"压裂方法是有效的体积压裂缝网形成的方法.      

A variational approach to analyze a natural fault with hydraulic fracture based on the strain energy density criterion

A simplified analytical model of the interaction between a hydraulic fracture and an existing natural fault is developed. The mechanical activation of the natural fault as a result of contact with a pressurized fracture is described for plane strain conditions and quasi-static fracture propagation approximation. Using a variational approach, the normal and shear stresses, as well as the boundaries of the open and sliding zones along the fault, are predicted for three stages of the fracturing process (fracture approaching, coalescence, and fluid penetration). An accumulated concentration of shear stress at the tip of the fault’s sliding zone is shown to create sufficient tensile stress to initiate a new tensile crack on the opposite side of the fault, provided either the differential in situ stress is low or the friction coefficient is sufficiently large. The results of direct numerical simulation of the fracture interaction fit the model predictions made from the strain energy density fracture criterion.     

一种XFEM断裂分析的裂尖单元新型改进函数

提出了一种适用于裂尖改进单元的新型改进函数, 基于三角变换的方法, 保留裂纹尖端场的应力奇异性和裂纹上、下表面的位移不连续性, 将常规扩展有限元法裂尖改进单元的4 项改进函数缩减为2 项, 裂尖改进单元的结点由常规的8 个改进自由度减少为4 个. 采用2 个正交的水平集函数表征材料内部裂纹面, 详细阐述了改进单元类型的判别方法, 给出一种改进单元的分区域积分方案. 最后, 若干断裂力学问题经典算例的数值计算结果表明:建议的裂尖改进函数具有较高的数值精度, 该方法是十分有效的.     

LC4—M材料复合型韧断主要影响因素的分析

通过分析LC4－M铝合金材料在不同复合比载荷下的断裂试验结果,参照常规断裂现象,修正了一般断裂试验中认定裂纹启裂方向的方法,结合不同复合比下裂尖附近应力三维、主应力方向的计算分析,得到：在裂端的钝化变形区域、应力三维度的极大值处,对应于裂纹的启裂位置,即使在高韧性材料中发生剪切断裂的情况下也是如此;裂纹的启裂方向在拉伸断裂时与启裂点最大拉应力方向有关,在剪切断裂时启裂点最大剪应力方向有关,引起两种形式断裂的主要因数和破坏机理有很大不同。     

THE STRESS FIELD AT MODE I CRACK TIP ON STEEL PIPE UNDER AXIAL TENSION

采用弹性断裂力学Westergaard 的方法,分析在轴向拉伸载荷作用下钢管的Ⅰ型裂纹尖端处的应力场. 基于Ⅰ型裂纹钢管应力场的特征,设定尖端处的应力艾雷函数,给出其应力场模型边界条件和应力场的解析函数,并利用裂纹尖端处的切平面研究裂纹尖端局部应力场,建立了钢管裂纹尖端应力场模型. 通过钢管与平板Ⅰ型裂纹应力场模型的对比,结果表明二者明显不同,钢管裂纹尖端处应力峰值影响范围仅与裂纹长度、拉应力相关.     

基于混合方法的二维水力压裂数值模拟

水力压裂在页岩气开采中被广泛使用,采用数值方法研究压裂机理具有重要意义.基于连续-非连续单元法(CDEM) 和中心型有限体积法(FVM),提出解决水力压裂流固耦合问题的二维混合数值计算模型.该混合模型中,使用CDEM 求解应力场和裂缝扩展过程,使用FVM 求解裂隙渗流场.应力场裂缝扩展和渗流场均使用显式迭代求解, 并通过相互之间数据交换实现流固耦合.通过与KGD 理论模型进行对比, 验证数值模型的正确性.通过与颗粒离散元数值结果进行对比,验证数值模型的有效性.通过计算复杂缝网压裂模型,研究水力压裂机理,并说明该数值模型在水力压裂模拟中具有很好的前景.      

横观各向同性油气藏水力压裂裂纹扩展规律研究

针对横观各向同性与各向同性油气藏水力压裂裂纹扩展的差异性,基于扩展有限元法建立水力压裂力学模型,通过ABAQUS子程序开发了各向同性和横观各向同性岩体的起裂判据。在各向同性岩体数值模拟结果与解析解以及现场压裂典型曲线对比吻合的基础上,得到了包含层理构造的横观各向同性岩体水力压裂过程中裂纹扩展规律。层理类岩体水力压裂的裂缝扩展方向由地应力状态、层理方向以及岩体与层理界面抗拉强度共同决定;水力压裂过程中,注水压力在裂纹尖端产生应力集中,层理面法向分量先达到界面抗拉强度时,裂纹沿层理方向开裂,反之裂纹沿垂直最小地应力的方向扩展;裂纹扩展速度随层理抗拉强度的增加而降低;由于地层的滤失,随压裂液的注入,裂纹长宽尺度增长速率降低。     

Investigation on critical triaxiality along unified yield loci at crack tip under mixed mode (I + II) fracture

Metallic components used in industries and day to day appliances often contain micro-cracks. In general, cracks occur in various orientations to the loading axis. The present paper discusses the criticality of stress triaxiality, a well-known ductile fracture parameter, on the yield loci at the crack tip. In the process, an old model of stress triaxiality has been generalized using unified strength theory to incorporate various convex and nonconvex failure criteria, including single shear, twin shear, etc. The new triaxiality model also reveals about the effect of intermediate principal stress at the crack tip for materials with and without strength difference. The crack initiation angles at the crack tip, obtained through the proposed model have been found to be in unison with those obtained through other fracture criteria.     

混凝土黏聚开裂模型若干进展

黏聚模型是用来描述混凝土断裂行为的基本模型, 首先介绍了混凝土的黏聚开裂模型的基本概念,总结了确定黏聚区的本构方程的各种方法,即直接单轴拉伸测试、J积分方法、R曲线法、柔度法和逆推法.然后介绍了黏聚模型在I型和复合型裂纹问题、疲劳断裂问题中的应用以及黏聚模型与混凝土尺寸效应的关系.最后对黏聚开裂模型与桥联模型、带状裂缝模型进行了比较和总结, 指出了该模型存在的问题, 并对其以后的发展方向提出了建议.      

页岩层理方位及强度对水力压裂的影响

为了研究页岩天然层理倾角及强度等对水力压裂裂纹扩展的影响,采用室内水力压裂实验,通过监测孔直接对裂纹扩展的实时监测和注水压力信息及试件压裂后的剖切,分析层理倾角、强度等对压裂裂纹扩展的影响。实验结果表明:水力压裂过程中,垂直最小地应力稳定扩展的主裂缝遇层理时,层理面与主裂缝初始扩展方向夹角越小,主裂缝越易沿着层理面方向扩展,层理面与主裂缝初始扩展方向夹角越大,主裂缝遇层理面时越易贯穿层理面沿原方向扩展;层理方位,地应力及基质抗拉强度不变,层理的抗拉强度远弱于基质抗拉强度时,主裂缝与层理面相遇后越易沿着层理面方向扩展,层理抗拉强度与基质抗拉强度越相近,主裂缝遇层理时越易贯穿层理沿原方向扩展;层理方位和强度不变,地应力及应力差越大,主裂缝遇层理后越易贯穿层理面沿原方向扩展。     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

Non-local theory solution for in-plane shear of through crack

A non-local theory of elasticity is applied to obtain the plane strain stress and displacement field for a through crack under in-plane shear by using Schmidt's method. Unlike the classical elasticity solution, a lattice parameter enters into the problem that make the stresses finite at crack tip. Both the angular variations of the circumferential stress and strain energy density function are examined to associate their stationary value with locations of possible fracture initiation. The former criterion predicted a crack initiation angle of 54° from the plane of shear for the non-local solution as compared with about 75° for the classical elasticity solution. The latter criterion based on energy density yields a crack initiation angle of 80° for a Poisson's ratio of 0.28. This is much closer to the value that is predicted by the classical crack tips solution of elasticity.     

A Zener–Stroh crack interacting with a coated inclusion with generalized Irwin plastic zone correction

The elastic–plastic fracture behavior of a Zener–Stroh crack interacting with a coated inclusion in composite materials has been investigated with crack tip plastic zone corrections. With the distributed dislocation method, the crack problem is formulated into a set of singular integral equations which are solved numerically. The plastic zone sizes at the both crack tips are determined by a generalized Irwin model where Von Mises stress yielding criterion is used. The stress intensity factor (SIF), the plastic zone size (PZS), the crack tip opening displacement (CTOD) and the effective stress intensity factor have been evaluated. In the numerical examples, the influence of the inclusion shear modulus, the coating-layer thickness and shear modulus, as well as the distance between the crack and inclusion, on the SIF, the PZS and the CTOD are discussed in detail. Numerical examples show that increasing the shear modulus or the thickness of the coating phase, the influence of the inclusion on the normalized SIF and the normalized PZS will be shielded.     

基于能量耗散的岩石脆性评价方法

脆性评价是页岩气钻井以及压裂工艺设计的重要环节.然而,现今对于脆性的定义十分模糊,且脆性评价的方法没有成熟的理论作为支撑.本文以岩石应力应变关系为基础,采用能量法来进行储层岩石脆性评价.根据断裂力学裂纹尖端的非线性场理论,以前端区作为尖端特殊的屏蔽构型,来描述储层岩石在断裂过程中能量的耗散. 通过计算本征内聚力与前端区稳态阻力比值来表征储层岩石的脆性指数. 研究表明,能量耗散相对越大,岩石的脆性指数越低,其塑性特征表现得越明显. 基于能量耗散的储层岩石脆性评价是分析岩石力学特性的重要内容,也为压裂选层提供了重要依据.      

纤维复合材料复合型断裂正应力强度因子比准则

对纤维增强复合材料的裂纹起裂及开裂方向准则进行了研究，提出了复合型断裂的正应力强度因子比准则（Normal Stress Intensity Factor Raito Criterion).此准则是一种综合考虑了正应力强度因子和剪应力强度因子对裂纹起裂的推动的准则，并且不需要预先确定材料的特征尺度，使用较方便，且预测结果是很好的。