Introduction
天然气水合物是包含含有在轿厢1客体分子氢键的水分子的笼形结构的结晶化合物。在海底和永久冻土地区的甲烷水合物（MHS）的大量非常有趣的未来能源资源，但可能会影响到全球的气候条件2。
2013年3月，日本石油，天然气和金属国家公司进行了世界上第一个海上生产测试，使用了“减压法”3,4东部南海海槽自然承载MH沉积物中提取气体。
气体水合物可以存储气体，例如甲烷1，氢5，CO 2 1,6，和臭氧7。因此，甲烷和氢气水合物研究作为势能储存和运输介质。为了减少二氧化碳释放到大气中，CO 2 seques 2排放使用二氧化碳在深海沉积物2水合物tration进行了研究6。臭氧是目前在水净化和食品杀菌使用。臭氧保鲜技术的研究已进行了，因为它是化学上不稳定7。在水合物的臭氧浓度是在臭氧化水或冰7比要高得多。
为了开发天然含MH沉积物和基于水合物技术的天然气产量，就必须了解天然气水合物的热性能。然而，气体水合物含沉积物的热性质数据和模型研究很少8。
的“减压方法”可用于通过减少低于水合物稳定的孔隙压力以解离氢的沉积物孔隙空间。在这个过程中，沉淀物的孔隙空间的分量，从水和氢水，氢，和气体发生变化。热性能“测量后者条件是困难的，因为氢的熔化热可能会影响测量结果。为了解决这个问题，村冈等在氢形成9进行热性质“测量在过冷的条件。
与此视频协议中，我们解释过冷合成砂水煤气氢样品的测定方法。
图1示出用于测量人工甲烷水合物轴承沉积物的热性能的实验装置。安装是如图参考9相同。该系统主要包括一高压容器中，压力和温度的控制，以及测量系统的热性能。高压容器用140毫米的内径为140毫米的高度由圆柱形不锈钢的;其与死体积除去内部容积是2110厘米3，而其压力限制为15兆帕。该transie nt的平面源（TPS）技术被用于测量热性能10。九与2.001毫米个人半径TPS探头放入容器内。九个探针9的布局示于图2中的参考9。在TPS探针被连接到热性能'分析器用电缆和在实验过程中手动切换。在TPS传感器，连接图，和设置在容器的细节示于图S1中 ，图2和的支持信息3中参考9。
图1：该实验装置用于测量所述人工甲烷水合物轴承沉积物的热性能的数字是从参考9修改。3956fig1large.jpg“目标=”_空白“>点击此处查看该图的放大版本。
使用了TPS方法测量每个样品的热性能。该方法的原则在参照图10进行说明。在该方法中，依赖于时间的温度上升，ΔTAVE，是
哪里
在等式1中，W 0是来自传感器的输出功率，r是传感器探头的半径，λ是样品的热导率，α为热扩散率，t是从电源供给开始的时间到传感器探头。D（τ）是一个无量纲时间相关的功能。τ （αT/ R）1/2给出。在等式2中，m是TPS探针的同心环的数量和I 0是一个修正贝塞尔函数。样品的热导率，热扩散率和比热同时由施加到作为电源被供给到传感器探头的温度上升反演分析确定。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
注：请查阅所有相关的材料安全数据表作为本研究采用高压可燃甲烷气体和大型高压容器。戴上头盔，护目镜和防护靴。如果温度控制系统停止，在容器内的压力与氢离解增加而增加。为了防止发生意外，强烈建议使用安全阀系统自动甲烷气体释放到大气中。安全阀系统可以在没有电源供给工作。
1.沙-水-甲烷气体样品9的制备
- 放置高压容器上的振动台上。
- 倾1.5升纯水中水瓶，4000克二氧化硅砂在砂瓶中。准确称取沙子和水的群众在沙和水的瓶子，分别。
- 倒入1升纯水中的高压容器从水瓶2,110 cm 3的内部容积直到水填满半内部容器。
- 打开振动台振动整个容器。设置振动速度和电源以分别为50赫兹和220瓦，。施加振动，直到步骤1.5的完成。通过振动容器移除在容器的底部的漏极线及烧结金属过滤器的残留空气。
- 倾从砂瓶子进入容器3300克硅砂在约1g秒的恒定速率-1使用而整个容器被振动，以确保均匀的填料的水表面附近保持的漏斗。
- 停止振动时水到达容器的边缘。
- 放置一个环作为在容器上，以防止水溢出的边缘暂时壁。
- 在50赫兹和220 W.再次振动该容器
- 当沙到达容器（高度140毫米）的边缘，关闭振动。
- 删除使用次临时墙，孔隙水Ë排水管线。倒入孔隙水回水瓶。
- 在50赫兹和300瓦1秒钟一次或两次振动容器包装的沙子，如果需要添加更多的沙子。
- 称取沙子和水的群众在沙和水瓶。计算沙子和水团从沙子和水的瓶子质量差异的船只。在这个实验中，沙子和水的容器中的群众分别3385克和823.6克。的水在容器内的质量被表示为W 总 。
- 覆盖高压容器用不锈钢盖，拧紧对角线相对的对的螺栓序列。
- 从振动台移动至高压容器，以用于实验中的表。
- 覆盖高压容器与热绝缘体，用于控制温度。
- 高压管道和冷却水流线连接到高压容器中。 开的输入和输出的气体管线的阀。在800毫升分钟-1的速度，直至没有过量的水排放入大气压力下的陷阱通风10升甲烷。砂放电被固定在容器的底部的烧结金属过滤器防止。残留的水保留在砂表面上因为亲水硅砂吸附的水分子。
- 权衡的水的质量在陷阱，W 陷阱 ，以确定容器中的气体体积。确定的残留水，W 资源利用公式群众，在血管W¯¯RES = W 总 - 含 陷阱 。在这种情况下，W 水库以及w 阱分别360.6克和463.0克。
- V 砂 / V 细胞 ，其中V 沙 T的音量-使用公式Ѱ= 1确定样品的孔隙率由砂质量与砂密度（ 即 ，ρS =2630公斤米-3），和V 细胞的比率来确定他沙是在容器的内部容积。样品的孔隙率Ѱ为0.39。
- 关闭输出天然气管线的阀门。注入甲烷，以增加甲烷的孔隙压力在容器到约1210兆帕，在室温（ 即 31.6℃）。
- 关闭输入气体管线的阀。
- 开始使用数据记录器在实验过程中记录在容器内的压力和温度。数据采样的时间间隔为5秒。总的实验时间约为3000分钟。
2. MH合成和热性能“过冷样品9测量
- 打开冷却器用于通过循环冷却剂从室温容器冷却至2.0℃。让我们从冷却器T中的冷却液循环Ó容器的底部，从那里到容器的盖子，并最后返回到冷却器。在容器中的温度变化率约为0.001℃，秒-1。
- 设置使用TPS分析软件测量参数。设置传感器类型传感器设计＃7577。设置的输出功率W 0至30毫瓦和测量时间为5秒。请注意，如果传感器类型或样品条件改变相应参数应该改变。设置参数以从1℃的提高温度至1.5℃。
- 计算过冷，ΔTSUP，用下面的方程式的程度：
SUPΔT= T EQ（P） - T的。 （3）
Ť 当量 （P）为MH的平衡温度为压力P的函数Ť 当量 （P）的是使用CSMGem软件1进行计算。0; P和T是通过使用压力和温度计，分别测量该容器中的压力和温度。 - 使用TPS分析器后ΔTSUP大于2℃同时测量热导率，热扩散率，和体积比热。
- 切换连接到每个测量后的热性能分析器TPS的探针。开关实验9在TPS探头和手动分析仪之间的电缆。连接图显示在图中的S2参考9。每个传感器的切换顺序是否定的。 6→2→7→5→1→9→4→3→8→6 ...。该序列是根据在传感器之间的距离，这是设置尽可能防止影响测量的余热。收集数据每3-5分钟。
- 重复测量，直到＆＃916;吨燮 达到200°再次℃。在该实验中，ΔŤSUP最初随着时间而增加。后ΔTSUP达到最大值，因为压力带MH的形成减小ΔTSUP逐渐减小到0℃。检查，如果ΔTSUP是之前使用等式3在TPS测量大于2℃。
- 确保该温度分布不受MH熔化。如果氢在测量过程中熔化时，温度不会增加，因为氢的熔化是吸热反应。在测量过程中检查的温度分布，并在结果部分进行了讨论。
- 使用TPS技术，所有的温度分布的数据进行热性能“分析。
3.样品9,11的饱和度变化的计算
注意：饱和为氢，水和气体的样品作为时间t的函数的程度使用的气体的状态方程来计算。所使用的计算的细节和方程先前所述11。
- 计算出的甲烷气体体积V 气体 ， 叔在时间t
其中Q是气体的容器中的初始体积V MH，叔- 1是氢的在时刻t的体积- 1，和R VHW是水和MH的体积比。
在公式5 中，n是MH的水合数（〜6），ρMH，ρ 水分别对应于氢和水的密度，以及w MH和w的 水表示氢和水，r的分子质量espectively。 - 计算量ΔM 吨 （摩尔）MH的从T形成- 1到T
其中，R是气体常数，P为甲烷气体的压力，和Z 吨 （T 气体 ，T，P 的气体 ，t）是甲烷的在时刻t的压缩系数。我们9和坂本等 11所使用的本尼迪克特-韦布-鲁宾（BWR）方程，Lee和凯斯勒改性，计算ž 牛逼 12，13。对于这个计算，公式（3-7.1） - 7参考13 - BWR的13式和李-凯斯勒常数（3-7.4）在表3中。 - 计算体积变化ΔV MH，
其中P 小号 是101325 Pa的基准压力，T s是273.15 K时的基准温度中，Z s是在P点S中的压缩系数和T 秒 （ŽS〜1）和V CH 4是在甲烷气体体积的比率MH的单位体积[米3米-3。使用的165.99为V CH4值[米3米-3。 - 计算体积V MH， 叔 MH的在时间t
- 计算水中v 水在压力容器的体积， 叔在时间t
其中V 水，1是水的初始体积。 - 重复使用公式计算。 4-9在时间t = 2，3，...，以确定在水，甲烷，和MH 11的饱和的变化。初始条件为t = 1， 即 ，V 气体 ，1 = Q。 t时刻的P和T是从数据日志9服用。计算结果示于下面的部分。
其中V 水，1是水的初始体积。 Subscription Required. Please recommend JoVE to your librarian.
Representative Results
图2a示出了不受MH熔化。ΔTc是温度的变化，由于热常数'测量的温度分布; 图2b示出了由氢熔化影响的温度分布。在图2b的配置文件不能通过方程1和2进行分析，因为这些方程通过假设稳定样品条件的。
图3a示出的压力，温度，和在容器中作为时间的函数过冷却度。后系统MH成核已达到的压力和温度平衡。 MH的形成是通过在时间t = 170分钟的急剧压力变化标志。双头箭头表示过冷却度大于2℃。热常数在此范围内进行测定。
图4示出的热常数的测定值的一个例子。实验条件T = 825分钟， 磷= 7.1兆帕，T = 2.4℃下，S H = 0.16，S G = 0.53，和S W = 0.31。 图4a示出的温度分布。在TPS分析软件记录在时间超过一个预定的时间间隔等间隔的200个数据点;从而，数据被选择为从200个数据点的分析。双头箭头表示在分析中所使用的数据的范围。时间范围的分析1和2分别是0-5秒，0.65-4.88秒。分析图1和2分别是不恰当的，适当的范围的例子。 图4b和使用每个分析范围在TPS技术得到图4c， 图4b示出的温度变化ΔTAVE（τ）和D（τ）与ΔTAVE（τ ）ΔT= C（T）。 （ΔTAVE之间的关系＃964;）和D（τ）的变化取决于分析范围图4c示出的温度T d比的时间t的平方根。从由TPS反演分析得到的线性拟合的温度数据的偏差是T d 设定 。分析1中的测量开始时的偏差是相当大的，如在图4c中 ，这表明在TPS传感器探头的绝缘层会影响测量结果所示。
表1列出了热常数在每个分析范围，如上所述。总到特征时间比由总分析时间限定（对于t = 2-4秒，总时间为4秒）由特性时间τ划分。需要注意的是，总以特性时间比使用TPS技术时应该小于1。这在参照
避免影响测量传感器探头，不应被用于在每个测量的开始的数据。 T d 设定的平均偏差被最小化， 如图4c所示 ，通过调整分析时间范围。总到特征时间比通过调整分析时间范围调整到统一。因此，我们通过分析从2不是1的热常数的值。
热导率，比热和热扩散率被示为时间分别在图5a，b和c，的函数。最后，我们总结了第t结果hermal性能和水合物饱和度。有关详细结果在二段给出。 4 REFFERENCE 9。
图2：温度分布为时间（ 一 ）不受MH熔化（过冷条件）及（b）受MH熔化功能，请点击这里查看该图的放大版本。
注意，这两个温度分布是从初步实验。测量时间比在实验更长，以便澄清熔融的热的影响。在初步实验中，测量时间t为40秒，并在输出功率W 0为20毫瓦（a）和50毫瓦（b）中。
<p系列=“jove_content”FO： - together.within页保留=“1”>
图3：（a）压力，温度，和在容器中作为时间的函数过冷却度。双头箭头表明超冷却程度大于2℃。热常数在此范围内进行测定。 （ 二 ）样品的氢，水和甲烷气体饱和度显示为的时间（从参考文献9转载）函数。 请点击此处查看该图的放大版本。
图4： 热常数的测定值的分析例 （a）使用的温度曲线。TPS测量方法。时间范围的分析1和2分别是0-5秒，0.65-4.88秒。 （b）与ΔTAVE（τ）的温度变化ΔTAVE（τ）和D（τ）的关系ΔT= C（T）。 （ 三 ）温T d比平方根时间t。 请点击此处查看该图的放大版本。
图5：（a）热导率λ作为时间的函数，（B）的比热ρC 数p作为时间的函数，以及（c）热扩散α为时间的函数。结果被转换为热性能如MH饱和的函数。转换后的结果和相关的讨论中报告参考9.数据显示范围T = 210-980分钟内的重叠。为了清楚，绘制数据代表三次测量从在该范围内的相同的传感器的平均值。这些数字已经从参考9修改。 请点击此处查看该图的放大版本。
| 分析范围，S | λ，W M - 1 的 K - 1 | ρCP，MJ米- 3 的 K - 1 | α，毫米2秒- 1 | 达尔通道AR。时间 | 意思是开发，℃ | |
| 分析1 | 0.00 - 5.00 | 2.12 | 0.938 | 2.26 | 2.11 | 0.01018 |
| 分析2 | 0.65 - 4.88 | 2.31 | 2.11 | 1.10 | 1.00 | 0.00061 |
表1： 对每个分析范围热常数分析图1和2分别是不恰当的，适当的范围的例子。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
估计MH的测量形成热的影响。 MH的形成热量从中文变化率的产品估计， 如图3b和形成H = 52.9千焦摩尔-1 MH 14的焓。因此，最高温度变化是0.00081℃，秒-1。这是比5秒的时间间隔期间1℃，1.5℃之间的TPS传感器的温度上升ΔTÇ低得多。详细的估计和讨论二段描述。 4参考9。
以下是至关重要的协议的步骤。第一步是保持样品过冷条件。第二步骤是通过保持在TPS传感器的温度上升ΔTÇ过冷却ΔTSUP的程度下进行热常数的测定值。
以确保测量不受温度漂移的影响，以下应被确认。首先，确保大宗温度变化比TPS传感器的温度升高ΔTÇ低得多。第二，要确保由于MH的生成热的温度变化比在TPS传感器的温度上升ΔTÇ低得多。
如果样品融化时，热导率和比热将由TPS技术发散到无限。在这种情况下，来自传感器改变输出功率或减少测量时间。
这种测量方法可应用于气体水合物-水-旅客气系统，其中包含氢气，CO 2和臭氧水合物的热特性，因为天然气水合物的特性低形成速率不是唯一的氢。这种方法的关键点是低鼠靶材料的相变的即因此，这种方法可以被应用到其它的材料具有低相变率。该测定方法也可以应用于从低浓度的THF溶液和四丁基溴化铵（TBAB）水合物形成四氢呋喃（THF）水合物，如果这些水合物的形成速率是超冷却的条件下足够慢。这里唯一的要求是要确保由于形成水合物的热的温度变化比在传感器的温度上升低得多，如上所述。另一方面，这种技术不能应用到水冰和化学计量THF溶液水合物的相变，因为在这些系统中的转变率是非常快，形成的热显著影响测量。
韦特等人 15测量包括砂，甲烷气，和MH样品的热导率。 Kumar 等。16测量用样品相同的部件的热扩散率。它们直接在砂孔隙形成MH通过在加压的甲烷气体环境中利用的水冰。所有的水冰的转化为氢。因此，他们测得的样品的热导率，直到氢的形成完全停止。这种方法具有这样的热特性的测量结果不受氢的形成或分解热，样品组合物是恒定的优点。但是，这种方法不能得到含有沙，水，甲烷和氢样品的热性能。黄和风扇测量的水合物支承砂样品17的热导率。它们在使用十二烷基硫酸钠（SDS）溶液，这有利于对MH形成砂孔形成氢。他们指出，燃气及水的可能留在毛孔沙子和气体显著影响测量结果。然而，他们没有报告笏的组成ER和气体。我们的测量协议具有赋予的热性能（热导率，热扩散率，和体积比热）和组合物中的MH-轴承沉积物包括沙，水，甲烷，和MH之间的关系的优点。
开发天然气水合物的大规模生产技术，需要对水合物形成的热常数，所提出的测量方法正是这么做的。
Subscription Required. Please recommend JoVE to your librarian.
Acknowledgments
这项研究是财政上MH21研究联盟甲烷水合物资源日本经济，贸易和工业部国家甲烷水合物开发计划的支持。笔者想感谢T.前川和S.转到他们与实验的援助。
经许可后转载自数字（村冈，M.，Susuki，N.，山口，H.，辻，T.，山本，Y.，能源燃料，29（3），2015年，1345年至1351年，2015年，DOI： 10.1021 / ef502350n）。版权所有（2015）美国化学学会。
Materials
| Name | Company | Catalog Number | Comments |
| TPS thermal probe, Hot disk sensor | Hot Disk AB Co., Sweden | #7577 | Kapton sensor type, sensor radius 2.001 mm |
| Hot disk thermal properties analyzer | Hot Disk AB Co., Sweden | TPS 2500 | |
| Toyoura standard silica sand | Toyoura Keiseki Kogyo Co., Ltd., Japan | N/A | |
| Methane gas, 99.9999% | Tokyo Gas Chemicals Co., Ltd., Japan | N/A | Grade 6 N, Volume 47 L, Charging pressure 14.7 MPa |
| Water Purification System, Elix Advantage 3 | Merck Millipore., U.S. | N/A | 5 MΩ cm (at 25 °C) resistivity |
| Vibrating table, Vivratory packer | Sinfonia Technology Co. Ltd., Japan | VGP-60 | |
| Chiller, Thermostatic Bath Circulator | THOMAS KAGAKU Co., Ltd., Japan | TRL-40SP | |
| Coorant, Aurora brine | Tokyo Fine Chemical Co.,Ltd., Japan | N/A | ethylene glycol 71 wt% |
| Temparature gage | Nitto Kouatsu., Japan | N/A | Pt 100, sheath-type platinum resistance temperature detector |
| Pressure gage | Kyowa Electronic Instruments., Japan | PG-200 KU | |
| Data logger | KEYENCE., Japan | NR-500 | |
| Mass flow controller | OVAL Co., Japan | F-221S-A-11-11A | Maximum flow 2,000 N ml/M, maximum design pressure 19.6 MPa |
References
- Sloan, E. D., Koh, C. A. Clathrate Hydrates of Natural Gases, 3rd ed. , 3rd ed, CRC Press. Boca Raton, FL. (2007).
- Hatzikiriakos, S. G., Englezos, P. The relationship between global warming and methane gas hydrates in the earth. Chem. Eng. Sci. 48 (23), 3963-3969 (1993).
- Yamamoto, K. Overview and introduction: pressure core-sampling and analyses in the 2012-2013 MH21 offshore test of gas production from methane hydrates in the eastern Nankai Trough. Mar. Petrol. Geol. 66 (Pt 2), 296 (2015).
- Fujii, T., et al. Geological setting and characterization of a methane hydrate reservoir distributed at the first offshore production test site on the Daini-Atsumi Knoll in the eastern Nankai Trough, Japan. Mar. Petrol. Geol. 66 (Pt 2), 310 (2015).
- Mao, W. L., et al. Hydrogen clusters in clathrate hydrate. Science. 297 (5590), 2247-2249 (2002).
- Lee, S., Liang, L., Riestenberg, D., West, O. R., Tsouris, C., Adams, E. CO2 hydrate composite for ocean carbon sequestration. Environ. Sci. Technol. 37 (16), 3701-3708 (2003).
- Muromachi, S., Ohmura, R., Takeya, S., Mori, H. Y. Clathrate Hydrates for Ozone Preservation. J. Phys. Chem. B. 114, 11430-11435 (2010).
- Waite, W. F., et al. Physical properties of hydrate-bearing sediments. Rev. Geophys. 47 (4), (2009).
- Muraoka, M., Susuki, N., Yamaguchi, H., Tsuji, T., Yamamoto, Y. Thermal properties of a supercooled synthetic sand-water-gas-methane hydrate sample. Energy Fuels. 29 (3), 1345-1351 (2015).
- Gustafsson, S. E. Transient plane source techniques for thermal conductivity and thermal diffusivity measurements of solid materials. Rev. Sci. Instrum. 62 (3), 797-804 (1991).
- Sakamoto, Y., Haneda, H., Kawamura, T., Aoki, K., Komai, T., Yamaguchi, T. Experimental Study on a New Enhanced Gas Recovery Method by Nitrogen Injection from a Methane Hydrate Reservoir. J. MMIJ. 123 (8), 386-393 (2007).
- Lee, B. I., Kesler, M. G. A generalized thermodynamic correlation based on three-parameter corresponding states. AIChE J. 21 (3), 510-527 (1975).
- Reid, R. C., Prausnitz, J. M., Poling, B. E. Chapter 3, Unit 3, 7. The properties of gases and liquids. , 4th ed, 47-49 (1987).
- Anderson, G. K. Enthalpy of dissociation and hydration number of methane hydrate from the Clapeyron equation. J. Chem. Thermodyn. 36 (12), 1119-1127 (2004).
- Waite, W. F., deMartin, B. J., Kirby, S. H., Pinkston, J., Ruppel, C. D. Thermal conductivity measurements in porous mixtures of methane hydrate and quartz sand. Geophys. Res. Lett. 29 (24), 82-1-82-4 (2002).
- Kumar, P., Turner, D., Sloan, E. D. Thermal diffusivity measurements of porous methane hydrate and hydrate-sediment mixtures. J. Geophys. Res. 109 (B1), (2004).
- Huang, D., Fan, S. Measuring and modeling thermal conductivity of gas hydrate-bearing sand. J. Geophys. Res. 110 (B1), (2005).
