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

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

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
计量
- 文章访问数: 313
- HTML全文浏览量: 29
- PDF下载量: 82
出版历程
- 收稿日期: 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. 引言

随着筑坝材料的最大粒径越来越大,高坝筑坝料的最大粒径已经超过800 mm^[1],常规土力学试验设备无法直接用于研究原筑坝料的变形特性^[2-3]。为了解决试验粒径限制的问题,通常将原级配堆石材料进行缩尺得到替代原级配材料,然后对缩尺料进行试验,进而确定原型坝的变形特性^[4]。然而,一些研究成果及监测资料都表明,缩尺前后粗粒料的力学特性有较大差异,缩尺后粗粒料的力学特性并不能完全反应原级配粗粒料的力学性质^[5]。因此,如何确定原级配筑坝材料力学参数,是准确预测坝体变形的一个重要因素。

对此,国内外学者进行大量的大尺寸三轴试验。Marachi等^[6]对不同最大粒径的堆石料进行三轴试验得出,大粒径较小粒径压缩性更大。Hunter等^[7]积累了大量的不同粒径堆石料坝体的变形监测结果,提出堆石料的压缩模量是随粒径的增大而减小的。Varadarajan等^[8-9]对堆石料和砂砾石料分别进行不同缩尺条件下的对比试验,得出堆石料模量随粒径增加而降低,而砂砾石料的模量则随粒径增大而增加。汪小刚^[10]根据堆石料在压缩中可破碎和砂砾石料不可破碎的特征,进行了石膏和钢材的不同尺度下的压缩试验,得出与Varadarajan相同的规律。

不同文献所反映的缩尺试验结果可能与缩尺方法、密度、侧限条件、荷载及母岩等有关,但笔者认为将筑坝材料整体上分为可破碎的堆石料和不可破碎的砂砾石料分开研究是理清目前尺度效应研究的途径,同时为避免各种因素的交叉影响,本文基于接触变形理论、室内多尺度试验,研究不同粒径尺度、不同试样尺度、不同侧限条件、不同干密度、不同破碎条件、不同荷载等级等多种组合影响下变形模量的演化规律和内在机理,以此为基础建立多尺度变形参数的统一修正模型,为解决坝体变形预测不准确提供新的思考方法与解决途径。

1. 颗粒材料尺度效应理论分析

1.1 基本理论

（1）基本假定

以理想弹性力学作为计算理论基础；接触表面的摩擦力可忽略不计。

（2）变形方程

点接触物体受力后其接触表面为椭圆；线接触物体受力后其接触表面为矩形。符合变形连续条件。

（3）物理方程

服从弹性虎克定律,接触表面上应力的变化规律与接触体的应变成线性关系,应变最大的接触面中心压应力最大。根据接触表面压应力分布规律求得表面接触压力所组成的合力应等于外加载荷。

1.2 弹性体接触变形分析

设两圆球体的半径分别为R₁,R₂,开始时在公切平面上的O点相互接触,如图1所示。

图 1 弹性体接触变形

Figure 1. Contact deformation of elastomer

下载: 全尺寸图片幻灯片

在两球的子午面截线上,与z₁和z₂相距距离r处的两点M和N两点的距离为

$Z_{1} + Z_{2} = \frac{1}{2} (\frac{1}{R_{1}} + \frac{1}{R_{2}})$ 。

(1)

设w₁表示球体1面上的点M由于局部变形所产生的沿Z₁轴方向的位移,w₂表示球2面上的点N由于局部变形所产生的沿Z₂轴方向的位移,两球的中心O₁、O₂彼此接近的距离为δ。由于对称性,由接触产生的压力q和位移w对于接触中心O都是轴对称的,计算出最终。由此可以计算出：

$q_{0} = \frac{3 P}{2 π a^{2}}$ ,

(2)

$a^{3} = \frac{3}{4} \cdot \frac{R_{1} R_{2}}{R_{1} + R_{2}} (\frac{1 - μ_{1}^{2}}{E_{1}} + \frac{1 - μ_{2}^{2}}{E_{2}}) p$ ,

(3)

