管道流动体系下天然气水合物生成模型的研究进展_周诗岽

基金项目：

国家自然科学基金（编号：５１１７６０１５）、中国石油科技创新基金项目（编号：２０１１Ｄ５００６０６０６）、江苏省研究生培养创新工程项目（编号：ＣＸＺＺ１３＿０７３５）、扬中市科技局资助项目（编号：ＹＥ２０１２０１５

）。　　作者简介：

周诗岽，１９７８年生，副教授；主要从事气体水合物方面的研究工作。地址：（２１３０１６）江苏省常州市钟楼区白云路常州大学石油工程学院。电话：（０５１９）８３２９０２８０。Ｅ－ｍａｉｌ：ｚｓｄ＠ｃｃｚｕ．ｅｄｕ．ｃｎ管道流动体系下天然气水合物生成模型的研究进展

周诗岽１，２　余益松２　王树立２　张国忠１　李泓

３１．

中国石油大学（华东）储运与建筑工程学院　２．常州大学石油工程学院油气储运技术省重点实验室３．

长江（扬中）电脱盐设备有限公司　　周诗岽等．管道流动体系下天然气水合物生成模型的研究进展．天然气工业，２０１４，３４（２）：９２－

９８．　　摘　要　管道流动体系下天然气水合物生成模型的建立对天然气水合物浆液的输送、

管输天然气水合物防治以及天然气水合物技术的应用都具有重要意义。为此，查阅了大量的国内外相关文献并进行了总结与分析，认识到目前对该类模型的研究较少，现有的模型也是在静态釜式反应器天然气水合物生成理论基础上拓展而来的，主要包括驱动力、成核速率、诱导时间、水合物生长等方面的模型，上述模型被广大研究者用于计算管道单个截面处天然气水合物的生成速率预测，具有较好的计算精度。但现有模型用于管道流动体系下天然气水合物生成特性的预测还不成熟，需要进一步开展管道流动体系下天然气水合物的生成机理、管道沿线温度变化、添加剂及其流动界面对气液传质的影响等研究，建立动力学、传热、传质三者相结合的管道流动体系下天然气水合物生成模型，以此来解决管道流动体系下天然气水合物生成预测的技术难题。

关键词　管道　天然气水合物　驱动力　诱导时间　模型　成核　流动体系　生成

ＤＯＩ：１０．３７８７／ｊ

．ｉｓｓｎ．１０００－０９７６．２０１４．０２．０１５Ｒｅｓｅａｒｃｈ　ｐｒｏｇ

ｒｅｓｓ　ｉｎ　ｔｈｅ　ｎａｔｕｒａｌ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｍｏｄｅｌ　ｕｎｄｅｒ　ｔｈｅ　ｐｉｐｅ　ｆｌｏｗ　ｓｙｓｔｅｍＺｈｏｕ　Ｓｈｉｄｏｎｇ１，２，Ｙｕ　Ｙｉｓｏｎｇ２，Ｗａｎｇ　Ｓｈｕｌｉ　２，Ｚｈａｎｇ　

Ｇｕｏｚｈｏｎｇ１，Ｌｉ　Ｈｏｎｇ３（１．Ｃｏｌｌｅｇｅ　ｏｆ　Ｐｉｐｅｌｉｎｅ　＆Ｃｉｖｉｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｃｈｉｎａ　Ｕｎｉｖｅｒｓｉｔｙ　

ｏｆ　Ｐｅｔｒｏｌｅｕｍ，Ｑｉｎｇｄａｏ，Ｓｈａｎｄｏｎｇ２６６５８０，Ｃｈｉｎａ；２．Ｊｉａｎｇｓｕ　Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　

ｏｆ　Ｏｉｌ　＆Ｇａｓ　Ｓｔｏｒａｇｅ　ａｎｄ　Ｔｒａｎｓｐｏｒｔａｔｉｏｎ　Ｔｅｃｈｎｏｌｏｇｙ∥Ｃｈａｎｇｚｈｏｕ　Ｕｎｉ－ｖｅｒｓｉｔｙ，Ｃｈａｎｇｚｈｏｕ，Ｊｉａｎｇｓｕ　２１３０１６，Ｃｈｉｎａ；３．Ｃｈａｎｇｊｉａｎｇ＜Ｙａｎｇｚｈｏｎｇ＞Ｅｌｅｃｔｒｉｃ　Ｄｅｓａｌｔｉｎｇ　Ｅｑｕｉｐ

ｍｅｎｔＣｏ．，Ｌｔｄ．，Ｙａｎｇｚｈｏｎｇ，Ｊｉａｎｇ

ｓｕ２１２２００，Ｃｈｉｎａ）ＮＡＴＵＲ．ＧＡＳ　ＩＮＤ．ＶＯＬＵＭＥ　３４，ＩＳＳＵＥ　２，ｐｐ．

９２－９８，２／２５／２０１４．（ＩＳＳＮ　１０００－０９７６；Ｉｎ　Ｃｈｉｎｅｓｅ）Ａｂｓｔｒａｃｔ：Ｔｈｅ　ｅｓｔａｂｌｉｓｈｍｅｎｔ　ｏｆ　ｔｈｅ　ｎａｔｕｒａｌ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｍｏｄｅｌ　ｕｎｄｅｒ　ｔｈｅ　ｓｙｓｔｅｍ　ｏｆ　ｐｉｐｅ　ｆｌｏｗ　ｉｓ　ｏｆ　ｓｉｇｎｉｆｉｃａｎｃｅ　ｔｏ　ｈｙ

ｄｒａｔｅｓｌｕｒｒｙ　

ｔｒａｎｓｐｏｒｔａｔｉｏｎ，ｐｉｐｅｌｉｎｅ　ｇａｓ　ｈｙｄｒａｔｅ　ｐｒｅｖｅｎｔｉｏｎ　ａｎｄ　ｃｏｎｔｒｏｌ，ａｎｄ　ｔｈｅ　ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ｈｙｄｒａｔｅ－ｒｅｌａｔｅｄ　ｔｅｃｈｎｏｌｏｇｉｅｓ．Ｉｎ　ｖｉｅｗ　ｏｆｔｈｉｓ，ａｃｃｏｒｄｉｎｇ　

ｔｏ　ａ　ｇｒｅａｔ　ｎｕｍｂｅｒ　ｏｆ　ｄｏｍｅｓｔｉｃ　ａｎｄ　ｆｏｒｅｉｇｎ　ｌｉｔｅｒａｔｕｒｅｓ，ｆｅｗ　ｓｔｕｄｉｅｓ　ｈａｖｅ　ｂｅｅｎ　ｆｏｕｎｄ　ａｂｏｕｔ　ｔｈｅ　ｓａｉｄ　ｍｏｄｅｌ　ａｎｄ　ｔｈｅ　ｅｘｉｓｔ－ｉｎｇ　

ｍｏｄｅｌｓ　ａｒｅ　ｍｏｓｔｌｙ　ｂａｓｅｄ　ｕｐｏｎ　ｔｈｅ　ｔｈｅｏｒｙ　ｏｆ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ｔｈｅ　ｓｔａｔｉｃ　ｒｅａｃｔｏｒ，ｓｕｃｈ　ａｓ　ｄｒｉｖｉｎｇ　ｆｏｒｃｅ　ｍｏｄｅｌ，ｎｕｃｌｅａｔｉｏｎ　ｒａｔｅｍｏｄｅｌ，ｉｎｄｕｃｔｉｏｎ　ｔｉｍｅ　ｍｏｄｅｌ　ａｎｄ　ｈｙｄｒａｔｅ　ｇｒｏｗｔｈ　ｍｏｄｅｌ，ｅｔｃ．Ｍａｎｙ　

ｓｃｈｏｌａｒｓ　ｕｓｅ　ｔｈｅｓｅ　ｍｏｄｅｌｓ　ｔｏ　ｃａｌｃｕｌａｔｅ　ｔｈｅ　ｒａｔｅ　ｏｆ　ｈｙｄｒａｔｅ　ｆｏｒｍａ－ｔｉｏｎ　ａｔ　ａ　ｓｉｎｇｌｅ　ｐｉｐｅ　ｃｒｏｓｓ　ｓｅｃｔｉｏｎ　ｗｉｔｈ　ａ　ｈｉｇｈ　ｃａｌｃｕｌａｔｉｏｎ　ａｃｃｕｒａｃｙ．Ｈｏｗｅｖｅｒ，ｔｈｅ　ｅｘｉｓｔｉｎｇ　ｍｏｄｅｌｓ　ｆｏｒ　ｐｒｅｄｉｃｔｉｎｇ　

ｔｈｅ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｕｎｄｅｒ　ｔｈｅ　ｓｙｓｔｅｍ　ｏｆ　ｐｉｐｅ　ｆｌｏｗ　ａｒｅ　ｎｏｔ　ｍａｔｕｒｅ　ｅｎｏｕｇｈ　ｓｏ　ｔｈａｔ　ｆｕｒｔｈｅｒ　ｓｔｕｄｉｅｓ　ａｒｅ　ｓｔｉｌｌ　ｎｅｅｄｅｄ　ｏｎ　ｔｈｅ　ｈｙ

ｄｒａｔｅ　ｆｏｒｍａ－ｔｉｏｎ　ｍｅｃｈａｎｉｓｍ　ｕｎｄｅｒ　ｔｈｅ　ｓｙｓｔｅｍ　ｏｆ　ｐｉｐｅ　ｆｌｏｗ　ａｎｄ　ｔｈｅ　ｉｍｐａｃｔｓ　ｏｆ　ｅｎｖｉｒｏｎｍｅｎｔａｌ　ｔｅｍｐｅｒａｔｕｒｅ　ｃｈａｎｇ

