客体分子所处晶穴结构不同对甲烷水合物稳定性的影响

采用分子动力学方法模拟了SⅠ型甲烷水合物受热分解微观过程,并对水合物分解过程中不同晶穴结构内客体分子对甲烷水合物稳定性的作用进行了研究.通过最终构象、均方位移和势能等性质的变化规律对分别缺失大晶穴和小晶穴中客体分子的2种水合物体系随模拟温度升高稳定性的变化进行了分析.模拟结果显示,随温度的上升,水合物稳定性逐渐下降直至彻底分解;而水合物分解速度与2种晶穴各自部分晶穴占有率相关,不能简单的通过整体晶穴占有率表示.对比相同注热过程中2种水合物体系分解状况,发现位于大晶穴内的客体分子对水合物稳定性影响更大,缺失大晶穴内客体分子的水合物更容易随温度升高而分解.     

微波场中乙烷及丙烷气体水合物的分解特性

研究了2.45GHz微波场中I型乙烷水合物及II型丙烷水合物的热激分解过程,基于晶体表面两步分解机制的动力学模型,结合传热传质分析了其分解特性.结果表明:水合物在微波场中的加热分解是一个与实际微波电磁场相互耦合的过程,微波体积加热的特点强化了水合物颗粒表层的传热传质过程,时间累积的热效应增大了水合物晶体破解速率;在120至540W入射功率下,乙烷、丙烷水合物气化速率分别达到0.109-0.400mol·min-·1L-1及0.090-0.222mol·min-1·L-1.在一定范围内增大微波功率可显著提高水合物分解速率,其中乙烷水合物一直处于功率主控区,丙烷水合物更早进入功率和分解动力机制共同控制区.     

13C固体核磁共振测定气体水合物结构实验研究

