發刊日期 |
2022年6月
|
---|---|
標題 | 數學模型論新冠清零及「與疫共存」之兩難困境 |
作者 | |
關鍵字 | |
檔案下載 | |
全文 |
先前《就新冠疫情數學建模談起》一文 (A) 反應機制(I)--- 針對一般疫情傳播狀況而言: $$X\,{\buildrel \beta Y\over \longrightarrow}\, Y\,{\buildrel \nu' \over \longrightarrow} Z,\quad Y\,{\buildrel \emptyset \over \longrightarrow}\, X,$$其中 $\nu =\nu'+\emptyset \gt\emptyset,$ 事實上, 目前台灣因COVID-19有效控制疫情, 死亡人數相對於總人口數甚微情況下, 可設定系統內總人口數為定值: $X+Y+Z=N_0$, 如此三聯立正微分方程可由此關係式降階為 $2\times 2$ 之正微分方程組, 因此可建立之疫情流行系統方程式如下: \begin{align*} \frac{dX}{dt}=\,&-\beta XY+\emptyset Y+m=f_1(X,Y),\qquad X(t=0)=X_i,\\ \frac{dY}{dt}=\,&\beta XY-\nu Y+m'=g_1(X,Y),\hskip 16pt\qquad Y(t=0)=Y_i\gt0, \end{align*}其中 $X, Y, Z$ 分別代表未受感染者 (Susceptibles)、 已被感染者(Infectious)、 自免疫系統去除者(Recovered) (如終生免疫者跟死亡者)。 而 $\beta ,\emptyset ,\nu ,\nu'$ 分別代表疾病傳播常數, 疫病復原回復後可再被感染 (Reinfection) 速率常數, 感染者總去除速率常數, 感染後得以免疫速率常數。 $m,m'$ 分別表因防疫破口, 無法有效阻斷, 所引入之未感染者速率及已被感染者速率。 切記: 所有瘟疫之結局皆是所有人口遠離流行疫情系統, 才能完全解封來過正常生活。 唯一差別在於疫情控制成功與否, 是Z將以不幸亡者身分離世、 抑或是活著恢復健康著離開此流行病學系統。 換言之, 防疫成功的目標指數是尋求在最少死亡人數且最多免疫者條件下, 脫離此疫情陰霾。 當 $\frac{dX}{dt}=\frac{dY}{dt}=0$ 條件下, 其平衡穩態滿足 \begin{align} \frac{dX}{dt}=\,&-\beta XY+\emptyset Y+m=0,\label{1}\\ \frac{dY}{dt}=\,&\beta XY-\nu Y+m'=0,\label{2} \end{align}由 $\eqref{1}+\eqref{2}$ 可得到 $Y_0=\dfrac{m+m'}{\nu -\emptyset}$, 代入 \eqref{2} 可得 $$X_0=\frac{\nu}{\beta} -m'\frac{\nu -\emptyset}{m+m'}=\frac{\nu m+\emptyset m'}{\beta (m+m')}.$$亦即此時平衡穩態 $SS$ (steady state) 是 $$(X_0,Y_0)= \Big(\frac{\nu m+\emptyset m'}{\beta (m+m')},\frac{m+m'}{\nu -\emptyset}\Big).$$
依據此系統條件來決定穩態 ($SS$) 特性及分析其穩定性 (Stability) 此系統方程組之特徵值(eigenvalues)可由如下行列式來決定 \begin{align} \det({\underline A}-\lambda {\underline I})=\left|\begin{array}{ccc} -\beta Y_0-\lambda &&\emptyset-\beta X_0\\ \beta Y_0&&\beta X_0-\nu -\lambda\end{array}\right|=\,&0,\nonumber\\ {\hbox{亦即 }} \lambda ^2+(\nu +\beta Y_0-\beta X_0 )\lambda +\beta Y_0 (\nu -\emptyset)=\,&0,\label{5} \end{align}由此可決定出其特徵值 $\lambda_1,\lambda_2$ 必須滿足上式, 並且以下式表之: $$\lambda^2-(\lambda_1+\lambda_2)\lambda+\lambda_1\lambda_2=(\lambda-\lambda_1)(\lambda-\lambda_2)=0.$$若系統可穩定收斂至穩態 $(X_0,Y_0)$ 下, 則兩特徵值皆必須滿足具有負實部之先決條件下, 亦即必須滿足 $$\hbox{(i)}\quad \beta(Y_0-X_0)+\nu \gt0 ,\qquad \hbox{(ii)}\quad \nu-\emptyset\gt0.$$事實上, 由於 $\nu \gt\emptyset$ 恆成立, 因此即使是不再清零的疫情系統中, 若仍要維持於長期疫情控制收斂穩定之下, 以達到疫情在完全可控制狀態下, 顯然其必要條件是滿足 $X_0 \lt Y_0+\frac{\nu}{\beta}$ (此不等式之確定, 亦即表示病毒溯源傳播疫調有其絕對之必要性。 任何時間未感染人數、 已感染人數與平衡 $X_0,Y_0$ 之動態值以及各病毒變異株相關之動力學參數 $\nu ,\beta$ 皆須確實全程掌握, 在此定義: $Y_0+\frac{\nu}\beta$ 為「關鍵」族群數目)。 若無法有效持續追蹤, 清零之可能目標即已無可能。 單由此點可知, 台灣在一開始即已在吹哨人提醒下, 當機立斷超前部署進行出入境管制戒備 (才有機會能設定清零 $Y_0=0$ 為目標), 而能在世界防疫戰中取得可掌可控之領先地位。 反觀世界一些國家早在疫情爆發初期, 即錯失良機導致國門失守, 甚至因未知傳播鏈過多, 即已然放棄疫調, 從此就被迫注定走入可能必須「與病毒共存」的道路上。 防疫成功與否, 在於是否能控制未被感染人口數目始終小於「關鍵」族群數目: $Y_0+\frac{\nu}\beta $, 以達到疫情可調可控之狀況下。 由上述訊息更可指出:現今若能與有效疫苗覆蓋率合併計算, 即可確實推論出是否疫情系統已在可掌可控的發展情況下, 漸次步出瘟疫陰霾? 再者, 此系統平衡穩態 $(X_0,Y_0)$ 先天存在著三個令人感到頗有意義的極限疫情情況, 在此分別點出來進行分析如下: (i) 第一個極限穩態: $\lim\limits_{m\to 0} (X_0,Y_0 )=(\frac{\emptyset}{\beta },\frac{m'}{\nu'})$, 代表即使採取更嚴格之邊境隔離檢疫管制, 若在人口中存在隱形帶原無症狀感染者存在, 則會發現系統中永遠無法達到零確診之情況, 此亦即是不少國家日日出現上萬人確診的可能原因所在。 (ii) 第二個極限穩態: $\lim\limits_{m'\to 0}(X_0,Y_0 )=(\frac{\nu}\beta ,\frac{m}{\nu'})$, 說明在實施完全有效入境邊境篩檢, 來剔除受到感染者, 卻不完全封城, 仍可自由出入境, 可能亦無以避免在人口群眾中始終有確診者存在。 台灣目前對此種可能發生之機會仍要格外小心注意。 (iii) 第三個極限穩態: $\lim\limits_{m\to 0,m'\to 0}(X_0,Y_0 )=(\frac{\nu}\beta ,0)$ 或 $(\frac{\emptyset}\beta ,0)$, 代表只有在完全鎖國封城之特殊條件下, 才能達到所謂不沾鍋的「完全零確診」結果。
上述三種極限穩態皆說明若無完全封關、 封城或是鎖國條件下, 「社區感染」或是「與病毒共存」的發生機率會是始終存在的。
時中部長所言實有所本。
一旦出現有破口未能有效堵絕, 未知傳播鏈隱形存在下, 如此即造成無法有效隔絕所有可能傳播感染機會, 亦代表疫情走上「與病毒共存」之機會即難以避免。
無可諱言地, 自始至終皆未曾完全封城或封關的台灣疫情而言, 所言之「與病毒共存」的危機事實上是處處存在、而且亦是在所難免。
但即使如此, 若能時時掌握疫調, 全民配合防疫生活確實執行, 疫情可掌可控仍是一可達到的目標。
此三種極限穩態更清楚詮釋前文 新近由於各種疫苗施打逐漸覆蓋, 究竟會對流行疫情系統內族群人口產生何種效應? 以下即是就「有效疫苗覆蓋」來探討疫情系統穩態平衡狀況, 更對何時才能「完全解封」, 「自由出入國際旅遊」做出之沙盤推演說明如下: (B) 反應機制(II)--- 針對有效疫苗被發展出來進行施打注射覆蓋(亦即因有效疫苗注射施打而免疫: $X\rightarrow Z$)之疫情傳播狀況(施行有效疫苗注射施打, 即可離開此疫情系統)而言: $$X\,{\buildrel \beta Y\over \longrightarrow}\, Y\,{\buildrel \nu' \over \longrightarrow} Z,\quad Y\,{\buildrel \emptyset \over \longrightarrow}\, X,\quad X\,{\buildrel \Psi \over \longrightarrow}\, Z,$$其中 $\nu =\nu'+\emptyset \gt\emptyset,$ 系統中所有人口數可為 $X+Y+Z=N_0$, 因此在可掌可控之疫情系統所建立之流行病學方程式下, 當 $\frac{dX}{dt}=\frac{dY}{dt}=0$ 之平衡穩態下, 須滿足: \begin{align*} \frac{dX}{dt}=\,&-\beta XY+\emptyset Y-\Psi X=0,\qquad X(t=0)=X_i,\\ \frac{dY}{dt}=\,&\beta XY-\nu Y=0,\hskip 50pt \qquad Y(t=0)=Y_i\gt0,\\ \frac{dZ}{dt}=\,&\nu' Y+\Psi X. \end{align*}由上述方程式可解出兩個穩態 $(X_0,Y_0 )$ 分別為 $SS_1(0,0)$, $SS_2(\frac \nu\beta ,\frac{\Psi\nu}{(\emptyset-\nu)\beta}$)。 但第二個穩態並不存在 ($Y_0 \lt 0$)。 因此 $SS_1(0,0)$ 為唯一存在之穩態。 代入線性方程組可得 \begin{align*} \frac{d}{dt} \Delta X=\,&(-\beta Y_0-\Psi)\Delta X+(\emptyset-\beta X_0 )\Delta Y,\tag*{($3'$)}\\ \frac{d}{dt} \Delta Y=\,&\beta Y_0 \Delta X+(\beta X_0-\nu )\Delta Y.\tag*{($4'$)} \end{align*}取定唯一存在之穩態 $SS_1(0,0)=(X_0,Y_0 )$ 代入可得 \begin{align*} \frac{d}{dt} \Delta X=&-\Psi\Delta X+\emptyset \Delta Y,\tag*{($3'$-1)}\\ \frac{d}{dt} \Delta Y=&-\nu \Delta Y.\tag*{($4'$-1)} \end{align*}$\det(\underline A-\lambda \underline I)=0$可得兩個特徵值分別皆為負數的 $-\Psi$ 與 $-\nu$。 解析解為 $$(\Delta X,\Delta Y)=\Big(\Delta X_i e^{-\Psi t}+\frac{\emptyset\Delta Y_i}{\Psi-\nu}(e^{-\nu t}-e^{-\Psi t}),\Delta Y_i e^{-\nu t}\Big),$$當達到平衡穩態時 $\lim\limits_{t\to \infty}\Delta X=\lim\limits_{t\to \infty}\Delta Y=0$。 換言之, 在有效疫苗注射覆蓋條件下, 平衡後任何時間皆維持在遠離此瘟疫的狀況, 至此不再有人存在此疫情系統中 (亦即完全解封正常生活)。 亦即是 $X(t)=X_0$, $Y(t)=Y_0=0$。 此結果即是前述的完全脫離瘟疫的平衡穩態 $SS_1(0, 0)$。 然而反觀, 另一完全與此反應機制背道而馳的例子則是在疫情爆發初期, 由於並無有效防疫方法及策略, 因此除了個人衛教自我防護及加強醫療備援外, 對政府下達之公共衛生之政策又應當如何因應? 以下即是所討論的狀況分析: (C) 反應機制(III)--- 針對有效疫苗未被發展出來, 因此無法產生有效施打預防疫苗注射覆蓋, 而且即使是受感染後康復, 仍會再次被感染狀況而言: $$X\,{\buildrel \beta Y\over \longrightarrow}\, Y\,{\buildrel \nu' \over \longrightarrow} Z,\quad Y\,{\buildrel \emptyset \over \longrightarrow}\, X,\quad Z\,{\buildrel \Psi' \over \longrightarrow}\, X,$$其中 $\nu =\nu'+\emptyset \gt\emptyset,$ 系統內總人口數為 $X+Y+Z=N_0$。 由此可推得其疫情系統微分方程組如下: \begin{align*} \frac{dX}{dt}=\,&-\beta XY+\emptyset Y+\Psi'X,\qquad X(t=0)=X_i,\\ \frac{dY}{dt}=\,&\beta XY-\nu Y,\hskip 56pt\qquad Y(t=0)=0,\\ \frac{dZ}{dt}=\,&\nu' Y-\Psi' Z,\hskip 61pt \qquad Z(t=0)=0. \end{align*}將總人口數 $N_0$ 代入, 原式可推導為 \begin{align*} \frac{dX}{dt}=\,&-\beta XY+\emptyset Y+\Psi'(N_0-X-Y)=0,\quad X(t=0)=X_i,\\ \frac{dY}{dt}=\,&\beta XY-\nu Y=0,\hskip 130pt\qquad Y(t=0)=0. \end{align*}由上述兩個方程式可推導解出兩個穩態 $(X_0,Y_0 )$ 分別為 $SS_1(N_0,0)$, $$SS_2\Big(\frac\nu\beta,\frac{\Psi'(N_0-\frac\nu\beta)}{\Psi'-\emptyset+\nu}\Big)=\Big(\frac\nu\beta,\frac{\Psi'(N_0-\frac\nu\beta)}{\Psi'+\nu'}\Big).$$對於 $SS_1(N_0,0)$, 即是台灣一直希望能將疫情控管在零確診 (清零) 的平衡穩態狀況下。 然而 $$SS_2\Big(\frac\nu\beta,\frac{\Psi'(N_0-\frac\nu\beta)}{\Psi'+\nu'}\Big)$$即是現今世界不少國家所提到之「與病毒共存」(例如:英國所提之「佛系防疫」)的平衡穩態。 而且在此平衡穩態下, 在整個人口族群中, $Y_0:Z_0$ 亦會呈現幾乎固定比例(亦即 $\frac{\Psi'}{\nu'}$)。 換言之, 當疫情開始爆發, 即已經無法有效掌握「零號病人」去向, 無法有效疫調追蹤隱形傳播鏈時, 此時就算後續疫情控制得宜, 可能「與病毒共存」亦是在所難免 (例如: 目前新聞報導提及不少國家, 日日皆有上萬人確診者, 即是實際發生之案例)。 以下即是證明滿足此兩個平衡穩態 $$SS_1(N_0,0),\qquad SS_2\Big(\frac\nu\beta,\frac{\Psi'(N_0-\frac\nu\beta)}{\Psi'+\nu'}\Big)$$之必要條件之推導過程。 首先再次代入線性方程組可得 \begin{align*} \frac{d}{dt} \Delta X=\,&(-\beta Y_0-\Psi')\Delta X+(\emptyset-\Psi'-\beta X_0 )\Delta Y,\tag*{($3''$)}\\ \frac{d}{dt} \Delta Y=\,&\beta Y_0 \Delta X+(\beta X_0-\nu )\Delta Y,\tag*{($4''$)} \end{align*}取定存在之穩態 $SS_1(N_0,0)=(X_0,Y_0 )$ 代入可得 \begin{align*} \frac{d}{dt} \Delta X=&-\Psi'\Delta X+(\emptyset-\Psi'-\beta N_0 )\Delta Y,\tag*{($3^*$-1)}\\ \frac{d}{dt} \Delta Y=&(\beta N_0-\nu )\Delta Y.\tag*{($4^*$-1)} \end{align*}上述兩 ODEs 可解得解析解為 \begin{align*} \Delta X=&\Delta X_i e^{-\Psi't}+(\emptyset-\Psi'-\beta N_0 )\Delta Y_i \Big(e^{-(\nu -\beta N_0 )t}-e^{-\Psi't} \Big),\\ \Delta Y=&\Delta Y_i e^{(\beta N_0-\nu)t}. \end{align*}若 $N_0 \lt \frac \nu \beta $, 當達到平衡穩態時, 可得 $\lim\limits_{t\to \infty}\Delta X=\lim\limits_{t\to \infty}\Delta Y=0$ 疫情可有清零 (零確診) 之穩定平衡點。 反之, 若 $N_0\gt\frac \nu\beta $, 則不會在平衡點 $SS_1(N_0,0)$ 穩定存在 (亦即疫情失控而遠離 $SS_1$)。 針對平衡穩態 $$SS_2\Big(\frac\nu\beta,\frac{\Psi'(N_0-\frac\nu\beta)}{\Psi'+\nu'}\Big)$$ 取定線性方程組 ($3^*$), ($4^*$) 之 $\det(\underline A-\lambda \underline I)=0$ 可得兩個特徵值 $\lambda_1,\lambda_2$ 必須滿足如下方程式: \begin{align*} \lambda ^2+\Big(\frac{\beta\Psi'(N_0-\frac\nu\beta)}{\Psi'+\nu'}+\Psi'\Big)\lambda+\beta\Psi'\Big(N_0-\frac v\beta\Big)=0,\\ {\hbox{亦即 }} \lambda ^2-(\lambda_1+\lambda_2)\lambda+\lambda_1\lambda_2=(\lambda-\lambda_1)(\lambda-\lambda_2)=0. \end{align*}若 $N_0\gt\frac\nu\beta$ 條件下, 可得到兩個具有負實部之特徵值, 因此疫情可達到趨於穩態之平衡點 $$SS_2\Big(\frac\nu\beta,\frac{\Psi'(N_0-\frac\nu\beta)}{\Psi'+\nu'}\Big),$$ (亦即「與病毒共存」之平衡狀態)。 由此可以說明何以在未知瘟疫爆發時, 第一時間就必須要做到最有效之隔離檢疫, 避免個體間產生不必要之接觸而形成傳播, 才能致使區隔出來之每個小單位區域中 $N_0 \lt \frac\nu\beta $, 如此疫情才有逐漸式微, 不易形成流行的機會。 事實上, $\frac\nu\beta$ 於流行病學上更稱之為「閥值人口」(threshold population)數目。 新近新聞指出丹麥成為歐盟首個解除疫情限制措施國家, 不少國家亦會跟進。 此作法亦即主動採取「與病毒共存」之策略。 取定此策略就代表爾後之現有醫療量能, 已經有足夠把握來應付 $Y_0:Z_0$ 幾乎呈現固定比例 $\frac{\Psi'}{\nu'}$ 之確診比例情況。
(D) 反應機制(I$^*$)---
針對前文 其中 $\nu =\nu'+\emptyset \gt\emptyset$。 若假設總人口數定值: $X+Y+Z=N_0$, 可建立出流行病學疫情系統方程式如下: \begin{align*} \frac{dX}{dt}=\,&-\beta XY+\emptyset Y=f(X,Y),\qquad X(t=0)=X_i,\\ \frac{dY}{dt}=\,&\beta XY-\nu Y=g(X,Y),\hskip 20pt \qquad Y(t=0)=Y_i\gt0. \end{align*}顯然在任何 $X_0$ 值條件下, 只要 $Y_0=0$ 之解, 皆可以滿足上述兩式為平衡穩態之結果。 其線性化 ODEs 亦是式 \eqref{3} 及 \eqref{4}, 若式 \eqref{5} 中設定 $Y_0=0$ 為解, 則可得 \begin{align} \lambda ^2+(\nu -\beta X_0 )\lambda =0, \end{align}其中兩個特徵值分別為 0 與 $\beta X_0-\nu $, 其相對應之特徵矩陣則分別為 $\left[\begin{array}{c}1\\ 0\end{array}\right]$ 與 $\left[\begin{array}{c}\emptyset-\beta X_0\\ -\beta X_0+\nu\end{array}\right]$。 若起始條件 $\left[\begin{array}{c}\Delta X(t=0)\\ \Delta Y(t=0)\end{array}\right]=\left[\begin{array}{c} \Delta X_i\\ \Delta Y_i\end{array}\right]=\left[\begin{array}{c}X_i-X_0\\ Y_i-Y_0\end{array}\right]$, 其矩陣解則可表為 $$\left[\begin{array}{c}\Delta X\\ \Delta Y\end{array}\right]=C_1\left[\begin{array}{c}1\\ 0\end{array}\right]+C_2 e^{(\beta X_0-\nu )t}\left[\begin{array}{c}\emptyset-\beta X_0\\ -\beta X_0+\nu\end{array}\right],$$其中 $C_1=X_i-X_0-\frac{Y_i-Y_0}{\nu -\beta X_0}(\emptyset-\beta X_0 )=C_1 (X_i,Y_i) $, 而且$C_2=\frac{Y_i-Y_0}{\nu -\beta X_0}$。 此式子亦指出在疫情爆發當下, 第一時間若能超前部署, 確實有效執行隔離分區檢疫, 致使起始條件即能滿足 $N_i=X_i\simeq X_0 \lt \frac\nu \beta $, 其中 $N_i$ 及 $X_i$ 分別代表起始系統內所有人口數及未被感染人口數(亦即起始時間 「零號確診病人」尚未出現, 並無感染現象發生)。 如此疫情才能逐漸壓制下來, 漸次走向「清零」。 除此之外, 更可得到關係式如下: (i) $\lim\limits_{t\to \infty}\Delta X=X_i-X_0-\dfrac{Y_i-Y_0}{\nu -\beta X_0}(\emptyset-\beta X_0)$ 或是 $\lim\limits_{t\to \infty}X=X_0 +C_1 (X_i,Y_i )=X_i-\dfrac{Y_i-Y_0}{\nu -\beta X_0}(\emptyset-\beta X_0)$, (ii) $\lim\limits_{t\to \infty}\Delta Y=0\quad\hbox{或是}\ \lim\limits_{t\to \infty}Y=Y_0$. 以上結果至少強調出三個疫情發展重點: (a)瘟疫發生時, 追蹤感染傳播鏈之疫調工作極為重要, 因為任何時間之 $X(t),Y(t)$ 以及各變異株傳染相關參數 $\beta ,\emptyset,\nu$ 皆主宰著疫情是否可有效控制達到可掌可控的情況下。 (b) 最為重要的是, 若在防疫無缺口來掌握疫情條件下, 更要確認在任何時間, $X(t) \lt \frac\nu \beta$ 此必要條件務必皆能滿足。 (c) 第一時間起始疫情條件 $(X_i,Y_i )$ ($Y_i$ 表起始時間之確診人口數目) 主宰著當確診數目歸零時, 仍存留於感染系統中, 可能具有被感染風險之人口數目 (「吹哨人」效應) (即 $\lim\limits_{t\to \infty}X=X_i-\frac{Y_i-Y_0}{\nu -\beta X_0}(\emptyset-\beta X_0)$)。 致謝感謝教育部教學實踐研究計畫(PEE1100793)及科技部計畫(MOST 109-2221-E-197-016-MY3)之經費補助, 更感謝藉由此計畫執行過程中能得力於國立宜蘭大學教學發展中心的夥伴們持續協助及鼓勵, 以及校方教師社群計畫支助協助能聯結跨校與人文社會科學專業教授們(玄奘大學社會工作學系王敏菱、 中國文化大學心理輔導學系陳柏霖、 玄奘大學藝術設計學院王振邦、 國立台東大學生命科學系黃祥恩、 國立東華大學教育與潛能開發學系張德勝等位教授)有腦力激盪實質交流, 方有如此人文關懷研究的進行。 最後亦致上於心底最誠摯的感謝: 能在1991年間於台大公共衛生學院林瑞雄院長手下擔任研究工作, 才能有機會對流行病學有初步的入門體會。更影響到後續論文能以創新的病毒疫情感染流行病學模型來理解微生物生產製程參考文獻本文作者陳博彥任教於國立宜蘭大學化學工程與材料工程學系 |