用分子动力学模拟甲烷水合物热激法分解

用分子动力学模拟方法研究甲烷水合物热激法分解，系统地研究注入340 K液态水的结构Ⅰ型甲烷水合物的分解机理.模拟显示水合物表层水分子与高温液态水分子接触获得热能，分子运动激烈，摆脱水分子间的氢键束缚，笼状结构被破坏.甲烷分子获得热能从笼中挣脱，向外体系扩散.热能通过分子碰撞从外层传递给内层水分子，水合物逐层分解.对比注入277K液态水体系模拟结果，得出热激法促进水合物分解. 关键词： 甲烷水合物 分子动力学模拟 热激法     

显微激光拉曼光谱测定甲烷水合物的水合指数

甲烷水合物是由甲烷气体分子与水分子在低温高压下形成的一种笼型结构化合物,广泛存在于海底陆架区和陆地冻土区,被认为是一种潜在的能源资源。在水合物的晶格中,水分子在氢键的作用下形成大小不同的笼子,甲烷分子可分别进入大笼(51262)和小笼(512)中。在自行研制的实验装置上,分别合成了一系列不同体系下的甲烷水合物,包括十二烷基硫酸钠(SDS)水溶液-甲烷体系、冰粉-甲烷体系以及冰粉-不同粒度砂-甲烷体系。对这些甲烷水合物样品进行了激光拉曼光谱分析,测定了其水合指数,笼占有率等结构参数。结果表明,这些甲烷水合物都为Ⅰ型结构,其水合指数和笼占有率基本不受沉积物粒径大小的影响。在3种体系中生成的水合物,大笼中甲烷分子基本占满,占有率大于97%;小笼中甲烷分子占有率为80%～86%,测得的水合指数为6.05～6.15。     

CH4水合物在SDS和CTAB溶液中的生成特性研究

为探究不同促进剂在甲烷水合物生成过程的微观作用机理,选取动力学促进剂十二烷基硫酸钠(SDS)和热力学促进剂十六烷基三甲基溴化铵(CTAB)作为添加剂,采用分子动力学方法研究其对甲烷水合物生成速率的影响.通过分析势能变化、均方位移、径向分布函数、分子簇生长速率,发现质量分数为0.9%SDS、1.2%SDS、1.2%CTAB、1.6%CTAB的溶液均可促进水合物生成.质量分数为1.2%的SDS溶液水合物生长速率最快,且SDS促进效果优于CTAB.通过分析甲烷分子密度分布云图,发现呈阴性的SDS分子头部基团吸附了大量甲烷分子,水分子受挤压向中间聚集;CTAB含氮的头部基团朝向均相溶液,包含在不稳定的水合物笼中,形成半笼型水合物.相比之下,CTAB溶液中水合物含气率更高.     

热力法开采地层天然气水合物的热动力评价

本文在推导热力作用下水合物地层温度分布基础上,发展了评价热力法开采天然气水合物的热效率(用于水合物分解的热量与输入总的热量之比)和能量效率(即输出能量与输入能量之比)的模型.模型分析表明,水合物地层热物性参数以及水合物饱和度决定了热力法开采的能量效率.在注入蒸汽开采初始条件下,能量效率可以达到7.0.     

SⅡ型R134a水合物稳定性的分子动力学模拟北大核心

SⅡ型R134a制冷剂笼状水合物是理想的蓄冷介质之一。基于分子动力学(MD)方法,建立晶体结构模型,在NPT系综下运算,通过观察最终构象,分析计算体系中粒子的均方位移(MSD)和径向分布函数(RDF),从微观上揭示了小分子和温度对SⅡ型水合物晶体结构稳定性的影响。结果表明:小分子对小晶穴(512)占有率从0.25变化到1.00,占有率为0.25和0.50时,水合物晶体结构分解严重,体系中分子混乱,晶体结构难以稳定存在,占有率在0.75以上时,水合物晶体结构清晰有序,能够稳定存在;当体系温度由260K增大到290K时,水合物结构的稳定性逐渐降低,在280K时晶体结构出现扭曲变形,此时有部分水合物开始分解。     

笼形水合物的科学与技术研究以及在能源和环境领域中的应用

