Акустояркостная термометрия биологических и биоподобных сред тема автореферата и диссертации по физике, 01.04.06 ВАК РФ

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

КРОТОВ Евгений Валерьевич

АКУСТОЯРКОСТНАЯ ТЕРМОМЕТРИЯ БИОЛОГИЧЕСКИХ И БИОПОДОБНЫХ СРЕД

01.04.06 - акустика

Автореферат

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

Нижний Новгород - 2003

Работа выполнена в Институте прикладной физики РАН (г. Нижний Новгород)

Научный руководитель: Официальные оппоненты:

Ведущая организация:

кандидат физико-математических наук А. Д. Мансфельд

доктор технических наук, профессор И. Я. Орлов

кандидат физико-математических наук Л. С. Долин

Институт радиотехники и электроники РАН (ИРЭ РАН)

Защита состоится « » ¿¿¿^/¿12003 г. в /& часов на заседании диссертационного совета Д 002.069.01 в Институте прикладной физики РАН (603950, г. Нижний Новгород, ул. Ульянова, 46).

С диссертацией можно ознакомиться в библиотеке Института прикладной физики РАН

Автореферат разослан « / ^ » ^ ¿Сс^^Р 2003 г.

Ученый секретарь диссертационного совета кандидат физико-математических наук ' ' А. И. Малеханов

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

Актуальность темы диссертации

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

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

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

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

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

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

РОС. НАЦИОНАЛЬНАЯ БИБЛИОТЕКА

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

Научная новизна

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

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

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

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

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

Практическая и научная ценность

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

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

Апробация работы. Публикации

Основные результаты диссертации изложены в 26 работах: в 13 статьях (из которых 2 опубликовано в отечественных и зарубежных рецензируемых научных журналах [9, 24], 1 статья в рецензируемом тематическом сборнике ИПФ РАН [19], 10 в тематических сборниках международных и всероссийских научных конференций и тематическом сборнике ИПФ РАН [7, 8, 13-16, 18, 21-23]), 13 докладах на международных, всероссийских и местных конференций [1—6, 10, 11,12, 17, 20, 25, 26].

Все основные результаты диссертации опубликованы.

