BIS-Based Model Reference Adaptive Control of Anesthesia Depth under Patient Variability and Time-Varying Delays

기술 요약: 환자 변동성 및 시변 지연 조건하에서의 BIS 기반 모델 참조 적응 제어(MRAC)를 이용한 마취 깊이 조절

문제 정의
전신 마취를 위한 약물 투여는 무의식과 진통 상태를 유지하기 위해 최면제(프로포폴)와 진통제(레미펜타닐)를 정밀하게 조절하는 것을 필요로 한다. 자동 폐쇄 루프 제어(closed-loop control)의 주요 과제는 환자별 약동학/약력학(PK/PD) 파라미터의 상당한 개별 차이와 약물 분포 과정에서 발생하는 생리적 시변 지연(time-varying delays)이다. PID 또는 모델 예측 제어(MPC)와 같은 기존의 제어 전략은 종종 오프라인 방식의 환자 맞춤형 파라미터 식별에 의존하거나 약물 분포가 즉각적이라고 가정한다. 이러한 한계는 환자의 생리적 상태가 변화하거나(예: 심박출량 변화) 지연 시간이 변동하는 실제 임상 환경에서의 강건성을 저하시킨다. 또한, 약물 농도와 BIS 출력 사이의 비선형 관계(약물 상호작용 항을 포함한 Hill 함수에 의해 결정됨)는 선형 피드백 제어기 설계의 복잡성을 가중시킨다.

방법론
저자들은 프로포폴과 레미펜타닐의 이중 약물 주입을 사용하여 마취 깊이를 조절하도록 설계된 모델 참조 적응 제어(MRAC) 프레임워크를 제안한다. 본 방법론은 다음의 구성 요소들을 통합한다:

1. 플랜트 모델링(Plant Modeling): 환자는 두 약물에 대해 표준 3-구획 PK 모델(Schnider)을 사용하며, 구획 간 전이 과정에 시변 지연($\tau(t)$)을 포함하도록 확장하여 모델링된다. 약력학(PD) 성분은 약물 간의 시너지 상호작용을 고려한 비선형 Hill 함수를 통해 효과 부위 농도(effect-site concentrations)와 BIS 출력을 연결한다.
2. 역 Hill 재구성(Inverse Hill Reconstruction): Hill 함수의 비선형성과 다중 입력 제어의 모호성을 해결하기 위해, 저자들은 단순화된 역 재구성 전략을 도입한다. 정규화된 레미펜타일 효과가 정규화된 프로포폴 효과에 비례한다($Y_r \approx \mu Y_p$)고 가정함으로써, 복잡한 다중 약물 상호작용을 단일 변수 비선형 관계로 축소한다. 그 후, 측정된 BIS 신호에 역 Hill 함수($H^{-1}_\mu$)를 적용하여 선형화된 프로포폴 효과 부위 농도($\hat{Y}_p$)를 추정하며, 이를 피드백 변수로 사용한다.
3. 적응 제어 법칙(Adaptive Control Law): 본 연구의 핵심인 MRAC 스킴은 상대 차수(relative degree)가 2이고 미지의 시변 지연이 존재하는 시스템을 위해 설계되었다.
   • 참조 모델(Reference Model): 2차 참조 모델($W_m(s)$)은 원하는 과도 응답 및 정상 상태 응답(목표 BIS 50)을 정의한다.
   • 제어 구조(Control Structure): 제어 입력은 입력과 추정된 출력으로부터 유도된 필터링된 회귀 변수(filtered regressors)를 사용하여 생성된다. 상대 차수를 처리하고 안정성을 보장하기 위해, 오차 동역학을 엄격한 양의 실수(strictly positive real, SPR)로 만들기 위한 1차 안정 필터 $Z(s) = s + p_0$가 도입된다.
   • 적응 및 강건성(Adaptation and Robustness): 제어 파라미터는 Lyapunov-Krasovskii 범함수를 통해 온라인으로 업데이트된다. 미지의 시변 지연 및 외부 교란(예: 수술적 자극)을 보상하기 위해, 추적 오차의 부호와 적분 항을 포함하는 특정 강건성 항이 포함된다.
