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

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

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

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

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

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

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

应用复变函数论，对反平面动态扩展裂纹问题进行了研究。通过自相似函数的方法可以获得若干问题的解析解。应用该法可以迅速地将所论问题转化为Riemann-Hilbert问题，并可以相当简单地得到问题的闭合解。通过叠加原理利用这些解，就可以求得任意复杂问题的解。     

星形裂纹的应力分析

利用复变函数的方法, 通过构造适当的保角映射研究了星形裂纹的平面弹性问题,给出了裂纹尖端I型与II型问题应力强度因子的解析解.并由此模拟出了经典的Griffith裂纹,共点均匀分布三裂纹,十字裂纹,对称八裂纹问题.      

压应力对压剪裂纹扩展的影响研究

本文针对压剪裂纹的启裂及扩展问题,研究了脆性材料中裂纹面压应力变化对其扩展的影响规律.借助双轴加载试验机可自由调节两个轴位移和力的优势,设计了一种单边对称双裂纹压剪试样.试验中,施加裂纹面压应力至不同的预设值后,使剪应力单调增大直至裂纹启裂及扩展,得到不同预设压应力下压剪裂纹启裂及扩展规律. 随着预设压应力增大,启裂角增大,裂纹扩展路径与初始裂纹的偏角也增大.但随着压应力增大,启裂角增大至一定值后趋于稳定,实验发现,依据裂纹是否闭合,压应力对压剪裂纹扩展的作用大致可分为两个阶段:闭合前,压应力对裂纹启裂载荷及启裂角、扩展路径均有影响,预设压应力增大,裂纹启裂载荷增大、启裂角增大,扩展路径愈来愈偏离初始裂纹方向;闭合后,压应力虽然增大,启裂角和临界压剪应力比始终恒定,压应力对临界剪力和扩展路径存在一定影响.研究表明,裂纹启裂角与启裂时的压剪应力比存在一定的对应关系.启裂时的压剪应力比增大,启裂角增大,启裂时的压剪应力比恒定,启裂角不变.      

非对称III型界面裂纹扩展的动态问题

通过复变函数论的方法,对非对称Ⅲ型界面裂纹扩展的动态问题进行了研究.采用自相似函数的方法可以轻易地将所论问题转化为Riemann-Hilbert问题,并求得了裂纹坐标原点分别受到变载荷$Pt/ x$, $Px^3 /t^2$作用下的解析解的一般表达式.通过Muskhelishvili方法可以相当简单地得到问题的闭合解. 利用这些解并采用叠加原理,可以求得任意复杂问题的解.      

非对称Ⅲ型界面裂纹扩展的动态问题

通过复变函数论的方法,对非对称Ⅲ型界面裂纹扩展的动态问题进行了研究.采用自相似函数的方法可以轻易地将所论问题转化为Riemann-Hilbert问题,并求得了裂纹坐标原点分别受到变载荷Pt/z,Px3/t2作用下的解析解的一般表达式.通过Muskhelishvili方法可以相当简单地得到问题的闭合解.利用这些解并采用叠加原理.可以求得任意复杂问题的解.     

修正的拉应力裂纹扩展准则及裂隙\=水压对裂纹扩展的影响

对脆性材料的第一主应力--拉应力裂纹扩展准则进行了补充和修正,修正的裂纹扩展准则能确定裂纹扩展步长.以平面斜置裂纹扩展为例,利用无网格Galerkin方法,对不含裂隙水压的二维裂纹扩展进行数值模拟,计算结果与试验结果一致,表明最大周向拉应力准则的正确性.在不同裂隙水压条件下,研究了二维裂纹初始破裂,并在给定水压下对二维裂纹扩展路径进行了数值模拟跟踪.结果表明裂隙水压对裂纹初始破裂方向、破裂步长、破裂载荷和裂隙岩体破裂强度有显著影响.有水压和无水压的扩展迹线不同,但后续的扩展趋势相同.     

页岩水力压裂中多簇裂缝扩展的全耦合模拟

水平井和水力压裂是页岩气开发中的关键技术。对水力压裂中多簇裂缝同时扩展的物理过程进行了数值模拟。采用扩展有限元法(XFEM)模拟岩石中裂缝沿着任意路径扩展,采用有限体积法(FVM)模拟裂缝中流体的流动,并且考虑井筒中流体流动以及在各簇裂缝间的流量动态分配。通过牛顿迭代对全耦合物理过程进行数值求解,重点研究了初始长度不同的两条裂缝的扩展过程,证明较大的射孔摩阻能促进两条裂缝的同时扩展。并通过算例证明了本方法的精度和有效性。     

