Rosenbaum Eilis J, English Niall J, Johnson J Karl, Shaw David W, Warzinski Robert P
National Energy Technology Laboratory, United States Department of Energy, PO Box 10940, Pittsburgh, Pennsylvania 15236, USA.
J Phys Chem B. 2007 Nov 22;111(46):13194-205. doi: 10.1021/jp074419o. Epub 2007 Oct 30.
A single-sided transient plane source technique has been used to determine the thermal conductivity and thermal diffusivity of a compacted methane hydrate sample over the temperature range of 261.5-277.4 K and at gas-phase pressures ranging from 3.8 to 14.2 MPa. The average thermal conductivity, 0.68 +/- 0.01 W/(m K), and thermal diffusivity, 2.04 x 10(-7) +/- 0.04 x 10(-7) m2/s, values are, respectively, higher and lower than previously reported values. Equilibrium molecular dynamics (MD) simulations of methane hydrate have also been performed in the NPT ensemble to estimate the thermal conductivity for methane compositions ranging from 80 to 100% of the maximum theoretical occupation, at 276 K and at pressures ranging from 0.1 to 100 MPa. Calculations were performed with three rigid potential models for water, namely, SPC/E, TIP4P-Ew, and TIP4P-FQ, the last of which includes the effects of polarizability. The thermal conductivities predicted from MD simulations were in reasonable agreement with experimental results, ranging from about 0.52 to 0.77 W/(m K) for the different potential models with the polarizable water model giving the best agreement with experiments. The MD simulation method was validated by comparing calculated and experimental thermal conductivity values for ice and liquid water. The simulations were in reasonable agreement with experimental data. The simulations predict a slight increase in the thermal conductivity with decreasing methane occupation of the hydrate cages. The thermal conductivity was found to be essentially independent of pressure in both simulations and experiments. Our experimental and simulation thermal conductivity results provide data to help predict gas hydrate stability in sediments for the purposes of production or estimating methane release into the environment due to gradual warming.
采用单面瞬态平面热源技术,在261.5 - 277.4 K的温度范围内以及3.8至14.2 MPa的气相压力下,测定了压实甲烷水合物样品的热导率和热扩散率。平均热导率为0.68±0.01 W/(m·K),热扩散率为2.04×10⁻⁷±0.04×10⁻⁷ m²/s,其值分别高于和低于先前报道的值。还在NPT系综中对甲烷水合物进行了平衡分子动力学(MD)模拟,以估计在276 K和0.1至100 MPa压力下,甲烷组成从最大理论占有率的80%到100%时的热导率。使用了三种水的刚性势模型进行计算,即SPC/E、TIP4P - Ew和TIP4P - FQ,其中最后一种模型考虑了极化率的影响。MD模拟预测的热导率与实验结果合理吻合,不同势模型的热导率范围约为0.52至0.77 W/(m·K),其中可极化水模型与实验结果吻合最佳。通过比较冰和液态水的计算热导率值与实验热导率值,对MD模拟方法进行了验证。模拟结果与实验数据合理吻合。模拟结果预测,随着水合物笼中甲烷占有率的降低,热导率会略有增加。在模拟和实验中均发现热导率基本与压力无关。我们的实验和模拟热导率结果提供了数据,有助于预测沉积物中气体水合物的稳定性,以便进行开采,或估计由于逐渐变暖导致甲烷释放到环境中的情况。