水合物开采模拟超重力试验装置的研发及应用

王路君; 王鹏; 朱斌; 王心博; 杨颂清; 陈云敏

doi:10.11779/CJGE20221300

水合物开采模拟超重力试验装置的研发及应用 English Version

王路君^{1, 2,},
王鹏^{1, 2, ,},
朱斌^{1, 2},
王心博²,
杨颂清²,
陈云敏^{1, 2}

1.
浙江大学超重力研究中心, 浙江杭州 310058
2.
浙江大学软弱土与环境土工教育部重点实验室, 浙江杭州 310058

基金项目:

国家自然科学基金项目 52127815

国家自然科学基金项目 51988101

国家自然科学基金项目 52078458

详细信息

作者简介:
王路君(1985— )，男，博士，副教授，从事水合物沉积物力学特性和土体多相多场耦合等方面的研究工作。E-mail: lujunwang@zju.edu.cn

通讯作者:
王鹏, E-mail: 12012041@zju.edu.cn

中图分类号: TU415
计量
- 文章访问数: 349
- HTML全文浏览量: 34
- PDF下载量: 84
出版历程
- 收稿日期: 2022-10-20
- 网络出版日期: 2023-03-13
- 刊出日期: 2024-01-31

Development and application of in-flight centrifuge apparatus for modelling hydrate exploitation

WANG Lujun^{1, 2,},
WANG Peng^{1, 2, ,},
ZHU Bin^{1, 2},
WANG Xinbo²,
YANG Songqing²,
CHEN Yunmin^{1, 2}

1.
Center for Hypergravity Experimental and Interdisciplinary Research, Zhejiang University, Hangzhou 310058, China
2.
Key Laboratory of Soft Soils and Geoenvironmental Engineering of the Ministry of Education, Zhejiang University, Hangzhou 310058, China

摘要

摘要: 天然气水合物是突破世界能源瓶颈的重要潜在能源，广泛分布于深海沉积物和陆域永冻土。水合物赋存在高压低温环境下，储层厚度达数十至百米，大尺度储层重力场显著影响其开采时的液气渗流、储层变形和传热过程，进而影响开采产气效率且可能诱发工程灾变。研发了一套能够在超重力场下模拟大尺度水合物储层长历时开采演变过程的离心模型试验装置，主要由高压釜系统、水合物制备模块、水压力控制模块、水合物开采模块、监测与响应模块组成。该装置能够在100g超重力下伺服调控高压釜内气液压力和温度环境，在模型尺度上有效再现水合物储层应力、温度和荷载条件。在系统介绍装置研发思路和技术特色的基础上，通过100g下水合物储层模型降压分解试验探析了原位等效重力场下储层不同深度处孔压演化规律及开采产气特性，再现了Nankai海槽原位试采产气速率峰值模式，为揭示水合物储层开采相变多相多场耦合机制和灾变防控提供了新的研究手段。
- 水合物 /
- 离心模型试验 /
- 降压开采 /
- 高压釜 /
- 产气速率 /
- 伺服控压
Abstract: The natural gas hydrate, widely distributed in deep sea sediment and terrestrial permafrost, is considered as an important potential energy to solve the increasingly serious energy crisis. Hydrate is formed under high pressure and low temperature conditions, where the reservoir is tens of meters in thickness. The gravity field of large-scale reservoir significantly affects the process of fluid seepage, skeleton deformation and heat transfer during exploitation, changing the gas production rate and further inducing possible engineering disasters. An in-flight centrifuge apparatus, which is composed of high-pressure vessel, hydrate synthesis module, water pressure control module, hydrate exploitation module and monitoring module, for modelling hydrate exploitation is developed. The apparatus can servo control the pressure and temperature conditions in the high-pressure vessel under 100g hyper gravity, and effectively reproduce the stress, temperature and loading conditions of hydrate reservoir in the model scale. A hydrate exploitation modelling test via depressurization at 100g is then carried out. The evolution of pore pressure at different depths and the gas production characteristics under the field scale stress gradient are analyzed, and the gas production characteristics of field production tests in Nankai trough are reproduced. The developed apparatus provides a novel experimental approach to investigate the multiphase coupling mechanics and disasters prevention during hydrate exploitation.
- hydrate /
- centrifuge modelling /
- depressurized exploitation /
- high-pressure vessel /
- gas production rate /
- pressure servo control