Результаты диссертации докладывались на IV Всероссийской конференции по биомеханике (Н. Новгород, 1998), I Всероссийской научно-технической конференции "Компьютерные технологии в науке, проектировании и производстве" (Н. Новгород, 1999), II международной научно-технической конференции «Проблемы и прикладные вопросы физики» (Саранск, 1999), 4-5-й Нижегородской сессии молодых ученых (Н. Новгород, 1999-2000), 1-й международной конференции молодых ученых, Young Scientists on Laser Optics (LO-YS) (С.-Петербург, Россия, 2000), 2-6-й Научной конференции по радиофизике, (Н. Новгород, 1998-2002), Biomedical Optics (BIOS) (San Jose, California, 2001-2002), XI сессии Российского Акустического Общества (Москва, 2001), International School for Young Scientists on Optics, Laser Physics & Biophysics, (SFM' 01) (Saratov, Russia, 2001).

Структура и объем диссертации

Материал диссертации изложен в 4 главах и заключении. Диссертация содержит 164 страниц, 48 рисунков и 3 таблицы, список литературы содержит 132 наименования.

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

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

Первая глава посвящена исследованию способов картирования поля внутренней температуры путем многоканального приема акустояркостного

сигнала и последующей реконструкцией данных с помощью алгоритма алгебраической реконструктивной томографии (APT).

В разделе 1.1 обоснована актуальность применения многоканальных сканирующих систем и поставлена задача к главе.

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

Приемная антенна многоканального акустояркостного термотомографа состоит из 12 пьезоприемников, помещенных в герметичный антенный блок. Пьезоприемники с помощью двигателя и системы передачи могут синхронно поворачиваться в диапазоне углов ± 45° и тем самым осуществляется сканирование исследуемой среды. Блок заполнен согласующей жидкостью для обеспечения акустического контакта с исследуемой средой и термостатирования пьезоприемников. Сигнал, регистрируемый каждым пьезоприемником через звукопрозрачное окно антенного блока, поступает на 12 независимых усилительных каналов, где происходит их усиление и квадратичное детектирование. Далее сигналы поступают через многоканальное АЦП в персональный компьютер, где с помощью специально разработанного программного обеспечения происходит их регистрация, накопление, сохранение и обработка. Компенсация осуществляется также в персональном компьютере путем запоминания уровня сигнала от эталона черного тела и его вычитанием из каждого последующего измерения.

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

Экспериментально измеренная чувствительность многоканального акустояркостного термотомографа составила 0.3±0.1 К при 5 с накопления сигнала.

В разделе 1.3 исследована возможность применения алгоритма алгебраической реконструктивной томографии (APT), предназначенного для реконструкции данных, полученных многоканальным акустояркостным термотомографом с механическим угловым сканированием; проведено численное моделирование работы многоканального акустотермотомографа с учетом параметров антенны и среды, а также геометрии ультразвукового «пучка» приема; исследована устойчивость АРТ-алгоритма к зашумлению сигнала. Показано, что АРТ-алгоритм наилучшим образом подходит для решения поставленной задачи, когда массивы данных неполны, а сами дан-

ные сильно зашумлены. Для увеличения качества реконструированного изображения разработана и исследована методика пространственной низкочастотной Фурье-фильтрации.

В разделе 1.4 экспериментально исследованы возможности пространственной локализации тепловых источников в однородных и неоднородных средах с поглощением с помощью многоканального акустояркостного термотомографа. Показано, что исследуемым методом можно локализовать температурные неоднородности с точностью определения координаты 2-4 мм.

В разделе 1.5 экспериментально исследовано применение многоканального акустояркостного термотомографа для регистрации температурной реакции человеческого организма на различные человеческие нагрузки in vivo. Показана возможность применения метода для регистрации температурной динамики икроножной мышцы человека при температурном воздействии на симметричную мышцу другой ноги, а также физиологической реакции печени на глюкозный тест.

В разделе 1.6 проведено обсуждение исследований и результатов первой главы.

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

В разделе 2.1 обоснована актуальность применения фокусированных антенн и поставлена задача к главе.

В разделе 2.2 исследованы возможности применения фокусирующей антенны в компенсационном акустояркостном термотомографе. Также предложено и исследовано совместное применение фокусирующей антенны и корреляционной обработки сигнала для повышения пространственного разрешения. Проведен численный сравнительный анализ возможностей локализации малых объектов с помощью компенсационного акустояркостного термотомографа с одиночной фокусирующей антенной и корреляционного акустояркостного термотомографа с двумя фокусирующими антеннами. Показано, что применение корреляционной обработки сигнала в акустояркостном термотомографе с фокусирующей антенной увеличивает пространственную разрешающую способность приборов. Исследовано влияние соотношения между диаметром фокусной «перетяжки» и пространственным радиусом корреляции на чувствительность прибора. Получено условие применимости корреляционного акустояркостного термотомографа с фоку-

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

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

В разделе 2.4 экспериментально исследованы параметры созданного акустояркостного термотомографа. Экспериментально измеренная чувствительность прибора по черному телу составила величину порядка 0.2 ± 0.05 К при 5 с накопления сигнала. Проведена прямая локализация модельных тепловых источников путем отображения данных сканирования в амплитудном или яркостном представлении. Экспериментально показано, что для однородных сред без поглощения применение компенсационного акустояркостного термотомографа с фокусирующей антенной позволяет получить продольное пространственное разрешение порядка 10 ± 2 мм и поперечное пространственное разрешение порядка 3 ± 1 мм.

В разделе 2.5 проведено обсуждение исследований и результатов второй главы.

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

В разделе 3.1 обоснована актуальность исследований по контролю внутренней температуры при взаимодействии оптического излучения с биологическими и модельными средами и поставлена задача к главе.

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

В разделе 3.3 экспериментально исследована температурная динамика образцов биологических тканей (говяжьей печени и куриной грудки) при непрерывном режиме воздействия лазерного излучения. С помощью контактного электротермометра получена зависимость изменения локальной термодинамической температуры образца говяжьей печени от мощности воздействия излучения. Экспериментально показана возможность локализации источника греющего излучения по регистрации динамики акустояр-

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

В разделе 3.4 проведено обсуждение исследований и результатов третьей главы.

Четвертая глава посвящена исследованию методик измерения физических параметров мягких биологических тканей: температуропроводности, акустического и оптического поглощения.

В разделе 4.1 обоснована актуальность разработки методик измерения физических параметров мягких биологических тканей и поставлена задача к главе.

В разделе 4.2 предложена и исследована методика измерения температуропроводности мягких биологических тканей по динамике остывания прогретой области путем регистрации изменения акустояркостного сигнала многоканальным ЛЯТ после выключения воздействия лазерного излучения. Измеренная температуропроводность образца говяжьей печени составила величину порядка 1.3 • 10"7± 0.2 • 10"7 м2/с. Этот результат хорошо согласуется с известными справочными данными составляющими (1.3 -5-1.5) • 1СГ7 м2/с.

В разделе 4.3 предложена и исследована методика измерения акустического поглощения мягких биологических тканей по данным регистрации акустояркостного сигнала многоканальным акустояркостным термотомографом, полученным при различных глубинах расположения источника греющего излучения. Измеренный коэффициент акустического поглощения образца говяжьей печени составил величину порядка 0.6 ± 0.1 см"1. Независимо коэффициент акустического поглощения был измерен с помощью активной ультразвуковой локации путем регистрации последовательности отраженных от границ образца ультразвуковых сигналов и составил величину порядка 0.62 см-1. Результаты хорошо согласуются с данными предложенной методики. Также восстановлено изменение термодинамической температуры прогреваемой области по известным данным акустоякркост-ной температуры, температуропроводности и акустического поглощения. Йри расчетах была использована модель среды с однородным акустическим поглощением, в которой термодинамическая температура зависит только от одной координаты. Результат хорошо согласуется с независимыми измерениями проведенными с помощью контактного электротермометра.

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

акустояркостного сигнала нескольких последовательно расположенных каналов многоканальным акустояркостным термотомографом. Измеренный коэффициент оптического поглощения образца говяжьей печени составил величину порядка 3.3±0.3 см"1. Независимо коэффициент оптического поглощения образца был измерен с помощью импульсного оптико-акустического метода, основанного на регистрации акустической волны, порождаемой быстрым тепловым расширением неоднородности, поглотившей импульсное лазерное излучение. Коэффициент оптического поглощения был измерен путем сравнения оптико-акустического сигнала образца с оптико-акустическим сигналом эталона и составил величину порядка 3.1 см"1. Результаты хорошо согласуются с данными предложенной методики. Показано, что предложенная методика может быть применима для измерения коэффициентов поглощения греющего излучения сильно рассеивающих мутных сред.

В разделе 4.5 проведено обсуждение исследований и результатов третьей главы.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Экспериментально показана возможность неинвазивного картирования поля температурных неоднородностей однородных и неоднородных биоподобных сред с поглощением по их собственному тепловому акустическому излучению методом многоканальной акустояркостной (АЯ) термометрии.

2. Экспериментально продемонстрирована возможность регистрации реакции человеческого организма на различные физиологические нагрузки in vivo многоканальным акустояркостным термотомографом.

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

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

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

СПИСОК РАБОТ ПО ТЕМЕ ДИССЕРТАЦИИ

1. Кротов Е.В., Ксенофонтов С.Ю., Мансфельд А.Д., Рейман A.M., Санин А.Г. Многоканальная акустическая термография биологических объектов // Тез. докл. на IV Всерос. конф. по биомеханике, Н. Новгород, 1998. С. 64.

2. Кротов Е.В., Ксенофонтов С.Ю., Мансфельд А.Д., Рейман A.M. Разработка и экспериментальные исследования возможностей акустотермо-графа // Труды научной конференции по радиофизике, посвященной 80-летию Нижегородского государственного университета им. Н.И. Лобачевского. Н. Новгород, ННГУ, 1998. С. 59.

3. Кротов Е.В., Мансфельд А.Д., Ксенофонтов С.Ю., Крашенинников А.Ю., Ершов К.В. Применение алгоритма алгебраической реконструктивной томографии для задач многолучевой акустотермометрии плоскослоистых сред // Тез. докл. I Всерос. н.-т. конф. "Компьютерные технологии в науке, проектировании и производстве", г. Н. Новгород, НГТУ, 1999. Т. X. С. 21.

4. Крашенинников А.Ю., Мансфельд А.Д., Рейман A.M., Кротов Е.В. Применение компьютерных методов для обработки сигналов компенсационного акусторадиометра с фокусированной антенной // Тез. докл. I Все-росс. н.-т. конф. "Компьютерные технологии в науке, проектировании и производстве". Н. Новгород, НГТУ, 1999. Т. 3. С. 19.

5. Крашенинников А.Ю., Мансфельд А.Д., Кротов Е.В., Ершов К.В. Исследование флуктуации коэффициента усиления радиотракта компенсационного радиометра // Тез. докл. на I Всерос. н.-т. конф. "Компьютерные технологии в науке, проектировании и производстве. Н. Новгород, НГТУ, 1999. Т. XV. С. 38.

6. Кротов Е.В., Крашенинников А.Ю., Мансфельд А.Д. Особенности применения фокусированных антенн для задач акустотепловидения // 4-я Нижегородская сессия молодых ученых: Сб. тез. докл. Н. Новгород, 1999. С. 95.

7. Крашенинников А.Ю., Мансфельд А.Д., Кротов Е.В., Ершов К.В. Применение фокусированной антенны для акустической термографии // Учебный эксперимент в высшей школе: Научно-методический журнал. Саранск, 2000. Т. 1.С. 16-21.

8. Кротов Е.В., Мансфельд А.Д., Рейман A.M., Ксенофонтов С.Ю. Многоканальная акустическая термотомография плоскослоистых сред // Учебный эксперимент в высшей школе: Научно-методический журнал. Саранск, 2000. Т. 2. С. 13-17.

9. Кротов Е.В., Ксенофонтов С.Ю., Мансфельд А.Д., Рейман A.M., Санин А.Г., Прудников М.Б. Экспериментальные исследования возможностей многоканальной акустической термотомографии // Изв. ВУЗов. Радиофизика. 1999. Т. XLII, № 5. С. 479^*84.

10. Krotov E.V., Yakovlev I.V., Reyman A.M., Kirillov A.G. The tomography of biological objects at centimeter depths by optoacoustic and thermoacous-tic methods // Technical Digest of Conference LO-YS'2000 (First International Conference for Young Scientists on Laser Optics), Saint-Petersburg, RUSSIA,

2000. C. 149.

11. Кротов E.B., Мансфельд А.Д., Жадобов M.B. Характеристики и возможности многоканального акустотермографа // Труды 4-й научной конференции по радиофизике. Н. Новгород: ННГУ, 2000. С. 163.

12. Кротов Е.В., Мансфельд А.Д. Исследования возможностей применения акустотермографа корреляционного типа с фокусированной антенной // 5-я сессия молодых ученых. Н. Новгород, 2000. С. 102.

13. Кротов Е.В. Пространственная локализация нагретых источников с помощью корреляционного акустотермографа с фокусированной антенной // Избранные труды конкурса молодых ученых Института прикладной физики РАН. Н. Новгород: ИПФ РАН, 2001. С. 77-82.

14. Krotov E.V., Xenophontov S.Yu., Mansfeld A.D., Reyman A.M., Sanin A.G. Experimental study of the potential of multichannel acoustic thermotomo-graphy//Proc. SPIE. 2001. Vol. 4256. P. 101-108.

15. Reyman A.M., Krotov E.V., Mansfeld A.D. Application of the correlation method to the acoustic thermograph with focused antenna // Proc. SPIE.

2001. Vol.4256. P. 109-118.

16. Кротов E.B., Мансфельд А.Д., Рейман A.M., Вилков A.B. Двумерная акустическая термография биологических объектов // Сборник Тр. XI Сессии РАО. М.: АКИН, 2001. Т. 3. С. 165-169.

17. Кротов Е.В., Вилков В.А., Мансфельд А.Д. Локализация нагретых объектов с помощью корреляционного акустотермографа с фокусированной антенной // Тр. 5-й научной конференции по радиофизике / ННГУ им. Лобачевского, Н. Новгород: ТАЛАМ, 2001. С. 152.

18. Krotov E.V., Reyman A.M., Vilkov V.A., Mansfeld A.D. Internal temperature mapping of biological objects by an acoustical brightness thermometer with focused antenna // Proc. SPIE. 2002. Vol. 4618. P. 31-37.

19. Кротов Е.В., Вилков В.А., Мансфельд А.Д., Рейман A.M., Акусто-яркостный термотомограф с фокусированной антенной // Методы акустической диагностики неоднородных сред: Сб. научных трудов. Н. Новгород: ИПФ РАН, 2002. С. 235-239.

20. Кротов Е.В., Вилков В.А., Мансфельд А.Д., Жадобов М.В., Рейман A.M. Двумерная акустояркостная термотомография биологических объектов с помощью акустотермографа с фокусированной антенной // Тр. 6-й научной конференции по радиофизике. Н. Новгород: ННГУ, 2002. С. 145.

21. Krotov E.V., Yakovlev I.V., Reyman A.M., Zhadobov M.V. ABT-controllable laser hyperthermia of biological objects // Proc. SPIE. 2002. Vol. 4622. P. 100-105.

22. Krotov E.V., Yakovlev I.V., Reyman A.M., Zhadobov M.V., Zharov V.P. Application of acoustical thermometry to noninvasive monitoring of internal temperature during laser hyperthermia // Proc. SPIE. 2002. Vol. 4618. P. 38^5.

23. Krotov E.V., Vilkov V.A., Mansfeld A.D., Reyman A.M. Experimental investigations of heated sources localization by acoustic brightness thermograph with focused antenna // Proc. SPIE. 2002. Vol. 4707. P. 288-293.

24. Krotov E.V., Zhadobov M.V., Reyman A.M., Volkov G.P., Zharov V.P. Detection of thermal acoustic radiation from laser-heated deep tissue II Appl. Phys. Lett. 2002. Vol. 81, № 21. P. 3918-3920.

25. Жадобов M.B., Кротов E.B., Рейман A.M. Лазерная гипертермия мягких биологических тканей с контролем температуры внутренних слоев методом акустояркостной термометрии // Тр. 6-й научной конференции по радиофизике. Н. Новгород: ННГУ, 2002. С. 164.

26. Жадобов М.В., Кротов Е.В., Рейман A.M. Альтернативные методы измерения оптического поглощения мягких биологических тканей // Тр. 6-й научной конференции по радиофизике. Н. Новгород: ННГУ, 2002. С. 165.

ОГЛАВЛЕНИЕ ДИССЕРТАЦИИ

ВВЕДЕНИЕ 4 Глава 1. МНОГОКАНАЛЬНАЯ АКУСТОЯРКОСТНАЯ ТЕРМОМЕТРИЯ БИОЛОГИЧЕСКИХ ОБЪЕКТОВ 24

1.1. Введение 24

1.2.Многоканальтный акустояркостный термотомограф с механическим угловым сканированием 28

1.3. АРТ алгоритм для восстановления профиля температур биологических объектов 39

1.4. Картирование поля внутренних температур биологических сред многоканальным акустояркостным термотомографом 55

1.5. Измерения акустояркостной температуры от vivo 65

1.6. Результаты и выводы 71 Глава 2. ПРИМЕНЕНИЕ ФОКУСИРУЮЩИХ АНТЕНН ДЛЯ ЗАДАЧ АКУСТОЯРКОСТНОЙ ТЕРМОМЕТРИИ 74

2.1. Введение 74

2.2. Локализация нагретых источников с помощью акустояркостного термотомографа с фокусирующей антенной работающем в компенсационном и корреляционном режиме приема 76

2.3. Основные узлы акустояркостного термотомографа с фокусирующей антенной 92

2.4. Локализация малых нагретых объектов акустояркостным термотомографом с фокусирующей антенной 104

2.5. Результаты и выводы 116 Глава 3. КОНТРОЛЬ ВНУТРЕННЕЙ ТЕМПЕРАТУРЫ МЯГКИХ БИОЛОГИЧЕСКИХ ТКАНЕЙ ПРИ ВОЗДЕЙСТВИИ ОПТИЧЕСКОГО ИЗЛУЧЕНИЯ 118

3.1. Введение 118

3.2. Импульсно периодический режим воздействия лазерного излучения 120

3.3. Контроль температуры при воздействии непрерывного лазерного излучения 127

3.4. Результаты и выводы 139 Глава 4. ИЗМЕРЕНИЕ ФИЗИЧЕСКИХ ПАРАМЕТРОВ МЯГКИХ БИОЛОГИЧЕСКИХ ТКАНЕЙ 140

4.1. Введение 140

4.2. Измерение температуропроводности мягких биологических

тканей 141

4.3. Измерение акустического поглощения мягких биологических тканей 143

4.4. Измерение оптического поглощения мягких биологических тканей 148

4.5. Результаты и выводы 152 ЗАКЛЮЧЕНИЕ 153 ЛИТЕРАТУРА 155

КРОТОВ Евгений Валерьевич

АКУСТОЯРКОСТНАЯ ТЕРМОМЕТРИЯ БИОЛОГИЧЕСКИХ И БИОПОДОБНЫХ СРЕД

Автореферат

Формат 60 х 90 '/и. Бумага офсетная № 1. Усл. печ. л. 1,0. Тираж 100 экз. Заказ № 51(2003).

Отпечатано в типографии Института прикладной физики РАН, 603950, г. Н. Новгород, ул. Ульянова, 46

7&Г7

7 89 7

ВВЕДЕНИЕ

Глава 1. МНОГОКАНАЛЬНАЯ АКУСТОЯРКОСТНАЯ

ТЕРМОМЕТРИЯ БИОЛОГИЧЕСКИХ ОБЪЕКТОВ

1.1. Введение

1.2.Многоканальтный акустояркостный термотомограф с механическим угловым сканированием

1.3. APT алгоритм для восстановления профиля температур биологических объектов

1.4. Картирование поля внутренних температур биологических сред многоканальным акустояркостным термотомографом

1.5. Измерения акустояркостной температуры in vivo 65 1.6 Результаты и выводы

Глава 2. ПРИМЕНЕНИЕ ФОКУСИРУЮЩИХ АНТЕНН ДЛЯ

ЗАДАЧ АКУСТОЯРКОСТ1ЮЙ ТЕРМОМЕТРИИ

2.1. Введение •

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

2.3. Основные узлы акустояркостного термотомографа с фокусирующей антенной

2.4. Локализация малых нагретых объектов акустояркостным термотомографом с фокусирующей антенной

2.5. Результаты и выводы

Глава 3. КОНТРОЛЬ ВНУТРЕННЕЙ ТЕМПЕРАТУРЫ МЯГКИХ БИОЛОГИЧЕСКИХ ТКАНЕЙ ПРИ ВОЗДЕЙСТВИИ ОПТИЧЕСКОГО ИЗЛУЧЕНИЯ

3.1. Введение

3.2. Импульсно периодический режим воздействия лазерного излучения

3.3. Контроль температуры при воздействии непрерывного лазерного излучения

3.4. Результаты и выводы

Глава 4. ИЗМЕРЕНИЕ ФИЗИЧЕСКИХ ПАРАМЕТРОВ МЯГКИХ БИОЛОГИЧЕСКИХ TKAI1ЕЙ

4.1. Введение

4.2. Измерение температуропроводности мягких биологических тканей

4.3. Измерение акустического поглощения мягких биологических тканей

4.4. Измерение оптического поглощения мягких биологических тканей

4.5. Результаты и выводы

Развитие современной медицины предполагает разработку более совершенных методов оценки состояния человеческого организма. Важную информацию о патологиях внутри биологических тканей несет в себе пространственное распределение температуры. Исследование распределения внутренней температуры биологических объектов имеет важное значение, не только для создания методов диагностики различных заболеваний, но и для понимания физиологических процессов, происходящих в организме, в том числе при воздействии на него различных факторов. Информация о температуре исследуемого органа облегчает раннюю диагностику при различных патологиях. Например, исследование температурного распределения при глюкозном тесте может дать информацию о наличии локальных патологических изменений организма [1.1].

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

Известные неинвазивные методы, такие, как пассивная микроволновая радиометрия [1.3], ядерный магнитный резонанс [1.4, 1.5] и температурный мониторинг с помощью активных ультразвуковых методов [1.6] обладают рядом недостатков. Ограничением первого метода является слабое пространственное разрешение. Второй и третий методы требуют калибровки по температуре, которую необходимо производить на каждом конкретном пациенте.

Одним из наиболее перспективных методов исследования температурных полей внутри биологических объектов является акустояркостная (АЯ) Термометрия, которая в последние годы оформилась как исследовательское направление и развивается сравнительно быстрыми темпами. АЯ Термометрия основана на регистрации акустического излучения, порождаемого тепловым движением атомов и молекул среды[1.7-1.9].

Начало исследований возможностей регистрации теплового акустического излучения относится к середине 20 - го века. Наличие шумового акустического излучения воды было теоретически описано Р.Х. Мелленом [1.7] в 1952 г. В его работе рассмотрена возможность измерения температуры океана с помощью регистрации точечным гидрофоном акустического шумового излучения. В своей работе Медлен описал шумовое акустическое излучение идеальной жидкости с точки зрения теории классической статистической механики.

Экспериментально спектр шумового акустического излучения был измерен впервые Д.Х.Эзроу [1.8] в 1962 г. В своей работе Эзроу основывался на процедуре измерения сопротивления излучения от тонкого диска нагруженного на воду.

Позднее в 1994 г. Ю.К. Барабаненковым и В.И.Пасечником [1.9] была предложена волновая статистическая теория теплового акустического излучения, построенная в рамках модели гидродинамических флуктуаций. Было выведено волновое стохастическое уравнение для флуктуаций давления со сторонними случайными механическими силами, рассчитаны с учетом аппаратных функций регистрирующего устройства тепловые флуктуации давления в безграничной среде и в акустическом резонаторе.

Дальнейшее развитие метода АЯ Термометрии заключалось в исследовании законов переноса акустического излучения в среде [1.10-1.14].

Метод определения глубинной температуры океана и других жидких и твердых сред по измерениям интенсивности их собственного теплового акустического шумового излучения на поверхности сред был впервые предложен В.И.Бабием [1.10] в 1974 г. В этой работе Бабий рассмотрел уравнение теплопереноса. Им была рассмотрена одномерная задача, когда плоская акустическая волна распространяется в плоскослоистой среде, т.е. в случае, когда распределение внутренней температуры и коэффициента поглощения являются функциями только одной координаты, направленной вглубь объекта. Изменение акустической энергии при прохождении волны сквозь элементарный объем среды обусловлено ее поглощением и излучением в соответствии с законами Кирхгофа. В своей работе Бабий показал, что интенсивность теплового акустического излучения на поверхности объекта можно определить через яркостный эквивалент температуры - акустояркостную температуру.

АЯ Температура — это интегральная характеристика равная термодинамической температуре абсолютно черного тела, создающего такой же поток акустического излучения, что и исследуемый объект [1.15].

В 1981 г. Т. Боуэн получил патент США на систему для пассивного дистанционного измерения температуры по акустическому тепловому излучению [1.11]. Эту систему он предложил использовать для контроля температуры в однородных мягких тканях биологических объектов по аналогии с уже известными системами, регистрирующими электромагнитное излучение [1.12, 1.13]. Также им была дана оценка пороговой чувствительности акустояркостного термометра (АЯТ).

Работы Боуэна [1.14] позволили поставить ряд задач медицинской диагностики решаемых в рамках метода АЯ Термометрии.

В 1985 г. появилась работа Ю.В.Гуляева, Э.Э.Годика и др. о возможностях акустотермографии при измерениях внутренней температуры биологических объектов по сравнению с СВЧ-радиометрией [1.16]. Целью работы было получение количественных оценок чувствительности, глубинности и пространственного разрешения акустотермографии. Единичный акустический приемник проигрывает в чувствительности радиометру приблизительно в 30 раз (из-за разницы на три порядка в полосе принимаемых частот), однако чувствительность AT можно повысить за счет использования матричного приемника. Чувствительность акустояркостного приемника увеличивается пропорционально y/~N , где N количество элементов матрицы, при условии, что все элементы матрицы направлены на одну и ту же 8 окрестность исследуемого пространства, а интенсивность теплового акустического излучения измеряемая каждым элементом матрицы суммируется. В этой работе было также показано, что использование метода акустотермометрии в ряде биомедицинских приложений может оказаться выгоднее, чем применение СВЧ — радиометрии, из-за более высокой проникающей и пространственной разрешающей способности акустотермометрии [1.17].

Позднее в ряде российских работ были исследованы различные схемы создания акустояркостных термометров, оценена пороговая чувствительность, от которой в основном зависит погрешность определения акустояркостной температуры [1.15, 1.19]. Также были исследовано влияние различных потерь [1.12], была дана аналогично

1.20], общая формула для пороговой чувствительности [1.15]. Уточнение формулы с учетом свойств пьезопреобразователя было проведено в работе [1.21].

Экспериментальная проверка физических принципов акустояр костной термометрии была проведена В.И.Пасечником [1.15,1.22], который использовал модуляционный прием теплового акустического излучения. В этих работах измерения проводились на модельных объектах: касторовом масле, глицерине, для которых характерна сильная температурная зависимость коэффициента поглощения. Пасечник рассматривал зависимость акустояркостной температуры образца от его поглощения, температуры и его толщины. В этих работах Пасечник ввел поправку в формулу Бабия [1.10], которая учитывает отражение от слоя при распространении акустической волны. Для экспериментального определения акустояркостной температуры независимо определялось поглощение в образцах. При этом экспериментальные значения акустояркостной температуры хорошо совпали с теоретическими.

Сообщения о создании первых экспериментальных образцов акустояркостных термометров появились в 1987 г. В работах [1.15, 1.19] использовался модуляционный прием сигнала. Устройство такого акустояркостного термометра подробно описано в работах [1.15,1. 23,1.24].

В работе [1.19] использовался акустоэлектронный модулятор: кристалл фотопроводящего сульфида кадмия, одним из свойств которого является изменение в затухании акустических волн кристалла при подаче на него светового излучения. Это свойство использовалось для подавления излучения исследуемой среды. В этой работе было проведено измерение акустояркостной температуры предварительно охлажденного пластилинового диска, находящегося в аквариуме с водой при комнатной температуре. Температура диска контролировалась двумя терморезисторами, вставленными в диск на глубине 5 мм и на глубине около 0,3 мм. Измерения проводились в течение 10 мин, за которые диск постепенно нагревался. Значения измеренной акустояркостной температуры находились между значениями измеряемых термодинамических температур диска на разных глубинах образца. Эти измерения показали, что регистрируемая акустояркостная температура определяется как некоторая усредненная температура исследуемого объекта.

