高压反应釜实验装置

『壹』 将物质进行高温高压处理以后，物质会发生那些变化可以用什么装置做高温高压试验

经过高温高压处理以后物质的性质肯定会有很多变化，比如内部结构啊，物相啊，熔点啊，沸点啊，导电性，导磁性，硬度啊等都可能会发生变化。
高温高压装置有：
1.活塞圆筒装置：提供的是密封的样品仓且样品仓很小，压力4GPa，温度2000℃
2.高温高压反应釜：压力是几百MPa，温度是几百℃
3.多面顶：压力几十GPa，温度两千多摄氏度

『贰』 惰性气体保护下在230摄氏度进行的有机合成所需要的仪器有哪些

高温高压反应釜该产品适用于进行化学、化工反应时，需要同时测定温度，添加惰性气体，压力，PH值， 电导，氧化还原电位的精密化学合成反应的实验。应用实例： 催化，高温高压合成，动力学检测等。
技术参数如下：
◆材质：釜体为纯钛材质/或316L不锈钢，操作台面为316L不锈钢，耐工作压力6MPa.工作温度：400℃以下。
◆加热方式：直接在反应釜大包围加热，简单方便，升温快，使用性强，用系列电热纤维炉/或电热环直接加热控温，加热功率：300W-500W。控温精度：±1℃
◆阀门管件为德国进龘口，接口均为标准配置。
◆搅拌速度，温度均为数字显示0-1500RPM，转速精度±1RPM
◆选用专龘利插入式针式K型传感器，工作温度：400℃以下 显示精度：±1℃标配含搅拌装置，纯钛材质反应釜本体，电热纤维炉，插入针式K型传感器等。
◆可充惰性气体，保护反应物不与空气直接接触。
◆可插入温度，压力，PH值，电导率传感器，保证准确测量相关数据。
◆技术亮点（1）：让繁琐的实验简单化，可实现在同一操作台面上不同环境下或同环境不同温度/不同压力的多元化串联实验，25ML/50ML/100ML/200ML250ML/500ML单独或并联、串连使用。
◆技术亮点（2）：可通过RS232通讯接口与计算机连接进行打印和显示用户需要的历史数值和实际曲线，通过软件实时监控整个系统。

『叁』 什么是高压反应釜

化学实验常常能带来精彩绝伦的视觉盛宴，但易燃、易爆、有毒介质的化学反应，却往往成为“不安分子”，增加了实验风险。化学实验中，美丽与危险并存，磁力/机械搅拌高压反应釜，兼具防止泄露、保障安全的可靠与控温精准、操作简单的高效，专为进行高温、高压下的化学反应而来。

1、体积小巧，整机结构简单可靠

力辰磁力/机械搅拌高压反应釜采用静密封/石墨密封结构，体积小巧、整机结构简单可靠，具有密封性能好、无任何泄露和污染的特点。
它是进行高温、高压下的化学反应的理想装置，更适合用于各种易燃易爆、剧毒、贵重介质及其它渗透力极强的化学介质进行搅拌反应。

『肆』 高压反应釜如何通氮气

用真空油泵连接高压反应釜进气阀。
用真空油泵连接高压反应釜进气阀，抽取釜内空气后，关闭进气阀派坦。卸下真空油泵管道。连接氮气管至进气阀，开启氮气调压阀及高压反应釜。
高压反应釜（磁力高压反应釜）是磁力传动装置应用于反应设备的典型创新，它从根本上解决了以前填料密封、机械密封无法克服胡羡携的轴封泄漏问题，无任何泄漏和污染，是国内目前进行高温、高压下的化学反应最为裤伏理想的装置，特别是进行易燃、易爆、有毒介质的化学反应，更加显示出它的优越性。

『伍』 高压反应釜的基本介绍

高压反应釜釜体材料主要采用 1cr18Ni10Ti 不锈钢制作，并可根据不同介质要求制作钛材（ TA2 ）、镍（ Ni6 ）及复合钢板，釜体结构有平盖、凸形盖以及带人孔的闭式反应釜体，釜盖上的开孔可根据用户要求进行设计；加热方式有夹套蒸气、夹套热油电加热等形式供用户订货时任意选配。 对有抛光要求的釜体内表面，可达到 以上的镜面抛光水平，对高粘度的物料加工成锥形底，便于放料、清洗。

