Academic literature on the topic 'Permafrost – Thermal conductivity'

Soil hydrothermal dynamics are a crucial parameter for understanding the internal physical conditions of the active layer in permafrost regions. It is very difficult to obtain data in permafrost regions, especially on the Qinghai–Tibet Plateau (QTP). Land surface modes (LSMs) provide an effective tool for soil hydrothermal dynamics. However, it is necessary to evaluate the simulation performance before using them. Here, we used two in situ sites along with the latest version of the Community Land Model (CLM5.0) to evaluate the simulated performance in the soil hydrothermal parameters of the model in permafrost regions on the QTP. Meanwhile, the effects of soil properties, thermal roughness length, and the freeze–thaw process on the simulation results were investigated. The results showed that CLM5.0 can capture the dynamic changes in soil hydrothermal changes well in permafrost regions on the QTP. Soil moisture and thermal conductivity were more sensitive to soil properties and the freeze–thaw process, while the thermal roughness length had a greater effect on soil temperature. Notably, although we improved the soil properties and thermal roughness length, there were still some errors, especially in the soil moisture and soil thermal conductivity. It may be caused by inappropriate hydrothermal parameterizations of the model, especially the soil thermal conductivity, hydraulic conductivity, unfrozen water scheme, and snow schemes. There is an urgent need for collaboration between experts in permafrost science, hydrological science, and modelers to develop the appropriate schemes for permafrost regions and enhance the LSMs.

Abstract. Seasonally variable thermal conductivity in active layers is one important factor that controls the thermal state of permafrost. The common assumption is that this conductivity is considerably lower in the thawed than in the frozen state, λt/λf < 1. Using a 9-year dataset from the Qinghai–Tibet Plateau (QTP) in conjunction with the GEOtop model, we demonstrate that the ratio λt/λf may approach or even exceed 1. This can happen in thick (> 1.5 m) active layers with strong seasonal total water content changes in the regions with summer-monsoon-dominated precipitation pattern. The conductivity ratio can be further increased by typical soil architectures that may lead to a dry interlayer. The unique pattern of soil hydraulic and thermal dynamics in the active layer can be one important contributor for the rapid permafrost warming at the study site. These findings suggest that, given the increase in air temperature and precipitation, soil hydraulic properties, particularly soil architecture in those thick active layers must be properly taken into account in permafrost models.

Abstract. Gravel (particle size ≥ 2 mm) is common in soil profiles of the Qinghai-Tibetan Plateau (QTP). It has different thermal and hydrological properties than other fine mineral soils (particle size < 2 mm), which may have significant impacts on the thermal and hydrological processes of soil. However, few models have considered gravel. In this study, we implemented the thermal and hydraulic properties of gravel into the Dynamic Organic Soil-Terrestrial Ecosystem Model to develop new schemes to simulate the dynamics of permafrost on the QTP. Results showed that: (1) the widely used Farouki thermal scheme always simulated higher thermal conductivity of frozen soils than unfrozen soils with the same soil water content; therefore it tends to overestimate permafrost thickness strongly; (2) there exists a soil moisture threshold, below which the new set of schemes with gravel simulated smaller thermal conductivity of frozen soils than unfrozen soils; (3) soil with gravel has higher hydraulic conductivity and poorer water retention capability; and simulations with gravel were usually drier than those without gravel; and (4) the new schemes simulated faster upward degradation than downward degradation; and the simulated permafrost thicknesses were sensitive to the fraction of gravel, the gravel size, the thickness of soil with gravel, and the subsurface drainage. To reduce the uncertainties in the projection of permafrost degradation on the QTP, more effort should be made to: (1) developing robust relationships between soil thermal and hydraulic properties and gravel characteristics based on laboratory work; and (2) compiling spatial datasets of the vertical distribution of gravel content based on measurements during drilling or the digging of soil pits.

