공생 SIR 모델에서 동시 감염 회복 속도의 영향

초록

본 논문은 네트워크 상에서 두 질병이 동시에 전파되는 대칭 SIR 공생 모델을 제시하고, 공동 감염자들의 회복률이 단일 감염자와 다를 때 전염병의 일시적 특성이 어떻게 변하는지를 분석한다. 평균장(mean‑field) 근사와 전이 대칭성을 이용해 차원 축소를 수행하고, 회복률 감소가 공동 감염 부담을 증가시키며 전염 기간과 피크를 연장한다는 정량적 결과를 도출한다.

상세 분석

논문은 먼저 기존의 2종 공생 접촉 과정(2SCP)을 확장하여, 두 질병 A와 B가 동시에 퍼지는 SIR 모델을 정의한다. 각 개인은 (X,Y)∈{s,i,r}² 형태의 9가지 상태를 가지며, 전염률 λ₁, λ₂와 회복률 μ를 기본 파라미터로 둔다. 공동 감염(i,i) 상태에서는 회복률을 μ̄<μ로 설정함으로써 ‘공생 회복’ 효과를 도입한다. 평균정도 ⟨k⟩를 갖는 무작위 네트워크에 대한 동질 평균장 방정식을 전개하고, 질병 자유 평형(DFE)에서의 선형화로부터 차례대로 전염 매트릭스 F와 전이 매트릭스 V를 구성한다. 차례대로 F와 V의 야코비안을 이용해 기본 재생산수 R₀=ρ(FV⁻¹)=max(R₀,A,R₀,B)임을 증명한다. 이는 두 질병이 독립적인 SIR 프로세스로 동작함을 의미하지만, 공동 감염 상태가 존재함으로써 비선형 단계에서 중요한 영향을 미친다.
특히 λ₁=λ₂=λ인 대칭 전염 상황에서는 (i,s)와 (s,i) 등 상태를 교환해도 시스템이 불변임을 보이고, 초기 조건이 대칭을 만족하면 해는 6차원 불변 부분공간으로 축소된다. 이 축소된 변수 a=(i,s)=(s,i), b=(r,s)=(s,r), d=(r,i)=(i,r), u=(s,s), c=(i,i), v=(r,r) 로 정의하고, 총 감염량 M(t)=a(t)+c(t)+d(t)와 총 감염자 수 J(t)=2a+2d+c 로 표현한다.
주요 정리는 두 가지이다. 첫째, c′(t)=2⟨k⟩λa(t)M(t)−μ̄c(t) 형태에서 μ̄가 감소하면 c(t)와 누적 부담 B(μ̄)=∫₀^∞c(t)dt가 단조 증가함을 보이는 단조성 정리와 그 증명(Prop. 1)을 제시한다. 둘째, c′≥−μ̄c를 이용해 c(t)≥c(t₀)e^{−μ̄(t−t₀)} 를 얻고, 이를 통해 전염 기간 하한 T_J(ε)≥t₀+μ̄^{-1}ln(c(t₀)/ε) 를 도출한다. 즉 μ̄가 작을수록 전염 기간이 역비례적으로 연장된다.
수치 실험에서는 ⟨k⟩=10, λ=0.25, μ=1을 기준으로 μ̄를 0.2~1 범위에서 변화시켜, μ̄ 감소가 (i) 공동 감염 비중, (ii) 전체 감염 피크, (iii) 전염 지속시간을 모두 증가시키는 것을 확인한다. 민감도 방정식 Σ′=Df·Σ+∂{μ̄}f 를 풀어 ∂{μ̄}J(t)=2σ_a+σ_c+2σ_d 가 전염 피크 근처에서 음수임을 보이며, 이는 회복률 감소가 피크를 상승시킨다는 1차 근사 결과와 일치한다.
결론적으로, 공동 감염 회복률 μ̄는 전염병 역학에서 무시할 수 없는 파라미터이며, 특히 고전파 전파 상황에서 전염 기간과 피크를 크게 확대한다는 점을 이론적·수치적으로 입증한다.