ｅｓ，ａｄｄｉｔｉｖｅｓ　ａｎｄ　ｔｈｅ　ｆｌｏｗ　ｉｎｔｅｒｆａｃｅｏｎ　ｔｈｅ　ｇａｓ－ｌｉｑｕｉｄ　ｍａｓｓ　ｔｒａｎｓｆｅｒ　ｉｎ　ｌｉｎｅ　ｐｉｐｅｓ，ｅｔｃ．Ｃｏｎｓｉｄｅｒｉｎｇ　

ｔｈｅ　ｄｙｎａｍｉｃｓ，ｈｅａｔ　ｔｒａｎｓｆｅｒ，ｍａｓｓ　ｔｒａｎｓｆｅｒ　ａｎｄ　ｏｔｈｅｒ　ｆａｃｔｏｒｓ，ｔｈｅ　ｈｙ－ｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｍｏｄｅｌ　ｕｎｄｅｒ　ｔｈｅ　ｓｙｓｔｅｍ　ｏｆ　ｐｉｐｅ　ｆｌｏｗ　ｉｓ　ｔｈｅ　ｋｅｙ　ｔｏ　ｔｈｅ　ｐｒｏｂｌｅｍ　ｏｆ　ｐｒｅｄｉｃｔｉｎｇ　

ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ｓｕｃｈ　ａ　ｃａｓｅ．Ｋｅｙｗｏｒｄｓ：ｐｉｐｅｌｉｎｅｓ，ｎａｔｕｒａｌ　ｇａｓ　ｈｙｄｒａｔｅ，ｄｒｉｖｉｎｇ　

ｆｏｒｃｅ，ｉｎｄｕｃｔｉｏｎ　ｔｉｍｅ，ｍｏｄｅｌ，ｎｕｃｌｅａｔｉｏｎ，ｆｌｏｗ　ｓｙｓｔｅｍ，ｆｏｒｍａｔｉｏｎ·２９·　　　　　　　　　　　　天　然　气　工　业　　　　　　　　　　　　　　　　　　２

０１４年２月

天然气水合物是由主客体分子相互作用而形成的

一种非化学计量的笼形络合物［

１］

，管道流动体系下水合物的生成规律及数学模拟是开展水合物资源的开采、管输水合物的防治以及水合物利用技术等研究的

基础［

２－

４］。管道流动体系下水合物的形成过程与静态釜式水合物反应器的中水合物的生成机理大体相同，都是一种由流体相向固体相转变并伴随着热量和物质转移的过程。因此可将其看作一种泛意义上的结晶动

力学过程［

５］

。当前对水合物动力学模型的研究大多是在釜式水合物反应器的基础上得到的，针对管道流动体系下水合物生成动力学模型的相关研究较少，然而，由于二者水合物生成过程大体相同，其相关模型可以相互借鉴。水合物的生成过程包括成核与生长２个阶

段［６］

，相应的动力学模型建立分别从驱动力、成核、生

长３个方向展开。

１　驱动力模型

物质运动状态的改变以及化学反应的发生都需要

在一定的势差驱动下才能进行，水合物的生成也不例外。通常，用来描述水合物生成驱动力的参数包括温差、过冷度、过饱和度、化学势差、浓度差、吉布斯自由

能等［

７－

９］。　　早期有关水合物驱动方面的研究相对较少。１

９世纪末，Ｍａｋｏｇｏｎ［１０］和Ｓｌｏａｎ［１１］

通过实验系统研究了水合物的驱动过程，并取得了突破性的进展。然而，由于当时对水合物生成过程中的驱动力不能形成统一的认识，这方面一直未形成统一的理论，Ｋ

ａｓｈｃｈｉｅｖ与Ｆｉｒｏｏｚａｂａｄｉ［１

２］

在前人研究的基础上，分别建立了气液相化学平衡且温度恒定、气液相化学平衡且压力恒定、气液相不满足化学平衡且温度恒定、气液相不满足化学平衡且压力恒定这４种不同条件下的水合物生成驱动力模型，并给出了具体的简化条件及求解方法。这一模型与数据拟合较好，为大多数学者所接受，管道流动体系下水合物的生成驱动模型大多也是在这基础上提出的。

Ｔａｌａｇ

ｈａｔ等［１

３－１６］在Ｋａｓｈｃｈｉｅｖ和Ｆｉｒｏｏｚａｂａｄｉ的基础上对其模型进行了简化，得到了模型（１），并在管道中双组分气体水合物生成研究过程中发展了Ｃｈｒｉｓ－

ｔｉａｎｓｅｎ等［１

７］

提出的以吉布斯自由能作为驱动力表征参数的驱动力模型（２

），且通过数值方法进行了求解。　　　Δｕ＝ｋＴｌｎφｅｘｐ

（ｐｅｘｐ，Ｔ）ｐｅｘ

ｐφｅｑ

（ｐｅｑ，Ｔ）ｐｅ［

］

ｑ＋（ｎｗυｗ－υｈ）

（ｐｅｘｐ－ｐｅ

ｑ）（１）　　　ΔＧｅｘｐ

＝ｎｗ（υｗ－υｈ）

（ｐｅｑ－ｐｅｘｐ）＋ＲＴ∑ｎｏｃ

ｋ

ｎｋｌｎｆｅ

ｑｋ

ｆｅｘ

ｐ（）

ｋ（２）式中Δｕ为过饱和度，Ｊ；Δ

Ｇｅｘｐ

为吉布斯自由能的变化量；ｐｅｑ为相平衡压力，ＭＰａ；ｐｅｘ

ｐ为实验条件下的压力，

ＭＰａ；Ｒ为气体常数，ｋＪ／（ｍｏｌ·Ｋ）；Ｔ为系统温度，Ｋ；

ｋ为波尔兹曼常数；ｎｗ为水合数；υｗ为水分子的体积，

ｍ３

；υｈ为单个水合物晶胞的体积与平衡压力下所含有

的气体分子数的比值，ｍ３；ｎｋ为气体分子的数量；ｆｅｑｋ为相平衡条件下ｋ组分气体的逸度系数；ｆｅｘｐ

ｋ为实验条件下ｋ组分气体的逸度系数；φｅｑ为相平衡条件下气

体的逸度系数；φｅｘｐ

为实验条件下气体的逸度系数。　　Ｓ

ａｒｓｈａｒ［１８－

２０］根据管流的实际情况对Ｋａｓｈｃｈｉｅｖ和Ｆｉｒｏｏｚａｂａｄｉ的成核模型进行了优化，

之后通过将模型与（７３％甲烷＋２７％丙烷）所组成的混合气体在１０ｍ长管道中，温度处于２７７．１５～２７８．１５Ｋ、

压力为２～３ＭＰ

ａ的实验数据进行了比较，结果表明得到的模型能够对此种条件下管道中水合物的生成速率以及水相组成进行有效的预测。

此外，将ＳＲＫ［２１］、Ｐｅｎｇ

－Ｒｏｂｉｎｓｏｎ［２２］、ＥＲ［２３］

、ＰＴ［２４］、ＶＰＴ［２

５］

方程应用到水合物生成驱动力计算中具有较高的计算精度。

２　微观成核模型

水合物的成核是指形成超过临界尺寸的稳定水合

物晶核的过程，当溶液处于过冷状态或过饱和状态时，就可能发生成核现象。成核现象的发生具有很大的随机性，然而，在驱动力较大的区域，其成核的可预测性增强。文献中曾基于釜式水合物合成装置而报道了多种水合物模型，其中有：成簇成核模型；界面成核模型；随机水合物成核与界面成簇模型；反应动力学机理模

型；双过程水合物成核模型等［

２６］

。以这些模型为基础，国外的很多研究者对不同操作条件下管道流动体系下水合物的生成过程进行了研究，并从成核速率及诱导时间这两方面建立了相关模型。２．１　成核速率模型

Ｎａｔａｒａｊ

ａｎ等［２７

］

在水合物成核速率研究过程中，采用了一种基于结晶过程的经验公式（３）

。这种模型由于缺乏基础数据的支持，是一种纯经验的数学模型，其准确性难以保证，在随后的应用过程中效果并不理想。

Ｊ＝ｋ０（

Ｓ－１）ｎ

（３）式中Ｊ为水合物的成核速率，１／（ｍ３

·ｓ）；ｋ０、

ｎ为常数；Ｓ为过饱和比率。

·

３９·第３４卷第２期　　　　　　　　　　　　　　　　　集　输　与　加　工　　　　　　　　　　　　

Ｋ

ａｓｈｃｈｉｅｖ与Ｆｉｒｏｏｚａｂａｄｉ［２８］

通过分析晶体的成核过程，认为水合物的成核首先发生在气液接触面而不

是液相的内部，基于此，提出了适用于单组分均相成核与三维非均相成核速率模型。研究了添加剂对水合物成核速率的影响，得到了与之相对应的模型（４～６）。由于该模型与实验数据吻合较好，更加接近于水合物的真实成核过程，且方程的求解较简单，因而被广大学者所接受。即

Ｊ＝Ａ（１＋ｋｎＣａ

［］）ｅΔ

ｕ／ｋＴｅｘｐ－４ｃ３υ２ｈσ３

ｅｆ２７ｋＴΔｕ（）２

（

４）Ｊ＝Ａ（１＋ｋｎＣａ

［］）Ｓ

ｅｘｐ－４ｃ３υ２ｈσ３

ｅｆ２７ｋＴ３ｌｎ２

［］Ｓ（５）Ｊ＝Ａ１＋ｋｎＣａ

［］ｅｘｐΔ

Ｓｅ

ΔＴ（）ｋＴｅｘｐ－４ｃ３υ２ｈσ３

ｅｆ２７ｋΔ

Ｓ２ｅ

ＴΔＴ（）２（６

）　　其中

Ｓ＝φ（ｐ，Ｔ）ｐ

φ（

ｐｅ，Ｔ）ｐ［］ｅｅｘｐ

（ｎｗυｗ－υｈ）（ｐ－ｐｅ）［］