智能高压反应釜
磁力传动装置是高压反应器实现物料混合的核心技术，是设计压力超过50bar的必须部件。专利公开后成为压混合设备的行业标准。
磁力传动装置将真正用于搅拌的部件密封在压力环境内，以此彻底解决了搅拌轴与反应器之间动密封的可靠性。由于磁力传动装置的基本制造材料均为铁磁性物质（铁基合金，镍基合金），磁导率极佳，磁场在穿过用于密封合金层时几乎没有损耗。由此外部驱动永磁体产生的旋转磁场可以有效被内部永磁体转子接受，进而带动反应器内的搅拌部件进行物料混合。

『陆』 水合物热物理参数的实验测定

自然界中水合物有99%是甲烷水合物，直接研究甲烷水合物的热物理参数有重要的实际意义。然而，这一工作长期以来困难重重，一方面是由于实验室合成甲烷水合物过程中存在“铠甲”效应(即外部的水合物生成后会形成一层厚厚的“壳”，导致生成不够彻底，生成的水合物中夹杂大量的气体、水和冰粒，并且合成过程非常缓慢);另一方面是实验技术和测试方法的局限性。青岛海洋地质研究所水合物实验室研制出一套沉积物中水合物分解过程中的热物理特性模拟实验装置，包括一个可编程控制变温实验箱一台，高压模拟实验装置一套，研制了高压热-TDR探针，购置了TDR仪和数据采集器，制作了高压和温度监测系统，并研制了计算机控制与数据采集系统一套。该实验装置的技术核心在于热-TDR探针的设计制作。TDR技术和热脉冲技术具有相对独立的探头，我们将二者有机结合，可以实现同时同地测量介质含水量、温度、容积热容量、热导率、热扩散系数等多项参数。不但避免了介质时空变异性的影响，还可以实现连续定位测定。

实验装置

如图75.13所示，模拟实验系统硬件部分包括可编程步入式变频高低温箱一台、高压模拟实验装置一套(包括其核心技术———耐高压热-TDR探针)、数据采集系统，软件部分我们自行设计编写了计算机控制与数据采集系统。

图75.13 实验装置简图

高压模拟实验设备主体部分是增压系统、两个高压釜体及插入反应体系中的热-TDR探针。气高压气瓶顶端有两个压力控制阀门，用于控制气瓶输出压力和釜体输入压力。阀门连接两个压力指示表，可以直接读出两处压力值，便于控制加压幅度。

高压釜体包含一个反应釜体和一个为搅拌釜体。两个高压釜体容积均为200cm³，最大工作压力30MPa。高压反应釜外层用不锈钢制作，采用自紧法螺纹密封，为保证螺纹密封效果，在连接部分采用两个O型密封圈进行密封。整个反应釜也是专门设计定做，通过测试，其密封效果可以保证实验顺利完成。搅拌釜体内装有聚四氟磁棒，下部是磁力搅拌器。反应釜体内部装有内筒(内筒用聚砜材料切割制成，聚砜具有力学性能优异，刚性大、耐磨、耐高压、热稳定性好等特点，适合在低温高压条件下作为水合物的反应容器材料)。容积为70cm³。热-TDR探针插入内筒所盛的反应物中发射热脉冲和测定反应体系温度、含水量等参数。压力表直接连接在气体管路上，便于采集数据和人工监控。

实验技术与方法

将沉积物装入模拟装置，采用逐渐升压的办法，测量压力对热物理参数的影响。当模拟装置内的压力达到预定的压力条件时，停止加压。室温下模拟装置放置一定的时间后，若压力没有发生变化即可开展水合物生成模拟实验(压力恒定48h)。启动监测装置，监测模拟装置内，温度、压力和TDR波形的变化。随着水合物逐渐生成，TDR波形逐渐发生变化，反射系数逐渐增加，相对距离缩短。

打开搅拌釜、反应釜进气阀门(阀门3、4、5)，打开抽真空口(阀门2)，其余阀门关闭，将系统抽真空。待系统负压稳定后，关闭抽真空口和抽真空机。打开除高压阀以外的所有阀门，通入实验所用的甲烷气清洗气路，重复3～4次。然后打开进气阀门(阀门1、3、4)，其余阀门关闭，开始向两个高压釜内加压。加至实验所需压力(4.0～7.0MPa)后关闭加压阀门稳定一段时间。打开磁力搅拌器直至搅拌釜内的甲烷气溶解在SDS溶液中达到饱和。打开搅拌釜和反应釜之间的阀门(阀门5)，使溶解了饱和甲烷气的SDS溶液流向反应釜，直至反应釜中的松散沉积物达到含水量饱和状态后关闭阀门5。开启控温箱开关，将温度设置为0.5℃。实验进入水合物合成阶段。水合物合成所需时间受多个条件影响，如水合物的“记忆效应”、温度“过冷度”、表面活性剂的添加等。水合物合成一般需要1、2d时间。水合物生成进度可以通过TDR波形图明显看出。

