63878 (Физический анализ магнитно-резонансных томографов)

Физический анализ магнитно-резонансных томографов

1 Физические основы магнитно-резонансной томографии

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

m₀ = s, (1)

где - гиромагнитное отношение; для протона гиромагнитное отношение равно 242,6 МГц /Тл.

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

Рисунок 1. Действие на спины ядер магнитных полей.

Если на образец вещества подействовать переменным магнитным полем H₁(t), направленным поперечно к полю В₀, то спины ядер при определенных условиях могут отклониться от направления В₀ и прецессировать вокруг него, как волчки, если их слегка толкнуть. Для этого частота переменного магнитного поля и индукция постоянного магнитного поля должны быть связаны равенством Лармора.

₀ = В₀ . (2)

Частота ₀ называется частотой ядерного магнитного резонанса (она называется также ларморовой). Отсюда следует, что при магнитной индукции В₀ = 1 Тл частота ЯМР протона будет равна

= 42,6 МГц.

Если же индукция будет равна 0, 12 Тл, то частота ЯМР для протонов составит 5 МГц. Как видим, эти частоты лежат в диапазоне коротких радиоволн, которые считаются безвредными. И только в очень сильных магнитных полях (до З Тл) частота ЯМР может быть достаточно большой – 120 МГц. Но электромагнитные колебания и такой частоты еще не причиняют заметного вреда, особенно учитывая, как увидим далее, их малое время воздействия.

Импульс H₁(t) называют высокочастотным, или радиочастотным (РЧ), импульсом. По его окончании спины ядер возвращаются в исходное состояние. Этот процесс называется релаксацией. При этом ядра излучают электромагнитные колебания, которые могут быть зарегистрированы с помощью специальных антенн – РЧ катушек. Для этого могут использоваться те же катушки, которые излучают поле H₁(t) или другие. Величина эхо-сигнала определяется макроскопической намагниченностью объекта, которая складывается из магнитных моментов отдельных ядер. Вектор макроскопической намагниченности также прецессирует вокруг направления В₀. Следует отметить, что прецессия магнитных моментов имеет место и при отсутствии РЧ поля Н₁, так как всегда есть какие либо флуктуации магнитного поля. Однако эта прецессия незначительна.

Выберем систему координат, в которой ось z направим вдоль поля В₀, а ось x – перпендикулярно плоскости РЧ катушек (рис.2). Вектор М₀ будет иметь

Рисунок 2. Составляющие вектора результирующей намагниченности.

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

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

Для скорости релаксации поперечной намагниченности особое значение имеют фазы спинов отдельных ядер. Во время действия РЧ импульса спины всех ядер прецессируют синхронно и синфазно. После окончания импульса начинается убывание поперечной намагниченности, называемое спадом свободной индукции (ССИ). При этом из-за влияния соседних атомов и молекул и неоднородности постоянного магнитного поля В₀ происходит расфазировка спинов, что особенно заметно в сложных соединениях, образующих жиры и мышечные ткани. По этой причине, даже при одинаковом изменении поперечных проекций магнитных моментов отдельных ядер, макроскопическая поперечная намагниченность убывает значительно быстрее, чем растет продольная. Это поясняет рис.3, где показаны взаимные положения трех элементарных магнитных моментов m₀ через некоторые равные интервалы времени. В четвертом случае элементарные моменты расположены под углами 120^о, поэтому суммарная поперечная намагниченность равна нулю.

Рисунок 3. Поперечная и продольная релаксации.

Продольная и поперечная намагниченности M_z и M_xy изменяются приблизительно по экспоненциальному закону но с разными постоянными времени Т₁ и Т₂: Т₁ > T₂. На рис.3 показан случай, когда угол отклонения вектора начальной намагниченности М₀ от оси z после воздействия РЧ импульса составил 90^о.

В табл. 1 приведены времена Т₁ и Т₂ для некоторых типов тканей головного и спинного мозга; - процентное содержание воды в соответствующих тканях. Как видно из таблицы, различные ткани мозга заметно отличаются своими постоянными Т₁ и Т₂. Эти постоянные имеют большие абсолютные значения в ликворах и крови. В белом и сером веществе мозга они сильно отличаются (Т₁>> T₂). Величина Т₁ с ростом магнитной индукции В₀ увеличивается.

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

Таблица 1. Постоянные релаксаций Т₁ и Т₂.

2. Уравнение Блоха

Уравнение Блоха является основой для анализа электромагнитных процессов, возникающих при ЯМР. Оно получено из феноменологических представлений (не физических) и хорошо описывает поведение макросистемы в магнитном поле. Это уравнение имеет вид

. (3)

Член отражает незатухающую прецессию (ротацию), где произведение пропорционально , т.е. 1/t; векторная сумма - поперечная намагниченность; Т₁ и Т₂ - постоянные времени продольной и поперечной релаксаций. Форма второго и третьего слагаемых уравнения Блоха говорит о том, что процесс релаксации предполагается экспоненциальным. Это допущение справедливо для жидкостных сред (ликворов), однако является весьма приближенным для жиров, серого и белого вещества мозга и совсем далеко от истины для твердых образований, у которых Т₁ и Т₂ очень малы.

Положим, что Т₁ и Т₂ весьма велики. Тогда вторым и третьим членами в уравнении (3) можно пренебречь. Допустим также, что Н = Н₀ и Н = kН₀. Тогда уравнение Блоха примет вид

. (4)

Начальные значения составляющих М обозначим как . Представим М в виде , где i, j, k – орты, и выполним перемножение векторов согласно правилу, которое записано в виде таблицы.

Сравнивая левые и правые части в (4), находим

. (5)

Решим систему (4.5), положив . Знак "минус" здесь необходим для правильного отражения действия градиентных полей, в чем убедимся далее. Дифференцируя первое уравнение системы (5), с учетом второго получаем

или .

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

.

Полное решение системы (5) будет иметь вид

,

, (6)

.

При учете в уравнении Блоха членов, содержащих Т₁ и Т₂ первое и второе уравнения системы (6) следует умножить на exp(-t/T₂), а третье уравнение примет вид

.

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

Поперечную намагниченность можно представить в компактной комплексной форме

или , (7)