Определение атомного номера вещества объектов по ослаблению пучков фотонов с энергиями до 10 МэВ тема автореферата и диссертации по физике, 01.04.16 ВАК РФ

Курилик, Александр Сергеевич АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2014 ГОД ЗАЩИТЫ
   
01.04.16 КОД ВАК РФ
Диссертация по физике на тему «Определение атомного номера вещества объектов по ослаблению пучков фотонов с энергиями до 10 МэВ»
 
Автореферат диссертации на тему "Определение атомного номера вещества объектов по ослаблению пучков фотонов с энергиями до 10 МэВ"

На правах рукописи УДК 539.12.04

Курилик Александр Сергеевич

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

Специальность 01.04.16 физика атомного ядра и элементарных частиц

АВТОРЕФЕРАТ диссертации на соискание учёной степени кандидата физико-математических паук

1 7 ИЮЛ 2014 -И И!0Л 2014 ■

Москва — 2014

005550525

005550525

Работа выполнена на кафедре общей ядерной физики физического факультета Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Московский государственный университет имени М.В. Ломоносова» и в Отделе электромагнитных процессов и взаимодействий атомных ядер Научно-исследовательского института ядерной физики имени Д.В. Скобельцына федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Московский государственный университет имени М.В. Ломоносова».

Научный руководитель:

Ишханов Борис Саркисович, заведующий кафедрой общей ядерной физики физического факультета МГУ имени М.В. Ломоносова, заведующий Отделом электромагнитных процессов н взаимодействий атомных ядер НИИЯФ имени Д.В. Скобельцына МГУ имени М.В. Ломоносова, доктор физико-математических наук, профессор.

Официальные оппоненты:

Гуров Юрий Борисович, доктор физико-математических наук, главный научный сотрудник кафедры «Физика элементарных частиц» Национального исследовательского ядерного университета «МИФИ».

Огородников Сергей Анатольевич, кандидат технических наук, директор по науке и технологиям ООО «Скантроник Системе».

Ведущая организация:

Федеральное государственное учреждение науки Институт ядерных исследований РАН г. Москва.

Защита состоится <10 » Оус^ХОтрЦ._2014 г. в 15 часов на заседании диссертационного совета Д 501.001.77 на базе Московского государственного университета имени М. В. Ломоносова по адресу: 119991, Москва, Ленинские горы, дом 1, строение 5, НИИЯФ МГУ, «19 корпус», ауд. 2-15.

С диссертацией можно ознакомиться в Фундаментальной библиотеке Московского государственного университета имени М.В. Ломоносова (Ломоносовский проспект 27).

Автореферат разослан «_о4_» ^/у-дьО!._2014 г.

Ученый секретарь Страхова С.И.

совета по защите диссертаций

на соискание ученой степени доктора наук, на соискание ученой степени кандидата наук Д 501.001.77,

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

Общая характеристика работы

Актуальность темы

Часто существует необходимость получить информацию о внутренней структуре и составе некоторого объекта. Например, обнаружить скрытые дефекты в продукции промышленного производства, проверить соответствие декларации и реального содержимого груза. Также важной задачей является контроль транспортных средств с целью выявления и пресечения незаконных перевозок — контрабанды. В последние годы особую актуальность приобрела проблема обеспечения безопасности и предотвращение терроризма [1,2]. Одним из основных элементов системы обеспечения безопасности являются технические средства, позволяющие осуществлять неразрушающпй контроль — исследовать, не вскрывая, грузы и транспортные средства.

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

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

Просвечивание излучением с двумя граничными энергиями позволяет восстановить распределение атомного номера Х{х,у) и количества вещества пО{х,у). Возможность измерения атомного номера 2 обусловлена различной энергетической зависимостью полных сечений взаимодействия фотонов с атомами различных химических элементов. Измерив прозрачности при двух разных энергиях и сопоставив с известными заранее зависимостями, можно сделать вывод о Z вещества исследуемого объекта. Т.к. в большинстве случаев материалы, из которых изготавливаются различные объекты, не состоят из атомов одного химического элемента с одним значением Z, а состоят из атомов различных химических элементов — вводится понятие эффективного атомного номера 7,ец.

7-источники давно и успешно применяются для дефектоскопии. Существуют как проекты досмотровых комплексов на основе источников 'ЩСэ или ®?Со [6,7], так и уже действующие 7-сканеры [8]. В нефтяной промышленности для контроля параметров смеси газа, нефти и воды применяется просвечивание 7-квантами от двух разных радиоактивных изотопов [9]. Обычно применяются изотопы, дающие 7-кванты небольших энергий. Рассматривается возможность использования изотопов ^Сэ и 27С0 для измерения Z объектов, расположенных в крупногабаритных контейнерах [10,11].

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

ных источников из-за небольших энергий и малых интенсивностей излучения. Для того, чтобы исследовать массивные объекты, используют ускорители электронов с тормозной мишенью, которые позволяют получать интенсивные потоки тормозных фотонов больших энергий. Существуют системы, просвечивающие досматриваемый объём и дающие в качестве результата двумерное изображение прозрачности К(х,у) [12-14]. Кроме стационарных, разработаны мобильные комплексы [15].

В работах С. Огородникова и В. Пструнина [16-20] был предложен метод распознавания веществ, основанный на анализе ослабления высокоэнергетичного тормозного излучения с двумя граничными энергиями, позволяющий различать вещества по эффективному атомному номеру и получать двумерные изображения Z{x,y) и пО(х,у). Реализован полностью функционирующий досмотровый комплекс, основанный на этих работах [21]. Проекты аналогичных систем описаны рядом авторов в работах [22-25]. Типичые значения граничных энергий, в предлагаемых проектах, составляют 6-9 и 3-6 для большей и меньшей энергий соответственно. Метод трёх энергий, упомянутый в работе С. Огородникова и В. Петрунина [16] и описанный в работе [26], является следующим шагом в интроскопии, позволяющим улучшить точность распознавания атомного номера неизвестных объектов с большим атомным номером. Также, в [27] описан проект системы с четырьмя ускорителями разных энергий, а в [28] — на основе одного ускорителя с переключением четырёх энергий.

Цели диссертационной работы

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

Основные результаты, полученные в диссертации

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

2) Впервые показано, что в области энергий фотонов от 0.5 до 10 МэВ наиболее характерные особенности и отличия полных сечений наблюдаются при энергиях 0.5, 1.5, 3.7 и 10 МэВ. Следовательно, для определения атомного номера 2 вещества объектов, посредством сравнения величин ослабления пучков монохроматических фотонов, целесообразно использование не более четырёх значений энергии фотонов.

