[論文レビュー] A reaction-diffusion system to better comprehend the unlockdown: Application of SEIR-type model with diffusion to the spatial spread of COVID-19 in France

Youcef Mammeri|arXiv (Cornell University)|May 6, 2020

COVID-19 epidemiological studies参考文献 24被引用数 47

ひとこと要約

本論文は、無症状者と移動人口を含む反応拡散型の SEIR 型モデルを開発し、France における COVID-19 の空間拡散を研究し、データからパラメータを推定し、R0 を計算し、解除ロックダウン戦略を分析する。

ABSTRACT

A reaction-diffusion model was developed describing the spread of the COVID-19 virus considering the mean daily movement of susceptible, exposed and asymptomatic individuals. The model was calibrated using data on the confirmed infection and death from France as well as their initial spatial distribution. First, the system of partial differential equations is studied, then the basic reproduction number, R0 is derived. Second, numerical simulations, based on a combination of level-set and finite differences, shown the spatial spread of COVID-19 from March 16 to June 16. Finally, scenarios of unlockdown are compared according to variation of distancing, or partially spatial lockdown.

研究の動機と目的

• COVID-19 の空間拡散の理解と移動性・距離確保の影響の理解を動機づける。
• 移動可能な感受性者、暴露者、無症状者を組み込んだ結合型 PDE/ODE SEIR 型フレームワークを開発する。
• フランスの確認済み症例と死亡データにモデルをキャリブレーションし、基本再生数 R0 を導出する。
• 数値シミュレーションを通じて空間伝播を分析し、解除戦略を評価する。

提案手法

• 感受性者 S、暴露 E、無症状 I_a 部分を拡散を伴う六区分 SEI_aI_sUR モデルを定式化する。
• lockdown/解除の影響をモデル化するために、時間変化する接触率 ω(t) と拡散 d(t) を組み込む。
• 拡散なしの次世代行列から R0 を導出し、指数成長の条件を確立する。
• 適切なパラメータ領域でグローバルに適切性を示し、感染症なしの挙動を解析する。
• 数値シミュレーションにはレベルセット法と有限差分法、Runge-Kutta 時間積分を用いる。
• Approximate Bayesian Computation を用いた非線形最小二乗法と準ニュートン法のステップで六つのパラメータ（rho, beta_e, beta_s, beta_a, p, mu）を較正する。

実験結果

リサーチクエスチョン

• RQ1 感受性者・暴露者・無症状者の空間拡散は France における COVID-19 の拡散にどのように影響するか？
• RQ2 lockdown と解除の政策は基本再生数と空間的分布にどのような影響を与えるか？
• RQ3 モデルは観測された感染・死亡推移を再現し、地図ベースの解除戦略を提供できるか？
• RQ4 データに最もよく適合するパラメータ値は何で、それらは暴露および無症状者からの伝播について何を意味するか？

主な発見

• モデルは較正の下で基本再生数 R0 = 3.425257 を得る。
• lockdown は実効再生数を約 3.42 から 0.38 に低下させる。
• 64% 以上の距離確保順守は、症候性感染の数を制御し、モデル領域で Reff を 1 未満に保つことができる。
• 解除戦略は空間的不均質を示し、特定の戦略下では東部地域がより影響を受け続ける。
• 介入なしでは、全国的に感染レベルが広く高くなると予測され、移動と地域的な制御の役割を強調する。
• 実効再生数 Reff(x,t) は S(x,t)/N(x,t) と接触率 ω(t) に依存し、その域内平均は流行の抑制を追跡する。