ｋＴ（７

）式中Ａ为动力学因子，１／（ｍ３

·ｓ）；ｃ为形状因子；ΔＳｅ为相变前后的熵变，Ｊ／Ｋ；ｐ为系统压力，ＭＰａ；ｐｅ为相

平衡压力，ＭＰａ；σｅｆ为比表面能的折算值，Ｊ／ｍ２

；ｋｎ为兰格缪尔吸附常数，ｍ３

；φ为逸度函数；

Ｃａ为溶液中的添加剂浓度，ｍ－３。

Ｔａｌａｇ

ｈａｔ等［２９］

采用小型循环流动管道对有、无抑制剂存在情况下单种气体组分的水合物生成情况进行了研究，将其所得实验数据与Ｋａｓｈｃｈｉｅｖ和Ｆｉｒｏ－ｏｚａｂａｄｉ所提出的模型相互印证，结果表明该模型准确性较高，适用于管道流动体系下水合物的成核过程。２．２　诱导时间模型

诱导时间是过饱和系统维持亚稳态特性的量度。

目前，有关水合物的诱导时间有２种定义：①从形成不稳定簇到超临界尺寸水合物核出现时所经历的时间；②基于溶液中由于水合物生成的体积变化与所消耗的气体分子数关系曲线而得到的从初始阶段到变化率转变处所经历的时间。

水合物生成动力学研究过程中，采用了多种装置及方法对诱导期进行了有效的测量，然而，有关水合物诱导时间的模型却较少。当前，水合物诱导时间的模型大多是在成核速率模型的基础上得到的。

Ｓｌｏａｎ和Ｆｌｅｙ

ｆｅｌ［３０］

在研究过程中发现水合物的诱导时间与客体分子的空隙占有率相关，只有当空隙占

有率为０．８１～０．８９时，

水合物生成才会表现出诱导时间的特性。Ｎａｔａｒａｊ

ａｎ［３１］

在随后推翻了这一结论，并根据水合物的成核速率模型建立了诱导时间模型，但由于该模型的建立依赖于成核速率模型，且其准确性在很大程度上取决于实验数据，因此，在实际应用过程中

很少被采用。Ｃｈｒｉｓｔｉａｎｓｅｎ等［３

２］

通过研究认为，相较于成核速率，诱导时间对生长速率更加敏感，并据此提出了一个经验模型，然而，这一模型与其他的实验数据拟合效果较差。

Ｋ

ａｓｈｃｈｉｅｖ与Ｆｉｒｏｏｚａｂａｄｉ［３３］

首次采用了水合物诱导时间的第二种定义，在水合物驱动力模型及成核速

率模型的基础上，针对水合物的成核类型而得到了不同的诱导时间（ｔｉ）模型，如式（８～１０）所示。　　瞬时成核：

ｔｉ＝

（ａｄ

Ｖｓ／ｃＮｃ）１

３　ｍＧ

（８

）　　过程成核＋体积内部成核：

ｔｉ＝（１＋３ｍ）ａｄ

ｃＧ３　

ｍ［］

Ｊ

１

１＋３　

（９

）　　过程成核＋界面成核：

ｔｉ＝（１＋３ｍ）ａｄ

ＶｓｃＡｓ

Ｇ３　

ｍ［］

Ｊ１

１＋３　

（１０

）　　当存在表面活性剂等添加剂时，其诱导时间（ｔｉ，ａ）

模型在上述模型的基础上进行了修正：

ｔｉ，ａ＝［（１＋ｋｇ

Ｃａ）３　ｍ（１＋ｋｎＣａ）］１

＋　ｔｉ（１１）式中ｔｉ为诱导时间，ｓ；ａｄ为水合物结晶的体积与溶液体积的比值；Ｖｓ为溶液的初始体积，ｍ

３

；Ｎｃ为微晶数量；ｍ为与晶核生长相关的一个常数，一般取值为０．５；

Ｇ为增长常数；Ａｓ为初始接触面面积，ｍ２

；ｋｇ为相关的吸附常数，ｍ

３

。　　Ｔａｌａｇ

ｈａｔ等［３４］

在对管道中气体水合物生成诱导时间进行研究的过程中，继续采用了Ｋａｓｈｃｈｉｅｖ和

Ｆｉｒｏｏｚａｂａｄｉ所提出的诱导时间模型，证明了该模型应用于动态体系下水合物诱导时间预测的有效性，并在此基础上得到了适用于管道流动体系下水合物诱导时间模型。

ｔｉ＝Ｋ

ｅ－Δ

ｕ／ｋＴ（１－ｅ

－Δｕ／ｋＴ）－

３ｍ１＋

３ｍｅｘｐ４ｃ３υ２ｈσ３

ｅ

ｆ２７（１＋３ｍ）ｋＴ（Δ

ｕ）［］

２

（１２

）式中Ｋ为动力学常数。

此外Ｄ

ａｖｉｅｓ等［３５］

采用多相流动装置进行了水合物的生成实验，将水合物的生成动力学与热力学结合

起来建立了相关模型，对长输油气管道中水合物的生

成堵塞情况有效地进行了预测。Ｎｉｃｈｏｌａｓ等［３

６］

遵循能量与质量的守恒关系，建立了水合物的沉积模型，在

·４９·　　　　　　　　　　　　天　然　气　工　业　　　　　　　　　　　　　　　　　　２

０１４年２月

精度允许的范围内，将其所得到的参数应用于流速为１

．８９Ｌ／ｍｉｎ与２．８３Ｌ／ｍｉｎ的单通道流动体系下进行水合物的沉积实验，并开展了一系列的水合物防堵研究。

３　宏观生长模型

当过饱和溶液中晶核的尺寸达到临界尺寸时，

体系的吉布斯自由能达到最大值，这时，体系为了维持稳定将自发地向吉布斯自由能较小的方向发展，从而进入水合物生长阶段。当前，流动状态下的水合物生长模型大多也是从静态水合物生长模型的基础上发展而来的，其生长过程大多通过气体消耗速率来衡量。　　Ｖｙ

ｓｎｉａｕｓｋａｓ和Ｂｉｓｈｎｏｉ［３７］

在对甲烷水合物与乙烷水合物生成过程进行研究的基础上，通过将水合物的生长过程与影响水合物生成的主要因素相关联，首

先提出了一个半经验模型。Ｋｕｍａｒ［３

８］

在随后的实验过程中，通过数据分析与比较，认为该模型描述水合物生长过程的精确度较低。即

ｒ＝Ａ１ａｓ

ｅｘｐγ

－ΔＥａ（）ＲＴｅｘｐ－ａΔ

Ｔ（）

ｂｐγ

（１３）式中ａｓ为气水界面的总表面积，

ｍ２；Ａ１为指前因子，ｍ３

／（ｍ２

·ｓ·Ｐａ）；ａ、ｂ、γ为常数；Δ

Ｅａ为活化能，ｋＪ／ｍｏｌ；ｐ为系统压力，

Ｐａ；ｒ为甲烷消耗速率，ｍ３

／ｓ；ΔＴ为过冷度，Ｋ。

刘陈伟等［３９］

对原油相、

水相、气相、水合物相四相管道流动体系中水合物的生成进行了模拟，在Ｖｙｓｎｉ－ａｕｓｋａｓ和Ｂｉｓｈｎｏｉ模型的基础上，针对具体的应用条件，采用了如下的模型计算水合物的形成速率：ＦＨ＝Ｋ

６Ａ２Ｅｗ０ｄ０

ＥｗＥ（）ｗ０２

３

ｅｘｐ

－ΔＥａ（

）ＲＴｅｘｐ－ａΔ

Ｔ（）

ｂｐγ

（１４

）式中ＦＨ为单位管长中气体的消耗速率，ｍ

３

／ｓ；Ｋ为综合预测常数，ｍ３／（ｍ２·ｓ·Ｐａ）；Ｅｗ０为水合物开始形成

时水的体积分数；Ｅｗ为管道中水相的体积分数；Ａ２为管道截面积，ｍ２

；ｄ０为水合物开始形成时的液滴直径，

ｍ。

Ｅｎｇ

ｌｅｚｏｓ等［４０］

在假设水合物的生长速率不随时间变化的基础上，综合考虑水合物的传质过程及反应

过程，提出了描述单组分及混合气体组分的水合物生长模型（１５）。该模型在搅拌式反应釜中与水合物生成数据吻合较好，然而，由于其生长速率不随时间变化的假设与实际不符，且流动情况下气液界面的不断更新使得其与静态釜式反应器传质过程存在着较大区别，使得该模型应用于流动体系下水合物的生成研究存在着较大的局限性。即

ｄ

ｎｄ［］ｔ＝ｋｒｋｄ

ｋｒ

＋ｋ

ｄ

Ａｐ

（ｆ－ｆｅｑ

）（１

５）式中Ａｐ为水合物颗粒壳体表面积，ｍ２

；ｆ为实验条件

下气体逸度，ＭＰａ；ｆｅｑ为三相平衡时的气体逸度，

ＭＰａ；ｋｒ为反应速率常数，ｍｏｌ／（ｍ２

·ＭＰａ·ｓ）；ｋｄ为传质系数，ｍ

ｏｌ／（ｍ２

·ＭＰａ·ｓ）；ｎ为水合物生长消耗气体分子摩尔数，ｍｏｌ；ｔ为时间，ｓ。　　孙长宇等

［４１］

在上述基础上，从管道中水合物生成

的实际情况出发，考虑液相中水合物的浓度，通过循环管路水合物实验测得的数据进行回归分析，建立了涵盖压力和体积分数的水合物生长模型。将该模型与实验数据的分析结果进行比较，发现其误差较小，精确度相对较高，但该模型的应用过于依赖其循环管路的水合物生成装置，其应用面较窄。即

ｄｎｄｔ

＝ｋｖ（１－υｈｙｄ）（ｆ－ｆｅｑ）

