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

 

Расчет характеристик ядерного реактора

1. Исходные данные

Исходные данные из [1] приведены в таблице 1.

Таблица 1. Исходные данные.

НаименованиеПринятое обозначениеЧисленное значениеМощность, МВт Тепловая: Электрическая: W Q 3070 930Обогащение, %q3,95Загрузка, т U-метm(U-мет)67,5Выгорание, МВт-сут / т U-метУВ4,2?104Время выдержки ТВЭЛов, сут.Тохл1020Массовое число цепочкиА107Коэффициент воспроизводстваКв0,752Заряд ядра нуклидов семейства изобарZ35

Тип реактора - ВВЭР-1000 (на тепловых нейтронах).

Химическая формула топлива - UО2.

2. Расчетная часть

.1 Определение удельного выгорания

Поскольку удельное выгорание в МВт·сут/ т U-мет. уже задано, пересчет в другие размерности труда не представляет. Единственное, что надо учесть, это измененный, по сравнению с природным, изотопный состав исходного материала. Поэтому требуется уточнение значений атомной и молярной масс урана и топлива с применением изотопных масс (в углеродной шкале) и заданной степени обогащения топлива U - 235.

Таблица 2.Содержание изотопов урана в природном и обогащенном топливе.

ИзотопМолярная масса,г/моль*Природный уран,масс %*Уран, обогащенный до 3,95 масс % по U-235, масс %U-238238,050899,273996,01875U-235235,04390,72043,95000U-234234,04100,00570,03125* - [1]

Массовый процент обогащенного урана был найден следующим образом:

В процессе обогащения содержание U-235 увеличилось от 0,7204 масс % до 3,95 масс%, то есть в 5,5 раз. Содержание U-234, который сопровождает U-235 в технологических схемах, также увеличилось в 5,5 раз по сравнению с природным.

%U-234(в топливе) = [q /% U-235(природный)]·U-234(природный) = [3,95/0,7204]·0,0057 = =0,03125 масс%.

Содержание U-238 в топливе было найдено как разность:

%U-238(в топливе) = 100% - %U-234(в топливе) - %U-235(в топливе) =100-0,03125-3,95= = 96,01875 масс%.

Находим молярную массу U-мет по формуле:

М(U-мет) = М(U-238)·%U-238(в топливе)+ М(U-235)·%U-235(в топливе)+ М(U-234)·%U-234(в топливе);

где М(U-238), М(U-235), М(U-234) - молярные массы изотопов U [4]

М(U-мет)=[238,0508·0,9601875]+[235,0439·0,0395]+[234,0410·0,0003125]=

=237,9308 г/моль.

Молярная масса топлива равна:

М(UО2) = 269,9206 г/моль,

где М(О) = 15,9949 г/моль.

Масса топлива в пересчете на UО2:

m(UО2) = [m(U-мет)/М(U-мет)]·М(UО2),

m(UО2) =[67,5/237,9308] ·269,9206= 76,575 т.

Масса U-235 в загрузке топлива с учетом обогащения:

m(U-235) = m(U-мет)·q = 67,5·0,0395 = 2,6663 т.

2.1.1Определение удельного выгорания в МВт·сут/т топлива

Удельное выгорание (УВ) топлива меньше удельного выгорания U-мет во столько раз, во сколько масса топлива UО2 больше массы загрузки U-мет:

УВ(U-235) = [m(U-мет)/ m(UО2)]·УВ(U-мет) = [67,5/ 76,575]·4,2·104 =

= 3,7·104 МВт·сут/т топлива.

2.1.2Определение удельного выгорания в МВт·сут/т U-235

УВ(U-235) = УВ(U-мет)/q = 3,7·104 /0,0395 = 9,37·105 МВт·сут/т U-235.

2.1.3Определение удельного выгорания в кг U-235 сгоревшего/ кг U-235 загруженного

Можно рассчитать как отношение удельного выгорания U-235 к предельному выгоранию:

УВ = УВ(U-235)/ПВ(U-235) = 9,37·105 /9,4·105 =

=0,9968 кг U-235 сгоревшего/ кг U-235 загруженного;

где ПВ=9,4·105 МВт·сут/т U-235 - предельная величина выгорания.

.1.4 Определение степени выгорания

Степень выгорания - это отношение реального удельного выгорания U-235 к предельному выгоранию:

СВ = 9,37·105 ·100/9,4·105 =99,68%.

.1.5 Определение fifa

Fifa (fissions per initial fussionablе atom) - число актов деления ядер на один первоначально загруженный делящийся атом.

Первоначально загружено ядер U-235:

N(U-235) = [m(U-235)·NA] /M(U-235) = [2,6663·106·6,02·1023] /235,0439 = 6,829·1027 ядер,

где NA=6,02·1023 моль -1 -постоянная Авогадро;

В реакторе при заданном удельном выгорании U-мет делится ядер:

,2·104·67,5·86400·3·1016 = 7,35·1027 ядер,

где 86400 с = 1 сут, 3·1016 дел/с = 1МВт

fifa = 7,35·1027 /6,829·1027 = 1,076

.1.6Определение fima

Fima (fissions per initial metal atom) - число актов деления ядер на один первоначально загруженный атом тяжелого металла.

Fima меньше fifa в n раз, где:

n = (N35 + N38 )/ N35 >> 1= (0,9601875+0,0395) / 0,0395 = 25,31 раз= fifa / n = 1,076 / 25,31 = 0,0425.

2.1.7Определение глубины выгорания

Глубина выгорания - доля выгоревших атомов топлива или тяжелого металла.

ГВ = fifa (1 + sЗ/ sД)35.

