Расчет быстродействующего канала гамма-датчика

Для того чтобы получить данные, инвариантные относительно общего светового потока, генерируемого при сцинтилляции, необходимо выработать также нормирующий сигнал U_z, пропорциональный сумме всех сигналов ФЭУ, т. е. ?q_i Определяя отношения U_x/U_z и U_y/U_z, можно вычислить координаты х и у независимо от энергии регистрируемых г-квантов и от статистического разброса амплитуд сигналов, вызванных одной и той же переданной энергией Ег. Неоднородность свойств ФЭУ и кристаллов, т. е. неоднородность чувствительности по площади СПЧД, может быть учтена при сложении сигналов или подбором напряжений питания ФЭУ (его усиления).

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

Центр системы координат можно выбрать произвольно. Если его расположить на краю системы ФЭУ и ограничиться только двумя суммирующими усилителями (по одному для U_х и U_у), то придется суммировать сигналы с выходов всех ФЭУ. Если же совместить его с центром симметрии системы, весовые коэффициенты для многих сигналов окажутся равными нулю, однако необходимы четыре усилителя (U_x⁺, U_x^-, U_y⁺, U_y^-) и последующее выделение модулей разностных сигналов (|U_x⁺- U_x^-| и | U_y⁺- U_y^-|).

Приведена структурная схема такого СПЧД, содержащего 19 ФЭУ (их условная нумерация показана там же). Сигналы с выходов ФЭУ через резисторы или конденсаторы, задающие весовые коэффициенты f_i, суммируются на входах усилителей У_х+, У_х-, У_у+ и У_у+; кроме того, сигналы с выходов всех ФЭУ суммируются на входе усилителя У_г. Дифференциальными усилителями, ДУ_x и ДУ_у выделяется разность усиленных сигналов, и выходные импульсы затем поступают для обработки непосредственно (при |U⁺|>|U^-|) или после инверсии (при |U+|<|U^-|). Выполнение того или иного условия определяют дискриминаторы Д_х и Д_у. При помощи каскадов аналогового деления амплитуд сигналов КАД_Х и КАДу осуществляется нормировка амплитуд выходных сигналов U_х и U_у, которые пропорциональны координатам сцинтилляции х и у.

Расчеты показывают, что оптимальные характеристики достигаются, когда край кристалла располагается у центров фотокатодов периферийных ФЭУ. Минимум позиционной нелинейности соответствует случаю, когда (h_сц+h_св)/d=0,45, где d - диаметр фотокатода ФЭУ, и приемлема удаленность сцинтилляций от плоскости фотокатода в пределах от 0,35d до 0,55d. Следовательно, h_сц =0,2d, a h_св =0,35d.

Позиционное разрешение рассмотренных СПЧД в значительной степени определяется флуктуациями числа первичных фотоэлектронов и относительно невелико. Особенно это сказывается при низких энергиях г-квантов (менее 0,1 МэВ). Система суммирования сигналов с весовыми коэффициентами, используемая в рассмотренных СПЧД, не оптимальна, поскольку значительный вклад в позиционный сигнал дают ФЭУ, воспринимающие малую долю световой вспышки и несущие малую информацию, т. е. вырабатывающие сигналы с амплитудой, подверженной большому статистическому разбросу. Позиционное разрешение монотонно ухудшается с увеличением значений h_сц+ h_св, при (h_сц + h_св)/d=0,45 оно близко к 0,3% в центре сцинтиллятора и доходит до 0,85% у его периферии. Сужение зоны амплитудного отбора сигнала U_z улучшает позиционное разрешение, но это приводит к снижению интегральной чувствительности детектора.

Дополнительное ухудшение позиционного разрешения может быть вызвано при многократных процессах взаимодействия г-квантов со сцинтиллятором (например, комптоновское рассеяние и затем фотопоглощение рассеянного кванта, так что в кристалле происходят две сцинтилляции на некотором расстоянии одна от другой). Чтобы снизить влияние этих эффектов, приходится ограничиваться энергиями г-квантов не превосходящими 0,4-0,5 МэВ.

