作成: 令和3年8月24日, 9月3日編集，創価大学理工学部 畝見達夫
SimEpidemic 個体ベース感染シミュレータ 技術文書
デルタ株の爆発的な感染拡大が全国に広まっている。 発症者のほとんどはワクチン未接種であり， 40〜50歳代に加え，20〜30歳代の患者にも重症化する場合がでてきている。
ワクチン接種の進捗に応じて感染者数は減少することが期待されるが， 本国および英米の現状からすると，減少に転じるには日本ではワクチン接種率が不十分であると推測される。
シミュレーションのパラメータを デルタ株下での感染対策#5 と同様の設定とする。 ただし，感染による獲得免疫の効果についてモデルを改良し， 最新の東京都の陽性患者数の推移と一致するよう， 試行的なシミュレーションを通してシナリオを再調整した。 また，前回のシミュレーションでは昨年12月22日からの64回の試行結果から， 今年1月1日から8月14日までの東京の陽性患者数推移との差が小さい4つの試行を抽出し， それらの継続として各16試行ずつ，合計で64回の試行を実行した。 今回は8月23日までの試行から4つを抽出して用いた。
より厳しい行動制限を8月30日から実施したと仮定し， 制限により集会頻度および移動頻度をそれまでの 20, 40, 60, 80, 100% とした場合について， 陽性患者数，感染者数，重症者数，市中不顕性感染者数の推移を調べる。 シミュレーションの実行は昨年12月22日からの始めているが，ここでは5月16日以降について示す。 下に陽性患者数の推移を示す。細い折れ線はパラメータとして設定した集会頻度の推移， 破線は平均±標準偏差である。
下に感染者数の推移を示す。灰色の細い線は1回目のワクチン接種率の推移であり， 値は右側の目盛りに相当する。
下に重症者数の推移を示す。
6月下旬以降で実データとシミュレーション結果の乖離が見られる。 これは，この頃から増加したデルタ株の重症化速度の上昇を適切にモデル化できていないためと考えられる。
下に市中不顕性感染者数，つまり，隔離されていない無症状の感染者数の推移を示す。
平均値のピークの日付と値は以下の通りである。
| 集会移動 | 日付 | ピーク値 | 東京都換算 | |
|---|---|---|---|---|
| 陽性患者数 | 20% | 9月6日 | 0.0445% | 6,218人 |
| 40% | 9月7日 | 0.0457% | 6,375人 | |
| 60% | 9月9日 | 0.0464% | 6,481人 | |
| 80% | 9月12日 | 0.0496% | 6,921人 | |
| 100% | 9月16日 | 0.0542% | 7,570人 | |
| 感染者数 | 20% | 9月5日 | 0.809% | 112,884人 |
| 40% | 9月7日 | 0.828% | 115,606人 | |
| 60% | 9月9日 | 0.836% | 116,652人 | |
| 80% | 9月15日 | 0.893% | 124,633人 | |
| 100% | 9月18日 | 0.967% | 135,051人 | |
| 重症者数 | 20% | 9月24日 | 0.00213% | 297人 |
| 40% | 9月24日 | 0.00213% | 297人 | |
| 60% | 9月23日 | 0.00216% | 301人 | |
| 80% | 9月29日 | 0.00226% | 315人 | |
| 100% | 9月30日 | 0.00240% | 335人 | |
| 市中不顕性感染者数 | 20% | 9月3日 | 0.356% | 49,657人 |
| 40% | 9月4日 | 0.361% | 50,376人 | |
| 60% | 9月6日 | 0.359% | 50,185人 | |
| 80% | 9月9日 | 0.381% | 53,171人 | |
| 100% | 9月14日 | 0.419% | 58,427人 |
行動制限の強さに応じて 9月初旬〜9月末にピークとなり，その後はワクチンの効果により減少すると推測される。 行動制限が弱いと、ピーク時 (9月中旬〜末) には， 陽性患者数 0.069%、感染者数 1.23%、重症者数 0.0031%、市中不顕性感染者数 0.53% に達する可能性がある。 免疫効果の経時的減衰により、感染者数が高い状態が持続する可能性が示唆され， 3回目の接種の実施についても考える必要があろう。