垂直裂缝在砂泥岩层理界面上拐折扩展特性研究

首先以大庆油田薄差油藏的砂泥岩储层力学参数为基础,根据相似原理,制备了砂泥岩相似材料试件。然后通过三点弯曲实验构建了Ⅰ型裂缝扩展,并采用数字散斑相关方法,获得了裂缝扩展过程中试件表面散斑场的变化规律。由此给出了界面位置处裂尖的应力奇异性阶数及应力强度因子,了解了材料强度差异对裂缝在界面扩展特性的影响。分析表明,当垂直裂缝从高强度材料一侧扩展时,界面的应力奇异性阶数大,界面的应力集中效应高,KⅡ/KⅠ比值小,裂缝拐折程度小,反之亦然。     

巴西圆盘泥岩试件裂纹扩展及应变演化实验研究

为了研究准静态加载条件下岩石试件巴西劈裂裂纹扩展规律,采用MTS试验机进行准静态加载,同时用高速摄像机记录裂纹扩展过程。采用白光数字散斑处理软件对摄像机记录的照片进行处理,得到试件裂纹扩展过程中应变场的演化情况。通过实验和分析可以看出,由于端部效应及加载方式的原因,因此裂纹起裂点在底部加载部位;泥岩试件表面裂纹的平均扩展速度为252m/s;岩石的非均质性即内部微缺陷、微裂纹使得泥岩试样的开裂并不是沿着中心直径方向,而是偏离一定的角度,初始偏离角度约为17°。裂纹扩展过程可以划分为三个阶段:泥岩试件宏观变形阶段(宏观无裂纹)、宏观裂纹稳定扩展阶段、宏观裂纹动态张裂阶段。同时,在裂纹扩展过程中,表面第一主应变场、水平位移场等变化明显,在开裂部位第一主应变最大。通过对圆盘泥岩试件裂纹扩展实验研究,可为研究岩石破裂及其演化规律提供依据。     

页岩气储层压裂数值模拟技术研究进展

页岩气储层水力压裂数值模拟既要考虑页岩储层岩石的特性,又要兼顾水平外分段压裂施工工艺,是一个非常棘手的力学难题.本文简述了页岩气储层岩石具有的地质力学特征和页岩气储层开发常用的水平井分段压裂技术;详述了扩展有限元、边界元、离散元在水力压裂裂缝模拟上的应用现状,指出了它们在处理裂缝问题的局限性和优越性,总结出边界元三维位移不连续法是模拟多裂缝扩展的有效方法.     

混凝土结构锈胀开裂的扩展有限元数值分析

依据非均匀锈胀理论提出钢筋锈胀作用的计算方法,应用扩展有限元法(XFEM)建立了钢筋锈胀保护层开裂的有限元模型.数值分析表明:采用XFEM与混凝土黏聚力模型能有效模拟混凝土开裂及裂纹扩展,避免了网格重剖分的问题;预裂纹的存在抑制了混凝土裂纹萌生,却加速了裂纹扩展贯通保护层,且萌生始于预裂纹尖端,而非钢筋-混凝土锈蚀层界面处;初始无损伤结构裂纹萌生位置对称分布于锈蚀层界面一定范围内,裂尖距交界面距离越大,单元受锈胀影响越小,最终贯通保护层主要是锈胀位移与锈蚀产物渗入裂缝产生作用力共同作用的结果,且裂纹扩展角趋于120°;提高混凝土等级和增大保护层厚度能有效延缓锈胀裂缝的产生与发展,有利于提高结构耐久性.     

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

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

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

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

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

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

水力压裂解析模型裂缝扩展参数敏感性分析

水力压裂形成复杂裂缝网络是致密储层油气开采的重要技术,掌握水压裂缝扩展机理是控制压裂行为和优化压裂效果的关键.水压裂缝动态扩展行为涉及储层岩体、注入压裂液、压裂实施工艺等方面,其中水力压裂扩展时间、压裂液流体动力粘度系数、压裂液流体注入流速、储层岩石剪切模量成为决定裂缝扩展长度和裂缝开度的重要因素.本研究采用KGD、P...     

基于单元破裂的岩石裂纹扩展模拟方法

传统离散元方法在处理破裂问题时, 采用界面上的准则进行判断, 裂纹只能沿着单元边界扩展. 当物理问题存在宏观或微观裂隙时, 在界面上应用准则具有其合理性; 而裂纹沿着单元边界扩展, 使得裂纹路径受网格影响较大, 扩展方向受到限制. 针对上述情况, 可以基于单元破裂的方式, 构建连续- 非连续单元法, 并应用于岩石裂纹扩展问题的模拟. 该方法在连续计算时, 将单元离散为具有物理意义的弹簧系统, 在局部坐标系下由弹簧特征长度、面积求解单元变形和应力, 通过更新局部坐标系和弹簧特征量, 可进一步计算块体大位移、大转动, 连续问题计算结果与有限元一致, 同时提高了计算效率. 在此基础上, 引入最大拉应力与莫尔—库伦的复合准则, 判断单元破裂状态和破裂方向, 并采用局部块体切割的方式, 在单元内形成初始裂纹. 裂纹两侧相应增加新的计算节点, 同时引入内聚力模型描述裂纹两侧的法向、切向作用与张开度及滑移变形之间的关系. 按此方式, 裂纹尖端处的扩展路径可穿过单元内部和单元边界, 在扩展方向的选取上更为准确. 最后, 通过三点弯曲梁、单切口平板拉伸、双切口试样等典型数值试验, 模拟裂纹在拉伸、压剪等各种应力状态下的扩展问题, 并对岩石单轴压缩试验的破坏过程进行模拟, 分析裂纹形成与应力—应变曲线各阶段之间的对应关系. 结果表明: 连续—非连续单元法通过单元内部破裂的方式, 可以显示模拟裂纹萌生、扩展、贯通直至形成宏观裂缝的过程.      

