CC BY
46
. NS
D =-------------—
So (Nb -1)
m ,,
Xi
JVS^V So
Nr zz(x *- x ,*)2
m
= 2Z±
NSo i=1 j=1
N — —2'
S0 (Xl - Xk )
NS -1
o
Выражение, стоящее в квадратных скобках, представляет собой исправленную дисперсию к-ого признака оО С учетом этого, количественная оценка меры компактности области, характеризующей референтное состояние в метрике Махаланобиса, будет равна удвоенной размерности пространства признаков
Д = Д = 2т.
Таким образом, критерий количественной оценки состояния некоторого объекта х может быть представлен в следующем виде:
1
2mNc
-о
Z dM (x> x ).
Выделение показателей, по совокупности которых следует производить интегральную оценку
состояния биосистем, трудно формализуемая задача. Ее решение определяется целями проводимых исследований, ограничениями, накладываемыми на условия их проведения, используемой измерительной аппаратурой, уровнем знаний об объекте исследования и т.д. Вместе с тем, можно определить общий подход к их формированию, следующий из системных свойств данных. Он следует из того, что описание свойств любой системы укладывается в определенную иерархическую структуру [5], каждый уровень которой соотносится с соответствующими методами получения данных. Как правило, на уровне популяций и организмов используются эпидемиологические, социологические, психологические методы; на уровне органов или тканей применяются физиологические; на клеточном уровне - микрометрические, цитогенетические; на внутриклеточном уровне - биохимические и биофизические методы измерений и т.п. Поэтому комплекс показателей, который в совокупности отражает свойства соответствующего элемента структурного описания, должен учитываться при его интегральной оценке.
СПИСОК ЛИТЕРАТУРЫ
1. Миронкина Ю.Н., Бобров А.Ф. Информационная технология статистического синтеза критериев и алгоритмов оценки функционального состояния человека в прикладных медикобиологических исследованиях // Информационные технологии. -1998. -№ 3. -С. 41—47.
2. Генкин А.А. Новая информационная технология анализа медицинских данных (программный комплекс ОМИС). — СПб.: Политехника, 1999. — 191 с.
3. Конрадов А.А. Статистические подходы к анализу многомерных гетерогенных биологических систем // Радиационная биология, радиоэкология. —1994. — Т 34. — Вып. 6. — С. 877-886.
4. Ту Дж., Гонсалес Р. Принципы распознавания образов. — М.: Мир, 1978. —416 с.
5. Клир Дж. Системология. Автоматизация решения системных задач / Пер. с англ. — М.: Радио и связь, 1990. — 554 с.
УДК 612.014.461
ИДЕНТИФИКАЦИЯ ПАРАМЕТРОВ МОДЕЛИ МОНИТОРИНГА ВОДНОГО БАЛАНСА ОРГАНИЗМА В РЕАЛЬНОМ ВРЕМЕНИ
А.Л. Тестов
Томский политехнический университет E-mail: [email protected]
Рассмотрены вопросы идентификации параметров математической модели мониторинга водного баланса организма в реальном времени, показано, что при соблюдении ряда условий методическая погрешность измерения величин жидкостных секторов с помощью данной модели не превысит Б...6 %, что позволит добиться более точных результатов по сравнению с аналогичными методиками.
Вода имеет ведущее биологическое значение в метаболических процессах и транспорте веществ организма. При различных патологиях происходят нарушения водного обмена, которые трудно поддаются коррекции, а их диагностика сложна. Владея методами оценки количества воды в организме и перераспределения ее по секторам в реальном времени, можно обнаружить скрытое накопление жидкости на ранних стадиях, что позволит предуп-
редить развитие осложнений и провести соответствующие действия до появления у больного патофизиологических и метаболических изменений. Одним из методов, позволяющих осуществить такую диагностику, является двухчастотная биоимпе-дансметрия [1], суть которой заключается в измерении параметров импеданса тела человека на разных частотах и оценки изменений водных объемов организма в соответствии измеренными парамет-
рами с помощью математического моделирования.
В работе [1] была разработана математическая модель оценки динамики объема водных секторов, которая может служить основой для построения автоматизированной биотехнической системы мониторинга динамики водного баланса организма. В этой статье мы рассмотрим метрологические аспекты идентификации параметров этой модели, для чего перепишем ее в следующем виде:
0,12-I -3 • Hf +
м3,
\п
3 - 0,7 •' 30Ш
ARз,
AR3,
AR3,
-1
= 0;
Н
V = V
,р,+1 кр1 Н
=
R3,
1
Rзo, Г
V = 1,05 •
Ь -1— V • (1 - н у
о„ R30,. хр1 1 ’
V = V + V
екж I кр1 ижг ?
V = V + V ,,
(1)
где Н - показатель гематокрита; АЛ300, АЛ30 - переменные и Л300, Л30 - постоянные составляющие активной части импеданса на частотах 300 и 30 кГц соответственно; Гр, Уж, Уиж, Уекж, ¥оео - объемы крови (КР), клеточной (КЖ), интерстициальной (ИЖ), внеклеточной (ВКЖ) жидкостей и общей воды организма (ОВО) соответственно; L - рост тела человека; о„ -удельная электропроводность плазмы; - количество измерений.
В системе ур. (1) погрешность идентификации трех величин определить достаточно просто: Гкр характеризуется погрешностью оценки показателя гематокрита, Г1жж и Гоео - алгебраической суммой абсолютных ошибок идентификации объемов крови и ИЖ, объемов КЖ и ВКЖ соответственно. Поэтому выражения определения этих параметров исключим из рассмотрения и остановимся подробнее на оставшихся уравнениях.
Для оценки ошибки идентификации показателя гематокрита воспользуемся математической моделью электропроводности крови на разных частотах, из которой было получено первое выражением в системе (1). Запишем:
1 - н + В • н
А = ■
(2)
[2]. Ограничиваясь линейной частью приращения функции А, запишем:
АА =
дА АН + дА
дН дВ
АВ,
(3)
где АА, АН, АВ - некоторые приращения (абсолютные погрешности) параметров А, Н, В.
Учитывая связь между абсолютной и относительной погрешностями: 8А=АА/А, 5В =АВ/В, 5Н=АН/Н, из выражения для предельной абсолютной погрешности функции А (3) можно получить выражение для определения предельной относительной погрешности показателя гематокрита 5Н:
5 А • (В + С)-5В • В
5Н = -
(4)
(1 - н у2
где Л=АЛ300/АЛ30, Б=оц/оп - отношение удельных электропроводностей цитоплазмы эритроцитов и плазмы крови.
Для определения чувствительности погрешности оценки А=ЛН,Б) используем общий прием, заключающийся в определении частных производных
В + 0,5 +1,5• ВС ’
где 5А, 5В-предельные относительные ошибки измерения отношения амплитуд переменных составляющих импеданса тела на разных частотах и отношения удельной электропроводности цитоплазмы эритроцитов и плазмы крови соответственно; С=(1-Н)/Н.
Для начала определим, в каком диапазоне могут изменяться измеряемые параметры В и Н. Изменение величины В напрямую зависит от удельной электропроводности цитоплазмы оц и плазмы оп крови. По понятным причинам каких-либо строгих величин для этих показателей не существует; даже в нормальном состоянии у разных людей они могут различаться -есть смысл говорить о средних и предельно допустимых значениях. Для оценочных расчетов возьмем (оп)шп равной 1,05 См/м, О)^ - 1,85 См/м [3, 4]. Принимая во внимание, что удельная электропроводность цитоплазмы всегда меньше таковой для плазмы, а также на основании данных ТЛ. Челидзе [3], примем (оц)га,п=0,6 См/м, (оц)тах=1,6 См/м. Таким образом, теоретически величина В может изменяться в пределах 0,37...1,0.
Гематокрит более важный показатель, чем оц и оп, широко используется в клинической практике, поэтому его нормы более четко определены. У здоровых людей гематокрит венозной и капиллярной крови равен 40.48 % (или 0,40.0,48) для мужчин и 36.42 % (или 0,36.0,42) для женщин [5].
Уменьшение гематокритной величины наблюдается при анемии, иногда до значительных цифр (20.25 %). Выраженное повышение (55.60 %) наблюдается при симптоматических эритроцито-зах, сопутствующих врожденным порокам сердца, легочной недостаточности, некоторым гемоглобинопатиям. Поэтому, допустимыми нормами для Н являются величины от 0,2 до 0,6.
Результаты расчетов относительной погрешности определения гематокрита представлены на рис. 1. Очевидно, что эта погрешность в физиологически значимом диапазоне не превышает 3.5 % при 5А=5В=1 %, а при отсутствии отклонений у пациента (нормальная величина показателя гема-токрита) ошибка составляет 1,5.2 %. Поэтому данную модель можно применять без опасения за искажение результата.
8H, %
H, o.e.
а) 8А = 8В = 2 %
sh, %
H, o.e.
б) 8А = 8В = 1 %
8H, %
H, o.e.
в) 8А = 1 %; 8В = 0,1 %
8H, %
H, o.e.
г) 8А = 0,1 %; 8В = 1 %.
Рис. 1. Изменение предельной относительной погрешности определения показателя гематокрита по измеренным значениям переменных составляющих активной части импеданса тела человека и удельной электропроводности цитоплазмы эритроцитов и плазмы крови. Сплошная кривая соответствует В = 0,4, точечная — В = 0,6, штриховая — В = 0,8
Кроме того, была произведена оценка степени влияния каждой составляющей на величину 8Н (рис. 1, в, г): анализ показал, что она весьма чувствительна к изменению параметра А и относительно слабо к отношению В. Точно определить это отношение во время мониторинга не представляется возможным (во-первых, это займет какое-то время, во-вторых, необходимо произвести забор крови, т.е. все лабораторные методы оценки ац и ап являются инвазивными). Если допустить, что В у разных людей будет заметно изменяться в силу индивидуальных особенностей, то эта погрешность будет оказывать влияние на абсолютную величину показателя гематокрита, но не на динамику его изменений. В численном выражении это выглядит так: отклонение отношения В от его среднего значения на 10 % может привести к увеличению погрешности 8Н до 5.6 % (при 8А = 1 %), что является вполне приемлемым результатом.
Теперь рассмотрим выражение из системы ур. (1), описывающее изменение объема клеточной жидкости тела человека. Для определения абсолютной и относительной погрешностей оценки количества КЖ воспользуемся тем же приемом, что и для гематокрита:
АУ =
дУ
дЬ
АЬ +
дУ
дУ„
дЛ3,
АЛ300 +
да
дУ„_
Аа +
дЛ30
(5)
где АУкж, А1, Аап _ некоторые приращения (абсолютные погрешности) параметров Укж, Ь, ап соответственно.
По аналогии с вышеизложенным для показателя гематокрита, заметим, что ошибка идентификации тем меньше, чем меньше каждая составляющая этой погрешности, и, что она (ошибка) изменяется в соответствии с вкладом каждой из этих составляющих в общую погрешность измерения. Поэтому также проанализируем возможные изменения относительной погрешности определения объема КЖ с целью определения факторов, влияющих на точность оценки данного параметра. Учитывая связь между абсолютной и относительной погрешностями, из выражения для предельной абсолютной погрешности (6) можно получить выражение для определения предельной относительной погрешности объема КЖ АУкж:
ЗУ = 23 Ь + 3а +
Л
Л
Л30 Л300
(6)
х3Л30 + л _ л
-3Л3,
где ЗЬ, 3ап, 3Я300, 3Я30 _ предельные относительные ошибки измерения Ь, ап, Я300, Я30.
Для начала выясним с какой точностью мы можем определить величины Ь и а„, так как напрямую от изменений этих параметров 3Укж не зависит, а
только от их относительных ошибок измерения. С хорошей достоверностью можно определить рост человека: обычно его измеряют с точностью до 1 см, т.е. ошибка составляет приблизительно 0,5 %. Несложно уменьшить эту ошибку до 0,2...0,3 %, измерив Ь с точностью 0,5 см. Но в клинических условиях часто нет возможности проводить подобные измерения, и чаще всего принимаются данные, предоставляемые самим пациентом, поэтому остановимся на 3Ь равной 0,5 %. С погрешностью определения удельной электропроводности цитоплазмы клеток все гораздо сложнее. В нормальном состоянии этот показатель в среднем равен 0,68 См/м, а у разных людей он колеблется в пределах 0,6...0,8 См/м [3, 4]. Таким образом, разброс значений у разных людей даже в нормальном состоянии составляет 10.15 %, что уж тут говорить о патологиях. Естественно, измерить удельную электропроводность цитоплазмы клеток во время мониторинга не представляется возможным, поэтому мы вынуждены принять эту величину равной среднему значению и по возможности вносить коррективы. Учитывая, что изменение электролитного состава клеток происходит относительно медленно, сделаем допущение, что во время оперативного контроля оно не повлияет на оценку динамики объема клеточной жидкости, а лишь на абсолютное его значение, поэтому для нашего исследования примем 3ап равной 0,01.
Очевидно, что выражение (6) существенно зависит от Я300 и Я30. Попытка проанализировать их возможные значения для минимизации погрешности идентификации показала, что при предельных относительных погрешностях величин, входящих в формулу (6), не превышающих 1 %, и 3Укж, равной 5 %, Я300 и Я30 должны равняться 210 и 630 Ом соответственно (коэффициент к = Я300/Я30 равен 0,33). Зависимость предельной относительной погрешности определения объема клеточной жидкости от коэффициента к представлена на рис. 2, из которого становится видно, что 3Укж <5 % при к<0,4.
%
Рис. 2. Динамика предельной относительной погрешности идентификации объема клеточной жидкости при различных значениях соотношения постоянных составляющих активной части импеданса тела человека на разных частотах к (3L =3ап =3Н300 =3Н30 = 1 °%)
Практические исследования свидетельствуют о других значениях к: в пределах 0,7.0,87 в различных условиях (пациенты, электроды, токопроводящая паста и др.), т.е. 3Укж может составить 10.20 %. Очевидно, что эта ошибка влияет на абсолютные значения 3Укж, носит систематический характер и не зависит от погрешности определения парамет-
к. о.е
ров, входящих в модель. Поэтому в дальнейшем для ее уменьшения будет необходимо ввести поправочный коэффициент в выражение, описывающее объем КЖ.
Чтобы оценить степень влияния погрешностей измерения импеданса тела человека на ошибку идентификации объема КЖ, будем считать, что поправки уже внесены и поэтому для расчетов примем Я300 = 210 Ом, Я30 = 630 Ом. Результаты этих расчетов представлены на рис. 3, из которого видно, что 8Укж не превышает 3.5 % (приемлемый результат) при 8Я300= 8Я30=1 %, а при увеличении 8Я300 и 8Я30 до 5 % ошибка идентификации составляет 10.12 %. Поэтому данную модель можно применять без опасения за искажение результата при точности измерения постоянных составляющих импеданса тела на разных частотах равной 1 %.
ЛУ =
дУ
дЬ
ЛЬ +
дУ
до„
Лст +
дУ
дН
ЛН +
дУ
дЯ3
ЛЯ30 +
дУ
дУ„,
ЛУ
28Ь + 8<уп + 8Я30 + 5 (1,5
Т?н 8Н+У)
(7)
(8)
где
5 =
(1 - Н )1,
У и
„р п
Я
Величина 8Уиж из (8) зависит как от значений измеряемых параметров, входящих в модель (1), так и от ошибок их регистрации. В таблице приведены результаты расчетов ошибки идентификации объема ИЖ при изменении значения одного из па-
раметров или его погрешности на 10 %. На основании этих данных можно сделать вывод, что на возрастание 8Уиж в большей степени чем остальные влияют 8Ь, 8ап, и 8Я30. Таким образом, для уменьшения результирующей ошибки необходимо их минимизировать.
Таблица. Относительные погрешности идентификации объема интерстициальной жидкости.
ЪКзоо, %
Рис. 3. Динамика предельной относительной погрешности идентификации объема клеточной жидкости при различных значениях погрешности измерения постоянной составляющей активной части импеданса тела человека на частоте 300 кГц 8 = 0,5 °%, 8о„ = 1 %). Сплошная кривая соответствует 8Н30 = 1 %, точечная — 8Н30 = 3 %, штриховая — 8Н30 = 5 %
Осталось проанализировать уравнение, описывающее изменение объема интерстициальной жидкости тела человека (1). Аналогично вышеизложенному для других уравнений модели, а также учитывая связь между абсолютной и относительной погрешностями, запишем выражения для определения предельной абсолютной и относительной погрешности объема ИЖ ЛУиж (7) и 8Уиж (8):
L, м Ы н, о.е. ЪН См/м Ъ°п R30, Ом Щ0 V,,, м3 ъv>р ЪVltж
1,8 0,01 0,5 0,01 0,68 0,01 500 0,01 0,005 0,01 0,055
1,62 0,01 0,5 0,01 0,68 0,01 500 0,01 0,005 0,01 0,059
1,98 0,052
1,8 0,005 0,5 0,01 0,68 0,01 500 0,01 0,005 0,01 0,043
0,02 0,079
1,8 0,01 0,45 0,01 0,68 0,01 500 0,01 0,005 0,01 0,057
0,55 0,053
1,8 0,01 0,5 0,005 0,68 0,01 500 0,01 0,005 0,01 0,053
0,02 0,058
1,8 0,01 0,5 0,01 0,612 0,01 500 0,01 0,005 0,01 0,053
0,748 0,057
1,8 0,01 0,5 0,01 0,68 0,005 500 0,01 0,005 0,01 0,049
0,02 0,067
1,8 0,01 0,5 0,01 0,68 0,01 450 0,01 0,005 0,01 0,053
550 0,057
1,8 0,01 0,5 0,01 0,68 0,01 500 0,005 0,005 0,01 0,049
0,02 0,067
1,8 0,01 0,5 0,01 0,68 0,01 500 0,01 0,0045 0,01 0,053
0,0055 0,057
1,8 0,01 0,5 0,01 0,68 0,01 500 0,01 0,005 0,005 0,054
0,02 0,057
5УМЖ, %
5Язо, %
Рис. 4. Динамика предельной относительной погрешности идентификации объема интерстициальной жидкости при различных значениях погрешности измерения постоянной составляющей активной части импеданса тела человека на частоте 30 кГц (^ = 0,5 °%, 8<г„ = 1,0 %, 8Н =8VKр = 2,0 %)
Обоснование для оценки погрешности определения величин Ь и стп, приведено выше, поэтому для настоящего случая также примем 8Ь равной 0,5 %, а 8оп — 1 %. Результаты расчетов относительной погрешности определения показателя гематокрита представлены на рис. 1, из которого видно, что при отсутствии отклонений у пациента (нормальная величина показателя гематокрита) 8Н = 2 %. Изменения объема крови в математической модели мониторинга водного баланса (1) зависит только от изменений показателя гематокрита (исключение составляют случаи острой кровопотери), следовательно, 8Укр можно принять равным 8Н, т.е. 2 %. В результа-
те получаем, что 8Уиж не превышает 5 % при 8Я30 = 1.1,5 %, а при увеличении 8Я30 до 5 % ошибка идентификации может составить 10 % (рис. 4).
Таким образом, на основании вышеизложенного можно сделать ряд заключений: относительная погрешность измеряемых составляющих импеданса (ЛЯ300/ЛЯ30, Я300, Я30) не должна превышать 1 %, а роста (Ь) — 0,5 %. Также необходимо сделать допущение, что разброс параметров оц и а„ от их сред-
него значения не превысит 1 % во время мониторинга. В этом случае предельные относительные погрешности идентификации объемов крови, ИЖ, КЖ не превысят 5.6 %, что говорит о высокой достоверности результатов, полученных с помощью данной математической модели (1), по сравнению с аналогичными разработками на основе данной методики (ошибка до 20.30 %).
СПИСОК ЛИТЕРАТУРЫ
1. Тестов А.Л. Математическая модель оценки динамики водного баланса в реальном времени // Современные средства и системы автоматизации — гарантии высокой эффективности производства: Матер. III научно-практ. конф. — Томск: Изд-во Том. ун-та, 2003. —С. 255-258.
2. Новицкий П.В. Зограф И.А. Оценка погрешностей результатов измерений. —2-е изд., перераб. и доп. — Л.: Энергоатомиздат, 1991. —304 с.: ил.
3. Челидзе Т.Л. Электропроводность цитоплазмы эритроцитов // Биофизика. —1980. —Т. 25. — № 6. — С. 27—30.
4. Уманский О.С. Биотехническая система экспресс-оценки группы гематологических параметров в одной микропробе: Дис. ... канд. тех. наук. —Томск, 1996. —183 с.
5. Гольдберг Д.И. Справочник по гематологии с атласом микрофотограмм. — Томск: Изд-во ТГУ, 1989. —468 с.