Abstract. Climate change projections still suffer from a limited representation of the permafrost–carbon feedback. Predicting the response of permafrost temperature to climate change requires accurate simulations of Arctic snow and soil properties. This study assesses the capacity of the coupled land surface and snow models ISBA-Crocus and ISBA-ES to simulate snow and soil properties at Bylot Island, a high Arctic site. Field measurements complemented with ERA-Interim reanalyses were used to drive the models and to evaluate simulation outputs. Snow height, density, temperature, thermal conductivity and thermal insulance are examined to determine the critical variables involved in the soil and snow thermal regime. Simulated soil properties are compared to measurements of thermal conductivity, temperature and water content. The simulated snow density profiles are unrealistic, which is most likely caused by the lack of representation in snow models of the upward water vapor fluxes generated by the strong temperature gradients within the snowpack. The resulting vertical profiles of thermal conductivity are inverted compared to observations, with high simulated values at the bottom of the snowpack. Still, ISBA-Crocus manages to successfully simulate the soil temperature in winter. Results are satisfactory in summer, but the temperature of the top soil could be better reproduced by adequately representing surface organic layers, i.e., mosses and litter, and in particular their water retention capacity. Transition periods (soil freezing and thawing) are the least well reproduced because the high basal snow thermal conductivity induces an excessively rapid heat transfer between the soil and the snow in simulations. Hence, global climate models should carefully consider Arctic snow thermal properties, and especially the thermal conductivity of the basal snow layer, to perform accurate predictions of the permafrost evolution under climate change.

Abstract. Thermal modeling is a powerful tool to infer the temperature regime of the ground in permafrost areas. We present a transient permafrost model, CryoGrid 2, that calculates ground temperatures according to conductive heat transfer in the soil and in the snow pack. CryoGrid 2 is forced by operational air temperature and snow depth products for potential permafrost areas in Southern Norway for the period 1958 to 2009 at 1 km spatial resolution. In total, an area of about 80 000 km2 is covered. The model results are validated against borehole temperatures, permafrost probability maps from "Bottom Temperature of Snow" measurements and inventories of landforms indicative of permafrost occurrence. The validation demonstrates that CryoGrid 2 can reproduce the observed lower permafrost limit to within 100 m at all validation sites, while the agreement between simulated and measured borehole temperatures is within 1 K for most sites. The number of grid cells with simulated permafrost does not change significantly between the 1960s the 1990s. In the 2000s, a significant reduction of about 40% of the area with average 2 m ground temperatures below 0 °C is found which mostly corresponds to degrading permafrost with still negative temperatures in deeper ground layers. The thermal conductivity of the snow is the largest source of uncertainty in CryoGrid 2 strongly affecting the simulated permafrost area. Finally, the prospects of employing CryoGrid 2 for an operational soil temperature product for Norway are discussed.

Abstract. Thermal modeling is a powerful tool to infer the temperature regime of the ground in permafrost areas. We present a transient permafrost model, CryoGrid 2, that calculates ground temperatures according to conductive heat transfer in the soil and in the snowpack. CryoGrid 2 is forced by operational air temperature and snow-depth products for potential permafrost areas in Southern Norway for the period 1958 to 2009 at 1 km2 spatial resolution. In total, an area of about 80 000 km2 is covered. The model results are validated against borehole temperatures, permafrost probability maps from "bottom temperature of snow" measurements and inventories of landforms indicative of permafrost occurrence. The validation demonstrates that CryoGrid 2 can reproduce the observed lower permafrost limit to within 100 m at all validation sites, while the agreement between simulated and measured borehole temperatures is within 1 K for most sites. The number of grid cells with simulated permafrost does not change significantly between the 1960s and 1990s. In the 2000s, a significant reduction of about 40% of the area with average 2 m ground temperatures below 0 °C is found, which mostly corresponds to degrading permafrost with still negative temperatures in deeper ground layers. The thermal conductivity of the snow is the largest source of uncertainty in CryoGrid 2, strongly affecting the simulated permafrost area. Finally, the prospects of employing CryoGrid 2 as an operational soil-temperature product for Norway are discussed.