EXPERIMENT STUDY OF I TYPE CRACK BLUNT EXTENSION INDUCED BY WEAK PLANE1)

天然裂缝或层理等弱界面的剪切破坏是体积压裂裂缝网络形成的关键因素。利用数字散斑相关技术,通过三点弯曲构建了I型裂缝经过弱界面时的扩展特性实验,获得了试件表面位移场和应变场的变化规律。实验结果表明：I型裂缝扩展至弱界面层后,裂缝发生短暂停滞扩展,裂缝尖端张开位移迅速增加,I型裂缝尖端钝化,弱界面剪切应变迅速增加,裂缝由I型裂缝转变为I-II复合型裂缝并转向扩展。     

基于XFEM-MBEM的嵌入式离散裂缝模型流固耦合数值模拟方法

离散缝网的表征与模拟是目前国内外研究的热点. 在非常规油气开发过程中, 由于地应力场的存在会对裂缝的流动属性产生显著影响, 若将裂缝视为静态对象, 与矿场数据会出现极大偏差, 因此要基于动态裂缝做更深入的研究. 本文针对致密油藏应力场?渗流场耦合力学问题, 提出了一种高效的混合数值离散化方法, 其中采用扩展有限元法 (XFEM) 求解岩石的弹性形变, 采用了混合边界元法 (MBEM) 精确计算基岩与裂缝间的非稳态窜流, 这两种数值格式是完全耦合的, 并对整体计算格式的时间项进行了全隐式求解, 可准确表征致密油藏开采过程中的裂缝变形及流体流动机理. 此外, 本文采用了嵌入式离散裂缝前处理算法显式表征大尺度水力压裂缝, 并考虑了支撑剂的作用; 采用了双孔有效应力原理和双重介质隐式裂缝表征方法, 可捕捉基质与小尺度天然裂缝的动态信息; 由此, 本文所提出的混合模型综合表征了基质?天然裂缝?水力压裂缝共同组成的致密油藏复杂渗流环境, 并通过几个实例论证了模型的准确性, 研究表明: 对致密油藏压裂水平井进行产能评价时, 应力场所引起渗流参数的改变及裂缝开度降低的影响不可忽略. 本文研究可为非常规油气资源的开发提供理论指导.      

A computationally efficient solution of the wave equation for the transient response of infinite reservoirs

In this paper, the transient response of an infinite reservoir is analyzed using the dual-reciprocity boundary element method. A vertical and an inclined-face rigid dam are analyzed under a transient loading. Sharan-type boundary-condition transmission is implemented in the formulation. The results are compared with the exact solution and those obtained by using the finite element method. It is seen that the application of the dual-reciprocity boundary element method is simpler and the results are in very good agreement with the exact solution and those obtained by using the finite element method.     

缝洞型油藏三维离散缝洞数值试井模型

缝洞型碳酸盐岩油藏发育着大尺度的溶洞和裂缝,非均质性极强,缝洞型碳酸盐岩油藏问题的研究成为了世界级难题之一.由于大尺度溶洞和裂缝对储层的流体流动起主导作用,因此,基于连续介质理论的双重介质或三重介质模型已不适合其中流体流动的描述.根据大型缝洞分布地质特征,探索性地提出了一种板块组合的复合架构离散缝洞模型描述该类油藏中的流体流动,将三维空间大裂缝用板块描述,溶洞用高渗透率和高孔隙度不规则多面体团块描述.将裂缝面用二维三角形单元离散, 溶洞和基质用三维四面体离散, 利用三维混合单元有限元法对建立的不定常渗流模型进行求解,得到了三维渗流条件下的试井理论曲线及压力场分布.通过对试井理论曲线特征的分析, 获得了各敏感参数对试井曲线的影响规律.通过对1口井的实际测试资料解释结果的分析,并与实际地震反射资料及生产实际资料的对比,发现本文所建立的模型可以较好地反映裂缝和溶洞的地质动态状况,并与实际生产状况具有较好的一致性.这一结果说明了所建模型的正确性以及测试资料分析结果的可靠性.      

考虑诱导缝多段压裂水平井非均质改造压力动态模型研究

