多簇压裂干扰应力变化规律及对裂纹扩展的影响

为研究水平井分段多簇压裂缝间的干扰应力及其对裂纹扩展的影响,在现有二维未考虑地应力的单裂缝干扰应力解析解的基础上,利用双平面复变函数保角变换得到了包含地应力项的三维干扰应力解析解。基于扩展有限元法建立三维多裂缝扩展力学模型,利用Python脚本二次开发平台实现了三维多裂缝水力压裂参数化建模,通过解析解与数值计算对比分析,得到如下结论。裂纹两侧裂纹面法向和走向干扰正应力分别为压应力和拉应力,均呈纺锤形,法向干扰应力影响范围大约为走向干扰应力的5倍;裂纹尖端裂纹面法向和走向干扰正应力分别为拉应力和压应力;裂纹尖端两侧存在干扰剪应力;考虑初始地应力对干扰应力解析解进行修正后的干扰应力值均变小;多簇压裂中裂缝间的干扰应力叠加,簇间距越小,叠加效果越强;多簇压裂的干扰应力使裂缝间裂纹面法向压应力增大,走向压应力减小,导致裂纹扩展注水压力升高,裂缝张开的宽度降低,不利于单裂缝的扩展;干扰应力使裂缝间应力差降低,甚至局部最小地应力方向发生改变,有利于形成复杂缝网。     

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

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

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

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

动态扩展裂纹的若干反平面问题的研究

采用复变函数论，对反平面条件下的动态裂纹扩展问题进行研究。通过自相似函数的方法可以获得解析解的一般表达式。应用该法可以很容易地将所讨论的问题转化为Riemann—Hilbert问题，并可以相当简单地得到问题的闭合解。文中分别对裂纹面受均布载荷、坐标原点受集中增加载荷、坐标原点受瞬时冲击载荷以及裂纹面受运动集中载荷Px／t作用下的动态裂纹扩展问题进行求解，得到了裂纹扩展位移、裂纹尖端的应力和动态应力强度因子的解析解。应用该解并通过叠加原理，就可以求得任意复杂问题的解。     

三维压裂缝网不稳定压力半解析求解方法

受地应力及压裂工艺影响, 大斜度井水力压裂缝网展布复杂, 缝网中存在不同倾斜方向、不同展布形态及不同贯穿程度的压裂缝. 本文通过将裂缝面离散为若干矩形微元实现裂缝形态有效表征, 将渗流过程划分为基质向裂缝流动及裂缝向井筒流动两阶段, 采用有限差分方法构建离散裂缝面内不稳定渗流数值解, 结合封闭边界面源函数及叠加原理构建基质内不稳定渗流解析解, 耦合裂缝内流动数值解与基质内流动解析解, 求解了三维压裂缝网不稳定压力. 基于积分中值定理提出了点源、特殊线源代替面源求解基质内渗流的求解方法, 分析了该方法的可行性及适用条件, 在保证模型精度的同时提升了计算效率. 研究表明, 在基质内采用点源函数面积分求解面源的方法可准确求解三维压裂缝网井底压力动态但计算效率极低, 基于积分中值定理的点源、特殊线源近似面源求解方法可大大提升计算效率, 且在裂缝微元划分较为精细(微元无因次边长小于0.15)时可取得较高精度, 基于该模型分析了裂缝导流能力、裂缝倾角、裂缝高度及裂缝段间距对压裂大斜度井典型试井曲线的影响.      

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

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

压-剪混合型定常扩展裂纹尖端的弹黏塑性场

假定黏性系数与塑性等效应变率的幂次成反比，考虑其黏性和裂纹面摩擦接触效应 建立了压-剪混合型定常扩展裂纹尖端弹黏塑性场的渐近方程，求得了裂纹尖端场不含应力、应变间 断的数值解. 并讨论了压-剪混合型裂纹数值解随各个参数的变化规律，计算结果 和分析表明，压-剪混合型裂纹尖端场是满塑性的，不含有弹性卸载区，黏性效应是研究扩展裂纹尖端场时的一个重要因素. 无论混合裂纹趋近I型还是趋近II型，静水压力随摩擦系数的增加都是增加的,裂纹面摩擦 效应是阻止裂纹扩展速度的因素，且摩擦作用越强，裂纹尖端场的韧性越高.     

三维切口/裂纹结构的扩展边界元法分析

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.      

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

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

复合射孔爆燃气体压裂裂缝起裂扩展研究

