УДК 519.218.28
DOI 10.17150/2500-2759.2021.31(2).221-229
МОДЕЛЬ КВАЗИСТАЦИОНАРНОЙ ЭПИДЕМИЧЕСКОЙ КИНЕТИКИ
А.В. Боровский1, А.Л. Галкин2
1 Байкальский государственный университет, г. Иркутск, Российская Федерация
2 Институт общей физики им. А.М. Прохорова РАН, г. Москва, Российская Федерация
Информация о статье
Дата поступления 12 апреля 2021 г.
Дата принятия к печати 21 июня 2021 г.
Дата онлайн-размещения 9 июля 2021 г.
Ключевые слова
Математическое моделирование эпидемий; теоретические модели эпидемий; эпидемическая кинетика; модель эпидемии с запаздыванием
Аннотация
Распространение инфекционных заболеваний негативно сказывается на социально-экономических показателях развития страны, ведет к ухудшению качества жизни граждан. В связи с этим представляется актуальным теоретическое исследование процесса пространственного распространения эпидемий. Целью данной работы является исследование модели эпидемической кинетики с запаздыванием по времени численными методами, в том числе асимптотических свойств эпидемического процесса в зависимости от источника заражения. В статье предложена новая теоретическая модель эпидемической кинетики, учитывающая скрытый инкубационный период заболевания в виде слагаемых с запаздыванием во времени. Модель учитывает четыре вида членов популяции — неинфицированных, не обладающих иммунитетом, активно инфицированных, выздоровевших и приобретших иммунитет, испытавших летальный исход. Разработанная модель эпидемической кинетики с запаздыванием учитывает возможность инфицирования членов популяции на интервале времени, равном инкубационному периоду заболевания, а также смертность только вне инкубационного периода. Предлагаемая модель эпидемической кинетики с запаздыванием позволяет системным образом изучать влияние различных эффектов на эпидемическую кривую, благодаря чему на практике возможно прогнозировать развитие эпидемической ситуации для введения необходимых противоэпидемических мер, что особенно актуально в условиях разразившейся мировой пандемии коронавируса.
MODEL OF QUASI-STATIONARY EPIDEMIC KINETICS
Andrei V. Borovsky1, Andrei L. Galkin2
1 Baikal State University, Irkutsk, the Russian Federation
2 Prokhorov General Physics Institute of the Russian Academy of Sciences, Moscow, the Russian Federation
Article info
Received April 12, 2021
Accepted June 21, 2021
Available online July 9, 2021
Keywords
Mathematical modeling of epidemics; theoretical models of epidemics; epidemic kinetics; lag epidemic model
Abstract
A new theoretical model of epidemic kinetics, which takes into account the latent incubation period of the disease in the form of time lagged terms, is viewed. The model takes into account four types of population members: the uninfected (non-immune), the actively infected, the recovered and acquired immunity, and the experienced a lethal outcome. The model considers the possibility of introducing anti-epidemic measures smoothly, as well as the presence of various types of infection of the uninfected contingent. Numerical calculations of the epidemic development show that the initial exponential growth of actively infected people after the introduction of quarantine measures is replaced by a decline in the epidemic curve within two — three weeks. Then, after three months, having a permanent source of infection, the epidemic enters a quasi-stationary mode of functioning. The quasi-stationary values statistics of actively infected individuals uniquely determines the size of the infection source. Calculations of the problem with a time-varying infection source describe the «second wave» of a separate intensity epidemic.
© Боровский А.В., Галкин А.Л., 2021
01 И 5<
а
л
V
п *
о
о
о
а
и ^
о ч
я ф
X X
о
п
о у
X
ч
ф
ч
2 О 2
Ы ^
Z
10 р
2 2
N N (О
Введение
Распространение ифекционных заболеваний оказывает большое влияние на жизнь людей. Они останавливают социально-экономическое развитие общества, отрицательно отражаются на демографических показателях, ухудшают качество жизни граждан. Все эти процессы нуждаются в изучении, причем в его основе должны лежать эпидемические модели. Их изучение приобретает особую актуальность в условиях разразившейся мировой пандемии коронавируса.
В медицинской литературе основоположниками теории эпидемий считают Даниила Бернулли, Уильяма Фарра, Петра Енько, Джона Браунли. Ссылки на их труды можно найти в Медицинской энциклопедии, Большой Советской Энциклопедии, а также в Википе-дии. Известны также труды их последователей [1—5]. Все работы делятся на два вида. Одни развивают теоретико-вероятностный подход для описания пространственного распространения эпидемии, другие используют балансный подход, применимый для населенного пункта. Эпидемический процесс во втором виде этих классических работ описывается с помощью системы обыкновенных дифференциальных уравнений эпидемической кинетики, решения которых характеризуют изменение численности в различных подгруппах популяции (восприимчивые, зараженные, иммунные, выздоровевшие и т.д.).
Математическая модель эпидемии гриппа Барояна — Рвачева [6] представляет собой уже систему нелинейных интегродифферен-циальных уравнений в частных производных с соответствующими граничными и начальными условиями. В интегральной части модели [7] исследователи пытаются учесть особенности вируса при развитии инфекционного процесса (мутации вируса), в том числе в период инкубации. Наряду с обычной характеристикой в виде средней частоты передачи возбудителя, модель отягощена такими понятиями, как функция развития периода инкубации и функция развития инфекционного периода. Модель оправдывает классификацию как эпидемическая динамика, поскольку описывает волнообразное распространение эпидемии между городами при учете пассажиропотоков.
Проблемы прогнозирования эпидемий, верификации теоретических моделей, описания локальных вспышек инфекции представлены в статьях [8-10]. Ход эпидемии коронавируса разбирался в интернет-статьях [11; 12]. Их авторы делают вывод, что эпидемия коронавируса не описывается известными
моделями эпидемической кинетики (SIR, SIS, SEIR, MSEIR). В работе [13] рассматривалось развитие эпидемии коронавируса в Австрии на основе дискретной модели.
В данной работе, следуя [14; 15], развивается новая модель эпидемической кинетики. Эта модель использует физическую аналогию с кинетикой заселенностей уровней атома [16]. Отличие рассматриваемой модели от известной модели SIR заключается в дополнительном учете запаздывающих слагаемых, а также наличии у известных слагаемых временной задержки. Модель эпидемической кинетики с запаздыванием учитывает возможность инфицирования членов популяции на интервале времени равном инкубационному периоду заболевания, а также смертность только вне инкубационного периода. В результате в уравнениях появляются переменные в момент времени t и t — т, где т — длительность скрытого инкубационного периода заболевания. Модель является открытой, т.е. может учитывать для населенного пункта внешний источник заражения, что еще больше приближает ее к реальной ситуации. Модель эпидкинетики позволяет рассчитать эпидемическую кривую — зависимость от времени числа активно инфицированных членов популяции.
Отметим, что предлагаемая модель эпидемической кинетики с запаздыванием в отличие от дискретной модели [13] позволяет системным образом изучать влияние различных эффектов на эпидемическую кривую. К их числу относятся: экспоненциальный рост активно инфицированных лиц на начальном этапе эпидемии, влияние противоэпидемических мер и их плавного введения на ограничение роста эпидемической кривой, влияние длительности инкубационного периода и времени течения болезни на эпидемическую кривую, учет наличия амбулаторных больных и группы лиц, бессимптомно переносящих заболевание, и др.
Целью данной работы является исследование модели эпидемической кинетики с запаздыванием по времени численными методами, в том числе асимптотических свойств эпидемического процесса в зависимости от источника заражения. Отметим также, что предлагаемая модель относится к изучению эпидемии в городе. Ее не следует применять в целом к стране, состоящей из совокупности городов и сельской местности. В этом случае ситуация усложняется наличием пассажиропотоков между населенными пунктами и зарубежными странами, что требует применения более сложных моделей.
Скорость заражения в населенном пункте
Анализ различных видов заражения, возможных в населенном пункте (город, поселок, вахтовое производство), проведенный авторами данной работы, позволяет записать следующее общее выражение для скорости изменения активно инфицированного контингента:
dN2 dt
= k^(t, t-T, t-(T + T + )) + k~2V2(t, t-T,
t-(t + t +))-k+N2 (t) - kcN2(t-t ) + A.
(1)
Здесь первых два слагаемых в правой части (1) описывают два основных вида заражения. Первый относится к заражению неин-фицированных лиц в небольших устойчивых коллективах. (Это заражение на производстве, в учреждениях, в организациях среднего, специального и высшего образования. Сюда же относится заражение в городском транспорте, в небольших магазинах, кафе и ресторанах и других предприятиях сферы обслуживания, спорта и культуры.) Второй вид заражения происходит в многоквартирных домах и осуществляется через общие инженерные системы (подъезд, лифт, мусоропровод, вентиляция, канализация). Функции ^ и описывают количество инфицированных членов популяции, участвующих в первом и втором виде заражения. Они различаются, так как в первом виде заражения участвуют граждане, которые выходят на улицы, а во втором случае добавляются амбулаторные больные. Коэффициенты k~ и к2- — частота заражения для первого и второго случаев, T — длительность инкубационного периода, t+ — время выздоровления.
Отрицательные слагаемые в (1) описывают скорость выздоровления и летального исхода. Слагаемое А отвечает постоянно действующему источнику заражения, например скорости заражения в медицинских учреждениях, а также в результате прибытия инфицированных лиц с внешним пассажиропотоком.
Авторы допускают, что частотные коэффициенты в выражении (1) должны быть подвергнуты медицинской верификации. Для этого нужно обладать необходимым объемом медицинской статистики. В данной статье для установления общих закономерностей ограничимся модельным подходом.
Численные исследования уравнений эпидемической кинетики с запаздыванием
Скорость заражения (1) содержит основные слагаемые. Если ограничиться их учетом
и не принимать во внимание сезонную миграцию населения, то уравнения эпидемической кинетики, учитывающие запаздывание, принимают вид [15]:
dN1 dt
(t) = -k - |N 2 (t)- N 2 (t - t )|- A,
N-(-)t) = k~ |N2(0-N2(t -т)|- к +N2(t)-dt
- kcN2(t - t) + A, dN
dN3(t) = k+N^t), dt
dN4 dt
(t) = kcN2(t - T).
(2)
Здесь Ы0 — численность популяции за вычетом доли, имеющей иммунитет к заболеванию. Ы1, И2, Ы3, Ы4 — численность здоровых, не имеющих иммунитета; инфицированных; выздоровевших и испытавших летальный исход членов популяции; все переменные берутся в момент времени t. Для заболевания COVID-19, вызванного корона-вирусом SARS-CoV-2 (в отличие, например, от гриппа), примем время выздоровления т + = 15 суток, тогда частота выздоровления есть обратная величина к+ = 1/15 сутки-1. Аналогично частота летального исхода кс = 1/50 сутки-1 . Длительность инкубационного периода т = 7 суток; частота заражения представлена выражением
k "(t) = p(t)A n(t) = p(t)A n0(t)
Z N (t) = N0.
Ni(t) N '
(3)
Используются следующие выражения для частоты контактов и вероятности заражения:
A n(t) =
An,n, 0 < t < t' An,n-(Ann - An,)t;—Г t' < t < t", A n,, t" < t
p(t) =
p , 0 < t < t'
' in'
P,n -(Pn - Pf) -t—71' < t < t"
(4)
Формулы (4) описывают плавное введение санитарных норм. Обозначения А п , А п, р, рг соответствуют начальному и конечному значениям частоты контактов и вероятности заражения.
п н
01 ш s<
а
л т
п *
о
о
о
а ^
о ч
я ф
X X
о
п
о у
X
ч
ф
ч
2 О 2
Ы ^
Z
ю р
2 2
N N «0
Постановку задачи (2) следует дополнить условием неотрицательности искомых функций и начальными условиями, например:
Ni(t) >0, . = 1, ...4, N(0) = ы0, N2(0) = 1, N3(0) = 0, N4(0) = 0. (5)
Для численного решения задачи эпид-кинетики (1)—(5) используется разностная схема, особенностью которой является запоминание предыстории к данному моменту t (например, с шагом 0,1 суток) уже полученных функций N|(s), s < t.
В приближении N1(t) « N0 (к-(t) в (3) не зависит от N0), второе уравнение в (2) становится независимым, остальные компоненты определяются через это решение
В приближении N0 « N0 на малых временах без учета запаздывания из второго уравнения системы (2)
dN2 ~dt
(t) = (k "(О)- k + )N 2(t) + Л
получаем аналитическое выражение
Л
N 2 (t) = N2(0)e
(k "(0) - k + )t
+ -
X
X e
24 ' (k (0)- k+) (k (0) - k+ )t -1"| (6)
Вспышка заболеваемости (6) при положительном инкременте к -(0) — к+ > 0 может гаситься только за счет запаздывающих слагаемых в (2), рис. 1. В расчетах принято Д п. = 30, Д п{ = 20 и р. = 0,01,
N
р{ = 0,005. Отметим, что волнистое поведение N2 справа от максимума обусловлено сложной структурой второго уравнения системы (2).
При более интенсивном развитии эпидемии (инкремент имеет более высокие положительные значения, Д п. = 34,5, Д п. = 23
' 1П ' ' f
и рп = 0,01, р( = 0,005) спад эпидемической кривой определяется спадом к вследствие принятия карантинных мер (рис. 2, 3).
Ход эпидемической кривой на рис. 2, 3 показывает крайне опасную ситуацию, возникающую из-за опоздания введения санитарных мер. При начале введения санитарных мер через 30 и 50 дней с момента начала эпидемии количество одновременно инфицированных достигает 1 400 и 25 000 членов популяции.
На рис. 4 представлена эпидемическая кривая, когда в населенном пункте действует постоянный источник заражения. Таким источником может быть инфицированный контингент, прибывающий в населенный пункт с внешним пассажиропотоком. Другим вариантом может быть заражение населения в медицинских учреждениях.
Сместим реализацию санитарных мер на 20-й день с начала эпидемии. Получим уменьшение инфицированного контингента до 3 500 чел. на 30-й день с последующим уменьшением. По статистике коронавирус-ной инфекции в РФ на коечный фонд в больницах попадает примерно пятая часть от числа инфицированных на данный момент. В рассмотренном варианте число госпитализированных составит 700 больных.
80
60
40 -
20
20
40
60
80
00
Рис. 1. Решение системы (2) без источника (А = 0) с ^(0) = Н0 = 6 • 105 (Иркутск), Г = 40, 7" = 70 в коэффициентах (4)
t
0
20
40
60
80
100
Рис. 3. Решение системы (2) без источника (Л = 0) с t' = 50, t" = 80 в коэффициентах (4)
и Дп. = 34,5, Дп = 23
in f
N
5 000 -
4 000
3 000
2 000
1 000 -
Рис. 4. Решение системы (2) с источником (Л = 10) с t' = 30, t" = 60 в коэффициентах (4)
и Дп. = 30, Дп = 20
in f
01 И 5<
а
л т
п *
о
о
о
а ^
о ч
я ф
X X
о
п
о у
X
ч
ф
ч 2 2 О 2
2
Ы ^
Z
ю
2
О 2
2
N N (О
0
t
Стационарный режим развития эпидемии
В стационарном случае:
dN2 dt
= k - |N2(t )-N 2 (t — t )|-k +H2(t )--kcN2(t-T) + Л = 0,
N(t) = N(t - T)
решение имеет вид
N =-
Л
(7)
(8)
2 (к+ + кс)
Стационарное значение (8) не зависит от Ы0, к-. Остальные компоненты определяются через (8):
N ~-At, N ~ -4+^, N 1 3 (к+ + кс) 4
(k + + kc) '
т.е. стационарным является только число инфицированных. Отметим, что соотношение (7) означает отсутствие изменений в количестве инфицированных в течение инкубационного периода. Поэтому решение (2) не выходит на стационарное решение (8), а только асимптотически к нему приближается. Источник A = const > 0 приводит к постоянному суточному приросту зараженных, который компенсируется убылью активно инфицированных в результате выздоровлений и летальных исходов. Так как k+ >> kc, то
N2 = Лт +.
(9)
Таким образом в стационарном режиме развития эпидемии количество активно инфицированных равно числу инфицируемых в сутки, умноженному на время выздоровле-
ния. Если А = 10, т + = 15, то Ы2 = 150 инфицированных.
Это обстоятельство может быть использовано для определения интенсивности источника постоянного заражения по наиболее достоверной части эпидемической статистики. Вертикальными прямыми на рис. 4, 5 отмечена область нарастания карантинных мер, горизонтальной прямой — целочисленное стационарное значение 115.
Сравнение рис. 2, 4, 5 показывает, что появление стационарного источника заражения в количестве 10 чел/сутки приведет к росту одновременно инфицированных с 1 500 до 5 400 членов популяции при начале введения санитарных мер на 30-й день эпидемии и до 3 500 инфицированных при начале введения санитарных мер на 20-й день эпидемии.
Рассмотрим задачу с переменным во времени источником
Л(0 = Л Л + ае
(10)
Два варианта параметров: Вариант 1: А = 10 чел/сутки; а = 1,5; t0 = 210 суток; Аt = 45 суток.
Вариант 2: А = 10 чел/сутки; а = 26; t0 = 210 суток; Аt = 60 суток.
Результаты расчета представлены на рис. 6 (сплошные графики с цифрами 1,2) в сравнении с действием постоянного источника А = 10 (точки). Пунктиром построены графики для двух вариантов (1,2, цифры штрихованные) по формуле
N2(t) =
A(t)
(k + + kc)
(11)
dN,
dt
(t) = k" {0,8|N2(t)- N2(t - г)| + 0,21N2(t)-
На больших в^менах расчеты с переменным источникоми с постоянным источником даюс одинакоыый +и,2птоти^<^с2сый предел
для N2(t). Такой же предел при малых
временах имеют пунктирные графики. Диапазон Ы >100 можно рассматривать как квазистационарный (изменения N2(t) за время инкубации мало, что не совсем выполняется для второго варианта). Временной ход N2(t) в квазистационарном диапазоне повторяет временной ход Л(Ы) с небольшим сдвигом. Это означает, что для решения обратной задачи определения источника Л(Ы) по статистике N2(t) в первом приближении можно использовать обращение соотношения (8):
л(ы) = ^(^а + кс).
Сравнение численного и аналитического стационарных значений дает дополнительную оценку погрешности методики — совпадение пяти знаков.
Отметим, что рассмотренные варианты расчета задачи с переменным во времени источником описывают вторую волну эпидемии различной интенсивности.
Дополнительное заражение
бессимптомными инфицированными
Введем в задачу модифицированную скорость заражения
N
- ^(Ы-(е ре+))} - к+ ^(е ре+ )|-
- аcN2(t - е) рЛ. (12)
В (12) учитывается заражение населения не только в инкубационный период заболевания, но также бессимптомными инфицированными в течение всего периода болезни. Скорость выздоровления взята на интервале заболевания. Результаты моделирования уравнений эпидемической кинетики с запаздыванием со скоростью (12) представлены на рис. 7.
Заключение
В работе представлено исследование новой эпидемической модели с запаздывающими членами. Системный анализ различных видов заражения приводит к скорости инфицирования (1). Выделение главных слагаемых в (1) позволяет редуцировать скорость заражения к видам (2) и (12). Численное решение системы нестационарных уравнений эпидемической кинетики представлено на рис. 1-7. Решения позволяют учесть плавное введение противоэпидемических мер, а именно: привыкание населения к ношению масок и введение частичных карантинных мер, связанное с уменьшением частоты контактов населения. Показано, что решения чувствительны к временной задержке при введении противо-
3 500 7 i
3 000 1
2 500 : г Н \
2 000 : \ Т "" \
1 500 /
1 000 Г | \ / \\
500 ])1 ................................
50
100
50
200
250
300
350
Рис. 6. Решение системы (2) со скоростями заражения (2) и (10), Г = 20, 7" = 50 в коэффициентах (4) и Д п.п = 30, = 20 с постоянным источником А = 10 (точки) и двумя вариантами переменного источника (10) (сплошные кривые 1,2)
Примечание: пунктиром построены графики для двух вариантов (1,2, цифры штрихованные)
по формуле (11).
01 И 5<
а
л т
п *
о
о
о
а ^
о ч
я ф
X X
о
п
о у
X
ф ^
п S
ч
ф
ч
2 2
о
м ^
2 Ы
^ Z
10 22
О 2
N ^
2
N (О
t
0
N2 new f = 20, t" = 50, A = 10
8 000 -
6 000 - /V
4 000 - ^ N2 /Гк Nv N4
2 000 /
20
40
60
80
100
Рис. 7. Решение системы (2) со скоростями заражения (2) и (12) с источником (Л = 10), t' = 20, Г' = 50 в коэффициентах (4) и Д пп = 30, = 20
эпидемических мер. Характерный вид эпидемической кривой для города с численностью населения 600 тыс. чел. показан на рис. 5-7. Прототипом является город Иркутск.
Примерно через 80 дней с момента начала введения противоэпидемических мер количество инфицированных граждан, пройдя через пик, стремится к стационарному значению (8). В дальнейшем эпидемический процесс развивается по квазистационарному сценарию, когда постоянный источник заражения в (1) приводит к постоянному числу инфицированных (8) в населенном пункте. Для подавления инфекции нужно вводить полный карантин, т.е. закрывать населенный пункт от внешнего пассажиропотока. Кроме того, необходимо вводить строгие санитарные меры по обслуживанию населения в медицинских учреждениях, препятствующие заражению граждан. Выход статистической кривой для активно инфицированных на плато1 не является признаком спада эпидемии, а свидетельствует о неподавленном квазистационарном источнике заражения. Если населенный пункт не закрывать, то эпидемия охватывает все регионы и повсюду развивается в квазистационарном режиме. Соответствующим математическим аппаратом для теоретического описания такой ситуации может являться уравнение
1 Коронавирус: статистика // Яндекс. URL: https://yandex.ru/covid19/stat?utm_source=main_
notif&geoId=11266.
Фоккера — Планка или статистическая теория с применением уравнений Колмогорова.
На нестационарном этапе развития эпидемии, который продолжается примерно 100 дней с момента начала эпидемии, существенное влияние на количественные характеристики числа переболевших оказывает бессимптомное инфицирование части населения (рис. 7). Например, на конец марта, по официальной статистике в г. Иркутске (население — 650 тыс. чел.), переболело 38 тыс. чел., что составляет 6 % населения. В то же время проводимая вакцинация показывает, что процент населения, обладающего иммунитетом, составляет 45 %. Вакцинировано на конец марта 4 % населения города. Это означает, что около 35 % населения переболели бессимптомно, либо обошлись самолечением.
Проведенные расчеты с переменным во времени квазистационарным источником заражения демонстрируют возможность описания второй волны эпидемии различной интенсивности в рамках рассматриваемой теории, а также подход к определению величины источника заражения в квазистационарном приближении по статистике активно инфицированного контингента. Более точное определение источника заражения А(^ может производиться путем решения обратной задачи А(^) = Н3^), где правая часть — данные статистики, а левая часть — решение системы (2).
СПИСОК ИСПОЛЬЗОВАННОМ ЛИТЕРАТУРЫ
1. Kermack W.O. A Contribution to the Mathematical Theory of Epidemics / W.O. Kermack, A.G. McKend-rick // Journal of the Royal Statistical Society. Series A: Statistics in Society. — 1957. — Vol. 120, iss. 1. — P. 48-60.
2. Bartlett M.S. Measles Periodicity and Community Size / M.S. Bartlett // Proceedings of The Royal. Society A: Mathematical, Physical and Engineering Sciences. — 1927. — Vol. 115, iss. 772. — P. 700-721.
t
0
3. Stochastic Two-Agent Epidemic Simulation Models for a Community of Families / L. Elveback, E. Ackerman, L. Gatewood, J. Fox // American Journal of Epidemiology. — 1971. — Vol. 93, № 4. — P. 267-280.
4. Бейли Н. Математика в биологии и медицине / Н. Бейли ; пер. с англ. Е.Г. Коваленко. — Москва : Мир, 1970. — 326 с.
5. Bartlett M.S. An Introduction to Stochastic Processes, with Special Reference to Methods and Applications / M.S. Bartlett. — Cambridge : At the University Press, 1978. — 388 p.
6. Бароян О.В. Моделирование и прогнозирование эпидемий гриппа для территории СССР / О.В. Бароян, Л.А. Рвачев, Ю.Г. Иванников. — Москва : Изд-во ин-та эпидемиологии и микробиологии им. Н.Ф. Гамалеи АМН СССР, 1977. — 546 с.
7. Боев Б.В. Прогнозно-аналитические модели эпидемий (оценка последствий техногенных аварий и природных катастроф) : лекция / Б.В. Боев. — URL: https://armscontrol.ru/course/lectures05a/bvb050324.pdf.
8. Современное состояние проблемы математического моделирования и прогнозирования эпидемического процесса / А.А. Лопатин, В.А. Сафронов, А.С. Раздорский, Е.В. Куклев // Проблемы особо опасных инфекций. — 2010. — № 3 (105). — С. 28-30.
9. Верификация математических моделей при чуме / А.А. Лопатин, Е.В. Куклев, В.А. Сафронов [и др.] // Проблемы особо опасных инфекций. — 2012. — № 3 (113). — С. 26-28.
10. Бачинский А.Г. Универсальная модель локальных эпидемий, вызываемых возбудителями особо опасных инфекций: моделирование новых инфекций / А.Г. Бачинский, Л.Ф. Низоленко // Проблемы особо опасных инфекций. — 2015. — № 1. — С. 51-53.
11. Тамм М.В. Коронавирусная инфекция в Москве: прогнозы и сценарии / М.В. Тамм // Фармакоэконо-мика. Современная фармакоэкономика и фармакоэпидемиология. — 2020. — Т. 13, № 1. — С. 43-51.
12. Иванов М.В. Математическое моделирование процесса пандемии / М.В. Иванов // Институт развития стратегических инициатив (ИРСИ). — 2020. — 30 апр. — URL: https://indsi.ru/2020/04/30.
13. Trigger S.A. Equation for Epidemic Spread with the Quarantine Measures: Application to COVID-19 / S.A. Trigger, E.B. Czerniawski. — DOI 10.1088/1402-4896/abb2e2 // Physica Scripta. — 2020. — Vol. 95, № 10. — P. 105001.
14. Боровский А.В. Модель эпидемии с запаздыванием / А.В. Боровский // System Analysis and Mathematical Modeling. — 2020. — Т. 2, № 4. — С. 53-63.
15. Боровский А.В. Математическая модель эпидемии со скрытым инкубационным периодом заражения / А.В. Боровский, А.Л. Галкин // Передовое развитие современной науки как драйвер роста экономики и социальной сферы : труды 2-й Всерос. науч.-практ. конф., Петрозаводск, 20 дек. 2020 г. — Петрозаводск, 2020. — С. 182-186.
16. Плазма многозарядных ионов. Элементарные процессы, кинетика и рентгеновские лазеры / А.В. Боровский, С.А. Запрягаев, О.И. Зацаринный, Н.Л. Манаков. — Санкт-Петербург : Химия : С.-Петерб. отд-ние, 1995. — 342 с.
Информация об авторах
Боровский Андрей Викторович — доктор физико-математических наук, профессор, кафедра математических методов и цифровых технологий, Байкальский государственный университет, г. Иркутск, Российская Федерация, e-mail: andrei-borovskii@mail.ru.
Галкин Андрей Леонидович — доктор физико-математических наук, Институт общей физики им. А.М. Прохорова Российской академии наук, г. Москва, Российская Федерация, e-mail: galkin@ kapella.gpi.ru.
Authors
Andrei V. Borovsky — D.Sc. in Physical and Mathematical Sciences, Professor, Department of Mathematical Methods and Digital Technologies, Baikal State University, Irkutsk, the Russian Federation, e-mail: an-drei-borovskii@mail.ru.
Andrei L. Galkin — D.Sc. in Physical and Mathematical Sciences, Prokhorov General Physics Institute of the Russian Academy of Sciences, Moscow, the Russian Federation, e-mail: galkin@kapella.gpi.ru.
Для цитирования
Боровский А.В. Модель квазистационарной эпидемической кинетики / А.В. Боровский, А.Л. Галкин. — DOI 10.17150/2500-2759.2021.31(2).221-229 // Известия Байкальского государственного университета. — 2021. — Т. 31, № 2. — С. 221-229.
For Citation
Borovsky A.V., Galkin A.L. Model of Quasi-Stationary Epidemic Kinetics. Izvestiya Baikal'skogo go-sudarstvennogo universiteta = Bulletin of Baikal State University, 2021, vol. 31, no. 2, pp. 221-229. DOI: 10.17150/2500-2759.2021.31(2).221-229. (In Russian).