基于两相麻坑模型估算含水合物沉积物下伏游离气饱和度的新方法及其在挪威Storegga滑坡区的应用

《Geosciences》：Two-Phase Pockmark Modeling and Gas Saturation Estimation Beneath Hydrate-Bearing Sediments: Insights from the Storegga Slide Zheng Su, Yifan Wu, Chao Yang and Nengyou Wu

本文针对含水合物沉积物下方游离气带中气体饱和度的定量估算难题，提出一种创新的两相麻坑模型。该研究以挪威Storegga滑坡区为应用案例，揭示了在麻坑形成前兆阶段，下伏游离气储层的有效气体饱和度阈值，为深入理解天然气水合物-麻坑系统的渗流动力学与地质灾变过程提供了关键定量约束。

　　

在广阔的大陆边缘海底，分布着无数被称为“麻坑”的碗状凹陷。它们并非地表的天然装饰，而是深部流体，尤其是甲烷气体，冲破层层沉积物束缚，向上渗漏并喷出海底后留下的“疤痕”。这些流体渗漏现象与海底天然气水合物的赋存密切相关。天然气水合物是一种在低温高压下由甲烷和水形成的冰状固体，其下方的游离气（Free Gas）是驱动渗流、形成麻坑的关键“动力源”。然而，含水合物沉积物下方的游离气带通常包含水和气体两相。一个关键的科学难题在于，如何仅仅通过观测到的海底麻坑形态及其下方气体烟囱的地震反射特征，来定量反演其下方游离气储层中的气体饱和度。这不仅是理解海底甲烷渗漏过程、评估其气候效应和地质灾害风险的核心，也对天然气水合物资源的勘探开发至关重要。遗憾的是，基于传统单一相模型的估算方法在此场景下存在局限。为此，研究人员在《Geosciences》期刊发表了一项研究，致力于解决这一难题。

为了开展研究，作者主要采用了理论建模与应用分析相结合的技术方法。首先，在Cathles等人（2010）提出的单相麻坑模型框架基础上，发展了适用于水合物环境下气-水两相系统的数值模型。该模型的核心是建立了麻坑深度、游离气带厚度与气体饱和度之间的解析关系。其次，将所建模型应用于挪威外海的Storegga滑坡区。这是一个地质特征清晰的研究区，其天然气水合物、麻坑及气体烟囱状地震异常均有详实的地震反射数据记录，为模型验证与应用提供了关键参数约束，如底模拟反射面（BSR）确定的封闭层深度、游离气带厚度和麻坑几何形态等。

3.1. 控制方程与解析公式

研究者推导了控制气体烟囱向上传播的核心方程。该方程描述了在浮力驱动下，气-水两相混合物在沉积物孔隙中向上运移的动力学过程。方程的关键在于引入了两相混合流体密度（ρ_f），并将其与孔隙水密度（ρ_w）的差值作为驱动压差（ΔP）的主要来源。通过求解该微分方程，可以将气体烟囱的生长高度（h_f）与时间（t）关联起来，从而估算烟囱从突破水合物封闭层到抵达海底所需的传播时间。

4. 应用：Storegga滑坡

模型在Storegga滑坡区的应用结果表明，在麻坑形成的初始（近阈值、喷发前）阶段，水合物层下方厚度为50-100米的游离气储层中，有效气体饱和度估计约为1.36%–1.58%。这一数值代表了麻坑形成和封闭层破裂所需的最低阈值条件，而非现今仍在持续渗流的状态。同时，模型预测的气体烟囱在初始阶段的传播时间尺度约为200年。这些估算为理解该区域古渗流事件的发生条件提供了定量依据。

5. 讨论

研究结论强调，本文开发的两相麻坑模型实现了对含水合物沉积物下伏游离气饱和度的定量估算，突破了传统单相模型的局限。模型应用揭示，在Storegga地区，只需相对较低的气体饱和度（约1.4%-1.6%）即可在特定地质时间内触发气体烟囱的向上传播和麻坑的形成。这一发现具有重要意义：首先，它提供了利用常见的地震和海底地形观测数据（麻坑深度、BSR深度、游离气带厚度）来间接估算深部气体饱和度的新工具，这对缺乏直接钻探数据的区域尤为重要。其次，该研究深化了对天然气水合物系统渗流动力学过程的理解，表明麻坑的形成是气体在毛细管封闭下逐渐累积、最终在达到阈值压力后发生突发性渗流的结果。最后，该模型和估算结果有助于评估海底甲烷渗漏的潜在通量及其对全球碳循环的影响，同时对海底工程设施的地质灾害风险评估具有参考价值。需要注意的是，该模型旨在约束渗流发生前的“近阈值”状态，而非描述持续渗流的动态过程。渗流活动可能因烟囱内次生水合物的形成或麻坑周围自生碳酸盐岩的胶结而停止，留下残留的气体烟囱构造。总之，这项研究为探索天然气水合物系统中的流体迁移机制及其环境与资源意义开辟了新的定量途径。