２

（１６

）式中ｋｖ为水合物的反应速率常数，ｍｏｌ／（ｓ·Ｐａ２

）；υｈｙ

ｄ为液相中水合物所占的体积分数。

美国科罗拉多矿业大学水合物研究中心针对水合物的生成开发出了ＯＬＧＡ模拟水合物的计算模块（ＣＭＳＨＹＫ）

，并将其成功地应用于管道多相流动过程中水合物的生长情况［４２］

。其主要模型如式（１７）。Ｋ

ｉｎｎａｒｉ等［４３］

运用该模型对管道中水合物生成分析发现，该模型中动力学控制水合物生长的因素太大，与实际过程不相符。即

－ｒｇ

ａｓ＝Ａｓｋ１ｅｘｐｋ

２

（）

ＴΔＴ（１７

）式中ｋ１为速率常数，ｋｇ／（ｍ２

·ｓ·Ｋ）；ｋ２为温度参数，Ｋ；ｒｇａｓ为气体消耗量，ｋｇ／ｓ。　　Ｓｋｏｖｂｏｒｇ等［４４］

基于水合物生成过程中的相平衡

原理，对Ｅｎｇ

ｌｅｚｏｓ等所建的模型进行了简化，建立了如下的生长模型。

ｄｎｄｔ

＝ｋＬＡｇ－ｌｃｗ０（ｘｉｎｔ－ｘｂ）

（１８）式中Ａｇ－ｌ为气液相界面面积，ｍ２；ｃｗ０为水分子初始浓

度，ｍｏｌ／ｍ３

；ｋＬ为液层质量传递系数，ｍ／ｓ；ｘｉｎｔ为实验

条件下界面气液相平衡时气体的摩尔分数；ｘｂ为实验条件下液相中水合物相与液相平衡时气体的摩尔分数。

Ｔａｌａｇ

ｈａｔ等［

１３－１６，

４５－

４６］

在前人研究的基础上通过小型循环管道中水合物的生成试验，发展了水合物抑制

剂体系下单组分气体形成水合物的宏观动力学生长模型，并通过进行一定的假设，建立了低剂量水合物抑制剂体系下水合物的生成模型，进而通过初始条件及边界条件对模型进行了求解，在此基础上针对管道中水

·

５９·第３４卷第２期　　　　　　　　　　　　　　　　　集　输　与　加　工　　　　　　　　　　　　

合物的生成规律提出了含水合物抑制剂体系下水合物的生长速率模型和气体消耗速率模型。

Ｇａ（ｔ，ｚ）＝εａｈＤｅｆｃＣｅ１＋ｋｇ

Ｃａ［ｅΔ

Ｇ（ｔ，ｚ）／ｋＴ－１

］（１

９）ｒｈ（

ｔ）＝∫

ｚ＝Ｌｐｉｐ

ｅｘｈ

ｋｂρｈ（ｔ，ｚ）Ｇａ（ｔ，ｚ）Ｊａ（ｔ，ｚ）Ｍｈ（ｔ，ｚ）

×ｔ３　

ｍａｉｎｔ（

ｚ）πｄ２ｐｉｐ

ｅ４

ｄｚ（２０）式中Ｇａ为生长函数；

ε为水合物晶体表面的黏附系数；ａｈ为水合物晶体的表面积，ｍ２

；Ｄｅｆｃ为扩散系数，ｍ２

／ｓ；Ｃｅ为气体成分的平衡溶解度，

ｍ－３；ΔＧ为吉布斯自由能的变化量；ｒｈ为气体消耗速率，ｍｏｌ／ｓ；Ｊａ为

固定成核速率，１／（ｍ３

·ｓ）；ｂ为形状系数；ｘｋｈ为水合

物相中各组分的气体含量；Ｌｐｉｐ

ｅ为管道的长度，ｍ；ａｉｎｔ为单位体积管道的气液接触面积，ｍ

２

；ｚ为衡量管道中位置的水平坐标；ｄｐｉｐ

ｅ为管道直径，ｍ。　　此外，

Ｆｉｄｅｌ－Ｄｕｆｏｕｒ等［４７］

研究了管道十二烷乳状液流动体系中水合物的结晶及堵塞情况，认为水滴周围的传质速率是气体消耗速率的主要控制因素。Ｓａｒ－

ｓｈａｒ等［４

８］

将水合物颗粒看成一个球体，在幂率增长模式的基础上，得到了水合物半径的增长模型。Ｇｒｅａｖｅｓ等［４

９］建立了水合物的生成速率模型，认为气体消耗速率与过冷度、比表面积、反应温度有关。Ｂａｌａｋｉｎ等［５

０］

基于质量守恒建立了一个计算模型，用于对管道湍流体系下水合物的生长情况进行预测，他们认为，管道流动体系下水合物的成核、扩散、生长均与流体的剪切强

度有关。仿真结果显示与数据拟合较好。Ｊａｓｓｉｍ［５１］

、Ｋｅｎｄｏｕｓｈ［５

２］

等也通过实验研究建立了管道中水合物的生长沉积模型。

可以看出，

管道中水合物生成过程的研究借鉴了静态釜式反应器水合物的生成理论。在晶体生成模型的基础上，Ｋａｓｈｃｈｉｅｖ和Ｆｉｒｏｏｚａｂａｄｉ分别从水合物的驱动力、成核和生长３方面系统地研究了管道水合物的生成过程，并建立了相关模型。这些模型由于准确性较高，被大多数研究者所采用并进行了一定程度的优化。然而，这一模型只是考虑了管道某一截面处水合物的生成情况，无法对管道的沿线生成情况进行有效预测。

４　结束语

当前，有关管道流动体系下水合物生成的模型较少，且大多是从静态水合物生成模型发展而来的。Ｋ

ａｓｈｃｈｉｅｖ和Ｆｉｒｏｏｚａｂａｄｉ通过对水合物的生成机理进行研究，分别在驱动力、成核、生长３方面取得了突破

性的进展，其所提出的水合物的驱动力模型、成核速率模型、诱导时间模型、气体消耗速率模型与管道流动体系下实验数据拟合较好，因而被广泛地应用于管道流动体系下水合物生成规律的研究。然而，这些模型的建立都是针对管道的某一截面提出的，不能对管道沿线水合物的生成情况进行有效的预测。管道流动体系下水合物的生成过程仍有许多问题尚待解决，这也将是今后水合物研究的重点，主要包括以下４个方面内容。

１

）管道流动体系下水合物的生成机理及其与釜式反应器中水合物生成过程的区别。

２）管道沿线温度的变化对水合物生成的影响。　　３）流动界面对气体的传质以及液体内部气体扩散的影响。

４）有无添加剂作用下管道沿线水合物的生成模型建立。

参　考　文　献

［１］ＭＯＲＡＤＩ　Ｓ，ＨＡＧＨＴＡＬＡＢ　Ａ，ＦＡＺＬＡｌ　Ａ．Ｐｒｅｄｉｃｔｉｏｎ　

ｏｆｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｃｏｎｄｉｔｉｏｎｓ　ｉｎ　ｔｈｅ　ｓｏｌｕｔｉｏｎｓ　ｃｏｎｔａｉｎｉｎｇ　ｅ－ｌｅｃｔｒｏｌｙｔｅ　ａｎｄ　ａｌｃｏｈｏｌ　ｉｎｈｉｂｉｔｏｒｓ　ａｎｄ　ｔｈｅｉｒ　ｍｉｘｔｕｒｅｓ　ｕｓｉｎｇ

ＵＮＩＱＵＡＣ－ＮＲＦ　ｍｏｄｅｌｓ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２０１３，３４９：６１－

６６．［２

］吴华，邹德永，何义忠，等．醇盐混合溶液中天然气水合物的形成［Ｊ］．石油与天然气化工，２０１１，４０（４）：３３２－３３６，３２３．ＷＵ　Ｈｕａ，ＺＯＵ　Ｄｅｙｏｎｇ，ＨＥ　Ｙｉｚｈｏｎｇ

，ｅｔ　ａｌ．Ｎａｔｕｒａｌ　ｇａｓｈｙｄｒａｔｅ　ｆｏｒｍｉｎｇ　ｃｏｎｄｉｔｉｏｎｓ　ｉｎ　ｒｉｃｈ　ｗａｔｅｒ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎ－ｇｉｎｅｅｒｉｎｇ　

ｏｆ　Ｏｉｌ　＆Ｇａｓ，２０１１，４０（４）：３３２－３３６，３２３．［３

］孙志高，刘成刚，周波．水合物储存气体促进技术实验研究［Ｊ］．石油与天然气化工，２０１１，４０（４）：３３７－

３３８．ＳＵＮ　Ｚｈｉｇａｏ，ＬＩＵ　Ｃｈｅｎｇｇａｎｇ，ＺＨＯＵ　Ｂｏ．Ｅｘｐｅｒｉｍｅｎｔａｌ　ｓｔｕｄｙｏｆ　ｇａｓ　ｓｔｏｒａｇｅ　ｉｎ　ｈｙ

ｄｒａｔｅｓ　ｗｉｔｈ　ｗｅｔ　ａｃｔｉｖｅ　ｃａｒｂｏｎ［Ｊ］．ＣｈｅｍｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ　ｏｆ　Ｏｉｌ　＆Ｇａｓ，２０１１，４０（４）：３３７－３３８．［４］张琦．天然气水合物的生产储运技术及现状［Ｊ］．石油与天然气化工，２０１３，４２（３）：２６１－

２６４．ＺＨＡＮＧ　Ｑｉ．Ｎａｔｕｒａｌ　ｇａｓ　ｈｙｄｒａｔｅ　ｓｔｏｒａｇｅ　ａｎｄ　ｔｒａｎｓｐｏｒｔａｔｉｏｎｔｅｃｈｎｏｌｏｇｙ　

