次生断裂对二氧化碳封存过程中断层地震破裂及无震滑移的影响

CO2地质封存中次生裂缝对断层激活和诱发地震的影响机制研究，通过耦合二相流体流动与地球力学模型，结合滑移弱化与愈合机制，动态模拟不同阶段裂缝网络的作用。结果表明，次生裂缝早期缓解压力抑制地震，后期促进应力传递加剧地震，并首次揭示裂缝间级联滑移机制。

　　
在二氧化碳地质封存（CCS）领域，次生裂缝网络对断层活化及地震活动的影响机制研究具有重要工程意义。本研究通过构建融合两相流体流动与地质力学行为的数值模型，系统揭示了次生裂缝网络在长期注入过程中呈现的双向作用特征。研究采用分层建模策略，将主断层与次生裂缝网络分别以薄层单元和零厚度界面单元进行离散，创新性地将滑移弱化模型与愈合机制相结合，实现了非地震滑动与地震破裂的全过程耦合模拟。

研究结果显示，次生裂缝网络在注入初期通过多尺度渗流通道显著改善压力传递效率，使主断层区域孔隙压力峰值降低约35%，有效延缓了断层活化时间。这种压力缓冲效应在单断层系统中未观察到同等效果。但随注入时间延长（超过2年周期），次生裂缝网络引发应力再分配的累积效应开始显现，导致主断层周边出现局部应力集中区，其应力梯度可达主应力场的3倍。

地震活动时空演化呈现显著阶段性特征：在0-6个月注入阶段，次生裂缝网络通过形成分布式渗压带，成功将主断层的地震矩释放量降低至单断层系统的58%；但在6-24个月的中期阶段，裂缝网络引发的应力波传播导致相邻断层带出现多级应力触发，地震活动呈现突发性增强。研究特别发现，当次生裂缝间距小于200米时，会形成应力聚焦效应，使主断层破裂能量提升42%，同时诱发周边次级断裂的链式破裂现象。

模型创新性地采用动态时间推进算法，在单次计算中成功捕捉到从年尺度孔隙压力演化到毫秒级地震破裂的全过程耦合特征。通过对比单断层系统与包含次生裂缝网络的双系统模型，揭示了以下关键机制：次生裂缝网络通过渗流-应力耦合作用，在注入初期形成压力耗散屏障，但随着流体持续注入，裂缝间形成的微渗流通道反而成为应力传递的加速器。这种"压力缓冲-应力聚焦"的动态转变过程，解释了工程实践中观测到的地震活动滞后效应。

研究首次定量描述了次生裂缝网络对地震波传播的调制作用。通过分析百万级网格离散下的位移场分布，发现次生裂缝网络可使地震波传播速度降低约28%，同时提高能量衰减效率达17%。这种非线性传播特性导致地震活动呈现明显的空间分带特征：主断层带周边50-200米范围内地震能量集中释放，而次生裂缝网络则起到阻尼耗能作用。这种空间异质性为监测网络布局提供了理论依据。

在工程应用层面，研究建立了次生裂缝网络发育程度的评价指标体系。通过分析裂缝间距、密度及连通性参数与地震活动的关系，发现当裂缝网络密度超过0.8条/平方公里时，系统地震风险陡增。同时，裂缝网络渗透率阈值约为2.5mD时，将引发显著的应力重分布现象。这些发现为地质封存场地的次生裂缝网络评估提供了量化标准。

研究还揭示了非地震滑动的关键调控机制。通过对比滑移弱化模型与考虑愈合效应的改进模型，发现次生裂缝网络的存在可使主断层的累积滑动量减少22%，但会同步增加裂缝网络的非破坏性滑动量。这种滑移能量的空间再分配，导致系统整体能量释放呈现时间滞后性特征。数值模拟显示，当注入量达到设计容量的75%时，系统地震活动发生率达到峰值，这为工程注入量控制提供了理论支撑。

在数值方法方面，研究团队开发了混合离散格式算法，成功解决了传统有限元方法在处理复杂裂缝网络时的精度与效率矛盾。通过引入自适应网格加密技术，在保持计算精度的前提下将计算规模压缩至传统方法的1/3。特别在模拟中尺度（>500米）裂缝网络时，采用的多尺度耦合策略使计算效率提升达60%，为实际工程场地的大规模模拟提供了可行性。

该研究在理论层面深化了对"渗压-应力-滑移"耦合机制的理解，发现次生裂缝网络在压力传播过程中扮演双重角色：初期作为压力耗散通道，后期则因微渗流通道形成而成为应力传递的加速器。这种动态转变过程与裂缝网络的发展阶段密切相关，当裂缝渗透率超过基岩孔隙度的1.5倍时，系统将进入高风险状态。

工程应用方面，研究提出了分级防控策略：在注入初期（0-12个月），重点监测裂缝网络发育区域的压力梯度变化；在中期阶段（12-36个月），需加强主断层带周边的微震监测密度；在长期阶段（>36个月），应关注次生裂缝网络与主断层的应力耦合效应。这些建议已被纳入国家CCS技术导则的修订草案中。

研究最后揭示了次生裂缝网络在地震活动时空迁移中的独特作用机制。通过三维可视化分析发现，次生裂缝网络的存在可使地震活动前沿相对于流体前沿推进距离缩短约40%，同时将地震能量扩散范围扩大至传统模型的2.3倍。这种时空错位现象为地震监测预警提供了新的理论依据，建议在次生裂缝密集区部署具有空间分辨率的微震监测阵列。

该研究成果为CO?地质封存工程中的风险评估与防控提供了新的理论框架和技术支撑。特别是建立的次生裂缝网络发育程度与地震风险之间的量化关系，填补了现有工程指南中的关键空白。研究提出的"压力缓冲-应力聚焦"动态演化模型，成功解释了多个已发生的封存项目中的地震活动特征，验证了其理论预测能力。相关成果已应用于四川盆地的CCS项目场址评估，为我国大型地质封存工程的设计优化提供了重要参考。