Protocole de mesure des propriétés thermiques d'un synthétique Sable-eau-gaz-méthane Hydrate Sample surfondue

Introduction

Les hydrates de gaz sont des composés cristallins qui comprennent des structures en cage de molécules d'eau liées à l' hydrogène contenant des molécules hôtes dans la cage ^1. De grandes quantités d'hydrates de méthane (MHS) dans les régions de plancher de l' océan et le pergélisol sont des ressources intéressantes énergétiques futurs , mais peuvent affecter les conditions climatiques mondiales ^2.

En Mars 2013, la Japan Oil, Gas, and Metals National Corporation menées premier test de production offshore du monde pour extraire le gaz de sédiments naturels MH-portant dans l'est de Nankai Trough par la "méthode de dépressurisation" ^3,4.

Les hydrates de gaz peuvent stocker des gaz tels que le méthane , ^{1, 5} un atome d' hydrogène, CO ₂ de ^1,6 et de l' ozone ^7. Par conséquent, le méthane et d'hydrogène hydrates sont étudiés comme le stockage d'énergie potentielle et les médias de transport. Pour réduire les émissions de CO ₂ dans l'atmosphère, CO ₂ sequestration utilisant le CO ₂ hydrates dans les sédiments océaniques profonds ont été étudiés ^6. L'ozone est actuellement utilisé dans la purification de l'eau et de la stérilisation des aliments. Les études sur la technologie de conservation de l' ozone ont été menées parce qu'il est chimiquement instable ^7. La concentration en ozone dans l' hydrate est beaucoup plus élevée que celle de l'eau ozonisée ou de la glace ^7.

Pour développer la production de gaz à partir des sédiments MH-portant naturelles et des technologies à base d'hydrates, il est impératif de comprendre les propriétés thermiques des hydrates de gaz. Cependant, les données sur les propriétés thermiques et études de modèles de sédiments d'hydrate portant gaz sont rares ^8.

La «méthode de dépressurisation» peut être utilisé pour dissocier MH dans l'espace poreux des sédiments en diminuant la pression de pore inférieure à la stabilité de l'hydrate. Dans ce procédé, les composants spatiaux interstitielles des sédiments changent de l'eau et de l'eau MH, MH, et le gaz. La mesure de Les propriétés thermiquesde cette dernière condition est difficile parce que la chaleur de fusion de MH peut affecter les mesures. Pour résoudre ce problème, Muraoka et al. Réalisée la mesure des propriétés thermiques dans des conditions surfusion pendant la formation MH ^9.

Avec ce protocole vidéo, nous expliquons la méthode de mesure de surfusion synthétique échantillon de sable eau-gaz-MH.

La figure 1 illustre le dispositif expérimental permettant de mesurer les propriétés thermiques du méthane artificiel hydrates de sédiments. La configuration est la même que celle indiquée dans la référence ^9. Le système comprend essentiellement un récipient à haute pression, de pression et de régulation de la température et des propriétés thermiques du système de mesure. Le récipient à haute pression est en acier inoxydable cylindrique avec un diamètre interne de 140 mm et une hauteur de 140 mm; son volume intérieur avec le volume mort est éliminé 2110 cm ³ et sa limite de pression est de 15 MPa. Le transie source plane nt (TPS) technique est utilisée pour mesurer les propriétés thermiques ^10. Neuf sondes TPS avec des rayons individuels de 2.001 mm sont placés à l'intérieur du navire. La disposition des neuf sondes ⁹ est représenté sur la figure 2 en référence ^9. Les sondes TPS sont connectées à l'analyseur des propriétés thermiques d'un câble et commutées manuellement pendant l'expérience. Les détails du capteur TPS, schéma de connexion et la configuration dans le récipient sont représentés sur les figures S1, 2 et 3 de l'information à l' appui en référence ^9.

Figure 1:. Le dispositif expérimental pour mesurer les propriétés thermiques de l'artificiel sédiments hydrate de méthane porteur Le chiffre est modifié de la référence ^9.3956fig1large.jpg "target =" _ blank "> S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.

La méthode TPS a été utilisé pour mesurer les propriétés thermiques de chaque échantillon. Les principes de la méthode sont décrits dans la référence ^10. Dans cette méthode, l'augmentation de la température en fonction du temps, AT _ave, est

où

