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

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

Иллюстрации

Показать все

Реферат

Настоящее изобретение, в общем, относится к области применения спинового ядерного магнитного резонанса, в частности к ядерной магнитно-резонансной спектроскопии (ЯМРС)/ ядерной магнитно-резонансной томографии (МРТ). Описан способ и устройство для генерации магнитно-резонансных данных, их получения и обработки.

С момента открытия спинового ядерного магнитного резонанса в плотном веществе Блохом [1] и Перселлом [1], независимо друг от друга, около 60 лет назад, он получил быстрое развитие и стал основным инструментом в исследованиях и инженерных разработках в области физики, химии, биологии, фармацевтики и пр. Особенно после новаторской работы Дамадьяна [3], Лотербура [3] и Менсфилда [4] в начале 1970-х развитие данного способа привнесло революционные изменения в технику диагностической визуализации в области медицины и здравоохранения.

В сущности, имеется две широкие области применения ядерного магнитного резонанса. Одна - ядерная магнитно-резонансная спектроскопия (спектрометр). Другая - ядерная магнитно-резонансная томография (сканер). В обеих областях требуется мощное магнитное поле В₀. В обеих используются одинаковые физические принципы, математические уравнения и большинство технологий получения и обработки данных, но их цели и результаты различны. Во избежание путаницы, в науке и промышленности используются следующие договоренности, в данном приложении используется сокращение ЯМРС для ядерной магнитно-резонансной спектроскопии (спектрометр), а термин МРТ используется для ядерной магнитно-резонансной томографии (сканер). ЯМРС часто используют в химических, физических и фармацевтических лабораториях для получения спиновых магнитно-резонансных частот, химических сдвигов и подробного спектра образца, тогда как МРТ используют в медицинских центрах и биологических лабораториях для получения 1D (одномерного), 2D (двумерного) и 3D (трехмерного) изображения ядерной спиновой плотности ρ, времени спин-решеточной релаксации T₁ и спин-спиновой релаксации Т₂ для тела человека или других образцов живых организмов.

Имеется две теоретические трактовки ядерного магнитного резонанса [5]. Одна основана на квантовой механике [5, 6] - всесторонняя и исчерпывающая, другая, основанная на полуклассическом электромагнетизме [5, 7] - феноменологическая. Оба эти описания взаимодополняющие. Описание квантовой механики можно количественно применить ко всем известным явлениям ядерного магнитного резонанса. Классические теории полезно применять для объяснения большинства экспериментов в области ядерного магнитного резонанса, за исключением некоторых нетипичных. Тем не менее, когда вопрос касается практического применения, классические теории превалируют. Классические уравнения Блоха в сочетании с импульсами поля высокой частоты (ВЧ) B₁, спиновым и градиентным эхом, пространственным кодированием и затуханием свободной индукции по большей части составляют сегодня так называемый импульсный магнитный резонанс. Современные приложения для ядерного магнитного резонанса фактически полностью теоретически основаны на классическом электромагнетизме [8].

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

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

на фиг. 1а изображен график уровней энергии спина - 1/2 в статическом однородном магнитном поле В₀;

на фиг. 1б изображены пространственные ориентации для спина - 1/2 в статическом однородном магнитном поле В₀;

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

на фиг. 3а и 3б схематически изображена приемная катушка согласно настоящему изобретению;

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

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

на фиг. 6 схематически изображены процессы получения, обработки и сокращения данных для МРТ согласно настоящему изобретению;

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

на фиг. 8а, 8б и 8в графически изображена коррекция функции корреляции согласно настоящему изобретению.

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

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

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

Согласно настоящему изобретению, случайный сигнал V_SR(t) в его оригинальном виде и без каких-либо манипуляций можно исследовать для получения обширных сведений о спектре S(U) спинового резонанса, спиновой плотности ρ и времени спиновой релаксации T₁ и Т₂.

