Abstract

In order to control the spread of the COVID-19 virus, this study proposes an ARCN-SUTS (air railway composite network susceptible-untested-tested-susceptible) model based on the correlation characteristics of the air railway composite network in mainland China. Furthermore, this study also puts forward a broken-edge decision-making strategy for the purpose of making decision about the edge efficiently broken and avoiding the second outbreak of the virus spread to minimize the economic losses for railway and civil aviation companies. Finally, simulation results demonstrate that the proposed strategy can effectively control the spread of the virus with minimal economic losses.

1. Introduction

Since December 2019, a typical human-to-human epidemic of novel coronavirus has broken out. At the beginning of the outbreak, the Chinese government attached great importance to this event. The Chinese government not only rapidly deployed various medical and rescue departments but also adopted effective joint prevention and control measures. Wuhan was locked down. Not only that, all cities, counties, and districts in China adopted prohibition measures to limit the flow of population to the greatest extent in order to achieve the purpose of curbing the spread of the epidemic. Till 24:00 on June 30, 2020, the total number of confirmed cases in the world was 10,185,374, and the cumulative death toll was 503,862 including 83,534 confirmed cases and 4,634 deaths in China [1, 2]. The average daily number of new cases in China at present has leveled off; however, the possibility of a rebound has not been ruled out [35]. From reference [6], it is known that the total passenger turnover in 2020 is 1,925,147 million person kilometers (pkm), where passenger turnover of Chinese civil aviation and railway is 145,747 million pkm in total accounting for 75.7% of the total passenger turnover. In order to guarantee the recovery of the social and economic order in China, it is necessary to restrict the traffic flow in civil aviation and railway transportation and find a feasible decision-making strategy to control the epidemic spread within a reasonable threshold through the combination transportation network among aviation and high-speed railway (HSR). It is necessary to make full use of the Chinese civil aviation and HSR transportation system and realize the maximization of transportation efficiency and economic benefits.

Studies on the spread of the virus in a complex network model have been going on for decades; references [7, 8] not only introduced classic virus propagation models and classifications but also made comparisons about their advantages and disadvantages. In particular, the worldwide spread of COVID-19 has inspired scholars to study the spreading models of the COVID-19 virus. References [911] introduced that the infected persons show some latent characteristics in the initial stage according to the characteristics of COVID-19 virus propagation. The explored nodes are added into the basic susceptible-infected-recovered (SIR) framework, and the immune loss area and death zone are inserted into the basic SIR framework to establish the susceptible-explored-infected-recovered-susceptible (SEIRS) model. References [12, 13] introduced that the threshold property of the basic reproduction number was introduced for the propagation model of age-dependent SIS (susceptible-infection-susceptible) diseases in heterogeneous networks. The introduction of a popular threshold can accurately grasp the control effect of the network. Reference [14] proved that no critical immune threshold was caused by the infinite connectivity fluctuation of scale-free networks. References [15, 16] analyzed the controllability of complex networks by combining graph theory and using the controllability theory of linear systems. Emergency decision-making on COVID-19 has become one of the hot issues in recent years. The aim of emergency decision-making is to reduce casualties and property losses[17,18]. Emergency decision-making could be viewed as a multiattribute decision-making problem due to many aspects taken into consideration [1921]. Reference [22] investigated a biobjective fuzzy emergency distribution center location model for reflecting the urgency and uncertainty of large-scale emergencies. Reference [23] developed a method for interactive multicriteria group decision-making with probabilistic linguistic term sets and applied it to the emergency assistance area selection of COVID-19 for Wuhan. Furthermore, it is also necessary to strengthen research on dynamic emergency decision-making for the COVID-19 spread based on traffic networks. It is known that a traffic network is a kind of complex network. The spread of the virus in the field of transportation is different from other biological viruses, computer viruses, and information dissemination models over social networks. The mobility of personnel in the field of the transportation network is strong, and the network is highly complex. From a global perspective, economic development has increasingly demanded air-rail combined transport, which has also led to a more complex transportation network. The associated research studies have been globally carried out. From the cooperation perspective, reference [24] proposed a model of aviation and railway intermodal transportation, which can obtain better-integrated transportation services and economic benefits with lower environmental costs. In response to the current situation of China’s aviation and railway [2527] analyzed the complex characteristics of the air railway composite network and constructed a mixed-integer programming model for the air railway composite network with the lowest intermodal cost. It is well known to all that the delay propagation of HSR train or flight over the network can be controlled by using complex network theory and propagation dynamics model during emergencies [2830]. Similarly, the combination of complex network theory and propagation model is applicable to control the spread of biological viruses over the transportation network. References [31, 32] improved the SEIR (susceptible-explored-infected-recovered) and SEI (susceptible-explored-infected) models based on the spreading characteristics of COVID-19 with system dynamics in epidemic areas. However, this study does not propose the strategy of network edge breaking so as to reduce the virus transmission rate and economic loss. Therefore, this study establishes the COVID-19 virus propagation dynamics model over the air railway composite network and studies the related broken-edge decision-making strategies to control the virus propagation among cities. Simultaneously, this study considers the epidemic threshold of the composite heterogeneous network of air and railway to ensure that the virus propagation among cities always lies within the epidemic threshold range over the network. This can control virus propagation over the air railway composite network. This study also studies the corresponding optimal strategy for the economic consideration of the aviation and railway corporations so as to minimize economic losses within the controllable virus propagation threshold range.