实验选择的热脉冲电源为12V直流电源。通过计算机直接控制热脉冲发射的时间和时长。由于实验采用的加热丝直径很小，加热时间过长容易导致加热丝绝缘层烧化;另外，水合物本身遇热容易分解。综合考虑上述各因素影响，加热时长一般掌握在4～8s内。另外，一个热脉冲发射过后需要一定的散热时间，待反应体系温度完全恢复到脉冲发射前的状态时再发射下一个脉冲。两个相邻的热脉冲之间发射间隔过短，反应体系内的余温会干扰实验结果;间隔过长则费时费电(刁少波等，2008)。

计算

采用了平行热线法和交叉热线法测定热物理特性。

1)交叉热线法。热导率计算公式为:

岩石矿物分析第四分册资源与环境调查分析技术

式中:λ为热导率;"₂、"₁为两次热脉冲的加热时间;T₂-T₁为温度的变化;Q为热源强度。

2)平行热线法。热扩散系数α计算公式为:

岩石矿物分析第四分册资源与环境调查分析技术

式中:t_m为达最高温度时的时间;t₀为热脉冲的加热时间;r为热电偶距线性热源的垂直距离。

容积热容量计算公式为:

岩石矿物分析第四分册资源与环境调查分析技术

式中:q为单位长度加热丝在单位时间内释放的热量;E_i(-x)为指数积分。

由λ=α·ρ_c计算出热导率。

『柒』 高压反应釜怎么做气密性实验，工作压力4.8MPa

常温加水达釜体约80％，然后置氮气（逐级升压，密闭后，观察釜内气压变化）最终至6.0MPa．同时用起泡液（肥皂水等），对各封闭，连接部件等检查．

『捌』 含水合物沉积物的纵波声速测定

李风光¹，陈光进¹，孙长宇¹，李清平²，郭绪强¹，杨兰英¹，潘恒¹

李风光（1985－），男，博士生，主要从事水合物研究，E-mail:[email protected]。

1.重质油国家重点实验室/中国石油大学，北京102249

2.中海石油研究中心，北京100027

摘要：自行设计搭建了天然气水合物纵波声速（V_P）测定实验装置，主要测定含水合物样品的声学性质。为使水合物在沉积物中分布均匀并能完全填满沉积物孔隙，采用四氢呋喃（THF）水溶液与甲烷气体在沉积物中生成水合物，测定了水合物生成过程中体系的声学性质变化，并分析了沉积物粒径及四氢呋喃水溶液浓度对声速的影响。实验结果表明：在水合物生成过程中，声速随着水合物饱和度的增加而增加，最后趋于定值；四氢呋喃水溶液的浓度越高，最后水合物样品的声速也越大，然而沉积物颗粒大小对声速结果几乎没有影响；波形的振幅变化表明在水合物生成过程中振幅随着水合物的生成先达到一最大值，随后受游离的甲烷气体的影响，振幅又会逐渐减小，最后趋于定值。根据水合物胶结模型的计算结果，模型计算值与实验值基本吻合。

关键词：水合物；沉积物；纵波声速；振幅；饱和度

Investigate on the P-Wave Velocity of Hydrate-Bearing Sand

Li Fengguang¹,Chen Guangjin^l,Sun Changyu¹,Li Qingping²,Guo Xuqiang¹,Yang Lanying¹,Pan Hen¹

1.State Key Laboratory of Heavy Oil Processing/China University of Petroleum,Beijing,102249,China

2.CNOOC Research Center,Beijing 100027,China

Abstract:An experimental apparatus was developed to measure P-wave velocity (V_P) o f gas-hydrate-bearing sediment.Tetrahydrofuran (THF) was added to quicken the hydrate formation in the porous media and to synthesize hydrate-bearing sediments with uniformdistribution.Methane acted as afree gas to participate in the hydrate formation.Five experimental runs were performed to examine the influence of sediment grain size and THF concentration on V_P.The P-wave velocity and the amplitude for the first arrival wave signal were collected in real time ring hydrate formation process.The experimental data showed that V_Pincreases monotonically with the increase of hydrate saturation in the sediment pore space and finally tends to be a constant value.This final V_Pvalue increases with the increase of initial THF content,but the effect of sand grain size on V_Pis inconclusive.The variations of amplitudefor the first arrival wave signal with elapsed time ring hydrate formation illustrates that the amplitude increases with the increase of hydrate saturation until it attains a maximum value and then decreases graally e to the effect of free methane gas penetrating into the hydrate-bearing sediment.The acoustic velocity of THF-hydrate filled sediment was also predicted based on the extended contact cement theory.The predicted results were close to the experimental data obtained in this work.

