《Physics and Chemistry of the Earth, Parts A/B/C》：Thermo-mechanical response of Sakesar limestone under cyclic steam injection: Insights from the Potwar Basin, Pakistan
编辑推荐：
本研究针对巴基斯坦波特瓦尔盆地Balkassar油田的始新统Sakesar组裂缝性碳酸盐岩储层，开展了循环蒸汽刺激（CSS）条件下的热-水-力（THM）耦合响应研究。通过集成实验力学测试、微CT成像和数值模拟，揭示了热循环导致岩石力学性能非线性退化（单轴抗压强度降低33%），但渗透率因裂缝网络连通性增强而翻倍（28升至56 mD）。THM模拟识别出井筒附近10米内为高滑移趋势区，表明热采在提升注入能力的同时也增加了裂缝剪切再活化的地质力学风险。该研究为波特瓦尔盆地及类似构造活动背景下碳酸盐岩储层的热采优化与安全设计提供了关键科学依据。
在全球能源需求持续增长和常规油气产量递减的背景下，如何有效动用成熟裂缝性碳酸盐岩储层中的剩余油成为行业焦点。热采强化（EOR）技术，特别是循环蒸汽刺激（CSS），通过降低原油粘度和改善储层渗流能力，被视为一种有效的开发方式。然而，注蒸汽过程中反复的加热-冷却循环会在储层岩石中引发复杂的热-力耦合过程，可能导致岩石强度劣化、裂缝系统演变，并诱发井筒稳定性风险。位于巴基斯坦波特瓦尔盆地的Balkassar油田就是一个典型的成熟裂缝性碳酸盐岩油田，其主力产层——始新统Sakesar组灰岩储层正面临产量递减的挑战。尽管该区域地质研究已有一定基础，但针对CSS条件下储层岩石热-力学行为，尤其是应力-应变响应、裂缝演化和渗透率变化规律的认识仍显不足，这制约了热采技术在该类储层中的安全高效应用。
为了填补这一知识空白，研究人员以Balkassar油田Sakesar组灰岩为研究对象，开展了一项集地质表征、岩石力学实验、微结构成像和数值模拟于一体的综合性研究。该研究旨在揭示循环热载荷下储层岩石的热-力学响应机理，量化热损伤对岩石力学性质和渗流能力的影响，并评估热采过程中的地质力学风险，从而为波特瓦尔盆地及类似碳酸盐岩储层的热采方案优化提供科学支撑。相关研究成果发表在《Physics and Chemistry of the Earth, Parts A/B/C》上。
为开展此项研究，作者主要应用了以下几项关键技术方法：首先，对取自Balkassar油田三口生产井（Balkassar-OXY-01, Well-7, Well-8）Sakesar组储层段（深度2467–2600米）的岩心样品进行了基线表征，包括岩石力学参数（单轴抗压强度UCS、巴西拉伸强度、杨氏模量等）和物性参数（孔隙度、渗透率）测试。核心实验是模拟现场CSS条件的室内热-力学循环测试，将岩样在25°C至180°C之间进行10次加热-冷却循环，并在每个循环后测量其力学性质和渗透率的变化。利用微计算机断层扫描（微CT）技术对热循环前后的岩样进行高分辨率（体素分辨率20微米）扫描，定量分析裂缝孔径、连通性的演变规律。最后，基于实验获取的参数，建立了耦合的热-水-力（THM）数值模型，模拟蒸汽注入过程中井筒附近的温度场、应力场重分布及裂缝再活化潜力（滑移趋势）。
4.1. 基线储层物性
研究首先明确了Sakesar灰岩的初始条件。平均孔隙度为21.1%，基线渗透率为28毫达西（mD），表明储层兼具基质孔隙和裂缝提供的渗流通道。单轴抗压强度（UCS）为47.5兆帕（MPa），巴西拉伸强度为4.9 MPa，静态杨氏模量为22.4吉帕（GPa），泊松比为0.24，这些参数表明该灰岩属于中等强度、具有一定脆性的裂缝性碳酸盐岩。应力-应变曲线呈现出典型的弹塑性响应，为评估热循环后的性能退化提供了基准。
4.2. 热循环下的力学性能退化
热-力学测试结果表明，随着热循环次数的增加，岩石的强度和刚度出现渐进式退化。经过10次循环后，单轴抗压强度（UCS）从47.5 MPa降至31.5 MPa（降幅33%），拉伸强度从4.9 MPa降至3.1 MPa，静态杨氏模量从22.4 GPa降至16.7 GPa。这种退化归因于热胀冷缩在矿物颗粒边界和微裂缝尖端产生的循环热应力，导致微裂纹萌生、扩展和汇聚。力学性能的退化在前期循环（如前5次）中较为显著，随后趋于缓和。在恒定法向刚度（CNS）条件下进行的节理剪切试验也显示，热循环后裂缝面的峰值剪切强度降低，而剪切膨胀量增加，表明热作用降低了裂缝的抗剪能力。
4.3. 裂缝孔径-渗透率关系
微CT成像分析揭示，渗透率的增强主要受裂缝网络连通性演化的控制，而不仅仅是单一裂缝孔径扩大的结果。经过10次热循环后，平均裂缝孔径从0.15毫米增至0.39毫米。更重要的是，图像三维重建显示，热循环促使原本孤立或连通性差的微裂缝相互连接，形成了更为连续的有效渗流网络。渗透率测量结果与此吻合，从基线28 mD不可逆地增加至56 mD。渗透率的增长在循环初期近似线性，但在后期（约第5次循环后）呈现非线性加速趋势，这与微CT观察到的裂缝网络发生“渗流”转变（即形成贯通流道）密切相关。加热和冷却过程中的渗透率滞后现象表明，裂缝在冷却后无法完全闭合，造成了渗透率的永久性提升。
4.4. 热应力重分布与滑移趋势
THM数值模拟结果清晰展示了蒸汽注入期间井筒周围强烈的热应力重分布。加热导致岩石热膨胀，在近井地带（约10米半径内）产生显著的热应力，改变了原地应力场。滑移趋势（Ts）分析表明，在加热阶段，井筒附近预存裂缝面上的剪应力（τ）增加，而有效正应力（σn）降低，导致Ts值显著升高，在部分区域接近甚至达到1.0（设定摩擦系数μ=0.6），这意味着这些裂缝处于剪切再活化的临界状态。高滑移趋势区与高温梯度区高度吻合，表明热诱导的剪切滑移是近井地带渗透率增强的重要机制之一，但也同时带来了裂缝失控扩展和井壁失稳的风险。
4.5. 多周期循环下储层物性演化
综合实验和模拟结果表明，力学性能退化和渗透率增强随热循环次数的演变并非线性。力学参数（如UCS、杨氏模量）的退化速率在经历前期快速下降后逐渐减缓，表现出某种“稳定化”趋势。而渗透率的提升则持续进行，甚至在力学退化减缓后仍能通过裂缝网络的进一步优化而继续改善。这种解耦现象说明，在一定的操作窗口内，通过热循环实现渗透率显著改善而将力学完整性维持在可接受水平是可行的。但需要指出，本研究限于10个循环，现场应用可能涉及更多周期，长期效应需进一步评估。
研究的讨论部分深入剖析了上述发现的工程意义。Sakesar灰岩在热循环下表现出的力学性能退化与渗透率增强的“权衡”关系，是热采设计必须面对的核心问题。THM模型识别出的近井高滑移趋势区，在Balkassar油田这种受喜马拉雅造山运动影响的挤压性构造背景下尤其需要警惕，因为地应力场本身就有利于剪切破裂。这意味着，热采在提升原油流动能力的同时，确实可能激活天然裂缝系统，对井筒完整性和断层密封性构成潜在威胁。因此，成功的现场应用依赖于精细的注入参数控制（如蒸汽温度、压力、注入速率、循环周期）以及实时监测（如微地震监测、井下温度和压力监测）。
综上所述，这项研究通过多学科交叉的方法，首次系统揭示了波特瓦尔盆地Sakesar组裂缝性碳酸盐岩在循环蒸汽刺激下的热-力学行为规律。核心结论是：循环热载荷会不可逆地改善储层渗透率，其主要机制是裂缝网络连通性的增强，而不仅仅是孔径扩大；但同时会导致岩石力学强度逐步退化，并在井筒附近产生高滑移趋势区，增加地质力学风险。该研究建立了一个将实验室实验、微观结构观测和数值模拟相结合的集成分析框架，为Balkassar油田及地质条件相似的裂缝性碳酸盐岩储层安全、高效地实施热采强化（EOR）策略提供了关键的理论依据和设计指南。研究成果强调，在追求注入能力提升的同时，必须对热采诱发的地质力学变化进行定量评估和有效管控，以实现提高采收率和保障长期生产安全之间的最佳平衡。
