Научная статья на тему 'Эквивалентные электрические модели биологических объектов'

Эквивалентные электрические модели биологических объектов Текст научной статьи по специальности «Физика»

CC BY
1385
217
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ЭЛЕКТРИЧЕСКИЙ ИМПЕДАНС / РЕОГРАФИЯ / ЭЛЕКТРИЧЕСКАЯ МОДЕЛЬ / БИОЛОГИЧЕСКИЕ ТКАНИ / ЭКВИВАЛЕНТНЫЕ СХЕМЫ ЗАМЕЩЕНИЯ / ELECTRICAL IMPEDANCE / RHEOGRAPHY / ELECTRIC MODEL / BIOLOGICAL TISSUES / EQUIVALENT SUBSTITUTION CIRCUITS

Аннотация научной статьи по физике, автор научной работы — Зуев А. Л., Мишланов В. Ю., Судаков А. И., Шакиров Н. В., Фролов А. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Зуев А. Л., Мишланов В. Ю., Судаков А. И., Шакиров Н. В., Фролов А. В.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Known electrical equivalent circuits are discussed, and modifications of these circuits are proposed. Examples of the analysis and description of the impedance spectrum of real biological objects are provided. The amplitude-frequency and phase-frequency characteristics of various circuits are analyzed. The frequency characteristic of the phase-shift angle is found to be rather informative. For three-element models consisting of two resistors and one capacitor, the necessary identity conditions are determined.

Текст научной работы на тему «Эквивалентные электрические модели биологических объектов»

УДК 531/534: [57+61]

ЭКВИВАЛЕНТНЫЕ ЭЛЕКТРИЧЕСКИЕ МОДЕЛИ БИОЛОГИЧЕСКИХ ОБЪЕКТОВ

А.Л. Зуев1,2, В.Ю. Мишланов3, А.И. Судаков1, Н.В. Шакиров1, А.В. Фролов4

1 Институт механики сплошных сред Уральского отделения Российской академии наук, Россия, 614013, Пермь, ул. Академика Королева, 1, e-mail: [email protected]

2 Пермский национальный исследовательский политехнический университет, Россия, 614990, Пермь, Комсомольский проспект, 29

3 Пермская государственная медицинская академия имени академика Е.А. Вагнера, Россия, 614990, Пермь, ул. Петропавловская, 26

4 Республиканский научно-практический центр «Кардиология», Республика Беларусь, 220036, Минск, ул. Р. Люксембург, 110

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

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

Введение

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

К настоящему времени разработано большое количество приборов и методов определения удельного сопротивления, проводимости и диэлектрической проницаемости. Было установлено, что величина биоэлектрического импеданса или адмиттанса существенно зависит от состояния тканей и кровообращения исследуемого участка тела человека. В [6, 8-10, 12] детально обсуждаются методы, приборы и схемы

© Зуев А.Л., Мишланов В.Ю., Судаков А.И., Шакиров Н.В., Фролов А.В., 2012

Зуев Андрей Леонидович, д.ф.-м.н., старший научный сотрудник ИМСС УрО РАН, Пермь

Мишланов Виталий Юрьевич, д.м.н., заведующий кафедрой пропедевтики внутренних болезней ПГМА,

Пермь

Судаков Андрей Иванович, старший научный сотрудник ИМСС УрО РАН, Пермь Шакиров Нагим Вагизович, старший научный сотрудник ИМСС УрО РАН, Пермь Фролов Александр Владимирович, д.б.н., заведующий лабораторией медицинских информационных технологий, Минск

09806267

100000 1000000

Частота, Г ц

а

Частота, Г ц б

Рис. 1. Спектральные характеристики сегмента тела человека (отведение «правая рука - правая нога») [14]

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

