Geophysical methods for the study of natural and human-induced changes in ground massifs of the permafrost zone

全文:

ВВЕДЕНИЕ

Температурный фактор определяет геомеханические и емкостные свойства трещиновато-пористых дисперсных и скальных грунтов криолитозоны. Глобальные и региональные изменения криолитозоны довольно медленны во времени и латерально неоднородны, что затрудняет их оценку. Однако, имеются инженерные объекты (АЭС, водохранилища, объекты ТЭК), которые создают локальную тепловую нагрузку на ММП, меняя их состояние вплоть до талого. Поскольку часть территории криолитозоны сейсмоактивна, очевидно, что и оценка сейсмического отклика при оттаивании ММП актуальна и должна проводиться с учетом геофизической изученности грунтового массива. Ранее был сделан вывод [9, 17, 18], что изменение фазового состава мерзлых пород, связанное с ростом температуры, должно вызывать существенное изменение их упругих свойств. Известно, что на многих гидротехнических сооружениях Западной Якутии [2, 3, 16, 17, 24] и ТЭЦ в г. Якутск [6] наблюдается развитие неблагоприятных инженерно-геокриологических процессов, влияющих на штатную эксплуатацию предприятий. Закономерности изменения упругих свойств мерзлых и талых грунтов могут быть установлены с помощью экспериментальных полевых и лабораторных исследований. Для изучения их состояния и внутреннего строения широко применяются геофизические методы. Были рассмотрены примеры с мерзлыми грунтами и расположенными на них инженерными объектами. На основе геофизических и инженерно-геологических данных, имеющихся для этих территорий, строились геосейсмические модели, исходя из представлений об эффективных сейсмогеокриологических моделях (ЭСМ), подробно рассмотренных в работе [13].

Были проанализированы закономерности изменения сейсмических свойств грунтов криолитозоны под действием техногенных факторов для решения задач сейсмического микрорайонирования (СМР). Рассмотрены примеры геофизических аномалий (тепловых, геоэлектрических, сейсмических) вблизи техногенных объектов, связанные с их отепляющим влиянием на мерзлые породы (водохранилище Сытыкан, каскад Вилюйских ГЭС и территория сооружений Главного корпуса и радиаторных охладителей Билибинской АЭС) [1–3, 16, 17, 24, 25].

Цель исследования — анализ ряда примеров геофизического контроля состояния мерзлых пород, подверженных отепляющему воздействию техногенных объектов.

ГЕОФИЗИЧЕСКИЙ КОНТРОЛЬ СОСТОЯНИЯ МАССИВОВ МЕРЗЛЫХ ГРУНТОВ НА ПРИМЕРЕ ГИДРОУЗЛОВ ЗАПАДНОЙ ЯКУТИИ

Известно, что особое место при анализе теплового воздействия инженерных сооружений на ММП занимают талики природного и техногенного происхождения. Им, в частности, отводится роль триггера в глобальных климатических эволюционных моделях деградации мерзлоты [28]. Талики являются ключевым агентом при региональном и глобальном изменении состояния ММП.

При эксплуатации многих гидротехнических сооружений в зоне распространения ММП (питьевые и гидротехнические водохранилища, хвостохранилища и другие промышленные объекты), существуют проблемы, обусловленные возникновением таликов техногенного происхождения [1–3, 26]. Рассматриваются примеры использования геофизических методов для обнаружения и мониторинга талых зон, приуроченных к гидротехническим сооружениям.

Детальные геофизические исследования [1, 26, 28] показали важную роль таликов природного и техногенного происхождения в ускорении таяния массивов мерзлых пород. Гидрогенные талики техногенного происхождения, с одной стороны, “головная боль” при эксплуатации многих гидротехнических сооружений в зоне ММП, а с другой — уникальные модельные объекты для геофизиков, гидрогеологов и геокриологов.

