基于近场动力学-有限元耦合方法的正交各向异性多孔弹性介质中流体驱动压裂的数值模拟

  • A+

本文推广了近场动力学-有限元耦合方法来模拟正交各向异性多孔弹性介质中的流体驱动压裂。为此,本文提出了一个键微模量和临界能量密度随键取向连续变化的正交各向异性非普通态近场动力学( NOSBPD )模型。推导了包含用于捕捉固体裂缝的正交各向异性NOSBPD模型和用于模拟具有各向异性渗透张量的流体流动的有限元模型( FEM )的耦合系统的控制方程。为了确定所提出的正交各向异性NOSBPD模型合适的离散参数,进行了m -收敛和Δ-趋同研究。为了验证模型,模拟了不同材料取向的紧凑拉伸试验中的断裂扩展。弹性和渗透率各向异性、注入速度以及多裂缝和天然裂缝的存在对正交各向异性孔隙弹性介质中流体驱动的裂缝扩展有显著影响。在不同的材料角度下,断裂似乎偏离初始缺口方向,并且这种趋势随着弹性模量和垂直于弯曲平面的断裂释放率的降低而变得更加突出。渗透率各向异性影响流体模式和裂缝扩展长度。较快的注入速度导致更宽的损伤区和更弥散的孔隙压力场。

基于近场动力学-有限元耦合方法的正交各向异性多孔弹性介质中流体驱动压裂的数值模拟
图1 正交各向异性多孔弹性介质中流体驱动压裂的半对称模型,用于三种情况:(a)工况I:单一初始切口;( b )工况II:两个平行的初始缺口;(c)工况III:含竖向天然裂缝
基于近场动力学-有限元耦合方法的正交各向异性多孔弹性介质中流体驱动压裂的数值模拟
图2 工况I中不同材料方向下半对称模型预测的裂纹形态:(a) t = 8.10 s时刻θ = 0 °;(b) t = 10.05 s时刻θ = 15 °;(c)在时间t = 10.03秒时,θ = 30 °;(d) t = 13.85 s时θ = 45 °;(e) t = 13.65 s时刻θ = 60 °;(f)在时间t = 13.00 s时,θ = 75 °;(g)在时间t = 13.00 s时,θ = 90
基于近场动力学-有限元耦合方法的正交各向异性多孔弹性介质中流体驱动压裂的数值模拟
图3 场景I中不同材料方向下半对称模型中的孔压分布
基于近场动力学-有限元耦合方法的正交各向异性多孔弹性介质中流体驱动压裂的数值模拟
图4 裂纹模式和相应的孔压分布:( a ) ( d ) t = 8.35 s时材料角θ = 0 °,( b ) ( e ) t = 15.00 s时θ = 45 °,( c ) ( f ) t = 10.00 s时θ = 90 °;( g ) ( j )在t = 24.00 s时θ = 0,( h ) ( k )在t = 24.00 s时θ = 45,( i ) ( l )在t = 29.50 s时θ = 90°
weinxin
我的微信公共号
我的微信公招扫一扫

发表评论

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: