Способы и базы данных для идентификации нуклидов

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

 

Настоящее изобретение относится к способам и базам данных для идентификации нуклидов.

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

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

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

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

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

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

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

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

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

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

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

Согласно предпочтительному осуществлению способ также содержит этап искусственной реконструкции измеренного первого спектра µ(Е).

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

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

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

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

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

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

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

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

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

Предпочтительно, чтобы модельные спектры из матриц откликов на поглощение были объединены в качестве векторов-столбцов в матрице F.

Также предпочтительно получать модельный спектр для каждой матрицы откликов на рассеяние.

Согласно еще одному предпочтительному осуществлению модельные спектры из матриц откликов на рассеяние объединяют в качестве векторов-столбцов в матрице М.

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

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

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

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

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

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

Краткое описание чертежей

На чертежах:

фиг.1 - иллюстрация осуществления настоящего изобретения;

фиг.2А - график, показывающий типичное реальное измерение европия Eu-152 при использовании NaI детектора;

фиг.2В - график, показывающий пример двух этапов процедуры моделирования;

фиг.3 - общий вид различных матриц откликов на поглощение;

фиг.4 - вид строк матрицы характеристик;

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

фиг.6 - вид реконструированного синтетического спектра в сравнении с измерением и полученным стандартным способом моделированием;

фиг.7 - результат из уравнения (6) для эталонного спектра европия;

фиг.8 - результат измерения для йода I-131 при использовании источника, заключенного в алюминиевый и кадмиевый кожух; и

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

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

- оценивание по меньшей мере одного источника излучения по измеренному спектру µ(Е) с помощью заданного гамма-спектроскопического прибора;

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

- определение на основании сравниваемых сценариев информации α, относящейся к материалу x, которым источник излучения окружен.

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

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

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

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

Матрица R откликов представляет собой поле m×n чисел, которое связывает ответ детектора в канале i (где i∈{1, …, m}) с энергией излучения. Последняя представляет собой энергию падающего излучения, которая возбуждает материал детектора. Падающее излучение, которое обычно является частицей с энергией E_Incident, может испытывать а) полное поглощение в материале детектора, при этом полная энергия E_Detected=E_Incident остается в материале детектора, b) частичное поглощение, тем самым количество E_Detected поглощенной энергии меньше E_Incident, или с) прохождение через детектор без поглощения энергии. Кроме того, не исключено отражение энергии, которая проходит через детектор, окружающим материалом, что влечет за собой отдельный цикл поглощения. Поэтому матрица R откликов является представлением физически поглощаемой энергии в детекторе. Иначе говоря, матрица R откликов представляет преобразование от исходных энергий излучения к спектральной картине, создаваемой детектором.

В этом осуществлении изобретения матрицы откликов моделируют с использованием кодов переноса излучения, описанных в Geant4. Geant4 является инструментальным программным средством для моделирования прохождения частиц через вещество и описано в S. Agostinelli et al., Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, Volume 506, Issue 3, 1 July 2003, Pages 250-303, а также в Allison, J., et al., IEEE Transactions on Nuclear Science 53, № 1 (2006), 270-278. В этих кодах подход Монте-Карло используется для воспроизведения траекторий частиц и волн. На основании геометрической и физической модели таким способом моделируют различные сценарии поглощения и рассеяния в детекторе.

Согласно этому осуществлению изобретения следующую процедуру выполняют для определения матрицы R откликов.

Сначала задают распределенный источник с моноэнергетическим линейчатым излучением, имеющим энергию E_i.

Начиная с канала i=1, моделируют излучение энергии Е₁=3 кэВ, для этой единственной энергии образуют модельный спектр R₁(E). Такая функция R₁(E₁) представляет собой гистограмму интервала [E₁=1 кэВ, E_N=3072 кэВ] энергий.

Затем энергию излучения повышают в соответствии с дискретными интервалами dE энергий и моделирование повторяют, чтобы образовать дополнительные спектральные функции для каждого дополнительного канала R_i(E), где i=1, 2, …, 1024, при этом энергия излучения Е₁₀₂₄ составляет 3072 кэВ.

Таким образом, при использовании этой процедуры мы получаем N=1024 функций R_i, где i=1, …, 1024, которые представляют отклики системы на энергии Е=1 кэВ, …, 3072 кэВ падающего излучения. Эти функции можно представить в виде матрицы с конкретными размерами N×N, где N=1024:

