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

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

<3?

РАДЧЕНКО СЕРГЕЙ ЕВГЕНЬЕВИЧ

Алгоритмы обработки измерительной информации для информационно-измерительных систем обнаружения азотосодержащих веществ на основе гамма-спектрометрии

Специальность: 05.11.16 — Информационно-измерительные и управляющие системы (промышленность)

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

2 4 ЯНВ 2013

НОВОСИБИРСК-2013

005048790

Работа выполнена в Федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Новосибирский государственный технический университет»

Научный руководитель: доктор технических наук, профессор

Вострецов Алексей Геннадьевич

Официальные оппоненты: Кирьянов Валерий Павлович, доктор технических

наук, Федеральное государственное бюджетное учреждение науки «Институт автоматики и электрометрии Сибиркого отделения РАН», в.н.с. Полосаткин Сергей Викторович, кандидат физико-математических наук, Федеральное государственное бюджетное учреждение науки «Институт Ядерной Физики Сибирского отделения РАН им. Г.И. Будкера», с.н.с.

Ведущая организация: Федерального государственное бюджетное

образовательное учреждение высшего

профессионального образования «Национальный исследовательский Томский политехнический университет»

Защита состоится «12» февраля 2013 г. в 10-00 часов на заседании диссертационного совета Д 212.173.05 при Новосибирском государственном техническом университете по адресу: 630092, Новосибирск, пр. Карла Маркса, 20.

С диссертацией можно ознакомиться в библиотеке Новосибирского государственного технического университета.

Автореферат разослан «/<3» января 2013 г.

Ученый секретарь диссертационного совета

Шпилевая Ольга Яковлевна

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

Большинство известных взрывчатых веществ характеризуются повышенным содержанием азота по сравнению с обычными веществами. Известно, что при облучении исследуемого объекта гамма-излучением с энергией квантов 9.17 МэВ, данное излучение резонансно поглощается ядрами атомов азота и рассеивается другими веществами. Сравнение спектров первичного и вторичного излучений позволяет определить наличие взрывчатых веществ в составе объекта. К достоинствам метода следует отнести высокую чувствительность, а также низкий уровень радиационного облучения окружающей среды и сканируемых объектов: на порядок ниже облучения, происходящего при использовании рентгеновского излучения, и на два порядка ниже, чем при облучении быстрыми нейтронами. Пространственное разрешение, получаемое при использовании метода, составляет 0.1-1.0 см, что делает его перспективным для обследования крупных грузов.

Известно несколько зарубежных работ и патентов, описывающих результаты работы указанного метода. Начало исследований, а также основные результаты были получены М. Голдбергом и Д. Вартски. В России аналогичные работы ведутся в Институте Ядерной Физики СО РАН им. Г.И. Будкера группой ученых под руководством профессора A.B. Бурдакова, в Национальном исследовательском томском политехническом университете основателем данного направления был профессор B.JI. Чахлов.

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

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

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

Для достижения поставленной цели были решены следующие основные задачи:

1. Разработка математической модели спектра вторичного гамма-излучения.

2. Представление спектра в ортонормированном базисе и определение интенсивности отдельных спектральных линий.

3. Выделение признаков различия гамма-спектров при наличии и отсуствии азота в исследуемом объекте и формулировка задачи обнаружения азотосодержащего вещества в анализируемой пробе.

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

5. Разработка алгоритмов совокупной оценки параметров спектральных линий гамма-излучения.

6. Исследование эффективности разработанных алгоритмов.

7. Оценка практической реализуемости разработанных алгоритмов в информационно-измерительных системах.

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

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

Основные результаты, выносимые на защиту:

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

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

3. Алгоритм совокупной оценки параметров спектральных линий на основе полных достаточных статистик.

Научная новизна работы состоит в следующем:

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

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

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

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

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

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

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

Личный вклад. Автором работы осуществлена формулировка задач, а также получены основные результаты, выносимые на защиту.

Внедрение результатов исследований. Результаты диссертационной работы были внедрены в НИР, выполненные по проектам №2.1.1./3465

"Исследование генерации и поглощения резонансного гамма-излучения в реакциях на ядрах углерода и азота" и №2.1.1./11431 "Исследование генерации и поглощения резонансного гамма-излучения в реакциях на ядрах углерода и азота" в рамках Аналитической ведомственной программы «Развитие научного потенциала высшей школы (2009-2011 г.)». Кроме того, полученные в диссертации результаты внедрены в Институте Ядерной Физики СО РАН им. Г.И. Будкера (г. Новосибирск) при создании экспериментальной установки для обнаружения азотосодержащих веществ и исследования резонансного поглощения гамма-квантов ядрами атомов азота.

Апробация работы. Основные положения диссертационной работы были представлены в виде докладов и обсуждались на следующих мероприятиях: Международный форум по стратегическим технологиям IFOST-2009 (Вьетнам, 2009), X международная конференция «Актуальные проблемы электронного приборостроения АПЭП-2010» (Новосибирск, 2010), «Научная сессия НГТУ» (Новосибирск, 2011), Всероссийская научно-техническая конференция «Современные проблемы радиоэлектроники» (Красноярск, 2011), XI международная конференция «Актуальные проблемы электронного приборостроения АПЭП-2012» (Новосибирск, 2012).

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

Структура диссертации. Диссертационная работа изложена на 156 страницах, состоит из введения, пяти глав, заключения, списка литературы из 60 наименований и приложения. Основной текст изложен на 128 страницах, содержит 40 рисунков.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

