非常规储层中常常发育大量天然裂缝,对水力压裂裂缝扩展形态影响显著。为了研究裂缝性储层中水力裂缝的扩展规律,往往要借助数值模拟。现有的水力压裂数值模拟方法主要有扩展有限元、离散元和边界元。分别介绍了上述3种方法的基础理论、优缺点以及应用情况。通过比较各种压裂数值模拟计算方法,在研究不同问题时要根据问题本身以及研究需求选用合适的数值方法,进而掌握水力压裂过程中的裂缝扩展规律,指导实际的水力压裂施工设计。 A large number of natural fractures were often developed in unconventional reservoirs; it had a significant effect on the fracture propagation morphology in hydraulic fracturing. Fractures were rich in unconventional reservoirs and these fractures will affect hydraulic fracturing significantly. Therefore, numerical simulation was an effective method to study fracture propagation in fractured reservoirs. Till now, the numerical simulation methods mainly included extended finite element method, distinct element method and boundary element method. In this paper, the basic theory, the advantages and disadvantages, and the applications of these three methods are all studied. By comparing these three methods, it is found that proper numerical methods are adopted in consideration of different research requirements. Furthermore, the fracture propagation law in hydraulic fracturing process can be grasped and the actual hydraulic fracturing operation design can be guided.
杨向同1,张杨1,范文同1,巴旦1,朴立文2,3,孙一流2,3
1中国石油塔里木油田分公司油气工程研究院,新疆 库尔勒
2油气资源与探测国家重点实验室(中国石油大学(北京)),北京
3中国石油大学(北京)石油工程学院,北京
作者简介:杨向同(1972-),男,高级工程师,主要从事高温高压井完井改造方面的工作。
收稿日期:2018年1月15日;录用日期:2018年2月25日;发布日期:2018年6月15日
非常规储层中常常发育大量天然裂缝,对水力压裂裂缝扩展形态影响显著。为了研究裂缝性储层中水力裂缝的扩展规律,往往要借助数值模拟。现有的水力压裂数值模拟方法主要有扩展有限元、离散元和边界元。分别介绍了上述3种方法的基础理论、优缺点以及应用情况。通过比较各种压裂数值模拟计算方法,在研究不同问题时要根据问题本身以及研究需求选用合适的数值方法,进而掌握水力压裂过程中的裂缝扩展规律,指导实际的水力压裂施工设计。
关键词 :裂缝性储层,水力压裂,扩展有限元,离散元,边界元
Copyright © 2018 by authors, Yangtze University and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
以页岩气、页岩油、致密气为主的非常规油气资源在新开发储层中占据了越来越大的比重,随着大规模储层改造技术日益成熟,水平井分段压裂技术在现场广泛应用,并取得了巨大的经济效益。与常规均质砂岩储层不同,这类储层中往往发育有大量的天然裂缝 [
地震技术、岩心切片分析技术及井下成像技术均已证明,非常规储层中发育有丰富的天然裂缝,与低胶结强度的沉积岩层理面共同形成了储层中的弱交界面。大量的室内岩石力学试验表明,储层非均质性及各向异性对水力裂缝的扩展及裂缝形态具有很大影响。
裂缝性储层各向异性的特点使得人工裂缝尖端应力场与常规储层扩展相比更为复杂。室内试验已证明,地应力分布、交界面强度、弱面方位角均是影响人工裂缝穿透天然裂缝行为的重要因素。为了判别水力压裂过程中,水力缝尖扩展与天然裂缝的相互影响,国内外学者提出了诸多判别准则及研究结果,如Warpinski等 [
裂缝性储层水力压裂数值模拟技术是模拟水平井分段压裂施工条件下,人工裂缝与天然裂缝相互影响形成的裂缝网络的逢高、缝宽等几何参数以及支撑剂分布情况。其目的是为了设计和优化施工程序,为后续油田整体开发提供理论依据。水力压裂模拟可以概括为4个过程的耦合:岩石在高压流体作用下发生断裂,形成流动通道;压裂液在裂缝内流动,将井底压力传递到缝尖;流体在孔隙介质中的渗流;支撑剂随压裂液流动分布。由于研究的问题具有复杂性,一般采用数值方法求解。目前,关于裂缝性储层中水力裂缝扩展数值模拟方面的研究,国内外学者已经做了大量的工作,取得了较大的研究进展。主要数值模拟方法包括:扩展有限元法,离散元法以及边界元法。
扩展有限单元法是有限单元法的一种变形。当研究目标区域含有裂缝时,有限元法在处理裂缝尖端的奇异性时准确性较差,同时在裂缝扩展过程中,有限元法需要实时变化整个区域的网格划分,导致计算速度缓慢。扩展有限元法通过引入带有不连续性质的局部加强函数,能够有效且准确地描述裂缝的不连续特征,同时在处理裂缝扩展问题时,不需要根据裂缝扩展路径重新构建整个区域的网格,在很大程度上提高了计算效率,同时还能更准确地反映裂缝的扩展规律(图1)。
图1. 扩展有限元法中的不同单元类型
王涛等 [
离散元方法的基本思想是将岩体分割为若干独立的单元,邻近的单元可以接触也可以分开。单元之间的相互作用力根据单元之间的重叠量,按力-位移定理求出,而单元的运动则根据该单元所受的不平衡力和不平衡力矩的大小按牛顿第二定律求出。按上述步骤循环计算,直到每一个单元均达到平衡状态。PFC颗粒流方法是将每一个独立单元视为规则球体的离散元方法。在颗粒流方法中,本构关系遵循力-位移关系,运动方程遵循牛顿第二定律。
离散元法 [
Hazzard等 [
边界元法分为直接边界元和间接边界元,目前常用的模拟裂缝扩展的边界元方法主要为间接边界元法中的位移不连续法。该方法是基于线弹性力学推导的,适用于处理包含裂缝等不连续区域的一种常用方法。尤其在处理水力压裂问题时,位移不连续法能够用简单、直接的形式描述复杂裂缝网络,同时考虑多裂缝间的相互干扰,因此被国内外的大量学者广泛采用。
边界元在一定程度上克服了由积分奇异性造成的困难,并应用到了很多工程及科学领域。其主要特点为: ① 将区域问题化简为边界问题,显著降低了需要求解的自由度数,提高了计算效率; ② 能够采用解析解与数值解相结合的方式,利用已知基本解作为边界积分方程的基函数,提高了计算精度; ③ 已知微分算子的基本解自动满足无限远处条件,因此边界元适合于处理无限域及半无限域问题。但边界元法求解以存在相应的微分算子的基本解为前提,一定程度上限制了其使用范围和发展前景。
Olson等 [
扩展有限元法在处理裂缝扩展问题时其裂缝扩展路径不受单元网格划分的影响,但是需要对整个研究目标进行网格划分,网格数目巨大,计算耗时长。
边界元法通过降维处理,能够将区域问题转化成边界问题,大大减少了计算量,但是边界元问题不能考虑研究对象内部的非均匀性和各向异性,同时随着计算次数的增多,边界元计算量逐渐增大,严重影响计算速度。
以PFC3D为代表的离散元软件在研究小尺寸试样的细观力学行为上具有一定的优势,但在模拟油田尺度的水力裂缝扩展时,由于涉及到单元众多,并不适用。
裂缝性储层中发育着大量的天然裂缝,显著影响着水力压裂裂缝扩展轨迹。而水力压裂数值模拟技术是研究裂缝扩展的有力工具,同时也是模拟水力压裂过程中裂缝扩展规律的重要方法。现有的水力压裂数值模拟方法主要是扩展有限元法、离散元法和边界元法,每种方法各有优缺点,在使用时要结合具体的研究对象和主要需要解决的问题进行筛选。
中国石油塔里木油田分公司项目(201015120220)。
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[
[编辑] 帅群