石油与天然气化工  2022, Vol. 51 Issue (2): 53-58
气主导盲管内水合物沉积特性实验研究
郭朋皞1 , 宋光春1 , 李玉星1 , 王武昌1 , 胡皓晨1 , Amadeu K. Sum2     
1. 山东省油气储运安全省级重点实验室·中国石油大学(华东);
2. 科罗拉多矿业学院
摘要目的 为了研究气主导盲管内水合物的沉积特性,使用专门设计的气主导盲管实验装置开展了一系列水合物管壁沉积实验。方法 实验中,竖直管内的气体在密度差的作用下产生自然对流,进而携带高温水蒸气在竖直管内运移并不断在低温的竖直管壁上冷凝,冷凝水最终与气体反应形成管壁水合物沉积层。结果 实验中竖直管内存在水合物沉积层生长最快、厚度最大且最容易形成堵塞的危险位置,管壁水合物沉积层存在厚度和表观体积减小的现象。根据水合物在形成过程中传质和传热条件的不同,竖直管内管壁水合物沉积层从上到下存在6种不同形态。结论 实验过程中,管壁水合物沉积层的生长速率会逐渐降低,这主要是由管壁水合物沉积层具有的“隔热”作用导致的。该研究结果可为深水流动安全保障盲管水合物防治技术的发展提供理论支撑。
关键词水合物    气主导    盲管    沉积    
Experimental investigation on the characteristics of hydrate deposition in gas-dominated deadlegs
Guo Penghao1 , Song Guangchun1 , Li Yuxing1 , Wang Wuchang1 , Hu Haochen1 , Amadeu K. Sum2     
1. Shandong Key Laboratory of Oil-Gas Storage and Transportation Safety, China University of Petroleum, Qingdao, Shandong, China;
2. Colorado School of Mines, Golden, Colorado, USA
Abstract: Objective In order to investigate the characteristics of hydrate deposition in gas-dominated deadlegs, a series of experiments on hydrate deposition were conducted in a specially designed gas-dominated deadleg setup. Methods During the experiments, natural convection happened on and the gas in the deadleg under the action of density difference, further induced the transportation of the high-temperature vapour inside the deadleg and the condensation of the vapour on the cold pipe wall. As a result, the condensed water eventually reacted with the gas and formed hydrate deposits on the pipe wall. Results In experiments, there existed a dangerous position in the deadleg where the hydrate deposits grew fastest, had the largest thickness and was most likely to form blockage. Moreover, the thickness and apparent volume of the hydrate deposits gradually decreased in the late stage of the experiments. According to the different mass and heat transfer conditions during the formation of hydrates, there existed six different hydrate deposit morphologies in the deadleg from top to bottom. Conclusions During the experiment, the growth rate of the hydrate deposits gradually decreased, which was mainly caused by the "insulation effect" of the hydrate deposits. The conclusions of this work can provide a guidance for the development of hydrate management strategy in deepwater deadlegs.
Key words: hydrate    gas-dominated    deadleg    deposition    

在深水油气资源的开发和输送过程中,天然气水合物极易在管道内形成并引发堵塞,阻碍正常生产[1-6]。近年来,越来越多的现场案例、理论及实验研究表明,沉积是水合物在油气管道内生成后导致管道堵塞的关键原因[7-9]。因此,有必要对管道内水合物的沉积特性展开深入研究[10-11]。盲管作为深水管输系统的重要组成部分[12],其内部的水合物沉积问题同样不容忽视。为了系统研究气主导盲管内水合物的沉积特性,美国科罗拉多矿业学院联合挪威国家石油公司设计制造了一系列气主导盲管实验装置[12],并开展了大量实验研究:首先,设计开发了盲管内水合物沉积层的检测方法[13];随后,系统研究了盲管内管壁温度、储水器内水温、盲管尺寸、盲管内水蒸气含量、复杂几何结构、物理约束及自然对流强度等因素对盲管内水合物沉积特性的影响[14-22];此外,G. C. Song等[23]还基于实验数据计算出了盲管内形成的水合物沉积层及沉积层内纯水合物的导热系数。

