НАН України: Прогноз розвитку епідемії COVID-19 в Україні на 13-26 жовтня 2021 року (інфографіка)
Робоча група з математичного моделювання проблем, пов’язаних з епідемією коронавірусу SARS-CoV-2 в Україні (базова установа – Інститут проблем математичних машин і систем НАН України) створена Розпорядженням Президії НАН України від 3 квітня 2020 р. №198.
Прогноз розвитку епідемії COVID-19 в Україні в період 13-26 жовтня 2021 р. (Прогноз РГ-53)
ЗМІСТ
Вступ
1. Загальнонаціональна епідемічна динаміка
1.1. Динаміка захворюваності, летальних випадків і тестування
1.2. Показники навантаження на лікарні
1.3. Спостережувана летальність, одужання та вікова структура інфікованих
2. Огляд демографічних даних
3. Аналіз груп регіонів
4. Регіональна епідемічна динаміка
5. Аналіз затримок оприлюднення даних
6. Прогноз розвитку епідемії в Україні з використанням статистичної моделі часових рядів Facebook Prophet
7. Порівняння з прогнозом від 28.09.2021 року
8. Прогноз розвитку епідемії в Україні з використанням компартментної моделі
Висновки
12.10.2021
Вступ
З початку квітня 2020 р. міжвідомча Робоча група (РГ) представників Національної академії наук України, Київського національного університету імені Тараса Шевченка, Вінницького національного технічного університету та Національної академії медичних наук України – з урахуванням світового досвіду математичного моделювання розвитку епідемії COVID-19, на основі статистичних даних про динаміку епідемії в Україні та країнах Європи – створювала і тестувала математичну модель SEIR–U. За результатами проведеного моделювання Робоча група підготувала документи «Прогноз розвитку епідемії коронавірусу SARS-CoV-2 в Україні».
Президія НАН України офіційно представила прогнози до державних органів. На ці прогнози надавав посилання Кабінет Міністрів України під час брифінгів.
У новому документі «Прогноз РГ-53»:
• здійснено аналіз первинних статистичних даних розвитку епідемії в Україні;
• порівняно статистичні дані з попереднім прогнозом «Прогноз РГ-52»;
• представлено прогноз розвитку епідемії на наступний період 13–26 жовтня 2021 року;
• проаналізовано динаміку поширення епідемії в регіонах України.
Робоча група використовує для аналізу такі джерела даних:
1. Дані щоденних звітів Міністерства охорони здоров’я (МОЗ) України. Дані про кількість нових виявлень, одужань, летальних випадків і підозр для кожної області України оприлюднюються на щоденних брифінгах МОЗ України, а також на сайті Ради національної безпеки і оборони (РНБО) України.
2. Первинні дані Центру громадського здоров’я Міністерства охорони здоров’я (ЦГЗ МОЗ) України. Опис структури даних та їх інтерпретацію наведено у документі «Прогноз РГ-9» від 26.06.2020 р.
3. Дані щодо обсягів ПЛР- та ІФА-тестування у кожній лабораторії України. Дані надає ЦГЗ МОЗ України, вони містять інформацію про загальну кількість проведених тестів, кількість позитивних тестів і ретестувань.
4. Дані Національної служби здоров’я України (НСЗУ) оновлюються щодня та містять дані про госпіталізації у кожній лікарні України.
5. Дані про загальну смертність з усіх причин і дані щодо мобільності населення з ресурсу «Економіка карантину».
6. Статистичні дані щодо країн світу: Worldometer, Financial Times, Ourworldindata.
1. Загальнонаціональна епідемічна динаміка
В останні два тижні захворюваність і далі швидко зростала за всіма ключовими показниками. Порівняно з попереднім двотижневим періодом, кількість нових випадків хвороби зросла у 1,9 раза, кількість позитивних результатів тестів ПЛР – у 2,1 раза, нових госпіталізацій – у 1,7 раза. При цьому темпи тижневого зростання в останні 1–3 дні демонструють певне зниження (особливо це стосується госпіталізацій, рис.4). Згладжені темпи тижневого зростання нових госпіталізацій знизилися до 25%, що є мінімальним значенням за останній місяць. Водночас, аналогічні очікування продовження зниження не виправдалися два тижні тому, коли було ще вираженіше зниження темпів зростання. Тоді темпи знову збільшились і стабілізувались: госпіталізації надалі демонстрували зростання з темпом близько 30% за тиждень, нові випадки хвороби – з темпом 40% за тиждень, а позитивні ПЛР-тести – з темпом 40-50% за тиждень.
За офіційною статистикою, на 12 жовтня 2021 року середня кількість осіб, що одужують, склала 5419 осіб на день. Середня за тиждень кількість нових інфікованих на 12.10.2021 р. становить 13175 на день для України, що майже вдвічі більше, ніж два тижні тому, і вп’ятеро більше, ніж місяць тому. Середня за тиждень кількість нових летальних випадків становить 264 на добу, що в 2,3 раза більше, ніж 2 тижні тому. Кількість активних інфікованих (які ще не одужали) також зросла більш ніж удвічі за два тижні та становить 205144. Нові госпіталізації зросли на 72% за два тижні – їхня кількість становить у середньому 3293 особи на день. Загалом темпи поширення лишилися майже незмінними за останні 2 тижні.
1.1. Динаміка захворюваності, летальних випадків і тестування
На рис.1 і 2 епідемічні показники наведено разом із даними за аналогічні дати 2020 року. Таке порівняння може бути доцільним для розуміння часових рамок перебігу нової хвилі.
Порівняно з минулорічними показниками у ті самі дати, позитивних тестів зараз удвічі більше, випадків хвороби – у 2,5 раза більше, а летальних випадків – втричі більше (рис.1, 4). Варто відзначити, що 7-денне середнє добового числа випадків хвороби склало 264, і це більше максимального значення за весь минулий рік. Крім того, за останні два тижні цей показник продемонстрував зростання у 2,3 раза.
Обсяги ПЛР-тестування (усереднені за 7 днів) вже перевищили максимальні торішні показники, але на 10% поступаються цьогорічним максимумам (рис.3).
1.2. Показники навантаження на лікарні
Кількість зайнятих ліжок хворими з підтвердженою коронавірусною інфекцією або підозрою на неї наблизилось до 30 тис., що більше за минулорічний максимум і становить близько 60% від цьогорічного максимуму (рис.5). При цьому кількість хворих, які потребують кисню, перевищує минулорічний максимум у 1,8 раза та сягнула 70% від цьогорічного максимуму, хоча ще два тижні тому їх було лише 42% від цьогорічного максимуму.
Водночас, ситуація у регіонах різна (див. дашборд РНБО, розділ «Ліжковий фонд під COVID-19»). Тоді як у Вінницькій, Івано-Франківській, Київській, Кіровоградській, Закарпатській, Полтавській областях та у Києві зайнятих ліжок менше 30% від максимумів попередніх двох хвиль, у Волинській, Донецькій, Запорізькій, Одеській, Харківській та Херсонській областях їхня кількість уже перевищила ці максимуми або впритул до них наблизилась. При цьому в деяких найскладніших регіонах уже розгорнуто безпрецедентно багато ліжкомісць.
На рис.6 разом з оперативними даними про нові госпіталізації показано динаміку у 2020 році. Кількість госпіталізацій за даними ЦГЗ уже вдвічі перевищує торішній показник для тих же дат і впритул наблизився до максимальних торішніх значень, поступаючись їм менш ніж на 5%.
1.3. Спостережувана летальність, одужання та вікова структура інфікованих
Спостережувана летальність для дат 4-5-тижневої давнини склала 2,5–3,0% проти значень 3,1–3,3% для дат 6-7-тижневої давнини (рис.7, таблиця 1).
Спостережувана летальність для найактуальніших дот ще зазнає ретроспективних змін, але найімовірніше, буде не менше 2,5% (рис.7). Отже, можливій захворюваності у 20 тис. щоденних виявлених хворих, очікуваній найближчим часом (див. прогнозну частину документу), відповідатимуть до 500 летальних випадків на день.
Середній вік виявлених хворих лишається низьким, а середній вік госпіталізованих хворих незначно зріс (рис.8).
Нинішня захворюваність серед дітей та молоді у вікових категоріях 0–9 і 10–19 років перевищила максимальні значення за весь час пандемії, хоча кількість госпіталізацій досі нижча від такого максимуму (рис.9, 10, таблиця 1).
Методологія. Пурпуровий графік на рис.7 показує очікувану спостережувану летальність з урахуванням статево-вікових коефіцієнтів летальності, що спостерігались у середньому в 2020 році. Значення спостережуваної летальності вище очікуваної на основі коефіцієнтів 2020 року в середньому означає, що спостережувана летальність у статево-вікових категоріях зросла, порівняно з 2020 роком. Пурпуровий графік коливається виключно за рахунок зміни статево-вікової структури виявлених хворих: якщо виявляють більше хворих похилого віку, то графік зростає, а якщо більше молоді – спадає.
2. Огляд демографічних даних
На початку поточного місяця надійшло оновлення даних Міністерства юстиції України щодо загальної смертності (ці дані доступні на сайті «Економіка карантину»), які дають відносно повну картину про число смертей до 25 вересня 2021 року (див. рис.11, де візуалізацію здійснено включно до 38-го тижня, який завершився 26 вересня 2021 року, та рис.12, де наведено динаміку оновлення даних).
Як і у попередні періоди високої захворюваності, надлишкова смертність над середнім рівнем 2015–2019 років виявилася вищою, ніж кількість летальних випадків при COVID-19. На основі візуалізації на рис.11 попередньо можна припустити, що перевищення становило приблизно два рази, хоча інші коливання загальної смертності є значними.
3. Аналіз груп регіонів
На рис.16 показано частку груп регіонів серед захворюваності, кольорову гаму яких задано на рис.17 (докладніше про групи регіонів див. у «Прогнозі РГ-38»).
На рис.13, 14 і 15 показано кількість нових хворих, позитивних тестів і госпіталізацій на 100 тис. населення.
В останні два тижні ситуація продовжила погіршуватись у всіх групах регіонів. Серед нових випадків хвороби найшвидше зростала частка південних регіонів, інцидентність у яких сягнула середнього рівня східних областей, а також групи регіонів «Центр та Північ», де захворюваність іще два тижні тому лишалась однією з найнижчих (рис.13, 16). Але тепер найменше випадків хвороби госпіталізують і виявляють у столиці та Київській області – за кількості позитивних тестів на рівні, близькому до середнього по країні, на тлі високого рівня тестування у столиці (рис.13–16).
4. Регіональна епідемічна динаміка
На картах рис.18–20 цифрами для регіонів позначено інцидентність, а також її робастна (очищена від шумових коливань) тижнева та двотижнева динаміка.
Попри незначне сповільнення зростання захворюваності на національному рівні, воно продовжилось у всіх, без винятків, регіонах України за кожним із таких оперативних показників, як нові госпіталізації (рис.18), нові позитивні тести (рис.19) і нові випадки хвороби (рис.20).
Найвищі темпи госпіталізації зі схожими показниками спостерігались у низці регіонів Сходу та Півдня країни, а також на Житомирщині, Сумщині й у Чернівецькій області, яка межує з Румунією, де захворюваність зараз вища, ніж в Україні (рис.18).
Крім того, у Чернівецькій області найбільша (після Луганської області) кількість летальних випадків у розрахунку на 100 тис. населення (рис. 22). Третю позицію за цим показником посідає Житомирська область, а далі розташовуються Донецька, Запорізька, Одеська, Сумська та Харківська області – з показником близько 6 летальних випадків на тиждень на 100 тис. населення.
Захворюваність у столиці лишається істотно нижчою, ніж загалом по країні, з огляду на дані про інтенсивність госпіталізацій та виявлення захворюваності, демонструючи також відносно слабкі показники зростання (особливо випадків хвороби) (рис.18, 20). Водночас, кількість позитивних тестів зростала несподівано швидко, що може свідчити про приховане зростання захворюваності (рис.19). Крім того, доволі швидко у столиці зростала кількість летальних випадків (рис.22).
5. Аналіз затримок оприлюднення даних
Рис.26 демонструє зміну кількості нових інфікованих, нових одужалих і нових летальних випадків за день згідно зі щоденними звітами МОЗ України для України загалом, що показані на момент публікації інформації (лівий графік), та ці ж дані, зведені до дат настання події (правий графік). Рис.26 (правий графік) демонструє нові випадки, показані на дати проведеного тестування, нові одужання, показані на дати одужання/виписки, та нові летальні випадки, показані на дату настання смерті. Дані на дату настання події є коректнішими і не мають нерегулярних стрибків, пов’язаних із затримками внесення даних до реєстру. Але дані на момент події змінюються ретроспективно, через це значення на кінцях інтервалу зазнаватимуть змін у майбутньому.\
Офіційні оприлюднені статистичні дані можна почасти пояснити на основі аналізу затримок оприлюднення інформації. Відомо, що щодня оголошувані нові випадки інфікування, нові смерті й нові одужання насправді не всі настають у попередній день. В окремих випадках можуть траплятися значні затримки між датою фактичного настання події (тестування, виписки, госпіталізації, смерті) і датою внесення цієї інформації до бази та її оприлюднення.
Хорошим індикатором реальної кількості нових виявлень на момент тестування може слугувати кількість ПЛР-тестів із позитивним результатом. Загальна кількість позитивних тестів стає відомою раніше, ніж інформацію про кожен випадок окремо буде внесено до реєстру. Графік на рис.28 показує співвідношення між кількістю позитивних тестів і кількістю нових виявлень на дату тестування.
Останнім часом кількість позитивних ПЛР-тестів непогано узгоджується із кількістю нових виявлень на день, порівняно з попередніми хвилями. Хоч, як і під час попередніх хвиль, нових випадків виявляють трохи менше за кількість позитивних тестів під час фази інтенсивного зростання.
Хоча загалом в Україні кількість позитивних ПЛР-тестів і нових виявлень добре узгоджується, але в Києві ситуація відрізняється. Рис.29 показує динаміку позитивних тестів і нових виявлень у Києві. Останнім часом кількість позитивних ПЛР-тестів значно перевищує кількість офіційних виявлень, що свідчить про великий відсоток тестів у приватних лабораторіях і/або значну частину тестувань не за місцем проживання. Позитивних тестів у Києві зараз роблять приблизно вдвічі більше за кількість офіційно виявлених випадків. Подібна ситуація спостерігалась і під час перебігу попередніх хвиль (рис.29).
Середні затримки стабільні та мінімальні за всю історію спостережень. Це свідчить про те, що ведення статистики суттєво поліпшилося, порівняно з минулим роком (рис.30). Структуру затримки в оприлюдненні нових випадків хвороби та нових летальних випадків показано на рис.31–34. У Києві ситуація із затримками результатів тестування поступово покращується, але залишається значно гіршою, ніж у середньому по країні (рис.33).
6. Прогноз розвитку епідемії в Україні з використанням статистичної моделі часових рядів Facebook Prophet
За допомогою методів статистичного аналізу та моделі часових рядів Facebook Prophet було досліджено динаміку щоденної кількості нових хворих для виявлення закономірностей поширення епідемії, для дослідження впливу свят і псевдосвят (аномальних дат на кшталт державних свят, теплих днів без опадів тощо), впливу тижневої та інших видів сезонної мінливості і виявлення їхнього характеру.
З огляду на те, що за щоденним приростом кількості нових хворих поточна хвиля сумарно по країні вже майже зрівнялася з весняною хвилею та стрімко наздоганяє найбільшу (весняну) хвилю, було вирішено повернутися до моделювання за всім рядом спостережень із початку першої великої хвилі – орієнтовно від 6 липня 2020 року. Такий підхід дозволяє більше врахувати вплив святкових днів, оскільки моделювання даних лише з кінця серпня 2021 року, як це робилося в цьому розділі у попередніх двох звітах, не враховувало такого впливу, оскільки єдиним державним святом був День Незалежності України, а значення тоді – суттєво меншими за нинішні. Водночас, найближчим часом очікуються 4 вихідні поспіль, зокрема святковий день 14 жовтня (День захисників і захисниць України), тому важливе вміння добре враховувати аномальний вплив таких днів. Для підвищення точності моделі додатково враховано несвяткові дні тижня, перенесені на інші дати для утворення кількох вихідних поспіль (наприклад, 15 жовтня 2021 року).
Отже, було побудовано модель, яка використовувалась у звітах протягом жовтня 2020 року – березня 2021 року, і за тим же рядом спостережень із 6 липня 2020 року. Модель ретельно враховує всі свята і псевдосвята (аномальні дати на кшталт державних свят і перенесених вихідних днів, теплих днів без опадів, усі дати зміни карантинних умов тощо). Ця модель за даними 06.07.2020 р. – 28.09.2021 р. дала прогноз на 13–26.10.2021 р. із сумарною відносною похибкою за останні 14 днів – 7,75% (рис.35, 36, таблиця 2).
Цікаво, що ця ж модель, але побудована без урахування останніх 7, а не 14 днів, тобто за даними 06.07.2020 р. – 05.10.2021 р. дала прогноз на 13–26.10.2021 р. із сумарною відносною похибкою за останні 7 днів – 1,51% і майже із такими ж значеннями, що доводить високу адекватність побудованої моделі станом на цей момент. Однак перша модель із 14-денним тестовим інтервалом, результати прогнозування за якою наведено на попередніх рисунках та у таблиці, видаються надійнішими.
Як і раніше, було побудовано моделі для 70 країн світу, зокрема для України, за всіма доступними даними за 2020–2021 роки і підготовлено прогноз на майбутні 7 днів. А потім за ними побудовано картограму приростів тренду країн Європи, де червоними колами показано збільшення кількості нових хворих, а синіми – зменшення; розмір кола пропорційний величині цього приросту.
Аналіз результатів моделювання приросту кількості нових підтверджених випадків захворювань показав таке:
- найближчим часом поточна хвиля в Україні, на жаль, наростатиме і далі;
- значення наростатимуть доволі стрімко, але нестабільно, саме тому попередні прогнози були неточними. Найбільша нестабільність очікується найближчим часом, за 4 вихідні поспіль (14-17.10.2021 р.);
- якщо щодо поточного тижня, через невизначеність показників під час свят, важко стверджувати, чи спостерігатиметься перевищення добового приросту кількості нових хворих у країні за весь попередній ряд спостережень, то на наступному тижні це точно відбуватиметься;
- аналіз прогнозів у країнах Європи показав, що у найвіддаленіших від України державах намітився спад кількості нових хворих, але у найближчих – цей показник, на жаль, зростає (особливо швидко – в Румунії), що може справляти додатковий негативний вплив на захворюваність у нашій країні;
- нестабільність динаміки зміни кількості нових хворих в Україні протягом цієї осені та невизначеність результатів спостережень під час 4-х вихідних поспіль (14–17 жовтня 2021 року), на жаль, не дає впевненості у підготовлених прогнозах, тому рекомендується використовувати їх обережно.
Обчислення за допомогою моделі Facebook Prophet і аналіз отриманих результатів виконали завідувач кафедри системного аналізу та інформаційних технологій (САІТ) Вінницького національного технічного університету (ВНТУ) доктор технічних наук, професор В.Б. Мокін і аспірант кафедри САІТ ВНТУ А.В. Лосенко.
7. Порівняння з прогнозом від 28.09.2021 року
Із рис.37 видно, що динаміка кількості нових виявлень дещо прискорилася, порівняно із середніми очікуваннями, і більше відповідає песимістичному сценарію розвитку епідемії.
8. Прогноз розвитку епідемії в Україні з використанням компартментної моделі
У таблиці 3 наведено середні прогнозні результати прогнозів для областей України на період до 19 та 26 жовтня 2021 року. Прогнози не враховують впливу зміни кількості тестувань, кількості підозр і зміни карантинних обмежень. При обчисленнях вважалося, що репродуктивне число та коефіцієнт летальності залишаються сталими впродовж прогностичного періоду. Репродуктивне число й коефіцієнт летальності обчислювалися з алгоритму калібрування математичної моделі, а для прогнозного сценарію використовувалися середні значення за останній тиждень.
Для врахування можливої зміни кількості контактів у часі було додатково розглянуто два сценарії для кожної області: з поступовим збільшенням рівня контактності на 25% та зменшенням на 25%. Зміна відбувалася поступово протягом тижня. У таблиці 2 наведено прогнозні середні рівні виявлення нових інфікованих на день і діапазон їхніх значень на 19 та 26 жовтня 2021 року. Наведено також оцінку середнього за тиждень значення репродуктивного числа для кожного регіону, що використовувалося для обчислення прогнозу.
Згідно із розглянутими сценаріями, регіонами з найбільшою очікуваною кількістю інфікованих є Дніпропетровська, Донецька, Житомирська, Запорізька, Львівська, Одеська, Сумська, Харківська, Херсонська області та м. Київ.
Для України загалом отримаємо такі прогнозні числові показники:
Репродуктивне число – 1.35 (середнє за останній тиждень, має тенденцію до зниження)
Середня кількість нових інфекцій за день на 19.10.21: [11399-17399] при середньому значенні 14897.
Середня кількість нових летальних випадків за день на 19.10.21: [309-367] при середньому значенні 338.
Середня кількість нових інфекцій за день на 26.10.21: [11927-31497] при середньому значенні 18640.
Середня кількість нових летальних випадків за день на 26.10.21: [376-468] при середньому значенні 422.
На рис.40 показано зміну в часі репродуктивного числа для України загалом, отриману в результаті калібрування математичної моделі на статистичних даних. Відтінками позначено оцінки репродуктивного числа для різних областей. На цей момент середньотижневе значення склало 1.35.
Таблиця 5 демонструє середні прогнозні значення нових виявлень компартментної моделі на період до 26 жовтня 2021 року.
Висновки
1. В останні два тижні захворюваність і далі швидко зростала за всіма ключовими показниками. Кількість нових випадків хвороби зросла вдвічі за 2 тижні та вп’ятеро за місяць. За два тижні кількість позитивних результатів тестів ПЛР зросла у 2,1 раза, а нових госпіталізацій – у 1,7 раза. Смертність виросла в 2,3 раза. За більшістю показників поточна хвиля зрівнялася та перевищила торішню осінню хвилю. Зокрема, смертність уже на 20% більша за пікові значення, а загальна кількість зайнятих ліжок у лікарнях – уже на 6%, і продовжує стрімко зростати. На цьому тижні очікуються рекордні значення кількості нових захворювань, але цьому, ймовірно, завадять довгі вихідні, що традиційно призводять до зменшення кількості офіційних виявлень. При збереженні поточних тенденцій можна очікувати досягнення показників весняної хвилі за 1–2 тижні.
2. Середня кількість нових виявлень в Україні збільшилася до 13175 осіб на день, середня кількість летальних випадків становить 264 смертей на день у середньому за останній тиждень. Згідно з прогнозними обчисленнями (в яких закладено поточну тенденцію до збільшення кількості нових виявлень), середня кількість нових випадків дорівнюватиме 14897 (SEIR_U), 14650 (Prophet) протягом тижня 13–19 жовтня 2021 року та 18640 (SEIR_U), 20922 (Prophet) – протягом тижня 20–26 жовтня 2021 року. Докладніше прогноз представлено у таблицях 2–5 і на рис.38–41.
3. У Києві ще спостерігається значна неузгодженість між кількістю позитивних ПЛР-тестів і нових виявлень, що свідчить про велику частку приватних лабораторій та тестування не за місцем проживання. Така специфіка даних у Києві ускладнює оцінку реальної епідемічної ситуації в столиці.
4. Спостережувана летальність дещо знизилась і перебувала у діапазоні від 2,5% до 3% для дат 4–5 тижнів тому.
5. Захворюваність серед дітей є найвищою за весь час пандемії, а інтенсивність їх госпіталізації наближається до абсолютних максимумів.
6. Найзагрозливішою, з погляду доступності обладнаного ліжкового фонду в лікарнях і показників смертності, є ситуація у Запорізькій, Одеській, Луганській, Чернівецькій, Житомирській та Харківській областях. У деяких областях кількість зайнятих ліжок уже перевищила максимуми обох попередніх хвиль.
Прес-служба НАН України