Показана необходимость использования не менее трёх энергий фотонов для однозначного определения 2.

3) Показано, что при использовании радиоактивных изотопов Се-137 и Со-60 в качестве источников монохроматических фотонов, отношение логарифмов прозрачности! для различения веществ с^<5 необходимо измерять в 100 — 1000 раз точнее, чем для различения веществ с 2 > 40.

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

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

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

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

Научная новизна

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

ного излучения с тремя граничными энергиями позволяет повысить точность измерения атомного номера вещества объектов.

Ценность научной работы

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

Практическая значимость

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

Аппробация работы

Статьи в реферируемых журнале«

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

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

• Измерения атомного номера неизвестных объектов на ускорителе электронов с перестраиваемой энергией / А.Н. Ермаков, Б.С. Ишханов, A.C. Курилик и др. // Известия РАН. Серия физическая. - 2012. — Т. 76, No 11. — С. 1337-1339.

• Измерения атомного номера и количества вещества объектов с помощью радиоактивных источников / С.С. Белышев, Б.С. Ишханов, A.C. Курплик, К.А. Стопани // Известия РАН. Серия физическая. - 2013. - Т. 77, No 7. - С. 1008-1011.

• Определение выходов ядерных реакций на основе анализа цепочек распадов / С.С. Белышев, К.А. Стопани, С.Ю. Трощиев и др. // ВМУ. Серия 3. Физика. Астрономия. - 2011. - No 4. - С. 42-47.

• Studying photonuclear reactions using the activation technique / S. Belyshev, A. Erma-kov, B. Ishkhanov et al. // Nuclear Instruments and Methods in Physics Research Section A. - 2014. - Vol. 745. - P. 133-137.

Доклады на конференциях

Материалы диссертации доложены лично автором и опубликованы в трудах следующих конференций и совещаний:

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

• Курплик А., Руденко Д., Стопани К. Много-пучковый метод сканирования объектов // Труды VII межвузовской научной школы молодых специалистов "Концентрированные потоки энергии в космической технике, электронике, экологии и медицине". - 2006. - С. 83-88.

• Многопучковый метод определения атомного номера / Б.С. Ишханов, A.C. Курилик, Д.С. Руденко и др. // Труды VIII межвузовской научной школы молодых специалистов "Концентрированные потоки энергии в космической технике, электронике, экологии и медицине". — Издательство МГУ Москва, 2007. — С. 160-164.

• Kurilik A., Rudenko D., Stopani К. 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.

• Определение атомного номера методом многопучкового сканирования объектов / С.С. Белышев, A.A. Кузнецов, A.C. Курилик, К.А. Стопани // Тезисы докладов 58 международного совещания по ядерной спектроскопии и структуре атомного ядра "ЯДРО 2008". - Saint-Petersburg, 2008. - С. 281.

• Многопучковый метод определения структуры и состава объектов / С.С. Белышев, A.A. Кузнецов, A.C. Курилик, К.А. Стопани // Сборник тезисов, материалы Четырнадцатой Всероссийской научной конференции студентов-физиков и молодых учёных "ВНКСФ-14". - Уфа, 2008. - С. 285.

• Кузнецов А., Куршшк А. Калибровка эффективности HPGe детектора по экспериментальным измерениям и моделированию Geant4 // Тезисы докладов 59 международного совещания по ядерной спектроскопии и структуре атомного ядра "Ядро-2009". - СПбГУ Санкт-Петербург, 2009. — С. 314.

• 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 Fcmto- and Nanotechnologies>, July 6-9, 2010, Saint-Petersburg, Russia. Book of Abstracts. — Saint-Petersburg, 2010.

- P. 418.

• 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.

Также материалы диссертации доложены и опубликованы в трудах следующего совещания:

• Автоматизация обработки данных гамма-активациопных экспериментов / С.С. Бе-лышев, А.А. Кузнецов, А.С. Курилик, К.А. Стопани // Тезисы докладов 58 международного совещания по ядерной спектроскопии и структуре атомного ядра "Ядро-2008". - Санкт-Петербург, 2008. - С. 282.

Материалы диссертации докладывались лично автором на научных конференциях "Ломоносовские чтения "в 200G, 2007, 2009, 2010, 2011 и 2013 годах, а также на выставке-конкурсе "Инновационный проект 2008".

Публикации

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

Личный вклад автора

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

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

Структура и объём диссертации

Диссертация состоит из введения, четырех глав и заключения. Полный объем диссертации 111 страниц текста с 70 рисунками и 12 таблицами. Список литературы содержит 106 наименований.

Содержание работы

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

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

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

ТСБ -рл -й

_1_1_|_|_|_I_

12 3 4 5 6 7 Энергия фотонов Е [МэВ|

12 3 4 5 6 7 Энергия фотонов Е [МэВ]

Рис. 1. Зависимости сечений основных процессов и полного сечения от энергии фотона Е для некоторых 2. Все сечения в работе взяты из [31].

зависимости сечений для некоторых атомных номеров 2. Зависимость сечения фотоэффекта от атомного номера вещества 2 и энергии фотона Е, имеет вид: арь ~ 2Ъ¡Е7!2. Сечение эффекта Комптона зависит от атомного номера 2 вещества и энергии фотона Е как: ас ~ 2/Е. Сечение образования пар <тпр зависит от атомного номера 2 вещества и энергии фотона Е как: егпр ~ 22 ■ 1п(.Е).

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

При прохождении через вещество фотоны выбывают из первичного пучка [32]. Ослабление узкого параллельного монохроматического пучка (п. 1.2.1) фотонов при прохождении через поглотитель описывается формулой:

ТУр) = N(0) ■ е~а{Е-г)""

(1)

где N(0) — количество фотонов в пучке до поглотителя, N(0) — количество фотонов, прошедших через поглотитель без взаимодействия, т.е. оставшихся в пучке, В — толщина поглотителя, <т(Е, 2) — здесь и далее, если не указано иное, под а понимается аы — полное эффективное сечение взаимодействия фотонов с веществом поглотителя, Е — начальная энергия фотонов, 2 — атомный номер вещества поглотителя, п — количество атомов вещества в см3.

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

(2)

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

ЛГ(£>) = ]Г ЛГ(Я{, о = 0) • е—

(3)

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

ос

N(0) = J ЛГ(£, О = 0) ■ ° ¿Е, (4)

о

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

С(Е>) = У УIУ(Т, Е) ■ п ° ■ ЩЕ, ЕЛер) ■ ЕЛер АЕ <Щ,ер, (5)

о о