где sЗ = 97,4 - сечение захвата U-235;

?sД = 583,5 - сечение деления U-235 [1]

ГВ = 1,076·(1+97,4 / 583,5) = 0,92208 или 92,21%.

.2 Определение скорости выгорания

Скорость выгорания - это удельная мощность, т.е. мощность, отнесенная к единице массы топлива, металла или U-235.

2.2.1 Определение скорости выгорания в МВт / т топлива

V = W / m(UО2) = 3070/ 76,575= 40,09 МВт / т топлива.

.2.2 Определение скорости выгорания в МВт / т U-мет

V = W / m(U-мет) = 3070 / 67,5 = 45,48 МВт / т U-мет.

.2.3 Определение скорости выгорания в МВт / т U-235

V = W / m(U-235) = 3070 / 2.6663 = 1151,408 МВт / т U-235.

.3 Определение времени кампании

УВU-мет = (W·tK) / mU-мет ;

Следовательно, tK = УВ U-мет ·m U-мет / W,

tK = 4.2·104 ·67,5 / 3070 = 923 сут.

.4 Определение суммарной b?- активности продуктов деления с учетом и без учета короткоживущих радионуклидов

Для реакторов на тепловых нейтронах и уранового топлива (любой степени обогащения) набор продуктов деления и их независимые выходы являются постоянными. Это дает возможность табулировать значения суммарной b-активности продуктов деления как функцию времени кампании при неизменной (как правило - единичной) удельной мощности реактора для неизменной (тоже, как правило - единичной) массы топлива любой степени обогащения.

Таблицы для определения b-активности часто имеют достаточно большой шаг по времени кампании и полученное время кампании может оказаться лежащим между двумя табулированными величинами. В этом случае необходимая величина активности может быть найдена путем графической линейной интерполяции. Если время кампании выходит за пределы табулированных величин, то возможно два варианта. Первый - поискать справочник, в котором будет интересующее значение времени кампании. Второй - использовать графическую нелинейную экстраполяцию, имея при этом в виду, что чем больше интервал времени, для которого проводится экстраполяция, тем больше шансов получить заметную погрешность в конечном результате. Далее полученная величина активности умножается на удельную мощность, на массу топлива и окончательный результат выдается в кюри и в беккерелях[1].

2.4.1Определение суммарной b?- активности продуктов деления без учета короткоживущих радионуклидов

2.4.1.1К концу облучения (Тохл=0)

Справочные данные для заданного реактора из [2] приведены в таблице 3.

Таблица 3. Полная активность смеси продуктов деления 235U, Ки / кг в зависимости от времени кампании tк и времени выдержки Tохл для единичной мощности реактора.

Тохл, сутtK, сут130180200360720?06737077177638599207207,479,6810,617,027,416510804,355,686,2310,317,6148

Рисунок 1. Зависимость активности продуктов деления от времени кампании.

По графику находим, что для tK = 923 сут. активность для единичной мощности реактора равняется 880,0 (Ku/кг) / (Вт/г).

Тогда удельная активность для полной мощности реактора:

AWуд = A0уд·V,

где V - удельная мощность (скорость выгорания), МВт / т U-мет.

AWуд = 880 (Ku/кг) / (Вт/г)*45,48 кВт / кг U-мет = 4,0*104 Ku / кг

Полная активность на всю загрузку, Ku:

A? = AWуд·m(U-мет)кг = 4,0*104 ·67,5·1000= 2,70·109 Ku.

Полная активность на всю загрузку, Бк:

AS = AS?Ku) ·3,7·1010 = 2,70·109 ·3,7·1010 ? 10,0·1019 Бк,

где 1Ku = 3,7·1010 Бк.

Итоговые результаты расчетов b?- активности продуктов деления для Тохл = 0 без учета короткоживущих радионуклидов представлены в таблице 4.

Таблица 4. Удельная и полная b?- активность ПД на конец кампании без учета короткоживущих радионуклидов.

ПроцессA0уд , (Ku/кг) / (Вт/г)АWуд, Ku / кгА??, KuА??, БкРабота реактора - облучение ТВЭЛов880,04,0*1042,70·10910,0·1019

2.4.1.2К концу выдержки (Тохл = 1020 суток)

По данным таблицы 3 строим зависимость активности ПД от времени охлаждения и находим искомую величину активности для единичной мощности реактора при tK = 923 сут, Тохл = 1020 дней.

Рисунок 2. Зависимость активности продуктов деления от времени охлаждения.

Из уравнения у=-0,0272х+47 следует, что при Тохл = 1020 суток, tK =720 дней, активность равна 19,233 (Ku/кг) / (Вт/г).

А для tK = ? активность равняется -0,0472*1020+199=150,836 (Ku/кг) / (Вт/г).

По полученным данным можно построить кривую, которая будет характеризовать изменение активности ПД во время выдержки.

Рисунок 3. Зависимость b?- активности долгоживущих продуктов деления от времени компании.

По графику находим, что для tK = 923 сут., величина активности ПД для единичной мощности реактора на конец времени выдержки будет равняться 120(Ku/кг) / (Вт/г).

Результаты расчетов в ?- активности продуктов деления без учета короткоживущих радионуклидов при Тохл = 1020 сут приведены в таблице 5.

Таблица 5.Удельная и полная в?- активность на конец выдержки.

ПроцессА°уд, (Ku/кг) / (Вт/г)АWуд,Ku / кгА??,KuА??,БкВыдержка ТВЭЛов1205,46·1033,684·10813,63·1018

2.4.2 Определение суммарной в- активности продуктов деления с учетом короткоживущих радионуклидов