Наиболее мощное воздействие на тепловой режим мерзлых пород связано с сооружением плотин и созданием искусственных водоемов различного назначения. Старейшая в России ГЭС в зоне распространения ММП – Вилюйская, существует с 1957 г., и на сегодня она входит в комплекс каскада Вилюйских ГЭС-1, 2, 3. Наше внимание было сосредоточено на изучении талых зон, связанных с береговым примыканием водоемов и плотин гидротехнических сооружений, геофизическими методами [20, 24–26].

Предназначение плотин, включая плотины в зоне многолетней мерзлоты, —сохранение запаса воды при ее контролируемом пропуске. Как показывает практика эксплуатации водохранилищ в зоне ММП такой контроль не всегда удается [20]. Разнообразные противофильтрационные меры (мерзлотные завесы, противофильтрационные мембраны) лишь частично и временно эффективны. В этой ситуации перед геофизическими методами исследования, наряду с задачей выявления и прослеживания таликовых зон, стоит задача опережающего прогноза возникновения и развития талика [15]. Одной из важнейших задач остается контроль гидротехнических сооружений, где центральное место отведено выявлению и мониторингу развития таликовых зон, как связанных с плотиной, так и с береговыми примыканиями [14, 20].

Решение таких задач невозможно без понимания особенностей геофизических полей, обусловленных особенностями геокриологической среды, где развивается фазовый переход “лед–вода” [19]. Ниже приведен пример результатов геофизических наблюдений, характеризующих деградацию мерзлой среды, вмещающей питьевой водоем на р. Сытыкан. Сытыканский гидроузел создан в 1972–1974 гг. для водоснабжения г. Удачный и Удачнинского горно-обогатительного комбината (УГОКа). В конструктивном отношении плотина гидроузла является насыпной грунтовой, мерзлого типа (проект), водонепроницаемость которой обеспечивается мерзлотной завесой, создаваемой системой термосифонов с принудительным воздушным и жидкостным охлаждением. В 1995 г. в основании плотины зафиксировано развитие фильтрационных таликов, а в последующие годы отмечена резкая активизация их развития с прогрессирующим оттаиванием мерзлого основания плотины (рис. 1). Фильтрация развивалась и в протаявших зонах правобережного примыкания (рис. 2). В 2005 г. по результатам термометрии был зафиксирован прорыв таликом цементационно-мерзлотной завесы.

 

Рис. 1. Геологический разрез (а), нулевые изотермы и динамика изменения контура талика (б), радарограмма (в) по оси Сытыканской плотины, 1997–1998 гг. [1]

 

Рис. 2. Сытыканский гидроузел: а - распространение отрицательных аномалий потенциала естественного электрического поля (ЕП) в верхнем бьефе водохранилища; отрицательным аномалиям соответствуют зоны инфильтрации [1]; б - площадной контроль обходной фильтрации в талой зоне правобережного примыкании по совместной интерпретации данных скважинной термометрии и электротомографии [2].

 

Пример мониторинга таликовой зоны in situ в правобережном примыкании Сытыканского водохранилища [19, 24, 26], сложенного неоднородными трещиноватыми карбонатными породами, разной степени глинистости, с использованием метода односкважинного радиопрофилирования на двух частотах (ОРВП) дан на рис. 3. Наряду со снижением эффективных сопротивлений (ρ) в талых зонах, относительная диэлектрическая проницаемость (ε) талых зон многократно выше, чем в мерзлых. Повторные радиоволновые просвечивания (см. рис. 3), дают возможность получить наглядное представление о внутреннем строении техногенного талика в зависимости от динамики его развития и времени года. Интенсивность этих аномалий максимальна в летний период, что отличает их от аномалий, приуроченных к глинистым прослоям. Наблюдается снижение интенсивности обоих параметров (ρ, ε) по мере удаления от береговой линии, наиболее заметное в проницаемых зонах. Слоистая структура проводящих, таликовых водонасыщенных зон, отмечаемая в марте, свидетельствует об их сохраняющейся пониженной активности и после зимнего периода.

 