В ИЯФ СО РАН им. Г.И. Будкера разработана и изготовлена установка, позволяющая проводить эксперименты по исследованию объектов на предмет наличия в них азотосодержащих веществ (рис. 1). Для генерации резонансных гамма-квантов используется мощный протонный ускоритель на основе тандема VITA. Он способен генерировать протонный пучок с энергией частиц около 2МэВ и током до 5мА в непрерывном режиме. Резонансные кванты получаются в реакции C13(p,y)N14. Реакция происходит в графитовой мишени, обогащенной изотопом С13, под действием протонов, ускоренных до энергии 1.75МэВ. При поглощении протона ядром С13 образуется возбужденное ядро N14, которое излучает гамма-квант с энергией 9.17МэВ. Эти гамма-кванты излучаются в форме конуса с углом раскрытия -80,7°.

Конический т

веерообразный /

пучок у

Ускоритель Мишень, ,3С , /

Материал, содержащий -

Детекторы детектор №4(у,р)С>3

гамма-излучения

Рис.1. Схема экспериментальной установки

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

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

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

Для описания формы функции интенсивности у'-го одиночного пика предложено использовать известное представление в виде гауссовой функции с

1 Г (е-£\)2]

единичной площадью /,(е -£,-)= ,— ехр<!--%—}, где е - энергия, £, -

у 1 л/2 ла [ 2сг \ 1

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

функция интенсивности одиночного у'-го пика описывается зависимостью

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

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

Рис.2. Модельная функция интенсивностей спектральных линий

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

энергия квантов. МэВ Рис. 3. Модель спектра гамма-излучения

Рассмотренная модель требует априорной информации о расположении пиков на энергетической оси, их ширине и интенсивности. Если при решении задачи обнаружения азотосодержащего вещества вся необходимая для построения соответствующего алгоритма информация является известной (известно положение обнаруживаемых пиков, ширина пика задается аппаратной функцией детектора, а знание интенсивности пика, как будет показано ниже, не требуется), то при решении задачи анализа спектра положение пиков и их интенсивности являются неизвестными. Поэтому в общем случае для анализа спектра необходима универсальная модель, позволяющая адекватно описывать наблюдаемый спектр. Для построения такой модели предложено представить вектор отсчетов функции интенсивностей спектральных линий в ортонормированном базисе, получаемом следующим образом. На первом этапе в качестве базисных принимается набор векторов _ (j = 0...Nj--1, Nj- -число базисных векторов, N -число

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

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

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

вектора отсчетов наблюдаемого спектра на пространство, натянутое на базисных ортонормированных векторах:

Л7-1

]=о

где Е - аппроксимирующий вектор отсчетов функции интенсивностей

спектральных линий, aJ■ = , Т7 - N -мерный вектор отсчетов

наблюдаемого спектра, (•,•) -скалярное произведение, у = -1, — N-

мерные базисные ортонормированные вектора, N^ - число базисных векторов.

Реальный спектр и результат аппроксимации функции интенсивностей спектральных линий по его N = 2000 отсчетам при использовании N^ - 76

базисных векторов приведены на рис. 4 (функция интенсивностей показана полужирной сплошной линией).

Рис. 4. Результат аппроксимации реального спектра

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

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

Результат применения этой процедуры к участку спектра 7.1 - 9.6 МэВ (рис. 5) приведен на рис. 6.

энергия квантов. МэВ

Рис.

Рис. 6.

Несмотря на то, что на рис. 5 визуально можно выделить не более четырех пиков, предложенная процедура позволила установить наличие 5 основных пиков (линии 9.17 МэВ (штриховая линия), 8.62 МэВ (точечная линия), 8.80 МэВ (сплошная линия), 8.98 МэВ (штрихпунктирная линия), 8.06 МэВ (точечная линия) и трех пиков single-escape.

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

Анализируемый объект облучается потоком гамма-квантов с широким энергетическим спектром. Если объект не содержит азота, то спектр излучения после прохождения объекта практически не изменится. Наличие атомов азота

5. Аппроксимация участка спектра 7.1 - 9.6 МэВ

ю

энергия квантов, МэВ Выделение вклада отдельных спектральных линий

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

Первый подход основан на использовании свойств спектральных линий 2.36 МэВ и 9.17 МэВ. Первый пик практически не поглощается атомами азота, в то время как второй подвергается резонансному поглощению. Степень поглощения зависит от концентрации азота. В отсутствие поглощения отношение интенсивностей этих пиков определяется генератором гамма квантов, является известным и равным К0. При наличии поглощения гамма-квантов отношение интенсивностей пиков 2.36 МэВ и 9.17 МэВ будет зависеть от концентрации атомов азота и превысит величину К0.

В качестве исходных данных для синтеза алгоритма обнаружения приняты два выборочных вектора - вектор х = {х,,г =1,...,и}, составленный из отсчетов спектра на выходе детектора в области 2.36 МэВ, и вектор У = »>'/>' = 1,•••,«}, составленный из отсчетов спектра в области энергетического пика 9.17 МэВ. Для описания функции интенсивности спектральных линий использовалась разработанная во второй главе модель, в соответствии с которой для энергетического пика 2.36 МэВ использована функция КБ/(Де,), а для энергетического пика 9.17 МэВ - функция 5У(Де,), где /(Де,), / = 1 отсчеты нормированной гауссовской функции с единичной площадью, 5" -неизвестная интенсивность пика 9.17 МэВ, К - коэффициент отношения интенсивностей пиков 2.36 МэВ и 9.17 МэВ, Де,- = е; -£■ - расстояние между /-м отсчетом спектра в области соответствующего пика и положением этого пика (2.36 МэВ или 9.17 МэВ) на оси энергий. Число отсчетов и выбрано таким образом, чтобы площадь под аппроксимирующей функцией /(Де,) между крайними отсчетами в пределах пика была не менее 99%.

Совместное распределение вероятностей векторов х и у имеет следующий вид:

У) =-~„-Ы-ехр{|>ё(га/(Де,•))]*,. + ¿[1ое(5/(Ае, ))]>>,} •

