제세동

계산물리 물리 전자기해석 마지막 수정: 2026-08-19 04:49:55

1. 개요[편집]

제세동 Defibrillation
목적	모든 재진입 파원을 동시에 소멸시키기
필수 모형	이중영역(bidomain) — 단일영역으로는 재현 불가
핵심 현상	가상전극 분극(VEP) · 개폐 여기(break excitation)
경쟁 가설	임계질량 · 취약성 상한(ULV)
에너지	체외 이상성 120~200 J · 삽입형 5~20 J
계산 병목	타원형 방정식 반복 풀이 (전체 시간의 대부분)

제세동(defibrillation)은 세동 중인 심근 전체에 짧고 강한 전기 충격을 가해 존재하는 모든 재진입 파원을 한꺼번에 소멸시키고, 정상 박동원이 다시 주도권을 잡을 수 있는 상태로 되돌리는 조작이다. 모형 쪽 언어로 옮기면 시공간 카오스 상태의 반응-확산계를 단 한 번의 공간 균일 섭동으로 균일 정지 상태의 흡인영역 안으로 밀어 넣는 문제다.

부정맥이 어떻게 생기고 유지되는지는 부정맥과 나선파에 있고, 이 문서는 그 반대 방향 — 끄는 쪽 — 만 다룬다. 특히 이 주제가 계산 심장학에서 특별한 이유는 하나다. 부정맥 시뮬레이션에서 잘 쓰던 단일영역 모형이 여기서는 원리적으로 틀린 답을 낸다.

2. 왜 단일영역으로는 안 되는가[편집]

부정맥에서 쓴 단일영역(monodomain) 모형은 세포 안팎을 하나의 유효 확산 매질로 뭉갠 것이다. 이중영역(bidomain) 모형은 세포내 공간과 세포외 공간을 각각 독립한 연속체로 두고, 두 전위 $\phi_i$ , $\phi_e$ 를 별도 미지수로 들고 간다. $V_m = \phi_i - \phi_e$ 로 놓으면

\chi\left(C_m\frac{\partial V_m}{\partial t} + I_{\rm ion}\right) = \nabla\cdot(\sigma_i \nabla V_m) + \nabla\cdot(\sigma_i \nabla \phi_e)

\nabla\cdot\big((\sigma_i + \sigma_e)\nabla\phi_e\big) = -\nabla\cdot(\sigma_i \nabla V_m) - I_{\rm e}

가 된다. 위는 시간 전진하는 포물형 방정식, 아래는 매 시각 풀어야 하는 타원형 방정식이다.

여기서 결정적인 조건이 나온다. 만약 두 전도도 텐서가 비례한다면( $\sigma_e = \lambda\,\sigma_i$ , 이른바 등이방성비), 아래 식을 대입해 $\phi_e$ 를 소거할 수 있고 계는 정확히 단일영역으로 축약된다. 즉 단일영역은 근사가 아니라 등이방성비라는 가정 아래의 정확한 축약이다.

그런데 실측은 그렇지 않다. 세포내 공간은 섬유 방향 대 횡방향 전도도비가 10 안팎인 데 비해, 세포외 공간은 2~3 정도로 훨씬 등방에 가깝다. 두 이방성 비가 다르다는 이 사실 하나가 모든 것을 바꾼다.

자기 자극에 의한 부정맥·전파 문제에서는 차이가 대체로 작다. 파면이 스스로 만드는 $\phi_e$ 는 작고, 단일영역이 몇 % 오차 안에서 같은 답을 낸다. 그래서 비용이 한 자릿수 이상 싼 단일영역을 쓰는 것이 국룰이다.
외부에서 전류를 주입할 때는 완전히 다르다. 등이방성비라면 균일 조직 내부에 아무리 센 외부 전기장을 걸어도 막전위는 경계에서만 변한다. 내부 어디에서도 $V_m$ 이 바뀌지 않으므로 충격이 조직 심부에 도달하지 못한다는 결론이 나오고, 이것은 실험과 정면으로 어긋난다. 이방성 비가 다르면 조직 내부 곳곳에서 막을 가로지르는 전류원이 생겨나며, 이것이 제세동이 작동하는 물리적 기반이다.