Биоэлектрический импеданс представляет собой соотношение между напряжением и током, и поэтому для его моделирования используются электротехнические элементы и эквивалентные схемы замещения, под которыми понимаются модельные электрические схемы, составленные из резисторов, конденсаторов и более сложных структурных элементов. Построение эквивалентной схемы исследуемого объекта является важным этапом изучения данного физического явления. Модель должна не только отражать основные свойства биообъекта в исследуемой частотной и температурной областях, но и предсказывать поведение объекта в более широком диапазоне частот и температур. Правильность интерпретации результатов электрических измерений усложняется тем, что эквивалентные схемы замещения, как правило, неизвестны, и поэтому исследуемый объект с электродами и измерительной аппаратурой представляет собой «черный ящик». Достоверность реографических формул, устанавливающих количественную связь между характеристиками реограммы и гемодинамическими параметрами, существенно зависит от соответствия электрофизической модели реальному объекту [3, 4, 12, 13]. Примеры анализа спектров импеданса для некоторых биологических сред, электролитических материалов и композитов приводятся в [1, 6, 7, 9].

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

Эквивалентные схемы замещения и их диаграммы

Полное электрическое сопротивление объекта можно записать в виде комплексной величины: Z = Z' — iZ”, где Z” - действительная (активная), Z” - мнимая (реактивная) составляющие импеданса. Наглядным средством отображения частотных свойств биоимпеданса могут служить дисперсионная зависимость полного сопротивления и угла фазового сдвига ф = агС£ (2" /Z') тока и напряжения от

частоты ю (рис. 1). Эта пара графиков называется диаграммой Боде.

Другим наглядным представлением импеданса является векторная диаграмма комплексного сопротивления (годограф, график Найквиста), показанная на рис. 2. Ромбиками обозначены данные экспериментальных измерений, проведенных на сегменте человека (отведение «запястье правой руки - щиколотка правой ноги»). Результаты, представленные на рис. 2, а, взяты из работы [14], а на рис. 2, б - из работы [5].

ей

о

X

ей

и

С

Л

н

ей

Г

К

X

45

40

35

30

25

20

15

10

5

200

225

250

275

300

325

Действительная часть импеданса, Ом

о

13

э

с

л

н

ей

Г

£

Е

К

X

80 -70 60 50 -40 -30 20 10 0

330

380

430

480

а

Действительная часть импеданса, Ом б

Рис. 2. Диаграммы Найквиста для сегмента тела человека (отведение «правая рука - правая нога»)

с

*

&

До

о.

я-,

б

Рис. 3. Трехэлементные модели

Использование биоэлектрического импеданса для оценки состава тела значительно возросло в последние годы. Для качественного и количественного моделирования электрических свойств отдельных частей человеческого тела широко используется трехэлементная модель (рис. 3, а), представляющая собой две параллельно включенные цепочки, одна из которых состоит из последовательно соединенных сопротивления Я1 и конденсатора Сі, а другая - только из активного сопротивления Я0. При низких частотах величина импеданса (7) определяется сопротивлением межклеточной жидкости Я0 = (7(ш ^ 0)). С ростом частоты

переменного тока реактивное сопротивление емкости (С1) уменьшается и все большая часть тока проходит внутри клеток. Используя правила расчета импедансов параллельного и последовательного соединения активного и реактивного элементов, для высоких частот находим Я1 = Я 7(ш ^ го)/(Я — 7(ш ^ ^)).

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

Я + (^0 + Я1)^0ю2 — іЯ^ С1Ю

7 =

1 + (Я0 + Я1)2 С12ю2

ф = —

Я0 С ш

1 + (Я0 + Я1) Я1С12 ш

2_2 •

(1)

(2)

Из формулы (2) следует, что фазовый угол при некоторой частоте ш = ш* достигает экстремума. Из условия экстремума фазового угла ф находим частоту, при которой отставание напряжения от тока максимальное:

Я - 7(<»)

*

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ш =

(3)

0

а