м

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

Н0 : К = К0 (азот отсутствует);

III: К> К0 (азот присутствует). ' '

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

1 V

1

1

2 я

хехр<--

(2)

./=1 ко ¡=Г

где X,- = 2^ + 1, % = 1 = 1,...,и.

Распределение (2) инвариантно относительно группы аддитивных преобразований С = {£: X, —» Х1 + Ъ1\У1 —► + , отражающих априорную

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

<*>(*. У) =

+1 - + 1 УЛ^) > С / 2;

/=1

О, в противном случае. Пороговая константа С определяется исходя из допустимого значения вероятности а принятия ложного решения о поглощении у-квантов и задается выражением:

где - квантиль порядка (1 - а) стандартного нормального распределения.

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

Совместное распределение выборочных векторов в данном случае имеет

вид

exp -SoX/CAe,.)] p(x, у) =-Ц-^--exp-j I>, log(S0/(Ae,))

ГМ ы

ехр^^ДДе,.)]

-exp^log^/iAe,))^

Псио Ul

;=i

и обладает достаточными статистиками Т(х,у) =

\ )=1 /=1 / 1=1

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

у(и,Т)=-

exp

1=1

и\{т-и)\

exp-! [/log

+ 7* log

<=i

•• (3)

Данное распределение характеризуется полезным Л = log мешающим ju = log

параметрами.

^ХДАе,-) . /=1

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

Н0:Л = 0;//е (-оо,оо) (азот отсутствует);

(4)

Я,: X > 0; /л е (-со, от) (азот присутствует).

В работе получен равномерно наиболее мощный несмещенный алгоритм решения задачи (4), решающая функция которого имеет следующий вид:

ГО, и < С(Т), (р{и,Т) = \у{Т), и = С(Г), [1, и>С{Т).

Здесь С(Г) - пороговая функция, у(Т) - параметр рандомизации, зависящие от достаточной для мешающего параметра ¡л статистики Т и заданного уровня а вероятности ложного принятия решения о наличии азотосодержащего вещества. Для их конкретизации при каждом фактическом значении статистики Т необходимо решить уравнение

J'<p(U,T)w(U | Я = 0,T)dU = а.

Исследование эффективности разработанных алгоритмов проводилось методом имитационного моделирования на ЭВМ. Зависимости вероятности правильного обнаружения от соотношения интенсивностей пиков для рассмотренных алгоритмов приведены ниже на рис. 7. Для каждого значения ^о / 5, моделировалось 5000 пар спектров, для каждой пары вычислялась решающая функция, вероятность правильного обнаружения алгоритма определялась по количеству положительных решений для выбранного значения 50 / . Уровень вероятности ложной тревоги задавался на уровне а = 0.01.

I 0.2

л _ -J7*____и_ ,

о -^-^---1

g I 1.05 1.1 1.15 1.2 1.25

с.

% соотношение интенсивностеи пиков Рис. 7. Зависимость вероятности правильного обнаружения от степени поглощения а) для алгоритма, основанного на сравнении интенсивностей двух пиков одного спектра; б) для алгоритма, основанного на методе контраста

Как видно из рисунков, оба алгоритма достаточно эффективны: при 20-процентном превышении пика в опорной выборке над пиком в рабочей выборке (степень поглощения около 20%) вероятность правильного обнаружения оказывается более 0.9, при этом контрастный алгоритм (линия «б») оказывается несколько эффективнее.

Эффективность работы алгоритмов также зависит и от набранной статистики при проведении измерений, объем которой определяется длительностью измерения. На рис. 8 приведены зависимости вероятности принятия правильного решения для контрастного алгоритма при разных интенсивностях опорного пика.

Рис. 8. Зависимость вероятности правильного обнаружения от степени поглощения при разных интенсивностях опорного пика

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

Форму отдельной спектральной линии можно аппроксимировать гауссовой кривой со средним значением, соответствующим положению спектральной линии на энергетической шкале, параметром <т2, определяющим ширину спектральной линии, а также параметром 5, определяющим интенсивность линии.

В окрестности одной спектральной линии, определяемой шириной

аппаратной функции детектора, формируется выборка х = (хП1,..,дсИ2), состоящая из п отсчетов, взятых в точках на оси энергий гамма-квантов, / = «|,..и2 (1 + я2-И| =п). Отсчеты х1 отстоят друг от друга по оси энергий на значение г, определяемое величиной младшего разряда АЦП, применяемого в системе: г,- =г'-г. Смещение линии относительно начала координат выражается параметром /0.

Интенсивность одного отсчета определяется выражением:

А- =-

('/-'о)

" 2а2

у/Ъгсг

где параметр 5 определяет интенсивность спектральной линии.

Совместная плотность распределения вероятностей определяется выражением "2

выборки

1

ехр

2а2

т,

Сг'оГ

2а2

жа

■ (5)

В работе показано, что распределение (5) обладает полными

»2 "2

достаточными статистиками 7|(х) = ^х, и T2{x)='^{t¡-t0)2x¡, на основе

которых с использованием теоремы Лемана-Шеффе получены эффективная оценка интенсивности и оценка параметра ширины спектрального пика:

¿('/"'о

I=71,

"2 ¡=Щ

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

Анализ погрешности полученных оценок проводился методом численного моделирования на ЭВМ. Для каждого значения оцениваемого параметра генерировались 1000 независимых выборок, по которым оценивались интенсивность исследуемой спектральной линии (параметр 5) и ее ширина (параметр а2) и вычислялись систематическая и случайная составляющие погрешности.

Результаты моделирования приведены на рис. 9.

ыо1

истинное значение интенсивности

а)

