На этой фотографии представлен общий
вид применяемого в указанных исследованиях спектрометра Лаборатории
радиационного контроля ЛРК-1 МИФИ. Спектрометр предназначен для
определения радионуклидного состава и активности радионуклидов в
пробах окружающей среды методами низкофоновой гамма-спектрометрии
с полупроводниковыми детекторами в стационарных лабораторных условиях.
Спектрометр располагается в полуподвальном проветриваемом помещении.
В состав спектрометра входят 7 индивидуальных
спектрометрических трактов (ИСТ), устройство накопления, обработки
и отображения результатов на базе многоканального анализатора импульсов
NOKIA LP-4900B с программным управлением, устройство обработки результатов
на базе персональной ЭВМ типа PC/AT.
Каждый ИСТ состоит из следующих основных элементов:
Блоки линейных усилителей и высоковольтного питания размещаются в крейтах "Вектор ВБЦ2-90", ORTEC 401, ORTEC 4001.
Собственно низкофоновая защита выполнена по модульному принципу и состоит из:
Применение внутренних модулей различных
конфигураций позволяет помещать в камеру измерительные контейнеры
различной геометрии, а также использовать детекторы разных конструкций.
Подъемный механизм обеспечивает доступ внутрь защиты при смене проб,
замене внутренних модулей и при извлечении детектора.
Вес защитной камеры в сборе с блоком детектирования не превышает
600 кг; габаритные размеры, включая пространство для открывания/сдвига
крышки: 0,6х0,8х1,3 м (WxDxH).
Для заливки азота используется метод вытеснения воздушным давлением: по трубке, введенной в горловину сосуда Дьюара, или заливной патрубок криостата детектора, подается азот из внешнего резервуара, в свободный объем которого насосом нагнетается воздух.
Для уменьшения вклада в фон установки гамма-излучения дочерних продуктов радионуклида 226Ra, внутренний объем защитной камеры постоянно продувается парами жидкого азота.
Исследуемые пробы помещаются в цилиндрические контейнеры объемом 145 см3 (диам. 95х30 мм) или 250 см3 (диам.95х50 мм), либо в контейнер с геометрией Маринелли объемом 500 см3 или 1000 см3. В первых двух случаях контейнеры устанавливаются в закрепленный непосредственно на крышке колпака криостата детектора фиксатор, воспроизводящий геометрию "проба-детектор".
Ниже представлены некоторые фоновые характеристики спектрометра.
Спектр фона для ORTEC GEM 30185-p+. |
|
Спектр фона для Canberra GL3830. |
||||||||||
|
Nuclide |
Energy, keV |
Count rate, 1/ksec |
|
|
Nuclide |
Energy, keV |
Count rate, 1/ksec |
||||
1. |
Ra-226 |
185.97 |
3.09 |
± |
0.45 |
|
1. |
Ra-226 |
185.97 |
|
< |
0.25 |
2. |
Th-232 |
238.63 |
2.36 |
± |
0.29 |
|
2. |
Th-232 |
238.63 |
|
< |
0.16 |
3. |
Th-232 |
240.99 |
|
< |
0.27 |
|
3. |
Th-232 |
240.99 |
|
< |
0.16 |
4. |
Ra-226 |
241.97 |
|
< |
0.29 |
|
4. |
Ra-226 |
241.97 |
|
< |
0.16 |
5. |
Ra-226 |
295.22 |
|
< |
0.31 |
|
5. |
Ra-226 |
295.22 |
|
< |
0.14 |
6. |
Th-232 |
338.33 |
0.44 |
± |
0.15 |
|
6. |
Th-232 |
338.33 |
|
< |
0.17 |
7. |
Ra-226 |
351.92 |
0.33 |
± |
0.21 |
|
7. |
Ra-226 |
351.92 |
|
< |
0.13 |
8. |
Th-232 |
463.01 |
|
< |
0.22 |
|
8. |
Th-232 |
463.01 |
|
< |
0.13 |
9. |
Annihil |
511.00 |
6.88 |
± |
0.35 |
|
9. |
Annihil |
511.00 |
4.52 |
± |
0.27 |
10. |
Pb(n,n') |
569.60 |
0.42 |
± |
0.10 |
|
10. |
Pb(n,n') |
569.60 |
|
< |
0.16 |
11. |
Th-232 |
583.19 |
0.87 |
± |
0.18 |
|
11. |
Th-232 |
583.19 |
0.14 |
± |
0.10 |
12. |
Ra-226 |
609.32 |
0.31 |
± |
0.13 |
|
12. |
Ra-226 |
609.32 |
|
< |
0.18 |
13. |
Cs-137 |
661.66 |
0.26 |
± |
0.15 |
|
13. |
Cs-137 |
661.66 |
|
< |
0.09 |
14. |
Cu(n,n') |
669.60 |
0.22 |
± |
0.16 |
|
14. |
Cu(n,n') |
669.60 |
|
< |
0.10 |
15. |
Th-232 |
727.33 |
0.21 |
± |
0.11 |
|
15. |
Th-232 |
727.33 |
|
< |
0.07 |
16. |
U-238 |
766.60 |
|
< |
0.17 |
|
16. |
U-238 |
766.60 |
|
< |
0.07 |
17. |
Th-232 |
794.94 |
0.11 |
± |
0.10 |
|
17. |
Th-232 |
794.94 |
|
< |
0.11 |
18. |
Pb(n,n') |
803.10 |
0.30 |
± |
0.12 |
|
18. |
Pb(n,n') |
803.10 |
|
< |
0.11 |
19. |
Fe + n |
847.00 |
|
< |
0.09 |
|
19. |
Fe + n |
847.00 |
|
< |
0.07 |
20. |
Th-232 |
860.56 |
0.16 |
± |
0.10 |
|
20. |
Th-232 |
860.56 |
|
< |
0.07 |
21. |
Th-232 |
911.21 |
0.32 |
± |
0.11 |
|
21. |
Th-232 |
911.21 |
|
< |
0.07 |
22. |
Cu(n,n') |
962.10 |
0.16 |
± |
0.10 |
|
22. |
Cu(n,n') |
962.10 |
|
< |
0.10 |
23. |
Th-232 |
964.76 |
|
< |
0.10 |
|
23. |
Th-232 |
964.76 |
|
< |
0.10 |
24. |
Th-232 |
968.97 |
0.13 |
± |
0.10 |
|
24. |
Th-232 |
968.97 |
|
< |
0.09 |
25. |
U-238 |
1001.00 |
0.22 |
± |
0.10 |
|
25. |
U-238 |
1001.00 |
|
< |
0.06 |
26. |
Cu(n,n') |
1115.50 |
|
< |
0.17 |
|
26. |
Cu(n,n') |
1115.50 |
|
< |
0.08 |
27. |
Ra-226 |
1120.31 |
|
< |
0.14 |
|
27. |
Ra-226 |
1120.31 |
|
< |
0.06 |
28. |
Co-60 |
1173.20 |
0.09 |
± |
0.09 |
|
28. |
Co-60 |
1173.20 |
|
< |
0.09 |
29. |
Co-60 |
1332.50 |
|
< |
0.11 |
|
29. |
Co-60 |
1332.50 |
|
< |
0.07 |
30. |
K-40 |
1460.80 |
0.57 |
± |
0.09 |
|
30. |
K-40 |
1460.80 |
0.12 |
± |
0.05 |
31. |
Ra-226 |
1764.55 |
0.18 |
± |
0.05 |
|
31. |
Ra-226 |
1764.55 |
|
< |
0.05 |
32. |
Ra-226 |
2204.17 |
0.08 |
± |
0.07 |
|
32. |
Ra-226 |
2204.17 |
|
< |
0.06 |
33. |
Ra-226 |
2447.84 |
|
< |
0.08 |
|
33. |
Ra-226 |
2447.84 |
|
< |
0.05 |
34. |
Th-232 |
2614.53 |
0.46 |
± |
0.07 |
|
34. |
Th-232 |
2614.53 |
0.08 |
± |
0.07 |
Непрерывный фон, отсчетов/кэВ/ксек |
Непрерывный фон, отсчетов/кэВ/ксек |
||||||||||
1. |
100 keV |
2.489 |
± |
0.015 |
1. |
100 keV |
1.170 |
± |
0.009 |
||
2. |
200 keV |
2.422 |
± |
0.015 |
2. |
200 keV |
1.253 |
± |
0.009 |
||
3. |
300 keV |
1.426 |
± |
0.012 |
3. |
300 keV |
0.735 |
± |
0.007 |
||
4. |
400 keV |
0.849 |
± |
0.009 |
|
4. |
400 keV |
0.461 |
± |
0.006 |
|
5. |
500 keV |
0.635 |
± |
0.008 |
5. |
500 keV |
0.406 |
± |
0.005 |
||
6. |
600 keV |
0.411 |
± |
0.006 |
6. |
600 keV |
0.250 |
± |
0.004 |
||
7. |
700 keV |
0.313 |
± |
0.005 |
7. |
700 keV |
0.206 |
± |
0.004 |
||
8. |
800 keV |
0.239 |
± |
0.005 |
8. |
800 keV |
0.153 |
± |
0.003 |
||
9. |
1000 keV |
0.165 |
± |
0.004 |
9. |
1000 keV |
0.105 |
± |
0.003 |
||
10. |
1200 keV |
0.121 |
± |
0.003 |
10. |
1200 keV |
0.083 |
± |
0.002 |
||
11. |
1500 keV |
0.086 |
± |
0.003 |
|
11. |
1500 keV |
0.054 |
± |
0.002 |
|
12. |
2000 keV |
0.053 |
± |
0.002 |
|
12. |
2000 keV |
0.037 |
± |
0.002 |
|
13. |
2500 keV |
0.035 |
± |
0.002 |
13. |
2500 keV |
0.028 |
± |
0.001 |
Программно-математическое обеспечение (ПМО) управления процессом измерений, разработанное для анализатора LP-4900B, позволяет:
Необходимое для получения представительного
спектра время измерения одной пробы, определяемое ее видом, активностью,
радионуклидным составом и решаемой задачей, обычно составляет величину
от 2 до 60 часов. Обеспечение сохранности информации при длительных
измерениях реализуется посредством записи на ГМД, помимо результирующего
спектра измерения, каждого из последовательно набранных парциальных
спектров (длительность одного парциального измерения составляет
1-4 часа). Такая избыточность, кроме того, позволяет оператору,
не влияя на текущие измерения и не прекращая их, выполнять с помощью
данного ПМО суммирование ранее измеренных спектров с приведением
их к общей энергетической калибровке, что эквивалентно периодической
коррекции коэффициента преобразования ИСТ.
Суммарный спектр, накопленный за все время измерений пробы (либо
полученный оператором ручным суммированием парциальных спектров
с их коррекцией или без таковой), снабженный информацией об имени
пробы, дате и времени измерения, записывается на ГМД 5,25" в файловом
формате MS DOS и передается для последующей обработки на PC/AT.
Спектр также может быть передан из анализатора в PC/AT с использованием
параллельной линии связи.
Для количественного анализа результатов спектрометрических измерений образцов разработано ПМО, состоящее из трех частей:
1. Интерактивная обработка спектров.
Целью обработки спектра является определение скоростей счета и других параметров пиков полного поглощения и оценка их погрешностей. Исходными данными для обработки служат измеренный спектр и информация об образце и условиях его измерения, характеристики спектрометра и библиотека данных о радионуклидах и их гамма-излучении. Обработка проводится интерактивным методом. Оператор, проводящий обработку, руководствуется представленными на экране изображением спектра, таблицей линий гамма-излучения с их характеристиками и подсказками программы и может управлять форматом изображения, выбирать пики или группы пиков для анализа, задавать или корректировать выбранные программой границы пиков и групп пиков, идентифицировать пики или корректировать идентификацию пиков программой. Используемые в программе алгоритмы позволяют определять положения, площади и ПШПВ одиночных пиков, проводить разрешение мультиплетов с определением положений и площадей составляющих, оценивать верхний предел площади для внесенных в таблицу, но не обнаруженных в спектре пиков. В ходе обработки выполняется автоматическая настройка энергетической шкалы спектрометра на данный спектр. Для всех определяемых величин проводится оценка их статистических и систематических погрешностей.
2. Определение радионуклидного состава образца.
Содержание в образце радионуклидов из заданного списка проводится по определенным при обработке спектра скоростям счета в пиках полного поглощения. После вычитания фоновых скоростей счета и внесения коррекции на плотность спектрометрического образца с помощью хранящихся в библиотеке данных значений чувствительности спектрометра (отношение скорости счета в пике к активности источника) для каждой линии гамма-излучения усреднением по всем линиям вычисляются активности радионуклидов или оценки их верхних пределов. На основании полученных результатов готовится протокол измерений, в котором перечислены характеристики пробы и приготовленного из нее спектрометрического образца, описаны условия измерений и приведены значения удельных активностей радионуклидов в пробе в единицах Бк/кг на дату отбора пробы. Если проба характеризуется определенной поверхностью, то приводятся также значения удельных активностей в единицах Бк/м2, а если объемом, то в единицах мБк/см3. Для анализа результатов измерений проб с послойным отбором по глубине почвы разработана программа, аппроксимирующая распределение удельной активности по глубине одной или двумя экспонентами с определением параметров распределения и запаса активности радионуклида в почве. Погрешности получаемых значений активностей радионуклидов определяются систематическими погрешностями используемых данных (характеристики применяемых при калибровке по эффективности и чувствительности образцовых источников, параметры для коррекции на самопоглощение в образце и др.) и статистикой спектра и составляют от 5% до 15%. Примеры протокола результатов обработки приведены на рис. 4,5.
3. Подготовка данных для анализа.
Для проведения описанных выше операций необходимы следующие данные:
Для получения этих величин и внесения
их в наборы данных для анализа разработан ряд вспомогательных программ.
3.1. Фоновые скорости счета определяются
из экспериментальных фоновых спектров, измеряемых периодически для
каждого ИСТ спектрометра, с помощью программы обработки и заносятся
в соответствующие наборы данных для каждого ИСТ спектрометра.
3.2. Для определения значений чувствительности
выполняется калибровка по эффективности для каждого сочетания "детектор
+ контейнер образца". Процедура калибровки состоит в определении
значений и погрешностей параметров аналитической зависимости эффективности
регистрации в пике полного поглощения от энергии гамма-квантов по
результатам измерений образцовых источников. Для радионуклидов из
списка, совпадающих с используемыми образцовыми источниками (137Cs,
134Cs, 226Ra, 232Th, 40K,
57Со, 54Mn, 22Na, 65Zn,
139Ce, 113Sn, 152Eu), значения
чувствительности определяются усреднением результатов измерений
образцовых источников, для остальных радионуклидов - умножением
определяемой из калибровочной кривой эффективности для каждой гамма-линии
на ее квантовый выход.
3.3. Зависимость ПШПВ пиков от энергии
определяется из тех же измерений с образцовыми источниками и аппроксимируется
линейной функцией.
3.4. Для коррекции зависимости эффективности
регистрации гамма-квантов от плотности источника используется аналитическая
функция, задающая зависимость корректирующего коэффициента от плотности
образца и энергии гамма-квантов. Параметры этой функции определяются
методом наименьших квадратов для каждой комбинации "детектор + контейнер
образца" из результатов обработки спектров источников 226Ra
c плотностью от 0.2 до 2.5 г/см3.
Полученные таким образом данные
хранятся в отдельных файлах для каждого набора "детектор+защита+контейнер"
и используются при обработке спектров и при определении радионуклидного
состава пробы.
Для проверки работоспособности методик
накопления и обработки данных было произведено количественное определение
содержания 137Cs, 232Th и 152Eu
в образцах, предоставленных НИКИЭТ в рамках программы "Метрологическая
экспертиза гамма-спектрометров на предприятиях 16-го ГНТУ Министерства
атомной энергии РФ". Отношения полученных оценок удельных активностей
к их паспортным значениям составили:
137Cs: 1.01 ± 0.05; 232Th: 0.95 ± 0.09; 152Eu: 1.02 ± 0.10
Пример 1: протокол обработки результатов спектрометрического анализа пробы. МОСКОВСКИЙ ИНЖЕНЕРНО-ФИЗИЧЕСКИЙ ИНСТИТУТ (ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ)
Результаты определения содержания радионуклидов в пробе на дату отбора
Пример 2: протокол обработки результатов спектрометрического анализа пробы. МОСКОВСКИЙ ИНЖЕНЕРНО-ФИЗИЧЕСКИЙ ИНСТИТУТ (ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ)
Результаты определения содержания радионуклидов в пробе на дату отбора
|