CC BY
8ISSN 0868-5886 НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2018, том 28, № 2, c. 20-24
^^ПРИБОРОСТРОЕНИЕ ФИЗИКО-ХИМИЧЕСКОЙ БИОЛОГИИ =
УДК 535.015 +57.087.1 + 004.418 © Д. А. Кравчук, И. Б. Старченко
МОДЕЛИРОВАНИЕ ПРОЦЕССА НАСЫЩЕНИЯ КИСЛОРОДОМ БИОЛОГИЧЕСКИХ ТКАНЕЙ С ПОМОЩЬЮ ОПТОАКУСТИЧЕСКОГО МЕТОДА
Работа посвящена моделированию оптоакустических сигналов от эритроцитов с учетом кислородонасыще-ния для лазерного излучения с длиной волны 700 нм. Смоделированы пространственные организации тканей, непересекающихся, случайно распределенных смесей оксигенированных и дизоксигенированных эритроцитов в двумерном пространстве. Представлена теоретическая модель для изучения влияния SO2 крови на оптоакустические сигналы. Для этого рассматривались смеси оксигенированных и дизоксигенирован-ных эритроцитов. Полученный сигнал от эритроцитов был рассчитан с использованием принципа линейной суперпозиции для сигналов, излучаемых отдельными эритроцитами. Было замечено, что амплитуда опто-акустического сигнала возрастала по мере уменьшения SO2 для оптического излучения 700 нм.
Кл. сл.: оптоакустический сигнал, оксигенация, дизоксигенация, кислородонасыщение, эритроциты, спектральная плотность мощности, лазер
ВВЕДЕНИЕ
Предположим, что короткий лазерный импульс облучает тканевую жидкую среду, которая поглощает свет под термоупругим расширением. В результате он излучает импульс акустического давления и может быть обнаружен ультразвуковым преобразователем. Обнаруженный сигнал полезен для исследований, а также для получения изображения тканевой среды — оптоакустического (ОА) изображения. Контраст ОА-изображения обычно лучше, чем обычное изображение обратного рассеяния ультразвука, поскольку разница оптического свойства поглощения намного сильнее, чем различие поглощения в различных областях ткани. В ОА локализованная информация может быть получена с помощью фокусированного ультразвукового преобразователя. Более того, более глубокие области тканей могут быть отображены с использованием ОА-методов (по сравнению с оптическим способом визуализации), поскольку он использует волны давления, рассеяние которых на два-три порядка меньше, чем у света, и, таким образом, распространяется на большее расстояние.
Молекулы гемоглобина, которые присутствуют в эритроцитах, ведут себя как доминирующие поглощающие свет хромофоры и формируют излучение ОА-сигнала из эритроцитов. Недавно были созданы пространственные карты оксигенирован-ных эритроцитов и дезоксигенированных эритроцитов для малой визуализации мозга животных с использованием ОА-метода [1].
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
Аналитическое выражение оптоакустического поля, генерируемого жидкой сферой при условии термического ограничения, можно записать в виде [2-10]
Pf ( q ) =
i^ßI0vsa'
CPr
[sin q - qcos q]
ikf (r-a)
(i - /3)
sin
- cos q + i/3vsmi
Вышеприведенное решение может быть получено путем решения волнового уравнения для поля давления, возникающего из-за поглощения света, в частотной области и согласования граничных условий на сферической границе [11-14]. В формуле л, СР, в и а — коэффициент поглощения, изобарическая удельная теплоемкость, коэффициент теплового расширения и радиус поглощающей сферы соответственно. Далее р = / pf
и V = уэ / vf — безразмерная плотность и скорость
звука (у5 — скорость звука в эритроцитах, Vf —
скорость звука в плазме крови) соответственно. Плотность и скорость звука в эритроцитах были взяты как рх = 1092 кг/м3 и = 1092 м/с соответственно. Числовые значения тех же самых коли-
х
X
2
честв окружающей среды (плазмы крови) были выбраны как Ру = 1005 кг/м3 и у/ = 1498 м/с [7-9].
Безразмерная частота определяется как д = = са ¡уш = кра, а — радиус поглощающей сферы (эритроцита). При этом к и ку — волновые числа акустической волны в поглощающей области и жидкой среде соответственно. При этом ш — частота модуляции падающего оптического излучения с интенсивностью 10.
МОДЕЛИРОВАНИЕ
Рассмотрим двумерную модель в области 100x100 мкм, которая была равномерно облучена оптическим излучением. Моделирование в 2D вычислительно легче, чем 3D, и достаточно для понимания происходящих процессов. Размер области исследования был взят 100 х 100 мкм и состоял из смесей оксигенированных и дезоксигенирован-ных эритроцитов при условии неперекрытия частиц и гематокрите Н = 45 %.
На рис. 1 показаны моделируемые реализации тканей, соответствующие уровням SO2 1 %
и 99 %. Моделирование образцов тканей проводилось по ранее разработанной методике [10-12].
Соотношения, при которых оксигенированные и дезоксигенированные эритроциты были смешаны, фиксировали насыщенность кислородом образца, как ^02) = N0 = (N0 +
Численные значения физических свойств эритроцита, используемые в моделировании, представлены в [7]. Коэффициент поглощения среды внутри эритроцита, насыщенного кислородом, оценивали путем умножения концентрации молекул гемоглобина внутри клетки с молярным коэффициентом экстинкции для оксигенированных молекул гемоглобина [8].
Спектральные плотности мощности для различных уровней S02 при длине волны лазера 700 нм изображены на рис. 2 в широком диапазоне частот. Спектральная плотность мощности уменьшалась по мере увеличения ^02). Эта тенденция отчетливо проявляется до 70 МГц, но после этого они перекрываются для некоторых случаев. Спектральная плотность мощности при ^02) 0 % на каждой частоте почти в 6 раз больше, чем при ^02) 100%, т. к. в этом случае отношение ^0 : № = 1: 6.
Рис. 1. Моделируемые реализации тканей, соответствующие уровням 802 1% и 99% при 45% гематокрите (красные круги представляют собой оксигенированные эритроциты, а синие круги обозначают дезоксигенированные эритроциты)
22
Д. А. КРАВЧУК, И. Б. СТАРЧЕНКО
Спектральная плотность при (SO2) = 1 %
Спектральная плотность при (SO2) = 25 %
a)
IZ
О
ct a)
■ 30
40 30 20 10 0
0 10 20 30 40 50 60 70 f, МГц
90 100
50
40 -■■
20 -■■
10 -
Спектральная плотность при (SO2) = 50 %
0 10 20 30 40 50 60 70 f, МГц
90 100
"0 10 20 30 40 50 60 70 f, МГц
a)
■ 30
90 100
50
40 -■■
20 -
10 -
Спектральная плотность при (SO2) = 99 %
0 10 20 30 40 50 60 70 f, МГц
90 100
Рис. 2. Спектральная плотность мощности ОАС при облучении оптическим излучением с длиной волны 700 нм.
Уровни SO2 в %: а — 1, б — 25, в — 50, г — 99
ЗАКЛЮЧЕНИЕ
Разработан алгоритм для моделирования реализации 2D тканей, состоящих из оксигенированных эритроцитов и дезоксигенированных эритроцитов, с учетом поглощения лазерного излучения форменными элементами крови. Уровень кислородной насыщенности оксигенированных эритроцитов принимался равным 100 % и 0 % — для клеток без кислорода, а также в долях 25 % и 50 %. Доля оксигенированных и дезоксигенированных фиксировала кислородное насыщение (802) образца крови как (802) = N = (N0 + N13), где N и ^ представляют собой количество клеток с кислородом и без кислорода. Результаты моделирования показали, что средняя амплитуда сигнала ОАС монотонно уменьшалась при повышении уровня 802 для лазерного излучения 700 нм. Проведенные исследования позволяют продолжить работу над созданием проточного цитомера [13-15].
СПИСОК ЛИТЕРАТУРЫ
Wang L. V. Prospects of photoacoustic tomography // Med. Phys. 2008. Vol. 35, no. 12. P. 5758-5767. Doi: 10.1118/1.3013698.
Diebold G.J. Photoacoustic monopole radiation: Waves from objects with symmetry in one, two and three dimensions // Photoacoustic Imaging and Spectroscopy. Edited by L.V. Wong. Taylor & Francis Group, LLC, London, 2009. P. 3-17.
Saha R.K., Kolios M.C. A simulation study on photoacoustic signals from red blood cells // J. Acoust. Soc. Am. 2011. Vol. 129, no. 5. P. 2935-2943. Doi: 10.1121/1.3570946.
Кравчук Д.А. Экспериментальные исследования и моделирование процесса генерации оптоакустических волн // Электронный научный журнал "Инженерный вестник Дона". 2017. Т. 45, № 2. URL: ivdon.ru/ru/magazine/archive/n2y2017/4234. Кравчук Д.А. Теоретические исследования генерации
б
а
в
г
оптоакустических волн в жидкости цилиндрическими поглотителями // Электронный научный журнал "Инженерный вестник Дона". 2017. Т. 46, № 3. URL: ivdon.ru/ru/magazine/archive/N3y2017/4350.
6. Кравчук Д.А. Аналитический результат генерации оптоакустических волн для сферических поглотителей в дальнем поле // Электронный научный журнал "Инженерный вестник Дона". 2017. Т. 47, № 4. URL: ivdon.ru/ru/magazine/archive/n4y2017/4436.
7. Shung K.K., Yuan Y.W., Fei D.Y., Tarbell J.M. Effect of flow disturbance on ultrasonic backscatter from blood // J. Acoust. Soc. Am. 1984. Vol. 75, no. 4. P. 1265-1272. Doi: 10.1121/1.390733.
8. Zhang H.F., Maslov K., Sivaramakrishnan M., Stoic a G., Wang L. V. Imaging of hemoglobin oxygen saturation variations in single vessels in vivo using photoacoustic microscopy // Appl. Phys. Lett. 2007. Vol. 90. P. 1-3.
9. Savery D., Cloutier G. Effects of red blood cell clustering and aisotropy on ultrasound blood backscatter: A Monte Carlo study // IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2005. Vol. 52. P. 94-103.
10. Кравчук Д.А., Старченко И.Б. Математическое моделирование оптикоакустического сигнала от сферических поглотителей на примере эритроцитов // Известия Юго-Западного государственного университета. Серия: Управление, вычислительная техника, информатика, медицинское приборостроение. 2017. Т. 7, № 3. C. 101-107.
11. Кравчук Д.А., Старченко И.Б. Математическое моделирование оптоакустического сигнала от агрегированных эритроцитов для оценки уровня агрегации // Научное приборостроение. 2018. Т. 28, № 1. C. 30-36. URL: http://213.170.69.26/mag/2018/abst1.php#abst4.
12. Кравчук Д.А. Система проточной лазерной диагности-
ки жидкостей при генерации оптоакустического сигнала на рассеивателях сферической формы // Качество и жизнь. Москва, 2017. № 4. С. 74-78.
13. Старченко И.Б., Кравчук Д.А., Кириченко И.А. Прототип оптоакустического лазерного цитомера // Медицинская техника. 2017. № 5. C. 4-7.
14. Кравчук Д.А. О методе моделирования оптоакустиче-ских сигналов от источников сферической формы на примере эритроцитов // Качество и жизнь. Москва, 2017. № 4. С. 78-80.
15. Starchenko I.B., Kravchuk D.A., Kirichenko I.A. An Op-toacoustic Laser Cytometer Prototype // Biomedical Engineering. 2018. Vol. 51, no. 5. P. 308-312.
ФГАОУ ВО "Южный федеральный университет" ИНЭП, г. Таганрог
Контакты: Кравчук Денис Александрович, [email protected]
Материал поступил в редакцию 15.03.2018
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2018, Vol. 28, No. 2, pp. 20-24
THE MODEL FOR DETERMINING OXYGEN SATURATION OF BIOLOGICAL TISSUES WITH THE HELP OF AN OPTOACOUSTIC METHOD
D. A. Kravchuk, I. B. Starchenko
Southwest State University, Taganrog, Russia
The work is devoted to modeling of optoacoustic signals from erythrocytes taking into account oxygen saturation for laser radiation with a wavelength of 700 nm. The spatial organization of tissues, disjoint randomly distributed mixtures of oxygenated and dysoxygenated erythrocytes in two-dimensional space is modeled. A theoretical model for studying the effect of SO2 on optoacoustic signals is presented. For this, mixtures of oxygenated and dysoxygenated erythrocytes were considered. The received signal from erythrocytes was calculated using the principle of linear superposition for signals emitted by individual red blood cells. It was observed that the amplitude of the optoacoustic signal increased with decreasing SO2 for optical radiation of 700 nm.
Keywords: optoacoustic signal, oxygenation, dysoxygenation, oxygen saturation, erythrocytes, power spectral density, laser
REFERENСES
1. Wang L.V. Prospects of photoacoustic tomography. Med. Phys, 2008, vol. 35, no. 12, pp. 5758-5767. Doi: 10.1118/1.3013698.
2. Diebold G.J. Photoacoustic monopole radiation: Waves from objects with symmetry in one, two and three dimensions. Photoacoustic Imaging and Spectroscopy, edited by L.V. Wong, Taylor & Francis Group, LLC, London, 2009, pp. 3-17.
3. Saha R.K., Kolios M.C. A simulation study on photoacoustic signals from red blood cells. J. Acoust. Soc. Am., 2011, vol. 129, no. 5, pp. 2935-2943.
Doi: 10.1121/1.3570946.
4. Kravchuk D.A. [The pilot studies and process modeling of generation of optoacoustic waves]. Elektronnyy nauchnyy zhurnal "Inzhenernyy vestnik Dona" [Online scientific magazine "Engineering Bulletin of Don"], 2017, vol. 45, no. 2. URL: ivdon.ru/ru/magazine/archive/n2y2017/4234. (In Russ.).
5. Kravchuk D.A. [Theoretical researches of generation of optoacoustic waves in liquid cylindrical absorbers]. Elektronnyy nauchnyy zhurnal "Inzhenernyy vestnik Dona" [Online scientific magazine "Engineering Bulletin of Don"], 2017, vol. 46, no. 3. (In Russ.).
URL: ivdon.ru/ru/magazine/archive/N3y2017/4350.
6. Kravchuk D.A. [Analytical result of generation of optoacoustic waves for spherical absorbers in a distant field]. Elektronnyy nauchnyy zhurnal "Inzhenernyy vestnik Dona" [Online scientific magazine "Engineering Bulletin of Don"], 2017, vol. 47, no. 4. (In Russ.).
URL: ivdon.ru/ru/magazine/archive/n4y2017/4436.
7. Shung K.K., Yuan Y.W., Fei D.Y., Tarbell J.M. Effect of flow disturbance on ultrasonic backscatter from blood. J. Acoust. Soc. Am., 1984, vol. 75, no. 4, pp. 1265-1272. Doi: 10.1121/1.390733.
8. Zhang H.F., Maslov K., Sivaramakrishnan M., Stoic a G., Wang L.V. Imaging of hemoglobin oxygen saturation variations in single vessels in vivo using photoacoustic mi-
croscopy. Appl. Phys. Lett., 2007, vol. 90, pp. 1-3.
9. Savery D., Cloutier G. Effect s of red blood cell clustering and aisotropy on ultrasound blood backscatter: A Monte Carlo study. IEEE Trans. Ultrason. Ferroelectr. Freq. Control., 2005, vol. 52, pp. 94-103.
10. Kravchuk D.A., Starchenko I.B. [Mathematical simulation of an optikoakustichesky signal from spherical absorbers on the example of erythrocytes]. Izvestiya Yugo-Zapadnogo gosudarstvennogo universiteta. Seriya: Upravlenie, vychislitel'naya technika, informatika, medi-zinskoe priborostroenie [News of the Southwest state university. Series: Control, ADP equipment, informatics, medical instrument making], 2017, vol. 7, no. 3, pp. 101— 107. (In Russ.).
11. Kravchuk D.A., Starchenko I.B. [Mathematical modeling of the optoacoustic signal from aggregated erythrocytes to assess the level of aggregation]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2018, vol. 28, no. 1, pp. 30—36. Doi: 10.18358/np-28-1-i3036. (In Russ.).
12. Kravchuk D.A. [The system of flowing laser diagnostic of liquids in case of generation of an optoaudible tone on dif-fusers of spherical shape]. Kachestvo i zhizn' [Quality and life], Moscow, 2017, no. 4, pp. 74—78. (In Russ.).
13. Starchenko I.B., Kravchuk D.A., Kirichenko I.A. [Prototype optoacoustic laser cytomeasure]. Medizinskaya technika [Medical equipment], 2017, no. 5, pp. 4—7. (In Russ.).
14. Kravchuk D.A. [About a method of simulation of op-toaudible tones from sources of spherical shape on the example of erythrocytes]. Kachestvo i zhizn' [Quality and life], Moscow, 2017, no. 4, pp. 78—80. (In Russ.).
15. Starchenko I.B., Kravchuk D.A., Kirichenko I.A. An Optoacoustic Laser Cytometer Prototype. Biomedical Engineering, 2018, vol. 51, no. 5, pp. 308—312.
Contacts: Kravchuk Denis Aleksandrovich, kravchukda@sfedu. ru
Article received in edition: 15.03.2018