е

о а:

3 |

1*Ш3 1.Ы03 1.2x10* 1.3х105

истинное значение интенсивности б)

Рис. 9. Зависимость погрешности оценки интенсивности спектральной линии от истинного значения интенсивности а) случайная составляющая; б) систематическая составляющая

Аналогичным образом было проведено моделирование работы алгоритма оценки ширины спектральной линии (рис. 10).

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

2 0054 §

& „

3 0.048 8

истинное значение параметра ширины

а)

истинное значение параметра ширины б)

Рис. 10. Зависимость погрешности оценки параметра ширины спектральной линии от истинного значения параметра а) случайная составляющая; б) систематическая составляющая

Оценка параметра ширины спектральной линии признана несостоятельной, поскольку с ростом объема выборки погрешность оценивания

не стремится к нулю. Однако в тех случаях, когда требуемая погрешность оценки находится на уровне >=5%, ее можно использовать для оценки ширины спектральной линии.

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

1. В состав досмотрового комплекса входит ЭВМ, осуществляющая управление комплексом. Модуль принятия решений (МПР) о наличии взрывоопасных грузов реализован в виде пользовательского приложения.

2. Независимо от конфигурации и организации системы управления комплексом, модуль принятия решений реализован в виде отдельного устройства.

Исходя из особенностей задачи обнаружения повышенных концентраций азота, сформулированы основные требования к программно-аппаратной реализации работы алгоритмов:

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

2. Возможность применения разных алгоритмов обнаружения.

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

4. Максимально возможный уровень защищенности от внешнего вмешательства в функционирование модуля.

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

6. Программная надежность, обеспечивающая максимально возможное время наработки на отказ.

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

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

В итоге структурная схема обработки данных выглядит следующим образом (рис. 12).

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

Рис. 12. Структурная схема обработки данных

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

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

Таким образом, структурная схема обработки данных принимает вид:

¡СУ;

Рис. 13. Структурная схема обработки данных при аппаратной реализации МПР

С целью определения элементной базы сформулированы следующие основные требования к МПР:

1. Возможность работы от стандартной сети.

2. Возможность автономной работы от аккумулятора.

3. Ведение журнала событий с привязкой ко времени.

4. Возможность переноса журнала на съемный носитель.

5. Расширенная индикация режимов работы и результатов обнаружения.

6. Возможность работы в реальном времени.

Согласно приведенным требованиям произведен выбор элементов основных составляющих устройства, в частности, микроконтроллера. Было предложено использование архитектуры ARM Cortex МЗ как одной из наиболее современных архитектур, оптимальных по соотношению цена — производительность. Показано, что применение микроконтроллера STM32F103RCT6 удовлетворяет всем требованиям к функциональности МПР.

В заключении сформулированы основные результаты работы, которые сводятся к следующему:

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

2. Разработаны и исследованы алгоритмы обнаружения азотосодержащих веществ на основе поглощения энергетического пика 9.17МэВ во вторичном гамма-излучении объекта. Алгоритмы имеют достаточно высокую эффективность, позволяя регистрировать малозаметные изменения интенсивности отдельных спектральных пиков.

3. На основе применения следствий теоремы Лемана-Шеффе для оценок на основе полных достаточных статистик получены алгоритмы оценивания параметров интенсивности и ширины спектральных линий вторичного энергетического спектра гамма-излучения после прохождения через исследуемый объект.

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

Приложение содержит исходный код программы по обработке гамма-спектров на ЭВМ, а также акты внедрения.

СПИСОК ПУБЛИКАЦИЙ ПО ТЕМЕ ДИССЕРТАЦИИ

Статьи в изданиях из Перечня ведущих российских рецензируемых научных журналов и изданий

1. С. Е. Радченко. Метод обнаружения поглощения гамма-квантов при прохождении через азотосодержащее вещество / А. Г. Вострецов, А. В. Бурдаков, С. Е. Радченко, А. С. Кузнецов, Ю. С. Суляев // Автометрия, 2010. Т. 46. №3. С. 22-29.

2. С. Е. Радченко. Контрастное обнаружение спектральных пиков при прохождении через вещество / С. Е. Радченко // Доклады Академии наук высшей школы Российской Федерации, 2011. №1. С. 127-136.

3. С. Е. Радченко. Аппроксимация энергетического спектра гамма-излучения в ортонормированном базисе / А. Г. Вострецов, С. Е. Радченко // Научный вестник Новосибирского государственного технического университета, 2011. № 4 (45). С. 3-11.

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

4. S.E. Radchenko. Statistical Algorithm for Detection of Gamma Ray Spectral Peak Absorption / A.G. Vostretsov, S.E. Radchenko // IFOST 2009. Proceedings of the 2009 International Forum on Strategic Technologies. October 2123, 2009, Ho Chi Minh City, Vietnam. Session 2. Information and Communication Technology- P.82-86 [Статистический алгоритм обнаружения поглощения спектральных пиков гамма-излучения].

5. С. Е. Радченко. Методы и алгоритмы обнаружения азотосодержащих веществ / А. Г. Вострецов, С. Е. Радченко // Материалы VIII международной конференции "Актуальные проблемы электронного приборостроения АПЭП - 2010" в 7 томах. Том 4. Новосибирск, 22-24 сентября, 2010. Новосибирск: Изд-во НГТУ, 2010. С. 162-166

6. С. Е. Радченко. Методика интерпретации энергетического спектра вторичного гамма-излучения на основе векторной аппроксимации / С.Е. Радченко // Материалы XI международной конференции "Актуальные проблемы электронного приборостроения АПЭП-2012" в 7 томах. Том 2. Новосибирск, 02-04 октября, 2012. Новосибирск: Изд-во НГТУ, 2012. С.37-41

