一种快速确定沉积物中气体水合物相平衡条件的优化方法

《Journal of Ocean Engineering and Science》：An Optimized Method for Rapidly Determining the Phase Equilibrium Conditions of Gas Hydrate within Sediments

　　

确定相平衡条件是解决科研和工程领域中与水合物相关问题的第一步。传统的步进加热等容法（SHI method）虽然常用且被接受，但耗时较长。在本研究中，研究人员对传统的步进加热等容法（SHI method）进行了优化，以快速确定沉积物中气体水合物的相平衡数据。与传统方法不同，优化方法提出了一种指数函数来表征每个温度步长下的水合物分解曲线（气体压力随时间变化）。校准后的关键参数可用于计算相平衡温度和压力。与传统方法相比，水合物分解曲线无需达到稳态。此外，研究人员采用优化方法测定了CH₄水合物在几种不同沉积物中的相平衡数据，并以传统方法作为参考。优化方法与传统方法在相平衡数据上高度一致，证实了优化方法能够有效确定相平衡数据。更重要的是，优化方法大幅节省了实验步骤的时间。

甲烷水合物（CH₄ hydrate）广泛存在于海洋大陆架和永久冻土区，是一种绿色能源，其开采常采用降压、加热或注入抑制剂等方法使水合物分解为甲烷气体和水。然而，水合物分解可能改变原位应力状态并削弱沉积物内聚力，引发海底滑坡和大规模沉降；同时，甲烷是强温室气体，若开采失控释放至大气，将加剧全球变暖。因此，系统理解气体水合物的形成与分解行为是未来商业开采的基础。此外，基于水合物的CO₂海底封存也被视为缓解全球变暖的潜在策略。相平衡条件是分析水合物形成/分解行为的第一步，对CH₄水合物储量评估、开采效率及底层稳定性评价至关重要。目前，最常用的相平衡测试方法是等容测量法，包括步进加热等容法（SHI method）和连续加热等容法（CHI method）。Tohidi等（2006）指出，SHI method比CHI method具有更高的精度和可重复性，因此被广泛采用。然而，SHI method耗时显著：每个温度步长（通常0.5 K）需约24小时甚至更长以达到平衡态，测定一条完整相平衡曲线约需一周时间，既不经济也不利于快速解决工程问题。因此，有必要开发一种有效降低测试时间的方法。

研究人员在本文中提出了一种快速测定沉积物中水合物相平衡条件的优化方法。该方法基于对传统SHI method的改进，通过推导指数函数表征水合物分解曲线，利用非稳态实验数据拟合关键参数，从而计算相平衡温度和压力，无需等待系统达到稳态。通过6个不同案例（包括细粒沉积物、砂质沉积物、不同水含量、含盐溶液及CO₂水合物）验证了该方法的合理性和有效性，并与传统方法进行对比。结果表明，优化方法所得相平衡数据与传统方法高度一致（R²>0.99），且显著节省实验时间（节省率27.5%~40%）。该研究发表在《Journal of Ocean Engineering and Science》，为快速获取水合物相平衡条件提供了实用工具，对水合物资源开发及工程应用具有重要价值。

研究人员采用的主要关键方法包括：基于步进加热等容法（SHI method）的优化策略，核心是推导并应用指数函数P(t)=P₀+ΔP[1-exp(-(t-t_i)/τ)]拟合每个温度步长下的压力-时间曲线，通过拟合得到参数ΔP和τ（或C_j和τ_j），进而利用极限计算相平衡压力P_eq=P₀+ΣΔT_j/C_j和温度T_eq=T₀+ΣΔT_j。实验装置采用结晶器（直径99.8 mm，高148.7 mm，最大压力20 MPa），通过水浴控温（-20℃~+100℃，精度±0.1℃）。样本来源：四种沉积物，其中沉积物(a)类似南海神狐海域细粒沉积物，沉积物(b)和(d)类似日本南海海槽砂质沉积物，沉积物(c)为添加膨胀土的砂质沉积物。水含量为5%~12.7%，盐溶液为0.15 mol/L和0.5 mol/L NaCl。

**Case 1: 细粒沉积物低水含量（ω=5%）**
通过优化方法对CH₄水合物在细粒沉积物（ω=5%）中设置8个温度步长，采用指数函数拟合压力-时间曲线，得到参数C_j和τ_j，拟合曲线与实测数据高度一致（R²=0.99）。与传统方法对比，优化方法所得相平衡数据几乎重合（R²=0.990），且测试时间从71.1小时缩短至46小时，节省35.3%。结果表明，优化方法能有效测定细粒沉积物中CH₄水合物的相平衡数据。

**Case 2: 细粒沉积物较高水含量（ω=10%）**
采用优化方法测定细粒沉积物（ω=10%）中CH₄水合物的相平衡数据，结果显示与传统方法完全一致（R²=0.996）。在低温段（<276 K）数据精确重合，高温段略有偏差，归因于拟合精度问题。优化方法节省时间27.5%。进一步比较Case 1和Case 2发现，较低水含量（ω=5%）的相平衡压力更高，这表明毛细管和吸附作用降低了孔隙水活性，使得水合物需要更苛刻的条件稳定。

**Case 3: 砂质沉积物（ω=10%）**
对砂质沉积物中CH₄水合物进行测试，优化方法与传统方法结果高度一致（R²=0.996），且节省时间40%。砂质沉积物中相平衡数据偏离纯水合物曲线，归因于小孔隙中的毛细效应。与细粒沉积物（ω=10%）相比，砂质沉积物的相平衡压力更低，说明细粒沉积物表面负电荷对水的吸附降低了水活性。

**Case 4 和 Case 5: 含盐溶液沉积物**
采用0.15 mol/L和0.5 mol/L NaCl溶液，优化方法测定CH₄水合物相平衡数据，与传统方法高度一致（R²分别为0.999和0.998）。盐度降低水活性，导致相平衡曲线向高压低温方向偏移，高盐浓度下偏移更显著。结果证实优化方法适用于含盐沉积物。

**Case 6: CO₂ 水合物**
对砂质沉积物中CO₂水合物测试，优化方法与传统方法结果一致（R²=0.998），表明该方法适用于CO₂水合物。

**时间消耗对比**
优化方法在6个案例中均显著节省时间，节省率27.5%~40%。当测量更多相平衡点时，时间节省更显著，对实际工程具有重要意义。

讨论部分指出，优化方法的核心是通过拟合获得准确参数C_j和τ_j，但短时间分解曲线可能导致参数不稳定，当前解决方案是选取足够长的曲线（通过多次拟合至参数稳定）。该方法仅适用于单组分气体系统，多组分系统需进一步研究。未来需在更多案例中验证并改进优化方法。

**结论**
本研究对传统的步进加热等容法（SHI method）进行了优化，以快速确定沉积物中气体水合物的相平衡数据。在优化方法中，推导了一个指数函数来表征每个温度步长下的水合物分解曲线（气体压力随时间变化），相应的关键参数可通过等容条件下步进加热的非稳态水合物分解测试进行校准，然后计算相平衡压力和温度。为验证优化方法的有效性，采用优化方法进行了6个案例的相平衡数据测定。优化方法得到的结果与传统方法高度一致，证明优化方法能够合理有效地确定沉积物中水合物的相平衡数据。更重要的是，与传统方法相比，优化方法可以节省大量时间。因此，优化方法高效且省时，但仍需在未来的研究中进一步考虑一些细节。