含水合物的粘土质粉砂在未排水的情况下会发生变形，这可能导致海底滑坡

《Journal of Ocean Engineering and Science》：Undrained deformation of hydrate-bearing clayey silt could generate submarine landslides

　　彭武|刘世静|舒穆利克·平克特|王云辉|黄磊|宋永晨|李阳辉
教育部海洋能源利用与节能重点实验室，大连理工大学，中国大连116024

**摘要**
随着全球能源勘探向深海环境扩展，基础设施部署越来越多地与含水合物的地层相交。海底滑坡是一种重大的地质灾害，其规模和能量远超大多数陆地斜坡失稳事件。然而，传统理论难以解释在看似稳定的深海浅层地层中广泛存在的层状条纹和弱层现象。通过低温高压原位CT三轴试验，我们证明了在非排水条件下，低渗透性含水合物沉积物（HBS）中可能会产生局部应变，从而导致渐进式的剪切带形成。在此过程中，局部结构不均匀性周围的应力集中成为微裂纹的核，这些微裂纹随后发展成宏观剪切带。这些结构改变了局部水合物的分解行为，产生了孔隙压力异常，并促进了弱层的传播。我们发现局部内部剪切先于宏观屈服发生，揭示了一种“由内而外”的破坏传播机制。这种级联的地质力学响应将微观沉积物行为与宏观斜坡不稳定性联系起来，为大规模海底滑坡的潜在触发机制提供了依据。

**1. 引言**
根据假设均匀变形的经典临界状态土力学理论，饱和沉积物在快速非排水剪切作用下应整体发生塑性变形而不改变体积，从而排除了局部剪切带形成的可能性。这一原理长期以来一直是评估海底斜坡稳定性的基础。然而，广泛的海洋勘探记录表明，在含天然气水合物的深海盆地中，从浅层沉积物失稳到大规模海底滑坡的现象普遍存在。这些自然观察结果清楚地表明，局部剪切变形和弱层传播不仅在自然条件下可能发生，而且非常普遍。理论预测与现场证据之间的这种明显矛盾揭示了我们对非排水条件下含水合物沉积物微观变形机制理解上的根本性差距。

海底滑坡作为深海地质灾害的关键表现形式，由于其触发机制和沉积物稳定性问题尚未得到充分理解，一直对地质科学界构成挑战[1]。它们的突然发生和广泛影响经常损坏海底基础设施，并对沿海社区构成灾难性风险[2][3][4]。自21世纪以来，来自深海钻探和海底观测网络的突破性数据表明，全球超过60%的活跃滑坡区与天然气水合物稳定区存在空间重叠，表明这两种现象之间存在深刻的动态联系[5]。此外，深部自由气体的积累或天然气水合物在稳定区底部的分解经常产生显著的孔隙流体超压。这种超压机制在上方高强度的水合物固结沉积物和下方含气地层之间形成了明显的力学弱区。在外部触发下，这一力学弱区通常会发展成基底滑移面，从而引发海底滑坡[6][7][8]。

天然气水合物（NGHs）在全球碳循环和能源资源开发中起着关键作用[9][10][11][12][13]。包括美国、日本和中国在内的主要经济体已经实施了国家勘探和开发战略，迄今为止已进行了多次成功的试生产活动[14][15][16]。然而，过于激进的开采策略可能会引起温度-压力条件的严重扰动，触发不受控制的水合物分解。尽管这种人为干扰受到自限性质的限制（需要持续的热量供应或压力下降——这些条件在大范围内不太可能自发发生），但它们仍可能突破这些热力学障碍，导致包括储层变形在内的级联地质力学失效[17]。此外，包括实验室实验[18,19]、数值模拟[2]和地质证据[20,21]在内的多项研究表明，这种分解可以引发海底滑坡[5,22]。

要解决这一知识瓶颈，需要深入理解微观力学行为。值得注意的是，全球近90%的天然气水合物存在于含水合物的粘土质粉砂沉积物中[23]，但其形成机制对沉积物力学行为的影响仍不完全清楚。这些细粒沉积物具有高孔隙率和低渗透性，在外部载荷作用下容易发生复杂的应变局部化[24,25]。因此，全面表征含水合物粘土质粉砂的力学行为对于揭示储层失稳机制和识别海底滑坡的潜在因素至关重要[26][27][28]。