能源与环境是人们越来越关注的问题,而笼形水合物同时在这两方面展示出了巨大的应用前景:(1)在广阔的海底大陆架中蕴藏着大量的以甲烷为主体的笼型水合物--天然气可燃冰;(2)将大气中日益增多的二氧化碳装入笼型水分子中并沉于海底以降减温室效应;(3)氢气可燃冰作为一种清洁的能源载体具有很高的能量密度、可重复再生、可快速允氢;(4)利用笼形水合物的形成过程,可以有效地分离气体、降低能耗、减轻污染,虽然笼形水合物在工程应用以及科学研究中有着重要的价值,仉其晶体结构,成键机制,温压相图,热化学与力学稳定性,合成与分解的反应动力学,声学弹性,与海底沉积的反应,以及扩散和输运性质等郜有待深入的研究.目前,科学家已经研发或者正在探讨将高压和低温装置与中子衍射技术以及激光光谱、热学测量、超声技术等有机的结合在一起,从而能够进行一系列的实验研究米解决诸多的基本科学问题.高压及低温环境下的中子衍射装置在氢气、甲烷以及二氧化碳气体水合物体系研究中显示了巨大的优势,并且在确定水合物晶体结构,气体分子在水合物笼格中的占据情况,以及气体分子在笼格中的分布状态等方面取得了巨大的进展.     

高盐溶液蒸发过程中的分子动力学模拟与研究

随着我国工业的快速发展,如何减量化处理工业排放的高盐废水已经成为了亟待解决的环境问题.在处理高盐废水的各种工艺中,热蒸发技术因具备脱盐效果好,灵活性好等特点得到了广泛应用.本文从微观角度出发,采用分子动力学模拟方法研究了LiCl、KCl、CaCl_(2)三种溶液在400 K和500 K两种温度加热条件下的蒸发过程,分析了Li^(+)、K^(+)、Ca^(2+),Cl^(-)四种盐离子对蒸发速率、水分子取向、氢键、溶液结构等性质的影响.研究结果表明温度的提升会对蒸发速率产生极大影响,温度的提高不会改变配位水化层的位置但会明显减少离子的配位水分子数,有利于提高活跃水分子的占比和蒸发速率.     

天然气水合物快速合成可视化实验研究

水合物法储运天然气技术在安全、高效储运天然气方面极具应用潜力，其中甲烷水合物合成动力学又对该技术的发展至关重要。为此，实验研究了十二烷基硫酸钠(SDS)浓度以及初始压力、初始温度对四氢呋喃(THF)+SDS组合溶液中甲烷水合物合成动力学的影响规律。研究表明，SDS的加入显著地改善了水合物合成速率和最终气体吸收量，但随着SDS浓度增加动力学促进效果逐渐减弱，即存在动力学效果最佳的SDS浓度。此外，SDS不但加速了水合物的生长，而且随SDS浓度增加水合物的生长模式也由优先气相转变为优先液相。最后发现，甲烷水合物最终气体吸收量随初始压力和初始温度升高而增加，合成速率随初始压力升高而增加随初始温度升高而减小。     

添加正辛醇的溴化锂水溶液汽液界面微观形态

采用分子动力学方法研究303 K时,添加正辛醇(1-octanol)的溴化锂水溶液汽液界面的微观形态、密度分布和结构性质.模拟结果表明: Li+、Br在汽液界面处解吸;正辛醇分子吸附在汽液界面并在界面处优势取向;醇分子的亲水基会尽可能多地趋向水分子,二者在界面处以氢键相互作用.同时,在非平衡条件下,采用分子动力学方法模拟了添加或未添加正辛醇的溴化锂水溶液吸收水蒸气的动态过程,发现在短时间内,部分水蒸气分子吸附在汽液界面处或到达溶液内部.     

温度相关的微波频率下氯化钠水溶液介电特性

针对NaCl水溶液(0.001–0.5 mol/L)介电特性, 实验调查了频率(200–6.25 GHz), 温度(293–353 K)、 浓度相关复介电常数. 结果表明: 频率增大的过程中虚部呈逐渐减小的趋势, 高温使离子扰动增大, 破坏了溶液内部水分子四面体结构和氢键构象而使介电常数实部减小. 与纯水相比, 溶液的损耗角正切在高温353 K低频区下降明显. 同时发现2.45和5.8 GHz的复介电常数随温度变化的温度窗效应, 温度窗效应导致微波加热时耗散功率的振荡变化, 温度分布不均匀现象在实验中得以证实. 关键词： 微波 复介电常数 介电特性 高温探头     