Key words:hydrate;sediments; P- wave Velocity;amplitude; saturation

0 引言

根据目前地质探测的结果，水合物主要分布在海底大陆架的沉积物及冻土带中^[1]。由于其储量巨大，作为一种潜在的能源资源，水合物的研究开展十分广泛^[2-4]。掌握含水合物沉积物的性质，如岩性、水合物饱和度、渗透性、密度、声速等物性，对水合物资源的估计以及将来的开采工作都有重要的意义^[5-10]。

受水合物存在条件的限制，对地层水合物样品进行原位性质测试异常困难，因而常常在实验室中人工合成水合物样品以测定其相关的物性^[8-9,11]。为了准确的测定水合物沉积物的物性，合成具有代表性的水合物样品至关重要。静态下甲烷气体在水中的溶解度很小，因而利用溶解的甲烷气合成水合物样品需要耗费相当长的时间，这给含水合物样品物性的研究带来很大的困难。Waite等^[8]测定了部分饱和水中含水合物沉积物样品的声学性质。研究表明沉积物中水合物饱和度相对有限，而且水合物生成极其缓慢。实验进行了1 400多h，沉积物中水合物的饱和度最高为70％，声波最后趋于稳定。在实验室条件下，合成水合物能完全填满沉积物孔隙的水合物样品相对较难。众所周知，四氢呋喃溶液跟水可以以任意比例混合，四氢呋喃的存在能加快水合物的生成速率，水合物在沉积物中可以均匀分布，因此实验室中经常用四氢呋喃替代甲烷合成水合物样品^[12,13]。

声速是一重要的地球物理性质参数，可以反映岩性、水合物丰度、矿藏分布等重要信息。实验室测定的声速数据可以为地震勘探的测井解释提供解释依据。勘探地层水合物最常用的方法是地震法，地震测井资料的解释需要对含水合物沉积物的物性准确了解，然而，含水合物沉积物样品的保真取样极其困难因而不能进行有效的物性测试。Priest等^[14]在实验室中合成了均匀分布的不同水合物饱和度的水合物样品，并测定了其纵波声速及横波声速，实验结果表明：甲烷水合物首先胶结砂粒，随后填充孔隙。Pearson等^[12]测定了四氢呋喃和水的混合物中生成的水合物样品的声速值。对Berea砂岩和Austin chalk样品，从水合物开始生成到水合物样品合成只需极短的时间，其声速变化都增长的很快，前者从2 500 m/s到4 500 m/s，后者从1 400 m/s到5 000 m/s；进一步的降温并没有降低样品的声速。然而，四氢呋喃生成Ⅱ型水合物，并且只占据水合物晶格中的大孔。这跟天然沉积物中的水合物区别较大，自然的水合物样品的晶格中大小孔主要被甲烷分子填充。本文利用实验室设计的水合物声速测定装置测定了含水合物沉积物样品的声学参数。四氢呋喃作为一种水合物的生成促进剂，并能使水合物在沉积物中分布较均匀，甲烷作为游离气体参加水合物生成反应。本文合成的水合物样品跟纯四氢呋喃水合物相比，四氢呋喃主要填充水合物晶格的大孔穴，而甲烷则填充小孔穴，跟天然水合物样品具有可比性。在实验过程中，V_p波形图通过软件记录下来，通过对波形图的分析可以得到声速、振幅等声学性质在水合物生成过程中的变化。

1 实验部分

1.1 实验装置

实验装置如图1所示。主要包括：高压反应釜、低温空气浴槽温控系统、天然气配气系统、温度压力测量系统、超声波声速测定系统5部分。

低温高压系统的核心部件是高压反应釜（图1b），釜的设计压力为32.0 MPa，由不锈钢材料加工而成，反应釜的容量为2.0 L（φ130×150），由江苏华安石油科研仪器公司制造。该反应釜的釜底和釜盖上安装有一对超声波探头，其中釜盖上的探头可以随滑杆上下移动，便于调节两探头间的样品测量距离。在高压反应釜侧壁开设有多个不同尺寸的孔径，便于安装温度和压力监测系

图1 声速测量实验装置示意图（a）和高压反应釜内部示意图（b）

1.甲烷气瓶；2,3,4,11.阀；5.高压反应釜；6.压电换能器；7.沉积物； 8.手柄；9.底部进气口； 10.Pt100; 12.空气浴；13.温度测量； 14.压力测量；15.超声波讯号发射接收仪；16.数字示波器；17.计算机