Ключевыми в данном изобретении являются две функции - автокорреляция R(t) и спектральная функция S(v) для случайных спиновых магнитно-резонансных излучений. Из R(t) и S(v) можно получить другие параметры для ЯМРС и МРТ с помощью уравнений (3, 6, 7 и 8), как будет описано ниже. Теорема Винера-Хинчина, уравнение (5), устанавливает связь между R(t) и S(v). Переход от R(t) к S(v) - прямое преобразование Фурье, а переход от S(v) к R(t) - обратное преобразование Фурье. Следовательно, R(t) либо S(v) можно сначала получить из случайного сигнала спинового излучения V(t). Существуют различные способы, а также коммерческое программное обеспечение для расчета R(t) и S(v). В описании настоящего изобретения спиновую функцию автокорреляции R(t) сначала получают из исходного необработанного шумового сигнала спинового резонанса V(t), затем из уравнения (5) получают S(v). Следовательно, случайными неспиновыми шумами, не относящимися к сигналу V(t), можно на практике пренебречь в операции авто/кросс корреляции. Если сначала получают S(v) из V(t), а затем R(t) из S(v) с помощью обратного преобразования Фурье, S(v) обычно загрязняется неспиновыми случайными шумами и прочими нежелательными излучениями. Такие загрязняющие компоненты в S(v) должны быть на практике удалены перед тем, как S(v) используется при получении R(t) и в обработке данных сигнала.

Для получения функции авто/кросс корреляции требуются два входных сигнала. Для этой цели, в реализации настоящего изобретения применены два набора приемных катушек. На фиг. 3а и 3б, как будет подробно описано ниже, схематически изображены эти катушки, изображая две возможные конфигурации. Эти два набора приемных катушек расположены вместе, окружая образец, подвергаемый ЯМРС/МРТ, они генерируют два необработанных напряжения V_a(t) и V_b(t) на двух наборах выводов катушек. Фактически V_a(t) и V_b(t) содержат шум сигналов спинового резонансного излучения V_SRa(t) и V_SRb(t), другой случайный электронный шум V_Na(t) и V_Nb(t), а напряжения поля ВЧ B₁ - это V_B1a(t) и V_B1b(t). V_N1a(t) и V_N1b(t) будут сокращены при кросс-корреляции по причине их статистической независимости (если не принимать во внимание взаимную индукцию). V_sRa(t)=V_sRb(t); V_B1а(t)=V_B1b(t), V_B1a(t) и V_B1b(1) не являются случайными. Они не могут быть исключены в операции корреляции, но их вклад в загрязненную функцию автокорреляции R'(t) можно устранить при «коррекции R'(t) для R(t)», как показано на фиг. 4, 6 и 7а - 7в и будет более подробно обсуждено ниже.

Имеется два возможных варианта реализации изобретения, как показано на фиг. 4 и 6 - для двойного набора приемных катушек, а на фиг. 7а - 7в для одиночного набора приемных катушек. Одиночный набор приемных катушек в настоящем описании изобретения такой же, как и при стандартной импульсной ЯМРС/МРТ. Варианты реализации изобретения с одиночным набором катушек (фиг. 7а - 7в) требуют меньших затрат на оборудование, но устранение неспиновых сигнальных шумов с корреляцией невозможно (фиг. 7а) или возможно только частично (фиг. 7б и 7в). Удаление этих неспиновых сигнальных шумов выполняют в функциональном блоке «Коррекция R'(t) для R(t)"» (фиг. 7а - 7в). Поперечное магнитное поле высокой частоты B₁ в данном изобретении необходимо для измерения времени релаксации Т₂, но если требуется только время релаксации T₁, измерение также можно выполнить без поля B₁. В любом случае данное поле высокой частоты B₁ должно быть непрерывным и постоянным. Применение данного поля B₁ изменяет состояние динамического равновесия спинового ансамбля, поскольку способствует возрастанию вероятности резонансного перехода P_B1. Это эффективный способ увеличить силу спинового сигнала. В большинстве случае желательно крайне слабое магнитное поле.

Основной результат способа ЯМРС/МРТ - спектр мощности спинового магнитно-резонансного излучения S(υ), спиновая плотность ρ, время спиновой релаксации T₁ и Т₂, а также их отображение. Все эти параметры имеют истинные значения, а не так называемые «взвешенное» ρ, «взвешенное» T₁ и «взвешенное» Т₂.

