Магнитно-резонансная томография тепловых эффектов в модельных средах тема автореферата и диссертации по физике, 01.04.03 ВАК РФ

Волков, Александр Александрович АВТОР
кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Москва МЕСТО ЗАЩИТЫ
2012 ГОД ЗАЩИТЫ
   
01.04.03 КОД ВАК РФ
Диссертация по физике на тему «Магнитно-резонансная томография тепловых эффектов в модельных средах»
 
Автореферат диссертации на тему "Магнитно-резонансная томография тепловых эффектов в модельных средах"

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

ВОЛКОВ Александр Александрович

МАГНИТНО-РЕЗОНАНСНАЯ ТОМОГРАФИЯ ТЕПЛОВЫХ ЭФФЕКТОВ В МОДЕЛЬНЫХ СРЕДАХ

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Л ,

Специальность 01.04.03 - радиофизика, 01.04.01 - приборы и методы экспериментальной физики

МОСКВА 2012

005045616

005045616

Работа выполнена в Институте общей физики имени A.M. Прохорова РАН и лаборатории магнитной томографии и спектроскопии факультета фундаментальной медицины МГУ имени М.В. Ломоносова

Научные руководители - доктор физико-математических наук,

главный научный сотрудник A.C. Прохоров

доктор физико-математических наук, профессор Ю.А. Пирогов

Официальные оппоненты - доктор физико-математических наук,

ведущий научный сотрудник A.B. Корженевский

кандидат физико-математических наук

Д.А. Куприянов

Ведущая организация - Физический институт

имени П.Н. Лебедева РАН

Защита состоится 10 мая 2012 года в 16 часов на заседании диссертационного совета Д 501.001.67 в Московском государственном университете имени М.В. Ломоносова по адресу: 119991 ГСП-1, Москва, Ленинские горы, МГУ, Физический факультет, аудитория имени Р.В. Хохлова.

С диссертацией можно ознакомиться в библиотеке физического факультета МГУ имени М.В. Ломоносова.

Автореферат разослан « 9 » апреля 2012 г.

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

кандидат физико-математических наук, доцент А.Ф. Королев

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

Актуальность темы.

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

В настоящей работе предпринята попытка на базе MP-томографов с магнитными полями 0.5 и 7 'Гл лаборатории магнитной томографии и спектроскопии факультета Фундаментальной медицины МТУ имени М.В. Ломоносова осуществить целенаправленное МРТ-термометрическое исследование физическими методами на модельных объектах и получить экспериментальный материал, доступный для количественного анализа.

Цели и задачи исследования:

1) создание аппаратуры и методик для проведения модельных МРТ-экспериментов в области термометрии;

2) оптимизация МРТ-методов аппаратных и компьютерных для наблюдения температурных полей;

3) разработка моделей описания наблюдаемых МРТ-методом тепловых полей;

4) наблюдение ЯМР-отклика водосодержащего раствора и оценка парциального вклада в ЯМР-сигнал молекул воды.

Защищаемые положения:

1) проведение физических модельных экспериментов является необходимой составляющей работ по продвижению МРТ-методов контроля температурных полей в гипертермию;

2) в модельных гипертермических экспериментах на МР-томографах возможно применение традиционных методов производства тепла с помощью металлосодержащих электронагревателей;

3) подбор параметров импульсных последовательностей и конфигурацией эксперимента позволяет осуществлять радикальное подавление и полное устранение артефактов;

4) наиболее удобными параметрами МРТ-контроля температурных полей на ядрах 'Н в области комнатных температур являются линейно зависящие от температуры время релаксации Г; и химический сдви;.

5) вклад в величину ЯМР-сигнала водосодержащих объектов вносят продукты автоионизации молекул Н20 - ионы Н30+ и ОН".

Новизна и практическая значимость результатов

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

Личный вклад автора

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

Апробация

Результаты диссертационных исследований докладывались на всероссийских и международных конференциях: 3-ем Евразийском конгрессе по медицинской физике и инженерии «Медицинская физика - 2010» (г. Москва, 21-25 июня 2010 г.); международных научных конференциях студентов, аспирантов и молодых учёных «Ломоносов-2010», «Ломоносов-2011» и «Ломоносов-2012» (г. Москва); 8-ой Всероссийской школе-семинаре «Физика и применение микроволн», (г. Звенигород, 23-28 мая 2011 г.); научно-практической конференции «Инновационный проект 2011» (г. Москва, 11 октября 2011 г.); и прошли апробацию на семинарах Института общей физики имени A.M. Прохорова в 2010-2011 г., на научном семинаре в томографическом центре университета штата Вашингтон (Сиэтл, 2008).

Публикации

Результаты диссертации отражены в 12 публикациях - 5 - в рецензируемых журналах, входящих в список ВАК, и 7 тезисах и трудах конференций.

Структура и объем работы

Диссертация состоит из введения, четырех глав, заключения и списка литературы. Работа изложена на 130 страницах, содержит 90 рисунков. Список цитируемой литературы включает 110 наименований

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

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

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

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

В первой части приводятся сведения о заложенных в основу МРТ принципах [1]. Описываются главные характеристики ЯМР - спин-решеточное и спин-спиновое времена релаксации, Тi и 7?, классические методы регистрации ЯМР, методы пространственно-селективного возбуждения спинов, методы построения МРТ-изображений. Обсуждаются сравнительные особенности стационарного и импульсного методов измерений на МР-томографах.

МРТ как метод исследования вещества - это практическая реализация идеи о том, что фундаментальный эффект ЯМР можно использовать для глубинного зондирования объектов и построения пространственных изображений их внутреннего строения. ЯМР-сигнал способны производить ядра с ненулевым спином, помещенные в постоянное магнитное поле б. В поле В спины ориентируются дискретным образом с интервалом энергий Е ~ В так, что при дополнительном наложении переменного электромагнитного поля с частотой со ~ В они резонансно поглощают/переизлучают электромагнитное излучение. Характеристикой ядер в отношении ЯМР является их магнитный момент //, типичные величины которого для наиболее используемых в медицине ядер составляют fj= 2.7-10"8 рад/(Тлс) для 'Н, и ¡л = 0.67 10 s рад/(Тл с) для |3С. В постоянном поле 1 Тл ЯМР на протонах реализуется на радиочастотах в области 40 МГц.

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

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

Важным качеством МРТ, отличающим его от других томографических методов и позволяющим без ограничений использовать МРТ для исследования живых тканей, является малость в условиях ЯМР энергетического возмущения АЕ исследуемой системы: А£~10'5 кТ в случае ЯМР на фоне АЕ ~10"3 кТ для ЭПР или ДЕ ~105 кТ для рентгеновского анализа.

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

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

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

МРТ-метод неинвазивен и в сравнении с контактными методами контроля тепла предполагает многократно расширенный охват исследуемого пространства. В той или иной степени, от температуры зависят все МРТ-регистрируемые параметры: времена Ti и 7"л намагниченность (протонная плотность), коэффициент диффузии и химический сдвиг. В последние годы опробовано множество вариантов использования этого факта в применении к МРТ-термометрии. По изменению коэффициента самодиффузии в отдельных случаях удается достигать термочувствительности 0.2 °С при объеме воксела 0.3 см3. Есть сообщения об измерении фазочувствительным методом температуры с шагом лучше 1°С с временным разрешением менее 10 секунд при объеме воксела 0.05 см2. Характерно, что подобные исследования проводятся в рамках биомедицинских работ, часто - на живых тканях и органах. О результатах судят по биологическому эффекту. Сообщения о физических модельных экспериментах на MP-томографах в области гипертермии в литературе отсутствуют.

Третья часть касается вопроса, который многие годы остается не решенным, хотя имеет для МРТ фундаментальное значение - о механизме ЯМР-активности воды, которая является одним из главных поставщиков протонов в биологических тканях. На ядрах 'Н проводится основная масса биомедицинских исследований (ПМР-спектроскопия), резонансная частота протона входит в названия фирменных ЯМР-установок: Bruker -200, 400, 1000. Принимается, что интегрально регистрируемый ЯМР-сигнал, или МРТ-контраст, исходит непосредственно от протонов и прямо отражает их количество в исследуемых объектах.

При исследовании водосодержащих сред считается естественным, что оба входящие в состав молекулы воды НгО протона 'Н ЯМР-активны. При этом игнорируется тот факт, что молекуле НгО присуще свойство спиновой изомерии: она способна с отношением 3:1 существовать в двух формах - в виде орто и пара молекул [3]. Спины 'Н в орто-молекуле расположены однонаправленно (складываются), а в пара-молекуле - разнонаправленно (взаимно компенсируются). Теоретически, орто-Н2О с суммарным спином 1 должна, как и отдельный протон со спином Vi, реагировать на внешнее магнитное поле и производить ЯМР-сигнал. Молекула пара-П20 с нулевым суммарным спином в ЯМР не активна. В приложении к МРТ последнее означает, что четверть находящейся в объекте воды остается невидимой. Так ли на самом деле - ответа нет.

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

Рис. 1. Магнитно-резонансные томографы, на которых проводились диссертационные

эксперименты: a) Bruker Tomikon S-50 с полем 0.5 Тл и б) Bruker BioSpec 70/30 с полем 7 Тл.

Вторая глава «МРТ-термометрия на модельных образцах с использованием традиционных методов нагрева» посвящена описанию, разработанных в диссертации МРТ-совмсстимых методов создания и регистрации температурных полей внугри модельных жидкостных образцов. С целью сравнения МР'Г-проявлений температурных полей в магнитных полях разной напряженности эксперименты проводились параллельно на двух представленных на рис. 1 MP-томографах: среднепольном медицинском томографоме Bruker Tomikon S50 с магнитным полем 0.5 Тл (частота магнитного резонанса 21, 2 МГц) и сильнопольном исследовательском томографоме Bruker BioSpec 70/30 с рабочим магнитным полем 7 Тл (частота 300 МГц). Проведение физических экспериментов на этих установках сопряжено с целым рядом специфических трудностей - ограниченностью рабочего пространства в сильном поле, ограниченностью выбора конструкционных материалов для любых аппаратурных приставок, невозможностью отключения магнитного поля на время монтажных работ, неприспособленностью узлов MP-томографов и рабочих протоколов для проведения термометрических измерений, дистанционностью измерений, невозможностью долговременной установки в МРТ нештатных устройств.

В первой части анализируются с точки зрения их пригодности для МРТ-термометрии термочувствительные МРТ-параметры, связанные с 1) изменением величины 7>релаксации (Г/-метод), сдвигом частоты резонанса протонов (Р-метод) и величиной коэффициента диффузии (/IDS-метод). Обсуждаются физические основы их

чувствительности к температуре, оценивается и сравнивается возможности создания на их основе измерительных МРТ-методик.

