This paper presents the results of numerical modeling of the permafrost thermal regime and thermobaric conditions of methane hydrates in the north of Western Siberia over the past 70 thousand years. The area of hydrate formation was determined and the rate of accumulation of hydrates was estimated in connection with the migration of fluid from the underlying gas-saturated layers under the conditions of cover glaciation. The estimates obtained for the change in hydrate saturation as a result of fluid migration during the 10 thousand-year glaciation period, depending on the permeability of the soil, are from 6 to 40% in the upper 350 m. Based on quantitative characteristics of the equilibrium and metastable states of methane hydrates, the conditions for the preservation of relict methane hydrates in permafrost under the paleoclimatic scenario were determined, taking into account periods of ice cover and transgression. It is shown that due to the effect of self-preservation at temperatures below –4°C, it is possible to preserve relict methane hydrates in the upper 200 m of soil under non-equilibrium conditions. The effect of lowering the temperature while the hydrates dissociate prevents the complete decomposition of the deposit and leads to an increase in the thickness of the frozen soil.