There are two main contributions in this study that are summarized as follows:(1)First, based on the susceptible-exposed-infected-susceptible (SEIS) model [33], the air railway composite network susceptible-untested-tested-susceptible (ARCN-SUTS) model is designed to describe the spread dynamics of the COVID-19 virus among city nodes.(2)Second, in order to avoid the second outbreak of the epidemic and minimize the economic loss over the network, the broken-edge decision-making strategy is proposed to control the spread of the virus in the composite network with less time and minimum economic losses.

The rest of this study is organized as follows. The construction and characteristic analyses of the composite network are introduced in Section 2. The ARCN-SUTS propagation dynamics model is developed and presented in Section 3 in detail. Section 4 presents a broken-edge decision-making strategy over the air railway composite network. Section 5 conducts simulations that demonstrate the effectiveness of the broken-edge decision-making strategy. Conclusions are drawn in the Section 6.

2. The Construction and Characteristic Analysis of the Composite Network

With the expansion of the scale of the air railway composite network, it owns higher complexity and lower controllability; therefore, the study of network complexity is particularly important. This section analyzes the characteristics of mainland China HSR subnetwork, aviation subnetwork, and air railway composite network. Data come from the websites [34, 35]. This study considers all nationwide cities at the prefecture level and above including 223 cities with HSR transport capacity and 155 cities with aviation transport capacity. Any city with HSR or civil aviation transport capacity is considered as a city node on the air railway composite network, and 261 city nodes in total over the composite network are obtained. In this composite network, if there is an HSR line between any two neighboring cities or a direct flight between any two cities, a connection between these two cities is created.

Figure 1 shows the network topology of the HSR subnetwork, aviation subnetwork, and air railway composite network in mainland China, respectively, where the dot represents the city, and the line indicates the connection between the two cities.

From the network topologies as shown in Figure 1, it is obvious that the HSR subnetwork has apparent community structure characteristics; however, the connection among the communities is relatively loose. While the aviation subnetwork has no obvious community structure, the air railway composite network effectively combines the advantages of the two transport networks, increases the cluster coefficient of the network, and reduces the average path length. Furthermore, the node degree numbers and ranking for the three network topologies mentioned above are shown in Figure 2.