Dans l' équation 1, W ₀ est la puissance de sortie du capteur, r est le rayon de la sonde de détection, λ est la conductivité thermique de l'échantillon, α est la diffusivité thermique, et t est le temps depuis le début de l'alimentation à la sonde de capteur. D (τ) est une fonction dépendant du temps sans dimension. τ (aT / r) ^1/2. Dans l' équation 2, m est le nombre d'anneaux concentriques de la sonde TPS et I ₀ est une fonction de Bessel modifiée. La conductivité thermique, la diffusivité thermique et la chaleur spécifique de l'échantillon sont simultanément déterminées par une analyse d'inversion appliquée à l'augmentation de la température que l'énergie est fournie à la sonde de détection.

Protocol

Remarque: S'il vous plaît consulter toutes les fiches de données de sécurité des matériaux pertinents comme cette étude utilise haute pression inflammable gaz de méthane et un grand navire à haute pression. Portez un casque, des lunettes de sécurité, et des bottes de sécurité. Si le système de régulation de température cesse, la pression dans la cuve augmente avec MH dissociation. Pour éviter tout accident, l'utilisation d'un système de soupape de sécurité est fortement recommandée pour libérer automatiquement le gaz de méthane dans l'atmosphère. Le système de soupape de sécurité peut fonctionner sans alimentation électrique.

1. Préparation des échantillons de gaz de sable à l' eau-méthane ⁹

1. Placer le récipient à haute pression sur la table vibrante.
2. Verser 1,5 L d'eau pure dans une bouteille d'eau et 4000 g de sable de silice dans une bouteille de sable. Peser les masses de sable et de l'eau dans les bouteilles de sable et d'eau, respectivement.
3. Versez 1 litre d'eau pure dans le récipient à haute pression avec un volume intérieur de 2110 cm ³ d'une bouteille d'eaujusqu'à ce que l'eau remplit la moitié de la cuve interne.
4. Tourner sur la table vibrante pour faire vibrer l'ensemble du récipient. Réglez le taux de vibration et de l'alimentation à 50 Hz et 220 W, respectivement. Appliquer la vibration jusqu'à la fin de l'étape 1.5. Retirer l'air résiduel dans la conduite d'évacuation et un filtre métallique fritté, au fond du récipient en faisant vibrer la cuve.
5. Verser 3,300 g de sable de silice d'une bouteille dans le récipient de sable à une vitesse constante d'environ 1 g s - ¹ en utilisant un entonnoir tenu près de la surface de l' eau alors que l'ensemble du récipient est vibré pour assurer l' emballage uniforme.
6. Arrêter la vibration lorsque l'eau atteint le bord du récipient.
7. Placer un anneau comme un mur temporaire sur le bord du récipient pour empêcher l'eau de se répandre.
8. Vibrer à nouveau le récipient à 50 Hz et 220 W.
9. Lorsque le sable atteint le bord de la cuve (hauteur 140 mm), éteindre la vibration.
10. Retirez le mur temporaire et de l'eau des pores en excès par eligne de fuite des e. Versez le dos excès d'eau des pores dans la bouteille d'eau.
11. Emballez le sable en faisant vibrer le récipient une ou deux fois à 50 Hz et 300 W pendant 1 sec et ajouter plus de sable si nécessaire.
12. Peser les masses de sable et de l'eau dans les bouteilles de sable et d'eau. Calculer les masses de sable et d'eau dans le récipient à partir des différences de masse dans les bouteilles de sable et d'eau. Dans cette expérience, les masses de sable et de l'eau dans la cuve étaient 3,385 g et 823,6 g, respectivement. La masse d'eau dans le récipient est notée en poids au _total.
13. Couvrir le récipient à haute pression avec un couvercle en acier inoxydable et serrer les boulons de paires diagonalement opposées dans l'ordre.
14. Déplacer le récipient à haute pression de la table vibrante à la table destinée à l'expérience.
15. Recouvrir la cuve haute pression avec l'isolant thermique pour contrôler la température.
16. Connecter les gazoducs à haute pression et le refroidissement des conduites d'écoulement de l'eau dans la cuve à haute pression. Ouvrir les vannes des conduites d'entrée et du gaz de sortie. Ventiler 10 L de méthane à un débit de 800 ml min ^-1 jusqu'à aucun rejet d'eau en excès dans le piège sous pression atmosphérique. La décharge de sable est empêchée par un filtre métallique fritté fixée sur le fond du récipient. L'eau résiduelle reste sur la surface du sable, car le sable de silice hydrophile adsorbe les molécules d'eau.
17. Peser la masse d'eau dans le piège, w _piège, pour déterminer le volume de gaz dans la cuve. Déterminer la masse d'eau résiduelle, _res w, dans le récipient en utilisant l'équation w _res = w _totale - w _piège. Dans ce cas, _res w et _piège w était 360,6 g et 463,0 g, respectivement.
18. Déterminer la porosité de l'échantillon en utilisant la formule Ѱ = 1 - V _{cellulaire sable} / V, où le _sable V est le volume de til sable déterminé par le rapport de la masse de sable à la densité du sable (c. -à ρ _s = 2.630 kg m ^-3) et les _cellules V est le volume interne du récipient. Ѱ la porosité de l'échantillon était de 0,39.
19. Fermer la vanne de la conduite de gaz de sortie. Injecter du méthane pour augmenter la pression des pores de méthane dans le récipient à environ 12,1 MPa à la température ambiante (c. -à- 31,6 ° C).
20. Fermer la vanne de la conduite de gaz d'entrée.
21. Commencer à enregistrer la pression et la température dans le récipient au cours de l'expérience en utilisant l'enregistreur de données. L'intervalle d'échantillonnage de données est 5 s. Le temps expérimentale totale est d'environ 3000 min.

