[논문 리뷰] Modeling, dynamics and optimal control of Ebola virus spread

Amira Rachah, D. F. Torres|arXiv (Cornell University)|2016. 03. 18.

Viral Infections and Outbreaks Research참고 문헌 18인용 수 4

한 줄 요약

이 연구는 실존 WHO 데이터를 통합하여 기저귀니아에서 에볼라 전파를 분석하기 위해 생명동역학과 유도 사망률을 고려한 수정된 SIR 모델을 개발한다. 최적 제어 이론을 적용하여 감염자 수와 백신 접종 비용을 최소화하며, 표적으로 삼은 백신 접종이 통제되지 않은 전파에 비해 피크 감염자 수를 65% 감소시키고 전염병 지속 기간을 15일 이상 단축시킨다.

ABSTRACT

We present a mathematical analysis of the early detection of Ebola virus. The propagation of the virus is analysed by using a Susceptible, Infected, Recovered (SIR) model. In order to provide useful predictions about the potential transmission of the virus, we analyse and simulate the SIR model with vital dynamics, by adding demographic effects and an induced death rate. Then, we compute the equilibria of the model. The numerical simulations confirm the theoretical analysis. Our study describes the 2015 detection of Ebola virus in Guinea, the parameters of the model being identified from the World Health Organization data. Finally, we consider an optimal control problem of the propagation of the Ebola virus, minimizing the number of infected individuals while taking into account the cost of vaccination.

연구 동기 및 목표

• 기본 인구 역학과 유도 사망률을 통합함으로써 에볼라 전파 모델링의 현실성 향상.
• 정확한 질병 전파 예측을 위한 생명동역학 하에서 SIR 모델의 평형점과 안정성 분석.
• 감염자 수를 최소화하면서 백신 접종 비용을 고려한 최적 제어 전략 설계 및 평가.
• 2015년 기니의 에볼라 유행 보고서를 바탕으로 한 데이터 기반 시뮬레이션을 통한 공중보건 정책에 실질적인 통찰 제공.

제안 방법

• 출생 및 사망률을 포함하고 감염자에게서 기인하는 유도 사망률까지 고려한 수정된 SIR 모델을 수립.
• 감염자 수를 줄이는 데 기여하는 백신 접종 노력에 해당하는 제어 변수를 도입.
• 감염자 수와 백신 접종 비용을 종합한 비용 기능을 최소화하기 위해 최적 제어 이론 적용.
• 두 가지 독립적인 수치적 방법을 사용하여 최적 제어 문제 해결: ACADO(국소 최적화) 및 BocopHJB(동적 프로그래밍을 통한 전역 최적화).
• ACADO와 BocopHJB의 해를 상호 교차 검증하여 전역 최소값으로 수렴함을 확인.
• 기니의 2015년 에볼라 유행에서 확보한 WHO 데이터를 사용하여 모델의 매개변수를 校정하고 예측의 타당성을 검증하기 위해 수치 시뮬레이션 수행.

실험 결과

연구 질문

• RQ1 생명동역학과 유도 사망률을 통합함으로써 에볼라 전파에 대한 SIR 모델의 안정성과 거동에 어떤 영향을 미치는가?
• RQ2 최적의 백신 접종 제어 전략이 피크 감염자 수와 전염병 지속 기간에 어떤 영향을 미치는가?
• RQ3 통제되지 않은 전파와 비교했을 때 최적 제어 전략은 감염 감소와 인구 생존율에 어떤 영향을 미치는가?
• RQ4 국소 최적화와 전역 최적화 방법을 모두 사용하여 최적 제어 해를 신뢰성 있게 계산할 수 있는가, 이는 해의 강건성 보장에 기여하는가?
• RQ5 모델의 예측은 2015년 기니의 실제 에볼라 유행 데이터와 어느 정도 일치하는가?

주요 결과

• 생명동역학과 유도 사망률을 통합한 모델은 기니에서의 에볼라 전파 역학을 더 현실적으로 표현한다.
• 수치적 시뮬레이션은 이론적 안정성 분석을 확인하며, 질병 없음 및 만성 감염 평형점이 시스템의 거동과 일치함을 보여준다.
• 최적 제어 하에서 피크 감염자 수는 통제되지 않은 경우 대비 16.3%로 감소하며, 통제되지 않은 경우 47.69%에 비해 크게 감소한다.
• 최적 백신 접종 제어 하에서 전염병 지속 기간은 60일 이상에서 46일로 단축된다.
• 제어 조건 하에서 감염 가능 인구가 더 빠르게 감소하여 캠페인 종료 시점에 34.78%에 도달하지만, 통제되지 않은 경우는 44.12%에 머무른다.
• 제어 조건 하에서 회복자 수가 더 빠르게 증가하여 빠른 면역 확보와 전파 감소 효과를 나타낸다.