Особенности биоимпедансных измерений на высоких частотах

можность оценить эффективность системы самой по себе и относительно однотипных систем.

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

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Кузнєцов Д.М. Інформаційно-вимірювальна система градуювання термоане-мометричних вимірювальних перетворювачів швидкості газа. Автореф. дисертації на здобуття наукового ступеня кандидата технічних наук. Донецький національний технічний університет. Донецьк, -2005.

2. Новопашенный Г.Н. Информационно-измерительные системы: Учебн. пособие для специальности „Информационно-измерительная техника” вузов. -М.: Высш. школа, 1977. -208 с.

3. Рабинович В.И., Цапенко М.И. Информационные характеристики средств измерения и контроля. - М.: Энергия, 1968. - 96 с.

4. Кузьмин И.В., Кедрус В.А. Основы теории информации и кодирования. - К.: Вища шк., 1997.

5. Касаткин А.С., Кузьмин И.В. Оценка эффективности автоматизированных систем контроля. - М.: Энергия, 1967. - 80 с.

6. Солодов А.В. Теория информации и ее применение к задачам измерения и контроля. М.: Наука, 1967. - 432 с.

7. Кузьмин И.В. Оценка эффективности и оптимизация автоматических систем контроля и управления. - М.: Советское радио, 1971. - 296 с.

8. Тарасенко Ф.П. Введение в курс теории информации. -Томск: Изд-во Томского университета, 1963. - 240 с.

Е.В. Мирошниченко, Е.С. Семенистая ОСОБЕННОСТИ БИОИМПЕДАНСНЫХ ИЗМЕРЕНИЙ НА ВЫСОКИХ ЧАСТОТАХ

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

• двухэлектродный метод;

• трехэлектродный метод;

• четырехэлектродный метод;

• пятиэлектродный метод.

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

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

Рис. 1

Этот метод корректно работает, если импеданс электрических контактов (Хс+ и Zc-) намного ниже, чем импеданс исследуемого тела (ИТ). На самом деле измеренный импеданс Z - это сумма ZИТ , ZC + и ZC_ . Контактный импеданс может быть больше чем импеданс ИТ, особенно на низких частотах, а также нестабилен и непредсказуем для математической коррекции. Таким образом, повторяемость и воспроизводимость измерений получается крайне низкой.

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

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

Четырехэлектродный метод, известен также как тетраполярный или метод Кельвина, используемый для измерения сопротивления материалов начиная с 1800 года. Он использует пару электродов для инжектирования тока в образец и другую пару электродов для измерения падения напряжения (рис 2). В принципе, так как ток не протекает через вольтметр (Zvmetr ), то инжектируемый ток полностью протекает через образец и падение напряжения на образце то же, которое регистрирует измеритель. Поэтому контактный импеданс никак не влияет на измерение.

Однако известно, что контактный импеданс влияет на измерение преимущественно на низких частотах, где его амплитуда может быть очень большой. Это происходит из-за того, что импеданс электрода приводит к появлению синфазного напряжения, которое не может быть полностью подавлено входным дифференциальным усилителем вольтметра. Например, предположим, что «земля» электрической цепи на рис. 2, соединена с электродом, принимающим ток (е4), тогда конечный импеданс этого электрода приведет к появлению падения напряжения (I  1е4), которое проявляется в виде синфазного напряжения на входе вольтметра. К сожалению, ситуация намного хуже так как внутренний импеданс электродов и паразитный импеданс вольтметра приведут к уменьшению коэффициента подавления синфазной помехи. Часть синфазного напряжения, присутствующего на образце, попав на дифференциальные входы усилителя, преобразуется в дифференциальное напряжение (рис. 3). Этот феномен был проанализирован [1] и было показано, что необходимо обеспечить входной импеданс максимально большим и сбалансированным.

Влияние паразитных емкостей на высоких частотах, показано на рис. 4. Некоторые ошибки возникают из-за следующих паразитных импе-дансов: ток, протекающий через образец уменьшается (шунтирование паразитным импедансом Z14), дифференциальное усиление уменьшается (шунтирование паразитным импедансом Z24), часть тока инжектируется на

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

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

Результат представленный на рис. 5. соответствует моделированию измерения биоимпеданса проводящего силикона, помещенного в солевой раствор (0,9 %). Параметры, использованные для моделирования совпадают с естественными, и, как можно увидеть, результаты моделирования достаточно близки к реальным измерениям. Несоответствие (30 %) между внутренними электродами должно быть хорошо видно на низких частотах, однако на этих частотах никаких ошибок не наблюдается.

с

Г%2

І-Е+01 1.Е*СВ »-Е+03 1.Е+04 !£+« 1.Е*0в

Іг^иепсу

2 ц

N 4

'"''ч

V

\

1 і

1.Е+01 1.Е+02 1.Е+03 1 Е-Л4 1.Е+4» 1Е+0Я