Thermal properties are important for featuring the water-heat transfer capacity of soil. They are also key to many processes in earth sciences, such as the land surface processes and ecological and geoenvironmental dynamics and their changes in permafrost regions. With loose and porous structures, the organic matter layer in soil strata substantially influences soil thermal conductivity. So far, thermal conductivity of mineral soils has been explored extensively and in depth, but there are only limited studies on that of organic soils. In this study, influences of soil temperature, soil moisture saturation (SMS), and soil organic matter (SOM) content on soil thermal conductivity were analyzed on the basis of laboratory experiments on the silt-organic soil mixtures of varied mixing ratios. Results show that soil thermal conductivity declines slowly with the lowering temperatures from 10 to 0°C; however, it increases and finally stabilizes when temperature further lowers from 0 to -10°C. It is important to note that thermal conductivity peaks in the temperature range of -2~0°C (silty and organic-poor soil) and -5~0°C (organic-rich soil), possibly due to phase changes of ice/water in warm permafrost. Under both thawed and frozen states, soil thermal conductivity is positively related with SMS. However, with rising SOM content, the growth rate of soil thermal conductivity with SMS slows gradually. Given the same SMS, soil thermal conductivity declines exponentially with increasing SOM content. Based on the experimental and theoretical analyses, a new empirical computational formula of soil thermal conductivity is established by taking into account of the SOM content, SMS, and soil temperature. The results may help better parameterize in simulating and predicting land surface processes and for optimizing frozen soil engineering designs and provide theoretical bases for exploring the dynamic mechanisms of environmental changes in cold regions under a changing climate.

Soil thermal conductivity is a dominant parameter of an unsteady heat-transfer process, which further influences the stability and sustainability of engineering applications in permafrost regions. In this work, a laboratory test for massive specimens is performed to reveal the distribution characteristics and the parameter-influencing mechanisms of soil thermal conductivity along the Qinghai–Tibet Engineering Corridor (QTEC). Based on the measurement data of 638 unfrozen and 860 frozen soil specimens, binary fitting, radial basis function (RBF) neural network and ternary fitting (for frozen soils) prediction models of soil thermal conductivity have been developed and compared. The results demonstrate that, (1) particle size and intrinsic heat-conducting capacity of the soil skeleton have a significant influence on the soil thermal conductivity, and the typical specimens in the QTEC can be classified as three clusters according to their thermal conductivity probability distribution and water-holding capacity; (2) dry density as well as water content sometimes does not have a strong positive correlation with thermal conductivity of natural soil samples, especially for multiple soil types and complex compositions; (3) both the RBF neural network method and ternary fitting method have favorable prediction accuracy and a wide application range. The maximum determination coefficient (R2) and quantitative proportion of relative error within ±10% ( P ± 10 % ) of each prediction model reaches up to 0.82, 0.88, 81.4% and 74.5%, respectively. Furthermore, because the ternary fitting method can only be used for frozen soils, the RBF neural network method is considered the optimal approach among all three prediction methods. This study can contribute to the construction and maintenance of engineering applications in permafrost regions.

Data representing the seasonal variation of thermal conductivity of the ground at depths within the seasonally active freezing/thawing zone are presented for a number of different soil conditions at four sites across Canada. An inexpensive probe apparatus suitable for routine field measurements is described.In all the cases examined, significant seasonal variations were confined to the first few decimetres. In addition to distinct seasonal differences associated with phase change, quite large changes occurred during the period when the soil was thawed in those cases where seasonal drying was possible. Below the seasonally active zone, thawed soil conductivities did not differ greatly among the three nonpermafrost sites in spite of soil composition ranging from marine clay to sandy silt. The data suggest that, even within a given soil layer, quite significant differences in thermal conductivity may be encountered in engineering structures such as embankments, presumably because of differences in drainage conditions. Key words: thermal conductivity, field measurements, phase relationships, drying, permafrost, clay, silt, peat.

