The unscheduled formation of gas hydrate plugs in oil and gas pipelines, which can lead to serious mechanical and personnel damage, is a problematic issue in the petroleum industry. Traditionally, thermodynamic inhibitors such as methanol have been used to control the formation of gas hydrates, but due to the large expenses and ecological risks associated with its use there is increased interest in the use of alternative hydrate inhibitors. They include kinetic inhibitors (KIs) and antiagglomerants (AAs) and as their names imply, function by interfering with the kinetics of hydrate formation and hydrate agglomeration. Recently, antifreeze proteins (AFPs) have shown to inhibit hydrates and have been proposed as hydrate inhibitors. Normally, AFPs function to protect the tissues of various organisms during freezing conditions. Initially they were found in polar fish, and were later recognized in insects, plants and microorganisms. AFPs are thought to function by lowering the freezing point of water through an adsorption-inhibition mechanism. This thesis has shown that antifreeze proteins (AFPs) are able to modify the crystal morphologies of structure II (sII) tetrahydrofuran (THF) similarly to the KI poly-N-vinylpyrrolidone (PVP) by adhering to the hydrate surface and inhibiting crystal growth. The AFPs were also tested on a high-pressure sII methane/ethane/propane hydrate and proved to have superior hydrate inhibition to PVP. Yet, the expense of purifying AFPs makes them impractical for industrial purposes, thus investigations into the use of cold-adapted bacteria as hydrate inhibitors proved that isolates capable of adsorbing to THF hydrate showed the most effective THF hydrate inhibition. These findings suggest a potential for the future development of biologically-based hydrate inhibitors. / Thesis (Master, Biology) -- Queen's University, 2009-09-01 10:04:00.72

The oil and gas industry is often plagued by the formation of clathrate hydrates in oil pipelines. While the industry originally had a heuristic of avoidance of clathrate hydrates they are moving to a heuristic of risk management. To successfully implement a risk management heuristic, time dependent phenomena of clathrate hydrate formation and flowline plugging must be known. The study of time dependent phenomena of formation and agglomeration are investigated using a TA Instruments AR-G2 rheometer with a pressure cell capable of operating at up to 13.8 MPa. Pressurized rheological experiments examine clathrate hydrates formed in situ. Both shear and oscillatory experiments have been conducted on the samples, giving flow and viscoelastic parameters. Shear experiments show sharp increases in viscosity upon clathrate hydrate formation indicating rapid aggregation. Transient oscillation experiments show a sharp increase in the elastic and loss moduli followed by a decrease in the loss moduli. Thus, both in situ clathrate hydrate formation and annealing are quantified. In addition these oscillatory measurements provided a novel technique for non-destructive investigation of clathrate hydrate aggregation over time.

As the oil industry moves from a heuristic of avoidance of hydrates to a heuristic of risk management time dependent phenomena of hydrate formation and plugging must be known. One of the key parameters to this process is the aggregation of hydrate particles, the fractal networks they form, and the effect these two parameters have on flow. Unfortunately the aggregation and fractal structure information is extremely difficult to acquire experimentally, for this reason a three-dimension discrete element method (3D-DEM) model has been implemented. The 3D-DEM model calculates detailed solutions to Newton's equations of motion for individual particles. In addition these particles are coupled with the surrounding fluid through computational fluid dynamics (CFD). This coupled 3D-DEM can be used to investigate what the effects of shear, suspending viscosity, attractive forces, and other relevant variables have on the structure, stresses, and positions of the hydrate particles over time. In addition, the effect on viscosity has been calculated using CFD and compared back to basic hard sphere theory.

