Search

Methane hydrates (MHs) have been considered as the future source of energy because of its vast resource volume and high energy density. Energy recovery from MH-bearing sediments has attracted intensifying research activities. Fundamentally, heat transfer, fluid flow through porous media, and the kinetics of hydrate reaction are the three key processes controlling the behavior of MH dissociation and the associated fluid production. Earlier studies have suggested that heterogeneous spatial distribution of SH is inevitable in MH-bearing samples synthesized in laboratory. In this paper, we extend our study to analyze numerically the simulation results from the two realizations of the samples (homogeneous and heterogeneous) to identify differences in the fluid production and to determine if they are sufficiently different. Additionally, we conduct a sensitivity analysis and a statistical analysis on the key transport and kinetic rate parameters that could affect hydrate dissociation and fluid production in the context of a heterogeneous hydrate-bearing sample, in an effort to provide insights that could lead to improved designs for laboratory experiments and (possibly) field applications. Our results suggest that the approximation of an artificial hydrate-bearing core with heterogeneous phase saturations by an assumption of uniform phase saturation distributions results in practically similar fluid production profile except for the very early stage with maximum 20.0% deviation in the water production. From the sensitivity and statistical analysis, we determine that gas production depends strongly on the kinetic rate constant, Kd0 and the composite thermal conductivity of the hydrate-bearing sediments, λθ; while, water production is very sensitive to Kd0 and the absolute permeability of the sandy medium, k. Understanding the effect of phase heterogeneity and the relative importance of key parameters on the production behavior of hydrate-bearing sediments could provide basis for novel production technologies that lead to enhanced gas production and energy efficiency in the energy recovery process.

Laboratory-created samples of methane hydrate (MH)-bearing media are a necessity because of the rarity and difficulty of obtaining naturally-occurring samples. The hypothesis that the inevitable heterogeneity in the phase saturations of the laboratory samples may lead to unreliable and non-repeatable results provided the impetus for this study, which aimed to determine the conditions under which maximum uniformity can be achieved. To that end, we designed four experiments involving different multi-stage cooling regimes (in terms of their duration and number of stages) to induce MH formation under excess-water conditions. In the absence of direct visualization capabilities, we analysed the experimental results by means of numerical simulation, which provided high-resolution predictions of the spatial distributions of the phase saturations in the cores and enabled the estimation of the parameters controlling the kinetic MH-formation behaviour through history-matching. Analysis of the numerical results indicated that, under the conditions of the experiments and with the design of the reactor, significant heterogeneities in phase saturation distributions were observed in all cases, leading to the conclusion that it is not possible to obtain cores with uniform phase saturation. Additionally, contrary to expectations, heterogeneities increased with the number of cooling stages and the duration of cooling, and this was attributed to imperfect insulation of the upper part of the reactor. A set of simulations involving perfect insulation of the reactor top confirmed the validity of this assumption: (a) predicting the formation of high-uniformity MH-bearing cores that became more homogeneous as the number of cooling stages and the length of the cooling period increased; and (b) providing important information for the improvement of the standard design of the experimental apparatus for the laboratory creation of MH-bearing cores using the excess water method.

Methane Hydrates (MHs) are a promising energy source abundantly available in nature. Understanding the complex processes of MH formation and dissociation is critical for the development of safe and efficient technologies for energy recovery. Many laboratory and numerical studies have investigated these processes using synthesized MH-bearing sediments. A near-universal issue encountered in these studies is the spatial heterogeneous hydrate distribution in the testing apparatus. In the absence of direct observations (e.g. using X-ray computed tomography) coupled with real time production data, the common assumption made in almost all numerical studies is a homogeneous distribution of the various phases. In an earlier study (Yin et al., 2018) that involved the numerical description of a set of experiments on MH-formation in sandy medium using the excess water method, we showed that spatially heterogeneous phase distribution is inevitable and significant. In the present study, we use as a starting point the results and observations at the end of the MH formation and seek to numerically reproduce the laboratory experiments of depressurization-induced dissociation of the spatially-heterogeneous MH distribution. This numerical study faithfully reproduces the geometry of the laboratory apparatus, the initial and boundary conditions of the system, and the parameters of the dissociation stimulus, capturing accurately all stages of the experimental process. Using inverse modelling (history-matching) that minimized deviations between the experimental observations and numerical predictions, we determined the values of all the important flow, thermal, and kinetic parameters that control the system behaviour, which yielded simulation results that were in excellent agreement with the measurements of key monitored variables, i.e. pressure, temperature, cumulative production of gas and water over time. We determined that at the onset of depressurization (when the pressure drop – the driving force of dissociation – is at its maximum), the rate of MH dissociation approaches that of an equilibrium reaction and is limited by the heat transfer from the system surroundings. As the effect of depressurization declines over time, the dissociation reaction becomes kinetically limited despite significant heat inflows from the boundaries, which lead to localized temperature increases in the reactor.