ABSTRACTBasal depth hoar that forms in Arctic snowpacks often has a low thermal conductivity, strongly contributing to the snowpack thermal insulance and impacting the permafrost thermal regime. At Ward Hunt Island (Canadian high Arctic, 83°05′N, 74°07′W) almost no depth hoar was observed in spring 2016 despite favorable thermal conditions. We hypothesize that depth hoar formation was impeded by the combination of two factors (1) strong winds in fall that formed hard dense wind slabs where water vapor transport was slow and (2) low soil moisture that led to rapid ground cooling with no zero-curtain period, which reduced soil temperature and the temperature gradient in the snowpack. Comparisons with detailed data from the subsequent winter at Ward Hunt and from Bylot Island (73°09′N, 80°00′W) and with data from Barrow and Alert indicate that both high wind speeds after snow onset and low soil moisture are necessary to prevent Arctic depth hoar formation. The role of convection to form depth hoar is discussed. A simple preliminary strategy to parameterize depth hoar thermal conductivity in snow schemes is proposed based on wind speed and soil moisture. Finally, warming-induced vegetation growth and soil moisture increase should reduce depth hoar thermal conductivity, potentially affecting permafrost temperature.

Le pergélisol est une composante majeure du système climatique terrestre. Avec le réchauffement du climat, la dégel du pergélisol profite à l'activité biochimique qui décompose davantage de matière organique dans les sols arctiques et la rejette dans l'atmosphère sous forme de gaz à effet de serre (CO2, CH4). Ce phénomène pourrait constituer une rétroaction climatique positive majeure. Prédire ces effets nécessite d'étudier l'évolution du régime thermique du pergélisol ainsi que des facteurs qui l'influencent. Le manteau neigeux, de par son pouvoir isolant, contrôle les échanges de chaleur entre le sol et l'atmosphère une grande partie de l'année. Le flux de chaleur à travers la neige dépend de la hauteur du manteau neigeux et de la conductivité thermique des couches de neige qui le constituent. Ces deux variables sont elles-même très dépendantes des conditions climatiques et de la présence de végétation. Nous réalisons ici le suivi des propriétés de la neige et du sol d'un site haut arctique de toundra herbacée (Île Bylot, 73N), et d'un site bas arctique à la frontière de la toundra arbustive et forestière (Umiujaq, 56N). Nous utilisons les données issues de stations de mesure automatiques complétées par des mesures manuelles. Une attention particulière est portée sur la conductivité thermique de la neige, car peu de données sont disponibles pour les régions arctiques. Le modèle numérique couplé ISBA-Crocus est ensuite utilisé pour simuler les propriétés de la neige et du sol des deux sites étudiés. Les résultats sont comparés aux mesures de terrain afin d'évaluer la capacité du modèle à simuler le régime thermique des sols arctiques.Nous avons pu caractériser les interactions atmosphère-neige-végétation qui façonnent la structure des manteaux neigeux arctiques. Le vent et la redistribution de neige qu'il induit sont des paramètres fondamentaux qui déterminent la hauteur et la conductivité thermique de la neige. Un couvert végétal haut et dense (arbustes, arbres) piège la neige soufflée et l'abrite du tassement éolien. De plus, la structure ligneuse des massifs arbustifs soutient la masse de neige et empêche son tassement. Cet abri procure à la neige une capacité d'isolation élevée qui retarde le gel du sol dès les premières accumulations. Le refroidissement atmosphérique se poursuivant, le manteau neigeux peu épais est soumis à un gradient thermique élevé qui provoque d'importants transferts de vapeur d'eau depuis le sol et les couches de neige basales, vers les couches supérieures et l'atmosphère. La croissance de givre de profondeur qui s'opère, favorisée à la fois par le gradient thermique élevé et la faible densité de la neige, aboutit à la formation de couches très isolantes en contact avec la surface du sol. Tant que le sol demeure relativement chaud, la croissance de givre de profondeur perdure. Finalement, des épisodes de fonte peuvent avoir lieu en automne durant la mise en place du manteau neigeux dans les régions arctiques. Le regel de la neige peut rapidement annuler ou même temporairement inverser l'effet isolant des interactions neige-végétation. Une surface de neige gelée ne subit pas l'effet du vent et empêche sa redistribution. La formation de croûtes de regel à forte conductivité thermique accélère le refroidissement du sol. Le manteau neigeux affecté par la fonte au début de l'hiver a donc une capacité d'isolation diminuée qui pourrait entraver le réchauffement des sols arctiques. Nos résultats de simulation montrent que ces différents effets ne sont pas correctement représentés dans les modèles de neige. Les erreurs dans les conductivités thermiques de la neige simulées sont particulièrement problématiques puisqu'elles interviennent lors de la période de gel du sol. Étant donné l'étendue des régions affectées par le pergélisol, ces erreurs sur la modélisation de la neige arctique pourraient significativement affecter les simulations climatiques et les projections de la hausse des températures globales
Permafrost is a major component of the Earth climatic system. Global warming provokes the degradation of permafrost which favors biogeochemical activity in Arctic soils. The decomposition of organic matter increases and results in the release of high amounts of greenhouse gases (CO2 and CH4) to the atmosphere. By amplifying the greenhouse effect induced by human activities, this phenomenon may constitute one of the strongest positive feedbacks on global warming. Predicting these effects requires to study the evolution of the permafrost thermal regime and the factors governing it. The snowpack, because of its insulating effect, modulates the heat fluxes between permafrost and atmosphere most of the year. The snow insulating capacity depends on snow height and thermal conductivity. These two variables are highly dependent on climatic conditions and on the presence of vegetation. Here we monitor the snow and soil physical properties at a high Arctic site typical of herbaceous tundra (Bylot Island, 73°N), and at a low Arctic site situated at the limit between shrub and forest tundra (Umiujaq, 56°N). We use data from automatic measurement stations and manual measurements. A special attention is given to the snow thermal conductivity because very few data are available for Arctic regions. Results are interpreted in relation to vegetation type and atmospheric conditions. The numerical coupled model ISBA-Crocus is then used to simulate snow and soil properties at our sites. Results are compared to field data in order to evaluate the model capacity to accurately simulate the permafrost thermal regime.We managed to describe atmosphere-snow-vegetation interactions that shape the structure of Arctic snowpacks. Wind and the snow redistribution it induces are fundamental parameters governing snow height and thermal conductivity. A high vegetation cover (i.e. shrubs and forest) traps blowing snow and shields it from wind compaction. Vegetation growth thus favors the formation of an insulating snowpack which slows down or even prevents soil freezing. Furthermore, the shrubs woody structure supports the snow mass and prevents the resulting compaction of bottom snow layers. Thus sheltered, snow in shrubs develops a high insulating capacity which delays soil freezing. Continued atmospheric cooling increases the thermal gradient in the snow, maintaining large water vapor transfers from the soil and the snow basal layers to upper layers and atmosphere. The growth of depth hoar, enhanced by the large thermal gradient and the low snow density, results in the formation of highly insulating snow layers thus constituting a positive feedback loop between soil temperature and snow insulation. As long as the soil stays relatively warm, depth hoar growth persists. Finally, if warm spells occur in autumn, they can trigger the partial melting of the early snowpack which can cancel or temporarily reverse the insulating effect of snow-vegetation interactions. A frozen snow surface prevents snow drifting and its redistribution. The presence of highly conductive refrozen layers facilitates soil cooling and reduces the thermal gradient. An early snowpack affected by melting is thus less insulative which could hamper Arctic soil warming. Simulation results show that these different effects are not correctly represented in snow models. Errors in the estimated snow thermal conductivities are particularly problematic as they highly affect the simulation of soil freezing. Given the area of permafrost-affected regions, these errors on Arctic snow modelling could significantly impact climate simulations and the global warming projections

