CC BY
55
УДК 612.15
М. З. Смирнов, А. Е. Пушкарева
ОСОБЕННОСТИ ИСПОЛЬЗОВАНИЯ УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ ДЛЯ МОДЕЛИРОВАНИЯ ВОЗДЕЙСТВИЯ ЛАЗЕРНОГО ИЗЛУЧЕНИЯ
НА КОЖУ
Проведено моделирование воздействия лазерного излучения на кожу. Для учета влияния кровотока в уравнение теплопроводности добавлены источники тепла, обусловленные потоком крови. Показано, что учет температурной зависимости скорости кровотока приводит к появлению максимума температуры и дальнейшей ее стабилизации.
Ключевые слова: математическое моделирование, уравнение теплопроводности, кровоток, кожа.
Введение. Воздействие лазерного излучения на биологическую ткань, например кожу, обычно сопровождается локальным повышением ее температуры. Нагрев, в свою очередь, вызывает ответные реакции организма. В частности, отмечаются сосудистые реакции, проявляющиеся даже на некотором расстоянии от области воздействия. Происходит местное расширение сосудов, снижается вязкость крови, что приводит к увеличению скорости ее течения.
В норме линейная скорость кровотока в капиллярах (диаметр d = 8 мкм) очень низкая и составляет примерно 0,1 мм/с, а в мелких венах (d = 0,75 мм) — 0,05 м/с. При локальном нагреве она может возрастать в несколько раз. Кроме того, скорость кровотока не мгновенно реагирует на изменения температуры. Время задержки довольно велико и составляет, согласно [1], 60—90 с. Такие изменения скорости течения крови влияют на результат воздействия лазерным излучением, особенно если результат зависит от степени термического повреждения биоткани, поскольку кровоток может быть дополнительным и достаточно эффективным механизмом отведения тепла от области обработки. Оценка влияния данного эффекта может быть проведена с помощью численного моделирования динамики процессов, происходящих при воздействии лазерного излучения на кожу.
В настоящей работе проведено моделирование лазерного нагрева кожи для трех случаев. В первом случае расчеты проводились без принятия во внимание кровотока, во втором — задавалась постоянная скорость течения крови, в третьем учитывалась зависимость температуры от скорости кровотока. Результаты, полученные для этих трех случаев, позволяют оценить влияние кровотока на лазерный нагрев кожи.
Модель. В расчетах была воспроизведена обработка кожи лазерным излучением с длиной волны 810 нм. Рассматривался импульсный режим с длительностью одного импульса 1 с и частотой повторения 0,8 Гц. Плотность энергии импульса 1,6 Дж/см , диаметр пятна на коже 7,5 см, длительность обработки 250 с. Считалось, что начальная температура кожи равна 37 °С и в процессе обработки поверхность кожи охлаждается воздушным потоком, температура которого равна 20 °С. Такой режим представляется характерным для лазерной обработки жировой ткани.
С учетом строения кожи была разработана ее модель, включающая в себя четыре слоя — эпидермис (толщина 70 мкм), базальный слой (15 мкм), дерма, или истинная кожа (3 мм), и подкожная жировая клетчатка (5 мм). Тип кожи — наиболее распространенный — 2-й по Фицпатрику. Кровеносные сосуды располагаются в дерме и жире. Скорость течения крови в сосудах разного типа различна и зависит от диаметра сосуда. Для рачетов были приняты следующие значения плотности потока крови: в дерме 15 мл/(мин-100 г), в подкожной
Особенности использования уравнения теплопроводности для моделирования воздействия 65
жировой клетчатке 5 мл/(мин-100 г) [1]. Скорость кровотока практически неизменна при температуре меньше 41,5 °С и линейно возрастает при дальнейшем нагреве.
Моделирование динамики оптических и тепловых процессов диффузии при воздействии оптического излучения на кожу проводилось на основе теории переноса излучения, а также теории теплопроводности. Распределение освещенности в среде рассчитывалось при помощи уравнения переноса излучения в диффузионном приближении [2]:
V2 ф5 (г)-3|да | 1г ф5 (г)+3| | 1г Е (г, 5о)-З^ V(E (г, ^о) = 0. (1)
где ф5 (г) — диффузная составляющая освещенности в точке г = (х, у, г), — коэффициент поглощения, 15 — коэффициент рассеяния, 11г = | +(1 -g) | — транспортный коэффициент затухания, Е (г, ,?о) — коллимированная освещенность в точке г, ,?о — направление распространения первичного луча, g — фактор анизотропии рассеяния.
На основе распределения освещенности внутри кожи была рассчитана объемная плотность источников тепла в среде:
0(г) = |аф(г)Ео/хр , (2)
где ф(г) — полная освещенность в точке г, отнесенная к единичной плотности мощности излучения на поверхности кожи, Ео — плотность энергии излучения, т р — длительность импульса.
Далее был проведен расчет распределения температуры в каждом слое (в эпидермисе и базальном слое), для чего использовалось уравнение теплопроводности:
рс ^^ ^(^Т(г, 0)+0(г), (3)
ОТ
где р — плотность ткани, с — теплоемкость, I — время, к = арс — теплопроводность, а — температуропроводность, Т — температура.
Данное уравнение применимо только для слоев, где нет кровотока — эпидермиса и ба-зального слоя. Для дермиса и жира следует дополнительно учитывать источники тепла, обусловленные потоком крови, Qb (, Т) . В этих слоях уравнение теплопроводности принимает следующий вид:
рсНзТ^^(^Т(г,0)+Q(г)+Qb (,Т) . (4)
Слагаемое Qb (,Т) можно записать следующим образом [2]:
Qb (^ Т)=рс [рь/(Ть -Т)], (5)
где рь — плотность крови, Ть — ее температура (36,5 °С), / — плотность потока крови в ткани, представляющая собой объем крови, переносимый в секунду в 1 грамме соответствующей ткани. В стационарных условиях плотность потока равна / и различна для дермиса и жира. Кроме того, стационарная плотность потока зависит от температуры, т.е. можно записать /о = /о(Т).
В настоящей работе рассматриваются нестационарные условия воздействия излучения на кожу. В этом случае / не равно /о(Т). Как уже было отмечено выше, при изменении температуры поток крови изменяется не мгновенно, а с задержкой £мау=бо—9о с. Для учета данной задержки используем следующее уравнение для потока крови:
±/=т/, (6)
д 1 ^е1ау
Таким образом, для расчета температуры кожи с учетом скорости кровотока и ее изменения вследствие нагрева, а также времени задержки этого изменения с использованием выражений (4)—(6) была получена система уравнений.
Результаты и их обсуждение. Как было указано выше, моделирование воздействия лазерного излучения на кожу проводилось для трех случаев. Для сравнения были получены графики зависимости температуры базального слоя (на глубине 83 мкм), дермы (на глубине 2 мм) и подкожной жировой клетчатки (на глубине 5 мм) от времени (рис. 1, 2), зависимости температуры от глубины для всех случаев (рис. 4) а также для третьего случая — зависимость плотности потока крови от времени (рис. 3).
На рис. 1, а приведены графики зависимости без учета кровотока, на рис. 1, б — с учетом постоянной скорости кровотока. а)
0
100
200
0
100
200
1, с
1, с Рис. 1
Временная зависимость температуры базального слоя (1) представляет собой чередование максимумов и минимумов. Такой характер объясняется импульсным режимом обработки, а также малым временем тепловой релаксации базальной мембраны. Зависимости температуры дермы (2) и жировой клетчатки (3) от времени не показывают подобных колебаний, поскольку происходит их сглаживание вследствие большей глубины указанных слоев.
Введение в расчет постоянной скорости кровотока дает снижение температуры всех слоев на несколько градусов, однако не изменяет характера временной зависимости температуры относительно первого случая (рис. 1). Нарастание температуры базального слоя происходит сначала достаточно резко, затем более медленно. В дерме скорость увеличения температуры меньше, а в подкожной жировой клетчатке начало обработки практически не сопровождается нагревом слоя.
Введение температурной зависимости скорости кровотока резко изменяет тепловую динамику кожи (рис. 2, здесь 1 — базальный слой, 2 — дерма, 3 — жировая клетчатка). Сначала нарастание температуры ничем не ограничено. После достижения температуры 41,5—42 °С начинает возрастать плотность потока крови (рис. 3, 1 — дерма, 2 — жировая клетчатка) и через некоторое время температура начинает снижаться, в дальнейшем приводя к ее стабилизации. Это приводит в свою очередь к уменьшению и стабилизации скорости кровотока. Изменение температуры в дерме и подкожной жировой клетчатке вследствие
Т, °С 65
55
45
35
100 Рис. 2
200
1, с
Особенности использования уравнения теплопроводности для моделирования воздействия 67
изменения потока крови оказывает влияние также на динамику температуры базального слоя (рис. 2). Максимальная температура данного слоя намного больше, чем дермы и жира, так как за счет содержащегося в нем меланина излучение поглощается сильнее. Поскольку поток крови в базальном слое отсутствует, не происходит непосредственного отведения тепла за счет кровотока.
Температурная зависимость потока крови также оказывает влияние и на распределение тепла по глубине (рис. 4, здесь 1 — без учета кровотока, 2 — скорость кровотока постоянна, 3 — скорость кровотока зависит от температуры). Видно, что характер охлаждения кожи при движении в глубь среды в третьем случае отличается от первых двух. Наибольшая скорость спада температуры в этом случае наблюдается в верхней части дермы, а не в подкожной жировой клетчатке. Интересно, что на глубине 8 мм температура одинакова и равна начальной температуре ткани 37 °С.
Заключение. Таким образом, на основе полученных результатов можно сделать вывод о том, что для наиболее достоверного расчета температуры в коже необходимо учитывать как скорость кровотока, так и ее зависимость от температуры, поскольку это явление оказывает большое влияние на результат воздействия лазерного излучения. Учет данной особенности живых биологических тканей осуществляется посредством введения в уравнение теплопроводности источников тепла, обусловленных потоком крови. В отличие от моделей, принимающих во внимание скорость кровотока постоянной или не учитывающих ее вообще, модель, в которой скорость кровотока рассматривается как функция от температуры, позволяет дать наиболее полную оценку динамики температуры в коже.
СПИСОК ЛИТЕРАТУРЫ
1. Sekins K. M., Emery A. F. Thermal Science for Physical Medicine. Therapeutic Heat and Cold / Ed. by J. F. Lehmann Baltimore: Williams & Wilkins, 1990. P. 62.
2. Star W. M. Diffusion theory of light transport // Optical-thermal response of laser-irradiated tissue / Ed. by A. J. Welch and M. J. C. van Gemert. NY: Plenum Press, 1995. P. 144.
Михаил Захарович Смирнов Александра Евгеньевна Пушкарева
Сведения об авторах канд. техн. наук, доцент; Palomar Medical Technologies Inc.; E-mail: msmirnov@palomarmedical.com
канд. техн. наук; Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кафедра лазерной техники и биомедицинской оптики; E-mail: alpushkareva@yandex.ru
Рекомендована кафедрой оптико-электронных приборов
Поступила в редакцию 29.05.08 г.