4. 안정성 분석(Stability Analysis): 이론적 안정성은 Lyapunov-Krasovskii 범함수를 사용하여 증명되며, 지연율의 변화가 1보다 작을 경우 모든 폐루프 신호가 유계(bounded)를 유지하고 추적 오차가 점근적으로 0으로 수렴함을 입증한다.

주요 기여
본 논문은 다섯 가지 구체적인 기여를 명시한다:

• 파라미터 독립성: 제안된 MRAC 방식은 기존의 PID나 이벤트 기반 제어기와 달리, 정확한 환자별 PK/PD 파라미터 식별이나 사전 튜닝을 요구하지 않는다.
• 명시적 지연 처리: 기존의 적응형 프레임워크(예: 분수 차수 MRAC)와 달리, 본 방법은 단순히 분수 차수 메모리 효과에 의존하는 것이 아니라, 생리적 시변 지연을 제어 설계 및 안정성 분석에 명시적으로 포함한다.
• 비선형 보상: 제안된 단순화된 역 Hill 재구성은 복잡한 다중 약물 상호작용에도 불구하고 비선형 BIS 신호로부터 선형 출력을 추정할 수 있게 하여, 적응형 출력 피드백 제어를 용이하게 한다.
• 강건성: 적응 법칙은 파라미터를 지속적으로 온라인 업데이트하여, 사전 제어기 재설정 없이도 환자 간 변동성과 수술 중 교란에 대해 안정성을 확보한다.
• 임상적 안전성: 제어기는 주입 속도를 임상적으로 허용되는 표준 범위 내로 조절하여, 과다 투여와 불충분한 마취를 모두 방지한다.

결과
제안된 방법은 다양한 생체 특징(연령, 체중, 키, 성별)과 무작위로 생성된 PK/PD 파라미터를 가진 12명의 가상 환자군을 대상으로 한 시뮬레이션을 통해 검증되었다.

• 유도 단계(Induction Phase): 제어기는 모든 환자에 대해 BIS를 기저치에서 목표치인 50까지 성공적으로 유도하였다. 목표 도달 시간(TT)은 104193초(평균 ~145초)였으며, 정착 시간(ST10)은 141360초(평균 ~263초)였다. BIS 최저치(BIS-NADIR)는 모든 경우에서 40 이상을 유지하여 과도한 마취 깊이를 피했다.
• 유지 단계(Maintenance Phase): 시뮬레이션된 수술적 교란(이중 단계 입력) 하에서도 제어기는 BIS를 목표 범위로 빠르게 복구하였다. 양(+)의 교란 및 음(-)의 교란으로부터의 회복 시간은 짧았으며(평균 ~0.6분), BIS는 40–60 범위 내에 머물렀고 오버슈트는 5% 미만, 편차는 ±2 이내로 매우 작았다.
• 약물 투여: 프로포폴과 레미펜타일의 주입 프로파일은 유도 및 유지 단계 모두에서 전형적인 임상 제한 범위 내에 머물렀다.
• 비교: 제안된 MRAC는 환자별 모델 파라미터나 오프라인 튜닝 없이도 경쟁력 있는 추적 성능을 달면서, 기존 전략들(분수 차수 MRAC, 최적화된 PID, 비선형 MPC)과 대등하거나 더 우수한 성능을 보였다.

의의 및 주장
저자들은 본 연구가 강건하고 개별화된 폐쇄 루프 마취 관리를 위한 유망한 솔루션을 제공한다고 주장한다. 본 연구의 주요 의의는 제어기가 환자별 파라미터 식별 없이도 효과적으로 작동할 수 있다는 점에 있으며, 이는 자동 마취 시스템의 임상 적용에 있어 주요 장벽을 해결한다. 시변 지연과 비선형 약물 상호작용을 명시적으로 모델링하고 보상함으로써, 제안된 MRAC 프레임워크는 다양한 생리적 조건에서도 안정성과 안전성을 유지할 수 있는 이론적 근거가 확실하고 확장 가능한 접근 방식을 제공한다. 본 연구는 이러한 시스템이 수동 방식이나 덜 적응적인 자동화 방식에 비해 환자 안전을 높이고, 약물 소비를 줄이며, 회복 시간을 개선할 수 있음을 시사한다.