三维离散元在煤层水平井井壁稳定的应用.pdfVIP

三维离散元在煤层水平井井壁稳定的应用 屈 平12 申瑞臣1 杨恒林1 张学磊2 付利1’2 摘要：为了阐明煤层中端割理和面割理等不连续面对井壁稳定的影响，利用三维离散元来分 析煤层水平井井壁稳定，实现了在三维空问中对不连续面的仿真，在充分、真实反映煤岩割理特性 的基础上，对煤层水平井进行了井壁稳定数值模拟分析。提出了虚拟节理和内外双层建模技术， 虚拟节理技术可实现节理在各个方向的随机分布和动态的断裂分析，内外双层技术可降低计算量 和消除边界的影响。模拟了沁水盆地3号煤层面割理和端割理的实际分布情况，最大程度反映了 煤层的真实微观特征。通过分析井眼的走向、井径、井底压力、割理密度对井壁稳定的影响，发现 当面割理的走向与井眼的走向相互垂直时，井壁最稳定；当两者相互平行时，井壁最不稳定；割理 密度越大，井壁的稳定性越差；井径越大，井壁越不稳定；当不考虑渗流时，井底压力越大井壁越 稳定。 关键词：3DEC煤层气离散元井壁稳定虚拟节理 由于煤岩的高度非连续特点，有限元等常规数值方法无法分析煤岩的井壁稳定，迫切需要 寻找新的数值模拟方法ll一…。而且煤岩的割理等微观结构对煤岩的力学行为影响较大，因此 在井壁稳定力学分析过程中必须充分、真实地考虑煤岩割理特性。李嗣贵‘1。把煤层视为节理 地层，张哲心1分析了井壁周围的诱导裂纹，都考虑了煤岩的非连续特点，但都没有考虑真实的 煤岩割理微观结构。笔者试图运用三维离散元，在充分、真实反映煤岩割理特性的基础上，针 对煤层水平井进行井壁稳定数值模拟分析。 1 三维离散元原理和力学模型 1．1基本原理 离散单元法不仅仅是一种数值计算方法，更是一种岩石力学的思维模式。它把岩体当作 各个离散的单元，单元之间的作用符合宏观牛顿定律，而单元内部适合微观的变形体弹塑性力 学。其实质是把结构面力学和岩石力学有机的统一，不连续介质和连续介质力学之间的有机 统一。离散单元法按空间维度，可分为二维和三维离散元；按研究对象可分为刚体和可变形体 离散元。离散单元法把岩体分为块体和不连续面(包括节理、断层等)两类，以块体为最小研 究单位即单元。不连续面一方面作为块体的边界，另一方面作为块体之问力作用的传媒。如 果是刚性块体，块体之间的作用力可以通过力和位移的关系求得，而个别块体的运动则完全根 据该单元所受的不平衡力和不平衡力矩的大小，遵循牛顿运动定律；如果是变形块体，块体之 间的作用力依然可以通过力和位移的关系求得，而个别块体的变形和运动则更加复杂。 作者简介：屈平，男，1982年生，2008年获中国石油勘探开发研究院硕士学位，现为中国石油勘探开发研究院在读博士 主要从事煤层气钻升技术研究。E—mail：qplol2@163．com。 一339— 1．2力学模型 1．2．1 不连续面的刚度模型 在离散单元中，块体与块体之间不需要遵守连续方程及变形协调方程，因此在运动过程 中，块体和块体之问可以分离、滑动和嵌入。此时，不连续面的刚度和强度将决定块体之间的 相互作用。 当块体嵌入的时候，刚度就起关键作用。刚度大则嵌入量就大，冈0度小则嵌入量小。在离 散元里，常用两个块体中问连接一弹簧来模拟刚度的作用。当接触刚度大，则单元之间的叠合 量大，反之亦然。 在三维空间里，不连续面为空间的任意形状的面，可以为曲面和平面，但一般简化为平面， 因此嵌入的时候接触区域不再是一个点而是一个面。其刚度模型为 f△F。=一矗。△u?／4。 { (1) l△F。=一惫。△“；A。 式中艮。和尼。——分别为节理面的法向刚度和切向刚度，MPa； A。——接触的面积，m2； △“?和△M；——分别为节理面法向和切向位移，m； △，。和△F。——分别为法向和切向作用力增量，N。 1．2．2块体运动和变形模型 对于刚性块体，运动完全遵循牛顿第二定律；对于可变形块体，块体一方面会运动，另一方 面会产生变形。因此刚性块体只需平动和转动两个方程就可以完全确定运动情