Hydrological processes exert strong control over biological and climatic processes in every ecosystem. They are particularly important in the boreal zone, where the average annual temperatures of the air and soil are relatively near the phase-change temperature of water (Chapter 4). Boreal hydrology is strongly controlled by processes related to freezing and thawing, particularly the presence or absence of permafrost. Flow in watersheds underlain by extensive permafrost is limited to the near-surface active layer and to small springs that connect the surface with the subpermafrost groundwater. Ice-rich permafrost, near the soil surface, impedes infiltration, resulting in soils that vary in moisture content from wet to saturated. Interior Alaska has a continental climate with relatively low precipitation (Chapter 4). Soils are typically aeolian or alluvial (Chapter 3). Consequently, in the absence of permafrost, infiltration is relatively high, yielding dry surface soils. In this way, discontinuous permafrost distribution magnifies the differences in soil moisture that might normally occur along topographic gradients. Hydrological processes in the boreal forest are unique due to highly organic soils with a porous organic mat on the surface, short thaw season, and warm summer and cold winter temperatures. The surface organic layer tends to be much thicker on north-facing slopes and in valley bottoms than on south-facing slopes and ridges, reflecting primarily the distribution of permafrost. Soils are cooler and wetter above permafrost, which retards decomposition, resulting in organic matter accumulation (Chapter 15). The markedly different material properties of the soil layers also influence hydrology. The highly porous near-surface soils allow rapid infiltration and, on hillsides, downslope drainage. The organic layer also has a relatively low thermal conductivity, resulting in slow thaw below thick organic layers. The thick organic layer limits the depth of thaw each summer to about 50–100 cm above permafrost (i.e., the active layer). As the active layer thaws, the hydraulic properties change. For example, the moisture-holding capacity increases, and additional subsurface layers become available for lateral flow. The mosaic of Alaskan vegetation depends not only on disturbance history (Chapter 7) but also on hydrology (Chapter 6).