где \У{Т, Е) — энергетический спектр фотонов тормозного излучения до просвечиваемого объекта, Т — верхняя граница спектра тормозного излучения, Е^ер — энергия, выделяемая фотоном при взаимодействии с детектором, ЩЕ, р) — отклик детектора — распределение вероятности того, что фотон с энергией Е выделит в детекторе энергию Едер при попадании в чувствительную область детектора.

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

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

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

а(Е,2) -п- £> = -Ь(А'). (7)

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

В случае пучков монохроматических фотонов удобно использовать отношение логарифмов прозрачностсй ЩЕ\, Е2,2\.

ШР Р ^ _ -а{Е2,2)-п-Р а(Е2,2)

П{Еи 2) - ЩкГ) - -а(Еи2)-п-0 - (8)

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

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

а(Еъ2) а(Е2,2) (У>

Как и в случае монохроматических фотонов, просвечивание пучком тормозного излучения (п. 1.3.2) с одной верхней границей Т позволяет получить только изображение прозрачности К(х, у) содержимого исследуемого объекта. В случае использования пучков с двумя разными верхними границами энергий появляется возможность разделять вещества по эффективному атомному номеру 2. В работе [19] предлагается просвечивать объекты пучками тормозного излучения с двумя верхними границами 8 и 4 МэВ. В работе [23] верхние границы энергии пучков тормозного излучения предлагается выбрать равными 9 и 5 МэВ. Однако, при использовании двух энергий, в некоторых случаях 2 объекта измеряется неоднозначно [19,20].

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

<т'(Е, 2) = а(Е, 2)/а(3.7, г). (10)

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

Далее эти зависимости были разделены на их значения при 2 = 82 (рис. 2 б). , <т(Е, 2)/а(3.7,2)

а{Е>2) = а(Е,82)/а(3.7,82У (И)

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

г = 99 -г - 92 -г = 82 -

2 = 74 -г = 43 2 = 26

3 4 5 6 7 Энергия фотонов Е [МэВ]

1.4 1.3 1.2 1.1

8.6 9.5

—1-1 1 1-1 1 к 1-1........ 1

__^

V \

1 = 99 - Х->-ч

г = 92-

2 = 82 - \\

•г - 74 - Х

1 = 48 ___

•г = 26

г = 13 ..............

2 = ! - 1 • «

Энергия фотонов Е (МэВ] б)

Рис. 2. Энсргстические зависимости полных сечений делённые на их значения при энергии фотона Е = 3.7 МэВ (а). Сечения делённые на их значения при энергии 3.7 МэВ и на сечение для свинца (6). Все сечения в работе взяты из [31).

В п. 1.4.2 рассматриваются вклады сечений отдельных процессов в полное течение и зависимости полного сечения от Е и 2. В п. 1.4.3 анализируются зависимости полных сечений от Е и 2.

В п. 1.4.4 и 1.4.5 сформулированы выводы и рекомендации. В области энергии от 0.5 до ~1.5 МэВ полные сечения доя всех 2 убывают. Причём, чем больше 2 — тем быстрее убывает сечснис, по причине возрастания вклада фотоэффекта при увеличении 2. Сечения для малых 2 ведут себя одинаково т.к. вклад фотоэффекта и рождения пар очень мал и полные сечения представляют собой сечение эффекта Комптона с масштабным коэффициентом ~ 2. В области энергий от ~1 до ~3.7 МэВ основным процессом для всех 2 является эффект Комптона. Наблюдается уменьшение влияния эффекта Комптона и возрастание влияния рождения пар. Убывание влияния эффекта Комптона сказывается, в основном, на сечениях для больших 2. Это приводит к уменьшению отличий поведений зависимостей полных сечений друг от друга для больших 2. Увеличение влияния рождения пар сказывается на сечениях для всех 2. Это приводит к тому, что энергетические зависимости сечений для малых 2 начинают отличаться друг от друга. В области энергий > 3.7 МэВ для больших 2 основным процессом является рождение пар. Это имеет два следствия. Первое — ссчсния для больших 2 возрастают в этой области энергий. Второе — характерное поведение полных сечении для больших 2 является почти одинаковым. Влияние рождения пар на ссчсния доя малых 2 приводит к тому, что энсргстические зависимости сечений убывают при увеличении Е по разному. Для того, чтобы отличать друг от друга различные вещества необходимо как минимум две разных энергии фотонов. Для того, чтобы отличать друг от друга вещества с малым 2 необходимо, чтобы одна из энергий была больше, чем ~2 МэВ. Для того, чтобы отличать друг от друга вещества с большим 2 необходимо, чтобы одна из энергий была меньше, чем ~1.5 МэВ. Причём, наибольшее отличие наблюдается при энергии 0.5 МэВ, что является нижней границей

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

Во второй главе описываются проведённые автором эксперименты по просвечиванию тестовых объектов фотонами от радиоактивных источников (раздел 2.1) и от ускорителя электронов с перестраиваемой энергией (раздел 2.2).

В п. 2.1.1 обосновывается выбор изотопов '||Cs и 27С0 в качестве источников гамма-квантов. При этом учитывались следующие характеристики — число и значения энергий испускаемых 7-квантов, период полураспада, доступность изотопа, химические и механические свойства изотопа. При энергии 7-квантов меньше нескольких сотен кэВ — требуется большое количество такого изотопа т.к. полное сечение в этой области энергий возрастает при уменьшении энергии фотона Е. Если в процессе распада испускаются 7-кванты нескольких энергий, причём в разных энергетических диапазонах — теряется одно из преимуществ радиоактивных источников — моноэнергетичпость испускаемых 7-квантов. Применение такого источника возможно при использовании 7-спектрометра в качестве детектора. При расмотрении характеристик широко применяемых и производимых в существенных количествах радиоактивных изотопов, сделан вывод о том, что без применения 7-спектрометров возможно использование изотопов 'j^Cs и 27С0, а также в некоторой степени 'ffBa.

В п. 2.1.2 анализируется зависимость отношений сечений от Z при энергиях гамма-квантов, характерных для выбранных источников. Если в качестве источников 7-квантов используются изотопы 1||Cs и ™Со, тогда

где oes — полное сечение взаимодействия с веществом фотонов с энергией Е = 662 кэВ, Осо — среднее значение сечений при энергиях 1173 и 1333 кэВ.

На рис. За показана зависимость отношения R(Z), полученная согласно (12). В области Z < 30 значения величины R почти не зависят от атомного номера Z. В области больших Z наблюдается заметное отличие значений R для разных Z.

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

R

Z

Рис. 3. а) — отношение Д(^) полных сечений взаимодействия фотонов с веществом при энергиях, соответствующих испускаемым источниками '^Сэ и ¡¡"Со. б) — ЛИ для веществ, отличающихся на = 5.