В 1992 г. Пасечником был получен патент [1.25] на акустояркостный термометр с немеханическим модулятором. В таком AT параллельно вторичной обмотке согласующего трансформатора через ключ было подсоединено активное сопротивление, термодинамическая температура которого совпадала с температурой пьезопреобразователя. При размыкании ключа, АЯТ принимал излучение от исследуемого объекта. При замкнутом ключе, модуль коэффициента электромеханического преобразования, снижается практически до нуля. Измеряемый в этом случае опорный электрический сигнал определяется температурой и емкостью пьезопреобразователя. На выходе такого АЯТ формируется напряжение, пропорциональное разности абсолютных температур тела и пьезопреобразователя.

При разработке первых акустотермометров в качестве приемных антенн использовались различные пьезопреобразователи. В работах В.И.Миргородского с соавторами [1.19] измерения проводились с помощью пьезопреобразователя из ниобата лития, акустически согласованного с исследуемой средой (водой) с помощью двух четвертьволновых слоев. Центральная частота пьезопреобразователя составляла 1,3 МГц, полоса пропускания около 0,7 МГц. В работе Т.Боуэна [1.12] использовался стандартный пьезопреобразователь с центральной частотой 2,25 МГц и полосой пропускания 1 МГц. В работах В.И.Пасечника регистраторами теплового акустического излучения служили: приемная антенна из пьезокерамики ЦТС-19-2 без согласующих слоев с центральной частотой 2 МГц и полосой пропускания - 0,7 МГц [1.22,1.24, 1.26] и приемная антенна из пьезокерамики ЦТС-19 с двумя четвертьволновыми согласующими слоями (f= 2 МГц, Д/= 1 МГц) [1.27], эта приемная антенна была изготовлена в ИПФ РАН под руководством А.Д. Мансфельда. В этих же работах было дано обоснование регистрации акустического излучения от биологических обектов именно на таких частотах приема.

