1. Introduction
Since the outbreak of the worldwide COVID-19 pandemic in 2020, the growth of local outbreaks has been controlled to some extent by the aggressive preventive and control measures implemented by the Chinese government and the concerted efforts of the whole nation [
1]. However, due to the emergence of various mutant strains of SARS-CoV-2 that are more infectious, with a significantly shorter transmission interval and higher viral load than the original coronavirus, the number of confirmed cases and deaths caused by coronavirus disease has again continued to rise, and there has even been a trend toward “breakthrough infection” among COVID-19 vaccine recipients in some countries. SARS-CoV-2 variants have been transmitted worldwide, and the need to prevent and control the spread of the epidemic has again been highlighted. Therefore, in the process of global cooperation to combat the epidemic, it is important to construct a more targeted SARS-CoV-2 transmission model in response to the changing virus transmission characteristics and the constantly improving epidemic prevention and control measures.
Traditional epidemic transmission methods are mainly mathematical, generally using differential equations to describe the relationship between changes in individual states, analyzing the characteristics of equilibrium solutions and the relationship between parameters, and obtaining the transmission pattern of infectious diseases in the population [
2]. Currently available simulation models of novel coronavirus pneumonia are mainly based on infectious disease dynamics models. Common infectious disease dynamics models are classified according to the type of infectious disease as SI, SIR, SIRS, SEIR models, etc. [
3]. They are used to divide the population into different groups according to their different health states, and differential equations are used to describe the dynamic evolution process of the number of individuals among different groups. For example, the SEIR model contains four groups of susceptible, exposed, infected, and recovered individuals, and the relationship between them is represented by the conversion rate in the differential equation [
4]. The traditional mathematical models of infectious diseases make some assumptions and simplifications on the complex process and factors of epidemic transmission, which can theoretically analyze and describe the macroscopic pattern of the spread and diffusion of novel coronavirus pneumonia, but at the same time, it also brings limitations and make it difficult to describe the complex process of epidemic transmission in detail [
5]. First of all, the assumption of uniform population mixing makes it difficult for the mathematical model of epidemic transmission to describe the differences of individual micro-behavior and interaction. Secondly, fewer parameters make it difficult to describe the wide variety of factors associated with the spread of an epidemic, especially the complex behavior of humans. Finally, the parameters assigned to the mean in the differential equation have difficulty describing the heterogeneity in the epidemic transmission process, such as the heterogeneous transmission capacity of infectious individuals and the heterogeneous time scale of disease course development [
6].
Based on the above limitations, researchers in China and other countries have improved various intervention assessment models for repeated outbreaks and to expand the mathematical model to predict the development of epidemic. J Huang et al. [
7] developed a global epidemic prediction system by improving the SIR model combined with temperature and humidity. In the system, assuming that there is no immune difference between individuals and the total population during the outbreak period is constant, the total population of each country is divided into the susceptible population (S), the infected population (I), and the cured and dead population (R). Leontitsis A et al. [
8] developed the SEAHIR model by improving the SEIR model and proposed the use of a six-segmented model for current or future epidemics. Teslya A et al. [
9] developed a COVID-19 regional transmission model to predict the incidence, maximum number of diagnoses, and time to maximum number of diagnoses. Ayele W et al. [
10] evaluated and predicted epidemics at different time points in a given region using the SEIR model. Yang Z et al. [
11] introduced the input/output flow on the basis of the SEIR model and divided the original propagation speed parameters into two flow directions based on the traffic data between provinces and cities, so that the population changes during the model deduction process are more realistic. Most of these models were developed by improving the existing mathematical models of epidemic transmission, which still have limitations. (1) Mathematical methods have insufficient research on the dynamic evolution of the infectious disease transmission process, without considering the mobility of the population; (2) the spread behavior is closely related to the individual behavior of the crowd, and mathematical models generally assume that the individual behavior attributes are the same, and It is not evenly mixed in the crowd and lacks the ability to describe individuals at the micro level, and it is difficult to reflect the influence of human heterogeneous social activities and flow patterns on the communication process [
12]; (3) with the increase of constraints, the complexity of the model increases, and the form of the solution is complex and cannot be analyzed [
13]; (4) it is useful to predict the spread of the original novel coronavirus over time without considering stochasticity, spatial effects, contact mode, interaction behavior, and adaptive changes as the environment changes [
14]. In fact, the previous control measures and prevention strategies for the original novel coronavirus are no longer fully suitable for the transmission characteristics of the current or even future mutated strains. The effectiveness of existing inactivated vaccines at blocking infection by mutated strains has also been compromised.
The interactions between human self-control and autonomous decision-making behavior and the spatial environment influence the transmission of the epidemic. With the continuous contact between individuals in the artificial society and the simulation of SARS-CoV-2 transmission based on individual contact, a macroscopic epidemic transmission will be presented and even a new emergence phenomenon may occur. A multi-agent system is an autonomous agent organization that interacts with each other in a shared environment [
15]. Agents are heterogeneous entities used to model different individuals [
16]. The environment is a virtual representation of the region where the agent is located. The interaction between agents and between agents and the environment is simulated through interactive behavior [
17]. Among them, the mobile interaction behavior of agents is often more suitable for simulating SARS-CoV-2 propagation dynamics. Therefore, in the present study, multi-agent modeling theory was combined with the current epidemic situation and virus transmission characteristics to construct a rapid SARS-CoV-2 transmission suppression assessment model under government control, mainly considering the role of the natural environment, individual attributes, and the interactions of people in the process of epidemic transmission. It is not only suitable for the original coronavirus but is also suitable for the variant strain of the coronavirus. It takes the human individual as the modeling object and abstracts the human being as an autonomous agent [
18]; the heterogeneous random behavior of individuals can be described in detail, and agents can be classified according to the heterogeneity of the population. Computer technology is used to build individual models of crowds using agents, define the structure of social relationships and contact behavior patterns between individuals, and combine mathematical models to describe intelligences’ behavior and interaction rules. At the same time, GIS technology is used to describe geospatial information to simulate the interaction between the agent and the spatial environment. To a certain extent, it reproduces the situational, behavioral, and evolutionary characteristics of epidemic spread and prevention and control and then carries out simulation to realize the emergence of group complex behavior.
In this model, the environmental model is refined, the study area is divided, and the population distribution and flow speed are affected by different building groups, climatic conditions, and road traffic. The reactive agent is used to establish an individual population model, and the attributes, behaviors, and interactions of agents are described randomly and heterogeneously. According to the external information input, agents make judgments and update their own decision-making and output behavior and then affect the environment. Intervention measures such as nucleic acid testing, vaccination, and tracking close contacts are introduced, and hospital agents are introduced to affect the time of treatment and the development of the disease course through distance and medical resources. At the same time, the government can sense the spread of the epidemic in real time and adjust the response level and prevention and control policies for major public health emergencies. Based on the dynamic feedback mechanisms between humans and the spatial environment, different intervention measures can be simulated for implementation at different times under government control, and different intervention strategies are provided to improve disease control methods and to analyze the effectiveness of SARS-CoV-2 prevention and control.
3. Multi-Agent Based Simulation for the Rapid Suppression of SARS-CoV-2 Propagation
3.1. Operational Flow of the Simulation System
The behaviors of human society in major public health emergencies follow simple local microscopic rules and complex overall macroscopic mechanisms with nonlinear characteristics. The use of computer modeling and simulation techniques to build a model with strong similarity to urban social systems is important for analyzing the spread, prevention, and control of SARS-CoV-2 in cities. The multi-agent-based government-controlled rapid SARS-CoV-2 propagation suppression assessment model developed in this study can be divided into three phases: the simulation preparation phase, the simulation operation phase, and the simulation data management phase. The operational flow of the model is shown in
Figure 5.
The simulation preparation stage involves the initial configuration of the spatial environment in the model, as well as information regarding the interactions between the agents by providing actionable data. The initialization of the spatial environment data involves converting the ArcGIS map data for Taiyuan city into ASCII format, which can be recognized by the model and used as the operating environment for various types of agents in the model. The agent data are initialized as a series of characteristics of agents according to the corresponding data, which are visualized, and this mainly involves marking government, hospital, university, station, and business-district agents in the model according to real geographic data and dividing the raster near the area to which they belong. Different intervention measures are then taken to prevent virus transmission when resident agents flow to this area. Based on the information described above, the healthy-person agents are distributed at different densities according to the proportions of the population in different municipal districts and important geographical locations, and the latent person agents are introduced into the areas with large flows near the station agents.
The simulation operation stage is crucial for model construction. In each simulation time step, as the environmental data are updated, various agents follow different behavior rules for interactive feedback and their own state transformation to form a self-organization and self-learning mechanism. The corresponding decisions and behaviors are produced, such as the number quarantined, close contacts, death statistics, number of cured people, and the corresponding intervention measures and the response level of the government are changed according to the data. The interaction feedback between agents in the model and the interaction feedback between the agents and space environment constitute a complex evolving spatiotemporal system.
The simulation data management phase ensures the stable operation of the system. The module control unit ensures the orderly operation of the various components of the model. Agent scheduling allocates resources to the agents running in each simulation time step. State control is mainly for controlling the state changes of the various types of agents and the environment in the model. Clock control is mainly responsible for changing the number of simulation steps, and environment control is mainly responsible for the import, change, and export of the environment in the model. These four control mechanisms work together to extract, store, and statistically visualize the data as the model runs. The data analysis module extracts the data from the model runs and then presents them statistically in graphical form.
3.2. Simulation Model Implementation
In the present study, six districts of Taiyuan City, Shanxi Province, were selected as the research area, and simulations were performed according to the real data. The main parameters in the model are as follows. The total population of the Taiyuan district was set to 3.02159 million [
26]. According to data from the seventh Taiyuan population census in 2020, the age composition of the population was set as 15.55% for 0–14 years old, 68.34% for 15–59 years old, and 16.11% for 60 years old and above. According to data regarding the national complete vaccination rate on 17 March 2022 published by the National Health Commission of the People’s Republic of China, the vaccination rate and mortality was set at 87.85% and the case fatality rate at 2.1%. The transmission rate was set at 76%, and the latency was set at 4.4 days [
27]. The mean time of death was set at 16.7 days [
28]. The average cure time was set at approximately two weeks for mild diseases and 3–6 weeks for severe and critical diseases [
29]. The wearing of masks could reduce the infection rate by 80% [
30]. Vaccination reduced the infection rate to 95% [
31]. According to the Novel Coronavirus Pneumonia Prevention and Control Program (Eighth Edition), the average number of nucleic acid detection tests in an isolation period was set at four. Novel coronavirus pneumonia is pneumonia caused by SARS-CoV-2. The number of medical staff and beds were set for 11 hospitals in Taiyuan city area with level 3a and above, which had the ability to treat the new type of pneumonia due to SARS-CoV-2. Each hospital provided 15% of the total number of beds and the total number of medical staff in accordance with the corresponding medical carrying capacity for the novel pneumonia due to SARS-CoV-2. According to the new management standard for setting up hospitals for treating the novel pneumonia due to SARS-CoV-2, the mild patients had to achieve a health care ratio of 1:2.5, with a bed protection ratio of 1:1, and the severe patients had to achieve a medical care ratio of 1:3, with a bed protection ratio of 1:6.
The operating interface of the model is shown in
Figure 6. During the operation of the model, real-time status data statistics were used for resident agents of different ages. According to the different initial population and infrastructure distributions in different municipal districts, the numbers of healthy people with different ages, quarantined people, close contacts, confirmed cases, deaths, and cured people were the outputs. The map part of the interface was the corresponding model space environment with the current operating state for virus transmission and prevention and control. Different landmarks represented the locations of different spatial agents, and the humanoid icons with different colors moving on the map represented resident agents in different states.
The model simulated the virus transmission process and evaluated the prevention and control effects of different intervention strategies under the control of the government based on the complex process of interaction between the change in the state of the resident agents and behavioral decisions, the level of government response to major public health emergencies, the adjustment of intervention policies of different organizations, and the allocation of medical resources, in order to provide a reference for government decision-making regarding urban virus prevention and control.
3.3. Model Validation
In order to verify the accuracy of the model and parameters, a comparison was made with the existing COVID-19 virus prevention and control simulation model [
24] in Taiyuan. The real data for the local confirmed cases in Taiyuan in April 2022 were compared with the simulated data, and the models were all adjusted to uniform parameters according to the real data and imported into the model for inversion. The results of the comparison are shown in
Figure 7.
As shown in
Figure 7b, with the advance of the simulation time step, the deviation between the cumulative confirmed number predicted by the existing COVID-19 virus prevention and control simulation model in Taiyuan and the real value gradually increases. The reason is that most of the transmission of COVID-19 is in crowded areas, and the flow willingness and speed of the crowd will also be affected to a certain extent according to the current level of prevention and control. Furthermore, the large-scale vaccination with COVID-19 vaccine, nucleic acid testing, and tracking close contacts have affected the traditional law of epidemic spread. However, the existing models do not take into account the willingness and velocity of resident agent flow under different levels of government control, and the target points for the flow of resident agents are randomly determined and not targeted. Existing models also do not divide the model space environment. The entire geographical space of Taiyuan is used as the simulation area, and the sparsely populated mountainous areas are not considered. The movement of people is very small, and the risk of virus transmission is small. At the same time, the state of the latent agent in the existing model is directly transformed into a confirmed agent, which does not take into account the importance of tracking close contacts for the spread of the epidemic. The intervention measures of large-scale vaccination of COVID-19 vaccine and nucleic acid testing have not been taken into account, so the simulation data deviates greatly from the real data.
Therefore, in the current model, the geospatial environment is narrowed down to the urban districts with dense population that are prone to outbreaks, and the municipal districts are divided into six regions according to the actual situation. In order to highlight human self-control and the interaction between autonomous decision-making behavior and the spatial environment, a new government agent is added to the current model, which adjusts the response level in real time based on the current epidemic data of the artificial society, and the resident agent also makes autonomous decisions based on the current response level, deciding whether to flow, as well as the flow speed and target area, and when the prevention and control level is higher, the flow speed of the crowd will be slowed down accordingly. At the same time, the hospital agents, stations, shopping areas, and university agents are added to the current model, and they are mapped to the space environment of the model according to the real geographical location and will adjust the corresponding personalized intervention measures according to the conditions of the municipal district where they are located and the current government agents adjusting the response level and prevention and control policies of major public health emergencies. The crowds around the geographical locations such as hospitals, stations, shopping areas, and universities are denser, and the flow will and speed are also different. In the current model, the visiting rate and the visiting time of the hospital agent will be adjusted according to the current level of prevention and control, and the distance between the agent and the hospital will also affect the visiting time. The resident agents contacted by the latent agents during the flow process will be transformed into close-contact agents, and the infection rate will be affected according to the self-protective measures taken before contact. After being tracked, the close-contact agent and the latent agent will be isolated and regularly tested for nucleic acid, so as to improve the accuracy of the model.
Based on the above improvements, as shown in
Table 3, four evaluation indicators are used to evaluate the errors of the two models, namely, the mean absolute error (MAE), the root mean squared error (RMSE), the mean absolute percentage error (MAPE), and the sum of square error (SSE) evaluates the predictive performance of different models in this study. These four indicators are all negative indicators, that is, the smaller the index value is and the smaller the deviation from the true value of the forecast value of the model, the stronger the forecasting ability of the model [
32]. The statistical results of each evaluation index of different models are shown in
Table 4. The MAE of the current model is reduced by 1.03; the RMSE is reduced by 1.45; the MAPE is reduced by 3.45%, and the SSE is reduced by 173. The simulated index numbers for the current model fitted well with the actual data, and they were significantly improved relative to the existing model, and the accuracy of the model is significantly improved, which was closer to the actual data, thereby indicating that the assessment model for the rapid suppression of SARS-CoV-2 transmission under government control is effective.
4. Simulation Scenarios and Analysis of Experimental Results
4.1. Government Response Level Adjustment
In the model, the government agent can make real-time adjustments to the response level and prevention and control policies for major public health emergencies based on the number of resident agents in different states. In order to simulate the spread of the virus and the effects of prevention and control under different response levels by the government, the parameters were set for the government agents in this scenario to simulate the spread of the virus under different response levels and with flows into the same number of latent agents from station agents or airport agents at four levels. According to the four emergency response levels defined by the National Health and Wellness Commission of the People’s Republic of China, the government agents responded at the different epidemic levels based on the relevant research findings and the specific parameters set in the emergency plan, as shown in
Table 5 [
33].
In particular, when the emergency response level of the government agent was at level IV, the prevention and control ability was the weakest. Intervention measures included a self-protection rate of only 20% and a medical treatment rate of 50%. There was no restriction on the flow of people, and the flow rate was the fastest. When the emergency response level was at level III, involving isolation, medical observation, patient treatment, and fixed-point medical treatment for close contacts, the flow of people slowed down, and the SARS-Cov-2 vaccine was researched. At level II, large-scale tracking and isolation, small-scale shutdown and suspension of classes, and the organization of a certain scale of vaccination were added. At level I, large-scale vaccination was conducted, the city was closed, and the population within the effective range of the city stopped moving.
Each experimental simulation was run 500 times, and the average values of repeated simulations were taken to calculate the experimental results. The experimental results are shown in
Figure 8.
According to the experimental data, when the emergency response was at level IV, the number of resident agents in different states continued to rise. When the medical resources were insufficient, the number of agents with diagnosed patient status increased rapidly, and it was difficult to completely cure them. Due to the lack of effective isolation and vaccination, the proportion of severely ill patients was much higher than that of mild-sickness patients. When the emergency response was at level III, the hospital agent adjusted the speed of treatment, strengthened the self-protection awareness of resident agents, slowed down the flow rate, and tracked a certain number of close contacts for effective isolation. The number of confirmed persons, latent persons, and close contacts decreased significantly compared with level IV, but it was difficult to rapidly curb the spread of the epidemic, and the proportion of severely ill patients was not effectively improved. When the emergency response reached level II, the transmission of the virus could be controlled, and the speed of medical treatment and self-protection awareness continued to improve. In addition, half of the population was vaccinated, which decreased the transmission rate rapidly, while the proportion of severely ill cases decreased, the speed of achieving complete cure increased significantly, and the death toll was only one person. When the emergency response reached level I, the closed places were shut down, and classes were suspended to rapidly control the spread of the virus. Large-scale vaccination formed an immune barrier in the population, so the number of latent persons and close contacts decreased significantly, and the confirmed cases were cured in 42 days. The experimental results showed that strong emergency response measures could quickly control the spread of the epidemic, but a strong response level inevitably incurred higher implementation costs and social costs. Therefore, timely evaluations and selecting appropriate prevention and control measures are essential for urban public health management.
4.2. Improvements in Disease Control Methods
The effects of intervention control measures under different government response levels were analyzed in this study. It is difficult to measure the effects of prevention and control measures on SARS-CoV-2 in a real artificial society when a small number of intervention measures are superimposed at the same time. Thus, in the next experiment, three different intervention strategies are designed, as shown in
Table 6.
Intervention strategy A involved implementing different intervention measures at different times, and intervention strategies B and C also included the first five intervention measures. In intervention strategy C, intervention measure 6 was added, comprising prior vaccination to better control the spread of SARS-CoV-2.
In the model, the willingness of the residents to flow was initially set at 0 (median willingness to flow), and the other parameters were consistent with those described in the implementation part of the simulation model. The corresponding prevention and control measures were implemented according to the different intervention strategies in different simulation time steps. The experimental results are shown in
Figure 8.
The awareness of residents about prevention and control was improved due to the effects of the comprehensive intervention measures. The prevention and control of individuals and organizational units was enhanced at all levels, and thus the number of confirmed persons, quarantined persons, close-contact persons, and latent persons decreased significantly and gradually tended towards zero. The number of cured persons increased significantly, and the time required for achieving a cure was shortened. The proportion of patients diagnosed as severe cases decreased. The deaths occurred in people aged over 65 years who had not been vaccinated, thereby indicating that elderly people with complex health conditions were at a high risk of death. The different intervention strategies also affected the increases in the numbers of newly confirmed persons, latent persons, and close-contact persons.
Figure 9a–c shows curves representing the changes in the number of resident agents in different states under the three intervention strategies over time. The experimental results obtained under intervention strategies C and B were lower than the total number under intervention strategy A, and the growth rate was higher under strategy A than the other two strategies, thereby indicating the importance of medical treatment, tracking, isolation, and limiting the timeliness of population movements to rapidly curb the spread of the epidemic. As shown in
Figure 9d, the proportion with mild disease was much higher under intervention strategy C than strategies A and B, thereby indicating that vaccination was the most effective measure for reducing the risk of severe disease. The simulation results showed that the number of agents in different states was effectively controlled under intervention strategy C, thereby indicating that taking multiple measures at the same time could effectively control the spread of SARS-CoV-2 and that vaccinating a large number of people could gradually build an immune barrier. The factors that affect the spread of SARS-CoV-2 are complex and diverse. Thus, it is difficult to effectively control the spread of SARS-CoV-2 by using a single intervention measure. Therefore, it is essential to implement multiple measures at the same time.
5. Conclusions
In the present study, a multi-agent modeling approach based on the current epidemic situation, virus transmission characteristics, real data, and the spatial environment in a realistic artificial society was used to construct a model of rapid SARS-CoV-2 transmission suppression assessment under government control, which can simulate various scenarios and hypotheses and which is suitable not only for traditional SARS-CoV-2 but also for its variant strains. It is another tool for epidemiological modeling and decision-making. The model can import a geospatial map and the corresponding coordinate positions, then adjust the epidemic parameters to be evaluated and simulated through the calculation of the mathematical model; by adding or deleting different types of agents in the agent model library, it can be used to simulate and predict different types of epidemic in different regions. Different from traditional model analysis methods that focus only on macroscopic patterns of epidemic spread, in the model, humans and environmental entities are abstracted as agents, thereby highlighting the effects of human self-control and autonomous behavior, population heterogeneity, and policy changes by environmental entities on suppressing the spread of the virus. This model combines the randomness and spatial effects of epidemic transmission to define crowd contact patterns and interaction behavior rules that allow the agents to make adaptive changes as the environment changes. In order to improve the prediction accuracy of the model, the environmental model was refined, and the intervention measures, such as vaccination, medical treatment, and close-contact tracing, were introduced in combination with the age condition and the self-protection of different reactive agents. In particular, the governmental agent can respond to the changing epidemic situation in real time, and multiple other environmental entities implement different prevention and control measures in a cooperative manner. The effectiveness of multiple combinations of interventions were analyzed and evaluated based on the dynamic interaction feedback mechanism between humans and the spatial environment to simulate virus transmission in complex urban systems and the changes under different interventions. By fitting real data and comparing it with the data from existing models, the accuracy of the model is significantly improved, and it can accurately predict the development of the epidemic in the city within a certain period of time.
We tested the proposed model in two different scenarios, and the simulation results indicated two main findings. First, an adjustment of the response levels and prevention and control policies for major public health emergencies by the government affected the spread of the virus, and stronger emergency response measures could more rapidly control the spread of the epidemic. Second, experimental comparisons of the improved disease control methods showed that implementing different combinations of intervention strategies in the same place at different times was significantly more effective for the prevention and control of SARS-CoV-2. In particular, vaccination was the most effective and economical method for containing the spread of the virus. Therefore, the agile adjustment of governmental prevention and control strategies according to the spread of the epidemic is crucial for controlling the spread of the epidemic, and a “one-size-fits-all” approach to prevention and control should not be adopted at the same time. The effectiveness of policies should be reasonably adjusted according to the actual situation, and governments should have an active role in promoting disease awareness in prevention and control efforts [
34].
The model should be improved by conducting further simulation experiments and considering epidemiological theories, as well as subdividing the population and geographical areas and adopting different prevention and control measures according to the epidemic situations in different geographical areas. It is also necessary to study the spread of the virus due to the cross-domain flow of the population, as well as the negative impacts of different prevention and control measures on social and economic development. More precise and normalized prevention and control measures can be explored to provide more effective decision support for ensuring urban health protection.