Figures 2(a)2(c) show ten city nodes with the highest degree of HSR subnetwork, aviation subnetwork, and air railway composite network, respectively. In the HSR subnetwork, Zhengzhou has the largest degree of 17 and 33 city nodes with the degree of 1 accounting for 14.8% of the HSR subnetwork. In the aviation subnetwork, Beijing has the largest degree of 121, and 4 city nodes have the degree of 1 accounting for 2.58% of the aviation subnetwork. The sum of degrees of all city nodes in the HSR subnetwork is 824 and that of all city nodes in the aviation subnetwork is 3,188. According to the 2019 National Statistical Yearbook, the annual passenger volume of the HSR is 2.29 billion and the annual passenger volume of civil aviation is 660 million [6, 36]. Therefore, it can be seen that the transit of civil aviation transportation is more flexible than that of HSR transportation when considering the direct connection among cities. However, HSR transportation has a stronger capability for sending passengers. Moreover, the probability of susceptibility to travel by HSR is slightly lower than that of travel by civil aviation with effective protective measures, but HSR has a greater impact on COVID-19 development than air [37, 38].

According to equation (1), the power-law fitting, P(k), of the degree distribution of the nodes over the network roves that the degree distributions of the HSR subnetwork, aviation subnetwork, and air railway composite network follow the different characteristics of the power-law distribution with different γ, respectively [39, 40]. The degree distribution analyses are shown in Figure 3.

As shown in Figure 3, it can be seen that the three networks follow the power-law distribution, and the index of the HSR subnetwork degree distribution, γ, is 1.5134, the index of the aviation subnetwork degree distribution is 0.422, and the index of the air railway composite network degree distribution is 0.7635. Therefore, they are scale-free and typical heterogeneous networks, which can simulate the virus propagation among cities. At present, although the spread of the COVID-19 virus within mainland China has slowed down, it cannot rule out undetected symptomatic and asymptomatic infections in some cities. The movement of the infected persons over the transportation network may lead other city nodes to be infected. According to reference [41], the prevalence threshold of viruses over heterogeneous networks is defined aswhere <k> is the average degree of the node and <k2> is the average degree of the second moment, that is,

According to equation (2), it can be seen that as the degree of nodes on the network increases, <k2>⟶∞, λc⟶0. When there is virus propagation in the air railway composite network and the average propagation rate on the network v >λc, the virus will exist on the network for a long time. This means that viruses can easily spread in scale-free networks. We believe that when a node’s degree ki= 0, the node is disconnected from other nodes, so the node can be regarded as out of the network. Therefore, the degree of each node over the air railway composite network is greater than 1 and is a positive integer. When the node’s degree k decreases, the prevalence threshold λc over the network will correspondingly increase. The traditional immunization strategies can avoid susceptible units to be infected in advance when the virus spreads; however, the susceptible units described in this study are city nodes and it is not possible to immunize the susceptibility status of city nodes in advance. The method for increasing the prevalence threshold of the composite network mentioned in this study can be adopted to avoid the further spread of the epidemic. The next section will introduce the ARCN-SUTS propagation dynamics model about COVID-19 virus propagating over the network.

3. ARCN-SUTS Propagation Dynamics Model

This study considers that the spread of the COVID-19 virus over the air railway composite network leads to the state change of city nodes as the population flow changes and establishes a virus propagation model.

Model assumptions for the state changes of city nodes are placed in this study.(1)This study considers the air railway composite network among city nodes within the scope of mainland China, regardless of the impact of the foreign entrance(2)The asymptomatic infection has a certain latent nature, and the detection rate for this group in cities is lower than that of the symptomatic group(3)The model is suitable for the stable epidemic situation, preventing the second outbreak of the epidemic over the air railway composite network and a more serious impact(4)assumed that no matter what scale of infection or asymptomatic infection occurs in a certain city node, it will immediately lead to state changes of city nodes and employ corresponding strategies

Population mobility is common among cities, especially among large cities. Generally, city nodes over a composite network are classified into three states including the state vulnerable to infection (susceptible, S), the state existing untested infected persons (untested, U), and a state existing tested infected persons (tested, T). Due to the population movement in a city and the spread of the COVID-19 virus from person to person, the status of the city nodes will change. As shown in the schematic diagram of the COVID-19 virus propagation model among cities in Figure 4, if there is no infected person in a city, it means that the city node is safely located in the susceptible state, S. If an untested infected person appears in a city, the city is located in an untested state, U, at this time. After using certain testing methods, the infected person is found in a city; then, the city is located in a tested state, T. Even if the infected city recovers, it may be infected again within a certain period of time; therefore, this city is still located in the state, S. In this case, this study proposes a “susceptible-untested-tested-susceptible” (ARCN-SUTS) model applied to the air railway composite network, which can better describe the spread of the virus among cities. The schematic diagram of the proposed propagation model is shown in Figure 4.

