COVID-19感染リスク数値化の検討

Research on guidelines for quantifying of covid-19 infection risk

 日本では2020年1月16日に初めて感染者が確認されて以降、新規感染者数などの情報はメディアなどを通じて毎日何度も伝えられている。
しかしながら、自分が生活している地域やこれから移動する地域にどれくらいの感染リスクがあるのかについての指標は公開されているものは少なく日常生活をするにあたり不安である。そこで、感染拡大から1年間で得られたCOVID-19感染症の知見や公開されている情報から相対的な感染リスクを数値化する事を試みた。

COVID-19の特徴

  • 感染経路と隔離期間

 COVID-19の感染経路はアメリカ疾病予防管理センタ(CDC)によると
close contact > airborne transmission > contact with contaminated surfaces
の順に大きいとしている。また、隔離が必要な期間については、症状のある人は症状が出てから少なくとも10日間経過し、かつ解熱剤を使用せず発熱がなくなってから24Hr以上経過していること、かつその他症状が改善するまでとしている。なお、臭覚味覚の異常は長期におよぶ事があるが感染には影響しないとして除くとしている。また、無症状の陽性者は症状がなければ陽性確認後10日間の隔離が必要としている。
https://www.cdc.gov/coronavirus/2019-ncov/prevent-getting-sick/how-covid-spreads.html

また、日本国内においては感染経路について花王のホームページにわかりやすくまとめられている。
これによると、飛沫(5μm以上)と飛沫核(5μm以下)による飛沫感染経路と手指へのウィルスの付着による接触感染経路があり、飛沫の大きさや飛距離は温度、湿度、空気の流れにより変化するとされている。
https://www.kao.com/jp/hygiene-science/expert/inactivation-methods-new-coronavirus/transmission-route/

  • ウィルス(SARS- CoV-2)の生存性

 感染経路において、ウィルス自体の生存性は感染に大きな影響を与える。例えば手指についたウィルスが感染力を維持しているかに影響するためである。
日本防菌防黴学会のホームページに掲載されている麻布大学 野田 衛氏の「新型コロナウイルスの基礎知識、集団予防および生存性・不活化」の報告にまとめられている。それによると、新型コロナウィルスの生存性に及ぼす温度の影響が示されており、温度が低い程ウィルスの生存性が長くなる事がわかっている。なお、ウィルスの生存性と湿度の関係については他のコロナウィルスでの報告から湿度が高いほど生存性が短くなる傾向があるがさらなる検証が必要とのことである。
https://www.saaaj.jp/covid/pdf/covid02.pdf

  • 飛沫の拡散性

 感染経路において、飛沫の飛ぶ距離、量などについても同様に重要な因子である。内閣官房によるCOVID-19 AI・シミュレーションプロジェクトの成果報告において、理化学研究所/神戸大学 坪倉 誠氏の「室内環境におけるウイルス飛沫感染の予測とその対策 #3」の報告にまとめられている。それによると湿度が低いほど前方の人に到達する飛沫の量が大きくなる事が報告されている。
http://ww.covid19-ai.jp/ja-jp/organization/riken/articles/article003

感染リスクのモデル化検討

 これまでの知見から感染リスクの数値化に向けて感染モデルの仮説を作成した。ここでモデルとしたのは自助努力ではどうにもならない外的因子とした。
リスクモデルは 基本再生産数×感染者数×温度×湿度×人の移動×抗体非保有者数となる。
ここで基本再生産数と抗体非保有者数は一定であると仮定し除外し
リスクモデル=感染者数×温度×湿度×人の移動度とした。

感染リスクモデル図

感染リスク因子の数値化

  • 基本再生産数(R)

 国立感染症研究所によるとR0とは, ある感染症に対して全く免疫を持たない集団の中で, 1人の感染者が平均して何名の二次感染者を発生させるかを推定した値であるとしている。そのため感染リスクを数値化するには重要な因子である。しかしながら、コロナ変異種の影響など現状不明確なものが多いため、このモデルでは単一因子として扱い因子から除外する。今後、コロナ変異種の基本再生産数などが従来比明確となった時点で補正を行う。

  • 感染者数データ(I)

国内における新規感染者数はNHKにより都道府県別まとめられたものが日々更新されている。本モデルではこれを参照する事とする。また、隔離が必要な期間を参考に過去10日間の平均値をウィルス保有者数として数値化する事とした。なお、都道府県別の数値は人口10万人当たりの感染者数に変換し、東京都における2020年6月1日から10月31日までの人口10万人当たりの新規感染者数の平均値が1となるように係数補正を行った。これは日本では2020年前半は検査数自体の変化が大きく比較する事が難しいためである。このようにして求めた指数を東京にあてはめると最大10程度の値を示した。
https://www3.nhk.or.jp/news/special/coronavirus/data/?utm_int=news_contents_news-closeup_002

  • 気温データ(T)

温度とウィルス生存性の関係は下図のように温度とウィルス生存度比をまとめた。

温度とウィルス生存比の関係

なお、東京都の2020年の平均気温16.6℃でウィルス生存度比が1となるよう更に係数補正を行っている。都道府県別の気象データは気象庁のホームページで公開されているデータを用いた。
https://www.data.jma.go.jp/obd/stats/etrn/index.php

  • 湿度データ(H)

湿度と飛沫拡散比の関係は下図のようまとめた。

湿度と飛沫拡散比の関係

東京都の2020年の平均湿度71.5%RH時の飛沫拡散比が1となるように係数補正をおこなった。なお温度データと同様に都道府県別の気象データは気象庁のホームページで公開されているデータを用いた。下図をみてわかるように湿度に対する飛沫拡散比の変化は最小湿度が概ね30%RHから100%RHで推移する事を加味すると0から4の間で推移する。これは湿度による変化が温度と比べて影響度が大きい事を示す。

  • 人の動き(M)

人の動きはスマートフォンなどの位置情報を利用した解析が公開されている。本解析ではGoogle社の公開しているコミュニティモビリティレポートを参照した。このレポートは日本では都道府県単位でのカテゴリー毎の日毎のデータが公開されている。本解析では小売店と娯楽施設、乗換駅、職場の数値を平均し基準日が1となるように補正した。なお、基準日は2020 年 1 月 3 日〜2 月 6 日の 5 週間における該当曜日の中央値となっている。なお、都道府県間のデータ補正は行わない事とした。
なお、上記により求めた東京都の数値は0.23から1.01の範囲であった。
https://www.google.com/covid19/mobility/

都道府県別感染リスク数値化

都道府県別の感染危険度は各因子を乗じた数字とした。各因子を単純に掛け算する事は少し雑破な印象も受けるが今後実態と合わせて各因子の寄与度の調整を検討する。
感染危険度=R×I×T×H×M×100
検証の意味合いで北海道、仙台、東京、大阪、福岡の感染危険度の推移を下図に示す。二度目の緊急事態宣言時(1月)の東京都の感染危険度が突出した数値となっている。また、地域間の危険度の比較はほぼ実感と合っているように思われる。

主要都市エビデンスグラフ

まとめ

 以上の経緯からCOVID-19感染危険度をマップ化し公開することとした。自分生活している地域の状況および旅行などで移動する際の参考としていただきたい。今後実感と合わない点があれば適宜改訂するとともに改訂履歴を表記することとする。