.4.2.1 К концу облучения (Тохл = 0).

Справочные данные для заданного реактора из [2] приведены в таблице 6.

Таблица 6.Удельная в- активность (в кюри) смеси продуктов деления U-235 на 1 кг металла в зависимости от времени кампании tK и времени выдержки Тохл при мощности реактора W = 1 Вт/ U-мет.

Тохл, сут.tK, сут.720?02307236872027,4165108017,6148

По данным таблицы 6, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 сут, Тохл= 0 сут (рис.4):

Рисунок 4. Зависимость в- активности смеси продуктов деления от времени кампании.

По данным таблицы 6, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 дня, Тохл= 0:

А°уд = 2348 Ku·г/кг·Вт.

Результаты расчетов b?- активности продуктов деления при Тохл = 0 представлены в таблице 7.

Таблица 7. Удельная и полная в?- активности продуктов деления на конец кампании.

ПроцессА°уд, (Ku/кг) / (Вт/г)АWуд, Ku / кгА??, KuА??, БкРабота реактора - облучение ТВЭЛов2,35·1031,07·1057,22·1092,67·1020

2.4.2.2К концу выдержки (Тохл = 1020 сут).

Значения b?- активности ТВЭЛов на конец выдержки с учетом короткоживущих радионуклидов совпадают с соответствующими значениями без их учета, т.к. за время охлаждения Тохл = 1020 сут. короткоживущие радионуклиды распадаются и не вносят вклад в значение активности.

.5 Определение суммарного г - эквивалента продуктов деления

Для оценки g - излучения продуктов деления воспользуемся величинами, характеризующими воздействие этого излучения, выраженными в г-экв Ra.

2.5.1Определение суммарного г - эквивалента продуктов деления без учета короткоживущих радионуклидов

2.5.1.1К концу облучения (Тохл = 0)

Справочные данные для заданного реактора из [2] приведены в таблице

Таблица 8.Удельный г-эквивалент долгоживущих продуктов деления 235U (в г-экв Ra) на 1 кг металла в зависимости от времени кампании tK и времени выдержки Тохл при мощности реактора W = 1 Вт/г 235U.

Тохл, сут.tK, сут.720?02402517201,382310801,0822,2

По данным таблицы 8, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 дня, Тохл= 0 (рис.5): А°уд = 242 г-экв Ra / кг.

Рис. 5 Зависимость г - эквивалента долгоживущих продуктов деления от времени кампании.

Результаты расчета г - эквивалента продуктов деления без учета короткоживущих радионуклидов для Тохл = 0 сут приведены в таблице 9.

Таблица 9. Удельный и полный г- эквивалент на конец кампании.

ПроцессА°уд,г-экв Ra / кгАWуд, г-экв Ra / кгА??, KuРабота реактора - облучение ТВЭЛов2421,10·1047,43·108

К концу выдержки (Тохл = 1020 сут)

По данным таблицы 8, аналогично пункту 2.4.1.2., находим искомую величину активности для единичной мощности реактора, tK = 923 сут,Тохл= 1020 сут (рис.6 и 7)

Рис.6. Зависимость г-эквивалента долгоживущих продуктов деления от времени выдержки.

Из уравнения у=-0,02722х+47 следует, что при Тохл = 1020 суток, tK =720 дней, активность равна 19,233 (г-экв Ra/кг).

А для tK = ? активность равняется -0,04722*1020+199=150,836 (г-экв Ra/кг).

Рис.7. Зависимость г-эквивалента долгоживущих продуктов деления от времени компании.

Результаты расчета г- эквивалента продуктов деления без учета короткоживущих радионуклидов для Тохл = 1020 сут приведены в таблице 10.

Таблица 10. Удельный и полный г- эквивалент на конец выдержки.

ПроцессА°уд, г -экв Ra /кг АWуд, г -экв Ra/ кг·А??,г -экв RaВыдержка ТВЭЛов38,001,73·1031,17·105

2.5.1Определение суммарного г - эквивалента продуктов деления с учетом короткоживущих радионуклидов

2.5.1.1К концу облучения (Тохл = 0).

Справочные данные для расчета заданного реактора из [2] приведены в таблице 11.

Таблица 11. Удельный г -эквивалент смеси продуктов деления 235U (в г-экв Ra) на 1 кг металла в зависимости от времени кампании tK и времени выдержки Тохл при мощности реактора W = 1 Вт/г 235U с учетом короткоживущих изотопов.

Тохл, сут.tK, сут.720?06326437201,3823,010801,0822,2

По данным таблицы 11, аналогично пункту 2.4.1.1. находим искомую величину активности для единичной мощности реактора, tK= 923 сут, Тохл= 0 сут (рис.8):

Рис. 8 Зависимость г-эквивалент смеси продуктов деления от времени кампании.

А°уд = 635 (г-экв Ra/кг)/(Вт/г).

Результаты расчета г - эквивалента продуктов деления с учетом короткоживущих радионуклидов для Тохл = 0 приведены в таблице 12.

Таблица 12. Удельный и полный г- эквивалент на конец кампании.

ПроцессА°уд,г -экв Ra / кг АWуд,г -экв Ra/ кг·А??,г -экв RaРабота реактора - облучение ТВЭЛов6352,89·10419,5·108

2.4.2.2К концу выдержки (Тохл = 1200 сут).

Значения ?г-эквивалента ТВЭЛов на конец выдержки с учетом короткоживущих радионуклидов совпадают с соответствующими значениями без их учета, т.к. за время охлаждения Тохл = 1020 сут короткоживущие радионуклиды распадаются и не вносят вклад в значение активности.