Медицинские температуры 30-45°С предельно высоки для их надежной МРТ-регистрации по той, прежде всего, причине, что в отношении ЯМР исследуемая система находится в режиме энергетического насыщения. При этом еще она имеет принципиально слабую связь с носителем температуры - атомно-молекулярным остовом ткани. Для достижения полезного для практического использования температурного разрешения (порядка десятых долей градуса) требуется регистрация температурных изменений, лежащих в пределах всего нескольких процентов абсолютной величины. На практике стационарные температурные картины, особенно при работе в клинических условиях, осложнены флуктуациями температуры в пределах нескольких градусов, обусловленных эффектами конвекции, перфузии, изменением уровня насыщения крови кислородом и т.п. Эти температурные искажения во время МРТ-сканирования ведут к систематической ошибке измерений. Для проведения корректной МРТ-термометрии требуются максимально короткие сроки сканирования. В диссертации с целью поиска путей сокращения времени сканирования и нивелирования тепловых флуктуационных эффектов осуществлена серия экспериментов по оптимизации параметров импульсных МР-последовательностей trueFISP, RARESE, PVM_PRESS для работы в режиме термометрии - рис. 2. Подбором параметров AV, RES и SEG достигался наиболее выгодный баланс времени и точности измерения контраста.

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

Рис. 2 7"/-взвешеиные изображения растительного масла, полученные при разных

наборах параметров импульсной последовательности true FISP (7 Тл). Обозначения на табличках: А V - количество накоплений сигнала, RES - пространственное разрешение, SEC -количество сегментов 7/. Контролируемые параметры: время релаксации Г/ и время измерения t.

Рис. 3 Модельный образец - блок из стеклянных трубок с навитой в центре 0.1 мм

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

Эксперименты проводились на модельных конструкциях, представляющих собой аксиально симметричные блоки стеклянных трубок с размещенными в центральной части проволочными нагревателями. В зависимости от материала нагревателя и МРТ-режимов прослеживалась пространственно-временная динамика артефактов, исходящих от материала нагревательного элемента и магнитного поля электрического тока (рис. 3). Проводилась оценка погрешностей измерений (искажений), вызванных присутствием в рабочей зоне проводящих сред - меди, стали, нихрома (рис. 4 и 5). Найдены способы снижения влияния артефактов путем варьирования ВАХ цепи нагревательного элемента, использования бинарной намотки провода, а также смещения зоны сканирования БОУ для подавления РЧ наводки от проводов, питающих нагревательный элемент и температурные датчики. Исследован вопрос о работоспособности температурных датчиков в условиях сильных магнитных полей, найдены условия их устойчивой работы. С учетом накопленных данных разработана МРТ-совместимая приставка, способная работать в условиях сильных магнитных полей и позволяющая создавать и регистрировать температурные поля внутри модельных жидкостных объемов (рис. 6).

Третья часть посвящена описанию термометрических МРТ-экспериментов, проведенных с помощью разработанной установки, по визуализации карт распределения и численной оценке времени ^/-релаксации и химического сдвига протонов на примере растительного масла. Осуществлена регистрация а) стационарного температурного градиента одновременно в пяти точках пространства с различньми температурами в рамках одного режима (ШдеИЗР и РУМ РИБЭЗ) и б) ступенчато изменяющейся температуры с одним температурным датчиком. Для количественного определения Т/ использована быстрая импульсная последовательность инверсии-восстановления 1гие_Р15Р (3 накопления, 0.1 см/пиксел, 13 шагов Т]: 4К-300) с автоматической перестройкой параметров последовательности в пределах одного режима (с возрастающими временами инверсии). При работе Р-методом использована последовательность РУМ_РЯЕЗЗ (32 накопления при зоне сканирования 10x10 мм). По совокупности измерений установлен линейный температурный рост обоих параметров, составивший 6.8 мс/'С (с точностью 3%) для 7>релаксации и 0.02 ррт/'С (с точностью 2.5%) для химического сдвига (рис. 7 и 8).

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

1.А

Рис. 4 Зависимость величины артефакта (количества искаженных вокселов) от силы

тока в проводе нагревательного элемента из 1) вольфрама, 2) меди и 3) нихрома. 7 Тл.

Рис. 5 Зависимость величины артефакта (количества искаженных вокселов) от силы

тока в 0.2 мм медной проволоке при 1) прямой и 2) встречной намотках. 7 Тл.

Рис. 6 Тгвзвешенные изображения деталей жидкостной кюветы-вставки (7 Тл).

Стрелкой отмечено расположение температурных датчиков. Контуром показана геометрия среза при сканировании в режимах П'ие РКР и РУК^РЯЕЭЗ.

fue. 7 Химический сдвиг линии резонанса протонов !Н в растительном масле при

изменении температуры. Импульсная последовательность PVM_PRESS. 7 Тл.

—.—I—•—I—■—i—•—i—■—i—-—i—-—i—>—i—•— 10 15 го ZS 30 35 40 <5 SO 55

T.C

Рис. 8 Температурные зависимости времени Г;-релаксации (последовательность

truc_FISP) и величины химического сдвига для растительного масла, полученные с помощью последовательностей и PVM_PRESS. 7 Тл.

В третьей главе «Регистрация температурного поля Г /-методом в условиях лазерной гипертермии» описываются эксперименты по МРТ-термометрии модельных объектов -растительного масла, сырого и термически коагулированного яичного белка, при генерации теплового поля путем световодной подачи внутрь исследуемых образцов лазерного излучения.