三轴试验已被广泛用于研究纯水合物样品和含水合物沉积物的刚度、强度和应力-应变行为[29,30]。然而，现有研究主要集中在粗粒甲烷水合物沉积物上[31][32][33][34]，而关于含水合物粘土质粉砂的宏观力学行为的研究较少。值得注意的是，罗等人[35]证明了水合物颗粒增加了沉积物的胶结性和破坏强度。刘等人[36]发现，随着有效围压和水合物饱和度的增加，峰值强度和残余强度也随之提高。土壤的宏观力学行为从根本上由其微观结构和破坏机制决定。X射线CT系统因其无损性和高空间分辨率而被越来越多地用于研究含水合物沉积物（HBS）的变形机制[37][38][39][40][41]。尽管在理解含水合物沉积物的微观变形机制方面取得了显著进展，但在含水合物粘土质粉砂中剪切带的形成和弱层传播的关键知识方面仍存在空白。因此，阐明含水合物粘土质粉砂的微观变形机制并建立储层失稳特性与级联灾害机制之间的联系仍然十分必要。

为此，本研究采用了低温高压原位CT三轴试验系统[42]，以准确模拟自然含水合物沉积物中的原位非排水剪切环境。利用这一技术能力，我们的工作将系统地表征含水合物粘土质粉砂在非排水单调剪切下的微观结构演变和宏观力学响应。特别关注于阐明应变局部化的形成机制以及这些低渗透性沉积物中剪切带的传播规律。这项研究将揭示直接观察这些微观力学过程可以揭示控制储层变形和破坏的内在机制。这种理解有望为揭示水合物分解在触发弱层形成中的作用提供坚实的实验基础，最终建立一个基于微观结构的海底斜坡稳定性预测模型新框架。

**2. 材料与方法**
2.1. 样品制备
将去离子水与干燥的南海粘土质粉砂按3:40的质量比混合（粒径0.2-380μm，平均粒径D50=11.43μm，粒密度Gs=2.73 g/cm3）[43]。然后将湿润的粘土质粉砂放入橡胶套管[44]中，制成直径=20 mm、长度=40 mm的致密样品。样品的初始孔隙率和初始含水量分别为45%和25%。随后将样品放入-4°C的冰箱中冷冻超过6小时，并让其自然解冻。

2.2. 实验步骤和CT扫描程序
传统的非排水三轴试验适用于水饱和土壤，模拟快速加载且无排水的情况。相比之下，天然气水合物储层（尤其是常见的低渗透性粘土质粉砂沉积物）由沉积物、水合物、水以及通常存在的自由气体组成的多相系统构成。在生产过程中，水合物分解产生的流体由于渗透性低而难以排出，导致孔隙压力升高，有效应力降低，从而削弱沉积物稳定性。此外，孔隙压力积累和由于排水受限导致的强度软化的主要力学后果与非排水试验的理论基础一致。因此，这里使用的非排水三轴试验通过模拟气体饱和的孔隙空间，有效模拟了这些储层在早期开发过程中的力学响应。虽然气体可压缩性允许一定的内部体积调整，但严格保持了“无质量交换”的条件，这对于地质灾害评估至关重要。