В работах [1.12,1.15,1.22,1.24,1.26,1.27] были проведены модельные экспериментальные измерения теплового акустического излучения.

В работе [1.12] с целью опробования метода акустотермографии зарегистрировалось изменение акустояркостной температуры акрилового образца, который предварительно выдерживался в воде при разных температурах (0°С, 20°С и 38°С).

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

В работах [1.11,1.22,1.28] было рассмотрено влияние рассеяния на формирование акустояркостной температуры. В работах [1.22,1.28] влияние рассеяния было рассмотрено теоретически и экспериментально. Вклад в акустояркостную температуру дают как волны, поглощенные в слое, так и рассеянные. Если все рассеянные волны будут поглощены в областях с одинаковой температурой, то коэффициент поглощения, будет определяться как сумма коэффициента поглощения в тканях и коэффициента рассеяния. Экспериментальные исследования влияния рассеяния на величину АЯ температуры (в том числе и на биологических образцах) показали, что при малых градиентах температуры в исследуемой среде рассеяние дает малый вклад в акустояркостную температуру.

Совместно с исследованиями биологических моделей проводились in vivo измерения теплового акустического излучения тела человека. В работах [1.29-1.35] были зарегистрированы изменения внутренней температуры при нагреве или охлаждении разных частей (бицепса, кисти, икроножной мышцы) тела человека. В работах [1.29-1.31] восстанавливали профиль температуры кисти руки человека. Также с помощью АЯ Термометрии было зарегистрировано повышение температуры бицепса и мышц предплечья после физической нагрузки [1.24]. (Следует отметить, что динамика нагрева бицепса при физической нагрузке подробно исследована с помощью СВЧ-радиометрии в работе [1.34]). При регистрации теплового акустического излучения кисти руки человека, погружаемой в воду при различной температуре, было показано [1.21,1.26], что акустояркостная температура ладони зависит от температуры воды. Это позволяет отслеживать процессы терморегуляции в живом организме. В работе [1.32] исследовали с помощью теплового акустического излучения термодинамические параметры организма человека. При этом охлаждали икроножную мышцу. С одной стороны ноги устанавливали AT, к другой подносили охладитель, и в течение 15 минут проводили измерения акустояркостной температуры мышцы. По результатам измерений рассчитали коэффициент температуропроводности, значения которого по порядку величины оказались близки к общепринятым данным.