Путем применения слабого высокочастотного поля B₁ для увеличения мощности сигналов спиновой эмиссии и применения авто/кросс корреляции для уменьшения неспиновых сигнальных шумовых эффектов, согласно данному изобретению, не только устранена необходимость сильного статического магнитного поля В₀, но также требуется меньше программного и аппаратного обеспечения. Следовательно, устройства для ЯМРС/МРТ с применением данного способа имеют меньшую стоимость, менее сложны, менее громоздки и более безопасны, более точны, имеют более высокое спектральное разрешение, проще в использовании, чем существующие импульсные ЯМРС/МРТ, таким образом, они предназначены для более широкого применения в науке, включая медицину и технику. Описание указанных импульсных ЯМРС/МРТ дается в терминах вектора намагниченности М исследуемого материала в статическом магнитном поле В₀, его прецессии относительно поля В₀ и сигналов затухания свободной индукции. Описание упомянутой теоретически основанной на квантовой теории непрерывной ЯМРС/МРТ согласно настоящему изобретению дано в терминах распределений спиновой заселенности, зеемановского расщепления спиновой энергии в магнитном поле В₀, вероятностей перехода между энергетическими уровнями Зеемана и случайного шума при спиновом магнитно-резонансном излучении.

Каждый атом имеет ядро. Каждое ядро состоит из протона (протонов) и нейтрона (нейтронов), за исключением ядра обыкновенного водорода ¹H, которое содержит один протон и ни одного нейтрона. Если ядро имеет один непарный протон или нейтрон, его механический спин относительно оси вызывает спиновой магнитный момент , где γ и h известны как спиновой гиромагнитный коэффициент и постоянная Планка соответственно. Символ I обозначает квантовое число спина. Наиболее часто используемое при МРТ ядро - протон ¹H; наиболее часто используемые при ЯМРС ядра - протон ¹H и углерод-13 ¹³С. Оба квантовых числа спина I для ¹H и ¹³С равны ¹/₂. Гиромагнитный коэффициент γ для ¹H и ¹³С соответственно равны 2,675×10⁸ рад с^-1 Тл^-1 и 6,729×10⁷ рад с^-1 Тл^-1.

Магнитные спины не имеют предпочтительной пространственной ориентации в среде при отсутствии магнитного поля. Если ансамбль ядерных магнитных спинов расположен в статическом магнитном поле В₀, направление которого обозначают осью z в декартовой системе координат, первоначальные вырожденные спиновые магнитные энергетические уровни разделяются на (2I+1) одинаково отстоящие друг от друга энергетические уровни Зеемана. Для случая ¹H или ¹³С I=¹/₂ и имеется только (2I+1)=2 уровня Зеемана. Высокий (верхний) энергетический уровень E_h и низкий (нижний) энергетический уровень E_h показаны на фиг. 1а. При этом имеет место сопутствующее разделение энергии и пространственное квантование спинов. Спины энергии E₁ выравнивают себя в положительном направлении В₀; спины энергии E_h выравнивают себя в отрицательном направлении В₀ (фиг. 1б). При достижении температурного равновесия спиновая плотность n_h на верхнем уровне и спиновая плотность n₁ на нижнем уровне удовлетворяют распределению Больцмана: , где h - постоянная Планка, a k_B и Т обозначают соответственно постоянную Больцмана и спиновую температуру. Хотя n_I>n_h всегда, n_I практически равно n_h, отношение (n_I-n_h)/(n_h+n_I) крайне мало, независимо от того, какую силу имеет поле В₀ в стандартных лабораторных условиях. Разность (n_I-n_h) зависит от мощности В₀. Чем больше В₀, тем больше (n_I-n_h), а чем больше (n_I-n_h), тем больше намагниченность М образца и выше частота Лармора. Это одна из главных причин, по которой в импульсной ЯМРС/МРТ стремятся использовать более высокое В₀. В настоящем изобретении не используют намагниченность М и разность (n_I-n_h), поэтому более сильное поле В₀ не означает обязательно лучшую производительность ЯМРС/МРТ.

Хотя при тепловом равновесии спиновые плотности на верхнем и на нижнем уровне n_h и n_I остаются устойчивыми пока поле В₀ и температура Т остаются неизменными, спины на верхнем энергетическом уровне постоянно, по причине спин-решеточных взаимодействий, переходят на нижний энергетический уровень с вероятностью перехода Р_B0 и наоборот, поэтому статистически dn_h/dt=dn_I/dt=0. Такое состояние называют «динамическим равновесием». Каждый спиновый переход с высокого на низкий уровень сопровождается эмиссией фотона с угловой скоростью (частотой Лармора) ω₀=γВ₀ (линейная частота v₀=ω₀/2π Герц). При каждой такой эмиссии фотона вырабатывается микроскопическое напряжение на выводах ближайшего устройства обнаружения (приемные катушки), расположенного вблизи исследуемого образца. Сложение всех таких эмиссий образует макроскопические сигналы спинового магнитного резонанса (излучения, обусловленного спиновым магнитным резонансом). Этот процесс является случайным. Полученные сигналы имеют форму стационарного случайного шума. Эти шумы в основном слабы, но их можно измерить с помощью современных электронных средств. Временная последовательность этих шумов образует на выводах приемных катушек постоянный стационарный случайный шумовой сигнал от упомянутого излучения, обусловленного спиновым магнитным резонансом.