实验程序采用基于微焦点X射线CT的气体水合物三轴装置（见补充图1、2、4）[45]。这种先进的方法能够在精确控制的压力和温度条件下对含水合物沉积物进行原位无损观察[17,46]。为了解决传统CT成像中甲烷水合物和孔隙流体之间的灰度模糊问题，我们在CT扫描过程中使用氙气替代客体分子，显著提高了相位区分的准确性[40,41,47]。尽管氙气和甲烷具有不同的分子间作用力，但现有文献表明，含氙气的水合物沉积物表现出与甲烷水合物高度相似的宏观力学行为、刚度和相变动力学，从而验证了其在本研究中的可靠性[47,48]。气体饱和的实验条件直接模拟了向上迁移的水合物稳定区底部的地质环境——这是商业开采的关键目标，因为那里的资源回收潜力更高[49]。实验程序如下：首先将冷冻的主样品放入压力室，确保所有连接处密封良好以防止泄漏。然后向压力室中充入氮气，并使用围压泵将围压提高到0.05 MPa，以避免损坏样品。接着将围压流体（氮气）的温度升高到18°C并保持6小时，以融化样品中的冰。通过监测样品内部温度和体积应变的稳定情况来确认冰的完全融化；读数的稳定表明固液相变已经完成。接下来，用氙气替换样品中的空气，并分别使用孔隙压力泵和围压泵将孔隙压力（氙气）和围压提高到0.60 MPa和0.65 MPa。需要注意的是，在水合物生成过程中，通过将压力泵设置为恒压模式来保持压力恒定。最后，将围压流体的温度降低到4°C以形成水合物（见补充图3）。当氙气的消耗停止时，认为水合物已经完全形成[50,51]。

含水合物粘土质粉砂样品在2.6 MPa的围压、0.6 MPa的初始孔隙压力和0.1%/min的应变率下进行非排水单阶段剪切，直到轴向应变达到21%，此时认为剪切过程完成。为了真实模拟实际海底水合物储层的含气、非饱和状态，剪切试验在严格定义的全局非排水条件下进行。与完全饱和的土壤不同，其中非排水边界规定了严格的恒定体积，而样品中的高压缩性自由气体改变了这一地质力学前提。在这里，全局非排水条件被定义为没有流体或气体交换的质量封闭系统。然而，由于气体的可压缩性，沉积物基质内部允许局部体积应变和微观多相调整。这种边界定义确保了我们的评估准确反映了深海储层固有的复杂封闭系统相互作用。在自然地质条件下，天然气水合物储层的应变率通常在10??至10?? %/s之间，远低于实验室三轴剪切率。为了模拟储层对开采引起的扰动的短期力学响应，本研究使用了0.1%/min的轴向应变率。这种加速的速率有效地捕捉了可靠的剪切强度和破坏行为，同时保持了水合物相的稳定性并减少了试验过程中的过度蠕变[52]。