ａｎｄ　ｉｔｓ　ｄｅｖｅｌｏｐｍｅｎｔ　ｓｔａｔｕｓ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉ－ｎｅｅｒｉｎｇ　ｏｆ　Ｏｉｌ　＆Ｇａｓ，２０１３，４２（３）：２６１－２６４．［５］ＺＡＲＥＮＥＺＨＡＤ　Ｂ，ＶＡＲＡＭＩＮＩＡＮ　Ｆ．Ａ　ｕｎｉｆｉｅｄ　ａｐｐ

ｒｏａｃｈｆｏｒ　ｄｅｓｃｒｉｐｔｉｏｎ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｋｉｎｅｔｉｃｓ　ｉｎ　ｔｈｅｐｒｅｓｅｎｃｅ　ｏｆ　ｋｉｎｅｔｉｃ　ｐｒｏｍｏｔｅｒｓ　ｉｎ　ｇａｓ　ｈｙ

ｄｒａｔｅ　ｃｏｎｖｅｒｔｅｒｓ［Ｊ］．Ｅｎｅｒｇｙ　Ｃｏｎｖｅｒｓｉｏｎ　ａｎｄ　Ｍａｎａｇｅｍｅｎｔ，２０１３，７３：１４４－１４９．［６］ＺＥＲＰＡ　Ｌ　Ｅ，ＲＡＯ　Ｉ，ＡＭＡＮ　Ｚ　Ｍ，ｅｔ　ａｌ．Ｍｕｌｔｉｐ

ｈａｓｅ　ｆｌｏｗｍｏｄｅｌｉｎｇ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅｓ　ｗｉｔｈ　ａ　ｓｉｍｐｌｅ　ｈｙｄｒｏｄｙｎａｍｉｃ　ｓｌｕｇｆｌｏｗ　ｍｏｄｅｌ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，２０１３，９９：２９８－

３０４．［７］ＹＡＮＧ　Ｍｉｎｇｊｕｎ，ＬＩＵ　Ｗｅｉｇｕｏ，ＳＯＮＧ　Ｙｏｎｇ

ｃｈｅｎ，ｅｔ　ａｌ．Ｅｆｆｅｃｔｓ　ｏｆ　ａｄｄｉｔｉｖｅ　ｍｉｘｔｕｒｅ（ＴＨＦ／ＳＤＳ）ｏｎ　ｔｈｅ　ｔｈｅｒｍｏｄｙ

－·６９·　　　　　　　　　　　　天　然　气　工　业　　　　　　　　　　　　　　　　　　２

０１４年２月

ｎａｍｉｃ　ａｎｄ　ｋｉｎｅｔｉｃ　ｐｒｏｐｅｒｔｉｅｓ　ｏｆ　ＣＯ２／Ｈ２ｈｙｄｒａｔｅ　ｉｎ　ｐｏｒｏｕｓｍｅｄｉａ［Ｊ］．Ｉｎｄｕｓｔｒｉａｌ　＆Ｅｎｇｉｎｅｅｒｉｎｇ　Ｃｈｅｍｉｓｔｒｙ　Ｒｅｓｅａｒｃｈ，２０１３，５２（１３）：４９１１－４９１８．

［８］ＭＯＴＴＡＨＥＤＩＮ　Ｍ，ＶＡＲＡＭＩＮＩＡＮ　Ｆ，ＭＡＦＡＫＨＥＲＩＫ．Ｍｏｄｅｌｉｎｇ　ｏｆ　ｍｅｔｈａｎｅ　ａｎｄ　ｅｔｈａｎｅ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｋｉ－ｎｅｔｉｃｓ　ｂａｓｅｄ　ｏｎ　ｎｏｎ－ｅｑｕｉｌｉｂｒｉｕｍ　ｔｈｅｒｍｏｄｙｎａｍｉｃｓ［Ｊ］．Ｊｏｕｒ－ｎａｌ　ｏｆ　Ｎｏｎ－Ｅｑｕｉｌｉｂｒｉｕｍ　Ｔｈｅｒｍｏｄｙａｎｍｉｃｓ，２０１１，３６（１）：３－２２．

［９］ＱＵ　Ｈ　Ｙ，ＭＵＮＫ　Ｔ，ＣＯＲＮＥＴＴ　Ｃ，ｅｔ　ａｌ．ｉｎｆｌｕｅｎｃｅ　ｏｆｔｅｍｐｅｒａｔｕｒｅ　ｏｎ　ｓｏｌｖｅｎｔ－ｍｅｄｉａｔｅｄ　ａｎｈｙｄｒａｔｅ－ｔｏ－ｈｙｄｒａｔｅｔｒａｎｓｆｏｒｍａｔｉｏｎ　ｋｉｎｅｔｉｃｓ［Ｊ］．Ｐｈａｒｍａｃｅｕｔｉｃａｌ　Ｒｅｓｅａｒｃｈ，２０１１，２８（２）：３６４－３７３．

［１０］ＭＡＫＯＧＯＮ　Ｙ　Ｆ，ＭＥＬＩＫＨＯＶ　Ｉ　Ｖ，ＫＯＺＬＯＶＳＫＡＹＡＥ　Ｄ，ｅｔ　ａｌ．Ｓｅｃｏｎｄａｒｙ　ｎｕｃｌｅａｔｉｏｎ　ｉｎ　ｔｈｅ　ｆｏｒｍａｔｉｏｎ　ｏｆ　ｍｅｔｈ－ａｎｅ　ｃｒｙｓｔａｌ　ｈｙｄｒａｔｅ［Ｊ］．Ｒｕｓｓｉａｎ　Ｊｏｕｒｎａｌ　ｏｆ　Ｐｈｙｓｉｃａｌ　Ｃｈｅｍ－ｉｓｔｒｙ　Ａ，２００７，８１（１０）：１６４５－１６４９．

［１１］ＳＬＯＡＮ　Ｅ　Ｄ．Ｃｌａｔｈｒａｔｅ　ｈｙｄｒａｔｅ　ｏｆ　ｎａｔｕｒａｌ　ｇａｓｅｓ［Ｍ］．ＮｅｗＹｏｒｋ：Ｍａｒｃｅｌ　Ｄｅｋｋｅｒ，１９９０．

［１２］ＫＡＳＨＣＨＩＥＶ　Ｄ，ＦＩＲＯＯＺＡＢＡＤＩ　Ａ．Ｄｒｉｖｉｎｇ　ｆｏｒｃｅ　ｆｏｒｃｒｙｓｔａｌｌｉｚａｔｉｏｎ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　ＣｒｙｓｔａｌＧｒｏｗｔｈ，２００２，２４１（２）：２２０－２３０．

［１３］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ．Ｅｖａｌｕａｔｉｏｎ　ｏｆ　ｖａｒｉｏｕｓ　ｔｙｐｅｓ　ｅｑｕａｔｉｏｎｓ　ｏｆｓｔａｔｅ　ｏｎ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｒａｔｅ　ｏｆ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｆｏｒ　ｂｉｎａｒｙ　ｇａｓｍｉｘｔｕｒｅｓ　ｉｎ　ｔｈｅ　ｐｒｅｓｅｎｃｅ　ｏｒ　ａｂｓｅｎｃｅ　ｏｆ　ｋｉｎｅｔｉｃ　ｈｙｄｒａｔｅ　ｉｎｈｉｂｉ－ｔｏｒｓ　ｉｎ　ａ　ｆｌｏｗ　ｍｉｎｉ－ｌｏｏｐ　ａｐｐａｒａｔｕｓ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２０１３，３４７：４５－５３．

［１４］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ．Ｅｘｐｅｒｉｍｅｎｔａｌ　ｉｎｖｅｓｔｉｇａｔｉｏｎ　ｏｆ　ｄｏｕｂｌｅｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ｔｈｅ　ｐｒｅｓｅｎｃｅ　ｏｆ　ｍｏｄｉｆｉｅｄ　ｓｔａｒｃｈａｓ　ａ　ｋｉｎｅｔｉｃ　ｉｎｈｉｂｉｔｏｒ　ｉｎ　ａ　ｆｌｏｗ　ｍｉｎｉ－ｌｏｏｐ　ａｐｐａｒａｔｕｓ［Ｊ］．Ｃａ－ｎａｄｉａｎ　Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１２，９０（２）：４２９－４３６．

［１５］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ．Ｅｎｈａｎｃｅｍｅｎｔ　ｏｆ　ｔｈｅ　ｐｅｒｆｏｒｍａｎｃｅ　ｏｆｋｉｎｅｔｉｃ　ｉｎｈｉｂｉｔｏｒｓ　ｉｎ　ｔｈｅ　ｐｒｅｓｅｎｃｅ　ｏｆ　ｐｏｌｙｅｔｈｙｌｅｎｅ　ｏｘｉｄｅ　ａｎｄｐｏｌｙｐｒｏｐｙｌｅｎｅ　ｏｘｉｄｅ　ｆｏｒ　ｂｉｎａｒｙ　ｍｉｘｔｕｒｅｓ　ｄｕｒｉｎｇ　ｇａｓ　ｈｙ－ｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ａ　ｆｌｏｗ　ｍｉｎｉ－ｌｏｏｐ　ａｐｐａｒａｔｕｓ［Ｊ］．Ｃａｎａｄｉ－ａｎ　Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１２，９０（１）：７９－８６．［１６］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ．Ｉｎｔｅｎｓｉｆｉｃａｔｉｏｎ　ｏｆ　ｔｈｅ　ｐｅｒｆｏｒｍａｎｃｅ　ｏｆｋｉｎｅｔｉｃ　ｉｎｈｉｂｉｔｏｒｓ　ｉｎ　ｔｈｅ　ｐｒｅｓｅｎｃｅ　ｏｆ　ｐｏｌｙｅｔｈｙｌｅｎｅ　ｏｘｉｄｅ　ａｎｄｐｏｌｙｐｒｏｐｙｌｅｎｅ　ｏｘｉｄｅ　ｆｏｒ　ｓｉｍｐｌｅ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ａｆｌｏｗ　ｍｉｎｉ－ｌｏｏｐ　ａｐｐａｒａｔｕｓ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２０１０，２８９（２）：１２９－１３４．