Указанные выше случайные спиновые резонансные излучения в магнитном поле В₀ происходят естественным образом и непрерывно, независимо от того, присутствует ли поперечное ВЧ магнитное поле B₁. Поперечное поле ВЧ B₁ при частоте Лармора стимулирует дополнительные случайные эмиссии с вероятностью перехода между двумя одинаковыми спиновыми зеемановскими энергетическими уровнями, что приводит к возрастанию спинового резонансного шумового сигнала . Так как имеется статистическая независимость между излучениями, вызванными В₀ и B₁, при условии, что имеются В₀ и B₁, общая вероятность спинового резонансного перехода и объединенный шумовой сигнал спинового магнитного резонансного излучения (фиг. 2). При правильном анализе, эта временная последовательность случайного шумового сигнала спинового резонансного излучения V_SR(t) дает богатую и подробную информацию о вероятностях спинового перехода, времени релаксации, спектре резонанса и прочем.

Вероятность перехода Р определяет скорость перехода между двумя энергетическими уровнями (фиг. 2). Размерность Р-1/с, обратное ей время релаксации Т (с), то есть Т=1/Р. Обратная величина - время спин-решеточной релаксации T₁, величина, обратная - время релаксации Т_B1 (T_B1 не является временем спин-спиновой релаксации Т₂), величина, обратная P_SR - время релаксации T_SR. Так как P_SR=Р_B0+P_B1, то 1/T_SR=1/T₁+1/T_B1 (фиг. 2).

Если спиновая плотность на верхнем энергетическом уровне - n_h, a соответствующая вероятность перехода - Р, количество переходов в секунду с верхнего уровня на нижний равняется n_h×P. Тогда мощность резонансного излучения W(t) может быть записана как W(t)=hv₀n_hP, где h и v₀ представляют постоянную Планка и частоту Лармора. W(t) пропорциональна квадрату шумового напряжения V(t). To есть увеличение Р приводит к увеличению резонансного шумового сигнала V(t).

Пока поле, указанное В₀, и температура Т остаются постоянными, вероятность и время спин-решеточной релаксации остаются постоянными. При воздействии поля B₁, как описано ниже, вероятность изменяется пропорционально квадрату B₁. По сравнению с высокочастотным полем B₁, статическое магнитное поле В₀ играет весьма умеренную роль в уменьшении или увеличении вероятности спинового перехода .

Вероятность перехода между зеемановскими уровнями, вызванного полем В1 с частотой Лармора, обозначают P_B1. При условии, что поле B₁ значительно слабее поля В₀ (это всегда имеет место в случае ЯМРС и МРТ), его воздействие на гамильтониан энергии спинового ансамбля можно рассматривать как возмущение. Тогда применяют стандартную для квантовой механики теорию возмущений, которая приводит к выражению, связывающему время спин-спиновой релаксации Т₂ и поле B₁ с вероятностью перехода, вызванного полем В1, Р_B1 [9, 10].

где γ и I - соответственно спиновые гиромагнитный коэффициент и спиновое квантовое число. Символ m - обозначает спиновое магнитное квантовое число. Для всех ядер с I=¹/₂ (тогда m=¹/₂), таких как протон (¹H), углерод-13 (¹³С) или фосфор-31 (³¹P), приведенное выше уравнение преобразуют в [9],

Если приложены статическое поле В₀ и поперечное высокочастотное поле с частотой Лармора, по причине статистической независимости, общая вероятность спинового резонансного перехода P_SR=Р_В0+P_B1. Поэтому в терминах времени релаксации

где T₁ и Т₂ - соответственно спин-решеточное и спин-спиновое время релаксации.