7. S.E. Radchenko. Method of Gamma-Spectrum Interpretation Using Vectorial Approximation / A.G. Vostretsov, S.E. Radchenko // APEIE 2012. Proceedings of the 2012 11th International Conference on Actual Problems of Electronic Instrument Engineering, Volume 1, October 2^1, 2012, Novosibirsk, Russia. P. 11-14 [Методика интерпретации спектра гамма-излучения на основе векторной аппроксимации].

Публикации в других изданиях

8. С. Е. Радченко. Совместная оценка параметров спектральных линий энергетического спектра гамма-излучения / С. Е. Радченко // Современные проблемы радиоэлектроники: сб. науч. тр. - Красноярск: Сиб. Федер. Ун-т, 2011. С. 174- 179.

Отпечатано в типографии Новосибирского государственного технического университета 630092, г.Новосибирск, пр. К. Маркса, 20, Тел./факс (383) 346-08-57 Формат 60 х 84/16. Объем 1.5 п.л. Тираж 100 экз. Заказ 85. Подписано в печать 25.12.2012 г.

0

Перечень используемых сокращений ^

Введение (Обнаружение азотосодержащих веществ гамма-спектрометрическим методом) ^

1. Аналитический обзор современного состояния проблемы контроля состава веществ применительно к задаче обнаружения азотосодержащих веществ ^ у

1.1 Средства обработки масс-спектров ^ у

1.2. Обнаружение взрывчатых веществ

1.2.1. Контактные методы обнаружения ВВ

1.2.2. Бесконтактные методы обнаружения ВВ

1.2.2.1. Методы на основе применения рентгеновского излучения

1.2.2.2. Электромагнитные методы

1.2.2.3. Методы на основе применения радаров

1.2.2.4. Методы на основе применения широкополосного терагерцового излучения

1.2.2.5. Методы на основе применения нейтронного излучения

1.2.2.6. Методы на основе применения гамма-излучения

1.3. Выводы

2. Модели гамма-излучения и формулировка задачи обнаружения азотосодержащих веществ ^

2.1. Выбор и обоснование модели спектра гамма-излучения 3 g

2.2. Представление спектра в ортонормированном полном базисе ^

2.2.1. Построение ортонормированного базиса ^

2.2.2. Представление спектра в ортонормированном базисе ^

2.2.3. Оптимизация параметров аппроксимации

2.2.4. Применение аппроксимации спектра в ортонормированном базисе

2.3. Признаки различия гамма-спектров различных веществ и формулировка задачи обнаружения повышенного содержания азота в анализируемой пробе

2.4. Выводы

3. Алгоритмы обнаружения азотосодержащих веществ в условиях априорной неопределенности спектра гамма-излучения

3.1. Алгоритм обнаружения азотосодержащих веществ на основе сравнения интенсивностей спектральных линий

3.2. Контрастный алгоритм обнаружения азотосодержащих веществ

3.3. Характеристики эффективности алгоритмов

3.4. Выводы

4. Оценка параметров спектральных линий'

4.1. Оценка интенсивности спектральных линий

4.2. Оценка ширины спектральных линий

4.3. Анализ погрешностей оценок

4.4. Выводы

5. Оценка практической реализуемости разработанных алгоритмов

5.1. Разработка программной реализации алгоритмов обнаружения повышенного содержания азота в объектах

5.1.1. Обнаружение повышенной концентрации азота по одному спектру ^

5.1.2. Контрастное обнаружение повышенного содержания азота

5.1.3. Аппроксимация спектра в ортонормированном полном базисе ^

5.2. Разработка модуля принятия решений о присутствии повышенных концентраций азота ^ ^

5.2.1. Оценка быстродействия устройства

5.2.2. Выбор и обоснование элементной базы щ

5.2.2.1. Выбор микроконтроллера

5.2.2.2. Питание устройства цу

5.2.2.3. Индикация режимов работы

5.2.2.4. Интерфейсы ввода-вывода | ^

5.3. Выводы 120 Заключение 122 Список литературы 124 Приложение . .

Перечень используемых сокращений

ВУ взрывное устройство

ВВ взрывчатое вещество

АЦП аналого-цифровой преобразователь

ЗУ запоминающее устройство

ФЭУ фотоэлектронный умножитель

ЭВМ электронная вычислительная машина

ИЦР ионный циклотронный резонанс

AMS Accelerator Mass Spectrometry

СВЧ сверхвысокочастотное излучение

ТНТ тринитротолуол

РНМ равномерно наиболее мощный

СКО среднеквадратическое отклонение

ОС операционная система

СГ система генерации

CP система регистрации

USB Universal Serial Bus

ПС плата спектрометра

МПР модуль принятия решений

CAN Controller Area Network

UART Universal Asynchronous Receiver-Transmitter

ARM Advanced RISC Machine

RISC Restricted Instruction Set Computer

SDIO Secure Digital Input/output

Telecommunications Industry Association / Electronic Industries

TIA/EIA

Alliance

Введение (Обнаружение азотосодержащих веществ гамма спектрометрическим методом)

Помимо использования нейтронного и электронного излучений, большие возможности для обнаружения ВВ представляет использование высокоэнергетического гамма-излучения. Как правило, в качестве первичного излучения выступает луч гамма-квантов с энергией около 9.2 МэВ [1-3].

Большинство известных взрывчатых веществ характеризуются повышенным содержанием азота по сравнению с обычными веществами. Известно, что при облучении исследуемого объекта гамма-излучением с энергией квантов 9.17МэВ, данное излучение резонансно поглощается ядрами атомов азота и рассеивается другими веществами [2]. Сравнение спектров резонансного и нерезонансного излучений позволит определить наличие взрывчатых веществ в составе объекта. К достоинствам данного метода необходимо отнести высокую чувствительность, а также низкий уровень радиационного облучения окружающей среды и сканируемых объектов - на один порядок ниже облучения, происходящего при использовании рентгеновского излучения, и на два порядка ниже, чем при облучении быстрыми нейтронами. Пространственное разрешение, получаемое при использовании метода, составляет 0.1-1.0 см, что делает его идеальным для обследования крупных грузов.