统、超声波发射接收系统及进气口和排水口等。热电阻Pt100的测量精度为±0.1 K，压力传感器测量精度为0.5％，量程为0～60 MPa；声速测量系统主要包含4部分：陶瓷压电换能器、超声波讯号发射接收仪、示波器和数据采集分析软件。其中换能器的发射主频为500k Hz～1.0 MHz，电压为400 V；美国进口5077PR数字式超声波讯号发射接收仪；Tektronix Technology公司TBS2012B型号数字示波器，采用了高速A/D数据采集卡。声速测量的误差大约±0.5％。

1.2 实验材料

沉积物材料选用松散的石英砂。实验前先将石英砂用去离子水清洗干净，然后用烘箱在393.2K下干燥12 h。然后将其筛分为20～40目， 40～60目，60～80目3种，3种目数的沉积物物性列于表1。

表1 沉积物物性表

甲烷气由北京北分气体工业有限公司提供，纯度为99.9％。四氢呋喃水溶液由纯度为99.8％的四氢呋喃和去离子水配制。四氢呋喃水溶液生成水合物的最佳的摩尔比为5.9%

1.3 实验步骤

沉积物中水合物饱和度和分布对含水合物沉积物样品的声学性质有很大影响^[9]。为合成沉积物中水合物均匀分布的水合物样品，采用以后步骤合成水合物样品：

1）先用去离子水清洗反应釜，然后将其擦干。并在反应器侧壁安装一支热电阻Pt100以测量水合物生成过程中温度的变化。