剪切结束后，停止轴向加载。在压力腔内，在原始垂直有效应力下进行单阶段连续减压分解。孔隙压力以0.1 MPa/min的快速恒定速率持续降低到0 MPa。这种特定的直接单阶段减压策略旨在模拟激进的商业开采情景或极端的井底压力下降。最近的研究强调了整合地质力学测量技术以评估不同减压策略下的储层稳定性的紧迫性，并提高结果解释的完整性[53]。为了解决这一关键概念差距，我们的加速提取模型利用先进的原位CT成像技术有效地可视化和评估与快速水合物分解相关的严重地质力学风险——如瞬态孔隙压力积累、严重的结构退化和随后的斜坡不稳定。虽然在现场通常更倾向于使用较慢的逐步减压方式以减少工程事故，但直接快速的单阶段减压可以放大瞬态地质力学响应，使其成为捕捉上限局部破坏机制和完成地质力学评估框架的理想方式。等待几个小时直到体积泵指示不再变化后，认为水合物已完全分解。有关详细的实验步骤和CT扫描程序，请参阅补充信息（文本S1）。3. 结果与讨论3.1. 微观相变机制和主要发生模式如图1所示，根据形态演变，水合物生长过程可以大致分为两个不同的阶段：下载：下载高分辨率图像（865KB）下载：下载全尺寸图像图1. 含水合物的粘土质粉砂的微观相变机制。(a) 水合物形成的相变过程。(b) 含水合物的粘土质粉砂形成过程中的微观机制图。初始生长阶段：水合物团簇在沉积物中迅速成核，表现为斑块状水合物。这些斑块具有相对均匀的分布和高内部密度。在这些团簇的边缘，树枝状水合物以扩散的方式出现，从斑块状水合物的中心向外径向扩展。随后，这些结构横向扩展并相互连接，合并成更大的水合物聚集体（补充图5）。如图1a所示，在形成过程中，水合物聚集体明显表现为附着在沉积物表面的薄膜。值得注意的是，随着水合物生长的持续，出现了独特的不规则气体流动通道——这是本实验中首次观察到的现象。后期稳定阶段：在这个阶段，树枝状水合物继续向外扩展，与斑块状水合物合并形成更大的、固结的团簇。这些团簇被气体通道空间隔离，保持没有水合物分布的独立状态——这种生长行为类似于冰晶的形成[54]。这种形态演变可以通过相场方法和Mullins–Sekerka不稳定性理论来解释，尽管该理论最初是为无约束相边界制定的，但其关于过冷驱动界面不稳定性的核心原理适用于水合物成核过程中的局部孔隙环境[55]。当相变的驱动力（过冷）增加时，扩散边界层的厚度超过一个临界值。因此，Mullins–Sekerka扰动的生长率从负值变为正值，触发生长界面的不稳定性。具体来说，增强的驱动力导致平面水合物界面不稳定，形成蜂窝状结构。随着生长的进行，这些蜂窝的侧壁发生周向不稳定性，导致不对称性。然后树枝状生长主要沿轴向进行，同时发展出侧向分支，最终形成复杂的、分枝的树枝状形态。控制这一独特现象的具体因素将在后续实验中进一步研究。除了形态演变外，在这个稳定阶段还建立了特征宽度约为0.24mm的明显气体流动通道。这些没有水合物积累的通道归因于两个主要机制：1. 热抑制：水合物形成的放热性质提高了团簇边缘的温度，从而减少了这些特定区域的局部过冷，抑制了进一步的水合物成核。2. 优先流动路径：在持续的气体供应下，气体流动优先沿着阻力最小的路径（即无水合物区域）迁移，建立了稳定的供应通道，在这些通道中由于高流速和缺乏成核位点，水合物生长被有效抑制[56]。3.2. HBS在剪切破坏过程中的微观结构变化从应力-应变曲线来看，HBS在轴向应变εa的21%范围内表现出应变软化行为[57]。这是由于水合物占据的大孔隙空间减少了水合物和土壤颗粒重新排列的可用空间。因此，这种空间限制促进了轴向加载下微裂纹-剪切带的共同演化，最终导致脆性破坏。根据斜率的变化，整个应力-应变过程分为三个不同的阶段：阶段I：线性弹性（轴向应变εa：0～4.5%）：偏应力（q）与εa之间的线性关系表明具有胡克弹性，斜率定义了初始刚度。宏观上，试样保持均匀的圆柱形几何形状，没有可见的外部损伤（图1，图2）。微观上（图1，图2），虽然斑块状水合物骨架相对稳定，但树枝状水合物沿着气体迁移路径垂直坍塌，导致气体通道明显变窄。这表明局部塑性屈服和微观结构损伤在宏观弹性范围内就开始了。下载：下载高分辨率图像（2MB）下载：下载全尺寸图像图2. 含水合物的粘土质粉砂在剪切过程中的宏观应力-应变响应和纵向截面发展。(a) 在有效围压（σc）为2 MPa的条件下，含水合物的粘土质粉砂在非排水三轴剪切试验中的宏观应力-应变响应。(b) 含水合物的粘土质粉砂局部高饱和区域的灰度观测数据。(c) 局部低饱和区域的纵向截面CT图像和微观结构变化。黄色和橙色线条分别代表低饱和水合物沉积物的区域1和2。横向突出的线条代表变形阶段沉积物中的气体通道，两条平行黄色线条之间的距离是气体通道的宽度。阶段II：塑性变形（轴向应变εa：4.5～9.2%）：在这个过渡阶段，渐进的剪切放大了微观尺度颗粒的相互作用——包括滑动、旋转和重叠。由摩擦阻力引起的能量耗散导致应力-应变曲线的斜率降低。关键的是，宏观变形模式从对称压缩转变为不对称变形。在微观尺度上（图2b3–4），斑块状水合物骨架开始发生差异性压缩和膨胀，而树枝状水合物继续向内坍塌，进一步减少孔隙连通性。下载：下载高分辨率图像（2MB）下载：下载全尺寸图像图3. 剪切带的形成和发展。(a) 使用粒子图像测速（PIV）分析从含水合物沉积物的增量位移场获得的亚毫米尺度上的相应局部变形。[63](b) 增量位移变化的散点图。(c) 剪切带内外的3D局部分析。下载：下载高分辨率图像（1MB）下载：下载全尺寸图像图4. 剪切后含水合物沉积物中的解离机制和类似断层的演化。(a) 3D重建和2D垂直CT截面，突出显示了异质解离过程。(b) 概念性多尺度机制图，说明了HBS试样中微观尺度剪切带演化与海底水合物储层中经常观察到的浅层沉积断层带渐进发展的形态相似性。阶段III：屈服破坏（轴向应变εa：9.2～21%）：在这个阶段，偏应力达到大约4.8 MPa的峰值，然后随着轴向应变的增加而减小，表现出明显的屈服行为。这一峰值标志着从塑性变形到破坏的临界转变。微观结构观察（图2b5–8）显示，不同的低饱和区域合并形成广泛的弱区。同时，高饱和水合物团簇作为刚性障碍物阻碍整体变形，迫使应变集中在这些形成的弱区。这种机制最终演变为贯穿的剪切破坏表面，随着破坏的进行，剪切带的倾斜角度增加。与Zhang等人[58]观察到的峰值强度差异可能源于有效形成应力和特定树枝状微观结构的变化[34]。下载：下载高分辨率图像（856KB）下载：下载全尺寸图像图5. 水合物解离与海底滑坡演化的跨尺度耦合机制。(a) 含水合物沉积层的准稳定相。(b) 由局部水合物解离引发的次级弱层的形成（红色虚线：新出现的弱层；蓝色箭头：地层沉降；灰色气泡：甲烷气体逸出）。(c) 在外部扰动下，弱层的渐进性合并，导致新形成的脱离表面（红色虚线）沿线的宏观应力集中，并触发大规模的坡度破坏。(d) 弱层形成的微观力学机制，强调颗粒重组和局部剪切界面处的摩擦滑移减弱（红色箭头指示剪切位移方向）。局部变形的力学分析：尽管整体上处于非排水条件，但由于气体可压缩性和剪切引起的膨胀，仍会发生局部体积变化。斑块状和树枝状水合物分布的异质性产生了非均匀的孔隙压力分布，驱动了局部流体交换[59]。为了解决这些细微特征，采用了基于ImageJ的[60]灰度统计分析，其中灰度值与密度相关（水合物 > 沉积物 > 孔隙）[46]。如图2b所示，阶段I–II显示出增加的平均灰度值，反映了以压实为主的变形。然而，在阶段III，灰度值减小，表明水合物结构内的颗粒破碎和膨胀。由于高饱和区域中的较大团簇在几何上受到限制，难以旋转或重新排列[61]，因此沉积物密度在局部增加，因为孔隙被破碎的水合物细粒填充。这种局部机制通过3D重建（图2c）得到了视觉验证。在15% εa时，离散的剪切带作为不连续的裂缝出现。关注低饱和区域（R-1），在阶段I–II期间该区域显著缩小（与气体通道坍塌的观察结果一致），将连接的侧支区域减少到其原始宽度的20–47%。在阶段III，R-2区域发生明显滑移，R-1通过位错驱动的重叠与R-2合并。这种在沉积物基质内的异质、类似流体的塑性变形反映了海底滑坡剖面中观察到的“阶梯式”剪切平面形态，为解离诱导的流体迁移对滑坡进展的级联效应提供了理论支持[62]。3.3. HBS中的应变局部化和剪切带演化特征在线性弹性阶段，沉积物的增量位移基本上是轴对称的。然而，当试样进入宏观塑性变形区（εa：6–9%）时，内部位移方向出现不连续性，标志着剪切带的初始形成。值得注意的是，这一内部特征比外部宏观观察（仅在15%应变时可见）提前了6–9%的应变幅度。随后，在塑性到屈服的过渡区（εa：9–12%）内，内部位移的方向和大小都变得明显不连续，形成了一个明显的过渡区。这种不连续的变形引起了相对滑移，产生了“位错和张力”现象，有效地创造了地质弱面。这种演化揭示了一种明显的“由内向外”的破坏传播机制。从机制上讲，这种“由内向外”的破坏传播本质上是由含水合物沉积物的内在空间异质性控制的。水合物胶结的非均匀分布自然产生了沉积物基质内不同刚度和强度的局部域。在非排水变形过程中，这些结构异质性作为强烈的内部应力集中器。因此，在机械薄弱区域（例如，如图2c所示，以较低水合物饱和度为特征的区域）内，局部屈服和微裂纹的形成提前发生。根据渐进破坏力学，随着轴向应变的积累，这些孤立的内部损伤区不会同时屈服；相反，它们沿着阻力最小的路径逐步传播、相互作用并合并。这种内部结构退化不断累积，最终形成一个完全穿透试样外部边界的宏观剪切带。因此，内部不连续性与外部宏观表现之间的时间滞后不仅仅是现象学特征，而是由异质性驱动的渐进破坏的直接地质力学后果。这种“由内向外”的机制对于早期识别海底滑坡的前兆信号具有关键意义。随着剪切的进行，剪切带的平坦度增加，而其厚度略有减小，最终导致剪切带的逐步破坏，最终剪切带的倾斜角度为57°，厚度为3.48毫米（大约是D50的34倍）。相比之下，Yoneda等人[41]的研究中，含水合物饱和度为62%的沉积物中剪切带的厚度仅为所用沙粒直径的2.3倍；Lu等人[64]的研究中，剪切带的出现滞后于水合物饱和度的增加，并且随着有效围压的增加而提前出现，剪切带的厚度随着有效围压的增加而减小。与之前的研究相比，本实验中的剪切带倾斜角度θ符合其规律。剪切带倾斜角度θ由Mohr-Coulomb破坏准则表示为θ = π / 4 + φ / 2，其中φ是材料的内摩擦角，φ为24°。从图3b可以看出，运动学分析显示，从3%到6%的轴向应变，增量位移向量分布在轴的两端。随着轴向应力的增加，这些点集中并均匀地对齐到水平轴的零点，而垂直位移可以忽略不计。（注意：0–3%轴向应变时的正垂直增量位移是由于加载头轻微初始倾斜造成的实验误差）。本质上，水平增量位移的比例随着剪切的增加而增加，宏观上表现为剪切带倾斜度的增加。这种运动学特征与由摩擦减弱引发的海底滑坡基底剪切带中的不稳定滑动具有类似的动态机制。通过比较剪切带内外的3D重建水合物分布（图3c），可以获得更深入的见解。空间相关性分析揭示了明显的变形异质性：在剪切带外部，沉积物经历均匀压实，同时气体迁移通道变窄；相反，在剪切带内部，变形由滑移控制，导致气体通道重新定向并表现出明显的不连续性。与含水合物的砂质沉积物中观察到的明显膨胀不同（例如，Higo等人[65]），这些粘土-粉砂沉积物中的剪切带没有明显的宏观体积膨胀。尽管从“位错和拉伸”机制理论上预期会有拉伸变形，但未观察到明显的膨胀归因于分辨率限制。然而，关键的是，该区域表现出整体变暗（灰度值降低），表明由滑移变形产生的拉伸微孔存在，但非常微小（<20μm）。剪切带内外观察到的变形异质性表明，在宏观上未排水（质量恒定）的沉积物中，体积调整是通过复杂的多物理耦合过程实现的。除了颗粒重排和过量孔隙气体的产生外，气体孔隙流体的可压缩性和潜在的应力诱导水合物熔化[66]是导致整体变形响应的两个关键因素。由于样品处于全局未排水条件下，这种高气体可压缩性理论上允许由局部孔隙压力梯度驱动的残余孔隙水的微观尺度重新分布。从机制上讲，局部剪切带内的膨胀行为将水分集中在该区域，这严重降低了局部有效正应力，促进了渐进的应变软化，并加速了沉积物的宏观强度退化。然而，由于固体水合物对孔隙的严重阻塞以及粘土粉砂基质的固有低渗透性，大规模流体迁移受到显著阻碍。我们承认，虽然原位CT成像有效地捕捉到了固体骨架的结构退化，但直接追踪这些微观流体重新分布仍然是一个限制。具体来说，剪切带内骨架的膨胀倾向压缩了其中的气体，为局部应变适应提供了必要的内部空间。同时，在剪切局部化产生的高应力集中区域，水合物可能发生“应力诱导熔化”——这是一种由于塑性变形导致局部平衡温度升高而引发的解离。这种相变不仅消耗能量并软化材料，还释放了填充孔隙空间的流体，进一步在动态反馈循环中改变局部体积和孔隙压力。此外，必须认识到这些微观力学和热力学反馈对施加的变形速率非常敏感。如方法论中详细说明的，相对较快的实验室应变速率促进了更窄、更集中的剪切带的发展，并由于强烈的未缓解的颗粒摩擦而放大了局部剪切加热。因此，这种快速的能量积累在剪切带内引发了更严重的局部应变软化和滑移破坏。相比之下，自然地质环境中极慢的应变速率将允许更广泛的流体扩散和热耗散。从概念上讲，这种长期的地质变形将导致更宽的剪切局部化区域、更分散的孔隙压力演变和较少的突然渗透性变化。因此，虽然我们的原位观察有效地捕捉到了沉积物中的微观尺度局部破坏机制，但在推断长期储层不稳定性时，必须仔细考虑这些依赖于速率的热力学和力学差异。因此，本研究中观察到的“未排水剪切下的体积变形”是气体可压缩性、水合物相变和骨架塑性变形的耦合结果，代表了比传统饱和流体土壤更复杂的行为。