3. 가상전극 분극[편집]

1989년 세풀베다·로스·윅스워가 이중영역 방정식을 유한차분으로 풀어 얻은 예측이 이 분야의 전환점이었다. 조직 표면 한 점에 음극 자극을 주면, 직관적으로는 그 주변이 둥글게 탈분극될 것 같지만 실제 답은 개뼈다귀(dog-bone) 모양이었다. 섬유 방향으로 길게 탈분극된 영역이 뻗고, 그 양옆 횡방향으로는 과분극된 영역이 나란히 생긴다. 전극이 하나뿐인데 조직 안에 양극과 음극이 동시에 만들어진 것처럼 보이므로 이것을 가상전극 분극(virtual electrode polarization, VEP)이라 부른다.

몇 년 뒤 전압 감응 염료를 쓴 광학 매핑이 이 개뼈다귀 패턴을 그대로 찍어 내면서, “수치 모형이 먼저 예측하고 실험이 나중에 확인한” 계산 생물학의 대표 사례가 됐다.¹

VEP가 왜 중요한가는 충격이 끝나는 순간에 드러난다.

충격 도중(make) 탈분극된 영역은 흥분하지만, 과분극된 영역은 오히려 흥분성이 강제로 회복된다 — 불활성화가 풀리고 막이 “재장전”된다.
충격이 끝나면(break) 탈분극 영역에 남은 전하가 옆의 재장전된 영역으로 흘러 들어가 새 파면을 만든다. 자극이 없어진 뒤에 흥분이 시작되는 이 현상이 개폐 여기(break excitation)이며, 고전적인 자극 이론에는 없던 것이다.

새로 태어난 파면은 탈분극 영역과 과분극 영역의 경계를 따라 생기므로 자유단을 가진 조각이다. 자유단은 감기고, 감기면 나선이 된다. 이것이 가상전극 유발 위상 특이점이며, 제세동 충격이 세동을 끄기도 하지만 만들기도 하는 이유를 하나의 기전으로 설명한다. 나선파에서 “어중간한 세기의 충격이 오히려 새 나선을 만든다”고 적힌 것의 정체가 이것이다.

4. 임계질량 대 취약성 상한[편집]

충격이 성공하는 조건에 대해 오래 경쟁한 두 가설이 있다.

임계질량 가설. 심근의 충분히 큰 비율(흔히 75~90%로 인용된다)을 동시에 탈분극시켜 불응으로 만들면, 살아남은 파면이 갈 곳을 잃어 소멸한다는 그림이다. 직관적이고, 세동 유지에 최소 조직 크기가 필요하다는 오래된 관찰과도 결이 맞는다. 문제는 이 가설이 성공한 충격 뒤에도 일부 심근이 여전히 흥분 가능하다는 관측을 설명하지 못한다는 것이다.

취약성 상한 가설(upper limit of vulnerability, ULV). 관점을 뒤집는다. 정상 리듬 중에 재분극 시기(심전도의 T파 구간)에 충격을 주면 오히려 세동이 유발되는데, 이 유발이 일어나는 세기에는 상한이 있다. 그보다 세게 때리면 세동이 유발되지 않는다. 이 상한이 ULV다. 그리고 세동 중에 준 충격이 실패하는 이유는 기존 파면을 못 지워서가 아니라 충격 자신이 만든 가상전극 분극으로 세동을 다시 켰기 때문이라는 것이 이 가설의 주장이다. 따라서 제세동 역치는 ULV와 거의 같아야 하고, 실제로 두 값이 잘 맞아떨어진다는 관측이 임상에서 반복 확인됐다.

VEP 기전은 두 가설 중 후자를 자연스럽게 지지한다. “지우기에 충분한 세기”가 아니라 “새로 만들지 않을 만큼 큰 세기”가 기준이라는 뒤집힌 논리는, 임상적으로도 실용적 결과를 낳았다. 삽입형 제세동기의 여유 마진을 확인할 때 매번 세동을 유발했다가 끄는 위험한 절차 대신, 정상 리듬 중 T파에 충격을 주어 ULV를 재는 방법이 대안으로 쓰인다.