结合线弹性断裂力学的裂缝尖端应力强度因子判据,建立了复合射孔爆燃气体压裂裂缝的起裂扩展模型,通过建立与多个变量相关的缝内气体压力分布函数,利用迭代法实现了模型的数值求解,获得了缝内气体压力分布随时间的动态变化规律,并分析了不同特征参数对裂缝起裂扩展与止裂过程的影响。实例计算结果表明:(1)随着裂缝扩展的进行,爆燃气体流动尖端与裂缝尖端经历了由重合到不重合再到重合的过程;(2)地应力越大,裂缝起裂扩展越困难,爆燃气体有效致裂作用时间越短,最终得到的裂缝扩展长度也越小;(3)初始裂缝越长,裂缝更容易起裂扩展,爆燃气体能量利用率越高,裂缝扩展更长;(4)岩石断裂韧性的改变对裂缝起裂、止裂和裂缝扩展长度没有明显的影响;(5)升压速率越小,爆燃气体有效致裂作用时间越长,最终裂缝扩展也更长,但对裂缝起裂压力与止裂压力几乎没有影响。     

压应力场中爆生裂纹分布与扩展特征实验分析

采用动静组合加载实验装置和数字激光焦散线实验系统,进行了0、3、6、9 MPa等4种压应力场中PMMA试件的爆破致裂实验,分析了沿静态主应力方向扩展的裂纹运动学和力学行为。实验结果表明:首先,静态竖向载荷在预制炮孔周围产生应力集中,在炮孔壁上下端部处出现最大拉应力;随后,在动态爆炸载荷的叠加作用下,裂纹优先在炮孔壁上最大拉应力位置处起裂,并沿最大主应力方向扩展;裂纹扩展过程中,静态竖向载荷越大,裂纹扩展速度越大,且裂纹尖端应力强度因子值越大。     

The dynamic behaviour of an oblique edge-crack under impact loading

When a stress wave (tensile or compressive) impinges on a crack existing in an elastic medium, reflection, refraction and diffraction-phenomena take place. A result of diffraction is the loading of the crack. While compressive stress-waves do not create any stress concentration at the tip of an existing crack, tensile stress-waves develop stresses at the tip which may cause a propagation of the crack. If the tensile pulse is weak the crack may propagate by steps under the action only of successive tensile stress-pulses, whereas intermediate compressive-stress pulses do not have any influence.A complete study of the phenomena of incubation, initiation and propagation of cracks in thin plates, when they are subjected to a compressive pulse, which is subsequently reflected from the free boundaries of the plate and changed to complicated wave-trains, was undertaken in this paper, based on the method of caustics. Interesting results that were derived from this experimental study are presented.     

Reverse shear on inclined planes at the tip of a sharp crack

Plane strain plastic yielding at a crack tip has been represented by edge dislocations with Burgers vectors parallel to symmetrical planes inclined at 70° and 45° to the plane of the crack. The plastic displacement and the stresses near the crack tip were calculated by a numerical method and the effect of a reduction in applied stress was determined. Removal of the whole or a part of the initial load produces reverse shear in regions of the slip band nearest the crack tip. The amount of reverse shear depends only on the reduction in the load and not on its initial value. The reverse shear is associated with the presence of negative dislocations and the stresses near the crack tip may become compressive even though the applied (remote) stress is still tensile. The degree and extent of compression depends on the reduction in applied stress and on its original value. It is argued that the residual compressive stresses produced under fluctuating loads may produce crack closure and crack arrest. The effect of residual plasticity in a slip band left behind a growing crack has been estimated. It is shown that after an overload the excess residual plasticity opposing crack opening rises to a maximum value when the crack tip has advanced some distance from the point where the overload was applied.     

A Study on Hydraulic Autofrettage of Thick-Walled Cylinders Incorporating Bauschinger Effect

Pressure vessels which are subjected to cyclic external or internal high pressure are used in many fields of industry and need to be sure of reliability and safety. To ensure of reliability and safety, thick-walled cylinder, such as a cannon or nuclear reactor, is autofrettaged to induce advantageous residual stresses. The compressive residual stress which was introduced by autofrettage process acts to offset the tensile residual stress induced by internal pressure. It increases operating pressure and restrains crack initiation and crack propagation. As the autofrettage level increases, the magnitude of compressive residual stress at the bore also increases. However, the Bauschinger effect reduces the compressive residual stresses with prior tensile plastic strain, and decreases the beneficial autofrettage effect. There are some differences between theoretical solution considering elastic-perfectly material behavior and real autofrettage process results. The purpose of the present paper is to predict the accurate residual stress of SNCM 8 high strength steel using the Kendall model which was adopted by ASME Code. The tensile and uniaxial Bauschinger effect tests of SNCM 8 were performed to evaluate Bauschinger effect factor(BEF), thereafter this constant was used in calculating the residual stress. The residual stress distribution which is considering the Bauschinger effect was profiled using Kendall model, and the results were compared with the analytical and Finite Element analysis. The results were found that residual stress incorporating the Bauschinger effect at bore was smaller than ideal calculation. These results should be considered in designing pressure vessels.     

高地应力岩体多孔爆破破岩机制