3.4 相变和沉积物变形：对海底滑坡触发机制的影响
水合物解离过程可以分为三个阶段：
（1）预解离阶段：在剪切后的立即状态下，树枝状和斑块状水合物的宏观形态基本保持完整。气体迁移继续沿着初始沉积结构建立的固有渗透率各向异性进行，保持了明确的流动路径。
（2）再生阶段：解离主要沿着垂直梯度（从底部到顶部）开始。然而，解离的吸热性质导致局部冷却，在气体通道附近引发热抑制解离[67]，从而减缓了进一步的破坏。同时，在沉积物顶部发生次级水合物形成，迁移的流体在孔隙空间内重新形成随机填充物——这种形态与主要生长纹理不同（补充图6）。
（3）目标解离阶段：主要解离前沿与剪切带垂直对齐，而次级前沿沿着局部低饱和区内的气体路径向内传播。关键的是，沉积物界面处的剪切诱导结构不连续性由于颗粒破碎和孔隙塌陷而演变为低渗透性滑移表面。如图5b所示，这些障碍物阻碍了孔隙压力的消散，形成了局部超压区。
值得注意的是，HBS中层次化剪切带的演变与海底水合物储层中经常观察到的浅层沉积断层带的渐进发展在运动学和形态学上具有明显的相似性（图4b）。在微观尺度上，渐进的水合物解离和颗粒破碎产生了低渗透性滑移表面。这一过程由解离引起的流体释放驱动，产生了局部超压，降低了有效应力，并促进了摩擦减弱。在宏观尺度上，这种结构退化与未固结海洋沉积物中连续弱界面的形成相类似，通过相邻层的水平偏移产生了明显的断层状分离区[68]。这种层次化的破坏序列表明了一种潜在的多尺度类比，控制着不稳定性：最初局部化的解离引发了微裂纹，这些微裂纹合并成传播的剪切带，从而在外部触发因素下为大规模破坏提供了结构上的预条件[4]。然而，必须认识到，将毫米级三轴未固结样品的局部变形机制外推到千米级的宏观断层带存在显著的不确定性。实验室测试中的固有边界条件和大规模效应限制了直接的定量放大。因此，这种多尺度机制图主要作为一个概念性假设。将这些三轴测试结果外推到自然沉积断层和滑坡触发机制的最终有效性需要通过大规模物理建模和原位地球物理监测来进一步确认。
基于前面的微观结构分析和解离验证的整合，图5描绘了驱动水合物相关海底滑坡的跨尺度地质力学耦合机制。这种宏观不稳定性的起源本质上根植于水合物的空间异质性。如在未排水剪切阶段所阐明的，应变局部化和剪切带优先在这些机械受损、低饱和度区域内开始并传播。关键的是，这些局部损伤区域为宏观尺度破坏在结构上进行了预条件。就像在外部扰动（例如，海洋变暖或地震活动）下一样，沉积物基质内部发生了热力学不稳定。从机制上讲，解离降低了沉积物的凝聚力，而气体和水的释放引发了过量孔隙压力的快速积累。这种超压显著降低了有效应力，使系统从稳定状态转变为临界状态。当积累的孔隙压力超过临界阈值时，流体驱动的膨胀在预先存在的剪切带区域内引发了微裂纹。这些微裂纹逐渐重组为相互连接的次级弱层（图5b），作为宏观破坏的前兆。
在持续的外部强迫下（例如，地震加载或快速沉积），这些局部弱层合并成穿透性的分离表面（图5c）。这一级联过程通过三种耦合机制协调了结构完整性的突然丧失：（1）水合物骨架的胶结破坏，（2）瞬态超压引起的局部流体化，以及（3）周围地层的渐进机械退化。这些效应的协同作用沿着发展的分离界面驱动了应力局部化，最终引发了大规模的海底滑坡。关键的是，由此产生的弱层（图5d）表现出与浅层沉积断层带类似的摩擦特性：明显的颗粒重组产生了结构上重新工作的颗粒界面，而沿着这些高度连接的路径的膨胀解离流体的集中迁移引发了显著的压强依赖性滑移减弱。我们假设这种严重的摩擦阻力降低本质上控制了宏观弱层的形成。鉴于广泛的地质文献已经确定这种结构受损层是大规模海底滑坡的主要触发因素[69,70]，我们的概念模型填补了一个关键的机制空白。因此，微观剪切带的渐进连通性不仅仅导致局部屈服；相反，它从根本上预准备了这些表面破坏的起源[71,72]。这一框架因此提供了一个从微观到宏观的地质力学联系假设，有助于理解海底斜坡不稳定性的规模和运动学。