［１７］ＣＨＲＩＳＴＩＡＮＳＥＮ　Ｌ　Ｉ，ＲＡＮＴＡＮＥＮ　Ｊ　Ｔ，ＶＯＮ　ＢＯＮＳ－ＤＯＲＦＦ　Ａ　Ｋ，ｅｔ　ａｌ．Ａ　ｎｏｖｅｌ　ｍｅｔｈｏｄ　ｏｆ　ｐｒｏｄｕｃｉｎｇ　ａ　ｍｉｃｒｏ－ｃｒｙｓｔａｌｌｉｎｅ　ｂｅｔａ－ｓｉｔｏｓｔｅｒｏｌ　ｓｕｓｐｅｎｓｉｏｎ　ｉｎ　ｏｉｌ［Ｊ］．ＥｕｒｏｐｅａｎＪｏｕｒｎａｌ　ｏｆ　Ｐｈａｒｍａｃｅｕｔｉｃａｌ　Ｓｃｉｅｎｃｅｓ，２００２，１５（３）：２６１－２６９．［１８］ＳＡＲＳＨＡＲ　Ｍ，ＦＡＴＨＩＫＡＬＡＪＡＨＩ　Ｊ，ＥＳＭＡＥＩＬＺＡ－ＤＥＨ　Ｆ．Ｋｉｎｅｔｉｃ　ｏｆ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｏｆ　ｐｒｏｐａｎｅ　ａｎｄ　ｉｔｓ

ｍｉｘｔｕｒｅ　ｗｉｔｈ　ｍｅｔｈａｎｅ　ｉｎ　ａ　ｃｉｒｃｕｌａｔｉｎｇ　ｆｌｏｗ　ｒｅａｃｔｏｒ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２０１０，２９８（１）：３８－４７．

［１９］ＳＡＲＳＨＡＲ　Ｍ，ＦＡＴＨＩＫＡＬＡＪＡＨＩ　Ｊ，ＥＳＭＡＥＩＬＺＡＤＥＨＦ．Ｅｘｐｅｒｉｍｅｎｔａｌ　ａｎｄ　ｔｈｅｏｒｅｔｉｃａｌ　ｓｔｕｄｙ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａ－

ｔｉｏｎ　ｉｎ　ａ　ｈｉｇｈ－ｐｒｅｓｓｕｒｅ　ｆｌｏｗ　ｌｏｏｐ［Ｊ］．Ｃａｎａｄｉａｎ　Ｊｏｕｒｎａｌ　ｏｆＣｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１０，８８（５）：７５１－７５７．

［２０］ＳＡＲＳＨＡＲ　Ｍ，ＦＡＴＨＩＫＡＬＡＪＡＨＩ　Ｊ，ＥＳＭＡＥＩＬＺＡ－ＤＥＨ　Ｆ．Ｉｎｄｕｃｔｉｏｎ　ｔｉｍｅ　ｏｆ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ａ　ｆｌｏｗ

ｌｏｏｐ［Ｊ］．Ｔｈｅｏｒｅｔｉｃａｌ　Ｆｏｕｎｄａｔｉｏｎｓ　ｏｆ　Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒ－

ｉｎｇ，２０１０，４４（２）：２０１－２０５．

［２１］ＩＬＬＢＥＩＧＩ　Ｍ，ＦＡＺＬＡＬＩ　Ａ，ＭＯＨＡＭＭＡＤＩ　Ａ　Ｈ．Ｔｈｅｒ－ｍｏｄｙｎａｍｉｃ　ｍｏｄｅｌ　ｆｏｒ　ｔｈｅ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｅｑｕｉｌｉｂｒｉｕｍ　ｃｏｎｄｉ－ｔｉｏｎｓ　ｏｆ　ｃｌａｔｈｒａｔｅ　ｈｙｄｒａｔｅｓ　ｏｆ　ｍｅｔｈａｎｅ　ｐｌｕｓ　ｗａｔｅｒ－ｓｏｌｕｂｌｅｏｒ－ｉｎｓｏｌｕｂｌｅ　ｈｙｄｒａｔｅ　ｆｏｒｍｅｒ［Ｊ］．Ｉｎｄｕｓｔｒｉａｌ　＆ＥｎｇｉｎｅｅｒｉｎｇＣｈｅｍｉｓｔｒｙ　Ｒｅｓｅａｒｃｈ，２０１１，５０（１５）：９４３７－９４５０．

［２２］ＨＡＧＨＴＡＬＡＢ　Ａ，ＺＡＲＥ　Ｍ，ＡＨＭＡＤＩ　Ａ　Ｎ．Ｐｒｅｄｉｃｔｉｏｎｏｆ　ｈｙｄｒａｔｅ　ｅｑｕｉｌｉｂｒｉｕｍ　ｃｏｎｄｉｔｉｏｎｓ　ｕｓｉｎｇ　ｅｌｅｃｔｒｏｌｙｔｅ　ｃｕｂｉｃｓｑｕａｒｅ－ｗｅｌｌ　ｅｑｕａｔｉｏｎ　ｏｆ　ｓｔａｔｅ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２０１２，３３３（１５）：７４－８６．

［２３］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ，ＥＳＭＡＥＩＬＺＡＤＥＨ　Ｆ，ＦＡ－ＴＨＩＫＡＬＡＪＡＨＩ　Ｊ．Ｅｆｆｅｃｔ　ｏｆ　ｖａｒｉｏｕｓ　ｔｙｐｅｓ　ｏｆ　ｅｑｕａｔｉｏｎｓ　ｏｆｓｔａｔｅ　ｏｎ　ｄｒｉｖｉｎｇ　ｆｏｒｃｅ　ｃａｌｃｕｌａｔｉｏｎ　ａｎｄ　ｇａｓ　ｃｏｎｓｕｍｐｔｉｏｎ　ｐｒｅ－ｄｉｃｔｉｏｎ　ｉｎ　ｓｉｍｐｌｅ　ａｎｄ　ｄｏｕｂｌｅ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｗｉｔｈ　ｏｒｗｉｔｈｏｕｔ　ｐｒｅｓｅｎｃｅ　ｏｆ　ｋｉｎｅｔｉｃ　ｉｎｈｉｂｉｔｏｒｓ　ｉｎ　ｂａｔｃｈ　ｓｙｓｔｅｍｓ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｃｏｍｍｕｎｉｃａｔｉｏｎ，２０１０，１９７（４）：５８４－６０５．

［２４］ＡＢＡＹ　Ｈ　Ｋ，ＳＶＡＲＴＡＡＳ　Ｔ　Ｍ，ＫＥ　Ｗ．Ｅｆｆｅｃｔ　ｏｆ　ｇａｓｃｏｍｐｏｓｉｔｉｏｎ　ｏｎ　ｓｌｌ　ｈｙｄｒａｔｅ　ｇｒｏｗｔｈ　ｋｉｎｅｔｉｃｓ［Ｊ］．Ｅｎｅｒｇｙ　＆

Ｆｕｅｌｓ，２０１１，２５（４）：１３３５－１３４１．

［２５］ＣＨＡＰＯＹ　Ａ，ＡＬＳＩＶＡＢＩ　Ｉ，ＧＨＯＬＩＮＥＺＨＡＤ　Ｊ．Ｃｌａｔｈ－ｒａｔｅ　ｈｙｄｒａｔｅ　ｅｑｕｉｌｉｂｒｉａ　ｉｎ　ｌｉｇｈｔ　ｏｌｅｆｉｎｓ　ａｎｄ　ｍｉｘｅｄ　ｍｅｔｈａｎｅ－ｏｌｅｆｉｎｓ　ｓｙｓｔｅｍｓ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２０１３，３３７（１５）：１５０－１５５．

［２６］陈光进，孙长宇，马庆兰．气体水合物科学与技术［Ｍ］．北京：化学工业出版社，２００７．

ＣＨＥＮ　Ｇｕａｎｊｉｎ，ＳＵＮ　Ｃｈａｎｇｙｕ，ＭＡ　Ｑｉｎｇｌａｎ．Ｓｃｉｅｎｃｅ　ａｎｄｔｅｃｈｎｏｌｏｇｙ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅ［Ｍ］．Ｂｅｉｊｉｎｇ：Ｃｈｅｍｉｃａｌ　ＩｎｄｕｓｔｒｙＰｒｅｓｓ，２００７．

［２７］ＮＡＴＡＲＡＪＡＮ　Ｖ，ＢＩＳＨＮＯＩ　Ｐ　Ｒ．Ｌａｎｇｍｕｉｒ　ｃｏｎｓｔａｎｔ　ｃｏｍｐｕ－ｔａｔｉｏｎｓ　ｆｏｒ　ｇａｓ　ｈｙｄｒａｔｅ　ｓｙｓｔｅｍｓ［Ｊ］．Ｉｎｄｕｓｔｒｉａｌ　ａｎｄ　ＥｎｇｉｎｅｅｒｉｎｇＣｈｅｍｉｓｔｒｙ　Ｒｅｓｅａｒｃｈ，１９９５，３４（４）：１４９４－１４９８．

［２８］ＫＡＳＨＣＨＩＥＶ　Ｄ，ＦＩＲＯＯＺＡＢＡＤＩ　Ａ．Ｎｕｃｌｅａｔｉｏｎ　ｏｆ　ｇａｓｈｙｄｒａｔｅｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｃｒｙｓｔａｌ　Ｇｒｏｗｔｈ，２００２，２４３（７）：４７６－４８９．