5. 역치와 에너지[편집]

제세동 역치(defibrillation threshold, DFT)라는 이름은 다소 부정확하다. 실제 관계는 계단 함수가 아니라 S자 용량-반응 곡선이다. 같은 세기의 충격이 어떤 때는 성공하고 어떤 때는 실패하므로, 정직한 기술량은 성공 확률 50%(또는 95%)에 대응하는 세기다. 임상에서 흔히 쓰는 단계적 감소 프로토콜로 얻은 “DFT”는 그 확률 분포에서 뽑은 표본 하나에 가깝고, 반복 측정하면 상당히 흔들린다.²

에너지 규모는 두 갈래다.

체외 제세동기. 과거 단상성 파형에서 200~~360 J를 쓰던 것이, 이상성(biphasic) 절단 지수 파형이 도입되며 120~~200 J로 내려왔고 첫 충격 성공률도 올라갔다.
삽입형 제세동기(ICD). 전극이 심장 안에 있어 훨씬 적게 든다. 역치가 10 J 안팎이고 기기 최대 출력이 30~40 J 규모라, 그 차이가 안전 마진이 된다.

파형이 왜 중요한가도 가상전극으로 읽힌다. 충격이 끝나는 시점에 막에 남은 잔류 분극이 개폐 여기의 씨앗이므로, 두 번째 위상에서 극성을 뒤집어 그 잔류 전하를 되돌려 놓으면 충격 직후 막이 정지 상태에 가까워지고 세동 재유발 확률이 떨어진다. 이상성 파형의 이점을 설명하는 이 “전하 되돌리기” 논리는, 방전 파형이 실제로는 축전기의 지수 방전을 잘라 쓴 것이라는 회로적 제약과 맞물려 파형 설계 문제(절단 시점, 두 위상의 길이 비)를 낳는다.

6. 저에너지 전략[편집]

강한 충격은 아프고, 심근을 손상시키며, 배터리를 먹는다. 그래서 더 적은 에너지로 같은 일을 하려는 연구가 오래 이어져 왔다.

항빈맥 페이싱(ATP). 충격이 아니라 페이싱이다. 재진입 회로에 흥분성 간극이 있으면 빈맥 주기보다 조금 짧은 주기로 몇 발 페이싱해 회로 안에 파를 하나 더 밀어 넣고, 그 파가 회로를 지운다. 통증이 없고 에너지가 수천 분의 일이다. 다만 간극이 있는 해부학적 재진입성 빈맥에만 통하고, 기능적 재진입인 세동에는 통하지 않는다. 회로 구조가 곧 치료 선택을 결정한다는 점이 부정맥의 분류와 그대로 이어진다.
고정 해제(unpinning). 나선 팁이 흉터나 혈관 같은 비흥분 장애물에 붙어 있으면 안정해져 잘 안 없어진다. 나선 주기에 맞춘 약한 주기 자극으로 팁을 표류시켜 장애물에서 떼어내면 그 뒤로는 훨씬 쉽게 소멸한다. 나선파의 공명 표류가 그대로 치료 아이디어가 된 사례.
저에너지 항세동 페이싱(LEAP). 조직 안의 자연스러운 이질성 — 혈관, 섬유 방향 불연속 — 하나하나가 약한 전기장에서도 가상전극이 된다는 점을 이용한다. 강한 충격 한 발 대신 약한 펄스를 여러 발 연달아 주면 이 분산된 가상전극들이 조직을 단계적으로 동기화한다. 2011년에 보고된 결과에서 필요한 에너지가 기존의 약 1/6 수준으로 떨어졌다. 여전히 실험·모형 단계의 접근이지만, “때리는 세기”가 아니라 “타이밍과 조직의 미세구조”가 지렛대라는 관점 전환이라는 점에서 계산 심장학이 임상에 던진 가장 구체적인 제안으로 꼽힌다.

7. 수치적으로 무엇이 어려운가[편집]

이중영역 시뮬레이션은 심장 계산 중에서도 유난히 비싸고 까다롭다.