HTML全文

0. 引言

离散单元法（discrete element method, DEM）^[1]是一种基于牛顿第二定律与颗粒间接触本构的数值仿真方法，能记录试样内部颗粒的位置与运动情况、捕捉颗粒间接触特征等微观信息，用于揭示砂土、砾土等颗粒状岩土工程材料力学特性的微观机理。

线性弹簧模型和Hertz-Mindlin模型是DEM常用的两类接触模型^[2]，均假设颗粒为弹性材料，模型参数包括颗粒的弹性参数。以Hertz-Mindlin接触模型为例，两个等直径球体的在法向接触力F_N的作用下，颗粒间接触面为圆形，其半径a可以由Hertz理论计算得到：

$a = \sqrt[\begin{subarray}{l} 3 \\ \end{subarray} ]{{\frac{{3(1 - {\nu _\text{P}}^2){F_\text{N}}R}}{{4{E_\text{p}}}}}}。$

(1)

式中：R为球体半径；E_p和ν_p分别为颗粒的杨氏模量和泊松比。

颗粒间接触的法向刚度K_N及切向刚度K_T计算公式如下：

$K_{\mathrm{N}}=\frac{2 E_{\mathrm{g}} a}{3\left(1-v_{\mathrm{P}}^2\right)} ,$

(2)

$K_{\mathrm{T}}=\frac{2 E_{\mathrm{P}} a}{\left(2-v_{\mathrm{P}}\right)\left(1+v_{\mathrm{P}}\right)}\left(1-\frac{F_{\mathrm{T}}}{\mu F_{\mathrm{N}}}\right)^{1 / 3} 。$

(3)

式中：µ为颗粒间摩擦系数；F_T为颗粒间切向接触力。

在外力作用下，颗粒间接触力增量与颗粒间位移增量通过接触刚度（式（2）、（3））关联。显然，颗粒的弹性参数是主要控制参数。

然而，学者们主要通过宏观力学试验的数值仿真结果反演得到接触模型参数，即选取多种参数组合模拟土单元体试验（如三轴试验、直剪试验等），与试验结果比对选取最优的模型参数。这类方法经验性较强，缺乏统一标准^[2]；此外，多个参数组合可能产生相近宏观模拟结果而微观特征存在较大差异，甚至某些参数组合还会导致模拟结果在微观尺度上失真^[3]。因此，有学者尝试开展颗粒尺度试验以标定接触模型参数，以期在宏微观尺度均得到合理的离散元仿真结果^[3-5]。

本文以玻璃珠为试验对象，尝试采用纳米压痕试验获得玻璃珠的弹性模量，并与玻璃块的弹性波速测试结果进行比对，验证纳米压痕试验可用于测量颗粒弹性模量。该方法可以用于测量真实砂土和砾石土颗粒的弹性模量。

1. 纳米压痕试验

1.1 测试材料

测试材料为直径2 mm的玻璃珠（GB-2-A）。利用X射线能谱仪测试得到玻璃珠能谱图如图 1所示，表明玻璃珠的化学元素组成主要为Na₂SiO₂∙2SiO₂。

图 1 玻璃珠X射线能谱图

Figure 1. X-ray energy spectra of glass beads

下载: 全尺寸图片幻灯片

1.2 纳米压痕仪试验原理

纳米压痕试验广泛用于测试金属、半导体、涂层、薄膜、生物材料等的杨氏模量和硬度等力学性质^[6-7]。本研究采用Agilent Technologies生产的G200纳米压痕仪（Nano Indenter G200），图 2是其内部结构示意图，将三角锥形状的金刚石压头压入所测材料表面，利用高精度传感器测量压头在加-卸载过程中的压入深度和压力。图 3是纳米压痕试验中典型的荷载-压入深度曲线，可用于计算得到所测材料的弹性模量。

图 2 G200纳米压痕仪内部结构示意图

Figure 2. Schematics of internal structures of Nano Indenter G200

下载: 全尺寸图片幻灯片

图 3 加-卸载过程的荷载-压入深度曲线

Figure 3. Load-depth curves during loading and unloading

下载: 全尺寸图片幻灯片

1.3 弹性模量计算方法