基于上述研究,本文将重点研究气主导盲管内水合物的沉积机理,水合物沉积层的形成过程、生长形态及生长速率。本研究结果可为深水流动安全保障盲管水合物防治技术的发展提供理论支撑。

1 实验设计
1.1 实验装置及材料

采用如图 1所示的实验装置开展盲管内的水合物沉积实验,装置的主体由一段竖直管及与竖直管相连通的储水器组成。其中,竖直管内径(d)2.4 cm,长28 cm。竖直管外部紧密缠绕有多圈铜制冷却盘管,冷却盘管与恒温水浴相连,用于控制竖直管的管壁温度。如图 1所示,竖直管缠绕有冷却盘管的部分被等分为五小节,从使用每节到第五小节底部的距离x与管径之比(x/d)来区分每小节的具体位置。与竖直管相连通的储水器底部放置有加热器,用于控制储水器内部的水温。同时,储水器中安装有磁力搅拌器,其搅拌作用可使储水器内的水均匀受热,并加快水分向竖直管内挥发。实验中,降低竖直管壁温度的同时提高储水器内的水温,水受热挥发产生的水蒸气会逐渐在过冷管壁上冷凝成液态水,并最终与管内气体反应生成管壁水合物沉积层。该实验装置的详细设计参数见参考文献[22]。

图 1     气主导盲管装置示意图

使用甲烷和乙烷的混合气及去离子水生成水合物,其中混合气中甲烷摩尔分数为74.2%,乙烷摩尔分数为25.8%。

1.2 实验步骤及工况

本实验主要包括管壁水合物沉积层的形成及管壁水合物沉积层的泄压加热分解两个主要步骤。其中,在分解水合物沉积层之前,需先将竖直管壁温度降至-10 ℃后泄压,以便在自保护效应的作用下(常压低温下水合物分解异常缓慢)拆卸竖直管,并使用管道内窥镜原位测量竖直管内壁水合物沉积层的厚度、分布和形貌[14, 24]

本实验涉及的主要变量为管壁温度、储水器水温及反应时间等。在相同的温压工况下,在不同反应时间进行实验,以此得到水合物沉积层随时间的演化情况。具体实验工况如表 1所列。

表 1    实验工况

2 实验结果及分析
2.1 沉积层的形成机理

以储水器内水温为60 ℃、管壁温度为4 ℃的工况(简称60 ℃/4 ℃工况)为例,对实验中气主导管壁水合物沉积层的形成机理进行说明。

图 2所示,实验开始前,整个实验系统内的温度稳定维持在25 ℃,无温度梯度。实验开始后,管壁温度降低为4 ℃,储水器内水温上升至60 ℃,实验系统内开始形成较大的温度梯度,并诱发气体产生密度梯度(经PVTSim软件计算可知,对于本实验所用的混合气,在压力为10 MPa、温度分别为4 ℃、30 ℃、60 ℃时,其密度分别为132.9 kg/m3、102.9 kg/m3和84.4 kg/m3)。接着,密度梯度会引发气体产生自然对流并携带水蒸气在盲管内运移。在这一过程中,高温水蒸气会在过冷的竖直管壁上冷凝形成小液滴,并逐渐铺展成液膜,储水器内的水位也会因此逐渐降低,如图 2所示。随后,管壁上的冷凝水会与混合气反应生成水合物沉积层并逐渐增厚,最终覆盖整个竖直管壁,甚至引发竖直管堵塞。

图 2     水合物沉积层形成机理示意图

2.2 沉积层生长过程及形态

以工况1~工况12为例,对实验中管壁水合物沉积层的生长过程及形态进行说明。图 3展示了工况1~工况12中管壁水合物沉积层的形态和厚度分布,图中白色部分即为管壁上形成的水合物沉积层。图 4则为工况1~工况12中管壁水合物沉积层沿管段的厚度分布曲线,图中使用短虚线标记了各工况下管壁水合物沉积层厚度最大的位置,蓝色虚线代表 3组重复实验,对比3种工况可知,本实验具有良好的可重复性。