4. 结论
总之，利用低温高压原位CT三轴成像，这项研究从根本上挑战了传统观点，即含水合物沉积物中的未排水剪切排除了异质变形的可能性。我们发现，宏观稳定性被一种关键的形态依赖机制所掩盖：树枝状水合物的过早崩塌，不同于坚固的斑块状骨架，是土壤强化和变形异质性的主要驱动因素。关键的是，我们的观察表明，这些微米级裂纹在宏观线性弹性范围内就引发了应变集中，证实了显著的结构退化在宏观屈服点之前就开始了。
关键的是，这项研究量化了一种独特的“由内而外”的破坏传播机制，其中内部剪切带的形成先于宏观表面表现，其时间滞后为3-6%。这种时间滞后表明，当前的监测基线可能从根本上低估了海底储层不稳定性的开始，因为关键损伤在沉积物核心中“地震学上”积累，然后才可能被外部检测到。
此外，我们建立了一个新的跨尺度耦合模型，其中剪切带的演变从根本上重新配置了水合物解离的热力学，触发了一系列地质力学反应：局部解离导致异常的孔隙压力积累；降低的有效正应力促进了沿剪切带的摩擦减弱；这些局部弱区合并将原始剪切带转变为次级的、穿透性的分离表面。我们假设这种“微裂纹到宏观滑坡”的连锁反应控制了大规模海底滑坡的核化。因此，检测这些初始“由内而外”的微裂纹的地球物理特征为深海地质灾害的早期预警和缓解提供了有希望的新途径。此外，必须明确指出，尽管这种三轴方法可以成功捕捉封闭系统中的微观力学耦合、局部孔隙压力响应和水合物分解在未排水条件下的变化，但它也有其固有的局限性。具体来说，它无法完全捕捉水合物颗粒和粉砂颗粒的更详细胶结过程，这是我们未来将研究的内容。

数据可用性
作者声明，支持本研究发现的数据可在主文章、补充信息和源数据文件中找到。源数据随本文提供。本研究中的CT图像堆栈可在Zenodo上获取（https://doi.org/10.5281/zenodo.14693486）。图像分割和灰度分析通过Fiji/ImageJ [73]进行（可在Fiji获取），体积渲染、横截面变化和交互式阈值处理由Avizo Software（https://www.thermofisher.cn/cn/zh/home/global/forms/industrial/avizo-software-trial.html）完成。