Уравнение (3) является основным уравнением в данном изобретении, которое устанавливает соотношение между априори известным полем B₁, величиной T_SR и временем релаксации T₁ и T₂. Если присутствуют как В₀, так и B₁, полный спиновый магнитно-резонансный сигнал . Шумовой сигнал содержит информацию о времени релаксации T₁; шумовой сигнал содержит информацию о времени релаксации Т₂. Полный непрерывный стационарный случайный шум V_SR(t), время релаксации которого T_SR, содержит информацию как о T₁, так и о T₂.

Предположим стационарный во времени сигнал V(t), который может представлять , или или , его функцию автокорреляции R(t) определяют как:

R(t) - четная функция вещественной переменной, R(t)=R(-t). Тогда функция автокорреляции, спектральная функция S(υ) и время релаксации Т сигнала V(t) могут быть строго определены с помощью двух математических теорем.

Теорема Винера-Хинчина, в которой говорится, что спектральная функция S(υ) временной функции V(t) есть преобразование Фурье ее функции автокорреляции R(t) [11, 12].

где j=(-1)^1/2, υ - частота. При наличии спектральной функции S(υ), спиновую резонансную частоту υ₀ можно найти как [13]

Теорема Борна-Вольфа, в которой говорится, что квадрат времени релаксации Т функции времени V(t) есть нормализованный второй момент квадрата ее функции автокорреляции R(t) [13, 14].

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

В дополнение к времени спиновой релаксации T₁ и Т₂, спиновая плотность ρ также является фундаментальным параметром в приложениях ЯМРС/МРТ. Спиновую плотность ρ можно получить из вероятности спинового перехода Р и спиновой резонансной спектральной функции S(υ) при резонансной частоте υ₀, так как S(v₀)=с×ρ×Р, где с - коэффициент калибровки. Обычно требуется относительная спиновая плотность ρ, которая дает следующее выражение:

Относительная спиновая численная плотность

Здесь вероятность Р и Т может представлять либо P_B1 и T₁, когда существует только В₀, или P_SR и T_SR, если имеется B₁. Если требуется абсолютная спиновая плотность ρ, необходимо вычислить коэффициент калибровки.

Приведенные выше семь уравнений (1/2, 3, 4, 5, 6, 7 и 8) создают основу анализа и предварительной обработки данных в настоящем изобретении. Одним словом, после получения шумового сигнала V_SR(t) спинового магнитно-резонансного излучения, получают его функцию R(t) автокорреляции или спектр мощности S(v). Спектральная функция спинового резонанса S(υ), резонансная частота v₀, спиновая численная плотность ρ, а также время релаксации T_SR могут быть получены с помощью уравнений (3, 4, 5, 6, 7 и 8). Без поля B₁ (B₁=0), V_SR(t)=V_B0(t), P_SR=P_B0, T_SR=T_B0=T₁; с полем B₁, V_SR(t)=V_B0(1)+V_B1(t), тогда согласно уравнению (3) T_SR зависит от T₁ и Т₂. Два набора измерений двух различных значений B₁ (одно B₁ может быть равно 0) дают два T_SR и два связанных уравнения (3), которые можно решить одновременно для T₁ и Т₂.

Параметры ρ, T₁ и Т₂, рассчитанные здесь согласно приведенным выше уравнениям и процедурам, представляют их «истинные» значения. Они отличаются от так называемых «взвешенных» T₁, Т₂ или ρ в импульсной ЯМРС/МРТ.

В действительности, уравнение (7) - не единственная формула для подсчета времени релаксации Т сигнала V(t) из функции автокорреляции R(t). Имеются и другие математические формулы для выполнения той же задачи. Например, Гудман [15] определяет время релаксации Т с помощью следующей формулы:

Оба уравнения (7) и (9) можно использовать для определения времени релаксации Т. Уравнение (9) проще, но уравнение (7) дает лучшее понимание физики. В данном описании изобретения расчеты времени релаксации основаны на уравнении (7).

1. Поперечное высокочастотное магнитное поле B₁

