模拟海水环境下MICP固化钙质砂的力学特性

李昊; 唐朝生; 刘博; 吕超; 程青; 施斌

doi:10.11779/CJGE202010019

模拟海水环境下MICP固化钙质砂的力学特性 English Version

南京大学地球科学与工程学院,江苏　南京　210023

基金项目:

国家杰出青年科学基金项目 41925012

国家自然科学基金项目 41902271

国家自然科学基金项目 41772280

国家自然科学基金项目 41572246

江苏省自然科学基金项目 BK20171228

江苏省自然科学基金项目 BK20170394

中央高校基本科研业务费专项资金项目 2018-2022

详细信息

作者简介:
李昊(1996—),男,硕士研究生,主要从事微生物地质工程研究工作。E-mail：lihao0531@smail.nju.edu.cn

通讯作者:
唐朝生, E-mail：tangchaosheng@nju.edu.cn

中图分类号: TU441
计量
- 文章访问数: 474
- HTML全文浏览量: 26
- PDF下载量: 386
出版历程
- 收稿日期: 2019-10-07
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-09-30

Mechanical behavior of MICP-cemented calcareous sand in simulated seawater environment

School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China

摘要

摘要: 钙质砂广泛分布于近海大陆架、海岸带及大洋岛礁,具有低强度和易破碎的特点,且长期受风浪侵蚀。为了改善钙质砂的力学特性,提出了基于微生物诱导碳酸钙沉积（MICP）固化技术,并在模拟海水环境中开展了一系列固化试验,测试了试样的无侧限抗压强度,同时与淡水环境下获得的试验结果进行了对比分析。此外,试验通过设置不同尿素浓度（0.25,0.5,1.0和1.5 mol/L）的胶结液,探究了尿素浓度对MICP固化钙质砂力学性能的影响及机理。研究结果表明：①MICP技术能够适用海洋环境,且对钙质砂的加固效果比淡水环境更佳,在海水环境中固化后试样的无侧限抗压强度相比淡水环境得到了成倍提高；②海水的弱碱性环境对提升脲酶菌活性和MICP固化效果具有积极作用；③尿素浓度对MICP的固化效果有重要影响,试样的无侧限抗压强度随尿素浓度的增加呈先增加后减小趋势,本次试验发现最优尿素浓度为1.0 mol/L；④MICP固化试样的无侧限抗压强度与微生物诱导生成的碳酸盐含量呈正相关关系。
- 海水环境 /
- MICP /
- 钙质砂 /
- 力学特性 /
- 碳酸盐含量 /
- 尿素浓度
Abstract: The calcareous sand is widely distributed in offshore continental shelves, coastal zones and oceanic reefs. It has the characteristics of low strength and fragility and has been eroded by wind and waves for long time. In order to improve the mechanical properties of the calcareous sand, a microbial-induced calcium carbonate precipitation (MICP) method is proposed, a series of cementing processes are reproduced in simulated seawater environment, and the unconfined compressive strength of the specimens is tested. At the same time, the results are compared with those obtained in fresh water environment. In addition, the effect and mechanism of urea concentration on the mechanical properties of MICP-treated calcareous sand are investigated by setting cementing solutions with different urea concentrations (0.25, 0.5, 1.0 and 1.5 mol/L). The results show that: (1) The MICP technology can be applied to marine environment, the reinforcement effect of the calcareous sand is better than that of fresh water environment, and the unconfined compressive strength of specimens after solidification in seawater is doubled compared with that in fresh water environment. (2) The alkaline environment of seawater has a positive effect on the activity of urea bacteria and the cementing effect of MICP. (3) The urea concentration has an important impact on the cementing effect of MICP, and the unconfined compressive strength of the specimen first increases and then decreases with the increase of urea concentration. It is found that the optimal urea concentration is 1.0 mol/L. (4) The unconfined compressive strength of MICP-treated specimens is positively correlated with the content of microbial-induced calcium carbonate.
- seawater environment /
- microbial-induced calcium carbonate precipitation /
- calcareous sand /
- mechanical behavior /
- carbonate content /
- urea concentration

HTML全文

0. 引言

天然气水合物是一种笼形晶状化合物（图 1），具有能量密度高、储量巨大等特点，其中97%以上分布在海底沉积物中^[1-2]。据估算，全球天然气水合物中含碳量相当于全球已探明传统化石燃料总碳量的两倍，被普遍认为是21世纪最具潜力的新型清洁能源之一^[2]。目前世界很多国家都开展了水合物勘探和试采工作，例如美国JIP计划^[3]、日本MHR & DP项目^[4]，挪威GANS项目^[5]。另外，俄罗斯、德国、韩国、印度等国家都相继开展了水合物相关的研究项目和开发计划^[6]。中国在水合物领域的研究起步较晚，但取得了重要成果。在南海北部陆坡圈定了6个水合物成矿远景区，19个成矿区带，25个有利区块和24个钻探目标区^[7-8]。特别是在神狐、东沙等海域钻探获取了水合物实物样品，发现了多种类型的高品位水合物储层，锁定了试采目标^[9]。并于2017年和2020年成功完成了两次水合物试采工作，取得了突破性进展^[10-11]。

图 1 水合物实物与分子结构图^[2]

Figure 1. Physical and molecular structural diagrams of hydrate

下载: 全尺寸图片幻灯片

目前水合物领域的研究主要集中在资源储量勘探评估、水合物开发技术和水合物地层稳定性评价3个方面。其中水合物地层稳定性评价是水合物安全可持续生产的前提条件，特别是考虑开采影响下水合物储层的动态响应及海底斜坡稳定性问题。通过海洋工程地质调查、室内模型试验和数值模拟研究发现，含水合物储层的海底斜坡失稳破坏，往往与水合物分解后产生的超孔隙压力有关^[12-14]。因此，近年来专家学者对与水合物分解相关的超孔压计算模型进行了深入研究。例如，Nixon等^[15]提出一个数学模型，通过假定水合物的分解比例，量化了不排水条件下水合物分解产生的超孔压。在此基础上，Xiong等^[16]考虑甲烷的溶解度对孔隙压力的影响，对模型进行了改进。Zhang等^[17]通过RK方程计算了真实甲烷气体的压缩系数，并将此系数引入到上述模型。Chen等^[18]利用此模型，通过考虑水合物分解过程中压力和温度的平衡，计算了水合物分解的最大超孔隙压力。此外，刘锋等^[19]、Grozic^[20]、Zhu等^[21]进行了相关研究。然而，这些研究中所采用的孔隙压力模型很少考虑水合物储层的热-水-化(THC)耦合作用，导致无法准确获取瞬态超孔隙压力及其演化特征，因此很难准确评价海底斜坡的实时动态稳定性。

本文考虑水合物注热开采的影响，发展了一种热-水-化（THC）耦合作用下海底斜坡动态稳定性简化评价方法，并自主开发了相应的数值计算程序，建立了水合物饱和度与超孔隙压力间的内在关系；进而获取了多物理场耦合作用下水合物分解引发的瞬态超孔压演化特征及储层的动态响应；并利用改进的极限平衡法对海底斜坡稳定性进行动态评价；最后以南海北部陆坡水合物试采区典型斜坡为目标，研究了水合物注热开采影响下水合物储层动态响应及实时稳定性演化。

1. 考虑THC耦合作用的超孔压确定方法

1.1 THC耦合方程

