Синергетика_Глава_3_10_Математические_методы_анализа_ритмограмм
.pdf3.10 Математические методы анализа ритмограмм
3.10.1. Введение
Значение клинических исследований вариабельности сердечного ритма (ВСР), выполненных начиная с 80-х годов прошлого века, значительно возросло, когда была показана прямая и независимая от других факторов связь ВСР со смертностью от острого инфаркта миокарда. Исследованиями, выполненными в эти годы, была установлена тесная корреляция нарушений автономной нервной регуляции с жизнеопасными желудочковыми аритмиями (ЖА) и сердечной смертью, включая внезапную. Риск последних резко возрастает при существенном увеличении симпатической и падении парасимпатической активности на фоне уменьшения мощности RR-спектра (равной величине дисперсии RR-интервалов). Эти параметры легко получить, анализируя энергетический спектр сердечного ритма. Более глубокая оценка состояния ВНС, взаимоотношения вегетативной нервной системы и ВСР, стала возможной в конце прошлого века в связи с широким внедрением в клиническую практику методов спектрального анализа сердечного ритма . Начало 90-х годов ознаменовалось существенным прогрессом метода, что было связано с внедрением новых цифровых многоканальных устройств 24часовой записи ЭКГ, а также с усовершенствованием методик математической обработки кардиограмм .
Сегодня к проблемам анализа причин возникновения и характера ВСР начали применять методы теории хаоса, нелинейных динамических систем и математического моделирования. Хотя эти методы еще не вышли из стен лабораторий, они являются многообещающими.
Области применения технологии вариабельности сердечного ритма - валеология, терапевтическая клиника в ее широком понимании, включая кардиологию, пульмонологию, нефрологию, неврологию, психиатрию, а также акушерство, перинатология, реаниматология и др. области клинической медицины. В развитых странах исследование ВСР включено в стандартные протоколы контроля над состоянием здоровья при профилактических осмотрах, наблюдении за беременными, ведении послеинфарктных больных.
3.10.2. Математические методы анализа вариабельности сердечного ритма. Ритмограмма
В середине прошлого века фундаментальные исследования в области физики и математики доказали, что многие явления окружающего нас мира, которые считались случайными и хаотическими подчиняются строгим закономерностям.
246
Это было настолько поразительным, что получило название «детерминированного хаоса». Явления детерминированного хаоса наблюдаются в сложных системах, имеющих большое количество функциональных связей. Одним из примеров такой сложной системы является сердечно-сосудистая система. На сердце суммируются все управляющие импульсы, поступающие от центральных структур головного мозга, симпатического и парасимпатического отделов вегетативной нервной системы и отражающие изменения внутренней и внешней среды. Работа сердца как исполнительного органа позволяет в определенной степени судить о работе всей системы в целом. Ритм сердца имеет сложную динамическую структуру, образованную суперпозицией многих управляющих воздействий, обусловленных, как принято считать, регулирующим влиянием центральной нервной системы, вегетативной нервной системы и нейрогуморальных воздействий. Такие широко распространенные методы исследования электрической оценки сердца как суточное мониторирование и электрокардиограмма высокого разрешения в настоящее время обращают особое внимание на анализ флуктуационных характеристик сердечного ритма. В последнее время для интерпретации сердечных ритмов применяется математический формализм фрактальной математики, использующий многие подходы нелинейной динамики.
Интерес к методу, исследующему вариабельность сердечного ритма, резко увеличился после сообщения о том, что ослабление колебаний частоты сердечных сокращений при остром инфаркте миокарда обладает высокой прогностической значимостью и является независимым указателем, свидетельствующим о большой вероятности наступления внезапной смерти [WV],[BFRS]. В этих работах было продемонстрировано, что при снижении общей мощности спектрального состава вариабельности сердечного ритма значительно возрастает опасность смертного исхода. Точно такая же закономерность прослеживалась и в отношении временных показателей. При значениях стандартного отклонения длительности RR интервалов ниже 50 мс смертность была в 5.3 раза выше, чем при значениях более 100 мс.
Для стандартизации получаемых данных и использования математического обеспечения в 1995 году специальная комиссия Европейского кардиологического общества и Американского общества электростимуляции и электрофизиологии выработали стандарты измерения, физиологической интерпретации и клинического применения вариабельности сердечного ритма. Отечественные методические рекомендации анализа вариабельности сердечного ритма при использовании различных электрокардиографических систем, подготовленные в соответствии с решением Комиссии по клиникодиагностическим приборам и аппаратам Комитета по новой медицинской
247
технике МЗ РФ (протокол N 4 от 11 апреля 2000 г), приведены в работе
[БИЧ].
Основной метод исследования вариабельности сердечного ритма |
|
основан на распознавании |
и обработке временных интервалов между RR |
зубцами (RR интервалы) |
электрокардиограммы, построении динамических |
рядов кардиоинтервалов и последующего анализа полученных данных различными математическими методами. Динамический ряд кардиоинтервалов может быть как стационарным так и нестационарным. Стационарные процессы характеризуются эргодичностью, когда усреднение случайного сигнала по времени совпадает с усреднением по множеству реализаций. Иными словами на любом участке времени мы должны получать одни и те же характеристики сердечного ритма. Нестационарные (или переходные) процессы имеют определенную тенденцию развития во времени и их характеристики зависят от в начала отсчета. Практически в каждой кардиограмме содержатся элементы нестационарности. Вначале предположим, что интервалограмма сердечного ритма является
стационарным случайным процессом и дадим |
соответствующую |
интерпретацию экспериментов. |
|
При анализе динамических рядов кардионтервалов следует различать кратковременные ('короткие') и долговременные ('длинные') записи. Под последними, как правило, понимают данные, получаемые при 24-х и 48 часовом мониторировании ЭКГ (Холтеровское мониторирование). К так называемым 'коротким' записям относят данные исследований, проводимых в течение минут, десятков минут или нескольких часов. Рабочей Группой анализ ВСР рекомендуется проводить на двух временных промежутках - 5- минутном и 24-часовом. Их предлагают называть кратковременным и долговременным. В кратковременных измерениях показатели ВСР ниже, чем в долговременных.
Во время изучения результатов кардиограммы в большом числе случаев приходится сталкиваться с различными артефактами записи. В большинстве случаев это достаточно хорошо известные шумы и наводки, определение источника которых не представляет каких-либо сложностей. Основную проблему при анализе кардиограммы представляют артефакты, имитирующие нарушения ритма сердца или псевдоаритмии. Для максимального устранения таких артефактов следует обратить внимание на математический алгоритм определения положения основных зубцов сердечного ритма. К сожалению, такие псевдоаритмии могут серьезно повлиять на результаты исследования и в конечном итоге на постановку правильного диагноза. В любом случае, при наличии сомнений врача в
248
истинности возникающих аритмий, необходимо их описание с указанием возможного артефактного характера и, если возможная псевдооаритмия может существенно повлиять на окончательный диагноз и прогноз больного, необходимо рекомендовать повторное обследование. Во многих случаях правильная оценка возможна только при одновременной интерпретации всех регистрируемых каналов электрокардиограммы.
Особое внимание при анализе RR ритма следует экстрасистолам – преждевременным сокращениям сердца, обусловленным возникновением импульса вне синусового узла. Экстрасистолы и другие аритмии, дефекты записи, ее зашумленность могут вносить изменения в оценку параметров вариабельности сердечного ритма. Конечно, предпочтительнее использовать короткие записи без экстрасистол, и шумов, однако, при некоторых обстоятельствах подобная избирательность, может приводить к необъективности. Необходимо помнить, что полученные результаты могут зависеть от наличия экстрасистол, поэтому, приводя результаты исследований следует указывать число и относительную длительность интерполированных или выброшенных из обработки RR интервалов.
Все множество методов исследования ВСР можно представить как:
1. статистические методы;
2. анализ гистограмм;
3. спектральные методы;
4. нелинейные методы.
Суть пространственно - временных методов вытекает из их названия. Анализу подвергается участок ЭКГ, который образован повторяющимися PQRST комплексами, возникающими в синусовом узле. На выбранном участке ЭКГ определяются длины RR-интервалов. Чтобы подчеркнуть их принадлежность к синусовым сердечным комплексам, их часто называют NNинтервалами. На рисунке S.1 представлен участок ЭКГ с последовательностью комплексов синусового происхождения. По этой последовательности на декартовой плоскости строят интервалограмму. В ней по оси абсцисс откладывают номера RR-интервалов, а по оси ординат - их длины, таким образом интервалограмма отображает зависимость длительности RR от номера цикла измерения. Интервалограмма дает хорошее представление об изменениях сердечного ритма и у специалиста со временем накапливается опыт качественной оценки ВСР на основе этого метода.
249
Рис. S1. Формирование интервалограммы. На выбранном участке ЭКГ определяются промежутки между соседними R - R пиками сердечного ритма. Такие промежутки называются RRn интервалами, которые зависят от номера удара n 0,1,2...N 1, где N - полное количество ударов сердца.
3.10.3. Статистические методы оценки вариабельности сердечного ритма
Временной анализ (time domain) относится к группе методов оценки вариабельности сердечного ритма, основанных на применении статистических программ к обработке выделенного количества RR интервалов между R зубцами кардиограммы с последующей физиологической и клинической интерпретацией полученных результатов. Пусть мы имеем динамический ряд интервалов RRi , где индекс i=0,1,2,3… N–1 характеризует множество событий, а величина N представляет собой суммарное количество измерений RR интервалов. Как правило, величина интервала между R зубцами кардиограммы измеряется в миллисекундах. При статистическом анализе вариабельности сердечного ритма используются следующие величины:
RRNN (мс) – средняя длительность RR интервалов. Показатель RRNN отражает конечный результат многочисленных регуляторных влияний на синусовый ритм сложившегося баланса между парасимпатическим и симпатическим отделами вегетативной нервной системы;
|
1 |
N 1 |
|
RRNN |
|
RRi |
(S.1). |
|
|||
|
N i 0 |
|
|
250
Заметим, что величина |
RRNN в точности совпадает со средним значением |
||
RR интервалов, индекс |
i изменяется в интервале |
i 0,1,...N 1, где |
N - |
полное количество ударов. Количество ударов |
сердца N 300 , |
что |
|
приблизительно соответствует продолжительности снятия кардиограммы, равной 5 мин.
|
1 |
|
N 1 |
|
||
RR |
|
|
RRi |
(S.2). |
||
|
|
|||||
|
N i 0 |
|
||||
Аналогичным образом введем среднее значение <RR2> |
||||||
|
|
|
|
1 |
N 1 |
|
RR2 |
|
RRi2 |
(S.3) |
|||
|
|
|||||
|
|
|
|
N i 0 |
|
|
и выборочную дисперсию величины RR интервалов |
||||||
( R)2 RR2 RR 2 |
(S.4). |
|||||
HR (ЧСС) – частота сердечных сокращений. Этот параметр, измеряемый в 1/мин иногда называемый как частота пульса, характеризует средний уровень функционирования системы кровообращения.
HR |
60 *1000 |
(S.5). |
|
RRNN |
|||
|
|
Отметим, что величина RRNN при подсчете (ЧСС) измеряется в мс.
SDNN (standart deviation of the NN interval) – стандартное отклонение величин нормальных RR-интервалов
|
1 |
N 1 |
2 |
1/ 2 |
|
SDNN |
|
RRi |
RRNN |
|
(S.6). |
|
|||||
|
N 1 i 0 |
|
|
|
|
Эта величина, представляющая собой стандартное отклонение σ=SDNN, является несмещенной оценкой дисперсии D=σ2 случайной величины RRi . Естественно, что при больших значениях N выполняется соотношение
SDNN |
( R)2 |
|
(S.7). |
|
251 |
|
|
Параметр SDNN является интегральным показателем, характеризующим вариабельность ритма сердца в целом, и зависит от влияния на синусовый узел симпатического и парасимпатического отделов ВНС. Увеличение или уменьшение этого показателя свидетельствует о смещении вегетативного баланса в сторону преобладания одного из отделов ВНС.
Отметим, что для холтеровского миниторирования, когда наблюдение за ритмом сердца ведется в течение большого интервала времени (например, в течение суток) вводятся величины SDNN-i (мс) средние значения стандартных отклонений за 5 минут, а также SDANN-i стандартное отклонение усредненных за 5 минут значений интервалов RR. Нормальные значения SDNN находятся в пределах 40-80 мс. Считается, что при анализе коротких записей рост SDNN указывает на усиление контура автономной регуляции. Уменьшение SDNN связано с усилением симпатической регуляции, которая подавляет активность автономного контура. Резкое снижение этой величины обусловлено значительным напряжением регуляторных систем, когда в процесс регуляции включаются высшие уровни управления, что ведет к почти полному подавлению активности автономного контура. Значения этого показателя могут зависеть от многих факторов, в частности, от основного заболевания.
По литературным данным, приведенным в работе [Сне] у больных артериальной гипертензией SDNN снижен - 38,1±4,1 мс, у больных ИБС со стабильной стенокардией напряжения показатель SDNN, был 32,8±3,9 мс, с острым коронарным синдромом - 39,3±9,8 мс. Среднее значение SDNN у здоровых людей составляет: до 25 лет - 70±10 мс, 26-40 лет 60,1±6 мс, старше 40 лет - у мужчин 60±8 мс, у женщин - 50±4 мс. По мнению ряда авторов, величина SDNN зависит также от возраста и частоты сердечных сокращений.
RMSSD – среднеквадратичная разностная характеристика (root mean sum successful deviation). Эта величина, измеряемая в мс, представляет собой
квадратный корень |
из среднего квадрата разностей величин |
||||||
последовательных пар RR интервалов. |
|
|
|
|
|||
|
1 |
N 1 |
|
|
2 |
1/ 2 |
|
RMSSD |
|
|
RR RR |
|
|
(S.8). |
|
|
|||||||
|
(N 1) |
i |
i 1 |
|
|
|
|
|
i 1 |
|
|
|
|
|
|
Сущность показателя RMSSD состоит в оценке степени различия двух соседних интервалов RR. То есть при гипотетическом отсутствии колебаний соседний RR интервалов (так называемый «вкопанный» ритм, возможном при фиксированной частоте сердечных сокращений, навязанной кардиостимулятором) или высокой степени вегетативной блокады сердца
252
данный показатель стремится к нулю. Чем больше между собой разница соседних интервалов (то есть чем выше синусовая аритмия), тем будут более высокие значения RMSSD. Естественно, что при расчете этого показателя не учитываются медленные составляющие сердечного ритма. Чем выше значения RMSSD, тем активнее звено парасимпатической регуляции. Считается, что показатель RMSSD, тестирует способность синусового узла к концентрации ритма сердца, регулируемой переходом роли основного водителя ритма к различным отделам этого узла, с неодинаковым уровнем синхронизации возбудимости и автоматизма. На фоне усиления симпатических влияний, отмечается уменьшение показателя RMSSD, то есть усиление концентрации ритма и наоборот, при активации тонуса вагуса, функция концентрации ритма сердца снижается [Мак2]. В норме значение этого показателя находится в переделах 20-50 мс. В работе [Сне] получено среднее значение RMSSD у здоровых людей до 25 лет 49±15,23 мс. У больных артериальной гипертензией отмечено снижение RMSSD - 25,2±3,3 мс .
3.10.4. Анализ гистограммы сердечного ритма
Под гистограммой понимается графическое изображение сгруппированных значений сердечных интервалов, при котором по оси абсцисс откладываются
значения RR интервалов [RRmin, RRmax], а |
по оси ординат - их количество |
|||||
(частота встречаемости), |
то есть число попаданий RR интервалов в данный |
|||||
интервал группировки. |
Сущность методов, исследующих закон |
|||||
распределения |
кардиоинтервалов, состоит |
в построении |
гистограммы, |
|||
показывающей |
кривую |
распределения |
кардиоинтервалов. |
Иногда эту |
||
гистограмму |
называют |
вариационной |
пульсограммой |
и |
целью |
|
исследования является непосредственное описание и интерпретации формы
этой кривой. |
Основные |
методы изучения |
такой |
гистограммы |
распределения принадлежат школе профессора Р.М.Баевского |
[БКК]. |
|||
Различают следующие типы гистограмм распределения ритма сердца:
нормальная гистограмма, близкая по виду к кривым Гаусса, типична для здоровых людей в состоянии покоя;
асимметричная - указывает на нарушение стационарности процесса, наблюдается при переходных состояниях.
Пусть для множества наблюдений случайной величины {RR0, RR1, …RRN-1} известны RRmin и Rmax – минимальное и максимальное значения выборки, а величина R=RRmax–Rmin представляет собой размах переменной
253
RRn. Наглядное описание представленных экспериментальных данных осуществляется с помощью группировки наблюдений в классы.
Под группировкой или классификацией мы будем понимать некоторое разбиение всего интервала R=RRmax–RRmin , содержащего N наблюдаемых результатов наблюдений {RR0, RR1, …RRN-1} на k интервалов. Для простоты мы будем считать, что все интервалы j=1,2,…k, которые мы будем называть интервалами группировки, имеют одинаковую длину и покрывают всю наблюдаемую область. Пусть n j - это число наблюдений в j интервале
группировки. Ясно, что суммарное количество наблюдений во всех интервалах группировки должно быть равно полному числу наблюдений
k
n j N .
j 1
Определим величину fj − относительную частоту попаданий наблюдений {RR0, RR1, …RRN-1} в j интервал группировки
f j n j N
Из определения видно, что попадания наблюдений в нормировки
для величин fj , имеющих смысл вероятности j интервал группировки, справедливо условие
k
f j 1 .
j 1
Графическое изображение зависимости частоты попадания элементов выборки от соответствующего интервала группировки называется гистограммой выборки. По определению высота каждого столбца гистограммы пропорциональна относительной частоте попадания наблюдений в данный интервал группировки.
Естественно, что величина интервала группировки , или количеcтво
интервалов k, существенно влияют на общий вид гистограммы. Обычно
предполагается, что величина k лежит в интервалах 5<k<20. Во многих
учебниках по статистике число интервалов разбиения k определяется с
помощью формулы Стерджесса [Пол]
k Int[1 log 2 (N)] Int[1 1,442ln( N)] Int[1 3,322lg( N)] . |
(4.7) |
В формуле (4.7), значок Int означает целую часть числа где log 2 , |
ln , log – |
логарифм по основанию 2, натуральный и десятичный логарифм. Из формулы (4.7) следует, что увеличение числа точек в сигнале N в два раза приведет к тому, что число интервалов разбиения k увеличится на единицу.
254
Если длина группировки мала, то влияние случайных колебаний начинает преобладать, так как каждый интервал содержит при этом лишь небольшое количество наблюдений. В этом случае значения гистограммы fj будут резко изменять свои значение при изменении аргумента j на единицу. Чем больше величина интервала группировки , тем более размазанными являются характерные черты гистограммы распределения.
В многолетней практике сложился традиционный подход к группировке кардиоинтервалов в диапазоне от 400 до 1300 мс с интервалом в 50 мс. Таким образом, выделяются 20 фиксированных диапазонов длительностей кардиоинтервалов, что позволяет сравнивать вариационные пульсограммы, полученные разными исследователями на разных группах исследований. При этом объем выборки, в которой производится группировка и построение вариационной пульсограммы, также стандартный - 5 минут. Другой способ построения вариационных пульсограмм заключается в том, чтобы вначале определить модальное значение кардиоинтервала, а затем, используя диапазоны по 50 мс, формировать гистограмму в обе стороны от моды. При использовании вариационной пульсометрии выделяются следующие основные параметры.
Рис. S.2. Ритмограмма RRn в зависимости от числа ударов n. Величина RRn измеряется в миллисекундах (ms).
Для построения гистограммы fj будем считать, что весь интервал изменения [RRmin; RRmax ] мы разобьем на интервалы группировки величиной ∆=50 мс.
255