Известно несколько зарубежных работ и патентов, описывающих результаты работы указанного метода [4-6]. Первые результаты были получены М. Голдбергом и Д. Вартски. В России аналогичные работы ведутся в Институте Ядерной Физики СО РАН им. Г.И. Будкера группой ученых под руководством профессора A.B. Бурдакова [3], в Национальном исследовательском томском политехническом университете основателем данного направления был профессор B.JI. Чахлов. В ИЯФ СО РАН разработана и изготовлена установка, позволяющая проводить эксперименты по исследованию объектов на предмет наличия взрывчатых веществ (рис. 1). Для генерации резонансных квантов используется мощный протонный ускоритель на основе тандема VITA [3]. Он способен генерировать протонный пучок с энергией частиц около 2МэВ и током до 5мА в непрерывном режиме. Резонансные кванты получаются в реакции C13(p,y)N14. Реакция происходит в графитовой мишени, обогащенной изотопом С13, под действием протонов, ускоренных до энергии 1.75МэВ. При поглощении протона ядром С13 образуется возбужденное ядро N14, которое излучает гамма-квант с энергией 9.17МэВ. Эти гамма-кванты излучаются в форме конуса с углом раскрытия -80,7°.

Исследуемый объект и сцинтилляционный детектор гамма-квантов находятся за коллиматорами, которые задерживают все гамма-кванты, угол вылета которых отличается от резонансного. В спектре гамма-квантов, регистрируемых детектором, кроме пика с энергией 9,17 МэВ присутствуют пики с меньшей энергией, так как возбужденное ядро азота имеет несколько уровней, и снятие возбуждения может происходить каскадным образом. Гамма-излучение возникает также при возбуждении еще пяти резонансов с каскадным распадом. Кроме того, в бомбардируемой протонами мишени присутствует изотоп С , дающий в аналогичной реакции интенсивную линию 2,36 МэВ.