Рис. 3. Примеры геоэлектрических разрезов: слева – эффективных сопротивлений (ρ); справа – относительной диэлектрической проницаемости (ε), по профилю скважин (4–5–6–7–10) в мерзлом массиве в крест береговой линии Сытыканского гидроузла (Зап. Якутия); ОРВП наблюдения: вверху – март 2001; в середине – август 2001; внизу – март 2002. Для мерзлых пород характерно высокое ρ и низкое ε; для талых — низкое ρ и высокое ε. Уровень дна около плотины 301 м [19, 26]. По горизонтальным осям расстояние от береговой линии, по вертикальным — абсолютные отметки глубины, [м].

 

По результатам электротомографии в сочетании с заливом электролита в пьезометрическую скважину (рис. 4) был выявлен участок повышения электропроводности (временное сечение на 90 мин в районе ПК100), и определены направление и скорость фильтрации, которая составила 288 м/ч, что гораздо выше ожидаемых на этом участке скоростей [2].

 

Рис. 4. Схематическое расположение наблюдательного профиля и пьезометрической (закачной) скважины, правобережного примыкания водохранилища Сытыкан, Зап. Якутия (а). Режимные (повторные) электротомографические измерения для выявления фильтрационных окон и скоростей фильтрации. Окно фильтрации на глубине 25 м в районе ПК100м на 90 мин после закачки (б) [2].

 

Схожая ситуация — развитие техногенного талика, наблюдается в правобережном примыкании старейшего в России гидротехнического сооружения на ММП Вилюйского водохранилища ВГЭС-1,2 начиная с конца 1996 г. Зарождение обходной фильтрации наблюдалось в прибортовой зоне водосбросного канала, который являлся локальным источником тепла для возникновения проницаемых талых зон в первичномерзлых породах берегового примыкания [2, 3, 20, 24].

Комплекс геофизических методов для изучения зон фильтрации в правобережном примыкании ВГЭС-1,2 включал термометрию, электротомографию, метод естественного поля, геотехнические исследования фильтрационных свойств грунтов, георадиолокационные исследования [2, 3, 24, 26]. Мониторинг межскважинного радиоволнового просвечивания (РВГИ), проводимый ООО “Радионда” и охватывающий период более 20 лет, позволил осуществить опережающий температурный прогноз развития межмерзлотного талика в правобережном примыкании ВГЭС-1,2 [3] (рис. 5).

 

Рис. 5. Геоэлектрические разрезы по результатам РВГИ: 1999–2022 гг. с диаграммами гамма каротажа и электрического сопротивления. Приведена термометрия по скважинам за 1999, 2013 и 2022 г. [3]

 

Эти исследования в последние годы были дополнены регистрацией микросейсм, связанных с работой турбин ВГЭС-1,2 [7]. Анализ результатов измерений амплитуд микросейсмического сигнала в скважинах (рис. 6) показал, что для высокоскоростных слоев грунтовых толщ, температура которых ниже нуля (при прочих равных условиях), амплитуда сигнала существенно меньше амплитуды на растепленных низкоскоростных грунтах (разница при этом может составлять 25 дБ и более). Причем наблюдения, сделанные в разное время года (при наличии и отсутствии воды в грунте), позволяют выявлять зоны с аномально высокой амплитудой сигнала в талой части скважины, что дает возможность связать выявленные участки с областями повышенной трещиноватости грунтового массива. Эти области можно объяснить результатом влияния фильтрационных потоков на грунты правобережного примыкания Вилюйской ГЭС-1,2. Отмечено [7], что большая амплитуда сигнала фиксируется в основном в высокочастотной части спектра (110–200 Гц). В то же время наблюдаемые аномалии в мерзлых грунтах меньше по амплитуде (0.7–5.0 дБ), но частотный диапазон немного шире (70–210 Гц). В границах этих областей с учетом данных термометрии можно прогнозировать возникновение негативных инженерно-геокриологических процессов несмотря на то, что температура грунтов пока еще находится ниже нулевой отметки. Таким образом, как и в случае с радиоволновыми исследованиями, можно говорить об опережающем температурном прогнозе развития талика. Результаты исследований свидетельствуют о сложной внутренней картине взаимодействия вод гидросооружения с трещиноватой мерзлой–талой толщей.

 