显微激光拉曼光谱原位观测甲烷水合物生成与分解的微观过程

采用显微激光拉曼光谱技术对高压透明毛细管中甲烷水合物的生成与分解的微观过程进行了原位观测,初步探讨了甲烷水合物笼型结构的变化规律.结果表明,甲烷水合物在生成过程中,甲烷分子的拉曼峰(2 917 cm-1)逐渐分裂为两个峰(2 905和2 915 cm-1),表明溶解态甲烷分子从单一的化学环境进入了两个有差异的化学环境中...     

Methane combustion in hydrate systems: Water-methane and water-methane-isopropanol

Kinetics of dissociation of synthetic and natural methane gas hydrates, and also double isopropanol-methane hydrate is investigated. Thermal fields of the sample surfaces are measured by means of thermal imaging in combustion of released methane with clathrate dissociation. The dissociation rates of natural hydrate and double hydrate with isopropanol are many times lower than those of synthetic methane hydrate. Methane combustion is accompanied by formation of a thin water film on the powder surface, which has a strong effect on the heat and mass transfer mechanisms. The experiments demonstrated partial self-preservation for methane hydrate and the absence of self-preservation for double isopropanol-methane hydrate. The experimentally observed dissociation rate of double isopropanol-methane hydrate is considerably lower than that of methane hydrate.     

Molecular dynamics simulations of PVP kinetic inhibitor in liquid water and hydrate/liquid water systems

The possible effects of PVP (poly(N-vinylpyrrolidone)) on the properties of liquid and water in clathrate hydrate has been investigated using NVT molecular dynamics simulations. A model for a monomer of the PVP polymer is immersed in three systems, liquid water, a unit cell of a hydrate in liquid water with a hydrate former and a third system where some of the liquid water molecules of this last system are replaced by a PVP monomer. Both molecular dynamics simulation and integral equation theory predict hydrogen bonding between the double bonded oxygen in the PVP ring and hydrogens in water. For the composite system, the PVP monomer has a preference for hydrogen bonding to hydrogens from the water molecules at the surface of the hydrate lattice. The simulations indicate that the PVP monomer tends to orient perpendicular to the hydrate surface. For the model systems in this study PVP may form hydrogen bonds with liquid water through the double bonded oxygen in the ring. When a hydrate crystal is immersed in the liquid water phase this hydrogen bonding is shifted towards the hydrate due to a more favourable Coulomb interaction involving hydrogens from more than one water molecule at the hydrate surface. The PVP monomer has a preference for perpendicular orientation with respect to the hydrate surface. A scheme is suggested for the characterization of kinetic hydrate inhibitors based on molecular dynamics simulations and on three basic properties. In addition to the energy between the active groups of the inhibitor and hydrate water another point of focus is the free energy changes in the interactions between the inhibitor and water as the charges are changed from zero to the original model charges. In particular the difference between this integral for the (hydrate water)–(PVP monomer) interaction and the (liquid water)–(PVP inhibitor) interaction should reflect the driving forces in freezing the inhibitor out from the liquid water phase and onto the hydrate surface. The third property in the characterization scheme is the diffusivities of groups connecting to the hydrate crystal, relative to the diffusivities of the hydrate crystal. Results are presented from simulations where a small cavity with a methane model as a guest is immersed in liquid water with free methane molecules at a temperature of 150K. Changes in structure, diffusivities and energy indicate a tendency towards a more solid–like structurde around the cavity.     

醇类抑制剂对甲烷水合物形成的影响

