赵建国,孙友宏,王海亮

(吉林大学建设工程学院,吉林长春 130026)

天然气水合物孔底冷冻取样器室内冷冻试验能量计算与分析

赵建国,孙友宏,王海亮

(吉林大学建设工程学院,吉林长春 130026)

针对以干冰作为冷冻剂、酒精作为载冷剂的天然气水合物孔底冷冻取样方法,提出干冰与酒精用量的计算方法。计算出在干冰双层保冷的方式下,干冰与酒精的用量。通过实验与计算机模拟分析,验证该计算方法的可行性。针对试验过程中出现的问题,提出对试验装置的改进,为取样器的设计提供参考。

天然气水合物;干冰;冷量计算;冷冻试验;计算机模拟

0 引言

FCS型天然气水合物孔底冷冻取样器,是利用外部冷源 (干冰)在孔底降低水合物岩心温度,采用主动式降温的方法降低水合物的临界分解压力,达到被动式降压来抑制水合物的分解,以此来获得高保真的水合物样品。

取样器冷源——干冰的保冷采用双层保冷的方法,保冷层内层为硬质聚氨酯泡沫塑料,外层为聚氨酯泡沫塑料。按照允许最大冷损失量的原则,设计内层保冷层的厚度为 15 mm,外层保冷层厚度为 3 mm。通过对比试验,可知采用聚氨酯泡沫塑料与硬质聚氨酯塑料的双层保冷方法可以大幅度提高干冰在井下的保存时间。但该保冷层厚度是在最大允许冷损失量的前提下设计的,干冰必定会有所损耗。

在孔内,只有岩心冷冻到一定的温度,才能实现被动式降压的效果。因此,在考虑最大冷损失量的前提下,必须对干冰的冷量进行进一步的分析计算。

利用干冰的低温特性 (-7815℃)和升华时强大的制冷能力 (636 kJ/kg),为试样的冷冻提供冷源。利用酒精的大的导热系数〔01185 W/(m·K)〕和低的凝固点 (-114℃)的特性,作为载冷剂。

如图 1所示,钻进过程中,酒精腔体 2下部阀体关闭。回次结束后,投球,冲洗液改变流向,将酒精腔体 2下部阀体打开,推动活塞 1,将酒精压入到干冰腔体 4中。由于酒精的催化作用,干冰 6迅速升华吸热,酒精被迅速降到极低的温度。随着反应的进行,被降温后的酒精沿着酒精通道 8进入到冷冻腔体 9中,通过铜质岩心管与岩心发生冷交换,最终将岩心冷冻到一定的温度。

2 干冰冷冻能量计算

2.1 理想状态干冰能量计算

根据能量守恒原理,设计的干冰量能够提供的冷量为将岩心、岩心管以及酒精冷冻一定温度的冷量。

图 1 取样器冷冻方式原理示意图

根据这一原理,首先不考虑保冷结构表面换热系数,忽略最大允许干冰冷损失量,计算出在理想状态下所需干冰用量。

干冰升华吸热总量等于将岩心、岩心管以及酒精冷却到一定温度所需冷量之和,即:

式中:V岩心、V岩心管、V冷冻腔——分别为岩心、岩心管和冷冻腔的体积;C岩心、C岩心管、C酒精——分别为岩心、岩心管和酒精的体积比热;qv——干冰体积升华热,按密度ρ=600 kg/m3计算,取 3181×108J/m3;K——岩心冷冻温度差。

干冰物理性质为:升华点 -7815℃,密度 600 kg/m3,升华热 6135×105J/kg,体积升华热 3181×108J/m3。

实验设计取心直径为 0105 m,取心长度 012 m,岩心管外径 0106 m,冷冻腔内径为 01074 m,冷冻温差 K取 30℃,岩心、岩心管和酒精物理性质以及冷冻腔体积见表 1。

表 1 岩心、岩心管和酒精物理性质以及冷冻腔体体积

将以上数据带入式 (3),经计算得所需干冰体积为 V=0117×10-3m3,也就是 0110 kg。

2.2 最大允许冷损失量计算

最大允许冷损失量是指以每平方米保冷层外表面积为单位的冷损失量,其单位为 W/m2。在实际井下钻进过程中,冷损失的主要原因为流动的冲洗液与保冷结构外表面存在着冷交换,因此,计算冷损失量首先需要确定保冷结构外表面换热系数α,α的取值按照下式计算:

式中:v——保冷结构表面流体流速,对于金刚石钻进,取最大上返流速 016 m/s。

最大允许冷损失量为:

经计算,允许最大冷损失量为 7617 W/m2。

2.3 冷冻过程中损失量计算