2.5Определение плотности потока нейтронов и количества образовавшихся Pu-239, 240, 241 с учетом всех путей их расходования

Расход плутония идёт в трех направлениях: радиоактивный распад, деление, захват. Каждое из направлений характеризуется своей константой, а следовательно характеризуется количественно. Общее количество наработанного плутония за вычетом количества, дожившего до конца кампании, дает общее число израсходованного плутония, которое в свою очередь распределяется пропорционально сечениям деления и захвата между реакциями деления и захвата. Плутоний, претерпевший деление, вносит свой вклад в создании мощности реактора, который не учитывался в расчете пункта 2.1, поэтому реальное количество выгоревшего U-235 будет ниже заявленного - его необходимо пересчитать.

Плутоний претерпевший захват нейтрона, вносит свой вклад в поставку более тяжелого изотопа плутония.

2.5.1Определение плотности потока нейтронов

рср = 2,2·1010 ·W / m(U-235) [1]

рср - средняя плотность потока нейтронов.

рср = 2,2·1010 ·3070/ 2,6663 = 2,533·1013 см-2·с-1.

.5.2Определение количества образовавшихся Pu-239, 240, 241 с учетом всех путей их расходования

Цель этого пункта - дать оценку вклада образующегося плутония-239 в работу реактора и определить количество плутония-239 и его более тяжелых изотопов, накопившихся к концу кампании. Это позволит исправить степень выгорания U-235, полученную в п.2.1 по чисто формальному признаку (через энерговыделение), и получить более объективную оценку глубины выгорания.[1]

Рассмотрим цепочку превращений Pu в реакторе:

Эволюция Pu-239 в реакторе на тепловых нейтронах зависит от мощности W, потока нейтронов р, времени кампании tK.

Пусть, W = const, то есть скорость поставки Pu-239 Р239 = const. При этом, р = const, если загрузка U-235 постоянна.

Данные для расчета эволюции Pu представлены в таблице 13.

Таблица 13.Данные для расчета эволюции плутония-239.

РадионуклидТ1/2,летl, с-1sд, барнsз, барнL*з,с-1L*д,с-1L,с-1Т°,с123456789Pu-2392,4·1049,16·10-13 744267,26,49·10-91,81·10-82,46·10-82,82·107Pu-2406,54·1033,36·10-120,032907,04·10-97,29·10-137,05·10-99,83·107Pu-2411,45·1011,52·10-910093688,94·10-92,45·10-83,50·10-81,98·107

, 5 - 1 барн = 10-24 см; [1]

6, 7 - L* = s·pср;

8 - L?= l + L*з +L*д ;

- Т°= ln2 / L.

Скорость деления ядер:

F = 3,1·1016·W = 3,1·1016·3070 = 9,21·1019 дел/с

Скорость поставки Pu:

P239 = F·KB·(1+?sз?/sд) = 9,21·1019·0,752·(1+ 97,4/583,5) = 8,082·1019 дел/с.

Общее количество наработанного за время кампании Pu-239:

N239 = P239 ·tK = 8,082·1019·923·86400 = 6,445·1027 ядер

m (Pu-239) = N239 ·M(Pu-239)/NA = 6,445·1027·239,0522/6,02·1023 = 255,94*104г =2559,4 кг ? 2,6 т.

где М (Pu-239) = 239,0522 г/моль - молярная масса Pu-239[1].

2.5.2.1Определение количество образовавшегося Pu-239

Количество Pu-239, дожившего до конца кампании, определим по формуле, тождественной случаю накопления дочернего радионуклида при вековом равновесии.

Вместо постоянной активности материнского нуклида используем постоянную скорость поставки Pu-239 при работе реактора с постоянной мощностью:

2,823*1027 шт.

1121 кг.

Общее количество израсходованного Pu-239:

m (Pu-239) - m (Pu-239,tK) = 2559,4 - 1121 = 1438,4 кг.

Рu-39 в реакторе расходуется по трем путям - реакция радиационного захвата, реакция деления и собственный радиоактивный распад. Последний путь расходования при работающем реакторе не вносит значительной вклад, т.к. период полураспада Pu-39 (T1/2 = 24 тыс. лет) значительно превосходит время кампании реактора (tк = 923 сут.), следовательно, определим количества Pu-39, израсходованные только на деление и захват, разделив полученное количество израсходованного Рu-39 через отношение сечений деления и захвата:

т.о. израсходовано на деление Pu-239 - 1057,87 кг

израсходовано на захват Pu-239(на образование Pu-240) - 380,53 кг.

2.5.2.2 Определение количество образовавшегося Pu-240

Количество Pu-240, дожившего до конца кампании:

N240,t =7,5899*1026 шт.

= 302,65 кг.

где М (Pu-240) = 240,0538 г/моль - молярная масса Pu-240.

Общее количество израсходованного Pu-240:

m (Pu-240) - m (Pu-240,tK) = 380,53 - 302,65 = 77,88 кг.

т.о. израсходовано на деление Pu-240 - 8·10-3 кг

израсходовано на захват Pu-240 - 77,872 кг;

то есть, почти весь Pu-240 расходуется на образование Pu-241.

Из полученных значений видно, что практически весь Рu - 40 расходуется на образование Рu - 41 по реакции захвата, т.к. его сечение захвата значительно превосходит сечение деления.

.5.2.3 Определение количество образовавшегося Pu-241

Количество Pu-241, дожившего до конца кампании:

N241(tK) = 8,603·1025 шт.

= 34,45 кг.

где М (Pu-241) = 241,0569 г/моль - молярная масса Pu-240.