$δ^{3} = \frac{9}{16} \cdot \frac{R_{1} + R_{2}}{R_{1} R_{2}} {(\frac{1 - μ_{1}^{2}}{E_{1}} + \frac{1 - μ_{2}^{2}}{E_{2}})}^{2} p^{2}$ 。

(4)

式中 q₀为接触中心的压力最大；a是接触圆半径；E₁,μ₁,E₂,μ₂分别是两个弹性体弹性常数；p为压力分布函数；δ为变形量；R₁,R₂分别是两个弹性体半径。

由上式得出,最大接触压应力与载荷不是线性关系,而是与载荷的立方根成正比,这时候因为随着载荷的增加,接触面积也在增大,其结果使接触面上的最大压应力的增长较载荷的增长为慢。

当弹性体与平面接触变形时,则可以将R₁=R,R₁=∞,E₁=E₂=E,μ₁=μ₂=0.3代入式（2）～（4）求解,得到下式：

$q_{0} = 0 .388 \sqrt[3]{\frac{P E^{2}}{R^{2}}}$ ,

(5)

$a = 1 .109 \sqrt[3]{\frac{P R}{E}}$ ,

(6)

$δ = 1 .231 \sqrt[3]{\frac{P^{2}}{E^{2} R}}$ 。

(7)

当计算两等直径弹性体接触时,则可将R₁=R,R₁=∞,E₁=E₂=E,μ₁=μ₂=0.3代入式（2）～（4）求解。

$a = 1 .109 \sqrt[3]{\frac{P R}{E}}$ ,

(8)

$δ = 1 .231 \sqrt[3]{\frac{2 P^{2}}{E^{2} R}}$ 。

(9)

1.3 粒径尺度对压缩模量影响研究

首先研究固定试样尺寸,改变材料的最大粒径情况下对压缩模量的影响,从材料的粒径改变揭示压缩模量的变化规律。推导假定为：在相同加载试样尺寸内,材料从单个半径为R的弹性体逐渐减小其半径,分别为R/2,R/4,R/6,…,R/1000,如图2所示。在此过程中密度保持不变,空间内的孔隙率保持不变,荷载等级保持不变,无侧限条件,受荷载P的作用下,逐步分析粒径逐渐缩小所引起的弹性模量的变化。

图 2 最大粒径改变示意

Figure 2. Schematic diagram of change of maximum particle size

下载: 全尺寸图片幻灯片

利用1.2节理论推导结果,计算图2中不同粒径下的变形量,进而计算更小的最大粒径下变形量,形成归一化压缩模量。此处定义归一化压缩模量含义为：设未进行缩尺时试件的压缩模量初始值为1,缩尺后压缩模量与初始压缩模量之比为归一化压缩模量。

为更直观表达级配组构条件下粒径尺度对压缩模量的影响,将初始最大粒径设为1000 mm,对应的压缩模量归一化为1,并化为无量刚量,则随着粒径的减小,建立如图3所示关系。

图 3 粒径尺度对压缩模量影响

Figure 3. Influences of particle size on compression modulus

下载: 全尺寸图片幻灯片

建立粒径尺度与归一化压缩模量的关系如下式所示：

$E_{gd} = 0 .5691+0 {.2116e}^{- \frac{ξ}{9 .9617}} +0 {.2374e}^{- \frac{ξ}{212 .3621}}$ ,

(10)

式中,E_gd为粒径尺度归一化模量,ξ为固定试样尺寸条件下的归一化粒径。

1.4 试件尺度对压缩模量影响研究

对于有级配组构情况下,首先确定组构规则,根据圆形几何组构规则,空隙部分首先用最大半径进行组构,组构后每个空间为三边形,在此三边空隙内继续组构,直至R/1000。组构弹性体的粒径符合指数递减规律,并形成悬浮型颗粒构架体系,同时构建最接近实际情况下的级配曲线,悬浮级配构建如图4所示,与之对应的级配曲线如图5。

图 4 D_max=0.6D时排列状态

Figure 4. Arrangement state of D_max = 0.6D

下载: 全尺寸图片幻灯片

图 5 D_max=0.6D的级配曲线

