CC BY
25
- © В.Н. Морозов, А.И. Каган,
В.Н. Татаринов, 2015
УДК 621.039
В.Н. Морозов, А.И. Каган, В.Н. Татаринов
ГЕОДИНАМИЧЕСКИЕ АСПЕКТЫ РАДИАЦИОННОЙ БЕЗОПАСНОСТИ НА АРХИПЕЛАГЕ НОВАЯ ЗЕМЛЯ
Представлены результаты моделирования напряженно-деформированного состояния пород в районе пролива Маточкин шар архипелага Новая Земля. Выполнен анализ возможного развития тектонического процесса и инфильтрации подземных вод на основании структурно-тектонической модели и расчета напряженно-деформированного состояния на локальных площадях на территории архипелага, потенциально пригодных для подземной изоляции радиоактивных отходов. Ключевые слова: напряженно-деформированное состояние, моделирование, метод конечных элементов, Нижнеканский массив.
Введение
На ядерном полигоне архипелага Новая Земля с 21 сентября 1955 г. по 24 октября 1990 г. (объявление моратория) было проведено 132 ядерных взрыва: 87 атмосферных (в том числе 84 воздушных), 1 наземный (7 сентября 1957 г.), 2 надводных (27 октября 1961 г. и 22 августа 1962 г.), 3 подводных (21 сентября 1955 г., 10 сентября 1957 г., 23 сентября 1961 г.) и 42 подземных. Подземные ядерные взрывы (ПЯВ) проводились в горизонтальных выработках и шахтных стволах на глубинах превышающих 400-500 м от поверхности [1]. При этом использовался подходящий рельеф и горно-геологические условия, обеспечивающие экологическую безопасность ПЯВ. Вместе с тем, геолого-тектонические условия района, сейсмотектонические процессы имели второстепенное значение при выборе мест заложения ПЯВ.
Архипелаг Новая Земля, являющийся продолжением Уральской горной системы, находится в условиях активного тектогенеза. Локальная сеть сейсмических станций, функционирующая последние 20 лет в этом районе позволила зафиксировать ряд землетрясений, в том числе и землетрясения М = 4,6 с эпицентром в районе пролива Маточкин шар. Причиной сейсмо-тектогенеза в этом районе, вероятнее всего, являются тектонические напряжения сжатия с векторами усилий широтного простирания. Как и в пределах Центрального Урала можно ожидать, что напряжения сжатия могут достигать величин 30-40 МПа. Новые возможности анализа напряженно-деформированного состояния (НДС) породных массивов в сочетании с космическими съемками и сейсмическими исследованиями позволяют выполнить оценку развития тектонического процесса, как с позиции экологической безопасности зон проведения ПЯВ, так и использования барьерных свойств многолетне-мерзлых пород (ММП), для создания объекта длительного хранения и окончательной изоляции РАО в геологических формациях.
ММП на территории архипелага Новая Земля имеют мощность до 300 м и можно допустить в течение последних 15-20 тыс. лет образование протяженных шлировых структур, которые генетически связаны с разломной тектоникой архипелага. Пластичные свойства льда, в сочетании с пленочной незамерзшей
водой в градиентных полях тектонических напряжений способствуют формированию сутуровидных структур льда, определяющих механизм возможной миграции радионуклидов в условиях возможного глобального потепления климата.
Принятие федерального закона «Об обращении с радиоактивными отходами...», а также предполагаемое расширение транспортных потоков по Северному морскому пути повышают актуальность вопросов обеспечения экологической безопасности в регионе, в том числе вопросов создания и безопасного размещения пунктов долговременного хранения и изоляции РАО. Перспектива дальнейшего использования полигона определяет целесообразность геоэкологического анализа устойчивости геологической среды северной части Южного острова с позиции геодинамики, включающего:
• Разработку структурно тектонической модели центральной части архипелага, в т.ч. палеотектоническую реконструкцию тектонического процесса;
• Моделирование НДС центральной части архипелага с использованием метода конечных элементов;
• Анализ НДС состояния центральной части архипелага с позиции дальнейшего развития тектонического процесса и прогноза формирования ледопо-родных структур в градиентных полях тектонических напряжений;
• Прогнозные оценки возможных путей фильтрации водных потоков при таянии ММП в условиях глобального потепления климата.
Целью настоящей работы является анализ возможного развития тектонического процесса в пределах архипелага Новая Земля на основании структурно-тектонической модели и расчета НДС локальных площадей на территории архипелага, потенциально пригодных для изоляции РАО.
Результаты исследований
Архипелаг Новая Земля входит в Пайхойско-Новоземельский складчатый комплекс, представляющий сложно построенное «складчатое сооружение» дугообразной формы. В пределах архипелага выделяются два крупных поднятия в южной и северной частях островов. Южно-Новоземельская и Северо-Ново-земельская антиформы разделяются Карамакульским прогибом. Структурами более мелкого порядка являются Пахтусовская и Литкенская антиформы [1].
Южно-Новоземельская антиформа (рис. 1) длиной 300 км с северо-востока разрезается глубинным Байдарицким разломом. Особенностью этой структуры являются широко развитые изоклинальные складки, образовавшиеся в результате горизонтального смещения надвинутой Карской зоны в западном направлении. Разломы, осложняющие структуры имеют субвертикальные плоскости смещения. Северо-Новоземельская антиформа, протяженностью 500-600 км, сформировалась в результате субширотного сжатия с доминированием изоклинальной складчатости и образованием надвиговых чешуй. Вдоль крутых разломов образуются широкие (10-15 км) шовные зоны, которые проходят по границам сопряженных Северо-Новоземельской антиформы и Карской впадины.
Геологические структуры Пайхойско-Новаземельского региона более молодые по сравнению с уральскими структурами, с которыми имеют определенное сходство. Время формирования Новаземельской складчатости датируется рубежом триаса и юры. На рубеже триаса и юры весь регион испытал складчатость и скучивание, сопровождающиеся внедрением гранитойдов с возрастом 210-180 млн лет [1]. Современная тектоническая структура сформировалась в период раннекиммерийского тектогенеза.
м- 60- 44= Ыг= сд
Рис. 1. Тектоническая схема района архипелага Новая Земля [1]: 1 - континентальная плита с чехлом верхнерифейских-нижнекембрийских пород; 2 - фрагменты океанического метабазитового комплекса и гранитоидов; 3, 4 - докембрийские комплексы: 3 - ме-зопротерозойский метабазит-ультрабазитовые (а) и мезоп1ротерозойские-неопротерозойские гранитоидные комплексы (б), 4 - неопротерозойские комплексы: русановский метабазитовый (а), митюшевский гранитоидный (б), раховский гранитоидный (в); 5 - байкальский орогенный пояс; 6 - кембрийско-каменноугольные шельфовые чехлы эпигренвильской и эпибайкальской плит: 7-12 - Новоземельский сегмент герцинского орогенного пояса: 7 - Центрально-Ново-земельский предпозднедевонский аллохтон, 8 - позднедевонский рифтогенный терригенно-вулканогенный комплекс, 9 - Восточно-Новоземельский позднепалеозойский аллохтон, 10 -верхнепермская-нижнетриасовая пестроцветная моласса, 11 - участки траппового магматизма пермо-триаса, 12 - массивы гранитоидов, преимущественно позднепермского возраста; 13-15 - надвиги герцинского орогенеза: 13 - Главный Новоземельский, 14 - Восточно-Но-воземельский, 15 - прочие; 16 - разломы Байдарацкой зоны сдвигов; 17 - сдвиги и надвиги прочие; 18 - границы стратиграфические
Рис. 2. Эпицентры землетрясений в арктической материковой окраине России
Анализ сейсмических явлений в регионе (рис. 2) позволяет выделить сейсмическую зону, которая трассирует географическое положение пролива Маточ-кин шар. Новоземельское землетрясение 01,08,1986 г. с М = 4,6 с эпицентром вблизи южного берега пролива Маточкин шар дает представление о тенденции современного развития тектонического процесса. Механизм очага землетрясения по данным о полярности Р-волн соответствует взбросу (взбросовой подвижке) в субмеридианной плоскости разрыва под углом 42° на глубине 15 км.
Если воспользоваться корреляционной связью между мощностью земной коры и магнитудой возможного землетрясения в виде [2]:
М = 0,087Н + 2,7, (1)
полученная для Монголо-Охотской зоны с мощностью земной коры до 50-55 км, то для архипелага Новая Земля получим величину
М = 6,18 при Р = 40 км. (2)
Далее, используя корреляционные связь между амплитудой землетрясения и вероятной длинной тектонического разлома, имеем:
. М - 4,76 . .
ь =--— = 473 км. (3)
0,003
Если исходить из потенциальной возможности реализации землетрясения в районе пролива Маточкин шар М = 5, получаем длину разлома равную 80 км.
В работе [2] приводится уравнение регрессии связи между возрастом пород фундамента и возможной максимальной магнитудой землетрясения:
1) М = 0,848 ■ Т + 5,49 (разлом Танлу),
2) М = 0,95 ■ Т + 2,58 (Моголо-Охотский район), где Т - миллиарды лет.
В пределах архипелага Новая Земля наиболее древними породами считаются нижнепротерозойские-нижнерифейские образования, метаморфизованные в эпидот-амфиболитовой фации породы района губы Северная Сульменева на западном побережье северного острова. Изотопное определение возраста (и-РЬ) по цирконам составило для амфиболитов и кристаллических сланцев
1550±80 и 1490±100 млн лет, для гранитов, инъецинирующих кристаллические сланцы - 1300±90 млн лет. В то же время возраст, определенный по амфиболу и биотиту (K-Ar) дает значение 645±50 и 587±27 млн лет, что, видимо, отражает возраст позднего метоморфизма.
По оптимистическому сценарию развития процесса тектогенеза в районе, положим значение возможной магнитуды землетрясения М = 5,96. Тогда, допуская возможность землетрясений с М = 6, существует вероятность пересечения острова тектоническим разломом длиной не менее 100 км. То есть, имеется опасность возникновения протяженных разломов и в большей степени это касается северной части Южного острова и южной Северного, т.е. районов прилегающих к проливу.
Потенциальная возможность возникновения тектонического разлома связана с высокой концентрацией напряжений, возникающих под действиями тектонических сил. В качестве модели локализации тектонических напряжений рассмотрим центральный район архипелага, разделенного проливом Маточ-кин шар, площадью около 10 км2 (рис. 3).
Превалирующим напряжениеми в регионе являются напряжения сжатия, ориентированные в субширотном направлении, поэтому в качестве рабочей модели мы приняли тектоническую схему, представленную на рис. 3, задав для определенности тектонические напряжения ст = 50 МПа и ст = 10 МПа (ось
1 1 хх уу v
Y направлена вдоль меридиана, ось X - перпендикулярно). Далее была использована программа расчета НДС, детально изложенная в работе [4].
На рис. 4 и 5 показана концентрация напряжений сдвига и интенсивности напряжений соответственно. Максимальные напряжения сдвига потенциально опасны с точки зрения формирования тектонических подвижек, т.е. образования тектонических разломов содержащих сдвиговую компоненту. Опасность возникновения разломов наиболее вероятна в зонах высокой концентрации т .
И в первую очередь это касается юго-западной окраины Северного острова, в зоне примыкающей к «излому» пролива Маточкин шар. Возможность форми-
Рис. 3. Структурно-тектоническая модель центральной части архипелага Новая Земля [1]: 1 - разломы по геологическим данным; 2 - границы блоков, выраженные различными типами; 3 - очаг землетрясения. Стрелки показывают направление действия тектонических напряжений
xy
Рис. 4. Распределение касательных напряжений т
Рис. 5. Распределение интенсивности напряжений
рования протяженного разлома с выходом в акваторию океана в этом районе наиболее вероятна.
Видно, что расчет НДС позволяет выделить зоны концентрации напряжений, потенциально опасные с позиции возникновения тектонических нарушений, так и зоны разрядки напряжений наиболее безопасные с позиции размещения хранилищ и могильников РАО. Исходя из представлений о наследственности развития тектогенеза в пределах архипелага Новая Земля, следует допустить непрерывный процесс «подкачки» в зоны их аномальной концентра-
ции с момента образования зоны многолетнее-мерзлых пород, т.е. в течение последних 8-10 тыс. лет. На это указывает современная сейсмичность района. Как известно, фильтрация через среду описывается законом Дарси:
<< = , (4)
¿х
где О - расход потока; а - площадь поперечного сечения потока; ду/дх - напорный градиент потока; К - коэффициент фильтрации. Закон Дарси можно записать в следующем виде:
^х , (5)
х ц ах
где V - величина скорости фильтрации через пористую среду; к - проницаемость среды; ц - динамическая вязкость жидкости или газа; Р - давление.
Коэффициент к имеет размерность скорости, и, если предположить, что выполняется условие полного насыщения, он зависит от геометрии пористого пространства (типа среды и характеристик пор), а также от удельного веса и вязкости жидкости. Он постоянен для данной жидкости, если пористая среда несжимаема и изотропна.
Для изотропной среды как среды, имеющей одну и ту же проницаемость во всех направлениях можно записать:
V =-ко ¿р,V =-ко ар,V =-ко ар
х ц ¿х' у ц ф' 2 ц ¿г ; или
V = -Ц V*Р), (6)
м-
где V - вектор скорости фильтрации. Применив операцию дивергенции к предыдущему уравнению и используя уравнение
^ ЭУх дУх
V • V = —х + —- + —2 = о
дх ду дг
получаем:
- ко V2 (Р) = о . (7)
ц
Таким образом, фильтрация сквозь изотропные среды может быть сведена к решению уравнения Лапласа с соответствующими граничными условиями. Если распределение Р известно, то скорости фильтрации могут быть получены из закона Дарси в форме (6).
Во многих случаях пористая среда анизотропна и проницаемость зависит от направления течения. В этом случае компоненты вектора скорости фильтрации можно записать в виде:
V = Кар. V = Кар. V = КОР
х ц ¿х у ц ¿у' 2 ц & ' (8)
Перейдя к векторной форме и вычислив дивергенцию, получим:
V•V = -кох др - V др - ког д2р ц дх2 ц ду2 ц дг2 .
Предыдущие уравнение может быть записано в таком виде:
д2Р д2Р д2Р .
—Т +—2 +—Т = 0 дх ду дz
Таким образом, мы снова получили уравнение Лапласа. Следовательно, истинный физический случай можно представить как фиктивный изотропный в преобразованных координатах.
Связь между компонентами деформации и смещения в плоской среде имеет вид:
Гди
е =
дх дК
ду
ди + ду
дх ду
(9)
В соответствии с законом Дарси будем исходить из предположения, что скорость фильтрации в градиентном поле тектонических напряжений пропорциональна градиенту действующих напряжений:
Vf =
или
V =
V
V
ив ди
дх
кв ™
ду
'Кх ' Г кВУ2и'
' кВУ2У _
(10)
в предположение некоторого постоянного коэффициента фильтрации к.
Скорость фильтрации, выраженная через компоненты напряжений, в условиях плосконапряженного состояния, имеет вид:
К = Ух + V = кв (ахх + ауу ) (11)
Для условий анизотропной проницаемости массива (кх ф ку ) .
V = V + У = в (кх а хх + ку а уу ) .
Тогда используя эти соотношения, представляется возможность расчета и прогнозирования развития гидрогеологических условий в районе заложения могильника РАО при изменение НДС породного массива и изменение климатических условий.
Обратим внимание, что вектор скоростей фильтрации ортогонален направлению действия сил, вызывающих компрессию вдоль соответствующих осей в выбранной системе координат. Уровень напряженного состояния породного массива в его локальных зонах определен величиной интенсивности напряжений. Исходя из предположений, что скорость фильтрации определяется этой
интегральной оценкой уровня действующих напряжений в элементе объема породного массива имеем:
Ц = дгаё (ст хх + ст ^^)
(12)
Это упрощает расчеты первого приближения, поскольку задача оценки напряженного состояния и фильтрации решаются раздельно. Зависимость скорости фильтрации от скорости развития деформационного процесса следует из ранее сделанных предположений. Если известна (из экспериментальных данных) скорость изменения тектонических напряжений в пространстве и времени то выражение (12) приводится к виду:
V (t) = кЕУ
дстх дt
с.
(12)
В этом случае скорость фильтрации, как функция координат и времени Щ(х, у, $ дает возможность оценок водопритоков в ослабленные зоны тектонической трещиноватости, как основу дальнейших расчетов скорости миграции и прогнозируемой миграции радионуклидов в процессах массопереноса.
Величина водопритока со скоростью фильтрации, зависящей от величины интенсивности напряжений можно записать в следующем виде:
О = БкУст1 (X)
К,
(13)
где 5 - площадь поперечного сечения (элементарная площадка); к - коэффициент фильтрации, при линейном изменении ст. на участке Ьх - Ьх , а следовательно рассчитать потенциальную опасность водопритока при растеплении ММП, используя экспериментальные данные для конкретного района исследований.
Рис. 6. Направление векторов скорости фильтрации при растеплении массива ММП
Фактически возможность раздельного решения упругой задачи и процесса фильтрации под действием локальных градиентов тектонических напряжений (т.е. без учета деформации порового скелета горных пород) существенно упрощает алгоритм расчета возможных направлений фильтрационных потоков в условиях ММП при их растеплении. Этот алгоритм расчета градиентов напряжений введен в программу расчета НДС массива. На рис. 6 представлены расчеты векторов
^ = grad ^
2 max
максимальных градиентов действующих тектонических напряжений.
В условиях отрицательных температур эти направления соответствуют перемещениям незамерзшей и пленочной воды в ослабленные зоны с формированием сегрегационных залежей льда. Как это видно из рис. 6, максимальные векторы градиентов действующих напряжений направлены внутрь тектонических нарушений.
Это значит:
При возникновение концентрации в вершине разлома (в зависимости от пространственного ориентирования), в условиях отрицательных температур, тектонические разломы (зоны тектонической нарушенности) являются аккумуляторами сегрегационного накопления льда.
При растепление породных массивов (предположительно за счет потепления климата) эти зоны, в соответствие с законом Дарси, являются «резервуарами» стока грунтовых вод.
При наличие градиента давления в породном массиве эти разломы собой каналы интенсивной фильтрации грунтовых вод, способствующие интенсивной миграции радионуклидов в геологической среде.
При растепление породных массивов ММП зоны концентрации тектонических напряжений опасны не только с позиции возникновения новых тектонических нарушений, но и способствуют возникновению ориентированных потоков воды в направление площадей, соответствующих более низким значениям интенсивности тектонических напряжений.
Эти общие фундаментальные закономерности связи между тектоническим напряжением и потенциальной угрозы массопереноса радионуклидов в процессе интенсивной фильтрации грунтовых вод, требуют определенной детализации с позиции разработки инженерных методов расчета радиационной угрозы для окружающей среды.
Потенциальная опасность пересечения разлом площади могильника РАО на территории островов, может быть оценена величиной физического риска возможности одномоментного эпизодического природного воздействия (землетрясение, оползень и др.) по методике изложенной в работах [5, 6].
Физическая уязвимость территории определяется как:
V (H ) = SH • S-1, (14)
где SH - площадь пораженная опасностью H, S( - площадь таксона.
Удельный риск поражения:
R (H) = P* (H) • V (H) = P' (H) • Sh • S-1 (15)
где P*(H) - вероятность события на территории SH.
Площадь архипелага Новая Земля составляет 83 тыс. км2. Эту площадь принимаем за величину SH. Площадь таксона - площадь района расположения возможного места подземной изоляции ВАО. Условно площадь таксона принимаем равной 800 тыс. км2.
Исходя из данных анализа тектонической нарушенности архипелага Новая Земля тектоническими разломами, принимаем их количество равным 1000, сформировавшихся в течение кайнозоя, т.е. за последние 20 млн лет. Получаем значение вероятности возникновения протяженного тектонического нарушения
P' (H ) = -1037 = 5 • 10-5
v ' 2•107
Подстановка P*(H) и принятых значений SH и S(, приводит к величине ожидаемого риска:
R (H) = 5 • 10-5 • 102 = 5 • 10-3 ,
как среднюю величину риска возможности возникновения тектонического разлома на площади порядка 800 тыс. км2. Ситуация осложняется следующим обстоятельством. Пролив Маточкин шар трассируется гипоцентрами неглубоких тектонических землетрясений (5-15 км - глубина гипоцентров) (см. рис. 3). Это касается и Южной окраины архипелага. При магнитуде одного из сильных тектонических землетрясений на Новой Земле M = 6 можно ожидать тектонический разлом L>100 км.
Если предположить, что это событие имеет периодичность один раз в 1000 лет, следует допустить, что вероятность формирования протяженного тектонического разлома в районах архипелага, примыкающих к проливу Маточкин шар, значительно выше, т.е. составляет
P * (H) = 10-3
При площади могильника РАО равной 8 км2, риск достигает величины
R2 (H ) = 10-3
Эти оценки, разумеется, не являются адекватными без расчета реальной угрозы распространения радионуклидов в окружающую среду, но они однозначно требуют постановки и решения соответствующей задачи, приведенной в последующем изложении.
Таким образом, архипелаг Новая Земля находится в зоне активных тектонических процессов. Рассмотрим потенциальную опасность размещения пункта долговременного хранения РАО в верхней части земной коры, размещенного в штольне на глубине 100 м над уровнем океана. Допустим, что в результате тектонической активизации возникает зона тектонического сдвига, захватывающая континентальную часть архипелага и прилегающих акваторий при ширине D, значительно превышающей геометрические размеры могильника.
На первой стадии формирования зоны сдвига, следует допустить формирования с образованием микропрослоев льда в его центральной части, т.е. формирование зоны сегрегации льда под действием напряжений сдвига. На второй, завершающей стадии, образуется зона сдвига (шириной D>200 м), захватывающая могильник РАО. Глубина проникновение зоны сдвига значительно ниже отметки дна прилегающей акватории. За счет высокой диспергации пород зоны сдвига морская вода постепенно ее заполняет с верхней отметкой, соот-
ветствующей уровню океана. На уровне порядка 1 км при нормальном температурном градиенте температура морской воды около 10 °С, в то время как на уровне моря температура воды соответствует температуре воды в океане (условно 5 °С).
Выше в зоне тектонической деструкции ММП имеем среднюю температуру порядка -2 °С. Эта физическая модель лежит в основе прогнозной модели, описывающей процесс растепления горных пород в ближней зоне могильника и опасность массопереноса радионуклидов в окружающую среду и в первую очередь в прибрежные воды архипелага Новая Земля.
Используя допустимые упрощения поставленной задачи, в окончательном виде сформулируем ее следующим образом.
Начальные условия
• Средняя температура мерзлой толщи на гулубине 0 - (-200 м) - -2 °С.
• На глубине (-200 м) - (-1000 м) начальная температура + 5 °С.
• На глубине свыше 1000 м температура + 10 °С.
• Температура РАО соответствует температуре окружающих пород - (-2 °С).
Эти условия соответствуют моменту времени t0 = 0, теплопроводность не
зависит от координаты Z, а температура от X.
Задача сводится к оценке возможности растепления сегрегированного льда в зоне могильника с последующим переносом радионуклидов в зону ниже отметки -200 м за время 103-4 лет. В связи с этим, следует допустить потенциальную угрозу изменения НДС зоны сдвига и возможность обратного отжима морской воды в океанические воды акватории. При этом недопустимым уровнем опасности является растепление ММП в зоне могильника до температуры 0 °С за время меньшее периода полураспада наиболее долгоживущих РАО.
Практическую значимость при решении данной задачи приобретает определение максимум глубины заложения могильника с известным периодом полураспада РАО в предположении описанной геодинамической ситуации, т.е. величины порядка H от уровня океана
Если принять во внимание, что максимальная скорость фильтрации вниз составляет 1-10-7 м/сек или 3 м/год, то в процессе фильтрации радионуклиды из могильника, находящегося выше уровня океана на 100 м, в самом неблагоприятном варианте достигнут уровня грунтовых вод менее чем через 35 лет. Разумеется, эта оценка соответствует пессимистическому сценарию развития событий, т.е. не учитывается период полураспада изотопов, не стационарность скорости вертикальной фильтрации, роль конвективного потока и др.
Опасность растепления зоны ММП даже мощностью порядка 500 м и расположение могильника и хранилища РАО в приповерхностной зоне, вероятность достижения радионуклидов уровня океана высока и составляет около 160 лет.
Расчеты возможного растепления горного массива над слоем океанической воды можно существенно упростить, если предположить, что температура на поверхности массива T и температура на уровне океана To(H) постоянны по времени, а также коэффициент температуропроводности одинаков по всей толще диспергированного массива (А = const). Уравнение температуропроводности:
ТТ"^ = 0 (16)
dt dx2 (16)
при этих условиях приводится к виду (одномерная температуропроводность)
дг
дх
= 0 (17)
Г х = ОТ = Т1 при \х = НТ = Т (н) •
Интегрирование с выписанными граничными условиями приводит:
Т = А ■ х + В или
Т (х ) = ( дтаёнТ )■ х + Т3 (18)
Следовательно на глубине заложения могильника Нм, температура в условиях стационарного потока тепла будет равна:
Тм = ( д^нТ )■ Нм + Т
Нулевая изотерма соответствует глубине
Н (0) = - Т
дг^нТ (19)
Рассмотрим вариант наиболее безопасного размещения РАО, т.е. допускаем, что штольня с РАО расположена выше уровня океана на величину + 400 м находится в ММП. Тектонический разлом с шириной значительно превышающей размеры ближней зоны могильника захватывает акваторий. На начальной стадии формирования дислокации, как это отмечалось выше, формируется система микротрещин, сопровождающаяся сегрегацией льда в центральной части разлома.
При образовании собственно сдвига за счет трения выделяется достаточное количество тепла для растепления породного массива в зоне сдвига, при этом контейнеры с РАО повреждены и начинается процесс миграции радионуклидов в окружающее пространство. Считаем известным следующие параметры:
Скорость выщелачивание или растворения радионуклидов в воде, а следовательно начальную концентрации радионуклидов в водном растворе С0.
Скорость инфильтрации грунтовых вод V по вертикальным трещинам, образовавшихся в результате тектонического сдвига.
Для «пессимистического» сценария принимаем, что скорость инфильтрации равна скорости переноса радионуклидов в зону разлома, заполненную морской водой до нулевой отметки.
Задача сводится к расчету радиоактивного загрязнения морской воды в зоне разлома, имеющего выход в акваторий.
Вместе с тем, в течение времени отрицательные температуры на поверхности разлома постепенно проникают на глубокие горизонты, возвращая температурные условия в исходное состояние. Таким образом, возникает необходимость рассчитать скорость продвижения 0 °С изотермы от поверхности к горизонту заложения могильника РАО.
По мере продвижения фронта отрицательных температур с поверхности охлаждение воды идет через образующиеся прослои льда. В работе [7] приведено решение системы дифференциальных уравнений
дТ д2Тс
— а
д " ^ , (20)
дТл = д2Т„
д л дг
2
где Тв, Тл - температура воздуха и льда, ав, ал - соответственно коэффициенты теплопроводности, £ - время, г - координаты. Время кристаллизации определяется
. Ь2 V ■ Ь
. —---
кр -т-
ав ' ал ав ' Т в (21)
где V - скрытая кристаллизация воды равная 0,33-109 Дж/м2.
Если воспользоваться этим решением, то значение скорости продвижения фазового перехода (изотермы 0 °С) при температуре воздуха -10 °С вглубь массива ориентировочно составляет 7-10-3 м/час. При заложение могильника или хранилища РАО на глубине 300 м от поверхности, РАО будут находиться при плюсовой температуре в интервале времени около 46 лет.
Таким образом, можно рассчитать время, в течение которого радионуклиды способны к миграции в рамках модели адвекции-диффузии в условиях положительных температур локального объема ММП.
Заключение
Расчеты напряженно-деформированного состояния центральной части архипелага Новая Земля в сочетание с результатами сейсмологических исследований дают основания предположить возможность возникновения новых тектонических разломных зон протяженностью более 100 км. Локальные зоны высокой концентрации тектонических напряжений являются источниками возможного возникновения и последующего распространения зон тектонической деструкции. Заполнение этих зон морской водой будет способствовать растеплению ММП с риском возможной миграции радионуклидов. Вместе с тем предварительные расчеты показывают, что радиационная безопасность мест подземной изоляции РАО зависит от их расположения относительно океана. При этом, чем ниже уровень действующих напряжений в структурно-тектонических, тем выше уровень их надежности.
Палеотектоническая реконструкция действующих тектонических напряжений не является достаточной, поскольку она постулирует наследственность тектонического процесса в его современной фазе. Возможности методов космической геодезии были бы хорошим дополнением в прогнозах развития локальной концентрации напряжений. Создание геодезического полигона и специализированные геолого-геофизические исследования (включающие натурные измерения напряжений в горных выработках и скважинах) способствовали бы более высокому уровню прогнозных оценок, в том числе и на длительный период времени, а также более обоснованной разработке критериев приемлемости РАО для предполагаемых к созданию объектов обращения, долговременного хранения и изоляции РАО на архипелаге Новая Земля.
СПИСОК ЛИТЕРАТУРЫ
1. Филатова Н.И., Хаин В.Е. Кратон Арктида и неопротерозойские-мезозойские ороген-ные пояса циркумполярного региона // Геотектоника. - 2010. - № 3. - С. 3-29.
2. Геология и полезные ископаемые России. Т. 5. Арктические и дальневосточные моря. Кн. 1. Арктические моря / Под ред. И.С. Грамберга, В.Л. Иванова, Ю.Е. Погребицкого. -СПб.: Изд-во ВСЕГЕИ, 2004. - 468 с.
3. Николаев В.В., Врублевский А.А., Ахмадулин В.А., Кузнецов В.Е. Геодинамика и сейсмическое районирование материковой части Дальнего Востока. - Владивосток: ДВО РАН, 2000. - 90 с.
4. Морозов В.Н., Колесников И.Ю., Татаринов В.Н. Моделирование уровней опасности напряженно-деформированного состояния в структурных блоках Нижнеканского гранито-идного массива (к выбору участков захоронения радиоактивных отходов) // Геоэкология. -2011. - № 6. - С. 524-542.
5. Морозов В.Н., Татаринов В.Н. Методика выбора участков земной коры для размещения экологически опасных отходов // Геоэкология. - 1996. - № 6. - С.109-120.
6. Морозов В.Н., Колесников И.Ю., Белов С.В., Татаринов В.Н. Напряженно-деформированное состояние Нижнеконского гранитоидного массива - района возможного захоронения радиоактивных отходов // Геоэкология. - 2008. - № 3. - С. 232-243.
7. Росман Г.И., Быховских Л.З., Самсонов Б.Г. Хранение и захоронение радиоактивных отходов. - М.: ВИМС, 2004. - 240 с. щгд
КОРОТКО ОБ АВТОРАХ_
Морозов Владислав Николаевич - доктор технических наук, главный научный сотрудник, e-mail: v.morozov@gcras.ru,
Каган Александр Иосифович - старший научный сотрудник, e-mail: a.kagan@gcras.ru, Татаринов Виктор Николаевич - доктор технических наук, зав. лабораторией, e-mail: v.tatarinov@gcras.ru, Геофизический центр РАН.
UDC 621.039
GEODYNAMIC ASPECTS OF RADIATION SAFETY AT ARHIPELAG NOVAYA ZEMLYA
Morozov V.N.1, Doctor of Technical Sciences, Chief Researcher, e-mail: v.morozov@gcras.ru, Kagan A.I.1, Senior Researcher, e-mail: a.kagan@gcras.ru,
Tatarinov V.N.1, Doctor of Technical Sciences, Head of Laboratory, e-mail: v.tatarinov@gcras.ru, 1 Geophysical Center of the Russian Academy of Sciences, 119296, Moscow, Russia.
The results of the modeling of stress-strain state rocks in the area Matochkin Shar of archipelago Novaya Zemlya. The analysis of the possible development of a tectonic process of infiltration and groundwater based on structural-tectonic model and calculation of the stress-strain state in local areas on the archipelago, potentially suitable for underground disposal of radioactive waste.
Key words: stress-strain state, modeling, finite element method, Nizhnekansky array.
REFERENCES
1. Filatova N.I., Khain V.E. Geotektonika. 2010, no 3, pp. 3-29.
2. Geologiya i poleznye iskopaemye Rossii, t. 5. Arkticheskie i dal'nevostochnye morya. Kn. 1. Arktich-eskie morya. Pod red. I.S. Gramberga, V.L. Ivanova, Yu.E. Pogrebitskogo (Geology and minerals of Russia, vol. 5. Arctic and Far East seas, book 1. Arctic seas. Gramberg I.S., Ivanov V.L., Pogrebitskiy Yu.E. (Eds.)), Saint-Petersburg, Izd-vo VSEGEI, 2004, 468 p.
3. Nikolaev V.V., Vrublevskii A.A., Akhmadulin V.A., Kuznetsov V.E. Geodinamika i seismicheskoe raion-irovanie materikovoi chasti Dalnego Vostoka (Geodynamic and seismic zoning of inland Far East), Vladivostok, DVO RAN, 2000, 90 p.
4. Morozov V.N., Kolesnikov l.Yu., Tatarinov V.N. Geoekologiya. 2011, no 6, pp. 524-542.
5. Morozov V.N., Tatarinov V.N. Geoekologiya. 1996, no 6. C.109-120.
6. Morozov V.N., Kolesnikov l.Yu., Belov S.V., Tatarinov V.N. Geoekologiya. 2008, no 3, pp. 232-243.
7. Rosman G.I., Bykhovskikh L.Z., Samsonov B.G. Khranenie i zakhoronenie radioaktivnykh otkhodov (Storage and disposal of nuclear waste), Moscow, VIMS, 2004, 240 p.