г

Рис. 4. Сплошная и пунктирная линии — сечения эффекта Комптона и фотоэффекта при энергии 662 кэВ, штриховая — среднее значение сечений эффекта Комптона при энергиях 1173 и 1333 кэВ, штрих-пунктирная — среднее значение сечений фотоэффекта при энергиях 1173 и 1333 кэВ.

На рис. 4 видно, что при используемых энергиях в области малых Z доминирует эффект Комптона. Сечение Комптоновского рассеяния линейно зависит от Z. Следствием этих двух причин является близость зависимости R(Z) к константе в области малых Z. С увеличением Z всё более существенный вклад вносит фотоэффект и при энергии 7-квантов 662 кэВ начинает доминировать у веществ с Z > 87. Влияние фотоэффекта проявляется в выраженном монотонном уменьшении значения отношения R(Z).

Для оценки возможности применения радиоактивных источников для измерения атомного номера неизвестных объектов был проведён эксперимент [33] по просвечиванию тестовых образцов пучками 7-квантов (п. 2.1.3). В ходе эксперимента просвечивались образцы разной толщины из разных веществ. Для измерений в качестве детектора 7-квантов использовался 7-спектрометр из сверхчистого германия Canberra GC3019 [34].

В п. 2.1.4 приводятся результаты эксперимента и выводы. Был получен набор спектров 7-квантов, зарегистрированных детектором. Анализ спектров, определение площадей пиков полного поглощения Р и подложки В проводились на основе методики и с помощью программного обеспечения, разработанных в НИИЯФ МГУ [35,36].

В результате сопоставления измеренных Я с известной зависимостью, были оценены наиболее вероятные значения атомных номеров Z. Затем, по измеренным прозрачностям К и атомным номерам Z объектов согласно (9), были оценены количества вещества nD. Соответствующие результаты для образцов из кадмия и свинца приведены на рис. 5 и С. Хорошее согласие измеренных и ожидаемых значений атомного номера Z и количества

о, барн/атом 35

30

25

20

15

10

5

0

л0теа5игес1' см"

х1023

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

_| х1023

2.0 ПОгеа1, см-2

Рис. 6. Ожидаемые и измеренные значения количества вещества пй образцов из свинца V и кадмия о. Штриховая линия — соответствует случаю пОтаг = ?гД.еа|.

вещества п£) свидетельствует о возможности, при использовании изотопов Сэ-137 и Со-60, измерять атомный номер 2 и количество вещества пО объектов с 2 > 40.

В п. 2.2.1 описывается эксперимент [37] на ускорителе с перестраиваемой энергией. Линейный ускоритель электронов работал в режиме импульсного переключения энергии пучка. Ускоренный пучок электронов падал на вольфрамовую тормозную мишень толщиной 1.5 мм. В тормозной мишени образуется тормозное излучение с максимальной энергией, равной энергии ускоренных электронов. При помощи коллиматора, установленного сразу после тормозной мишени, вырезался тонкий пучок фотонов тормозного излучения. За коллиматором была расположена подвижная платформа с закреплёнными на ней тестовыми образцами из различных материалов. На расстоянии ~ 768 см от тормозной мишени располагался детектор излучения, прошедшего через образцы, с помощью которого измерялась суммарная доза С. Для уменьшения влияния рассеянного излучения, перед детектором располагался ещё один коллиматор, а сам детектор был окружен свинцовой защитой. Измерения проводились при четырёх разных значениях энергии ускорителя Т кз 7.3, 6.2, 5.5, 3.5 МэВ.

В п. 2.2.2 приводятся результаты эксперимента и выводы Для каждого объекта при четырёх фиксированных энергиях ускорителя Т измерялось значение дозы С излучения, прошедшего через исследуемый объект. Далее вычислялась прозрачность объекта К(Т, 2, п£>) = С(Т, 2, пО)/С(Т, 2,0). Затем набор прозрачностей, полученных при всех четырёх энергиях, используя МНК-подобный алгоритм сопоставлялся, с предварительно полученной таблицей прозрачностей, в результате чего получалось наиболее вероятное

Када«й р-Х-.

т =4..... л :

Рис. 7. Отклонения измеренных значений атомных номеров от ожидаемых: при различных значениях прозрачности К, рассчитанной для энергии фотонов 3.7 МэВ.

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

В третьей главе описывается выполненное автором моделирование отдельных элементов интроскопической системы, проведённое с использованием СЕАМТ4 [38], представляющего собой пакет программ для моделирования прохождения частиц через вещество на основе методов Монте-Карло. СЕА1МТ4 включает в себя инструменты для описания геометрии, используемых веществ и материалов, силовых полей, содержит множество физических моделей процессов взаимодействия частиц с веществом.

В разделе 3.1 рассматривается схема установки типовой системы предназначенной для просвечивания крупногабаритных транспортных средств. Поперечные размеры области, в которой осуществляется перемещение исследуемых объектов, определяются их максимальными шириной и высотой. Толщина и материал стенок определяют минимальную требуемую проникающую способность тормозного излучения. Максимальная высота может достигать 5.4 м, ширина — 3.6 м. У контейнеров-цистерн толщина стенки из нержавеющей стали может достигать 25 мм.

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

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

В п. 3.2.1 кратко рассматривается взаимодействие электронов с веществом. При взаимодействии с веществом электронов с энергиями (Те < 1 МэВ) основными являются ионизационные потери [30]. При больших энергиях увеличиваются радиационные потери. Электрон, испытывающий ускорение а в поле ядра, излучает энергию. Мощность излучения пропорциональна X"1. По этой причине, для изготовления тормозных мишеней-конверторов используют материалы с большим

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

В п. 3.2.3 обосновывается выбор оптимальной толщины тормозной мишени. При энергии электронов 3 МэВ интенсивность тормозного излучения в области 1.5-3 МэВ снижается при увеличении угла вылета фотонов. Толщина 2 мм и более является нецелесообразной, т.к. начинает уменьшаться интенсивность спектра в области верхней границы при углах вылета ~30°. При энергии электронов 10 МэВ интенсивность излучения существенно уменьшается при увеличении угла вылета фотонов. Особенно это заметно для 0.5 мм мишени, т.е. использование мишени толщиной 0.5 мм и менее нецелесообразно. В результате, был сделан вывод, что оптимальная толщина тормозной мишени из вольфрама составляет 1—1.5 мм.