Общее количество израсходованного Pu-241:

m (Pu-241) - m (Pu-241,tK) = 77,872 - 34,45 = 43,43 кг.

т.о. израсходовано на деление Pu-241 - 31,81 кг.

израсходовано на захват Pu-241 - 11,61 кг.

Результаты расчета эволюции изотопов Pu приведены в таблице 14.

Таблица 14. Эволюция изотопов плутония во время кампании

РадионуклидНаработано за время кампании, кгИзрасходовано на деление, кгИзрасходовано на захват, кгДожило до конца кампании, кгPu-2392559,401057,87380,531121,00Pu-240380,538·10-377,872302,65Pu-24177,8731,8111,6134,45

2.5.3 Определение количества реально сгоревшего U-235

С учетом полученных данных уточним количество реально сгоревшего U-35.

В реакторе при заданном удельном выгорании U-мет делится ядер:

,2·104·67,5·86400·3,1·1016=7,593·1027 ядер,

где 86400 с = 1 сут, 3,1·1016 = 1МВт,

Тогда массу сгоревшего U-235 равна:

m(U-235) = [N(U-235) · M(U-235)] / NA = [7,593·1027·235,0439] /6,02·1023 = 2964,6 кг,

где NA=6,02·1023 моль -1 -постоянная Авогадро;

M(U-235)=235,0439 г/моль - молярная масса U-235.

Учитывая, что часть Pu-239, расходуемая на деление, идет на замену U-235, получим, что реально сгорает U-235:

=

= 2964,6 - 1057,87 - 0,008 - 31,81 = 1874,912 кг.

кг U-235 / кг топл.

Тогда степень выгорания U-235:

СВреал. = mU-35,сгор/ mU-35,загр. = 1874,912/2666,3 = 70,32%,

где ПВ=9,4·105МВт·сут/т U-235 - предельная величина выгорания.

Уточним значения fifa, fima, глубины выгорания:

в реакторе делится ядер U-235:

N = 1874,912 ·1000·6,02·1023 / 235,0439 = 4,802·1027 шт.

fifa = 4,802·1027 / 6,829·1027 = 0,703

где 6,829·1027 шт.- число, первоначально загруженных ядер U-235.

fima = 0,703/25,31 = 0,0278

ГВ = 0,703·100%·(1+97,4 / 583,5) = 82,04%.

2.6Определение весового содержания продуктов деления (ПД)

Примем, что на 1 кг сгоревшего U-235 образуется 1 кг ПД. Тогда масса ПД с учетом того, что на 1 кг загруженного U-235 сгорает 703г U-235 (из ftfa) и что масса загрузки U-235 составляет 2666,3 кг, равна 1874,41 г.

24,48 гПД/кгUО2.

43,92 кг ПД / т U-мет.

.7 Определение коэффициента очистки топлива

В общем случае коэффициенты очистки и времена выдержки облученного топлива зависят от многих факторов, таких как тип реактора, тип топливного цикла, дальнейшее использование топлива, степень выгорания и ряда других. Тем не менее есть общие "моменты": выдержка облученного топлива желательна для полного завершения эволюции цепочки U-238 (n,g, 2b)... или Th-232(n,g, 2b)...;распада U-237 в урановом топливе; возможно более заметного снижения удельной активности облученного материала с целью избежать в последующем заметных радиационно химических осложнений при переработке топлива. Очистка дает возможность либо получить новое ядерное топливо для повторного использования, либо направить его (топливо) на другие цели (оборона, использование в народном хозяйстве). Реактор на тепловых нейтронах, гетерогенный, на природном (или слабообогащенном) топливе, с невысокой плотностью потока нейтронов и умеренным временем кампании (т.е. с малой степенью выгорания). Время выдержки лимитируется главным образом необходимостью исключить в будущем работу с ураном-237. Это ведет к временам выдержки порядка 80-100 дней. Коэффициент очистки определяется из основной посылки: активность топлива после радиохимической переработки не должна превышать активности "первичного" урана (на один a?- 2b распада ближайших дочерних, что дает удельную активность по b-2,44·10+4 Бк/г).[1] При активности облученного ядерного топлива (Тохл =1020 сут), равной АWуд=5,46 Ku/г U-мет (п. 2.4.1.2.) коэффициент очистки равен:оч = 5,46·3,7·1010 / 2,44·104 = 8,28·106.

Определение суммарной концентрации ПД после растворения ТВЭЛов при соотношении Т:Ж = 1:10. Ранее найдено, что на 1 кг UО2 образуется 24,48 г ПД (п. 2.7.1.). После растворения ТВЭЛов концентрация ПД составит:

кг ПД/ т жидк.фазы;

где т.

Определение изменения активности для временных в абсолютных единицах (Бк) для цепочки с заданным А.

Цепочку с А = 107 выписываем из [1]:

Таблица 15. Характеристики радионуклидов цепочки с А = 107.

Радионуклидhнез, долиТ1/2l, с-1Мо1071,07E-033,5c1,98E-01Tc1071,82E-0429,0с2,39E-02Ru1071,69E-044,2мин2,75E-03Rh1071,10E-0621,7мин5,32E-04mPd1074,17E-1022,0c3,15E-02Pd1071,19E-106,50E+06лет3,38E-15mAg107044,3с1,56E-02

пояснение к таблице:

з - независимый выход, доли;

T1/2 - период полураспада, с;

л - постоянная распада, с-1;

Примечания: ветви, для которых приход или уход в соседние цепочки с большим(меньшим) А составляет меньше 5% и для которых л больше Л*з не учитываются.

