CC BY-NC
24Развитие взрывных процессов в криолитозоне в связи с формированием Ямальского кратера
Хименков Александр Николаевич
кандидат геолого-минералогических наук ведущий научный сотрудник, Институт геоэкологии РАН 101000, Россия, г. Москва, Уланский ппроезд, 13, стр. 2
Станиловская Юлия Викторовна
младший научный сотрудник, Институт геоэкологии им. Е.М Сергеева Российской академии наук
101000, Россия, г. Москва, Уланский переулок, 13
Сергеев Дмитрий Олегович
кандидат геолого-минералогических наук заведующий, Институт геоэкологии Российской академии наук 101000, Россия, г. Москва, Уланский проезд, 13, стр. 2
И [email protected] Власов Александр Николаевич
доктор технических наук Директор, Институт прикладной механики Российской академии наук 125 040, Россия, г. Москва, ул. Ленинградский Проспект, 7
И [email protected] Волков-Богородский Дмитрий Борисович
кандидат физико-математических наук старший научный сотрудник, Институт прикладной механики Российской академии наук
125040, Россия, г. Москва, Ленинградский проспект, 7
И [email protected] Мерзляков Владимир Павлович
кандидат технических наук ведущий научный сотрудник, Институт геоэкологии Российской академии наук 101000, Россия, г. Москва, Уланский проезд, 13, стр. 2
Типенко Геннадий Сергеевич
кандидат физико-математических наук ведущий научный сотрудник, Институт геоэкологии Российской академии наук 101000, Россия, г. Москва, Уланский проезд, 13, стр. 2
Статья из рубрики "Многолетнемерзлые породы и подземные льды Арктики"
Аннотация. Аннотация. Предметом исследования являются естественные взрывные процессы в криолитозоны. Объектом исследования является Ямальский кратер. Авторами подробно рассматриваются группы естественных взрывных процессов на территории распространения многолетнемерзлых пород. Предложен новый подход в теории формирования воронок газового выброса. Он базируется на механизме миграции газовых флюидов из зоны диссоциации газогидратов в вышележащую толщу многолетнемёрзлых пород. Зона диссоциации формируется на участках локального прогрева за счёт поступления тепла сверху. Показана стадийность процессов подготовки взрывного процесса сформировавшего Ямальский кратер. Проведён расчёт давления выброса мёрзлой кровли при формировании Ямальского кратера. Выделены, и охарактеризованы 4 фазы подготовки естественных взрывных процессов. Рассмотрены различные сценарии подготовки естественных взрывных процессов в криолитозоне. Создание феноменологической модели процессов тепломассопереноса и напряженно-деформированного состояния в мёрзлых породах содержащих газовые гидраты при изменении температуры и давления. Экспериментальное и теоретическое исследование процессов разложения газогидратов Разработка подходов к оценке обеспечения устойчивости природно-технических систем Арктики, включающих мёрзлые породы содержащие газогидраты в условиях меняющегося климата и растущей техногенной нагрузки. Основными выводами проведённого исследования являются:взрывы гидролакколитов и образование воронок газового выброса относятся к одной группе процессов - физическим взрывам естественного происхождения;подготовка взрывных процессов в многолетнемёрзлвх породах, проходит через несколько стадий;при подготовке взрывов естественного происхождения в многолетнемёрзлых породах, большую роль играет миграция газовых флюидов
Ключевые слова: пластические деформации, стадийность, Ямальский кратер, воронка газового выброса, флюиды, диссоциация, газовые гидраты, криолитозона, физический взрыв, зона транзита
DOI: 10.7256/2453-8922.2017.4.25094
Дата направления в редакцию: 25-12-2017
Дата рецензирования: 25-12-2017
Работа выполнена при поддержке гранта РФФИ № 17-05-00294.
Введение
В области распространения многолетнемёрзлых пород (ММП) есть группа процессов известных давно, но не привлекавших к себе особого внимания, поэтому слабо изученных. К ним относятся взрывы гидролакколитов и наледных бугров. В последние годы к ним прибавились сходные по внешним проявлениям образования связанные с газовыми выбросами.
Прежде всего, следует определить, можно ли объединить отмеченные процессы в одну группу. Мы считаем, что их все можно отнести к одной группе физических взрывов. Рассмотрим насколько понятие взрыв применимо к процессам, наблюдаемым в
криолитозоне. Согласно ГОСТ Р 22.0.08-96, взрыв — это «процесс выделения энергии за короткий промежуток времени, связанный с мгновенным физико-химическим изменением состояния вещества, приводящим к возникновению скачка давления или ударной волны, сопровождающейся образованием сжатых газов или паров, способных производить
работу»
Взрывные процессы, наблюдаемые в криолитозоне, относятся к физическим взрывам, связанным с возникновением внутреннего давления в полости, заполненной жидкостью или газом, превышающим предельно допустимые значения прочностных характеристик вмещающей среды. Физические взрывы чаще всего связаны с неконтролируемым
высвобождением потенциальной энергии сжатых газов из замкнутых объемов К ним же следует отнести взрывы, формирующиеся при промерзании водного потока, или при всестороннем промерзании воды в замкнутом объеме, а так же взрывы сжатого газа содержащегося в грунтовом массиве вне зависимости от его генезиса. Основной причиной их образования является не химическая реакция, а физический процесс, обусловленный высвобождением внутренний энергии жидкости или газа за короткое время. Физические взрывы не сопровождаются химическими реакциями с выделением энергии, они являются следствием перехода запасенной потенциальной энергии в кинетическую энергию движения газов и жидкостей. После того когда замкнутая полость с содержащимся газообразным или жидким веществом при повышенном давлении разрушается, происходит выброс этого вещества с генерацией волн давления. При этом потенциальная энергия сжатого газа или жидкости переходит в энергию ударной волны.
Подготовка взрывов в естественных условиях сопровождается комплексом сменяющих друг друга процессов и образований. Предлагаемая статья посвящена выявлению общих закономерностей развития естественных взрывных процессов в криолитозоне.
Взрывы гидролакколитов
В геокриологии взрывные процессы в естественных условиях известны давно. Связаны они в большинстве случаев с взрывами гидролакколитов и наледных бугров. Образование и взрывы гидролакколитов обусловлено, как правило, промерзанием несквозных подозерных таликов и возникновением напорных подземных вод. В некоторых случаях такие взрывы наблюдались
и в сезонноталом слое. В 1930-х годах В.И. Андреев-2!, проводивший исследования на Ямале, сообщал, что, по рассказам ненцев, в Байдарацкой тундре зимою со стороны бугра раздавался оглушительный треск (взрыв), после чего здесь находили надземную наледь. В связи с этим автор предположил, что подо льдом находится полость, наполненная водой, а, быть может, отчасти и газом, которые во время глубокого
растрескивания бугров и выбрасываются на поверхность. О.И. Баженова-3-, рассматривая механизмы функционирования озерно-флювиальных систем Южного Забайкалья, приводит данные о взрывах гидролакколитов, сформировавшихся в озерных котловинах. Они имеют овальную форму, высоту 1.5-2.5 м, длину 50-100 м. Взрывы происходят не только зимой и весной, но и летом. Часто взрывы бугров сопровождаются сильным звуком, напоминающим орудийный выстрел, который слышен в радиусе 7-10 км. При этом возникают новые формы рельефа, представленные воронками — кратерами взрыва. Диаметр воронок изменяется от 1 до 15-25 м, глубина составляет 2-5 м. Во время взрыва из воронок выбрасывается большой объем льда и грунта. В днище пади Арангот 24 мая 1964 г. в 14 ч по местному времени А.Н. Скляревская наблюдала взрыв крупного гидролакколита (размеры 30х50 м при высоте 2.5 м, лед залегал с глубины
0.35 м), который был слышен на расстоянии до 6 км. Во время взрыва из центральной части бугра было выброшено большое количество песчано-дресвяного материала с глыбами льда, размер которых достигал 2х1.5х0.7 м. В результате взрыва образовалась
продолговатой формы воронка длиной 15 шириной 2-4 и глубиной около 2 м [4]. Н.С. Богомолов 20 июля 1964 г. наблюдал взрыв гидролакколита в долине р. Урейки (левый приток Акши). Ширина основания бугра 18-20 м, высота 0.8-1 м. Гидролакколит разорвался с грохотом, после которого последовало шумное извержение громадного количества воды в виде столба размером 3х2.5х1.7 м, с дебитом 12 м3/с,
продолжавшееся 15 мин [4]. Размеры гидролакколитов изменяются в больших пределах: диаметр основания от 20 до 250 м, а высота — от 2 до 70 м. Кровля бугров обычно состоит из переслаивающихся льдистых суглинков и супесей мощностью 2-10 м, под
которыми залегает ледяное ядро мощностью в несколько десятков метров (до 60 м) I33!. Растущие гидролакколиты подстилаются талыми водоносными горизонтами, или водяными линзами мощностью до 2 м, обладающими гидравлическим напором (рис. 1).
Взрывы наледных бугров
Для Забайкалья имеются наблюдения взрывов наледных бугров. Так, 28 марта 1927 г. произошел взрыв наледного бугра в долине р. Онон, общий объем льдогрунтовой массы, выброшенной взрывом, составил 508 м3. Самая большая глыба льда с прослоями песка,
гравия и галечника имела толщину 2 м, ширину 6-9 и длину 18 м I19!. Взрывом родникового бугра пучения вблизи д. Бырца 27 июля 1938 г. лед, песок и галечник были подняты на высоту 8-12 м и отброшены на расстояние 15 м. Фонтан воды высотой 2 м функционировал около 2 ч. В результате взрыва образовалась воронка диаметром 4.5 м
и глубиной 4 м ¿¿П.
Общее для взрывов гидролакколитов и наледных бугров
Следует отметить, что непосредственно взрыву гидролакколита или наледного бугра предшествует длительный процесс подготовки. Всестороннее промерзание талика в основание гидролакколита или сужение водного потока в основании наледного бугра приволит к возрастанию гидравлического напора. Усиление давления на мёрзлую кровлю приводит к её пластической деформации выражающейся в формировании бугра пучения. После того как предел прочности мёрзлой кровли будет преодолён произойдёт её разрыв, сопровождающийся извержениями фонтанов воды или жидкой грязи в несколько метров высотой и разбросом обломков на десятки метров вокруг. Судя по приведенным Дж. Маккеем данным, давление в водяной линзе (по высоте столба воды)
составляет около 2,5 МПа "З^. Несмотря на то, что взрывы гидролакколитов связаны с промерзанием водоносного горизонта, следует отметить, что в данном процессе всегда присутствует газовая составляющая. Это связано с тем, что при льдообразовании всегда выделяется газ, растворенный в воде.
Рис. 1. Сопоставление размеров структурных элементов растущего гидролакколита (пинго)
на Севере Аляски и Ямальского кратера
Воронки газового выброса
В 2014 г. на Ямале, в 30 км южнее Бованенковского нефтегазоконденсатного месторождения (НГКМ) была обнаружена круглая воронка (рис. 2), глубиной около 50. Её диаметр в верхней части около 30 м, с глубины 8м, 15м (рис.3). Непосредственно к воронке примыкает концентрический вал из выброшенной породы (рис.13 ). Никаких следов хозяйственной деятельности поблизости не обнаружено. Позже на севере Западной Сибири было обнаружено еще несколько воронок, близких по строению. Большинство исследователей (Богоявленский, М.О. Лейбман, А.И. Кизяков, В.В. Оленченко, М.И. Эпов, И.Н. Ельцов и др.), изучавших данные образования, сходятся во мнении, что их формирование обусловлено выбросом га зов. Источник газа различен,
диссоциация гидратосодержащих мерзлых пород [31, 14, 29, 18 и Др-1, или поступающих с
больших глубин из подмерзлотных горизонтов J5—61 Это дало основание В.И. Богоявленскому первому выделить данные образования как воронки газового выброса
Хотя это понятие закрепилось и широко используется, оно довольно неопределённо и позволяет широкое толкование. К выбросам, например, относят кратковременные или длительные поступления любых загрязняющих веществ или избыточного тепла из
техногенных источников ПШ. Морфология воронок и разброс пород вокруг, позволяют отнести процесс их формирования к физическому взрыву. В данной публикации не рассматриваются различные гипотезы происхождения газа, участвующего в образовании воронок. Это отдельная тема. Мы исходим из того, что подготовка физического взрыва в горной породе будь то вулкан, гейзер, газовый выброс в шахте или в многолетнемёрзлых грунтах, однотипна. В силу различных обстоятельств в породе формируется текучая субстанция (вода, газ, газонасыщенная лава), которая под воздействием под воздействием давления перераспределяется и накапливается в определённой зоне. В данной зоне повышается давление, вызывающее пластические деформации во вмещающей породе. После превышения предела прочности происходит физический взрыв (газовый выброс)с образованием специфических форм.
При формировании гидролакколитов и воронок газового выброса наблюдаются схожие процессы и явления. В обоих случаях в мёрзлых породах возникает повышенное давление, приводящее к росту многолетних бугров, и выбросов ледогрунтового
материала. Для гидролакколитов зафиксированы значения свыше 2 МПа^33^, для воронок газового выброса, наши расчёты (смотри ниже) и расчёты других авторов показали близкие к этому значения. Тем не менее, несмотря на сходство между этими типами криогенных образований существует значительные различия. Рассмотренные выше материалы показывают, что для гидролакколитов горизонтальные размеры связанных с ними таликов на много превосходят вертикальные. Бугры высотой 10-15м и
шириной от 30 до 300 м имеют площадь водосбора (талика) до 700 м в поперечнике I32!. В случае воронок газового выброса горизонтальные размеры полости в разы меньше вертикальных. Данные морфологические различия (рис. 1) связаны с различными движущими силами. В случае гидролакколитов это преимущественно вода, перемещающаяся в талой породе. В воронке газового выброса газ, движущийся в мерзлой льдистой породе.
Анализ имеющихся материалов
К настоящему времени наибольший объём информации о строении воронок газового выброса получен при изучении воронки, расположенной в 30 километрах к югу от Бованенковского ГКМ. Морфология данного образования, расширенная коническая верхняя часть и суженная цилиндрическая нижняя часть очень похожа на вулканический кратер. Как будет показано ниже, процессы формирующие воронку так же сходны с вулканическими процессами. Поэтому для рассматриваемого данного объекта мы используем название Ямальский кратер, предложенное новосибирскими исследователями
Рис. 2. Ямальский кратер. Фото В.В. Оленченко, 2014 г.
Форма Ямальского кратера - правильный эллиптический цилиндр, немного расширенный у поверхности (рис. 3). Глубина сезонно талого слоя по периметру равна 0,8 м. Материал, выброшенный из воронки, распределён равномерно во всех направлениях. Наибольшая концентрация обломков была обнаружена по окружности на расстоянии 140 м от центра воронки, а отдельные обломки разлетались на расстояние более 180 м. Разрез кратера сложен практически чистым льдом. В выброшенном материале
обнаружены камни - окатанные песчаники, сланцы I29!.
Рис. 3. Принципиальная схема строения Ямальского кратера 1 - вода озера на дне кратера
Особенностью геологического строения участка расположения Ямальского кратера является наличие горизонта газогидратов на глубине 60-80 м, выделенного по геофизическим данным ПШ.
Мы считаем, что основной причиной формирования Ямальского кратера является физический взрыв газа, образовавшегося при диссоциации газогидратов данного горизонта. Между выделением газа и взрывом проходит несколько этапов, со своим набором процессов и геологических образований. Для правильного понимания механизмов образования Ямальского кратера необходимо провести анализ строения слагающих его пород.
Формирование каверн и гротов в зоне разложения газогидратов
В нижней части кратера наблюдается серия каверн и гротов (рис. 4) размером от десятков сантиметров. По мере вытаивания, мелкие каверны сливаются, образуя большие гроты с гладкими стенками размером в несколько метров (рис. 6, 8). По примерным оценкам объём газогидрата метана сформировавшего Ямальский кратер
составляет около 12,5 м3 гидрата [15]. Такой незначительный объём вполне соответствует зоне небольших каверн и пустот, наблюдаемых в основании кратера. Аналогичные образования обнаружены в зонах выбросов газа (за счёт разложения газогидратов) в окрестностях острова Беннетта (Новосибирские острова) и в Охотском море. Здесь, в районе газового фонтана, обнаружены характерные формы: провальные воронки, ямы, поддонные ниши и пещеры, образующих систему сообщающихся гротов
Вероятно развитие системы замкнутых и соединённых между собой полостей, является типичным при диссоциации гидратонасыщенных пород.
Рис. 4. Каверны и гроты в зоне разложения газогидратов
(июль 2014). Фото В. А. Пушкарева
Течение флюидов газов сквозь мёрзлую толщу
Новым явлением, которое было обнаружено при исследовании воронок связанных с диссоциацией газогидратов, является формирование газовых флюидов, и их миграция сквозь льдистую мёрзлую толщу. На фотографии, сделанной в ноябре 2014 года в нижней части Ямальского кратера (рис. 5) хорошо видны потоки газовых флюидов. Нижние их части совпадают с зоной каверн и гротов и соответствуют зоне залегания
газогидратного слоя Верхние части газовых флюидов контактируют с залегающим выше, газонасыщенным ледогрунтом (рис. 6).
Рис. 5. Флюиды в нижней части Ямальского кратера (фото В.А. Пушкарёва)
Рис. 6. Связь газовых флюидов с ячеистыми образованиями на стенках Ямальского
кратера (фото В.А. Пушкарёва)
Формирование штока газонасыщенного ледогрунта
Характерным элементом поверхности вертикальных стенок Ямальского кратера (рис. 3) являются многочисленные круглые ячейки ( рис. 6, 7, 8, 9). Они распределены неравномерно, местами образую изометричные скопления, иногда вытягиваются в вертикально ориентиррованные цепочки. Размер ячеек составляет от 2-3 см до 30-40см. В нижней части воронки они контактируют с зоной газовых флюидов (рис. 6). В верхней части кратера содержание газовых включений в ледогрунте увеличивается, здесь поверхность кратера становится похожей на сжатые соты (рис. 9). Наблюдаемые на стенках воронки ячейки представляют собой отпечатки газовых пузырей содержащихся в мёрзлой породе.
Другим элементом строения вертикальных стенок Ямальского кратера, являются многочисленные пластические и разрывные деформации (рис. 7, 8, 9, 10). По всей протяжённости видны признаки течения льда и ледогрунта.
Ещё одним структурным элементом рассматриваемого образования, является
вертикальная ориентировка слоёв ледогрунта —22]. Слой с вертикальной ориентировкой формирует краевые части ледогрунтового тела образуя своеобразную оболочку в виде слоистого стакана. Механизм формирования вертикальной слоистости пока не выяснен. Это может быть результатом деформации изначально горизонтальной
слоистой ледогрунтовой толщи I22!, или формированием слоистой контактной зоны между движущейся под давлением льдистой породы и окружающего массива мёрзлых пород
[25]. В обоих случаях слой с вертикальной слоистостью, в стенках кратера сформирован в результате возникших избыточных давлений и деформаций. Данный слой соответствует границам кратера, он выделялся при исследованиях проводимых в 2014 году, к июлю 2015 года он вытаял и на его месте обнажились породы с субгоризонтальной
слоистостью [22].
Рассмотренные материалы позволяют предположить, что над зоной диссоциации газогидратов в льдистой мёрзлой толще, под давлением газовых флюидов сформировался газонасыщенный ледогрунтовый шток цилиндрической формы. В силу неоднородности возникающих давлений здесь происходило интенсивная перестройка первоначального строения, сопровождающаяся пластическими и разрывными
деформациями мёрзлого субстрата [25].
Рис. 9. Верхняя часть газонасышенного ледогрунта Газонасыщенный лёд над зоной разложения газогидратов июль 2014 (Фото В.В. Оленченко)
Пластические деформации кровли мёрзлых пород
Формирование газонасыщенного ледогрунтового штока вызвало выдавливание перекрывающих пород. Поверхностный слой многолетнемерзлых пород, имея низкие отрицательные температуры, является непроницаемым экраном для поступающих снизу газовых флюидов. Вследствие этого, на границе экрана формируется область повышенного давления. Это приводит к пластическим деформациям многолетнемёрзлых пород и формированию бугра пучения на поверхности. На рис 10. показано строение мёрзлых пород верхней части кратера, образовавшегося после взрыва разрушившего
виден фрагмент мёрзлой породы,
Рис. 10. Деформация кровли слоистых суглинков над областью диссоциации газовых гидратов, июль 2015.(Фото В.И. Богоявленского)
Стадия взрыва
Для лучшего понимания процессов, происходящих при взрывном образовании воронок газового выброса, воспользуемся предложенной М. А. Алидибировым Г1998"1 моделью вулканического взрыва как двухстадийного процесса. На стадии подготовки, в жерле вулкана формируется пористое газонасыщенное магматическое тело. На стадии развития взрыва, вызываемой резкой декомпрессией, происходит фрагментация (дробление) этого тела, высвобождение потенциальной энергии сжатых газов и выброс продуктов дробления (пирокластики) в атмосферу. Взрывная волна разгрузки, постепенно продвигаясь вглубь массива, слой за слоем разрушает пористый материал. Расширение газа в сторону области низкого давления приведет к выбросу продуктов дробления в
атмосферу. При быстром сбросе внешнего давления процесс газоотделения не играет существенной роли. Время релаксации напряжений оказываются значительно больше
чем время разрушения. При этом, само разрушение магмы носит хрупкий характер -Ш.
Аналогичные процессы наблюдаются и при формировании Ямальского кратера. В результате проникновения газовых флюидов, образующихся при диссоциации газогидратов, формируется газонасыщенное ледогрунтовое тело, перекрытое газонепроницаемой деформируемой кровлей. Давление в газовых включениях значительно превосходит атмосферное. Продвигаясь вглубь массива, взрывная волна слой за слоем будет выбрасывать ледогрунтовый материал. Причём чем ниже будет зона разрушения, тем энергичнее будет дробиться материал и дальше выбрасываться обломки. Данный процесс захватит только растеплённые участки мёрзлых пород, находящихся над областью разложения газогидратов. Окружающие многолетнемёрзлые низкотемпературные породы в рассматриваемые процессы включены не будут. Процесс послойного разрушения происходи в течении нескольких секунд. В результате в земле образуется круглая полость глубиной около 60 м. В нижней часть полости располагается зона каверн и гротов, соответствующая зоне диссоциации газогидратов. В верхней части Ямальского кратера формируется конусообразное расширение, соответствующее скалывающим напряжениям при первичном взрыве.
Стадийность формирования Ямальского кратера
Изучение структуры пород, слагающих Ямальский кратер, позволило выделить пос ле до в а те льно с ть с фо рмиро в а в ших е го проце с сов .
I стадия
Формирование озера над многолетнее мёрзлыми породами, содержащими пластовую залежь подземного льда. Под озером формируется чаша протаивания и зона повышенной, по сравнению с окружающими породами, температурой.
II с та дия
Зона повышения температуры в многолетнемёрзлых породах достигает слоя газогидратов. Начинается процесс диссоциации газогидратов с выделением газа и переохлаждённой воды. Газововодяные флюиды под воздействием давления начинают мигрировать в наименее прочную область, каковой является зона повышенной температуры под озером и насыщать первичный ледогрунтовый субстрат газом.
III с та дия
Разложение метастабильного газового гидрата вызывает рост давления, которое в значительно превышает пластовое. В результате происходит движение флюидов, направленное вверх и сопровождающееся пластическими и разрывными деформациями. В мёрзлой породе, представленной слоистой ледогрунтовой пластовой залежью за счёт насыщения флюидами, постепенно формируется шток, который состоит из слоистого газонасыщенного льда, с многочисленными ячейками различных размеров. Высота образовавшегося штока около 40 м, диаметр около 15 м. Формирующийся ледовогазовый шток выпучивает вышележащий слой. По мере обмеления водоема и начала промерзания талика на нижней границе талика за счёт миграции формируется инфильтрационно-сегрегационные льды. Переход из субаквального состояния в субаэральное понижает температуру верхнего слоя грунта, что увеличивает его прочность. При этом образуется мёрзлый газонепроницаемый экран, деформация
которого формирует бугор пучения.
IV стадия После того, как предел прочности кровли будет преодолён, происходит разрыв мёрзлой кровли и резкий сброс давления. Расширение газа в верхнем слое газонасыщенного ледогрунта приводит к взрыву. Механизм взрыва описан выше. В первую очередь будет разрушен сформировавшийся бугор пучения. Его фрагменты в виде концентрического вала расположены непосредственно вблизи кратера (рис. 2, 3). Последовательность смены стадий развития Ямальского кратера приведена на рис 11.
Рис. 11. Стадии развития Ямальского кратера (I, II, III, IV) -t25-
Обозначениям — покровный горизонт; 2 — слой льда между талыми и мёрзлыми породами; 3 — мерзлый газонасыщенный ледогрунт со следами пластических деформаций ; 4 —инфильтрационно-сегрегационный лёд;
5 — газонепроницаемая кровля ММП; 6 — зона повышения температуры в ММП под озером; 7 — ММП вне отепляющего воздействия озера; 8 — слой гидратосодержащих ММП; 9 — зона разуплотнения в слое гидратосодержащих ММП, примыкающая к кратеру; 10 — направление движения флюидов; 11 — газовые флюиды; 12 — гроты и каверны в нижней части кратера; 13 — озеро; 14 — кратер, сформировавшийся после выброса газонасыщенного ледогрунта; 15 — талик (оформление Д. Волкова)
Факторы, повлиявшие на формирование Ямальского кратера
Климатические изменения
Следует отметить данные по температурам многолетнемёрзлых пород, имеющиеся для Ямала показали недостаточность возмущающих воздействий климатических изменений для возможности разложения газовых гидратов на глубинах, соответствующих глубине воронки (около 60 м). Температуры пород в скважинах геокриологического стационара Марре - Сале ВСЕГИНГЕО показывают, что с 2001 по 20015 практически не изменились и
колеблются в диапазоне -5-6°С t12!. На глубинах от 20 до 60м наблюдается интервал безградиентного распределения температур (рис. 15). Вероятно, это связано с тем, что в засолённых грунтах широко распространённых на Ямале фазовые переходы происходят в области более низких температур. Незначительные колебания температурного поля будут сглаживаться фазовыми переходами. Приведённые данные показывают недостаточность возмущающих воздействий климатических изменений для возможности разложения газовых гидратов на глубинах, соответствующих глубине воронки (около 60 м) (рис. 12). Что позволяет усомниться в утверждениях некоторых авторов о ведущей роли резкого повышения летних температур в течении нескольких лет, предшествующих образованию
Ямальского кратера 117, 14]. Климатические изменения последних десятилетий могли значительно изменить температуру многолетнемёрзлых пород только в верхнем 10-20 метровом горизонте.
а
-18 -12 -в
г* с
-ь ч
\ Л *
Л } 1
VII - гоо - VII- 2002
X- 2008-гтъа* VIII- П1 1 гооэ т1п
Рис. 12. Температура пород в скважинах геокриологического стационара Марре-Сале ВСЕГИНГЕО 1121
Повышение температуры многолетнемёрзлых пород под воздействием
поверхностного водоёма
Поверхностные водоёмы являются наиболее сильным отепляющим фактором в криолитозоне. Поэтому необходимо рассмотреть материалы по формированию температурного поля в многолетнемёрзлых породах под поверхностными водоёмами.
Существуют лишь единичные измерения температурного поля в субаквальных озёрных осадках криолитозоны. Рассмотрим данные по исследованию озера Иллисарвик, расположенного в дельте р. Маккензи. В озере размером 300 х 600 м с максимальной глубиной 4,5 м. Среднегодовая температура мёрзлых пород в районе озера находится в диапазоне - -8 - -10°С. Мощность мерзлоты - 400 - 600м. В глубоководной части озера кровля мёрзлых пород залегает на глубине 20 - 30 м. Температура мёрзлых пород на глубине около 60 м ниже дна озера составляет около - 2°С , а на глубине 90 м около -3°С 130! (рис. 13).
Рис. 13. Температурный профиль, в субаквальных осадках центральной, глубоководной части озера Иллисарвик (глубина 4,5 м) i30!
Для севера Западной Сибири, где были обнаружены кратеры газового выброса, данных о температурах субаквальных осадков нет. Отепляющее влияние озера на многолетнемёрзлые породы можно оценить анализируя температуры в скважинах пробуренных на хасыреях - днищах осушившихся озёр (рис. 14). Температуры, замеренные на поверхности хасыреев, на глубине б0-80м находятся в диапазоне: -3 - -2,5°С (рис.14а ) , для III морской террасы и --1,8 - -1,б°С для поймы (рис. 14б)[1б]. Температуры многолетнемёрзлых пород вне зоны озёр на глубине б0-80м, находятся для III морской террасы и поймы в диапазоне -5,5 - -5°С
б
Рис. 14. Температура криогенной толщи под хасыреем на территории Бованенковского ГКМ
а) III морская терраса; б) пойма
Расчёты стационарного состояния температурного поля под дном озера радиусом 45 м, проведённые Г.С. Типенко, показали, что на глубине 60-80 м температура пород находится в диапазоне -2 - -2,5°С, что соответствует результатам натурных измерений.
Рис. 15
Рис. 15. Стационарное температурное поле под дном озера радиусом 45 м.
Важно выяснить, не только, степень воздействия озёр на температуру грунтов, но и морфологию формирующейся под озером зоны отепления. На трёхмерной диаграмме сейсмического зондирования видно (рис. 16), что под озёрами формируются зоны
пониженных скоростей в виде перевёрнутых конусов I34!. Пониженные скорости объясняются тем, что повышение температуры приводит к увеличению количества незамёрзшей воды, газовых пузырьков, а для засолённых грунтов появление свободной минерализованной воды (криопэгов). Интересно, что для небольших озёр иногда конусов не видно, но на срезах, проведённых через 100 м, выделяются округлые пятна. Это является свидетельством того, что даже небольшие озёра оказывают значительное отепляющее воздействие на многолетнемёрзлые породы.
Рис. 16. Трёхмерное распределение областей пониженных скоростей под озерами [34]
Рассмотренные материалы свидетельствуют, что в низкотемпературных мёрзлых породах под озёрами формируются зоны растепления, в которых температуры пород находятся в диапазоне -1 - -3°С. Мощность данных зон может превышать 100м. При повышении температуры в многолетнемёрзлых породах и льдах происходит уменьшение прочности, появление локальных дефектов, связанных с объёмными температурными напряжениями. Повышение температуры под озером изменяет структуру и свойства многолетнемёрзлых пород и льдов, формируя благоприятные условия для проникновения в них газовых флюидов.
Зона растепления, достигает глубины залегания гидратосодержащих мёрзлых пород , вызывая их диссоциацию и формируя область миграции газовых флюидов. При этом, диссоциация газогидратов и миграция флюидов происходит при отрицательной те мпе ра туре .
Исследования экспедициии МГУ под руководством В.З. Хилимонюк показали-24!, что на этой территории наиболее низкие температуры пород находятся в диапазоне от -1°С до -5°С. Наиболее низкие температуры (-4 - -5°С) формируются на водораздельных поверхностях III морской террасы в пределах ландшафтов кустарничково-моховой тундры. А под озШрами и полосами стока, вследствие отепляющего влияния водных покровов с открытым зеркалом воды наблюдаются наиболее высокие среднегодовые температуры 0 - -1°С. На территории развиты несквозные талики. Георадарное зондирование выявило несквозной талик мощностью 2-3 м, сформировавшийся под дном
современного термокарстового озера i24!. Таким образом, поверхностные условия, такие как наличие поверхностных водоёмов, подозёрных таликов благоприятствуют локальному прогреву ММП. Близкое залегание пластовых льдов на склонах, способствует развитию оползней и формированию в их тыловых частях озёр. Наличие пластовых залежей приводит к протаиванию и формированию небольших, но глубоких незамерзающих озёр.
Давление газов
Разложение газового гидрата вызывает рост давления, которое может значительно превышать не только гидростатическое, но и пластовое ^^ При температурах в диапазоне 268, 15 Т К (-5°С) -f 273, 15 Т К (0°С) равновесное давление в системе газ —
вода (лед) — гидрат находится в диапазоне 2,57 - 2,17 МПа [27]. В области локального прогрева ММП под озёрами температуры многолетнемёрзлых пород находятся в диапазоне -1 - -3°С. При этих температурах, в мёрзлых породах и льдах развивается система дефектов, пор, ослабленных зон, а поровое давление соответствует
гидростатическому. На глубине 60-80 м, где выделяется гидратосодержащий слой гидростатическое давление составляет около 0,7 МПа. Из приведённых данных видно, что на этой глубине, в зоне разложения газогидратов, создаются давления, в три раза превышающие гидростатическое. Это обусловит фильтрацию газовых флюидов вверх.
При разложении газогидратов соблюдается определенное равновесие между давлением и температурой. В случае неполного разложения при изменении термобарических условий происходит восстановление температуры или давления до равновесного
значения [26]. В опытах Р. Ханта, моделирующих условия разложения газогидратов в керне, показано, что после каждого сброса давления оно возвращалось к значениям,
соответствующим пластовому I23!. В случае Ямальского кратера давление, возникающее
при диссоциации газогидратов, будет поддерживаться на всем протяжении формирования газонасыщенного льда. В зоне диссоциации давление, вследствие оттока газовых флюидов будет падать, но оно постоянно будет восстанавливаться до тех пор, пока здесь сохраняются термобарические условия или существуют сами газогидраты.
Оценим минимальное давление необходимое для взрыва, обусловившего формирование Ямальского кратера.
Предполагается, что газ под давлением выталкивает тело, составленное из конической и цилиндрической частей (рис. 17). Необходимая величина давления зависит от глубины выброса.
Рис. 17. Схема соотношения конической и цилиндрической частей Ямальского кратера
Конические стенки верхней части кратера несут следы разрыва породы под углом к горизонтали 400- 600, в то время, как почти вертикальные стенки нижней части воронки
имеют вид гладкого мёрзлого грунта [6]. Поэтому представим приближенно верхнюю часть выброшенного тела в виде перевёрнутого усеченного конуса, нижний радиус которого равен , а верхний -^2 , высота равна Н .
Отрыв происходит под действием растягивающих напряжений С , распределенных по боковой поверхности конуса. Для того, чтобы «сорвать» верхнюю часть выброшенного тела давлением с глубины Н следует преодолеть вертикальную силу сцепления боковой поверхности с грунтом Р , силу веса конуса (3 и силу атмосферного давления,
распределённого по внешней поверхности площадью на глубине #равно
ж К
2 . Т. о. критическое давление
, где
. (1)
Сила
Р0 = РП7ГЩ
2 , где - предельно допустимое растягивающее напряжение,
Зададим величины, входящие в вышеприведенные формулы, равными:
. Д1 = 12.5т
р=1600—^
= 0.4МРа 8 =
Тогда
^2 — 18.5т Я = 8т F = 6Л9xl08N, G = 9.59xl07N, Р0 = 1.07х108Г^р - 8.53х108М. Критическое
= 1.74МРа
давление "Л Наши данные близки к результатам расчётов В.И.
Богоявленского и И.А. Гарагаша, определивших, что для разрушения мёрзлой кровли Ямальского кратера достаточно давления даже в 1,25 МПа
Проведём расчёт давлений, определяющих выброс пород на разных глубинах. При глубине Н'уН необходимо учесть вес цилиндрического тела Ст и силу среза по боковым вертикальным стенкам Р . Считая приближенно цилиндрическое тело круговым, примем его радиус — 12.5т / а высоту, равной Н — Н . В таком случае Р — -Р + Сг . Сила ^ — gp7гií Г^2тгД(/г — Нгде Т^ _ предельно допустимое касательное
напряжение на срез. Тогда критическое давление равно
(2)
где
(3)
т+ — 1МРа ь и — т™ р"=1.7бМРа
Задавая значение ' í j-j.vj.j- / величину 1%—г1 —Шт., получим ^10 Зависимость критического давления от глубины представлена в таблице 1 и рис 18. Таблица 1. Зависимость критического давления от глубины.
Глубина, 'V'1'1- 8 18 28 38 60
п р,,МРа Давление, Л Ь/ 1,74 3,49 5,25 7,01 10,88
Рис. 18. Зависимость критического давления от глубины
Как видно из проведённых расчёта, уже на глубине более 20м давления, создающиеся при разложении газогидратов, в диапазоне названных выше температур, не приведёт к выбросу кровли, тем более на глубине 60м. Рассмотренный выше процесс последовательного послойного разрушения газонасыщенных мерзлых пород позволяет преодолеть это несоответствие.
Сценарии возможного развития взрывных процессов в многолетнемёрзлых
грунтах
Рассмотренный пример формирования Ямальского кратера, представляет лишь один из в а риа нто в ра з в ития в о ро но к г а з ов ого в ыброс а .
Различие природных условий, обуславливают различные сценарии подготовки взрывных процессов в многолетнемёрзлых грунтах.
При незначительном и кратковременном повышении температуры, диссоциация газогидратов может быстро остановиться. В результате в многолетнемёрзлых породах будет сформирована зона повышенной пористости за счет небольших газовых полостей.
Если озеро будет достаточно глубоким, а сформировавшийся под ним талик мощным, то высвободившийся при диссоциации газогидратов газ, будет выделяться со дна водоёма.
Это наблюдается во многих озерах Ямала Выделение газа длится от нескольких месяцев до нескольких лет, но без взрыва, поскольку есть возможность свободного его выхода.
Если гидратосодержащие многолетнемёрзлые породы залегают на незначительной глубине (первые десятки метров), то давление газов, возникшее при разложении газогидратов, может вызвать взрыв, минуя стадию миграции флюидов. При этом, бугор пучения предшествующий взрыву, будет слабо выражен и даже отсутствовать.
Несмотря на разнообразие естественных взрывных процессов в криолитозоне, разнообразия сценариев и условий их подготовки можно выделить некоторые общие, присущие для всех черты. Прежде всего обязательно наличие флюида (fluidus — текучий ). В качестве флюида может выступать любое вещество, поведение которого при деформации может быть описано законами механики жидкостей. В геологических процессах флюидами могут выступать газы, водные растворы, нефть, илы, магма, и даже
"твёрдые" вещества В нашем случае мы рассматриваем в роли флюидов воду, газ и их сочетание.
Анализ материалов по строению гидролакколитов, наледных бугров пучения, воронок газового выброса позволил выработать общую схему развития взрывных процессов в криолитозоне (рис. 19) .
Рис. 19. Общая структура взрывных образований в ММП
1 - зона формирования флюидов; 2 - зона транзита флюидов; 3 - зона накопления
Стадийность в подготовке взрывных процессов
* *
• я
Л 1 1
» *
флюидов; 4 - зона пластических деформаций
Разнообразие условий реализации взрывных процессов выражается в многообразии сценариев их развития.
Первой фазой подготовки взрывного процесса является формирование области возникновения флюида. Для деятельного слоя, это непромёрзшая зона. Флюиды также формируются на участках многолетнемёрзлых пород подвергшихся протаиванию. При промерзании таких участков водонасыщенные талые породы оказываются в замкнутом состоянии при постоянно увеличивающемся давлении. Для массива гидратосодержащих мёрзлых пород это области повышения температуры, где возможна диссоциация газогидратов. Для области распространения криопэгов это зоны, где возникает дополнительное давление (возможно при понижении температуры и частичном промерзании криопэга), обуславливающее их миграцию и перераспределение в многолетнемёрзлых породах. Как, уже отмечалось выше, существует предположение, что взрывные процессы в криолитозоне могут происходить за счёт проникновения подмерзлотных газовых флюидов. Хотя таких выходов, в местах воронок газового выброса, не обнаружено, следует рассмотреть возможность и условия их развития. В начальной фазе на границе мёрзлых и талых пород должно сформироваться зона повышенного давления газа или повышение температуры.
Вторая фаза связана с перераспределением флюида в область меньшего давления. При этом формируется зона транзита. Для промерзающих талых пород поток грунтовых вод, направлен в сторону растущих бугров пучения(сезонных или многолетних). Для области распространения криопэгов зона транзита обеспечивает отжимание засолённых вод в области с меньшим давлением. Для гидратосодержащих пород, это участки высокотемпературных многолетнемёрзлых пород в качестве флюида выступает выделившийся при диссоциации газогидратов газ, возможно, какую-то роль играет переохлаждённая вода. Участки проникновения флюдов, из под слоя многолетнемёрзлых пород, приурочены к зонам повышенных тепловых потоков.
Третья фаза связана с формированием зоны накопления флюида. Чаще всего эта зона локализуется на небольшой глубине, под кровлей низкотемпературных или льдонасыщенных многолетнемёрзлых пород. Данный слой является экраном, препятствующим выход флюида на поверхность. Линзы воды прослеживаются в основании бугров пучения. Наиболее благоприятными зонами для накопления газа являются льдистые и низкотемпературные мёрзлые породы.
Четвертая фаза связана с возникновение пластических деформаций газо- и водонепроницаемой кровли, при повышении давления в области накопления флюида. Внешне области пластических деформаций проявляются в виде бугров пучения. После того как давление преодолеет прочность кровли происходит взрыв.
Выделенные фазы отражают общую последовательность событий, подготовки естественных взрывных процессов, но не пространственное расположение зон, соответствующих разным фазам. Например, при промерзании внутригрунтового водного потока 1, 2 и 3 фазы располагаются в пределах одной зоны. При промерзании талика 1 и 2 фазы расположены в одних границах. При быстром нарастании давления 4 фаза может быть слабовыраженной.
Возникает вопрос, как соотносятся наши взгляды о роли флюидов в естественных взрывных процессах в криолитозоне с представлением о инъекционном механизмом льдообразования? Согласно П. А. Шумскому, инъекционные льды формируются при
промерзании воды или водонасыщенного грунта внедрившегося под напором в горную породу [Шумский, 1955]. То есть, инъекционное льдообразование включает в себя только один механизм напорных процессов, реализующийся в том случае, когда давление в водоносном горизонте превышает предел прочности породы и происходит её гидроразрыв. Это лишь один процесс, завершающий серию процессов обеспечивающих напорную миграцию грунтовых вод. При этом газовая составляющая и её перераспределение, механизмом инъекционного льдообразования не рассматриваются.
Заключение
Первичной причиной формирования Ямальского кратера является диссоциация газогидратов, залегающих в толще многолетнемёрзлых пород.
Данный процесс обусловлен локальным прогревом пород под воздействием поверхностного водоёма. Под озером формируется область теплового воздействия, имеющая двучленное строение. Непосредственно под озером расположена талая зона с температурами выше 0°С, ниже расположена зона высокотемпературных мерзлых пород ^ - -1 - -3°С). После того, как зона высокотемпературных мерзлых пород достигает горизонта гидратосодержащих мёрзлых пород, запускается механизм диссоциации газогидратов.
Изучение Ямальского кратера позволило выявить стадийность его развития. Каждой стадии соответствует определённый набор процессов преобразования многолетнемёрзлых пород и определённые элементы строения кратера. В результате диссоциации газовых гидратов формируется зона повышенного давления, превышающая гидростатическое. Вследствие этого газовые флюиды проникают в вышележащие слои. Возникающие пластические деформации кровли приводят к формированию бугра пучения. В последующем, когда пластические деформации мёрзлой кровли, достигнут предела прочности, происходит её разрыв, выброс газонасыщенного ледогрунта, и формирования кратера.
Диссоциация газогидратов осуществляется в локальных зонах повышения температуры многолетнемёрзлых пород под водоёмами, и на других участках с благоприятными поверхностными условиями, а так же на территориях техногенных воздействий.
В зависимости от соотношения различных факторов сценарии развития воронок газового выброса и, соответственно, их морфология, могут различаться.
Анализ материалов по строению Ямальского кратера показал, что в его образовании большую роль играет миграция газовых флюидов сквозь мёрзлые породы. В настоящее время данное явление не изучено.
При всей близости процессов сопровождающих взрывы гидролакколитов, наледных бугров и формирование воронок газового выброса имеются значительные различия в причинах их возникновения и истории развития. Эти различия должны учитываться при разработке методики исследования данных образований.
Воронки газового выброса в зоне распространения многолетнемёрзлых пород представляют опасность при строительстве и эксплуатации инженерных сооружений. По мере освоения Арктики тепловое воздействие на многолетнемёрзлые породы будет возрастать, опасность взрывных процессов для инженерных сооружений так же увеличится. Тем не менее, данная группа процессов не только не учитывается при расчётах, выборе проектных решений и прогнозах взаимодействия инженерного
сооружения с многолетнемёрзлыми грунтами, но даже не включена в группу опасных геологических процессов.
Библиография
1. Алидибиров М. А. Механизм фрагментации сильновязкой магмы при вулканических взрывах (экспериментальное исследование) : автореферат диссертации на соискание ученой степени д-ра физ.-мат. наук: 04.00.22. М., 1998.-44 с.
2. Андреев В.И. Гидролакколиты (булгунняхи) в Западно-Сибирских тундрах // Известия Государственного географического общества.1936. Т. 68, вып. 2. С. 186210.
3. Баженова О.И. Современная динамика озерно-флювиальных систем Онон-Торейской высокой равнины (Южное Забайкалье) // Вестник Томского государственного университета. 2013. № 371. С. 171-177.
4. Богомолов Н.С., Скляревская А.Н. О взрывах гидролакколитов в южной части Читинской области // Наледи Сибири. М.: Наука, 1969. С. 127-130. Богоявленский В.И. Угроза катастрофических выбросов газа из криолитозоны Арктики. Воронки Ямала и Таймыра // Бурение и Нефть. № 9, 2 0 1 4. С. 12-17.
5. Богоявленский В.И. Выбросы газа и нефти на суше и акваториях Арктики и Мирового океана // БУРЕНИЕ И НЕФТЬ, ИЮНЬ, 2015, С. 4-10.
6. Богоявленский В.И., Гарагаш И.А. Обоснование процесса образования кратеров газового выброса в Арктике математическим моделированием // Арктика: экология и экономика. № 3 (19), 2015. С. 12-17.
7. Васильева З.А. Джафаров Д.С., Аметова Т.А. Косвенные техногенные признаки индикации газогидратов в криолитозоне // Криосфера Земли. Т. XV. № 1, 2011. С. 61-67.
8. Горная энциклопедия / Гл. ред.: Е. А. Козловский. М.: Советская энциклопедия. Т.5, 1991. 541 с.
9. ГОСТ Р 22.0.08-96. Безопасность в чрезвычайных ситуациях. Техногенные чрезвычайные ситуации. Взрывы. Термины и определения. М.: ИПК Издательство стандартов, 1996.
10. Гражданская защита. Энциклопедия. Т.1 Ф-И / Под общ. Ред. С.К. Шойгу; МЧС России. М.: Московская типография № 2, 2006. 551 с.
11. Девисилов В.А., Дроздова Т.И., Тимофеева С.С. Теория горения и взрыва. Практикум: учебное пособие. М.: Форум, 2012. 352 с.
12. Дубровин В.А., Крицук Л.Н., Полякова Е.И. Температура, состав и возраст отложений шельфа карского моря в районе геокриологического стационара Марре Сале // Криосфера Земли, 2015, т. XIX, № 4. С. 3-16
13. Дядин Ю.А., Гущин А.Л. Газовые гидраты // Соросовский Образовательный журнал. Вып. 3, 1998. С. 55-64.
14. Кизяков А.И., Сонюшкин А.В., Лейбман М.О., Зимин М.В., Хомутов А.В. Геоморфологические условия образования воронки газового выброса и динамика этой формы на центральном Ямале // Криосфера Земли. Т. XIX. № 2, 2015. С. 15-25.
15. Кожина Л.Ю., Микляева Е.С., Перлова Е.В., Синицкий А.И., Ткачева Е.В., Черкасов В.А. Опасные современные проявления криоактивности - основные результаты исследования ямальского кратера// Научный вестник Ямало-Ненецкого автономного округа, № 2/87. Салехард, 2015. С.19-28.
16. Криосфера нефтегазоконденсатных месторождений полуострова Ямал. Т. 2. Криосфера Бованенковского нефтегазоконденсатного месторождения / под общ.
ред. Ю.Б. Баду, Н.А. Гафарова, Е.Е. Подборного. М.: ООО Газпром экспо, 2013. 424 с.
17. Лейбман М.О., Плеханов А.В. Ямальская воронка газового выброса // Холод'ОК. № 2(12), 2014. С. 9-15.
18. Оленченко В.В., Синицкий А.И., Антонов Е.Ю, Ельцов И.Н., Кушнаренко О.Н., Плотников А.Е., Потапов В.В, Эпов М.И. Результаты геофизических исследований территории геологического новообразования «Ямальский кратер» // Криосфера Земли, 2015, т. XIX, № 4. С. 94-106.
19. Петров В.Г. Наледи на Амурско-Якутской магистрали. Л.: Изд-во АН СССР, 1930. 177 с.
20. Сизов О.В. Дистанционный анализ последствий поверхностных газопроявлений на севере Западной Сибири // Геоматика. Вып. 1, 2015. С. 53-68.
21. Стругов А.С. Взрыв гидролакколита (Читинская область) // Природа, 1955. №6. С. 117.
22. Стрелецкая И.Д., Лейбман М.О., Кизяков А.И., Облогов Г.Е., Васильев А.А., Хомутов А.В., Дворников Ю.А.Подземные льды и их роль в формировании воронки газового выброса на полуострове Ямал // Вестник Московского университета. Серия 5: География, 2017. Т. 1, № 2. С. 91-99.
23. Хант Дж. Геохимия и геология нефти и газа. М.: Мир, 1982. 704 с.
24. Хилимонюк В.З., Оспенников Е.Н., Булдович С.Н., Гунар А.Ю., Горшков Е.И. Геокриологические условия территории расположения ямальского кратера // Пятая конференция геокриологов России. Т. 2. М.: Университетская книга, 2016. С. 245255.
25. Хименков А.Н., Сергеев Д.О., Станиловская Ю.В., Власов А.Н., Волков -Богородский Д.Б. Газовые выбросы в криолитозоне, как новый вид геокриологических опасностей // Геориск. №3, 2017. С. 58-65.
26. Чувилин Е.М., Буханов Б.А., Гребёнкин С.И, Дорошин В.В., Иоспа А.В. Экспериментальное изучение прочности мёрзлых гидратосодержащих грунтов в условиях проявления самоконсервации порового гидрата // Пятая конференция геокриологов России. Т. 3. М.: Университетская книга, 2016. С. 180-186.
27. Шостак Н.А. Моделирование образования и диссоциации гидратов при разработке и эксплуатации нефтяных и газовых месторождений: автореф. дис... канд. техн. наук. Краснодар: Кубан. гос. технол. ун-т., 2015. 23 с.
28. Шумский П.А. Основы структурного ледоведения М.: Изд-во АН СССР, 1955. 492 с.
29. Эпов М.И., Ельцов И.Н., Оленченко В.В., Потапов В.В., Кушнаренко О. Н., Плотников А.Е., Синицкий А.И. Бермудский треугольник Ямала // Наука из первых рук. Вып. 5(59), 2014. С. 14-23.
30. Burges M., Judge A.S., Taylor A., Allen D.V. Ground temperature studies of permafrost growth at a drained lake site, Mackenzie Delta (MD). In Proceedings of the 4th Canadian Permafrost Conferences, Calgary, AB, Canada, 2-6 March 1982; pp. 3-11.
31. Leibman M.O., Kizyakov A.I., Plekhanov A.V., Streletskaya I.D. New permafrost feature — deep crater in central Yamal (West Siberia, Russia) as a response to local climate fl uctuations // Geography environment. 2014. V. 7, № 4. P. 68-80.
32. Mackay J.R. Pingos of the Tuktoyaktuk Peninsula Area, Northwest Territories // Geogr. Phis. Quart. 1979. V. 33, № 1. P. 3-61.
33. Mackay J.R. Pingo Growth and collapse, Tuktoyaktuk Peninsula Area, Western Arctic Coast, Canada: along-term fi eld study // Geographie physique et Quaternaire. 1998.
Vol. 52, No. 3, p. 271-323. 34. Majorowicz J., Osadetz K., Safanda J. Models of Talik, Permafrost and Gas Hydrate Histories-Beaufort Mackenzie Basin, Canada // Energies, 2015, vol. 8, no. 7, p. 67386764.