Температурное поле, создаваемое лазерным источником, в общем случае, по времени неустойчиво. При подаче световой энергии внутрь образца температура в зоне нагрева растет, а при выключении источника излучения падает. Для корректной регистрации МРТ-картины требуется хотя бы минутное поддержание стационарной картины температурного поля. Отработка нужного режима выполнялась в 2 этапа: калибровочные термические измерения вне магнитного поля и собственно МРТ-исследование (в присутствие поля). Цель первого этапа - регистрация и анализ пространственно-временного поведения профиля температурного поля и подбор параметров нагрева, оптимальных для использования в МРТ. Цель МРТ-измерепий визуализация зоны нагрева в реальном времени.

В первой части рассматривается пирометрический метод регистрации температурного поля. Обосновывается подбор веществ для исследования, описывается установка для пирометрических измерений: передвижной блок с закрепленным на нем пирометром с длиной фокусировки 10 мм (рис. 9). Блок движется прямолинейно над

размещенным под ним образцом (фантомом), в глубину которого до центра введен волоконный световод. Пироприемник регистрирует температуру поверхности по линии над световодом. В фантом вводится ИК излучение с длиной волны 1.07 мкм от лазерной установки «Модуль ГФ» с дискретными уровнями мощности 2, 4 и 6 Вт. Проводится линейное сканирование поверхностной зоны нагрева, являющейся проекцией на поверхность глубинного теплового процесса, вызванного лазерным нагревом. По результатам сканирования строятся пространственно-временные зависимости пиросигнала (поверхностной температуры), которые закладываются в основу дальнейшего моделирования объемного температурного поля внутри фантома.

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

аТ 1 Си ---«\7"Г :<•/, ... (1)

где а = х/(ср)> X ~ коэффициент температуропроводности, с - теплоемкость, р - плотность вещества, с/ - функция источника.

Без источника тепла и при теплопроводности ткани, не зависящей от координаты и

температуры, распределение температуры передается гауссианом вида

д

^ . (2)

где А — площадь под колоколообразной кривой (характеризует полное количество неравновесной тепловой энергии в системе). При остывании без подачи тепла гауссиан монотонно расплывается со временем с сохранением площади А и падением температуры вблизи максимума по известному закону Т ~ А / При наличии источника цА) задача

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

1 !

1

фантом

Рис. 9 Схема пирометрических измерений при лазерном нагреве.

МР-томография

Нелинейный выброс

Теоретическая модель диффузии тепла 6Т! д!

Рис, 10 Температурный выброс при включении процесса подачи тепла.

заданием процесса нагрева линейной зависимостью А = Q + кТ, оставляющей уравнение теплопроводности (1) линейным. Из решения уравнения и подгонки параметров гауссиана под данные калибровочных пирометрических измерений найдены временные промежутки относительной устойчивости температурного поля, т.е. режимы, удобные для сканирования МРТ-методом.

В третьей части излагаются результаты эксперимента по МРТ ткани в режиме релаксации Т/ при тепловом воздействии на ткань. Исследовались пространственно-временные изменения МРТ-контраста в зависимости от температуры. Были осуществлены три режима сканирования - метод спин-эхо (SE), метод инверсии-восстановления (IR) и быстрый метод true_FISP, и построены временные зависимости Т/ для исследуемых образцов. Из сравнения этих данных с калибровочными пирометрическими данными временные зависимости Т/ пересчитывались в температурные. Зарегистрирован выброс Т/ в начальный момент времени и дальнейший его медленный температурой рост (Рис. 10). Найденный эффект описан ненулевой добавкой Q в формуле подачи тепла. Это означает резкий рост температуры в центре зоны нагрева. Практическим выводом для гипертермии может быть вывод об опасности неконтролируемо быстрого включения тепла с целью нагрева ткани даже в режиме совсем малой мощности.

Четвертая глава посвящена поиску ответа на вопрос о происхождении сигнала ЯМР в жидкой воде. Предложен эксперимент, который дает возможность выявить, индивидуальный протон ли 'Н со спином V% или о/отго-молекула Н2О со спином 1 резонируют в магнитном поле.

В первой■ части обосновывается схема проведения эксперимента. Предложено с погрешностью не хуже нескольких процентов измерить интегральный ЯМР-сигнал S от образцов жидкой воды, содержащих разное количество п молекул HjO для того, чтобы определить, пропорционально чему изменяется S от образцов - числу протонов 2п или числу оршо-молекул 1.5л. Различить два случая предложено путем введения в образцы НгО дозированных количеств тяжелой воды D2O, опираясь на тот факт, что концентрации протонов и опто-молекул в растворах H2O-D2O разного состава ведут себя по-разному.

Благодаря обратимости реакции

H20+D20^2HD0 (3)

