1. Introduction
Global warming has become one of the main issues in the public domain, and its importance is steadily increasing. The temperature increase can be in good part ascribed to the increase in anthropogenic carbon emissions, as a result of mankind’s massive dependence on carbon-based fossil fuels for the majority of energy production. In order to prevent or at least decelerate the advancement of global climate changes, the carbon footprint of human activities needs to be reduced [
1]. Combustion of cleaner fuels, increased use of renewable energies and a focus on CO
2 capture technologies can all contribute to this ultimate goal [
2].
Natural gas use in power generation provides low capital costs and favorable heat rates, for this reason is an attractive fuel [
3]. Moreover, the decarbonisation in the energy sector indicates a process of progressive decrease over time of the H/C ratio in the different sources of energy supply. Methane, with a ratio of 4, presents the highest value among fossil energy sources.
In addition to conventional gas and oil wells, unconventional deposits such as shale gas, tight gas, or coal bed methane deposits have been studied as huge gas resources. Nevertheless, the largest amount of methane worldwide is trapped in the permafrost and oceanic sediments, as natural gas hydrates (NGH).
Gas hydrates are solid crystalline compounds, in which small gas molecules are trapped inside cages of hydrogen-bonded water molecules. They form when the guest gas (i.e., CH
4) and water coexist at sufficiently low temperatures and high pressures. There are important results in literature about experimental and molecular simulations, which describe in detail the mechanisms of gas hydrate nucleation and crystal growth [
4,
5] thermal and mechanical properties [
6,
7], mass and heat transport processes [
8,
9,
10].
One volume of solid methane hydrate contains 164 volumes of gaseous methane at standard temperature and pressure (STP), which is equivalent to room temperature at 18 MPa [
11]. The amount of organic carbon contained in NGH reserves around the globe is estimated to be twice the amount contained in all currently recoverable worldwide conventional hydrocarbon resources [
12]. Numerical estimates of world NGH deposits range from the first published number of 3053 × 10
15 Sm
3 [
13] to more conservative values of 15 × 10
15 Sm
3 [
14], 2.5 × 10
15 Sm
3 [
15] and 120 × 10
15 Sm
3 [
16].
NGH exist in those areas where both biogenic or thermogenic gas sources occur under large, cold masses providing high pressure values and low temperatures. These conditions occur beneath quite thick frozen soils such as arctic permafrost. Several areas in the Arctic show potential for having gas hydrate accumulations, in North America, in Russia, in Norway.
Another primary hydrate stability zone are regions in the subsea floor sediment, where temperatures and pore pressures are suitable for stable NGH [
17]. NGH marine deposits are generally found along the continental slopes and shelves with abundance of organic matter supply and with water depths exceeding 500 m. In the deepest ocean regions, below around 2000 to 3000 m, hydrates are scarcely found because of the limited availability of methane that in the ocean is produced by microbial decomposition of organic matter sinking down from the sunlit zone near the surface. Conversely, no NGHs are found in marginal seas because of the insufficient pressure for their stability [
18]. Recent programs of gas hydrate field drilling have assessed gas hydrate accumulations in the marine sands in the Cascadia Margin, Gulf of Mexico and Nankai Trough, and also other types in offshore India, Malaysia and Korea [
12].
The natural gas production from oceanic and permafrost sediments is currently being performed using methods such as injection of hydrate inhibitors, depressurization and thermal stimulation.
The chemical inhibition method seeks to displace the NGH local conditions of pressure and temperature outside the thermodynamic conditions of hydrate stability zone through the injection of a liquid inhibitor from the surface down to NGH bearing layers [
19].
The depressurization method dissociates methane hydrate by decreasing the pressure in the wellbore drilled inside the NGH sediments [
20]. When pressure decreases, NGH becomes thermodynamically unstable and dissociates absorbing sensible heat from the surrounding sediments as well as geothermal heat flow. Depressurization method has been extensively investigated at laboratory scales [
21].
The thermal stimulation method decomposes NGH by raising temperature beyond the hydrate stability zone. Most of the thermal dissociation tests have been performed via the injection of heat transfer fluids, such as steam, into the hydrate sediment [
22]. Other proposed thermal stimulation techniques are microwave heating, electro-magnetic heating, and in-situ combustion [
19]. These methods are considered more promising than the conventional heated fluid injection methods, thanks to the direct energy contribution necessary for the dissociation at the hydrate zone [
17].
The in-situ combustion process consists in locating a point heat source in the hydrate region, via direct combustion of a liquid fuel and an oxidant [
23,
24].
Recovering methane from hydrate reservoirs can be also combined with simultaneous permanent CO
2 sequestration [
23,
25,
26]. The CO
2 injection into the NGH sediments causes the gaseous CH
4 release and the formation of CO
2 hydrate, even if global conditions of pressure and temperature remain within the CH
4 hydrate stability zone. This gives the opportunity to simultaneously recover an energy resource while entrapping a greenhouse gas.
Many experimental and numerical studies have been performed to understand the CH
4-CO
2 replacement [
27]. Baig et al. in [
28] stated that conversion of CH
4 hydrates in CO
2 hydrates occurs according to two mechanisms: new CO
2 hydrate formation starting from gaseous CO
2 and free water surrounding hydrates and a direct conversion in the solid state.
The solid-state exchange is slow because of mass transport limitations. The formation of new hydrates is a faster process since it releases heat, assisting the in situ dissociation of the methane hydrate. This is faster, because heat transport is very rapid in these systems [
29].
There is a sound body of knowledge in literature about energy potential of NGH reservoirs, with studies in laboratory settings and some recent field trials, deeply reviewed in [
12]. According to a hydrate research roadmap by the U.S. National Methane Hydrate R&D program [
30], the key research topics to be investigated are the development technologies of sampling, characterization, and production, the investigation on the role of gas hydrate in the global climate change. Some production technologies have been recently proposed, such as the use of horizontal wellbore [
31].
The present paper aims at addressing the issue of a production technology development, proposing a membrane-based CH
4-CO
2 replacement process for methane recovery and CO
2 storage in NGH. Production control during NGH exploitation is considered by the NGH scientific community a key priority which involves appropriate drilling and hazard monitoring technologies [
32]. Nevertheless, to the best of our knowledge, studies on the CH
4-CO
2 replacement process, which take into account the separation efficiency and propose a process to recover pure methane, are not present in literature.
In the present paper, an integrated process with a membrane section is proposed and analyzed from an energy and environmental point of view. The use of mechanical devices such as membranes [
33] is proposed to separate the CO
2 and CH
4 which have not taken part to the replacement process. In this way, pure methane can be recovered and CO
2 can be recirculated in the sediment.
The methane recovery and the CH4-CO2 replacement in the hydrate phase is studied experimentally in a 1 L lab scale reactor. The CH4 hydrates are formed in a matrix resembling geological properties of natural seafloor and are subsequently dissociated by a CO2 stream supplied in the upper part of the reactor. Experimental results will be used as input for energy and environmental evaluations to prove the sustainability of the process.
3. Experimental Results
In this section, results of the experimental tests carried out in laboratory are presented and discussed.
3.1. Hydrate Formation
An experimental test with a hydrate saturation equal to 10% was carried out to obtained experimental data useful for energy and environmental analysis of the membrane based integrated process. A 10% hydrate saturation was chosen because a good solution to well simulate the hydrate deposit conditions and effectively control temperature is to maintain a low sand pore saturation. In
Table 3 the experimental conditions and the obtained results are showed.
Test 1 was carried out at 10% hydrate saturation with an initial pressure and temperature respectively of 42.0 bar and 1.0 °C. To reach the initial pressure equal to 42.0 bar, 0.2307 moles were injected. The equilibrium pressure achieved at the end of the hydrate formation phase was 28.3 bar. At the end of the test, the dissociation pressure, which was reached inside the reactor because of the hydrate dissociation, was equal to 16.5 bar corresponding to CH
4 moles in the hydrate phase equal to 0.0769. Assuming an ideal hydration number equal to 6 and considering the negligible methane dissolved in water at system conditions (0.0012 mol kg
H2O−1), a total of 0.948 moles of water are required in the hydrate formation.
Figure 6 shows temperature profiles when the methane hydrate start forming in Test 1, carried out with saturation of 10%.
Figure 6 shows the temperature vs. time diagram, in which the temperature increase due to the exothermic nature of the formation is visible. Temperature increase due to hydrate formation is bigger in the lower part of the reactor (corresponding to thermocouples T11 and T16), while in the upper part (T07), temperature starts increasing previously but it then decreases sharper. The thermocouple T02 is positioned in the headspace, where only gas is present and there is no sand, temperature remains quite constant, suggesting no hydrate formation in this area.
The temperature increase in the middle part of the diagram shows the exothermic nature of the hydrate formation process, and this is particularly noticeable in a small scale reactor like this. In order to avoid that thermodynamic conditions strays from the equilibrium curve, it is necessary to cool the reactor and then reducing the peak amplitude, achieving the previous temperature value as soon as possible. Monitoring temperature is also necessary to obtain the same value both at the beginning and at the ending of the test.
CH
4 hydrate formation starts in the T07 region and, and then it moves towards the T11 and T16 regions. This is consistent with other experimental results. Zhou et al. [
38], in a 60 L reactor, found that hydrate formation starts from the gas-liquid interface in the top head space and then the hydrate formation front improves downward, thanks to the larger free water availability in the lower part of the reactor. From this comparison, it appears that temporal and spatial distribution of the CH
4 hydrate formation does not depend on the reactor’s size.
Figure 7 reports pressure profile during the formation and dissociation process for Test 1.
The hydrate formation causes a pressure decrease due to the methane adsorption inside the clathrate solid structure. Once pressure is stabilized, the reactor was opened to eject all gaseous CH
4 present inside. Then the vessel was closed again and pressure started increasing, due to the hydrate dissociation and the releasing of methane from solid to gaseous phase. In both situations the pressure variation happens in the first hours, then it assumes an asymptotic trend until reaching the equilibrium value.
Figure 8 shows the pressure-temperature graph, in which the theoretical thermodynamic conditions and the experimental one are compared.
The distance between the theoretical line and experimental points is greater in the upper area, corresponding to the hydrate formation beginning, or when reaction has not involved yet the whole reactor internal volume. The temperature increase due to the formation reaction moves the experimental conditions towards the equilibrium line and during hydrate formation there is an overlap between T07, T11 and T16 points and the equilibrium graph until the hydrate formation is completed.
3.2. CO2 Replacement
The CO
2 replacement was carried out by depressurization of the internal environment; the substitution was then started by bringing the reactor from the value of 28.3 bar to a pressure of 20.4 bar. In
Table 4 values of the experimental parameters are shown.
CO
2 injection started and continued for 1620 s. After that in the reactor 0.119 CO
2 moles are present. At the end of the CO
2 injection, the gas composition is 75.4% CO
2 and 24.6% CH
4.
Figure 9 shows the temperature trends during the replacement test, starting from internal reactor conditions immediately after the methane hydrate formation.
The test starts at an average internal temperature equal to 1.5 °C. Over time, the first temperature peak occurs in the higher part of the reactor at T07, but it is the lowest temperature value reached, with an increase of about 2 °C. The higher temperature value is reached in the lower area at T16, which is the second peak in temporal order. This suggests that the zone involved firstly with an exothermic process is the lower part where the CO2 injection point is located. This is also an area of the reactor rich in water, so the higher temperature increase with respect to the other thermocouples could suggest new hydrate formation rather than a replacement process. A second phase follows in the upper part of the sediment, near the gaseous headspace. The last peak occurs in T11, after few minutes, in the central area of the reactor.
The comparison of the obtained temperature peaks in the three different areas of the sediment suggests that a CH4-CO2 replacement took place in the central part of the reactor (T07 and T11) where the formation of CO2 hydrates and the simultaneous dissociation of CH4 hydrates avoided the formation of high temperature peaks.
Figure 10 shows the trend of the pressure during the test. In the left side of the graph, is visible the pressure drop made to move below the stability conditions for methane hydrates, while remaining in the CO
2 hydrate formation area. The constant pressure drop due to the formation of CO
2 hydrates, ends with the reactor opening, carried out to leave only the hydrates formed and the gas contained therein inside. The pressure increase in the right side represents the dissociation process of hydrates and the release of CH
4 and CO
2 contained in them.
Once the CO2 fluxing was completed, the reactor was closed and a sample of gas was withdrawn to make the gas chromatographic analysis. In this way the concentration of CO2 in the reactor was determined. A second sample was collected and analyzed at the end of the dissociation phase to evaluate the new concentrations of the two gases. Knowing pressure and temperature conditions, it was therefore determined both the number of moles of CO2 present before the starting of the replacement phase, and the number of moles stored in form of hydrate.
The CH
4-CO
2 hydrate dissociation causes a pressure increase from 1 to 6.05 bar (see
Table 3); the gas sample analysis proves that 78.47% of entrapped gas into hydrates is CO
2, while the remaining part is CH
4. The total gas moles which have formed hydrates are 0.043 mol (0.034 mol of CO
2 and 0.009 mol of CH
4).
Two parameters are calculated for the described process, in order to quantify the CO
2 sequestration and the CH
4-CO
2 exchange. Capture efficiency is the ratio of the CO
2 moles stored in the system through CO
2 hydrate formation to the total injected moles (Equation (3)):
The exchange factor is the ratio of the CO
2 moles sequestered in the hydrate phase to the CH
4 moles produced from the hydrate phase and it gives an extent of the replacement process (Equation (4)). In this case, it is equal to 0.54.
At the end of the process, 0.063 CH4 moles were released from the system and 0.034 CO2 moles were captured in the hydrate phase.
4. Energy and Environmental Evaluations
An energy and environmental evaluation was carried out considering that the replacement process is integrated with a membrane-based separation process, used for the purification of the methane released by the replacement process and for the recirculation of the CO2 separated.
The evaluation is based both on the experimental data obtained in laboratory and on theoretical calculations, in which the energy content associated to the recovered methane is compared to the energy consumption spent for its recovery and for CO2 sequestration.
On the membrane separation section, the following assumptions were made: (i) the composition of the gas flow sent to the membrane separation section is exclusively formed by CO2 and CH4; (ii) with a separation efficiency of 0.98, the retentate gas has a composition in volume of 3% CO2 and 97% CH4; (iii) the permeate is pure CO2, which is recompressed and recirculated in the vessel to take part in the replacement process.
Since the methane heat of combustion is 50,000 kJ kg−1, the energy content of the recovered methane (0.063 mol) is 50.40 kJ. Nevertheless, the CH4-CO2 replacement process presents energy consumption due to: (i) CO2 compression for its injection in the sediment, (ii) compression of the separated CO2 to be recirculated in the sediment. These contributions were evaluated considering the following assumptions. The compression work was calculated considering the compression as a multi-stage adiabatic process with a compression efficiency of 85% and inter-cooling. The energy consumption for the cooling operations, in which the environment can be used as a heat sink, was considered zero. CO2 is assumed to be compressed from atmospheric pressure to 42 bar, with an initial temperature of 298 K (25 °C), through a 3-stage compression, with inter-cooling until 293 K (20 °C) with the environment as a heat sink. The work of compression for the injected CO2 is equal to 1.91 kJ (361.93 kJ kgCO2−1).
After passing the membrane separation section at 3 bar, the permeated CO
2 is recirculated to be re-injected in the hydrate-bearing sediment inside the reactor. The permeated CO
2 is compressed to 42 bar starting from 3 bar and an initial temperature of 277 K, through a 2-stage compression with inter-cooling until 293 K (20 °C) with the environment as heat sink. The work of compression for the recirculation the permeated CO
2 is equal to 0.48 kJ (252.81 kJ kg
CO2−1). For the described contributions, considering a global efficiency of 0.46, the primary energy consumption is equal to 2.39 kJ, equivalent to the combustion of 0.006 moles of CH
4. The primary energy consumption for the injection of a CO
2 mass unit (kg
CO2) is equal to 1336.4 kJ kg
CO2−1, which corresponds to the combustion of 1.671 moles of CH
4. Results are summarized in
Figure 11. The ratio between the spent energy and the stored energy in the recovered CH
4 is 4.75%.
The carbon footprint of the CH
4-CO
2 exchange process is calculated as a balance of the CO
2 produced in the process and the CO
2 stored in system. The CO
2 produced in the process is due to: (i) combustion of the recovered CH
4; (ii) CO
2 related to the primary energy consumption; (iii) CO
2 from the membrane gas retentate. As the CH
4-CO
2 carbon ratio is 1:1, the stoichiometric combustion of 0.063 recovered CH
4 moles produces the same number of CO
2 moles. On the same basis, the primary energy consumption associated to the compression of the injected CO
2, brings to the production of 0.007 moles of CO
2. CO
2 from the membrane separation is 0.001 mol (
Figure 12). The amount of CO
2 moles (0.075 moles) is calculated adding the moles stored in hydrates to the moles dissolved in water.
Also in this second case we can speak of permanent storage, as the thermodynamic conditions of the environment will lead to the formation of new hydrates. The carbon footprint therefore results in an estimated negative value, equal to 0.004 moles sequestrated, thus proving that the exchange process acts as carbon storage sink with consequent beneficial environmental effects.
5. Conclusions
The present paper proposes a process for pure methane recovery and carbon storage in natural gas hydrate sediments based on membrane separation. The CH4-CO2 exchange process was studied in a simulated hydrate-bearing sediment in a 1 L lab reactor, where the CH4 hydrates are formed in a quartz sand matrix, partially saturated with water. The gaseous flow, which results from the CH4-CO2 exchange process, is a CH4-CO2 mixture. The separation of the two components is proposed to be carried out by a membrane-based separation section.
Test of artificial reproduction of natural gas hydrates with a hydrate saturation of 10% show that the hydrate formation goes downward, to the lower part of the reactor. From comparison with literature data, temporal and spatial distribution of the CH4 hydrate formation is not affected by the reactor’s size. The CO2 replacement shows a capture factor of 0.28 and an exchange factor equal to 0.54.
Energy and environmental evaluations were carried out to assess the viability of the process. The energy analysis takes into account the process energy costs, related to the compression and heating of the injected CO2 and the compression of the CO2 separated in the membrane and recirculated in the sediment. Results show that the spent energy constitutes 4.75% of the energy stored in the recovered methane.
The carbon footprint of the CH4-CO2 exchange process is calculated as a balance of the CO2 produced in the process and the CO2 stored in system. Results provide an estimated negative value, equal to 0.004 moles sequestered, thus proving the environmental benefit of the exchange process. Coupling the CH4-CO2 exchange process in the methane hydrate reservoirs with a membrane based separation process to obtain pure gaseous methane is proved to be viable.
Further experimental studies to improve exchange and capture efficiency are planned. Geological and mechanical properties of the reproduced hydrates reservoir will be monitored, especially during the dissociation and exchange process.