Pipelines carrying warm contents through permafrost areas present a key challenge to the development of petroleum reserves in the Arctic region. Evaluations have to be made during the pipeline design to avoid signification thaw settlement and large-scale permafrost damage. Current practice in estimating pipeline thaw settlement is primarily based on the assumption of a complete consolidation of the thaw-unstable layer. This often results in an over-conservative estimate and costly over-design. An integrated three-dimensional (3D) finite element (FE) model, developed recently as part of J P Kenny’s in-house simulator tools, has been used successfully to simulate the thermal and mechanical interaction between the pipeline and surrounding permafrost and gives a better estimate of the heat transfer, thaw settlement and pipeline deformation over the pipeline’s service life. Given the lack of experimental data and the complex nature of this problem, the sensitivity study on the major variables, presented in this paper and based on the 3D FE model, offers an in-depth insight into the problem and provides better-based guidelines for proper design of pipelines in permafrost. A number of thermal parameters are studied in this paper and compared based on their impact on the final deformation of the pipeline, so as to identify key parameters in the pipeline-permafrost thaw settlement processes. These parameters include the hydrocarbon content temperature and the convection coefficient, inside the pipe; the thickness and thermal conductivity of the insulation layer, on the pipe exterior; and the initial temperature, soil void ratio and the surface condition of the permafrost/soil, as boundary conditions. Results of current study improve the understanding of the pipeline-permafrost interaction from the heat transfer perspective and provide better guidance to the pipeline design in the permafrost environment.

Abstract Permafrost thawing caused by the hot crude pipeline is a major threat to the safe operation of buried pipelines in permafrost zone. In this paper, the process of thawing and consolidation of frozen soil is considered, and a three-dimensional (3D) finite element model of buried pipelines in permafrost zone is established using ABAQUS. The calculation of thaw settlement displacement of frozen soil based on moisture-heat-stress coupled was carried out, and the deformation and stress of buried pipelines were analyzed. The effects of ground temperature, oil temperature, thermal conductivity of insulation material and soil distribution along the pipeline on the vertical displacement and longitudinal stress of buried pipelines in frozen soil were studied. Research results show that in thaw-unstable soil, the vertical displacement and stress of the pipeline increase significantly with the increase of the average ground temperature, and change on ground temperature amplitude has a little effect on the vertical displacement and longitudinal stress of the pipeline in thaw settlement zone. It is 1/3 of the vertical displacement of the pipeline without a heat insulating layer. When the thermal conductivity of the insulation material is less than 0.4 W/m °C, the vertical displacement of the pipeline in the thawing zone can be further reduced by reducing the thermal conductivity of the insulation material. When clay and sand appear alternately along the pipeline, the vertical displacement and longitudinal stress of the pipeline can be reduced by reducing the length of clay section. This study has certain reference value for optimizing the design parameters of buried pipelines in permafrost zone and reducing the impact of differential thaw settlement of frozen soils on the safe operation of pipelines.