Figure 5. Grading curves of D_max = 0.6D

下载: 全尺寸图片幻灯片

为从材料的试样尺度角度揭示加载试验的面积改变对压缩模量的影响,逐渐增大试样尺度,分别为2R,4R,6R,…,1000R,求解在不同荷载加载面积下的压缩模量,如图6所示。

图 6 试样尺度逐渐增大示意

Figure 6. Increase of sample size

下载: 全尺寸图片幻灯片

根据不同试样尺度条件下的归一化压缩模量,引入径径比的概念,既试样尺寸与最大粒径之比,建立径径比与压缩模量的关系如图7所示。可得出,当径径比大于5时,归一化压缩模量变稳定,从理论角度论证了压缩试验可靠性的必要条件。

图 7 压缩模量与径径比的关系

Figure 7. Relationship between compression modulus and diameter diameter ratio

下载: 全尺寸图片幻灯片

2. 变形特性的外部影响因素研究

2.1 侧限条件对压缩模量的影响

对侧向不允许变形的压缩弹性体,取一个单元进行分析其弹性阶段受力变形。该单元受到三向应力 $σ_{x}$ , $σ_{y}$ , $σ_{z}$ 作用,在z方向的压力作用下,试样中的竖向有效应力为 $σ_{z}$ ,又因为试样的受力条件属于轴对称问题,故水平向有效应力 $σ_{x}$ = $σ_{y}$ ,由于土体的侧向变形被限制住,故 $ε_{x}$ = $ε_{y}$ = 0,于是：

$ε_{x} = \frac{σ_{x}}{E_{s}} - μ \frac{σ_{y}}{E_{s}} - μ \frac{σ_{z}}{E_{s}} = 0$ ,

(11)

式中,E,E_s分别为无侧线和完全侧限条件下的模量,μ为组构材料的泊松比。μ取0.3时,建立完全侧向与无侧限模量的关系为

$E_{s} = 1 .346 E$ 。

(12)

由此可得完全侧向与无侧限模量随粒径整体的增大倍数,在μ取0.3时增大约35%。

2.2 干密度与压缩模量的关系研究

本研究推导成立的边界条件是从原型尺度到实验室尺度均保持同样的干密度或孔隙比,而在实际室内试验中,由于粒径的缩小通常表现为筑坝原型干密度大,而室内实验试件的干密度小,且很难使压缩试件的干密度达到现场的大小,这就为实现同密度或孔隙率条件下的室内试验造成一定的障碍。为解决这一障碍,找出室内试验密度或孔隙率与压缩模量的关系,采用此规律修正后即便在密度未达到原型密度下亦可进行实验并进行回溯。

对若干个工程的材料均进行不同干密度下的压缩试验,建立如图8所示趋势图。

图 8 密度与压缩模量关系

Figure 8. Relationship between density and compression modulus

下载: 全尺寸图片幻灯片

建立如下式所示的关系。

$f (ρ_{d}) = 2 .331+ \frac{0 .702 - 2 .331}{1+ {(\frac{ρ_{d}}{1 .996})}^{29 .891}}$ ,

(13)

式中, $f (ρ_{d})$ 为受密度影响下试样的压缩模量, $ρ_{d}$ 为干密度。

2.3 颗粒边缘破碎影响分析

为分析堆石料破碎与粒径间的关系,在单一粒径情况下假定大小颗粒的粒径在特定比例下,某一应力状态中发生破碎,保持该应力状态不变,不断减小最大粒径。计算基本参数为：大粒径与小粒径之比为20,在缩小大粒径过程中,保持该比例不变。

图9不同粒径下的试样破碎变形计算。由式（2）～（4）可计算小粒径上的最大应力。在粒径不断缩小的过程中,单个小颗粒承受的外部应力会逐渐减小,化为等效破碎变形量,建立破碎量与粒径变化如图10所示。

图 9 不同粒径下试样破碎变形计算

Figure 9. Calculation of crushing deformation of samples with different particle sizes

下载: 全尺寸图片幻灯片

图 10 粒径与破碎变形量的关系

Figure 10. Relationship between particle size and crushing deformation

下载: 全尺寸图片幻灯片

