CC BY
57
сигнала и дрейфа ЦЛМ в частотной области, такая замена возможна. В случае же анализа длительности ЦЛМ во временной области и при исследовании тонкой структуры дрейфа частот ЦЛМ указанная замена не позволит получать воспроизводимые результаты, что накладывает существенные ограничения на перенос алгоритмов обработки данных ЭЭГ с использованием вейвлета Morlet на системы реального времени при анализе электромагнитных процессов, происходящих в головном мозге.
Литература
1. Канайкин А.М., Попов А.А., Рощина К.А. [и др.]. Обнаружение артефактов в сигнале электроэнцефалограммы с помощью вейвлет-преобразования // Электроника и связь. 2011. № 4. С. 126-130.
2. Туровский Я.А., Кургалин С.Д., Максимов А.В. Выбор анализирующих вейвлетов для системы с параллельной обработкой биомедицинских данных // Вестн. Воронеж. гос. ун-та: Сер. Системный анализ и информационные технологии. 2011. № 2. С. 74-79.
3. Лазоренко О.В., Лазоренко С.В., Черногор Л.Ф. Вейв-лет-анализ модельных сигналов с особенностями // Радиофизи-
ка и радиоастрономия. 2007. Т. 12. № 2. С. 182-204.
4. Torrence C., Compo P.G., A Practical Guide to Wavelet Analysis, Bulletin of the American Meteorological Society, 1998, Vol. 79, no. 1, pp. 61-78.
5. Зенков Л.Р. Клиническая электроэнцефалография (с элементами эпилептологии). М.: МЕДпресс-информ, 2011. 356 с.
References
1. Kanaikin А.М., Popov А.А., Roshchina К.А., Chertov O.R., Shashkov V.A., Elektronika i svyaz. Tematich. vyp. «Elektronika i nanotekhnologii» [Electronics and Communications. Special issue «Electronics and nanotechnologies»], 2011, no. 4, pp. 126-130.
2. Turovski Ya.A., Kurgalin S.D., Maksimov A.V., Vestnik Voronezhskogo gos. univ., Ser. Sisnemny analyz i inform. tekhno-logii [Proc. of Voronezh State Univ., Series. System analysis and IT], 2011, no. 2, pp. 74-79.
3. Lazorenko O.V., Lazorenko S.V., Chernogor L.F., Radiofizika i radioastronomia [Telecommunications and Radio Engineering], 2007, Vol. 12, no. 2, pp. 182-204.
4. Torrence C., Compo P.G., Bulletin of the American Meteorological Society, 1998, Vol. 79, no. 1, pp. 61-78.
5. Zenkov L.R. Klinicheskaya elektroentsephalografiya s elementami epileptologii [Clinical electroencephalography with elements of epileptology], Moscow, MEDPress-Inform, 2011, 356 p.
УДК 519.711.3; 629.78.07-051:612.014.4; 613.693
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ДИНАМИКИ ГИПОКСИЧЕСКИХ СОСТОЯНИЙ ЧЕЛОВЕКА
(Работа поддержана РФФИ, грант № 12-08-01273-а)
А.В. Богомолов, д.т.н., профессор, ведущий научный сотрудник; Ю.А. Кукушкин, д.т.н., профессор, ведущий научный сотрудник; М.В. Дворников, д.м.н., профессор, начальник отдела (Научно-исследовательский испытательный центр (авиационно-космической медицины и военной эргономики) 4 Центрального научно-исследовательского института Минобороны России, Петровско-Разумовская аллея, 12а, г. Москва, 127083, Россия, [email protected], [email protected], [email protected])
Рассматривается математическое обеспечение моделирования динамики гипоксических состояний человека в интересах рискометрии его здоровья и работоспособности в условиях низкого барометрического давления (высотные полеты, альпинизм, космонавтика и др.). Математическая модель дозы гипоксии в виде интеграла интенсивности гипоксического воздействия в пределах времени действия этого фактора дает возможность перейти в динамическую область исследования эффектов воздействия гипоксической газовой среды (оценивать ее влияние на состояние человека в реальном времени), позволяя рассчитать резервное время сохранения работоспособности человека в условиях гипоксической гипоксии и время восстановления работоспособности после прекращения воздействия гипоксической газовой среды. С помощью разработанной математической модели обеспечивается расчет адекватных оценок вероятностей рассматриваемых неблагоприятных эффектов при любых значениях скоростей и профилей изменения парциального давления кислорода в трахее человека, кратности и времени воздействия гипоксической среды. Верификация математической модели осуществлена по результатам вычислительного эксперимента с использованием данных реального авиационного инцидента с разгерметизацией салона самолета Ту-154. Полученные результаты обеспечивают поддержку принятия решений при выборе оптимальной траектории полета летательных аппаратов в случаях разгерметизации кабины или салона, при оказании помощи раненым и пораженным в чрезвычайных ситуациях в условиях высокогорья и в ряде других практических задач, связанных с пребыванием человека в условиях гипоксической гипоксии.
Ключевые слова: рискометрия здоровья, моделирование гипоксических состояний, опасные состояния человека.
MATHEMATICAL MODELING OF THE HUMAN HYPOXIC CONDITIONS DYNAMICS Bogomolov A V., Ph.D., professor, leading researcher; Kukushkin Yu.A., Ph.D., professor, leading researcher;
Dvornikov M. V., Ph.D., professor, head of department (Research and Testing Center (Aerospace Medicine and Military Ergonomics) 4 Central Research Institute of the Ministry of Defense of Russia, Petrovsky-Razumovsky alley, 12a, Moscow, 127083, Russia, [email protected], [email protected], [email protected])
Abstract. The article describes the simulation software of the human hypoxic conditions dynamics according to health risk metrics and human performance in terms of low barometric pressure (high altitude flights, mountaineering, space, etc.). The hypoxia dose mathematical model as the integral intensity of hypoxic impact within the time of this factor activity allows passing into the dynamic area of study the effects of hypoxic gas environment (to size up its impact on the human condition in real time). It allows calculating the buffer time of human efficiency saving in conditions of hypoxic hypoxia and recovery time after cut-off of. The developed model provides calculation of probabilities adequate valuation of adverse effects with any velocity profiles and oxygen partial pressure in a trachea, multiplicities and exposure time of hypoxic outward influence.
A mathematical model verification based on the results of computer simulation using data from real aircraft cabin depressuriza-tion incident with the Tu-154.
These results provide decision support when choosing the optimal flight path of aircrafts in the case of cabin depressurizati on, in assisting the injured and affected in emergencies in the high altitude and in a number of other problems associated with the presence of human beings in hypoxic hypoxia.
Keywords: health risk, hypoxic conditions modeling, dangerous condition of the person.
Ситуации, связанные с быстротекущим воздействием на человека среды с малым содержанием кислорода в окружающем воздухе, приводят к развитию гипоксических состояний (состояний кислородного голодания как всего организма в целом, так и отдельных органов и тканей), обусловливающих высокий риск потери сознания человеком и даже его гибели [1-5].
В интересах контроля работоспособности человека в условиях гипоксии необходимы прогностические оценки состояния человека с учетом влияния изменяемой гипоксической среды на его организм. Известные математические модели [1, 4] позволяют рассчитать такие оценки в статических условиях, то есть при сохранении постоянства интенсивности гипоксического воздействия в течение анализируемого интервала времени, что существенно ограничивает область их применения. Преодолеть указанный недостаток можно за счет разработки новых теоретических подходов к нормированию воздействия гипоксической газовой среды на организм человека и моделирования гипоксических состояний, базирующихся на дозо-вом принципе.
Математическая модель динамики гипоксических состояний
Аналитически доза гипоксии определяется вы-
ражением
D(t) = J[plo"Pm(t)]dt, где D(t) -
экспозиционная доза гипоксического воздействия; Р™5оо =37,8 кПа - парциальное давление кислорода в трахее при значении барометрического давления 44 кПа (соответствует высоте над уровнем моря 6,5 км); рт(() - текущее значение парциального давления кислорода в трахее; Т - интервал интегрирования.
Парциальное давление кислорода в трахее связано с барометрическим давлением зависимостью рт = F (р - рН0), где ^ - относительная концентрация кислорода в атмосфере (^=0,21); р - барометрическое давление; рНг0 - парциальное давление насыщенного водяного пара при температуре 310 °К (37 °С), равное 6,27 кПа (47 мм рт. ст.).
При численном интегрировании выражение для дозы гипоксии можно записать в виде
п
О = Е [рт - рт (^) ] &, где
1=1
Разность р™ш - рт (() определена как интенсивность гипоксического воздействия (V).
В высотной физиологии [1, 4] получены эмпирические результаты, позволяющие рассчитать оценки минимального (), среднего () и
максимального (ТПС ) резервного времени потери сознания при пребывании человека на различных высотах (при разном барометрическом давлении).
По этим данным методом Левенберга-Марк-вардта синтезированы регрессионные уравнения, описывающие зависимости времени потери сознания от интенсивности гипоксического воздействия. Эти уравнения имеют следующий вид:
Г™, = 186,15 U~
> грП ? 1 П
Статистическая значимость синтезированных уравнений по критерию Фишера соответствует уровню значимости р<0,001.
С учетом результатов структурной идентификации законов распределения резервного времени потери сознания получены аналитические выражения для расчета оценок вероятностей потери сознания в зависимости от времени t пребывания в гипоксической газовой среде, интенсивность которой постоянна, то есть
w; (t) =
0,5 - Ф
3,224х-
0,5 + Ф
t- Тпс (U)
1 i mean V /
Тпс (U) - Тпс. (U)
' mean V / i min V /
еСЛИ ( U- Т^ап (U) ) < 0,
ti- Тпс (U)
i i mean
3,224х
Тпт (и) - ТпС (и)
, / 1 теап V /
если ( Г- ТГтеап{Р) ) > 0 , где Ф - функция Лапласа.
Параметры этой модели определены при соблюдении условия постоянства барометрического давления в течение всего времени пребывания человека в газовой среде, что затрудняет ее использование применительно к динамическим условиям. Однако дозовый подход позволяет преодолеть
™ = 710.37 TJ •
эти трудности за счет того, что, определив дозу гипоксии при изменении рт и разделив ее значение на время измерения, можно вычислить вели-
п
чину АС = Д / 2 -Аt, эквивалентную получен-
1=1
ной дозе гипоксии при соблюдении условия постоянства парциального давления кислорода в трахее (при полете на постоянной высоте), а следовательно, и однозначно связанных с ним барометрическим давлением и высотой над уровнем моря. Тогда эквивалентная дозе величина парциального давления кислорода в трахее будет равна
рт = Арт + рт .
Эта величина используется для динамических расчетов параметр°в Т™еш, т Опреде-
лив значения этих параметров и зная текущее время ti=iAt, можно рассчитать оценки вероятностей потери сознания.
Коэффициенты выражений, аппроксимирующих зависимости 7ПП, С! и ТПСх от и, равны
дозам, соответствующим указанным значениям времени переносимости гипоксического воздействия. С учетом этих значений получена модель, позволяющая определять вероятность потери сознания в зависимости от дозы гипоксии в произвольный момент времени и
п
£[/(/.) Д/- 367,23
3,224х ---- ,
181,01
п
если ( ££/(/.) А/-367,23 ) < 0 , /=1 I
А/-367,23
3,224х—--- ,
343,14
п
если ( Д/-367,23 ) > 0 ,
/=1 /
п
где 2 и )Аt - значение дозы гипоксии в момент
-=1
времени ti; At - интервал времени, через который определяется доза; п - количество измерений доп
зы; а = 1 + 0,0004542 и (^ ) Аt - 367,23.
-=1
Методология дозового подхода и реализующее ее математическое обеспечение позволяют рассчитать время сохранения работоспособности человека в условиях гипоксической гипоксии и ее восстановления после потери сознания, обусловленной воздействием гипоксии. Временем начала восстановления работоспособности человека является момент изменения знака дозы гипоксии, а временем полного восстановления сознания - момент, в который наступает равенство по модулю положительных и отрицательных доз гипоксии.
Предложенный метод обеспечивает расчет адекватных оценок вероятностей рассматривае-
мых неблагоприятных эффектов при любых значениях скорости и профиля изменения парциального давления кислорода в трахее, времени и кратности воздействия гипоксической среды.
Верификация математической модели
на основе вычислительного эксперимента
Для верификации разработанной математической модели проведен вычислительный эксперимент с использованием данных реального авиационного инцидента, произошедшего на самолете Ту-154 в 1998 году. Фабула происшествия изложена в публикации «Новой газеты» (№ 22 (494) 1998 г. ). Кроме того, были использованы результаты расшифровки параметров изменения высоты полета (над уровнем моря) и высоты в салоне самолета (определяется парциальным давлением кислорода в салоне) в период от 1 700-й до 3 000-й секунды полета по данным средств объективного контроля, предоставленные специалистами Межгосударственного авиационного комитета.
В период от 1 750-й до 1 800-й секунды полета имело место превышение избыточного давления в гермокабине (более 0,56 кгс/см2), которое на 1 830-й секунде (на 30-й минуте) сменилось резким (за 100 секунд) увеличением высоты в салоне (с полной его разгерметизацией).
Экипаж изменил эшелон полета с 10 000 м до 8 000 м, а затем до 7 500 м. Полет на этой высоте продолжался в течение 510 секунд (8,5 мин.). Затем самолет снизился до высоты 4 000 м и после устранения разгерметизации кабины к 3 800-й секунде (к 64-й минуте) полет был продолжен на высоте 9 000 м и завершен на аэродроме назначения к 5 700-й секунде (95-й минуте) полета.
Результаты расследования авиационного инцидента показали, что одним из членов экипажа была допущена ошибка: вместо переключения на резервный регулятор соотношения давления был включен тумблер клапана выпускного (типа Кл. 6894) для принудительной разгерметизации салона. Эта ошибка была распознана лишь через 12-14 минут. Экипаж все это время в соответствии с руководством по летной эксплуатации находился в кислородных масках.
Второй ошибкой экипажа явилось то, что самолет снизился лишь до высоты 8 000-7 500 м, а не до 4 000 м, как положено в таких ситуациях. При этом пассажиры (52 человека) и стюарды (6 человек) находились в разгерметизированном салоне без дополнительного кислородного обеспечения (точка кислородного питания была только в помещении для стюардов).
Как свидетельствуют пассажиры и (неофициально) стюарды, практически все пассажиры в той или иной степени выраженности потеряли сознание. К счастью, за этот относительно непродолжительный период никто из них не погиб.
0,5 - ф
0,5 + ф
Традиционный вариант расчета риска потери сознания у пассажиров при постоянно изменяющейся высоте в салоне в течение 25 минут пребывания на высотах от 3 000 до 8 700 м не позволил дать даже ориентировочную оценку вероятности сохранения сознания, хотя суммарное время пребывания на высоте более 7 500 м приближалось к максимально возможным значениям резервного времени.
Использование предлагаемого подхода при ретроспективной оценке инцидента позволило по предложенным математическим моделям рассчитать дозу гипоксического воздействия и получить количественную оценку риска потери сознания.
Результаты расчета вероятности потери сознания пассажирами (92 %) хорошо согласуются с материалами расследования авиационного инцидента. Из 58 человек, находившихся в салоне без дополнительного кислородного обеспечения, только один сохранил относительную работоспособность и сумел воспользоваться кислородной маской. То есть фактически 98 % людей, находившихся в пассажирском салоне воздушного судна, потеряли сознание.
Результаты позволяют сделать вывод о высокой достоверности получаемых при использовании разработанных на основе дозового принципа нормирования воздействия на человека гипокси-ческой газовой среды математических моделей оценок риска снижения и утраты работоспособности человека в условиях выраженных степеней острого гипоксического состояния, прогнозных значений резервного времени наступления потери сознания и потребного времени устранения дефицита кислорода при снижении воздушного судна до высоты, входящей в зону полной компенсации, что весьма важно для разработки систем жизнеобеспечения и систем обеспечения жизнедеятельности летательных аппаратов [1, 4, 5].
В заключение следует отметить, что математическая модель дозы гипоксии в виде интеграла интенсивности гипоксического воздействия в пределах времени действия этого фактора позволяет переходить в динамическую область исследования эффектов воздействия гипоксической газовой среды и оценивать ее влияние на состояние человека
при изменении интенсивности гипоксического воздействия во времени. Разработанный на основе математической модели метод расчета риска потери работоспособности человеком в условиях низкого барометрического давления обеспечивает поддержку принятия решений при выборе оптимальной траектории полета летательных аппаратов в случае разгерметизации кабины или салона, при оказании помощи пораженным в чрезвычайных ситуациях в условиях высокогорья и в ряде других практических задач, связанных с пребыванием человека в условиях гипоксической гипоксии.
Литература
1. Авиационная медицина: руководство; [под ред. Н.М. Рудного, П.В. Васильева, С.А. Гозулова]. М.: Медицина, 1986. 578 с.
2. Кукушкин Ю.А., Богомолов А.В, Дворников М.В., Кисляков Ю.Ю., Рыженков С.П. Расчет риска потери работоспособности человеком в условиях низкого барометрического давления // Полет. 2012. № 11. С. 37-45.
3. Кукушкин Ю.А., Богомолов А.В., Дворников М.В., Степанов В.К., Сухолитко В.А. Методика определения потенциальной ненадежности действий и резервного времени сохранения работоспособности летчика в высотном полете // Проблемы безопасности полетов. 2002. № 11. С. 22-27.
4. Ушаков И.Б., Черняков И.Н., Шишов А.А. Физиология высотного полета. Воронеж: Истоки, 2008. 147 с.
5. Федоров М.В., Богомолов А.В., Цыганок Г.В., Айвазян С.А. Технология планирования многофакторных экспериментальных исследований и построения эмпирических моделей комбинированных воздействий факторов на операторов эрга-тических систем // Информационно-измерительные и управляющие системы. 2010. № 5. С. 53-61.
References
1. Rudny N.M., Vasilyev P.V., Gozulov S.A., Aviatsionnaya meditsina: rukovodstvo [Aviation Medicine: manual], Moscow, Meditsina, 1986, 578 р.
2. Kukushkin Yu.A., Bogomolov A.V., Dvornikov M.V., Kislyakov Yu.Yu., Ryzhenkov S.P., Polyot [Flying], 2012, no. 11, pp. 37-45.
3. Kukushkin Yu.A., Bogomolov A.V., Dvornikov M.V., Stepanov V.K., Sukholitko V.A., Problemy bezopasnosti polyotov [Problems of flight safety], 2002, no. 11, pp. 22-27.
4. Ushakov I.B., Chernyakov I.N., Shishov A.A., Fiziologiya vysotnogo polyota [Physiology of high-altitude flight], Voronezh, Istoki, 2008, 147 p.
5. Fedorov M.V., Bogomolov A.V., Tsyganok G.V., Ayvazyan S.A., Informatsionno-izmeritelnye i upravlyayushchie sistemy [Information-measuring and control systems], 2010, no. 5, pp. 53-61.
УДК 62-503.5
ПРИМЕНИМОСТЬ ПРОГРАММНЫХ КОМПЛЕКСОВ ДЛЯ РАБОТ СО СНИМКАМИ К ЗАДАЧЕ АНАЛИЗА РЕНТГЕНОГРАММ
М.В. Овечкин, аспирант; А.И. Сердюк, д.т.н., профессор, директор АКИ (Оренбургский государственный университет, просп. Победы, 13, г. Оренбург, 460018, Россия,
[email protected], [email protected])
Разнообразные системы диагностики технических объектов и управления технологическими процессами широко распространены в различных областях промышленности. Контроль и диагностика в большинстве систем оптическо-