Abstract Natural gas hydrate is widely distributed in the permafrost and marine deposits, and is regarded as an energy resource with great potential. The low-frequency electric heating assisted depressurization (LF-EHAD) has been proven to be an efficient method for exploiting hydrate sediments, which involves complex multi-physics processes, i.e. current conduction, multiphase flow, chemical reaction and heat transfer. The physical properties vary greatly in different hydrate sediments, which may profoundly affect the hydrate decomposition in the LF-EHAD process. In order to evaluate the influence of hydrate-bearing sediment properties on the gas production behavior and energy utilization efficiency of the LF-EHAD method, a geological model was first established based on the data of hydrate sediments in the Shenhu Area. Then, the influence of permeability, porosity, thermal conductivity, specific heat capacity, hydrate saturation and hydrate-bearing layer (HBL) thickness on gas production behavior is comprehensively analyzed by numerical simulation method. Finally, the energy efficiency ratio under different sediment properties is compared. Results indicate that higher gas production is obtained in the high-permeability hydrate sediments during depressurization. However, after the electric heating is implemented, the gas production first increases and then tends to be insensitive as the permeability decreases. With the increasing of porosity, the gas production during depressurization decreases due to the low effective permeability; while in the electric heating stage, this effect is reversed. High thermal conductivity is beneficial to enhance the heat conduction, thus promoting the hydrate decomposition. During depressurization, the gas production is enhanced with the increase of specific heat capacity. However, more heat is consumed to increase the reservoir temperature during electric heating, thereby reducing the gas production. High hydrate saturation is not conducive to depressurization because of the low effective permeability. After electric heating, the gas production increases significantly. High HBL thickness results in a higher gas production during depressurization, while in the electric heating stage, the gas production first increases and then remains unchanged with the increase of thickness, due to the limited heat supply. The comparison results of energy efficiency suggest that electric heating is more advantageous for hydrate sediments with low permeability, high porosity, high thermal conductivity, low specific heat capacity, high hydrate saturation and high HBL thickness. The findings in this work can provide a useful reference for evaluating the application of the LF-EHAD method in gas hydrate sediments.

Gas hydrates are anticipated to be a promising energy source with global distribution underneath the oceans and in permafrost regions. Hydrate dissociation referring to the phase transition of the solid gas hydrate into gas and liquid water, can occur due to altered environmental conditions like increase in temperature, decrease in pressure, or the injection of hydrate inhibitors. Numerical modeling work suggested that energy efficiency of the thermal simulation and inhibitor injection to dissociate gas hydrate in geological formations are quite low, and possibly negative. Depressurization is expected to be effective in many situations; however, the hydrate bearing sediments should be well enough confined by low permeability strata for the depressurization to be effective. In current study, numerical analyses were performed using the open-source reservoir simulator HYDRATERESSIM (HRS) to investigate methane gas hydrate dissociation behavior and production by depressurization, where the system pressure is lowered below the hydrate three-phase (gas-MH-water) stability pressure. Two-dimensional (2D) axisymmetric radial coordinates are used for the numerical simulations with local grid refinement. This paper presents the approach and results of numerical investigation conducted for a hypothetical hydrate-bearing reservoir in an offshore deep water site. A comprehensive parametric sensitivity study was undertaken to assess the effects of (1) various reservoir properties including sediment permeability (absolute permeability and relative permeability function), hydrate distribution, pore compressibility, porosity formation heterogeneity and thermal conductivity; and (2) operational condition on hydrate dissociation and production. Gas productivity and advance of hydrate dissociation front with time after the commencement of depressurization are discussed. This numerical study highlights the need and importance to obtain an accurate evaluation of reservoir properties from different data source, including seismic data, logging while drilling (LWD) data, core samples, field and laboratory testing data and drilling information, on gas hydrate production prediction. It can provide useful insights on the planning and design of a field production well test.