Выделяя в формуле (1) активную (2') и реактивную (2") составляющие импеданса, а затем исключая частоту ш, получаем в комплексной плоскости (2', 2") уравнение окружности:

Z -

2

+

( z ")2 =

(4)

Уравнение (4) описывает полуокружность с центром, лежащим на оси 2' в точке, отстоящей от начала координат на расстоянии (Я0 + 2(да))/2, и радиусом,

равным (Я0 -2(да))/2. Если в эквивалентной схеме, изображенной на рис. 3, б,

сопротивления Я2 и Яда выбрать равными: Я2 = АЯ = (2(ш ^ 0) - Яда),

Яда = (2 (ш ^ да)), то годографом этой схемы будет полуокружность, также

описываемая уравнением (4). Экстремум тангенса фазового угла для данной модели достигается при частоте

*

ш =

C2(Ro - Z(даЖ/ZK)'

(5)

Из сравнения формул (3) и (5) следует необходимое условие тождественности трехэлементных моделей, приведенных на рис. 3:

Cl

1 -

Rl

Ro

(б)

Таким образом, зная экстремальную частоту ш *, начальное Z(0) и конечное Z(да) значения импеданса, можно легко найти параметры (R0, R, R1, R2, C1, C2) трехэлементных эквивалентных схем замещения. Диаграммы Боде и Найквиста этих трехпараметрических моделей совпадут, если емкости выбрать удовлетворяющими соотношению (6).

Диаграммы Найквиста для многих биообъектов представляют собой «сдавленную полуокружность» (окружность с центром, лежащим ниже оси Z ') [6, 14] или результат перекрытия нескольких полуокружностей [1, 10]. Смещение центра окружности с горизонтальной оси легко достигается, если в эквивалентных схемах, изображенных на рис. 3, конденсаторы заменить на элемент постоянной фазы «Constant Phase Element - CPE» (рис. 4), импеданс которого задается формулой [11]

Z = A (ішт)

(7)

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

Я<с

7сре

г>

А Я

а

б

Рис. 4. Эквивалентные модели Коула-Коула

2

а

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

На примере модели Коула-Коула, представленной на рис. 4, б,

продемонстрируем схему построения диаграмм Боде и Найквиста и результаты сравнения экспериментальных и расчетных данных.

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

(R - я,)

Z = я. +■

1 + (raT)a cos |aK|-і (шт)

\a . I К

sin | a— 2

1 + 2(®T)a cos| aK-j + (QT)2a

(8)

tg ф =

(T®)a sin |a2 J( Ro- я»)

R +(T®)a cosI aK |(Ro + R,) + R»(T®)2a

(9)

(ш* т)

2 a

ш Т) =

Rl

R,

(10)

На рис. 1 представлены результаты сравнения экспериментальных [14] (отмечены ромбиками) и расчетных данных (показаны квадратиками), выполненных по формулам (8)-(10) при следующих параметрах: R0 = 305,3 Ом; R<X) = 211,4 Ом; а = 0,89; т = 4,810-6 с.

На этом же рисунке треугольниками показана аппроксимация, выполненная с использованием модели, изображенной на рис. 6, а. При этом значения параметров модели принимались равными: R1 = 973 Ом; R2 = 227 Ом; R3 = 307 Ом; С1 = 2,217-10-9 Ф; С2 = 410-9 Ф. Так как последняя модель содержит два конденсатора, то можно ввести в рассмотрение два времени релаксации (см. приложение) и найти их отношение:

^ = RC1 = 2,110-6 с, ^2 = R2C2 = 9 105 с, XJХ2 = 2,3.

(11)

Заметим, что для данного биологического объекта времена релаксации оказались одного порядка и сопоставимыми с постоянной времени т, входящей в модель Коула-Коула. Для того чтобы построить годограф Найквиста, выделим в формуле (8) действительную и мнимую части:

(Ro - R,)

Z' - R =

1 + (шт)** cos [aK

1 + 2 (шт^ cos |aK j + (шт)2a

(12)

Z =-

(шт)a sin[ aK

1 + 2(шт^ cos|aK | + (шт)2

(13)

и найдем их отношение

Z' - R>. t [ к|

-----------+ ctg| a— I =

Z" I 2 і

1

(шт)a sin[aK

(14)

Подставляя угловую частоту, выраженную из уравнения (14), в соотношение (13), после несложных преобразований получим уравнение окружности с центром, лежащим ниже действительной оси:

) Т --Т0,8Ы]

^Z , - Ro + Я,

\

2

AR

1

2

К

(15)

sin a

2

Если в качестве электрической схемы выбрать модель, представленную на рис. 4, а, а сопротивление Я1 принять равным Я0Я,/(Я0 - Я,), то ее импедансные

характеристики (модуль импеданса, угол сдвига фаз и экстремальная частота) совпадут с результатами вычислений, выполненных по формулам (8)-(10). Заметим, что уравнение окружности (4) и условие экстремальности фазового угла (5) непосредственно следуют из соотношений (15) и (10) при a = 1.

Импеданс многих биообъектов не описывается трех- или четырехэлементными схемами, рассмотренными выше; их годограф Найквиста отличается от окружности, а диаграмма фазового угла характеризуется несколькими экстремумами. Это хорошо иллюстрируют экспериментальные данные реограммы легкого [10], приведенные на рис. 5. Следует отметить два четко выраженных экстремума на диаграмме фазового угла (см. рис. 5, а) и пересечения окружностей, имеющих различные диаметры и места расположения их центров (см. рис. 5, б).

В работе [10] для расчета импеданса использовалась пятиэлементная модель, приведенная на рис. 6, б. Значения параметров, входящих в эту модель, принимались равными: Я1 = 49,7 Ом; Я2 = 97,9 Ом; Я» = 49,8 Ом; C1 = 1,44 10-7 Ф; C2 = 1,47 10-9 Ф. Вычислим для этой модели времена релаксации и их отношение:

X1 = R1C1 = 71,57-10-7 с, X2 = Я2 С2 = 1,4410-7 с, V X2 = 50.

(16)

Хорошая аппроксимация этих экспериментальных данных также достигается, если использовать модель, приведенную на рис. 6, а (Я1 = 74,6 Ом; Я2 = 609 Ом; Я3 = 194 Ом; С\ = 6,610-9 Ф; С2 = 8,810-8 Ф). Для этой модели соответствующие времена релаксации, обозначенные как L1 и L2, оказались равными

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

L1 = 706,6 10-7 с, L2 = 17,73 10-7 с,

L1 / L2 = 40.

(17)

10

100 1000 10000 100000 1000000

а

Действительная часть импеданса, Ом б

л,

-I--------------И

С,

ь

&

С2

I

я

а б

Рис. 6. Эквивалентные пятиэлементные модели

3000 3 2800 ° 2600

2400

2200

2000

К 1800

1600

1400

10

100 1000 10000 100000 1000000

10

100 1000 Частота, Г ц а

10000 100000

о

3

Он

и

сб

и

к

«

о

ч

о

и

-15-

-25 х

Частота, Гц

б

Рис. 7. Частотные зависимости реограммы человеческой крови

А -I------------------У

С1

о

с2

—00-

а

а б

Рис. 8. Эквивалентные шестиэлементные и семиэлементные модели

Частотные зависимости биологических растворов и крови [3] отличаются от вышерассмотренных тем, что дисперсия импеданса наблюдается при значительно меньших частотах (рис. 7).

Такое поведение биологических тканей можно достаточно хорошо описать эквивалентными моделями, приведенными на рис. 8. Эти эквивалентные схемы обобщают соответствующие модели, приведенные на рис. 6. Как частный случай из модели, приведенной на рис. 8, б, при стремлении Я4 к бесконечности вытекает шестиэлементная эквивалентная схема, предлагаемая в работе [2] для описания импеданса крови.

Результаты вычислений, выполненных с использованием шестиэлементной модели

(значения параметров: Я\ = 4148,7 Ом; Я2 = 8041,7 Ом; Яз = 3972,2 Ом; С = 1,18 10 6 Ф;

_8

С2 = 7,57-10“ Ф; С3 = 1 Ф), показаны на рис. 7 треугольниками. Аппроксимация реограммы человеческой крови семипараметрической моделью (Я\ = 170 Ом; Я2 = 170 Ом; Я3 = 2840 Ом; Я4 = 3610 Ом; С1 = 3,110_6 Ф; С2 = 3,1-10_6 Ф; С3 = 1,46-10-6 Ф) показана квадратиками.

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

Обсуждение результатов

Если диаграмма Найквиста исследуемой биоткани представляет собой окружность (или кривую, близкую к ней) с центром, находящимся на действительной оси, то такую среду можно удовлетворительно моделировать эквивалентной схемой замещения, состоящей из двух сопротивлений и одного конденсатора (рис. 3). При этом значения сопротивлений находятся из условий R0 = (Z (ш ^ 0)) и R^ = (Z (ш ^ го)) ,

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

Если центр окружности на диаграмме Найквиста находится ниже действительной оси, то спектральные характеристики объекта достаточно хорошо описываются моделью Коула-Коула или пятиэлементной электрической схемой. Из рис. 1 видно, что кривые Z = Z(ш) и ф = ф(ш) удовлетворительно описываются как моделью Коула-Коула, так и пятипараметрической электрической схемой, представленной на рис. 6, а. Следует отметить, что в модель Коула-Коула не входит в явном виде элемент, моделирующий емкостные свойства среды. Однако если задать а = 1, то элемент CPE вырождается в конденсатор с емкостью C, а параметр т становится равным произведению AR C. Следовательно, в данном случае постоянная времени совпадает с временем релаксации.

Поскольку модели, представленные на рис. 6, содержат по два конденсатора, то можно ввести в рассмотрение несколько времен релаксации. В том случае, когда отношение времен релаксации меньше 10, на частотной зависимости фазового угла имеется только один минимум (см. формулу (11) и рис. 1, б).

Если времена релаксации электрической схемы отличаются в десятки раз, то кривая зависимости угла сдвига фаз от частоты имеет два минимума, а диаграмма Найквиста состоит из двух пересекающихся окружностей (см. формулу (17) и рис. 5).

Из рис. 5 видно, что обе модели предсказывают два экстремума на зависимости угла сдвига фаз от частоты. Из полученных выше значений времен релаксации можно предположить, что необходимым условием нескольких экстремумов на частотной зависимости угла сдвига фаз является отличие времен релаксации в десятки раз (см. формулы (16) и (17)).

Из вышеприведенных рисунков (см. рис. 1, б и рис. 5, а) видно, что спектральные характеристики сегментов тела человека слабо изменяются вплоть до частот порядка 1000 Г ц.

Выводы

Спектральные характеристики отдельных сегментов тела человека очень слабо изменяются вплоть до частот порядка 1000 Гц. Дисперсия же импеданса крови человека наблюдается при значительно меньших частотах.

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

Сформулировано условие тождественности трехпараметрических электрических схем: одна из них представляет собой параллельную цепь R2C2 с добавочным резистором R«, а вторая - последовательную цепь R1C1, зашунтированную резистором R0.

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

Приложение

На примере параллельной схемы замещения вычислим импеданс.

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

U + CU = J, Я

(I)

где и_ напряжение; J _ ток, J = Q , Q _ электрический заряд.

Найдем решения уравнения (I) для двух конкретных изменений силы тока. В первом случае задается постоянный электрический заряд емкости Q = Q0 (т.е. J = 0) и находится изменение напряжения во времени

U=texP (тJ

(II)

где X - параметр, характеризующий скорость приближения напряжения к равновесному состоянию, X = R C. Уравнение (II) описывает релаксационный процесс, и поэтому параметр X часто называют временем релаксации.

Во втором случае в правую часть уравнения (I) подставляется периодическая функция J = J0 exp (mt) и находится периодическое решение

R (1 - mRC) . .

U = — -------J0 exp(mt).

1 + (шЯС)

Разделив напряжение на ток, находим изменение импеданса 2 в зависимости от

частоты:

Z=

Я (1 - ішХ)

1 + (шХ)2

(III)

Параметр X имеет размерность времени. Размерность времени имеет и параметр т, входящий в выражение для импеданса элемента постоянной фазы (см. формулу (7)).

Покажем, что в случае а = 1 эти параметры совпадают. Если в модели, приведенной на рис. 4, б, сопротивление Яда выбрать равным нулю, то из формулы (8) находим

Я (1 -/шт)

2 = , , .Л (^)

1 + (шт)

Из сравнения формул (III) и (IV) видно, что сопротивление Я0 = Я, а постоянная времени т равна произведению ЯС (времени релаксации). Из формул (III) и (IV)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

, - 1

следует, что параметры X и т определяют характеристическую частоту ш =

ЯС

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

Благодарности

Работа поддержана Российским фондом фундаментальных исследований (проекты РФФИ-Урал № 11-08-96005 и № 11-04-96017).

Список литературы

1. Буянова Е.С., Емельянова Ю.В. Импедансная спектроскопия электролитических материалов: учебное пособие. - Екатеринбург: Из-во УрГУ, 2008. - 70 с.

2. Ефремов А.В., Ибрагимов Р.Р., Манвелидзе Р.А. и др. Устройство для измерения активной и емкостной составляющих импеданса биологических тканей: патент РФ № 2196504, 2003.

3. Зуев А.Л., Мишланов В.Ю., Судаков А.И., Шакиров Н.В. Экспериментальное моделирование реографической диагностики биологических жидкостей // Российский журнал биомеханики. -2010. - Т. 14, № 3. - C. 68-78.

4. Кубарев А.М., Борисов В.И. Пульсация крови в артериальной системе и ее влияние на электрическое сопротивление тела // Нижегородский медицинский журнал. - 2008. - № 4. - С. 35-41.

5. Можаев В.А. Векторные диаграммы комплексного сопротивления и проводимости для оценки активной клеточной массы биообъекта // Материалы X научно-практической конф. «Диагностика и лечение нарушений регуляции сердечно-сосудистой системы». - М., 2008. - С. 92-98.

6. Николаев Д.В., Смирнов А.В., Бобринская И.Г., Руднев С.Г. Биоимпедансный анализ состава тела. -М.: Наука, 2009. - 404 с.

7. Поклонский Н.А., Горбачук Н.И. Основы импедансной спектроскопии композитов. - Минск: Изд-во БГУ, 2005. - 130 с.

8. Цветков А.А. Биоимпедансные методы контроля системной гемодинамики. - М.: Слово, 2010. -330 с.

9. Шван Х.П., Фостер К.Р. Воздействие высокочастотных полей на биологические системы: Электрические свойства и биофизические механизмы // ТИИЭР. - 1980. - Т. 68, № 1. - С. 121-132.

10. Ackmann J.J. Complex bioelectric impedance measurement system for the frequency range from 5 Hz to

1 MHz // Annals of Biomedical Engineering. - 1993. - Vol. 21 - P. 135-146.

11. Cole K.S., Cole R.H. Dispersion and absorption in dielectrics. I. Alternating current characteristics // J. Chem. Phys. - 1941. - Vol. 9. - P. 341-351.

12. Dai T., Adler А. Blood impedance characterization from pulsatile measurements // Proceedings of the Canadian Conference on Electrical and Computer Engineering, CCECE 2006, May 7-10, 2006, Ottawa Congress Centre. - Ottawa, 2006. - P. 983-986.

13. Macdonald J.R. Impedance spectroscopy // Annals of Biomedical Engineering. - 1992. - Vol. 20. -P. 289-305.

14. Nordbotten B.J., Martinsen O.G., Grimnes S. Methods for calculating phase angle from measured whole body bioimpedance modulus // J. Phys.: Conf. Ser. - 2010. - Vol. 224. - P. 253-268.

EQUIVALENT ELECTRICAL MODELS OF BIOLOGICAL OBJECTS

A.L. Zuev, V.Yu. Mishlanov, A.I. Sudakov, N.V. Shakirov (Perm, Russia),

A.V. Frolov (Minsk, Belarus)

Known electrical equivalent circuits are discussed, and modifications of these circuits are proposed. Examples of the analysis and description of the impedance spectrum of real biological objects are provided. The amplitude-frequency and phase-frequency characteristics of various circuits are analyzed. The frequency characteristic of the phase-shift angle is found to be rather informative. For three-element models consisting of two resistors and one capacitor, the necessary identity conditions are determined.

Key words: electrical impedance, rheography, electric model, biological tissues, equivalent substitution circuits.

Получено 02 февраля 2012

i Надоели баннеры? Вы всегда можете отключить рекламу.