Биологические
науки/11.Биоинженерия и
биоинформатика
Асп. Индюхин А.А., к.т.н. Коржук Н.Л., к.б.н. Индюхин А.Ф.
Тульский государственный университет, Россия
Метод динамического анализа ЭЭГ в неврологической
диагностике
Анализ электроэнцефалограмм (ЭЭГ) в неврологии
длительное время был ориентирован прежде всего на выявление грубой патологии. Только
к концу 20 века возобладало очевидное мнение об эффективности ЭЭГ как
инструмента исследования здорового мозга и центральной нервной системы (ЦНС). В
то же время многие вопросы ЭЭГ-диагностики негрубых нарушений ЦНС и сенсорных
нарушений остаются по-прежнему актуальными [1]. Эффективным диагностическим
средством зарекомендовали себя функции когерентности. Их использование
опирается на теорию функциональных систем (ФС) П.К. Анохина. Управляющие
элементы ФС сосредоточены в коре головного мозга, и их функционирование
сопровождается вовлечением в синхронную, когерентную деятельность различных
нейронных ансамблей. Функции когерентности (ФК) дают численную оценку степени
взаимодействия различных участков коры головного мозга в частотной области. Расчет
ФК базируется на спектральных преобразованиях, прежде всего на быстром
преобразовании Фурье. Оно может быть осуществлено только после того, как
регистрация определенной эпохи ЭЭГ закончена, и имеется определенное количество
точек исходного сигнала. Еще более затруднительным моментом является
необходимость визуального анализа эпохи для оценки качества и «типичности»
записи для конкретного пациента. Сам по себе расчет ФК математически очень
сложен, и может быть выполнен только с помощью ЭВМ. Обычно опция расчета ФК
входит в пакет программного обеспечения компьютерных электроэнцефалографов, и
для ее осуществления необходима достаточно высокая квалификация пользователя.
В нашей работе мы предлагаем полную
альтернативу расчету ФК – новый метод анализа ЭЭГ, позволяющий в реальном
времени оценить уровень взаимодействия тех участков коры, с которых
регистрируется ЭЭГ. Метод получил название динамический
анализ ЭЭГ. Он состоит в том, что сигнал ЭЭГ подается на цифровой
самонастраивающийся фильтр, который выделяет сигнал доминирующего ритма.
Настройка осуществляется известным из радиотехники методом ФАПЧ (фазовая
автоподстройка частоты). На вход регулятора подаются входной и выходной сигнал
фильтра (рис. 1).
Рис. 1. Блок-схема самонастраивающегося фильтра
Регулятор измеряет разницу фаз сигналов и вырабатывает
управляющий сигнал, изменяющий резонансную частоту полосового фильтра. По
окончании настройки доминирующий ритм проходит сквозь фильтр без фазовых и
амплитудных искажений. В результате в любой момент времени известно мгновенное
значение сигнала доминирующего ритма и его частота.
Уравнения полосового
фильтра:
UФ
= К1(UР)UВ + К2(UР)UВz-2 - К3(UР)UФz-1 - К4(UР)UФz-2;
где: UФ – выходной сигнал полосового фильтра;
UВ –
входной сигнал полосового фильтра;
UР –
выходной сигнал регулятора;
К1, К2, К3,
К4 – коэффициенты рекуррентного уравнения;
ωР – резонансная
частота полосового фильтра;
ωП - полоса
пропускания;
τ – шаг дискретизации.
Применение динамического анализа (ДА) позволяет,
не проводя спектральных преобразований, определить частоту доминирующего ритма.
Проведенные исследования показали, что частота, амплитуда и длина веретена
доминирующего ритма (в ЭЭГ чаще всего это альфа-ритм) являются ценными
диагностическими показателями в онтогенетическом плане [1].
Другим приложением ДА оказался сравнительно
простой алгоритм автоматической диагностики эпилепсии [2]. При прохождении
через полосовой фильтр комплекса «острая – медленная волна» резонансная частота
фильтра изменяется по визуально сходному закону, что позволяет корреляционным
методом обнаружить искомый комплекс.
Систему ДА можно использовать в компенсационном
режиме: вычитая из исходного сигнала ЭЭГ сигнал доминирующего ритма, визуализировать
те низкоамплитудные составляющие, которые незаметны на фоне более
высокоамплитудной активности. В частности, таким образом можно регистрировать
единичные реализации когнитивных вызванных потенциалов Р300 [3]. Поставив перед
испытуемым задачу «проговаривать» слова, соответствующие подаваемому стимулу,
можно проследить формирование команды в моторном центре речи, т. е. определить
не только электрофизиологические, но и психофизиологические характеристики пациента.
Упомянутые ФК доказали свою диагностическую
эффективность и прогностическую ценность в клинике критических состояний –
комы, тяжелых черепно-мозговых травм, опухолей мозга, у недоношенных детей.
Меньшее распространение получили исследования здоровых детей и детей с
неврологическими нарушениями.
Для оценки ФК используются их средние уровни в
стандартных частотных диапазонах и частоты, на которых ФК достигает максимума.
Аналогичные показатели предложены и в рассматриваемом способе ДА: уровни и
частоты синхронизации [4].
Рассмотрим структуру системы, осуществляющей расчет
частот и уровней синхронизации (рис. 2).
Рис. 2. Структурная схема устройства определения
уровня k1 и частоты f1 синхронизации
для одного частотного диапазона. ФД1 – фильтры диапазона, выходные сигналы
которых подаются на блок произведений; ОС1 - осреднитель; ФП1– фильтр
переменной частоты настройки.
Выходные
сигналы усилителя ЭЭГ попарно поступают на одинаковые фильтры диапазонов ФДi, после которых сигналы идут на блоки произведений, на
выходе которых при наличии в сигналах общих синусоидальных составляющих будет
присутствовать некоторый уровень постоянного сигнала. Количество анализируемых
диапазонов может соответствовать стандартным шести (дельта, тета-1, тета-2,
альфа, бета-1, бета-2), либо быть каким-либо другим.
Рассмотрены четыре группы детей в возрасте от 4
до 10 лет: «Норма» - 29 человек, «Минимальная мозговая дисфункция» - 33 человека,
«Нарушения зрения» - 19 человек, «Нарушения слуха» - 22 человека. Для всех были
получены уровни и частоты синхронизации в диапазонах дельта, тета, альфа, и бета по 120 парам отведений
(для 16 каналов регистрации). Проведено сравнение групп с неврологическими
нарушениями с группой «Норма» по всем 960 параметром и точным методом Фишера
получена достоверность отличий. Результаты приведены в таблице.
Как видно из таблицы, число отличий по частоте
синхронизации превосходит аналогичный показатель по уровню синхронизации.
На рис. 3 показана топография
выявленных отличий в тета-диапазоне. Несмотря на некоторые совпадения между
выявленными парами, видно, что идентификация конкретного пациента по группам
нарушений не вызовет больших трудностей.
Таблица
Количество маркеров, по которым группы
детей с неврологическими нарушениями достоверно (р < 0,05) отличаются от
здоровых по уровню и частоте синхронизации
Тип неврологического нарушения |
Нарушение зрения |
Нарушение слуха |
Минимальная мозговая дисфункция |
|||||||||
Частотный диапазон |
δ |
θ |
α |
β |
δ |
θ |
α |
β |
δ |
θ |
α |
β |
Уровень синхронизации |
5 |
6 |
6 |
8 |
5 |
5 |
2 |
6 |
0 |
8 |
6 |
12 |
Частота синхронизации |
1 |
16 |
32 |
4 |
1 |
10 |
24 |
10 |
9 |
18 |
4 |
2 |
а)
б)
с)
Рис. 3. Распределение маркеров достоверных отличий по
частоте синхронизации в тета-диапазоне в группах с неврологическими нарушениями
по отношению к норме; а) минимальная мозговая дисфункция; б) нарушения зрения;
в) нарушения слуха.
С использованием
всех перечисленных возможностей ДА разрабатывается высокотехнологичный
диагностический комплекс экспресс-оценки неврологического и
психофизиологического состояния пациента в виде портативного прибора с
минимальным временем обследования и выдачей заключения на миниатюрный экран. Комплекс
соответствует по своей реализации востребованному на современном этапе процессу
проведения комплексного изучения состояния здоровья и динамики
психофизиологического развития подрастающего поколения.
Литература:
1.
Жеребцова,
В.А. Системный анализ механизмов организации высших психических функций в
онтогенезе [Текст] / Жеребцова Валентина Александровна. Дисс. … докт. биол.
наук. – Тула, 2004.
2.
Мишина,
Л.М. Подсистема электроэнцефалографической диагностики эпилепсии [Текст] / Л.М.
Мишина, А.А. Индюхин, А.Ф. Индюхин, Н.С. Тархов // Медицинская техника, 2012. -
№ 1 (271). – С. 22 – 25.
3.
Жеребцова,
В.А. Способ диагностики длиннолатентного вызванного потенциала мозга и устройство
для его осуществления / В.А. Жеребцова., А.Ф. Индюхин, Э.М. Соколов, А.А. Хадарцев, В.П. Васильев, В.И.
Морозов // Патент РФ № 2240036 от 20.11.2004 г.
4.
Коржук, Н.Л. Способ
электроэнцефалографической диагностики неврологических нарушений [Текст] / Н.Л.
Коржук, А.А. Индюхин, А.Ф. Индюхин, В.В. Савельев // Известия ЮФУ. Технические науки. Тематический
выпуск: «Медицинские информационные системы». Таганрог: Изд-во ТТИ ЮФУ, 2010. -
№ 8 (109). – С. 121 – 127.