Приложение 1.1. Рекомендации по реализации консервативного подхода при моделировании миграции радионуклидов в системе барьеров безопасности и вмещающей геологической среде в рамках оценки долговременной безопасности пунктов приповерхностного захоронения радиоактивных отходов
Приложение N 1.1
к руководству по безопасности
при использовании атомной энергии
"Оценка долговременной безопасности
пунктов приповерхностного захоронения
радиоактивных отходов", утвержденному
приказом Федеральной
службы по экологическому,
технологическому и атомному надзору
от 14 декабря 2016 г. N 531
Рекомендации
по реализации консервативного подхода при моделировании
миграции радионуклидов в системе барьеров безопасности
и вмещающей геологической среде в рамках оценки
долговременной безопасности пунктов приповерхностного
захоронения радиоактивных отходов
1. Вариантом реализации консервативного подхода в рамках оценки долговременной безопасности ППЗРО является использование одномерного камерного моделирования для расчета миграции радионуклидов в системе барьеров безопасности и вмещающей геологической среде.
2. Создание камерной модели, описывающей систему захоронения РАО.
2.1. Система захоронения РАО рассматривается как совокупность камер (блоков), между которыми происходит перенос радионуклидов (камерная модель).
2.2. Камерная модель системы захоронения РАО имеет блочную структуру, включающую:
область захоронения РАО;
инженерные барьеры безопасности ППЗРО;
элементы ближней зоны системы захоронения РАО;
элементы дальней зоны системы захоронения РАО;
элементы окружающей среды (биосфера).
Пример камерной модели, состоящей из трех блоков, приведен на рисунке 1.
2.3. Элементы ППЗРО, в которых не происходит перенос (миграция) радионуклидов, не учитываются в камерной модели.
2.4. Перенос радионуклидов между камерами описывается системой обыкновенных дифференциальных уравнений первого порядка:
где:
i и j - камеры;
N и M - количество (активность) радионуклидов N и M в камере (N - дочерний радионуклид в цепочке распада), Бк;
S(t) - внешний источник радионуклида N, Бк/год;
Для камерной модели, приведенной на рисунке 1, система линейных дифференциальных уравнений имеет вид:
где:
S(t) характеризует поступление радионуклида в первую камеру модели;
|
|
|
\/
|
S(t) - источник РАО
|
||||
1. Инженерный барьер ППЗРО
|
--------------->
|
2. Ненасыщенная зона
|
--------------->
|
3. Водоносный горизонт
|
Рисунок 1. Трехкамерная модель системы захоронения РАО
2.5. В камерной модели предполагается так называемое полное смешивание радионуклида в соответствующей камере (мгновенное достижение равномерной концентрации радионуклида в камере).
2.6. Камерная модель системы захоронения РАО создается с учетом рекомендаций по использованию упрощенной модели переноса радионуклидов из ППЗРО, приведенных в п. 82 настоящего Руководства по безопасности.
3. Определение коэффициентов переноса радионуклидов между камерами.
3.1. Вертикальный перенос радионуклидов в упаковках РАО, инженерных барьерах безопасности, ненасыщенной зоне.
Перенос радионуклидов в упаковках РАО, инженерных барьерах безопасности и ненасыщенной зоне описывается коэффициентом
где:
qin - скорость инфильтрации в среде (упаковках РАО/инженерных барьерах безопасности/ненасыщенной зоне), м/год;
d - глубина (толщина) среды, через которую происходит перенос радионуклидов, м.
Безразмерный коэффициент задержки R определяется следующим образом:
где:
Kd - коэффициент межфазного распределения радионуклида в среде, м3/кг;
3.2. Конвективный перенос.
Конвективный поток
где:
vi - скорость конвективного потока, м/год;
Ci - объемная активность (концентрация) радионуклида в поровой воде камеры i, Бк/м3.
Скорость конвективного потока vi определяется следующим образом:
где:
ki - коэффициент фильтрации в камере i, м/год,
Коэффициент конвективного переноса между камерами i и j определяется следующим образом:
где:
vi - скорость конвективного потока в камере i, м/год;
Ri - коэффициент задержки для материала камеры i;
Li - длина камеры i, м.
3.3. Диффузионный/дисперсионный перенос.
Диффузионный/дисперсионный поток
где:
Di - коэффициент дисперсии в камере i, м2/год, объединяющий процессы молекулярной диффузии в жидкости и дисперсии, возникающей при движении жидкости в поровой среде. Коэффициент дисперсии Di определяется следующим образом:
где:
De - коэффициент молекулярной диффузии в жидкости, м2/год;
Диффузионный/дисперсионный поток моделируется посредством задания прямого и обратного потоков с коэффициентами переноса между камерами i и j:
4. Реализация консервативного подхода.
Из параметров, от которых зависят коэффициенты переноса (7), (10), (11), наибольшей неопределенностью характеризуются
В задачах оценки долговременной безопасности пунктов приповерхностного захоронения радиоактивных отходов, в которых значения расчетных величин (например, концентраций радионуклидов в подземных водах) сравниваются с контрольными уровнями, рекомендуется проведение серии расчетов для различных допустимых значений исходных параметров, в том числе