Имитационное моделирование распространения эпидемий на основе агентного подхода

15

Размещено на http://www.allbest.ru/

01.00.00 Физико-математические науки

Имитационное моделирование распространения эпидемий на основе агентного подхода

Утакаева Ирина Хайрлыевна

к.ф.-м.н., доцент

SPIN-код: 9185-4692, AuthorID: 589215

Сегодня инфекционные болезни остаются одной из ведущих причин преждевременной смерти людей на Земле. Агентное моделирование может сыграть важную роль в прогнозировании распространения болезни и в оценке мер по локализации. Целью работы является построение имитационной мультиагентной модели распространения эпидемии для формирования мер по эффективному снижению уровня заболеваемости. Использование имитационного мультиагентного подхода в моделировании эпидемий обусловлено тем, что подход позволяет рассматривать много факторов влияющих на эпидемический процесс, дает возможность проводить численные эксперименты.

Процессы пространственного распространения и временного изменения этих двух групп эпидемий автор называет инфекционной динамикой. Обычно трудно реализуемые пространственные составляющие динамики в предлагаемых моделях берёт на себя топология предфрактального графа, которая наращивается объёмными графами - затравками, а динамика наращения предфрактального графа, называемая его распознаванием, отвечает за временную составляющую процесса. Под агентом понимается элементарный участник исследования. Агент активен, находится в некотором состоянии, которое может меняться при влиянии факторов. К свойствам агента отнесятся характеристики, формирующие уровень иммунитета: рост, вес, пол, доход, семейное положение, образование, география

Today, infectious diseases remain a leading cause of premature deaths in the world. Agent-based modeling can play an important role in predicting the spread of disease and to assess the containment measures. The aim is to construct a multi-agent simulation model for the formation of epidemic measures to reduce effectively their incidence. Using the multi-agent simulation approach to modeling of epidemics due to the fact that the approach allows us to consider a number of factors influencing the epidemic process, makes it possible to carry out numerical experiments. The processes of the spatial distribution and temporal variation of these two groups of epidemics of infectious the author calls dynamics. Usually hard-implemented spatial components of the dynamics in the proposed model can be substituted by predfractal topology of the graph, which is built up by voluminous graphs - primers, and the dynamics of compounding prefractal graph, called its recognition, is responsible for the timing of the process component. Under the term of agent, we consider an elementary study participant. An agent is active; it is in a state that may change under the influence of factors. The properties of the agent are attributed characteristics that form the level of immunity: height, weight, gender, income, marital status, education, geography

Ключевые слова: МОДЕЛЬ, АГЕНТ, ЭПИДЕМИЯ, РАСПРОСТРАНЕНИЕ, ГРАФ

Keywords: MODEL, AGENT, EPIDEMIC, SPREAD, GRAPH

Doi: 10.21515/1990-4665-121-085

По данным Всемирной организации здравоохранения ежегодно на земном шаре переносят инфекционные заболевания свыше 1 млрд. человек. Целью предупреждающих мероприятий является воздействие на источник инфекции, чтобы уменьшить обсеменение внешней среды, локализовать распространение микробов, повысить устойчивость населения к заболеваниям. Основными возбудителями инфекционных болезней являются вирусы, бактерии и простейшие. Заболевший становится источником возбудителей болезни. Он может заразить окружающих при контакте или путем обсеменения возбудителями различных объектов внешней среды. Особенно опасны больные, которые своевременно не обращаются за медицинской помощью. Выявить инфекционного больного - задача медицинских работников.

Механизм заражения имеет эпидемиологическое значение и положен в основу классификации инфекционных болезней: кишечные, инфекции дыхательных путей, кровяные инфекции, инфекции наружных покровов, инфекции с различным механизмом передачи. Различают воздушный, водный, пищевой и контактно-бытовой пути передачи инфекции. Есть немало моделей, описывающих процесс распространения инфекционных заболеваний, но они не учитывают демографическую структуру мегаполиса и систему контактов между жителями. Имитация эпидемиологической ситуации позволяет проверить эффективность. Зная пути передачи инфекции, власти могут воздействовать на сети с помощью закрытия заведений, введение карантина, изоляция и лечение.