1.106Ru: Л*з = sз·pср = 2,533·1013 ·0,15·10-24 = 3,8·10-12 с-1, так как для Ru л= =2,75·10-3 больше Л*з, то реакцию захвата для него не учитываем.

2.107Pd: Л*з = sз·pср = 2,533·1013 ·10,1·10-24 = 2,58·10-10 с-1, так как для Pd л меньше Л*з, то реакцию захвата учитываем и необходимо рассчитать обобщенную константу распада: Л = л + Л*з = 2,58·10-10 с-1, т.к. она нам пригодится для дальнейших расчетов.

3.10646Pdстаб: рассчитаем Л*з для стабильного изотопа палладия из более легкой цепочки: Л*з = sз·pср = 2,533·1013 ·0,24·10-24 = 6,08·10-12 с-1, она нам также пригодится в дальнейших расчетах.

За время кампании (tK = 923 сут)

Накопление Mo

Родоначальник цепочки - Mo является короткоживущим. Его активность уже через 35 секунд достигнет своей предельной величины, равной скорости его поставки, и после этого меняться не будет.

АMo max = PMo = F·hMo = 9,21·1019 ·1,07·10-3 = 9,86·1016 Бк

lg АMo max = 16,994;

Расчетная формула:

расчетный интервал времени 1-35с

Результаты расчета приведены в таблице 16.

Таблица 16.Эволюция Mo за время кампании.

t,секN,штA,БкlgtlgA18,94E+161,77E+160,00016,248104,29E+178,49E+161,00016,929154,72E+179,35E+161,17616,971204,88E+179,67E+161,30116,985254,94E+179,78E+161,39816,991304,94E+179,83E+161,47716,992354,94E+179,85E+161,54416,993797472004,94E+179,85E+167,90216,994

пояснение к таблице:- время от начала кампании, с;- число ядер данного изотопа;

A - активность данного изотопа, А=N*л, Бк;

(примечание: последнее значение времени - время кампании);

Накопление Tc

Так как периоды полураспада молибдена и технеция отличаются не сильно, то накопление технеция рассчитываем по двучленной формуле.

Расчетная формула:

Расчетный интервал времени, учитывая что Т1/2(Tc) = 29 с выбираем от 36с до 290 c. К этому времени скорость поставки технеция сравнится со скоростью его распада и в дальнейшем меняться не будет. Этот интервал разбиваем на 10 точек.

Таблица 17. Изменение активности Tc за время кампании.

t,секN,штA,БкlgtlgA362,55E+186,08E+161,55616,78563,41E+188,15E+161,74816,91763,95E+189,44E+161,88116,97964,28E+181,02E+171,98217,011264,56E+181,09E+172,10017,041464,66E+181,11E+172,16417,051864,76E+181,14E+172,27017,062004,78E+181,14E+172,30117,062464,81E+181,15E+172,39117,062904,82E+181,15E+172,46217,06797472004,82E+181,15E+177,90217,06

Накопление Ru

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

Расчетная формула:

Расчетный интервал времени - от 300 с до 2500 с.

Результаты расчета приведены в таблице 18.

Таблица 18.Эволюция Ru за время кампании.

t,секN,штA,БкlgtlgA3002,44E+196,70E+162,47716,836003,74E+191,03E+172,77817,018004,17E+191,15E+172,90317,0610004,42E+191,22E+173,00017,0813004,61E+191,27E+173,11417,1015004,67E+191,29E+173,17617,1118004,72E+191,30E+173,25517,1120004,74E+191,30E+173,30117,1122004,75E+191,31E+173,34217,1224004,75E+191,31E+173,38017,1225004,75E+191,31E+173,39817,12797472004,76E+191,31E+177,90217,12

Накопление Rh

Расчетная формула:

Расчетный интервал времени - от 2550 с до 13020 с (217 мин).

Результаты расчета приведены в таблице 19.

Таблица 19.Эволюция Rh за время кампании.

t,секN,штA,БкlgtlgA25501,68E+208,92E+163,40716,9530001,84E+209,81E+163,47716,9935001,99E+201,06E+173,54417,0245002,18E+201,16E+173,65317,0760002,34E+201,24E+173,77817,0980002,42E+201,29E+173,90317,1190002,44E+201,30E+173,95417,11100002,45E+201,30E+174,00017,11115002,46E+201,31E+174,06117,12130202,46E+201,31E+174,11517,12797472002,46E+201,31E+177,90217,12

Накопление Pdm

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

Расчётная формула:

Расчетный интервал времени - от 1 с до 220 с.

Результаты расчета приведены в таблице 20.

Таблица 20.Эволюция Pdm за время кампании.

t,секN,штA,БкlgtlgA13,78E+101,19E+090,0009,08205,70E+111,80E+101,30110,25551,00E+123,16E+101,74010,50751,10E+123,48E+101,87510,541001,17E+123,68E+102,00010,571251,20E+123,77E+102,09710,581501,21E+123,81E+102,17610,581751,21E+123,83E+102,24310,582001,22E+123,83E+102,30110,582201,22E+123,84E+102,34210,58797472001,22E+123,84E+107,90210,58

Накопление Pd

Палладий накапливается за счет своей цепочки и за счет прихода из более легкой цепочки в результате реакции радиационного захвата стабильным 106-м палладием. Поэтому накопление палладия рассчитываем по двум цепочкам: по своей и по 106-й. По своей цепочке рассчитаем его по одночленной формуле, просуммировав все независимые выходы предшественников, а так как палладий расходуется по двум направлениям: распад и захват, то используем обобщенную константу распада:

