МЕДИЦИНСКАЯ КИБЕРНЕТИКА И ИНФОРМАТИКА
УДК 612.087: 616.151.4-073.213: 612.13: 612.117.2
РАЗРАБОТКА МАТЕМАТИЧЕСКИХ МОДЕЛЕЙ ДЛЯ ОЦЕНКИ ФУНКЦИОНАЛЬНОГО СОСТОЯНИЯ СЕРДЕЧНО-СОСУДИСТОЙ СИСТЕМЫ
Эльшан Асиф оглы Рзаев1, Сакит Рауф оглы Расулов2*, Аббас Гейдар оглы Рзаев3
'Анкарский университет, Анкара, Турция;
2Азербайджанская государственная нефтяная академия, г. Баку, Азербайджан;
3Институт систем управления национальной академии наук Азербайджана, г. Баку, Азербайджан
Реферат DOI: 10.17750/KMJ2015-681
Цель. Разработка математических моделей гемодинамики (с учётом аномалии в гемореологии), позволяющих диагностировать функциональное состояние сосудов/сердечно-сосудистой системы.
Методы. При разработке математических моделей использованы методы эхокардиографии, математического моделирования, седиментации и законы реологии, механики и физиологии человека.
Результаты. Получены следующие математические модели: определения вязкости коллоидной дисперсной системы крови с учётом концентрации дисперсной фазы (клеток крови) и структурообразования крови; скорости стеснённого отстоя клеток крови в зависимости от стоксовой скорости отстоя и концентрации клеток; упруго-вязкостного свойства крови в виде зависимости скорости сдвига от вязкости крови, упругости (модуль Юнга) и напряжения сдвига; фильтрации крови в сосудах (модифицированная форма закона Дарси) с учётом напряжения и скорости сдвига разрушенной и неразрушенной структуры коллоидной системы крови. Показано, что при нарушении скорости кровотока происходят динамический отстой клеток крови и образование структуры (тромба), движение которой не подчиняется закону Ньютона. Введены новые показатели для диагностики функционального состояния сосудов и оценки степени тяжести сосудистой недостаточности.
Вывод. Разработанные модели гемореологии позволяют адекватно оценить диагностический показатель кровотока в сердечно-сосудистой системе организма человека.
Ключевые слова: вязкость, напряжение сдвига, скорость сдвига, кровь, плазма, капилляр.
DEVELOPING MATHEMATICAL MODELS FOR CARDIOVASCULAR SYSTEM FUNCTIONAL ASSESSMENTS
E.A. Rzaev1, S.R. Rasulov 2, A.G. Rzaev3
'Ankara University, Ankara, Turkey;
2Azerbaijan State Oil Academy, Baku, Azerbaijan;
3Institute of Control System of the Azerbaijan National Academy of Sciences, Baku, Azerbaijan
Aim. Development of mathematical models of circulation (considering anomaly in hemorheology) allowing to diagnose functional condition of vessels/cardiovascular system.
Methods. Echocardiography, mathematical modeling, sedimentation and rheology laws, human mechanics and physiology methods were used for developing mathematical models.
Results. The following mathematical models were obtained: for determination of colloid dispersive blood system viscosity, considering concentration of dispersive phase (blood cells) and blood structure formation; velocity of inconvenient blood cells sedimentation depending on flow velocity of sediment and cell concentration; parameters of blood elasticity and viscosity as a connection between the velocity change and blood viscosity, Youngs elasticity and change tension; blood filtration in vessels (modified form of Darcys law) considering tension and changes of destroyed and undestroyed colloid blood system structure velocity. It was shown that impairments of blood flow velocity leads to blood cells sedimentation and thrombus structure formation, which is not moving according to Newton's law. New indicators for diagnosing functional condition of vessels and estimating the severity of vascular insufficiency are introduced.
Conclusion. Developed hemorheologic models allow to adequately estimate human cardiovascular bloodflow.
Keywords: viscosity, pressure change, blood, serum, microvasculature.
Успех лечения при заболеваниях сердечно-сосудистой системы (ССС) зависит от правильной оценки функционального состояния последней, определяемой точностью измерения, учёта взаимосвязи и взаимного влияния важнейших показателей исследуемой системы и квалификации врача. К числу самых распространённых методов определения функционального состояния сердца и отдельных элементов ССС относятся электронная стетоскопия, электрокардиография (ЭКГ), эхокардиография (ЭхоКГ) [2, 4, 5] и метод измерения артериального давления. При электронной стетоскопии улавливаются сердечные тоны и шумы. При ЭКГ регистрируют электрические поля, образующиеся при работе сердца; конечным результатом становится графическое представление раз-
Адрес для переписки: rasulovsakit@gmail.com
ности потенциалов, возникающих в результате работы сердца. ЭхоКГ представляет собой один из важнейших технических методов определения состояния ССС, при котором посредством ультразвука получают информацию о форме, работе сердца и сосудов. Однако эти методы носят локальный характер и не могут определить влияние гемореологии и гемодинамики на состояние ССС.
Из вышеизложенного следует, что для повышения качества диагностики функционального состояния ССС необходимо разрабатывать математические модели, отражающие влияние изменения гемореологии, аномальной вязкости на гемодинамику и позволяющие качественно оценить сосудистую недостаточность.
Существует классическая гидродинамическая математическая модель Навье-Стокса [1], отражающая взаимосвязь между градиентом
давления, скоростью потока, плотностью и вязкостью жидкости:
р & = АР + цА2и + (1 ц + ц)АЛ™ (1),
где р, и — соответственно плотность, кг/м3 и линейная скорость потока крови (м/с); ц, цу — коэффициенты динамической и объёмной вязкости крови (Па-с); АР—градиент давлений (Па/м), V — скорость отстоя клеточных элементов.
Если к системе уравнений Навье-Стокса добавить уравнение неразрывности потока
( V = ^ V = = ^ V
•>11 ■'2 2 ■'пп
то математически движение вязкой жидкости будет полностью описано, так как уравнение будет включать значения всех факторов, влияющих на движение вязкой жидкости, где /. (1=1, п) — площадь поперечного сечения сосудов (артерий, артериол, капилляров, венул и вен, артериовеноз-ных анастомозов).
Однако ввиду сложности уравнения (1) его решение возможно только для некоторых частных случаев.
Существует также математическая модель Пу-азейля-Гагена в виде зависимости объёмной скорости движения крови (О) по сосудам от разности давлений в начале и конце сосуда (АР), радиуса (г) и длины сосуда (I) и вязкости жидкости (ц):
0= АР
^ 8 ц1
Это уравнение также не может быть использовано для оценки функционального состояния ССС, так как оно выведено для труб постоянного диаметра, а сосуды имеют разные диаметры (аорта — 25 мм, артерия — 4 мм, капилляр — 8 мкм, ве-нула — 20 мкм, вена — 5 мм, полая вена — 30 мм). Причём, несмотря на минимальный диаметр отдельного капилляра (8 мкм), их общий просвет приблизительно в 600-800 раз больше просвета аорты. Следовательно, расчётная скорость кровотока в капиллярах составляет около 0,06 см/с, поэтому основная проблема структурирования крови и образования тромба возникает именно в системе капилляров. Для устранения данного недостатка и достижения поставленной цели с успехом можно применить классический закон Дарси:
и= К Ега±Р (2),
где К — проницаемость системы капилляров (м2).
В последнее время появилась работа по модульному моделированию ССС человека [3], где используется модель артериального дерева, созданная на основе модели, описывающей течение крови в артериальном сосудистом русле человека. При этом с целью учёта эластичности стенок артерий применяют закон Гука:
м
R2„
(3),
где Е — модуль упругости Юнга; — толщина стенки артерии; R, Rg — соответственно радиусы артерий при текущем и нулевом давлениях; Р — текущее давление (напряжение) в артериях.
Постановка задачи.
Анализ вышеприведённых математических моделей, используемых для оценки функционального состояния сосудов, показывает, что существующие математические модели не учитывают влияния главного фактора — аномальной вязкости крови — на объёмную скорость последней. Следовательно, разработка математических моделей кровотока с учётом вязкости крови с целью диагностики функционального состояния сосудов является актуальной задачей.
Решение задачи.
I вариант. Любое описание циркуляции жидкостей внутри сосудов отличается по сравнению с технологическими жидкостями тем, что сосуды не являются ригидными (твёрдыми) трубками с постоянным поперечным сечением. Кроме того, при протекании крови в капиллярах происходит набухание эритроцитов, при котором они приобретают приблизительно сферическую форму. При этом математическую модель, описывающую гемореологию, можно представить комбинацией пружин и поршней.
Натяжение пружины (эритроцита) пропорционально её деформации, а сила, действующая на поршень, пропорциональна скорости деформации. Следовательно, пружины и поршни в модели представляют соответственно упругие и вязкие свойства крови.
Соединение основных элементов механической модели — пружины и поршня в параллельную цепь называется элементом Кельвина-Фойг-та (рис. 1), который описывается следующей формулой:
х=Еу+р^ или^+|у-^- = 0 (4),
где т — напряжение сдвига или сокращение оболочки эритроцита (Па); у — деформация; ёу/ёт — скорость деформации (1/с); Е — модуль упругости Юнга (Па); ц — модуль вязкости (Па-с).
Деформация у в таком теле (крови) под действием постоянной нагрузки т0 развивается во времени. Скорость её снижается, так как на упругий элемент Гука (пружины) приходится всё большее усилие. Когда скорость деформации уменьшается до нуля, деформация достигает максимального значения (см. рис. 1). При условии постоянного напряжения математическая модель Кельвина-Фойгта (1) принимает следующий вид:
Для решения этого дифференциального уравнения разделим переменные
<Ъ{ йг
Рис. 1. Модель вязкоупругого элемента Кельвина-Фойгта (а) и зависимость деформации при т=т0 (б) и при т=0 (в) от времени
0,5 (р
Рис. 2. Зависимость вязкости дисперсной системы от объёмной доли твёрдых сферических частиц размером: 1 — 0,1 мкм; 2 — 0,5 мкм; 3 — 1 мкм; 4 — 1,5 мкм
Окончательное решение приводит к следующему уравнению:
, Е ч
1 - ехр(--()
И
(5)
£
В работе [3] приведено множество эмпирических и полуэмпирических моделей, определяющих зависимость эффективной вязкости дисперсной системы (то есть крови) от напряжения сдвига и концентрации частиц. Следует отметить, что описанные в литературе реологические модели, в том числе и модели для вычисления вязкости, не являются приближениями физических законов, а представляют собой эмпирические или полуэмпирические приближения, описывающие кривые течения крови в определённом интервале скоростей сдвига. Исходя из этого, эффективную вязкость крови с объёмным содержанием твёрдых частиц клеток крови ф предлагают определить по эмпирической формуле, являющейся модифицированной формой уравнения Эйнштейна:
1+2,54 + —ф ехр( ———-4 1(Ф-Ф„,)2
(6),
ёмная доля частиц (эритроцитов, лейкоцитов и тромбоцитов); фт — концентрация частиц при их максимальной упаковке; ц — эффективная вязкость крови (Па-с).
Как видно из формулы (6) и рис. 2, эффективная вязкость (це) существенно зависит от концентрации клеток в крови. При определённой температуре (~36,6 °С) повышение концентрации клеток от 33 до 53% обусловливает утроенное значение эффективной вязкости.
Застой и ацидоз в области микроциркуляции придают крови вязко-эластические качества геля. Кровь как неньютоновская жидкость обладает вязкостью с многочисленными аномалиями, зависящей одновременно от концентрации клеток в крови, протеинов и жиров в плазме, эластичности эритроцитов, скопления всех клеточных элементов и образования структуры, водородного показателя (рН) крови, температуры и др.
Реологический характер плазменных макромолекул и клеток крови определяется состоянием равновесия между силой диффузии и тенденцией их стеснённого осаждения по закону [7]:
0=1,50,(1-9)
щ -(Зщ, + 2йе)(1 + ^)2 + 2(цф + + (Зцф+2цс)(1 + Л-)3
(7),
где т — коэффициент, который определяется экспериментально; в нашем случае т=2,2; ц — динамическая вязкость плазмы крови (Па^с); ф — объ-
где и^ — стоксовая скорость отстоя, определяемая по формуле: _ , , ,
и, =--—•
9ц
. гб'2У/з
где к = — ; р, р — соответственно плотность клеток и плазмы (кг/м ); /ис, ц — вязкость дисперсном
Рис. 3. Характерная кривая течения структурированной крови
среды плазмы и фазы клеток и макромолекул крови (Па-с); г — радиус частицы (м).
Таким образом, формулы (5), (6) и (7) могут быть использованы для оценки состояния процесса кровотока в капиллярах.
II вариант. Как отмечено выше, кровь с аномальной вязкостью по своим свойствам отличается от обычной, и её описание не подчиняется законам течения ньютоновских жидкостей.
Характерная кривая течения структурированной крови приведена на рис. 3. Кривая течения имеет два асимптотических состояния: ньютоновскую прямолинейную зависимость и зависимость, описываемую уравнением Бингама.
На кривой 2 участок ОА соответствует малым скоростям течения (явление ползучести), при которых системе наносится незначительный сдвиг, так как разрушения, связанные с течением, успевают тиксотропно восстанавливаться, то есть течение системы происходит без разрушения её структуры. Это так называемое явление ползучести течения крови протекает при максимальной вязкости кровяной системы.
Далее с ростом напряжения т или перепада давлений по длине сосудов (АР) начинается постепенное разрушение временных контактов, в том числе водородных связей, между элементами структуры (эритроцитами и плазмой) и образование других. В результате возникает динамическое равновесие, скорость деформации (течение) резко возрастает, и реологическая кривая выходит за участок ВС кривой 2, соответствующий пластической вязкости (ц*) кровяной системы. Экстраполяция кривой ВС по оси 0-т даёт значение тр, которое количественно характеризует сдвиговую прочность структуры и является предельным динамическим напряжением сдвига, соответствующим предельному напряжению разрушения структуры. Участок пластического течения ВС описывается уравнением Шведова-Бингама:
где т0 — напряжение сдвига, соответствующее началу разрушения структуры; ц* — пластическая вязкость системы; у — скорость сдвига структурированной крови.
Для крови аномальной вязкости закон фильтрации (протекания) Дарси отклоняется от клас-684
сической формы и может быть записан в нелинейной форме:
к дР
\х(Т,^ас1Р) дх
(9)
Из уравнения Ньютона получим:
с учётом (8)
Определив из выражения (8) =т-т0 окончательно получим выражение для эффективной вязкости крови в виде:
(10),
из которого следует, что с ростом т значение ц уменьшается и в пределе при т » т0 величина ц ^ ц* соответствует кровяной системе с полностью разрушенной структурой (участок СD) т=тр. Таким образом, вязкость структурированной системы в процессе течения под действием возрастающего напряжения сдвига или перепада давлений по длине сосудов изменяется от ц, отвечающего неразрушенной структуре, до ц^, характерного для полностью разрушенной структуры. Подставляя (10) в формулу (9), получим уравнение Дарси для структурированной кровяной системы:
где х — координата поперечного сечения капилляра.
При т >> т0 это выражение переходит в обычное уравнение Дарси [формулу (10)].
В случае аномальной кровяной системы для описания зависимости напряжения сдвига от градиента давлений может быть использована следующая формула:
1п^ = а [_Ша№_ т0 \igradP) 0
gradP =4^2,
где
% — пиковая скорость кровотока, определяе-
мая с помощью ЭхоКГ [6].
Очевидно, что показатель п в зависимости от температуры и свойств капилляров характеризует кривизну участка ВС (см. рис. 3).
Таким образом, скорость кровотока (и) и отстоя клеточных элементов (V), эффективная вязкость (¡1), напряжение (х) и скорость сдвига (у ) структурированной крови служат диагностическими параметрами состояния процесса кровообращения. Математические модели (5), (6), (7) и (11) позволяют более реально, чем существующие формулы, количественно определить значения этих параметров.
Для диагностики функционального состояния сосудов и ССС вводится новый показатель — подвижность крови в сосудах:
г= к (1 - V) = 2 и ц (1 т ) АР
и по нему определяют степень тяжести сосудистой недостаточности:
D„ - D„
С==
D
100% ,
где Q — объёмная скорость кровотока; l — общая длина сосудов; АР — перепад давления по длине сосудов; DH и DT — соответственно подвижность крови в нормальных (здоровых) и оцениваемых сосудах. При СТ >70% состояние сосудов считают тяжёлым.
ВЫВОД
Разработанные модели гемореологии позволяют адекватно оценить диагностический показатель кровотока в сердечно-сосудистой системе организма человека. Показано, что при нарушении скорости кровотока происходят динамический отстой клеток крови и образование структуры (тромба), в этом случае движение крови не подчиняется закону Ньютона. Введены новые показатели для диагностики функционального состояния кровеносных сосудов и оценки степени тяжести их недостаточности.
ЛИТЕРАТУРА
1. Алёхин М.Н., Сидоренко Б.А. Современные подходы к эхокардиографической оценке систолической функции сердца // Кардиология. — 2007. — Т. 47, №7. — С. 4-12. [Alekhin M.N., Sidorenko B.A. Contemporary approaches to echocardiographical assessment of systolic function of the heart. Kardiologiya. 2007; 47 (7): 4-12. (In Russ.)]
2. Вилкенсхоф У., Крук И. Справочник по эхокар-диографии, пер. с нем. — М.: Медицинская литература, 2008. — 229 с. [Wilkenshoff U., Kruck I. Handbuch der Echokardiografie. Georg Thieme Verlag KG. 2008: 308 p. Russ Ed.: Wilkenshoff U., Kruck I. Spravochnik po ekhokardiografii. Moscow: Meditsinskaya literatura. 2008: 229 p. (In Russ.)]
3. Киселёв И.Н., Семисалов Б.В., Бибердорф Э.А. и др. Модульное моделирование сердечно-сосудистой системы человека // Матем. биол. и биоин-форм. — 2012. — Т. 7, №2. — С. 703-736. [Kiselev I.N., Semisalov B.V., Biberdorf E.A. Modular modeling of the human cardiovascular system. Matematicheskaya biologiya i bioinformatika. 2012; 7 (2): 703-736. (In Russ.)]
4. Матвеенко В.Е, Кирсанов Е.А. Вязкость и структура дисперсных систем // Вестн. Моск. унта. Сер. 2: Химия. — 2011. — Т. 52, №4. — С. 243-276. [Matveenko V.N., Kirsanov E.A. The viscosity and structure of dispersed systems. Vestnik Moskovskogo universiteta. Seriya 2: Khimiya. 2011; 52 (4): 243-276. (In Russ.)]
5. Млюк В.Г., Млюк С.Э. Основные принципы гемодинамики и ультразвукового исследования сосудов. Руководство по ультразвуковой диагностике / Под ред. В В. Митькова. —. М.: Видар, 1997. — Т. 4. — С. 158-220. [Mlyuk V.G., Mlyuk S.E. Osnovnye printsipy gemodinamiki i ul'trazvukovogo issledovaniya sosudov. Rukovodstvo po ul'trazvukovoy diagnostike. (Main principles of hemocirculation and vessel ultrasonography. Guidelines on unltrasonography.) Ed. by V.V. Mit'kov. Moscow: Vidar. 1997: 158-220. (In Russ.)]
6. Райдинг Э. Эхокардиография. Практическое руководство. 3-е изд. Перевод с англ. — М.: МЕД-пресс-информ, 2013. — 280 c. [Ryding A. Essential Echocardiography. Elsevier Ltd. 2013: 272 p. Russ Ed.: Ryding A. Ekhokardiografiya. Prakticheskoe rukovodstvo. 3rd ed. Moscow: MEDpress-inform. 2013: 280 p. (In Russ.)]
7. Kelbaliev G.I., Rzaev A.G., Guseinov A.F., Kasumov A.A. Sedimentation of particles from a consentrated dispersed flow. — Plenum Publishing Corporation, 1992. — P. 1064-1067.