Приложение 1.1. Рекомендации по реализации консервативного подхода при моделировании миграции радионуклидов в системе барьеров безопасности и вмещающей геологической среде в рамках оценки долговременной безопасности пунктов приповерхностного захоронения радиоактивных отходов

Приложение N 1.1
к руководству по безопасности
при использовании атомной энергии
"Оценка долговременной безопасности
пунктов приповерхностного захоронения
радиоактивных отходов", утвержденному
приказом Федеральной
службы по экологическому,
технологическому и атомному надзору
от 14 декабря 2016 г. N 531

Рекомендации
по реализации консервативного подхода при моделировании
миграции радионуклидов в системе барьеров безопасности
и вмещающей геологической среде в рамках оценки
долговременной безопасности пунктов приповерхностного
захоронения радиоактивных отходов

1. Вариантом реализации консервативного подхода в рамках оценки долговременной безопасности ППЗРО является использование одномерного камерного моделирования для расчета миграции радионуклидов в системе барьеров безопасности и вмещающей геологической среде.

2. Создание камерной модели, описывающей систему захоронения РАО.

2.1. Система захоронения РАО рассматривается как совокупность камер (блоков), между которыми происходит перенос радионуклидов (камерная модель).

2.2. Камерная модель системы захоронения РАО имеет блочную структуру, включающую:

область захоронения РАО;

инженерные барьеры безопасности ППЗРО;

элементы ближней зоны системы захоронения РАО;

элементы дальней зоны системы захоронения РАО;

элементы окружающей среды (биосфера).

Пример камерной модели, состоящей из трех блоков, приведен на рисунке 1.

2.3. Элементы ППЗРО, в которых не происходит перенос (миграция) радионуклидов, не учитываются в камерной модели.

2.4. Перенос радионуклидов между камерами описывается системой обыкновенных дифференциальных уравнений первого порядка:

(1),

где:

i и j - камеры;

N и M - количество (активность) радионуклидов N и M в камере (N - дочерний радионуклид в цепочке распада), Бк;

S(t) - внешний источник радионуклида N, Бк/год;

, - постоянные распада для радионуклидов N и M, год-1;

и - коэффициенты переноса, характеризующие процессы поступления и выхода радионуклида N из камер i и j, год-1.

Для камерной модели, приведенной на рисунке 1, система линейных дифференциальных уравнений имеет вид:

(2),

где:

S(t) характеризует поступление радионуклида в первую камеру модели; - постоянная радиоактивного распада, год-1.

|
|
|
\/
S(t) - источник РАО
1. Инженерный барьер ППЗРО
--------------->
2. Ненасыщенная зона
--------------->
3. Водоносный горизонт

Рисунок 1. Трехкамерная модель системы захоронения РАО

2.5. В камерной модели предполагается так называемое полное смешивание радионуклида в соответствующей камере (мгновенное достижение равномерной концентрации радионуклида в камере).

2.6. Камерная модель системы захоронения РАО создается с учетом рекомендаций по использованию упрощенной модели переноса радионуклидов из ППЗРО, приведенных в п. 82 настоящего Руководства по безопасности.

3. Определение коэффициентов переноса радионуклидов между камерами.

3.1. Вертикальный перенос радионуклидов в упаковках РАО, инженерных барьерах безопасности, ненасыщенной зоне.

Перенос радионуклидов в упаковках РАО, инженерных барьерах безопасности и ненасыщенной зоне описывается коэффициентом :

, (3),

где:

qin - скорость инфильтрации в среде (упаковках РАО/инженерных барьерах безопасности/ненасыщенной зоне), м/год;

- активная пористость среды;

d - глубина (толщина) среды, через которую происходит перенос радионуклидов, м.

Безразмерный коэффициент задержки R определяется следующим образом:

, (4),

где:

Kd - коэффициент межфазного распределения радионуклида в среде, м3/кг;

- плотность среды, кг/м3.

3.2. Конвективный перенос.

Конвективный поток (Бк/год) радионуклида из камеры i в камеру j определяется следующим образом:

, (5),

где:

vi - скорость конвективного потока, м/год;

- общая поверхность камер i и j, м2;

Ci - объемная активность (концентрация) радионуклида в поровой воде камеры i, Бк/м3.

Скорость конвективного потока vi определяется следующим образом:

, (6),

где:

ki - коэффициент фильтрации в камере i, м/год,

- градиент напора.

Коэффициент конвективного переноса между камерами i и j определяется следующим образом:

, (7),

где:

vi - скорость конвективного потока в камере i, м/год;

Ri - коэффициент задержки для материала камеры i;

Li - длина камеры i, м.

3.3. Диффузионный/дисперсионный перенос.

Диффузионный/дисперсионный поток радионуклидов из камеры i в камеру j определяется следующим образом:

, (8),

где:

Di - коэффициент дисперсии в камере i, м2/год, объединяющий процессы молекулярной диффузии в жидкости и дисперсии, возникающей при движении жидкости в поровой среде. Коэффициент дисперсии Di определяется следующим образом:

, (9),

где:

De - коэффициент молекулярной диффузии в жидкости, м2/год;

- полная пористость в камере i;

- продольная дисперсия, м;

- длина, которая задается как расстояние между центрами камер или меньшая из длин двух камер в направлении потока, м.

Диффузионный/дисперсионный поток моделируется посредством задания прямого и обратного потоков с коэффициентами переноса между камерами i и j:

, (10),

. (11)

4. Реализация консервативного подхода.

Из параметров, от которых зависят коэффициенты переноса (7), (10), (11), наибольшей неопределенностью характеризуются и Kd, значения которых для схожих условий могут изменяться в широких диапазонах. Значения этих параметров, с одной стороны, влияют на расчетную скорость переноса радионуклидов в рассматриваемых средах, а с другой - на размытие фронта движения радионуклидов и их концентрацию.

В задачах оценки долговременной безопасности пунктов приповерхностного захоронения радиоактивных отходов, в которых значения расчетных величин (например, концентраций радионуклидов в подземных водах) сравниваются с контрольными уровнями, рекомендуется проведение серии расчетов для различных допустимых значений исходных параметров, в том числе и Kd. В ходе таких расчетов для условий конкретных сценариев определяются результаты, соответствующие наиболее неблагоприятным последствиям, которые в итоге принимаются в качестве консервативной оценки.