Introduction
ガスハイドレートは、ケージ1のゲスト分子を含む水素結合水分子のケージ構造を含む結晶性化合物です。海底や永久凍土地域におけるメタンハイドレート（MHS）、大量の興味深い将来のエネルギー資源であるが、地球規模の気候条件2に影響を与えることがあります。
2013年3月には、日本の石油、ガス、および金属鉱物資源機構（株）は、「減圧法」3,4を用いて、東部南海トラフにおける天然MH-ベアリング堆積物からガスを抽出するために、世界初の海洋産出試験を実施しました。
ガスハイドレートは、メタン1、水素5、CO 2 1,6、およびオゾン7のようなガスを保存することができます。したがって、メタンと水素水和物は潜在的なエネルギー貯蔵や輸送媒体として検討されています。大気中に放出されたCO 2排出量、CO 2 sequesを低減するために、深海堆積物中のCO 2ハイドレートが6を検討されている使用してtration。オゾンは、現在、浄水、食品殺菌に使用されます。それは7化学的に不安定であるため、オゾン保全技術の研究が行われています。水和物中のオゾン濃度は、オゾン水や氷7よりもはるかに高いです。
自然MH-ベアリング堆積物や水和物ベースの技術からのガス生産を開発するためには、ガスハイドレートの熱特性を理解することが不可欠です。しかし、ガスハイドレート堆積物の熱特性データとモデル研究は、8が不足しています。
「減圧法は、「水和物の安定性以下の細孔の圧力を減少させることによって堆積物の細孔空間内MHを解離するために使用することができます。このプロセスでは、堆積物の孔空間成分は、水、MH、およびガスに水から、およびMH変わります。熱的特性」の測定MHの融解熱は、測定に影響を与える可能性があるため、後者の条件では困難です。この問題を解決するために、村岡らは MH形成9の間に過冷却状態での熱的特性」の測定を行いました。
このビデオプロトコルでは、我々は過冷却合成砂 - 水 - ガス-MHサンプルの測定方法を説明します。
図1は、人工メタンハイドレート堆積物の熱特性を測定するための実験を示しています。参照9に示すように、セットアップは同じです。システムは、主に、高圧容器、圧力及び温度制御、及び測定システムの熱特性を含みます。高圧容器を140 [mm]の内径および140ミリメートルの高さを有する円筒形ステンレス鋼で構成されています。削除デッドボリュームとその内部容積は2110センチメートル3であり、その圧力限界は15 MPaです。 transie NT面光源（TPS）技術は、熱特性10を測定するために使用されます。 2.001ミリメートルの個々の半径を有するナインTPSプローブは、容器の内側に配置されています。 9プローブ9のレイアウトは、文献9の図2に示されています。 TPSプローブは、ケーブルで熱特性」アナライザに接続し、実験中に手動で切り替えられます。容器内のTPSセンサ接続図、および設定の詳細は、図S1、2、および基準9における支援情報の3に示されています。
図1：人工メタンハイドレート堆積物の熱特性を測定するための実験装置の図を参照9から変更されています。3956fig1large.jpg "ターゲット=" _空白 ">この図の拡大版をご覧になるにはこちらをクリックしてください。
TPSの方法は、各サンプルの熱特性を測定しました。この方法の原理は文献10に記載されています。この方法では、時間依存性の温度上昇は、ΔTaveが 、あります
どこ
0は、センサからの出力電力であり、Wは 、式1において、Rはセンサプローブの半径であり、λは、試料の熱伝導率であり、αは熱拡散率であり、tは、電力供給を開始してからの時間でありますセンサプローブに。D（τ）は無次元時間依存する関数です。τ EM>（αT/ R）1/2で与えられます。式2において、mは TPSプローブの同心リングの数であり、I 0は、修正されたベッセル関数です。試料の熱伝導率、熱拡散率、比熱を同時に電力がセンサプローブに供給されるように、温度上昇に印加される反転分析によって決定されます。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
注：この研究は、高圧可燃性のメタンガスと大高圧容器を使用するように、関連するすべての物質安全データシートを参照してください。ヘルメット、安全メガネ、安全靴を着用してください。温度制御システムが停止した場合、容器内の圧力は、MHの分解と共に増加します。事故を防止するために、安全弁システムの使用が強く、自動的に大気中にメタンガスを解放することをお勧めします。安全バルブシステムは、電力供給なしで動作することができます。
砂-水-メタンガスサンプル9の調製
- 振動台上に高圧容器を置きます。
- 砂ボトルに水のボトルと4000グラムのシリカ砂の中に、純水1.5 Lを注ぎます。正確にはそれぞれ、砂と水のボトルに砂と水の質量を量ります。
- 水のボトルから2,110 cm 3程度の容積を有する高圧容器内の純水1Lに注ぎ水が半分内側容器を満たすまで。
- 容器全体を振動させる振動台をオンにします。それぞれ、50 Hzから220 Wに振動速度と電源を設定します。ステップ1.5が完了するまで振動を適用します。容器を振動させることにより容器の底にドレイン線と焼結金属フィルター内の残留空気を除去します。
- 約1g秒の一定の速度で容器に砂ボトルから3,300 gシリカ砂を流動-1容器全体を均一な充填を確実にするために振動させながら水面付近に保持された漏斗を使用して。
- 水は容器の縁に到達したときの振動を停止します。
- こぼれるから水を防ぐために、容器の縁に一時的な壁としてリングを置きます。
- 50 Hzから220 W.で再び容器を振動させます
- 砂が容器（高さ140ミリメートル）の縁に到達すると、振動をオフにします。
- 目を使用して、一時的な壁と過剰間隙水を除去電子ドレンライン。バックウォーターボトルに過剰間隙水を注ぎます。
- 50 Hzから1秒間300 Wで一度か二度、容器を振動させることにより砂を梱包し、必要に応じてより多くの砂を追加します。
- 砂と水のボトルに砂と水の質量を秤量します。砂と水のボトルでの質量差から容器に砂と水の質量を計算します。この実験では、容器内の砂と水の質量はそれぞれ、3385グラムと823.6グラムでした。容器内の水の質量を合計 Wとして示されます。
- ステンレス製の蓋付きの高圧容器を覆い、順番に対角線上の対のボルトを締めます。
- 実験のために意図されたテーブルに振動台から高圧容器を移動します。
- 温度を制御するための断熱材を有する高圧容器を覆います。
- 高圧力容器に高圧配管と冷却水の流れラインを接続します。 入力と出力のガスパイプラインのバルブを開きます。分800ミリリットルの割合で10 Lのメタンを換気-1余分な水は、大気圧下で罠に排出しなくなるまで。砂放電容器の底部に固定された焼結金属フィルタによって阻止されます。親水性シリカ砂は、水分子を吸着するので、残留水は、砂表面上に残ります。
- 容器内のガス量を決定するために、 トラップ wは 、トラップ中の水の質量を量ります。 トラップ ワット -残留水の質量、 ワット 解像度 、 総 ワット のres = ワット方程式を用いて容器内を決定します。この場合には、 ワット 解像度とワット トラップはそれぞれ360.6グラムと463.0グラムでした。
- V 砂はトンの容積であるV 砂 / V セル 、 -式Ѱ= 1を使用してサンプルの空隙率を決定します砂の密度に砂の質量の比によって決定される、彼は砂（ すなわち 、ρsは 2,630キロメートル-3 =）、およびV 細胞が容器の内部容積です。サンプルの空隙率Ѱは0.39でした。
- 出力ガスラインのバルブを閉じます。室温（ すなわち 、31.6℃）で約12.1メガパスカルに容器内のメタンの間隙水圧を高めるために、メタンを注入します。
- 入力ガスラインのバルブを閉じます。
- データロガーを使用して実験中に容器内の圧力と温度の記録を開始します。データのサンプリング間隔は5秒です。総実験時間は約3000分です。
2. MH合成と過冷却サンプル9の熱的性質」の測定
- 冷却液を循環させることにより、室温から2.0℃に容器を冷却するための冷却装置の電源をオンにします。チラートンからクーラント循環をしてみましょうそこから容器の蓋に、容器の底部O、および最終的に戻ってチラーへ。容器内の温度変化率は約0.001°C秒-1でした。
- TPSアナライザソフトウェアを使用して、測定パラメータを設定します。デザイン＃7577をセンサにセンサタイプを設定します。 5秒にmWの30への出力電力W 0と測定時間を設定します。センサの種類やサンプルの状態が変化した場合に適切なパラメータを変更する必要があることに注意してください。 1.5℃まで1℃から温度を上昇させるためにパラメータを設定します。
- 次の式で、SUPΔT、過冷却度を計算します。
ΔTSUP = T EQ（P） - T。 （3）
T 当量 （P）は、圧力Pの関数としてMHの平衡温度でありますTの 当量 （P）CSMGemソフトウェア1を使用して計算されます。0、PおよびTはそれぞれ、圧力および温度計を用いて測定し、容器内の圧力および温度です。 - ΔTの SUPが 2℃以上であるの後に同時にTPS分析器を用いて熱伝導率、熱拡散率、及び容積比熱を測定します。
- 各測定後の熱特性アナライザに接続されたTPSプローブを切り替えます。実験9時に手動TPSプローブとアナライザの間にケーブルを切り替えます。接続図を参照し9の図S2に示します。各センサ用のスイッチングシーケンスはnoです。 6→2→7→5→1→9→4→3→8→6 ...。シーケンスは、測定に影響を与えるから残留熱を防ぐために可能な限り設定されているセンサとの間の距離に基づいています。データごとに3-5分を収集します。
- ＃＆まで測定を繰り返し916; T（商標） 2°Cが再び達します。この実験では、ΔTの SUPは、最初は時間と共に増加します。 ΔTの SUPが最大値に達した後、圧力がMHの形成を減少するため、ΔTの SUPを徐々に0℃に低下します。 ΔTの SUPは、前式（3）を使用して、TPS測定に2°C以上であるかどうかをチェックします。
- 温度プロファイルは、MHの溶融によって影響されないことを確認してください。 MHは、測定中に溶けた場合はMHの溶融は吸熱反応であるため、温度が増加することはありません。測定中の温度プロファイルを確認し、結果のセクションで説明されています。
- TPSの技術を使用して、すべての温度プロファイルデータのための熱的特性」解析を行います。
サンプル9,11の彩度を変更する3.計算
注意：時間tの関数として試料中MHに対する飽和、水、及びガスの程度は、気体の状態方程式を用いて計算されます。使用される計算の詳細と方程式は、以前11に記載されています。
- 時刻tにおけるメタンガスの体積V ガス 、Tを計算します
Qは、容器内のガスの初期体積であり、V MH、T - 1は、時刻tにおけるMHの体積である- 1、及びR VHW水とMHの体積比です。
式5、nはMHの水和数（〜6）、ρMHであり、ρ 水それぞれMHと水の密度に対応し、 ワット MHとワット 水は MHと水、Rの分子量を表しますespectively。 - トンに1 - トンから形成されたMHの量ΔMtの （モル）を計算します
Rは気体定数であり、Pは、メタンガスの圧力であり、 およびZ T（T ガス 、T、P ガス 、t）は時刻tにおけるメタンの圧縮係数です。 Z tの 12、13を計算するため、リーとケスラーによって変更された私たち9と坂本ら 11は 、ベネディクト・ウェッブ・ルビン（BWR）の式を使用しています。この計算には、式（3から7.1） -参照13の7 - （3から7.4）BWR式13とリー・ケスラー定数のが表3で使用されています。 - 体積変化ΔV MHを計算し 、
ここで、P sの 101325 Paでの基準圧力であり、T sは 273.15 Kの基準温度であり、Z Sは P sで圧縮係数であり、T S（Z 秒 〜1）、およびV CH4はメタンガス体積の比でありますMHの単位体積[Nmで3メートル -3]。 165.99 [Nmで3メートル -3]のV CH4値を使用します。 - 容積V MH、時刻tにおけるMHのトンを計算します
- 水Vの 水の量を計算し、時刻tでの圧力容器でT
ここで、V 水は、1水の初期量です。 - 式を使用して計算を繰り返します。時刻tにおける4-9 = 2、3、...は、水、メタン、およびMH 11の彩度の変化を決定します。初期条件は、= 1、 すなわち 、V ガス 、1 = Q Tです。時刻tにおけるPおよびTは、データログ9から取り出されます。計算結果は、次のセクションで示されています。
ここで、V 水は、1水の初期量です。 Subscription Required. Please recommend JoVE to your librarian.
Representative Results
図2aは、MHの溶融に影響されない温度プロファイルを示している。ΔTcは 、熱定数'測定に温度変化である。 図2bは、MHの溶融に影響された温度プロファイルを示しています。これらの式は、安定したサンプルの条件を仮定することによって導出されるので、 図2b中のプロファイルは式1および2を介して分析することができません。
図3aは、圧力、温度、及び時間の関数として容器に過冷却の程度を示します。システムは、圧力と温度平衡に達した後にMHが核形成します。 MHの形成は、時間t = 170分での急激な圧力変化によって特徴付けられます。双頭矢印は、過冷却度が2℃以上であることを示しています。熱定数は、この範囲内で測定しました。
図4は、熱定数」の測定値の一例を示しています。実験条件は、P = 825分であったT= 7.1 MPaで、T = 2.4°C、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及び4cは 、各解析範囲内のTPSの技術を用いて得た。 図4Bは、ΔT の AVEと温度変化ΔTの AVE（τ）およびD（τ）を示します（τ ）=ΔTC（T）。 ΔTaveとの関係（＆＃964;）およびD（τ）は、解析範囲によって異なり、図4Cは、時間tの平方根対温度T dを示しています。 TPS反転分析によって得られた線形フィットからの温度データの偏差は、T Dです。 TPSセンサプローブの絶縁層は、測定に影響を与えることを示唆している図4cに示すように、測定の開始時に分析1の偏差は、かなり大きいです。
前述したように、表1は、各解析範囲内の熱定数を示します。特徴的な時間比率の合計は、特徴時間τで割った総分析時間（t = 2-4秒間、合計時間は4秒である）で定義されます。 TPSの技術を使用した場合の特性の時間割合の合計が1未満であるべきであることに留意されたいです。このことは、文献に記載されています
測定値に影響を与えるセンサプローブを回避するために、各測定の開始時にデータが使用されるべきではありません。分析時間範囲を調整することにより、 図4cに示すように、T dの平均偏差が最小化されます。特性時間比合計分析時間範囲を調整することにより、統一に調整されます。したがって、我々は分析2ない1からの熱定数の値を採用しました。
熱伝導率、比熱、熱拡散率は、それぞれ、 図5A、B、およびCにおける時間の関数として示されています。最後に、我々は、Tの結果を要約しますhermalプロパティおよびハイドレート飽和。その結果に関する詳細は章に記載されています。 refference 9の4。
図2：時間の関数としての温度プロファイル（a）は、（過冷却状態）を溶融MHの影響を受けないと、（b）MH溶融によって影響を受ける。 この図の拡大版をご覧になるにはこちらをクリックしてください。
両方の温度プロファイルは、予備実験からのものであることに注意してください。測定時間は、融解熱の影響を明らかにするために実験でよりも長いです。予備実験では、計測時間tが 40秒であり、出力電力W 0は、20 MW（a）および50 MW（B）でした。
<キープtogether.withinページ= "1">：FO Pクラス= "jove_content"
図3（a） の圧力、温度、及び時間の関数として容器に過冷却度。双頭矢印は、過冷却度が2℃以上であることを示しています。熱定数は、この範囲内で測定しました。 （b）は MH、水、およびサンプルのメタンガス飽和が（参照9から再版）時間の関数として示されている。 この図の拡大版をご覧になるにはこちらをクリックしてください。
図4： 熱定数'測定の分析例 （a）の温度プロファイル使用して。TPS測定方法。時間は分析1の範囲であり、2は、それぞれ0~5秒、0.65から4.88秒、です。 （b）は ΔTaveの持つ温度変化ΔT アベニュー （τ）とD（τ）との関係（τ）=ΔTC（t）とします 。 （c）の平方根時間tと温度T dを。 この図の拡大版をご覧になるにはこちらをクリックしてください。
図5：時間の関数としての（a）の熱伝導率λ、時間の関数としての（b）の比熱ρCpを 、時間の関数としての（c）の熱拡散率α。結果はMH飽和の関数としての熱特性に変換しました。変換結果と関連する議論は文献で報告されています。 9.データは、範囲トン = 210から980分以内にオーバーラップを示しています。明確にするために、プロットされたデータは、この範囲内の同一のセンサからの3回の測定の平均を表します。これらの数字は、参照9から変更されている。 この図の拡大版をご覧になるにはこちらをクリックしてください。
| 分析範囲の | λ、Wのメートル- 1 K - 1 | ρCpを 、MJメートル- 3 K - 1 | α、ミリメートル2秒- 1 | Chの総AR。時間 | デベロッパーを意味する。、°C | |
| 分析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の生成熱の影響を推定しました。 図3b及び形成H = 52.9 kJのモル-1 MH 14のためのエンタルピーに示すように、MHの形成熱はS Hの変化率の生成物から推定しました。従って、最高温度変化が0.00081°C秒-1でした。これは、5秒の時間間隔の間に1℃から1.5℃の間TPSセンサの温度上昇ΔTcよりはるかに低かったです。詳細評価と議論をSECに説明されています。参照9の4。
以下は、重要なプロトコルステップです。最初のステップは、サンプルの過冷却状態を維持しています。第二段階は、過冷却のΔTSUP度以下TPSセンサーの温度上昇ΔTcを維持することによって熱定数」の測定を行っています。
測定は、温度ドリフトの影響を受けないことを確実にするために、次のように確認されるべきです。まず、バルク温度変化がTPSセンサの温度上昇ΔTcよりはるかに低いことを保証します。第二に、MHの形成熱による温度変化がTPSセンサーの温度上昇ΔTcよりもはるかに低いことを確認してください。
試料が融解した場合、熱伝導率および比熱は、TPS技術によって無限大に発散します。このような場合には、センサからの出力を変更したり、測定時間を減少させます。
ガスハイドレートの特性低い生成速度はMHに固有ではないため、この測定方法は、水素、CO 2、及びオゾン水和物を含有するガス水和物水ゲストガスシステムの熱的特性に適用することができます。この方法で重要な点は、低ラットでありますターゲット材料の相転移のE。したがって、この方法は、低い相転移の速度で他の材料にも適用することができます。これらの水和物の生成速度は、過冷却条件下で十分に遅い場合は、この測定方法は、低濃度のTHF溶液とテトラブチルアンモニウムブロマイド（TBAB）水和物から形成されたテトラヒドロフラン（THF）水和物にも適用することができます。ここでの唯一の要件は、上記のように水和物の生成熱による温度変化は、センサの温度上昇よりもはるかに低いことを保証することです。これらのシステムの遷移速度は非常に速く、形成熱が大きく、測定に影響を与えるため、一方で、この技術は、水氷および化学量論THF溶液、水和物の相転移には適用できません。
ウェイトら 15は、砂、メタンガス、およびMHを含む試料の熱伝導率を測定しました。クマーら。16は同じ構成要素を有するサンプルを使用して熱拡散率を測定しました。これらは、加圧されたメタンガス雰囲気下で氷を使用して、砂の細孔内に直接MHを形成しました。氷の全ては、MHに変換しました。 MH形成が完全に停止するまで、このように、それらは、試料の熱伝導率を測定しました。この方法は、熱的特性の測定はMHの形成または解離熱により、試料の組成が一定であることが影響されないという利点を有します。しかし、この方法は、砂、水、メタン、およびMHを含む試料の熱特性を与えることはできません。黄及びファンハイドレートサンド試料17の熱伝導率を測定しました。これらは、MHの形成を促進し、ドデシル硫酸ナトリウム（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).