图 3     工况1~工况12管壁水合物沉积层的形态和厚度分布

图 4     工况1~工况12管壁水合物沉积层沿管段的厚度分布曲线

图 3图 4可知,实验开始初期,管壁上形成的水合物沉积层较薄,且其厚度沿管段的分布较为均匀。随后,已形成的管壁水合物沉积层因生长而逐渐增厚,同时,其厚度沿管段的分布也开始逐渐产生差异,变得不再均匀。其中,以60 ℃/4 ℃工况为例,由图 4(a)可知,管段x/d为7.3~8.5的位置处,管壁水合物沉积层生长速率最快,沉积层厚度最大,最容易形成管段水合物堵塞。因此,称x/d为7.3~8.5的位置为60 ℃/4 ℃工况下的危险位置。由图 4可知,当竖直管壁温度发生变化时,管段的危险位置也会发生变化,其中60 ℃/8 ℃和60 ℃/12 ℃工况下,危险位置的x/d分别为7.3和7.7~8.5。

图 4发现:靠近竖直管底部形成的水合物沉积层,其厚度会随时间的增加而逐渐减小,如红色虚线框部分所示。具体原因分析如下:实验中形成的水合物沉积层具有疏松多孔的结构,未转化的冷凝水就填充在这些孔穴中;当某些孔穴内未转化冷凝水的含量达到一定程度后,孔穴就会在重力作用下坍塌,造成管壁水合物沉积层对应位置处厚度的减小。

图 4还可以看出,实验中管壁水合物沉积层厚度逐渐减小的现象主要发生在竖直管底部,其原因是竖直管底部温度和水蒸气体积含量较高,水蒸气冷凝速率快,未转化冷凝水在孔穴内的积累速度也较快。因此,水合物沉积层更容易发生坍塌并导致厚度减小;而竖直管上部,温度和水蒸气体积含量均较低,冷凝水不易在孔穴内累积,沉积层也就不易因坍塌而导致厚度减小。

60 ℃/4 ℃工况中管壁水合物沉积层厚度减小的现象,相对于60 ℃/8 ℃工况和60 ℃/12 ℃工况发生的时间更晚,且发生的位置更低,这主要是因为60 ℃/4 ℃工况中管壁温度更低,而管壁水合物沉积层在更大过冷度下的生长速率和冷凝水消耗速率更快,这就大大降低了沉积层内未转化冷凝水的积累速率,进而在一定程度上抑制了沉积层厚度减小的现象。

图 3图 5可知,竖直管内形成的管壁水合物沉积层基本沿竖直管中轴线呈对称分布和生长。此外,由图 3还看出,从竖直管顶部到竖直管底部,管壁水合物沉积层的形态逐渐由粗糙变为平滑(竖直管最下端沉积层的粗糙形态是由冷凝水结冰导致的)。根据实验观察,在管壁水合物沉积层形态由粗糙到平滑的转变过程中,管壁水合物沉积层形态还可以大致细分为6种,如图 6所示:最上端的针型结构、上端接近最大厚度处的鳞片状结构、最大厚度处的平滑结构、最大厚度处下部的有孔结构、下端接近最大厚度处的平滑结构、最下端的绒毛结构。

图 5     从竖直管顶部观察到的管壁水合物沉积层形态

图 6     管壁水合物沉积层形态由粗糙到平滑的转变

2.3 沉积层生长速率

以工况1~工况4为例,对实验过程中管壁水合物沉积层生长速率的变化进行说明。图 7为工况1~工况4实验中气体消耗量及水合物沉积层体积的变化曲线。由图 7可知,实验过程中管壁水合物沉积层的生长速率逐渐降低,这一规律同样可以从管壁水合物沉积层的体积变化曲线中得出。造成管壁水合物沉积层生长速率逐渐降低的原因,主要是系统压力的逐渐降低和沉积层厚度的逐渐增加。

图 7     工况1~工况4气体消耗量及水合物沉积层体积变化曲线