Ускоритель Мишень, 13(

Материал, содержащий азот

Пучок протонов 1,75 МэВ

Конический веерообразны! пучок

13(Р,у)№4

Детекторы гамма-излуч ения детектор

4(у,р)С13

Рис. 1. Схема экспериментальной установки

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

Для формирования спектров с достаточной для анализа статистикой за малые промежутки времени (десятки миллисекунд) используется модуль спектрометра, также разработанный в Институте Ядерной Физики СО РАН им. Г.И. Будкера. В основу его работы положена оригинальная методика обработки наложенных событий [7]. С помощью дискретного фильтра форма исходного сигнала детектора преобразуется в короткий сигнал гауссовой формы, амплитуда которого пропорциональна энергии гамма-кванта, а паразитный пьедестал близок к нулю. Для формирования по отклику сцинтиллятора импульсной характеристики данного фильтра используется следующая последовательность действий:

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

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

3. с помощью обратного преобразования Фурье формируется импульсная характеристика фильтра.

Величина ошибки при обработке таким фильтром двух наложенных событий не превышает 1.3% от амплитуды, а характерная длительность выходного сигнала фильтра намного меньше характерной длительности сигнала детектора, что позволяет сделать вывод об эффективности данной методики обработки наложенных событий. Она реализована в цифровом узле спектрометра, применяемого для регистрации спектров гамма-излучения [7].

Данный узел включает в себя:

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

2. формирующий фильтр, преобразующий сигнал детектора в короткий сигнал гауссовой формы

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

4. два модуля автоинкрементных ЗУ, работающих в режиме набора статистики поочередно

5. набор управляющих регистров

6. буфер обмена данными и контроллер канала связи

Узел отбора полезных событий сравнивает текущее значение амплитуды сигнала с заданным оператором опорным уровнем, фиксирует факт наличия события с энергией, превышающей этот порог, и на последующем временном интервале, несколько большем полуширины выходного сигнала фильтра, выделяет его максимальное значение. В процессе выполнения этой операции используется интерполирующий полином, позволяющий по последовательности текущих амплитудных значений фильтра вычислить истинное значение его максимума с погрешностью менее 0.1%. Далее полученный отсчет умножается на корректирующий коэффициент, величина которого зависит от текущего значения температуры сцинтиллятора, и передается блоку автоинкрементных ЗУ. В этом ЗУ формируется распределение, связывающее друг с другом энергию и количество событий на интервале набора статистики. С целью исключения мертвого времени и совмещения процедур накопления и считывания информации используется два модуля автоинкрементных ЗУ, работающих поочередно. Пока в одном из:них в течение отведенного ему кванта времени идет процесс накопления данных, второе «переписывает» накопленный информационный массив в промежуточный буфер обмена для его последующей передачи контроллером канала связи серверу диагностики. Контроллер канала связи, помимо пересылки данных серверу, обеспечивает доступ последнего к регистрам управления анализатора, к ЗУ опорных констант, хранящему массивы коэффициентов формирующего фильтра корректора интенсивности вспышек, к регистрам, содержащим результаты измерений температуры сцинтиллятора и напряжения высоковольтного источника питания ФЭУ. Поддержка контроллером канала связи этих процедур позволяет программными средствами адаптировать его шкалу и характеристики узла обработки данных к требованиям эксперимента и параметрам сцинтиллятора. Таким образом, предложенная методика регистрации гамма-спектров обеспечивает надежную и качественную регистрацию спектров гамма-излучения.

Однако, как правило, подобные усилия, схемотехнические и алгоритмические усовершенствования становятся завершающим этапом обработки информации. Изучение ряда работ по гамма- и масс-спектрометрии [1-6, 8-15] приводит к выводу о том, основное внимание авторов уделяется развитию аппаратной части информационно-измерительных систем, а также совершенствованию аппаратуры и методов регистрации происходящих событий. Известно несколько работ, призванных повысить качество получаемых данных [16, 17, 18]. Но при этом собственно интерпретация результатов зачастую сопровождается достаточно простыми операциями, которые сводятся, как правило, к подсчету и сравнению количества событий в определенных областях спектров. В то же время чувствительность метода и достоверность обнаружения веществ во многом определяется характеристиками алгоритмов обработки информации, поступающей от первичных измерительных преобразователей, которая, в силу специфики применяемых детекторов высокоэнергетического излучения, носит статистический характер. Существующие алгоритмы не учитывают возможные изменения фоновой и полезной составляющих спектра и не обеспечивают автоматическую подстройку под данные изменения с поддержкой постоянного уровня ложного обнаружения и максимизацией вероятности правильного обнаружения. Поэтому тема диссертационной работы, посвященной созданию эффективных методов и алгоритмов обработки первичной измерительной информации в системах обнаружения азотосодержащих веществ, является актуальной.

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

Для достижения поставленной цели были решены следующие задачи:

1. Разработка математической модели спектра вторичного гамма-излучения.

2. Представление спектра в ортонормированном базисе и определение интенсивности отдельных спектральных линий.

3. Выделение признаков различия гамма-спектров при наличии и отсуствии азота в исследуемом объекте и формулировка задачи обнаружения азотосодержащего вещества в анализируемой пробе.

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

5. Разработка алгоритмов совокупной оценки параметров спектральных линий гамма-излучения.

6. Исследование эффективности разработанных алгоритмов.

7. Оценка практической реализуемости разработанных алгоритмов в информационно-измерительных системах.

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

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

Научная новизна работы состоит в следующем:

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

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

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

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

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

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

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

Личный вклад. Автором работы осуществлена формулировка задач, а также получены основные результаты, выносимые на защиту.

Внедрение результатов исследований. Результаты диссертационной работы были внедрены в НИР, выполненные по проектам №2.1.1./3465 "Исследование генерации и поглощения резонансного гамма-излучения в реакциях на ядрах углерода и азота" и №2.1.1./11431 "Исследование генерации и поглощения резонансного гамма-излучения в реакциях на ядрах углерода и азота" в рамках Аналитической ведомственной программы «Развитие научного потенциала высшей школы (2009-2011 г.)». Также полученные в диссертации результаты внедрены в Институте Ядерной Физики СО РАН им. Г.И. Будкера (г. Новосибирск) при создании экспериментальной установки для обнаружения азотосодержащих веществ и исследования резонансного поглощения гамма-квантов ядрами атомов азота.

Апробация работы. Основные положения диссертационной работы были представлены в виде докладов и обсуждались на следующих мероприятиях: Международный форум по стратегическим технологиям ПЮ8Т-2009 (Вьетнам, 2009), X международная конференция «Актуальные проблемы электронного приборостроения АПЭП-2010» (Новосибирск, 2010), «Научная сессия НГТУ» (Новосибирск, 2011), Всероссийская научно-техническая конференция «Современные проблемы радиоэлектроники» (Красноярск, 2011), XI международная конференция «Актуальные проблемы электронного приборостроения АПЭП-2012» (Новосибирск, 2012).

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

Структура диссертации. Диссертационная работа изложена на 156 страницах, состоит из введения, пяти глав, заключения, списка литературы из 60 наименований и приложения. Основной текст изложен на 128 страницах, содержит 40 рисунков.

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

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

3. Разработаны и исследованы алгоритмы обнаружения азотосодержащих веществ на основе поглощения энергетического пика 9.17МэВ во вторичном гамма-излучении объекта. а. Алгоритмы имеют достаточно высокую эффективность, позволяя регистрировать малозаметные изменения интенсивности отдельных спектральных пиков. b. Простота решающей функции алгоритмов существенно упрощает их практическую реализацию, одновременно повышая быстродействие конечного устройства. c. Совместно с аппаратной частью эксперимента реализация работы алгоритмов предоставляет возможности создания автоматизированного досмотрового комплекса, не требующего специальных навыков обслуживающего персонала, что исключает влияние человеческого фактора при принятии решений. Операция получения спектра не требует отбора частиц и паров в пробоотборник, что также повышает устойчивость данного метода обнаружения к различным неблагоприятным факторам. Осуществлена проверка работы алгоритмов имитационным моделированием на ЭВМ, а также применением алгоритмов к реальным экспериментальным данным.

4. На основе применения следствий теоремы Лемана-Шеффе для оценок на основе полных достаточных статистик получены алгоритмы оценивания параметров интенсивности и ширины спектральных линий вторичного энергетического спектра гамма-излучения после прохождения через исследуемый объект. Моделирование алгоритмов оценивания показало их эффективность.

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

6. Заключение

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

1. Canadian Patent 2,086,054 A1, December 22,1992.

2. International Patent WO 2008/050327 A2, May 2, 2008.

3. US Patent US 2010/0038550 Al, February 18,2010.

4. International Patent WO 88/01788 Al, 10 March, 1988.

5. European Patent 0 218 240 A2, October 7, 1986.

6. Хильченко А.Д. Аппаратная инфраструктура измерительных и управляющих систем плазменных установок ИЯФ СО РАН. Диссертация на соискание ученой степени доктора технических наук. Новосибирск, 2010, С. 139-144.

7. Accelerator Mass Spectrometers. Dedicated and Multielement Systems. -High Voltage Engineering Europa B.V.

8. Казгов M.A., Калинин Б.А., Волобуев П.В. Интерпретация масс-спектров метастабильных ионов. Аналитика и контроль, 2009. Том 13, №4, С. 199-202.

9. Kuznetsov А.V., Osetrov O.I. Detection of Improvised Explosives (IE) and Explosive Devices (IED). Detection and Disposal of Improvised Explosives. -Springer, 2006. P.7-25.

10. US Patent US 2003/0193019 Al, October 16, 2003.

11. Pesente S., Lunardon M, Nebbia G., Viesti G. Detection of Improvised Explosives Devices (IED) by Using Tagged Neutron Beams. Detection and Disposal of Improvised Explosives. - Springer, 2006. P.69-85.

12. Reber E., Blackwood L., Edwards A. Idaho Explosives Detection System: Development and Enhancements. Sens Imaging, 2007. №8, P. 121-130.

13. International Patent WO 2007/084160 A2, July 26,2007.

14. Vakhtin D.N., Gorshkov I.Yu., Evsenin A.V., Kuznetsov A.V., Osetrov O.I. Senna Portable Sensor for Exposives Detection Based on Nanosecond Neutron Analysis. - Detection and Disposal of Improvised Explosives. - Springer, 2006. P.87-96.

15. Томилов A.B., Калинин Б.А., Александров O.E, Селезнёв В.Д.

16. Математическая обработка масс-спектра с не полностью разрешенными пиками. Аналитика и контроль, 2008. Том 12, №3-4, С. 107-112.

17. Перегудов О.Н., Покровский В.А., Рогульский Ю.В., Бугай А.Н. Метод определения параметров формы пика секторного масс-спектрометра. -Масс-спектрометрия, 2007. №4(1), С. 43-48.

18. Гребенщиков О.А., Залесский В.Б. Оценка динамики дрейфа спектральных линий и компьютерная коррекция развертки масс-спектра. -Журнал прикладной масс-спектроскопии, 2004. Том 71, №3, С. 407-411.

19. De Hoffmann Е., Stroobant V. Mass Spectrometry. Principles and Applications. Third edition // John Wiley & Sons, 2007.

20. Разников B.B., Разникова M.O. Информационно-аналитическая масс-спектрометрия. -М.: Наука, 1991.-248 с. ISBN 5-02-001547-4.

21. Гилл Ф., Мюррей У., Райт М. Практическая оптимизация. М.: Мир, 1985.-509 с.

22. А. Г. Вострецов, А. В. Бурдаков, С. Е. Радченко, А. С. Кузнецов, Ю. С. Суляев. Метод обнаружения поглощения гамма-квантов припрохождении через азотосодержащее вещество. Автометрия, 2010. - Т. 46. -№ 3. - С. 22 - 29.

23. Vostretsov A.G., Radchenko S.E. Statistical Algorithm for Detection of Gamma Ray Spectral Peak Absorption. Proceedings of the 2009 international forum on strategic technologies. Vietnam: Ho Chi Minh, - 2009.

24. Васюков B.H. Введение в теорию сигналов: Учебное пособие. НГТУ - Новосибирск, 2001.

25. Фрэнке JI. Теория сигналов. М., «Сов. Радио», 1974.

26. Беклемишев Д. В. Курс аналитической геометрии и линейной алгебры. М., «Наука», 1987. — 319 с.

27. Радченко С.Е. Совместная оценка параметров спектральных линий энергетического спектра гамма-излучения. Современные проблемы радиоэлектроники: сб. науч. тр., 2011. - 563с.

28. V. J. Zeps, Е. G. Adelberger, A. Garcia, С. A. Gossett, Н. Е. Swanson, W. Haeberli, P. A. Quin, and J. Sromicki. Parity mixing of the 0+-0"/=l doublet in 14N.-Physical Review C, 1995. V.51,N.3.-PP. 1494-1520.

29. Корн Г., Корн Т. Алгебра матриц и матричное исчисление. Справочник по математике — 4-е издание. — М.: Наука, 1978. — С. 392— 394.

30. Беклемишев Д.В. Дополнительные главы линейной алгебры. М.: Наука, 1983.-С. 196.

31. Онучнн А.П. Экспериментальные методы ядерной физики: учеб. пособие. Новосибирск: Изд-во НГТУ, 2010. — 220 с.

32. Гельфанд И.М. Лекции по линейной алгебре. — М.: Наука, 1971.

33. Богданович В.А., Вострецов А.Г. Теорйя устойчивого обнаружения, различения и оценивания сигналов. М.: Физматлит, 2004. - 320 с. ISBN 59221-0505-8.

34. Lehmann E.L. Testing Statistical Hypotheses. N. Y: Wiley, 1959. (Русский перевод: Леман Э. Проверка статистических гипотез: Пер. с англ. Ю.В.Прохорова. - М.: Наука, 1979. - 408 е.).

35. Кендалл М.Дж., Стыоарт А. Теория распределений. М., «Наука», 1966.-588 с.

36. Левин Б.Р. Теоретические основы статистической радиотехники. Т.1 — М., «Сов. Радио», 1969. 751 с.

37. Гельфанд И.М., Шилов Г.Е. Обобщенные функции и действия над ними. — М., «Государственное издательство физико-математической литературы», 1959. 470 с.

38. Боровков A.A. Математическая статистика. Оценка параметров. Проверка гипотез. М., «Наука», 1984. - 472 с.

39. Сосулин Ю.Г. Теоретические основы радиолокации и радионавигации: Учеб. пособие для вузов. М.: Радио и связь, 1992. - 304 с.

40. Сквайре Дж. Практическая физика. М., «Мир», 1966. - 246 с.