该曲线采用函数表示为

$δ = 1 {.023e}^{- \frac{D}{42 .848}}$ 。

(14)

2.4 荷载等级对压缩模量影响研究

为了解荷载等级对压缩模量的影响关系,在有填充条件下,分别计算半径为R/n的弹性体,在有填充作用下分别受荷载2p,3p,4p,…,10p作用下的压缩模量。压缩模量的增大倍数仅与荷载大小直接相关,绘制压缩模量与荷载等级趋势关系如图11所示。并可用下式：

图 11 荷载等级与压缩模量的倍数关系

Figure 11. Relationship between load grade and compression modulus

下载: 全尺寸图片幻灯片

$f (p) = 2.7354 - 1.9981 e^{- \frac{P_{i}}{6.4893}}$ ,

(15)

式中, $f (p)$ 为荷载影响下的压缩模量放大倍数,p_i为荷载等级倍数,为当前加压等级与初始加压等级的比值。

3. 堆石料多尺度统一修正模型的建立

3.1 模型的提出

多尺度统一修正模型在于将缩尺材料的压缩模量求解为筑坝原型材料的压缩模量,在此过程中综合考虑缩尺过程、荷载等级、干密度和侧限状态、破碎效应的影响,对此提出涵盖此参数的统一修正模型公式为

$F (ξ, p, ρ_{d}, μ, D) = f_{1} (ξ) \cdot f_{2} (p_{i}) \cdot f_{3} (ρ_{d}) \cdot f_{4} (μ) \cdot f_{5} (D)$ 。

(16)

式中 $f_{1} (ξ)$ 为粒径影响下的的回溯函数； $f_{2} (p_{i})$ 为荷载等级影响下的修正函数； $f_{3} (ρ_{d})$ 为干密度影响下的修正函数； $f_{4} (μ)$ 为泊松比影响下的侧限修正函数。 $f_{5} (D)$ 为粒径影响下的破碎修正函数。根据回溯过程,建立粒径尺度和试样尺度影响下归一化压缩模量的关系,如图12所示。

图 12 粒径影响下的的回溯函数图

Figure 12. Backtracking function under influence of particle size

下载: 全尺寸图片幻灯片

具体函数关系为

$\begin{matrix} f_{1} (ξ) = \frac{A_{1} + \frac{A_{2} - A_{1}}{1 + {(ξ_{1} / A_{3})}^{A_{4}}}}{A_{1} + \frac{A_{2} - A_{1}}{1 + {(ξ_{0} / A_{3})}^{A_{4}}}} \cdot \frac{B_{1} + B_{2} e^{- \frac{ξ_{2}}{B_{3}}} + B_{4} e^{- \frac{ξ_{2}}{B_{5}}}}{B_{1} + B_{2} e^{- \frac{ξ_{1}}{B_{3}}} + B_{4} e^{- \frac{ξ_{1}}{B_{5}}}} \cdot \\ \frac{A_{1} + \frac{A_{2} - A_{1}}{1 + {(ξ_{3} / A_{3})}^{A_{4}}}}{A_{1} + \frac{A_{2} - A_{1}}{1 + {(ξ_{2} / A_{3})}^{A_{4}}}} \end{matrix}$ 。

(17)

式中 A_i表示固定材料最大粒径下,试样尺寸对归一化压缩模量的影响系数；B_i表示固定试样尺寸条件下材料粒径变化与归一化压缩模量的系数；ξ为转换过程的粒径。

3.2 基于外部影响因素的统一修正

（1）泊松比影响下的侧限修正

侧限条件是筑坝料变形应力状态的具体描述,侧限所引起的压缩模量改变可通过泊松比表示,泊松比对压缩模量的影响关系,可用下式表示：

$f_{4} (μ) = 1 - \frac{2 μ^{2}}{1 - μ}$ ,

(18)

式中,μ为筑坝料的泊松比,可通过压缩试验获取。

（2）干密度影响下的修正

由2.4节所提出的干密度对变形模量的影响模型,可得到在未达到原型密度下的回溯。如下式所示：