It is known that COVID-19 infection can be divided into two types: symptomatic infection and asymptomatic infection. The susceptible city node S can change into the state U1 where the symptomatic untested infection persons exist. At this time, the nodes in state U1 begin to infect other nodes over the network with propagation rate αi. The city constantly checks the infected person through some testing methods. If the city node locates in the state U1, it can change from the state U1 to the state T1 with the detection rate, βi, where T1 is the state with the tested symptomatic infection. Once the infected persons are found in the city, they will be immediately treated until all the infected persons are cured. Then, the city returns to the susceptible state S from state T1 with the rate γi. Similarly, the susceptible city node S can also change into the state U2 where asymptomatic untested infection persons exist. The nodes in state U2 begin to infect other nodes over the network with the propagation rate ai. If the city node is located in the state U2, it can change from the state U2 to the state T2 with the testing rate, bi, where the state T2 is the tested asymptomatic infection. Then, the city returns to the susceptible state S from state T2 with rate ci. U1 and U2 are hidden states, and T1 and T2 are dominant states.

The propagation dynamic model of ARCN-SUTS is shown as follows:whereαi is the symptomatic infection rate of each city, 0 < αi < 1, iϵ {1, 2,…, num}, and num is the number of city nodes.βi is the symptomatic infection testing rate in each city associated with the virus testing level of city i, 0 < βi < 1, iϵ {1, 2,…,num}.γi represents the recovery rate of a city for the city changing from testing symptomatic infection state into a susceptible state associated with the medical treatment level of city i, 0 < γi < 1, iϵ {1, 2,…,num}.ai refers to the asymptomatic infection rate of each city, 0 < ai < 1, iϵ {1, 2,…,num}.bi is the testing rate for asymptomatic infection in each city associated with the testing level of city i, 0 < bi < 1, iϵ {1, 2,…,num}.ci is the recovery rate for the city changing from the testing asymptomatic infection state to a susceptible state associated with the medical treatment level of city i, 0 < ci < 1, iϵ {1, 2,…,num}.τi is the time delay. It takes some time for city i to recover from infection states, T1 or T2. τi refers to the shortest time for all infected persons in city i to be cured.αi and ai represent the infectious rate for symptomatic infection and asymptomatic infection, respectively. These parameters also indicate city state change due to the population flow over the transport network. The mobility of viruses among city nodes can be regarded as the mobility of viruses carried by all travelers. The greater the external passenger flow, the higher the external propagation capacity of the node city.

This section uses the SIR virus propagation model among people and the real data for COVID-19 to obtain the parameters for SIR and then considers other factors to obtain the parameters more accurately mentioned above. This section will describe the procedures to obtain the corresponding parameters in detail. It is known that there are 5 types of states among people, namely, susceptible, explored, asymptomatic infected, symptomatic infected, and recovered. Among them, explored, asymptomatic, and symptomatic infected persons are classified as infection status.

Model assumptions for the propagation model among people are placed as follows:(1)In reality, patients in the incubation stage are classified as infections(2)The model only considers the patients who are tested, treated, and cured in the same city and does not consider the patients transited from other cities(3)The whole curing period for a city starts from the first patient’s treated to the last patient cured(4)The model does not consider reinfection cases

Figure 5 depicts the SIR propagation model of the COVID-19 virus in people.S, I, and R are defined as the susceptible, the infectious, and the recovered.mi is the infection rate of the virus in city i, 0 < mi < 1, iϵ{1, 2,…,num}. Ni is the cure rate of infected persons in city i, 0 < ni < 1, iϵ{1, 2,…,num}.

Therefore, the dynamic model for COVID-19 virus SIR propagation among people is established as follows:

According to the daily number of diagnosed and recovered people in each city, the number changes of the infected persons in all cities are obtained. Figure 6 shows the analysis of the changes in the number of infected persons in some cities.