В разделе 3.3 обосновывается выбор оптимального материала и размера элемента детектора. Для регистрации фотонов применяются различные типы детекторов. Для регистрации фотонов в диапазоне до 10 МэВ можно использовать полупроводниковые [39,40] и неорганические сцинтилляционные детекторы [41]. Однако, для регистрации тормозного излучения, получаемого с помощью импульсных электронных ускорителей применение полупроводниковых детекторов не оправдано. Оптимальным является использование неорганических сцинтилляторов соединённых с фотодиодами. В инспекционных системах обычно используются кристаллы СсГ\¥04 размером 3 х 3 х 20 мм3 [19], 5 х 5 х 30 мм3 [24].

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

и гамма-излучения. В данном случае основным ограничивающим фактором является поглощение фотонов люминесценции в длинном кристалле сцинтиллятора.

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

Показано что наибольшей вероятностью регистрации, благодаря наименьшей радиационной длине, обладает РЬ\¥04; этот же сцинтиллятор, благодаря наименьшему радиусу Мольера, имеет самые компактные поперечные размеры области энерговыделения. Обычно в интроскопических системах использовались кристаллы Сс1\\Ю4, обладающего большим световым выходом при близких прочих характеристиках. Однако, в связи с созданием детектров Большого Адронного Коллайдера, технология производства РЬ\У04 в больших количествах была существенно улучшена, что привело одновременно к улучшению характеристик и снижению стоимости. Если выбрать поперечный размер детектирующего элемента равным ~ 5 мм, полезный сигнал "расползается" только по двум соседним пикселям. Оптимальной глубиной кристалла РЬ'\\Ю4 можно считать величину 30-70 мм. При меньшей глубине - значительная часть фотонов будет пролетать сквозь кристалл. Большая глубина кристалла не целесообразна, т.к. часть кристалла сверх оптимальной толщины почти не увеличивает вероятность регистрации фотонов.

В разделе 3.4 описывается моделирование отклика детектора. Оно проводилось для детекторов из Р\\Ю размером 5 х 5 х 50 мм3. Кристаллы предполагается размещать с шагом 5 мм. Для оценки сигнала, получаемого на детекторе, было проведено моделирование взаимодействия фотонов с детектором. На рис. 8 показан отклик детектора (т.е. распреде-

1.0е+О0 1.0е-01 ™ 1.9е-02 1.0е-03

О 2 4 6 8 10

Выделенная энергии Еаер [НэВ]

Рис. 8. Отклик детектора при некоторых энергиях фотона. Детектор - Р\¥0 размером 5 х 5 х 50 мм3. Моделирование выполнялось для N = 5 • 10® фотонов с энергиями в диапазоне 0-10 МэВ с шагом 100 кэВ. Видны пики полного поглощения, одиночного и двойного вылета, а также т.н. комптоновская подложка.

ление вероятности выделенной энергии) для различных энергий фотонов. Знание отклика детектора необходимо для того, чтобы корректно моделировать спектр возможных значений энерговыделения при получении методом Монте-Карло оценок точности измерения 2.

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

В разделе 4.1 показывается возможность однозначного и устойчивого измерения £ с помощью монохроматических фотонов трёх энергий: = 0.5 МэВ, Е2 = 1.5 МэВ,£з = 3.7 МэВ, с использованием простого отношения

Д32 = 1п(АГ3) -ЩК,) = а(Е3,г) -сг{Еиг)

Д21 ь2^) <?2(Е2,г)

В разделе 4.2 рассматриваются возможности типовой инспекционной системы, использующей две энергии тормозного излучения на примере системы описанной в работе [24]. Граничные энергии тормозного излучения в указанной работе равны 9 и 5 МэВ.

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

Рис. 9. Поиск наилучшего соответствия измеренной пары значений прозрачностей при двух различных энергиях определённому атомному номеру. К\ и К2 прозрачности при двух различных энергиях. Кы и К2„ — измереные значения прозрачности объекта. 2 и Z'

— параметрические кривые, соответствующие веществам из различных диапазонов, г 1 —г4

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

ному 2 соседней ближайшей линией будет соответствующая не 2 ± 1, а соответствующая веществу из другого диапазона. В некоторых случаях величина отличия по 2 у соседних линий может достигать 98. Это приводит к тому, что ближайшими к измеренной точке будут линии соответствующие различным диапазонам атомных номеров 2 (рис. 9). Назовём область в которой линии соответствующие разным 2 пересекаются и накладываются областью неоднозначности.

Границы области неоднозначности в диапазоне прозрачностей 10~4 < К(Т = 9 МэВ) < < 1 показаны на рис. 10. Видно, что области неоднозначности соответствует существенная

Рис. 10. Области, в которых атомный номер определяется однозначно, расположены справа вверху и слева внизу (оставлены белыми). Область неоднозначности расположена между кривыми (выделена серым) и составляет 46% площади графика.

часть всех возможных комбинаций К(Т, 2, пИ). Таким образом можно утверждать, что в 54% всех возможных веществ и толщин атомный номер можно измерить однозначно, в пределах статистической погрешности. В 46% случаев измеренное значение 2 с вероятностью 50% будет близко к истинному и с вероятностью 50% будет соответствовать совершенно другому значению из другого диапазона веществ.

В разделе 4.3 показаны преимущества системы на основе трёх энергий тормозного излучения 9, 6 и 3 МэВ. Для трёх энергий, совокупность всех возможных комбинаций К(Т, 2, пО) в пределах 10~4 < К(Т = 9 МэВ) < 1 и 1 < Я ^ 99 представляет собой семейство параметрических кривых в трёхмерном пространстве прозрачностей. При использовании трёх энергий линии, соответствующие разным 2 не накладываются и не пересекаются. Две "конкурирующих" группы кривых, которые, при использовании двух энергий, накладываются одна на другую в двухмерном пространстве прозрачностей, при использовании трёх энергий отстоят друг от друга на некотором расстоянии в трёхмерном пространстве прозрачностей. Это позволяет, при низком уровне статистического разброса значений измеренных прозрачностей, избежать неоднозначности при измерении атомного номера.

В разделе 4.4 перечислены различные факторы влияющие на точность измерения 2\ Неоднозначность методов основанных на просвечивании объектов тормозным излучением

20

о —

1е-04

11

1е-02 К (Т = 9 МэВ)

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