В работе [1.36] изучали корреляцию между показаниями акустотермометра, регистрирующего сигнал, излучаемый печенью, и уровнем сахара в крови при некоторых физиологических воздействиях. Испытуемому давали стандартные дозы глюкозы и измеряли акустояркостную температуру печени. Проведенный совместный анализ показал наличие корреляции между акустояркостной температурой и количеством сахара в крови.

Следует отметить, что при измерениях теплового акустического излучения тела человека требуется достаточная точность для определения его акустояркостной температуры. Можно привести несколько цифр, показывающих пределы применимости акустояркостных термометров. Например, перегрев злокачественной опухоли относительно здоровой ткани составляет в среднем 0.3-0.5 К реже 1.5 К [1.37], область нагрева из за влияния кровотока и диффузионных свойств тканей составляет величину порядка 1см. Это значит, что акустояркостная температура такой опухоли порядка 0.1-0.3 К. При гипертермии в онкологии необходимо нагреть опухоль приблизительно до 42°С и поддерживать температуру неизменной с точностью 0.3-0.5 К. При физиологических реакциях организма температурные изменения составляют около 1 К. Таким образом, чем выше точность AT, тем шире круг задач, в которых может быть использован метод АЯТ.

Исследование механизмов формирования акустического излучения и возможностей его регистрации позволили выделить преимущества АЯ Термометрии по сравнению с другими методами, а также поставить целый комплекс решаемых задач. Однако следует отметить, что вопросы повышения чувствительности АЯ Термометров до сих пор остаются актуальными.

Благодаря хорошей проникающей способности ультразвуковых волн в мегагерцовом диапазоне частот, АЯ Термометрия позволяет регистрировать излучение от структур, находящихся на глубине в несколько сантиметров от поверхности биологического объекта [1.17]. В то же время в таком диапазоне частот длины волн составляют миллиметры и доли миллиметров, что позволяет реализовать высокую направленность акустической антенны при ее сравнительно малых размерах, и следовательно, высокую разрешающую способность.

Большая проникающая способность и высокая направленность антенн позволяет поставить перед АЯ Термометрией задачи локализации нагретых источников внутри исследуемых объектов и картирования температурных полей.

Задача восстановления профиля объекта по результатам измерений является обратной. Различают одномерные, двумерные и трехмерные обратные задачи. Зависимость распределения внутренней температуры только от глубины объекта, т.е. исследуемую среду можно представить в виде параллельных поверхности объекта слоев, например как в [1.10], то восстановление температуры является одномерной обратной задачей[1.38]. Двумерная обратная задача - это восстановление распределения внутренней температуры в некоторой области на плоскости[1.39-1.43]. Восстановление объемного пространственного распределения внутренней температуры - трехмерная обратная задача[1.44,1.45].

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

В работах [1.11-1.14] рассмотрена одномерная обратная задача только для однородной среды, в которой принималось, что вариации плотности среды и скорости звука меньше 5 %, а вариации коэффициента поглощения на фиксированной частоте не более 15 %. В этих работах показано, что при изменении распределения внутренней температуры объекта меняется частотная зависимость акустояркостной температуры. Это дает возможность отличать одно распределение внутренней температуры от другого. Таким образом, в работах [1.10,1.11] сформулирована одномерная обратная задача акустотермографии, но не указано способов ее решения.

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

Вариант решения трехмерной обратной задачи как суммы двумерных задач был предложен в работах [1.44,1.45].

Коротко рассмотрим математические алгоритмы для решения обратных задач. Методы реконструкции внутренних характеристик исследуемых, в том числе, и биологических объектов описаны в работах [1.46,1.47,1.48]. Например, в рентгеновской вычислительной томографии регистрируют поглощение рентгеновского излучения в организме человека. Источники излучения находятся в одной плоскости в разных точках вне организма. Лучи проходят сквозь исследуемый объект под разными углами. Приемники, расположенные в той же плоскости вне организма, измеряют интенсивность выходящего из него излучения. Восстанавливают распределение коэффициента поглощения рентгеновского излучения в тканях организма. Разработано множество восстанавливающих искомые характеристики алгоритмов, однако вид интегрального уравнения в этих задачах отличается от интегрального уравнения акустотермографии. Поэтому выводы вычислительной томографии напрямую не могут использоваться в акустотермографии.

В 1998 г. И.П.Боровиковым и др. [1.49] рассмотрена возможность при решении обратных задач акустотермографии использовать априорную информацию о том, что восстанавливаемое распределение внутренней температуры должно удовлетворять уравнению теплопереноса с учетом кровотока. В работах [1.17,1.50] авторами рассмотрен регуляризирующий алгоритм (по теории регуляризации - вариант локальной регуляризации), позволяющий восстанавливать пространственное распределение внутренней температуры объекта. Однако решение обратной задачи таким методом предполагает приближение абсолютно тонкого акустического луча, однако такое приближение при поставленных в работах начальных условиях, а именно: размерах апертуры 3-5 длин волн и глубине проникновения порядка 100 длин волн ♦ физически нереализуемо. По оценкам соотношение между шириной пучка и поперечными размерами исследуемого пространства должно составлять в дальней зоне порядка 10-15%.

В работах [1.51,1.52] предлагалось для получения набора акустояркостных температур при решении обратной задачи в плоскослоистой среде использовать повороты акустического приемника на разные углы. Интегральное уравнение для акустояркостной температуры, получаемое в этом случае, аналогично уравнению переноса. Заметим, что в данном случае двумерная обратная задача сводится к одномерной. В работе [1.53] был рассмотрен алгоритм локализации нагретых областей небольшого размера. Источником теплового акустического излучения служила длинная тонкая трубка (диаметр трубки меньше диаметра пьезопреобразователя) с нагретой жидкостью. Трубку располагали в аквариуме с водой на некоторой глубине параллельно поверхности воды. Снимали показания двух датчиков, пересекающиеся аппаратные функции которых расположены в плоскости, перпендикулярной оси трубки. Датчики перемещали по поверхности воды в той же плоскости. Таким образом, решали задачу локализации на плоскости практически точечного источника излучения, находящегося в непоглощающей среде. Если оба датчика показывают отличную от нуля акустояркостную температуру, то в области пересечения аппаратных функций существует тепловой источник. Заметим, что поставленная задача локализации на плоскости точечного источника отличается от задачи восстановления распределения температуры в области.

Как уже сказано выше, глубинную температуру биологических объектов можно измерять с помощью неинвазивного метода СВЧ-радиометрии, основанного на измерении шумового электромагнитного излучения [1.54]. Измеряют радиояркостную температуру [1.55] - аналог акустояркостной. Восстановление внутреннего распределения температуры по известным значениям является некорректной обратной задачей. Для решения одномерной обратной задачи СВЧ-радиометрии предложено использовать метод регуляризации по А.Н.Тихонову (К.П.Гайкович и др. [1.56,1.57]). Авторы восстанавливают температурные профили в слоистых модельных средах и в тканях организма человека по значениям радиояркостной температуры, измеренным на разных частотах. В частности, в работе [1.57] восстановлен температурный профиль, имеющий характерный максимум в слое, где расположена опухоль. Отметим, что мультиспектральное зондирование с большой детальностью практически трудно реализуемо, поскольку количество приемников сигнала определяется количеством выбранных для измерений частот. Как результат, важный для дальнейшего анализа, отметим, что решена одномерная обратная задача для СВЧ-радиометрии - задача, в которой интегральное уравнение аналогично уравнению переноса для акустотермографии.

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