As shown in Figure 6, it takes about 15 days from the first tested infection to the peak number for each city, and the number gradually decreases with the increase in isolation measures and continuous treatments. The corresponding fitted parameters mi and ni by using the least-squares method are shown in Table 1.

According to references [4244], the degree of the city nodes over the composite network, k, and total passenger flows, f, can determine the relationships among αi and mi, and γi and ni. The relationships are described in equation (6), where fi is the flow of passengers passing through the city i. The corresponding αi and γi values for the selected cities are shown in Table 2.

Next, as reported, the relationships between ai, αi, and bi, βi are described as

The testing rate βi of the city i is related to city size, which is characterized by city GDP (data coming from website [6]). βi is calculated according to equation (8).where Dn is the number of the current daily testing in China, pop is the current population of China, GDPi is GDP of city i, and is the average GDP of all cities in China. Therefore, βi will float according to the city size. As reported, the cure rate of asymptomatic infection is slightly higher than that of symptomatic infection; therefore, in this study, γi = (1 − r)ci, r is a random number, and 0 < r < 0.1.

4. Broken-Edge Decision-Making Strategy

Based on the analysis of the degree of the air railway composite network node in Section 2, it is known that the greater the degree k is, the smaller the epidemic threshold λc is. In order to avoid the second outbreak of the epidemic, this study puts forward a broken-edge decision-making strategy and analyzes the controllability of the epidemic spread over the air railway composite network. Among them, the schematic diagram of the broken-edge decision-making strategy is shown in Figure 7.

As illustrated in Figure 7, during the epidemic duration, when a symptomatic or asymptomatic infection is tested in the city i over the air railway composite network, the state of city i changes into state T1 or state T2; then, cut off the edges are connected with city i according to the strategy. If all the infected persons in city i recover, simultaneously, the broken-edge cities j are not located in states T1 or T2, then the state of city i returns to state S and reconnects the broken edges among cities i and j. The broken-edge strategy over the air railway composite network can ensure network controllability; however, it is inevitable to cause economic losses for aviation and railway companies. Therefore, it is necessary to minimize economic losses on the premise of virus propagation controlled. The objective function and constraints are shown as follows:where eij is the unit freight rate from city i to city j, is the total transport volume from city i to city j, and xij is a binary variable. If the connected edge between city i and city j is broken, xij = 1; otherwise, it is 0.

The ergodic optimization algorithm of the broken-edge decision-making strategy is adopted to minimize the economic loss over the air railway composite network, and the strategy algorithm is illustrated in Table 3.

Based on the broken-edge decision-making strategy proposed in this study, it can effectively cut off the connection between the outbreak city i and some susceptible cities j. In Section 5, the simulations are conducted to analyze the effectiveness of the proposed strategy.

5. Simulation Analyses

The air railway composite network with 261 city nodes proposed in Section 2 is considered, and cities with different states are randomly set. The initial numbers of cities for different states are shown in Table 4.

The state transition of each city for the next moment is determined by using the ARCN-SUTS propagation dynamics model (3). The transmission rate and recovered rate in the model (3) are obtained by using fitting model (4) with actual data. Then, the experiment carries out the broken-edge decision-making strategy mentioned in Section 4, and Monte Carlo runs are conducted 32,769 times.

Under the same initial conditions, Figure 8 demonstrates the number change of cities in different states over the air railway composite network with and without the broken-edge decision-making strategy, respectively. As seen from Figures 8(a) and 8(b), without adopting the broken-edge decision-making strategy, the epidemics will continue to spread across the entire network and cannot be completely controlled. When using the broken-edge decision-making strategy, all cities can be moved to the state S after approximately 256 days. It further demonstrates that the broken-edge decision-making immunization strategy can effectively shorten the recovery time of all cities to susceptible status and thus control the spread of the epidemic by the simulation comparisons.