Рассмотренные СПЧД используют в приборах, получивших название гамма-камер. Такую аппаратуру в последние годы выпускают многие фирмы, например «Nuclear Chicago», «Ohio Nuclear» (США), «Toshiba» (Япония) и др. В типичной гамма-ка- мере [42], описанной в литературе, были использованы монокристалл NaI(Tl) с D = 293 мм и h_cц=12,7 мм, световод с h_cв = 38 мм и 19 ФЭУ диаметром 80 мм.

СПЧД аналогового типа с временным преобразованием сигналов содержат, как и рассмотренные выше устройства, монокристалл NaI(Tl), поверхность которого перекрыта мозаикой ФЭУ с диаметром d и расстоянием между центрами фотокатодов l_Х и l_У. Однако для формирования координатных сигналов используется суммирование сигналов ФЭУ на отводах линий задержек, специальное формирование биполярных сигналов и выделение временных интервалов между моментом действия г-кванта и моментами пересечения сигналами нулевого уровня [17, 11, 6, 10].

В таких СПЧД сигналы от ФЭУ поступают на отводы двух линий задержек ЛЗ_X и ЛЗ_У причем точки подключения ФЭУ к линиям задержек зависят от координат центров фотокатодов (ФЭУ, центры фотокатодов которых имеют одну и ту же координату, подсоединены к одному отводу). Время распространения сигнала по линии задержки между двумя отводами Дt_x или Дt_y такое, что Дt_xv = l_x и Дt_yv = ly, где v -- скорость распространения сигнала. Следовательно, время распространения сигналов от точки, где возникла сцинтилляция, к краям линий, пропорционально расстояниям от этой точки до краев кристалла, а разница во времени пропорциональна координатам x и y.

Сигналы с выходов линий задержек ЛЗ_X и ЛЗ_У поступают на каскады БФ, осуществляющие биполярное формирование, а затем после усиления каскадами У_х+, У_х-, У_у+, У_у- дискриминаторы Д_х+, Д_х-_, Д_у+, Д_у- выделяют моменты t_x+, t_x-, t_y+, t_y- пересечения сигналами нулевого уровня. Кроме того, сигналы, непосредственно снимаемые со всех ФЭУ, суммируются на усилителе У_z, отбираются по амплитуде каскадом Д_г и выделяется момент t₀ регистрации г-кванта в сцинтилляторе. В работе [10] показано, что при оптимальном формировании интервал |t_x+ -- t_x- | пропорционален координате |х| события, а интервал \t_y+--t_y- | -- координате |у|. Для этого каскадами время-амплитудного преобразования ВАП_Х+, ВАП_Х-, ВАПу+ , ВАП_у- интервалы (t_x+ --t₀) и другие преобразуются в пропорциональную им амплитуду сигналов, а дифференциальными усилителями ДУ_Х и ДУ_У выделяются разностные сигналы U_х и U_у, амплитуда которых пропорциональна координатам.

Каскады БФ снабжают средством регулировки (в достаточно широких пределах) параметров импульсной переходной характеристики для приближения формы сигналов к оптимальной. Это позволяет получить удовлетворительные позиционные линейность и разрешение в широком диапазоне значений h_сц и h_cв. Важно, что при выделении координатных сигналов момент пересечения нулевой линии не зависит от амплитуды импульсов; это позволяет исключить нормирование сигналов по амплитуде импульса U_z. Однако временное преобразование требует дополнительного времени (обычно несколько микросекунд) и устройство должно содержать режектор наложений импульсов для того, чтобы уменьшить искажения при высокой частоте сигналов.

На базе таких СПЧД выполнено несколько сцинтилляционных гамма-камер, отличающихся размером монокристалла, числом и способом размещения ФЭУ. Типичное такое устройство содержит кристалл NaI(Tl) с D = 387 мм и h_cц=12,7 мм и 30 ФЭУ с d=76 мм и l_x=1_У=87 мм [11]. Позиционное разрешение СПЧД с временным формированием сигнала оказывается примерно в 1,5 раза лучше, чем детекторов с линейным суммированием импульсов ФЭУ. Данный метод широко применяется в однофотонный эмиссионный гамма-томограф „ЭФАТОМ“ и ПЭТ.