9ге

Рис. 5

Результаты моделирования, приведенные на рис. 6, соответствуют измерениям на гипотетической живой ткани (чистое удельное сопротивление тканей, р = 7000  см ). Ошибка на высоких частотах больше чем в пре-

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

frequency

N

\

\

А

1.Е*С1 1.E+Q2 1.Е+СЗ 1x404 1-Е-К VE+M

frequency

Рис. 6

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

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

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

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

• инжектируемая электрическая энергия (в случае большого тока) может вызвать различные виды повреждения живых тканей.

В живых тканях, соотношение напряжение/ток линейно до тех пор, пока прикладываемое электрическое поле не превышает порога нелиней -

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

Так как локальные электрические поля пропорциональны плотности тока (Е = J  р ), то порог нелинейности для данной ткани обычно вычисляется как порог плотности тока. Человеческая кожа может считаться линейным объектом при плотностях тока меньших, чем 100мкА/см2 на частоте 100Гц.

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

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

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

Пятиэлектродный метод

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

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

Рис. 7 показывает расположение электродов, требуемую коммутацию для использования этого метода. В период времени ^ электрод 1 переключается для инжекции тока. Это подразумевает то, что приложенный ток

будет протекать от электрода 1 до электрода 4. 2тг - результат этого измерения, а 2т2 - результат измерения, когда подключается электрод 5 и ток протекает от этого электрода к электроду 4.

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

Т *

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

1 л I

*

— = — = к а .

2 V

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

В данном примере можно определить две различные константы: к1 для фрагмента, составленного электродами 1, 2, 3 и 4 и к2 для фрагмента составленного электродами 5, 2, 3 и 4. Взаимосвязь между этими двумя значениями описывается следующим выражением:

т=-2-=21-,

к, 22

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

Когда выполнено четырехэлектродное измерение и синфазное напряжение, вызванное 2Е4 рассматривается как единственный источник ошибки, то измеряемый импеданс 2т будет:

V+ 12Е4

=-СМЩ= V +—^Е4— = 2 +—^Е4—,

т I I СМЯЯТ СМЯЯТ

где СМЯЯ - эффективный коэффициент понижения синфазной помехи, определенный в [6]. Если два фрагмента одновременно использующие измеряемый электрод 2Е4 используются при измерении одного и того же образца (одинаковой проводимости), тогда:

2 _ 2Е4

21 = 1 =-----СМЩ- ^ СМШ-1 = 2т1 _12т2 ;

22 2______2е4 Т 2е4

т2 СМШТ

7=7 2Е4 = 7 2т1 _1  2т2 7 ^7 = 1

21 2т1 СМЖТ 2т1 2Е4

Таким образом, имеется возможность скомпенсировать влияние синфазного напряжения, вызванного 2Е4 . Необходимым условием являет-

ся то, что в обоих случаях измерения импеданса CMRRT и ZE4 сохранят свои значения неизменными.

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

Результаты, приведенные на рис. 8, соответствуют опыту, в котором использовался раствора №С1 0.9 % с проводимостью 70£- ^ и они не имеют каких-либо заметных всплесков на частоте до 1 МГц. Как можно видеть, пятиэлектродный метод дает значительно лучшие результаты, особенно для мнимой части комплексного импеданса. Однако он не устранил ошибки на очень низкой частоте.

Известно, что импеданс электродов зависит от многих факторов и один из них это инжектируемый ток. Это значит, что импеданс электродов может быть нелинейным [11,12], поэтому можно считать, что разность плотности тока между измерениями возникает из-за отличия ZE4.

Ьт

Рис. 8

В случае, когда 7Е4 различен в двух измерениях (7Е41, 7Е4_2).

7 _ 7Е4_1

71 _ т1 СМШТ

-1 _я_ 72

7т2

-Е4 2

7і _Я-

СМШ-1 _

СМШТ

Я' 7т2 _ 7т1

Я- 7Е4_2 _ 7Е4_1 7т1•7Е4 2 _ 7т2-7Е4 1

Я-7

Е4 2

_7

Е 4 _ 1

что эквивалентно уравнению 5 в случае если 7і

Е4_1

Е4_2

чЕ4*

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

„ - 7т1 -Р- 7Е4_1 _ 7т2 - 7Е4 1 , 7т1 •

71 _ Я •-----------------------------------_ Я •

Р_7

т2

Я-Р- 7Е4_1 _ 7Е4_1 Я

Яр_ 1

71 _ 7С1 +е _ е_ Я

(7т1 7т2 ) + е ;

Я_ 1

^ 7т1 'Р_ 7т2 7т1 _ 7т2 ^

Яр_ 1

Я_ 1

8=1-

Если принять, что ошибка для четырехэлектродного метода:

ZE4

84 электр » ситТ ’

тогда ошибка пятиэлектродного метода:

8 =

1- 1