为了准确模拟致密油藏水平井大规模压裂形成复杂裂缝网络系统和非均质储层井底压力变化,建立考虑诱导缝矩形非均质储层多段压裂水平井不稳定渗流数学模型,耦合裂缝模型与储层模型得到有限导流裂缝拉普拉斯空间井底压力解,对两种非均质储层模型分别利用数值解、边界元和已有模型验证其准确性.基于压力导数曲线特征进行流动阶段划分和参数敏感性分析,得到以下结果:和常规压裂水平井井底压力导数曲线相比较,理想模式下,考虑诱导缝影响时特有的流动阶段是综合线性流阶段、诱导缝向压裂裂缝“补充”阶段、储层线性流动阶段和拟边界控制流阶段.诱导缝条数的增加加剧了综合线性流阶段的持续时间,降低了流体渗流阻力,早期阶段压力曲线越低;当诱导缝与压裂裂缝导流能力一定时,裂缝导流能力越大,线性流持续时间越长;当所有压裂裂缝不在一个区域时,沿井筒方向两端区域低渗透率弱化了低渗区域诱导缝流体向压裂裂缝“补充”阶段,因此,沿井筒方向两端区域渗透率越低,早期阶段压力曲线越高;当所有压裂裂缝在一个区域时,渗透率变化只影响径向流阶段之后压力曲线形态,外区渗透率越低,早期径向流阶段之后压力曲线越高.通过实例验证,表明该模型和方法的实用性和准确性.     

沥青混合料三点弯曲断裂的扩展有限元模拟

利用随机骨料生成和投放算法建立沥青混合料两相异质模型,借助ABAQUS软件开展了分别在跨中、偏中20mm和40mm不同位置包含预切口的沥青混合料三点弯曲试件断裂行为的扩展有限元模拟,分析了预切口位置对裂纹扩展路径、载荷位移曲线等的影响规律,并与实验结果进行了比较。结果表明,预切口偏离跨中越远,裂纹扩展方向偏离预切口延长线的角度越大,试件的峰值载荷越大,对试件承载能力的削弱越弱。实验和数值模拟结果在定性上的一致性表明,利用扩展有限元方法模拟沥青混合料断裂行为是可行的。     

页岩气藏压裂水平井试井分析

页岩气藏资源丰富,开发潜力巨大,已成为目前研究的热点.与常规气藏相比,页岩气藏运移机制复杂,流动模式呈非线性,有必要考虑页岩气的吸附解吸,天然微裂缝的应力敏感性,人工裂缝内的非达西流等非线性因素对压裂水平井压力响应的影响. 基于双重介质和离散裂缝混合模型,分别采用Langmuir等温吸附方程描述吸附解吸,渗透率指数模型描述应力敏感,Forchheimer方程描述非达西效应,建立页岩气藏压裂水平井数值试井模型. 运用伽辽金有限元法对模型进行求解.根据试井特征曲线,划分流动阶段,着重分析非线性因素对压力响应的影响.结果表明:页岩气藏压裂水平井存在压裂裂缝线性流、压裂裂缝径向流、地层线性流、系统径向流及封闭边界影响5 种流动阶段.吸附解吸的影响发生窜流之后,Langmuir吸附体积增大,拟压力导数曲线凹槽更加明显,系统径向流出现时间与压力波传播到边界时间均延迟;天然裂缝系统的应力敏感性主要影响试井曲线的晚期段,拟压力和拟压力导数曲线均表现为上翘,应力敏感效应越强,上翘幅度越大;高速非达西效应对早期段影响较大,非达西效应越强,拟压力降幅度越大,试井曲线上翘.与解析解的对比以及矿场实例验证了模型的正确性与适用性.      

水力压裂中裂纹扩展的数值模拟

水力压裂是在高压粘滞流体或清水作用下地层内裂缝起裂与扩展的过程。由于包含岩石断裂和流-固耦合等复杂问题,对该过程的数值模拟具有相当大的挑战性。本文建立基于有限元与离散元混合方法的裂纹模型,模拟岩石裂纹扩展,实现了连续向非连续的转化;建立双重介质流动模型,裂隙流作为孔隙渗流的压力边界,孔隙渗流反作用裂隙的压力求解,处理了流体在基岩与人工裂缝中的协调流动;将裂纹模型与流体流动模式进行结合,建立断裂-应力-渗流耦合形式的力学模型,进一步分析了水力压裂的基本过程,综合多种数值计算方法,编写程序,在验证岩体裂纹模型与双重介质流动模型有效性的基础上,对压裂过程进行复现,将模拟结果与文献结果进行了对比,并讨论了所构建模型的优缺点。     

广义扩展有限元法及其在裂纹扩展分析中的应用

结合广义有限元法(GFEM)和扩展有限元法(XFEM)的特点,提出了一种新的数值方法——广义扩展有限元法(GXFEM)。阐述了广义扩展有限元法的基本原理,对相关公式进行推导,探讨数值实施中需注意的重要问题,给出利用广义扩展有限元法进行断裂分析时应力强度因子的计算方法,编写了广义扩展有限元法程序。通过算例进行了应力强度因子的计算,模拟了结构裂纹的扩展过程。算例结果表明,利用广义扩展有限元法计算裂纹扩展问题,不需要进行过密的网格划分,且网格在裂纹扩展后无需重新剖分,具有相当高的计算精度。