마취 중 심박동 리듬 변동 현상의 모델링과 정량화

초록

본 논문은 마취 깊이에 따라 즉시 심박동(Instantaneous Heart Rate, IHR)이 보이는 “리듬‑비리듬 전이” 현상을 수학적으로 모델링하고, 다중탭퍼 동기압축 변환(Multitaper Synchrosqueezing Transform)을 이용해 시간‑변화 전력 스펙트럼을 추정한다. 이를 기반으로 비리듬‑리듬 비율(NRR) 지수를 정의하고, 임상 데이터베이스에서 기존 마취 깊이 지표와 비교·검증하였다. 결과는 NRR 지수가 근육 반응 및 통증 자극을 기존 지표보다 조기에 포착함을 보여주며, 실시간 구현을 위한 블렌딩 보간법도 제안한다.

상세 분석

이 연구는 마취 중 발생하는 IHR의 규칙적·불규칙적 변동을 “리듬‑비리듬 전이 현상”이라 명명하고, 이를 정량화하기 위해 두 가지 핵심 기술을 도입한다. 첫째, 생리학적 특성을 반영한 적응 조화 모델(adaptive harmonic model)을 제시한다. 이 모델은 IHR을 기본 주파수 성분과 그 고조파들의 가변 진폭·위상으로 표현함으로써, 깊은 마취에서는 주파수가 고정된 강한 조화 성분이 지배하고, 얕은 마취에서는 비조화적(노이즈‑유사) 성분이 증가한다는 가정을 수학적으로 정립한다. 모델의 존재와 유일성을 보장하기 위해 신호가 충분히 스무스하고, 고조파 간 최소 주파수 간격이 존재한다는 조건을 명시한다.

둘째, 시간‑주파수 분석 도구로 다중탭퍼 동기압축 변환(Multitaper Synchrosqueezing Transform, MT‑SST)을 활용한다. 기존의 단일 윈도우 STFT나 CWT에 비해 MT‑SST는 스펙트럼 누설을 최소화하고, 시간 해상도와 주파수 해상도를 동시에 향상시킨다. 특히, 동기압축 단계에서 에너지 집중을 재배치함으로써 각 조화 성분의 순간 전력을 정확히 추출한다. 이를 통해 얻은 시간‑변화 전력 스펙트럼(tvPS)에서 리듬 성분 전력( Rhythmic Power )과 비리듬 성분 전력( Non‑rhythmic Power )을 구분하고, 두 전력의 비율을 로그 스케일로 변환한 것이 NRR 지수이다.

임상 검증에서는 100명 이상(정확한 수치는 논문에 명시되지 않음)의 마취 환자 데이터를 이용해 NRR 지수를 BIS, HRV, PSI 등 기존 지표와 비교하였다. 통계적으로 NRR는 근육 반응(예: 움직임, 피부 절개에 대한 반응)과의 상관관계가 가장 높았으며, 특히 통증 자극 시 NRR가 급격히 상승하는 패턴을 보였다. 이는 비리듬 성분이 교감신경 활성화와 직접 연관될 가능성을 시사한다. 또한, NRR와 흡입 마취제(Sevoflurane) 농도 사이의 상관관계 분석 결과, NRR가 농도 변화에 민감하게 반응함을 확인하였다.

실시간 적용을 위해 저자는 B‑스플라인 기반 블렌딩 보간 연산자를 설계하였다. 이 연산자는 비균등하게 샘플링된 R‑R 간격 데이터를 연속적인 시간축으로 매끄럽게 보간하고, 실시간 MT‑SST 계산에 필요한 고해상도 윈도우를 제공한다. 구현 결과, 1 kHz 샘플링 기준에서 200 ms 지연 이하의 실시간 NRR 추출이 가능함을 보고하였다.

한계점으로는 (1) IHR 추출 과정에서 R‑peak 검출 오류가 NRR에 직접적인 영향을 미칠 수 있음, (2) 현재 연구는 주로 흡입 마취제와 절개 자극에 초점을 맞추었으며, 다른 마취제나 환자군(소아, 고령자)에서의 일반화 검증이 부족함, (3) 실시간 구현 시 연산 부하가 높은 편이어서 임베디드 시스템 적용을 위한 최적화가 필요함을 들 수 있다.

전반적으로, 적응 조화 모델과 MT‑SST를 결합한 NRR 지수는 마취 깊이와 통증 자극을 기존 지표보다 민감하게 포착할 수 있는 새로운 바이오마커로서, 향후 마취 모니터링 시스템에 통합될 가능성을 보여준다.