В разделе 4.5 проведено модельное сравнение точности измерения 2 одинаковых объектов с помощю систем с двумя (Т = 9 и 5 МэВ) и тремя (Т = 9, 6 и 3 МэВ) граничными энергиями тормозного излучения. Рассмотрена возможность измерения атомного номера 2 независимо в каждом пикселе, т.е. без применения алгоритмов обработки изображения. Частота переключения энергий ускорителя - 200 и 300 Гц соответственно. Расстояние от тормозной мишени до линейки детекторов - Им. Детектор — линейка РЬ\¥С>4 размером 5 х 4 х 50 мм3 установленных с шагом 5 мм. Суммарное количество фотонов, попадающих на один детектирующий элемент (пиксель) в отсутствие поглотителя за одно измерение составляет ~ 5 ■ 106. Число испытаний 10000. В качестве критерия точности использовалось среднеквадратичное отклонение 2е„ измеренных значений 2^ от истинного атомного номера 2тел\ вещества, из которого изготовлен просвечиваемый объект. Чем меньше 2е1Г тем лучше. Результат для поглотителя из урана толщиной 1 см показан на рис. 11.

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

В разделе 4.7 проведён выбор оптимальных значений энергий для случаев двух и трёх граничных энергий. Рассмотрено влияние значений энергий на точность измерения 2. В качестве тестового объекта был выбран куб со стороной 10 см из урана. Ослабление тормозного спектра с верхней границей 9 МэВ при прохождении слоя урана толщиной 10 см составляет ~ 104. Измеренные в каждом пикселе значения прозрачностей усредняются по площади боковой грани куба, что уменьшает статистический разброс измеренных прозрачностей в 20 раз, т.к. площадь боковой грани куба в 400 раз больше площади одного пиксела. Среднее суммарное количество фотонов, попадающих на боковую грань, в отсутствие поглотителя за одно измерение составляет ~ 2 • 109.

В п. 4.7.1 проведена оптимизация для случая тормозного излучения с двумя граничными энергиями. Оптимальными значениями энергий являются: Тц = 10 МэВ и 71, = 5.5 МэВ. Однако, следует заметить, что выбор оптимальных значений граничных энергий тормозного излучения зависит от значения большей энергии, которое определяется требуемой проникающей способностью интроскопической установки и предельно допустимой мощностью дозы. Например, при значении большей энергии равном 6 МэВ

^real " 92 | Zerr - 32.1

Zreal * 92 | Zerr - 1.2

19 20 30 40 50 60 70 80 90 100 Измеренный атомный номер Znsr

10 20 30 40 50 60 70 80 90 100 Измеренный атомный номер г

Рис. 11. Сравнение точности измерения 2 в каждом пикселе для метода, использующего две энергии (Т = 9 и 5 МэВ) — слева, и для метода, использующего три энергии (Т = 9, 6 и 3 МэВ) — справа. Поглотитель — уран толщиной 1 см. Видно, что сгенерированные точки для метода, использующего две энергии (Г = 9 и 5 МэВ) попадают в область, где линии, соответствующие теоретическим значениям прозрачностей, накладываются и пересекаются, что приводит к неоднозначности. В ~ 50% случаев результатом измерения будет значение 2 ~ 43. Т.е. абсолютная ошибка по 2 может составлять ~ 49. Среднеквадратичные отклонения измеренных 2таг от истинного 21ел\ составляют 32.1 (Т = 9 и 5 МэВ) и 1.2 (Г = 9, 6 и 3 МэВ).

— оптимальным значением меньшей энергии будет 3.5 МэВ. Данными параметрами обладает импульсный линейный ускоритель описанный в [42] для инсиекционно-досмотрового комплекса описанного в [21].

