Simulation of Spatiotemporal Relationship between COVID-19 Propagation and Regional Economic Development in China

: Coronavirus disease 2019 (COVID-19) propagation in 2019 posed serious threats and challenges to human public health and safety. Currently, COVID-19 is still not effectively controlled in certain countries and regions. This study combines the traditional susceptible-exposed-infectious-recovered (SEIR) model, system dynamics (SD) model, and cellular automata (CA) model to construct a spatiotemporal dynamics pandemic model (SDPM). The SDPM is used to dynamically depict the spatiotemporal diffusion and outbreak of COVID-19 through research on the relationship between epidemic spread and regional development. The results show that: (1) There is a positive correlation between regional development and epidemic spread. The more developed the regional economy, especially in areas with short-range population migration from Hubei Province, the more severe the epidemic spread; and (2) The spatial isolation and control measures adopted by the government can effectively prevent the COVID-19 spread. The results explore the relationship between COVID-19 spread and regional economic development by studying the spatial and temporal transmission characteristics of COVID-19, and provide a scientiﬁc reference for the government to formulate reasonable response measures.


Introduction
In December 2019, the first patient with coronavirus disease 2019 (COVID-19) was diagnosed in Wuhan [1,2]. Wuhan is the capital of Hubei Province in China, and is well known as the "Nine Province Thoroughfares". Wuhan is the center of the Triangle of Central China and is one of the most important aviation, railway, and highway hubs in China. By the end of 2019, the permanent population of Wuhan reached 11.212 million [3]. Owing to the dual effects of the Spring Festival travel and Wuhan's particular economic and geographical situation, COVID-19 spread rapidly to other provinces in China and certain other countries and regions around the world. Wuhan announced the prohibition of people from leaving Wuhan on 23 January 2020, and lifted the spatial control measures on 8 April 2020 [4]. As of 0:00 on 8 April 2020, a total of 81,802 cumulative confirmed cases, 77,279 cumulative cured cases, and 3333 cumulative deaths had been reported according to the records from 31 provinces (autonomous regions and municipalities) and the Xinjiang Production and Construction Corps [5]. The cumulative deaths due to COVID-19 have far exceeded those due to severe acute respiratory syndrome (SARS) in 2003 [6].
COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2) [7]. The rapid propagation characteristic of COVID-19 poses a serious threat to human life and health [8]. Public health issues have become an important aspect of social security. Acute infectious diseases are widely and quickly spread by simultaneously infecting and affecting multiple individuals. COVID-19 transmission characteristics and development of social security. Acute infectious diseases are widely and quickly spread by simultaneously infecting and affecting multiple individuals. COVID-19 transmission characteristics and development trends have attracted wide attention and research by domestic and foreign scholars [9][10][11][12][13]. Scholars have used a variety of methods to construct different epidemic simulation models to seek efficient and rapid preventive measures. Simulations can visually illustrate the rules for epidemic diffusion and enable the exploration of COVID-19 trends under various policy interventions. The simulation results can also provide a quantitative scientific basis and reference for the government, enabling those in power to take reasonable and effective preventive measures as soon as possible.
The susceptible-exposed-infectious-recovered (SEIR) model is a common infectious disease model that usually expresses a complex infectious disease system through multiple equations [14]. In the SEIR model, people are divided into four classes, namely: susceptible (referring to healthy people, who have a certain probability of being infected after contact with an infected person), exposed (referring to people who are infected by the infectious disease, but temporarily cannot propagate the disease to other people), infectious (referring to patients with infectious disease, who can transmit the disease to healthy people), and recovered (referring to patients who are cured) [15,16]. The relationships between these four classes are shown in Figure 1.
There are two main problems with the traditional SEIR model: (1) It is assumed by default that convalescent patients carry antibodies for the infectious disease after being cured. That is, there will be no second infection in recovered patients. However, according to the relevant reports, certain COVID-19 survivors have been reinfected; that is, they can be diagnosed as COVID-19 patients again; and (2) The default-exposed person is partially converted into an infected patient, but the exposed person is not contagious. However, COVID-19 has also exhibited infectiousness during the incubation period. Therefore, the traditional SEIR model does not represent the transmission characteristics of the novel coronavirus.
Based on the above considerations, this study builds an extended SEIR-SD model based on the classical SEIR model and system dynamics (SD) model, which is used to simulate and predict the macro-quantity spread trends for COVID-19. The cellular automata (CA) model is used to construct an SEIR-CA model to simulate patients' spatiotemporal distribution. The CA model has been widely used in simulation studies of epidemics [17][18][19]. The CA model can easily express the spread relationship between patients and surrounding persons. Through simulation research on patients' distribution in time and space, the epidemic spread is explored at the micro-level. The spatiotemporal dynamics pandemic model (SDPM), which is based on the SEIR-SD-CA models, can reveal the propagation process of infectious diseases more intuitively and effectively. The model can also explore the laws of COVID-19 spread, both in macro and microscales of time and space. Section 2 introduces and builds the extended research methods. Section 3 describes the data used in this study, as well as presenting the simulation results and discussion. The conclusions are drawn in Section 4.

