PFC离散元数值模拟深度解析——围压卸载条件下砂岩损伤及渗透演化颗粒流模拟

一、引言与痛点直击

在深部地下工程开采过程中，围岩常处于高应力状态。随着开采进行，围压卸载导致岩体应力状态发生显著变化，进而引发岩体损伤和渗透特性改变。围压卸载条件下砂岩的损伤演化及渗透规律研究，对于深部工程稳定性评价和突水灾害预测具有重要意义。

图1 砂岩试样

▎ 围压卸载是深部工程岩体的典型受力状态

▎ 岩体损伤与渗透演化存在密切耦合关系

▎ PFC离散元方法可模拟岩体破裂和流体运移

▎ 流固耦合算法是实现损伤-渗透模拟的关键

二、研究概况与实验方法

本研究采用室内三轴试验与PFC离散元数值模拟相结合的方法，系统研究了围压卸载条件下砂岩的损伤演化规律和渗透特性变化。试验考虑了不同围压水平（10MPa、20MPa、30MPa）对砂岩力学行为和渗透特性的影响。

▶ 2.1 试样与测试设备

试验采用标准砂岩试样，通过孔隙度测试和纵横波速度测试获取岩样的基本物理力学参数。

图2 孔隙度和纵横波速度测试装置

孔隙度测试采用氦气孔隙度仪，纵横波速度测试采用超声波脉冲透射法。这些基础参数为后续数值模拟的细观参数标定提供了依据。

▶ 2.2 试验装置

试验采用MTS815岩石力学试验系统，配备围压加载装置和渗透压控制系统。试验过程中通过声发射监测系统实时记录岩体内部的微破裂活动。

图3 试验装置示意图

━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━

三、PFC离散元数值模拟

为了深入理解围压卸载条件下砂岩的损伤演化及渗透规律，本研究采用PFC2D离散元软件建立了流固耦合数值模型。PFC通过圆形颗粒和粘结模型模拟岩石材料，能够很好地再现岩体的破裂过程和流体运移。

▶ 3.1 流固耦合算法

数值模拟采用PFC2D内置的流固耦合算法，通过管域模型描述颗粒间的流体通道，实现应力场与渗流场的耦合计算。

图4 流固耦合算法模型

管域模型将颗粒间的孔隙抽象为管道网络，流体在管道中流动遵循达西定律。通过计算每个时步的渗透压分布，实现流体与固体颗粒的相互作用。

图5 PFC2D渗流计算过程渗透压分布

PFC2D渗流计算过程展示了渗透压从初始状态到稳定分布的演化过程。经过多次计算循环，渗透压分布逐渐趋于稳定，形成线性分布。

▶ 3.2 模型验证

数值模拟结果与室内试验结果进行了对比验证。通过调整细观参数，使数值模拟的应力-应变曲线与试验曲线吻合良好。

图6 数值模拟与试验应力-应变曲线对比

从对比结果可以看出，数值模拟能够较好地再现砂岩在围压卸载条件下的应力-应变关系，包括弹性阶段、屈服阶段和破坏阶段。不同围压下的模拟曲线与试验曲线趋势一致，验证了模型的可靠性。

四、损伤演化分析

通过PFC数值模拟，详细分析了围压卸载条件下砂岩的损伤演化过程，包括裂纹萌生、扩展和贯通的整个过程。

▶ 4.1 声发射特征

数值模拟中的声发射事件数反映了岩体内部微破裂的活跃程度。通过监测声发射事件数的变化，可以识别岩体损伤的不同阶段。

图7 偏应力与AE事件数关系曲线

偏应力与AE事件数关系曲线显示，在弹性阶段（I区）声发射事件数较少；进入屈服阶段（II区）后，声发射事件数开始增加；在破坏阶段（III区），声发射事件数急剧增加，表明岩体内部微破裂大量产生；在残余阶段（IV区），声发射事件数逐渐减少。

▶ 4.2 裂纹演化规律

裂纹数是表征岩体损伤程度的重要指标。数值模拟追踪了不同围压下裂纹数随轴向应变的演化规律。

图8 裂纹数与轴向应变关系曲线

裂纹数与轴向应变关系曲线显示，在卸荷点之前裂纹数增长缓慢；卸荷之后，裂纹数快速增长。围压越低，裂纹数增长越快，表明低围压条件下岩体更容易发生损伤。

图9 砂岩试样损伤演化过程

砂岩试样损伤演化图展示了从初始状态到最终破坏的裂纹分布变化。从图中可以看出，裂纹首先在某些薄弱区域萌生，然后逐渐扩展、连接，最终形成贯穿性破裂带。

▶ 4.3 接触力链分析

接触力链是颗粒材料内部力传递的主要途径。通过分析接触力链的演化，可以揭示岩体损伤的力学机制。

