Явление ядерного магнитного резонанса (ЯМР) с момента своего открытия было положено в основу функционирования большого количества приборов, предназначенных для анализа состава и структуры веществ и материалов. Основной тенденцией развития этого направления было использование всё более сильных магнитов для увеличения разрешающей способности ЯМР-спектроскопии. Применение магнитов на основе сверхпроводящих соленоидов вывело ядерный магнитный резонанс в число базовых аналитических технологий 21 столетия, однако стоимость таких приборов оказывается чрезвычайно высокой.
Развитие микроэлектроники и появление новых магнитожёстких материалов в последние годы выявило новые возможности для создания малогабаритных и достаточно дешёвых ЯМР-анализаторов, пригодных для широкого применения в самых разных областях хозяйственной деятельности человека, включая научные исследования. Это уже позволило, наряду c уменьшением размеров и массы, приблизить возможности низкополевых приборов к ЯМР-спектрометрам в сильных полях по целому ряду показателей. В настоящее время работы, направленные на миниатюризацию ЯМР приборов контроля, активно продолжаются [1–3]. В настоящей работе представлен вариант генератора последовательностей радиочастотных импульсов, основанный на применении алгоритма прямого цифрового синтеза и современной интегральной схемотехники.
Особенности импульсного ЯМР и требования к аппаратуре
Основными узлами ЯМР-релаксометра являются: источник постоянного магнитного поля H0, генератор радиочастотных импульсов; датчик, состоящий из катушки, включенной в состав колебательного контура, который используется как для создания импульсного магнитного поля H1 в анализируемом образце, так и для приёма слабых сигналов ЯМР; усилитель и детектор сигналов ЯМР. Имеется также устройство управления, как правило взаимодействующее с компьютером.
Сформулируем требования для радиочастотного передатчика ЯМР-релаксометра [4, 5], в состав которого входит генератор последовательности радиочастотных импульсов. Основное требование – способность создавать короткие радиочастотные импульсы высокой мощности в микросекундном диапазоне длительностей. Требуемая мощность зависит от объёма и свойств анализируемого образца и задаётся выходным усилителем передатчика. При этом время нарастания и спада огибающей импульса должно быть мало по сравнению с длительностью импульса. Частота несущей радиочастотного импульса должна быть равна резонансной частоте ядер атомов исследуемого вещества, которая определяется как частота прецессии Лармора fл.
При анализе образцов в ЯМР-релаксометрах используются импульсные последовательности, которые накладывают требования на длительность и амплитуду импульсов, поворачивающих вектор ядерной намагниченности M на заданный угол (90 °, 180 ° и др.) в зависимости от измеряемого показателя. Этот угол прямо пропорционален амплитуде ВЧ поля H1 и длительности импульса t, то есть, по существу, он определяется площадью импульса.
Обычно управляют длительностью импульса, так как её проще оперативно изменять. Следовательно, генератор передатчика должен иметь возможность формировать радиочастотные импульсы разной длительности с точной и управляемой выдержкой времени. А для методов на основе спинового эха требуется высокая точность временных интервалов и между импульсами.
Также в ряде методик [6, 7] требуются изменения фазы радиочастотных колебаний относительно опорной (начальной). То есть генератор передатчика ЯМР-релаксометра должен иметь возможность управления фазой радиочастотных колебаний.
В современных условиях передатчик ЯМР-релаксометра может быть построен полностью на основе аналоговой схемотехники или на основе цифровой и аналоговой схемотехники. Во втором случае на цифровую часть приходится генератор (синтезатор) последовательностей радиочастотных импульсов, а на аналоговую – фильтры и силовые выходные каскады. Такой подход позволяет исключить из конструкции некоторые аналоговые узлы, в частности для задачи разных длительностей импульса и для коммутации несущей сигнала. При использовании современных программируемых больших интегральных схем (БИС) это несомненно ведёт к уменьшению габаритов цифровой части конечного изделия, особенно если данная БИС будет основой и для других цифровых узлов прибора. Также существуют проблемы, связанные с уменьшением размеров аналоговых силовых узлов: передатчика и приёмника ЯМР-релаксометра, но их решение вынесено за рамки настоящей работы для будущих соисканий.
К цифровому генератору применяется ещё одно важное требование. Это наличие встроенного программируемого устройства управления, задающего параметры и последовательность импульсов. При этом устройство должно быть синхронизовано с опорной частотой генератора, а программирование должно осуществляться перед началом сеанса, так как управление непосредственно в ходе сеанса измерений может вносить дополнительную погрешность интервалов времени в импульсной последовательности.
Основные аспекты прямого цифрового синтеза
Рассматриваемый радиочастотный генератор может быть основан на алгоритме прямого цифрового синтеза (ПЦС) [8]. Классическая структурная схема такого генератора показана на рис. 1.
Рис. 1. Классическая структурная схема генератора на основе ПЦС
В ядро генератора входит аккумулятор фазы и преобразователь фаза-амплитуда. Аккумулятор фазы при каждом импульсе от опорного тактового генератора увеличивает своё значение на величину условного приращения фазы TUNE. Преобразователь фаза-амплитуда формирует на выходе значение, соответствующее входному коду фазы. Аккумулятор фазы должен переполняться при фазе 2π, то есть реальное значение фазы φ связано с кодированным значением фазы φcode на выходе аккумулятора выражением
, (1)
где Nаcc – двоичная разрядность аккумулятора фазы. Исходя из этого преобразователь фаза-амплитуда можно реализовать на основе таблицы соответствия между кодированной амплитудой и кодированной фазой. Данное соответствие устанавливается формулой
(2)
где N – двоичная разрядность выхода преобразователя фаза-амплитуда и соответственно цифро-аналогового преобразователя ЦАП. При этом в реальной таблице обычно кодируют только первый квадрант (от 0 до π/2), остальные квадранты получают путём манипуляции со знаком и косвенным индексом обращения к таблице. Управление частотой генератора возможно при помощи изменения параметра приращения фазы TUNE. Частота на выходе генератора зависит от TUNE в соответствии с выражением
(3)
где fт – опорная тактовая частота. То есть частоту можно задавать с шагом (fт/2Nасс) в пределах от 0 до fт/2, в соответствии с теоремой Котельникова.
Цифровой синтез программируемых последовательностей РЧ импульсов в генераторе передатчика ЯМР-релаксометра
На основе ранее приведенных требований к генератору передатчика ЯМР-релаксометра предложена усложнённая схема ПЦС, ориентированная на генерацию последовательностей радиочастотных импульсов различной длительности и фазы с возможностью программирования. Данная схема приведена на рис. 2.
Рис. 2. Структурная схема генератора последовательностей РЧ импульсов на основе ПЦС
Помимо этого, для каждого импульса МПУ формирует код фазы синтезируемого радиочастотного импульса. Если требуется формировать фазу из набора 0, π/2, π, 3π/2, то для кода сдвига фазы достаточно только 2 бита, которые будут складываться с учётом переноса с двумя старшими битами кода фазы φcode.
Реализация генератора на основе алгоритма прямого цифрового синтеза с использованием ПЛИС
Существуют современные примеры конструкций генераторов на основе ПЦС для ЯМР-релаксометров. Обычно они основаны на использовании специализированной микросхемы ПЦС, что требует присутствия внешнего быстродействующего управляющего устройства, например микроконтроллера или ПЛИС, для запуска и прекращения синтеза импульса [9, 10]. Так как управление запуском осуществляется командами через внешнюю шину, а микропроцессорное управляющее устройство может выполнять и другие подпрограммы, то не исключены ошибки формирования интервалов времени в последовательности импульсов. Другим вариантом является использование схемотехники на основе цифрового сигнального процессора (ЦСП, DSP) [11]. Однако при использовании ЦСП могут возникнуть указанные выше проблемы в том случае, если от него требуется параллельное выполнение нескольких операций по обработки потоков данных. В случае ЯМР-релаксометра, помимо синтеза передаваемого сигнала, требуется анализ принимаемого сигнала по определённым алгоритмам и поддержка внешних быстродействующих интерфейсов.
Возможность распараллеливания процессов обработки данных при условии их полной синхронизации даёт применение ПЛИС, тем более, что в рамках конструкции ЯМР-релаксометра на одной ПЛИС могут быть реализованы цифровые узлы генератора, приёмника и высокоскоростного интерфейса [12]. Недостатком конструкции в указанной работе является невозможность оперативной перенастройки параметров импульсной последовательности, так как она описывается VHDL кодом, загружаемым в процессе конфигурирования ПЛИС.
В рамках настоящей работы предлагается реализация генератора с ПЦС на ПЛИС Altera, лишённого данного недостатка. При этом для аппаратного описания каждого блока применён язык VerilogHDL[13], что позволяет легко перенести разработанное техническое решение на ПЛИС других производителей или заказать производство серии микросхем на базовом микрокристалле с требуемой конфигурацией.
Схема верхнего уровня, сформированная в среде автоматизированного проектирования структур на ПЛИС Altera Quartus II Lite, приведена на рис. 3.
Рис. 3. Схема верхнего уровня генератора с ПЦС, сформированная в среде автоматизированного проектирования структур на ПЛИС Altera Quartus II Lite
Частота опорного тактового генератора выбрана 160 МГц. Блок аккумулятора фазы (PH_ACC_A) выполнен с разрядностью 10 двоичных разрядов. Через шину TUNE [8:0] задаётся код приращения фазы в каждом такте. В соответствии с выражением (3) это позволяет производить настройку частоты с шагом 0,156 МГц.
Преобразователь фаза-амплитуда реализован в виде блока P_to_A_code, таблица перекодирования в его составе синтезирована на основе выражения (2). Указанный блок имеет два восьмиразрядных выхода: B0 и B. Выход B0 запускается после запуска аккумулятора фазы, как правило, при первом импульсе, и с него снимается опорный цифровой сигнал с нулевой фазой. Выход B запускается только при появлении сигнала Imp, его фаза определяется кодом на входе Ph [1:0]. Сигнал с выхода B подаётся на вход ЦАП через шину DATA_DAC, которая синхронизируется сигналом strob_DAC.
Счётчик условного времени (блок h_cnt) описан как классический двоичный счётчик с тактовой частотой fT = 160 МГц и количеством разрядов 40. Таким образом, квант условного времени tq может быть рассчитан по формуле (4) и составит 6,25*10-9с (6,25 нс), а максимальная продолжительность времени программируемой последовательности tc (рассчитывается по формуле (5)) составит 6871,94 с.
(4)
(5)
где Nh_cnt – разрядность счётчика.
В предложенном варианте МПУ (блок MPU4) позволяет производить программирование последовательности импульсов с возможностью выбора фазы несущей для каждого радиочастотного импульса из ряда 0, π/2, π, 3π/2. Количество импульсов в последовательности определяется объёмом памяти конкретной ПЛИС. Запуск и остановка со сбросом аккумулятора фазы осуществляется через линию S_P, а установка фазы через порт SetPh[1:0], запуск синтеза импульсов производится через линию Imp. МПУ имеет интерфейс программирования, реализованный в виде шин PrgA(адрес), PrgD(данные), WE(разрешение записи), который позволяет программировать последовательности импульсов, записывая их параметры в соответствующие регистры.
Перед созданием макетного образца с указанными параметрами экономически целесообразно произвести компьютерное моделирование выполнения генератором заданных функций. Для этих целей использована система моделирования ModelSim [14], совместимая со средой Quartus II Lite с возможностью визуализации результатов. Результаты моделирования представлены на рис. 4 и 5.
Рис. 4. Последовательность радиочастотных импульсов разной ширины, полученная при моделировании генератора в среде ModelSim
Рис. 5. Пример генерации последовательности радиочастотных импульсов (вверху) с частотой 20 МГц и фазой 0 (а), π/2 (б), π (в), 3π/2 (г) относительно опорной несущей (внизу)
На рис. 4 представлена последовательность из нескольких импульсов разной ширины. На рис. 5 представлены импульсы с разной фазой относительно опорного сигнала, а именно 0, π/2, π, 3π/2. Стоит отметить, что масштабы временных интервалов оптимизированы для визуализации возможностей предлагаемого генератора и могут несколько отличаться от применяемых в практических методиках ЯМР-спектроскопии.
Заключение
В настоящей работе предложена структурная схема генератора последовательностей радиочастотных импульсов для ЯМР-релаксометра, реализованная с использованием ПЛИС. Методами компьютерного моделирования показана возможность генерации радиочастотных импульсов различной длительности с достаточно точной выдержкой по времени и возможность задавать фазу несущей для каждого конкретного радиочастотного импульса, что требуется в многоимпульсных методиках измерения времён ЯМР-релаксации. Частота несущей также может быть изменена без переконфигурирования ПЛИС, то есть генератор может быть использован для приборов, регистрирующих времена релаксации ЯМР разного типа атомных ядер в образце.
Возможно применить предложенную структуру генератора последовательностей радиочастотных импульсов при проектировании цифровой части приёмника и передатчика ЯМР-релаксометра на одной микросхеме ПЛИС, что позволит миниатюризировать данный прибор.
Библиографическая ссылка
Лапшенков Е.М., Волков В.Я. МИНИАТЮРИЗАЦИЯ ГЕНЕРАТОРА ПОСЛЕДОВАТЕЛЬНОСТЕЙ РАДИОЧАСТОТНЫХ ИМПУЛЬСОВ ЯМР-РЕЛАКСОМЕТРА // Современные наукоемкие технологии. – 2017. – № 9. – С. 42-47;URL: https://top-technologies.ru/ru/article/view?id=36798 (дата обращения: 21.11.2024).