равновесный раствор H2O-D2O с необходимостью содержит смешанные молекулы HDO. При малых начальных концентрациях НгО практически все молекулы НгО разрушаются и идут на образование HDO. Однако при соотношении 50/50 вода НгО отдает на строительство HDO только половину своих молекул, сама оставаясь в растворе в заметном количестве - до четверти объема раствора. Добавление D2O в образцы Н2О количества протонов не изменяет и ведет лишь к перераспределению протонов между молекулами НгО и HDO. Количество протонов с коэффициентом 2 линейно отслеживает увеличение исходной концентрации НгО. В то же время, D2O разрушает молекулы Н2О, как орто, так и пара, и уменьшает их количества. Разрушение молекул Н2О молекулами D2O порождает два конкурирующих процесса - уменьшение числа ЯМР-активных орто-молекул, уменьшающее S, и появление из разрушенных орто и пара-молекул ЯМР-активных молекул HDO, увеличивающее S. Если считать, что ЯМР-сигнал образуется протонами вне зависимости от того в состав, каких молекул они входят, то можно ожидать, что в процессе увеличения концентрации Н2О в присутствие тяжелой воды D20 ЯМР-сигнал растет линейно. Если, в другом варианте полагать, что ЯМР-активны только ортпо-молекулы, и неактивны nape-молекулы, то сигнал S должен расти медленнее и, главное, нелинейно. Заключение о линейности или нелинейности роста S подверглось в диссертации экспериментальной проверки.

Во второй части описывается эксперимент по сравнению числа ЯМР-активных протонов раствора HDO 50% и с опорным веществом - дистиллированной водой. На 0.5 Тл MP-томографе выполнено прецизионное сравнительное измерение интегральной амплитуды ЯМР-сигнала от образцов чистой воды и 50/50 ее смеси с тяжелой водой с помощью двух методов - ЯМР спектроскопии и МРТ-сканирования по протонной плотности. В качестве исследуемых образцов использовались три пробы 50% HDO, различающиеся по концентрации на величину ошибки меньше 0.1% и две пробы дистиллированной воды. Тяжелая вода D2O с молярной массой 20.04 грамм/моль содержится в обычной воде с молярной массой 18.04 грамм/моль в пропорции 1: 6400. Приготовление раствора с равным количеством молекул Н2О проводилось разбавлением исходного 90%-пого раствора D2O.

Окончательный результат — отрицательный: с погрешностью не хуже 5%, проявлений спиновой изомерии не зарегистрировано. Это свидетельствует о том, что в жидкой воде существуют неспаренные спины 'Н. Предположительно, они образуются в результате автоионизации молекулы Н20 при распаде её на ионы Н30+ и ОН".

Выводы

1. Аппаратурно реализованы режимы внутреннего нагрева модельных образцов (яичного белка, растительного масла и воды) лазерным излучением с Х=1.07 мкм при одновременном пирометрическом контроле температуры поверхности. Сравнением Г;-взвешенных МРТ-изображений температурных полей с данными пирометрических измерений в рамках диффузионной модели распространения тепла найдена общая закономерность роста времени релаксации Г/ с характерным значением 6,5 мс/°С, совпадающим с литературными данными. Обнаружен нелинейный характер нагрева в центре зоны. Высказано предположение о возможности возникновения в таких условиях опасных для медицинского применения неконтролируемых температурных выбросов.

2. Из сравнительных измерений Г;-релаксометрических карт модельных образцов в магнитных полях 0.5 и 7 Тл установлено десятикратное увеличение быстродействия МРТ-методик в сильном поле для достижения одинаковых по качеству контрастов.

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

4. Выявлены пространственно-временные закономерности формирования артефактов МР-изображений, исходящих от металлических проводников в условиях протекания в них токов (0.1 - 2 А). Реализованы методики уменьшения артефактов за счет компенсации индуцируемого токами магнитного поля путем бинарной намотки проволоки. Предложен способ нагрева подачей тока синфазно со считывающим импульсом РЧ катушки МР-юм01рафа. Предложены и реализованы программные методы улучшения качества МР-изображений и увеличения точности измерений с помощью оптимизации импульсных последовательностей КАКВ ЗЕ, 1гие_1Т5Р и РУМ_РЯЕ88.

5. Впервые температурные зависимости химического сдвига исследованы в модельных экспериментах. В магнитном поле 7 Тл для растительного масла на фоне артефактов выделена линейная Температурная зависимость резонансной частоты протонов с разрешением 0.02 ррт/°С, близким к результатам идеальных МРТ-измерений.

6. Предложен и реализован эксперимент по установлению источника ЯМР-сишала в жидкой воде. Методом сравнительных измерений интегральной амплитуды ЯМР-сигнала от

чистой воды и 50/50 ее смеси с тяжелой водой в магнитном поле 0.5 Тл в томографическом и спектроскопическом режимах с погрешностью 5% проявлений спиновой изомерии не зарегистрировано. Высказано предположение о том, что источниками сигнала ЯМР в воде являются продукты самоионизации молекулы Н20 - ионы НзО+ и ОН', обладающие неспаренными протонами.

Публикации по теме диссертации

1. A.A. Волков. В.Н. Никифоров, Ю.А. Пирогов, A.B. Иванов, A.C. Прохоров, Регистрация температурного поля методом магнитно-резонансной томографии, Медицинская физика, 20И, №1 (49), стр. 75-81.

2. A.A. Волков. П.О. Капралов, В.И. Тихонов, A.A. Волков, О происхождении сигнала ЯМР в жидкой воде, Известия РАН, серия физическая, 2011, том 75, №12, стр. 1776-177.

3. A.A. Волков. С.К. Какагельдыев, A.C. Прохоров, Ю.А. Пирогов, Традиционные методы нагрева в применении к магнитно-резонансной термометрии, Радиоэлектроника: электронный журнал, 2012, URL: http://jre.cplire.ru/jre/janl2/12/text.html.

4. A.A. Волков. С.К. Какагельдыев, A.C. Прохоров, Ю.А. Пирогов, МРТ-методы регистрации температурных полей в применении к гипертермии. Медицинская физика, 2011,№1 (51).