В п. 4.7.2 проведена оптимизация для случая тормозного излучения с тремя граничными энергиями. Оптимальными значениями энергий являются: Тн = 10 МэВ, Тм = 5.5 МэВ и Т[, = 2 МэВ. На рис. 12 приведены распределения вероятностей измеренных значений 2 уранового куба со стороной 10 см при некоторых энергиях тормозных пучков.

10 20 30 40 50 60 70 80 90 100 Измеренный атомный номер гт5Г

Т = 9,5 МэВ [24]

2геа\ = 92 I 1етг =14-7

10 20 30 40 50 60 70 80 90 100 Измеренный атомный номер

0 10 20 30 40 50 60 70 80 90 100 Измеренный атомный номер гт5Г

Т = 9,6,3 МэВ [26]

10 20 30 40 50 60 70 80 90 100 Измеренный атомный номер 7т5г

Т = 10,5.5,2 МэВ

Т= 10,5.5 МэВ

Рис. 12. Сравнение распределений вероятностей измеренных значений ^ десятисантиметрового уранового куба при некоторых энергиях тормозных пучков. Среднее суммарное количество фотонов, попадающих на боковую грань, в отсутствие поглотителя за одно измерение составляет ~ 2 ■ 109. Метод двух энергий — слева. Метод трёх энергий — справа. Верхний ряд — до оптимизации. Нижний ряд — после оптимизации.

В разделе 4.8 приводится описание модельного сравнения оптимизированных систем с двумя и с тремя граничными энергиями в широком динамическом диапазоне прозрачно-стей различных веществ при разных средних мощностях электронного ускорителя. Оценка точности измерения 2 проводилась независимо для каждого пиксела без какой либо специальной обработки изображений прозрачностей. Сравнение было проведено при суммарном количестве 1.5 • 105, 1.5 ■ 107 и 1.5 • 10° тормозных фотонов, приходящихся на один пиксель детектора за одно измерение. Показано преимущество систем с тремя энергиями.

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

Список литературы

[1] Bjorkholm P. Detection of weapons of mass destruction // Port Technology International. - PT19-5/1.

[2] Medalia J. — Detection of Nuclear Weapons and Materials: Science, Technologies, Observations. — Congressional Research Service, 2010.

[3] U.S. Congress, Office of Technology Assessment. — Technology Against Terrorism: The Federal Effort, 1991.

[4] U.S. Department of Justice.— Survey of Commercially Available Explosives Detection Technologies and Equipment, 1998.

[5] Report on the workshop on the role of the nuclear physics research community in combating terrorism. — DOE/SC-OOG2, 2002.

[6] Snell M. P. Gamma ray technology: The practical container inspection alternative // Port Technology International. — 2002. — Mar. — no. 16. — P. 83-88.

[7] Rapiscan security products inc. — URL: http://www.rapiscan.com.

[8] Science applications international corporation. — URL: http://www.saic.com.

[9] Volumetric fraction measurement in oil-water-gas multiphase flow with dual energy gamma-ray system / Li Dong-hui, Wu Ying-xiang, Li Zhi-biao, Zhong Xing-fu // Journal of Zhejiang University SCIENCE. — 2005. — no. 6A(12). — P. 1405-1411.

[10] Study on Material Discrimination by Atomic Number Using Dual Energy Gamma-Rays / Y.M. Gil, Y.S. Lee, H.S. Lee et al.— 2005.— URL: http://psl.postech.ac.kr/ publication/dom_conference_env/kns_2005_f_ymgil.pdf.

[11] Zhang W., Li X., Xu Z. Feasibility study on application of ct using gamina-ray with doubleenergy to container security inspection // Procedia Engineering. — No. 7. — 2010. — P. 203-208.

[12] Bjorkholm P. Cargo screening: selection of modality // Port Technology International.— PT17-08/1.

[13] Bio-Imaging Research, Inc.— URL: http://www.bio-imaging.com.

[14] Bjorkholm P. Highly integrated x-ray system // Port Technology International. — PT20-35/1.

[15] L-3 communications security and detection systems. — URL: http://www.L-3com.com/ xray.

[16] Novikov V., Ogorodnikov S., Petrunin V. Dual energy method of material recognition in high energy introscopy systems // Problems of atomic science and technology. — 1999. — no. 4. - P. 93-95.

[IT] Ogorodnikov S., Petrunin V., Vorogushin M. Application of high-penetrating introscopy systems for recognition of materials // Proceedings of БРАС 2000, Vienna, Austria. — 2000. — P. 2583-2585. — URL: http: //accelconf. web. cern. ch/accelconf /е00/ PAPERS/WEP6B12.pdf.

[18] Vorogushin M., Ogorodnikov S., Petrunin V. Experiments on material recognition for 8 mev customs inspection system for trucks and large-scale containers // Proceedings of XX International Linac Conference, Monterey, California. — 2000. — P. 642-644. — URL: http://www.slac.Stanford.edu/econf/C000821/TUE13.pdf.

[19] Ogorodnikov S., Petrunin V. Processing of interlaced images in 4-10 mev dual energy customs system for material recognition // Physical review special topics - accelerators and beams. — 2002. — Vol. 5.

[20] Material discrimination technology for cargo inspection with pulse-to-pulse linear electron accelerator / S. Ogorodnikov, R. Apevalov, M. Arlychev et al. // Proceedings of 4th International Particle Accelerator Conference (IPAC 2013) : Shanghai, China.— 2013.— P. 3699-3701.— URL: http://accelconf.web.cern.ch/AccelConf/IPAC2013/papers/ thpwa033.pdf.

[21] ООО "Скантроник Системе". Реализованные проекты.— 2013.— URL: http:// scantronicsystems.com/proj ects.

[22] Bowser G., Husemann R. Technology advances help solve cargo inspection challenges // Port Technology International. — 1999. — no. 9. — P. 79-82.

[23] Bjorkholm P. Wmd detection // Cargo Security International. — 06/07 2005.

[24] Bjorkholm P. Dual energy detection of weapons of mass destruction // Port Technology International. — PT22-6/4.

[25] Zhiqiang C., Xuewu W. Cargo x-ray imaging technology for material discrimination // Port Technology International. — PT30-41/2.

[26] Многопучковый метод сканирования объектов / B.C. Ишханов, А.С. Курилик, Д.С. Руденко и др. // Известия РАН. Серия физическая.— 2008.— Т. 72, № 6.— С. 908-911.

[27] Комплексы для инспекции элементного состава грузов на основе линейных ускорителей электронов / А.А. Завадцев, Д.А. Завадцев, В.И. Каминский и др. — Москва : Национальный исследовательский ядерный университет МИФИ, 2011.— 148 с.

[28] Ishkhanov В., Shvedunov V., Pakhomov N., Varzar S. Multi-energy cargo inspection system based on an electron accelerator. — 2006.— US Patent App. 12/088,707; CN Patent 101,375,153; EP Patent App. EP20,060,851,714; WO Patent App. PCT/US2006/038,495.

[29] Yao W. et al. Review of particle physics // Journal of Physics G: Nuclear and Particle Physics. - 2006. - Vol. 33. — 1232 p.

[30] Частицы и атомные ядра. Практикум. Учебное пособие с грифом УМО по классическому университетскому образованию. Издание пятое, переработанное. Под общей редакцией B.C. Ишханова / О.И. Василенко, В.К. Гришин, Н.Г. Гончарова и др. — Издательство Московского университета, 2004. — С. 392.

[31] The livermore evaluated photon data library.— 1997.— URL: http://wMW-nds.iaea. org/epdl97/.

[32] Прайс В. Регистрация ядерного излучения: пер. с англ. — М. : Издательство иностранной литературы, 1960. — 464 с.

[33] Измерения атомного номера и количества вещества объектов с помощью радиоактивных источников / С.С. Белышев, B.C. Ишханов, A.C. Курилик, К.А. Стопани // Известия РАН. Серия физическая. — 2013. - Т. 77, № 7. — С. 1008-1011.

[34] Canberra industries inc. — URL: http://www.canberra.com.

[35] Определение выходов ядерных реакций на основе анализа цепочек распадов / С.С. Белышев, К.А. Стопани, С.Ю. Трощисв и др. // ВМУ. Серия 3. Физика. Астрономия. —

2011.- № 4.-С. 42-47.

[36] Studying photonuclear reactions using the activation technique / S.S. Belyshev, A.N. Er-makov, B.S. Ishkhanov et al. // Nuclear Instruments and Methods in Physics Research Section A. — 2014. — Vol. 745. — P. 133-137.

[37] Измерения атомного номера неизвестных объектов на ускорителе электронов с перестраиваемой энергией / А.Н. Ермаков, B.C. Ишханов, A.C. Курилик, В.И. Швсду-нов // Известия РАН. Серия физическая. - 2012. - Т. 76, № 11. - С. 1337-1339.

[38] Agostinelli S. et al. Geant4-a simulation toolkit // Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. - 2003. — Vol. 506, no. 3. — P. 250-303.

[39] Lutz G. Semiconductor Radiation Detectors. — Springer, 2007. — 353 p.

[40] Чилингаров А. Г. Координатные полупроводниковые детекторы в физике элементарных частиц // Физика элементарных частиц и атомного ядра. — 1992. — Т. 23.

[41] Inorganic Scintillators for Detector Systems / P. Lecoq, A. Annenkov, A. Gektin et al.— Springer, 2006. - 251 p.

[42] Ишханов Б., Шведунов В. Исследования и разработка ускорителей электронов в НИ-ИЯФ МГУ // Вестник Московского Университета. Серия 3. Физика. Астрономия. —

2012. - № 6. - С. 9-24.

Заказ № 101-Р/06/2014 Подписано в печать 26.06.14 Тираж 100 экз. Усл. пл. 1,2

ООО "Цифровичок", г. Москва, Большой Чудов пер., д.5

тел. (495)649-83-30 \ jj www.cfr.ru ; e-mail: zakpark@cfr.ru

 
Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Курилик, Александр Сергеевич, Москва

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ "МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ имени М.В. ЛОМОНОСОВА" ФИЗИЧЕСКИЙ ФАКУЛЬТЕТ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ ЯДЕРНОЙ ФИЗИКИ

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

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

04201460802 УДК 539.12.04

Курилик Александр Сергеевич

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

Специальность 01.04.16 физика атомного ядра и элементарных частиц

ДИССЕРТАЦИЯ на соискание учёной степени кандидата физико-математических наук

Научный руководитель д. ф.-м. н., профессор Б.С. Ишханов

Москва — 2014

Оглавление

Введение 4

Актуальность темы диссертации............................. 4

Задачи интроскопии................................. 4

Существующие решения............................................................5

Цели диссертационной работы............................................................11

Основные результаты, полученные в диссертации......................................11

Научная новизна............................................................................12

Ценность научной работы ................................................................13

Практическая значимость ................................................................13

Аппробация работы........................................................................13

Публикации ................................................................................16

Личный вклад автора......................................................................16

Структура и объём диссертации..........................................................16

Содержание диссертации..................................................................17

1. Взаимодействие фотонов с веществом 19

1.1. Сечения процессов взаимодействия фотонов с веществом......................19

1.2. Ослабление потока фотонов........................................................23

1.2.1. Монохроматические пучки фотонов........................................23

1.2.2. Пучки фотонов с непрерывным спектром ................................25

1.3. Принцип измерения Z и пБ..........................................................26

1.3.1. Монохроматические пучки фотонов........................................26

1.3.2. Тормозное излучение........................................................28

1.4. Анализ сечений ......................................................................29

1.4.1. Нормированные сечения....................................................29

1.4.2. Вклады сечений основных процессов в полное сечение..................35

1.4.3. Полные сечения..............................................................40

1.4.4. Выводы........................................................................40

1.4.5. Рекомендации................................................................41

1.4.6. Соотнешение вероятностей процессов рассеяния и поглощения фотонов 42

2. Экспериментальные измерения атомного номера объектов 44

2.1. Радиоактивные источники Сз-137 и Со-60........................................44

2.1.1. Выбор источников гамма-квантов..........................................44

2.1.2. Соотношение Я в случае Сз-137 и Со-60..................................45

2.1.3. Описание эксперимента......................................................46

2.1.4. Результаты эксперимента и выводы........................................48

2.2. Тормозное излучение................................................................53

2.2.1. Описание эксперимента......................................................53

2.2.2. Результаты эксперимента и выводы........................................55

3. Моделирование на СЕА1ЧТ4 58

3.1. Схема установки для интроскопии крупногабаритных объектов................58

3.2. Тормозное излучение................................................................59

3.2.1. Взаимодействие электронов с веществом..................................59

3.2.2. Моделирование тормозного излучения....................................60

3.2.3. Выбор оптимальной толщины тормозной мишени........................62

3.3. Выбор типа и размера детектора..................................................62

3.4. Отклик детектора....................................................................66

4. Сравнение методов и оптимизация их параметров 69

4.1. Пучки монохроматических фотонов трёх энергий................................69

4.2. Тормозное излучение с двумя граничными энергиями..........................71

4.3. Тормозное излучение с тремя граничными энергиями ..........................77

4.4. Статистические флуктуации........................................................80

4.5. Измерение Z в каждом пикселе....................................................81

4.6. Уменьшение уровня статистических флуктуаций................................89

4.7. Выбор оптимальных энергий........................................................89

4.7.1. Тормозное излучение с двумя граничными энергиями..................90

4.7.2. Тормозное излучение с тремя граничными энергиями..................91

4.8. Сравнение в широком динамическом диапазоне..................................96

Заключение 101

Литература 103

Введение

Актуальность темы диссертации

Задачи интроскопии

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

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

В разных странах, в силу географических особенностей, также существуют различные требования, предъявляемые к таможенным системам. В США основной целью применения досмотровых комплексов является обнаружение оружия массового поражения и его составляющих, радиоактивных и взрывчатых веществ, оружия и наркотиков. Считается, что наибольшую опасность представляет ядерное оружие массового поражения (WMD — nuclear weapons of mass destruction) на основе высоко обогащённого урана (HEU) и оружей-

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

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

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

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

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

Существующие решения

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

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

• рентгеновские установки с верхней границей энергии фотонов до сотен кэВ

• радиоактивные источники 7-квантов с энергиями до нескольких МэВ

• источники тормозного излучения на основе ускорителей электронов

• импульсные и непрерывные источники тепловых и быстрых нейтронов

• неионизирующее электромагнитное излучение, в том числе т.н. терагерцовое излучение (THz source)

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

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

• взаимодействие фотонов с атомными электронами: фотоэффект, когерентное и некогерентное рассеяние, рождение е~/е+-пар

• взаимодействие фотонов с атомными ядрами: рождение е~/е+-пар, резонансное рассеяние (флюоресценция), ГДР, (7,n), (7,f)

• взаимодействие нейтронов с атомными ядрами: (11,7), (n,n), (n,f)

• взаимодействие электромагнитного излучения с молекулами и атомами

• испускание 7-квантов атомными ядрами при переходе из возбуждённых состояний

• испускание нейтронов деления

По типу и энергии регистрируемых частиц:

• фотоны

• нейтроны

• заряженные частицы

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

• формирование 2D (2-dimensional) / 3D (З-dimensional) изображения распределения линейного коэффициента ослабления исследуемого объекта

• формирование 2D / 3D изображения распределения эффективного атомного номера и количества вещества исследуемого объекта

• формирование 2D / 3D изображения распределения изотопного состава исследуемого объекта

• получение информации о наличии в объекте того или иного вещества

• получение информации о наличии в объекте делящихся и радиоактивных материалов

Радиоактивные и делящиеся материалы при своём распаде испускают излучение и частицы, которые могут быть зарегистрированы детекторами, расположенными вне досматриваемого объёма. Такие методы, основанные на регистрации продуктов самопроизвольного радиоактивного распада, называются пассивными. Эти системы позволяют об-

наруживать даже незначительные количества радиоактивных веществ [6]. Однако, существенное количество делящихся и радиоактивных веществ может остаться незамеченным, если они будут перевозиться в специальном защищённом контейнере. В то же время, достаточно высока вероятность ложного срабатывания, т.к. многие обычные грузы содержат радиоактивные изотопы.

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

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

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

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

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

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

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

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

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

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

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

изотопов ^Св и 27Со для измерения 2 объектов, расположенных в крупногабаритных контейнерах [26,27].

Типовые активности источников, используемых в гамма-сканерах, составляют ^Св ~ 7.4- Ю10 Бк, источника 2