Накопление палладия за счёт прихода из 106-й цепочки рассчитаем следующим образом: мы должны учесть реакцию радиационного захвата 106-м стабильным палладием, а т.к. его накопление определяется в основном эволюцией долгоживущего 106Ru, следовательно и рутений мы должны учесть тоже. Ниобий, молибден и технеций являются короткоживущими, поэтому мы приплюсуем их независимые выходы к рутению 106-му. Родием 106-м мы можем пренебречь в силу его незначительного независимого выхода и малого периода полураспада (29,9с), он получается связующим звеном между 106-м рутением и стабильным 106-м палладием, и практически не оказывает никакого влияния на эволюцию цепочки. Накопление палладия 107-го рассчитаем по трёхчленной формуле:

где

ЛPd-107 - обобщенная константа распада Pd-107.

Л*з,Pd-106 - модифицированная константа радиационного захвата Pd-106.

Расчетная формула:

Расчетный интервал времени от 4ч до 11ч, т.к. Pd долгоживущий, то ограничимся 4-мя точками.

Результаты расчета приведены в таблице 21.

Таблица 21. Накопление Pd за время компании.

t,секN,штA,БкlgtlgA130201,71E+215,78E+064,126,765000006,55E+222,21E+085,708,3510000001,31E+234,43E+086,008,65100000001,31E+244,43E+097,009,65797472001,06E+253,58E+107,9010,55

Накопление Agm

Т.к. у метастабильного серебра нет собственного независимого выхода и единственный путь его поставки - палладий 107-й, то сопоставив их периоды полураспада можно с уверенностью говорить о вековом равновесии, которое наступит уже через 443с. Исходя из этого, расчеты проводить не целесообразно, т.к. повторятся расчёты материнского радионуклида.

За время выдержки (Тохл = 1020 сут).

К моменту выгрузки радионуклиды цепочки А = 107 будут иметь характеристики, представленные в таблице 22.

Таблица 22. Исходные данные для расчета эволюции цепочки с массовым числом А = 107 при выдержке ТВЭЛов.

Радионуклид?, с-1Начальное число ядер,N0Мо1071,98E-014,94E+17Tc1072,39E-024,82E+18Ru1072,75E-034,76E+19Rh1075,32E-042,46E+20mPd1073,15E-021,22E+12Pd1073,38E-151,06E+25

Эволюцию цепочки во время выдержки рассчитаем по уравнению Бейтмана, в соответствии с которым для цепочки вида:

¾® N2 ¾® N3 ?¾®???? Nj ?¾®???? Ni ,

Если принять, что при t =0 N1 = N01, a N2 = N3 = ...Nj...= Ni = 0, то справедливо

Распад Mo

Mo - короткоживущий радионуклид. Его распад считаем по простой формуле:

Расчетный интервал времени от 1с до 125с.

Результаты расчета приведены в таблице 23.

Таблица 23. Распад Mo за время выдержки.

t,секN,штA,БкlgtlgA14,05E+178,02E+16016,90106,82E+161,35E+16116,13209,42E+151,86E+151,3015,27401,80E+143,55E+131,6013,55603,42E+126,78E+111,7811,83806,52E+101,29E+101,9010,111001,24E+092,46E+082,008,391258,81E+061,74E+062,106,24

Распад Tc

Т.к. периоды полураспада Мо и Тс соизмеримы,то распад Тс считаем по двучленной формуле:

Расчетный интервал времени от 1 с до 1000 с.

Результаты расчета приведены в таблице 24.

Таблица 24. Распад Тс за время выдержки.

t,секN,штA,БкlgtlgA14,79E+181,14E+17017,057501,63E+183,90E+161,7016,591758,96E+172,14E+161,8816,3301004,93E+171,18E+162,0016,0721501,49E+173,56E+152,1815,5512004,52E+161,08E+152,3015,0335003,48E+138,31E+112,7011,92010002,25E+085,37E+063,006,730

Распад Ru

Т.к. периоды полураспада Тс и Ru соизмеримы, то распад Ru считаем по двучленной формуле:

Расчетный интервал времени от 1 с до 2520 с.

Результаты расчета приведены в таблице 25.

Таблица 25. Распад Ru за время выдержки.

t,секN,штA,БкlgtlgA14,76E+191,309E+17017,1173002,35E+196,463E+162,4816,8106001,03E+192,833E+162,7816,4529004,51E+181,240E+162,9516,09412001,98E+185,445E+153,0815,73615008,66E+172,382E+153,1815,37718003,80E+171,045E+153,2615,01921001,66E+174,565E+143,3214,65925005,54E+161,524E+143,4014,18330001,40E+163,850E+133,4813,58590009,56E+082,629E+063,956,420

Распад Rh

Расчетная формула:

Расчетный интервал времени от 1 с до 45000 с.

Результаты расчета приведены в таблице 26.

Таблица 26. Распад Rh за время выдержки.

t,секN,штA,БкlgtlgA12,46E+201,31E+17017,1220001,07E+205,71E+163,3016,7630006,31E+193,36E+163,4816,5350002,18E+191,16E+163,7016,0670007,52E+184,00E+153,8515,60100001,52E+188,11E+144,0014,91110008,96E+174,76E+144,0414,68120005,26E+172,80E+144,0814,45130203,06E+171,63E+144,1214,21450001,25E+106,65E+064,656,82

Распад Pdm

Расчетная формула:

Расчетный интервал времени от 1 с до 220 с.

Результаты расчета приведены в таблице 27.

Таблица 27.Эволюция Pdm за время выдержки.