$$\text{R}=\left(\begin{array}{cccc}{R}_1(E_1)& R_2(E_1)& \mathrm{...}& R_1(E_1)\\ R_1(E_2)& R_2(E_2)& \mathrm{...}& R_1(E_2)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ R_1(E_N)& R_2(E_N)& \mathrm{...}& R_N(E_N)\end{array}\right)$$ . $$$$ (1)

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

После того как энергии Е излучения и интенсивности ω_i становятся известными, матрицу R откликов можно использовать для образования спектра δ:

$$\delta (E)={\displaystyle \sum _i{\omega }_i{R}_i(E)}$$ . $$$$ (2)

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

ζ(Е)=δ(Е)⊗ρ(Е). $$$$ (3)

В качестве варианта столбцы R могут быть свернуты до применения уравнения (2).

На фиг.2В показано моделирование, основанное как на несвернутых данных, так и на спектре, образованном с включением собственной разрешающей способности ρ детектора. В частности, на фиг.2В представлен пример двух состояний моделирования: верхний спектр, обозначенный (а), является спектром европия Eu-152, моделированного с использованием (идеального) дельтаобразного отклика детектора наряду с данными об энергии и интенсивности из базы данных. Этот искусственный спектр был получен в соответствии с уравнением (2). Такой же спектр, обозначенный (b), теперь включает в себя собственную разрешающую способность детектора в соответствии с уравнением (3). Следует отметить, что все же имеются различия между свернутым модельным спектром (b) на фиг.2В и реальным измерением, отображенным на фиг.2А, которое показывает типичное измерение европия Eu-152 с использованием детектора из NaI (иодида натрия), где в этом спектре показаны 512 каналов, при этом каждый канал отражает энергию 3 кэВ, в результате чего охватывается диапазон энергий от 0 кэВ до 1536 кэВ.

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

Следуя процедуре моделирования, изложенной выше, используя поглотители, создают набор откликов А₀, А₁, А₂, А₃, А₄, …, А₂₁, А₂₂ на поглощение, которые охватывают весь диапазон ослабления от 100%-ного пропускания до пропускания меньше 1%.

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

В качестве примера показана третья матрица поглощения:

$$A_3=\left(\begin{array}{cccc}{A}_{31}(E_1)& A_{32}(E_2)& \mathrm{...}& A_{3N}(E_N)\\ A_{31}(E_2)& A_{32}(E_2)& \mathrm{...}& A_{3N}(E_2)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ A_{31}(E_N)& A_{32}(E_N)& \mathrm{...}& A_{3N}(E_N)\end{array}\right)$$ . $$$$ (4)

Заметим, что двойное индексирование применено для различения вида ослабления и столбцов.

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

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

Кроме того, в осуществлении используют линейный каталог, содержащий все необходимые положения ε_i пика и интенсивности E_i пика, так что формулу (2) можно использовать для построения искусственных спектров для всех нуклидов в базе данных. Данные линейного каталога основаны на файле оцененных данных по структуре ядра (ENSDF) (дополнительные сведения см., например, в G. Audi et al., The 1995 update to the atomic mass evaluation, Nucl. Phys., 595; 409-480, 1995).

Применяя формулы (2) и (3), получаем модельные спектры для каждой матрицы откликов на поглощение, обозначенные ξ_А0, ξ_А1, …, ξ_А22:

$${\xi }_{Ax}=\left({\displaystyle \sum _i{\omega }_i{A}_{xi}}\right)\otimes \rho (E)$$ , $$$$ (5)

при этом индекс x обозначает число матриц поглощения. Модельные спектры объединяют в качестве векторов-столбцов в новой матрице:

$$\tilde{F}=({\xi }_{A0},{\xi }_{A1},\mathrm{...},{\xi }_{A22})$$ . $$$$ (6)

В качестве варианта вместо использования моделирования можно использовать набор спектров, составленных из измерений. Однако преимущество этого осуществления заключается в том, что отсутствует необходимость в предшествующих калибровочных измерениях, поскольку эта часть композиции поглощений заменяется моделированием. На следующем этапе все содержимое матрицы рассеяния добавляют к матрице $$\tilde{F}$$ характеристик

$$F=\left(\tilde{F},M_1,M_2,\mathrm{...},M_N\right)$$ , $$$$ (7)

где M₁, …, M_N представляют собой векторы-столбцы матрицы М рассеяния. Результирующая матрица F характеристик имеет размер (N+23)×N. Для минимизации количества строк, необходимых для описания рассеяния, эффекты рассеяния учитывают путем выбора 10 сортированных вкладов рассеяния из матрицы М со следующими энергиями: 50 кэВ, 100 кэВ, 200 кэВ, …, 900 кэВ, 1000 кэВ.