3.3 Двухдетекторный однофотонный эмиссионный компьютерный томограф „ЭФАТОМ“

3.3.1 Описание

Двухдетекторный однофотонный эмиссионный компьютерный томограф „ЭФАТОМ“ предназначен для медицинских диагностических исследований внутренних органов и систем человека на основе визуализации распределения фармацевтических препаратов, меченных гамма-излучающими радионуклидами (радиофармпрепаратов). Результатом работы томографа являются одноканальные цифровые изображения, представляющие собой проекции трехмерного распределения радиофармпрепарата на плоскость детектора [1,2]. В НИИЭФА им. Д.В. Ефремова в период 2003?2007 гг. был разработан и изготовлен опытный образец гамма-томографа [14,27]. Технические и клинические испытания гамма-томографа были проведены в 2007?2008 гг. в клинической больнице № 83 (Москва).

3.3.2 Устройство блока детектирования томографа

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

Излучаемые исследуемым объектом гамма-кванты, прошедшие через коллиматор, взаимодействуют со сцинтилляционным кристаллом, в результате чего возникают сцинтилляционные световые вспышки. Интенсивность вспышки пропорциональна энергии, которую гамма-квант оставил в сцинтилляторе. Ампилтуды импульсов с анода ФЭУ зависят от интенсивности световой вспышки и от ее положения. Эту зависимость называют амплитудно-пространственной характеристикой (АПХ. Далее каждую сцинтилляцию регистрируют сразу несколько ФЭУ. Сигнал с ФЭУ, прошедший через предусилитель, поступает на отдельный аналого-цифровой преобразователь. В результате электронная схема оцифровывает входные сигналы с ФЭУ, вычитает базовую линию, проверяет попадание сигнала во временное окно и передает данные по высокочастотной шине данных в компьютер сбора. Компьютер сбора принимает цифровые данные, которые сгруппированы в пакеты по 56 байт, с платы электроники, проводит передает вычисленные координаты в центральный компьютер по протоколу UDP через интерфейс Ethernet.

3.3.3 Вычисление координат событий

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

3.3.4 Энжеровский алгоритм с коррекцией нелинейности

Для определения координат по классическому методу Энжера используются сигналы со всех ФЭУ, и вычисление производится по форме взвешенного среднего [5]. Координаты центров ФЭУ суммируются с весами, пропорциональными амплитудам сигналов ФЭУ:

где X, Y -- оценка координаты события, Xci, Yci -- координаты центров ФЭУ, Ui = Ui(x, y) -- АПХ соответствующего ФЭУ, N -- общее число ФЭУ. Реальная АПХ нелинейная, поэтому вычисленные координаты имеют нелинейные искажения. Для проведения коррекции нелинейности необходимо составить таблицу соответствия энжеровских координат реальным координатам события. Эта таблица строится по информации, полученной в процессе калибровки детектирующей головки. Калибровка заключается в прохождении коллимированного точечного источника по координатной сетке, перекрывающей всю площадь кристалла, и определении отклика всех ФЭУ в каждой точке соответственно. Шаг сетки, как показал проведенный расчет и экспериментальные результаты, должен быть не более 4 mm. Основной недостаток классического энжеровского алгоритма -- это ухудшение собственного пространственного разрешения с увеличением количества ФЭУ детектирующей головки. Это связано с „шумом“ ФЭУ, находящихся далеко от точки, в которой произошло сцинтилляционное событие, в результате чего точка на картине размывается и ухудшается пространственное разрешение. Также классический алгоритм не позволяет регистрировать события за границами области центров крайних ФЭУ [5].

3.3.5 Кластерный энжеровский алгоритм

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

где X, Y -- оценка координаты события, Xci, Yci -- координаты центров ФЭУ, Ui -- амплитуды сигналов ФЭУ, Nc -- число ФЭУ в кластере. Для реализации этого алгоритма необходимо предварительно рассчитать таблицы перехода от кластерных энферовских координат к реальным координатам для каждого кластера, т. е. найти реальные координаты сцинтилляционного события, соответствующие кластерным энжеровским координатам. Расчет этих таблиц проводится с помощью минимизации квадратичного функционала вида:

f (x, y) =(XA ? XA(x, y))² + (YA ? YA(x, y))²,

где x, y -- реальные координаты, в которых произошло сцинтилляционное событие, XA, YA -- табличные энжеровские значения, XA(x, y), YA(x, y) -- 2D-сплайны зависимостей энжеровских кластерных координат от реальных координат. Приемлемая точность достигается при шаге по энжеровским координатам -- 0.08 mm. Достоинства описанного алгоритма заключатся в хороших результатах по пространственному разрешению и нелинейности, а также в высоком быстродействии.

3.3.6 Амплитудно-пространственные характеристики

Для построения таблиц соответствия энжеровских координат реальным были проведены измерения АПХ для всех ФЭУ в сборке детектирующей головки. Измерения производились с помощью коллимированного точечного источника гамма-квантов Co-57 с энергией 122 keV. Источник позиционировался в точке с заданной координатой на поверхности детектора, в которой определялись средние отклики от каждого ФЭУ. Далее источник передвигался по одной из координат на один шаг, равный 4 mm.

ЗАКЛЮЧЕНИЕ

1. Как видно использование отражателей не улучшает показателей длительности импульса на 50% амплитуды t_mеаn, нс, а в некоторых экспериментах даже ухудшает его. Следовательно, можно сделать вывод, что не целесообразно использовать отражатели в установках данного типа.

2. Использование в экспериментах коллиматора значительно (приблизительно в 1.8 раза) улучшает показатели по длительности импульса на 50% амплитуды t_mеаn, нс. Так же использование коллиматора обеспечивает защиту детектора от внешнего излучения, вследствие чего увеличивается динамический диапазон, что позволяет получать более точные измерения.

3. Проводились исследования зависимости влияния на ВАХ площади окна коллимационной головки для сцинтиллятора ПС-Н2. Показатели по длительности импульса на 50% амплитуды t_mеаn, нс лучше там, где использовалась коллимационная головка с меньшей площадью окна пропускания.

4. Без использования коллиматора видна зависимость длительности импульса на 50% амплитуды t_mеаn, нс от различных объемов сцинтиллятора. Можно сделать вывод, что чем меньше сцинтиллятор, тем лучше показатели длительности импульса на 50% амплитуды. Но уменьшать объем сцинтиллятора бесконечно нельзя, так как тогда на него будет попадать меньше частиц и, следовательно, придется увеличивать время проведения эксперимента.

5. Проводились испытания по влиянию на ВАХ формы окна коллимационной головки. Были использованы два коллиматора с круглым окном (площадь окна составляет S=254,34мм²) и квадратным окном (площадь окна составляет S=283,24мм²). Из исследований видно, что ВАХ при использовании круглого окна лучше, чем при использовании квадратного окна.

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

7. Использование многоканальных схем поможет улучшить производительность контроля объектов. Так же появится возможность сопоставления наличия дефекта в контролируемом образце с местонахождением данного дефекта.

СПИСОК ЛИТЕРАТУРЫ

1. Anger H. O. -- IEEE Trans., 1966, v. NS--13, JVfe 3, p. 380--386.

Anger H. O. Scintillation camera. -- Rev. Scient. Instrum., 1958, v. 29, № 1, p. 27--32.

2. Backer R. G., Scrimger J. W. -- Phys. Med. and Biol., 1967, v. 12, № 1, p. 51--56.

3. Charpak G., Dick L., Feuvrais L. -- Ibid., 1962, v. 15, № 4, p. 323--326.

4. Grawley J. C. W., Veall N. -- In: Medical radioisotope scintigraphy, 1972, v. 1, Vienna, IAEA, 1973, p. 105--112.

5. Grosshog G. -- Nucl. Instrum, and Methods, 1975, v. 127, No 1, p. 183--1I88.

6. Hiramoto Т., Tanaka E., Nohara N. -- J. Nucl. Med., 1971, v. 12, No 4, p. 161--168.

7. http://www.greenstar.ru/articles/obzor_gama_cpektorv.pdf

8. Lee L. J., Allred J. C., Goodman C. -- Ibid., 1974, v. 119, № 1, p. 29--33.

9. Miraldi F., Yon E. Т., Rejaly A. -- In: Medical radioisotope scintigraphy, 1972, v. 1, Vienna, IAEA, 1973, p. 121--132.

10. Tanaka E., Hiramoto Т., Nohara N. -- J. Nucl. Med., 1970, v. 11, No 9, p. 542--549.

11. Tanaka E., Nohara N., Kumano N. e. a. -- In: Medical radioisotope scintigraphy, 1972, v. 1, Vienna, IAEA, 1973, p. 169--180.

12. А.Н. Майоров, С.В. Мамиконян, Л.И. Косарев, В.Г. Фирстов, Радиоизотопная дефектоскопия (методы и аппаратура), М., Атомиздат, 1976, 208с.

13. Балдин Б. Ю., Вертоградов J1. С., Крумштейн В. В. и Др* Координатные детекторы на жидких сцинтилляторах. -- Приборы и техника эксперимента, 1978, № 3, с. 78--83.

14. Бондаренко В. Г., Григорьев В. А., Каплин В. А. Эффективная скорость распространения светового сигнала в сцинтилляцион- ном координатно-чувствительном детекторе. -- Приборы и техн. эксперимента, 1979, № 3, с. 75--77.

15.Бондаренко В. Г., Григорьев В. А., Макляев Е. Ф. и др. Сцинтилляционный координатно-чувствительный детектор с высоким временным разрешением. -- Приборы и техника эксперимента, 1976 № 1, с. 44--46.

16.Вяземский В. О., Ломоносов И. И., Писаревский А. И. и др. Сцинтилляционный метод в радиометрии. М.: Атомиздат, 1961.

17.Горн Л. С., Костылев В. А., Наркевич Б. Я. и др. Приборы для радиоизотопной диагностики в медицине. М., Атомиздат, 1978.

18. Калашникова В. И., Козодаев М. С. Детекторы элементарных частиц. М.: Наука, 1966.

19.Колюбин В.А., Косарев Л.И. и др., Методы и аппаратура томографического, интроскопического и радиометрического контроля качества продукции, Вопросы атомной науки и техники, серия техническая физика и автоматизация, выпуск 62, 2007, с. 69-81.

20. Кузелев, Е.А. Жуковский, В.Н. и др., Современное оборудование для радиографии Химическое и нефтегазовое машиностроение, №4, 2004, С.44-48.

21. Льюис В., Уэлс Ф. Миллимикросекундная импульсная техника. М.: Изд-во иностр. лит., 1956.

22. Марков А. А. Электронные методы экспериментальной физики. Ч. 2. М.: Изд. МИФИ, 1956.

23. Матвеев В. В., Соколов А. Д., Фотоумножители в сцинтилляционных счетчиках. М.: Госатомиздат, 1962.

24. Неразрушающий контроль, справочник в 7 томах под ред. В.В. Клюева, книга 2 «Радиационный контроль», М., «Машиностроение», 2003.

25.Современная аппаратура для неразрушающего контроля качества сварных соединений трубопроводов радиографическим методом, Тяжелое машиностроение, 9, 2005, с. 8-10.

26. Сотников С. К., Ефимов Б. В., Цитович А. П. Метод стабилизации тракта усиления сцинтилляционного счетчика, -- ПТЭ, 1965, № 1, с. 100-- 104.

27.Технические средства и методики автоматизированного неразрушающего контроля сварных соединений на основе цифровых систем регистрации рентгеновского и гамма излучения, Тяжелое машиностроение, 4, 2005, с. 4-6.

Размещено на Allbest.ru