Аналогично импульсной ЯМРС/МРТ, в настоящем изобретении используется поперечное (в плоскости x-y) магнитное поле высокой частоты B₁, генерируемое набором передающих катушек. Однако, в отличие от поля высокой частоты B₁ в импульсной ЯМРС/МРТ, в настоящем изобретении используется непрерывное (не импульсное) и крайне слабое поле B₁ для возбуждения спинового магнитно-резонансного излучения. Это широкополосное (гораздо шире относительно диапазона указанных линий спинового резонансного излучения) высокочастотное магнитное поле. По причине непрерывного режима работы некоторые из высокочастотных магнитных полей B₁ неизбежно перехватываются приемными катушками и генерируют добавочное напряжение U_B1(t) наряду с сигналом спинового магнитно-резонансного излучения V_SR(t). Затем U_B1 и V_SR вместе подаются с приемных катушек на любые электронные устройства. U_B1 является нежелательным загрязнением (помехой) для V_SR. Соответственно, его следует подавлять до минимально возможного уровня. Для практического устранения этого загрязнения можно применить три независимых способа: (1) механический, путем изготовления приемной катушки (катушек) специальной конструкции, установки и настройки. В одном варианте расположения, приемную катушку располагают перпендикулярно (90 градусов) к набору передающих катушек для устранения перекрестных помех между ними; (2) электронный, с помощью применения некоторой цепи компенсации; и (3) численный, путем применения корректирующих методов для окончательного устранения влияния U_B1. Подробное описание третьего способа приведено ниже.

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

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