Рассмотрим ситуацию, когда в городе началось распространение эпидемии и власти должны принять меры по предотвращению. Проводить массовую вакцинацию, ввести карантин, назначить антибиотики? Если отправить в город запасы противовирусных препаратов, то страны-доноры не смогут защитить население от инфекции. От выбора властей зависит не только жизнь людей, но экономическое и социальное благополучие страны. Чтобы выбрать оптимальную стратегию управления, нужно смоделировать сценарии развития событий. Имитация эпидемиологической ситуации позволяет проверить на модели эффективность мероприятий в борьбе со вспышками инфекционных заболеваний. Моделирование путей перемещения жителей города даёт динамическую картину социальной сети - аналогичную цепь контактов использует возбудитель инфекции, распространяясь в популяции. Техника социометрии впервые предложена Дж. Морено. Термин «социальная сеть» введен в 1954г. социологом Дж. Барнсом. В настоящее время ощущается дефицит систематического изложения методов и алгоритмов сетевого анализа, пригодных для современных прикладных исследований. Агентное моделирование может сыграть важную роль в прогнозировании распространения болезни, а также в оценке мер по локализации. Целью работы является построение имитационной мультиагентной модели распространения эпидемии для формирования мер по эффективному снижению уровня заболеваемости. Использование имитационного мультиагентного подхода в моделировании эпидемий обусловлено тем, что подход позволяет рассматривать много факторов влияющих на эпидемический процесс, дает возможность проводить численные эксперименты. Под социальной сетью на качественном уровне понимается структура, состоящая из множества агентов (субъектов индивидуальных или коллективных, например индивидов, семей, групп и организаций) и определённого на нём множества отношений (совокупности связей между агентами, например, знакомство, дружба, сотрудничество). Формально социальная сеть представляет граф , в котором - множество вершин (агентов) и - множество рёбер, соответствующих взаимодействию агентов. Итак, в эпидемиологии, под агентом будем понимать элементарного участника исследования. Агент активен, находится в некотором состоянии, которое может меняться при влиянии факторов. К свойствам агента можно отнести следующие характеристики, которые формируют в совокупности уровень иммунитета: рост h, вес w, пол s, доход d, семейное положение q, образование е, география g.

В модели целый искусственный город агентов, обладающих одинаковым набором свойств, но разных по значениям, т.е. некоторая БД агентов. Эти данные необходимы не только для моделирования на существующих эпидемиях, но мы должны обладать такого рода БД и для внезапных появлений новых вирусов. Обозначим через вершины отдельно взятых членов семьи, состоящей из человек, которые проживают вместе. Соединим вершины () и () ребром в случае, когда между членами семьи и имеется тесный бытовой контакт, достаточный для заражения исследуемым инфекционным заболеванием. В случае заражения одного из членов семьи, от него инфекционным заболеванием, передающимся воздушно-капельным путём, могут заразиться остальные члены семьи. Модель контактов между членами семьи, состоящей из 4 человек, будет представлять граф из 5 человек 5-вершинный полный граф :

15

Размещено на http://www.allbest.ru/

15

Размещено на http://www.allbest.ru/

Рис. 1. Модели контактов и

В графе-модели увеличится число вершин и рёбер, но останется удивительно стабильной его структура, состоящая из разномасштабных самоподобных частей с разным числом вершин в этих частях. Осталось назвать части затравками, а графовую модель - предфрактальным графом.

Таким образом, на этапе модель, описывающая контакты в семье, будет представлять собой -вершинный полный граф . На этапе модель, описывающая контакты семей, проживающих на одной лестничной площадке, представляется графом , который можно получить из графа , применяя операцию ЗВЗ к каждой её вершине.

Причем операцию ЗВЗ можно производить случайным образом графами или . Продолжая процесс, можно описать контакты многоэтажного дома, жилого квартала, города и т.д.

эпидемия имитационный мультиагентный

Рис. 2. Модель бытовых контактов людей, живущих на одной лестничной площадке в квартирах под номерами №1-№4.

Рис. 3 - Модель более сложных связей жителей одного подъезда

Повторяя этот процесс при , структура контактов будет представлять собой фрактальный граф , порождённый двумя полными затравками и . Если семьи состоят из другого числа людей, то затравки будут других типов. Всякий предфрактальный граф , из траектории построения фрактального графа является структурой контактов на -м этапе. Процесс построения фрактального графа можно обобщить на случай, когда операция ЗВЗ производится множеством затравок , . Полученная модель будет более адекватной, если дополнительно учесть ряд факторов:

1. Важной особенностью является то, что контакты человека с членами семьи оказываются более частыми, тесными и продолжительными, чем с соседом, а контакты с жителями соседних подъездов ещё слабее и т.д. Эти характеристики с течением времени меняются пропорционально, поэтому логично всем рёбрам предфрактального графа , в зависимости от их ранга определить вес по правилу , где - ранг ребра, - номер ребра -го ранга; - коэффициент пропорциональности, влияющий на изменение веса ребра.