Oliver-Pharr方法^[8]是广为应用的纳米压痕试验分析方法，定义压头与材料间的接触刚度S^*为荷载P与压入深度h关系曲线的初始卸载斜率：

${S^ * }{\text{ = }}\frac{{{\text{d}}P}}{{{\text{d}}h}}|{h_{\max }} 。$

(4)

式中：h_max为加载过程的最大压入深度。

本研究采用Berkovich压头，利用下式计算压头的有效接触深度h_c以及等效弹性模量E_r：

${h_\text{c}} = h - 0.75\frac{P}{{{S^ * }}} \text{，}$

(5)

${E_\text{r}} = \frac{{\sqrt {\text{π }} {S^ * }}}{{2{\beta _\text{s}}\sqrt {{A_\text{c}}} }} 。$

(6)

式中：β_s为形状校正因子，对于Berkovich压头取β_s = 1.05；A_c为压头在样品表面的接触面积，对于Berkovich压头取A_c = 24.56h_c²。

进一步，可利用下式计算得到所测材料的杨氏模量E_p：

${E_\text{p}} = (1 - {\nu _\text{P}}^2){\left[ {\frac{1}{{{E_\text{r}}}} - \frac{{1 - \nu _\text{i}^2}}{{{E_\text{i}}}}} \right]^{ - 1}} 。$

(7)

式中：E_i， $\nu _\text{i}^{}$ 为压头的弹性模量与泊松比，Berkovich压头材料为金刚石，取E_i = 1140 GPa， $\nu _\text{i}^{}$ = 0.07； ${\nu _\text{P}}$ 为材料泊松比。材料的硬度H则可根据峰值荷载P_max以及接触面积A_c计算：

$H = \frac{{{P_{\max }}}}{{{A_\text{c}}}} 。$

(8)

测试时将玻璃珠固定于载物台上，选择玻璃珠顶部区域，尽可能保证压头垂直压入玻璃珠样品。由于测试结果具有一定离散性，选取不同区域测试10组以上。图 4是纳米压痕受测颗粒与压入点的照片。

图 4 纳米压痕受测颗粒表面与压入点

Figure 4. Particle surface and indentation point

下载: 全尺寸图片幻灯片

玻璃珠的多组荷载-压入深度曲线如图 5所示。荷载-压入深度曲线的斜率随压入深度增加而增大，表明法向接触刚度随压入深度与接触力增加而增大，这与经典的Hertz接触模型中法向接触刚度与接触力的关系相似^[9-10]。本文取压入深度为1500 nm时的荷载-压入深度曲线用于计算杨氏模量。图 6是颗粒杨氏模量E_p和硬度H的频率统计图。E_p主要分布于50~60 GPa，平均值为59.7 GPa。

图 5 玻璃珠荷载-压入深度曲线

Figure 5. Load-indentation depth curves of tested glass beads

下载: 全尺寸图片幻灯片

图 6 玻璃珠杨氏模量与硬度的频率分布图

Figure 6. Frequency distribution of Young's modulus and hardness

下载: 全尺寸图片幻灯片

2. 测试结果验证

为验证纳米压痕试验的准确性，采用Pundit PL-200超声波仪测试玻璃试块的P波波速V_p和S波波速V_s，用下式计算剪切模量G_P和侧限模量M_P：

$G_{\mathrm{P}}=\rho V_{\mathrm{s}}^2,$

(9)

$M_{\mathrm{P}}=\rho V_{\mathrm{p}}^2 \quad 。$

(10)

式中：ρ为玻璃块的密度（2533 kg/m³）。

玻璃的泊松比 ${\nu _P}$ 由下式计算得到：

$v_{\mathrm{P}}=\frac{1}{2}\left(1-\frac{1}{M_{\mathrm{P}} / G_{\mathrm{P}}-1}\right) 。$

(11)

玻璃的杨氏模量E_p由下式得到：

$E_{\mathrm{P}}=2\left(1+v_{\mathrm{P}}\right) G_{\mathrm{P}} 。$

(12)

本文利用超声波仪测试玻璃试块的V_p和V_s。在3个不同方向采用160 kHz的P波信号与50 kHz的S波信号作为入射波，接收到的波形如图 7所示。P波和S波分别通过约9.8 μs（54 mm）与16 μs（49 mm）到达接收探头，得到V_p = 5506 m/s，V_s = 3063 m/s。进一步计算得到 ${\nu _P}$ = 0.276，E_p = 60.7 GPa。可见，波速测试得到的E_p与纳米压痕试验得到的结果相近。