5. A.A. Волков. Н.В. Анисимов, В.Н. Никифоров, Ю.А. Пирогов, A.C. Прохоров, Поиск ЯМР-сигнала от спиновых изомеров воды в смеси H2O-D2O, Биофизика, 2012, в печати.

6. A.A. Волков. Ю.А. Пирогов, В.Н. Никифоров, Магнитно-резонансная термометрия методом Т ¡-релаксометрии, 3-й Евразийский конгресс по медицинской физике и инженерии «Медицинская физика-2010», 21-25 июня, Москва, 2010 г. Сборник материалов.

7. A.A. Волков, Ю.А. Пирогов, С.К. Какагельдыев, Магнитно-резонансная термометрия, Международная научная конференция студентов, аспирантов и молодых учёных «Ломоносов-2011», 11-15 апреля, Москва, 2011 г. Сборник материалов.

8. A.A. Волков. В.Н. Никифоров, Ю.А. Пирогов, Магнитно-резонансная термометрия в модельных средах. Международная научная конференция студентов, аспирантов и молодых учёных «Ломоносов-2010», 12-15 апреля, Москва, 2010 г. Сборник материалов.

9. A.A. Волков. Ю.А. Пирогов, В.Н. Никифоров, Магнитно-резонансная термометрия модельных объектов, 8 Всероссийская школа-семинар «Физика и применение микроволн», 23-28 мая, Звенигород, 2011 г. Сборник тезисов.

10. A.A. Волков. C.K. Какагельдыев, A.C. Прохоров, B.H. Никифоров, Ю.А. Пирогов, Традиционные методы нагрева в применении к магнитно-резонансной гипертермии, Научно-практическая конференция «Инновационный проект 2011», 11 октября, Москва, 2011 г. Сборник материалов.

11. A.A. Волков. П.О. Капралов, С. К. Какагельдыев, Ю.А. Пирогов, Традиционные методы нагрева в применении к магнитно-резонансной гипертермии, 8-ая Зимняя молодежная школа-конференция «Магнитный резонанс и его приложения Spinus 2011», 27 ноября - 3 декабря 2011 г. Санкт-Петербург. Сборник материалов.

12. С.К. Какагельдыев, A.A. Волков. Ю.А. Пирогов, Методы регистрации температурных полей в применении к магнитно-резонансной термометрии, Международная научная конференция студентов, аспирантов и молодых учёных «Ломоносов-2012», 9-13 апреля, Москва, 2012 г. Сборник материалов.

Цитированная литература

1. R.A. Jones, К. Kvaemess, P.A. Rinck, Т.Е. Southon. Магнитный резонанс в медицине, ред. П.А. Ринк, Издание 3, 1993.

2. V. Rieke, K.B. Pauly. MR Thermometry, J. Magn. Reson. Imaging 27, 376 (2008).

3. N. Bloembergen, E.M. Purcell, R.Y. Pound, Relaxation Effects in Nuclear Magnetic Resonance Absorption, Phys. Rev. 73, N 7, 679 (1948).

Подписано в печать:

09.04.2012

Заказ №7144 Тираж - 100 экз. Печать трафаретная. Типография «11-й ФОРМАТ» ИНН 7726330900 115230, Москва, Варшавское ш., 36 (499) 788-78-56 www.autoreferat.ru

 
Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Волков, Александр Александрович, Москва

ИНСТИТУТ ОБЩЕЙ ФИЗИКИ имени A.M. ПРОХОРОВА РАН

ЛАБОРАТРИЯ МАГНИТНОЙ ТОМОГРАФИИ И СПЕКТРОСКОПИИ ФАКУЛЬТЕТА ФУНДАМЕНТАЛЬНОЙ

МЕДИЦИНЫ МГУ имени М.В. ЛОМОНОСОВА

61 12-1/792

ВОЛКОВ Александр Александрович

МАГНИТНО-РЕЗОНАНСНАЯ ТОМОГРАФИЯ ТЕПЛОВЫХ ЭФФЕКТОВ В МОДЕЛЬНЫХ СРЕДАХ

Диссертация на соискание ученой степени кандидата физико-математических наук

Специальность 01.04.03 радиофизика и 01.04.01 - приборы и методы экспериментальной физики

МОСКВА 2012

ОГЛАВЛЕНИЕ

Введение..........................................................................................................................4

Глава 1. Проблема применения МРТ в гипертермии. Обзор литературы..............7

Часть 1. Основные принципы магнитно-резонансной томографии........................7

1.1.1 МРТ-метод........................................................................................................7

1.1.2 ЯМР-сигнал......................................................................................................9

1.1.3. Спад свободной индукции...........................................................................12

1.1.4. ЯМР-релаксация...........................................................................................13

1.1.5. Импульсные последовательности...............................................................18

1.1.6. Пространственное детектирование избранных спинов............................25

1.1.7. Построение изображения.............................................................................26

Часть 2. МРТ как неинвазивный способ регистрации температурных полей.....27

1.2.1 Современное состояние термотерапии (гипертермии).............................28

1.2.2. Работа с температурными полями в термотерапии..................................30

1.2.3. Проблема контроля температуры в термотерапии...................................34

1.2.4. МРТ как инструмент мониторинга температурных полей......................36

1.2.5. МРТ-методы регистрации температурных полей.....................................38

Часть 3. Проблема регистрации спиновых изомеров воды...................................46