在本工作中,甲烷水合物的生长动力学是通过甲醇、乙醇、乙二醇三种不同醇类抑制剂存在下的分子动力学模拟研究的.模拟结果发现,三种醇类都可作为甲烷水合物的抑制剂,醇类分子中的亲水性羟基极大地破坏了水合物笼的结构,并且羟基可以与局部的液态水分子形成氢键,从而增加了形成水合物笼型结构的难度,导致甲烷水合物的生长速率降低.对于甲醇分子,甲醇分子的亲水性羟基与水分子形成氢键从而破坏了水分子结构,而亲油性甲基对周围的水分子具有簇效应,两者都会降低水合物生长速率;对于乙二醇和乙醇分子,它们只含有羟基,特别是乙二醇分子含有两个羟基,其对H₂O分子有很强的吸附作用,导致水合物生长速率降低.在抑制效果方面,甲醇分子最优,乙二醇稍微优于乙醇.     

Theoretical investigation of structures and compositions of double neon-methane clathrate hydrates,depending on gas phase composition and pressure

The region of existence of neon clathrate hydrates is an actual problem of hydrate chemistry. The current work presents theoretical study of the equilibrium formation conditions of pure neon clathrate hydrates and double clathrate hydrates of neon-methane mixture. The structures and properties of double clathrate hydrates were described within the scope of the previously developed molecular clathrate hydrate model that takes into account the influence of guest molecules on the host lattice, interaction of guest molecules between themselves, and the possibility of multiple filling of host lattice cages by guest molecules. The model makes it possible to find an equilibrium state and thermodynamic properties of clathrate hydrates at given values of p and T. In the present work, we considered the properties of double clathrate hydrates in the range of pressures from 0 to 4 kbar at 250 K. The results of modeling have shown that the mass fraction of neon in double clathrate hydrate of Ne and CH₄ mixture of cubic structure I (sI) can reach 26%, and 22.5% in double hydrate of cubic structure II (sII) even at a low methane concentration (1%) in gas phase, at high pressure. It is shown that in double clathrate hydrates of the Ne and CH₄ mixture at high pressures, phase transition sII-sI can occur.     

Raman Spectroscopy of Nitrogen Clathrate Hydrates

在253 K和16 MPa的压力下,于实验室内合成了氮气水合物,用显微共焦拉曼光谱对其N-N和O-H键伸缩振动的光谱特征进行了研究．结果表明,氮气水合物中的N-N和O-H键的拉曼峰分别为2322.4和3092.1 cm^-1,与天然的空气水合物中的数据十分接近．另外,还测定了液氮和溶解于水中的氮分子中N-N键的拉曼峰值,分别为2326.6和2325.0 cm^-1．氮气笼型水合物分解的拉曼谱图表明,氮分子同时进入水合物的大笼和小笼中,但由于氮分子在大、小笼中的环境氛围十分接近,其拉曼位移相差不大,故拉曼谱图只能显示N-N键伸缩振动一个峰．     

Influence of gravitational potential on the thermodynamic stability of pure and mixed clathrate hydrates

Clathrate hydrates are an ice-like material consisting of gas molecules confined within cavities in a crystalline water lattice. Phase equilibria of clathrate hydrates systems was described using the statistical mechanical theory of van der Waals and Platteeuw. This theory makes use of the fractional occupancy of cavities within the clathrate hydrate lattice in the determination of chemical equilibria. Classical density functional theory with intermolecular interactions restricted to the first hydration shell was employed to determine the fractional occupancy. In addition to the external field describing the gas-water interactions, the effect of a gravitational field was introduced. The results of the calculations show that although the gravitational potential term may be orders of magnitude smaller than the thermal kinetic energy of the gas species or the hydrogen-bond energy holding the clathrate lattice together, it can nevertheless influence the phase equilibrium of the clathrate hydrate system to some degree. The effect of the magnitudes of both the gravitational potential and the local gravitational field are considered too.     

Structure and dynamics of hydrogen molecules in the novel clathrate hydrate by high pressure neutron diffraction

The D2 clathrate hydrate crystal structure was determined as a function of temperature and pressure by neutron diffraction for the first time. The hydrogen occupancy in the (32+X)H2.136H(2)O, x=0-16 clathrate can be reversibly varied by changing the large (hexakaidecahedral) cage occupancy between two and four molecules, while remaining single occupancy of the small (dodecahedral) cage. Above 130-160 K, the guest D2 molecules were found in the delocalized state, rotating around the centers of the cages. Decrease of temperature results in rotation freezing followed by a complete localization below 50 K.