2. Из теории эпидемий известно, что каждый человек обладает определённым уровнем иммунитета от инфекционного заболевания: часть людей делает прививки от инфекционных заболеваний и становятся неподверженными инфицированию и обладают максимальным уровнем иммунитета; люди с ослабленной иммунитетом более подвержены инфицированию. Чтобы учесть этот факт, присвоим вершинам графа , представляющего модель бытовых контактов, веса , . Под весом вершины будем понимать коэффициент, пропорциональный степени защищённости от инфекционного заболевания. Если человек входит в контакт с инфицированным, в математической интерпретации он заражается заболеванием с некоторой вероятностью =, . Возможны случаи:

· - человек не подвержен заболеванию, ему сделана прививка или уже переболел определённым заболеванием и больше не подвержен;

· - человек обязательно будет инфицирован при осуществлении бытовых контактов, достаточных для инфицирования;

· - человек подвержен заболеванию и при осуществлении бытовых контактов, возможно заражение.

Для социальных сетей ключевым показателем в процессе инфицирования является «эпидемический порог» (порог перколяции)--критическая вероятность заражения соседа, при превышении которой «инфекция» распространяется по всей сети. Если социальную сеть представить случайным предфрактальным графом, то инфекция с вероятностью заражения, выше эпидемического порога, экспоненциально размножается; вирус с вероятностью заражения, ниже порога перколяции, экспоненциально «вымирает». Предположим, что если человек, обладающий вероятностью заражения исследуемым инфекционным заболеванием имеет некоторый бытовой контакт с инфицированным, то становится инфицированным, если же , контактируемый остаётся здоровым. Для простоты исследования с учётом общей случайности заражения веса , вершинам графа приписываются произвольным образом генератором случайных чисел.

В теории перколяции рассматриваются процессы протекания слева-направо либо снизу-вверх, определим понятие «протекания» («просачиваемости») на предфрактальном графе со взвешенными вершинами, которые необходимо принять за исток и сток. Будем считать, что существует протекание между различными фиксированными затравками одного ранга и , представляющими исток и сток, если существует маршрут из вершин , соединяющий две вершины, где и , , .

Алгоритм поиска порога перколяции предфрактального графа;

1 ;

2 Осуществляется поиск всевозможных маршрутов, соединяющих исток и сток при ; если существует такой маршрут, то возвращаемся к началу шага 2 при , в противном случае и переходим к шагу 3;

3 . Если , полагаем и алгоритм прекращает свою работу, в противном случае возвращаемся к шагу 2.

Полученное значение обозначим через - порог перколяции.

Чтобы выявить связи, которые необходимо блокировать для исключения дальнейшего развития эпидемии, необходимо задать ранг рассматриваемой затравки и указать процент заражённых узлов, при котором возможно применение мер по карантину. Как только в затравке указанного ранга достигается процент инфицированных, производится блокирование всех рёбер, инцидентных инфицированным вершинам.

Предположим, что задана модель в виде предфрактального графа, взвешенного по всем вершинам и рёбрам. Рассмотрим структуру распространения инфекционного заболевания, если задан первичный источник инфекции. С математической точки зрения на этапе структура распространения инфекции представляет собой -вершинную звезду . В терминологии теории графов вершина, обозначающая первичный источник инфекции, является корнем.

15

Размещено на http://www.allbest.ru/

Рис.4. Структура распространения инфекции в виде дерева

Предположим, что на этапах распространения инфекции каждый «новый» источник в состоянии породить ещё источников инфекции. Процесс распространения инфекции в смысле замещения одного источника инфекции источниками, соответствует процессу порождения фрактального графа , порождённого множеством затравок , каждая из которых представляет собой -вершинную звезду степени . Ставя в соответствие математическое описание словесному, дальнейший переход к следующему уровню распространения будем производить по следующим правилам:

1 В построении графа участвует множество затравок , , что соответствует произвольному числу заражаемых индивидов уровня ;

2 Если вершина не является висячей, то она не замещается;

3 Замещаемая затравкой вершина выбирается из подмножества висячих вершин, а само е ребро становится инцидентным центру звезды;

4 Если какая-либо висячая вершина оказалась не замещённой затравкой, то она называется «замороженной» и по отношению к ней операция ЗВЗ не применяется ни на каком из последующих этапов.

Таким образом, из графа , применяя к висячим вершинам операцию ЗВЗ , получим граф , который является структурой распространения инфекции на следующем этапе . Всякий предфрактальный граф , из траектории фрактального графа оказывается структурой распространения инфекции на -м этапе. На рис. 5 предфрактальный граф , порождённый множеством различных затравок-звёзд. Произведя L переходов, получим предфрактальное корневоедерево , являющееся деревом распространения инфекционного заболевания и обозначим .

Рис. 5 - Конструкция в виде корневого дерева