采用高功率1H去偶结合魔角旋转13C固体核磁共振技术,在低温常压条件下对合成的乙烷和丙烷气体水合物进行了测试,获得了两种纯气体水合物的13C核磁共振谱图,初步建立了固体核磁共振波谱法测定天然气水合物的实验方法.实验表明:乙烷水合物的13C核磁共振谱图中仅有一条谱线(δ7.7),结构类型为sI,且乙烷分子仅填充在大笼中(...     

气体分子对甲烷水合物稳定性的影响

通过B3LYP方法, 在6-31G(d,p)水平下, 分别优化了结构I型甲烷水合物十二面体和十四面体晶穴结构. 结果表明, CH4分子使晶穴的相互作用能降低, 增强了晶穴的稳定性. 计算了晶穴中甲烷分子C—H键的对称伸缩振动频率, 计算结果与实验值相符合. 研究发现CH4分子影响晶穴中氧原子的电荷分布, 从而增强了氢键的稳定性. 通过分子动力学方法研究水合物晶胞中气体分子的占有率对水合物稳定性的影响, 进一步说明气体分子对水合物晶穴稳定性的重要作用.     

过氧化氢水溶液作用下甲烷水合物分解特性的分子动力学模拟研究

采用分子动力学模拟方法研究过氧化氢水(HP)溶液作用下结构I型(SI)甲烷水合物晶体分解特性. 系统分析甲烷水合物在过氧化氢水溶液作用下由晶态向液态转变过程的机理, 对比相同摩尔浓度乙二醇(EG)溶液作用下甲烷水合物分解变化规律, 得出HP与水合物热力学抑制剂EG一样对甲烷水合物分解具有促进作用, 为HP溶液促进甲烷水合物分解实验研究提供参考.     

南海神狐海域及祁连山冻土区天然气水合物的拉曼光谱特征

采用显微激光拉曼光谱对我国在南海神狐海域及祁连山冻土区首次钻获的天然气水合物实物样品进行了详细的研究, 探讨了其笼型结构特征及其气体组成. 结果表明, 南海神狐海域天然气水合物样品是典型的I型结构(sI)水合物, 气体组分主要是甲烷, 占99%以上; 水合物大笼的甲烷占有率大于99%, 小笼为86%, 水合指数为5.99. 祁连山冻土区天然气水合物气体组分相对复杂, 主要成分除甲烷外(70%左右), 还有相当数量的乙烷、丙烷及丁烷等烃类气体, 从拉曼谱图上可初步判断其为II型结构(sII)水合物; 水合物的小、大笼的甲烷占有率的比值(θ_S/θ_L)为26.38, 远远大于南海神弧海域水合物的0.87, 这主要是由于祁连山水合物气体组分中的大分子(乙烷、丙烷及丁烷等)优先占据水合物的大笼, 大大减少了大笼中甲烷分子的数量.     

天然气水合物抑制剂的制备及抑制性能评价

采用溶液聚合法合成了聚乙烯己内酰胺PVCap及共聚物,利用FT-IR和~1H-NMR表征了它们的结构。将动力学抑制剂PVCap、乙烯基己内酰胺-乙二醇共聚物P(VCL-APEG)、乙烯基己内酰胺一含酯单体共聚物P(VCL-A)、聚乙烯吡咯烷酮PVP分别与甲醇进行复配,通过THF方法,研究不同抑制剂对天然气水合物的抑制性能,探讨了动力学抑制剂添加量、复合抑制剂中甲醇添加量等对抑制性能的影响。实验结果表明,甲醇的加入提高了动力学抑制剂的抑制能力,当1%的PVCap、PVP、P(VCL-APEG)、P(VCL-A)分别与5%甲醇复配时,其诱导时间分别可达887min、639min、420min、3300min。     

短链烷烃二元混合物在分子筛上吸附分离的分子模拟

用巨正则是系综Monte Carlo(GCMC)与构型偏倚(CBMC)相结合的方法模拟了 MFI分子筛对甲烷-丙烷、乙烷-丙烷体系(300K，345kPa)的吸附平衡，模拟结果与 文献实验结果相吻合，分别模拟了FER，ISV，MEL，MFI，MOR，TON等六种分子筛对 甲烷－丙烷、乙烷－丙烷体系（300K，345kPa）的吸附，得出甲烷－丙烷体系中分 子筛对较长链烷烃的选择性大小顺序（气相乙烷摩尔分数为0.5时）为ISV＞MEI＞ MEL＞FER＞TON＞MOR,对乙烷－丙烷体系选择性大小顺序（气相乙烷摩尔分数为0. 5时）为ISV＞MOR＞MFI＞FER＞MEL＞TON. MOR型分子筛对两个不同体系的吸附行为 表现出明显的不同，两个体中ISV的吸附量均最大，MFI，MEL，FER次之，此三种分 子筛具有相拟的吸附量，MOR和TON型分子筛吸附量较低。     

烷烃混合物在Cu-BTC中的吸附与分离

用巨正则系综Monte Carlo (GCMC)和构型导向Monte Carlo (CBMC)相结合的方法模拟了298 K下甲烷-乙烷-丙烷体系以及正丁烷-异丁烷体系在1,3,5-苯三甲酸铜(II) (Cu-BTC)中的吸附行为. 结果表明, Cu-BTC对丙烷以及异丁烷的吸附分离都有较好的选择性. 通过我们发展的“材料剖面成像”方法研究了烷烃混合物在Cu-BTC中不同压力下的吸附位点, 从而进一步分析了烷烃混合物在Cu-BTC中的分离性能. 结果发现, 在吸附过程中主要存在着两种效应, 即能量效应和尺寸效应的竞争. 在甲烷-乙烷-丙烷体系中, 较高压力下, 由于尺寸效应的影响, 丙烷主要吸附在主孔道中, 而对甲烷和乙烷组分, 能量效应占主导地位, 从而导致乙烷主要吸附在四面体孔内, 甲烷则主要吸附在三角形孔窗外. 在正丁烷-异丁烷体系中, 能量效应起主导作用, 从而使异丁烷主要吸附在四面体孔内, 而正丁烷主要吸附在主孔道中.     

用分子动力学模拟天然气水合物的抑制效应

用分子动力学模拟方法确定了结构H型(SH)天然气水合物的稳定晶体生长面为(001), 系统研究了277 K时三种动力学抑制剂对此晶面的影响. 模拟显示抑制剂中的氧与表面水分子形成氢键, 从而破坏原有的稳定结构, 造成水合物笼型结构坍塌, 达到抑制水合物形成的效果. 比较三种不同动力学抑制剂对SH的抑制效果得出: PVCap＞PEO＞PVP. 在此基础上研究了PVCap对天然气水合物结构I型(SI), 结构II型(SII)和SH三种不同晶型的抑制效应. 模拟发现抑制效果的次序为: SH＞SI＞SII.     

Kinetic studies of gas hydrate formation with low-dosage hydrate inhibitors

Pipeline blockage by gas hydrates is a serious problem in the petroleum industry. Low-dosage inhibitors have been developed for its cost-effective and environmentally acceptable characteristics. In a 1.072-L reactor with methane, ethane and propane gas mixture under the pressure of about 8.5 MPa at 4 °C, hydrate formation was investigated with low-dosage hydrate inhibitors PVP and GHI1, the change of the compressibility factor and gas composition in the gas phase was analyzed, the gas contents in hydrates were compared with PVP and GHI1 added, and the inhibition mechanism of GHI1 was discussed. The results show that PVP and GHI1 could effectively inhibit the growth of gas hydrates but not nucleation. Under the experimental condition with PVP added, methane and ethane occupied the small cavities of the hydrate crystal unit and the ability of ethane entering into hydrate cavities was weaker than that of methane. GHI1 could effectively inhibit molecules which could more readily form hydrates. The ether and hydroxy group of diethylene glycol monobutyl ether have the responsibility for stronger inhibition ability of GHI1 than PVP.     

Effective control of gas hydrate dissociation above the melting point of ice

Direct measurements of the dissociation behaviors of pure methane and ethane hydrates trapped in sintered tetrahydrofuran hydrate through a temperature ramping method showed that the tetrahydrofuran hydrate controls dissociation of the gas hydrates under thermodynamic instability at temperatures above the melting point of ice.     

Gas hydrate nucleation and cage formation at a water/methane interface

Nucleation of gas hydrates remains a poorly understood phenomenon, despite its importance as a critical step in understanding the performance and mode of action of low dosage hydrate inhibitors. We present here a detailed analysis of the structural and mechanistic processes by which gas hydrates nucleate in a molecular dynamics simulation of dissolved methane at a methane/water interface. It was found that hydrate initially nucleates into a phase consistent with none of the common bulk crystal structures, but containing structural units of all of them. The process of water cage formation has been found to correlate strongly with the collective arrangement of methane molecules.     

Mechanisms for thermal conduction in hydrogen hydrate

Extensive equilibrium molecular dynamics simulations have been performed to investigate thermal conduction mechanisms via the Green-Kubo approach for (type II) hydrogen hydrate, at 0.05 kbar and between 30 and 250 K, for both lightly filled H(2) hydrates (1s4l) and for more densely filled H(2) systems (2s4l), in which four H(2) molecules are present in the large cavities, with respective single- and double-occupation of the small cages. The TIP4P water model was used in conjunction with a fully atomistic hydrogen potential along with long-range Ewald electrostatics. It was found that substantially less damping in guest-host energy transfer is present in hydrogen hydrate as is observed in common type I clathrates (e.g., methane hydrate), but more akin in to previous results for type II and H methane hydrate polymorphs. This gives rise to larger thermal conductivities relative to common type I hydrates, and also larger than type II and H methane hydrate polymorphs, and a more crystal-like temperature dependence of the thermal conductivity.     

On the thermodynamic stability and structural transition of clathrate hydrates

Gas mixtures of methane and ethane form structure II clathrate hydrates despite the fact that each of pure methane and pure ethane gases forms the structure I hydrate. Optimization of the interaction potential parameters for methane and ethane is attempted so as to reproduce the dissociation pressures of each simple hydrate containing either methane or ethane alone. An account for the structural transitions between type I and type II hydrates upon changing the mole fraction of the gas mixture is given on the basis of the van der Waals and Platteeuw theory with these optimized potentials. Cage occupancies of the two kinds of hydrates are also calculated as functions of the mole fraction at the dissociation pressure and at a fixed pressure well above the dissociation pressure.     

Kinetic inhibitor of hydrate crystallization

We present the results of a combined theoretical/experimental study into a new class of kinetic inhibitor of gas hydrate formation. The inhibitors are based on quaternary ammonium zwitterions, and were identified from a computational screen. Molecular dynamics simulations were used to characterize the effect of the inhibitor on the interface between a type II hydrate and natural gas. These simulations show that the inhibitor is bifunctional, with the hydrophobic end being compatible with the water structure present at the hydrate interface, while the negatively charged functional group promotes a long ranged water structure that is inconsistent with the hydrate phase; the sulfonate-induced structure was found to propagate strongly over several solvation shells. The compound was subsequently synthesized and used in an experimental study of both THF and ethane hydrate formation, and was shown to have an activity that was comparable with an existing commercial kinetic inhibitor: PVP.     

常见客体分子对笼型水合物晶格常数的影响

Natural gas hydrates are considered as ideal alternative energy resources for the future, and the relevant basic and applied research has become more attractive in recent years. The influence of guest molecules on the hydrate crystal lattice parameters is of great significances to the understanding of hydrate structural characteristics, hydrate formation/decomposition mechanisms, and phase stability behaviors. In this study, we test a series of artificial hydrate samples containing different guest molecules (e.g. methane, ethane, propane, iso-butane, carbon dioxide, tetrahydrofuran, methane + 2, 2-dimethylbutane, and methane + methyl cyclohexane) by a low-temperature powder X-ray diffraction (PXRD). Results show that PXRD effectively elucidates structural characteristics of the natural gas hydrate samples, including crystal lattice parameters and structure types. The relationships between guest molecule sizes and crystal lattice parameters reveal that different guest molecules have different controlling behaviors on the hydrate types and crystal lattice constants. First, a positive correlation between the lattice constants and the van der Waals diameters of homologous hydrocarbon gases was observed in the single-guest-component hydrates. Small hydrocarbon homologous gases, such as methane and ethane, tended to form sI hydrates, whereas relatively larger molecules, such as propane and iso-butane, generated sⅡ hydrates. The hydrate crystal lattice constants increased with increasing guest molecule size. The types of hydrates composed of oxygen-containing guest molecules (such as CO₂ and THF) were also controlled by the van der Waals diameters. However, no positive correlation between the lattice constants and the van der Waals diameters of guest molecules in hydrocarbon hydrates was observed for CO₂ hydrate and THF hydrate, probably due to the special interactions between the guest oxygen atoms and hydrate "cages". Furthermore, the influences of the macromolecules and auxiliary small molecules on the lengths of the different crystal axes of the sH hydrates showed inverse trends. Compared to the methane + 2, 2-dimethylbutane hydrate sample, the length of the a-axis direction of the methane + methyl cyclohexane hydrate sample was slightly smaller, whereas the length of the c-axis direction was slightly longer. The crystal a-axis length of the sH hydrate sample formed with nitrogen molecules was slightly longer, whereas the c-axis was shorter than that of the methane + 2, 2-dimethylbutane hydrate sample at the same temperature.     

Poly(styrene/pentafluorostyrene)-block-poly(vinyl alcohol/vinylpyrrolidone) amphiphilic block copolymers for kinetic gas hydrate inhibitors: Synthesis,micellization behavior,and methane hydrate kinetic inhibition

Amphiphilic block copolymers of short poly(styrene) (PS) or poly(2,3,4,5,6-pentafluorostyrene) (PPFS) segments with comparatively longer poly(vinyl acetate) or poly(vinylpyrrolidone) (PVP) segments are synthesized using a 2-cyanopropan-2-yl N-methyl-N-(pyridin-4-yl)dithiocarbamate switchable reversible addition–fragmentation chain transfer (RAFT) agent toward application as kinetic gas hydrate inhibitors (KHIs). Polymerization conditions are optimized to provide water-soluble block copolymers by first polymerizing more activated monomers such as S and PFS to form a defined macro chain-transfer agent (linear degree of polymerization with conversion, comparatively low dispersity) followed by chain extensions with less activated monomers VAc or VP by switching to the deprotonated form of the RAFT agent. The critical micelle concentrations of these amphiphilic block copolymers (after VAc unit hydrolysis to vinyl alcohol units) are measured using zeta surface potential measurements to estimate physical behavior once mixed with the hydrates. A PS-poly(vinyl alcohol) block copolymer improved inhibition to 49% compared to the pure methane–water system with no KHIs. This inhibition was further reduced by 27% by substituting the PS with a more hydrophobic PPFS. A block copolymer of PS–PVP exhibited 20% greater inhibition than the PVP homopolymer and substituting PS with a more hydrophobic PPFS resulted in a 35% further decreased in methane KHI. © 2018 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2018 , 56, 2445–2457, 56, 2445–2457     

P-T stability conditions of methane hydrate in sediment from South China Sea

For reasonable assessment and safe exploitation of marine gas hydrate resource, it is important to determine the stability conditions of gas hydrates in marine sediment. In this paper, the seafloor water sample and sediment sample (saturated with pore water) from Shenhu Area of South China Sea were used to synthesize methane hydrates, and the stability conditions of methane hydrates were investigated by multi-step heating dissociation method. Preliminary experimental results show that the dissociation temperature of methane hydrate both in seafloor water and marine sediment, under any given pressure, is depressed by approximately -1.4 K relative to the pure water system. This phenomenon indicates that hydrate stability in marine sediment is mainly affected by pore water ions.     

Activation energy of methyl radical decay in methane hydrate

The thermal stability of gamma-ray-induced methyl radicals in methane hydrate was studied using the ESR method at atmospheric pressure and 210-260 K. The methyl radical decay proceeded with the second-order reaction, and ethane molecules were generated from the dimerization process. The methyl radical decay proceeds by two different temperature-dependent processes, that is, the respective activation energies of these processes are 20.0 +/- 1.6 kJ/mol for the lower temperature region of 210-230 K and 54.8 +/- 5.7 kJ/mol for the higher temperature region of 235-260 K. The former agrees well with the enthalpy change of methane hydrate dissociation into ice and gaseous methane, while the latter agrees well with the enthalpy change into liquid water and gaseous methane. The present findings reveal that methane hydrates dissociate into liquid (supercooled) water and gaseous methane in the temperature range of 235-260 K.