결합된 포물형 + 타원형. 매 시간 스텝마다 $\phi_e$ 에 대한 큰 타원형 방정식을 풀어야 한다. 실행 시간의 대부분이 여기 들어가고, 그래서 전처리기 선택이 성능의 전부다. 계수가 강하게 이방적이고 위치에 따라 주축이 회전하므로 기하 다중격자가 잘 안 듣고, 대수 다중격자법을 전처리기로 얹은 켤레기울기법이 사실상 표준이다.
특이한 타원형 계. 경계가 전부 절연(노이만) 조건이면 $\phi_e$ 는 상수만큼의 불확정성을 가지므로 계 행렬이 특이하다. 평균 0 조건을 걸거나 한 점을 접지해야 하고, 우변이 영공간에 직교해야 한다 — 물리적으로는 주입한 전류의 총합이 0이어야 한다는 조건이다. 전극 하나만 놓고 돌리면 여기서 조용히 발산한다.
충격 구간의 시간 해상도. 평소 $\Delta t$ 가 10~50 µs 규모라면, 충격 중에는 막전위가 수십 µs 안에 수백 mV 움직이므로 1 µs 수준까지 내려야 한다. 충격 전후로 스텝을 바꾸는 적응 전략이 필수.
이온 모형이 외삽 영역에 들어간다. 세포 모형의 전류 방정식은 대략 $-100$ ~ $+50$ mV 범위의 실측에 맞춘 것인데, 충격 중 $V_m$ 은 $+300$ mV나 $-400$ mV까지 간다. 그 영역에서 모형은 검증된 적이 없고, 지수 항이 그대로 폭주해 수치적으로 터지기도 한다. 실무에서는 전기천공(electroporation) 전류 항을 추가해 극단 전위에서 막이 새는 물리를 넣어 준다 — 이것 자체가 검증하기 어려운 추가 가정이라, 제세동 시뮬레이션 결과의 불확실성에서 큰 몫을 차지한다.³
도구. openCARP, Chaste 등 이중영역을 지원하는 오픈소스 코드가 있고, 대규모 실행은 대부분 MPI + GPU 컴퓨팅 혼합이다. 반응항은 격자점별로 완전히 독립이라 잘 흩어지지만, 타원형 풀이의 전역 통신이 확장성의 천장을 만든다. 전형적인 영역 분할법 문제.
검증. 표면 광학 매핑과 심전도가 비교 대상이지만, 둘 다 표면 또는 원거리 관측이라 심부의 3차원 필라멘트를 직접 확인해 주지 못한다. 모형이 만든 예쁜 그림 중 어디까지가 검증된 것인지 구분하는 것이 이 분야의 상시 숙제이며, 검증 및 확인의 정신이 특히 절실한 자리다.

8. 관련 문서[편집]

이 계보가 유쾌한 이유는 이론적 예측이 “직관과 반대라서” 한동안 의심받았다는 데 있다. 전극 하나 꽂았는데 옆에 과분극이 생긴다는 말을 실험가들이 곧이곧대로 믿을 리 없었고, 실제로 개뼈다귀가 눈으로 보이기 전까지는 수치 아티팩트 취급을 받았다. 지금은 그 반대로, 개뼈다귀가 안 나오는 이중영역 코드는 이방성 비를 잘못 넣은 것으로 의심받는다. ↩
“역치”라는 이름이 붙었는데 실은 분포의 대표값이라는 상황은 불응기에서도 똑같이 나온다. 흥분성 계에서 “역치”라 불리는 양은 대체로 결정론적 상수가 아니라 프로토콜과 잡음에 의존하는 통계량이다. 보고서에 숫자 하나만 적고 프로토콜을 안 적으면 그 숫자는 재현이 안 된다. ↩
그래서 제세동 시뮬레이션의 결과를 볼 때는 “충격을 몇 J로 주었나”보다 “그 코드가 $V_m = +200$ mV 에서 무엇을 하도록 되어 있나”를 먼저 물어보는 것이 맞다. 논문 본문에 없고 부록 어딘가에 한 줄로 적혀 있거나, 아예 안 적혀 있는 경우가 흔하다. 검증되지 않은 외삽이 결론을 좌우하는 전형적 구조라, 불확실성 정량화 관점에서 보면 손대야 할 곳이 명확한 편이다. ↩

제세동