Одно из решений задачи предложено в первой главе настоящей работы. Учитывая возможности АЯ Термометрии, локализацию нагретых источников внутри исследуемой среды можно проводить, используя многоканальную схему[1.58] с многоракурсным сканированием исследуемой среды, которая, как будет показано в первой главе, выгодно отличается от предыдущих схем, как по чувствительности, так и по общему времени сканирования. Для повышения чувствительности впервые применялась компенсационная схема приема акустического сигнала. Такая схема приема обеспечивает по сравнению с используемыми ранее модуляционными схемами почти двукратный выигрыш в чувствительности прибора.

Дтя восстановления профиля температур по полученным данным впервые использовался один из прямых томографических алгоритмов — алгоритм алгебраической реконструктивной томографии (APT). Этот алгоритм, как будет показано в первой главе, наилучшим образом подходит в случае, когда массивы данных неполны, а сами данные сильно зашумлены.

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

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

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

Однако следует отметить, что вопросы повышения чувствительности остаются актуальными. Это связано с тем, что АЯ термометр регистрирует акустояркостную температуру[1.15], которая как правило значительно ниже термодинамической температуры исследуемого объекта.

Одним из возможных способов повышения чувствительности приборов является применение фокусированных антенн. Площадь сечения «пучка» фокусированной антенны изменяется в зависимости от расстояния до антенны и достигает своего минимума в фокусе. При этом происходит перераспределение принятой антенной интенсивности сигнала, пришедшего из различных точек пространства, и максимальный сигнал приходит из фокусной «перетяжки». Чувствительность прибора в этом случае, увеличивается за счет коэффициента усиления фокусированной антенны.

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

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

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

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

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

Описываемые в литературе теоретические исследования корреляционных свойств теплового акустического излучения можно разделить на две части: исследование собственно корреляционных свойств теплового излучения в среде и вопросы регистрации коррелированных сигналов в ультразвуковом диапазоне конкретными физическими приемниками. Реальные изучаемые объекты являются неоднородными, как по температуре, так и по поглощению. Например, в ряде биомедицинских приложений, в исследуемой среде (полупространстве) существует неизменное во времени распределение внутренней температуры и коэффициента поглощения ультразвука. Задача расчета корреляционной функции для неоднородной среды требует своего решения. При этом пока не рассчитаны и коррелированные сигналы, принимаемые конкретными преобразователями с учетом их аппаратных функций.

Коррелированный сигнал можно измерить с помощью двух пьезопреобразователей, на выходе которых устанавливают перемножитель. Это дает возможность получить произведение давлений, регистрируемых на двух пьезоприемниках. В главе 1 рассмотрены измерения интенсивности теплового акустического излучения, т.е. измерения среднего квадрата давления на пьезопреобразователе. В то же время по аналогии с радиолокацией в работе Р.Хессемера и др. [1.59] (патент США, 1983 г.) рассматривали возможность получения распределения внутренней температуры с использованием корреляционного приема сигнала несколькими датчиками.

В ряде более поздних работ разных авторов [1.60-1.62] эта идея получила дальнейшее развитие. В данных работах теоретически обсуждали возможность регистрации корреляционной функции второго порядка теплового акустического излучения парами датчиков. При этом предлагали использовать полученные данные для восстановления внутренней температуры. С этой же целью также обсуждали возможность использовать корреляционный момент четвертого порядка [1.63].

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

Рассматривая современный уровень исследований по аналогичным проблемам в смежных с акустикой областях, отметим, что теоретические исследования корреляционных свойств теплового излучения проведены для электромагнитного излучения [1.66-1.67]. В работе [1.66] была рассчитана пространственно-временная корреляционная функция излучения абсолютно черного тела, а в работе [1.67] изучен прием с помощью антенн теплового электромагнитного излучения в однородной среде. Хотя теория корреляционных свойств теплового электромагнитного излучения достаточно хорошо разработана, но тепловое акустическое излучение по физике явления существенно отличается от электромагнитного. Отметим наиболее существенные отличия измерений теплового акустического излучения от регистрации шумового электромагнитного излучения, проводимого, например, в радиоастрономии [1.20].

В радиоастрономии измеряют шумовое излучение от нагретых источников, расположенных в среде с абсолютной термодинамической температурой, близкой к нулю. В акустотермографии температура источников теплового акустического излучения из биологического объекта (например, из тела человека) ненамного (максимум на 15 - 20 градусов, т.е. на 5-7 %) отличается от температуры окружающей среды. При прочих равных условиях это увеличивает погрешность акустических измерений.

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

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

Напомним, что в ряде работ [1.60,1.61] предлагали восстанавливать внутреннюю температуру с помощью разновидностей корреляционного приема теплового излучения. Варианты конструкции пассивного акустического термотомографа обсуждали теоретически в ряде работ В.И. Миргородский, С.В.Пешин и др. [1.60,1.61].

В работах [1.60,1.61] для восстановления внутренней температуры предлагали измерять корреляционную функцию теплового акустического излучения четвертого порядка. Для проверки возможностей этого метода был проведен модельный эксперимент [1.60]: с помощью четырех микрофонов регистрировали шумовой звуковой сигнал, излучаемый громкоговорителем размером около 50 см. Характерный размер восстановленного пространственного распределения интенсивности излучения (по уровню 0,5 от максимума) был приблизительно в 2,5 раза больше (120 см). Подобное увеличение размеров области шумового источника получили и в численных экспериментах. Нам представляется, что потенциально подобный метод позволяет £ локализовать в пространстве достаточно интенсивные небольшие шумовые акустические источники.

Дтя оценки возможного пространственного разрешения локализации источника некогерентного звукового сигнала В.И. Миргородским и др. в работах [1.68] был проведен еще один модельный эксперимент. Измерения корреляционного момента четвертого порядка звукового шумового сигнала со средней частотой ~ 1 КГц проводили пятью микрофонами, расположенными в вершинах двух тетраэдров с общим основанием. Шумовой источник размером около 1 мм был локализован в области с линейньш размером менее 40 см.

Отметим, что в вопросе применения корреляционного приема теплового акустического излучения имеются разные точки зрения. В ряде работ [1.60,1.61] предложено использовать корреляционный прием теплового акустического излучения для восстановления внутренней температуры в теле человека, а в работах [1.64,1.65,1.69] рассмотрена возможность восстанавливать таким методом распределение коэффициента поглощения в биологических объектах. Однако результаты работ предполагают наличие малых по сравнению с апертурой приемных антенн объектов.

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

