一种基于PFC离散元的页岩水力压裂数值模拟方法技术

本发明专利技术公开了一种基于PFC离散元的页岩水力压裂数值模拟方法，涉及页岩气藏水力压裂技术领域；基于PFC离散元方法探索水力压裂裂缝扩展机理，将渗流与离散元进行耦合，模拟页岩储层水力压裂裂缝起裂机制，进而指导页岩储层压裂设计；明确页岩水力压裂起裂机理和裂缝最终展布形态，揭示水力裂缝与页岩交互作用问题，实现页岩储层水力压裂最终裂缝展布形态的精准预测；从微观尺度精确表征水力裂缝起裂和扩展机理，准确捕捉水力裂缝宏观展布形态；精确再现水力裂缝与页岩交互作用，使得到的模拟结果与实际情况更加吻合。结果与实际情况更加吻合。结果与实际情况更加吻合。

【技术实现步骤摘要】
一种基于PFC离散元的页岩水力压裂数值模拟方法


[0001]本专利技术涉及页岩气藏水力压裂
，具体涉及一种基于PFC离散元的页岩水力压裂数值模拟方法。

技术介绍

[0002]我国页岩气储量巨大，探明可采页岩气资源量高达31.6万亿m3，资源潜力巨大。页岩气藏是一个“自生、自储”系统，页岩储层物性极差，多为纳米孔喉，孔隙度常为1％
‑
7％，基质渗透率较小。国内外多采用水平井分段多簇压裂技术对页岩储层进行整体压裂，从而改善储层渗流能力，增大储层改造体积，增加页岩气产量。水力压裂作为有效开采页岩气的核心技术之一，广泛应用于石油天然气等行业，通常向页岩储层注入足够压力的压裂液，页岩断裂产生大量的水力裂缝，形成气液高速渗流通道，以改造页岩气储层，达到增产的目的。从水力压裂技术被开发至今，最受关注的是水力压裂作用下岩石的破裂机理和水力裂缝产生的方向和大小。尽管室内实验能够探索裂缝起裂规律，但难以精确表征岩石力学行为和水力压裂过程中的流动
‑
位移
‑
应力耦合问题。同时，数值模拟具有成本低、可重复及试验周期较短的优势，因而得到了更广泛的应用，大量工程实际的检验也验证了数值模拟的准确性，应用前景十分广阔。
[0003]由于页岩储层的复杂性，在解决水力压裂问题时，常规的有限元法、块体离散元法等数值模拟方法表现出一定的局限性。近些年来，基于离散元的数值模拟方法在处理岩土工程的问题上得到了广泛的应用。离散元方法与其他数值模拟方法对比，其最大的优势就是离散元数值模拟不受物体变形量的影响，针对多相介质耦合作用下的力学行为和岩体的破裂等问题更具有仿真性，能够直观的看出也岩体的变化过程和结果。PFC离散元法主要是通过圆盘与圆盘之间的相互作用来模拟页岩中颗粒单元之间的力学行为。PFC程序是在基于离散元理论的基础上通过显式差分的运算法则综合研发得到的，PFC离散元法与普通连续介质的数值模拟方法对比，最大的优点就是能够充分体现出岩体的实体结构及其内部颗粒单元之间的相互作用，这种以颗粒单元为基本研究对象，在针对页岩颗粒单元力学行为的观测与研究时更为直观，是一种运用非连续的计算方法处理复杂变形问题的基本理论。正是由于PFC离散元是以颗粒单元为基本的计算单元，以牛顿第二定律为基本的力学关系，因此在PFC离散元数值模拟中对颗粒单元性质的描述更能表现出研究对象的本质性。

技术实现思路