The first project in this work used classical molecular dynamics to study the ice recrystallization inhibition potential of a series of carbohydrates and alcochols, using the hydration index, partial molar volumes and isothermal compressibilities as parameters for measuring their cryogenic efficacy. Unfortunately, after 8 months of testing, this work demonstrates that the accuracy and precision of the density extracted from simulations is not sufficient in providing accurate partial molar volumes. As a result, this work clearly demonstrates that current classical molecular dynamics technology cannot probe the volumetric properties of interest with sufficient accuracy to aid in the research and development of novel cryoprotectants.The second project in this work used molecular dynamics simulations to evaluate the Gibbs free energy change of substituting CO2 in sI clathrate hydrates by N2,CH4, SO2 and H2S flue gas impurities under conditions proposed for CO2 capture (273 K, 10 bar). Our results demonstrate that CO2 substitutions by N2 in the small sI cages were thermodynamically favored. This substitution is problematic in terms of efficient CO2 capture, since the small cages make up 25% of the sI clathrate cages, therefore a significant amount of energy could be spent on removing N2 from the flue gas rather than CO2. The thermodynamics of CO2 substitution by CH4, SO2 and H2S in sI clathrate hydrates was also examined. The substitution of CO2 by these gases in both the small and large cages were determined to be favorable. This suggests that these gases may also disrupt the CO2 capture by sI clathrate hydrates if they are present in large concentrations in the combustion flue stream. Similar substitution thermodynamics at 200 K and 10 bar were also studied. With one exception, we found that the substitution free energies do not significantly change and do not alter the sign of thermodynamics. Thus, using a lower capture temperature does not significantly change the substitution free energies and their implications for CO2 capture by sI clathrate hydrates.

When ice is compressed alongside a gas, crystalline 'host - guest' inclusion compounds known as gas clathrate hydrates form. These compounds are of interest not only for their environmental and possible technological impact as gas storage and separation materials, but also for their ability to probe networks not readily adopted by the pure `host' water molecules, and to study the interactions between water and gas molecules. Despite the pressure dependent crystal structures being fully determined for a large variety of `guest' gas species there is still relatively little known about the crystal structures in small guest gas systems such as H2 hydrate. The majority of structural studies have been done with x-ray diffraction and report a number of conflicting structures or hydrogen contents for the four known stable phases (sII, C0, C1 and C2). As this is a very hydrogen rich system the most ideal method to study the structure is neutron diffraction, which is able to fully determine the location of the hydrogen atoms within the structure and would allow a direct measurement of any hydrogen ordering within the host structure and the H2 content. In this work the phase diagram of the deuterated analogue of the H2-H2O system is explored at low pressures (below 0.3 GPa) with neutron diffraction. In the pressure/temperature region where the sII phase is known to be stable, two metastable phases were observed between the formation of sII from ice Ih and that this transition sequence occurred in line with Ostwald's Rule of Stages. One of these metastable phases was the C0 phase known to be stable in the H2-H2O system above 0.5 GPa, and the other is a new structure not previously observed in this system and is dubbed in this work as C-1 . Prior to this work the C0 phase has been reported with various structures that were determined with x-ray diffraction, and here the crystal structure and H2 content at low pressure are determined with neutron diffraction. The C0 phase was found to form a similar host structure to those of the previous studies with spiral guest sites but is best described with highly mobile H2 guests and a higher symmetry space group which make it the same structure as the spiral hydrate structure (s-Sp) recently observed in the CO2 hydrate system. In addition to this structure being determined at pressure a sample of C0 was also recovered to ambient pressure at low temperature and its structure/H2 content is presented as it was warmed to decomposition. The crystal structure of the C-1 phase was determined to be similar to ice Ih and a sample was recovered to ambient pressure to study its decomposition behaviour. Evidence for a similar structure in the helium hydrate system at low pressure is also reported here. This work was then extended to higher pressures with the recent developments of a hydrogen-compatible gas loader and large-volume diamond anvil cells. Several test experiments on gas-loaded Paris-Edinburgh presses are described on systems that are similar to hydrogen-water like urea-hydrogen and neon-water. And a further preliminary high pressure study on the deuterated analogue of the H2- H2O system in a diamond anvil cell between 3.6 and 28 GPa shows decomposition behaviour as pressure was increased.