含水合物沉积物可视为由多组分（土骨架、水、水合物）和多相（固、液、气）组成的动态系统（图 2）。根据地质调查^[22]，海底水合物大多储存于相对封闭的沉积地层中，其形成由沉积物沉积、孔隙流体流动及底部甲烷供给等因素决定，且这类储层中一般缺乏明显的横向变化，因此只需考虑竖向的流体运移和能量迁移。基于此，采用Xu等^[23]提出的一维简化THC耦合模型描述海底水合物地层的赋存特性。该模型通过求解质量、能量和甲烷守恒方程，给出了水合物储层边界的计算公式，并推导了初始地层温度–压力分布和水合物聚集的表达式。

图 2 水合物沉积物多组分多相孔隙介质示意图

Figure 2. Schematic diagram of multi-component multi-phase porous sediments containing gas hydrate

下载: 全尺寸图片幻灯片

首先，根据质量守恒给出一维简化情况下流体运移的控制方程：

，

$- \frac{{k{\rho _{\text{l}}}}}{{{\mu _{\text{l}}}}}\left( {\frac{{{\text{d}}P}}{{{\text{d}}z}} + {\rho _{\text{l}}}g} \right) = {q_{\text{f}}}\text{，}$

(1)

式中，k为渗透率，ρ_l，μ_l分别为液体相的密度和动态黏度，P为地层压力，z为空间坐标，g为重力加速度，q_f为总质量通量。

其次，水合物储存系统的能量迁移可通过能量守恒方程描述，其简化模式为

，

$- \lambda \frac{{{\text{d}}T}}{{{\text{d}}z}}{\text{ + }}{q_{\text{f}}}{C_{{\text{pl}}}}T{\text{ = }}{q_{\text{e}}}\text{，}$

(2)

式中，λ为有效导热系数，T为地层温度，C_pl为水合物的比热容，q_e为总能量通量。

最后，储存系统中的甲烷运移可由甲烷介质的质量守恒方程描述：

，

$- \phi {K_{\text{m}}}{S_{\text{l}}}\frac{{{\text{d}}{M_{\text{l}}}}}{{{\text{d}}z}}{\text{ + }}{q_{\text{f}}}{M_{\text{l}}}{\text{ = }}{q_{\text{m}}}\text{，}$

(3)

式中，ϕ为沉积物孔隙度，K_m为扩散系数，S_l为液体相饱和度，M_l为液体相中甲烷浓度，q_m为与深度相关的甲烷通量。

以上3组基本控制方程间的耦合作用，其具体的建立过程为对式（1），（2）进行积分，可得到关于压力剖面P和地层温度剖面T的函数表达式：

，

$P - {P_0} = - \left( {\frac{{{q_{\text{f}}}{\mu _{\text{l}}}}}{{k{\rho _{\text{l}}}}}{q_{\text{f}}} + {\rho _{\text{l}}}g} \right)z\text{，}$

(4)

$T= \begin{cases}-\left(\xi-T_0\right) \mathrm{e}^{\frac{q_{\mathrm{e}}}{\lambda \cdot \xi} z}+\xi & \left(q_{\mathrm{f}} \neq 0\right) \\ -\frac{q_{\mathrm{e}}}{\lambda} z+T_0 & \left(q_{\mathrm{f}}=0\right)\end{cases}\text{，}$

(5)

式中，ξ为参数，ξ=q_e/(q_fC_pl)，P₀，T₀为海底面处的压力和温度。联立式（4），（5），可以得到任意深度z处的压力与温度之间的关系式。

水合物的形成需要适当的温压条件，地层温度剖面和水合物相平衡曲线共同确定了水合物的理论稳定区间，如图 3所示。Sloan等^[24]通过试验测定了甲烷水合物的相平衡曲线，表达为