В третьей главе проведены исследования взаимодействия оптического излучения с мягкими биологическими тканями методами АЯТ. Актуальность таких исследований связана с тем, что в настоящее время активно развиваются различные научные направления для решения круга задач медицинской диагностики, основанные как на чисто активных или пассивных методах, так на их комбинации. Так, например, получивший в настоящее время мощное развитие такой метод как оптоакустическая томография[1.70] основан на активном оптическом импульсном облучении исследуемого объекта и пассивной регистрации акустического импульса, генерируемого в результате быстрого расширения поглощающей оптической неоднородности. Преимущество этого метода по сравнению с оптическими в большей порядка нескольких сантиметров глубине диагностического зондирования. Такой метод позволяет проводить диагностику внутренних заболеваний и воспалительных процессов, в том числе залегающих на глубинах до нескольких сантиметров, которые сопровождаются изменением оптических свойств относительно здоровых тканей [1.71]. Однако влияние импульсно - периодического лазерного излучения на свойства биологических тканей, обусловленное нагревом поглощающей области, накладывает ограничения на временные и энергетические параметры воздействия. В этом случае возможности акустояркостной термометрии позволяют проводить контроль изменения температуры внутри биологических тканей.

Другим интересным аспектом для медицинской практики является совмещение акустояркостной термометрии и использование непрерывных лазеров, которое открывает новые возможности для современной медицинской диагностики. Во — первых важное значение имеет применение методов АЯТ для контроля внутренней температуры при различных гипертермических процедурах в том числе лазерной ♦ гипертермии [1.72]. При гипертермии нагревают опухоль до некоторой оптимальной температуры, чаще всего до 42° С. Совместно с нагревом проводят радиотерапию или химиотерапию пораженной ткани. При этом важной проблемой является безболезненный неинвазивньй контроль глубинной температуры. Хотя точность таких измерений ниже, чем точность традиционных инвазивных методов, неинвазивный контроль желателен, если с его помощью можно измерить температуру на глубине 3-8 см, в объеме около 1 см с точностью 0,3-0,5 К. Во-вторых, применение непрерывного лазерного излучения позволяет «подсвечивать» область интересную для АЯТ диагностики, т.е. повышать температуру больной ткани, поскольку в большинстве случаев оптическое поглощение больной и здоровой ткани различно из-за образования паразитных сосудов, изменения структуры самой ткани и т.д.

Впервые применение метода АЯТ для исследования эффекта теплового воздействия при лазерной гипертермии было описано в работе [1.73] В этих работах исследовались образцы свиная и говяжья печень. Исследуемый образец представлял собой прямоугольный параллелепипед с характерным линейным размером в несколько сантиметров. Одна из граней образца находилась в акустическом контакте с AT. Нагрев осуществлялся лазерным излучением с длиной волны 1.02 мкм, вводимым внутрь биоткани через волоконный световод. Нагрев проводился около 5 мин, мощность лазера изменялась от 1.6 до 5 Вт. В работе было отмечено повышение температуры в образце на ~ 100 К. Регистрация акустояркостной температуры после выключения нагрева проводилась в течение 20 - 30 мин с точностью до ~ 5 К.

В третьей главе показана возможность применения методов АЯТ для контроля температуры при оптическом воздействии импульсно - периодического и излучения Nd:YAG лазера с длиной волны 1.064мкм, непрерывном воздействии излучения Nd:YAG лазера на однородные и неоднородные биологические ткани. Также показано, что с помощью применения многоканального АЯ термометра возможно проводить поперечную локализацию лазерного воздействия, а также контроль распространения тепловой волны, порождаемой лазерным воздействием.

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

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

Положения, выносимые на защиту:

1) Картирование поля внутренних температур можно проводить методами АЯ Термометрии используя многоканальный прием акустического излучения с последующей реконструкцией данных с помощью томографических алгоритмов. При этом локализация температурных неоднородностей (объектов) биологических тканей возможна при разнице термодинамических температур порядка 2-3 К и характерных размерах этих объектов порядка 1 см. Точность восстановления томографического изображения, т.е. разница между истинным пространственным распределением неоднородностей и восстановленным изображением может быть порядка 2-4 мм. Меньшие размеры объектов повлекут за собой уменьшение их АЯ Температуры регистрируемой АЯ Термотомографом.

2) Для картирования поля температур с малыми объектами возможно использование фокусированных антенн. Локализация малых объектов возможна за счет перераспределения чувствительности приема вдоль главной оси антенны.

3) Методы АЯ Термометрии могут быть применимы для контроля различных гипертермических процедур, в том числе при лазерной гипертермии.

4) С помощью совмещения методов АЯ Термометрии и лазерной гипертермии возможно проводить измерения температуропроводности, акустического и оптического поглощения, мягких биологических тканей.

Материал диссертации изложен в 4 главах и заключении. Диссертация содержит 164 страниц, 48 рисунков и 3 таблицы, список литературы содержит 132 наименований. Результаты работы докладывались на:

- IV Всероссийской конференции по биомеханике г. Н.Новгород, 1998;

- Научной конференции по радиофизике, посвященной 80 - летию ИНГУ, г. Н.Новгород, 1998;

I Всероссийской научно - технической конференции "Компьютерные технологии в науке, проектировании и производстве", г. Н.Новгород, 1999;

II международной научно - технической конференции «Проблемы и прикладные вопросы физики»;

3 - ей Научной конференции по радиофизике, г. Н.Новгород, 1999;

- 4-ой Нижегородской сессии молодых ученых, г. Н.Новгород, 1999;

- 1 - ой международной конференции молодых ученых, Young Scientists on Laser Optics (LO-YS), г. С.- Петербург, Россия, 2000;

4-ой Научной конференции по радиофизике, г. Н.Новгород, 2000;

5-ой Нижегородской сессии молодых ученых, г. Н.Новгород, 2000;

- Biomedical Optics (BIOS), San Jose, 2001;

- XI сессии Российского Акустического Общества, Москва, 2001;

- International School for Young Scientists on Optics, Laser Physics & Biophysics, (SFM' 01), Saratov, Russia, 2001;

- Biomedical Optics (BIOS), San Jose, 2002;

- 5-ой Научной конференции по радиофизике, г. Н.Новгород, 2001;

- 6-ой Научной конференции по радиофизике, г. Н.Новгород, 2002.

Основные результаты диссертационной работы опубликованы в 26 печатных работах, приведенных в списке литературы.

Основные результаты работы следующие:

1. Экспериментально показана возможность неинвазивного картирования поля температурных неоднородностей однородных и неоднородных биоподобных сред с поглощением по их собственному тепловому акустическому излучению методом многоканальной акустояркостной (АЯ) термометрии.

2. Экспериментально продемонстрирована возможность регистрации реакции человеческого организма на различные физиологические нагрузки in vivo многоканальным акустояркостным термотомографом.

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

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

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

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

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

ЗАКЛЮЧЕНИЕ

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

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

1. Акустотермометр работает следующим образом. Акустическое излучение регистрируется приемной антенной (ПА), далее усиливается УС (усилитель высокой частоты) и детектируется квадратичным детектором. После детектирования сигнал поступает на интегратор и проходит на регистрирующее большого количества приемных датчиков в ограниченном 29

2. Структурная схема канала многоканального акустотермографа. 30

3. Антенный блок многоканального акустотермотомографа 32