Extension of SEIR-SD Model
The SD model mainly integrates system theory, cybernetics, and computer simulations. It dynamically explores the interaction and influence by constructing complex causal relationships between various elements. Through a combination of quantitative There are two main problems with the traditional SEIR model: (1) It is assumed by default that convalescent patients carry antibodies for the infectious disease after being cured. That is, there will be no second infection in recovered patients. However, according to the relevant reports, certain COVID-19 survivors have been reinfected; that is, they can be diagnosed as COVID-19 patients again; and (2) The default-exposed person is partially converted into an infected patient, but the exposed person is not contagious. However, COVID-19 has also exhibited infectiousness during the incubation period. Therefore, the traditional SEIR model does not represent the transmission characteristics of the novel coronavirus.
Based on the above considerations, this study builds an extended SEIR-SD model based on the classical SEIR model and system dynamics (SD) model, which is used to simulate and predict the macro-quantity spread trends for COVID-19. The cellular automata (CA) model is used to construct an SEIR-CA model to simulate patients' spatiotemporal distribution. The CA model has been widely used in simulation studies of epidemics [17][18][19]. The CA model can easily express the spread relationship between patients and surrounding persons. Through simulation research on patients' distribution in time and space, the epidemic spread is explored at the micro-level. The spatiotemporal dynamics pandemic model (SDPM), which is based on the SEIR-SD-CA models, can reveal the propagation process of infectious diseases more intuitively and effectively. The model can also explore the laws of COVID-19 spread, both in macro and microscales of time and space. Section 2 introduces and builds the extended research methods. Section 3 describes the data used in this study, as well as presenting the simulation results and discussion. The conclusions are drawn in Section 4.

Extension of SEIR-SD Model
The SD model mainly integrates system theory, cybernetics, and computer simulations. It dynamically explores the interaction and influence by constructing complex causal relationships between various elements. Through a combination of quantitative and qualitative methods, the SD model can be used to simulate the behavior and function of complex systems. In this work, an extended SEIR-SD model based on the principles of the traditional SEIR model and SD model was constructed. The relevant formulas in the traditional SEIR model can be expressed as Equations (1)- (5). The number of susceptible persons, S t+1 , at time t + 1 can be expressed as: where S t reflects the number of susceptible persons (i.e., healthy people); N t represents the total population in the study area; I t is the number of infectious patients at time t; γ indicates the average number of persons contacted by infectious patients each day; and β is the infection probability while a susceptible person is in contact with infectious patients. The number of exposed persons, E t+1 , at time t + 1 can be expressed as follows: where E t is the number of exposed persons at time t; and σ is the probability of an exposed person transforming into an infectious patient. The number of infectious persons, I t+1 , at time t + 1 can be expressed as: where ω is the probability of an infectious patient transforming into a recovered person. The number of recovered persons, R t+1 , at time t + 1 can be calculated as: where R t is the number of recovered persons at time t.
The total population N t in the study area at time t can be written as: There are two disadvantages of the traditional SEIR model, namely: (1) By default, it is assumed that the survivors carry antibodies for this type of infectious disease after being cured, which means the survivors will no longer be infected again; and (2) The default is that the exposed person is partially transformed into an infectious patient, but does not have infectiousness.
In this study, the SEIR model was extended as follows: (1) After infected persons are cured, they are automatically converted into healthy people and do not carry COVID-19 antibodies. Some recovered patients may still be reinfected and become exposed people again, and other patients become susceptible people, which means there is still the possibility of being infected; (2) Exposed persons are also infectious, and healthy people are also at risk of being infected after contact with exposed persons; and (3) Based on the original four classifications, death is added to the extended model. Some infected patients eventually die because they are incurable. The expanded relationship is shown in Figure 2. The number of susceptible persons, , at time t + 1 can be extended as: The extended SEIR model relational expressions are expressed as Equations (6)- (11). The number of susceptible persons, S t+1 , at time t + 1 can be extended as: where γ 1 indicates the average number of persons in contact with infectious patients each day; γ 2 indicates the average number of persons contacted by exposed persons each day; β 1 is the infection probability while susceptible people contact infectious patients; β 2 is the probability of infected people contacting exposed persons; and ϕ is the probability of recovered persons transforming into infected patients. The number of exposed persons, E t+1 , at time t + 1 can be extended as follows: The number of infectious persons, I t+1 , at time t + 1 can be extended as: where µ is the death probability of infected patients. The number of recovered persons, R t+1 , at time t + 1 can be calculated as: The number of deaths, D t+1 , at time t + 1 can be expressed as: The total population, N t , in the study area at time t can be written as: According to Equations (6)- (11), an SEIR-SD-CA causality loop diagram is constructed for COVID-19 ( Figure 3), which is used to conduct dynamic simulation research on the epidemic situation.