Глава 2. МРТ-термометрия модельных образцов с использованием традиционных методов нагрева...................................................................................47

Часть 1. Анализ термочувствительных параметров МРТ-метода на базе литературных данных................................................................................................48

2.1.1. Метод химического сдвига (Р-метод)........................................................48

2.1.2. Диффузионный метод (D-метод)................................................................52

2.1.3. Метод Т]-релаксации....................................................................................53

Часть 2. Описание экспериментальной установки.................................................55

2.2.1. Конфигурация МРТ аппаратуры.................................................................55

2.2.2. Описание импульсных последовательностей............................................57

2.2.3. Особенности МРТ-реализации фазового кодирования............................64

2.2.4. Определение Тг и Т2-релаксационных характеристик............................65

2.2.5. Установка для создания и регистрации квазистационарного температурного градиента.....................................................................................66

2.2.6. Установка для создания и регистрации стационарных температурных полей различной интенсивности...............................................72

2.2.7. Установка для исследования артефактов МР-изображений....................72

2.2.8. Технические требования к устройствам на совместимость с МРТ-аппаратурой.............................................................................................................75

2.2.9. Артефакты МРТ-изображений от металосодержащих материалов........78

Часть 3. МРТ-термометрия модельных образцов Тг и Р-методами....................83

2.3.1 Термометрия стационарных температурных полей различных интенсивностей.......................................................................................................87

2.3.2 МРТ-регистрация стационарного температурного градиента..................90

Глава 3. Регистрация температурного поля в условиях лазерной гипертермии... 99

3.1.1. Установка для создания и регистрации температурного поля на основе медицинского лазера и пироприемника..................................................99

3.1.2. Выбор образцов..........................................................................................100

Часть 2. Теоретическая модель диффузии тепла..................................................102

Часть 3. МРТ мониторинг лазерного нагрева.......................................................105

Глава 4. Поиск ЯМР-сигнала от спиновых изомеров воды в смеси Н20-020. ..110

Часть 1. Фазовая диаграмма раствора Н20Л320...................................................110

Часть 2. Описание эксперимента...........................................................................114

4.2.1 Приготовление образцов............................................................................114

4.2.2. Эксперимент...............................................................................................115

4.2.3. Результаты...................................................................................................116

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

Выводы........................................................................................................................ 121

Список литературы..................................................................................................122

Благодарности...........................................................................................................132

Введение

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

В настоящей работе предпринята попытка на базе MP-томографов с магнитными полями 0.5 и 7 Тл лаборатории магнитной томографии и спектроскопии факультета фундаментальной медицины МГУ имени М.В. Ломоносова осуществить целенаправленное МРТ-термометрическое исследование физическими методами на модельных объектах и получить экспериментальный материал, доступный для количественного анализа.

Цели и задачи исследования

- Создание аппаратуры и методик для проведения модельных МРТ-экспериментов в области термометрии.

- Оптимизация МРТ-методов аппаратных и компьютерных для наблюдения температурных полей.

- Разработка моделей описания наблюдаемых МРТ-методом тепловых полей.

- Наблюдение ЯМР-отклика водосодержащего раствора и оценка парциального вклада в ЯМР-сигнал молекул воды.

Защищаемые положения

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

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

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

- Наиболее удобными параметрами МРТ-контроля температурных полей на ядрах 'Н в области комнатных температур являются линейно зависящие от температуры время релаксации Т] и химический сдвиг.

- Вклад в величину ЯМР-сигнала водосодержащих объектов вносят продукты автоионизации молекул Н20 - ионы Н30+ и ОН".

Новизна и практическая значимость результатов.

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

Личный вклад автора

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

Апробация

Результаты диссертационных исследований докладывались на всероссийских и международных конференциях: 3-ем Евразийском конгрессе по медицинской физике и инженерии «Медицинская физика - 2010» (г. Москва, 21-25 июня 2010 г.); международных научных конференциях студентов, аспирантов и молодых учёных «Ломоносов-2010», «Ломоносов-2011» и «Ломоносов-2012» (г. Москва); 8-ой Всероссийской школе-семинаре «Физика и применение микроволн», (г. Звенигород, 23-28 мая 2011 г.); научно-практической конференции «Инновационный проект 2011» (г. Москва, 11 октября 2011 г.) и прошли апробацию на семинарах Института общей физики имени A.M. Прохорова в 2010-2011 г., на научном семинаре в томографическом центре университета штата Вашингтон (Сиэтл, 2008).

Публикации

Результаты диссертации отражены в 12 публикациях - 5 - в рецензируемых журналах, входящих в список ВАК, и 7 тезисах и трудах конференций.

Структура и объем работы

Диссертация состоит из введения, четырех глав, заключения и списка литературы. Работа изложена на 130 страницах, содержит 80 рисунков. Список цитируемой литературы включает 110 наименований.

Глава 1. Проблема применения МРТ в гипертермии. Обзор литературы.

Часть 1. Основные принципы магнитно-резонансной томографии. 1.1.1 МРТ-метод.

МРТ как метод исследования вещества - это практическая реализация идеи о том, что фундаментальный эффект ЯМР можно использовать для глубинного зондирования объектов и построения пространственных изображений их внутреннего строения [1-9]. ЯМР-сигнал способны производить ядра с ненулевым спином, помещенные в постоянное магнитное поле В. В поле В спины ориентируются дискретным образом с интервалом энергий Е ~ В так, что при дополнительном наложении переменного электромагнитного поля с частотой со ~ В они резонансно поглощают/переизлучают электромагнитное излучение. Характеристикой ядер в отношении ЯМР является их магнитный момент ц, типичные величины которого для наиболее используемых в медицине ядер составляют ц= 2.7-10"8 рад/(Тл с) для 'Н, и ц = 0.67-10"8 рад/(Тл с) для 13С. В постоянном поле 1 Тл ЯМР на протонах реализуется на радиочастотах в области 40 МГц.