На фиг.4 показаны элементы матрицы F характеристик для Eu-152, нуклида с многочисленными линиями, который широко известен как тестовый источник. Показаны строки матрицы F характеристик, при этом каждая строка имеет длину 1024. F₁ отражает спектр почти без экранирования. F₁₀, F₁₅, F₂₀ представляют примеры различного экранирования, и F₂₄, F₂₅, F₂₆ представляют три различных следа рассеяния. В этом случае матрица F характеристик представляет собой матрицу 26×1024. Графическое представление полной матрицы F можно обнаружить на фиг.5, на которой представлена так называемая водопадная диаграмма матрицы F характеристик. Первые десять столбцов отражают рассеяние, последующие 22 столбца - варианты различного экранирования европия. Поглощение изменяет спектр путем снижения интенсивностей линий с низкими энергиями. Как можно видеть из срезов F₁₀ и F₁₅, если поглощение является достаточно сильным, линии удаляются из спектра вследствие этого эффекта. Следовательно, этот баланс интенсивностей является различным для всех 23 различных случаев поглощения. Дополнительное рассеяние также показано на фиг.4. Рассеяние действует как большой сгусток непрерывного излучения, расположенный также в области низких энергий и зависящий от энергии рассеиваемого излучения.

Теперь будет описано, каким образом матрицу F характеристик используют при реальных измерениях излучения, например в сценариях контроля. При этом задача для матрицы F характеристик заключается в получении новой информации об измеренном спектре µ(Е). Для этого находят наиболее вероятное сочетание строк матрицы характеристик, которое обеспечивает наилучшее соответствие измеренному спектру µ(Е). Математически должна быть решена задача инверсии. Несколько методов можно использовать для оценивания этого наиболее вероятного сочетания: стандартные методы наименьших квадратов, инверсию методом наименьших квадратов с неотрицательными решениями или генетические алгоритмы.

В одном осуществлении изобретения реализуется инверсия методом наименьших квадратов с неотрицательными решениями. Принимается, что в своей основе оператор $$\widehat{L}$$ должен быть оператором метода наименьших квадратов с неотрицательными решениями, которым решается задача инверсии матрицы F характеристик:

$$\alpha =\widehat{L}[F]\mu$$ , $$$$ (8)

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

Дальнейшие детали источника, представленные значениями α, можно видеть на фиг.7, на которой результат для вектора α показан применительно к эталонному спектру европия Eu-152. На фиг.7 можно видеть, что характеристика небольшого поглощения и три эффекта интенсивного поглощения обнаруживаются вплоть до характеристики под номером 20. Рассеяние образует характеристики от 24 до 34, на которых проявляется сочетание различных вкладов рассеивателя. Из этого результата можно вывести, что а) имеется смесь поглотителей, при этом с высокой вероятностью сильным поглотителем является кадмий или свинец. Для лучшего понимания этой интерпретации данные из таблицы на фиг.3 следует сравнивать с соответствующими коэффициентами α из фиг.7. Первые 23 строки отражают влияния материалов A_i. Строки с 24 по 34 соответствуют влиянию рассеяния от 50 кэВ до 1000 кэВ, а что касается строк 27, 28, 29, то в них проявляется максимальный вклад, при этом рассеяние определяется сложным составом рассеиваемых фотонов в диапазоне энергий от 300 кэВ до 500 кэВ.

Затем синтетический (модельный) спектр s(E) реконструируют на основании параметра α:

$$s(E)=F\alpha$$ . $$$$ (9)

Этому синтетическому спектру s(E) следует придать наилучшее соответствие измерению при использовании математических базисных векторов, охваченных F. Конечно, наиболее важная проверка этого моделирования заключается в непосредственном сравнении с первоначально измеренным спектром µ(Е). Это показано на фиг.6, на которой для сравнения приведены реконструкция s(E), начальное измерение µ(Е) и стандартное моделирование st(E). Обычно способом согласно осуществлению изобретения измерение воспроизводится с большей детальностью. Пиковые отношения интенсивностей разрешаются намного лучше, чем при стандартном моделировании.