冷冻腔内部为 3 mm厚的硬质聚氨酯泡沫绝热层对低温酒精进行保冷,但仍会有部分冷量的散失,需要对这部分散失的冷量进行计算。根据《工业设备及管道绝热工程设计规范》(GB/T 50264-97),圆筒形单层绝热结构冷损失量应按下式计算:

式中:Q——以每平方米绝热层外表面积表示的热损失量,W/m2;T0-Tα——温度差 K,取 50 ℃;D1——保冷层外径 ,取 01079 m;D0——保冷层内径,取 01073 m;λ——硬质聚氨酯导热系数,取01025 W/(m·K);α——表面传热系数,按式 (4)计算 ,v取 0 m/s。

将已知量代入式(6)中可算得冷冻腔体内冷损失量为 143 W/m2。

2.4 考虑冷损失量后干冰量计算

以 600 m深钻孔为例,考虑下钻时间 1 h,钻进时间 015 h,提钻时间也为 1 h,则冷冻介质干冰需保存时间为 115 h。设计所需干冰冷量为最大允许损失冷量、理论计算需干冰冷量之与冷冻过程中冷损失量之和,即:

式中:D0——保冷层内径 ,01049 m;D1——干冰腔保冷层外径,01079 m;D2——冷冻腔保冷层外径;t1——干冰保存时间,1×3600 s;t1——提钻时间,1×3600 s;Q1——干冰腔保冷层外表面单位面积冷损失量,7617 W/m2;Q2——冷冻腔保冷层外表面单位面积冷损失量,143 W/m2;H2——冷冻腔高度;qv——干冰体积升华热 ,3181 ×108J/m3;V理论——理论干冰计算量,0117×10-3m3。

将以上数据代入式 (12)可算出干冰腔体高度H1=0154 m。

五是，致歉者对自己的错误行为说出一些补救措施，以减少对双方的不良影响。例如，“我打烂了你妈妈的盆栽，真是对不住，我会赔偿的。”

实际所需要干冰体积为:

将 D0与 H代入 (11)可以算出 V实际=110×10-3m3,即 016 kg。

3 室内冷冻试验

3.1 实验装置结构原理

室内冷冻试验采用图 2所示冷冻装置。整个试验装置主要由酒精腔、干冰腔上接手、干冰腔、干冰腔下接手、冷冻腔组成。干冰被放置在岩心管上方的双层保温干冰腔内。实验开始后,向干冰腔倒入酒精。酒精与干冰热交换后,成为低温液体,沿着接手通道流入冷冻腔中,即岩心管和衬管环状间隙。然后通过岩心管向内部模拟岩心传递冷量,最后完成对岩心的冷冻。

图 2 冷冻试验装置结构图

天然气水合物主要赋存于泥质砂岩、细砂岩、粉砂岩、石灰岩等沉积层中,以细脉状充填于沉积物或岩石的裂隙中。天然气水合物与冰有明显的相似性,其性质比较见表 2。因此,可以通过冷冻含冰的沉积层样品来模拟实际天然气水合物地层。试验所用冷冻样品为砂土样。

表 2 天然气水合物和冰的性质比较

经过前面冷冻能量计算可知,要使模拟岩心达到 50℃的温降,需要干冰 016 kg,酒精 114 L。试验开始之前首先测量岩心样品的初始温度为 1615℃,酒精的初始温度为 1118℃。试验中每隔 2 min记录一次温度数据,待样品温度恢复至零度时停止记录数据。

3.2 实验结果与分析

图 3 岩样和酒精温度变化曲线图

图 4 冷冻岩样实物照片

由图 3分析可知,酒精在 36 min的时候温度降至最低为 -2617℃左右,因而确定酒精与干冰充分反应的时间为 36 min。而模拟岩心在这段时间内降温速度最快,由 1615℃降到 -413℃,平均温降约为 016℃/min。之后温降逐渐趋于缓慢,大约 60 min的时候降到最低的 -14182℃,然后温度开始上升。岩心的温度降幅为 3018℃。

同时可以看到,岩心的降温曲线出现轻微的波动,这是由于实际实验过程中,酒精的注入速率会影响干冰升华的剧烈程度,从而导致酒精温度变化的波动,而酒精温度的变化直接影响到岩心冷冻的速率。

模拟岩心在 60 min以后升温过程中,温度要低于载冷剂酒精 3～4℃,这是因为载冷剂酒精与模拟岩心温度达到平衡之后,停止向岩心传递冷量。但由于受外界环境的影响,使低温酒精通过冷冻腔外壁大量散失冷量,而此时模拟岩心由于周围低温酒精的保冷作用,冷量的散失要相对滞后,所以载冷剂酒精的温度在 60 min后高于土样温度。

4 冷冻试验有限元模拟分析

4.1 模型的建立

