Pinnelli S. R. Prasad*a,
Burla Sai Kiranab and
Kandadai Sowjanyaa
aGas Hydrate Division, CSIR–National Geophysical Research Institute (CSIR–NGRI), Hyderabad – 500 007, India. E-mail: psrprasad@ngri.res.in; Fax: +91 40 2717 1564; Tel: +91 40 2701 2710
bAcademy of Scientific and Innovative Research (AcSIR), CSIR–NGRI Campus, Hyderabad – 500 007, India
First published on 7th May 2020
Methane hydrates are promising materials for storage and transportation of natural gas; however, the slow kinetics and inefficient water to hydrate conversions impede its broad scale utilisation. The purpose of the present study is to demonstrate rapid (2–3 h) and efficient methane hydrate conversions by utilising the water molecules confined in the intra- and inter-granular space of silica powders. All the experiments were conducted with amorphous silica (10 g) powders of 2–30 μm; 10–20 nm grain size, to mimic the hydrate formations in fine sand and clay dominated environments under moderate methane pressure (7–8 MPa). Encasing of methane molecules in hydrate cages was confirmed by Raman spectroscopic (ex situ) and thermodynamic phase boundary measurements. The present studies reveal that the water to hydrate conversion is relatively slower in 10–20 nm grain size silica, although the nucleation event is rapid in both silicas. The process of hydrate conversion is vastly diffusion-controlled, and this was distinctly observed during the hydrate growth in nanosize silica.
On the other hand, the clathrate hydrates themselves have many interesting applications, such as fuel gas storage and transportation,4,5 gas separation from mixtures (flue gas),6 water purification,7 and carbon dioxide capture, storage/sequestration and utility.8–12 Typically, about 160–180 volumes (STP) of guest molecules can be encased in the hydrate cages. Their thermodynamic stability conditions are much milder (∼8 MPa and ∼10 °C), unlike liquefied (∼−162 °C and 0.1 MPa) or compressed (∼25 MPa and 25 °C) forms for natural gas.1 Slower formation kinetics and lower water to hydrate conversions are the prime obstacles for its industrial applications.13 Despite inherent risks and associated energy-intensive process in CNG and LNG technologies; they are the preferred methodologies for fuel gas storage & transportation. However, recent experiments on GH formation in porous medium signify that the process of hydrate nucleation and growth could be rapid.14–17 Thus the hydrate formation and dissociation in porous media or sediments should be understood, which is dependent on the medium property. Generally, gas hydrates occur in the form of segregated nodules, lenses, pellets or sheets in fine particle sediments, also as interstitial pore-filling between particles in coarse-grained sediments.18 Natural sediments hosting gas hydrates are generally characterised by mineral particle, organic debris, pore water, mud and clay, etc., so it is complicated to specify the impact of each factor. Because the deposition of sediment particles with different size controls the pore space, the amount of pore water and free gas, etc., and thus the particle size plays an important role in the hydrate morphology and distribution.19,20 Recent gas hydrate explorations in the Indian offshore region have indicated their presence in clay, fine-grained and coarse-grained sand dominated geological locations covering Krishna–Godavari (KG), Andaman and Mahanadi basins.21,22
The silica is a dominant constituent in oceanic sedimentary strata, and its morphology and grain size are also important factors in the hydrate formation mechanism.23–25 Chong et al.,26 examined methane hydrate formation in silica sand of different sizes namely sand-silt (0.063 mm) to granular pebbles (3 mm) and have reported that formation extending over 70 h. On the other hand, the addition of a small amount of silica to H2O–CH4 hydrate forming system has proven to be promoting methane hydrate formation kinetics.24,25,27 Further, series of earlier experiments from our group using hollow silica, comprising of inner void surrounded by a thinner solid shell, have shown several advantages over silica sand.28–31 In particular, the hollow silica helps in achieving rapid and efficient methane conversions with high confined water (4–6 g H2O/g SiO2).28–31 Further, it was established that the hydrate crystal growth morphology in low density (hollow) silica, that we used earlier, distinctly differs with that of silica sand or gels or capillary.32–35 In the former, the hydrate growth, predominantly, occurs on the top of the matrix bed, and the water subsequently is transported through the capillary action to the top of the bed.32 Whereas, in silica sand bed reactors, the hydrate growth reported to occur on the top layer of the sediment bed and subsequently spreads to the deeper parts; thereby the water confined in the interstitials were converted to hydrates.33–35 The solid silica samples of suitable size, are thus more appropriate, in mimicking the natural gas hydrate mechanism and therefore we choose two SiO2 samples having grain sizes 2–30 μm; 10–20 nm, for the present study with different water to SiO2 ratios, representative of typical fine-grained and clay-rich geological environments. The density of these solid silica samples is high, and they completely sink in water in the water column. Therefore intergranular spaces are saturated at lower water content compared to low-density hollow silica.28–31 Methane/carbon dioxide adsorption measurements were conducted earlier using carbonaceous material and reported enhancement in gas storage capacity; however, the studies were inconclusive about the mechanism, namely, was it due to hydrate formation or due to increased adsorption under wet conditions.36–38 Particularly enhanced isothermal (∼275 K) CO2 gas adsorption in ordered mesoporous carbon matrix with higher pre-adsorbed water molecules at 3.6 MPa is very attractive. It is 2.26 time higher than corresponding dry carbon matrix.38 In another report Casco et al.,39 have reported the formation of methane hydrates in the water confined in macro-, meso- and micropores of activated carbon. The authors have shown that the water in macropores was converted to hydrates at 3–4 MPa, while that in meso- and micropores require higher pressure (6–8 MPa). It is worth noticing that the experimental pressure and temperature employed in those studies are well within hydrate stable region. Liu and Liang40 recently reported the methane hydrate formation in a silica gel matrix, particularly, below ice melting temperature. The gas uptake was 0.150 mol mol−1, and the water to hydrate conversion was nearly 100% at 253 K and 6.0 MPa. Faster formation kinetics of 160 min (@263 K) was reported in this system, with gas uptake ratio of 52.96 mol min−1 m−3. Thus, indicating the surface area and morphological properties of the host matrix also play a critical role in the hydrate formation kinetics. On the other hand, the kinetics and water to hydrate conversion in consolidated sediments also critically depend on several parameters such as permeability, degree of water saturation and salinity.41
The present study aims at providing better insight into the mechanism of hydrate formation in SiO2–H2O–CH4 system, by conducting experiments by varying the water in matrix pores. In this study, we used solid silicas, wherein the confined water molecules occupy intergranular space. The driving force (pressure difference between the experimental and the equilibrium pressures at the hydrate nucleation temperature) was kept constant. Additionally, morphological and physical properties of the host sediments could influence the kinetics of hydrate formation and dissociation.38–43 The formation kinetics and ability of confined water into hydrates under the gas-rich environment is studied.
The experimental procedure followed for gas hydrate synthesis has already been described earlier,25,26 and a schematic layout is shown in Fig. 1. Briefly, the central part was an SS-316 cylindrical vessel, which can withstand gas pressures up to 15 MPa, and the volume of the vessel was 250 mL. A cold fluid (water + glycol mixture) was circulated around the vessel with the help of a circulator to bring and maintain the temperature inside the cell at the desired level. A platinum resistance thermometer (Pt100) was inserted into the vessel to measure temperature with an accuracy of ±0.5 K. Pressure in the vessel was measured with a pressure transducer (WIKA, type A-10 for pressure range 0 to 25 MPa with ±0.5% accuracy).
The reactor vessel was first filled with 10 g of silica powder and appropriated amount of water was added drop by drop and the vessel was pressurised (@300 K) with methane gas. The atmospheric gases in the experimental cell were diluted by purging with methane gas before the experiments, and the gas was filled to the desired level using the Teledyne ISCO syringe pump to pressure and temperature outside of the hydrate stability zone. Then, the reactor was isolated from the ISCO pump/gas tank by closing the gas inlet valve. Subsequently, a cold fluid from the chiller was circulated to bring down the temperature of the reactor, and the hydrate formation was detected by a sharp pressure drop at a particular temperature. Experimental set-up for probing the gas hydrate formation/dissociation under isochoric conditions is less expensive and simpler to implement. Thus, we adopted this methodology in our laboratory. The temperature and pressure were logged for every 60 seconds of the time interval. Insignificant head-pressure drop in the reactor observed over a longer duration indicates saturation in hydrate conversion. The bath temperature was slowly increased during the warming cycle. The heating rates during the hydrate formation (cooling cycle) and dissociation (warming cycle) regions were kept as 1–2° h−1.
Run no. | Rw | Nanoshel silica (sample #1) | Aldrich silica (sample #2) | ||||
---|---|---|---|---|---|---|---|
nW (mol) | nCH4 (mol) | % H2O conversion | nW (mol) | nCH4 (mol) | % H2O conversion | ||
Unsaturated samples | |||||||
1 | 0.5 | 0.278 | 0.0735 | 152.02 | 0.278 | 0.0875 | 180.98 |
2 | 0.5 | 0.278 | 0.0705 | 145.82 | 0.278 | 0.0850 | 175.81 |
3 | 0.5 | — | — | — | 0.278 | 0.0860 | 177.88 |
Average value | 0.0720 | 148.92 | 0.0862 | 178.22 | |||
Standard deviation | 0.0021 | 4.38 | 0.0013 | 2.60 | |||
4 | 1.0 | 0.556 | 0.1250 | 129.27 | 0.556 | 0.1700 | 175.81 |
5 | 1.0 | 0.556 | 0.1125 | 116.34 | 0.556 | 0.1725 | 178.39 |
6 | 1.0 | 0.556 | 0.0950 | 98.25 | 0.556 | 0.1635 | 169.09 |
Average value | 0.1108 | 114.62 | 0.1687 | 174.43 | |||
Standard deviation | 0.0150 | 15.58 | 0.0046 | 4.80 | |||
7 | 1.5 | — | — | — | 0.834 | 0.2430 | 167.54 |
8 | 1.5 | — | — | — | 0.834 | 0.2400 | 165.47 |
9 | 1.5 | 0.834 | 0.1425 | 98.25 | 0.834 | 0.2450 | 168.92 |
10 | 1.5 | 0.834 | 0.1525 | 105.14 | 0.834 | 0.2350 | 162.02 |
11 | 1.5 | 0.834 | 0.1625 | 112.04 | 0.834 | 0.2250 | 155.13 |
Average value | 0.1525 | 105.14 | 0.2376 | 163.82 | |||
Standard deviation | 0.01 | 6.89 | 0.0079 | 5.51 | |||
![]() |
|||||||
Saturated samples | |||||||
12 | 2.0 | 1.11 | 0.1875 | 97.13 | 1.11 | 0.1400 | 72.52 |
13 | 2.0 | 1.11 | 0.1925 | 99.72 | 1.11 | 0.1900 | 98.42 |
14 | 2.0 | 1.11 | 0.1950 | 101.01 | 1.11 | 0.1600 | 82.88 |
15 | 2.0 | 1.11 | 0.1950 | 101.01 | 1.11 | 0.1375 | 71.22 |
Average value | 0.1925 | 99.72 | 0.1569 | 81.26 | |||
Standard deviation | 0.0035 | 1.829 | 0.0243 | 12.57 | |||
![]() |
|||||||
Oversaturated samples | |||||||
16 | 2.2 | — | — | — | 1.223 | 0.1550 | 72.87 |
17 | 2.2 | — | — | — | 1.223 | 0.1450 | 68.17 |
18 | 2.2 | — | — | — | 1.223 | 0.1600 | 75.22 |
Average value | 0.1533 | 72.09 | |||||
Standard deviation | 0.0076 | 3.59 | |||||
19 | 2.5 | 1.389 | 0.0525 | 21.73 | — | — | — |
20 | 2.5 | 1.389 | 0.0550 | 22.76 | — | — | — |
21 | 2.5 | 1.389 | 0.0530 | 21.94 | — | — | — |
Average value | 0.0535 | 22.14 | |||||
Standard deviation | 0.0013 | 0.54 |
We calculated the water conversion to hydrates by considering the total water added to the SiO2 matrix, amount of CH4 consumed during the hydrate conversion process and assumed the hydration number as 5.75.28–31 However, in literature, different authors have used different hydration numbers in the range of 5.75 to 6.1. This variation has been mainly attributed due to the cage occupancy factors.1 Throughout this study, we have used the ideal hydration number (5.75) assuming full occupancy of the cages, even though this is about 3% less than the estimated hydration number.44 As seen from Table 1, the water to hydrate conversion factor exceeds 100% for lower Rw ratios. Such anomalous variations in dry or partially wet porous systems were reported earlier by several authors,38,39,45 and precise mechanism remained resolved. However, the hydrate conversion factor quickly reaches to 100% in coarse-grained SiO2 matrix (sample #1) as Rw is 1.5–2.0. Whereas there is a significant contrast for the water conversion factor in fine-grained SiO2 matrix (sample #2), such as it is more than 100% at Rw = 1.5 and very much less than 100% at Rw = 2.0. Casco et al.39 have also estimated the stoichiometry of the methane hydrates synthesised in the confined nano space, of activated carbons, from the adsorbed methane gas and water at different Rw ratios. They have also suggested the stoichiometry to be CH4·(3.6–4.1)H2O, particularly, in undersaturated conditions. Further, they estimate the stoichiometry very close to natural one (i.e., 5.75) in oversaturated samples. Zhou et al.,38 also reported similar observations in CO2 hydrates synthesised in the mesoporous carbon matrix. They also reported a progressive increase in isothermal gas adsorption as a function of increasing Rw ratio in the range 0–2.35. In contrast, the gas adsorption capacity decreases drastically at higher Rw (= 3.55) values. Furthermore, estimated stoichiometry (2.65–3.26 for Rw = 1.55–2.35) in this system also showed more significant deviations from the ideal value of 5.75. Similarly, Zhang et al.45 have also investigated the role of pre-adsorbed water in ordered mesoporous silica system in CO2 gas adsorption. They reported the highest storage capacity (12.8 mmol g−1) at Rw = 2.48, which is 1.86 times higher than the dry (Rw = 0.0) samples. The gas adsorption capacity decreased considerably at higher Rw (= 3.07) value. Thus, these experiments indicate the necessity for critical water content in porous systems for higher gas storage capacity. The main reason for higher gas intake was advocated as hydrate formation because the temperature and pressure conditions of isothermal adsorption measurements were well within hydrate stability regions. At higher water content (higher than the critical value) the gas storage capacity decreases significantly, which could be due to plugging in the gas circulation paths.38,45 On the other hand, the authors45 have also reported a decrease in gas adsorption, compare to dry conditions when Rw = 1.08 and 1.52. Integrated experimental studies using high-pressure differential calorimetry (HP-DSC) along with standard hydrate formation/dissociation measurements could provide useful information to quantify gas uptake due to adsorption and hydrate conversion.
Further hydrate formation was very sluggish at higher Rw ratios. The pressure drop was negligibly small in these systems upon increasing Rw to 3.0 and 2.5 respectively, and this water was enough to convert the samples into a thick paste (see ESI Fig. 5†). Therefore, we reduced Rw to 2.5 and 2.2, respectively, in sample 1 and 2 and conducted another set of experiments. The hydrate conversion factor for sample #1 decreased significantly to about 22.14%. On the other hand, in sample #2, the hydrate conversion is comparable for Rw = 2.0 and 2.2. However, the hydrate formation rate is considerably slower (see the following paragraph). Prolonged gas consumption during the warming cycle is evident in Fig. 3. A plot of CH4 gas content verses the temperature is reconstructed from the logged p, T data using the procedure described earlier.27 The behaviour could be attributed to the apparent blockages in the gas circulation paths due to ice formation during the cooling cycle and was eventually opened up during the warming cycle.45,46 Subsequently considerable amount of gas was adsorbed, resulting in more hydrate formation before the onset of hydrate melting. Such behaviour was observed only in water-saturated systems with Rw = 2.0 and 2.2.
![]() | ||
Fig. 3 A trajectory showing the variations in the amount of methane gas at different temperatures during the cooling (black) and warming (red) cycles for SiO2 (sample #2) with Rw = 2.2. |
The hydrate formation is a phase transformation, involving solid hydrates, from starting water and gas phases. To estimate the gas consumption and the kinetics in hydrate conversion, one needs time information on its onset and induction time. It is well known that the hydrate formation is an exothermic process accompanied by a small degree of increase in the local temperature and also associated with a decrease in the gas pressure. However, it has not always been the case in an isochoric process. Ideally, spectroscopic methods are suitable for detecting a state change such as hydrate conversion. Nevertheless, a distinct change in pressure decrease in conjunction with constant and/or temperature spike over sometime within the hydrate stability zone is considered as the onset point for hydrate nucleation.47 Fig. 4 shows the CH4 gas uptake per mole H2O in sample #1 (a) and sample #2 (b) in 10 h after the hydrate onset. These gas uptake kinetics plots are the average of repeat measurements tabulated in Table 1 with different Rw values. Initially, the gas uptake is rapid and eventually reaches a plateau, indicating the saturation in hydrate conversion. It is worth noting that the maximum methane gas uptake in sI hydrate, consisting of 8 empty cages formed from 46H2O molecules, is 0.173 (shown as a grey line in Fig. 3 and 4). The formation kinetics in both the systems is rapid, and 85–90% of gas uptake has occurred within 2–3 h. Further, in sample #1 (see Fig. 3a) the maximum gas uptake approached 100% of the permissible limit for Rw = 1.5–2.0 and this is an indication of pore water saturation state. The gas uptake understandably is more than 100%, and this mainly could be due to gas adsorption on wet/dry SiO2. Bagherzadeh et al.48 have shown through molecular simulations with a limited amount of water that a meniscus is formed between silica surfaces separated by 4–6 nm. Further, their studies on water and gas molecular distribution indicated that the water number density in the layers adjacent to the silica is higher, and these layers are more structured and less mobile compared with water layers far from the surface. Additionally, the concentrations of the gases are significantly higher at the liquid and silica interfaces than in other locations in the gas phase. In other words, higher free-gas adsorption was evidenced on a wet surface. Further, Bagherzadeh et al.48 speculated that “enhanced concentration of gases at the interface along with the extended contact area (curved meniscus compared with flat interface in the absence of silica surfaces) between water and guest molecules at the meniscus as well as lesser mobility of water molecules near the silica surface may provide a mechanism for the heterogeneous nucleation of the clathrate hydrate in a water-wetting porous medium”. Therefore, presences of the mineral body such as silica facilitate in the creation of nucleation sites and kinetically promote clathrate hydrates,49 However, they will not alter the space group of hydrates. Indeed the experiments conducted H2O–CH4 system with a small amount of silica supports such a mechanism.24,25,27 However, depending on the extent of hydrophobic or hydrophilic nature of the matrix materials, the lattice constant of the hydrates can show a small perturbation.50 The experimental studies reported in current research also support such conjecture.31,47 Nesterov et al.,47 investigated the role of certain oxides such as silicon, aluminium, zinc, titanium and manganese as active promoters for natural gas hydrates. Striking observations from their studies are; the presence of such oxides along with the surface functional groups, particularly, carbon-containing groups (e.g., carbonates) causes a decrease in the induction time of hydrate formation. The size of inducer particles should be much larger than the size of the hydrate nucleus. Induction time of gas hydrate formation decreases with an increase in the size of the promoter particles. From Fig. 3a, it is seen that the gas uptake, for sample #1, in the initial 0.5–1.0 h duration is only 40–50% and later on the gas uptake became slower for all Rw values. The time taken for 90% of the maximum (t90) gas uptake varies in the range of 3.6–2.3 h for Rw = 0.5 to 2.0. On the other hand, the gas uptake in nanosized SiO2 (sample #2, Fig. 4b) is comparatively slower. In fact, the gas uptake is 93, 77 and 67% of the maximum for under-saturated systems (Rw = 0.5, 1.0, 1.5) in 3.6 h duration; while it is ∼50% for saturated system with Rw = 2.0 or 2.2. The hydrate growth times in our experiments were restricted to ∼20 h; however, the hydrate growth may take place beyond these time scales. Our results are confirmatory to the literature data that the growth in fine-grained sediments is generally slower. This fact is further illustrated in Fig. 5 through the long-time behaviour of gas intake in two samples. The maximum hydrate conversion in micro-grained SiO2 (sample #1) is >98% (for saturated) and ∼20% (for oversaturated) on longer time scales; while the hydrate growth in nanosized SiO2 is slower but eventually reached higher hydrate conversions.
![]() | ||
Fig. 5 The long duration kinetics plot of methane uptake during the hydrate formation in both silica samples. The corresponding Rw values are shown along the curves. |
As said, the formation of gas hydrates is a phase transformation to solid-state from liquid (H2O) and gaseous guest molecules. The process of methane hydrate formation can be fitted using the well-known classical crystallisation Avrami model for phase change.51–53 The following equations govern the model:
α = 1 − exp(−k(t)n) | (1) |
S. no | Rw | Nanoshel silica (sample #1) | Aldrich silica (sample #2) | ||
---|---|---|---|---|---|
n | k (mmol per mol H2O per min) | n | k (mmol per mol H2O per min) | ||
1 | 0.5 | 1.26 | 8.0 × 10−7 | 1.11 | 16.1 × 10−7 |
2 | 1.0 | 1.26 | 9.2 × 10−7 | 1.22 | 9.3 × 10−7 |
3 | 1.5 | 1.64 | 3.2 × 10−7 | 1.48 | 2.8 × 10−7 |
4 | 2.0 | 1.30 | 8.5 × 10−7 | 1.51 | 2.1 × 10−7 |
5 | 2.2 | — | — | 1.23 | 4.4 × 10−7 |
6 | 2.5 | 1.08 | 4.9 × 10−7 | — | — |
From the previous studies, it is well established that the Raman spectroscopy has unambiguous spectral signatures for the methane molecules in free/dissolved state or encased in clathrate cages.44,54 The methane molecules encaged in the clathrate phase have distinguishable Raman shifts at 2915 and 2904 cm−1 corresponding to methane in two different cages, whereas, in the gas/dissolved phase it shows a single mode at 2918 cm−1. Such diverse signatures are useful in establishing the molecular configuration of methane molecules. We collected the Raman signatures of the hydrates grown in these two silica matrices and the detailed procedure adopted for such measurements is described earlier.44 From Fig. 6, we can realise that the recorded spectrum in the C–H stretching mode region on both the samples is comparable to the methane hydrate features. The Raman shifts at 2915 and 2904 cm−1, full width at the half-maximum intensity and their relative intensity ratio are comparable with methane molecules encased in sI cages.44,54
Contemporary literature states that the hydrate grain size is one of the parameters that portray phase boundary (PB) points of methane hydrates. In particular, methane hydrates formed in narrow pores (∼10 nm) are less stable than bulk hydrates and its PB point shifts to inhibition side. Furthermore, a progressive shift towards the bulk PB point has been reported as a function of increasing hydrate grain size.14 On the other hand, pure hydrates of 100 μm or bigger have no effect and their phase boundary points coincide well with bulk hydrates.14,55,56 Observed phase boundary behaviour, as shown in Fig. 7, in two silica samples investigated presently, agrees with this scenario. In all the experiments, the hydrate samples were dissociated with a heating rate of 1–2° h−1 and the reported values are the averages of repeat measurements. In under-saturate systems (Rw = 0.5) the PB point (279.8 ± 0.5 K, & 7.4 MPa) is lower by about 3 K in nanosized silica (sample #2); while that for sample #1 (283.8 ± 0.5 K, & 7.7 MPa) closely matches with bulk hydrates. The PB points shifts progressively towards 282.8 ± 0.5 K (7.7 MPa) as Rw ratio is increased to 2.2. We conjecture this could be due to the formation of thicker water layer (and therefore hydrate) within the nanosized silica. However, no such systematic shift is evident in micron-sized silica (sample #1). Still, the PB point has shifted by about 1.5–2.0 K to the right of bulk hydrate, i.e., towards the promotion side, and reasons could be multifold. Generally, the location of the thermal sensor, particularly when heat transfer is non-uniform, and faster hydrate dissociation rate with insufficient equilibration time. However, the primary goal of these measurements is to ensure the hydrate formation in water + silica system, and there is a measurable shift in their PB points. Tohidi et al.,57 have demonstrated that insufficient equilibration time in step heating could cause measurable drift in the dissociation temperature. Further, Park et al.,58 have reported that the PB points of methane hydrates, within the silica matrix (dry water), move to the promotion side by a small increment and this has been mainly attributed due to heat transfer efficiency of the matrix material.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ra01754j |
This journal is © The Royal Society of Chemistry 2020 |