$f_{3} (ρ_{d}) = D_{1} + \frac{D_{2} - D_{1}}{1 + {(ρ_{d} / D_{3})}^{D_{4}}}$ ,

(19)

式中,D_i(i=1,2,3,4,)为待定系数, $ρ_{d}$ 为干密度。

（3）破碎效应影响修正

筑坝材料整体上可分为可压缩破碎的堆石料和不破碎的砾石料,当采用堆石料筑坝时,往往表现出粒径越大压缩模量越低的情况,由此引入破碎效应的修正,如下式。当采用砂砾石料时不发生材料颗粒破碎,可不考虑破碎效应。

$f_{5} (D) = α e^{D / r}$ ,

(20)

式中,D为最大粒径, $α$ , $r$ 为待定系数。

（4）荷载等级影响下的修正函数

对于筑坝料,其压缩模量与试验施加的荷载应力大小有关,表现为荷载应力水平越大压缩模量越大,呈非线性正相关,具体如下式表示：

$f_{2} (p_{i}) = C_{1} + C_{2} e^{C_{3} ρ_{d}}$ ,

(21)

式中,C_i为受荷载等级影响下的压缩模量关系的非线性表达系数。

综合以上各影响因素,选取实际试验中所涉及到的情况进行选用。

3.3 多尺度统一修正模型的应用

以河口村大坝为例,坝体为混凝土面板堆石坝,最大坝高122.5 m,面板堆石坝设计主要考虑的因素是大坝的应力变形状态,而大坝应力变形分析的主要手段是有限元数值分析。将室内缩尺三轴试验所得E-B模型参数进行的坝体沉降计算与经统一修正模型修正后的模型参数计算结果进行对比,并与工程监测值相校核。

主要考虑主堆石料计算参数的影响,主堆石区填筑材料为级配碎石,室内三轴试验所得主应力差与轴向应变如图13所示。根据该三轴试验结果计算主堆石材料的E-B模型参数如表1所示。采用该参数计算所得最大竖向沉降竣工时61.2 cm,蓄水期68.1 cm。

图 13 主堆石料三轴试验

Figure 13. Triaxial tests on main rockfill

下载: 全尺寸图片幻灯片

表 1 修正前后河口村主堆石材料E-B模型参数

Table 1. 　E-B model parameters of main rockfill materials in Hekou Village

分类	干重度/(kN·m^-3)	K	n	R_f	K_ur	c/kPa	φ₀/(°)	Δφ/(°)	K_b	m
修正前	20.7	1660	0.21	0.85	3320	0	54	10.6	380	0.14
修正后	20.7	1312	0.25	0.82	3320	0	54	10.6	581	0.14

下载: 导出CSV

| 显示表格

本文的模型最要考虑变形参数的影响,E-B模型中与变形相关的过程量为E_i,根据三轴试验过程和材料属性,采用统一修正模型进行尺度及各影响因素的影响系数计算,如表2所示。

表 2 河口村主堆石材料修正计算

Table 2. Correction calculation of main rockfill materials in Hekou Village

项目	f₁(ξ)	f₂(μ)	f₃(ρ_d)	f₄(D)	f₅(p)	综合影响F(·)
修正值	1.19	1.0	1.47	0.56	1.0	0.84

下载: 导出CSV

| 显示表格

采用表2中修正后的E-B模型参数,计算总体沉降量对比如表3所示。

表 3 河口村坝体沉降计算对比

Table 3. Comparison of settlement calculation of Hekou Village Dam

项目	计算最大沉降/cm		提升准确度/%
项目	竣工	蓄水	提升准确度/%
室内缩尺试验	61.2	68.1
监测数据	78.9	109.7	17.43~28.89
统一修正模型	75.6	99.7

下载: 导出CSV

| 显示表格

4. 结论

针对筑坝料的变形特性的尺寸效应,以接触变形理论为基础,从粒径尺度、试件尺度及外部影响因素角度,结合试验的方法进行了研究,主要得到以下结论：

（1）基于弹性力学接触变形理论,建立了粒径尺度与归一化压缩模量的指数相关模型,论证了试件尺度中径径比大于5时压缩试验的变形可靠性。