Наблюдение ЯМР возможно для всех ядер второй и третьей групп таблицы Менделеева. Ядерный спин изотопов с четным числом протонов Z и

19 1А

четным числом нейтронов N равен нулю. На таких ядрах, типа С О, ЯМР принципиально не наблюдаем. Для наблюдения доступны группы ядер с четным Z и нечетным N или, наоборот, с нечетным Z при четном N.

Явление ЯМР имеет описание на двух языках - классическом и квантовом. При строгих последовательных рассмотрениях использует один из них, но при рассмотрениях общего характера - оба одновременно, без оговорок и переходов. С позиций квантовой механики угловой и магнитный моменты квантуются, т.е. имеют только дискретные значения. Величина проекции магнитного момента Zp. = yhml/27t. Для протонов она составляет ± уЬ/2л. В отсутствие внешнего магнитного поля спиновые состояния вырождены по энергии. При приложении к системе спинов внешнего магнитного поля В0 вырождение снимается, и ядра с направленными "по полю" и "против поля" магнитными моментами оказываются в состояниях с разной энергией. Для всей системы ядер возникает возможность энергетических переходов с одного уровня на другой. Работа двухуровневого перехода составляет физическую основу спектроскопии ядерного магнитного резонанса, и далее - МРТ. Следует специально отметить, что МРТ для медицинских целей естественным образом ориентирован на использование его при комнатной температуре (хотя сам источник постоянного магнитного поля находится при температуре жидкого гелия, 4.2 К ).

Важным качеством MPT, отличающим его от других томографических методов и позволяющим без ограничений использовать МРТ для исследования живых тканей, является малость в условиях ЯМР энергетического возмущения АЕ исследуемой системы: АЕ ~10"5 кТ в случае ЯМР на фоне АЕ -10° кТ для ЭПР или АЕ ~105 кТ для рентгеновского анализа. Метод ЯМР обладает высокой динамической чувствительностью, т.е. способен к регистрации нестационарных картин в широком спектральном

о

диапазоне 1 -10" с.

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

л л 1-1 1 ¿г

и О. Базовый изотоп кислорода в составе воды О в ЯМР не активен. Подавляющее число исследований проводится на ядрах 'Н (протонно-магнитная спектроскопия, ПМР).

МРТ-метод за последние 15 лет стал одним из ведущих методов неинвазивной диагностики [1]. Он является мощным инструментом изучения состояния биологической ткани. С помощью МРТ исследуют объекты по целому ряду параметров, получают трехмерные изображения внутренних структур и проводят оценку наиболее значимых характеристик. В мире в настоящее время работают около 10.000 магнитно-резонансных томографов. Впитав в себе весь опыт многолетней работы физиков, инженеров и конструкторов, МРТ-технология являет собой сгусток достижений современной науки и техники и дает ученым безграничные возможности изучения вещества. Устройство современных томографов настолько сложно, а методики многопараметричны, что работа на MP-томографах даже над стандартными задачами предполагает проведение целого комплекса предварительных (пристрелочных) исследований. Ученые во всем мире совершенствуют методики работы с томографами и открывают все новые технические возможности реализации своих задач. Несмотря на то что, в основу работы всех МРТ-установок заложено всего одно фундаментальное явление - ядерно-магнитный резонанс - число вариаций его использования в МРТ очень велико, а условия работы на разных установках сильно отличаться. Суммируя, можно сказать, что МРТ - очень молодой и сам по себе мало изученный экспериментальный метод исследования вещества. В

физической энциклопедии выпуска 1998 года, например, описание МРТ еще отсутствует.

1.1.2 ЯМР-сигнал.

В ЯМР томографии регистрация сигнала происходит от резонирующих

19 1"* 14. 1 q Л'?

ядер - водорода Н , Н , углерода С J, азота N , фтора F , натрия Na , фосфора Р31. Чаще всего, однако, используются протоны водорода Н1 по причине их высокого естественного содержания в биологических тканях [1-6]. Магнитный резонанс для ядер обусловлен наличием у этих частиц спиновых магнитных моментов. Последовательное описание магнитного резонанса возможно только на основе квантовой механики. Ряд важных понятий вытекает, однако, из электродинамики макроскопической системы зарядов.

Магнитный момент вводится для систем зарядов, движущихся стационарно в ограниченном объеме, как:

М = l/2cXnenrn х vn, (1.1)

где еп - заряд, гп - радиус-вектор заряда с номером n, vn - скорость.

Если все заряды и их массы одинаковы, то (1.1) можно переписать как:

M=^Inmrnxvn= yG, (1.2)

где G = Цпгп х pn - механический момент импульса системы, у = е/2шс -гиромагнитное отношение, константа пропорциональности между механическим и магнитным моментами.

В квантовомеханической интерпретации в трехмерной системе координат с направлением магнитного поля вдоль оси Z, взаимодействие магнитного момента с полем описывается гамильтонианом:

K = gNPNHIz. (1-3)

Гамильтониан имеет два собственных значения (два уровня энергии) и две собственные волновые функци