［２９］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ．Ｅｆｆｅｃｔ　ｏｆ　ｖａｒｉｏｕｓ　ｔｙｐｅｓ　ｏｆ　ｅｑｕａｔｉｏｎｓ　ｏｆｓｔａｔｅ　ｆｏｒ　ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｓｉｍｐｌｅ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｗｉｔｈｏｒ　ｗｉｔｈｏｕｔ　ｔｈｅ　ｐｒｅｓｅｎｃｅ　ｏｆ　ｋｉｎｅｔｉｃ　ｉｎｈｉｂｉｔｏｒｓ　ｉｎ　ａ　ｆｌｏｗ　ｍｉｎ－ｉ－ｌｏｏｐ　ａｐｐａｒａｔｕｓ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，２００９，２８６（１）：３３－４２．

［３０］ＳＬＯＡＮ　Ｅ　Ｄ，ＦＬＥＹＦＥＬ　Ｆ．Ｈｙｄｒａｔｅ　ｄｉｓｓｏｃｉａｔｉｏｎ　ｅｎｔｈａｌｐｙａｎｄ　ｇｕｅｓｔ　ｓｉｚｅ［Ｊ］．Ｆｌｕｉｄ　Ｐｈａｓｅ　Ｅｑｕｉｌｉｂｒｉａ，１９９２，７６（２）：１２３－１４０．

［３１］ＮＡＴＡＲＡＪＡＮ　Ｖ，ＢＩＳＨＮＯＩ　Ｐ　Ｒ，ＫＡＬＯＧＥＲＡＫＩＳ　Ｎ．Ｉｎｄｕｃｔｉｏｎ　ｐｈｅｎｏｍｅｎａ　ｉｎ　ｇａｓ　ｈｙｄｒａｔｅ　ｎｕｃｌｅａｔｉｏｎ［Ｊ］．Ｃｈｅｍｉ－ｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，１９９４，４９（１３）：２０７５－２０８７．

·

７

９

·

第３４卷第２期　　　　　　　　　　　　　　　　　集　输　与　加　工　　　　　　　　　　　　

［３２］ＣＨＲＩＳＴＩＡＮＳＥＮ　Ｒ　Ｌ，ＢＡＮＳＡＬ　Ｖ，ＳＬＯＡＮ　ＪＲ　Ｅ　Ｄ．Ａ－ｖｏｉｄｉｎｇ　ｈｙｄｒａｔｅｓ　ｉｎ　ｔｈｅ　ｐｅｔｒｏｌｅｕｍ　ｉｎｄｕｓｔｒｙ：Ｋｉｎｅｔｉｃｓ　ｏｆｆｏｒｍａｔｉｏｎ［Ｃ］∥Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆ　ｔｈｅ　Ｕｎｉｖｅｒｓｉｔｙ　ｏｆ　ＴｕｌｓａＣｅｎｔｅｎｎｉａｌ　Ｐｅｔｒｏｌｅｕｍ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｙｍｐｏｓｉｕｍ，２９－３１Ａｕ－ｇｕｓｔ　１９９４．Ｔｕｌｓａ：Ｔｕｌｓａ　ｏｆ　Ｕｎｉｖｅｒｓｉｔｙ，２００６．

［３３］ＫＡＳＨＣＨＩＥＶ　Ｄ，ＦＩＲＯＯＺＡＢＡＤＩ　Ａ．Ｉｎｄｕｃｔｉｏｎ　ｔｉｍｅ　ｉｎｃｒｙｓｔａｌｌｉｚａｔｉｏｎ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　ＣｒｙｓｔａｌＧｒｏｗｔｈ，２００３，２５（３）：４９９－５１５．

［３４］ＴＡＬＡＧＨＡＴ　Ｍ　Ｒ．Ｐｒｅｄｉｃｔｉｏｎ　ｏｆ　ｉｎｄｕｃｔｉｏｎ　ｔｉｍｅ　ｆｏｒ　ｎａｔｕ－ｒａｌ　ｇａｓ　ｃｏｍｐｏｎｅｎｔｓ　ｄｕｒｉｎｇ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｉｎ　ｔｈｅｐｒｅｓｅｎｃｅ　ｏｆ　ｋｉｎｅｔｉｃ　ｈｙｄｒａｔｅ　ｉｎｈｉｂｉｔｏｒｓ　ｉｎ　ａ　ｆｌｏｗ　ｍｉｎｉ－ｌｏｏｐａｐｐａｒａｔｕｓ［Ｊ］．Ｃａｎａｄｉａｎ　Ｊｏｕｒｎａｌ　ｏｆ　Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１３，９１（４）：７９０－７９７．

［３５］ＤＡＶＩＥＳ　Ｓ　Ｒ，ＨＥＳＴＥＲ　Ｋ　Ｃ，ＬＡＣＨＡＮＣＥ　Ｊ　Ｗ．Ｓｔｕｄｉｅｓｏｆ　ｈｙｄｒａｔｅ　ｎｕｃｌｅａｔｉｏｎ　ｗｉｔｈ　ｈｉｇｈ　ｐｒｅｓｓｕｒｅ　ｄｉｆｆｅｒｅｎｔｉａｌ　ｓｃａｎ－ｎｉｎｇ　ｃａｌｏｒｉｍｅｔｒｙ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，２００９，６４（２）：３７０－３７５．

［３６］ＮＩＣＨＯＬＡＳ　Ｊ　Ｗ，ＫＯＨ　Ｃ　Ａ，ＳＬＯＡＮ　Ｅ　Ｄ．ａ　ｐｒｅｌｉｍｉｎａｒｙａｐｐｒｏａｃｈ　ｔｏ　ｍｏｄｅｌｉｎｇ　ｇａｓ　ｈｙｄｒａｔｅ／ｉｃｅ　ｄｅｐｏｓｉｔｉｏｎ　ｆｒｏｍ　ｄｉｓ－ｓｏｌｖｅｄ　ｗａｔｅｒ　ｉｎ　ａ　ｌｉｑｕｉｄ　ｃｏｎｄｅｎｓａｔｅ　ｓｙｓｔｅｍ［Ｊ］．ＡＩｃｈＥＪｏｕｒｎａｌ，２００９，５５（７）：１８８９－１８９７．

［３７］ＶＹＳＮＩＡＵＳＫＡＳ　Ａ，ＢＩＳＨＮＯＩ　Ｐ　Ｒ．Ｋｉｎｅｔｉｃｓ　ｏｆ　ｅｔｈａｎｅｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，１９８５，４０（２）：２９９－３０３．

［３８］ＫＵＭＡＲ　Ａ　Ｊ，ＭＡＩＮＩ　Ｂ，ＢＩＳＨＮＯＩ　Ｐ　Ｒ．Ｅｘｐｅｒｉｍｅｎｔａｌｄｅｔｅｒｍｉｎａｔｉｏｎ　ｏｆ　ｐｅｒｍｅａｂｉｌｉｔｙ　ｉｎ　ｔｈｅ　ｐｒｅｓｅｎｃｅ　ｏｆ　ｈｙｄｒａｔｅｓａｎｄ　ｉｔｓ　ｅｆｆｅｃｔ　ｏｎ　ｔｈｅ　ｄｉｓｓｏｃｉａｔｉｏｎ　ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ　ｏｆ　ｇａｓ　ｈｙ－ｄｒａｔｅｓ　ｉｎ　ｐｏｒｏｕｓ　ｍｅｄｉａ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｐｅｔｒｏｌｅｕｍ　Ｓｃｉｅｎｃｅａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１０，７０（２）：１０９－１１７．

［３９］刘陈伟，李明忠，梁晨，等．含水合物的油包水体系流动数值模拟［Ｊ］．应用力学学报，２０１３，３０（４）：５７４－５８０．

ＬＩＵ　Ｃｈｅｎｗｅｉ，ＬＩ　Ｍｉｎｇｚｈｏｎｇ，ＬＩＡＮＧ　Ｃｈｅｎ，ｅｔ　ａｌ．Ｎｕ－

ｍｅｒｉｃａｌ　ｓｉｍｕｌａｔｉｏｎ　ｏｆ　ｗａｔｅｒ／ｏｉｌ　ｅｍｕｌｓｉｏｎ　ｆｌｏｗｉｎｇ　ｌａｗ　ｂｙｃｏｎｓｉｄｅｒｉｎｇ　ｇａｓ　ｈｙｄｒａｔｅ　ｔｒａｎｓｉｔｉｏｎ［Ｊ］．Ｃｈｉｎｅｓｅ　Ｊｏｕｒｎａｌ　ｏｆＡｐｐｌｉｅｄ　Ｍｅｃｈａｎｉｃｓ，２０１３，３０（４）：５７４－５８０．

［４０］ＥＮＧＬＥＺＯＳ　Ｐ，ＫＡＬＯＧＥＲＡＫＩＳ　Ｎ，ＤＨＯＬＡＢＨＡＩ　ＰＤ，ｅｔ　ａｌ．Ｋｉｎｅｔｉｃｓ　ｏｆ　ｆｏｒｍａｔｉｏｎ　ｏｆ　ｍｅｔｈａｎｅ　ａｎｄ　ｅｔｈａｎｅ　ｇａｓｈｙｄｒａｔｅｓ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，１９８７，４２

（１１）：２６４７－２６５８．

［４１］孙长宇，陈光进，郭天民．水合物成核动力学研究现状［Ｊ］．石油学报，２００１，２２（４）：８２－８６．

ＳＵＮ　Ｃｈａｎｇｙｕ，ＣＨＥＮ　Ｇｕａｎｊｉｎ，ＧＵＯ　Ｔｉａｎｍｉｎ．Ｔｈｅ　ｓｔａ－ｔｕｓ　ｏｆ　ｔｈｅ　ｋｉｎｅｔｉｃｓ　ｏｆ　ｈｙｄｒａｔｅ　ｎｕｃｌｅａｔｉｏｎ［Ｊ］．Ａｃｔａ　ＰｅｔｒｏｌｅｉＳｉｎｉｃａ，２００１，２２（４）：８２－８６．