Rules of the SEIR-CA Model
The CA model originates from complex science [20] and is a bottom-up microscopic dynamic model that is discrete in both time and space [19]. The CA model has the ability to simulate the spatiotemporal evolution of complex systems and reveals macroscopic

Rules of the SEIR-CA Model
The CA model originates from complex science [20] and is a bottom-up microscopic dynamic model that is discrete in both time and space [19]. The CA model has the ability to simulate the spatiotemporal evolution of complex systems and reveals macroscopic complex phenomena through the evolutionary behavior of each cell on the microscopic unit. There are five basic elements in the CA model: cell, cell state, cell space, neighborhood, and transition rule.
(1) The cell is the most basic and smallest unit in the CA model and is the object of model research.
(2) The cell state is a series of data used to represent the cell's current attributes. Each cell has only one state at a certain time.
(3) The cell space is the collection of spatial network points where the cells are located. In this study, the cell space refers to the scope of the study area.
(4) The neighborhood consists of the other cells within a certain range around the central cell and is an important component of the CA model. Common neighborhoods include the Von Neumann neighborhood and Moore neighborhood. The Von Neumann neighborhood includes four basic adjacent cells, and the Moore neighborhood includes eight adjacent cells [21].
(5) The transition rule is the core element of the CA model. Based on the current cell state and the current neighborhood cell state, the next moment cell state is determined by the transition function. The traditional CA transition function is: where S t+1 represents the cell state at time t + 1; S t is the cell state at time t; and N t refers to the neighborhood cell state at time t.
In this study, based on COVID-19 propagation characteristics, the CA transition rule is extended as: where Pop t represents the population distribution at time t; GDP t is the economic development at time t; Migw t indicates the population migration from Wuhan to the cell at time t; Migh t is the population migration from Hubei Province to the cell at time t; and Ran t represents sudden random conditions at time t.
The above relationship can be expressed as: where ∑ 8 i=1 N i,t is the influence probability of the neighborhood cells at time t; P pop,t is the influence probability of the population distribution at time t; P GDP,t is the probability of impact on economic development at time t; P Migw,t is the influence probability of population migration from Wuhan to the target cell at time t; P Migh,t represents the probability of population migration from Hubei Province to the target cell at time t; and P Ran,t indicates the probability of random events at time t.

Model Implementation
The SDPM was constructed and implemented based on the NetLogo platform and Logo language. NetLogo was first developed by Northwestern University in 1999 [22] and is a free and open-source software package that has been continuously updated. NetLogo can build an easy-to-operate graphical user interface and implement module functions based on a specific programming language, which is highly suitable for building complex and dynamic models.
The SDPM built on NetLogo has a friendly interface that can be conveniently operated (as shown in Figure 4). In this model, various simulation scenarios can be automatically simulated and run by inputting the initial operating parameters. Moreover, changes in Land 2021, 10, 599 6 of 15 the epidemic distribution can be dynamically monitored in real time, as well as dynamic changes in macro data. For map display data, the model outputs data in standard ASCII format, which can be easily imported into ArcGIS and other software for further data analysis and mining. For macro statistics, the model outputs data in a standard CSV format. The output of the standard data results provides a basis for further research.
Logo language. NetLogo was first developed by Northwestern University in 1999 [22] and is a free and open-source software package that has been continuously updated. NetLogo can build an easy-to-operate graphical user interface and implement module functions based on a specific programming language, which is highly suitable for building complex and dynamic models.
The SDPM built on NetLogo has a friendly interface that can be conveniently operated (as shown in Figure 4). In this model, various simulation scenarios can be automatically simulated and run by inputting the initial operating parameters. Moreover, changes in the epidemic distribution can be dynamically monitored in real time, as well as dynamic changes in macro data. For map display data, the model outputs data in standard ASCII format, which can be easily imported into ArcGIS and other software for further data analysis and mining. For macro statistics, the model outputs data in a standard CSV format. The output of the standard data results provides a basis for further research.