[0004]为解决现有技术存在的不足，本专利技术提出了一种基于PFC离散元的页岩水力压裂数值模拟方法，基于PFC离散元方法探索水力压裂裂缝扩展机理，将渗流与离散元进行耦合，模拟页岩储层水力压裂裂缝起裂机制，进而指导页岩储层压裂设计；明确页岩水力压裂起裂机理和裂缝最终展布形态，揭示水力裂缝与页岩交互作用问题，实现页岩储层水力压裂最终裂缝展布形态的精准预测。
[0005]采用的技术方案为：
[0006]一种基于PFC离散元的页岩水力压裂数值模拟方法，包括以下步骤：
[0007]步骤一：建立外部墙体单元，在墙体内生成颗粒单元，在中心位置构建圆形区域，删除颗粒单元，建立注液井眼；
[0008]根据单轴压缩实验结果、巴西劈裂实验结果与单轴压缩模拟结果、巴西劈裂模拟结果对比分析输出的应力应变曲线，标定页岩细观参数；
[0009]赋予颗粒单元基本物性参数，通过平行粘结单元粘结相邻的颗粒单元，生成页岩储层离散元数值模型；
[0010]所述页岩细观参数包括所有颗粒单元基本物性参数及平行粘结单元参数；
[0011]步骤二：设定相邻颗粒单元围成的区域为流体域，设定颗粒单元接触垂直处为流体管道，在页岩储层离散元数值模型上建立网络流动模型；
[0012]步骤三：在页岩边界墙体施加围压，设定井眼注入流体的速度，进行渗流
‑
位移
‑
应力耦合的数值模拟计算；
[0013]步骤四：实时监测颗粒单元间粘结状态，判断粘结单元是否发生破坏，当粘结单元所承受应力大于粘结单元的法向强度或切向强度时，粘结单元为破坏状态，视为产生微裂缝；计算裂缝数量及水力裂缝扩展形态；计算结束后，导出并保存页岩储层水力裂缝扩展形态模拟结果。
[0014]所述流体域的变量计算如下：
[0015]在二维模型中，通过立方定律得到平行板通道的流量，如式(1)所示：
[0016][0017]其中，w为通道的开度，单位m；ΔP为相邻流体域单元之间的压力差，单位Pa；L为通道的长度，单位m；μ为流体的粘度，Pa
·
s；
[0018]流体域单元中每一个时步的压力变化，如式(2)所示：
[0019][0020]其中，K
f
为流体体积模数，单位Pa；V
d
为流体域单元的表观体积，单位m3；Δt为时间步长，单位s。
[0021]所述颗粒单元，粘结后的法向作用力、切向作用力和弯矩通过平行粘结模型计算；
[0022]当法向接触力为拉力时，如式(3)所示：
[0023]w＝w0+ng
ꢀꢀꢀ
(3)
[0024]其中，w0为颗粒间接触力为0时对应的通道开度，称为残余开度，n为距离压缩因子，g为相邻颗粒表面的法向距离；
[0025]当法向接触力为压力时,如式(4)所示：
[0026][0027]其中，F0为通道开度是时颗粒单元间接触力的大小，F为颗粒单元之间的接触力，根据颗粒单元之间的叠合量，按力
‑
位移定理得出。
[0028]所述力
‑
位移定理对于“颗粒单元
‑
颗粒单元”接触，单位法向量n
i
的定义，如式(5)所示：
[0029][0030]其中，和分别为颗粒单元1和颗粒单元2的位置坐标，d为两球球心之间的距离；在“颗粒单元
‑
墙体”接触中，单位法向量n
i
则定义为球心距离墙体最大直线的方向；
[0031]“颗粒单元
‑
颗粒单元”或“颗粒单元
‑
墙体”之间的“重叠”量，如式(6)所示：
[0032][0033]颗粒单元
‑
颗粒单元或颗粒单元
‑
墙体接触产生的接触力，分解为法向接触力和切向接触力如式(7)所示：
[0034][0035]其中，法向接触力计算，如式(8)所示：
[0036][0037]其中，K
n
为接触处的法向刚度，是一个割线模量，使总的法向位移与法向接触力发生联系；
[0038]接触力矢量F
i
产生的切向接触力以增量的形式计算，当颗粒单元间的接触形成时，总的切向接触力先初始化为零，之后颗粒单元间的每个相对位移增量引起的弹性切向增量会累加到当前值上，从而产生切向接触力，通过更新每个时间步长Δt内的单位法向量n
i
接触中心点x
C
的坐标来对接触运动，新的时间本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于PFC离散元的页岩水力压裂数值模拟方法，其特征在于：包括以下步骤：步骤一：建立外部墙体单元，在墙体内生成颗粒单元，在中心位置构建圆形区域，删除颗粒单元，建立注液井眼；根据单轴压缩实验结果、巴西劈裂实验结果与单轴压缩模拟结果、巴西劈裂模拟结果对比分析输出的应力应变曲线，标定页岩细观参数；赋予颗粒单元基本物性参数，通过平行粘结单元粘结相邻的颗粒单元，生成页岩储层离散元数值模型；步骤二：设定相邻颗粒单元围成的区域为流体域，设定颗粒单元接触垂直处为流体管道，在页岩储层离散元数值模型上建立网络流动模型；步骤三：在页岩边界墙体施加围压，设定井眼注入流体的速度，进行渗流
‑
位移
‑
应力耦合的数值模拟计算；步骤四：实时监测颗粒单元间粘结状态，判断粘结单元是否发生破坏，当粘结单元所承受应力大于粘结单元的法向强度或切向强度时，粘结单元为破坏状态，视为产生微裂缝；计算裂缝数量及水力裂缝扩展形态；计算结束后，导出并保存页岩储层水力裂缝扩展形态模拟结果。2.如权利要求1所述的基于PFC离散元的页岩水力压裂数值模拟方法，其特征在于：所述页岩细观参数包括所有颗粒单元基本物性参数及平行粘结单元参数。3.如权利要求1所述的基于PFC离散元的页岩水力压裂数值模拟方法，其特征在于：所述流体域的变量计算如下：在二维模型中，通过立方定律得到平行板通道的流量，如式(1)所示：其中，w为通道的开度，单位m；ΔP为相邻流体域单元之间的压力差，单位Pa；L为通道的长度，单位m；μ为流体的粘度，单位Pa
·
s；流体域单元中每一个时步的压力变化，如式(2)所示：其中，K
f
为流体体积模数，单位Pa；V
d
为流体域单元的表观体积，单位m3；Δt为时间步长，单位s。4.如权利要求1所述的基于PFC离散元的页岩水力压裂数值模拟方法，其特征在于：所述颗粒单元，粘结后的法向作用力、切向作用力和弯矩通过平行粘结模型计算；当法向接触力为拉力时，如式(3)所示：w＝w0+ng
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)其中，w0为颗粒间接触力为0时对应的通道开度，称为残余开...

【专利技术属性】
技术研发人员：张永利，孙爽，程瑶，
申请(专利权)人：辽宁工程技术大学，
类型：发明
国别省市：