（2）揭示了与尺度效应相关的5个外部因素,侧限条件、干密度、边缘破碎及荷载等级等多种组合对变形特性的影响规律和内在机理,并分别建立了各因素与变形的相关模式。

（3）基于试验回溯现场过程,建立了综合考虑缩尺过程、侧限条件、干密度、边缘破碎及荷载等级因素下的多尺度变形参数统一修正模型。

（4）通过对大坝堆石料变形E-B模型计算参数的修正及变形计算,验证了多尺度变形参数统一修正模型的适用性与可靠性,较直接采用缩尺试验结果准确度提升了约20%。

图 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

施引文献(49)

期刊类型引用(21)

1.	孟珂，王笑梅，杜晓冉，张晓曼，罗娟. 罕见多矿物晶体共生标本的综合鉴定. 矿产综合利用. 2025(01): 200-205 . 百度学术
2.	刘勇，张志康，魏建平，徐向宇，郜英俊. 柔性刀具冲击破煤能量演化及关键参数. 煤炭学报. 2025(02): 965-974 . 百度学术
3.	李红丽. 非均质岩石单轴压缩下损伤演化规律数值模拟研究. 有色矿冶. 2024(01): 43-48 . 百度学术
4.	鞠明和，陶泽军，蔚立元，姜礼杰，郑彦龙，邹春江. 钢粒子迟滞重复冲击破岩硬岩损伤破裂特征研究. 岩土力学. 2024(04): 1242-1255 . 百度学术
5.	吴泽兵，袁若飞，张文溪，王刚，胡诗尧. PDC混合布齿钻头破碎非均质花岗岩数值模拟. 天然气工业. 2024(05): 105-117 . 百度学术
6.	周元，吕威帆，王颖轶. 基于块体离散元法的盾构掘进围岩与管片变形模拟研究. 都市快轨交通. 2024(03): 125-134 . 百度学术
7.	裴书锋，郝文锋，王营利，王一汀，曾凤娟. 双江口水电站花岗岩单轴压缩微观破坏机制研究. 西北水电. 2024(04): 62-68 . 百度学术
8.	张国桥，孙鹏，吴祥业，王婧雅，郭文斌，田宇航. 基于PFC-GBM非均质模型的砂岩裂纹演化细观规律研究. 中国矿业. 2024(09): 158-169 . 百度学术
9.	冯龙飞，王双明，王晓东，解嘉豪，窦林名. 煤单轴峰后动态冲击破坏特征及差异机制模拟研究. 煤炭学报. 2024(S2): 714-730 . 百度学术
10.	马文强，王酒婷. 花岗岩受压宏-细观破坏特征及能量演化规律. 信阳师范学院学报(自然科学版). 2023(02): 314-320 . 百度学术
11.	张涛，蔚立元，苏海健，高亚楠，贺虎，魏江波. 基于多级力链网络分析的花岗岩压缩特性的矿物尺寸效应研究. 岩石力学与工程学报. 2023(08): 1988-2003 . 百度学术
12.	乔世范，刘钰，王刚，张细宝，张海凤，董常瑞，谭晶仁，檀俊坤. 考虑岩石细观结构的TBM滚刀破岩过程数值研究. 中国安全生产科学技术. 2023(07): 106-112 . 百度学术
13.	向衍斌. 煤系岩石单轴压缩损伤破坏演化规律与表征. 煤矿安全. 2023(09): 88-95 . 百度学术
14.	赵光明，高宇，吴旭坤. 岩石变刚度实验条件下力学与声发射特性. 安徽理工大学学报(自然科学版). 2023(06): 63-72 . 百度学术
15.	王桂林，王润秋，孙帆. 块体离散元颗粒模型细观参数标定方法及花岗岩细观演化模拟. 长江科学院院报. 2022(01): 86-93 . 百度学术
16.	张涛，蔚立元，鞠明和，李明，苏海健，季浩奇. 基于PFC3D-GBM的晶体–单元体尺寸比对花岗岩动态拉伸特性影响分析. 岩石力学与工程学报. 2022(03): 468-478 . 百度学术
17.	兰恒星，包含，孙巍锋，刘世杰. 岩体多尺度异质性及其力学行为. 工程地质学报. 2022(01): 37-52 . 百度学术
18.	李博，梁秦源，周宇，赵程，伍法权. 基于CT-GBM重构法的花岗岩裂纹扩展规律研究. 岩石力学与工程学报. 2022(06): 1114-1125 . 百度学术
19.	Tongzhao Zhang，Hongguang Ji，Xiaobo Su，Shuang You，Daolu Quan，Zhou Zhang，Jinzhe Li. Evaluation and classification of rock heterogeneity based on acoustic emission detection. International Journal of Minerals, Metallurgy and Materials. 2022(12): 2117-2125 . 必应学术
20.	郑强强，徐颖，胡浩，钱佳威，宗琦，谢平. 单轴荷载作用下砂岩的破裂与速度结构层析成像. 岩土工程学报. 2021(06): 1069-1077 . 本站查看
21.	李博，朱强，张丰收，赵程，伍法权. 基于矿物晶体模型的非均质性岩石双裂纹扩展规律研究. 岩石力学与工程学报. 2021(06): 1119-1131 . 百度学术

