다중균주 전염병 모델에서 일시적 교차 면역을 포함한 새로운 혼돈 끌개

초록

본 연구는 일시적 교차 면역을 고려한 다중균주 전염병 모델에서, 2차 감염 시 감염성 감소(‘역 ADE’) 구간에서도 결정론적 혼돈 끌개가 존재함을 수학적으로 입증한다. 기존에는 2차 감염 시 감염성이 크게 증가하는 ADE 조건에서만 혼돈이 보고되었으나, 본 논문은 보다 현실적인 역 ADE 파라미터 영역에서도 광범위한 혼돈 영역을 발견한다. 이는 특히 뎅기열과 같이 2차 감염 시 입원율이 높아 실제 전파 기여도가 감소하는 경우에 적용 가능하며, 다중균주 전염병 모델 전반에 걸친 새로운 분석 가능성을 제시한다.

상세 분석

논문은 기본적인 SIR‑type 다중균주 모델에 일시적 교차 면역(temporary cross‑immunity)과 2차 감염 시 감염성 변화를 매개변수 β₂로 도입한다. β₁은 1차 감염 시 전파율, β₂는 2차 감염 시 전파율이며, β₂/β₁>1이면 전통적인 ADE, β₂/β₁<1이면 역 ADE라 정의한다. 교차 면역 기간 τ는 감염 후 일정 기간 동안 다른 균주에 대한 감수성을 완전히 차단한다는 가정이다. 모델은 4차 비선형 미분방정식 체계로 구성되며, 파라미터 공간을 전역적으로 탐색하기 위해 연속적인 파라미터 스위핑과 Poincaré 섹션, 최대 Lyapunov 지수 계산을 수행한다.

주요 결과는 β₂/β₁<1인 역 ADE 구간에서도 최대 Lyapunov 지수가 양수(≈0.020.15)인 영역이 존재한다는 점이다. 이는 전통적으로 혼돈이 발생하려면 β₂/β₁≫1이어야 한다는 기존 연구와는 대조적이다. 특히 τ가 3090일 사이일 때 혼돈 영역이 크게 확대되며, τ가 짧아질수록(≤15일) 혼돈이 사라지고 주기적 궤도로 수렴한다. 이는 교차 면역이 일정 기간 동안 감염자 풀을 감소시켜 비선형 피드백을 강화하고, 감염성 감소 효과와 결합될 때 복잡한 동역학을 야기한다는 메커니즘을 시사한다.

또한, 분기 다이어그램을 통해 주기‑2, 주기‑4, 그리고 혼돈 구간이 교차하는 복합적인 구조를 확인하였다. Hopf 분기점은 β₂/β₁≈0.8에서 발생하고, 그 이후 연속적인 period‑doubling cascade가 진행되어 혼돈으로 이어진다. 모델의 민감도 분석에서는 감염 지속 기간(γ⁻¹)과 회복률이 혼돈 발생에 큰 영향을 미치며, 특히 회복률이 감소하면 혼돈 영역이 확대된다.

이러한 수학적 결과는 뎅기열과 같은 실제 질환에 직접적인 함의를 가진다. 2차 감염 시 입원율이 높아 실제 전파 기여도가 감소하는 상황을 역 ADE 파라미터로 표현할 수 있으며, 일시적 교차 면역은 실제로 2~3개월 정도 지속되는 것으로 알려져 있다. 따라서 모델이 제시하는 혼돈 영역은 실제 역학 데이터와 일치할 가능성이 높으며, 기존에 혼돈을 설명하기 위해 가정했던 과도한 ADE 가정을 완화할 수 있다.