Разработан программный продукт в среде FreePascal для имитационного моделирования с применением агентного подхода, который позволяет моделировать эпидемиологическую ситуацию, рассчитать эпидемический порог, порог перколяции, выбрать меры по карантину, чтобы разрушить сеть распространения.

СПИСОК ЛИТЕРАТУРЫ

Жиганова Л.П. Биотерроризм и агротерроризм-реальная угроза безопасности общества.

Бароян О.В., Рвачев Л.А. Математика и эпидемиология. - Москва, «Знание», 1977г.

Бейли Н. Математика в биологии и медицине. - Москва, «МИР», 1970 .

Боев Б.В. Современные этапы математического моделирования процессов развития и распространения инфекционных заболеваний // Эпидемиологическая кибернетика: модели, информация, эксперименты. Москва, 1991.

Воробьев А.А., Боев Б.В., Бондаренко В.М., Гинцбург А. Л. Проблема биотерроризма в современных условиях// ЖМЭИ.-2002.-№3.

Горелик А.Л., Скрипкин В.А.Методы распознавания. Москва,2004.

Емеличев В.А. и др. Лекции по теории графов. Москва, 1990.

Кочкаров А.М. Распознавание фрактальных графов: Алгоритмический подход. Нижний Архыз, 1998г.

Мандельброт Б. Фракталы, случай и финансы. Москва,2004.

Мандельброт Б. Фрактальная геометрия природы. Москва, 2002.

Онищенко Г.Г., Сандахчиев Л.С., Нетесов С.В., Щелкунов С.В. Биотерроризм как национальная и глобальная угроза //ЖМЭИ. -2000 №6.

Сапрыкин В. Биотерроризм: страшная реальность наших дней//Южный регион. Компетентно.

Сорокин Ю.Д. Биологический терроризм - мифы и реальность// Проблемы местного самоуправления

Тарасевич Ю.Ю. Перколяция: теория, приложения, алгоритмы. Москва,2002.

Утакаева И.Х., Кочкаров Р.А. К вопросу об алгоритмах распознавания предфрактального графа . НТВ СПбГПУ. Санкт-Петербург, 2010.

Шредер М. Фракталы, хаос, степенные законы. Москва,2005.

http://tehnic.by

http://armscontrol.ru/course

https://publications.hse.ru/articles/97267128.

References

1. Zhiganova L.P. Bioterrorizm i agroterrorizm-real'naja ugroza bezopasnosti obshhestva.

2. Barojan O.V., Rvachev L.A. Matematika i jepidemiologija. - Moskva, «Znanie», 1977g.

3. Bejli N. Matematika v biologii i medicine. - Moskva, «MIR», 1970 .

4. Boev B.V. Sovremennye jetapy matematicheskogo modelirovanija processov razvitija i rasprostranenija infekcionnyh zabolevanij // Jepidemiologicheskaja kibernetika: modeli, informacija, jeksperimenty. Moskva, 1991.

5. Vorob'ev A.A., Boev B.V., Bondarenko V.M., Gincburg A. L. Problema bioterrorizma v sovremennyh uslovijah// ZhMJeI.-2002.-№3.

6. Gorelik A.L., Skripkin V.A.Metody raspoznavanija. Moskva,2004.

7. Emelichev V.A. i dr. Lekcii po teorii grafov. Moskva, 1990.

8. Kochkarov A.M. Raspoznavanie fraktal'nyh grafov: Algoritmicheskij podhod. Nizhnij Arhyz, 1998g.

9. Mandel'brot B. Fraktaly, sluchaj i finansy. Moskva,2004.

10. Mandel'brot B. Fraktal'naja geometrija prirody. Moskva, 2002.

11. Onishhenko G.G., Sandahchiev L.S., Netesov S.V., Shhelkunov S.V. Bioterrorizm kak nacional'naja i global'naja ugroza //ZhMJeI. -2000 №6.

12. Saprykin V. Bioterrorizm: strashnaja real'nost' nashih dnej//Juzhnyj region. Kompetentno.

13. Sorokin Ju.D. Biologicheskij terrorizm - mify i real'nost'// Problemy mestnogo samoupravlenija

14. Tarasevich Ju.Ju. Perkoljacija: teorija, prilozhenija, algoritmy. Moskva,2002.

15. Utakaeva I.H., Kochkarov R.A. K voprosu ob algoritmah raspoznavanija predfraktal'nogo grafa . NTV SPbGPU. Sankt-Peterburg, 2010.

16. Shreder M. Fraktaly, haos, stepennye zakony. Moskva,2005.

17. http://tehnic.by

18. http://armscontrol.ru/course

19. https://publications.hse.ru/articles/97267128.

????????? ?? Allbest.ru