(Ь-1)- 84электр .

Из формулы 15 следует, что импеданс электрода 4 должен быть постоянным (Д — 1) и 1 должна быть как можно больше. В противном случае можно получить ошибки большие, чем при использовании более простого четырехэлектродного метода. К сожалению, необходимые условия противоположны друг другу; большое значение 1 подразумевает различное распределение тока, а это ведет к большому значению Д.

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

Было проведено измерение силикона при помощи пятиэлектродного метода (рис. 9) как и в прошлом эксперименте. Первый электрод, используемый для инжекции тока 1+1 и второй 1+2. Эта конфигурация привела к 1 ~ 1.5.

Как и в прошлом эксперименте, результат (рис. 9) показывает, что пятиэлектродный метод снижает ошибку на низкой частоте, но не способен устранить ее полностью.

I 20

II

г о

I

=- -20

N

-40

-40

Ттецчепсу

1000 йгечиепсу

Рис. 9

Далее измерялся импеданс электрода 4 (электрод, принимающий ток I) относительно инжектируемого электрода (1+1 или 1+2). Эти значения импеданса не совпадают со значениями из предыдущего измерения, но их разность совпадает с учетом величины вычисленной ранее ошибки.

Frequency

Рис. 10

Несмотря на то, что разность между обоими измерениями мала (рис. 10), растущая ошибка может достигать значения 0,1 и это означает, что возникнет значительная ошибка на низкой частоте. Как видно из рис. 11, абсолютное значение полученной ошибки хорошо совпадает с ее ожидаемыми значениями. Следовательно, этот эксперимент подтверждает, что изменение импеданса электрода служит причиной неспособности пятиэлектродного метода устранять ошибку на низкой частоте.

Frequency

Рис. 11

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

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

Другие недостатки пятиэлектродного метода были обнаружены при анализе эффекта добавления случайной ошибки, такой, например, как шум между двумя измерениями:

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

• Четырехэлектродный метод является наилучшим для мягких тканей, однако он наиболее чувствительный к ошибкам.

• Полоса частот для измерений импеданса определяется источником ошибок.

• Ошибки измерений зависят от уровня инжектируемого тока во время измерения.

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Pethig R. Dielectric Properties of Body Tissues. Clin. Phys. Physiol. Meas. 8 (suppl A). 1987. Р. 5 - 12.

2. Brown B.H. Tissue Impedance Methods. Prog. Med. Enviro. Phys. 1983. Р. 85 -110.

3. Barber D.C. and Brown B.H. Applied Potential Tomography (Review Article). J. Phys. E:Sci. Instrum 17.1984. Р. 723 - 733.

Выводы.

4. Plonsey R. (1984) Quantitive Formulations of Electrophysiological Sources of Potential Fields in Volume Conductors. IEEE Trans. Biomed. Eng. Р. 868 - 872. 1984.

5. PlonseyR. Bioelectric Phenomena. McGraw Hill, New York. 1969.

6. Pethig R. Dielectric Properties of Body Tissues. Clin. Phys. Physiol. Meas. 8 (suppl A). Р. 5 - 12. 1987.

7. Duck F.A. Physical Properties of Tissue - A Comprehensive Reference Book. Academic Press, London. Р. 167 - 223. 1990.

8. Brown B.H. Tissue Impedance Methods. Prog. Med. Enviro. Phys. Р. 85 - 110. 1983.

9. Barber D.C., Brown B.H. Progress in Electrical Impedance Tomography. Inverse Problems in Partial Differential Equations. D. Colton, R. Ewing, and W. Run-dell(eds). Society for Industrial and Applied Mathemetics, Philadelphia. 151 -164. 1990.

10. Henderson R.P., Webster J.G. An Impedance Camera for Spatially Specific Measurements of the Thorax. IEEE Trans. Biomed. Eng. Р. 250 - 254. 1978.

11. Nopp P., Rapp E., Pfutzner H., Nakesch H. andRuhsam Ch. Dielectric Properties of Lung Tissue as a Function of Air Content. Phys. Med. Biol. Р. 699 - 716. 1993.

12. Visser K.R. Electric Properties of Flowing Blood and Impedance Cardiography. Ann. Biomed. Eng. Р. 463 - 473. 1989.

А.С. Рябоконь МЕТОДЫ СНИЖЕНИЯ ЭНЕРГОПОТРЕБЛЕНИЯ НОСИМЫХ КАРДИОМОНИТОРОВ

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

Данная задача реализуется с использованием GPS-модуля и GSM-модуля, поддерживающего технологию передачи данных GPRS. Подсистема связи и привязки на местности должна обеспечивать:

- минимальные затраты дополнительных ресурсов ядра системы;

- время работы с автономным питанием - не менее 4-6 часов;

- режимы передачи данных - циклический и по запросу.

Применение собственного процессора GSM- или GPS-модуля для

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