深部岩体爆破破岩是爆炸荷载与高地应力共同作用的结果。基于一些简化假设，建立了一个高地应力岩体双孔爆破计算模型，采用光滑粒子流体力学-有限元方法耦合数值模拟方法，研究了高地应力作用下炮孔间裂纹的传播及贯通过程，分析了炮孔周围应力场动态演化过程与分布特征。研究结果表明:爆破引起的岩体开裂主要是环向动拉应力所致，地应力对岩体的压缩降低了炮孔周围环向动拉应力、缩短了环向动拉应力的作用时间，因而对爆炸致裂起抑制作用；静水地应力条件下多孔爆破时，垂直于炮孔连线方向传播的爆生裂纹更易受到地应力的抑制；对于高地应力岩体爆破，炮孔间的裂纹扩展长度随地应力水平的提高而减小，裂纹主要沿最大主应力方向扩展，因此沿最大主应力方向布置炮孔、缩短炮孔间距有利于炮孔间裂纹的连接贯通，形成良好的爆破开挖面。     

Plastic yielding on inclined slip-planes at a crack tip

Plastic yield at crack tips on singular slip-planes, inclined to the crack plane, has been studied under plane-strain conditions for combined tension, hydrostatic stress, and in-plane shear. The singular integral equation, which represents the equilibrium condition of edge dislocations on the slip-planes, is transformed into a Fredholm integral equation in order to avoid difficulties that occur with its numerical solution. Results are presented for the slip-band length, the plastic crack-tip opening displacement, stress fields, and crack-opening contours. A series expansion of the results obtained numerically confirms approximate analytical expressions given by J.R. Rice (1974), up to the third-order in the applied stresses. The results of finite element methods agree with values of the crack-tip opening displacement obtained here to within 10 per cent. Ahead of the crack tip, the principal tensile stresses exceed the principal shear stresses by a factor of 10, approximately.     

混凝土拉伸断裂的细观数值分析

根据混凝土试件拉伸和三点弯曲的物理模型,用梁-颗粒模型BPM 2D(B eam-Particle M ode l)模拟了混凝土拉伸和三点弯曲试件微裂纹的萌生、扩展直至试件宏观破坏的全过程。在梁-颗粒模型中用三种类型梁单元形成混凝土细观数值模型,每种类型梁单元的力学性质均按韦伯(W e ibu ll)分布随机赋值以模拟混凝土细观结构的非均匀性。数值模拟结果给出了混凝土拉伸应力-应变曲线和三点弯曲载荷-位移曲线,以及混凝土试件破坏过程最大应力分布图和裂纹扩展图。数值模拟结果显示混凝土破坏过程实际上就是微裂纹萌生、扩展、贯通,直到宏观裂纹产生导致混凝土失稳断裂的过程。通过对数值模拟结果的分析,揭示出混凝土在拉伸条件下裂纹尖端的拉应力集中是裂纹扩展的动力,混凝土组成材料力学性质的非均匀性是造成裂纹扩展路径曲折的重要原因。     

Retardation of a crack with connections between the faces using an induced thermoelastic stress field

This paper considers local temperature variations near the tip of a crack in the presence of regions in which the crack faces interact. It is assumed that these regions are adjacent to the crack tip and are comparable in size to the crack size. The problem of local temperature variations consists of delay or retardation of crack growth. For a crack with connections between the crack faces subjected to external tensile loads, an induced thermoelastic stress field, and the stresses at the connections preventing crack opening, the boundary-value problem of the equilibrium of the crack reduces to a system of nonlinear singular integrodifferential equations with a Cauchy kernel. The normal and tangential stresses at the connections are found by solving this system of equations. The stress intensity factors are calculated. The energy characteristics of cracks with tip regions are considered. The limiting equilibrium condition for cracks with tip regions is formulated using the criterion of limiting stretching of the connections.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 133–143, January–February, 2005     

Analysis of crack extension in anisotropic materials based on local normal stress

The normal stress ratio theory is applied to predict crack extension behavior in center-notched unidirectional graphite-epoxy of arbitrary fiber axis orientation, subjected to arbitrary far-field planar loading. The theory is applied within analytical solutions for two infinite plate geometries: a plate with a sharp center crack, and a plate with an elliptical center flaw. A critical analytical case is identified suggesting that application of the theory within a stress solution modelling crack tip shape may increase the accuracy of crack growth direction predictions. Crack extension direction, location of crack extension, and critical stress predictions of the theory are compared to those obtained from experiments on specimens subjected to tensile, shear, and mixed-mode far-field loading. The comparison shows that, applied within each analytical solution, the normal stress ratio theory provides verifiable predictions of crack growth behavior. By modelling actual notch tip shape, the elliptical notch solution is able to provide accurate qualitative predictions of the origin of crack extension along the periphery of a cut notch tip in a way that the sharp crack analysis cannot. The sharp notch solution appears to provide slightly more accurate crack growth direction predictions, however. Also, in predicting critical applied far-field stresses, the sharp crack solution appears to exhibit a stronger ability to model subtle experimental trends.