Рис. 6. Результаты измерения амплитуд микросейсмического сигнала и температуры грунтов в скважинe NB-5 в 2019 и 2022 г.: 1 — насыпной грунт; 2 — долерит темно-серый, крепкий, сильнотрещиноватый, местами разрушенный до крупного щебня и глыб; 3 — переслаивание алевролитов, песчаников и доломита; 4 — известняк доломитизированный, местами разрушенный до супеси с щебнем и дресвой [7].

 

Рассмотренные примеры лишь часть комплексных геофизических работ на гидротехнических объектах Западной Якутии. Тем не менее они помогают понять пространственно-временные закономерности развития активных талых зон за относительно малый интервал времени. Одной из важных особенностей возникновения межмерзлотного активного талика в правобережных примыканиях ВГЭС-1,2 и Сытыкана следует считать наличие локального источника тепла, инициирующего его зарождение (рис. 2, 3) в береговом примыкании. На наш взгляд, таким источником является поток воды через обходной канал, греющий его дно и стенки. В сочетании с гидронапором он способствует прогрессивному развитию талика. Не следует исключать и наличие изначально воздушно-проницаемых мерзлых зон в трещиноватой среде берегового примыкания, связанных с тепловой историей последних.

ТЕХНОГЕННЫЕ ИЗМЕНЕНИЯ СЕЙСМИЧЕСКИХ СВОЙСТВ СКАЛЬНЫХ ММП ЧУКОТКИ НА ПРИМЕРЕ БИЛИБИНСКАЯ АЭС (БиАЭС)

Известно, что упругие свойства скальных грунтов зависят от состава, текстуры и структуры, трещиноватости, льдистости и теплового состояния породы [4, 5, 8, 12, 17]. В исследованиях авторов [17] были рассмотрены вопросы, связанные с особенностью проведения СМР на ММП. Было показано, что наличие выветрелой трещиноватой зоны в мерзлых скальных грунтах и ее обводнение при их оттаивании приводят к заметному изменению сейсмических свойств массива, и как следствие, к различному проявлению сейсмического эффекта. Были проанализированы закономерности изменения сейсмических свойств изначально мерзлых скальных грунтов на промплощадне БиАЭС. За время эксплуатации станции с 1976 г. произошли существенные изменения мерзлотных условий из-за растепления пород под основными сооружениями БиАЭС. При этом в основании тепловыделяющих сооружений (Радиаторный охладитель и Главный корпус) образовалась чаша оттаивания глубиной 30–40 м (1991 г.) (рис. 7).

 

Рис. 7. Пример изменения первоначального состояния мерзлых скальных грунтов по разрезу на промплощадке БиАЭС за время ее эксплуатации с 1976 г. Графики изменения температуры в скв. 8 (предположительно лето, 1991 г.), скв. 10 (предположительно зима, 1991 г.) и скв. 12 (режимные измерения в течение 2012 г.).

 