图 7 玻璃块P波与S波信号

Figure 7. Signals of P and S waves of glass block

下载: 全尺寸图片幻灯片

3. 离散元数值仿真

本研究利用纳米压痕试验所得E_p用于离散元数值仿真，模拟土体在小应变（γ < 10^-5）扭剪条件下的力学特性。颗粒之间的接触采用Hertz-Mindlin接触模型，颗粒与墙体间接触采用线性接触模型，模型参数见表 1，其中E_p取自纳米压痕试验。

表 1 离散元数值仿真接触模型参数

Table 1. Parameters of contact model in DEM simulation

E_p/ GPa	${\nu _P}$	墙体刚度/ (kN·m^-1)	β	μ	颗粒-墙体摩擦系数(两端/侧墙)
60	0.276	5×10⁶	0.1	0.25	10⁵ / 0.001

下载: 导出CSV

| 显示表格

试样采用半径扩大法生成颗粒：首先随机生成半径为目标半径一半的颗粒，使颗粒间不重叠，之后再将颗粒半径扩大1倍。试样生成后，对试样进行30 kPa等向预固结，随后对墙体进行伺服控制达到指定固结有效应力（σ'为50，100，200，400 kPa），并通过设定不同的颗粒间摩擦系数µ得到不同初始孔隙比的试样。在剪切阶段，固定试样底部墙体，在试样顶部墙体施加扭矩实现扭剪运动，扭剪速度为0.02 rad/s。本研究仅对颗粒间接触设置黏滞阻尼β = 0.1。扭剪过程中剪应力τ与剪应变γ的计算公式为

$\tau=R_{\mathrm{s}}\left(k^{\prime} \frac{T}{J}\right),$

(13)

$\gamma=R_{\mathrm{s}}\left(k^{\prime} \frac{\theta}{l}\right)。$

(14)

式中：J为试样截面的极惯性矩；T为扭矩；R_s为试样半径；θ为试样扭剪角度；l为试样高度；k'取0.8。

剪切模量G可通过加载阶段的τ-γ曲线得到：

$G=\frac{\tau}{\gamma} 。$

(15)

本研究利用共振柱所测的G_max数据验证离散元模拟结果（试验结果来自文献[11]）。试验材料为直径2 mm的玻璃珠（即GB-2-A）；试样直径为38 mm、高度为81 mm，通过干沉积法制备得到。图 8对比了离散元数值仿真与共振柱试验得到的小应变剪切模量G_max。G_max随孔隙比e增加而减小，模拟结果的G_max-e关系与试验结果基本重合。

图 8 各有效应力条件下小应变剪切模量与孔隙比关系

Figure 8. G_max-e relationships under various effective confining pressures of test and simulated results

下载: 全尺寸图片幻灯片

4. 讨论

文献主要采用颗粒间压缩试验用于标定颗粒的弹性参数。例如，Reddy等^[3]利用真实砂土颗粒间的压缩试验所得的力-变形曲线反演出颗粒杨氏模量，用于离散元数值仿真。然而这类方法在实际应用中存在一些局限：首先，颗粒间压缩试验并非常规试验，且试验结果的解读存在一定主观性；其次，颗粒间的力-变形曲线受到颗粒接触点局部形态的影响，而接触点的颗粒局部形态量化表征亦非易事；此外，这类方法需要基于特定的颗粒接触模型，可能无法全面考虑颗粒形状、颗粒表面特征对力-变形曲线的影响。本研究采用的纳米压痕试验是较为常见的成熟技术手段。对于常见的砂土、砾土颗粒，纳米尺度的压入深度所得结果受颗粒形状影响有限，亦能反映颗粒表面特性的影响。本研究利用形状规则的球形材料验证了利用纳米压痕试验标定颗粒弹性参数的可能性，具有进一步用于天然砂土、砾土颗粒的潜力，并有待进一步研究。

5. 结论