t,секN,штA,БкlgtlgA11,18E+123,72E+10010,571403,46E+111,09E+103,3010,037809,82E+103,09E+093,489,4901202,78E+108,77E+083,708,9431607,90E+092,49E+083,858,3962002,24E+097,06E+074,007,8492406,35E+082,00E+074,047,3012801,80E+085,68E+064,086,7543601,45E+074,57E+054,125,660

Распад Pd

Расчетная формула:

Так как палладий гораздо более долгоживущий, чем его предшественники, то расчет ведем по простейшей формуле, при этом даже не учитывая начальные числа ядер его предшественников, т.к. они не внесут особого вклада в эволюцию палладия.

За небольшое по сравнению с периодом полураспада палладия время выдержки его активность не изменится и для построения графика можно ограничиться 4-мя точками.

Результаты расчета приведены в таблице 28.

Таблица 28.Эволюция Pd в целом за время выдержки.

t,секN,штA,БкlgtlgA5000001,03E+253,47E+105,7010,54110000001,03E+253,47E+106,0010,541100000001,03E+253,47E+107,0010,541881280001,03E+253,47E+107,9510,541

Графические зависимости lg A-lg tK, lg A-lg Tохл для радионуклидов цепочки А= 107 по данным таблиц 16 - 21, 23-28 представлены на рисунке 9.

Рисунок 9. Графические зависимости lg A-lg tK, lg A-lg Tохл для радионуклидов цепочки А= 107.

3. Определение содержания стабильного (конечного) члена цепочки

к концу компании:

N?=F*знез*tk=9,21·1019*0,00142*923*86400=10,43*1024 шт,

где знез-общий независимый выход цепочки.

F-скорость деления чдер.

к концу выдержки:

Nост= ?Nчл.цеп.

В этой формуле не будем учитывать те радионуклиды цепочки, у которых время уменьшения их активности в 10 раз относительно первоначальной не сравнимо с временем выдержки - 1020 сут, значит нет смысла для них рассчитывать числа оставшихся ядер на момент окончания выдержки, т.е. для Mo, Tc, Ru, Rh, Pdm, Agm. Для оставшегося Pd107 с Т1/2=6,5·106 лет рассчитаем число оставшихся ядер цепочки по формуле Бейтмана за время охлаждения.

Число оставшихся ядер для Pd107 рассчитываем по формуле:

NостPd= =1,03E+25шт.

Nстаб= N?- Nост = 10,43*1024 - 1,03E+25=0,013*1025 шт- число оставшихся стабильных ядер цепочки на момент окончания выдержки.

mстаб=[Nстаб · MAg] / NA=0,013*1025*107,8682/6,02*1023=23,294 г- масса оставшихся стабильных ядер цепочки на момент окончания выдержки.

3.1Эволюция радионуклидов с предложенным Z = 35 в координатах lg A-lg Tохл

Из [3] выписываем изменение активности радионуклидов элемента с Z = 35 (Br) во время охлаждения для реактора ВВЭР-1000 с обогащением 3,95%.

В [3] приведен только 1 радионуклид Br, у которого Т1/2=35,30 ч.

Для него из [3] выписываем значения удельной активности в зависимости от времени. Справочные данные приведены в таблице 29.

Таблица29.Изменение активности, расп/(с·т-мет) ПД во время охлаждения.

радионуклидТ1/260 сут90 сут120 сут180 сут1 г2,3 г3 г10 летBr-8235,30 ч8,08Е+1-------

Для уточнения данных воспользуемся источником [1], из которого выпишем все изотопы брома.

По периодам полураспада данных радионуклидов можно сказать, что время снижения их активности в 10 раз относительно первоначальной не сравнимо с временем охлаждения - 1020 сут., поэтому нет смысла строить для них графическую зависимость, т.к. они полностью распадутся.

Самый долгоживущий из них радионуклид с периодом полураспада равным 35,30 часа уже ко времени равном 90 суток распадется полностью. (Мы можем говорить об этом, опираясь на данные, выписанные из [3])

Эволюция стабильных изотопов брома, их 2 шт., будет определяться эволюцией цепочек, которым они принадлежат.

3.2 Определение скорости образования гремучего газа при растворении ТВЭЛов в л/мин·кг топл. и общего количества выделившегося газа за 6 часов

Время растворения ТВЭЛов:

сек;

Суммарная в -активность продуктов деления на конец времени охлаждения (п. 2.4.2.):

Бк;

Начальный выход водорода (с.41[1]):

шт/100 эв;

Объем водорода за 6 часов растворения (с.41[1]):

л;

где: 2*106 - средняя энергия в -частиц, эВ;

Поскольку гремучий газ - это смесь Н2 и О2 в соотношении 2:1, то:

л;

Скорость образования гремучего газа (с.41[1]):

л/мин*кг;

топливо ядерный реактор изотоп

Список используемой литературы

1.Особенности эволюции радионуклидов в природных и техногенных системах/ Учебное пособие / Н.Д. Бетенеков, Ю.В. Егоров, Т.А. Недобух, В.Д. Пузако. Екатеринбург: ГОУ УГТУ-УПИ, 2001, 72с.

2.Гусев Н.Г. и др. Гамма-излучение радиоактивных изотопов и продуктов деления.-М:НФМЛ,1958.-208 с.

.Колобашкин В.М. и др. Радиационные характеристики облученного ядерного топлива: Справочник.-М.:Энергоатомиздат.1983.-384 с.

.Курс лекций по курсу Радиохимии III часть, 2006г. Н.Д. Бетенеков.

Расчет характеристик ядерного реактора 1. Исходные данные Исходные данные из [1] приведены

Больше работ по теме: