«Определение атомного номера вещества объектов по ослаблению пучков фотонов с энергиями до 10 МэВ ...»

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ

УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ

"МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

имени М.В. ЛОМОНОСОВА"

ФИЗИЧЕСКИЙ ФАКУЛЬТЕТ

НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ЯДЕРНОЙ ФИЗИКИ

имени Д.В. СКОБЕЛЬЦЫНА

На правах рукописи

УДК 539.12.04 Курилик Александр Сергеевич Определение атомного номера вещества объектов по ослаблению пучков фотонов с энергиями до 10 МэВ Специальность 01.04.16 физика атомного ядра и элементарных частиц

ДИССЕРТАЦИЯ

на соискание учёной степени кандидата физико-математических наук

Научный руководитель д. ф.-м. н., профессор Б.С. Ишханов Москва — Оглавление Введение Актуальность темы диссертации............................. Задачи интроскопии................................. Существующие решения.............................. Цели диссертационной работы.............................. Основные результаты, полученные в диссертации................... Научная новизна...................................... Ценность научной работы................................ Практическая значимость................................ Аппробация работы.................................... Публикации........................................ Личный вклад автора................................... Структура и объём диссертации............................. Содержание диссертации................................. 1. Взаимодействие фотонов с веществом 1.1. Сечения процессов взаимодействия фотонов с веществом........... 1.2. Ослабление потока фотонов............................ 1.2.1. Монохроматические пучки фотонов.................... 1.2.2. Пучки фотонов с непрерывным спектром................ 1.3. Принцип измерения Z и nD............................. 1.3.1. Монохроматические пучки фотонов.................... 1.3.2. Тормозное излучение............................ 1.4. Анализ сечений................................... 1.4.1. Нормированные сечения.......................... 1.4.2. Вклады сечений основных процессов в полное сечение......... 1.4.3. Полные сечения............................... 1.4.4. Выводы.................................... 1.4.5. Рекомендации................................ 1.4.6. Соотнешение вероятностей процессов рассеяния и поглощения фотонов 2. Экспериментальные измерения атомного номера объектов 2.1. Радиоактивные источники Cs-137 и Co-60.................... 2.1.1. Выбор источников гамма-квантов..................... 2.1.2. Соотношение R в случае Cs-137 и Co-60................. 2.1.3. Описание эксперимента........................... 2.1.4. Результаты эксперимента и выводы.................... 2.2. Тормозное излучение................................ 2.2.1. Описание эксперимента........................... 2.2.2. Результаты эксперимента и выводы.................... 3. Моделирование на GEANT4 3.1. Схема установки для интроскопии крупногабаритных объектов........ 3.2. Тормозное излучение................................ 3.2.1. Взаимодействие электронов с веществом................. 3.2.2. Моделирование тормозного излучения.................. 3.2.3. Выбор оптимальной толщины тормозной мишени............ 3.3. Выбор типа и размера детектора......................... 3.4. Отклик детектора.................................. 4. Сравнение методов и оптимизация их параметров 4.1. Пучки монохроматических фотонов трёх энергий................ 4.2. Тормозное излучение с двумя граничными энергиями............. 4.3. Тормозное излучение с тремя граничными энергиями............. 4.6. Уменьшение уровня статистических флуктуаций................ 4.7.1. Тормозное излучение с двумя граничными энергиями......... 4.7.2. Тормозное излучение с тремя граничными энергиями......... 4.8. Сравнение в широком динамическом диапазоне................. Введение Актуальность темы диссертации Задачи интроскопии Часто существует необходимость получить информацию о внутренней структуре некоторого объекта. Например, обнаружить скрытые дефекты в продукции промышленного производства, проверить соответствие декларации и реального содержимого груза. В последние годы особую актуальность приобрела проблема обеспечения безопасности и предотвращение терроризма в морских и речных портах, аэропортах, железнодорожных вокзалах. Также важной задачей является контроль содержимого грузов и транспортных средств с целью выявления и пресечения незаконных перевозок — контрабанды. Одним из основных элементов системы обеспечения безопасности являются технические средства, позволяющие осуществлять неразрушающий контроль — исследовать, не вскрывая, грузы и транспортные средства.

В зависимости от целей применения, к интроскопическим системам предъявляются различные требования. Например, при таможенном досмотре основными, в некотором смысле взаимоисключающими требованиями являются: высокая пропускная способность системы и максимально возможное снижение мощности дозы, создаваемой установкой в процессе функционирования. Аналогичные требования предъявляются к медицинским установкам. К дефектоскопическим системам подобных жёстких требований может не предъявляться, так как основной целью является обнаружение брака. В исследовательских интроскопических установках основной ограничивающий фактор — радиационные повреждения, возникающие в исследуемом объекте.

В разных странах, в силу географических особенностей, также существуют различные требования, предъявляемые к таможенным системам. В США основной целью применения досмотровых комплексов является обнаружение оружия массового поражения и его составляющих, радиоактивных и взрывчатых веществ, оружия и наркотиков. Считается, что наибольшую опасность представляет ядерное оружие массового поражения (WMD — nuclear weapons of mass destruction) на основе высоко обогащённого урана (HEU) и оружейного плутония (WGPu) [1, 2], а также устройства распыляющие радиоактивные вещества (RDD — radiological dispersal device). Считается желательным обнаружение 100 см3, что соответствует кубу со стороной 4.6 см, вещества с большим. Веществом с большим (High- material) считается вещество с атомным номером > 72. В России основные усилия сконцентрированы на обнаружении оружия 26 и взрывчатых веществ 5.

Результатом работы досмотрового комплекса в общем случае должна быть информация о пространственном распределении атомного номера (,, ) и количества вещества (,, ) в досматриваемом объёме. Для этого весь досматриваемый объём разбивается на небольшие элементарные объёмы — т.н. воксели (Voxel — Volumetric Picture Element).

Однако получение трёхмерных изображений содержимого крупногабаритных грузов и транспортных средств за ограниченное время досмотра является трудновыполнимой задачей. Наиболее распространённым способом является получение одной двумерной проекции трёхмерного содержимого при просвечивании пучком фотонов. Результатом является изображение распределения прозрачности (, ), состоящее из пикселей (Pixel — Picture Element). На основе этой информации делается предположение об объектах, содержащихся в досматриваемом объёме. Затем проводится сравнение с декларацией и делается вывод о наличии несоответствий.

Важной характеристикой досмотрового комплекса является пропускная способность — количество досматриваемых грузов или транспортных средств в единицу времени. В первую очередь, она зависит от скорости сканирования. Обычно скорость перемещения исследуемого объекта относительно сканирующей системы составляет от 0.1 до 0.8 м/с.

Системы, просвечивающие грузы и транспорт пучками фотонов и/или нейтронов имеют ещё одну важную характеристику — расстояние между источником излучения и детектором. Это расстояние определяется поперечными размерами грузов и транспортных средств. Обычно это расстояние составляет от 3 до 10 м.

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

Существующие решения Информацию о внутренней структуре исследуемого объекта можно получить с помощью проникающего излучения и частиц. Применяемые и перспективные технологии, их преимущества, недостатки и области применения рассмотрены в [2–5]. Существующие и разрабатываемые методы можно классифицировать по различным признакам.

По типу источника проникающего излучения:

• рентгеновские установки с верхней границей энергии фотонов до сотен кэВ • радиоактивные источники -квантов с энергиями до нескольких МэВ • источники тормозного излучения на основе ускорителей электронов • импульсные и непрерывные источники тепловых и быстрых нейтронов • неионизирующее электромагнитное излучение, в том числе т.н. терагерцовое излучение (THz source) • пассивные методы, основанные на регистрации излучений, испускаемых при распаде нестабильных изотопов, содержащихся в исследуемом объекте По взаимодействиям и процессам, используемым для идентификации веществ и изотопов:

• взаимодействие фотонов с атомными электронами: фотоэффект, когерентное и некогерентное рассеяние, рождение /+ -пар • взаимодействие фотонов с атомными ядрами: рождение /+ -пар, резонансное рассеяние (флюоресценция), ГДР, (,n), (,f) • взаимодействие нейтронов с атомными ядрами: (n,), (n,n), (n,f) • взаимодействие электромагнитного излучения с молекулами и атомами • испускание -квантов атомными ядрами при переходе из возбуждённых состояний • испускание нейтронов деления По типу и энергии регистрируемых частиц:

• фотоны • нейтроны • заряженные частицы По функциональным возможностям системы:

• формирование 2D (2-dimensional) / 3D (3-dimensional) изображения распределения линейного коэффициента ослабления исследуемого объекта • формирование 2D / 3D изображения распределения эффективного атомного номера и количества вещества исследуемого объекта • формирование 2D / 3D изображения распределения изотопного состава исследуемого объекта • получение информации о наличии в объекте того или иного вещества • получение информации о наличии в объекте делящихся и радиоактивных материалов Радиоактивные и делящиеся материалы при своём распаде испускают излучение и частицы, которые могут быть зарегистрированы детекторами, расположенными вне досматриваемого объёма. Такие методы, основанные на регистрации продуктов самопроизвольного радиоактивного распада, называются пассивными. Эти системы позволяют обнаруживать даже незначительные количества радиоактивных веществ [6]. Однако, существенное количество делящихся и радиоактивных веществ может остаться незамеченным, если они будут перевозиться в специальном защищённом контейнере. В то же время, достаточно высока вероятность ложного срабатывания, т.к. многие обычные грузы содержат радиоактивные изотопы.

Т.н. активные методы предполагают облучение досматриваемого объёма и регистрацию излучения, прошедшего без взаимодействия, рассеянного излучения, вторичных частиц и продуктов ядерных реакций.

Заряженные частицы испытывают кулоновское рассеяние, вследствие чего обладают низкой проникающей способностью и поэтому не подходят для использования в качестве первичных. Возможно использование заряженных частиц высоких энергий. Однако получение таких частиц с помощью технических средств весьма расточительно. Существуют проекты обходящие это ограничение — использующие высокоэнергетичные мюоны космического излучения — т.н. мюонная томография [7].

Методы, использующие нейтронное излучение [8, 9] могут обеспечить высокую чувствительность и избирательность в определении веществ и изотопов. Однако эти методы требуют относительно большего времени, затрачиваемого на обследование одного объекта. Пространственное разрешение изображений, получаемых при просвечивании пучками нейтронов невысокое, 2 см из-за сложности создания необходимых потоков нейтронов и их регистрации. Кроме того, срок службы нейтронного источника невелик — порядка 2000 часов.

Ядерная резонансная флюоресценция является весьма перспективной технологией, позволяющей быстро и точно определять пространственное расположение и изотопный состав содержимого инспектируемого объекта [10,11]. Однако, данное направление требует дополнительных исследований.

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

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

Также возможно комбинирование в одной системе различных типов первичных частиц, детекторов и анализируемых физических процессов. Например, нейтронов с энергией 14 МэВ и фотонов от радиоактивных источников или ускорителей [12, 13].

Задача просвечивания объектов небольшого размера решена давно. Исследуемый объм просвечивается пучками фотонов, получаемых с помощью рентгеновских трубок. Просвечивание объекта с одной стороны с использованием пучка одной граничной энергии позволяет получить только изображение прозрачности (, ). Из набора двумерных проекций, полученных при просвечивании с разных сторон, можно восстановить трёхмерное распределение линейного коэффициента ослабления (,, ). Просвечивание излучением с двумя граничными энергиями позволяет восстановить распределение атомного номера (, ) и количества вещества (, ). Возможность измерения атомного номера обусловлена различной энергетической зависимостью полных сечений взаимодействия фотонов с атомами различных химических элементов. Измерив прозрачности при двух разных энергиях и сопоставив с известными заранее зависимостями, можно сделать вывод о вещества исследуемого объекта. При просвечивании объекта с разных сторон рентгеновскими пучками двух разных энергий можно восстановить наиболее полную информацию об объёмном распределении (,, ) и (,, ) [14,15]. Т.к. в большинстве случаев материалы, из которых изготавливаются различные объекты, не состоят из атомов одного химического элемента с одним значением, а состоят из атомов различных химических элементов — вводится понятие эффективного атомного номера.

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

Рентгеновское излучение также применимо для исследования приповерхностного слоя крупногабаритных объектов на предмет скрытых вложений в двойных стенках, методом так называемого обратного рассеяния [16, 17]. Существующие в настоящее время системы позволяют различать не только форму объектов, но и материалы, из которых они изготовлены [18, 19].

-источники давно и успешно применяются для дефектоскопии. Существуют как проекты досмотровых комплексов на основе источников 55 Cs или 27 Co [17, 20], так и уже действующие -сканеры [21]. Также некоторыми авторами было предложено анализируя двумерные проекции, полученные под несколькими разными углами, восстанавливать трёхмерную информацию о внутреннем содержимом объекта [22–24]. В нефтяной промышленности для контроля параметров смеси газа, нефти и воды применяется просвечивание -квантами от двух разных радиоактивных изотопов [25]. Обычно применяются изотопы, дающие -кванты небольших энергий. Рассматривается возможность использования изотопов 55 Cs и 27 Co для измерения объектов, расположенных в крупногабаритных контейнерах [26, 27].

Типовые активности источников, используемых в гамма-сканерах, составляют 55 Cs 7.4 · 1010 Бк, источника 60 Co до 18.5 · 1010 Бк. Срок службы цезиевого источника рассчитан на 15 лет, кобальтового — на 5 лет, самого досмотрового комплекса — на 30 лет.

Радиоактивные изотопы при распаде могут испускать -кванты с энергиями до нескольких МэВ [28]. Что существенно больше максимальной энергии фотонов серийно выпускаемых рентгеновских трубок не превышающей 450 кэВ [29].

При энергиях фотонов от 150 кэВ до 3.7 МэВ полные сечения взаимодействия () для всех монотонно убывают с увеличением энергии фотонов. Следовательно, поток фотонов с большей энергией меньше ослабляется при прохождении через поглотитель.

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

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

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

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

Крупногабаритные объекты обладают большой поглощающей способностью, что делает невозможным применение рентгеновских трубок и радиоактивных источников изза небольших энергий и малых интенсивностей излучения. Для того, чтобы исследовать массивные объекты, используют ускорители электронов с тормозной мишенью, которые позволяют получать интенсивные потоки тормозных фотонов больших энергий. Существуют системы, просвечивающие досматриваемый объём и дающие в качестве результата двумерное изображение прозрачности (, ) [31–33]. Кроме стационарных, разработаны мобильные комплексы [34]. Разработаны методы, позволяющие восстанавливать пространственное распределение плотности в инспектируемом объёме [22, 35, 36].

В работах С. Огородникова и В. Петрунина [37–41] был предложен метод распознавания веществ, основанный на анализе ослабления высокоэнергетичного тормозного излучения с двумя граничными энергиями, позволяющий различать вещества по эффективному атомному номеру и получать двумерные изображения (, ) и (, ). Реализован полностью функционирующий досмотровый комплекс, основанный на этих работах [42]. Проекты аналогичных систем описаны рядом авторов в работах [43–46]. Типичые значения граничных энергий, в предлагаемых проектах, составляют 6-9 и 3-6 для большей и меньшей энергий соответственно. Метод трёх энергий, упомянутый в работе С. Огородникова и В. Петрунина [37] и описанный в работе [47], является следующим шагом в интроскопии, позволяющим улучшить точность распознавания атомного номера неизвестных объектов с большим атомным номером. Также, в [48] описан проект системы с четырьмя ускорителями разных энергий, а в [49] — на основе одного ускорителя с переключением четырёх энергий.

Следует учитывать, что спектр фотонов, получаемых с помощью ускорителя, не моноэнергетический, что является причиной значительных сложностей [50, 51] при измерении сечений различных процессов взаимодействия фотонов с атомами и атомными ядрами.

Следовательно, те же проблемы возникнут и при решении другой обратной задачи — измерении атомного номера вещества исследуемого неизвестного объекта. Кроме того, для измерения необходима высокая стабильность от импульса к импульсу тока, энергии и положения электронного пучка, что приводит к увеличению стоимости как самого ускорителя, так и его обслуживания.

Важной характеристикой возможности просвечивать объекты является проникающая способность — максимальная толщина слоя железа, за которой инспекционная система может обнаружить наличие скрытого объекта. Для систем на основе ускорителей, обычно указывается величина проникающаей способности 40 см железа [31]. Для -сканеров проникающая способность указывается от 5 [21] до 19 [17] см железа. Проникающая способность зависит от активности источника, расстояния от источника до детектора, высоты и ширины элемента детектора, скорости движения досматриваемого объёма и, при изменении указанных параметров, может варьироваться в достаточно широких пределах.

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

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

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

Основные результаты, полученные в диссертации 1) В результате анализа поведения сечений фотоэффекта, когерентного и некогерентного рассеяния, эффекта рождения электрон-позитронных пар, а также полных сечений взаимодействия фотонов с энергиями от 0.5 до 10 МэВ с атомами и атомными ядрами в зависимости от энергии фотона и от атомного номера химического элемента, вкладов сечений отдельных процессов в полные сечения было установлено, что для определения элементов с меньше железа необходимо, чтобы одна из энергий фотонов была >1.5 МэВ, а для определения элементов с больше вольфрама необходимо, чтобы одна из энергий фотонов была 40.

4) Впервые проведено экспериментальное измерение атомного номера вещества тестовых образцов посредством сравнения величин ослабления пучков тормозного излучения при четырёх энергиях ускорителя электронов с перестраиваемой энергией — 7.3, 6.2, 5.5 и 3.5 МэВ. Согласие измеренных и ожидаемых значений подтверждает необходимость использования как минимум трёх энергий ускорителя электронов.

5) Впервые определены области атомных номеров и величин ослабления 1/ объектов, в пределах которых, при просвечивании тормозным излучением при двух энергиях ускорителя электронов, атомный номер измеряется однозначно. Определены величины отличий между измеренными и истинными значениями в области, в которой измеряется неоднозначно. Показано, что причиной возникновения неоднозначности является пересечение параметрических кривых, соответствующих различным в двухмерном пространстве прозрачностей.

6) Впервые показана возможность однозначного измерения атомного номера объектов при просвечивании тормозным излучением при трёх энергиях ускорителя электронов. Проведено сравнение методов, использующих тормозное излучение с двумя и тремя энергиями. Показано, что использование трёх энергий позволяет в ряде случаев существенно улучшить результаты измерения.

7) Впервые получены оптимальные значения энергий ускорителя электронов, для методов с двумя и тремя энергиями, при которых достигается наилучшая точность измерения атомных номеров > 74, посредством сопоставления измеренных величин ослабления тормозного излучения. Впервые выполнено сравнение точностей методов измерения атомных номеров в широком диапазоне и радиоскопических прозрачностей при различных параметрах пучков тормозного излучения. Показано преимущество метода использующего три энергии.

Научная новизна Научная новизна представленной работы обусловлена тем, что на основе оценок повышения точности измерения атомного номера, полученных с использованием разработанного комплекса программ, впервые показаны принципиальные преимущества метода, использующего три и более энергий по сравнению с методом, использующим только две энергии. Сравнение результатов моделирования и проведённых экспериментов показывает, что использование тормозного излучения с тремя и более граничными энергиями позволяет повысить точность измерения атомного номера вещества объектов.

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

Практическая значимость Практическая значимость данной работы, в первую очередь, связана с развитием метода исследования объектов посредством просвечивания пучками тормозного или гаммаизлучения. Показаны области применимости методов измерения атомного номера и количества вещества объектов, как при использовании источников тормозных фотонов с несколькими граничными энергиями, так и при использовании изотопов Cs-137 и Co-60.

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

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

• Многопучковый метод сканирования объектов / Б.С. Ишханов, А.С. Курилик, Д.С. Руденко и др. // Известия РАН. Серия физическая. — 2008. — Т. 72, No 6.

— С. 908–911. [47] • Измерения атомного номера неизвестных объектов на ускорителе электронов с перестраиваемой энергией / А.Н. Ермаков, Б.С. Ишханов, А.С. Курилик и др. // Известия РАН. Серия физическая. — 2012. — Т. 76, No 11. — С. 1337–1339. [52] • Измерения атомного номера и количества вещества объектов с помощью радиоактивных источников / С.С. Белышев, Б.С. Ишханов, А.С. Курилик, К.А. Стопани // Известия РАН. Серия физическая. — 2013. — Т. 77, No 7. — С. 1008–1011. [53] • Определение выходов ядерных реакций на основе анализа цепочек распадов / С.С. Белышев, К.А. Стопани, С.Ю. Трощиев и др. // ВМУ. Серия 3. Физика. Астрономия. — 2011. — No 4. — С. 42–47. [54] • Studying photonuclear reactions using the activation technique / S. Belyshev, A. Ermakov, B. Ishkhanov et al. // Nuclear Instruments and Methods in Physics Research Section A. — 2014. — Vol. 745. — P. 133–137. [55] Доклады на конференциях Материалы диссертации доложены лично автором и опубликованы в трудах следующих конференций и совещаний:

• Курилик А., Руденко Д., Стопани К. Разработка эффективных методов сканирования объектов // Труды шестой Межвузовской научной школы молодых специалистов "Концентрированные потоки энергии в космической технике, электронике, экологии и медицине". — Москва, 2005. — С. 102–106. [56] • Курилик А., Руденко Д., Стопани К. Много-пучковый метод сканирования объектов // Труды VII межвузовской научной школы молодых специалистов "Концентрированные потоки энергии в космической технике, электронике, экологии и медицине".

— 2006. — С. 83–88. [57] • Многопучковый метод определения атомного номера / Б.С. Ишханов, А.С. Курилик, Д.С. Руденко и др. // Труды VIII межвузовской научной школы молодых специалистов "Концентрированные потоки энергии в космической технике, электронике, экологии и медицине". — Издательство МГУ Москва, 2007. — С. 160–164. [58] • Kurilik A., Rudenko D., Stopani K. Recognition of structure and composition of objects by scanning beams of gamma quantas // LVII International conference on nuclear physics "NUCLEUS 2007". Book of abstracts. — Saint-Petersburg, 2007. — P. 295. [59] • Определение атомного номера методом многопучкового сканирования объектов / С.С. Белышев, А.А. Кузнецов, А.С. Курилик, К.А. Стопани // Тезисы докладов международного совещания по ядерной спектроскопии и структуре атомного ядра "ЯДРО 2008". — Saint-Petersburg, 2008. — С. 281. [60] • Многопучковый метод определения структуры и состава объектов / С.С. Белышев, А.А. Кузнецов, А.С. Курилик, К.А. Стопани // Сборник тезисов, материалы Четырнадцатой Всероссийской научной конференции студентов-физиков и молодых учёных "ВНКСФ-14". — Уфа, 2008. — С. 285. [61] • Кузнецов А., Курилик А. Калибровка эффективности HPGe детектора по экспериментальным измерениям и моделированию Geant4 // Тезисы докладов 59 международного совещания по ядерной спектроскопии и структуре атомного ядра "ЯдроСПбГУ Санкт-Петербург, 2009. — С. 314. [62] • Methods of detection of hazardous materials using an electron accelerator / S.S. Belyshev, A.A. Kuznetsov, A.S. Kurilik, K.A. Stopani // LX International Conference on Nuclear Physics NUCLEUS 2010 Methods of Nuclear Physics for Femto- and Nanotechnologies, July 6 – 9, 2010, Saint-Petersburg, Russia. Book of Abstracts. — Saint-Petersburg, 2010.

• Kurilik A. Optimal photon energies to measure atomic numbers of objects // LXIII International Conference NUCLEUS 2013 "Fundamental Problems of Nuclear Physics and Atomic Power Engineering October 8 – 12, 2013, Moscow, Russia. Book of Abstracts.

— Saint-Petersburg, 2013. — P. 254. [64] Также материалы диссертации доложены и опубликованы в трудах следующего совещания:

• Автоматизация обработки данных гамма-активационных экспериментов / С.С. Белышев, А.А. Кузнецов, А.С. Курилик, К.А. Стопани // Тезисы докладов 58 международного совещания по ядерной спектроскопии и структуре атомного ядра "ЯдроСанкт-Петербург, 2008. — С. 282. [65] Материалы диссертации докладывались лично автором на научных конференциях "Ломоносовские чтения"в 2006, 2007, 2009, 2010, 2011 и 2013 годах, а также на выставкеконкурсе "Инновационный проект 2008".

Публикации Основные результаты по теме диссертации изложены в 15 печатных работах, 5 из которых изданы в журналах, рекомендованных ВАК.

Личный вклад автора Автором был разработан комплекс программ, позволяющий моделировать, обрабатывать и анализировать процесс генерации, прохождения через объекты и регистрации фотонов тормозного и гамма-излучения. Автором были выполнены сравнение и анализ сечений различных процессов взаимодействия фотонов с атомами различных веществ. С помощью пакета программ GEANT4 автором было выполнено моделирование углового и энергетического распределения тормозного излучения при различных энергиях электронного пучка и толщины вольфрамовой мишени; моделирование отклика и эффективности линейки сцинтилляционных детекторов, дозиметра БДКС-96, HPGe-детектора Canberra GC3019. В рамках работ по калибровке HPGe-детектора Canberra GC3019, автором была создана модель этого детектора в GEANT4 и были получены оценки неизвестных параметров детектора. Автором был проведен эксперимент по просвечиванию тестовых объектов пучками тормозных фотонов при четырёх значениях граничной энергии. Автором было осуществлено планирование эксперимента, проведение измерений и обработка полученных данных. Автором был проведен эксперимент по просвечиванию тестовых объектов гамма-квантами от радиоактивных источников. Автором было осуществлено планирование эксперимента, проведение измерений и обработка полученных данных. С помощью разработанного автором комплекса программ, были показаны области применимости методов, использующих, как ускорители электронов с перестраиваемой энергией, так и различные радиоактивные источники. Автором показано преимущество использования как минимум трёх энергий ускорителя по сравнению с использованием только двух энергий. С помощью разработанного автором комплекса программ, были получены оценки значений энергий ускорителя, при которых достигается наибольшая возможная точность измерения атомного номера вещества объектов.

Структура и объём диссертации Диссертация состоит из введения, четырех глав и заключения. Полный объем диссертации 111 страниц текста с 70 рисунками и 12 таблицами. Список литературы содержит 106 наименований. Текст диссертации подготовлен с помощью системы L TEX [66]. Список литературы оформлен с помощью BibTEX [67]. Использован стиль gost2008 [68], соответствующий ГОСТ Р 7.0.5-2008 введённому в действие 01.01.2009.

Содержание диссертации Во введении даётся краткий обзор методов изучения внутренней структуры и состава неизвестных объектов. Описываются цели интроскопии вообще и цели применения досмотровых комплексов в частности. Обосновывается актуальность исследований в этой области. Приводится классификация различных технологий и методов, используемых в интроскопии, их возможности и особенности.

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

В первой главе рассматриваются основные процессы взаимодействия фотонов с энергией до 10 МэВ с атомами вещества, зависимости сечений от энергии фотона и от атомного номера вещества. Описывается ослабление пучка фотонов при прохождении через слой вещества, как в случае монохроматических пучков, так и пучков с непрерывным спектром.

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

Во второй главе описываются проведённые автором эксперименты по просвечиванию тестовых объектов фотонами от радиоактивных источников и от ускорителя электронов с перестраиваемой энергией. В первой части главы обосновывается выбор изотопов 55 Cs и 27 Co в качестве источников гамма-квантов. Анализируется зависимость отношений сечений от при энергиях гамма-квантов, характерных для выбранных источников. Приводится описание эксперимента и его результаты. Делается вывод о возможносях данного метода. В следующей части главы описывается эксперимент с использованием ускорителя с перестраиваемой энергией. Приводятся его результаты и делаются выводы.

В третьей главе описывается проведённое автором моделирование отдельных элементов интроскопической системы, выполненное с использованием пакета программ GEANT4.

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

Четвёртая глава посвящена описанию методики оптимизации интроскопических систем. В качестве основных параметров системы, влияющих на точность измерения, выбраны количество и значения энергий монохроматических пучков фотонов или граничных энергий тормозного излучения. Показывается потенциальная возможность однозначного и устойчивого измерения с помощью трёх пучков монохроматических фотонов. Подробно рассматриваются возможности типовой инспекционной системы, использующей две энергии тормозного излучения. Показаны преимущества системы на основе трёх энергий тормозного излучения. Перечислены различные факторы влияющие на точность измерения. Проведено модельное сравнение точности измерения одинаковых объектов с помощю систем с двумя и тремя граничными энергиями тормозного излучения. Сформулирован критерий, на основе которого проводилась оптимизация значений граничных энергий тормозного излучения. Проведён выбор оптимальных значений энергий для случаев двух и трёх граничных энергий. Далее приводится описание модельного сравнения оптимизированных систем с двумя и с тремя граничными энергиями в широком динамическом диапазоне прозрачностей при равных средних мощностях электронного ускорителя.

Показано преимущество систем с тремя энергиями.

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

Глава 1.

Взаимодействие фотонов с веществом 1.1. Сечения процессов взаимодействия фотонов с веществом Рассмотрим основные процессы взаимодействия фотонов с веществом и зависимости вероятностей этих процессов от энергии налетающего фотона и атомного номера вещества, с которым происходит взаимодействие. Как известно [69–71], фотоны с энергией до 10 МэВ взаимодействуют с атомными ядрами и атомами вещества посредством трёх основных физических процессов: фотоэлектрического поглощения, комптоновского рассеяния и образования пар. Кроме того, с меньшей вероятностью могут происходить релеевское рассеяние и фотоядерные взаимодействия.

Фотоэффект — поглощение фотона атомом, при котором один из атомных электронов покидает атом. Энергия фотона должна быть больше энергии связи электрона в атоме.

Зависимость сечения фотоэффекта от атомного номера вещества и энергии фотона, имеет вид:

Зависимости сечения фотоэффекта ph от энергии фотона и от атомного номера приведены на рис. 1. Скачки в сечениях возникают при превышении энергией фотона потенциалов ионизации атомных оболочек. При увеличении растёт энергия связи электронов К-оболочки — для фотоэффекта на этой оболочке нужны фотоны всё больших энергий.

Релеевское или когерентное рассеяние — упругое рассеяние фотона на атоме без ионизации и возбуждения. Энергия рассеянного фотона равна энергии первичного фотона.

Зависимости сечения когерентного рассеяния R от энергии фотона и от атомного номера приведены на рис. 2.

Рис. 1. Зависимости сечения фотоэффекта ph от энергии фотона для разных и от атомного номера при разных. Видно, что вероятность фотоэффекта убывает с увеличением и возрастает с увеличением, за исключением скачков в области малых у веществ с большим, что связано с эффектами на соответствующих атомных оболочках.

Все сечения в данной работе взяты из [72].

Рис. 2. Зависимости сечения релеевского рассеяния R от энергии фотона и атомного номера. Видно, что вероятность когерентного рассеяния убывает с увеличением и возрастает с увеличением.

Эффект Комптона — некогерентное рассеяние фотона на атомном электроне. Энергия рассеянного фотона меньше энергии первичного фотона. Электрон, на котором произошло рассеяние, покидает атом. Сечение эффекта Комптона зависит от атомного номера вещества и энергии фотона как:

Зависимости сечения некогерентного рассеяния C от энергии фотона и от атомного номера приведены на рис. 3.

Рис. 3. Зависимости сечения комптоновского рассеяния C от энергии фотона и атомного номера. Видно что вероятность некогерентного рассеяния убывает с увеличением и возрастает с увеличением.

Образование электрон-позитронных пар может происходить в поле атомного ядра (np ) и в поле атомных электронов (ep ). В процессе образования пар первичный фотон поглощается и его энергия переходит в массу и кинетическую энергию позитрона и электрона.

Рождение пары возможно если энергия фотона больше 2. Сечение образования пар np зависит от атомного номера вещества и энергии фотона как:

Зависимости суммы сечений рождения пар в поле ядра и атомных электронов p от энергии фотона и от атомного номера приведены на рис. 4.

Под действием фотонов могут также происходить фотоядерные взаимодействия (nuc ), вероятность которых, однако, относительно невелика [73].

Полная вероятность взаимодействия фотона с атомом вещества является суммой вероятностей отдельных процессов:

Рис. 4. Зависимости сечения образования пар p от энергии фотона и атомного номера. Приведены зависимости суммы сечений рождения пар в поле ядра и атомных электронов p = np + ep. Видно что вероятность рождения пар возрастает с увеличением энергии фотона и атомного номера.

Рис. 5. Зависимости суммарного сечения tot от энергии фотона и атомного номера.

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

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

1.2. Ослабление потока фотонов 1.2.1. Монохроматические пучки фотонов При прохождении через вещество фотоны могут либо поглощаться (фотоэффект, рождение пар), либо рассеиваться (релеевское и комптоновское рассеяние). В результате этих процессов фотоны выбывают из пучка. Количество остающихся в пучке фотонов уменьшается [73, 74].

Рис. 6. Схема расположения источника фотонов, поглотителя, коллиматоров и детектора при измерении ослабления монохроматического пучка фотонов [73]. — толщина поглотителя, (0) — число фотонов в пучке до поглотителя, () — число фотонов, оставшихся в пучке после поглотителя.

Ослабление узкого параллельного монохроматического пучка фотонов при прохождении через поглотитель (рис. 6) описывается формулой:

где (0) — количество фотонов в пучке до поглотителя, () — количество фотонов, прошедших через поглотитель без взаимодействия, т.е. оставшихся в пучке, — толщина поглотителя, (, ) — здесь и далее, если не указано иное, под понимается tot — полное эффективное сечение взаимодействия фотонов с веществом поглотителя, — начальная энергия фотонов, — атомный номер вещества поглотителя, — количество атомов вещества в см3. Произведение величин и представляет собой количество атомов вещества приходящееся на см2 поперечного сечения поглотителя. В некоторых случаях это произведение удобно представлять в виде одной величины. Сечение в [барн/атом] зависит от энергии фотона и атомного номера и не зависит от плотности вещества. Количество атомов вещества в единице объёма удобно вычислять следующим образом:

где — плотность вещества поглотителя [г/см3 ], ат — масса одного атома [г], — атомная масса вещества поглотителя [а.е.м.]. Значения а.е.м. [г] и других констант взяты из [75].

Прозрачность поглотителя определяется как доля фотонов оставшихся в первичном пучке после прохождения поглотителя:

Из определения прозрачности следует, что эта безразмерная величина меняется в пределах от 0 до 1. Следовательно величина ослабления:

может меняться от 1 до, в пределе,. Величину ослабления также можно понимать как необходимое количество (0) фотонов до поглотителя для того, чтобы после поглотителя в среднем оставался один непровзаимодействовавший фотон, т.е. при () = 1.

Рис. 7. Энергетические зависимости ослабления узкого монохроматического пучка фотонов поглотителями разной толщины из воды, железа и свинца. Величины ослабления получены согласно (3).

На рис. 7 приведены энергетические зависимости величин ослабления узкого монохроматического пучка фотонов при прохождении через различные поглотители. Видно, что в области низких энергий величины ослабления существенно превышают значения при больших энергиях. Например, железный поглотитель толщиной 40 см ослабляет поток фотонов с энергией 10 МэВ в 104 раз, а с энергией 1 МэВ в 108 раз. Кроме того, чем больше толщина поглотителя, тем более выражено относительное отличие величин ослабления при разных энергиях фотонов.

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

Ослабление узкого коллимированного пучка фотонов с непрерывным энергетическим спектром при прохождении через слой поглотителя описывается соотношением (5), подобным (4):

где — энергия фотонов, (, = 0) — количество фотонов с энергией до поглотителя, () — суммарное количество фотонов прошедших без взаимодействия, (, ) — полное сечение взаимодействия фотонов с веществом поглотителя, — количество атомов вещества в единице объёма поглотителя, — толщина слоя поглотителя.

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

где (, ) — энергетический спектр фотонов тормозного излучения до просвечиваемого объекта, — верхняя граница спектра тормозного излучения, dep — энергия, выделяемая фотоном при взаимодействии с детектором, (, dep ) — отклик детектора — распределение вероятности того, что фотон с энергией выделит в детекторе энергию dep при попадании в чувствительную область детектора.

В таком случае, прозрачность объекта, при энергии ускорителя, определяется как отношение зарегистрированной дозы прошедшего через объект излучения (,, ) к дозе измеренной без объекта (,, 0):

где: - количество атомов вещества приходящееся на см2 поперечного сечения просвечиваемого объекта.

1.3. Принцип измерения Z и nD 1.3.1. Монохроматические пучки фотонов При просвечивании неизвестного объекта монохроматическим пучком фотонов с энергией можно измерить величину прозрачности объекта для фотонов с энергией. Из (2) следует что:

Видно, что одинаковой прозрачностью могут обладать объекты из разных веществ, с разными плотностями и толщинами. Необходимым условием равенства прозрачностей двух разных объектов 1 = 2 является равенство произведения (, ) · · :

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

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

Таким образом в левой части (8) остаётся произведение двух величин и. Величина сечения, в свою очередь, зависит от известной энергии фотонов и неизвестного значения атомного номера. Табличные данные по сечениям различных процессов доступны для от 1 до 99. Следовательно, несмотря на то, что значение неизвестно, можно получить ограничения на значения :

увеличением. Поэтому подобная оценка ограничений на значения является правомерной. Однако оценка (9) является единственной полезной информацией, которую можно извлечь из одного измеренного значения прозрачности. При этом одному измеренному значению будут соответствовать 99 пар значений и.

Оценить значения двух неизвестных величин и можно, сопоставляя значения двух прозрачностей 1 и 2, измеренных при двух различающихся энергиях фотонов 1 и 2. В случае пучков монохроматических фотонов удобно использовать отношение логарифмов прозрачностей (1, 2, ):

Это отношение зависит только от отношения сечений при двух энергиях 1 и 2 фотонов. Измерив экспериментально прозрачности 1 и 2 неизвестного объекта при энергиях фотонов 1 и 2, можно получить измеренное значение отношения Теоретическая зависимость отношения (10) от атомного номера может быть получена как отношение известных полных сечений 1 () и 2 () при энергиях фотонов 1 и 2 в виде Сопоставляя измеренное значение с заранее известной зависимостью (1, 2, ), можно оценить наиболее вероятное значение атомного номера неизвестного объекта.

Необходимым условием однозначности измерения являеется монотонность зависимости (). Иначе одному значению будут соответствовать два или более значения.

Более подробно, на примере изотопов 55 Cs и 27 Co, это вопрос рассмотрен в п. 2.1.2.

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

1.3.2. Тормозное излучение Как и в случае монохроматических фотонов, просвечивание пучком тормозного излучения с одной верхней границей позволяет получить только изображение прозрачности (, ) содержимого исследуемого объекта. На рис. 8 приведён пример того, что объекты из совершенно разных веществ могут иметь одинаковую прозрачность. Т.е. по одной прозрачности нельзя сделать никакого вывода об атомном номере вещества объекта.

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

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

В работе [40] предлагается просвечивать объекты пучками тормозного излучения с двумя верхними границами 8 и 4 МэВ. Перед началом измерений составляется таблица соотношений логарифмов прозрачностей:

с которой потом сравниваются значения отношений логарифмов измеренных прозрачностей. На основании сравнения делается вывод о наиболее вероятном значении неизвестного объекта. Следует заметить, что соотношение (12), в отличие от (10), зависит от В работе [44] верхние границы энергии пучков тормозного излучения предлагается сделать равными 9 и 5 МэВ. Создаётся таблица, содержащая соотношения энерговыделения ( = 9 МэВ) в детекторе после прохождения слоя поглотителя 9 МэВ спектром к ( = 5 МэВ) — энерговыделению 5 МэВ спектра, для различных веществ и разных толщин поглотителя.

имея измеренные ( = 9 МэВ) и ( = 5 МэВ), путём сравнения с табличными данными, можно оценить вещества поглотителя.

Однако, метод двух энергий имеет характерную особенность, заключающуюся в том, что в некоторых случаях объекта измеряется неоднозначно [40, 47].

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

1.4.1. Нормированные сечения Для того, чтобы оценить различия характерного поведения энергетических зависимостей сечений фотоэффекта () для разных веществ, разделим эти сечения на их значения при энергии фотонов 3.7 МэВ, т.е. проведём следующую нормировку:

Нормировка необходима, поскольку абсолютные значения сечений различных веществ отличаются. Кроме того, количество вещества в законе ослабления потока фотонов входит как масштабный коэффициент при сечении. Таким образом, эта нормировка позволяет исключить из рассмотрения масштабный фактор и сконцентрироваться на отличиях характерных особенностей сечений присущих различным веществам. Выбор энергии при Рис. 9. Энергетические зависимости сечений фотоэффекта, отнормированные на их значения при энергии фотона = 3.7 МэВ: (, ) = (, )/(3.7, ) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для железа: (, ) = = [(, )/(3.7, )]/[(, 26)/(3.7, 26))] (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для свинца: (, ) = [(, )/(3.7, ))]/[(, 82)/ /(3.7, 82)] (в). Видно, что максимальные отличия (при сделанных нормировках) достигают 25% при энергии 1 МэВ и 20% при энергии 10 МэВ. Наибольшее отличие наблюдается между зависимостями для = 1 и 99. Зависимости для больших отличаются друг от друга существенно слабее, чем зависимости для малых.

Рис. 10. Зависимости сечений фотоэффекта от атомного номера, отнормированные на их значения при = 26: (, )/(, 26) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/((, 26)/(3.7, 26)) (б).

Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82:

((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Видно, что при энергии 500 кэВ зависимости б) и в) не являются монотонными. Наибольшая крутизна зависимостей наблюдается при энергиях 1 МэВ и 10 МэВ, при этом зависимости можно считать монотонными.

Рис. 11. Энергетические зависимости сечений релеевского рассеяния, отнормированные на их значения при энергии фотона 3.7 МэВ: (, )/(3.7, ) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/ /((, 26)/(3.7, 26)) (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82: ((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Видно, что максимальные отличия (при сделанных нормировках) достигают 8% при энергии 0.5 МэВ и < 1% при энергии 10 МэВ. Наибольшее отличие наблюдается между зависимостями для = 1 и 99. Зависимости для больших отличаются друг от друга сильнее, чем зависимости для малых.

Рис. 12. Зависимости сечений релеевского рассеяния от атомного номера, отнормированные на их значения при = 26: (, )/(, 26) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/((, 26)/ /(3.7, 26)) (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82: ((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Наибольшая крутизна зависимостей наблюдается при энергиях 0.5 МэВ и 1 МэВ, при этом зависимости можно считать монотонными. Зависимости при энергиях 2.5 МэВ отличаются друг от друга очень слабо.

Рис. 13. Энергетические зависимости сечений комптон-эффекта, отнормированные на их значения при энергии фотона 3.7 МэВ: (, )/(3.7, ) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/((, 26)/ /(3.7, 26)) (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82: ((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Видно, что максимальные отличия (при сделанных нормировках) достигают 2% при энергии 0.5 МэВ и < 0.3% при энергии 10 МэВ. Наибольшее отличие наблюдается между зависимостями для = 1 и 99.

Рис. 14. Зависимости сечений комптон-эффекта от атомного номера, отнормированные на их значения при = 26: (, )/(, 26) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/((, 26)/ /(3.7, 26)) (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82: ((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Наибольшая крутизна зависимостей наблюдается при энергиях 0.5 МэВ и 1 МэВ, при этом зависимости можно считать монотонными. Зависимости при энергиях 2.5 МэВ отличаются друг от друга очень слабо.

Рис. 15. Энергетические зависимости сечений рождения пар, отнормированные на их значения при энергии фотона 3.7 МэВ: (, )/(3.7, ) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/((, 26)/ /(3.7, 26)) (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82: ((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Видно, что максимальные отличия, в зависимости от нормировки, достигают 80-140% при энергии 1.5 МэВ и 60-70% при энергии 10 МэВ. Наибольшее отличие наблюдается между зависимостями для = 1 и 99. Зависимости для больших отличаются друг от друга слабее, чем зависимости для малых.

Рис. 16. Зависимости сечений рождения пар от атомного номера, отнормированные на их значения при = 26: (, )/(, 26) (а). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 26: ((, )/(3.7, ))/((, 26)/(3.7, 26)) (б). Сечения, отнормированные на их значения при энергии 3.7 МэВ и на сечения для = 82: ((, )/(3.7, ))/((, 82)/(3.7, 82)) (в). Наибольшая крутизна зависимостей наблюдается при энергиях 1.1 МэВ и 2.5 МэВ, при этом зависимости можно считать монотонными. Зависимости при энергиях 5 МэВ отличаются друг от друга слабо, за исключением диапазона малых.

которой осуществляется нормировка равной 3.7 МэВ обусловлен тем, что до этой энергии полные сечения для всех убывают. При энергии фотонов большей чем 3.7 МэВ сечения для малых продолжают также убывать, а для больших, напротив, начинают возрастать. Т.е. при энергии 3.7 МэВ имеется характерная особенность — сечения, соответствующие веществам с большим имеют минимум при этой энергии. Результат проделанной нормировки показан на рис. 9 а. Видно что отличия незначительны.

Далее, для того чтобы оценить отличия характерного поведения зависимостей сечений фотоэффекта от атомного номера () при разных энергиях фотона, разделим эти зависимости на их значения при = 26:

Выбор значения атомного номера, при котором происходит нормировка равным 26 обусловлен тем, что железо является наиболее распространённым конструкционным материалом. Результат нормировки приведен на рис. 10 а. Такая нормировка показывает характерные отличия сечений для разных веществ от сечения для железа. Видно, что отличия также незначительны, за исключением области малых энергий.

Проведём двойную нормировку. Сначала разделим все сечения на их значения при энергии фотонов 3.7 МэВ. Затем разделим их на сечения для железа ( = 26):

Полученные зависимости (рис. 9 б и рис. 10 б) показывают отличия характерного поведения сечений взаимодействия для разных веществ от характерного поведения сечения для железа.

На рис. 9 в и рис. 10 в показаны сечения, отнормированные при энергии фотонов 3.7 МэВ приведённые по отношению к сечению свинца ( = 82):

Выбор свинца для нормировки обусловлен тем, что свинец является распространённым материалом с большим атомным номером. Видно, что характерное поведение сечений фотоэффекта для веществ с большим отличается слабо.

На рис. 11 и 12 показаны аналогичные зависимости нормированных сечений когерентного рассеяния, на рис. 13 и 14 — некогерентного рассеяния, на рис. 15 и 16 — рождения электрон-позитронных пар.

Из всех рассмотренных процессов наибольшие относительные отличия в диапазоне энергий фотонов от 500 кэВ до 10 МэВ наблюдаются между сечениями рождения пар и, при сделанных нормировках, достигают 140 210%. Несколько меньшие относительные отличия наблюдаются между сечениями фотоэффекта 45%. Максимальные относительные отличия сечений релеевского рассеяния существенно меньше 8%. Относительными отличиями сечений эффекта Комптона можно пренебречь из-за малой величины этих отличий 2%.

1.4.2. Вклады сечений основных процессов в полное сечение Рассмотрим совместно сечения основных процессов и их вклад в полное сечение. На рис. 17 приведены энергетические зависимости сечений для некоторых атомных номеров. На рис. 18 приведены зависимости этих же процессов от атомного номера при некоторых энергиях фотонов. Графики сечений выполнены в разных масштабах, т.к.

абсолютные значения сечений сильно отличаются у разных веществ. На рис. 19 и рис. приведены вклады различных процессов в полное сечение взаимодействия. Все графики вкладов выполнены в одинаковом масштабе для удобства сравнения.

При анализе зависимостей сечений от атомного номера можно заметить следующее:

• Для веществ с малым атомным номером основным процессом взаимодействия в рассматриваемом диапазоне энергий является эффект Комптона.

• С увеличением атомного номера в области низких энергий наибольший вклад в полное сечение начинает давать фотоэффект. Это приводит к тому, что энергетическая зависимость полного сечения у веществ с большим спадает быстрее чем у веществ с малым.

• В области больших энергий с увеличением атомного номера всё больший вклад в полное сечение даёт процесс рождения пар. В отличие от прочих процессов, вероятность рождения пар возрастает при увеличении энергии фотона. Эти два фактора приводят к тому, что энергетическая зависимость суммарного сечения в области больших энергий спадает всё медленнее при увеличении, а затем для больших начинает возрастать.

• В области энергий от 1 до 4 МэВ наибольший вклад в полное сечение взаимодействия продолжает давать эффект Комптона.

При рассмотрении зависимостей сечений от атомного номера можно заметить следующее:

• В области низких энергий у веществ с малым доминирующим процессом является эффект Комптона, а у веществ с большим — фотоэффект.

• В средней области энергий для всех основным процессом является некогерентное рассеяние.

• С ростом энергии у малых основным процессом взаимодействия остаётся некогерентное рассеяние, а для больших на первое место выходит образование электронпозитронных пар.

Рис. 17. Зависимости сечений основных процессов и полного сечения от энергии фотона для некоторых. Зависимости вкладов основных процессов в полное сечение взаимодействия от энергии фотона для некоторых, выполненные в одинаковом масштабе, приведены на рис. 19.

Рис. 18. Зависимости сечений основных процессов и полного сечения от атомного номера при некоторых энергиях фотона. Зависимости вкладов основных процессов в полное сечение взаимодействия от атомного номера при некоторых энергиях фотона, выполненные в одинаковом масштабе, приведены на рис. 20.

Рис. 19. Зависимости вкладов сечений основных процессов в полное сечение от энергии фотона для разных. Видно следующее:

• У веществ с малым атомным номером во всём рассматриваемом диапазоне энергий основным доминирующим эффектом является некогерентное рассеяние.

• При увеличении атомного номера до 26 всё более существенный вклад в области больших энергий даёт рождение пар, а в области малых энергий существенных изменений не наблюдается. Таким образом для того чтобы отличать между собой вещества с необходимо просвечивать объекты фотонами с энергиями 2 МэВ.

• Для веществ с 74 величина вклада рождения пар при 5 МэВ меняется очень слабо.

Рис. 20. Зависимости вкладов сечений основных процессов в полное сечение от атомного номера при некоторых энергиях фотона. Видно следующее:

• При низких энергиях у веществ с небольшим доминирует некогерентное рассеяние, а у веществ с большим доминирует фотоэффект.

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

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

1.4.3. Полные сечения Рассмотрим отличия характерного вида суммарных сечений различных веществ. Аналогично проделанному в п. 1.4.1, разделим зависимости полных сечений на значения этих сечений при энергии фотонов 3.7 МэВ (14). Графики полученных отношений приведены на рис. 21 а. Далее разделим зависимости сечений от на их значения при = 26 (15).

Полученные зависимости приведены на рис. 22 а. На рис. 21 б, в и рис. 22 б, в полные сечения приведены относительно сечений для железа (16) и свинца (17) соответственно.

Рис. 21. Энергетические зависимости полных сечений делённые на их значения при энергии фотона = 3.7 МэВ: () = ()/(3.7) (а). Сечения делённые на их значения при энергии 3.7 МэВ и на сечение для железа: (, ) = [(, )/(3.7, )]/[(, 26)/ /(3.7, 26)] (б). Сечения делённые на их значения при энергии 3.7 МэВ и на сечение для свинца: (, ) = [(, )/(3.7, )]/[(, 82)/(3.7, 82)] (в). Видно, что максимальные отличия, в рассмотренных случаях, достигают 60-100% при энергии фотонов 0.5 МэВ, 40% при энергии 1.5 МэВ и 60-80% при энергии 10 МэВ. Сечения у веществ с большим отличаются друг от друга значительно слабее, чем у веществ с малым. Причём в области энергий 2 МэВ сечения для больших отличаются друг от друга сильнее чем в области 2 МэВ.

1.4.4. Выводы В результате проведённого анализа сечений можно сделать следующие выводы:

• Рассматриваемую область энергий 0.5 10 МэВ можно условно разбить на несколько диапазонов, в которых характерное поведение зависимостей полных сечений для разных заметно отличается.

• В области энергий от 0.5 до 1.5 МэВ полные сечения для всех убывают. Причём, чем больше — тем быстрее убывает сечение, по причине возрастания вклада фотоэффекта при увеличении. Сечения для малых ведут себя одинаково т.к. вклад фотоэффекта Рис. 22. Зависимости полных сечений взаимодействия фотонов с веществом от атомного номера делённые на их значения при = 26: () = ()/(26) (а). Сечения делённые на их значения при энергии 3.7 МэВ и на сечения для железа: (, ) = [(, )/ /(3.7, )]/[(, 26)/(3.7, 26)] (б). Сечения делённые на их значения при энергии 3.7 МэВ и на сечения для свинца: (, ) = [(, )/(3.7, )]/[(, 82)/(3.7, 82)] (в). Видно, что в области малых наибольшая крутизна зависимостей наблюдается при больших энергиях 7.5 и 10 МэВ. В области больших наибольшая крутизна зависимостей наблюдается при энергиях 0.5 и 1 МэВ, а при энергиях 2.5 МэВ зависимости отличаются друг от друга слабо. Кроме того, видно, что полученные зависимости б) и в) не являются монотонными.

тона с масштабным коэффициентом.

• В области энергий от 1 до 3.7 МэВ основным процессом для всех является эффект Комптона. Наблюдается уменьшение влияния эффекта Комптона и возрастание влияния рождения пар. Убывание влияния эффекта Комптона сказывается, в основном, на сечениях для больших. Это приводит к уменьшению отличий поведений зависимостей полных сечений друг от друга для больших. Увеличение влияния рождения пар сказывается на сечениях для всех. Это приводит к тому, что энергетические зависимости сечений для малых начинают отличаться друг от друга.

• В области энергий 3.7 МэВ для больших основным процессом является рождение пар. Это имеет два следствия. Первое — сечения для больших возрастают в этой области энергий. Второе — характерное поведение полных сечений для больших является почти одинаковым. Влияние рождения пар на сечения для малых приводит к тому, что энергетические зависимости сечений убывают при увеличении по разному.

1.4.5. Рекомендации На основе выполненного исследования были выработаны следующие рекомендации:

• Для того, чтобы отличать друг от друга различные вещества необходимо как минимум две разных энергии фотонов.

• Для того, чтобы отличать друг от друга вещества с малым необходимо, чтобы одна из энергий была больше, чем 2 МэВ. Для того, чтобы отличать друг от друга вещества с большим необходимо, чтобы одна из энергий была меньше, чем 1.5 МэВ. Причём, наибольшее отличие наблюдается при энергии 0.5 МэВ. При энергии 1.5 МэВ характерное поведение полных сечений для всех оказывается практически одинаковым. Локальный минимум в сечениях для малых расположен при энергии 10 МэВ, а для больших при энергии 3.7 МэВ.

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

• Наиболее оптимальными значениями энергий фотонов можно считать 0.5, 1.5, 3.7 и 10 МэВ. Использование большего числа энергий нецелесообразно, т.к., при наличии ограничивающих факторов (суммарной средней мощности ускорителя, предельной активности источников, времени просвечивания, предельной мощности дозы и т.п.) затраты на промежуточные области энергии, уменьшат число фотонов в областях энергий, соответствующих наиболее характерным изменениям сечений.

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

1.4.6. Соотнешение вероятностей процессов рассеяния и поглощения фотонов Отдельно следует рассмотреть соотношение вероятностей процессов рассеяния и поглощения фотонов. Процессами рассеяния являются когерентное (релеевское) и некогерентное (эффект Комптона) рассеяние. К процессам поглощения фотонов, кроме фотоэффекта, можно причислить и эффект образования электрон-позитронной пары.

Построим зависимости от энергии фотона и атомного номера вещества величин вкладов сечений процессов рассеяния и поглощения в суммарные сечения (рис. 23), т.е.

отнормируем сечения рассеяния и поглощения на полное сечение:

Рис. 23. Вклады сечений процессов рассеяния и поглощения в суммарное сечение взаимодействия фотонов с веществом рассчитаны по формулам (18). Вклады процессов рассеяния изображены сплошными линиями. Вклады процессов поглощения изображены точками. Видно что при энергии 1.5 МэВ для всех, фотоны в основном будут рассеиваться.

А при энергиях 0.5 МэВ и 5 МэВ всё большее количество фотонов будет поглощаться.

Таким образом, если некие неизвестные объекты облучать фотонами с энергиями, например, 0.5 и 1.5 МэВ, то появляется возможность различать вещества, сравнивая количество поглощённых и рассеявшихся фотонов. Эта возможность довольно давно реализована для улучшения точности распознавания содержимого небольшой ручной клади и обнаружения запрещённых вложений в двойных стенках контейнеров, грузовиков и автомобилей [16]. Также, рядом авторов предлагается использование этой методики для обнаружения мин [18], [19]. Использование больших энергий представляется нецелесообразным, т.к. проникающая способность такого излучения больше, чем у фотонов с небольшой энергией и, следовательно, рассеявшиеся фотоны будут возникать на разной глубине что приведёт к неопределённостям и затруднениям при сопоставлении и анализе величин рассеянного излучения.

Глава 2.

Экспериментальные измерения атомного номера объектов 2.1. Радиоактивные источники Cs-137 и Co- 2.1.1. Выбор источников гамма-квантов При выборе изотопов важно учитывать следующие характеристики — число и значения энергий испускаемых -квантов, период полураспада, доступность изотопа, химические и механические свойства изотопа.

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

Предположим, что 1/2 = 1 год, а требуемый срок эксплуатации источника до плановой замены — 10 лет - тогда за эти 10 лет активность такого источника уменьшится в 103 раз, что вряд ли является приемлемым. Если период полураспада большой - то требуется большое количество изотопа. Следствием чего будет самопоглощение в самом источнике.

Например, для получения одной и той же активности, количества вещества изотопа с 1/2 = 100 лет потребуется в 10 раз больше, чем количества вещества изотопа с 1/2 = = 10 лет.

В табл. 1 приведены некоторые из широко применяемых и производимых в существенных количествах радиоактивных изотопов. Более полную информацию о доступности различных изотопов можно найти в [77,78]. Видно, что без применения -спектрометров возможно использование изотопов 55 Cs и 27 Co, а также в некоторой степени 56 Ba.

Таблица 1. Радиоактивные изотопы — кандидаты для использования в досмотровых комплексах. Приведены изотопы с периодами полураспада 3 100 лет, энергиями -квантов > 100 кэВ и вероятностями испускания -квантов > 15%. Характеристики изотопов взяты из [28, 76].

2.1.2. Соотношение R в случае Cs-137 и Co- где Cs — полное сечение взаимодействия с веществом фотонов с энергией = 662 кэВ, Co — среднее значение сечений при энергиях 1173 и 1333 кэВ, Можно считать, что энергии 1173 и 1333 кэВ расположены близко друг к другу, по сравнению с энергией 662 кэВ. На рис. 24а показана зависимость отношения (), полученная согласно (19). В области 30 значения величины почти не зависят от атомного номера. В области больших наблюдается заметное отличие значений для разных Однозначность и устойчивость измерения определяются отличием для разных веществ. Если нескольким разным соответствует одно и тоже значение, то нет никакой возможности выбрать из них наиболее вероятное. Возникает неоднозначность. На рис. 24б показано с какой точностью должна быть измерена величина для того, чтобы можно было различить двух образцов различающиеся по на 5 единиц. Видно, что при прочих равных условиях, вещества с большим различаются на два порядка лучше, чем вещества с малым.

На рис. 25 видно, что при используемых энергиях в области малых доминирует эффект Комптона. Сечение Комптоновского рассеяния линейно зависит от. Следствием этих двух причин является близость зависимости () к константе в области малых. С Рис. 24. a) — отношение () полных сечений взаимодействия фотонов с веществом при энергиях, соответствующих испускаемым источниками 55 Cs и 27 Co. Значения величины получены согласно (19). Значения сечений взяты из [72]. б) — для веществ, отличающихся на = 5.

увеличением всё более существенный вклад вносит фотоэффект и при энергии -квантов 662 кэВ начинает доминировать у веществ с 87. Влияние фотоэффекта проявляется в выраженном монотонном уменьшении значения отношения ().

Для того, чтобы уменьшить неопределённости в области малых необходимо использовать -источник с энергиями меньшими, чем у 55 Cs что позволит воспользоваться тем, что при меньших энергиях вклад фотоэффекта будет более существенным. Однако, полное сечение также будет расти с уменьшением энергии -квантов. Это потребует использования источника с большей активностью чем, у источника 137 Cs. Однако, источник большей активности неизбежно имеет большую массу и размер. При достаточно больших размерах начинает сказываться самопоглощение и рассеяние в самом источнике. Количество вылетающих из источника -квантов перестаёт линейно зависеть от количества вещества источника, кроме того, источник начинает греться.

2.1.3. Описание эксперимента Для оценки возможности применения радиоактивных источников для измерения атомного номера неизвестных объектов был проведён эксперимент [53] по просвечиванию тестовых образцов пучками -квантов. Использовались калиброванные радиоактивные источники: 137 Cs с активностью 92 кБк испускающий -кванты с одной энергией 662 кэВ Рис. 25. Зависимости сечений фотоэффекта и эффекта Комптона от атомного номера для -квантов испускаемых, 55 Cs и 27 Co. Значения сечений взяты из [72]. Сплошная и пунктирная линии — сечения эффекта Комптона и фотоэффекта при энергии 662 кэВ, штриховая — среднее значение сечений эффекта Комптона при энергиях и 1333 кэВ, штрих-пунктирная — среднее значение сечений фотоэффекта при энергиях 1173 и 1333 кэВ.

и 27 Co с активностью 39 кБк испускающий за один акт распада два -кванта с энергиями 1173 и 1333 кэВ. В ходе эксперимента просвечивались образцы разной толщины из разных веществ. Параметры образцов приведены в табл. 2.

Таблица 2. Образцы использовавшиеся в эксперименте. - атомный номер, - толщина.

Для измерений в качестве детектор -квантов использовался -спектрометр из сверхчистого германия Canberra GC3019 [79]. С целью уменьшения фона детектор располагается в свинцовой защите показанной на рис. 26 в которой предусмотрена камера для исследуемых образцов. Внутрь этой камеры (рис. 28) были помещены держатель для калиброванных источников и свинцовый коллиматор. Схема расположения показана на рис. 27.

Схема устройства детектора показана на рис. 29. Для каждого источника проводились изРис. 26. Детектор -излучения Canberra GC3019 с криостатом в свинцовой защите, установленный в НИИЯФ МГУ.

мерения без образцов для получения величины (0) и с каждым образцом для получения величин (). Время одного измерения составляло 10 мин.

Рис. 27. Схема расположения источника -квантов, свинцового коллиматора, просвечиваемого объекта и детектора -излучения.

2.1.4. Результаты эксперимента и выводы В результате проведённой серии измерений был получен набор спектров -квантов, зарегистрированных детектором. Анализ спектров, определение площадей пиков полного поглощения и подложки проводились на основе методики и с помощью программного Рис. 28. На фотографии изображены свинцовый коллиматор, просвечиваемый образец и детектор -квантов Canberra GC3019 расположенные в измерительной камере (вид сверху).

Рис. 29. Схема устройства детектора -излучения Canberra GC3019 [79]. Приведены оценённые значения характерных размеров кристалла, + контакта, конструкционных элементов. Производитель приводит далеко не все важные параметры детектора. Оценка параметров проводилась в процессе калибровки эффективности детектора [62].

обеспечения, разработанных в НИИЯФ МГУ [54, 55, 65]. На рис. 30 приведены спектры, полученные на источниках 55 Cs и 27 Co в отсутствие образцов.

Рис. 30. Спектр зарегистрированных энергий -квантов провзаимодействовавших с детектором. Источники — 55 Cs и 27 Co. dep — зарегистрированное значение энергии выделившейся в детекторе. — число событий.

Рассмотрим процесс формирования спектра зарегистированных энергий -квантов [80].

При попадании -кванта в область детектора возможны три варианта:

• -квант проходит сквозь детектор, не провзаимодействовав с ним никакого сигнала детектор не формирует, • -квант взаимодействует с детектором, вся энергия -кванта и вторичных частиц передаётся на ионизацию чувствительной области детектора (dep = ) формируется сигнал, соответствующий пику полного поглощения на единицу увеличивается площадь пика полного поглощения, • -квант взаимодействует с детектором, однако он рассеявшись и/или вторичные частицы покидают детектор, унося часть энергии (0 < dep < ) формируется сигнал, дающий вклад в подложку спектра на единицу увеличивается площадь подложки.

При расчёте прозрачности объекта учитывались только площади пиков полного поглощения. Прозрачность для источника 55 Cs вычисляется как Cs = Cs ()/Cs (0). Для источника 27 Co учитывалась сумма обоих пиков Co = (1173 кэВ) + (1333 кэВ). Тогда Co = Co ()/Co (0). Без образцов количество регистрируемых за 10 мин -квантов в пике полного поглощения составляло 8.5 · 104 для источника 137 Cs и 6.5 · 104 для источника 27 Co.

Экспериментальное значение отношения логарифмов измеренных прозрачностей для источников 55 Cs и 27 Co можно получить следующим образом Рис. 31. Отношения логарифмов экспериментально измеренных прозрачностей объектов (табл. 2), полученных согласно (20): — оргстекло, — алюминий и дюраль, — железо, — кадмий, — свинец. Сплошная линия — зависимость отношения полных сечений взаимодействия от атомного номера, полученная согласно (19). Значения сечений взяты из [72].

Рис. 32. Измеренные значения атомного номера образцов из свинца и кадмия.

- известная толщина тестовых образцов. Штриховая линия — ожидаемое значение для свинца, пунктирная — для кадмия.

На рис. 31 показаны отношения логарифмов измеренных прозрачностей тестовых объектов в сравнении с зависимостью отношений сечений от атомного номера. Видно что Рис. 33. Ожидаемые и измеренные значения количества вещества образцов из свинца и кадмия. Штриховая линия — соответствует случаю msr = real.

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

Можно только сделать вывод что 50. Значения для образцов из свинца заметно отличается от значений для остальных образцов. Можно уверенно сделать вывод что 83. В целом наблюдается согласие измеренных значений с теоретической зависимостью.

Сопоставив измеренные с известной зависимостью, можно оценить наиболее вероятные значения атомных номеров. Затем, по измеренным прозрачностям и атомным номерам объектов согласно (11), можно оценить количество вещества. Соответствующие результаты для образцов из кадмия и свинца приведены на рис. 32 и 33.

Хорошее согласие измеренных и ожидаемых значений атомного номера и количества вещества свидетельствует о потенциальной возможности, при использовании изотопов Cs-137 и Co-60, относительно просто измерять атомный номер и количество вещества объектов с > 40. Для определения атомного номера объектов с < 40 необходимо существенно увеличить активность источников.

2.2. Тормозное излучение 2.2.1. Описание эксперимента Общая схема эксперимента [52] и расположение основных элементов показана на рис. 34. Линейный ускоритель электронов [81] работает в режиме импульсного переключения энергии пучка. Длительность импульса составляет 3 мкс, расстояние между импульсами 10 мс. Значение энергии электронного пучка может изменяться от 3 до 8 МэВ.

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

Рис. 34. Общая схема эксперимента и расположение основных элементов установки. У - электронный ускоритель, Э - пучок ускоренных электронов, М - тормозная мишень из вольфрама, К1 - свинцовый коллиматор с размерами 380.52 см3, - пучок фотонов, О исследуемые образцы на подвижной платформе, К2 - свинцовый коллиматор с размерами 82 2 2.5 см3, Д - детектор, З - свинцовая защита детектора, - расстояние между тормозной мишенью и детектором 768 ± 1см.

При помощи коллиматора, установленного сразу после тормозной мишени, вырезается тонкий пучок. За коллиматором расположена подвижная платформа с закреплёнными на ней тестовыми образцами из различных материалов (рис. 35).

Для эксперимента были изготовлены тестовые образцы из разных материалов (табл. 3).

Поперечный размер всех образцов был одинаковым и составлял 4 4 см.

На расстоянии 768 см от тормозной мишени располагался детектор излучения, прошедшего через образцы, с помощью которого измерялась суммарная доза в мЗв за 10 с. В качестве детектора использовался “Блок детектирования БДКС-96 (дозиметр - радиометр ДКС-96г)” производства ООО НПП “Доза” [82]. Детектирующий элемент - тканеэквивалентный пластмассовый сцинтиллятор, 45 20 мм. Диапазон энергии регистрируемого излучения 0,015 10 МэВ. Для уменьшения влияния рассеянного излучения, перед детектором расположен ещё один коллиматор, а сам детектор окружён свинцовой защитой.

Рис. 35. Тестовые образцы на подвижной платформе. Al - эталонный образец из алюминия толщиной 20 см, О - “неизвестные” образцы, П - подвижная платформа, В - видеокамера для контроля положения образцов, К1 - выходное отверстие коллиматора.

Таблица 3. Таблица. Тестовые образцы. - атомный номер, - толщина образца.

Измерения проводились при четырёх разных значениях энергии ускорителя 7.3, 6.2, 5.5, 3.5 МэВ. При оценке эффективного атомного номера тестового образца msr учитывались все четыре значения прозрачностей полученные при четырёх разных энергиях ускорителя.

Контроль стабильности дозы осуществлялся при отсутствии образцов на пути пучка.

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

Предварительно, с использованием пакета GEANT4 [83], было выполнено моделирование характеристик тормозного излучения при разных энергиях электронного пучка.

Также было выполнено моделирование отклика детектора. Результаты моделирования использовались при построения таблицы прозрачностей (,, ), со значениями которой проводится сравнение [40] экспериментально измеренных прозрачностей и для построения зависимости прозрачности образца из алюминия толщиной 20 см от значения верхней границы тормозного спектра (рис. 36), при помощи которой, по измеренному значению прозрачности, методом ослабления [84], оценивалось значение энергии ускорителя.

Рис. 36. Зависимость прозрачности поглотителя из алюминия толщиной 20 см от энергии ускорителя.

2.2.2. Результаты эксперимента и выводы В результате эксперимента для каждого объекта при четырёх фиксированных энергиях ускорителя измерялось значение дозы излучения, прошедшего через исследуемый объект. Далее вычислялась прозрачность объекта (,, ) = (,, )/(,, 0).

Затем набор прозрачностей, полученных при всех четырёх энергиях, используя МНКподобный алгоритм сопоставлялся, с предварительно полученной таблицей прозрачностей, в результате чего получалось наиболее вероятное значение атомного номера исследуемого объекта. На рис. 37 приведены результаты измерения атомного номера тестовых объектов при различных значениях прозрачности, рассчитанной для энергии фотонов 3.7 МэВ.

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

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

Рис. 37. Измеренные значения атомного номера и ожидаемые значения (прямые линии) тестовых образцов: — свинец, = 82; — кадмий, = 48; — чугун, 26;

• — дюраль, 13; — оргстекло, 4 при различных значениях прозрачности, рассчитанной для энергии фотонов 3.7 МэВ. Параметры образцов приведены в табл. 3.

Видно, что измеренные значения в целом совпадают с ожидаемыми.

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

Глава 3.

Моделирование на GEANT Моделирование тормозного излучения, прохождения частиц через вещество и отклика детектора проводились с помощью GEANT4 [83, 85], представляющего собой пакет программ для моделирования прохождения частиц через вещество на основе методов МонтеКарло. Включает в себя инструменты для описания геометрии, используемых веществ и материалов, силовых полей. Содержит множество физических моделей процессов взаимодействия частиц с веществом.

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

Перечислим некоторые характеристики грузовых вагонов и контейнеров. Максимальная длина у наиболее распространённых железнодорожных вагонов не превышает 25 м.

Максимальная высота < 5.4 м, ширина < 3.6 м. У контейнеров-цистерн толщина стенки из нержавеющей стали может достигать 25 мм. Соответственно суммарная толщина стальных стенок может достигать 5 см.

Типовая схема установки изображена на рис. 38. Пучок электронов из ускорителя попадает на тормозную мишень из вольфрама, в которой происходит генерация тормозного излучения. С помощью щелевого коллиматора формируются тонкий веерообразный пучок тормозного излучения. Далее этот пучок проходит через исследуемое транспортное средство и попадает на детектор. В досмотровых комплексах линейка детекторов обычно имеет Г-образную форму из-за простоты изготовления и размещения. В медицинских системах линейка детекторов имеет обычно круговую форму. В данной работе, для упрощения расчётов при моделировании, форма детекторной линейки выбрана в виде дуги окружности с углом 60°и с центром, совпадающим с тормозной мишенью ускорителя.

Сигнал от детектора идёт на ЭВМ для обработки и анализа.

Рис. 38. Схема установки для интроскопии крупногабаритных объектов, используемая для моделирования.

3.2. Тормозное излучение 3.2.1. Взаимодействие электронов с веществом Как известно, при взаимодействии с веществом электронов с энергиями ( 1 МэВ) основными являются ионизационные потери [69, 70] — возбуждение и ионизация атомов.