本研究采用纳米压痕试验测试了2 mm玻璃珠的杨氏模量，与通过弹性波速测试得到的杨氏模量基本相同，证明纳米压痕试验可以用于测量球形颗粒的杨氏模量。进一步地，利用纳米压痕测得的杨氏模量用于离散元数值仿真，成功模拟了扭剪模式下的玻璃珠小应变剪切特性。该方法相比于文献中基于颗粒间压缩试验反演杨氏模量的方法具有多方面优势，具有进一步应用于天然砂土、砾土颗粒的潜力。

图 1 装置原理图

Figure 1. Schematic diagram of apparatus

下载: 全尺寸图片幻灯片

图 2 装置系统组成

Figure 2. System composition of apparatus

下载: 全尺寸图片幻灯片

图 3 高压釜系统整体结构

Figure 3. High-pressure vessel system

下载: 全尺寸图片幻灯片

图 4 回压泵-回压阀伺服调控开采井压力

Figure 4. Pump-valve combined control of pressure in vertical well

下载: 全尺寸图片幻灯片

图 5 测点空间分布和模型剖示图

Figure 5. Distribution of measuring points and model diagram

下载: 全尺寸图片幻灯片

图 6 ZJU400离心机及装置搭载于吊篮

Figure 6. ZJU400 centrifuge at Zhejiang University and apparatus loaded on centrifuge

下载: 全尺寸图片幻灯片

图 7 水合物储层模型制备流程

Figure 7. Process of hydrate reservoir model preparation

下载: 全尺寸图片幻灯片

图 8 水合物制备过程中孔压时程曲线

Figure 8. Evolution of pressure during hydrate synthesis

下载: 全尺寸图片幻灯片

图 9 水合物制备过程中温度时程曲线

Figure 9. Evolution of temperature during hydrate synthesis

下载: 全尺寸图片幻灯片

图 10 水合物制备过程中T₃-P₃测点的温压路径

Figure 10. Temperature-pressure response throughout hydrate synthesis process of T₃-P₃

下载: 全尺寸图片幻灯片

图 11 100g下水合物降压分解时模型的温压路径

Figure 11. Temperature-pressure response of model during hydrate dissociation by depressurization at 100g

下载: 全尺寸图片幻灯片

图 12 100g下水合物降压分解时模型孔压变化

Figure 12. Evolution of excess pore pressure at 100g

下载: 全尺寸图片幻灯片

图 13 100g和1g下水合物分解归一化产气量

Figure 13. Gas production rates during hydrate dissociation by depressurization at 100g

下载: 全尺寸图片幻灯片

图 14 100g下水合物降压分解产气过程

Figure 14. Gas production during hydrate dissociation by depressurization at 100g

下载: 全尺寸图片幻灯片

参考文献(32)