Asymmetrical dynamic propagation problems on mode Ⅲ interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on mode Ⅲ interface crack are studied. The universal representations of analytical solutions are obtained by the approaches of serf-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.     

Asymmetrical dynamic propagation problems on mode III interface crack

By the application of the theory of complex functions, asymmetrical dynamic propagation problems on modeⅢinterface crack are studied. The universal representations of analytical solutions are obtained by the approaches of self-similar function. The problems researched can be facilely transformed into Riemann-Hilbert problems and analytical solution to an asymmetrical propagation crack under the condition of point loads and unit-step loads, respectively, is acquired. After those solutions were used by superposition theorem, the solutions of arbitrarily complex problems could be attained.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Dynamic propagation problems concerning surfaces of asymmetrical mode III crack subjected to moving loads

With the theory of complex functions, dynamic propagation problems concerning surfaces of asymmetrical mode III crack subjected to moving loads are investigated. General representations of analytical solutions are obtained with self-similar functions. The problems can be easily converted into Riemann-Hilbert problems using this technique. Analytical solutions to stress, displacement and dynamic stress intensity factor under constant and unit-step moving loads on the surfaces of asymmetrical extension crack, respectively, are obtained. By applying these solutions, together with the superposition principle, solutions of discretionarily intricate problems can be found. Project supported by the Post-Doctoral Science Foundation of China (No. 2005038199) and the Natural Science Foundation of Heilongjiang Province of China (No. ZJG04-08)     

DYNAMIC PROPAGATION PROBLEM ON DUGDALE MODEL OF MODE Ⅲ INTERFACE CRACK

By the theory of complex functions, the dynamic propagation problem on Dugdale model of mode Ⅲ interface crack for nonlinear characters of materials was studied. The general expressions of analytical solutions are obtained by the methods of self-similar functions. The problems dealt with can be easily transformed into Riemann-Hilbert problems and their closed solutions are attained rather simply by this approach. After those solutions were utilized by superposition theorem, the solutions of arbitrarily complex problems could be obtained.     

基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。     

应变率对含裂隙红砂岩裂纹扩展模式及破碎特征的影响

以含不同倾角预制裂纹的长方形板状红砂岩为研究对象,采用分离式霍普金森压杆沿试样宽度方向施加冲击荷载,使用高速摄像机记录裂纹扩展过程,获得试样的裂纹路径特征以及动态压缩强度和动态弹性模量,利用筛分统计法分析试样碎块分布特征,结合分形理论定量描述试样破碎程度及特点,探讨中应变率条件下含裂隙试样裂纹扩展模式与动态力学性质和破碎程度的相互关系。研究结果表明,应变率较高时试样会更多地出现远场裂纹和离层裂纹,并且相比相关低应变率实验结果,中应变率范围内试样破坏模式及裂纹分布情况随应变率的变化规律是不同的。随着应变率的提高,试样大体上从1条拉伸裂纹的临界破坏演变成X形剪切裂纹为主的复杂裂隙网络,并且不同角度预制裂隙对于这种裂纹扩展模式的演变有重要影响。在预制裂纹倾角一定的情况下,岩样动态压缩强度和动态弹性模量表现出明显的应变率效应,不同角度预制裂纹对于试样的应变率敏感性有显著影响。随裂纹倾角的增大,试样的动态强度、动态弹性模量和分形维数表现出的变化趋势具有一定的相似性,大体呈现先减小后增大的趋势,裂纹倾角为45°的试样的动态压缩强度、动态弹性模量和分形维数均为最小。随应变率的升高,不同预制裂纹倾角的试样碎块分布更加分散,应变率越高,预制裂纹倾角对于岩石冲击破碎程度、分形维数的影响越显著。

     

Stress-strain field near crack tip and calculation of critical stress of crack propagation in a pure bending beam of rectangular section with one-sided mode i crack

This paper studies the stress-strain field near crack tip in a pure bending beam of rectangular section with one-sided mode I crack by the analytic method of Ref. [1], then it gives the stress and strain components at the crack tip when the crack propagates and further it obtains the formulas of calculating the elastic deformed area width, the deformed intensity area width and the equation groups of calculating the critical stress of crack propagation, last the equation group of calculating critical stress of crack propagation is verified by calculating instance. The maximum error is 0.18%. First Received May 7, 1994.     

Analysis of stress intensity factor in orthotropic bi-material mixed interface crack

Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With consideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial differential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit,undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship between the bi-material stress intensity factor and the ratio of bi-material Young's modulus are given in the numerical analysis.