［４２］ＺＥＲＰＡ　Ｌ　Ｅ，ＳＬＯＡＮ　Ｅ　Ｄ，ＳＵＭ　Ａ　Ｋ，ｅｔ　ａｌ．Ｏｖｅｒｖｉｅｗ　ｏｆＣＳＭＨｙＫ：Ａ　ｔｒａｎｓｉｅｎｔ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ　ｍｏｄｅｌ［Ｊ］．Ｊｏｕｒ－ｎａｌ　ｏｆ　Ｐｅｔｒｏｌｅｕｍ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ，２０１２，９８／９９：１２２－１２９．

［４３］ＫＩＮＮＡＲＩ　Ｋ，ＬＡＢＥＳ－ＣＡＲＲＩＥＲ　Ｃ，ＨＡＢＥＴＩＮＯＶＡ　Ｅ，ｅｔ

ａｌ．Ｒｅｄｕｃｅｄ　ｃｈｅｍｉｃａｌ　ｉｎｊｅｃｔｉｏｎ　ｓｔｒａｔｅｇｙ　ｆｏｒ　ｈｙｄｒａｔｅ　ｃｏｎｔｒｏｌ　ｏｆｓｕｂｓｅａ　ｔｅｍｐｌａｔｅｓ　ａｎｄ　ｓｐｏｏｌｓ［Ｃ］∥５ｔｈ　Ｎｏｒｔｈ　Ａｍｅｒｉｃａｎ　Ｃｏｎｆｅｒ－ｅｎｃｅ　ｏｎ　Ｍｕｌｔｉｐｈａｓｅ　Ｔｅｃｈｎｏｌｏｇｙ，３１Ｍａｙ　２００６－２Ｊｕｎｅ　２００６，Ｂａｎｆｆ　Ｃａｎａｄａ．Ｃｒａｎｆｉｅｌｄ　Ｂｅｄｆｏｒｄｓｈｉｒｅ：ＢＨＲ　Ｇｒｏｕｐ　Ｌｉｍｉｔｅｄ，２００６．

［４４］ＳＫＯＶＢＯＲＧ　Ｐ，ＲＡＳＭＵＳＳＥＮ　Ｐ　Ａ．Ｍａｓｓ　ｔｒａｎｓｐｏｒｔ　ｌｉｍ－ｉｔｅｄ　ｍｏｄｅｌ　ｆｏｒ　ｔｈｅ　ｇｒｏｗｔｈ　ｏｆ　ｍｅｔｈａｎｅ　ａｎｄ　ｅｔｈａｎｅ　ｇａｓ　ｈｙ－ｄｒａｔｅｓ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，１９９４，４９（８）：１１３１－１１４３．

［４５］周厚安，汪波，金洪，等．川渝气田天然气水合物防治技术研究与应用进展［Ｊ］．石油与天然气化工，２０１２，４１（３）：３００－３０３，３５９．

ＺＨＯＵ　Ｈｏｕ′ａｎ，ＷＡＮＧ　Ｂｏ，ＪＩＮ　Ｈｏｎｇ，ｅｔ　ａｌ．Ｎａｔｕｒａｌ　ｇａｓｈｙｄｒａｔｅ　ｐｒｅｖｅｎｔｉｏｎ　ｔｅｃｈｎｏｌｏｇｙ　ａｎｄ　ｉｔｓ　ａｐｐｌｉｃａｔｉｏｎ　ｐｒｏｇｒｅｓｓ　ｉｎＳｉｃｈｕａｎ　ａｎｄ　Ｃｈｏｎｇｑｉｎｇ　ｇａｓ　ｆｉｅｌｄ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　ｏｆＯｉｌ　＆Ｇａｓ，２０１２，４１（３）：３００－３０３，３５９．

［４６］戴兴学，杜建伟，唐翠萍，等．化学类添加剂抑制天然气水合物形成的实验研究［Ｊ］．石油与天然气化工，２０１１，４０（１）：１１－１４．

ＤＡＩ　Ｘｉｎｇｘｕｅ，ＤＵ　Ｊｉａｎｗｅｉ，ＴＡＮＧ　Ｃｕｉｐｉｎｇ，ｅｔ　ａｌ．Ｅｘｐｅｒ－ｉｍｅｎｔａｌ　ｓｔｕｄｙ　ｏｎ　ｃｈｅｍｉｃａｌ　ａｄｄｉｔｉｖｅｓ　ｆｏｒ　ｉｎｈｉｂｉｔｉｎｇ　ｎａｔｕｒａｌｇａｓ　ｈｙｄｒａｔｅ　ｆｏｒｍａｔｉｏｎ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　ｏｆ　Ｏｉｌ　＆Ｇａｓ，２０１１，４０（１）：１１－１４．

［４７］ＦＩＤＥＬ－ＤＵＦＯＵＲ　Ａ，ＧＲＵＹ　Ｆ，ＨＥＲＲＩ　Ｊ　Ｍ．Ｒｈｅｏｌｏｇｙ　ｏｆｍｅｔｈａｎｅ　ｈｙｄｒａｔｅ　ｓｌｕｒｒｉｅｓ　ｄｕｒｉｎｇ　ｔｈｅｉｒ　ｃｒｙｓｔａｌｌｉｚａｔｉｏｎ　ｉｎ　ａｗａｔｅｒ　ｉｎ　ｄｏｄｅｃａｎｅ　ｅｍｕｌｓｉｏｎ　ｕｎｄｅｒ　ｆｌｏｗｉｎｇ［Ｊ］．ＣｈｅｍｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，２００６，６１（２）：５０５－５１５．

［４８］ＳＡＲＳＨＡＲ　Ｍ，ＦＡＴＨＩＫＡＬＡＪＡＨＩ　Ｊ，ＥＳＭＡＥＩＬＺＡ－ＤＥＨ　Ｆ．Ｓｔｕｄｙ　ｏｆ　ｃａｐｔｕｒｉｎｇ　ｅｍｉｔｔｅｄ　ＣＯ２ｉｎ　ｔｈｅ　ｆｏｒｍ　ｏｆｈｙｄｒａｔｅｓ　ｉｎ　ａ　ｔｕｂｕｌａｒ　ｒｅａｃｔｏｒ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ

Ｃｏｍｍｕｎｉｃａｔｉｏｎ，２００９，１９６（１１）：１３４８－１３６５．

［４９］ＧＲＥＡＶＥＳ　Ｄ，ＢＯＸＡＬＬ　Ｊ，ＭＵＬＬＩＧＡＮ　Ｊ．Ｈｙｄｒａｔｅ　ｆｏｒｍａ－ｔｉｏｎ　ｆｒｏｍ　ｈｉｇｈ　ｗａｔｅｒ　ｃｏｎｔｅｎｔ－ｃｒｕｄｅ　ｏｉｌ　ｅｍｕｌｓｉｏｎｓ［Ｊ］．Ｃｈｅｍｉ－ｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，２００８，６３（１８）：４５７０－４５７９．

［５０］ＢＡＬＡＫＩＮ　Ｂ　Ｖ，ＨＯＦＦＭＡＮＮ　Ａ　Ｃ，ＫＯＳＩＮＳＫＩ　Ｐ．Ｅｘ－ｐｅｒｉｍｅｎｔａｌ　ｓｔｕｄｙ　ａｎｄ　ｃｏｍｐｕｔａｔｉｏｎａｌ　ｆｌｕｉｄ　ｄｙｎａｍｉｃｓ　ｍｏｄｅｌ－ｉｎｇ　ｏｆ　ｄｅｐｏｓｉｔｉｏｎ　ｏｆ　ｈｙｄｒａｔｅ　ｐａｒｔｉｃｌｅｓ　ｉｎ　ａ　ｐｉｐｅｌｉｎｅ　ｗｉｔｈｔｕｒｂｕｌｅｎｔ　ｗａｔｅｒ　ｆｌｏｗ［Ｊ］．Ｃｈｅｍｉｃａｌ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ｓｃｉｅｎｃｅ，２０１１，６６（４）：７５５－７６５．

［５１］ＪＡＳＳＩＭ　Ｅ，ＡＢＤＩ　Ｍ　Ａ，ＭＵＺＹＣＨＫＡ　Ｙ．Ａ　ｎｅｗ　ａｐ－ｐｒｏａｃｈ　ｔｏ　ｉｎｖｅｓｔｉｇａｔｅ　ｈｙｄｒａｔｅ　ｄｅｐｏｓｉｔｉｏｎ　ｉｎ　ｇａｓ－ｄｏｍｉｎａｔｅｄｆｌｏｗｌｉｎｅｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｎａｔｕｒａｌ　Ｇａｓ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｅｎｇｉ－

ｎｅｅｒｉｎｇ，２０１０，２（４）：１６３－１７７．

［５２］ＫＥＮＤＯＵＳＨ　Ａ　Ａ，ＪＯＵＤＩ　Ｋ　Ａ，ＪＡＳＳＩＭ　Ｎ　Ａ．Ｔｈｅｇｒｏｗｔｈ　ｒａｔｅ　ｏｆ　ｇａｓ　ｈｙｄｒａｔｅ　ｆｒｏｍ　ｒｅｆｒｉｇｅｒａｎｔ　Ｒ１２［Ｊ］．Ｅｘ－

ｐｅｒｉｍｅｎｔａｌ　Ｔｈｅｒｍａｌ　ａｎｄ　Ｆｌｕｉｄ　Ｓｃｉｅｎｃｅ，２００６，３０（７）：６４３－

６５１．

（修改回稿日期　２０１３－１２－０６　编辑　何　明）

·

８

９

·　　　　　　　　　　　　天　然　气　工　业　　　　　　　　　　　　　　　　　　２０１４年２月

本文来源：https://www.bwwdw.com/article/6qy1.html