图10 配位数与接触力链分布

配位数与轴向应变关系曲线显示，随着损伤发展，配位数逐渐降低，表明颗粒间的接触数目减少。图中嵌入的接触力链云图展示了不同应变阶段的力链分布，可以清晰看到力链从均匀分布向局部集中的转变过程。

五、能量演化分析

能量分析是理解岩体损伤机制的重要手段。通过追踪总输入能量、弹性应变能和耗散能的演化，可以定量描述岩体损伤过程中的能量转换。

▶ 5.1 能量演化曲线

数值模拟记录了加载过程中各类能量的演化。总输入能量是外部做功的累积，弹性应变能是储存在岩体内部的能量，耗散能是用于产生裂纹和摩擦消耗的能量。

图11 能量演化曲线

能量演化曲线显示，在弹性阶段，大部分输入能量转化为弹性应变能储存起来；进入屈服阶段后，耗散能开始明显增加；在破坏阶段，耗散能急剧增加，弹性应变能快速释放。图中A、B、C三点分别对应屈服起始点、峰值点和残余阶段。

▶ 5.2 能量分布特征

不同围压下，能量的分布特征存在显著差异。高围压条件下，岩体能够储存更多的弹性应变能，破坏时释放的能量也更大。

图12 能量分布对比

能量分布对比图显示，随着围压增加，总输入能量和弹性应变能均呈线性增加趋势。弹性应变能占比保持在85-90%之间，耗散能占比约为10-15%。这表明围压越高，岩体储能能力越强，但也意味着破坏时可能释放更多的能量。

六、渗透演化分析

渗透演化是围压卸载条件下砂岩损伤的重要表现。通过流固耦合数值模拟，详细分析了渗透率随损伤发展的演化规律。

▶ 6.1 裂纹与渗流的关系

裂纹的产生和扩展为流体运移提供了通道。通过分析裂纹区域面积比（CAR）和渗流区域面积比（KAR）的关系，可以建立损伤与渗透的定量联系。

图13 应力云图演化

应力云图演化展示了围压卸载过程中试样内部应力的重新分布。随着围压卸载，应力集中区域逐渐发展，最终导致裂纹萌生和扩展。

图14 裂纹与渗流演化过程

裂纹与渗流演化过程图展示了裂纹扩展与渗透率分布的对应关系。红色区域表示高渗透率区域，与裂纹分布高度吻合，表明裂纹是流体运移的主要通道。

图15 裂纹与渗流区域面积比演化

裂纹与渗流区域面积比演化曲线显示，CAR和KAR都随轴向应变增加而增加，但KAR的增长速度略低于CAR，表明并非所有裂纹都参与渗流。

图16 渗流与裂纹区域面积比关系

渗流与裂纹区域面积比关系曲线显示，KAR与CAR呈非线性关系。当CAR较小时，KAR增长较慢；当CAR超过一定阈值后，KAR快速增加，表明渗流网络开始贯通。

▶ 6.2 渗透率演化规律

渗透率是表征岩体渗流能力的重要参数。数值模拟追踪了不同围压下渗透率随轴向应变的演化规律。

图17 不同围压下渗透率云图

不同围压下渗透率云图显示，10MPa围压条件下，试样内部形成了明显的优势渗流通道，渗透率分布不均匀；随着围压增加，渗透率分布趋于均匀，优势通道效应减弱。

图18 渗透率分布特征

渗透率分布特征进一步展示了围压对渗流模式的影响。高围压抑制了裂纹的扩展，使得渗流通道更加分散。

图19 渗透率演化规律

渗透率演化规律曲线显示，在初始阶段，渗透率略有下降，这是由于岩体压密所致；随着裂纹萌生和扩展，渗透率开始增加；进入破坏阶段后，渗透率急剧增加，增幅可达1-2个数量级。

七、工程应用建议

基于本研究的试验和数值模拟结果，为深部工程围岩稳定性评价和突水灾害预测提出以下建议。

通过本研究，我们深入理解了围压卸载条件下砂岩的损伤演化及渗透规律。PFC离散元流固耦合数值模拟方法，为研究这类复杂问题提供了有效的技术手段。研究成果可为深部工程围岩稳定性评价和突水灾害预测提供科学依据。如果您在科研或工程中遇到相关问题，欢迎咨询交流！

如有其他需要，欢迎关注我的咸鱼号:Orangetu

本推送内容基于公开文献进行梳理与解读，著作权归原作者及出版机构所有。内容仅代表个人对原文的理解与归纳，旨在学术交流与知识传播，不作商业用途。如需引用或进一步探讨，请务必参阅原始文献，若涉及版权事宜请联系我们处理。

感谢您的关注，麻烦多多对推文点赞、收藏和转发，并衷心希望您多多指教！帮助CAE数值仿真模拟进步提升！