Ключевые слова дисперсия интервала QT, дисперсия интервала QTa, модели генеза электрокардиосигналов, дисперсия интервала QTa | Key words QT-interval dispersion, QT-interval dispersion, models of the electrocardiosignal genesis, QTa-interval dispersion |
Аннотация Вопросы генеза, измерения и информативности дисперсии интервала QTa (от начала комплекса QRS до вершины зубца T) рассмотрены в терминах биофизических моделей электрической активности сердца. | Annotation Problems of genesis, measurement, and informativity of the QTa-interval (from the QRS onset to the T wave apex) have been considered in terms of biophysical models of the cardiac electrical activity. |
Автор Баум, О. В., Попов, Л. А., Волошин, В. И., Муромцева, Г. А. | Номера и рубрики ВА-N29 от 15/12/2002, стр. 10-18 /.. Оригинальные исследования |
Презумпция информативности дисперсии QT-интервала (QTD) как предиктора возможного развития желудочковых аритмий поддерживается отсутствием общепринятых доказательств, опровергающих гипотезу о связи пространственной негомогенности процесса реполяризации в миокарде желудочков с распределением QT-интервала по длительностям, наблюдаемым в условиях реальных электрокардиографических измерений в разных отведениях. Нет, впрочем, единого биофизически обоснованного мнения и в поддержку этой гипотезы. Широкий спектр мнений по данной проблеме [1-9] объясняется во многом сложностями идентификации окончания зубца T и, следовательно, измерения длительности QT-интервала [8, 10-14].
По итогам теоретического рассмотрения проблем QTD и результатам проведенных компьютерных экспериментов на моделях электрической активности сердца [8, 14, 15] «тайну» QT дисперсии можно трактовать следующим образом. С одной стороны, наблюдаемое различие измеряемых длительностей QT интервала от отведения к отведению, tQT(l), можно объяснить в общем случае только используемыми алгоритмами определения конца зубца T, способами проведения изолинии, введением порогов идентификации в связи с принятыми правилами кодирования или с реальными значениями отношения сигнал/шум, а также ошибками измерения.
С другой стороны, как это ни парадоксально, именно измерения длительности QT-интервала с некоторым порогом делают, в принципе, возможным использование такого параметра как QTD в качестве показателя существующей негомогенности длительности трансмембранного потенциала (ТМП) на поверхности эпикарда и эндокарда [8]. В то же время, результаты модельного картирования длительности QT-интервала на поверхности виртуальной сферы отображения методом SSPM (Standard surfase potential method) [16] не позволяют на данный момент однозначно трактовать (вероятно, в силу упомянутого субъективизма измерений QT) наблюдаемую связь между распределением длительностей ТМП на поверхности электрически активного миокарда модели и распределением tQT(y, q), где y и q - угловые координаты поверхности сферы отображения, в терминах количественных моделей.
Наряду с измерениями tQT, которые имеют самостоятельное значение для определения синдрома удлинения QT-интервала [17, 18], и вычислением QTD, в качестве параметров, характеризующих пространственную негомогенность процесса реполяризации желудочков, используются и другие, менее изученные индексы [19-21], среди которых следует отметить дисперсию интервала QTa.
Интервал QTa измеряется между началом комплекса QRS и временем достижения зубцом T своего пикового значения в каждом из отведений выбранной системы. На множестве измеренных значений tQTa(l) ищутся максимальное и минимальное значения длительности интервала QTa. И, аналогично вычислению QTD, определяется QTa-дисперсия как QTaD = tQTamax - tQTamin .
Этот индекс, внешне привлекательный кажущейся простотой процедуры определения вершины зубца T, также как и некоторые другие новые индексы, начинают использовать во многих исследованиях [22-25].
Интересно отметить, что одна из первых публикаций [26], где для удобства измерений интервал QTa использовался в качестве характеристики длительности электрической систолы, относится, как отмечается в классической работе H.C.Bazett [27], еще к началу прошлого века.
Учитывая сложности и субъективизм измерений длительности QT-интервала, а также имеющиеся сомнения по поводу наличия электрофизиологического смысла у параметра QTD [8], представляет интерес исследовать дисперсию QTa-интервала, QTaD, в качестве претендента на роль более объективного и воспроизводимого аналога QT-дисперсии или в качестве ее более информативной (по результатам сравнительной оценки) альтернативы.
Цель данной работы: 1) рассмотреть и проанализировать процесс формирования длительности QTa-интервала в различных отведениях и QTa-дисперсии, используя векторкардиографическую интерпретацию модели эквивалентного генератора сердца в виде одиночного диполя; 2) получить экспериментальные оценки величин tQTa(l) и QTaD в стандартной и ортогональной системах отведений для выборки практически здоровых субъектов; 3) провести с помощью биофизической модели электрической активности сердца компьютерные эксперименты по оценке влияния эпи-эндокардиального распределения длительности ТМП на распределение tQTa(l) и величину QTaD.
ПРЕДВАРИТЕЛЬНЫЙ АНАЛИЗ
Рассмотрим сначала, как и для QT [8], процесс формирования QTa-интервала в разных отведениях в терминах ВКГ-петель, соответствующих дипольной модели эквивалентного генератора сердца. При этом будем помнить, конечно, об ограничениях такой модели, связанных с отсутствием в сигнале вклада от мультипольных компонент выше первого порядка [8].
На рис. 1 схематически изображены в качестве примера проекция ВКГ-петли зубца T на некоторую плоскость, а также оси нескольких произвольных отведений (l1, l2, l3, l4). Вершина зубца T в любом отведении соответствует на петле точке, в которой нормаль к оси этого отведения касается петли (например, точка Ta 4 для отведения l 4 ). Конец моментного вектора ВКГ, обежав сначала петлю QRS за время tQRS и «потоптавшись» на месте за время почти изоэлектричного в норме сегмента ST, двинется далее по петле T в соответствие с направлением ее обхода и достигнет точки Ta 4 за время, которое от начала комплекса QRS как раз и составляет QTa-интервал в отведении l4 , tQTa(l4). В этой точке зубец T в l4 достигнет своего максимального значения.
Рис. 1. Схематическое изображение проекции петли T на некоторую плоскость; l1 - l4 -произвольные отведения, оси которых лежат в данной плоскости. Ta i - точка, в которой нормаль к отведению l i касается петли T; эта точка во временной области соответствует вершине зубца T в отведении l i. Дисперсия интервала QTa в системе отведений l 1 - l 4, QTaD, определяется, главным образом, временем, за которое конец моментного вектора T пройдет вдоль петли расстояние между крайними точками Ta (в данном случае - от Ta 4 до Ta 1).
Скорость движения моментной точки по петле T в общем случае не является константой и определяется скоростями изменения сигналов в двух отведениях, выбранных для задания петли T. Множество точек Ta(l) для всех рассматриваемых в этой плоскости отведений определяет множество интервалов QTa(l), а значения tQTa min и tQTa max , соответствующие первой и последней по ходу петли точкам Ta, показывают размах (дисперсию) QTa-интервала для отведений данной плоскости. Оси большинства систем отведений (стандартной, ортогональной) не лежат в одной плоскости, и это надо иметь в виду, желая получить наглядную картину генеза всего множества интервалов QTa в рассматриваемой системе отведений.
Практически же, для определения дисперсии QTa-интервала значения tQTa(l) измеряют непосредственно по записям, желательно синхронным, скалярной ЭКГ и далее рассчитывают их дисперсию как разность между tQTa max и tQTa min .
Таким образом, мы видим, что, в отличие от QTD, которая теоретически должна равняться нулю, за исключением тех случаев, когда в число измеряемых отведений попадают отведения с «необычной» проекцией конечной части петли T на ось отведения [8], разброс значений tQTa(l), и, следовательно, величина дисперсии QTaD имеют четкий физический смысл, являются объективными параметрами измерений и должны определяться формой петли, проекцией «касательных» точек петли на оси отведений и скоростью перемещения конца моментного вектора по петле T.
Для вытянутых T-петель величина QTaD ожидается меньшей, чем для широких округлых петель. Пространственные повороты сердца, приводящие к сложным поворотам интегральных векторов QRS [28] и, следовательно, к изменениям T, должны приводить к перераспределению точек Ta на петле T, т.е. к изменению распределения tQTa(l) и, возможно, к изменению величины QTaD. Эти рассуждения справедливы для дипольного приближения, а для близких к сердцу отведений свой вклад в формирование и изменение tQTa(l) могут вносить также высшие мультипольные компоненты эквивалентного генератора.
В силу априори кажущейся относительной простоты процедуры определения максимального значения амплитуды зубца T и соответствующей этому максимуму временной точки, погрешность измерения tQTa(l) и QTaD должна быть меньше, а воспроизводимость результатов выше, чем при измерении параметров tQT(l) и QTD. Меньше должен быть и разброс между ручными и автоматическими измерениями. Увеличения погрешности измерения tQTa и QTaD можно ожидать только в тех случаях, когда в одном или в нескольких из рассматриваемых отведений наблюдаются зубцы T с плоской вершиной или низкоамплитудные, на уровне шумов, зубцы T.
Похоже, однако, что параметр tQTa, как временной интервал до вершины зубца T, имеет свой прямой смысл только в том случае, когда петля T проектируется целиком на одну из полуосей отведения, положительную или отрицательную, т.е. когда зубец T монофазен. Проблема двуфазных зубцов T, т.е. выбор вершины T, может возникать в отведениях, ось которых почти нормальна интегральному вектору зубца T.
Что касается диагностической ценности параметров tQTa(l) и QTaD, то следует отметить, что пространственная форма петли T и скорость ометания ее моментным вектором определяются электрофизиологическим состоянием сердца и распределением форм кривой ТМП по поверхности электрически активного миокарда. Поэтому такие параметры зубца T, как распределение tQTa(l) и характеристика его неоднородности - QTaD, априори оцениваемые к тому же как достаточно стабильно воспроизводимые и помехоустойчивые, могут, в принципе, оказаться потенциально информативными, конечно, наряду и в совокупности с другими параметрами процесса реполяризации.
Перейдем теперь к экспериментальной оценке tQTa(l) и QTaD для практически здоровых лиц (с целью получения статистических оценок этих параметров для нормы), а также к проверке некоторых гипотез на биофизических моделях.
МАТЕРИАЛ И МЕТОДЫ
Экспериментальная оценка параметров QTa-интервала выполнена с помощью разработанного в Институте теоретической и экспериментальной биофизики РАН пакета компьютерных программ «Уран» для синхронной регистрации и автоматической обработки ЭКГ в 12 стандартных и 3 ортогональных отведениях. Пакет адаптирован, в частности, к компьютерному электрокардиографу ИДК ЭКГ 12-1.1 (фирма «Геолинк») и в такой версии эксплуатируется в Отделе массовых обследований ГосНИЦентра профилактической медицины МЗ РФ.
В качестве материала к данной работе из общей базы данных были использованы записи 12 общепринятых и 3 ортогональных по Франку отведений 250 практически здоровых лиц (220 мужчин и 30 женщин; возраст, соответственно, 37±10 и 37±7 лет). На этой же выборке ранее [8, 14] были получены для нормы оценки длительности и дисперсии QT-интервала для стандартной и ортогональной систем отведений.
Автоматическая и ручная обработка сигналов и определение параметров QTa-интервала велись по методикам, сходным с описанными в [8, 14]. Ручная разметка вершины зубца T осуществлялась при помощи курсоров на экране компьютера по крупномасштабной ЭКГ (КМ ЭКГ). Для КМ ЭКГ использовались масштабы по усилению и скорости развертки, соответственно в диапазонах 60-100 мм/мВ и 150-200 мм/с. Для обеих исследуемых систем отведений (стандартной и ортогональной по Франку) в следующих друг за другом 3-4 кардиоциклах каждой записи определялось распределение по отведениям длительностей и корректированных согласно формуле H.C.Bazett [27] длительностей QTa-интервала, соответственно, tQTa(l) и tQTac(l), а также дисперсия QTaD и корректированная дисперсия QTaDc (по всем отведениям каждой из систем).
В компьютерных экспериментах использовались биофизическая модель электрической активности сердца и система моделирования, кратко описанные в [8]. Параметрами модели являются электрофизиологические и анатомо-физиологические характеристики сердца: геометрия желудочков и специализированной проводящей системы, скорости распространения процесса деполяризации по миокарду, ножкам пучка Гиса и волокнам Пуркинье, форма ТМП на эпикарде и эндокарде желудочков, ориентация собственных анатомических осей сердца относительно исходной системы координат и некоторые другие биофизические характеристики миокарда.
Примеры исходных модельных ЭКГ для одного из возможных наборов среднефизиологических значений параметров модели, принятых за условную норму, приведены в [8] (стандартные отведения) и [14] (X, Y, Z отведения системы Франка). Множество модельных ЭКГ обрабатывалось, как и множество реальных ЭКГ, алгоритмами идентификации и измерения зубцов и интервалов кардиосигнала с целью получения и сравнительного анализа значений tQTa(l) и QTaD.
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
1. Табл. 1 и табл. 2 содержат результаты ручных (по КМ ЭКГ) измерений средних значений длительности и корректированной длительности QTa-интервалов в каждом из 12 отведений стандартной системы (табл. 1) и в 3 отведениях ортогональной системы Франка (табл. 2) для мужчин и для женщин.
Таблица 1. Распределение tQTa(l) и tQTac(l) в стандартных отведениях.
Отведения |
tQTa(l) ± s, мс |
tQTac(l) ± s, мс1/2 |
||
М |
Ж |
М |
Ж |
|
I |
306 ± 24 |
314 ± 24 |
320 ± 19 |
331 ± 17 |
II |
312 ± 22 |
321 ± 25 |
326 ± 18 |
337 ± 17 |
III |
314 ± 31 |
319 ± 31 |
329 ± 30 |
335 ± 26 |
aVR |
308 ± 23 |
315 ± 26 |
323 ± 19 |
332 ± 18 |
aVL |
300 ± 30 |
307 ± 28 |
315 ± 27 |
322 ± 21 |
aVF |
315 ± 22 |
324 ± 25 |
330 ± 20 |
341 ± 16 |
V1 |
303 ± 38 |
316 ± 40 |
316 ± 34 |
332 ± 34 |
V2 |
294 ± 21 |
305 ± 22 |
307 ± 19 |
321 ± 22 |
V3 |
297 ± 21 |
315 ± 29 |
311 ± 17 |
331 ± 21 |
V4 |
302 ± 23 |
316 ± 25 |
316 ± 18 |
333 ± 17 |
V5 |
309 ± 24 |
318 ± 26 |
323 ± 18 |
335 ± 17 |
V6 |
311 ± 23 |
320 ± 25 |
325 ± 18 |
336 ± 16 |
Таблица 2. Распределение tQTa(l) и tQTac(l) в ортогональной системе Франка.
tQTa(l) ± s, мс |
tQTac(l) ± s, мс1/2 |
|||
Отведения |
М |
Ж |
М |
Ж |
X |
298 ± 24 |
307 ± 23 |
313 ± 19 |
329 ± 17 |
Y |
302 ± 22 |
308 ± 22 |
316 ± 20 |
330 ± 17 |
Z |
286 ± 21 |
297 ± 29 |
300 ± 19 |
318 ± 26 |
Длительность QTa-интервала у женщин во всех отведениях оказалась в среднем больше чем у мужчин. Отношение длительностей tQTa/tQT в каждом отведении близко к 0,7. Отведения, в которых достигались максимальное и минимальное значения длительности tQTa у конкретного пациента не совпадали, как правило, с отведениями, в которых таковыми являлись длительности tQT.
2. Автоматические измерения длительности QTa-интервала (которые в данном исследовании хронологически предшествовали измерениям ручным) выявили значительное количество случаев большой вариабельности измеряемых значений дисперсии QTa-интервала, QTaD, от цикла к циклу. Это создавало впечатление, как выяснилось далее - неоправданное, наличия у этих пациентов повышенного значения временной дисперсии QTaD.
С помощью КМ ЭКГ были проанализированы случаи (их оказалось около 50 из общего числа 250 наборов ЭКГ), когда разница длительностей интервала QTa, измеренных автоматически в любых двух из обработанных для данного пациента кардиоциклов, составляла более 30 мс и не укладывалась при этом в рамки формулы H.C.Bazett. Большинство рассмотренных случаев, как и следовало ожидать, было представлено плоскими «шероховатыми» (рис. 2) или двуфазными с почти равными амплитудами (рис. 3) зубцами T. ВКГ-картина этих случаев соответствует таким формам петель T и их ориентации, когда некоторая слабо изогнутая часть петли или ее интегральный вектор оказываются почти ортогональными к оси какого-либо отведения.
Рис. 2. Пример зубца T с «плоской» вершиной (Пациент 381900). Петля T во фронтальной плоскости (a) и ЭКГ в отведении aVF (b). Слабо изогнутая часть петли T почти ортогональна к оси aVF, что может приводить (при незначительных случайных колебаниях изолинии) к увеличению погрешности во время измерения величины tQTa(aVF).
Рис. 3. Двуфазный зубец T с почти равными амплитудами обеих фаз в III отведении (пациент 382273). Интегральный вектор петли T во фронтальной плоскости почти ортогонален оси III отведения. Возможен триггерный эффект при определении величины tQTa(III), например, в последовательности кардиоциклов, что может привести к необъективному завышению оценки временной дисперсии QTa-интервала.
Понятно, что в этих случаях небольшая (например, дыхательная [29]) ротация сердца или слабые колебания изолинии могут приводить к скачкообразным цикл-цикл изменениям длительности QTa-интервала в отдельных отведениях и, как следствие, величины QTa-дисперсии. Эти изменения не имеют, естественно, диагностического значения и являются вторичными по отношению к первичным изменениям QTaD за счет нарушений процесса реполяризации.
На рис. 4 изображены схемы двух систем представления осей отведений для стандартной системы во фронтальной (F) и горизонтальной (H) плоскостях: шестиосевая система Бейли для отведений от конечностей (a) и шестиосевая система для грудных отведений (b). Принимая во внимание статистические характеристики ориентации максимального вектора петли T [30-31] (сектора ATS F,H на рис. 4 a,b соответствуют 95-процентным доверительным границам) и вид осредненных петель T [32] для нормальной ЭКГ, можно ожидать, что наиболее часто «плохими» для определения вершины T отведениями в норме должны оказываться III стандартное отведение во фронтальной плоскости и V 1 в горизонтальной (т.е. отведения, оси которых попадают в сектора BLS F,H , ортогональные секторам ATS F,H ), что и подтверждается результатами наших измерений. На рис. 4 показаны также F и H проекции петли T для одного из пациентов выборки (наибольшая погрешность цикл-цикл измерений tQTa наблюдалась в отведении III). При выраженных горизонтальной или вертикальной позиции сердца проблемы измерения tQTa, как это можно легко представить по рис. 4, могут возникать, соответственно, в отведениях aVF или aVL (в F), а также в отведениях V 2 или V 6 (в H). Чтобы избежать больших погрешностей при оценке средних величин QTaD, определение tQTa в данной работе велось нами далее по результатам ручных измерений. Рекомендации по объективизации измерений в автоматическом режиме представлены ниже.
Рис. 4. Схемы шестиосевой системы Бейли для отведений от конечностей (a) и шестиосевой системы грудных отведений (b). Сектор ATS соответствует 95-процентным доверительным границам ориентации проекции максимального вектора петли T в F и H плоскостях [2]. В сектора BLS, ортогональные по своему положению сектору ATS, попадают оси тех отведений, в которых в норме могут наблюдаться максимальные погрешности определения положения вершины зубца T. Показаны проекции петли T одного из пациентов выборки на F и H плоскости.
3. В табл. 3 приведены средние значения дисперсии и корректированной дисперсии QTa-интервала в стандартной и ортогональной системах отведений для мужчин и женщин на фоне соответствующих значений дисперсии QT-интервала из [8, 14].
Таблица 3. Сравнительная оценка QTD и QTaD, QTDc и QTaDc.
Отведения |
QTD ± s, мс |
QTaD ± s, мс |
QTDc ± s, мс1/2 |
QTaDc ± s, мс1/2 |
|
М |
Станд. |
56 ± 17 |
50 ± 18 |
59 ± 17 |
52 ± 18 |
Ортог. |
18 ± 13 |
20 ± 10 |
18 ± 9 |
21 ± 11 |
|
Ж |
Станд. |
54 ± 15 |
51 ± 17 |
58 ± 17 |
53 ± 18 |
Ортог. |
16 ± 11 |
23 ± 22 |
17 ± 12 |
24 ± 22 |
Средние значения дисперсии и корректированной дисперсии QTa-интервала, рассчитанные как QTaD = tQTa(l)max - tQTa(l)min по всем 12 отведениям стандартной системы и по 3 отведениям системы Франка, у мужчин и женщин практически не различаются. Сравнение QTaD для стандартной и ортогональной систем показывает, что QTaD 3 значительно меньше чем QTaD 12 . Величина QTa-дисперсии имеет тот же порядок что и величина QT-дисперсии, причем в среднем QTaD 3 несколько больше чем QTD 3 , а QTaD 12 несколько меньше чем QTD 12 .
У конкретных пациентов разность D = QTaD - QTD может быть как положительной, так и отрицательной. В табл. 4 приведено распределение значений D для 213 мужчин.
Таблица 4. Гистограмма распределения разности D = QTaD -QTD (213 мужчин).
Δ, мс |
-50 |
-40 |
-30 |
-20 |
-10 |
0 |
10 |
20 |
30 |
40 |
50 |
n |
5 |
8 |
21 |
27 |
22 |
73 |
27 |
15 |
11 |
3 |
1 |
4. В качестве примера формирования величины QTaD в стандартной и ортогональной системах отведений приведем типичный результат измерения корректированных значений tQTac(l) во фронтальной (F) и горизонтальной (H) плоскостях для одного конкретного пациента из общей выборки (пац. Ю., муж., 40 лет):
tQTac(I, II, III, aVR, aVL, aVF) = 368, 363, 380, 361, 330, 368 мс;
tQTac(V 1 - V 6 ) = 391, 332, 344, 351, 342, 347 мс;
tQTac(X, Y, Z) = 347, 364, 336 мс.
ВКГ-петли в F и H плоскостях для этого примера представлены на рис. 5. Петли построены по следующим парам сигналов: X-Y и I-aVF (в F плоскости), X-Z и V 6 -V 2 (в H плоскости). Для удобства сравнения результатов петли двух систем были промасштабированы по максимальным значениям сигналов соответствующих пар отведений.
Рис. 5. ВКГ-петли зубца T в F и H плоскостях для типичного значения величины QTaD в норме (пац. Ю. из общей выборки, муж., 40 лет). QTaD 12 = 58 мс, QTaD 3 = 27 мс; корректированные значения QTaD: QTaDc 12 = 61 мс, QTaDc 3 = 28 мс.
Посмотрим сначала дисперсию значений tQTac у этого пациента раздельно в двух плоскостях для сигналов стандартной и ортогональной систем:
QTaDc(I - aVF) = 50 мс, QTaDc(V 1 - V 6 ) = 59 мс,
QTaDc(X, Y) = 17 мс, QTaDc(Y, Z) = 28 мс.
В результате (объединяя измерения tQTac в двух плоскостях и находя новые значения tQTac max и tQTac min ) получим: QTaDc 12 = 61 мс, QTaDc 3 = 28 мс. Меньшее значение QTaD 3 по сравнению с QTaD 12 в данном случае определяется, в частности, более узкой, как видно из рис. 5, формой петли, построенной на сигналах ортогональной системы отведений, а также тем, что «касательные» точки некоторых отведений стандартной системы оказываются на петле левее и правее (по ходу петли), чем касательные точки трех отведений ортогональной системы.
Понятно, что величина QTaD определяется формой пространственной петли T, ее ориентацией в пространстве относительно осей отведений и скоростью перемещения конца моментного вектора T вдоль петли. QTaD зависит также, но в меньшей степени, от проекции начальной части петли QRS на оси отведений (фазовый сдвиг t 0 QRS согласно [8]). Изменения петли, включая патологические, будут приводить к перераспределению положения на ней точек, в которых нормали к осям отведений касаются петли T (см. рис. 1), и, следовательно, к увеличению или уменьшению QTa-дисперсии.
Основной вопрос - как связаны, если такая связь существует, возможные изменения негомогенности процесса реполяризации в желудочках и изменения величины QTaD. Рассмотрим некоторые аспекты этого вопроса в рамках полученных результатов модельных экспериментов.
5. Результаты компьютерного эксперимента (для тех же параметров модели, которые были использованы в работах [8, 14]) по исследованию распределения tQTa(l) в стандартной и ортогональной системах отведений для разных вариантов распределения длительностей ТМП на поверхности s электрически активного миокарда желудочков сведены в табл. 5.
Таблица 5. Изменения распределения tQTa(l) в модели электрической активности сердца для разных вариантов распределения двух форм ТМП (длинного - Long и короткого - Short) на эндокарде и эпикарде моделируемых желудочков. Рассчитанные значения QTaD 12 и QTaD 3 (результаты измерения указаны в мс ).
|
Варианты распределения (ТМП endo / ТМП epi) |
|||
Long / Short |
Short / Long |
Long / Long |
Short / Short |
|
I |
288 |
290 |
343 |
259 |
II |
300 |
295 |
348 |
288 |
III |
355 |
348 |
355 |
295 |
AVR |
300 |
295 |
343 |
286 |
AVL |
288 |
281 |
* |
* |
AVF |
312 |
305 |
348 |
288 |
V1 |
319 |
346 |
319 |
262 |
V2 |
302 |
307 |
343 |
281 |
V3 |
295 |
288 |
353 |
293 |
V4 |
293 |
288 |
312 |
252 |
V5 |
305 |
290 |
324 |
264 |
V6 |
312 |
295 |
324 |
281 |
QTaD12 |
67 |
67 |
43 |
43 |
X |
288 |
283 |
** |
252 |
Y |
312 |
314 |
** |
288 |
Z |
295 |
305 |
** |
288 |
QTaD3 |
24 |
31 |
** |
36 |
Были заданы два типа ТМП, длинный (Long) и короткий (Short), с длительностями соответственно около 300 мс и 250 мс. Распределение ТМП(s) для условной нормы выбрано таким образом, что на эндокарде генерируется более длинный ТМП, чем на эпикарде. Рассмотрено 4 варианта распределения (ТМП endo / ТМП epi ): L/S, S/L, L/L, S/S.
Результаты, полученные с помощью модели электрической активности сердца для исходных значений ее параметров, соответствующих критериям условной нормы, хорошо согласуются со средними величинами характеристик QTa-интервала, полученными для нормальной ЭКГ в экспериментах с группой практически здоровых лиц (табл. 1-3). Это касается величин длительностей QTa-интервала и их наиболее часто встречающегося распределения по отведениям, а также величины QTaD. На рис. 6 показаны модельные петли зубца T, построенные в двух плоскостях по сигналам X, Y, Z отведений.
Рис. 6. Модельные петли зубца T («Пациент» 00023) в F и H плоскостях (a), исходные ЭКГ в отведениях X, Y, Z (b).
В компьютерных экспериментах по изменению распределения длинного и короткого ТМП по поверхности эпикарда и эндокарда наблюдались характерные для аналогичных электрофизиологических экспериментов [33] изменения амплитуды и полярности зубца T.
Рассмотренные варианты распределения (ТМП endo / ТМП epi ) задают фактически трансмуральную разность (дисперсию) длительностей ТМП, D tТМП. Условно нормальное значение DtТМП (вариант L/S, длительности ТМП близки, соответственно, к значениям 300 мс и 250 мс) оказывается равным 50 мс. Анализ табл. 5 позволяет отметить при этом некоторые особенности полученных результатов:
·
При изменении знака DtТМП (варианты L/S и S/L) изменения величин tQTa составили 1-5 % от длительности обобщенного QT-интервала, т.е. находились в пределах ошибки измерения. Величина QTaD 12 практически не изменилась, а изменения QTaD 3 составили около 25 %.·
Значения tQTa(l) для двух вариантов гомогенного распределения ТМП на эндокарде и эпикарде (варианты L/L и S/S) различаются во всех отведениях на величину tТМПL - tТМПS = 50 мс, что не удивительно. Более интересным наблюдением является тот факт, что в обоих случаях при имеющем место равенстве нулю величины D tТМП (т.е. тогда, когда на эпикарде и эндокарде генерировался один и тот же ТМП) величина измеренного значения QTaD оставалась равной некоторому постоянному значению. Эта «остаточная» величина QTaD зависит, в частности, как показывают эксперименты, от геометрии модели.Вращение модели желудочков вокруг анатомических осей сердца приводит к изменению распределения tQTa(l) и измеренных значений QTaD, а также к перераспределению отведений, в которых достигаются максимальное и минимальное значения tQTa, что находится в соответствии с наблюдаемыми картинами перераспределения касательных точек на модельной петле T.
В то же время, изменения tQTa(l) и QTaD, наблюдаемые в рассмотренных вариантах распределения ТМП(s), не поддаются аналитическому описанию. Однако перспективность использования биофизических моделей генеза электрокардиосигналов для выяснения механизмов, связывающих распределение ТМП(s) с регистрируемым распределением tQTa(l) и величиной QTaD, очевидна. В возникающем на практике распределении tQTa(l) существенную роль, наряду с ТМП(s), играют, несомненно, и позиционные факторы, включая геометрию миокарда, которая должна оказывать влияние на распределение по величине элементов матрицы k li [8], связывающей в прямой задаче ТМП(s) с регистрируемыми потенциалами. Необходим многофакторный компьютерный эксперимент.
6. Возвращаясь к практическим измерениям QTa-интервала, рассмотрим, опираясь на полученные результаты и не вдаваясь в технические подробности, вопрос о повышении объективности и точности автоматического измерения tQTa(l) и QTaD.
Как было показано, наибольшую погрешность в определение вершины зубца T вносят отведения, в которых T двуфазен с почти равными амплитудами или имеет плоскую протяженную вершину. Погрешность может достигать нескольких десятков мс. Но таких отведений обычно встречается не более одного в каждой из двух плоскостей, фронтальной и горизонтальной, в которых лежат отведения стандартной и ортогональной систем. К тому же, чаще всего, это отведения с низкоамплитудным плохо идентифицируемым зубцом T.
Поэтому такие отведения могут быть без потери информации исключены из рассмотрения при расчете величины QTaD. Определять «плохие» отведения достаточно удобно в полуавтоматическом режиме с помощью компьютерной программы, которая выводит на экран схему шестиосевой системы Бейли во фронтальной плоскости и схему осей шести грудных отведений в горизонтальной плоскости (рис. 4), а также соответствующие проекции петли T. В автоматическом режиме ось отведения-кандидата на исключение из рассмотрения определяется по признаку ее ортогональности интегральному вектору петли T, а плоские вершины T аппроксимируются параболой, что повышает воспроизводимость результатов измерений от цикла к циклу [34].
ВЫВОДЫ
1. В отличие от дисперсии QT-интервала (QTD), которая теоретически должна равняться нулю, а в реальных измерениях с порогом зависит от многих субъективных факторов, разброс по отведениям значений длительности интервала QTa, tQTa(l), и дисперсия QTaD имеют четкий физический смысл, являются объективными параметрами измерений и в рамках векторкардиографической модели трактуются в терминах положения на петле T точек, в которых касаются петли нормали к осям отведений.
2. Величина QTaD зависит, главным образом, от формы пространственной петли T, ее ориентации в пространстве относительно осей отведений и от скорости перемещения конца моментного вектора T вдоль петли. Изменения формы и ориентации петли T, включая патологические, приводят к перераспределению положения на ней «касательных» точек и, следовательно, к увеличению или уменьшению QTa-дисперсии.
3. Диагностическая ценность параметров tQTa(l) и QTaD основывается на моделях современной электрокардиологии, согласно которым пространственная форма петли T и скорость ометания ее моментным вектором определяются электрофизиологическим состоянием сердца и распределением форм ТМП на поверхности электрически активного миокарда. Поэтому такие параметры реполяризации желудочков, как распределение tQTa(l) и характеристика его неоднородности - QTaD, являются, в принципе, потенциально информативными по отношению к нарушениям процесса реполяризации.
Использование биофизических моделей генеза электрокардиосигналов в режиме прямой задачи позволяет проследить в общих чертах связь между распределением ТМП в сердце и параметрами реполяризационной части кардиоцикла.
4. Необходимость совершенствования компьютерных методов практического измерения реполяризационных параметров кардиоцикла определяется, в частности, важностью проблемы поиска неинвазивных предикторов возникновения нарушений сердечного ритма.
Возможность повышения объективности и точности автоматического измерения tQTa(l) и QTaD от цикла к циклу позволяет оценить эти параметры как более стабильно воспроизводимые и помехоустойчивые по сравнению с tQT(l) и QTD.
Учитывая также преимущества QTaD перед QTD по части физико-физиологической содержательности, можно рекомендовать вычеркнуть QTD из числа номинантов в категории «Самый информативный показатель негомогенности реполяризации», не подвергая в то же время сомнению диагностическую ценность обобщенной длительности QT-интервала.
5. Что касается прогностической значимости величины QTaD, как в пространственной, так и во временной областях, для определения риска развития тяжелых нарушений ритма, то эта гипотеза продолжает оставаться в фокусе экспериментальных исследований, находя свои частные подтверждения или опровержения. Вопрос теоретической оценки существования и механизмов такой значимости остается, на наш взгляд, открытым и требует дальнейшей углубленной разработки биофизической модели, связывающей распределение ТМП в сердце с регистрируемыми параметрами процесса реполяризации на торсе и, главное, разработки, если это в принципе возможно, корректных методов решения обратной задачи, т.е. оценки локальной негомогенности длительности ТМП по неинвазивным измерениям экстраклеточного электрического поля сердца.
Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и правительства Московской области (проект № 01-01-97003).
ЛИТЕРАТУРА
1. Day C.P., McComb J.M., Campbell R.W.F. - QT dispersion: an indication of arrhythmia risk in patients with long QT intervals - Br. Heart J., 1990, 63, 342-344.
2. Higman P.D., Campbell R.W.F. QT dispersion. - Br. Heart J., 1994, 71, 508-510.
3. Lux R.L., Fuller M.S., MacLeod R.S. et al. QT Interval Dispersion: Dispersion of Ventricular Repolarization or Dispersion of QT Interval? - J. Electrocardiol, 1997, 30 (suppl.), 176-180.
4. Kautzner J., Malik M. QT Interval Dispersion and its Clinical Utility - Pacing Clin. Electrophysiol, 1997, 20, 2625-2640.
5. Malik M., Camm A.J. Mystery of QTc Interval Dispersion - Am. J. Cardiol., 1997, 79 (6), 785-787.
6. Никитин Ю.П., Кузнецов А.А. Дисперсия интервала QT (обзор). - Кардиология, 1998, 5, 58-63.
7. Antzelevitch C., Shimizu W., Yan G.X., Sicouri S. Cellular Basis for QT Dispersion - J. Electrocardiol, 1998, 31 (suppl.), 168-175.
8. Баум О.В., Попов Л.А., Волошин В.И., Муромцева Г.А. QT-дисперсия: Модели и измерения. - Вестник Аритмологии, 2000, 20, 6-17.
9. Lee K.W., Kligfield P., Dower G.E., Okin P.M. QT dispersion, T-wave projection, and heterogeneity of repolarization in patients with coronary artery disease. - Am. J. Cardiol., 2001, 87, 2, 148-151.
10. McLaughlin N.B., Campbell R.W.F., Murray A. Comparison of Automatic QT Measurement Techniques in the Normal 12 Lead Electrocardiogram - Br. Heart J., 1995, 74, 84-89.
11. Hill J.A., Friedman P.L. Measurement of QT interval and QT dispersion. - Lancet, 1997, 349, 894-895.
12. Batchvarov V., Yi G., Guo X., Savelieva I., Camm A.J., Malik M. QT Interval and QT Dispersion Measured with the Threshold Method Depend on Threshold Level - Pacing Clin. Electrophysiol., 1998, 21, 2372-2375.
13. Kors J.A., van Herpen G. Measurement Error as a Source of QT Dispersion: a Computerized Analysis - Heart, 1998, 80 (5), 453-458.
14. Баум О.В., Попов Л.А., Волошин В.И., Муромцева Г.А. QT-дисперсия в ортогональных и псевдоортогональных отведениях. - Вестник Аритмологии, 2002, 26, 49-56.
15. Baum O.V., Popov L.A., Voloshin V.I., Muromtseva G.A. QT-dispersion: Models and measurements. - Electrocardiology ‘2000 (L. De Ambroggi Ed.), Roma: Casa Editrice Scientifica Internazionale, 2001.
16. Baum O.V., Popov L.A., Voloshin V.I. QT-mapping and QT-dispersion on the standard spherical surface: A model study. - XXVIII International Congress on Electrocardiology. Abstract book. San-Paolo, 2001, p. 155.
17. Surawicz B., Knoebel S. Long QT. Good, Bad or Independent? Review. - JACC, 1984, 4, 2, 398-413.
18. Остроумова О.Д. Удлинение интервала QT. - Русский медицинский журнал, 2001, 9, 18, 750-753.
19. Zareba W., Moss A.J., le Cessie S. Dispersion of ventricular repolarization and sudden death in ischaemic heart disease. - Am. J. Cardiol., 1994, 74, 550-553.
20. Zabel M., Portnoy S., Franz M.R. Electrocardiographic indexes of dispersion of ventricular repolarization: an isolated heart validation study. - J. Am. Coll. Cardiol., 1995, 25, 746-752.
21. Dilaveris P., Gialafos E., Pantazis A. et al. The spatial QRS-T angle as a marker of ventricular repolarization in hypertension. - J. Hum. Hypertens., 2001, 15, 63-70.
22. Savelieva I., Yap Y.G., Yi G. et al. Relation of ventricular repolarization to cardiac cycle length in normal subjects, hypertrophic cardiomyopathy, and patients with myocardial infarction. - Clin. cardiol., 1999, 22, 649-654.
23. Haapalahti P., Makijarvi M., Korhonen P. et al. Magnetocardiographic QT dispersion during cardiovascular autonomic function tests. - Basic Res. Cardiol., 2000, 95, 424-430.
24. Kuroda N., Ohnishi Y., Adachi K., Yokoyama M. Relationship between the QT indices and the microvolt-level T wave alternans in cardiomyopathy. - Jpn. Circ. J., 2001, 65, 974-978.
25. Benatar A., Cools F., Decraene T. et al. The T wave as a marker of dispersion of ventricular repolarization in premature infants before and while on treatment with the I(Kr) channel blocker cisapride. - Cardiol. Young, 2002, 12, 32-36.
26. Buchanan F. Trans. Oxford Univ. Junior Scientific Club, 1909; Science Progress, 1910, II, V, XVII, 60 (Öčň. ďî [27]).
27. Bazett H.C. An analysis of the time-relation of electrocardiogramms. - Heart, 1920, 7, 353-367.
28. Баум О.В., Попов Л.А., Волошин В.И., Кочеткова О.Д. Исследование на компьютерных моделях влияния анатомической ориентации сердца на результаты измерения его электрической оси. - Измерительная техника, 2001. №6, 71-75.
29. Рутткай-Недецкий И. Влияние дыхания и расположения сердца в грудной клетке на электрическое поле сердца. - Теоретические основы электрокардиологии (под ред. К.В. Нельсона и Д.Б. Гезеловица), М.: Наука, 1979, 152-167.
30. Chou T.C., Helm R.A., Kaplan S. Clnical vectorcardiography. N.Y. - L.: Grune and A
31. Титомир Л.И., Рутткай-Недецкий И. Анализ ортогональной электрокардиограммы. М.: Наука, 1990.
32. Groeben J. Observations on the variability of the electrocardiographic waveform in normal men and women. Palo Alto: Centr. Calif. Heart Assoc., 1977.
33. Baum O.V. Thermal Regimes of Heart, Myocardium Repolarization and ST-T Genesis. A Model Study. - Electrocardiology ‘97. Singapore etc.: World Scientific, 1998, 75-80.
34. Баум О.В., Попов ЛА., Волошин В.И. Алгоритмы автоматического измерения QTaD в стандартной и ортогональной системах отведений. - Измерительная Техника (в печати).