2）室温下将含饱和四氢呋喃水溶液的沉积物装入反应釜中，如图1(b）所示。

3）沉积物装入反应釜后，通过旋转手柄将沉积物样品压实。作用在手柄上的有效应力为500 k Pa，压紧后的样品长度在50mm左右。

4）连接好管线后在室温下抽真空2 min，以排除空气对实验造成的影响。空气浴温度设为278.2 K，开始水合物的生成实验。开始先通入0.5 MPa甲烷气检验装置的气密性，在通入甲烷气时同时开启声速测量软件开始采集波形图。

5）气密性检验完好后，继续从反应釜底部注入甲烷气。进气过程在1 min左右，进气结束后反应釜内压力约为12 MPa，关闭管线截止阀，水合物样品将快速合成。从进气开始连续记录温度、压力及波形图等数据，直至反应结束。

每一组水合物样品合成实验都按照上面的实验步骤进行，然后对波形图处理分析得到样品的声学参数。

2 结果与讨论

合成5个水合物样品的沉积物粒径和THF水溶液的摩尔比及声速数据均列入表2。

表2 合成的水合物样品的最终声速数据

根据表2所示，实验采用了3.0％和5.9% 2种THF水溶液的摩尔比合成水合物样品，水溶液在沉积物中均为100％饱和状态。对于饱和THF水溶液沉积物的声速差别不大，均在1 750 m/s左右；在水合物生成过程中，声速随着水合物饱和度的增加而不断增大，最后趋于一定值。THF水溶液在常压下277.6 K时以1:17(THF:H₂O）的比例合成Ⅱ型水合物^[13]，然而水合物的饱和度很难精确计算。尽管如此，根据水合物样品声速测量结果可以看出，THF浓度高时合成的水合物样品的声速要高于THF浓度低时的声速值，说明高浓度的THF溶液合成的水合物饱和度较高，因而声速较高。为进一步了解水合物生成过程中温度、压力、声速及振幅等各参数的变化情况，以样品5为例详细说明。

2.1 水合物生成过程中声速测定

图2为样品5(60～80目沉积物＋5.9%THF水溶液＋甲烷体系）在水合物合成过程中温度、压力随时间的变化关系图。从图中可以看出，在水合物生成初期，甲烷气被大量消耗，在4.2 min内压力从13.52降到11.29 MPa，说明消耗的甲烷气生成大量水合物。水合物生成为放热反应，受此影响，温度在70s内升高了25.5 K（从278.1到303.6 K）。7 h后，温度和压力变化趋于稳定。从20 h到30 h，温度和压力的降低是受降温的影响。对于水合物样品合成后温度变化对声速结果的影响，从图3声速随时间的变化可以看出，温度降低后声速没有变化。

图2 水合物生成过程中温度、压力随时间的变化

图3 水合物生成过程中声速随时间的变化

图4 图3中A,B,C,D时刻的波形图

图4为图3中在某一具体时刻（A、B、C和D）的波形图。在时刻A（图4a）为甲烷气没有完全注入前，声速为1 856 m/s，声速比含饱和T H F溶液的声速略高，说明已有少量水合物开始生成；此后，随着水合物的开始生成，在0.6 h后的B点，甲烷气体已经通入结束，声速V_P增加到3 078 m/s（图4b）；在1.3 h后的C点，声速继续增加到3 585 m/s（图4c）；反应进行到7 h的D 点，声速最后趋于定值3 827 m/s（图4d）。

2.2 声速振幅测量

图5为水合物样品5的波形振幅随反应时间的变化关系。从图中可以看出，水合物生成过程中波形的演化规律。振幅的变化与声速的变化并不相同，在水合物快速生成的前7 h，振幅随水合物饱和度的增加而不断增大；在声速趋于一定值的过程中，振幅开始不断降低最后稳定于一定值，这种变化规律与声波穿过水合物样品时的衰减具有一定的关系。根据Priest等^[15]研究结果表明，声波能量的衰减主要包括：①几何形状散射衰减；②扩散衰减；③本征衰减。几何形状散射衰减取决于实验所选材料的几何形状，而扩散衰减与沉积物粒径大小与声波的波长均有关系，这两者均可以通过选择合适的沉积物粒径的大小来避免，因而对声波信号的衰减主要受本征衰减（它主要取决于沉积物的材料、岩性、孔隙填充介质和饱和度等）及声波的频率等影响^[16]。

对于松散的沉积物，P波信号受衰减的影响无法通过沉积物，因而无法测量到波形信号。而对于饱和水的沉积物，低频的波形信号可以通过，因此可以采集到波形信号，但是波形很微弱，振幅值也较小。水合物生成后，开始生成的絮状水合物在溶液中仅仅改变了溶液的弹性性质；随着水合物的大量生成，水合物开始胶结沉积物颗粒，含水合物的沉积物的刚性随着水合物量的增加而不断增大，高频能量的信号受衰减减弱也能通过水合物样品，因而声速振幅不断增加。在图5中开始时的振幅出现了短暂的降低，这是受甲烷进气的影响^[8]。

图5 水合物生成过程中振幅随时间的变化关系

声速增加到最大值后趋于一定值，波形的振幅达到最大值后反而开始呈现不断降低的趋势， 80 h后趋于稳定。这是由于受反应釜内游离的甲烷气的影响。THF水合物生成的为Ⅱ型水合物， THF只能占据水合物晶格的大孔穴，而甲烷分子较小，可以填充在水合物晶格中的小孔，这为甲烷气体分子在水合物中的渗透提供了可能；同时，在甲烷气进气过程中，溶液中会有大量的微小气泡，这些微小气泡部分在水合物生成过程中被消耗，部分存在于沉积物孔隙中有助于裂纹形成。生成的大量的水合物虽然能填满沉积物的孔隙，但水合物与沉积物颗粒胶结成岩的过程中，受水合物的体积膨胀作用，水合物样品中也会不可避免产生微小的裂纹，这些裂纹会使顶部的自由气更容易扩散渗透到反应釜的底部。从而使波形信号的振幅在达到最大值后开始呈降低的趋势，这也可以从图2的温度、压力的变化曲线上找到证据，声速在7h时达到最大值，体系的压力从11.2降低到实验结束时的10.8MPa，游离的甲烷气体的影响，使声波信号中的高频信号无法通过，因而波形的振幅降低。

2.3 水合物样品的声速

合成的5个水合物样品的具体参数均已列入表2中（沉积物的粒径、THF水溶液的浓度及测定的声速）。饱和THF水溶液沉积物的声速在1 706～1 782 m/s，生成水合物后样品的声速为3 295～3 984 m/s；对20～40目的样品1、2及40～60目的样品3、4,THF水溶液的浓度为5.9％时生成水合物后样品的声速要比浓度为3.0％的声速值高。对于样品1、3和5比较，看不出沉积物粒径对于声速的影响。

2.4 模型计算

许多学者尝试建立沉积物中水合物饱和度与纵波声速的关系。Dvorkin等^[17]提出了4种可能的水合物分布模型：1）水合物存在于流体中，只改变流体的性质；2）水合物作为部分沉积物的骨架，影响骨架的性质；3）水合物只胶结于沉积物颗粒接触处，水合物作为胶结剂；4）水合物不只在沉积物颗粒接触处胶结，而且水合物还包裹在沉积物颗粒的表面。自然环境中的水合物样品，水合物在沉积物孔隙中的分布仍然是一个颇有争议的问题。加拿大的麦肯齐三角洲（Mallik 2L-38）的测井数据表明，水合物主要分布在沉积物的孔隙中，没有胶结沉积物颗粒^[18]；然而，实验室中的测量结果显示在富含游离气条件下生成的水合物多胶结沉积物颗粒^[8-9]。Priest等^[19]分别测量了利用过量水和过量气的方法合成水合物样品的声速，测量结果表明水合物对声速的影响取决于合成的水合物在沉积物孔隙中存在的状态。实验室测量的结果提供了一种对于一给定地质环境判断水合物分布的模式。

对于水合物分布模式1），生成的水合物悬浮在流体中，水合物只影响流体的体积模量，而没有影响流体的剪切流量，但对沉积物的弹性模量没有贡献；对于水合物分布模型2），生成的水合物作为部分的沉积物骨架，虽然作为沉积物的一种组分，由于水合物的物性参数跟沉积物的岩性相差较大，因而水合物的存在对声速的影响也很微弱；相反，对水合物分布模型3）和4），水合物作为沉积物颗粒间的胶结剂，水合物的胶结作用大大地增加了水合物样品的刚性，因而使水合物样品的声速增加。Dvorkin等^[20]的胶结理论可适用于水合物饱和度较高时的情况，或者沉积物孔隙中完全填满水合物。在本文的模型计算中， THF水溶液的在摩尔比为5.9％时刚好完全生成水合物，因此假设沉积物孔隙中水合物饱和度为100％，含THF水合物沉积物样品的体积模量K和剪切模量G由下式（1）计算^[17]：

南海天然气水合物富集规律与开采基础研究专集

式中：K_h和G_h分别为水合物的体积模量和剪切模量（GPa）；φ_c为松散沉积物的临界孔隙度（φ_c=0.4）；n为单颗粒平均接触数，8.5；参数S_n和S_τ根据Dvorkin和Nur等^[21]提供的计算方法求得：

南海天然气水合物富集规律与开采基础研究专集

南海天然气水合物富集规律与开采基础研究专集

式中：v_q、v_h分别为沉积物及水合物的泊松比。计算方法分别为v_q＝0.5(K_q-2/3G_q)/(K_q+ 1/3G_q） (K_q和G_。分别为石英砂的体积模量和剪切模量）、v_h=0.5(K_h-2/3G_h)/(K_h+1/3G_h），π为常数，3.14；参数α跟水合物的分布相关。对模型3），α＝2[（φ_c－φ）/(3nl－φ_c)]^0.25（φ为孔隙率，当THF水合物完全填满沉积物孔隙时，φ值为0）；对模型4），a=[2（φ_c－φ）/3(1－φ_c)]^0.5。

含水合物沉积物的弹性模量求得后，水合物样品的声速可由式（11）计算：

南海天然气水合物富集规律与开采基础研究专集

其中ρ为水合物样品的体积密度，可以式（12）计算求得：

南海天然气水合物富集规律与开采基础研究专集

模型计算过程中用到的参数列于表3中。在模型计算中，由于THF+CH₄水合物的弹性模量数据无法得到，且该实验合成水合物与纯甲烷水合物差别较小，因此计算中采用纯甲烷水合物的弹性模量参数代替合成的THF水合物。本文只对THF水溶液在浓度为5.9％化学计量比生成的水合物样品2、4及5进行了预测，结果列于表4中。当THF水溶液的浓度小于5.9％时，合成的水合物样品的饱和度很难确定，因而没有对其进行预测。根据表4中实验值与计算值的比较可以看出，胶结理论可以较好的对含水合物沉积物样品进行预测。对样品2和5，实验测量值与模型3）的计算值基本吻合。而对样品4，实验值与模型4）的计算值更接近。同样可以看出，水合物在沉积物中非常有明显的胶结作用。

表3 模型计算中的参数列表

表4 模型计算结果跟实验值对比

3 结论

1）在实验室搭建的水合物声速测定实验装置上，利用不同沉积物粒径和不同的THF水溶液浓度合成了5个的水合物样品，并测量了水合物生成过程中水合物样品的声学性质变化。

2）沉积物孔隙中的填充物对沉积物样品的声速具有重要的影响。对含饱和THF水溶液的沉积物，声速为1 706～1 782 m/s；在THF水合物生成过程中，水合物样品的声速随饱和度的增加而增大，最后趋于一定值；水合物生成反应完全结束后，含水合物沉积物的声速为3 295～3 984 m/s。声速值的大小在很大程度上取决于沉积物孔隙中水合物的饱和度。

3）振幅的变化趋势与声速变化不同，在水合物生成过程中，振幅先随着水合物饱和度的增加而不断增大，达到一个最大值后开始降低，最后趋于稳定，这主要是受游离的甲烷气体在水合物样品中不断渗透的影响结果；同时，水合物胶结模型的计算值与实验值基本一致。

参考文献

[1]Clennell M B,Hovland M,Booth J S,et al.Formation of Natural Gas Hydrates in Marine Sediments 1.Conceptual Model of Gas Hydrate Growth Conditioned by Host Sediment Properties[J].J Geophys Res,1999,104(B10):22985-23003.

[2]Sloan E D.Fundamental Principles and Applications of Natural Gas Hydrates[J].Nature,2003,426(20):353-359.

[3]Klauda J B,Sandler S I.Global Distribution of Methane Hydrate in Ocean Sediment[J].Energy Fuels,2005,19 (2):459-470.

[4]Reagan M T,Moridis G J.Large-scale Simulation of Methane Hydrate Dissociation Along the West Spitsbergen Margin[J].Geophys Res Lett,2009,36:doi:10.1029/2009GL041332.

[5]Pearson C F,Halleck P M,Mc Gulre P L,et al.Natural Gas Hydrate Deposits:A Review of in Situ Properties[J].J Phys Chem,1983,87:4180-4185.

[6]Buffet B A,Zatsepina O Y.Formation of Gas Hydrate from Dissolved Gas in Natural Porous Media[J].Marine Geology,2000,164:69-77.

[7]Carcione J M,Gei D.Gas-Hydrate Concentration Estimated from P-and S-Wave Velocities at the Mallik 2L-38 Research Well,Mackenzie Delta,Canada[J].J Applied Geophys,2004,56:73-78.

[8]Waite W F,Winters W J,Mason D H.Methane Hydrate Formation in Partially Water-Saturated Ottawa Sand[J].Am Mineral,2004,89:1202-1207.

[9]Winters W J,Waite W F,Mason D H,et al.Methane Gas Hydrate Effect on Sediment Acoustic and Strength Properties[J].J Petrol Sci Eng,2007,57:127-135.

[10]Westbrook G K,Chand S,Rossi G,et al.Estimation of Gas Hydrate Concentrationfrom Multi-Component Seismic Data at Sites on the Continental Margins of NW Svalbard and the Storegga Region of Norway[J].Marine Petroleum Geology,2008,25:744-758.

[11]Winters W J,Pecher I A,Waite W F,et al.Physical Properties and Rock Physics Models of Sediment Containing Natural and Laboratory-Formed Methane Gas Hydrate[J].Am Mineral,2004,89:1221-1227.

[12]Pearson C,Murphy J,Hermes R.Acoustic and Resistivity Measurements on Rock Samples Containing Tetrahydrofuran Hydrates:Laboratory Analogues to Natural Gas Hydrate Deposits[J].J Geophys Res,1986,91 (B14):14132-14138.

[13]Sabase Y,Nagashima K.Growth Mode Transition of Tetrahydrofuran Clathrate Hydrate s in the Guest/Host Concentration Boundary layer[J].J Phys Chem B,2009,113:15304-15311.

[14]Priest JA,Best AI,Clayton CR LALaboratory Investiga-tioni nto the Seismic Velocities of Methane Gas Hydrate-Bearing Sand[J].J Geophys Res,2005,110,B0,doi:10.1029/2004JB003259.

[15]Priest J A,Best A I,Clayton C R I.Atenuation of Seismic Waves in Methane Gas Hydrate-Bearing Sand[J].Geophysical Journal International,2006,164:149-159,doi:10.1111/j.1365-246X.2005.02831.x.

[16]Pride S R,Berryman J G,Harris J M.Seismic Atenuation Due to Wave-Inced Flow[J].J Geophys Res,2004,109:B01201,doi:10.1029/2003JB002639.

[17]Dvorkin J,Helgerud M B,Waite W F,et al.Introction to Physical Properties and Elasticity Models[C]//In Mas M D.Natural Gas Hydrate in Oceanic and Permafrost Environments.Kluwer:Dordrecht,Netherlands,2000:245-260

[18]Winters W J,Dallimore S R,Collett T S,et al.Physical Properties of Sedimentsfrom the JAPEX/JNOC/GSC Malik 2L-38 Gas Hydrate Research Well[J].Bull.Geol Surv Can,Rep 1999,544:95-100.

[19]Priest J A,Rees E V L,Clayton C R I.Influence of Gas Hydrate Morphology on the Seismic Velocities of Sands[J].J Geophys Res,2009,114: Bl 1205,doi:10.1029/2009JB006284.

[20]Dvorkin J,Berryman J,Nur A.Elastic Moli of Cemented Sphere Packs[J].Mechanics of Materials,1999,31:461-469

[21]Dvorkin J,Nur A.Elasticity of High-Porosity Sandstones:Theory for Two North Sea Data Sets[J].Geophysic,1996,61:1363-1370.