The first project in this work used classical molecular dynamics to study the ice recrystallization inhibition potential of a series of carbohydrates and alcochols, using the hydration index, partial molar volumes and isothermal compressibilities as parameters for measuring their cryogenic efficacy. Unfortunately, after 8 months of testing, this work demonstrates that the accuracy and precision of the density extracted from simulations is not sufficient in providing accurate partial molar volumes. As a result, this work clearly demonstrates that current classical molecular dynamics technology cannot probe the volumetric properties of interest with sufficient accuracy to aid in the research and development of novel cryoprotectants.The second project in this work used molecular dynamics simulations to evaluate the Gibbs free energy change of substituting CO2 in sI clathrate hydrates by N2,CH4, SO2 and H2S flue gas impurities under conditions proposed for CO2 capture (273 K, 10 bar). Our results demonstrate that CO2 substitutions by N2 in the small sI cages were thermodynamically favored. This substitution is problematic in terms of efficient CO2 capture, since the small cages make up 25% of the sI clathrate cages, therefore a significant amount of energy could be spent on removing N2 from the flue gas rather than CO2. The thermodynamics of CO2 substitution by CH4, SO2 and H2S in sI clathrate hydrates was also examined. The substitution of CO2 by these gases in both the small and large cages were determined to be favorable. This suggests that these gases may also disrupt the CO2 capture by sI clathrate hydrates if they are present in large concentrations in the combustion flue stream. Similar substitution thermodynamics at 200 K and 10 bar were also studied. With one exception, we found that the substitution free energies do not significantly change and do not alter the sign of thermodynamics. Thus, using a lower capture temperature does not significantly change the substitution free energies and their implications for CO2 capture by sI clathrate hydrates.

Clathrate hydrates can be reasonable choices for high-density hydrogen storage into compact host media, which is an essential task for hydrogen-based future society. However, conventional storage scheme where aqueous solution is frozen with hydrogen gas was impractically slow for practical use. Here we propose a much faster scheme where hydrogen gas was directly charged into hydrogen-free, crystalline hydrate powders. The storage kinetics was observed in situ by nuclear magnetic resonance (NMR) spectroscopy in a pressurized tube cell. At pressures up to 20 MPa the storage was complete within 80 minutes, as observed by growth of stored-hydrogen peak into the hydrate. Since the rate-determining step of current storage scheme is body diffusion of hydrogen within the crystalline hydrate media, we have measured the diffusion coefficient of hydrogen molecules using the pulsed field gradient NMR method. The results show that at temperatures down to 250 K the stored hydrogen is highly mobile, so that the powdered hydrate media should work well even in cold environments. Compared with more prevailing hydrogen storage media such as metal hydrides, the clathrate hydrate could offer even more advantages: It is free from hydrogen embrittlement, more chemically durable, more environmentally benign, as well as economically quite affordable.

The first project in this work used classical molecular dynamics to study the ice recrystallization inhibition potential of a series of carbohydrates and alcochols, using the hydration index, partial molar volumes and isothermal compressibilities as parameters for measuring their cryogenic efficacy. Unfortunately, after 8 months of testing, this work demonstrates that the accuracy and precision of the density extracted from simulations is not sufficient in providing accurate partial molar volumes. As a result, this work clearly demonstrates that current classical molecular dynamics technology cannot probe the volumetric properties of interest with sufficient accuracy to aid in the research and development of novel cryoprotectants.The second project in this work used molecular dynamics simulations to evaluate the Gibbs free energy change of substituting CO2 in sI clathrate hydrates by N2,CH4, SO2 and H2S flue gas impurities under conditions proposed for CO2 capture (273 K, 10 bar). Our results demonstrate that CO2 substitutions by N2 in the small sI cages were thermodynamically favored. This substitution is problematic in terms of efficient CO2 capture, since the small cages make up 25% of the sI clathrate cages, therefore a significant amount of energy could be spent on removing N2 from the flue gas rather than CO2. The thermodynamics of CO2 substitution by CH4, SO2 and H2S in sI clathrate hydrates was also examined. The substitution of CO2 by these gases in both the small and large cages were determined to be favorable. This suggests that these gases may also disrupt the CO2 capture by sI clathrate hydrates if they are present in large concentrations in the combustion flue stream. Similar substitution thermodynamics at 200 K and 10 bar were also studied. With one exception, we found that the substitution free energies do not significantly change and do not alter the sign of thermodynamics. Thus, using a lower capture temperature does not significantly change the substitution free energies and their implications for CO2 capture by sI clathrate hydrates.

No description available.

