A practical multi-field coupling distinct element method for methane hydrate bearing sediments
-
摘要: 水合物储量丰富且污染小,是全球关注的未来绿色能源。而水合物的大规模开采将致使水合物沉积物（即能源土）力学性质劣化,从而引发出砂、海底滑坡等工程和地质灾害,制约水合物的安全高效开采。针对目前大多数多场耦合分析框架无法准确描述颗粒运移等大变形问题,基于离散元商业计算软件Particle Flow Code (PFC)和含水合物多孔介质的多相流分析程序TOUGH+HYDRATE (T+H),构建了一种实用型水合物开采多场耦合分析框架T+H+PFC。通过将力学计算模块PFC所得孔隙率参数和T+H计算所得温度、压力、盐浓度信息交换、更新,实现了T+H与PFC的耦合计算。然后,通过一维固结和一维热传导试验验证了该耦合程序的有效性。最后,模拟分析了一维降压、升温开采下能源土动态响应的与水合物弱化过程。结果表明,耦合模拟与试验结果较为吻合,T+H+PFC在多场耦合下水合物开采引起的地层变形、出砂、海底滑坡等问题分析中具有良好的应用前景。Abstract: The methane hydrate (MH) has been attracting extensive attention as a promising green energy source because of the abundant reserve and more environmental friendliness. However, large-scale exploitation induces the weakening of MH bearing sediments (MHBS) structure, which can result in engineering problems such as sand jamming, and environmental disasters such as landslide, and this will constrain efficient and safe exploitation of MH. Given that several multi-field coupling simulators fail to accurately describe large deformation of soils such as particle migration, a practical multi-field method T+H+PFC for MH exploitation analysis is established based on the commercial distinct element method (DEM), software particle flow code (PFC) and numerical code TOUGH+HYDRATE (T+H) for the multiphase flow analysis of hydrate-bearing geologic systems. The coupling computation can be achieved by exchanging the information between the porosity obtained by DEM and the temperature, pressure and salinity obtained by T+H. Then, the method is validated by conducting numerical modeling of one-dimensional (1D) consolidation and 1D heat conduction tests. Finally, focus is on the dynamic responses of MHBS and weakening of MH under the conditions of 1D dissociation tests with depressurization and heat stimulating methods. The numerical results are found to be in good agreement with the experimental data available, which demonstrates that the T+H+PFC method has a great promise in the MH decomposition-induced multi-field coupling analysis, such as deformation of seabed, sand production and submarine landslide.
-
0. 引言
近年来，中国大力发展南海地区基础设施建设，钙质砂作为沿海地区的主要岩土体受到广泛使用。有别于陆源石英砂，钙质砂作为难溶类碳酸盐（多为碳酸钙）物质的集合体，具有形状不规则、内孔隙高、颗粒易破碎、棱角度高等性质[1-3]。国内外学者对于钙质砂的力学特性开展了大量研究：张晨阳等[4]分析了细粒含量对钙质砂力学特性的影响，认为钙质砂峰值偏应力与剪胀性随细粒含量的增加而逐渐减小；黄宏翔等[5]得出在相同条件下，钙质砂的抗剪强度和残余强度均大于石英砂；Rui等[6]发现不规则颗粒形状的钙质砂表现出较慢的累积孔隙压力速率，显示出更高的抗液化性。对于钙质砂的颗粒破碎特性，研究表明在低围压情况下就会发生颗粒破碎，颗粒形状对其力学性质的影响较大，在高围压情况下钙质砂颗粒破碎情况加剧，颗粒破碎对其力学性质影响愈加显著[7]。Coop等[8]指出钙质砂的颗粒破碎程度受轴向应力大小及初始颗粒形状影响；汪轶群等[9]在三轴试验中发现随着初始围压的增大，钙质砂颗粒破碎程度加大，土样整体剪胀趋势减小，且破碎引起的能量耗散增加。
上述研究为以钙质砂的颗粒特性为切入点，对其力学性质进行了研究。但主要以钙质砂作为纯净砂进行分析，而仅考虑某一粒径、某一种类砂土多与实际工程情况不符，存在一定的局限性。中国南海地区广泛分布源于岛礁中的碳酸盐土体和沿岸地区的非碳酸盐类土所结合的碳酸盐混合土[10]，其力学性质势必由两种土体共同决定。李小刚等[11]通过三轴试验发现，混合砂峰值强度和相对破碎率随石英砂含量的增加而减小，黏聚力及峰值内摩擦角随碳酸盐含量的增加而增加。Lade等[12]指出土体中粉土、粉砂等细颗粒的含量是海底斜坡、土坝和路堤等产生液化破坏的主要原因。含细颗粒土体的试验同样表明细颗粒含量、粒径分布和颗粒形状特征是影响砂土抗剪强度和液化反应的主要因素[13-15]，且细颗粒含量与粉土特性对砂土静态液化行为的影响存在耦合关系[16]。此外，砂土中不同颗粒特征的粗颗粒和细颗粒相互作用，对混合土体的整体力学性质也产生显著影响[17-18]。
综上所述，细颗粒对钙质混合土体力学性质影响显著，而现有研究中多数仅考虑钙质细粒含量对钙质砂力学特性的影响，在对钙质混合砂的研究中，细颗粒特征不同对混合砂的影响程度亦不同。基于此，本文采用角形石英粉和圆形玻璃珠作为细颗粒，分别与钙质砂混合，对钙质混合砂进行不同轴向压力条件下的环剪试验，探究不同细颗粒含量及形状对钙质混合砂峰值强度、残余强度和轴向变形的影响，对钙质混合砂的颗粒破碎特征进行研究分析，进一步探讨颗粒破碎特征对钙质混合砂性质的影响规律，从而为南海地区工程建设提供参考。
1. 试验材料及方案
1.1 试验材料
试验所用钙质砂取自南海区域某岛礁附近。试验前先将砂样烘干，采用标准筛进行筛分，得到粒径范围为0.1~1 mm的钙质砂颗粒。玻璃珠因具有与土体相似的矿物学特征[19]，多用于模拟砂土，其形状接近理想的球形；石英粉则更具棱角和不规则性。因此，采用粒径小于0.075 mm的石英粉及玻璃珠模拟混合砂中不同形状的石英细颗粒。3类基本砂样的物理指标参数见表 1，颗粒级配曲线见图 1。按试样总质量的0%，10%，20%，30%，40%，将石英粉、玻璃珠细颗粒分别与钙质砂均匀混合得到9种试样：含石英粉试样，由J10，J20，J30，J40表示；含玻璃珠试样，由Y10，Y20，Y30，Y40表示；0表示纯钙质砂试样。
表 1 基本砂样的物性指标Table 1. Physical properties of basic sand samples试样 圆度 宽长比 D10/
μmD50/
μmD60/
μm不均匀系数Cu 曲率系数Cc 颗粒相对质量密度 钙质砂 0.517 0.541 374.5 637.7 643.5 1.72 1.28 2.75 玻璃珠 0.924 0.893 35.1 50.8 58.6 1.67 0.92 2.64 石英粉 0.635 0.746 27.2 45.7 51.3 1.89 1.05 2.66 注：①圆度为单颗粒投影面积与等周长圆面积之比；宽长比为单颗粒最小Feret直径与最大Feret直径之比。②圆度、宽长比取采样（1000粒）均值。 ![]()
图 1 3种砂样颗粒级配曲线Figure 1. Grain-size distribution curves of three kinds of sand particles为定量化表示土体形态特征，多通过颗粒形状扫描仪器，根据颗粒轮廓图像建立二维形状参数（圆度、宽长比、粗糙度等），能够较好地表述颗粒形状与其力学性质之间的关系[20-22]。对于二元颗粒混合物，参考Yang等[23]提出的组合形状参数概念，认为二元混合物组合形状参数由各颗粒材料的形状特征值及其对应的质量占比决定，组合圆度、组合宽长比越接近于1，表明颗粒整体越接近于圆形、越规整。如钙质混合砂J30的组合圆度、组合宽长比计算分别为
Sc−J30=Sc−J×Fc+Sc−G×(1−Fc)=0.635×30%+0.517×70%=0.552, (1) Ar−J30=Ar−J×Fc+Ar−G×(1−Fc)=0.746×30%+0.541×70%=0.603 (2) 图像颗粒分析仪如图 2所示，扫描基本砂样轮廓图如图 3所示，混合砂样特性指标如表 2所示。
表 2 混合砂样特性指标Table 2. Characteristics of mixed sand samples砂样 相对密实度Dr/% 最大孔隙比emax 最小孔隙比emin 组合圆度Sc 组合宽长比Ar 0 60 1.109 0.812 0.517 0.541 J10 60 1.010 0.707 0.529 0.562 J20 60 0.893 0.553 0.541 0.582 J30 60 0.794 0.467 0.552 0.603 J40 60 0.743 0.430 0.564 0.623 Y10 60 0.977 0.724 0.558 0.541 Y20 60 0.856 0.605 0.598 0.576 Y30 60 0.744 0.512 0.639 0.611 Y40 60 0.685 0.491 0.680 0.647 1.2 试验方案
试验采用HJ-1型环剪仪，其试样尺寸为外径100 mm，内径60 mm，高20 mm。天然状态下土体的相对密实度在20%~80%，且天然密实度多为中密状态[24]，故考虑实际工程情况，控制钙质混合砂的相对密实度Dr =60%（中密），进行不同轴向压力条件下的环剪试验。试验步骤如下：①采用漏斗法将试样装入容器。首先，将环剪仪的上、下剪切盒固定在一起，然后将试样缓慢倒入漏斗之中，通过漏斗将砂样以零高度的落差缓慢倒入盒中，最后将试样抚平，使其高度一致。②采用浸水法饱和试样。试样装好后倒入蒸馏水，使其完全浸没，浸泡4 h确保试样饱和度达到要求。浸泡完成后将仪器组装完成。③逐级施加100，200，300 kPa的轴向压力，待固结完成后，以1°/min的剪切速率进行环向剪切试验。
2. 试验结果分析
2.1 剪切应力-剪切位移曲线分析
试验结束后，得到不同细颗粒含量及形状条件下的剪切应力-剪切位移曲线，如图 4所示。由图 4（a）可得，纯钙质砂的剪切应力随着剪切位移的增加而逐渐增加，呈近线性增长，随后，剪切应力的增长趋势变缓。当轴向压力为100，200，300 kPa时，剪切应力分别在剪切位移为10.5，11.9，13.6 mm时达到峰值强度，其值分别为100.4，177.9，251 kPa；剪切位移分别在28，30，31 mm处得到土样的残余强度，其值分别为88.23，159.4，213.51 kPa。由各应力-位移图线得：各组试样的抗剪强度均随轴向压力的增加而增大；当剪切位移相同时，轴向压力越大，试样的剪切应力越大；同时，随轴向压力的增加，各组试样峰值强度出现延后性，且各组试样的应力-位移曲线均呈现软化特性。
对不同细颗粒含量及形状条件下钙质混合砂的峰值强度进行统计并绘制图 5。如图 5（a）所示，随着细颗粒含量的增加，混合砂峰值强度逐渐降低。与纯钙质砂相比，细颗粒含量为10%时，混合砂的峰值强度略小于纯钙质砂，其中，含石英粉混合砂峰值强度最大降幅为12%，含玻璃珠混合砂峰值强度最大降幅为16%，此时细颗粒填充土体骨架的孔隙部分，土体骨架仍由钙质砂颗粒组成；细颗粒含量为20%时，含石英粉混合砂峰值强度最大降幅为17%，含玻璃珠混合砂峰值强度最大降幅为29%；细颗粒含量为30%时，峰值强度持续降低，含石英粉混合砂峰值强度最大降幅为22%，含玻璃珠混合砂峰值强度最大降幅扩大至36%；当细颗粒含量增加至40%，峰值强度发生大幅度降低，含石英粉混合砂峰值强度最大降幅为30%，含玻璃珠混合砂峰值强度最大降幅则高达47%，表明此时土体骨架由钙质砂颗粒承担转变为钙质砂、细颗粒共同承担。此外，如图 5（b），（c）所示，峰值强度受细颗粒形状影响，随组合圆度和组合宽长比增加而逐渐减小，含玻璃珠混合砂出现更大降幅。结合剪切应力-剪切位移曲线（图 4）中含玻璃珠混合砂的曲线波动更为明显，表明细颗粒形状的改变影响了骨架颗粒间的接触状态以及力链传递，进而降低土体骨架的稳定性，圆形细颗粒的加入使得土体骨架的稳定性更差。各组试样峰值强度及下降幅度如表 3所示。
![]()
图 5 不同细颗粒含量及形状下的峰值强度Figure 5. Peak intensities at different contents and shapes of fine particles表 3 不同轴向压力下试样的峰值强度及下降幅度Table 3. Peak strengths and decrease amplitudes of samples under different axial pressures轴压 类别 0 J10 J20 J30 J40 Y10 Y20 Y30 Y40 100 kPa 峰值强度/kPa 100.4 93.6 90.7 78.0 70.2 84.8 71.7 64.3 53.6 降幅/% — 7 10 22 30 16 29 36 47 200 kPa 峰值强度/kPa 177.9 156.5 148.7 142.3 139.4 148.7 133.6 120.9 105.3 降幅/% — 12 16 20 22 16 25 32 41 300 kPa 峰值强度/kPa 251.0 241.3 208.1 198.4 188.6 232.0 201.3 185.2 158.9 降幅/% — 4 17 21 25 8 20 26 37 2.2 不同细粒含量及形状下的强度指标特征
图 6为混合砂内摩擦角随组合形状参数的变化规律。纯钙质砂的内摩擦角正切值为0.74（φ=36.4°），混合砂内摩擦角正切值随细颗粒含量、组合圆度及组合宽长比增加均表现线性减小趋势。在细颗粒含量为40%时，组合宽长比和组合圆度达到最大值，混合砂的整体颗粒形状趋于规整，含石英粉混合砂内摩擦角正切值达到最小值0.57（φ=29.7°），其降幅为23%；含玻璃珠混合砂的内摩擦角正切值为0.53（φ= 27.8°），其降幅为28%，此时表现出更小的抗剪强度。原因在于，细颗粒的填充降低了土体骨架的稳定性，在粗颗粒间起到润滑的作用，降低了颗粒间的摩擦，因而内摩擦角逐渐减小；玻璃珠因颗粒形状圆滑，随着细颗粒含量的增加，混合砂的组合形状参数逐渐接近圆形颗粒的形状参数，颗粒之间运动产生的逐渐摩擦力较小，而石英粉的颗粒形状较不规则，在颗粒间的相对运动中，不利于滚动，需要较多的力进行滑动，故含玻璃珠混合砂的内摩擦角呈现更大的降幅。
![]()
图 6 不同细颗粒含量及形状下的内摩擦角Figure 6. Angles of internal friction at different contents and shapesof fine particles钙质砂为无黏性土，但由于具有不规则的颗粒形状，颗粒之间的嵌固咬合作用产生了不可忽视的表观黏聚力[25]。图 7为混合砂表观黏聚力随颗粒含量及组合形状参数的变化规律。试验发现，随着细颗粒含量的增加，混合砂表观黏聚力由纯钙质砂时的28 kPa开始逐渐降低，含石英粉混合砂的表观黏聚力降幅较小，最终降至17.2 kPa，而含玻璃珠的混合砂表观黏聚力呈现出较大的下降趋势，最终仅为1.2 kPa。其原因是，随着细颗粒含量的增加，钙质砂颗粒所形成的骨架逐渐被细颗粒填充，减少了钙质砂颗粒之间的接触，进而降低了钙质砂颗粒之间的嵌固咬合作用；当细颗粒含量为40%时，土体骨架由钙质砂及细颗粒共同承担，由于外形更接近于球形的玻璃珠细颗粒并不具有咬合能力，导致含玻璃珠混合砂的表观黏聚力大幅下降。
![]()
图 7 不同细颗粒含量及形状下的表观黏聚力Figure 7. Apparent cohesions at different contents and shapes of fine particles2.3 不同细粒含量及形状下的残余强度特征
为定量化表述钙质混合砂的残余强度特性，采用软化系数IB分析混合砂在剪切试验后的应变软化程度。IB越大，说明砂样的强度在峰值后衰减程度越大，软化越明显，试样的变形稳定性越低。其公式为[26]
IB=σdp−σdrσdp (3) 式中：σdp为剪切破坏时的峰值强度（kPa）；σdr为剪切破坏后残余强度（kPa）。
图 8为不同细颗粒含量和形状条件下混合砂的软化系数。随着轴向压力的增大，各组混合砂的软化系数均呈现先减小后增加的趋势，轴向压力为200 kPa时，各组试样软化系数最小，试样相对稳定。细颗粒含量为10%时软化系数达到最大值，混合砂的强度在峰值后变化最大，软化程度明显。细颗粒含量超过10%后，混合砂的软化系数随着细颗粒含量的增加而减小，且在细颗粒含量为40%时软化系数达到最小值。含石英粉混合砂的软化系数明显小于含玻璃珠混合砂，表明含石英粉混合砂的变形稳定性较高，剪切应力在峰值强度后的衰减较少。混合砂中玻璃珠含量为10%时，软化最为明显，软化系数达到最大值，其数值均在0.14以上，其原因主要是由于细颗粒含量的增加，充填在钙质砂颗粒所组成骨架的孔隙中，而玻璃珠因其颗粒形状较为圆滑，在颗粒运动中起到了润滑作用，使得含玻璃珠的钙质砂混合砂受到较小的力时就能发生颗粒运动，因此试样的变形稳定性较低，在强度达到峰值后有着较大的衰减程度。与张晨阳等[4]研究结果类似，在细粒含量为10%时稳定性最差，后续随着细颗粒含量的不断增加，钙质砂的软化系数逐渐降低，试样稳定性逐渐提高。
![]()
图 8 不同轴向压力条件下试样的脆性指标Figure 8. Brittleness indexes of samples at different axial pressures2.4 不同细粒含量及形状下的变形特征
钙质混合砂在环剪过程中的轴向位移变化可在一定程度上反映出其剪胀特征。图 9为轴向压力100 kPa时环剪试验的轴向位移变形结果。从图 9中可以看出，纯钙质砂在剪切开始时出现小幅度剪缩，随着剪切位移不断增加，出现明显的剪胀现象，最大剪胀量为-0.438 mm。由图 9（a）可知，当含石英粉的混合砂中细颗粒含量分别为10%，20%，30%，40%时，最大剪缩量为0.181，0.203，0.26，0.536 mm。表明随着细颗粒含量的增加，混合砂的变形特征逐渐由剪胀状态转变为剪缩状态。由图 9（b）可知，当含玻璃珠混合砂中细颗粒含量为10%，20%，30%，40%时，最大剪缩量为分别为0.261，0.331，0.421，0.604 mm。表明同一围压下随着细颗粒含量的增加，试样的剪胀性逐渐减小，剪缩状态逐渐明显，当细颗粒含量为10%~30%时，两种不同细颗粒的混合砂的轴向变形都较为接近，相差幅度较小，当细颗粒含量增加至40%时，其轴向变形的增长趋势明显增加，且含玻璃珠的混合砂的剪缩量比含石英粉的混合砂具有明显的增长趋势，且剪缩特征更加明显。
为研究不同轴向压力条件下混合砂的轴向变形情况，对细颗粒含量为20%和40%的混合砂的轴向变形进行汇总分析，如图 10所示。由图 10（a）可知，细颗粒含量为20%的混合砂在轴向压力为100 kPa时，轴向位移先逐渐增大，随后发生小幅度的减小，最后趋于稳定，两种混合砂的最大剪缩量分别为0.203，0.262 mm；当轴向压力为200 kPa时，混合砂均呈现出明显的剪缩特征，剪缩量随剪切位移增加而增加，两种混合砂的最大剪缩量分别为0.41，0.662 mm；轴向压力为300 kPa时，混合砂均呈现出较大的剪缩特征，随着剪切位移的增加，两种混合砂的最大剪缩量分别为0.723，0.834 mm。如图 10（b）所示，细颗粒含量为40%时混合砂均表现出明显的剪缩特征。随着剪切位移的增加，剪缩量呈现出不断上升的趋势，且随着轴向压力的增加，两种混合砂的最大剪缩量分别达到0.536，0.604，0.718，0.843，0.974，1.181 mm。表明随着轴向压力不断增加，混合砂的最大剪缩量逐渐增加，呈现出明显的剪缩特征，是因为随着轴向压力不断增加，试样在剪切过程中发生较大程度的颗粒破碎现象，大颗粒逐渐减少，致使混合砂的剪缩现象显著。
![]()
图 10 细颗粒含量为20%和40%的钙质混合砂的轴向变形−剪切位移关系曲线Figure 10. Axial deformation-shear displacement curves for calcareous mixed sand with fine particles contents of 20% and 40%3. 颗粒破碎特征
钙质混合砂特殊的剪切特征主要是由于其颗粒特征和颗粒破碎等多重因素共同作用的结果，为了对不同细颗粒含量及形状条件下混合砂的剪切过程进行细观分析，需建立细颗粒含量及形状与颗粒破碎之间的联系。对剪切后试样进行颗粒级配试验，得到剪切前后各组试样的颗粒级配曲线，采用Hardin[27]所提出的相对破碎指标Br进行颗粒破碎的定量化评价，相对颗粒破碎率的计算原理如图 11所示，其计算公式如下：
Br=BtBp0。 (4) 式中：Bp0为初始破碎势，即试样试验前级配曲线的积分；Bt为总破碎势，即试验前后的整体破碎势之差。
图 12为不同细颗粒含量及形状条件下钙质混合砂的破碎关系图。由图 12可知，各组混合砂的相对破碎率均小于3%。在轴向压力为100，200，300 kPa条件下，试样的最大破碎率分别为0.9%，1.3%，2.4%。表明随轴向压力的不断增加，更多的钙质砂颗粒发生破碎，混合砂的颗粒破碎度呈现出不断上升的趋势。在细颗粒含量0~40%的增长过程中，当轴向应力为100 kPa时，两种混合砂的相对破碎率由0.9%分别减小至0.41%和0.29%；当轴向应力为200 kPa时，两种混合砂的相对破碎率由1.3%分别减小至0.87%和0.78%；当轴向应力为300 kPa时，两种混合砂的相对破碎率由2.4%分别减小至1.88%和1.78%，该现象是因为细颗粒含量的增加降低了钙质砂之间的咬合程度，颗粒破碎率也随之减小。与含玻璃珠混合砂的相对破碎率相比，含石英粉混合砂的相对破碎率更高，二者差值在0.07%~0.16%，其主要原因是石英粉的颗粒形状与玻璃珠相比更为不规则，棱角度较高，对于钙质砂颗粒之间的咬合降低程度较小，致使含石英粉的钙质混合砂有着较高的颗粒破碎率。
![]()
图 12 不同细颗粒含量及形状条件下钙质混合砂的破碎关系图Figure 12. Crushing relationship of calcareous mixed sand under different contents and shapes of fine particles图 13为钙质混合砂的相对破碎率随组合形状参数变化趋势图。由图 13可知，混合砂的相对破碎率随着细颗粒含量的增加不断减小，与组合圆度及组合宽长比均表现出线性负相关关系，表明当混合砂整体的组合形状参数越偏离球形颗粒，颗粒轮廓呈现枝状时，其破碎程度越高。细颗粒含量的增加使得混合砂整体的组合形状参数向球形颗粒靠近，包围在钙质砂周围细颗粒增加，减少了钙质砂颗粒之间的接触，起到润滑作用，使得颗粒破碎程度不断降低。玻璃珠含量的不断增加，表现出更为明显的润滑作用，且混合砂整体的组合形状参数更加趋近球形颗粒的形状参数，含玻璃珠混合砂的颗粒破碎程度明显较小。
![]()
图 13 相对破碎率与组合形状参数关系图Figure 13. Relationship between relative crushing rate and composite shape parameters4. 结论
通过对不同细颗粒形状及含量条件下的钙质混合砂进行环剪试验，分析细颗粒形状及含量间的差异对钙质混合砂力学特性及破碎特性的影响，得到以下5点结论。
（1）钙质混合砂的应力-位移曲线均呈现软化特性，含圆形玻璃珠混合砂的峰值强度明显小于含角形石英粉混合砂的峰值强度，随着细颗粒含量不断增加，混合砂的峰值强度逐渐降低，且含玻璃珠混合砂峰值强度有着更大的降低趋势。
（2）钙质混合砂的强度指标随细颗粒含量、组合圆度、组合宽长比增加均呈现减小趋势。由于玻璃珠表面圆滑，致使含玻璃珠混合砂表现出更低的内摩擦角和表观黏聚力。
（3）钙质混合砂的软化系数随轴向压力的增加先减小后增大，在200 kPa时达到最小值，并随细颗粒含量增加表现出先增大后减小的变化趋势，当细颗粒含量为10%时，软化系数最大，当细颗粒含量超过10%后，混合砂的软化系数逐渐减小，且含石英粉试样的软化系数明显低于含玻璃珠试样。
（4）纯钙质砂表现出明显的剪胀，且随细颗粒含量和轴向压力增加，混合砂的剪缩量逐渐增大，当细颗粒含量达到40%时呈现出明显的剪缩，含玻璃珠混合砂的剪缩增量更大。
（5）在不同轴压下，钙质混合砂的相对破碎率与组合圆度及组合宽长比均呈线性负相关，且轴压的增加会造成相对破碎率的增加。随细颗粒含量的增加，混合砂组合形状参数趋于球形形状参数，颗粒破碎程度减小，含玻璃珠混合砂具有更小的相对破碎率。
-
![]()
图 5 一维固结模拟结果与理论对比
Figure 5. Comparison between simulated results and theoretical solution of 1D consolidation tests
![]()
图 6 一维热传导模拟与理论结果对比
Figure 6. Comparison between simulated results and analytical solutions of 1D heat conduction tests
![]()
图 8 一维降压开采试验和T+H模拟结果
Figure 8. Comparison of 1D dissociation test results obtained from experiments and T+H simulations
![]()
图 10 一维降压下开采产气量、变形、胶结破坏、压力演化
Figure 10. Evolution of gas production, vertical strain, percentage of bond breakage and pressure for 1D depressurization tests
![]()
图 11 胶结破坏与变形关系曲线
Figure 11. Relationship between strain and percentage of bond breakage during hydrate dissociation
![]()
图 12 一维升温开采产气量、变形、胶结破坏演化
Figure 12. Evolution of gas production, vertical strain and percentage of bond breakage from 1D heat stimulation tests
![]()
图 13 一维升温开采中分解前沿面与t1/2关系
Figure 13. Evolution of dissociation front location against t1/2 from 1D heat stimulation tests
表 1 水合物分解模型参数
Table 1 Parameters for MH dissociation
模型网格 网格尺寸/mm 初始压力/MPa 温度/℃ 盐度/% 开采压力/MPa 土体参数 Δx Δy Δz 渗透系数/m2 土体密度/(kg·m-3) 土体比热容/(J·kg-1·K-1) 土体湿导热系数/(W·m-1·℃-1) 1×1×10 45 45 15 10.5 12.5 0 4 1.5×10-12 2650 1000 3.3 
下载: 导出CSV
-
[1] MILKOV A V. Global estimates of hydrate-bound gas in marine sediments: how much is really out there?[J]. Earth-Science Reviews, 2004, 66(3/4): 183-197.
[2] WAITE W F, SANTAMARINA J C, CORTES D D, et al. Physical properties of hydrate-bearing sediments[J]. Reviews of Geophysics, 2009, 47(4): 465-484.
[3] SOGA K, NG M Y A, LEE S L, et al. Characterisation and Engineering Properties of Methane Hydrate Soils[M]. London: Taylor and Francis, 2006: 2591-2642.
[4] SONG Y C, YANG L, ZHAO J F, et al. The status of natural gas hydrate research in China: A review[J]. Renewable and Sustainable Energy Reviews, 2014, 31: 778-791. doi: 10.1016/j.rser.2013.12.025
[5] YAN C, REN X, CHENG Y, et al. Geomechanical issues in the exploitation of natural gas hydrate[J]. Gondwana Research, 2020, 81: 403-422. doi: 10.1016/j.gr.2019.11.014
[6] MORIDIS G J, COLLETT T S, DALLIMORE S R, et al. Numerical studies of gas production from several CH4 hydrate zones at the Mallik site, Mackenzie Delta, Canada[J]. Journal of Petroleum Science and Engineering, 2004, 43(3/4): 219-238.
[7] TANG L G, LI X S, FENG Z P, et al. Control mechanisms for gas hydrate production by depressurization in different scale hydrate reservoirs[J]. Energy and Fuels, 2007, 21(1): 227-233. doi: 10.1021/ef0601869
[8] MORIDIS G J, REAGAN M T. Estimating the upper limit of gas production from Class 2 hydrate accumulations in the permafrost: 2. Alternative well designs and sensitivity analysis[J]. Journal of Petroleum Science and Engineering, 2011, 76(3/4): 124-137.
[9] OYAMA H, KONNO Y, MASUDA Y, et al. Dependence of depressurization-induced dissociation of methane hydrate bearing laboratory cores on heat transfer[J]. Energy & Fuels, 2009, 23(10): 4995-5002.
[10] YIN Z Y, MORIDIS G, CHONG Z R, et al. Numerical analysis of experiments on thermally induced dissociation of methane hydrates in porous media[J]. Industrial & Engineering Chemistry Research, 2018, 57(17): 5776-5791.
[11] MIYAZAKI K, MASUI A, SAKAMOTO Y, et al. Triaxial compressive properties of artificial methane-hydrate-bearing sediment[J]. Journal of Geophysical Research, 2011, 116(B6): B06102.
[12] HYODO M, YONEDA J, YOSHIMOTO N, et al. Mechanical and dissociation properties of methane hydrate-bearing sand in deep seabed[J]. Soils and Foundations, 2013, 53(2): 299-314. doi: 10.1016/j.sandf.2013.02.010
[13] HYODO M, LI Y H, YONEDA J, et al. Effects of dissociation on the shear strength and deformation behavior of methane hydrate-bearing sediments[J]. Marine and Petroleum Geology, 2014, 51(2): 52-62.
[14] SONG Y C, ZHU Y M, LIU W G, et al. Experimental research on the mechanical properties of methane hydrate-bearing sediments during hydrate dissociation[J]. Marine and Petroleum Geology, 2014, 51(2): 70-78.
[15] KIMOTO S, OKA F, FUSHITA T, et al. A chemo-thermo- mechanically coupled numerical simulation of the subsurface ground deformations due to methane hydrate dissociation[J]. Computers and Geotechnics, 2007, 34(4): 216-228. doi: 10.1016/j.compgeo.2007.02.006
[16] RUTQVIST J, MORIDIS G J. Numerical studies on the geomechanical stability of hydrate-bearing sediments[J]. SPE Journal, 2009, 14(2): 267-282. doi: 10.2118/126129-PA
[17] UCHIDA S. Numerical Investigation of Geomechanical Behaviour of Hydrate Bearing Sediments[D]. Cambridge: University of Cambridge, 2012.
[18] GUPTA S, WOHLMUTH B, HELMIG R. Multi-rate time stepping schemes for hydro-geomechanical model for subsurface methane hydrate reservoirs[J]. Advances in Water Resources, 2016, 91: 78-87. doi: 10.1016/j.advwatres.2016.02.013
[19] ZHOU M L, SOGA K, YAMAMOTO K, et al. Geomechanical responses during depressurization of hydrate- bearing sediment formation over a long methane gas production period[J]. Geomechanics for Energy and the Environment, 2020, 23: 100111. doi: 10.1016/j.gete.2018.12.002
[20] 蒋明镜, 付昌, 贺洁, 等. 不同开采方法下深海能源土的离散元模拟[J]. 岩土力学, 2015, 36(增刊2): 639-647. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S2093.htm JIANG Ming-jing, FU Chang, HE Jie, et al. Distinct element simulations of exploitation of methane hydrate bearing sediments with different methods[J]. Rock and Soil Mechanics, 2015, 36(S2): 639-647. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2015S2093.htm
[21] CUNDALL P A, STRACK O D L. A discrete numerical model for granular assemblies[J]. Géotechnique, 1979, 29(1): 47-65. doi: 10.1680/geot.1979.29.1.47
[22] JIANG M J, SUN R H, ARROYO M, et al. Salinity effects on the mechanical behaviour of methane hydrate bearing sediments: a DEM investigation[J]. Computers and Geotechnics, 2021, 133: 104067. doi: 10.1016/j.compgeo.2021.104067
[23] MORIDIS G. User's manual for the hydrate v1. 5 option of TOUGH+ v1. 5: A code for the simulation of system behavior in hydrate-bearing geologic media[M]. Oak, Ridge US Lawrence Berkeley National Laboratory, 2014.
[24] 蒋明镜. 现代土力学研究的新视野——宏微观土力学[J]. 岩土工程学报, 2019, 41(2): 195-254. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201902002.htm JIANG Ming-jing. New paradigm for modern soil mechanics: Geomechanics from micro to macro[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(2): 195-254. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201902002.htm
[25] ITASCA PFC 6.0 (Particle Flow Code) Documentation[M]. Minmeapolics: ITASCA Inc, 2019.
[26] KIM H C, BISHNOI P R, HEIDEMANN R A, et al. Kinetics of methane hydrate decomposition[J]. Chemical Engineering Science, 1987, 42(7): 1645-1653. doi: 10.1016/0009-2509(87)80169-0
[27] LIU X, FLEMINGS P B P. Dynamic multiphase flow model of hydrate formation in marine sediments[J]. Journal of Geophysical Research: Solid Earth, 2007, 112(B3): B03101.
[28] ZHANG A, JIANG M J, THORNTON C. A coupled CFD-DEM method with moving mesh for simulating undrained triaxial tests on granular soils[J]. Granular Matter, 2020, 22(1): 1-13. doi: 10.1007/s10035-019-0969-4
[29] HU L T, WINTERFELD P H, FAKCHAROENPHOL P, et al. A novel fully-coupled flow and geomechanics model in enhanced geothermal reservoirs[J]. Journal of Petroleum Science and Engineering, 2013, 107: 1-11. doi: 10.1016/j.petrol.2013.04.005
[30] JAEGER J C, COOK NGW, ZIMMERMAN R W, Fundamentals of Rock Mechanics[M]. 4th eds. Malden: Blackwell Publishing, 2007.
[31] WANG Y, KOU X, FENG J C, et al. Sediment deformation and strain evaluation during methane hydrate dissociation in a novel experimental apparatus[J]. Applied Energy, 2020, 262: 114397. doi: 10.1016/j.apenergy.2019.114397
[32] 蒋明镜, 张望城. 一种考虑流体状态方程的土体CFD-DEM耦合数值方法[J]. 岩土工程学报, 2014, 36(5): 793-801. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201405002.htm JIANG Ming-jing, ZHANG Wang-cheng. Coupled CFD-DEM method for soils incorporating equation of state for liquid[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(5): 793-801. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201405002.htm
[33] CHEN F, DRUMM E C, GUIOCHON G. Coupled discrete element and finite volume solution of two classical soil mechanics problems[J]. Computers and Geotechnics, 2011, 38(5): 638-647. doi: 10.1016/j.compgeo.2011.03.009
[34] JIANG M J, KONRAD J M, LEROUEIL S. An efficient technique for generating homogeneous specimens for DEM studies[J]. Computers and Geotechnics, 2003, 30(7): 579-597. doi: 10.1016/S0266-352X(03)00064-8
[35] JIANG M J, SHEN Z F, WANG J F. A novel three- dimensional contact model for granulates incorporating rolling and twisting resistances[J]. Computers and Geotechnics, 2015, 65: 147-163. doi: 10.1016/j.compgeo.2014.12.011
[36] WEIBULL W. The phenomenon of rupture in solids[J]. Proceedings of Royal Swedish Institute of Engineering Research, 1939, 153: 1-55.
-
期刊类型引用(7)
1. 刘志，杨阳. 黄原胶联合微生物加固改善混合砂蓄水特性试验研究. 土木与环境工程学报(中英文). 2025(03): 49-57 . 
百度学术
2. 梁越，冉裕星，许彬，张鑫强，何慧汝. 细颗粒含量影响渗流侵蚀规律的细观机理研究. 岩土工程学报. 2025(05): 1099-1106 . 本站查看
3. 高霞，梅雯艳，吴强，崔祥龙. 球度和扁平度对含瓦斯水合物煤体宏细观力学性质影响研究. 中国安全生产科学技术. 2025(05): 86-94 . 百度学术
4. 蔡国庆，刁显锋，杨芮，王北辰，高帅，刘韬. 基于CFD-DEM的流-固耦合数值建模方法研究进展. 哈尔滨工业大学学报. 2024(01): 17-32 . 百度学术
5. 程建毅，郭晓军，李泳. 泥石流物源土体标度分布参数与粘聚力的关系. 山地学报. 2024(03): 401-410 . 百度学术
6. 孙增春，刘汉龙，肖杨. 砂-粉混合料的分数阶塑性本构模型. 岩土工程学报. 2024(08): 1596-1604 . 本站查看
7. 孙建强. 颗粒形状对黄土干密度与抗剪强度的影响. 黑龙江工程学院学报. 2024(06): 8-15+23 . 百度学术
其他类型引用(5)
计量
- 文章访问数: 367
- HTML全文浏览量: 23
- PDF下载量: 301
- 被引次数: 12