其他类型引用(28)

资源附件(1)

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

图(14)

计量

文章访问数: 313
HTML全文浏览量: 29
PDF下载量: 82
被引次数: 49

0. 引言
1. 颗粒材料尺度效应理论分析
1.1 基本理论
1.2 弹性体接触变形分析
1.3 粒径尺度对压缩模量影响研究
1.4 试件尺度对压缩模量影响研究
2. 变形特性的外部影响因素研究
2.1 侧限条件对压缩模量的影响
2.2 干密度与压缩模量的关系研究
2.3 颗粒边缘破碎影响分析
2.4 荷载等级对压缩模量影响研究
3. 堆石料多尺度统一修正模型的建立
3.1 模型的提出
3.2 基于外部影响因素的统一修正
3.3 多尺度统一修正模型的应用
4. 结论

0. 引言
1. 颗粒材料尺度效应理论分析
1.1 基本理论
1.2 弹性体接触变形分析
1.3 粒径尺度对压缩模量影响研究
1.4 试件尺度对压缩模量影响研究
2. 变形特性的外部影响因素研究
2.1 侧限条件对压缩模量的影响
2.2 干密度与压缩模量的关系研究
2.3 颗粒边缘破碎影响分析
2.4 荷载等级对压缩模量影响研究
3. 堆石料多尺度统一修正模型的建立
3.1 模型的提出
3.2 基于外部影响因素的统一修正
3.3 多尺度统一修正模型的应用
4. 结论

参考文献(32)

施引文献

资源附件(1)

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

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

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

计量

出版历程

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

0. 引言

1. 颗粒材料尺度效应理论分析

1.1 基本理论

1.2 弹性体接触变形分析

1.3 粒径尺度对压缩模量影响研究

1.4 试件尺度对压缩模量影响研究

2. 变形特性的外部影响因素研究

2.1 侧限条件对压缩模量的影响

2.2 干密度与压缩模量的关系研究

2.3 颗粒边缘破碎影响分析

2.4 荷载等级对压缩模量影响研究

3. 堆石料多尺度统一修正模型的建立

3.1 模型的提出

3.2 基于外部影响因素的统一修正

3.3 多尺度统一修正模型的应用

4. 结论

期刊类型引用(21)

其他类型引用(28)

其他相关附件

PDF格式

计量

出版历程

目录

0. 引言

1. 颗粒材料尺度效应理论分析

1.1 基本理论

1.2 弹性体接触变形分析

1.3 粒径尺度对压缩模量影响研究

1.4 试件尺度对压缩模量影响研究

2. 变形特性的外部影响因素研究

2.1 侧限条件对压缩模量的影响

2.2 干密度与压缩模量的关系研究

2.3 颗粒边缘破碎影响分析

2.4 荷载等级对压缩模量影响研究

3. 堆石料多尺度统一修正模型的建立

3.1 模型的提出

3.2 基于外部影响因素的统一修正

3.3 多尺度统一修正模型的应用

4. 结论

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

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