Вследствие аварии на атомной станции «Фукусима» постоянно наблюдаются сильные выбросы радиации в океан посредством переноса радионуклидов грунтовыми водами. Для предотвращения распространения радионуклидов в грунтовых водах было принято решение оградить здания турбин и реакторов ледяной стеной многометровой глубины. Эта ледяная стена должна ограничить как приток грунтовых вод в зону реакторов, так и отток из нее грунтовых вод, загрязненных радионуклидами. С этой целью по периметру атомной станции, как показано ниже на рисунке, будут размещены вертикальные охлаждающие устройства.
Схема расположения охлаждающих элементов по периметру атомной станции «Фукусима»
Зеленым цветом показан водонепроницаемый слой скальных пород.
Желтым цветом показан слой насыпного грунта, в котором имеет место движение грунтовых вод.
R/B и T/B – здания реактора и турбин соответственно.
Синими стрелками указаны потоки воды.
Для обоснования принятых проектных решений необходимо выполнить компьютерное моделирование искусственной заморозки грунта на заданной площадке в специализированном программном обеспечении Frost 3D.
При выполнении компьютерного моделирования, исходя из геометрических размеров рассматриваемой области и геолого-литологического строения грунтов, были приняты следующие размеры расчетной области: длина — 450 метров, ширина — 210 метров, высота – 30 метров.
На основании информации, полученной по инженерно-геологическим скважинам, воспроизводилось геолого-литологическое строение грунта в рассматриваемой области (рисунок ниже).
При создании трехмерной геометрической модели принималось во внимание наличие фундаментов и оснований зданий реакторов и турбин.
Глубина залегания фундаментов зданий реакторов и турбин составляет порядка 10 метров. При выполнении теплотехнического расчета предполагалось, что температура данных фундаментов и оснований соответствует температуре внутри здания.
Расположение фундаментов зданий реакторов и турбин в программе Frost 3D
Для замораживания грунта по предварительным оценкам будет использовано около 30 км охлаждающих элементов. Эти элементы представляют собой набор вертикально установленных охлаждающих труб. Длина каждой трубы – 30 метров. Расстояние между трубами – 1 метр. Таким образом, будет установлено 1073 охлаждающие трубы по периметру атомной станции. Для их охлаждения предполагается использовать 14 холодильных машин мощностью по 400 кВт каждая.
Таким образом, исходя из геолого-литологического строения грунтов, расположения фундаментов и оснований реакторов, турбин и других сооружений, а также расположения вертикальных охлаждающих труб, в программе Frost 3D была создана трехмерная геометрическая модель атомной станции «Фукусима» (рисунки ниже).
Трёхмерная модель грунтов и охлаждающих устройств в программе Frost 3D
Трёхмерная геометрическая модель объектов инфраструктуры АЭС «Фукусима» в программе Frost 3D
При выполнении теплотехнического расчета принимались следующие теплофизические свойства грунтов:
Таблица 1 — Теплофизические свойства грунтов возле АЭС «Фукусима»
Теплопроводность грунта в талом состоянии [Вт/м*К] |
Теплопроводность грунта в мерзлом состоянии, [Вт/м*К] |
Теплоемкость грунта в талом состоянии, [МДж/м3*К] |
Теплоемкость грунта в мерзлом состоянии, [МДж/м3*К] |
|
---|---|---|---|---|
Насыпной грунт |
2.20 |
3.40 |
2.78 |
2.03 |
Скалистые породы |
2.00 |
2.16 |
2.40 |
1.95 |
Проведенные геологические изыскания показали, что в насыпном грунте фильтрация составляет порядка 0.1 м/день. В скалистых породах фильтрация на порядок меньше.
На поверхности расчетной области задавалось изменение температуры от времени. При этом использовалась среднемесячная максимальная температура с целью моделирования наиболее неблагоприятных для промерзания грунта климатических условий.
Задание зависимости температуры от времени в программе Frost 3D
На нижней и боковых границах расчетной области был задан нулевой тепловой поток. При этом эти границы области моделирования расположены достаточно далеко от интересующего участка, чтобы не оказывать влияние на расчет тепловых процессов для рассматриваемого периода времени.
Для моделирования тепловых полей трехмерная геометрия рассматриваемого участка дискретизировалась гексаэдрической сеткой, содержащей 17 828 087 узлов.
Гексаэдрическая расчётная сетка, созданная в программе Frost 3D
Моделируемый период времени составлял два года, начиная с даты 01.09.2014.
Время, затраченное программой Frost 3D на просчет этой модели на персональном компьютере с центральным процессором Intel Core i7 3770 – 3.4 ГГц и объемом оперативной памяти 16 ГБ составило три часа. Время расчета этой же модели в программе Frost 3D на графическом процессоре Nvidia Tesla K20c составило менее 6 минут.
Ниже представлены результаты моделирования распределения температуры в рассматриваемом участке для различных моментов времени.
01.12.2014 – Распределение температуры через три месяца работы охлаждающих устройств
Для лучшей наглядности результатов моделирования, рассмотрим распределение температуры на более детализированном участке моделируемой области. На рисунке ниже этот участок обозначен красным прямоугольником.
Участок расчётной области для детального анализа распределения тепловых полей и незамерзшей воды
Более отчетливо фронт промерзания можно наблюдать по распределению количества замерзшей воды в порах грунта.
Распределение количества замерзшей воды в сечении XZ на 10.12.2014. Синим цветом обозначен грунт с незамерзшей влагой.
Распределение количества замерзшей воды на сечениях через полгода показывают полную проморозку грунта в местах нахождения охлаждающих устройств.
Дальнейшая работа охлаждающих устройств в режиме максимальной мощности больше не является необходимой, и охлаждающие элементы можно переключить в режим поддержания температуры.
Кроме этого, как показывают нижеприведенные распределения количества замерзшей воды еще через полгода работы, дальнейшая заморозка практически не повлияет на изменение положения фронта промерзания:
31.08.2016 – Распределение температуры через 2 года работы охлаждающих устройств
Таким образом, современные компьютерные технологии позволяют прогнозировать распределение температур и количество незамерзшей влаги в грунте вокруг АЭС «Фукусима». Это позволит проверить правильность проектных решений при построении ледяной стены по периметру АЭС «Фукусима».