2. MH Synthèse et la mesure de l'échantillon surfondue ⁹ de Propriétés thermiques

1. Mettre en marche le refroidisseur pour refroidir le récipient de la température ambiante à 2.0 ° C en faisant circuler le liquide de refroidissement. Laissez circuler le liquide de refroidissement du refroidisseur to le fond du récipient, à partir de là vers le couvercle de la cuve, et enfin vers le refroidisseur. Le taux de variation de la température dans la cuve était d' environ 0.001 ° C sec ^-1.
2. Définissez les paramètres de mesure à l'aide du logiciel de l'analyseur TPS. Réglez le type de capteur à la conception # 7577. Régler la puissance de sortie W ₀ à 30 mW et le temps de mesure de 5 secondes. Notez que les paramètres appropriés doivent être modifiés si le type de capteur ou conditions échantillons changent. Réglez les paramètres pour augmenter la température de 1 ° C à 1,5 ° C.
3. Calculer le degré de surfusion, AT _sup, avec l'équation suivante:
   AT _sup = T _eq (P) - T. (3)
   _Teq (P) est la température d'équilibre de MH en fonction de la pression P. _Teq (P) est calculée à l' aide du logiciel CSMGem ^1.0, P et T sont la pression et la température dans le récipient de mesure à l'aide de jauges de pression et de température, respectivement.
4. Mesurer simultanément la conductivité thermique, diffusivité thermique, et la chaleur spécifique volumétrique en utilisant l'analyseur TPS après AT _sup est supérieur à 2 ° C.
5. Mettez la sonde TPS reliée aux propriétés analyseur thermique après chaque mesure. Mettez les câbles entre les sondes TPS et l'analyseur manuellement pendant l'expérience ^9. Le schéma de raccordement est illustré à la figure S2 en référence ^9. La séquence de commutation de chaque capteur est négative. 6 → 2 → 7 → 5 → 1 → 9 → 4 → 3 → 8 → 6 .... La séquence est basée sur la distance entre les capteurs, qui est fixé dans la mesure du possible d'éviter que la chaleur résiduelle d'affecter les mesures. Recueillir des données toutes les 3-5 min.
6. Répéter les mesures jusqu'à ce que & #916; T _sup atteint 2 ° C à nouveau. Dans cette expérience, Δ T _sup augmente initialement avec le temps. Après _sup AT atteint la valeur maximale, diminue progressivement _sup AT à 0 ° C , car la pression diminue avec la formation de MH. Vérifiez si AT _sup est supérieur à 2 ° C avant les mesures TPS utilisant l' équation 3.
7. Assurez-vous que le profil de température est pas affectée par MH fusion. Si MH fond lors de la mesure, la température ne sera pas augmenter à cause de la fusion MH est une réaction endothermique. Vérifiez le profil de température pendant les mesures, et est discuté dans la section des résultats.
8. Effectuer l'analyse des propriétés thermiques pour toutes les données de profil de température en utilisant la technique TPS.

3. Calcul de la Saturation Changement de l'échantillon ^9,11

Remarque:Le degré de saturation pour MH, d' eau et de gaz dans l'échantillon en fonction du temps t est calculée en utilisant l'équation d'état du gaz. Les détails de calcul et les équations utilisées sont décrites précédemment ^11.