首先对于模拟岩心的传冷是轴对称的,径向的温度梯度远远大于轴向的温度梯度,所以我们将实际的三维岩心简化为二维的模型;其次铜质的岩心管导热率与模拟岩心和酒精相比要大得多,其对实验结果的影响非常微小,所以在模拟过程中忽略岩心管对传热的影响。

4.2 降温过程模拟分析

计算结果可用温度分布云图和具体温度值显示,我们分析所得冻结 60 min之后的温度梯度图见图 5,温度变化曲线图见图 6。

图 5 降温过程温度分布图

图 6 降温过程温度变化曲线图

由图 6分析可知,岩心的温度在前 20 min降低很快,降温幅度达到了 22℃,后 40 min温度降低较慢,降温幅度为 10℃。岩心中心点在冷冻至第 12 min时温度降至 0℃,岩样计时结束时的温度为 -141633℃。

4.3 升温过程的模拟分析

计算结果可用温度分布云图和具体温度值显示,我们模拟得温度升到 0℃的温度梯度图见图 7,温度变化曲线图见图 8。

由图 8分析可知,岩心升温过程的温度曲线较平缓,接近于直线,说明在升温过程中岩心向外散热比较均匀,且岩心周围的低温酒精对岩心也有一定的保冷作用。

5 结论与改进

(1)通过模拟分析与试验结果对比可知,模拟岩心体积中心点温度变化基本一致。说明建立的二维岩心模型可以真实地反映岩心冷冻温度的变化。从最后获得的冷冻样品可以看出,采用该冷冻方法可以成功地获取冷冻岩心。

(2)实验结果与理论计算结果存在很大误差,理论计算的岩心温降为 50℃,而实验过程中模拟岩心温度只降到 -1418℃,其温度降与理论计算相差近 20℃。经过分析可能存在以下几点原因:

图 7 升温过程温度分布图

图 8 升温过程温度分布图

①酒精注入干冰腔后,与干冰反应。与此同时,一部分酒精还未来得及与干冰充分反应便通过干冰腔下接手的酒精通道进入干冰腔内部,等到干冰充分反应升华之后,冷冻腔内已充满酒精。与干冰充分反应的酒精则留在干冰腔体内部,无法进入冷冻腔。

②由于受到尺寸的限制,酒精通道直径只有 10 mm,大大限制了干冰腔体与冷冻腔体内酒精的冷传递速度,导致没有持续的冷源补给冷冻腔体内酒精,这样也就无法对岩心持续冷冻。

所以,在今后的试验中,需要对实验装置进一步改进,首先在酒精进入冷冻腔之前,使其能够与干冰充分反应,其次在干冰腔与冷冻腔之间增加类似于热管的传热装置,加快干冰腔与冷冻腔之间的热传导,保证对冷冻腔冷源的持续补给。

[1]郭威,孙友宏.天然气水合物孔底冷冻取样方法的室内试验研究[J].探矿工程 (岩土钻掘工程),2009,36(5).

[2] 郭威,孙友宏,等.干冰升华式孔底冷冻取样器冷源保冷方法的实验研究[J].探矿工程 (岩土钻掘工程),2009,36(9).

[3]王海亮.FCS-108型天然气水合物空地冷冻取样器的研制及实验研究[D].吉林大学,2009.

[4] GB/T 50264-97,工业设备及管道绝热工程设计规范[S].

Energy Calculation and Analysis on Laboratory Freezing Exper iment of Freezing Sampler for Gas Hydrate

Z HAO Jian-guo,SUN You-hong,WANG Hai-liang(College of Construction Engineering,Jilin University,Changchun Jilin 130026,China)

According to freezing sampling method of gas hydrate with dry ice as refrigerant and alcohol as cooling agent,calculation was given for dry ice and alcohol consumption,and the consumption of dry ice and alcoholwith the double-layer insulation were calculated out.The feasibility of the calculationmethodwas verified through experiments and computer sim- ulations.Considering the problems in the exper iment,some improvementwas suggested to the experimental device to pro-vide a reference to sampler design.

gas hydrate;dry ice;cold calculation;freezing experiment;computer simulation

1 取样器冷冻方式原理

P634.5

A

1672-7428(2010)01-0008-05

2009-11-16

科技部国际合作与交流专项“天然气水合物孔底冷冻取样方法及取样器的研究”(项目编号:2007DFR60100);中国地质调查局项目“青藏高原冻土带天然气水合物调查评价”(项目编号:1212010818055)

赵建国 (1983-),男 (汉族),河南孟州人,吉林大学硕士研究生,地质工程专业,从事天然气水合物勘探技术的研究,吉林省长春市西民主大街 6号,,jgzhao07@mails.jlu.edu.cn。