В настоящем способе непрерывной ЯМРС/МРТ применяют два типа приемных (обнаруживающих) катушек. Один состоит из двух приемных катушек для генерации двух одинаковых шумовых сигналов спинового резонансного излучения V_SRa(t) и V_SRb(t). Другой тип состоит из одиночной приемной катушки, такой же, как и при импульсной ЯМРС/МРТ, используемой для генерации шумового сигнала спинового резонансного излучения V_SR(t). На фиг. 3а и 3б приведены эти две приемные катушки. Их можно установить так, чтобы они с двух сторон окружали исследуемый образец (фиг. 3а), или намотать совместно вокруг исследуемого образца (фиг. 3б). С двух выводов катушки 10 или пары катушек 10 сигнал V_a(t) подают на электронные устройства 14 и 62; с двух выводов катушки 12 или пары катушек 12 сигнал V_b(t) подают на электронные устройства 16 и 64 (см. фиг. 4 и фиг. 6). Необработанные напряжения V_a(t) и V_b(t) представляют собой сумму шумового сигнала спинового резонансного излучения V_SR(t), U_B1(t) и V_n(t). Здесь V_n(t) представляет все типы неспиновых случайных шумов, испускаемых катушками (а затем и последующими электронными устройствами V_a(t)=V_SRa(t)+U_B1(t)+V_na(t) и V_b(t)=V_SRb(t)+U_B1b(t)+V_nb(t). V_SRa(t)=V_SRb(t) и U_B1a(t)=U_B1b(t), но V_na(t)≠V_nb(t). Однако эти три типа сигналов V_SR(t), U_B1(t) и V_n(t) статистически взаимно независимы. Кроме того, V_na(t) статистически независим от V_nb(t), если не принимать во внимание взаимную индуктивность между этими двумя катушками.

3. Описание непрерывной точной ЯМРС и МРТ

Основные параметры в приложениях ЯМРС или МРТ - линейный профиль спинового магнитно-резонансного излучения (спектр мощности) S(v), спиновая плотность ρ, время спин-решеточной (продольной) релаксации T₁ и время спин-спиновой (поперечной) релаксации Т₂. Другие параметры, необходимые для некоторых специальных видов ЯМРС/МРТ, можно получить из данных измерений.

В общем случае, приведенные выше параметры являются функциями положения x, y и z в объеме образца, что требует, таким образом, 1-мерного, 2-мерного или 3-мерного сканирования. Образцы для ЯМРС-приложений обычно однородны по своему объему, и эти параметры предаются без изменений по всему объему.

3-1. Ядерная магнитно-резонансная спектроскопия (ЯМРС)

Задачи ЯМРС-приложений требуют, в основном, получения ρ, T₁, Т₂ и подробного спинового резонансного спектра высокого разрешения для однородного исследуемого образца. Если определяют пространственное распределение этих параметров, оно становится предметом магнитно-резонансной спектроскопии (МРС).

На фиг.4 приведена блок-схема процесса генерации, получения и обработки данных в ЯМРС-приложениях. Блок 20 включает в себя механизм для ЯМРС, такой как магниты для поля B₀, высокочастотную предающую катушку для поперечного (в плоскости x-y) поля B₁ (блок 22) и наборы приемных катушек. Исследуемый образец помещают в статическое однородное поле B₀, в котором нет градиентного поля, так как требуется сканирование или спиновая локализация. В блоке 20 на фиг.4 имеются две приемные катушки, как описано выше и показано на фиг.3а и 3б. После того как образец поместили в магнитное поле В₀, в образце естественным образом происходят спиновые магнитно-резонансные излучения, в результате чего на выводах каждой из двух катушек генерируются два необработанных сигнала V_a(t) и V_b(t). Для электронных устройств 14, V_a=V_SRa+V_na+U_B1a; для электронных устройств 16, V_b=V_SRb+V_nb+U_B1b. После раздельного прохождения через электронные устройства 14 и 16 (электронные устройства 14 и 16 идентичны), V_a(t) и V_b(t) встречаются в авто/кросс-корреляторе 24 для корреляции. Коррелятор 24 работает как автокоррелятор для V_sr и U_B1, возвращая необработанную функцию автокорреляции R′(t)=R_SR(t) от V_sr+R_B1(t) от U_B1, так как V_SRa=V_SRb и U_B1a=U_B1b. Для статистически независимых V_na(t) и V_nb(t) коррелятор 24 работает как кросс-коррелятор, возвращая функцию корреляции R_n(t)=0. Таким образом, необработанная R′(t)=R_SR(t)+R_a1(t)+R_n(t)=R_SR(t)+R_B1(t). Для ЯМРС приложений требуется только R_SR(t). R_B1(t) и остаточную R_N(1) необходимо исключить из R′(t). Эту задачу выполняет блок коррекции 26 на фиг.4. Исправление 26 на вход представляет собой зашумленную необработанную функцию корреляции R′(t), после коррекции на выходе получают спиновую функцию автокорреляции R(t)=R_SR(t). Это тот случай, когда присутствуют как статическое поле B₀, так и ВЧ поле B₁. Если ВЧ поле B₁ не применяют, U_B1(t)=0 и R_B1(t) отсутствует, следовательно, блок коррекции 26 возможно становится необязательным.

Если получили R(t), теорема Винера-Хинчина (5) и (6) дает спектр спинового резонанса S(v) и частоту спинового резонанса ν₀. Уравнение (7) дает точное значение для времени релаксации T из R(t), а уравнение (8) дает спиновую численную плотность ρ (блок 28).

Если требуются только S(v), ρ и Т₁, не требуется применять в описанной выше процедуре ВЧ поле B₁. Без B₁ из одного набора измерений можно получить S(v), ρ и T₁. С другой стороны, если требуется найти время релаксации Т₂, приведенная выше процедура может быть повторена дважды для генерации двух R(t) для двух различных значений B₁ (одно из двух значений B₁ можно взять равным 0). Два R(t) дают два T_SR согласно уравнению (7). Используя два этих значения R_SR, время релаксации T₁ и Т₂ можно получить путем совместного решения двух уравнений (3), одного для первого B₁, а второго - для второго B₁.

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

Когда применяется высокочастотное поле В₁, его необходимое напряжение может быть выбрано согласно гиромагнитному соотношению у спинов, предположительным временам T₁ и Т₂ релаксации спинов и уравнению (3). Для обычно встречающихся в ЯМР времен T₁ и Т₂ релаксации требуется очень малое поле B₁.

3-2. Ядерная магнитно-резонансная томография (МРТ)

Задача МРТ-приложений состоит в получении пространственного распределения, то есть 1-мерного, 2-мерного или 3-мерного изображения спиновой плотности ρ, спин-решеточного времени релаксации T₁ и спин-спинового времени релаксации Т₂ для образцов, таких как ткани и тело человека и тому подобное. Спиновая резонансная частота - обычно известный параметр. В таком случае необходимо иметь в наличии специальный прибор или устройство для спиновой пространственной локализации.

Как это принципиально показано на фиг. 5а для одномерного изображения, в изучаемом образце такое устройство магнитным образом выделяет похожий на стержень объем 50 в сектор 52 высотой z и толщиной Δz. Вдоль оси объема y устройство создает монотонно возрастающее (или убывающее), направленное по оси z частотно-кодируемое магнитное поле B_e(у) (54 на фиг.5б), так что каждому спиновому пространственному элементу с координатой у в данном узком объеме сопоставляют уникальное магнитное поле B_e=B_e(y), при этом соответствующая спиновая резонансная частота ν_e(y)=γ×[B₀+B_e(y)]/2π. Таким образом, все без исключения простра