As shown in Figure 9, among the 32,769 simulations, the number of epidemic-controlled cases accounts for 90.78% in total, where the number of epidemics controlled within 100 days cases accounts for 72.69%, and the number of epidemics uncontrolled cases accounts for 9.22%. Therefore, this model effectively controls the COVID-19 epidemics. The ergodic algorithm was further used to calculate the shortest time for all city nodes recovered to the susceptible state S in the air railway composite network with the broken-edge decision-making strategy.

Under the control strategy with the shortest epidemic duration, Figure 10 demonstrates the number changes of cities in various states. It is obvious that the shortest time for all city nodes recovered to state S is 14 days. The broken-edge strategies for some city nodes are listed in Table 5.

According to the data analysis in Table 5, in the initial stage of the outbreak, the appropriate number of edges was cut off on the nodes of cities in states T1 and T2 to quickly curb the outbreak. When the virus epidemic becomes stable, the number of broken edges will be gradually reduced, and finally the air railway composite network will be stable and safe. The average times to control the epidemics for different cases are obtained through the broken-edge decision-making strategy as shown in Table 6. Figure 9 shows that the broken-edge decision-making strategy has a 90.78% probability of controlling the duration of the epidemic within 100 days, and the average duration of the epidemic is 42.0358 days as shown in Table 6.

The shortest epidemic duration with this strategy is 14 days as shown in Figure 10. The minimum time is far less than the average duration for controlling epidemics within 0∼100 days as shown in Table 6, so it can be concluded that the broken-edge decision-making strategy control effect is better. Furthermore, on the premise of controlling the average propagation rate on the network, , that is,  < λc, the objective and constraint functions (8) will be taken into account, and the total economic losses are minimized for aviation and railway companies by using the broken-edge decision-making strategy.

Partial areas of minimum economic losses for aviation and railway companies with broken-edge decision-making strategies are shown in Figure 11. The red dot represents the least economic losses over the whole time frame of the epidemic spread simulation, and it is the global optimization point. Figure 12 demonstrates the variation of the city node number in each state with minimum economic loss broken-edge decision-making strategy. As illustrated in Figure 11, the total minimum loss is 944.9 million RMB yuan. Correspondingly, it takes 21 days to completely control the epidemics. From another perspective, the total loss is 2.55 billion RMB yuan with the shortest epidemic duration, and it takes 14 days to completely control the epidemics.

Table 7 lists the average economic losses for aviation and railway companies during different time frames. According to Table 7, the minimum economic loss happens within 100 days under the epidemics controlled, which is far less than that of other time frames. As analyzed in Tables 5 and 7, the average losses within the short time frame are less than those within the long time frame. Therefore, it is recommended to cut off the multiple edges in the initial stage within a reasonable range. As the epidemic becomes stable, it can reduce the number of broken edges in order to control the spread of the COVID-19 virus in the early stage and less economic losses for airlines and railway companies.

6. Conclusions

To control the spread of the epidemics over the air railway composite network and minimize the economic loss of airlines and railway companies, this study proposes an efficient decision-making strategy after the outbreak of COVID-19. The main works of this study are summarized as follows. First, this study considers an air railway composite network composed of prefecture-level cities and above in mainland China and analyzes the degree and degree distribution of the networks, which verify the scale-free characteristics of the networks. The traffic flow among city nodes can be more accurately simulated over the air railway composite network. Second, this study establishes the ARCN-SUTS model based on the propagation characteristics of the COVID-19 virus among city nodes over the air railway composite network, which more objectively and rationally describes the variations of city states. Third, the broken-edge decision-making strategy is proposed to control the spread of the virus in the composite network with less time and minimize economic losses. It effectively solves the problem of dynamic emergency decision-making for the spread of the COVID-19 epidemic based on traffic networks. The effect of the strategy is introduced from the economic benefits and the time benefit assessment for controlling COVID-19. Furthermore, the corresponding simulation results show the effectiveness of the broken-edge decision-making strategy for controlling the spread of COVID-19 among cities over the air railway composite network. However, this study has not conducted a comprehensive study on the control methods of epidemic transmission between urban nodes through transportation networks. Therefore, it is not yet possible to effectively make a comparative analysis among multiple control methods. The focus of future research will enhance the development of control strategies in this area.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.