L'hydrate de méthane est un cristal non-stœchiométrique dans lequel les molécules d'eau forment des cages liées par liaison hydrogène qui piégent des molécules de méthane. Des ressources abondantes en hydrate de méthane peuvent être trouvées sur Terre, en particulier dans les roches poreuses minérales (par exemple, l'argile, le permafrost, les fonds marins, etc.). Pour cette raison, la compréhension de la thermodynamique et de la cinétique de formation de l'hydrate de méthane confiné dans des milieux poreux suscite beaucoup d'attention. Dans cette thèse, nous combinons la modélisation moléculaire et des approches théoriques pour déterminer la thermodynamique et la cinétique de formation de l'hydrate de méthane confiné dans des milieux poreux. Tout d'abord, l'état de l'art en matière de thermodynamique et de cinétique de formation de l'hydrate de méthane est présenté. Deuxièmement, différentes stratégies de simulation moléculaire, y compris des calculs d'énergie libre utilisant l'approche de la molécule d'Einstein, la méthode de coexistence directe et la technique textit{hyperparallel tempering}, sont utilisées pour évaluer la stabilité de l'hydrate de méthane à différentes températures et pressions. Troisièmement, parmi ces stratégies, la méthode de coexistence directe est choisie pour déterminer le déplacement du point de fusion lors du confinement dans des pores, $Delta T_m = T_m^{pore} - T_m^{bulk} $ où $ T_m^{pore}$ et $T_m^{bulk}$ sont les températures de fusion d'hydrate de méthane non confiné et confiné. Nous avons constaté que le confinement diminue la température de fusion, $T_m^{pore} < T_m^{bulk} $. Le changement de température de fusion en utilisant la méthode de la coexistence directe est cohérent avec l'équation de Gibbs-Thompson qui prédit que le décalage de la température de fusion dépend linéairement de l'inverse de la taille des pores, $Delta T_m/T_m^{bulk} sim k_{GT}/ D_p$. La validité quantitative de cette équation thermodynamique classique pour décrire de tels effets de confinement et de surface est également abordée. Les tensions de surface des interfaces hydrate-substrat et eau-substrat sont déterminées à l'aide de la dynamique moléculaire pour valider quantitativement l'équation de Gibbs-Thompson. Des simulations de dynamique moléculaire sont également effectuées pour déterminer les propriétés thermodynamiques importantes de l'hydrate de méthane non confiné et confiné: (a) conductivité thermique $lambda$ en utilisant le formalisme de Green-Kubo et la fonction d'autocorrélation du flux thermique; (b) expansion thermique $alpha_P$ et compressibilité isotherme $kappa_T$. Enfin, des conclusions et perspectives pour des travaux futurs sont présentées. / Methane hydrate is a non-stoichiometric crystal in which water molecules form hydrogen-bonded cages that entrap methane molecules. Abundant methane hydrate resources can be found on Earth, especially trapped in mineral porous rocks (e.g., clay, permafrost, seafloor, etc.). For this reason, understanding the thermodynamics and formation kinetics of methane hydrate confined in porous media is receiving a great deal of attention. In this thesis, we combine computer modeling and theoretical approaches to determine the thermodynamics and formation kinetics of methane hydrate confined in porous media. First, the state-of-the-art on the thermodynamics and formation kinetics of methane hydrate is presented. Second, different molecular simulation strategies, including free energy calculations using the Einstein molecule approach, the direct coexistence method, and the hyperparallel tempering technique, are used to assess the phase stability of bulk methane hydrate at various temperatures and pressures. Third, among these strategies, the direct coexistence method is chosen to determine the shift in melting point upon confinement in pores, $Delta T_m = T_{m}^{pore} - T_{m}^{bulk}$ where $T_m^{pore}$ and $T_m^{bulk}$ are the melting temperatures of bulk and confined methane hydrate. We found that confinement decreases the melting temperature, $T_m^{pore}<T_m^{bulk}$. The shift in melting temperature using the direct coexistence method is consistent with the Gibbs-Thompson equation which predicts that the shift in melting temperature linearly depends on the reciprocal of pore width, i.e., $Delta T_m/T_m^{bulk} sim k_{GB}/D_p$. The quantitative validity of this classical thermodynamic equation to describe such confinement and surface effects is also addressed. The surface tensions of methane hydrate-substrate and liquid water-substrate interfaces are determined using molecular dynamics to quantitatively validate the Gibbs-Thompson equation. Molecular dynamics simulations are also performed to determine important thermodynamic properties of bulk and confined methane hydrate: (a) thermal conductivity $lambda$ using the Green-Kubo formalism and the autocorrelation function of the heat-flux and (b) the thermal expansion $alpha_P$ and isothermal compressibility $kappa_T$. Finally, some conclusions and perspectives for future work are given.