[1]	SLOAN E D Jr. Fundamental principles and applications of natural gas hydrates[J]. Nature, 2003, 426(6964): 353-359. doi: 10.1038/nature02135
[2]	BOSWELL R, COLLETT T S. Current perspectives on gas hydrate resources[J]. Energy & Environmental Science, 2011, 4(4): 1206-1215.
[3]	RICE W. Hydrogen production from methane hydrate with sequestering of carbon dioxide[J]. International Journal of Hydrogen Energy, 2006, 31(14): 1955-1963. doi: 10.1016/j.ijhydene.2006.01.017
[4]	LI X S, YANG B, ZHANG Y, et al. Experimental investigation into gas production from methane hydrate in sediment by depressurization in a novel pilot-scale hydrate simulator[J]. Applied Energy, 2012, 93: 722-732. doi: 10.1016/j.apenergy.2012.01.009
[5]	TANG L G, XIAO R, HUANG C, et al. Experimental investigation of production behavior of gas hydrate under thermal stimulation in unconsolidated sediment[J]. Energy & Fuels, 2005, 19(6): 2402-2407.
[6]	ZHANG L X, YANG L, WANG J Q, et al. Enhanced CH₄ recovery and CO₂ storage via thermal stimulation in the CH₄/CO₂ replacement of methane hydrate[J]. Chemical Engineering Journal, 2017, 308: 40-49. doi: 10.1016/j.cej.2016.09.047
[7]	ZHOU J Z, ZHOU Y A, YANG Z J, et al. Dissociation-induced deformation of hydrate-bearing silty sand during depressurization under constant effective stress[J]. Geophysical Research Letters, 2021, 48(14): e2021GL092860. doi: 10.1029/2021GL092860
[8]	袁思敏, 王路君, 朱斌, 等. 考虑固相分解的含水合物沉积物体积应变分析模型[J]. 岩土工程学报, 2022, 44(6): 1044-1052. doi: 10.11779/CJGE202206008 YUAN Simin, WANG Lujun, ZHU Bin, et al. Volumetric strain analysis model for gas hydrate-bearing sediment considering effects of hydrate dissociation[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(6): 1044-1052. (in Chinese) doi: 10.11779/CJGE202206008
[9]	WU Z R, LI Y H, SUN X, et al. Experimental study on the effect of methane hydrate decomposition on gas phase permeability of clayey sediments[J]. Applied Energy, 2018, 230: 1304-1310. doi: 10.1016/j.apenergy.2018.09.053
[10]	蒋明镜, 陈意茹, 卢国文. 一种实用型深海能源土多场耦合离散元数值方法[J]. 岩土工程学报, 2021, 43(8): 1391-1398. doi: 10.11779/CJGE202108003 JIANG Mingjing, CHEN Yiru, LU Guowen. A practical multi-field coupling distinct element method for methane hydrate bearing sediments[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(8): 1391-1398. (in Chinese) doi: 10.11779/CJGE202108003
[11]	SANTAMARINA J C, DAI S, TERZARIOL M, et al. Hydro-bio-geomechanical properties of hydrate-bearing sediments from Nankai Trough[J]. Marine and Petroleum Geology, 2015, 66: 434-450. doi: 10.1016/j.marpetgeo.2015.02.033
[12]	CHONG Z R, YANG S H B, BABU P, et al. Review of natural gas hydrates as an energy resource: prospects and challenges[J]. Applied Energy, 2016, 162: 1633-1652. doi: 10.1016/j.apenergy.2014.12.061
[13]	YE J L, QIN X W, XIE W W, et al. The second natural gas hydrate production test in the South China Sea[J]. China Geology, 2020, 3(2): 197-209. doi: 10.31035/cg2020043
[14]	WANG Y, FENG J C, LI X S, et al. Large scale experimental investigation on influences of reservoir temperature and production pressure on gas production from methane hydrate in sandy sediment[J]. Energy & Fuels, 2016, 30(4): 2760-2770.
[15]	李桂琴, 李刚, 陈朝阳, 等. 多孔介质中甲烷水合物不同分解方法实验研究[J]. 化工进展, 2013, 32(6): 1230-1235. LI Guiqin, LI Gang, CHEN Zhaoyang, et al. Experimental investigation on different methods of dissociation of methane hydrate in porous sediment[J]. Chemical Industry and Engineering Progress, 2013, 32(6): 1230-1235. (in Chinese)
[16]	李彦龙, 刘昌岭, 廖华林, 等. 泥质粉砂沉积物—天然气水合物混合体系的力学特性[J]. 天然气工业, 2020, 40(8): 159-168. LI Yanlong, LIU Changling, LIAO Hualin, et al. Mechanical properties of the mixed system of clayey-silt sediments and natural gas hydrates[J]. Natural Gas Industry, 2020, 40(8): 159-168. (in Chinese)
[17]	WU Z R, LIU W G, ZHENG J N, et al. Effect of methane hydrate dissociation and reformation on the permeability of clayey sediments[J]. Applied Energy, 2020, 261: 114479. doi: 10.1016/j.apenergy.2019.114479
[18]	MCCALLUM S D, RIESTENBERG D E, ZATSEPINA O Y, et al. Effect of pressure vessel size on the formation of gas hydrates[J]. Journal of Petroleum Science and Engineering, 2007, 56(1/2/3): 54-64.
[19]	NAGAO J. Development of methane hydrate production method[J]. Synthesiology, 2012, 5(2): 89-97. doi: 10.5571/synth.5.89
[20]	ZHANG X H, LU X B, SHI Y H, et al. Centrifuge experimental study on instability of seabed stratum caused by gas hydrate dissociation[J]. Ocean Engineering, 2015, 105: 1-9. doi: 10.1016/j.oceaneng.2015.06.006
[21]	CHEN Y M, TANG Y, LING D S, et al. Hypergravity experiments on multiphase media evolution[J]. Science China Technological Sciences, 2022, 65(12): 2791-2808. doi: 10.1007/s11431-022-2125-x
[22]	KWON T H, OH T M, CHOO Y W, et al. Geomechanical and thermal responses of hydrate-bearing sediments subjected to thermal stimulation: physical modeling using a geotechnical centrifuge[J]. Energy & Fuels, 2013, 27(8): 4507-4522.
[23]	NG C W W, BAGHBANREZVAN S, LAU S Y, et al. Effects of hydrate dissociation on vertical casing–sediment interaction in carbon dioxide hydrate–bearing sand: novel In-flight centrifuge modeling[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2022, 148(3): 04021199. doi: 10.1061/(ASCE)GT.1943-5606.0002732
[24]	ZHU B, WANG L J, YANG S Q, et al. Pressure-control Temperature-control Hypergravity Experimental Device for Simulating Deep-sea Seabed Responses: US20210263005[P]. 2021-08-26.
[25]	朱斌, 王路君, 杨颂清, 等. 天然气水合物降压开采超重力模拟系统: CN108490151A[P]. 2020-07-03. ZHU Bin, WANG Lujun, YANG Songqing, et al. Supergravity simulation system for depressurization production of natural gas hydrates: CN108490151A[P]. 2020-07-03. (in Chinese)
[26]	朱斌, 王路君, 杨颂清, 等. 超重力条件下的天然气水合物热激法开采模拟装置: CN108386164A[P]. 2019-07-16. ZHU Bin, WANG Lujun, YANG Songqing, et al. Natural Gas Hydrate Heat Shock Method Mining Simulation Device under Supergravity Condition: CN108386164A[P]. 2019-07-16. (in Chinese)
[27]	杨颂清. 海底水合物开采多场相互作用数值分析与超重力装置研发[D]. 杭州: 浙江大学, 2020. YANG Songqing. Numerical Analyses on Multi-field Interactions in Soils during Marine Gas Hydrates Exploitation and Design of Centrifuge Modelling Device[D]. Hangzhou: Zhejiang University, 2020. (in Chinese)
[28]	王路君, 朱斌, 王心博, 等. 可模拟真实地层中天然气水合物储层的模型制备装置及方法: CN113072990A[P]. 2022-03-01. WANG Lujun, ZHU Bin, WANG Xinbo, et al. Model Preparation Device and Method Capable of Simulating Natural Gas Hydrate Reservoir in Real Stratum: CN113072990A[P]. 2022-03-01. (in Chinese)
[29]	CHEN Y, KONG L, ZHOU Y, et al. Development of a large geotechnical centrifuge at Zhejiang university[M]//Physical Modelling in Geotechnics. Zurich: CRC Press, 2010.
[30]	KNEAFSEY T J, TOMUTSA L, MORIDIS G J, et al. Methane hydrate formation and dissociation in a partially saturated core-scale sand sample[J]. Journal of Petroleum Science and Engineering, 2007, 56(1/2/3): 108-126.
[31]	ITO T, KOMATSU Y, FUJII T, et al. Lithological features of hydrate-bearing sediments and their relationship with gas hydrate saturation in the eastern Nankai Trough, Japan[J]. Marine and Petroleum Geology, 2015, 66: 368-378. doi: 10.1016/j.marpetgeo.2015.02.022
[32]	YAMAMOTO K, WANG X X, TAMAKI M, et al. The second offshore production of methane hydrate in the Nankai Trough and gas production behavior from a heterogeneous methane hydrate reservoir[J]. RSC Advances, 2019, 9(45): 25987-26013. doi: 10.1039/C9RA00755E

施引文献

资源附件(1)

其他相关附件
- PDF格式
  11-202402-G22-1300审稿意见与作者答复点击下载(325KB)

图(14)

计量

文章访问数: 349
HTML全文浏览量: 34
PDF下载量: 84
被引次数: 0

0. 引言
1. 纳米压痕试验
1.1 测试材料
1.2 纳米压痕仪试验原理
1.3 弹性模量计算方法
2. 测试结果验证
3. 离散元数值仿真
4. 讨论
5. 结论

水合物开采模拟超重力试验装置的研发及应用 English Version

作者简介: 王路君(1985— )，男，博士，副教授，从事水合物沉积物力学特性和土体多相多场耦合等方面的研究工作。E-mail: lujunwang@zju.edu.cn

通讯作者: 王鹏, E-mail: 12012041@zju.edu.cn

计量

出版历程