本实验为定容实验,故系统压力会随着实验的进行而逐渐降低,即水合物生成的主要驱动力过冷度会逐渐降低。因此,管壁水合物沉积层的生长速率会逐渐降低。根据实验过程中系统压力的变化曲线可知,各工况实验过程中由水合物沉积层形成和生长所导致的系统压降超过0.6 MPa。这种情况下,由PVTSim软件计算出的实验过程中的最大过冷度降仅为0.4 ℃。由此可以看出,实验过程中系统压力或过冷度降低对管壁水合物沉积层生长速率的影响很小,几乎可以忽略不计。

水合物沉积层逐渐增厚对沉积层生长速率降低的影响分析如下:类似于管壁上形成的蜡沉积层[25],管壁上形成的水合物沉积层同样具有“隔热”作用[26],会大大削弱管壁和管道内部间的热传递。因此,管壁水合物沉积层的生长速率会随着沉积层的形成和增厚而逐渐降低。当管壁水合物沉积层的厚度生长到一定程度时,在水合物沉积层的“隔热”作用下,管壁水合物沉积层内表面(即生长前沿)的温度会逐渐由管壁温度升高至实验压力下对应的水合物相平衡温度,此时,管壁水合物沉积层就会最终停止生长。

3 结论

利用气主导盲管实验装置开展了一系列水合物管壁沉积实验,获得的主要结论如下:

(1) 竖直管内的气体会在密度差的作用下产生自然对流,竖直管内的高温水蒸气会随着气体的自然对流在竖直管内运移,并不断在低温的竖直管壁上冷凝,竖直管壁上的冷凝水最终会与气体反应形成管壁水合物沉积层。

(2) 竖直管内存在水合物沉积层生长最快、厚度最大且最容易形成堵塞的危险位置。管壁水合物沉积层存在厚度和表观体积减小的现象。根据水合物在形成过程中传质和传热条件的不同,竖直管内管壁水合物沉积层从上到下存在6种不同形态。

(3) 管壁水合物沉积层的生长速率会逐渐降低,这主要是由管壁水合物沉积层具有的“隔热”作用导致的。

