전염병 동역학의 보존법칙과 ‘역학적 모멘텀’ 새로운 통합 개념

초록

저자들은 감염 가능성을 가중한 유행병 ‘모멘텀’을 정의하고, 이를 통해 모든 전염병 모델이 동일한 위상평면 상에서 보존량을 따라 움직인다는 수학적 법칙을 발견했다. 모멘텀은 관측된 발병률과 세대 간격 분포만으로 계산 가능하며, 전염성(R₀)과 사전 면역 수준을 동일 시계열에서 별도로 추정할 수 있게 한다. 1918년 인플루엔자 사례를 재분석해 이 방법의 실용성을 보여준다.

상세 분석

이 논문은 기존 SIR·SEIR 등 다양한 전염병 모델을 하나의 통합 프레임워크로 묶는 ‘역학적 모멘텀( epidemic momentum)’이라는 새로운 변수 도입으로 시작한다. 모멘텀 Y는 현재 감염자 수에 그들의 남은 전염 가능성을 가중한 가중합으로 정의되며, 수식적으로는 Y(τ)=∫₀^∞ι(τ−α)R_α/R₀ dα 혹은 Y(τ)=ι(τ)−∫₀^∞ι(τ−α)g(α)dα 로 표현된다. 여기서 ι는 발병률, g(α) 는 세대 간격 확률밀도, R_α는 감염 연령 α에서의 감소된 재생산수이다. 이 정의는 Renewal 방정식 체계와 완전히 호환되며, SIR·SEIR에서는 Y가 바로 유병률과 동일함을 보여준다.

모멘텀을 도입함으로써 얻어지는 가장 중요한 결과는 ‘보존량 C(x,y)’이다. Susceptible 비율 X와 모멘텀 Y 사이의 위상관계 dY/dX=−1+ĥx/X 를 적분하면
 Y(x)=y_i+(x_i−x)−ĥx·ln(x_i/x)
이라는 식이 나오고, 이는 모든 모델이 동일한 곡선(등위선) 위를 따라 움직인다는 것을 의미한다. 여기서 ĥx=1/R₀이며, C(x,y)=Y+ĥx·V(x/ĥx) (V(u)=u−1−ln u) 로 정의되는 보존량은 초기 조건에 따라 고정된다. 따라서 전염병 진행 과정은 이 보존량에 의해 강제되며, 모델 복잡도는 단지 시간축의 재스케일링(τ→T(τ))에만 영향을 미친다.

보존량을 이용하면 사전 면역 수준 z₋=1−x₋ (x₋는 역학적 궤적의 왼쪽 교점)와 최종 감염 규모 z₊=x₋−x₊ (x₊는 오른쪽 교점)를 직접 계산할 수 있다. 기존의 최종 규모 공식은 사전 면역을 외부 입력으로 가정했지만, 여기서는 R₀만 알면 x₋, x₊를 역으로 구해 z₋와 z₊를 추정한다. 또한, 발병 초기와 말기의 지수 성장·감소율 λ₋, λ₊는 Laplace 변환 L