Определение составляющих импеданса биообъекта Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

УДК 681.335

DOI: 10.20310/1810-0198-2017-22-1-104-108

ОПРЕДЕЛЕНИЕ СОСТАВЛЯЮЩИХ ИМПЕДАНСА БИООБЪЕКТА

© Е.И. Глинкин, Е.А. Рухлова

Тамбовский государственный технический университет 392000, Российская Федерация, г. Тамбов, ул. Советская, 106 E-mail: glinkinei@rambler.ru

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

Ключевые слова: способ определения; постоянная времени; установившееся значение потенциала; импульсная динамическая характеристика; нормируемые параметры; погрешности

Предлагаемое изобретение относится к медицине и может быть использовано для оценки функционального состояния организма.

Известен способ бесконтактного измерения удельного электросопротивления [1], согласно которому измеряемый образец размещают на катушке индуктивности колебательного контура и измеряют изменение добротности контура, при этом индуктивность контура измеряют в диапазоне 135-155 МГц с помощью катушек Гельмгольца, а удельное электросопротивление образца определяют по формуле.

Данный способ обладает низкой точностью из-за изменения конструктивных параметров.

По способу измерения резистивной и емкостной составляющих комплексного сопротивления [2] измеряемое сопротивление периодически подключают вначале к эталонному источнику напряжения на априорно заданное время t\ , затем измеряемое сопротивление закорачивают, мгновенное значение падения напряжения и 1 на сопротивлении измеряют в конце временного промежутка tl . Через априорно заданное время 12 = Ь после закорачивания измеряют значение падения напряжения и 2 . Резистивную и емкостную составляющие вычисляют по формулам:

(

Rx = Ro

ln(l - U) _U2_

lnU2) U i

Л

-l

Cx = —

'-UO-^

Ro

lnl l

U o J К U i

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

Наиболее близким к заявленному техническому решению является способ определения составляющих импеданса биообъекта [3], заключающийся в том, что на биообъект подается через электроды импульс стабилизированного тока определенной полярности (например, положительной) и амплитудой /0. Вследствие емкостного характера реактивной составляющей импеданса биообъекта происходит переходной процесс нарастания напряжения на биообъекте, которое измеряется в фиксированные два момента времени tl и 12 после начала импульса тока, получая, соответственно, значения напряжения и 1 и и 2 . Измерение в момент времени 12 производится, когда емкость тканей биообъекта заряжена полностью и переходной процесс закончился. Величина стабилизированного тока /о выбирается такой, чтобы за время действия импульса тока произошел полный заряд емкости тканей биообъекта. Тогда напряжение на биообъекте пропорционально величине активной составляющей импеданса биообъекта. Активное сопротивление К биообъекта определяется по формуле (при параллельной схеме замещения биообъекта) К = и2110. Эквивалентная емкость С тканей биообъекта вычисляется с помощью выражения

C=i Ufi- R 1-

Недостатками прототипа являются: низкая точность из-за наличия динамической и методической погрешности и низкая оперативность, вызванные необходимостью ожидания установившегося режима ВАХ.

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

Данная техническая задача решается за счет того, что в способе определения составляющих импеданса биологического объекта [4], заключающемся в подаче

j

на биообъект импульса стабилизированного тока, измерении напряжения на биообъекте в фиксированные два момента времени после начала импульса тока, в отличие от прототипа, дополнительно измеряют амплитуду стабилизированного тока 10, моменты времени фиксации напряжения представляют собой t\ и 12 , причем 12 = 2 ^ ; а в качестве составляющих импеданса биообъекта определяют активное сопротивление К и эквивалентную емкость С тканей биообъекта, которые рассчитывают по формулам.

Сущность предлагаемого способа [4] поясняется на рис. 1-4.

Для определения составляющих импеданса биологического объекта на тело пациента в месте измерения сопротивления накладывают измерительные электроды, прикладывают напряжение на измерительную ячейку, состоящую из последовательно включенных измеряемого комплексного и эталонного сопротивлений (рис. 1). После включения напряжения питания на биообъект подают через электроды импульс стабилизированного тока, измеряют его амплитуду /0 (рис. 2б).

В моменты времени t\ и 12 , причем t2 = 2^ , фиксируют (рис. 2а) значения падений напряжения Ц 1 и Ц 2 соответственно, на эталонном сопротивлении ^ . По измеренным значениям напряжения и времени находят активное сопротивление и эквивалентную емкость тканей биообъекта через установившееся значение потенциала Е с постоянной времени Т.

Экспериментальная зависимость и(г) = и динамического процесса (рис. 2а) от импульса стабилизированного тока (рис. 2б) изменяется по экспоненциальному закону:

и = Е ■ (1 - е т )

(1)

Зависимость (1) связывает между собой измеряемое значение амплитуды и напряжения за время г исследования до установившегося значения Е потенциала с постоянной времени Т.