We analyse numerically an earlier experimental study that involved the formation of methane hydrates by the excess water method in a small reactor filled with a sandy porous medium, and seek to address questions about the type of the hydration reaction and the phase heterogeneity in the resulting hydrate-bearing sand. Using a fine discretization describing the reactor assembly, the experimental process is faithfully replicated numerically. The multi-stage process of hydrate formation is subdivided in 7 steps. The experimental data from the continuously-monitored pressure and temperature during each step are used for comparison against the numerical predictions, the identification of the dominant processes and the determination of the associated parameters through a history-matching process that minimizes deviations between observations and simulation results. The results of this first-ever study on this subject demonstrate unequivocally that the hydration reaction is a kinetic (as opposed to an equilibrium) process, and that the spatial distributions of the various phases (aqueous, gas and hydrate) at the end of the formation process are strongly heterogeneous. This has serious implications in simulation studies of hydrate dissociation that assume uniform initial phase saturation distributions. The history-matching process indicates that (a) the system behaviour is sensitive to some flow parameters (porosity and irreducible water saturation) only during the first water injection, (b) it is insensitive to the sand intrinsic permeability during all steps of the study, and (c) thermal processes dominate after the first water injection, yielding estimates of the thermal properties of the sand and of time-variable key parameters of the kinetic reaction.

Numerical simulation is essential for the prediction and evaluation of hydrocarbon reservoir performance. Numerical simulators developed for the description of the behavior of hydrates under production and the corresponding flow of fluids and heat accounting for all known processes are powerful, but they need validation through comparison to field or experimental data in order to instill confidence in their predictions. In this study, we analyze by means of numerical simulation the results of an experiment of methane hydrate dissociation by thermal stimulation in unconsolidated porous media heated through the vessel walls. The physics captured by the model include multicomponent heat and mass transfer, multiphase flow through porous media, and the phase behavior of the CH4 + H2O system involved in methane hydrate formation and dissociation. The set of governing equations consists of the mass and energy conservation equations coupled with constitutive relationships, i.e., the dissolution of gas in H2O, relative permeability and capillary pressure models, composite thermal conductivity models, and methane hydrate phase equilibria. The model geometry describes accurately the hydrate reactor used in a recent experimental study investigating methane hydrate dissociation behavior [Chong et. al. Appl. Energy 2016, 177, 409-421]. The cumulative gas production is estimated and validated against three tests of experimental data involving different boundary temperatures, showing a good agreement between observations and numerical predictions. The predicted evolution of the spatial distributions of different phases over time shows that hydrate dissociation progresses inward from the reactor boundary to the center, methane gas accumulates to the top of the reactor because of buoyancy, and water migrates down to the bottom of the reactor because of gravity. A sharp hydrate dissociation front is predicted, and the estimated location of hydrate dissociation front suggests a linear relationship with the square root of time. A sensitivity analysis on the thermal conductivity of sand under fully saturated conditions is conducted to elucidate its effect on the gas production behavior. In addition, the energy efficiency ratio computed from the simulation of this boundary-wall heating technique varies from 14.0 to 16.2. Deviations between observations and predictions of the evolution of the temperature profile are attributed to initial heterogeneous distribution of the hydrate phase in the hydrate reactor.