Ниже показываются несколько примерных случаев, в которых применяется способ. На фиг.8 показано измерение йода I-131, при этом источник был заключен внутри алюминиевого и кадмиевого кожуха, а на фиг.9 показано измерение йода I-131, когда в противоположность реконструкции и стандартному моделированию источник был размещен в модели тела человека. Йод I-131 представляет собой медицинский изотоп, широко используемый для радиологической диагностики и лечения. Здесь обсуждаются два случая: для фиг.8 радиоактивный источник располагали в цилиндрической оболочке, образованной из алюминиевых и кадмиевых слоев. Из спектра материала можно видеть, что поглотитель непосредственно вытекает из проведенного анализа. Кроме того, только небольшая степень рассеяния связана с источником. В отличие от этого на фиг.9 показан I-131 в таком же алюминиевом контейнере, но дополнительно внутри модели тела человека. Фактически, рассеяние излучения сильно возрастает, что можно видеть из области низких энергий. Реконструкция дает спектр α материала с ясно выраженным рассеянием. Алюминиевый поглотитель все еще обнаруживается, но также возникает дополнительное небольшое поглощение от материала модели. Поэтому на основании информации α о материале, которая приводит к приведенным выше результатам, можно легко сделать вывод, что медицинские источники на самом деле были заключены в теле человека, создающем высокую степень рассеяния, или в стальном или алюминиевом цилиндре.

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

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

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

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

1. Способ идентификации по меньшей мере одного источника излучения, в частности по меньшей мере одного нуклида, который по меньшей мере частично окружен заданным материалом, при этом способ содержит следующие этапы, на которых:
- оценивают по меньшей мере один источник излучения по измеренному спектру (µ(E)) посредством заданного гамма-спектроскопического прибора;
- сравнивают измеренный спектр со спектрами поглощения и спектрами рассеяния из базы данных, содержащей данные, представляющие спектры для множества сценариев поглощения и для множества сценариев рассеяния, относящихся к упомянутому по меньшей мере одному оцениваемому источнику;
- определяют на основе сравниваемых сценариев информацию (α), относящуюся к материалу (x), которым по меньшей мере частично окружен источник излучения.

2. Способ по п.1, причем сценарий поглощения представляет вид и конфигурацию поглощающего материала (x), а сценарий рассеяния представляет вид и конфигурацию рассеивающего материала.

3. Способ по п.1, причем спектр поглощения описывают вектором (ξ_Ax), вычисляемым на основе матрицы (А_х) откликов, представляющей соответствующий поглощающий материал (x), и функции (ρ) отклика, представляющей гамма-спектроскопический прибор.

4. Способ по п.2, причем спектр поглощения описывают вектором (ξ_Ax), вычисляемым на основе матрицы (А_х) откликов, представляющей соответствующий поглощающий материал (x), и функции (ρ) отклика, представляющей гамма-спектроскопический прибор.

5. Способ по любому из пп.1-4, причем стандартная процедура идентификации нуклидов базируется на одном из способа поиска пика и способа сопоставления с образом.

6. Способ по п.5, причем данные, содержащиеся в базе данных, представляющие спектры, находятся в виде матриц (F) характеристик, при этом матрица (F) характеристик образована множеством векторов (ξ_A0, ξ_A1, …, ξ_A22), описывающих спектры поглощения, и строками матрицы (М) рассеяния, представляющей рассеивающий материал.

7. Способ по п.6, причем информацию (α) определяют находя одно сочетание строк матрицы (F) характеристик, которое обеспечивает наилучшее соответствие измеренному первому спектру (µ(Е)).

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

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

10. Способ по любому из пп.1-4 или 6-9, причем данные, представляющие спектры, получают с помощью процедур моделирования.

11. Способ по п.5, причем данные, представляющие спектры, получают с помощью процедур моделирования.

12. Способ по любому из пп.1-4, 6-9 или 11, причем гамма-спектроскопический прибор представляет собой сцинтилляционный детектор излучения, в частности ручной детектор.

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

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

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

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

17. Способ по п.15, причем для по меньшей мере одной матрицы (М) откликов на рассеяние получают модельный спектр.

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

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

20. Способ по п.18 или 19, причем в ходе процедуры моделирования моделируют множество матриц откликов на рассеяние.

21. Способ по п.20, причем множество сценариев поглощения и рассеяния моделируют на основании физической и геометрический модели гамма-спектроскопического прибора.

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

23. Способ по любому из пп.1-4, причем данные, содержащиеся в базе данных, представляющие спектры, находятся в виде матриц (F) характеристик, при этом матрица (F) характеристик образована множеством векторов-столбцов (ξ_A0, ξ_A1, …, ξ_А22) и столбцами матрицы (М) рассеяния, представляющей рассеивающий материал.

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