Параметры Е и Т однозначно определяют динамическую характеристику эксперимента по зависимости (1).

Регистрация параметров Е и Т организована по двум измеренным значениям амплитуды иь и2 напряжения в два момента времени гь г2 из системы уравнений по формуле (1) для первого и второго измерений:

и

и1 = Е ■ (1 - е т ),

- Ъ,

Щ = Е ■ (1 - е т ).

Рис. 1. Схема замещения измерительной ячейки

а)

б)

в)

Рис. 2. Экспериментальная зависимость динамического процесса: а - напряжения; б - импульса тока; в - тока

и запишем отношение:

и-у

t 1п(1 - -2 ) _Е_

^ = 1п(1 -

Е

Решение в явном виде получено при кратном от-t2 „

ношении — = 2 после приведения к общему знамена-

Выразим из уравнений системы ^ и г2:

ч =-т ■ 1п(1 - Цю

t2 =-Т ■ 1п(1 - Ц2)

1п(1 - Ц2) = 2 ■ 1п(1 - Цт) .

(2) Е Е

Проэкспоненциируем данное уравнение и выразим параметр Е:

IJr

1» 1Я

Рис. 3. Динамические характеристики: 1 - исследуемая; 2 - эквивалентная

Рис. 4. Погрешность

E —-

U,2

2 U - U 2

(3)

Для нахождения Т подставим выражение (3) в первое уравнение системы (2):

т — ■ -t'

ln(U2 -1) U )

(4)

где и1 - напряжение на биообъекте в момент времени tl ; и 2 - напряжение на биообъекте в момент времени 12 .

С помощью параметров Е и Т определяют значение активного сопротивления:

R = E/I0 .

Это обусловлено значением тока:

r dU CE

I — C-—-e т ■

dt T T

(5)

CE _E _ 1—1 н — RC ~ R ~10'

где начальный ток 1Н тождественен амплитуде стабилизированного тока 10. Эквивалентная емкость тканей биообъекта в свою очередь определяется как:

T

C — -, R

т. к. Т = Я ■ С .

Адекватность и эффективность способа представлены ниже.

1. Адекватность предлагаемого способа физике эксперимента доказывает математическое моделирование исследуемой и (Г) ИДХ 1 относительно эквивалента 2 экспериментальной иэ (t) ИДХ. По полученным значениям К и С определяется значение периода Т (согласно формулы Т = R•C), строятся исследуемая 1 и эквивалентная 2 ИДХ (рис. 3).

Затем проводится оценка адекватности полученных зависимостей по формуле определения относительной погрешности:

U э (t) - Ui (t)| -100%

и э (О

ее оценка представлена на рис. 4. При этом погрешность 8 отклонения и, (t) относительно иэ (t) не

превышает 1,5-10-13 %.

2. Повышение точности за счет методической и динамической погрешности приведем на примере активного сопротивления:

R — U — . 1

El 1 -

e т

1 н • e т

где RH = const - информативный параметр ИДХ сопротивления.

Эффективность по точности определяется нелинейностью п сопротивления R относительно постоянного сопротивления RH предлагаемого способа:

_R__ V Ян

ет

-1

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

3. Динамическая погрешность 8 определяется нелинейностью п:

Si:

Я0 - R

Ro

1 -

R_

Ro

= 1 -л|,

В начальный момент времени t = 0, когда e = 1:

Таблица 1

Оперативность предлагаемого способа

Результаты Эквивалент Прототип Предлагаемый способ е % Епрот.' е '% Епр.сп.

Е, мВ 120 114 120 5 6,4-10-12

T.10-7, с 42 21,1 42 49,7 6,7-10-12

R, кОм 34,3 32,6 34,3 5 6,5-10-12

С, пФ 123 65 123 47,1 2,7-10-13

Т- е- Si =

и также растет по экспоненте (рис. 2в) с

увеличением времени Г, в то время как мгновенное значение и ИДХ стремится по асимптоте к установившемуся потенциалу (рис. 2а).

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

4. Повышение оперативности предлагаемого способа оценивается эффективностью времени измерения Г. В предлагаемом способе t < Т измерения не превышает постоянную времени, а для прототипа в 3-5 раз больше ^ = (3 - 5)Т для погрешности (5-1)% определения установившегося потенциала Е.

Из эффективности, ц = ■(3—5)Т = (3 - 5) следует, t Т

что оперативность предлагаемого способа в 3-5 раз выше известных способов.

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

Анализ табл. 1 показывает, что точность предлагаемого метода на несколько порядков выше за счет учета динамической погрешности и устранения методической погрешности.

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

ФОРМУЛА ИЗОБРЕТЕНИЯ

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

времени после начала импульса тока, отличающийся тем, что дополнительно измеряют амплитуду стабилизированного тока 10 , моменты времени фиксации напряжения представляют собой t\ и 12, причем 12 = 2^; а в качестве составляющих импеданса биообъекта определяют активное сопротивление К и эквивалентную емкость С тканей биообъекта, которые рассчитывают по следующим формулам:

л = Е, 10

где Е - установившееся значение потенциала с постоянной времени Т,

E = ■

U2

T = -

2Ui - U 2 ti

ln(— -1) Ui

где и! и и2 , соответственно, напряжение на биообъ-

Т

екте в моменты времени ^ и 12 ; при этом С = — .

К

СПИСОК ЛИТЕРАТУРЫ

А. с. № 1642410 СССР. Способ бесконтактного измерения удельного сопротивления. МПК5 G 01 R 27/02. 1991. Бюл. № 14. Патент № 2003123 РФ. Способ измерения резистивной и емкостной составляющих комплексного сопротивления. G 01 R 27/26. 1993. Бюл. № 41.

А. с. № 1397024 СССР. Способ определения составляющих импеданса биообъекта. А61 В5/05. 1988. Бюл. № 19. Патент № 2509531 РФ. Способ определения составляющих импеданса биообъекта / А.В. Наумова, Е.И. Глинкин. А 61 В 5/053. 2014. Бюл. № 8.

Поступила в редакцию 25 ноября 2016 г.

Глинкин Евгений Иванович, Тамбовский государственный технический университет, г. Тамбов, Российская Федерация, доктор технических наук, профессор, профессор кафедры биомедицинской техники, заслуженный изобретатель Российской Федерации, e-mail: glinkinei@rambler.ru

Рухлова Евгения Алексеевна, Тамбовский государственный технический университет, г. Тамбов, Российская Федерация, магистрант по направлению подготовки «Биотехнические системы и технологии», кафедра биомедицинской техники, e-mail: eugeniya.ruhlova2011@yandex.ru

UDC 681.335

DOI: 10.20310/1810-0198-2017-22-1-104-108

DETERMINING THE IMPEDANCE COMPONENTS OF A BIOLOGICAL OBJECT

© E.I. Glinkin, E.A. Rukhlova

Tambov State Technical University 106 Sovetskaya St., Tambov, Russian Federation, 392000 E-mail: glinkinei@rambler.ru

A method of detecting the resistance and the equivalent capacity of biological tissue by dynamic impulse response with normalized parameters is proposed. These parameters are calculated by the voltage amplitude which changes in two time-points.

Key words: method of determining; time constant; steady value of potential; pulse dynamic characteristic; normalized parameters; error

REFERENCES

1. Sposob beskontaktnogo izmereniya udel'nogo soprotivleniya [The method of nonintrusive measurement of resistivity constant]. Author's license no. 1642410 SSSR, MPK5 G 01 R 27/02, 1991. (In Russian).

2. Sposob izmereniya rezistivnoy i emkostnoy sostavlyayushchikh kompleksnogo soprotivleniya [The way of resistance measurement and capacitive impedance]. Patent no. 2003123 RF, G 01 R 27/26, 1993. (In Russian).

3. Sposob opredeleniya sostavlyayushchikh impedansa bioob"ekta [The way of determining the impedance components of a biological object]. Author's license no. 1397024 SSSR, A61 B5/05, 1988. (In Russian).

4. Naumova A.V., Glinkin E.I. Sposob opredeleniya sostavlyayushchikh impedansa bioob"ekta [The way of determining the impedance components of a biological object]. Patent no. 2509531 RF, A 61 B 5/053, 2014. (In Russian).

Received 25 November 2016

Glinkin Evgeniy Ivanovich, Tambov State Technical University, Tambov, Russian Federation, Doctor of Technics, Professor, Professor of Biomedical Technics Department, Honored Inventor of Russian Federation, e-mail: glinkinei@rambler.ru Rukhlova Evgeniya Alekseevna, Tambov State Technical University, Tambov, Russian Federation, Master's Degree Student on Training Direction "Biotechnical Systems and Technologies", Biomedical Technics Department, e-mail: eugeniya.ruhlova2011 @yandex.ru

Информация для цитирования:

Глинкин Е.И., Рухлова Е.А. Определение составляющих импеданса биообъекта // Вестник Тамбовского университета. Серия Естественные и технические науки. Тамбов, 2017. Т. 22. Вып. 1. С. 104-108. DOI: 10.20310/1810-0198-2017-22-1-104-108

Glinkin E.I., Rukhlova E.A. Opredelenie sostavlyayushchikh impedansa bioob"ekta [Determining the impedance components of a biological object]. Vestnik Tambovskogo universiteta. Seriya Estestvennye i tekhnicheskie nauki — Tambov University Reports. Series: Natural and Technical Sciences, 2017, vol. 22, no. 1, pp. 104-108. DOI: 10.20310/1810-0198-2017-22-1-104-108 (In Russian).