1. Calculer le V _gaz de volume de gaz méthane, _t au temps t
   
   où Q est le volume initial de gaz dans la cuve, V _{MH, t - 1} est le volume de MH à l' instant t - 1 et R _VHW est le rapport volumique de l' eau et MH.
   
   Dans l' équation 5, n est le nombre de MH d'hydratation (~ 6), et ρ ρ _{MH eau} correspond à la densité de MH et de l' eau, respectivement, et w _MH et de l' _eau w représentent la masse moléculaire MH et de l' eau, respectively.
2. Calculer la quantité ΔM _t (mole) de MH formé à partir de t - 1 à t
   
   où R est la constante des gaz, P est la pression du gaz de méthane, et _{Zt (gaz} T, _T, P _{gaz, t)} est le coefficient de méthane de compression à l' instant t. Nous ⁹ et Sakamoto et al. ¹¹ avons utilisé l'équation Benedict-Webb-Rubin (BWR), tel que modifié par Lee et Kesler, pour calculer Z _t ^{12, 13.} Pour ce calcul, les formules (3-7.1) - (3-7.4) de l'équation BWR ¹³ et les constantes de Lee-Kesler sont utilisées dans les tableaux 3 - 7 de référence ^13.
3. Calculer le changement de volume Δ V _MH, t - 1 à t
   
   où P _s est la pression de 101325 Pa de référence, T _s est la température de référence de 273,15 K, Z _s est le coefficient de compression au niveau de P _s et _de T (Z _s ~ 1), et V _CH4 est le rapport du volume de gaz méthane dans le unité de volume MH [Nm ³ m ^-3]. Utilisez une valeur V _CH4 de 165,99 [Nm ³ m ^-3].
4. Calculer le volume V _{MH, t} de MH à l' instant t
   
5. Calculer le volume _d'eau de l' eau de V, _T dans le récipient sous pression à l' instant t
   où l' _eau de V ₁ est le volume initial d'eau.
6. Répéter les calculs en utilisant les équations. 9.4 à l' instant t = 2, 3, ... afin de déterminer la variation de la saturation de l' eau, de méthane et ¹¹ MH. La condition initiale est t = 1, c. -à- _gaz V, ₁ = Q. Le P et T à l' instant t sont prises à partir des journaux de données ^9. Les résultats du calcul sont présentés dans la section suivante.

Representative Results

La figure 2a représente le profil de température qui ne soit pas affectée par MH fusion. _C AT est la variation de température due à la mesure des constantes thermiques. La figure 2b montre le profil de température qui est affectée par MH fusion. Le profil de la figure 2b ne peut être analysé par les équations 1 et 2 parce que ces équations sont dérivées en supposant que les conditions d'échantillonnage stables.

La figure 3a montre la pression, la température et le degré de super - refroidissement dans le récipient en fonction du temps. nucléation MH après que le système a atteint la pression et de l'équilibre de température. La formation de MH est marquée par un changement de pression drastique au temps t = 170 min. Les doubles flèches montrent que le degré de surfusion est supérieur à 2 ° C. Les constantes thermiques ont été mesurées dans cette plage. t. Le calcul de la saturation est décrite dans la section 3. La saturation est définie comme S _{i, t} = V _{i, t} / (V _cellule - V _sable), où i désigne le MH, l' eau et les composants de gaz de méthane à l' instant t. A t = 170 min, MH a commencé à former et S _MH a augmenté de manière significative. Entre 170 et 2500 min, _MH S a augmenté de 0 à 0,32, alors que S S _eau et de _gaz a diminué de 0,43 à 0,18 et de 0,56 à 0,50, respectivement. Après 2500 minutes, le MH, l'eau et la saturation en gaz étaient à peu près constante.

