《Petroleum Science》：Ultrasonic monitoring of shear rupture processes in hydrate-bearing sediments via discrete element method-based simulation
编辑推荐：
水合物沉积物在剪切变形过程中的声学响应对于解释海洋水合物储层的不稳定过程至关重要，然而在力学非均质系统中，潜在机制仍受到不足的约束。现有研究主要关注静态刚度-速度关系或具有均匀水合物分布的沉积物，尚不清楚层状水合物非均质性如何控制剪切诱导的声学行为和微观结构演
水合物沉积物在剪切变形过程中的声学响应对于解释海洋水合物储层的不稳定过程至关重要，然而在力学非均质系统中，潜在机制仍受到不足的约束。现有研究主要关注静态刚度-速度关系或具有均匀水合物分布的沉积物，尚不清楚层状水合物非均质性如何控制剪切诱导的声学行为和微观结构演化。在此，研究人员开发了一种基于离散元方法（DEM）的超声脉冲测速方法，以研究具有层状水合物分布的沉积物在剪切变形下的耦合力学、声学和微观结构响应。结果识别出两种由子层水合物分布控制的失效路径。当一个子层无水合物时，分布式微裂纹逐渐发展，导致信号振幅持续衰减和波形逐渐退化。相比之下，当所有子层都存在水合物时，声学响应经历三个阶段：微损伤引起的振幅减小和波形展宽、以相位滞后为特征的剪切局部化起始、以及界面主导的宏观剪切破坏。相应地，纵波速度（P-wave velocity）呈现出三阶段演化，包括初始压实期间的短暂增加、与微裂纹成核相关的持续衰减、以及临界波动阶段的异常恢复，反映了剪切带内的局部愈合。时频分析进一步表明，高频能量聚集与微剪切带成核相关，而低频不连续性源于层间界面效应，提供了沉积物非均质性和损伤演化的频谱指标。通过将声学、力学和微观结构响应整合到一个一致的框架中，本研究阐明了层状水合物非均质性如何控制剪切诱导的声学特征，并为监测水合物层的不稳定性提供了新的见解。
# 论文解读：基于离散元方法模拟的水合物沉积物剪切破裂过程超声监测
## 研究背景与问题
海洋天然气水合物储层在开采或地质扰动下可能发生剪切失稳，其声学响应是监测储层稳定性的关键手段。然而，天然水合物沉积物往往呈现层状非均质分布，现有研究多集中于均匀水合物分布或静态刚度-速度关系，未能揭示层状非均质性如何控制剪切诱导的声学行为与微观结构演化。这一认知空白限制了超声监测技术在复杂储层中的应用。为此，研究人员基于离散元方法（DEM）开发了超声脉冲测速模拟框架，旨在阐明层状水合物沉积物在剪切变形下的力学-声学-微观结构耦合机制，为预测失稳提供理论依据。该论文发表在《Petroleum Science》。
## 主要关键技术方法
研究人员采用离散元方法（DEM）建立圆柱形沉积物试样（半径1.5 mm，高度6 mm），通过“砂+水合物”方法随机生成砂粒（粒径0.1–0.4 mm）和水合物颗粒（半径0.03 mm）。砂粒间采用线性滚动阻力接触模型，砂粒与水合物、水合物间采用线性平行粘结模型。通过各向同性压缩至有效围压1 MPa固结后，制备四组层状水合物分布试样（上下子层饱和度不同）。同步进行三轴剪切（轴向应变速率1×10
?9 m/s至应变15%）与超声脉冲测速（频率0.5 MHz，每1%应变记录一次），监测宏观应力、应变、孔隙率及微观颗粒滑动、裂纹扩展。方法验证通过对比模拟纵波速度（V
p）与实验数据完成。
## 研究结果
### 3.1 波形特征（Waveform characteristics）
通过分析不同轴向应变（ε=0%–15%）下的输出波形，发现层状水合物分布控制两种失效路径。当存在无水合物子层时，波形振幅持续衰减（ε=12%时超过30%），无显著相移，反映分布式微裂纹扩展。当所有子层均含水合物时，声学响应经历三阶段：ε<6%时振幅减小和波形展宽（微损伤）；ε=9%时出现0.5–1.0 μs相位滞后（剪切局部化起始）；ε=12%–15%时相位延迟>2.5 μs且振幅衰减>50%（界面主导宏观剪切破坏）。
### 3.2 压缩：应力、裂纹、纵波速度（Compression: stress, crack, P-wave velocity）
纵波速度（V
p）呈现三阶段演化：初始压实阶段（ε<2%）V
p短暂增加最多7%（骨架密实化）；加速衰减阶段（2%≤ε≤9%）V
p持续下降0.5%–1.1%/应变（微裂纹成核）；临界波动阶段（ε>9%）V
p异常恢复（如从1200增至1400 m/s），对应剪切带内局部锁定和颗粒再胶结，形成损伤自抑制效应。
### 3.3 沉积物剪切破裂特征（Sediment shear rupture characteristics）
通过孔隙率、体积应变率和裂纹分布图分析，高饱和度子层呈现收缩，低饱和度子层呈现膨胀。剪切裂纹主要集中于高饱和度区域（下部），拉伸裂纹集中于低饱和度区域（上部）。层状沉积物体积应变率波动更大，更难预测失稳行为。
### 3.4 波谱特征（Wave spectrum characterization）
时频分析显示，高频能量聚集与微剪切带成核对应，其时空分布与体积应变率异常峰吻合。低频不连续性源于层间界面效应，可作为识别层状非均质性的特征指标。不同水合物分布下，高能带结构（横向 vs. 纵向）反映局部剪切带或弥散裂缝网的发育。
## 讨论与结论
讨论部分指出，层状水合物分布通过控制裂纹成核与剪切局部化路径，主导声学响应特征。异常V
p恢复与应力强化、裂纹增长率降低共同构成损伤自抑制现象，为工程中压力调控提升储层安全提供了可能。结论总结：本研究建立了基于DEM的超声脉冲测速框架，揭示层状水合物沉积物剪切破坏的两种模式（渐进式分布式损伤 vs. 三阶段局部化破坏），明确了纵波速度三阶段演化与时频特征对损伤演化及非均质性的指示作用。该成果为水合物储层稳定性监测提供了新见解。