$T_{\mathrm{diss}}=\left\{\begin{array}{ll} 9.6349 \ln P+197.65 & (P \geqslant 2576.74 \mathrm{kPa}) \\ 0.0158 P+232.7 & (P<2576.74 \mathrm{kPa}) \end{array},\right.$

(6)

图 3 水合物理论稳定区域与生成区域的空间关系示意图^[26]

Figure 3. Schematic diagram of spatial relationship between hydrate stability zone and hydrate formation zone

下载: 全尺寸图片幻灯片

式中，T_diss为水合物分解的相平衡温度。

将式（4）~（6）联立，可得到水合物理论稳定区域的上边界z_st和下边界z_sb：

${z_{{\text{st}}}} = 0 \text{，}$

(7)

$z_{\mathrm{sb}}=\left\{\begin{array}{ll} -\frac{\lambda}{q_{\mathrm{f}} C_{\mathrm{pl}}} \ln \left(\frac{q_{\mathrm{e}}-q_{\mathrm{f}} C_{\mathrm{pl}} T_0}{q_{\mathrm{e}}-q_{\mathrm{f}} C_{\mathrm{pl}} T_{\mathrm{diss}}}\right) & \left(q_{\mathrm{f}} \neq 0\right) \\ -\frac{\lambda}{q_{\mathrm{e}}}\left(T_{\mathrm{diss}}-T_0\right) & \left(q_{\mathrm{f}}=0\right) \end{array},\right.$

(8)

式中，上边界z_st和下边界z_sb见图 3右侧标注。另外，水合物的生成不仅需要适宜的温度和压力，还取决于孔隙中甲烷的可用量。只有当孔隙中甲烷含量大于或等于甲烷溶解度时，才可能在孔隙中形成水合物，如图 3曲线（c）和曲线（d）所示。在水合物生成区域内部，多余的甲烷会与水结合形成水合物，因此孔隙甲烷含量等于甲烷溶解度，可由下式描述：

${M_{{\text{sl}}}}(T) = {M_{{\text{sl,diss}}}}({T_{{\text{diss}}}},P,S) \cdot \exp \left( {\frac{{T - {T_{{\text{diss}}}}}}{\beta }} \right)\text{，}$

(9)

式中，M_sl为任意温度、压力及孔隙中盐浓度影响下的甲烷溶解度，M_{sl, diss}(T_diss, P, S)为甲烷溶解度关于相平衡温度T_diss、压力P和盐度S的函数，由试验测定^[25]，β=14.4 ℃为经验常数。

而在水合物生成区域上部，孔隙中甲烷含量小于其溶解度，甲烷介质以溶质赋存于孔隙中，可由式（3）描述。为确定水合物生成区域的上边界z_ot和下边界z_ob(图 3右侧标注)，在0→z_ot区间内对式（3）进行积分，并辅以式(4)，(5)，(9)可得

$z_{\mathrm{ot}}= \begin{cases}-\frac{\phi K_{\mathrm{m}}}{q_{\mathrm{f}}} \ln \left(\frac{q_{\mathrm{mt}}-q_{\mathrm{f}} M_0}{q_{\mathrm{mt}}-q_{\mathrm{f}} M_{\mathrm{sl}}\left(z_{\mathrm{ot}}\right)}\right) & \left(q_{\mathrm{f}} \neq 0\right) \\ -\frac{\phi K_{\mathrm{m}}}{q_{\mathrm{mt}}}\left(M_{\mathrm{sl}}\left(z_{\mathrm{ot}}\right)-M_0\right) & \left(q_{\mathrm{f}}=0\right)\end{cases} \text{，}$

(10)

${z_{{\text{ob}}}}: - \phi {K_{\text{m}}}{\left( {\frac{{{\text{d}}{M_{{\text{sl}}}}}}{{{\text{d}}z}}} \right)_{z = {z_{{\text{ob}}}}}}{\text{ + }}{q_{\text{f}}}{M_{{\text{sl}}}}({z_{{\text{ob}}}}){\text{ = }}{q_{\text{m}}} \text{，}$

(11)

式中，M₀为海底面处的甲烷含量，q_mt为z_ot处的甲烷通量。

综合考虑式（7），（8），（10），（11），同时满足水合物理论稳定区域条件和水合物生成区域条件即可得出水合物真实赋存区域，其上边界THZ和下边界BHZ如图 3所示。

进一步，考虑水合物形成的时间尺度和空间分布与甲烷供应有关，可给出确定整个水合物储层初始水合物饱和度的表达式^[23]。甲烷水合物填充孔隙度ϕ到任意平均饱和度S_h所需的时间t₀近似为

${t_0} \approx {S_{\text{h}}}{\left( {\frac{{\partial {S_{\text{h}}}}}{{\partial t}}} \right)^{ - 1}}\text{，}$

(12)

式中，t为时间，∂S_h/∂t为水合物的积累速率，

$\frac{{\partial {S_{\text{h}}}}}{{\partial t}} = \frac{1}{{\phi ({\rho _{\text{h}}}{M_{\text{h}}} - {\rho _{\text{l}}}{M_{{\text{sl}}}})}}\left( {\phi {K_{\text{m}}}\frac{{{{\text{d}}^2}{M_{{\text{sl}}}}}}{{{\text{d}}{z^2}}} - {q_{\text{f}}}\frac{{{\text{d}}{M_{{\text{sl}}}}}}{{{\text{d}}z}}} \right) \text{，}$

(13)

其中，ρ_h，M_h为水合物密度和水合物中甲烷的质量分数。综上，由式（1）~（13）构建了一套完整的控制方程，可用于计算地层初始状态下的温度、压力剖面和水合物饱和度剖面。

1.2 超孔压计算

水合物储层注热开采会破坏地层中温度的动态平衡，一旦地层中的温度T超过水合物的相平衡温度T_diss，水合物会产生分解。地层温度剖面T的演化可以用一维热传导方程描述：

$\frac{{\partial T(z,t)}}{{\partial t}} = \kappa \frac{{{\partial ^2}T}}{{\partial {z^2}}} + s(z,t)\text{，}$

(14)

式中，κ为地层有效热扩散系数，假设为常量，这是因为水合物的导热系数接近于水（相差 < 10%）^[27]，s(z，t)为源项，用以考虑可能的热量增加和消耗。需要注意的是，在地层升温过程中（见图 4（a），A→B），一旦温度超出水合物的稳定区域，水合物会发生分解。这其中包含水合物的相变过程，因此需要消耗一定的潜热。这意味着当地层温度T_{i, j}（下标i代表时间步，j代表某地层深度）达到T_diss后，随着热量的增加，温度会保持T_diss不变，这会延迟温度在地层中的传播（见图 4（b））。

图 4 注热开采导致的地层温度变化路径示意图

Figure 4. Schematic diagram of change path of stratum temperature caused by hydrate heat injection exploitation

下载: 全尺寸图片幻灯片

为了计算潜热的消耗，采用热积分法^[28]进行简化。具体实施过程：首先利用有限差分法求解式(14)，计算在任意时间步i时地层深度j的温度T_{i, j}。在水合物分解区域将此温度修正为T_diss，进而会产生一个温度增量dQ_{i, j}=C_ph(T_ij-T_diss)，此温度增量会根据时间步和深度更新。将每一时间步所产生的热量累计求和，可以计算出地层深度j增加的总增量∑dQ_{i, j}。一旦累积的总热量超过水合物的分解潜热，则意味着相变完成，地层温度重新按照式（14）升高。在此过程中，位于深度j处，时刻i时的水合物饱和度S_{h i, j}为

${S_{{\text{h}} i,j}} = {S_{{\text{h}} j}} \times \frac{{{L_{\text{h}}} - \Sigma {\text{d}}{Q_{i,j}}}}{{{L_{\text{h}}}}}\text{，}$

(15)

式中，S_{h, j}为深度j处的初始饱和度，L_h为水合物的分解潜热。

根据式（15）可得出在地层温度变化过程中水合物的分解量，

$\Delta {S_{{\text{h }}i{\text{,}}j}} = \Sigma ({S_{{\text{h }}i{\text{,}}j}} - {S_{{\text{h }}i - 1{\text{,}}j}})。$

(16)

最后，利用气体非逃逸方法计算温度升高过程中瞬态超孔隙压力及其演化。此方法由Nixon等^[15]首先提出，后续不少学者对该方法进行了完善^[16-21]。这里在前人工作的基础上，充分考虑水合物赋存系统的THC耦合作用，获取水合物实时分解量(式(16))，并引入到该方法中，实现升温后地层中的真实超孔隙压力演化。其中做了如下假设：①依据中国南海水合物试采区的地质调查结果，其水合物储层沉积物以细粒为主，渗透系数较小^[9]，因此忽略水合物分解后气体的渗流；②水合物赋存区的孔隙水在水合物分解前已处于甲烷饱和状态，因此忽略分解后生成甲烷在孔隙水中的溶解；③相对于气体的压缩量其余介质的压缩量很小，忽略不计。

基于以上假设，单位水合物分解后沉积物的体积变化等于生成气体的体积∆V_unit，即

$\Delta {V_{{\text{unit}}}} = {V_{\text{g}}} + {V_{{\text{wh}}}} - 1\text{，}$

(17)

式中，V_g为水合物分解产生的气体体积，V_g=164.6P_atmT_nat/(T_atmP_nat)，P，T分别为地层压力和温度，下标atm和nat分别为标准条件和地层真实条件，V_wh为水合物分解产生水的体积。故沉积物总体积变化为

$\Delta V = \phi \Delta {S_{{\text{h}}i,j}}V(\frac{{164.6{p_{{\text{atm}}}}{T_{{\text{nat}}}}}}{{{T_{{\text{atm}}}}{p_{{\text{nat}}}}}} + {V_{{\text{wh}}\theta }} - 1)。$

(18)

根据有效应力原理，超孔隙压力与沉积物变形有关，

$\begin{aligned} u_{\mathrm{e} i, j} & =-E_{\mathrm{s}} \Delta \varepsilon_{\mathrm{v}}=E_{\mathrm{s}} \frac{\Delta V}{V} \\ & =E_{\mathrm{s}} \phi \Delta S_{\mathrm{h} i, j}\left(\frac{164.6 p_{\text {atm }} T_{\text {nat }}}{T_{\mathrm{atm}} p_{\text {nat }}}+V_{\mathrm{wh}}-1\right), \end{aligned}$

(19)

式中，u_{e i, j}为时刻i时深度j处的超孔隙压力，ε_v为沉积物体应变，E_s为压缩模量，其余参数的物理意义同前。综上，可以得到地层升温后考虑THC耦合的超孔隙压力及其演化特征。

需要指出的是，在超孔压计算中做了一些合理性假设，忽略了一部分气体渗流和其余介质的压缩量，故计算出的超孔压会出现大于实际值的情况。另外，水合物分解产生的孔压会增大地层压力P_nat，进而反向影响水合物的分解，使水合物分解减缓，上述超孔压计算方法为单向耦合过程，即不考虑这种超孔压的反向耦合效应，使得海底斜坡的稳定性评价偏于安全。

1.3 方法验证

为了检验上述超孔压计算方法的正确性，复刻了Reagan等^[29]基于TOUGH+HYDRATE的模拟案例。TOUGH+HYDRATE为目前关于水合物分解模拟较为普遍的商业软件之一，并得到了广泛应用^[30-31]。此案例描述了浅层水合物赋存的典型场景，其中水深为320 m，海底温度为0.4℃，地温梯度为3℃/100m，模拟了海底水合物层100 a内，以每年0.03℃上升工况下水合物饱和度的动态响应。

对比文献[29]中设定升温工况下的计算结果，可见不同时刻（0，10，100 a）水合物饱和度的演化曲线具有较好的一致性，如图 5所示。根据式（19）可知，超孔压计算结果的合理性取决于水合物饱和度变化量∆S_{hi, j}的正确性，上述饱和度演化曲线的一致性很好的检验了本文提出的超孔压计算方法的正确性。

图 5 水合物饱和度剖面演化曲线对比图

Figure 5. Comparison of evolution curves of hydrate saturation profile

下载: 全尺寸图片幻灯片

2. 海底斜坡动态稳定性评价方法

极限平衡方法(LEM)被广泛用于海底斜坡的稳定性评价，尤其对于长度远大于沉积层厚度的无限坡。海洋工程地质调查表明，大多数水合物储层与海底面平行分布，贮藏于一定厚度的不透水覆盖层下方^[10]，因此将海底斜坡概化为由上覆层、水合物储层和下卧层组成的无限斜坡（图 6（a））。注热开采条件下，水平井附近的水合物开始分解，将水合物的分解前缘视为潜在滑动面，并以此对海底斜坡进行稳定性评价，如图 6所示。

图 6 水合物开采影响下典型海底斜坡稳定性评价示意图^[32]

Figure 6. Schematic diagram of stability evaluation of typical submarine slope under influences of hydrate exploitation

下载: 全尺寸图片幻灯片

根据极限平衡方法，海底斜坡的稳定性安全系数等于滑动面上的抗剪强度和剪应力的比值，

${F_{\text{s}}} = \frac{{\int {{{\vec \tau }_{\text{r}}}{\text{d}}l} }}{{\int {{{\vec \tau }_{\text{d}}}{\text{d}}l} }}\text{，}$

(20)

式中， ${\vec \tau _{\text{r}}}$ ， ${\vec \tau _{\text{d}}}$ 分别为滑动面处沉积物的抗剪强度和剪应力，dl为单元滑动面长度，通过建立垂直于l轴的力系平衡，可以导出：

${\vec \tau _{\text{r}}} = c' + {\sigma '_0}\tan \varphi '{\text{ = }}c' + (\sigma - {u_{\text{e}}})\tan \varphi ' \text{，}$

(21)

${\vec \tau _{\text{d}}} = \gamma 'H\sin \alpha \cos \alpha \text{，}$

(22)

式中，φ′为有效内摩擦角， ${\sigma '_0}$ 为上覆有效应力 ${\sigma '_0}$ = $\sigma - {u_{\text{e}}}$ ， $\sigma$ 为总上覆压力 $\sigma$ = $\gamma '$ Hcos² $\alpha$ ， $\gamma '$ 为沉积物有效重度，H为滑面深度，c′为沉积物有效黏聚力。水合物分解会引起沉积物胶结强度弱化，进而导致黏聚力降低^{[31, 33]}，这里假设黏聚力随着水合物的分解呈线性折减：

${c'_i} = - ({c'_1} - {c'_2})\frac{{\Delta {S_{{{\text{h}}^{}}i,j}}}}{{{S_{{{\text{h}}^{}}i,j}}}}{\text{ + }}{c'_1} \text{，}$

(23)

式中， ${c'_i}$ 为i时刻的黏聚力， ${c'_1}$ ， ${c'_2}$ 分别为初始黏聚力和最终黏聚力。

将式（21），（22）代入式（20）中，可得出海底斜坡不同深度、不同时刻的瞬态稳定性安全系数F_{s i, j}为

${F_{{{\text{s}}_{}}i,j}} = \frac{{{{c'}_i}}}{{\gamma 'H\sin \alpha \cos \alpha }} + \left( {1 - \frac{{{u_{{{\text{e}}_{}}i,j}}}}{\sigma }} \right)\frac{{\tan \varphi '}}{{\tan \alpha '}} 。$

(24)

3. 注热开采下海底斜坡稳定性实例分析

3.1 研究区域

选取南海北部陆坡白云凹陷作为研究区，区域整体位于珠江口盆地珠二坳陷内，走向近东西，面积大于2.0×10⁴ km²，是南海北部最大的一个深水凹陷，其北部与番禺低隆起相接，南端是珠江口盆地的南部隆起带^[34]。区域内发育有大量现代海底峡谷（见图 7），沉积物失稳频发。中国地质调查局于2019年10月—2020年4月在此区域进行了第二次天然气水合物试采，实现了连续产气30 d，总产气量86.14×10⁴ m³，大大提高了日产气量和产气总量^[11]。此次试采首次提出了利用水平井开采水合物的工作思路，可有效克服直井开采大压差导致二次水合物和冰生成的风险。由于复杂地质环境影响，开采引起的水合物分解可能在海底斜坡内部产生潜在滑动面，即使在坡度很小的情况下也会引发斜坡破坏。特别是考虑到海底浅表层沉积物强度低^[35-36]，易在外部作用干扰下发生土体失稳，进而可能导致海底滑坡等地质灾害，威胁海洋开采平台、锚固基础和海底管线等水下设施的安全^[37-39]。这里主要探讨水平井注热开采工况下水合物储层的超孔压演化过程及海底斜坡动态稳定性评价。

图 7 研究区位置及区域海底地形条件^{[11, 40]}

Figure 7. Location of study area and topography of regional seabed

下载: 全尺寸图片幻灯片

为服务于第二次水合物试采井的轨迹优选，在试采水平井靶点位置附近进行了先导井测试工作，并利用测井得到的P波速度和S波速度预测了水合物饱和度S_h和游离气饱和度S_g分布（图 8）^[41]。分析可见，该区域水合物储层位于海底面以下深度207.8~253.4 m，厚度为45.6 m，平均孔隙度为37.3%，平均水合物饱和度为31.0%^[41]，全井段岩性主要是泥质粉砂，渗透系数较小。另外，海洋地质调查显示此区域水深1000~1500 m，平均地形坡度为3°，海底面温度为3.3~3.7℃，地热梯度为45~67 ℃/km^[9]。基于开采井原位测试数据，构建二维海底斜坡模型，由上覆层、水合物层和下卧层组成，如图 6（a）所示。为达到最大传热效率，将水平井设置在水合物储层中间位置，即H=230.6 m处。设定由于注热影响，地层温度以0.1 K/d的速度匀速增加，总升高温度为10 K，当温度达到设定值后保持不变。模型所用具体参数见表 1。

图 8 先导井测井响应柱状图^[41]

Figure 8. Histogram logging response of pilot well

下载: 全尺寸图片幻灯片

表 1 模型参数^{[9, 11, 26, 41-42]}

Table 1. Model parameters

参数		取值	单位
坡体参数	z (深度)	—	m
	α (坡度)	3	°
	H(水平井位置)	230.6	m
地层参数	T₀(海底温度)	3.5	℃
	G/(地热梯度)	56	K/km
	S_h (初始饱和度)	31	%
	k (渗透率)	2.38	mD
	κ (有效热扩散系数)	4.5×10⁻⁷	m²/s
	λ (有效导热系数)	3.66	W/(m·K)
	c₁′ (初始黏聚力)	800	kPa
	c₂′ (最终黏聚力)	30	kPa
	φ′ (内摩擦角)	8.03	°
	ρ (体密度)	1650	kg/m³
	ϕ (孔隙度)	0.37	m³/m³
通用参数	ρ_l/(液相密度)	1024	kg/m³
	ρ_h (水合物密度)	930	kg/m³
	C_pl (液体比热)	4.18×10³	J/(kg·K)
	C_ph/(水合物比热)	2.16×10³	J/(kg·K)
	μ_l (液体动态黏度)	8.87×10⁻⁴	kg/(m·s)
	S (液体盐度)	3.5	%
	L_h/(水合物潜热)	56552-16.8/T	J/kg

下载: 导出CSV

| 显示表格

3.2 水合物储层动态响应分析

图 9为模拟工况下地层温度剖面的演化曲线，其中THZ和BHZ分别表示水合物储层的上顶面和下底面。初始地温曲线（0 a）和水合物相平衡曲线（T_diss）所围成的阴影部分为水合物的理论稳定区间（下底面为z_sb）。可以看到，由于水平井中热量的注入，水合物储层中部的温度开始升高，并逐渐向上下扩散，一旦温度超过水合物的相平衡温度，水合物开始分解（图 10）。可以观察到，在曲线T_diss附近地层的温度出现局部的降低，这是因为水合物分解是吸热过程，消耗了一定量的潜热，潜热的消耗会造成温度传导的延迟。

图 9 地层温度剖面演化曲线

Figure 9. Evolution curves of stratum temperature profile

下载: 全尺寸图片幻灯片

图 10 升温过程中水合物储层动态响应

Figure 10. Dynamic responses of hydrate reservoir during heating process

下载: 全尺寸图片幻灯片

图 10（a）为水合物饱和度的演化曲线。

可以看到，随着热量的注入储层中部水合物最先开始分解，饱和度逐渐降为0，在地层中形成一条分解带，且随着热量的传递，上下两侧的水合物也开始分解，分解带不断延伸。图 10（b）为水合物分解前缘的演化曲线，可以更加直观地看到水合物的分解及传播过程，其中红色曲线代表水合物分解前缘，黑色曲线为水合物储层的上下界面。可以看到，分解前缘最先产生于储层中部，并不断向两侧移动，直到与水合物储层上下界面重合，意味着水合物全部分解；在储层中部，水合物分解前缘的演化非常迅速，而在储层两侧则相对缓慢。这是因为储层中部靠近水平井热源，此处的水合物分解反应更加剧烈，而如上所述，水合物分解会导致地层中热量的延迟传导，故储层两侧的水合物分解需要更长的时间。

3.3 海底斜坡动态稳定性分析

图 11显示海底斜坡安全系数沿深度的变化曲线。可以看到，初始状态下海底斜坡处于稳定状态，且由于水合物的胶结作用，增加了沉积物的强度，使得初始状态下水合物储层处的安全系数明显高于地层其它位置。但随着水合物分解，海底斜坡的安全系数迅速减小，并达到临界状态。一方面是由于胶结作用的丧失，沉积物强度迅速降低，另一方面水合物分解产生的超孔隙压力减小了地层的有效应力，随着水合物分解区域的不断延伸，不稳定区域也不断扩大。

图 11 海底斜坡安全系数沿深度的变化

Figure 11. Variation of safety factor of submarine slope along depth

下载: 全尺寸图片幻灯片

图 12为海底斜坡最小安全系数的演化曲线。可以看到，经过96 d的升温后，地层中的温度达到分解条件，水合物开始分解，此时海底斜坡安全系数开始下降，在开采358 d后安全系数降为1.0，意味着海底斜坡达到临界状态。此时为了后期的持续安全开采，需要调整开采方案。例如：调整施工参数，配合降压法，既可以有效降低地层中的孔隙压力，又可以提高开采效率。

图 12 海底斜坡最小安全系数随时间的变化曲线

Figure 12. Variation of minimum safety factor of submarine slope with time

下载: 全尺寸图片幻灯片

4. 结论

（1）考虑注热开采影响下海底水合物储层的THC耦合作用，推导了水合物分解引起的瞬态超孔隙压力求解格式，并基于极限平衡分析理论建立了海底斜坡动态稳定性简化评价方法。

（2）水平井注热开采影响下，水合物分解前缘最先形成于储层中部，并不断向上、下两侧延伸；由于水合物潜热消耗的影响，导致热量延迟传导，使储层上下两侧的水合物分解需要更长时间。

（3）水合物试采区典型海底斜坡初步分析表明，由于水合物的胶结作用，初始状态下水合物储层处的安全系数偏高；升温96 d后水合物开始分解，海底斜坡安全系数迅速降低，并于358 d后到达临界状态；为保证水合物的持续安全开采，提出调整开采方案。

（4）本文基于无限坡理论分析了水合物注热开采影响下海底斜坡的动态稳定性问题，而实际含水合物层的海底斜坡失稳破坏或滑坡模式，可能表现为弧形、盘式或拉裂刺穿等破坏形式^[13]，未来应加强不同破坏模式的海底斜坡动态稳定性评价方法的深入研究。

衷心感谢中国科学院武汉岩土力学研究所汪稔研究员为本次研究提供的钙质砂。

图 1 无侧限抗压强度试验砂柱模型及试样

Figure 1. Molds and specimens for unconfined compression strength tests

下载: 全尺寸图片幻灯片

图 2 5组试样的应力-应变曲线

Figure 2. Stress-strain curves of five sand samples

下载: 全尺寸图片幻灯片

图 3 模拟海水环境下不同尿素浓度的无侧限抗压强度

Figure 3. Unconfined compressive strengths under seawater with different urea concentrations

下载: 全尺寸图片幻灯片

图 4 无侧限抗压试验后0.25与1.0 mol/L尿素浓度试样对比

Figure 4. Comparison of unconfined compressive strengths under urea concentrations of 0.25 and 1.0 mol/L

下载: 全尺寸图片幻灯片

图 5 不同胶结环境无侧限抗压强度对比

Figure 5. Comparison of unconfined compressive strengths in different cementing environments

下载: 全尺寸图片幻灯片

图 6 不同试样的碳酸盐含量与无侧限抗压强度

Figure 6. Carbonate contents and unconfined compressive strengths of different specimens

下载: 全尺寸图片幻灯片

图 7 不同碳酸盐含量对应的无侧限抗压强度

Figure 7. Unconfined compressive strengths corresponding to different carbonate contents

下载: 全尺寸图片幻灯片

表 1 不同试样碳酸盐沉淀含量

Table 1 Carbonate precipitation content in different samples

碳酸盐含量/%	淡水环境	海水环境
碳酸盐含量/%	S0	S1	S2	S3	S4
平行样1	17.25	6.97	23.57	27.95	27.51
平行样2	19.36	10.11	24.26	26.13	25.42
平行样3	17.94	—	22.39	26.84	26.28
平均值	18.18	8.54	23.41	26.97	26.40
注：S1-平行样3未完成固化,无碳酸盐含量及无侧限抗压强度数据。

下载: 导出CSV

参考文献(63)

[1]	ALBA J L, AUDIBERT J M. Pile design in calcareous and carbonaceous granular materials, an historic overview[C]//The Second International Conference on Engineering for Calcareous Sediments, 1999, Manama: 29-44.
[2]	陈海洋, 汪稔, 李建国, 等. 钙质砂颗粒的形状分析[J]. 岩土力学, 2005, 26(9): 1389-1392. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200509007.htm CHEN Hai-yang, WANG Ren, LI Jian-guo, et al. Grain shape analysis of calcareous soil[J]. Rock and Soil Mechanics, 2005, 26(9): 1389-1392. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200509007.htm
[3]	蒋明镜, 吴迪, 曹培, 等. 基于SEM图片的钙质砂连通孔隙分析[J]. 岩土工程学报, 2017, 39(增刊1): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2017S1002.htm JIANG Ming-jing, WU Di, CAO Pei, et al. Connected inner pore analysis of calcareous sands using SEM[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(S1): 1-5. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2017S1002.htm
[4]	刘崇权, 汪稔. 钙质砂物理力学性质初探[J]. 岩土力学, 1998, 19(1): 32-37, 44. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX199801005.htm LIU Chong-quan, WANG Ren. Preliminary research on physical and mechanical properties of calcareous sand[J]. Rock and Soil Mechanics, 1998, 19(1): 32-37, 44. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX199801005.htm
[5]	李彦彬, 李飒, 刘小龙, 等. 颗粒破碎对钙质砂压缩特性影响的试验研究[J]. 工程地质学报, 2020, 28(2): 352-359. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202002016.htm LI Yan-bin, LI Sa, LIU Xiao-long, et al. Study on effect of particle breakage on compression properties of calcareous sands with oedometertests[J]. Journal of Engineering Geology, 2020, 28(2): 352-359. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202002016.htm
[6]	孙吉主, 黄明利, 汪稔. 内孔隙与各向异性对钙质砂液化特性的影响[J]. 岩土力学, 2002, 23(2): 166-169. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200202007.htm SUN Ji-zhu, HUANG Ming-li, WANG Ren. Influence of inner pore and anisotropy on liquefaction characteristics of calcareous sand[J]. Rock and Soil Mechanics, 2002, 23(2): 166-169. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX200202007.htm
[7]	孙吉主, 汪稔. 钙质砂的颗粒破碎和剪胀特性的围压效应[J]. 岩石力学与工程学报, 2004, 23(4): 641-644. doi: 10.3321/j.issn:1000-6915.2004.04.021 SUN Ji-zhu, WANG Ren. Influence of confining pressure on particle breakage and shear expansion of calcareous sand[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(4): 641-644. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.04.021
[8]	张丙树, 顾凯, 李金文, 等. 钙质砂破碎过程及其微观机制试验研究[J]. 工程地质学报, 2019, 28(4): 725-733. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202004006.htm ZHANG Bing-shu, GU Kai, LI Jin-wen, et al. Study on crushing process of calcareous sand and its microscopic mechanism[J]. Journal of Engineering Geology, 2019, 28(4): 725-733. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ202004006.htm
[9]	OHNO S, OCHIA H, YASUFUKU N. Estimation of pile settlement in calcareous sans[C]//The Second International Conference on Engineering for Calcareous Sediments, 1999, Manama: 1-6.
[10]	虞海珍, 汪稔. 钙质砂动强度试验研究[J]. 岩土力学, 1999, 20(4): 6-11. doi: 10.3969/j.issn.1000-7598.1999.04.002 YU Hai-zhen, WANG Ren. The cyclic strength test research on calcareous sand[J]. Rock and Soil Mechanics, 1999, 20(4): 6-11. (in Chinese) doi: 10.3969/j.issn.1000-7598.1999.04.002
[11]	朱长歧, 周斌, 刘海峰. 天然胶结钙质土强度及微观结构研究[J]. 岩土力学, 2014, 35(6): 1655-1663. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201406022.htm ZHU Chang-qi, ZHOU Bin, LIU Hai-feng. Investigation on strength and microstracture of naturally cemented calcareous soil[J]. Rock and Soil Mechanics, 2014, 35(6): 1655-1663. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201406022.htm
[12]	张家铭. 钙质砂基本力学性质及颗粒破碎影响研究[D]. 武汉: 中国科学院研究生院(武汉岩土力学研究所), 2004. ZHANG Jia-ming. Study on the Fundamental Mechanical Characteristics of Calcareous Sand and the Influence of Particle Breakage[D]. Wuhan: Institute of Rock & Soil Mechanics, Chinese Academy of Sciences, 2004. (in Chinese)
[13]	吴京平, 褚瑶. 颗粒破碎对钙质砂变形及强度特性的影响[J]. 岩土工程学报, 1997, 19(5): 51-57. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC705.007.htm WU Jing-ping, CHU Yao. Influence of particle breakage on deformation and strength properties of calcareous sands[J]. Chinese Journal of Geotechnical Engineering, 1997, 19(5): 51-57. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC705.007.htm
[14]	王新志. 南沙群岛珊瑚礁工程地质特性及大型工程建设可行性研究[D]. 武汉: 中国科学院研究生院(武汉岩土力学研究所), 2008. WANG Xin-zhi. Study on Engineering Geological Properties of Coral Reefs and Feasibility of Large Project Construction on Nanshaislands[D]. Wuhan: Institute of Rock & Soil Mechanics, Chinese Academy of Sciences, 2008. (in Chinese)
[15]	孙宗勋. 南沙群岛珊瑚砂工程性质研究[J]. 热带海洋学报, 2000, 19(2): 1-8. doi: 10.3969/j.issn.1009-5470.2000.02.001 SUN Zong-xun. Engineering properties of coral sands in Nanshaislands[J]. Journal of Tropical Oceanography, 2000, 19(2): 1-8. (in Chinese) doi: 10.3969/j.issn.1009-5470.2000.02.001
[16]	WANG X Z, JIAO Y Y, WANG R, et al. Engineering characteristics of the calcareous sand in Nansha Islands, South China Sea[J]. Engineering Geology, 2011, 120(1/2/3/4): 40-47.
[17]	MURFF J D. Pile Capacity in calcareous sands: state if the art[J]. Journal of Waterway, Port, Coastal and Ocean Engineering, 1987, 113(5): 490-507.
[18]	MCCLELLAND B. Calcareous sediments: an engineering enigma[C]//Proc Int Cong On Calcareous Sediments, 1988, Perth: 777-784.
[19]	樊恒辉, 高建恩, 吴普特. 土壤固化剂研究现状与展望[J]. 西北农林科技大学学报(自然科学版), 2006, 34(2): 141-146. https://www.cnki.com.cn/Article/CJFDTOTAL-XBNY200602029.htm FAN Heng-hui, GAO Jian-en, WU Pu-te. Prospect of researches on soil stabilizer[J]. Journal of Northwest A & F University (Natural Science Edition), 2006, 34(2): 141-146. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XBNY200602029.htm
[20]	宋云, 李培中, 郝润琴. 中国土壤固化/稳定化技术应用现状及建议[J]. 环境保护, 2015, 43(15): 28-33. https://www.cnki.com.cn/Article/CJFDTOTAL-HJBU201515013.htm SONG Yun, LI Pei-zhong, HAO Run-qin. Analysis on the application status and advice of solidification/stabilization in China[J]. Environmental Protection,,2015, 43(15): 28-33. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HJBU201515013.htm
[21]	王灵秀, 王雅明, 郝庆军, 等. 水泥行业熟料生产CO2排放调查研究[J]. 中国建材科技, 2010(增刊2): 96-99. https://www.cnki.com.cn/Article/CJFDTOTAL-JCKJ2010S2022.htm WANG Ling-xiu, WANG Ya-ming, HAO Qing-jun, et al. Survey and research on CO2 emissions in clinker production of cement industry[J]. China Building Materials Science & Technology, 2010(S2): 96-99. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JCKJ2010S2022.htm
[22]	钱春香, 王安辉, 王欣. 微生物灌浆加固土体研究进展[J]. 岩土力学, 2015, 36(6): 1537-1548. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201506003.htm QIAN Chun-xiang, WANG An-hui, WANG Xin. Advances of soil improvement with bio-grouting[J]. Rock and Soil Mechanics, 2015, 36(6): 1537-1548. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201506003.htm
[23]	荣辉, 钱春香. 微生物水泥基材料特征[J]. 科学通报, 2012, 57(9): 770-775. https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201209013.htm RONG Hui, QIAN Chun-xiang. Characterization of microbe cementitious materials[J]. Chinese Science Bulletin, 2012, 57(9): 770-775. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201209013.htm
[24]	刘璐, 沈扬, 刘汉龙, 等. 微生物胶结在防治堤坝破坏中的应用研究[J]. 岩土力学, 2016, 37(12): 3410-3416. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201612009.htm LIU Lu, SHEN Yang, LIU Han-long, et al. Application of bio-cement in erosion control of levees[J]. Rock and Soil Mechanics, 2016, 37(12): 3410-3416. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201612009.htm
[25]	崔明娟, 郑俊杰, 赖汉江. 颗粒粒径对微生物固化砂土强度影响的试验研究[J]. 岩土力学, 2016, 37(增刊2): 397-402. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2016S2051.htm CUI Ming-juan, ZHENG Jun-jie, LAI Han-jiang. Experimental study of effect of particle size on strength of bio-cemented sand[J]. Rock and Soil Mechanics, 2016, 37(S2): 397-402. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2016S2051.htm
[26]	邵光辉, 冯建挺, 赵志峰, 等. 微生物砂浆防护粉土坡面的强度与抗侵蚀性影响因素分析[J]. 农业工程学报, 2017, 33(11): 141-147. https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201711017.htm SHAO Guang-hui, FENG Jian-ting, ZHAO Zhi-feng, et al. Influence factor analysis related to strength and anti-erosion stability of silt slope with microbial mortar protective covering[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(11): 141-147. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-NYGU201711017.htm
[27]	DEJONG J T, SOGA K, KAVAZANJIAN E, et al. Biogeochemical processes and geotechnical applications: Progress, opportunities and challenges[J]. Géotechnique, 2013, 63(4): 287-301.
[28]	MARTINEZ B C, DEJONG J T, GINN J, et al. Experimental optimization of microbial-induced carbonate precipitation for soil improvement[J]. Journal of Geotechnical & Geo- environmental Engineering, 2013, 139(4): 587-598.
[29]	QABANY A A, SOGA K, SANTAMARINA C. Factors affecting efficiency of microbially induced calcite precipitation[J]. Journal of Geotechnical &Geoenvironmental Engineering, 2012, 138(8): 992-1001.
[30]	何稼, 楚剑, 刘汉龙, 等. 微生物岩土技术的研究进展[J]. 岩土工程学报, 2016, 38(4): 643-653. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201604010.htm HE Jia, CHU Jian, LIU Han-long, et al. Research advances in biogeotechnologies[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(4): 643-653. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201604010.htm
[31]	李明东, LIN Li, 张振东, 等. 微生物矿化碳酸钙改良土体的进展、展望与工程应用技术设计[J]. 土木工程学报, 2016, 49(10): 80-87. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201610013.htm LI Ming-dong, LIN Li, ZHANG Zhen-dong, et al. Review, outlook and application technology design on soil improvement by microbial induced calcium carbonate precipitation[J]. China Civil Engineering Journal, 2016, 49(10): 80-87. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201610013.htm
[32]	FERRIS F G, HALLBECK L, KENNEDY C B, et al. Geochemistry of acidic Rio Tinto headwaters and role of bacteria in solid phase metal partitioning[J]. Chemical Geology, 2004, 212(3/4): 291-300.
[33]	DEJONG J T, MORTENSEN B M, MARTINEZ B C, et al. Bio-mediated soil improvement[J]. Ecological Engineering, 2010, 36(2): 197-210.
[34]	DEJONG J T, FRITZGES M B, NÜSSLEIN K. Microbially induced cementation to control sand response to undrained shear[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(11): 1381-1392.
[35]	PAASSEN V L A. Bio-mediated ground improvement: from laboratory experiment to pilot applications[C]//Geo-Frontiers Congress 2011, 2011, Dallas: 4099-4108.
[36]	CHU J, STABNIKOV V, IVANOV V. Microbially induced calcium carbonate precipitation on surface or in the bulk of soil[J]. Geomicrobiology Journal, 2012, 29(6): 544-549.
[37]	GOMEZ M G, MARTINEZ B C, DEJONG J T, et al. Field-scale bio-cementation tests to improve sands[J]. Proceedings of the Institution of Civil Engineers Ground Improvement, 2014, 168(3): 1-11.
[38]	JIANG N J, SOGA K. The applicability of microbially induced calcite precipitation (MICP) for internal erosion control in gravel-sand mixtures[J]. Géotechnique, 2016, 67(1): 42-55.
[39]	JIANG N J, SOGA K, KUO M. Microbially induced carbonate precipitation for seepage-induced internal erosion control in sand-clay mixtures[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 143(3): 04016100.
[40]	LIU L, LIU H, XIAO Y, et al. Biocementation of calcareous sand using soluble calcium derived from calcareous sand[J]. Bulletin of Engineering Geology and the Environment, 2018, 77: 1781-1791.
[41]	方祥位, 申春妮, 楚剑, 等. 微生物沉积碳酸钙固化珊瑚砂的试验研究[J]. 岩土力学, 2015, 36(10): 2773-2779. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510005.htm FANG Xiang-wei, SHEN Chun-ni, CHU Jian, et al. An experimental study of coral sand enhanced through microbially-induced precipitation of calcium carbonate[J]. Rock and Soil Mechanics, 2015, 36(10): 2773-2779. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510005.htm
[42]	李捷, 方祥位, 申春妮, 等. 含水率对珊瑚砂微生物固化体力学特性影响研究[J]. 工业建筑, 2016, 46(12): 93-97. https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ201612019.htm LI Jie, FANG Xiang-wei, SHEN Chun-ni, et al. Influence of moisture content on mechanical properties of biocemented coral sand columns[J]. Industrial Construction, 2016, 46(12): 93-97. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GYJZ201612019.htm
[43]	方祥位, 李晶鑫, 李捷, 等. 珊瑚砂微生物固化体单轴损伤本构模型[J]. 地下空间与工程学报, 2018, 14(5): 93-98. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201805012.htm FANG Xiang-wei, LI Jing-xin, LI Jie, et al. Damage constitutive model of biocemented coral sand columnsunder unconfined compression[J]. Chinese Journal of Underground Space and Engineering, 2018, 14(5): 93-98. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201805012.htm
[44]	方祥位, 李晶鑫, 李捷, 等. 珊瑚砂微生物固化体三轴压缩试验及损伤本构模型研究[J]. 岩土力学, 2018, 39(增刊1): 10-17. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S1002.htm FANG Xiang-wei, LI Jing-xin, LI Jie, et al. Study of triaxial compression test and damage constitutive model of biocemented coral sand columns[J]. Rock and Soil Mechanics, 2018, 39(S1): 10-17. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S1002.htm
[45]	刘汉龙, 肖鹏, 肖杨, 等. MICP胶结钙质砂动力特性试验研究[J]. 岩土工程学报, 2018, 40(1): 38-45. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201801003.htm LIU Han-long, XIAO Peng, XIAO Yang, et al. Dynamic behaviors of MICP-treated calcareous sand in cyclic tests[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(1): 38-45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201801003.htm
[46]	马瑞男, 郭红仙, 程晓辉, 等. 微生物拌和加固钙质砂渗透特性试验研究[J]. 岩土力学, 2018, 39(增刊2): 217-223. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S2031.htm MA Rui-nan, GUO Hong-xian, CHENG Xiao-hui, et al. Permeability experiment study of calcareous sand treated by microbially induced carbonate precipitation using mixing methods[J]. Rock and Soil Mechanics, 2018, 39(S2): 217-223. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2018S2031.htm
[47]	CHENG L, SHAHIN M, CORDRUWISCH R, et al. Soil stabilisation by microbial-induced calcite precipitation (MICP): investigation into some physical and environmental aspects[C]//7th International Congress on Environmental Geotechnics, 2014, Melbourne.
[48]	欧益希, 方祥位, 张楠, 等. 溶液盐度对微生物固化珊瑚砂的影响[J]. 后勤工程学院学报, 2016, 32(1): 78-82. https://www.cnki.com.cn/Article/CJFDTOTAL-HQGC201601015.htm OU Yi-xi, FANG Xiang-wei, ZHANG Nan, et al. Influence of solution salinity on microbial biocementation of coral sand[J]. Journal of Logistical Engineering University, 2016, 32(1): 78-82. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-HQGC201601015.htm
[49]	彭劼, 田艳梅, 杨建贵. 海水环境下MICP加固珊瑚砂试验[J]. 水利水电科技进展, 2019, 39(1): 62-66. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201901011.htm PENG Jie, TIAN Yan-mei, YANG Jian-gui. Experiments of coral sand reinforcement using MICP in seawater environment[J]. Advances in Science and Technology of Water Resources, 2019, 39(1): 62-66. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201901011.htm
[50]	吴洋, 练继建, 闫玥, 等. 巴氏芽孢八叠球菌及相关微生物的生物矿化的分子机理与应用[J]. 中国生物工程杂志, 2017, 37(8): 96-103. https://www.cnki.com.cn/Article/CJFDTOTAL-SWGJ201708014.htm WU Yang, LIAN Ji-jian, YAN Yue, et al. Mechanism and applications of bio-mineralization induced by sporosarcinapasteurii and related mic- roorganisms[J]. China Biotechnology, 2017, 37(8): 96-103. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWGJ201708014.htm
[51]	谢约翰, 唐朝生, 尹黎阳, 等. 纤维加筋微生物固化砂土的力学特性[J]. 岩土工程学报, 2019, 41(4): 675-682. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201904014.htm XIE Yue-han, TANG Chao-sheng, YIN Li-yang, et al. Mechanical behavior of microbial-induced calcite precipitation (MICP)-treated soil with fiber reinforcement[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(4): 675-682. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201904014.htm
[52]	尹黎阳. 基于微生物诱导碳酸钙沉淀技术(MICP)固化钙质砂的力学性能研究[D]. 南京: 南京大学, 2019. YIN Li-yang. Experimental Study on the Mechanical Behavior of Fiber Reinforcement Calcareous Sand Cemented by MICP Technology[D]. Nanjing: Nanjing University, 2019. (in Chinese)
[53]	CHOI S G, PARK S S, WU S, et al. Methods for calcium carbonate content measurement of biocemented soils[J]. Journal of Materials in Civil Engineering, 2017, 29(11): 06017015.
[54]	赵茜. 微生物诱导碳酸钙沉淀(MICP)固化土壤实验研究[D]. 北京: 中国地质大学, 2014. ZHAO Qian. Experimental Study on Soil Improvement Using Microbial Induced Calcite Precipitation (MICP)[D]. Beijing: China University of Geosciences, 2014. (in Chinese)
[55]	欧益希, 方祥位, 申春妮, 等. 颗粒粒径对微生物固化珊瑚砂的影响[J]. 水利与建筑工程学报, 2016, 14(2): 35-39. https://www.cnki.com.cn/Article/CJFDTOTAL-FSJS201602007.htm OU Yi-xi, FANG Xiang-wei, SHEN Chun-ni, et al. Influence of particle sizes of coral sand on bio-cementation[J]. Journal of Water Resources and Architectural Engineering, 2016, 14(2): 35-39. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-FSJS201602007.htm
[56]	BARIYANGA J. Concentration Gradient[M]//BINDER M D, HIROKAWA N, WINDHORST U, ed. Encyclopedia of Neuroscience, Berlin: Springer, 2009.
[57]	CHENG L, SHAHIN M A, CORD-RUWISCH R. Bio-cementation of sandy soil using microbially induced carbonate precipitation for marine environments[J]. Géotechnique, 2014, 64(12): 1010-1013.
[58]	尹黎阳, 唐朝生, 谢约翰, 等. 微生物矿化作用改善岩土材料性能的影响因素[J]. 岩土力学, 2019, 40(7): 2525-2546. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201907007.htm YIN Li-yang, TANG Chao-sheng, XIE Yue-han, et al. Factors affecting improvement in engineering properties of geomaterials by microbial-induced calcite precipitation[J]. Rock and Soil Mechanics, 2019, 40(7): 2525-2546. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201907007.htm
[59]	STOCKS-FISCHER S, GALINAT J K, BANG S. Microbiological precipitation of CaCO₃[J]. Soil Biology & Biochemistry, 1999, 31(11): 1563-1571.
[60]	黄琰, 罗学刚, 杜菲. 微生物诱导方解石沉积加固的影响因素[J]. 西南科技大学学报, 2009, 24(3): 87-93. https://www.cnki.com.cn/Article/CJFDTOTAL-XNGX200903017.htm HUANG Yan, LUO Xue-gang, DU Fei. Studies on the effect factor of microbiologically induced calcite precipitation[J]. Journal of Southwest University of Science and Technology, 2009, 24(3): 87-93. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XNGX200903017.htm
[61]	FERRIS F G, FYFE W S, BEVERIDGE T J. Bacteria as nucleation sites for authigenic minerals in a metal- contaminated lake sediment[J]. Chemical Geology, 1987, 63(3): 225-232.
[62]	FENG K, MONTOYA B M. Influence of confinement and cementation level on the behavior of microbial-induced calcite precipitated sands under monotonic drained loading[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2016, 142(1): 04015057.
[63]	CHOI S G, JIAN C, BROWN R C, et al. Sustainable biocement production via microbially-induced calcium carbonate precipitation: use of limestone and acetic acid derived from pyrolysis of lignocellulosic biomass[J]. Acs Sustainable Chemistry, 2017, 5(6): 5183-5190.

施引文献(2)

期刊类型引用(2)

1.	刘姝，李文杰，李莉佳，罗永江，陶瑞，李晓璇，杨亚会. 天然气水合物开采方法研究现状及展望. 钻探工程. 2024(05): 12-23 . 百度学术
2.	于倩男，唐慧敏，李承龙，梁爽，梁家修，陈志静，张琨. 天然气水合物注热分解渗流特征及数值模拟. 东北石油大学学报. 2023(06): 38-54+127-128 . 百度学术