参考文献
[1]
SLOAN E D. Clathrate Hydrates of Natural Gases[M]. New York: Marcel Dekker, 2007: 1-5.
[2]
宋光春, 李玉星, 王武昌, 等. 油气混输管道中天然气水合物的形成和堵塞过程研究[J]. 石油与天然气化工, 2017, 46(2): 38-43. DOI:10.3969/j.issn.1007-3426.2017.02.008
[3]
陈兵, 张曙旋, 郭焕焕. 含杂质CO2管道输送水合物形成规律研究[J]. 石油与天然气化工, 2021, 50(2): 42-47. DOI:10.3969/j.issn.1007-3426.2021.02.007
[4]
杨健, 冯莹莹, 张本健, 等. 超高压含硫气井井筒内天然气水合物解堵技术[J]. 天然气工业, 2020, 40(9): 64-69. DOI:10.3787/j.issn.1000-0976.2020.09.008
[5]
宫敬, 史博会, 吕晓方. 多相混输管道水合物生成及其浆液输送[J]. 中国石油大学学报(自然科学版), 2013, 37(5): 163-167. DOI:10.3969/j.issn.1673-5005.2013.05.023
[6]
宫敬, 史博会, 陈玉川, 等. 含天然气水合物的海底多相管输及其堵塞风险管控[J]. 天然气工业, 2020, 40(12): 133-142. DOI:10.3787/j.issn.1000-0976.2020.12.015
[7]
李长俊, 黄婷, 贾文龙. 深水天然气水合物及其管道输送技术[J]. 科学通报, 2016, 61(22): 2449-2462.
[8]
ZHAO J F, WANG B, SUM A K. Dynamics of hydrate formation and deposition under pseudo multiphase flow[J]. AIChE Journal, 2017, 63(9): 4136-4146. DOI:10.1002/aic.15722
[9]
ZHOU S D, YAN H Y, SU D, et al. Investigation on the kinetics of carbon dioxide hydrate formation using flow loop testing[J]. Journal of Natural Gas Science and Engineering, 2018, 49: 385-392. DOI:10.1016/j.jngse.2017.10.022
[10]
张保勇, 于洋, 靳凯, 等. 松散沉积物中降压幅度和饱和度对天然气水合物分解过程的影响[J]. 天然气工业, 2020, 40(8): 133-140. DOI:10.3787/j.issn.1000-0976.2020.08.010
[11]
李彦龙, 刘昌岭, 廖华林, 等. 泥质粉砂沉积物-天然气水合物混合体系的力学特性[J]. 天然气工业, 2020, 40(8): 159-168. DOI:10.3787/j.issn.1000-0976.2020.08.013
[12]
KINNARI K J, ASKVIK K M, LI X Y, et al. Hydrate management of deadlegs in oil and gas production systems: background and development of experimental systems[J]. Energy & Fuels, 2017, 31: 11783-11792.
[13]
SA J H, LEE B R, ZHANG X W, et al. Hydrate management in deadlegs: detection of hydrate deposition using permittivity probe[J]. Energy & Fuels, 2018, 32(2): 1693-1702.
[14]
SA J H, LEE B R, ZHANG X W, et al. Hydrate management in deadlegs: hydrate deposition characterization in a 1-in. vertical pipe system[J]. Energy & Fuels, 2017, 31(12): 13536-13544.
[15]
SA J H, ZHANG X W, SUM A K. Hydrate management in deadlegs: effect of driving force on hydrate deposition[J]. Fuel, 2020, 279: 118481. DOI:10.1016/j.fuel.2020.118481
[16]
ZHANG X W, LEE B R, SA J H, et al. Hydrate management in deadlegs: effect of wall temperature on hydrate deposition[J]. Energy & Fuels, 2018, 32: 3254-3262.
[17]
ZHANG X W, LEE B R, SA J H, et al. Hydrate management in deadlegs: effect of header temperature on hydrate deposition[J]. Energy & Fuels, 2017, 31: 11802-11810.
[18]
ZHANG X W, LEE B R, SA J H, et al. Hydrate management in deadlegs: effect of pipe size on hydrate deposition[J]. Energy & Fuels, 2020, 34: 1422-1431.
[19]
ZHANG X W, SA J H, SUM A K. Hydrate management in deadlegs: effect of water vapor content on hydrate deposition[J]. Fuel, 2020, 273: 117714. DOI:10.1016/j.fuel.2020.117714
[20]
SA J H, ZHANG X W, SUM A K. Hydrate management in deadlegs: hydrate deposition in pipes with complex geometry[J]. Fuel, 2020, 269: 117440. DOI:10.1016/j.fuel.2020.117440
[21]
SA J H, ZHANG X W, LI X Y, et al. Hydrate management in deadlegs: limiting hydrate deposition with physical restriction[J]. Fuel, 2020, 270: 117506. DOI:10.1016/j.fuel.2020.117506
[22]
SONG G C, LI Y X, SUM A K. Hydrate management in deadlegs: effect of natural convection on hydrate deposition[J]. Energy & Fuels, 2020, 34: 14094-14100.
[23]
SONG G C, Li Y X, SUM A K. Hydrate management in deadlegs: thermal conductivity of hydrate deposits[J]. Energy & Fuels, 2021, 35: 3112-3118.
[24]
宋光春. 管道内水合物颗粒聚集和膜生长沉积特性研究[D]. 青岛: 中国石油大学(华东), 2019.
[25]
SINGH P, VENKATESAN R, FOGLER S H. Formation and aging of incipient thin film wax-oil gels[J]. AIChE Journal, 2000, 46: 1059-1074. DOI:10.1002/aic.690460517
[26]
LINGELEM M N, MAJEED A I, STANGE E. Industrial experience in evaluation of hydrate formation, inhibition, and dissociation in pipeline design and operation[J]. Annals of the New York Academy of Science, 1994, 715: 75-93. DOI:10.1111/j.1749-6632.1994.tb38825.x