La figure 4 montre un exemple d'les mesures des constantes thermiques. Les conditions expérimentales étaient t = 825 min, P= 7,1 MPa T = 2,4 ° C, S _h = 0,16, S _g = 0,53 et S _w = 0,31. La figure 4a représente le profil de température. Le logiciel d'analyse 200 TPS enregistre les points de données régulièrement espacés dans le temps sur un intervalle de temps prédéfini; Ainsi, les données sont sélectionnées pour une analyse à partir des points de données 200. Les doubles flèches indiquent la plage de données utilisées dans l'analyse. Le temps varie d'analyses 1 et 2 sont 0-5 sec et 0,65 à 4,88 sec, respectivement. Analyses 1 et 2 sont des exemples de plages inappropriées et appropriées, respectivement. Les figures 4b et 4c ont été obtenus en utilisant la technique TPS dans chaque plage d'analyse. La figure 4b montre _l'ave changement de température de AT (de τ) et D (τ) avec AT _ave (τ ) = AT _c (t). La relation entre AT _ave (&# 964;.) Et D (τ) varie en fonction de l'analyse plage Figure 4c montre la température T _d vs la racine carrée du temps t. L'écart entre les données de température de l'ajustement linéaire obtenu par l'analyse d'inversion de TPS est _Td. La déviation de l' analyse 1 au début de la mesure est assez grande, comme le montre la figure 4c, ce qui suggère que la couche isolante de la sonde de détection de TPS affecte les mesures.

Le tableau 1 énumère les constantes thermiques dans chaque plage d'analyse tel que mentionné ci - dessus. Le rapport total de temps caractéristique est définie par le temps total d'analyse (pour t = 2 à 4 secondes, le temps total est de 4 s) , divisé par caractéristique τ du temps. A noter que le rapport total de temps caractéristique doit être inférieure à 1 lors de l'utilisation de la technique de TPS. Ceci est décrit en référence T _d.

Afin d'éviter la sonde de détection affectant les mesures, données au début de chaque mesure ne doivent pas être utilisés. L'écart moyen de _Td est minimisé, comme le montre la figure 4c, en ajustant l'intervalle de temps d'analyse. Le rapport total de temps caractéristique est ajustée à l'unité en ajustant l'intervalle de temps d'analyse. Par conséquent, nous avons adopté les valeurs des constantes thermiques de l'analyse 2 pas 1.

Conductivité thermique, la chaleur spécifique et de la diffusivité thermique sont présentées en fonction du temps dans les figures 5a, b et c, respectivement. Enfin, nous résumons les résultats du tpropriétés Hermal et la saturation de l'hydrate. Les détails concernant les résultats sont donnés à la Sec. 4 refference ^9.

Figure 2:. Les profils de température en fonction du temps (a) pas affectée par MH de fusion (conditions de surfusion) et (b) affectée par MH fusion S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.

A noter que les deux profils de température sont des expériences préliminaires. Le temps de mesure est plus longue que dans l'expérience afin de clarifier l'effet de la chaleur de fusion. Dans les expériences préliminaires, le temps de mesure t était de 40 s et la puissance de sortie W ₀ a 20 mW (a) et 50 mW (b).

<p class = "jove_content" fo: keep-together.within-page = "1">
Figure 3: (a) la pression, la température et le degré de super - refroidissement dans le récipient en fonction du temps. Les doubles flèches montrent que le degré de super refroidissement est supérieur à 2 ° C. Les constantes thermiques ont été mesurées dans cette plage. (B) Le MH, l' eau et le méthane saturation en gaz de l'échantillon sont présentées en fonction du temps (réimprimé de la référence ^9). S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 4: Exemple d'analyse des mesures des constantes thermiques (a) le profil de température en utilisant le.méthode de mesure TPS. Le temps varie d'analyses 1 et 2 sont 0-5 sec et 0,65 à 4,88 sec, respectivement. (B) Relation entre _l'ave changement de température de AT (de τ) et D (τ) avec AT _ave (τ) = AT _c (t). (C) Température T _d vs le carré du temps de racine t. S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.

Figure 5: (a) la conductivité thermique λ en fonction du temps, (b) ρC chaleur spécifique _p en fonction du temps, et (c) diffusivité thermique α en fonction du temps.Les résultats ont été convertis en des propriétés thermiques en fonction de la saturation MH. Les résultats convertis et discussion pertinents sont présentés dans Réf. 9. Les données montrent un chevauchement dans l'intervalle t = 210-980 min. Pour plus de clarté, les données tracées représentent la moyenne de trois mesures à partir du même capteur dans cette plage. Ces chiffres ont été modifiés de la référence ^9. S'il vous plaît cliquer ici pour voir une version plus grande de cette figure.

	gamme d'analyse, s	λ, W m - 1 K - 1	ρC p, MJ m - 3 K - 1	α, mm 2 s - 1	Total à Char. Temps	Moyenne Dev., ° C
Analyse 1	0,00 à 5,00	2.12	0,938	2.26	2.11	0,01018
Analyse 2	0,65 à 4,88	2.31	2.11	1.10	1.00	0.00061

Tableau 1:. Constantes thermiques pour chaque gamme d'analyse Analyses 1 et 2 sont des exemples de plages inappropriées et appropriées, respectivement.

Discussion

L'effet de la chaleur de formation de MH sur mesure a été estimé. La chaleur de formation de MH a été estimé à partir des produits de la vitesse de changement S _h comme indiqué sur la figure 3b et l'enthalpie de formation H = 52,9 kJ mol ^-1 MH ^14. Par conséquent, la variation de la température maximale était de 0,00081 ° C sec ^-1. Cela a été beaucoup plus faible que l'augmentation de température AT _c du capteur TPS comprise entre 1 ° C et 1,5 ° C pendant l'intervalle de temps de 5 sec. estimation détaillée et la discussion sont décrits dans Sec. 4 de référence ^9.

Voici les étapes critiques de protocole. Première étape est de maintenir les conditions échantillon de surfusion. La seconde étape effectue les mesures de constantes thermiques en maintenant l'augmentation de température AT _c du capteur TPS au- dessous du degré de surfusion AT _sup.

Faire en sorte que la mesure ne soit pas influencée par la variation de température, les conditions suivantes doivent être confirmées. Tout d' abord, assurez -vous que le changement de température de la masse est beaucoup plus faible que l'augmentation de température AT _c du capteur TPS. D' autre part, veiller à ce que la variation de température due à la chaleur de formation de MH est beaucoup plus faible que l'augmentation de température AT _c du capteur de TPS.

Si un échantillon fond, la conductivité thermique et la chaleur spécifique divergent à l'infini par la technique de TPS. Dans de tels cas, modifier la puissance de sortie du capteur ou de diminuer le temps de mesure.

Cette méthode de mesure peut être appliquée aux propriétés thermiques du système de gaz d' hydrates de gaz-eau-hôtes, qui contient de l' hydrogène, CO _2, et hydrates d'ozone, car le taux de faible formation caractéristique de l' hydrate de gaz ne sont pas uniques à MH. Le point clé de cette méthode est le faible rate de transition de phase du matériau cible. Par conséquent, cette méthode peut être appliquée à d'autres matériaux ayant un taux de transition de phase faible. Cette méthode de mesure peut également être appliquée à du tétrahydrofuranne (THF), l'hydrate formé à partir d'une solution de THF à faible concentration et le bromure de tétra butyl (TBAB) d'hydrate si le taux de ces hydrates de formation est suffisamment lente dans des conditions de surfusion. La seule condition est d'assurer que la variation de température due à la chaleur de formation de l'hydrate est nettement inférieure à l'augmentation de la température dans le capteur, comme mentionné ci-dessus. D'autre part, cette technique ne peut pas être appliquée à l'eau glacée et stoechiométrique THF solution-hydrate de transition de phase, car le taux de transition dans ces systèmes est très rapide et la chaleur de formation affecte de manière significative les mesures.

Waite et al. ¹⁵ ont mesuré la conductivité thermique des échantillons comprenant du sable, du gaz méthane, et MH. Kumar et al.¹⁶ mesurer la diffusivité thermique en utilisant des échantillons avec les mêmes composants. Ils ont formé MH directement dans les pores de sable en utilisant la glace d'eau dans une atmosphère de gaz méthane sous pression. Tout de la glace d'eau a été converti en MH. Ainsi, ils ont mesuré la conductivité thermique de l'échantillon jusqu'à la formation MH complètement arrêtée. Cette méthode présente l'avantage que les mesures des propriétés thermiques ne sont pas affectées par la formation ou la dissociation thermique des MH et que la composition de l'échantillon est constante. Cependant, cette méthode ne peut pas donner les propriétés thermiques des échantillons comprenant le sable, l'eau, le méthane, et MH. Huang et Fan ont mesuré la conductivité thermique d'un hydrate de support de l' échantillon de sable ^17. Ils ont formé MH dans les pores de sable à l'aide de dodécylsulfate de sodium (SDS) solution, ce qui a facilité la formation MH. Ils ont noté que le gaz et l'eau est probablement resté dans les pores de sable et le gaz affectés de manière significative les mesures. Cependant, ils ne signalent pas la composition du water et du gaz. Notre protocole de mesure a l'avantage de donner la relation entre les propriétés thermiques (conductivité thermique, diffusivité thermique et la chaleur spécifique volumétrique) et la composition du sédiment MH-roulement comprenant le sable, l'eau, le méthane, et MH.

Pour développer des technologies de l'hydrate de gaz de production de masse, les constantes thermiques de la formation d'hydrates sont nécessaires, et la méthode de mesure proposée fait exactement cela.

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