Study Area and Data
In this study, we used COVID-19-related national statistics data (excluding Hong Kong, Macao, and Taiwan) from 23 January 2020 (when Wuhan announced a measure to forbid departure from Wuhan) [23]

Study Area and Data
In this study, we used COVID-19-related national statistics data (excluding Hong Kong, Macao, and Taiwan) from 23 January 2020 (when Wuhan announced a measure to forbid departure from Wuhan) [23] to 8 April 2020 (when the government lifted the travel limit measures) [24]. The data included cumulative confirmed cases, existing confirmed cases, daily increased confirmed cases, cumulative recovered cases, daily increased recovered cases, cumulative deaths, daily increased deaths, cumulative suspected cases, existing suspected cases, daily increased suspected cases, cumulative close contacts, existing medical observers, and daily medical observers ( Figure 5). These data were obtained from the National Health Commission of the People's Republic of China.
In addition, we obtained the migration rates from Wuhan to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan), and the migration rates from Hubei Province to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan). The data were spatially vectorized and obtained from the Baidu epidemic real-time big data report (https://voice.baidu.com/act/newpneumonia/newpneumonia, accessed on 8 April 2020). In this study, the average level of population migration from Wuhan to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan) before Wuhan's closure of all channels was obtained, as well as the average level of population migration from Hubei Province to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan) before and after the Spring Festival travel season, respectively ( Figure 6).
Land 2021, 10, 599 7 of 15 limit measures) [24]. The data included cumulative confirmed cases, existing confirmed cases, daily increased confirmed cases, cumulative recovered cases, daily increased recovered cases, cumulative deaths, daily increased deaths, cumulative suspected cases, existing suspected cases, daily increased suspected cases, cumulative close contacts, existing medical observers, and daily medical observers ( Figure 5). These data were obtained from the National Health Commission of the People's Republic of China. In addition, we obtained the migration rates from Wuhan to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan), and the migration rates from Hubei Province to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan). The data were spatially vectorized and obtained from the Baidu epidemic real-time big data report (https://voice.baidu.com/act/newpneumonia/newpneumonia, accessed on 8 April 2020). In this study, the average level of population migration from Wuhan to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan) before Wuhan's closure of all channels was obtained, as well as the average level of population migration from The spatial location data of daily COVID-19 patients were collected in this study from Tencent Kandian (https://ncov.html5.qq.com/community, accessed on 20 March 2020) and Nandu Media (https://m.mp.oeeee.com/h5/pages/v20/nCovcase/, accessed on 20 March 2020). The national GDP and population spatial distribution data were obtained from the Resource and Environment Data Cloud Platform (http://www.resdc.cn/Default. aspx, accessed in 2015) [25]. Detailed data are shown in Table 1.
In this study, all spatial data were unified into the Universal Transverse Mercator (UTM) coordinate system. Moreover, to eliminate the influence of different dimensions, the migration rate, population distribution data, and GDP distribution data were normalized according to the following normalization equation [26,27]: Hubei Province to other prefecture-level cities (excluding Hong Kong, Macao, and Taiwan) before and after the Spring Festival travel season, respectively ( Figure 6). The spatial location data of daily COVID-19 patients were collected in this study from Tencent Kandian (https://ncov.html5.qq.com/community, accessed on 20 March 2020) and Nandu Media (https://m.mp.oeeee.com/h5/pages/v20/nCovcase/, accessed on 20 March 2020). The national GDP and population spatial distribution data were obtained from the Resource and Environment Data Cloud Platform (http://www.resdc.cn/Default.aspx, accessed in 2015) [25]. Detailed data are shown in Table 1.

Simulation Results of COVID-19 Propagation in Macro Quantity
In the SEIR-SD model, the time was set to 1 d (one day). The values of the parameters were determined by references [28][29][30][31] and the law of historical data development. The specific values are presented in Table 2. Based on the SEIR-SD model, we simulated the trend of cumulative confirmed cases, cumulative cured cases, and cumulative death cases with respect to the time from 23 January 2020 to 8 April 2020 (Figure 7). Through a visual comparison, it can be seen that the simulation results are relatively close to the real data. Simultaneously, we quantitatively calculated the simulation precision (as shown in Figure 8). It can be seen from Figure 8 that, in addition to some deviations in the simulation effect of the cumulative recovered cases, the cumulative confirmed cases and cumulative death cases both achieved good simulation results. It can be seen from the figures that the epidemic has been basically under control since March in China. The number of confirmed cases and death cases did not increase significantly. Moreover, with the continuous investment in medical staff, medical facilities, and construction of cabin hospitals, the cure rate has been increasing, while the death rate has gradually decreased.

Scenario Simulations for COVID-19 Spatial Propagation
In this study, three different simulation scenarios were considered: (1) according to the current development trend, the government adopted a certain degree of interference measures, such as traffic control; (2) the epidemic spread freely, that is, the government did not take any spatial control measures; and (3) the authorities took measures to further strengthen urban management and control of Wuhan and even the entire Hubei Province.
The cell state is the person's state who is in the cell, namely: susceptible, exposed, infectious, and recovered. This study only considers the research on the COVID-19 spatial spread; that is, for the time being, only the transition from a healthy person to a confirmed patient is considered, whereas the healing cell state and death cell state are not taken into account. Additionally, owing to the limitation of data acquisition, we only obtained the spatial location data for certain confirmed patients from 6 February 2020 to 20 March 2020. Therefore, only some sample data were used to approximately estimate the COVID-19 spatial spread trend.

Scenario 1: Current Development Trend
Based on the space control measures taken by the government at this stage, the epidemic spread trend in the country (excluding Hong Kong, Macao, and Taiwan) centered on Hubei Province was explored. The simulation results are presented in Figure 9. It can be seen from the figure that the epidemic is centered on Hubei Province, and gradually spread to surrounding provinces, especially in areas close to Hubei Province, such as Beijing, the Pearl River Delta, and the Yangtze River Delta. The spatial traffic control measures implemented by the government have significantly prevented the spread of the epidemic to more cities and regions, and effectively limited the speed and scope of It can be seen from the figures that the epidemic has been basically under control since March in China. The number of confirmed cases and death cases did not increase significantly. Moreover, with the continuous investment in medical staff, medical facilities, and construction of cabin hospitals, the cure rate has been increasing, while the death rate has gradually decreased.

Scenario Simulations for COVID-19 Spatial Propagation
In this study, three different simulation scenarios were considered: (1) according to the current development trend, the government adopted a certain degree of interference measures, such as traffic control; (2) the epidemic spread freely, that is, the government did not take any spatial control measures; and (3) the authorities took measures to further strengthen urban management and control of Wuhan and even the entire Hubei Province.
The cell state is the person's state who is in the cell, namely: susceptible, exposed, infectious, and recovered. This study only considers the research on the COVID-19 spatial spread; that is, for the time being, only the transition from a healthy person to a confirmed patient is considered, whereas the healing cell state and death cell state are not taken into account. Additionally, owing to the limitation of data acquisition, we only obtained the spatial location data for certain confirmed patients from 6 February 2020 to 20 March 2020. Therefore, only some sample data were used to approximately estimate the COVID-19 spatial spread trend.

Scenario 1: Current Development Trend
Based on the space control measures taken by the government at this stage, the epidemic spread trend in the country (excluding Hong Kong, Macao, and Taiwan) centered on Hubei Province was explored. The simulation results are presented in Figure 9. It can be seen from the figure that the epidemic is centered on Hubei Province, and gradually spread to surrounding provinces, especially in areas close to Hubei Province, such as Beijing, the Pearl River Delta, and the Yangtze River Delta. The spatial traffic control measures implemented by the government have significantly prevented the spread of the epidemic to more cities and regions, and effectively limited the speed and scope of COVID-19 propagation. Owing to the limitation of data acquisition, the real data we obtained at each time point are sample points data. At the same time, different time points cannot guarantee that the sample points are the same, such that accurate calculations cannot be performed. However, we can see through the visual comparison method that the spatial distribution of simulated data and real data on 20 March 2020 is essentially the same. Assuming that at the beginning of the epidemic, the government takes space isolation control measures to prohibit unnecessary movement of people, the epidemic spread is shown in Figure 11.
Under these control measures, the epidemic is controlled to a great extent. The epidemic mainly propagates in neighboring provinces and economically developed areas around Hubei Province, but does not spread in other provinces and regions across the country.

Discussion
In this study, we calculated the macro-statistical data for the partial spatial distribution of confirmed patients in the above three scenarios. Figure 12 shows the results of the comparison. It is obvious that, from 6 February 2020 to 20 March 2020, the number of confirmed patients increased to 38,684 according to the current trend. Without any government intervention, which means that personnel can naturally flow between provinces and cities, the number increased to 85,007, which is more than twice the number under Scenario 1. If in the early stage of the epidemic the authorities took traffic control measures to prohibit the free movement of people in Wuhan and Hubei Province, the number of confirmed patients would be 28,258, which would greatly control the epidemic spread.

Scenario 2: Free Diffusion
Under this scenario, assuming that the authorities do not take any effective control measures, people flow freely. Population migration was calculated based on the migration rate prior to the closure of Wuhan. The epidemic spread effect is shown in Figure 10.
It can be seen from Figure 10 that, in the context of no government intervention, the epidemic spreads rapidly across the country. Combining these results with Figure 6 (the migration rate from Wuhan and Hubei Province to other prefecture-level cities), it is observed that in the provinces and cities that experienced short-range migration from Wuhan, there gradually appeared confirmed patients and spread to other prefecture-level cities, especially the Xinjiang Uygur Autonomous Region, Inner Mongolia Autonomous Region, and Three Northeastern Provinces.

Scenario 3: Tighter Control Measures
Assuming that at the beginning of the epidemic, the government takes space isolation control measures to prohibit unnecessary movement of people, the epidemic spread is shown in Figure 11.
Under these control measures, the epidemic is controlled to a great extent. The epidemic mainly propagates in neighboring provinces and economically developed areas around Hubei Province, but does not spread in other provinces and regions across the country.

Discussion
In this study, we calculated the macro-statistical data for the partial spatial distribution of confirmed patients in the above three scenarios. Figure 12 shows the results of the comparison. It is obvious that, from 6 February 2020 to 20 March 2020, the number of confirmed patients increased to 38,684 according to the current trend. Without any government intervention, which means that personnel can naturally flow between provinces and cities, the number increased to 85,007, which is more than twice the number under Scenario 1. If in the early stage of the epidemic the authorities took traffic control measures to prohibit the free movement of people in Wuhan and Hubei Province, the number of confirmed patients would be 28,258, which would greatly control the epidemic spread.

Conclusions
Based on the SD model and the CA model, the SDPM is proposed by extending the traditional SEIR model. This model can dynamically simulate the macroscopic quantitative change characteristics of COVID-19 over time, such as cumulative confirmed patients, cumulative cured cases, and cumulative deaths. Meanwhile, the model can also describe the spread of COVID-19 in space. The model realized the multi-angle and mul-

Conclusions
Based on the SD model and the CA model, the SDPM is proposed by extending the traditional SEIR model. This model can dynamically simulate the macroscopic quantitative change characteristics of COVID-19 over time, such as cumulative confirmed patients, cumulative cured cases, and cumulative deaths. Meanwhile, the model can also describe the spread of COVID-19 in space. The model realized the multi-angle and multi-dimensional real-time dynamic tracking and description of COVID-19 both in time and space and at macro and micro scales.
By exploring the relationship between the epidemic transmission characteristics and regional development power in time and space, it was found that the more developed the region, the more serious the epidemic spread. It was also shown that the greater the population migration from Wuhan and Hubei Province, the more severe were the effects of the epidemic. At the same time, in this study, three different simulation scenarios regarding the government's spatial traffic control were established. The results showed that the government's spatial traffic control measures, the prohibition of unnecessary movement of people, and the timeliness of control measure implementation will largely curb the spread of the epidemic to other regions.
COVID-19 propagation is complex and variable, with many influencing factors. In this study, due to the limitations of data acquisition, only certain confirmed patients' spatial distribution data were used for simulation research, and the acquisition time was only from 6 February 2020 to 20 March 2020. If it were possible to acquire longer time series and more large-scale spatial data of confirmed patients, the results would be more accurate and effective for obtaining epidemic spread characteristics.