14 / 19
А.А. Грешилов
14
Для идентификации взрыва по двум измеренным изотопам при-
меняют объединение двух видов деления урана 235 (нейтронами
спектра деления и нейтронами с энергией 14 МэВ) в один вид деле-
ния и объединение двух видов деления плутония 239 (нейтронами
спектра деления и нейтронами с энергией 14 МэВ) в один вид деле-
ния, что приводит к сокращению числа идентифицируемых видов
деления (вместо четырех рассматривают два).
При этом удельную активность рассчитывают по формуле (13):
для делящегося материала урана 235 с вектором независимых
выходов
235
235
235
14
U
U
U
1
1
1
f
c
c
,
где
235
U
f
и
235
14
U
— независимые выходы элементов изобарных це-
почек при делении урана 235 нейтронами спектра деления и нейтро-
нами с энергией 14 МэВ соответственно;
1
c
— параметр, учитываю-
щий доли независимых выходов
235
U
f
и
235
14
U
в их сумме,
1
0, 1
c
;
для делящегося материала плутония 239 с вектором
независимых выходов
239
239
239
14
Pu
Pu
Pu
2
2
(1 )
f
c
c
,
где
239
Pu
f
и
239
14
Pu
— независимые выходы элементов изобарных це-
почек при делении плутония 239 нейтронами спектра деления и
нейтронами с энергией 14 МэВ соответственно;
2
c
— параметр, учи-
тывающий доли независимых выходов
239
Pu
f
и
239
14
Pu
в их сумме,
2
0, 1
c
.
Задавая двумерную сетку по
1
c
и
2
c
с шагом
1
c
и
2
c
соответ-
ственно и находя минимум (17) для разных
1
c
и
2
c
, за
истинные
вклады источников
235
U
N
и
239
Pu
N
в суммарную активность изото-
пов ксенона принимают те, при которых сумма квадратов невязок
системы (13) м и н и м а л ь н а.
Проведено
имитационное моделирование реализации предла-
гаемого способа
на персональном компьютере с процессором Intel
Celeron 2,40 ГГц с объемом оперативной памяти 768 Мбайт в мате-
матическом пакете Matlab 7.0.
Имитировали ситуацию отбора пробы через 6 дней после взрыва
и измерения активностей пяти изотопов (
85
Kr
m
,
131
Xe
m
,
133
Xe
m
,
133
Xe
,
135
Xe
). Результаты рассчитаны при условии, что момент сепа-
рации предполагается известным: 3 ч после события. Значения изме-