Территория исследования сложена трещиноватыми скальными породами верхнетриасового возраста, представленными песчаниками, алевролитами и аргиллитами. На рис. 7 приведен геологический разрез по профилю III–III с расположением режимных термометрических скважин (ВНИИОСП, https://niiosp.ru/ru/) на участке I очереди АЭС (скв. 7, 8, 10 и 12), последние две скважины – под Главным корпусом АЭС. В скв. 3, расположенной возле Радиаторного охладителя, максимальная температура в 1991 г. составляла +4°С, а в 2012 г. уже +8°С. Чаша оттаивания по данным термометрических скважин, расположенных вблизи Главного корпуса и сооружений Радиаторного охладителя, развивалась во времени, и за 20 лет глубина ее увеличилась примерно на 20 м, т.е. в среднем 1 м в год. Самые последние доступные результаты термометрии относятся к 2012 г. В термометрической скв. 8 глубина оттаивания на 1991 г. составляла 40 м от поверхности под сооружениями Радиаторного охладителя и 30–35 м (скв. 12) под Главным корпусом БиАЭС, а в 2012 г. в скв. 12 глубина оттайки достигла 55 м (см. рис. 7). Необходимо отметить неоднородность интенсивности прогрева пород под разными сооружениями БиАЭС, а именно, снижение температуры от Главного корпуса реакторного отделения до зоны Радиаторного охладителя, что и определило латеральную неоднородность оттаивания скального основания¹.

В табл. 1 приведено описание расчетных геосейсмических моделей:

 

Таблица 1. Расчетные геосейсмические модели для грунтовых комплексов при изменении мерзлотных условий под основанием тепловыделяющих сооружений БиАЭС за время ее эксплуатации

Модель 1. Многолетнемерзлая толща пород						Модель 2. Оттаявшая толща пород к 1991 г. на общую глубину 35 м
№		Vp	Vs	r	H, м	№		Vp	Vs	r	H, м
1	Насыпные грунты талые в летний период	0.87	0.49	2.1	4	1	Насыпные, дресвяно-щебенистые грунты	0.87	0.49	2.1	4
2	Многолетнемерзлые триасовые породы: сильнотрещиноватые	3.5	1.95	2.65	31	2	Оттаявшие сильнотрещиноватые триасовые породы	1.95	1.08	2.3	6
3	Многолетнемерзлые триасовые породы: трещиноватые	3.9	2.15	2.65	40	3	Оттаявшие трещиноватые триасовые породы	2.6	1.5	2.65	25
4	Многолетнемерзлые триасовые породы: слаботрещиноватые	4.1	2.3	2.7	∞	4	Многолетнемерзлые слаботрещиноватые триасовые породы	3.7	2.15	2.65	40
						5	Многолетнемерзлые породы	4.1	2.3	2.7	∞
Модель 3. Оттаявшая толща пород к 2012 г. на общую глубину 55 м						Модель 4. Оттаявшая толща пород к 2023 г. на общую глубину 65 м (прогнозная)
№		Vp	Vs	r	H, м	№		Vp	Vs	r	H, м
1	Насыпные, дресвяно-щебенистые грунты	0.87	0.49	2.1	4	1	Насыпные, дресвяно-щебенистые грунты	0.87	0.49	2.1	4
2	Оттаявшие сильнотрещиноватые триасовые породы	1.95	1.08	2.3	6	2	Оттаявшие сильнотрещиноватые триасовые породы	1.95	1.08	2.3	6
3	Оттаявшие трещиноватые триасовые породы	2.6	1.5	2.65	45	3	Оттаявшие трещиноватые триасовые породы	2.6	1.5	2.65	55
4	Многолетнемерзлые слаботрещиноватые триасовые породы	3.7	2.15	2.65	20	4	Многолетнемерзлые слаботрещиноватые триасовые породы	3.7	2.15	2.65	10
5	Многолетнемерзлые породы	4.1	2.3	2.7	∞	5	Многолетнемерзлые породы	4.1	2.3	2.7	∞

 

• модель 1 для грунтов в мерзлом состоянии;
• модель 2 для отложений в зоне отепляющего техногенного воздействия основных сооружений БиАЭС (Радиаторного охладителя и Главнонго корпуса) при оттаивании грунтов до глубины 35 м (1991 г.);
• модель 3 для отложений в зоне техногенного воздействия основных сооружений БиАЭС при оттаивании грунтов до глубины 55 м (2012 г.);
• модель 4 для отложений в зоне техногенного воздействия основных сооружений БиАЭС при прогнозном оттаивании к 2023 г. до 65 м;

При построении геосейсмических моделей 1–4 основной информацией послужили результаты геолого-геофизических исследований, в том числе скважинной сейсморазведки.

Так как задачей исследований являлось уточнение сейсмических условий площадки БиАЭС с учетом теплового воздействия ее сооружений на многолетнемерзлые скальные грунты, то были проанализированы изменения, произошедшие в мерзлом массиве под БиАЭС за несколько десятков лет ее эксплуатации. Были рассмотрены сейсмические свойства грунтовых комплексов, где произошло оттаивание ММП за 36 лет существования станции c 1976 г. по 2012 г. и при прогнозном оттаивании к 2023 г. (модели 2–4), а также участок с сохранившимися в летний период мерзлыми грунтами (модель 1).

По результатам метода МСЖ были получены значения приращения сейсмической интенсивности для участков оттаявших грунтов ΔI = +0.28 ÷ +0.38 балла при среднем значении ΔI = +0.34 балла.

Возможность частотной избирательности толщи грунтов оценивалась путем расчета амплитудно-частотных характеристик (АЧХ) с помощью компьютерных программ. Для расчетов АЧХ грунтов использовалась программа COEF-9, основанная на матричном методе, разработанная в ИФЗ РАН им. О.Ю. Шмидта Л.И. Ратниковой, А.Л. Левшиным [10]. По АЧХ в диапазоне периодов от 0.05 до 0.20 с (частотный состав местных землетрясений) в среднем приращение сейсмической интенсивности на поверхности оттаявших пород до глубины 55–65 м в зоне техногенного влияния основных сооружений относительно участка распространения многолетнемерзлых скальных грунтов оценивается в ∆I = +0.25 балла.

Следующая серия расчетов была проведена по известной программе NERA [21], которая позволила оценить изменение уровня сейсмических колебаний на поверхности оттаявших грунтов, по сравнению с ММП на площадке БиАЭС. Для расчетов была использована акселерограмма, синтезированная для поверхности мерзлого скального грунта по программе PSEQGN [27] с характеристиками местного землетрясения по данным ИФЗ РАН². На рис. 8 приведены фрагменты расчетных акселерограмм и соответствующие им спектры реакции (Sa).

 

Рис. 8. Примеры расчетных акселерограмм и соответствующих им спектров реакции Sa (5%) · 10^-2 в долях g от предполагаемого местного землетрясения: М1 (модель 1) - на поверхности толщи мерзлых скальных пород площадки БиАЭС (А_max = 9.1 см/с²); М2 (модель 2) - на поверхности оттаявшей толщи скальных пород при общей глубине оттаивания 35 м в 1991 г. (А_max = 10.0 см/с²); М4 (модель 4-прогнозная) - на поверхности при прогнозной общей глубине оттаивания 65 м к 2023 г. (А_max = 12.0 см/с²).

 

При этом приращения ∆I по расчетному методу СМР (по значениям Amax) на поверхности оттаявшей толщи при разных глубинах оттаивания относительно поверхности многолетнемерзлой толщи (модель 1) составили:

∆I = 3.3lg (10.01 / 9.1) ≈ +0.16 баллов (модель 2);

∆I = 3.3lg (11.02 / 9.1) ≈ +0.27 баллов (модель 3);

∆I = 3.3lg (12.0 / 9.1) ≈ +0.40 баллов (модель 4-прогнозная).

По результатам исследований было установлено, что при оттаивании ММП в области отепляющего действия Главного корпуса и Радиаторного охладителя скорости продольных и поперечных волн в зонах трещиноватости скальных пород триасового возраста уменьшаются в 1.40–1.55 раза.

За время эксплуатации станции в основании тепловыделяющих сооружений образовалась чаша оттаивания, максимальная глубина которой к 2012 г. составила 55 м. Значения приращений сейсмической интенсивности (∆I) на поверхности оттаявшей толщи трещиноватых скальных пород под главным корпусом БиАЭС, рассчитанное по значениям Amax в одномерном варианте и диапазоне периодов 0.05 ÷ 0.20 с от возможных местных землетрясений, составило относительно кровли многолетнемерзлой толщи ∆I = +0.27 балла и может возрасти до ∆I = +0.40 баллов (при прогнозной глубине оттаивания 65 м к 2023 г.).

Результаты расчетов ∆I по разным методам СМР показали близкие значения (по МСЖ +0.28 ÷ +0.38 балла, с учетом АЧХ +0.25 баллов, по значениям Amax +0.27 ÷ +0.40 балла).

В дальнейшем для уточнения сейсмической реакции латерально неоднородно оттаивающей толщи необходимо использование двух- и трехмерных сейсмологических моделей, как предлагается в работе [22].

ЗАКЛЮЧЕНИЕ

Опыт комплексного геофизического изучения массивов мерзлых пород в зоне гидротехнических сооружений важен, как при оценке упругих характеристик среды и сейсмической опасности этой зоны [11, 18], так и для понимания ее внутреннего строения. Накопленный исследователями за многие годы опыт изучения толщ ММП дает возможность лучше понимать механизмы их эволюции на современном этапе и оценивать деградацию мерзлых толщ в условиях глобального изменения климата.

При исследованиях авторов в 2022–2023 г. [16–18] и анализе материалов прошлых лет [4, 5, 13, 18, 23] рассмотрены закономерности изменения электрических, тепловых и сейсмических свойств различных грунтов, находящихся в мерзлом и талом состоянии. Исследования показывают, что образование талых зон в массиве ММП существенно изменяет в нем характеристики значения удельного электрического сопротивления (ρ), диэлектрическую проницаемость (ε) и сейсмического (упругие модули, коэффициенты затухания) полей. Меняются сейсмические условия и характеристики сейсмических воздействий (пиковые значения ускорений (PGA-Аmax), спектры реакции-отклика (Sa) и кривых динамичности (β)), которые необходимо учитывать при СМР.

Геофизический мониторинг системы “водоем – мерзлый массив – климат”, наряду с мониторингом ММП в естественных условиях, в случае гидротехнических сооружений дает детальное представление о развитии межмерзлотных таликов техногенной природы. Эта информация важна как для разработки прогностических моделей эволюции ММП, подвергающейся техногенному воздействию вблизи крупных инженерных объектов, так и моделей эволюции криосферы в условиях изменения климата.

¹ Отчет. Билибинская АЭС. I очередь. Реконструкция. Проект. Промплощадка. Технический отчет об инженерно- геологических изысканий по созданию режимной сети. БКИ МО “Атомэнергопроект”, М., 1992.

² Итоговый отчет по теме «Обоснование и уточнение оценок сейсмичности площадки размещения Билибинской АЭС», ИФЗ РАН, М., 2002.

作者简介

B. Trifonov

Sergeev Institute of Environmental Geoscience of the Russian Academy of Sciences

编辑信件的主要联系方式.
Email: igelab@mail.ru
俄罗斯联邦, Bldg. 2, 13, Ulansky All., Moscow, 101000

S. Milanovsky

Sergeev Institute of Environmental Geoscience of the Russian Academy of Sciences; Sсhmidt Institute of Physics of the Earth of the Russian Academy of Sciences

Email: svetmil@mail.ru
俄罗斯联邦, Bldg. 2, 13, Ulansky All., Moscow, 101000; Bldg. 1, 10, B. Gruzinskaya St., Moscow, 123242

V. Nesynov

Sergeev Institute of Environmental Geoscience of the Russian Academy of Sciences

Email: igelab@mail.ru
俄罗斯联邦, Bldg. 2, 13, Ulansky All., Moscow, 101000

参考

补充文件