Location Optimization Using a Hierarchical Location-Allocation Model for Trauma Centers in Shenzhen, China

Abstract

Trauma is considered a “modern civilized sickness”, and its occurrence substantially affects all of society, as well as individuals. The implementation of trauma emergency systems in cities with young, prosperous, and highly mobile populations is necessary and significant. A complete trauma emergency system includes both low-level trauma centers that offer basic emergency services and high-level trauma centers that offer comprehensive services. GIS and operational research methods were used to solve the location problem associated with these centers. This study analyzed the spatial distribution characteristics of trauma demands and the candidate locations of trauma centers based on a spatial analysis and presented a hierarchical location-allocation model for low- and high-level trauma centers in Shenzhen. The response, coverage, treatment and cost capacities of the trauma center locations were considered, and an ant colony optimization was used to calculate the optimal solution. The objectives of this study were to optimize trauma center locations, improve the allocation of medical trauma resources and reduce the rate of deaths and disabilities due to trauma.

1. Introduction

Trauma, which is a severe injury, is typically caused by a violent attack or an accident, including injuries from falls, crushing injuries, gunshot wounds, and traffic accidents, among others. There are more than five million trauma-related deaths each year, accounting for 9% of all deaths [1,2]. In China, the number of trauma-related deaths is more than 75,000 annually, and the average age of patients is between 18 and 40 years [3,4,5,6]. For younger adults, trauma is the primary cause of death.

Trauma occurs suddenly, and experts have shown that trauma emergencies should be treated in the “Golden Hour” to reduce the rate of death and disability [7,8,9,10,11]. A trauma emergency includes three stages: pre-hospital emergency, emergency disposal and intensive care unit (ICU) observation [12]. In the pre-hospital emergency stage, firefighters, police, and medical personnel should arrive on the scene of an incident in a timely manner and offer basic life support [13]. In remote or underserved areas, helicopters are required to transport trauma patients [14]. In the emergency disposal stage, patients can receive more comprehensive service from trauma centers. For a severe trauma patient who has damaged multiple organs or multiple systems, the nearest trauma center should rapidly evaluate and stabilize the patient and then prepare the patient for transport to a high-level trauma center [15]. In the ICU observation stage, less severe cases require only basic observation, which can be provided in low-level trauma centers, whereas severe cases require intensive ICU observation that is provided only at high-level trauma centers [16]. However, trauma emergency systems have not matured in China, and the variability in medical environments as well as the confusion related to the allocation of medical resources affect the ability of trauma patients to receive emergency services in time, leading to high trauma-related mortality. Thus, it is crucial to optimize the locations of trauma centers to improve the trauma emergency system.

According to the description of a trauma emergency system, optimally locating trauma centers can be considered a complex problem [17,18,19]. On the one hand, trauma emergency systems are nested hierarchy systems that offer different types of services for trauma patients of different levels [20]. Thus, the problem is a hierarchical location problem that includes optimizing the locations of low-level and high-level trauma centers. On the other hand, it is a multi-objective location problem for high-level trauma centers, which must consider the principles of fairness, efficiency, economy, applicability and so on [21]. First, the locations of high-level trauma centers should significantly affect the capacity to respond to demands (i.e., the ability to provide timely emergency services and reduce deaths from accidents). Second, the locations of high-level trauma centers should cover a wide area, encompassing trauma hot spots as much as possible, while also minimizing the number of blind spots. Additional coverage ensures efficient accident control. Third, general hospitals that are candidates for high-level trauma centers should have strong medical care capacities to offer comprehensive treatment. Fourth, the costs of establishing high-level trauma centers should be considered because high-level trauma centers require substantial investments of labor and resources. Therefore, while meeting the maximum demands, the number of high-level trauma centers should be minimized.

Several models have been established to solve different location problems of emergency facilities. In general, these models can be categorized as basic location models and comprehensive location models. A basic location model considers only a single-objective problem. For example, Lu [22] presented a generalized weighted vertex P-center Model aimed at minimizing the maximum demand-weighted travel time between facilities and demands. The model plays an important role in locating urgent relief distribution centers. Shariat-Mohaymany et al. [23] presented a modified Location Set Covering Model for emergency service vehicles such as ambulances. The model calculates the minimum number of required ambulances for each demand that can access a certain area within a specified coverage time. Indriasari et al. [24] presented a Maximal Service Area Model that is aimed at maximizing the total service area of a specified number of facilities to optimally locate emergency facilities. The model was applied to solve the location problem of fire stations in South Jakarta, Indonesia. Comprehensive location models consider multi-objective problems that include the relationships and trade-offs among response time, coverage, cost, accessibility and capacity in location planning. For example, Beheshtifar and Alimoahmmadi [25] presented a multi-objective location-allocation model that considers four objectives: minimizing the total travel cost, inequity in clinic accessibility, land-use incompatibility in the study area, and the costs of land acquisition and facility establishment at optimum locations for new clinics. Baray and Cliquet [26] presented a hierarchical location-allocation model that combines the Maximum Covering Model and the P-center Model to determine the optimum locations for three levels of maternity hospitals in France. The model locates different levels of maternity hospitals in a manner that covers the entire territory and maximizes perinatal care accessibility. Widener and Horner [27] presented a hierarchical capacitated-median model based on GIS and spatial optimization strategies to determine relief distribution points. The model aimed to minimize the travel costs of different demand levels associated with an appropriate facility that provided various relief services.

In this study, a hierarchical location-allocation model is presented for trauma centers in Shenzhen. The Location Set Covering Model was used to calculate the smallest number of low-level trauma centers that could meet the demands of the covered area in the optimal time. A multi-objective model that included response, coverage, treatment and cost capacities was used to solve the location-allocation problem of high-level trauma centers. The weight of each factor was obtained using an analytical hierarchy process from a previous study [28]. Additionally, the model resolution step relied on ant colony optimization (ACO). This approach has been used to integrate interdisciplinary knowledge, such as the combination of geographical knowledge and medical knowledge, and study spatial epidemiology [29,30,31,32,33]. Similarly, the study incorporated the characteristics of trauma emergencies into the location decisions of trauma centers. The goals of the study were to improve the trauma emergency system in Shenzhen, to improve the allocation of medical trauma resources, to reduce the rates of death and disability due to trauma and to provide reference findings for planning departments.

2. Materials and Methods

2.1. Study Area

Shenzhen, formerly known as Baoan County, is one of the Special Economic Zones located in southern Guangdong Province, China, and it lies due north of Hong Kong. Since 2011, the city has evolved into 10 administration districts: Futian, Luohu, Yantian, Nanshan, Baoan, Guangming, Longhua, Longgang, Pingshan and Dapeng (Figure 1).

The beautiful environment, booming economy, preferential policies and comprehensive facilities of Shenzhen attract many young people to the area annually. According to the Shenzhen Statistical Yearbooks [34], in 2013, the population was more than 15 million, including approximately 10 million residents between the ages of 20 and 50. Extensive studies [15] have shown that the highest incidence rate of trauma occurs in individuals between 20 and 50 years of age; thus, because the high proportion of the population in Shenzhen is within this age range, the area has a high potential for trauma cases. In addition, a complex road network, insufficient medical resources and limitations in service quality increase the level of difficulty in obtaining emergency treatment. Consequently, optimizing the locations of trauma centers in Shenzhen is extremely urgent and significant.

Shenzhen currently has no trauma centers. In cities, trauma centers, especially high-level centers, typically belong to large- or medium-sized general hospitals because existing medical resources can be reused and patients can move from pre-hospital care to trauma care seamlessly. Therefore, the problem of locating trauma centers can be considered a combination of selecting of the most appropriate hospitals as trauma centers and the determination of their locations.

2.2. Data Description

2.2.1. Trauma Data

The trauma data from 2013 provided by the Shenzhen Center for Health Information totaled 50,241 records. These data included the call time, response time, age, gender, condition, diagnosis, address, and other information. To acquire the locations of trauma events, address geocoding [35] was used to convert the trauma attribute data into spatial data. Then, a spatial statistical analysis of the trauma data was conducted in each district. Figure 2 shows that more trauma cases occurred in Baoan and Longgang, and fewer occurred in Yantian and Dapeng.

Table 1. Numbers of trauma cases in each district in Shenzhen.

District	Male Patients (Persons)	Female Patients (Persons)	Total Patients (Persons)	Density (Persons/Square Kilometer)
Shenzhen	36,110	14,171	50,281	25.80
Baoan	11,827	4031	15,858	42.05
Longgang	6781	2754	9535	24.62
Longhua	4222	1526	5748	32.71
Futian	2957	1434	4391	60.53
Nanshan	3070	1371	4441	26.19
Luohu	2205	1196	3401	43.07
Guangming	2747	885	3632	23.39
Pingshan	1340	505	1845	11.05
Yantian	538	283	821	11.13
Dapeng	423	186	609	2.09

shows that the most trauma cases (15,858) occurred in Baoan, while the fewest trauma cases occurred in Dapeng (609). The highest density of trauma cases occurred in Futian (60.53 per square kilometer), while the smallest density occurred in Dapeng (2.09 per square kilometer).

In addition, a spatio-temporal change analysis was performed using the trauma data. Figure 3 shows that the mean centers that describe the central tendency of trauma data from January to December of 2013 were extremely close. These centers were all located in southeastern Longhua. Meanwhile, the standard deviational ellipses that describe the dispersion characteristics of trauma data were also extremely similar: the long axis runs from northwest to southeast and short axis runs from northeast to southwest. In other words, the monthly distribution characteristics of trauma data are similar, and optimizing the locations of trauma centers based on the data is reasonable and reliable.

2.2.2. Hospital Data

Hospitals in China are divided into three levels ranging from the lowest level (Level I) to the highest level (Level III). Each level has three classes (A, B, and C). When the level and class of a hospital are higher on these scales, then the medical facilities, technical level and health management capabilities of the hospital are better. In 2013, Shenzhen had 143 hospitals, including 16 Level III hospitals, 29 Level II hospitals, 41 Level I hospitals and 57 hospitals that did not achieve any of the grading standards.

Because the hospital that a trauma center relies on must be a comprehensive hospital, candidate hospitals that were selected for trauma centers are shown in Figure 4, in which the red symbols represent Level III hospitals, orange symbols represent Level II hospitals, and blue symbols represent Level I hospitals. Figure 4 and

Table 2. Numbers of candidate hospitals in each district in Shenzhen.

District	Level III (No.)	Level II (No.)	Level I (No.)
Shenzhen	10	16	33
Baoan	1	4	4
Longgang	2	3	13
Longhua	0	2	2
Futian	3	0	6
Nanshan	1	1	1
Luohu	1	2	2
Guangming	1	2	0
Pingshan	0	1	2
Yantian	0	1	1
Dapeng	1	0	2

show that there are large regional disparities in medical service resources among the districts in Shenzhen. Most of the high-level hospitals are in Futian and Luohu, and fewer hospitals are located in Yantian, Pingshan, and Dapeng in eastern Shenzhen.

2.2.3. Road Data

According to their functions within the road network, city roads were divided into seven levels in Shenzhen: highway, expressway, arterial road, secondary trunk road, branch road, residential street and relief road. The average vehicle speed on urban roads was determined based on the code for the design of urban road alignment (

Table 3. Design of vehicle speed on urban roads in Shenzhen.

Road Level	Road Name	Specified Speed (km/h)	Design Speed (km/h)
1	highway	60–120	120
2	expressway	60–80	80
3	arterial road	60–30	60
4	secondary trunk road	50–20	50
5	branch road	40–20	40
6	residential street	0–20	20
9	relief road	0–10	10

) [36], and the fastest speed of the range for each level was used as the design speed because of the priority given to ambulances [37].

In the process of optimizing the locations of trauma centers, travel time is an important parameter. Compared with using Euclidean distance divided by average speed to obtain the travel time between facilities and demands, using the network distance divided by the corresponding speed is more reasonable and accurate.

2.3. Methodology

2.3.1. Location-Allocation Models

A location-allocation model is a method used to determine the optimal locations of facilities and to spatially allocate facilities based on the demands. Classical research regarding location-allocation models includes the P-center Model [38], the P-median Model [38] and the Covering Model. The Covering Model includes the Location Set Covering Model [39] and the Maximal Covering Location Model [40]. The major difference between these models is associated with the criteria or objective function used to solve location problems. However, as the complexity of location problems increases, the basic location models noted above—which consider only a single objective—cannot solve the problems; thus, the concept of a comprehensive location model was proposed. Comprehensive location models are usually based on a basic location model and are established according to different location objectives and the hierarchical attributes of facilities.

For example, urban shelters are hierarchical facilities because they are divided into different levels according to the types of services they offer, including emergency shelters, temporary shelters and long-term shelters. At the first level are emergency shelters, which offer basic services for sudden disastrous accidents. At the next level are temporary shelters, which offer relatively comprehensive services to solve crowding problems experienced by the emergency shelters and to mitigate the effects of secondary disasters. At the third level are long-term shelters, which offer the most comprehensive services for victims whose housing was seriously damaged. Low-level facilities offer low-level services, and high-level facilities offer high-level services. Because of the different requirements, a fixed location-allocation model cannot solve this type of problem, and personalized location-allocation models of each facility level constitute a hierarchical location-allocation model [41].

In general, a hierarchical facility system can be separated into two types: nested hierarchy systems and non-nested hierarchy systems. In a nested hierarchy system, the high-level facilities offer both high- and low-level services. In a non-nested hierarchy system, each facility level offers different services. A trauma center system is a nested hierarchy system in which low-level facilities serve as first responders to offer basic emergency services, while high-level facilities offer not only basic emergency services but also services that low-level facilities cannot offer for trauma patients in serious condition. Therefore, the primary purpose of a low-level trauma center is to cover all demands within a certain timeframe to guarantee basic emergency services and control over the situation “in time”, where “in time” reflects the emergency standard of the “platinum 10 min”, which is the maximum time in should take trauma emergency services to arrive at a trauma scene. The location selection of high-level trauma centers requires the consideration of not only the minimum total weighted travel time from the low-level to the high-level trauma centers but also the area coverage capability, medical treatment capability and construction costs.

In this study, we use the Location Set Covering Model to solve the location-allocation problem of low-level trauma centers, which is used to locate emergency service facilities such as fire stations and hospitals. The aim of this model is to determine the minimum cost required to construct the facilities and to determine their locations so that each demand is covered by at least one facility in a specified threshold (distance or time). When the cost of constructing each facility is equal, the aim is to minimize the number of eligible facilities and to determine their locations.

The objective function is described as follows:

$$\mathrm{Minimize}z={\displaystyle \sum _{j\in J}{c}_j{x}_j}$$

(1)

$$\mathrm{Subject}\mathrm{to}{\displaystyle \sum _{j\in N_i}{x}_j}\ge 1\forall i\in I$$

(2)

$$x_j\in \{0,1\}\forall j\in J$$

(3)

where I and J denote the set of trauma demands and candidate low-level trauma centers, respectively; $c_j$ denotes the construction cost when a low-level trauma center is located at j; $x_j=1$ denotes a low-level trauma center that is located at j, and $x_j=0$ otherwise; $N_i=\{j\in {J|t}_{ij}\le T\}$ denotes the set of low-level trauma centers that can cover the trauma demand i within a specified threshold; T denotes the “platinum 10 min”; and $t_{ij}$ denotes the travel time between i and j.

Equation (1) minimizes the total construction costs of the low-level trauma centers; Equation (2) ensures that each trauma demand is covered by at least one low-level trauma center; Equation (3) ensures that the value of the variable is 0 or 1.

Meanwhile, we use a multi-objective model that considers the response, coverage, treatment and cost capacities of the candidate trauma centers to solve the location-allocation problem of high-level trauma centers. Compared to the basic location model, the considerations of the multi-objective model are more comprehensive and specific. The objective function is as follows:

$$\mathrm{Minimize}{z}_1={\displaystyle \sum _{i\in I}{\displaystyle \sum _{j\in J}{w}_i{t}_{ij}{y}_{ij}}}$$

(4)

$$\mathrm{Maximize}{z}_2={\displaystyle \sum _{i\in I}{w}_i{y}_i}$$

(5)

$$\mathrm{Maximize}{z}_3={\displaystyle \sum _{j\in J}{k}_j{x}_j}$$

(6)

$$\mathrm{Minimize}{z}_4={\displaystyle \sum _{j\in J}{c}_j{x}_j}$$

(7)

$$\mathrm{Subject}\mathrm{to}{\displaystyle \sum _{j\in J}{x}_j}=p$$

(8)

$${\displaystyle \sum _{j\in J}{y}_{ij}}=1\forall i\in I$$

(9)

$$y_{ij}-x_j\le 0\forall i\in I,j\in J$$

(10)

$$x_j\in \{0,1\}\forall j\in J$$

(11)

$$y_{ij}\in \{0,1\}\forall i\in I,j\in J$$

(12)

$$y_i\in \{0,1\}\forall i\in I$$

(13)

where I and J denote the set of trauma demands and candidate high-level trauma centers, respectively; $w_i$ denotes the trauma demands at node i; $t_{ij}$ denotes the travel time between i and j; $y_{ij}=1$ if the trauma demand i is assigned to a high-level trauma center j, and $y_{ij}=0$ otherwise; $y_i=1$ if a trauma demand i is covered by a high-level trauma center, and $y_i=0$ otherwise; $k_j$ denotes the medical quality of high-level trauma center j; $c_j$ denotes the construction cost when a high-level trauma center is located at j; $x_j=1$ denotes that a high-level trauma center is located at j, and $x_j=0$ otherwise; and p denotes the number of high-level trauma centers to be located.

Equation (4) minimizes the total weighted travel time between high-level trauma centers and trauma demands; Equation (5) maximizes the coverage of the trauma demands; Equation (6) maximizes the total medical quality of the high-level trauma centers; Equation (7) minimizes the total construction cost of the high-level trauma centers; Equation (8) ensures that the desired number of high-level trauma centers is determined; Equation (9) ensures that each trauma demand is covered by a high-level trauma center; Equation (10) ensures that only the selected high-level trauma centers can provide service for trauma demands; and Equations (11)–(13) ensure that the value of the variable is 0 or 1.

To reduce the number of contradictory objectives, a weighted method is often used to convert the multiple objectives to a single objective. Combining the above objective functions, the final objective function of high-level trauma centers is as follows:

$$\mathrm{Maximize}z=-{\gamma }_1{\displaystyle \sum _{i\in I}{\displaystyle \sum _{j\in J}{w}_i{t}_{ij}{y}_{ij}}}+{\gamma }_2{\displaystyle \sum _{i\in I}{w}_i{y}_i}+{\gamma }_3{\displaystyle \sum _{j\in J}{k}_j{x}_j}-{\gamma }_4{\displaystyle \sum _{j\in J}{c}_j{x}_j}$$

(14)

where ${\gamma }_1$, ${\gamma }_2$, ${\gamma }_3$, and ${\gamma }_4$ are the weight coefficients of each objective function, and the result includes constructing the judgment matrix, inviting an expert to evaluate the function and performing a consistency check [28].

2.3.2. Ant Colony Optimization

Finding the optimal locations of trauma centers is a typical NP problem that requires solving a large number of polynomials. Compared to exhaustive methods, heuristic methods such as ACO can save a great deal of computational time. ACO, which was introduced by Dorigo et al. [42], is a popular intelligent heuristic algorithm that was inspired by the foraging behavior of ants. This type of ant behavior is a collective behavior in which each ant deposits a chemical substance (pheromone) on a path from the anthill to food; the more pheromone that is deposited, the greater the probability that other ants in the community will follow the path. The intensity of the pheromone weakens over time. Eventually, the path with the highest pheromone intensity will be the shortest.

ACO has often been used to solve the Travelling Salesman Problem (TSP), but it has also been applied to solve the Maximum Independent Set Problem (MISP), Knapsack Problem (KP), and Quadratic Assignment Problem (QAP), among others [43,44,45,46]. ACO can also be applied to solve location-allocation problems. Figure 5 shows the flowchart of the ACO algorithm used to solve the location-allocation problem involving high-level trauma centers in this study. The ACO algorithm process is as follows (Figure 5).

Step 1: Initialization of the algorithm (including iternum (maximum iterations), m (number of ants), $\alpha$ (narrowing parameter), $\beta$ (evaporation rate), ρ (volatilization parameter), Q (pheromone intensity), z_min (minimum z), and so on).

Next, import $N_j$ (the set of candidate high-level trauma centers), $N_i$ (the set of trauma demands), $\mathrm{A}=\left[t_{ij}\right]$ (the travel time matrix between trauma demand i and high-level trauma center j), $N_j^i$ (the set of trauma demands that can be addressed from high-level trauma center j within the “Golden Hour”), and Zj (the matrix of candidate high-level trauma centers, which calculated by the response matrix, coverage matrix, treatment matrix and cost matrix of the candidates).

Then, set $\mathrm{iter}=0$ (initial number), $\mathrm{t}={\mathrm{t}}_0=0$ (initial time), $V_j^k\left(\mathrm{t}\right)=\varnothing$ (the set of selected high-level trauma centers at time t is null), $N_j^k\left(\mathrm{t}\right)=N_j$ (the set of unselected high-level trauma centers at time t is universal), $V_i^k\left(\mathrm{t}\right)=\varnothing$ (the set of selected trauma demands at time t is null), $N_i^k\left(\mathrm{t}\right)=N_i$ (the set of unselected trauma demands at time t is universal), $Z_k=0$ (initial value), and $BestN=\varnothing$ (the optimal locations of high-level trauma centers).

Step 2: Start internal iterations, and when the iteration number reaches iternum, go to Step 11.

Step 3: Randomly generate f candidate high-level trauma center IDs for ant k to visit. Simultaneously, update the set of selected/unselected high-level trauma centers and the selected/unselected trauma demands.

Step 4: If the number of selected high-level trauma centers reaches p (the number of required facilities), then go to Step 7.

Step 5: If the set of unselected high-level trauma centers is null, then go to Step 10.

Step 6: Select the high-level trauma center with the highest probability for ant k to visit according to Equation (15). Simultaneously, update t to t+1, and update the selected/unselected high-level trauma centers and selected/unselected trauma demands. Then, go to Step 4.

$$p_j^k(t)=\{\begin{array}{c}\frac{{\tau }_j^{\alpha }(iter){\eta }_j^{\beta }(t)}{{\displaystyle {\sum }_{j\in N_j^k}{\tau }_j^{\alpha }(iter){\eta }_j^{\beta }(t)}},j\in N_j^k(t)\\ 0,j\notin N_j^k(t)\end{array}$$

(15)

In Equation (15), $p_j^k\left(t\right)$ denotes the probability that ant k will visit high-level trauma center j at time t, ${\tau }_j\left(iter\right)$ denotes the pheromone intensity of high-level trauma center j, ${\eta }_j\left(t\right)$ denotes the heuristic function in Equation (15), $N_j^k$ denotes the set of high-level trauma centers that ant k has not visited at time t, α denotes the narrowing parameter, β denotes the pheromone evaporation rate, and iter denotes the number of iterations.

$${\eta }_j(t)=\frac{\left|N_i^k(t)\cap N_j^i\right|}{{\displaystyle {\sum }_{i\in N_i^k(t)}{z}_i}},j\in N_j^k(t)$$

(16)

In Equation (16), $N_i^k\left(t\right)$ denotes the set of trauma demands that ant k has not visited at time t, $N_j^i$ denotes the set of trauma demands that can be covered by high-level trauma center j, and ${\sum }_{i\in N_i^k\left(t\right)}{z}_i$ denotes the number of trauma demands that ant k has not visited at time t.

Step 7: Calculate z. If z < z_min, then z = z_min, and BestN = $N_j^k$.

Step 8: Start the next ant’s visit, k = k + 1; if k ≤ m, then go to Step 3.

Step 9: Update the pheromone intensity according to Equation (16) and then go to Step 2.

$${\tau }_j(iter+1)=(1-\rho ){\tau }_j(iter)+\Delta {\tau }_j$$

(17)

In Equation (17), $\rho$ denotes the volatilization parameter, iter denotes the number of iterations, and ${\Delta \tau }_j$ denotes the sum of the pheromone intensities at high-level trauma center j.

$$\Delta {\tau }_j={\displaystyle {\sum }_{k=1}^m\Delta {\tau }_j^k}$$

(18)

In Equation (18), ${\Delta \tau }_j^k$ denotes the pheromone intensity of ant k at high-level trauma center j.

$$\Delta {\tau }_j^k=\{\begin{array}{c}\frac{Q}{L_K},j\in V_j^k\\ 0,j\notin V_j^k\end{array}$$

(19)

In Equation (19), Q denotes the pheromone intensity, $L_k$ denotes the number of high-level trauma centers that ant k has visited, and $V_j^k$ denotes the set of high-level trauma centers that ant k has visited at time t.

Step 10: Print “Error: The task cannot be solved”, then go to Step 12.

Step 11: Complete the task and print the values of BestN and z.

Step 12: End.

3. Results and Discussion

3.1. Low-Level Trauma Centers

A Location Set Covering Model was used to solve the location problem for low-level trauma centers. In this problem, the demand data were the trauma demands, the facility data were the candidate comprehensive hospitals and the specified travel-time threshold was 10 min. ArcGIS Network Analyst has location-allocation models that include a Location Set Covering Model. Figure 6 shows the location results of low-level trauma centers using the Location Set Covering Model from the ArcGIS Network Analyst Toolbox. The model calculates that at least 46 low-level trauma centers can ensure that each accessible demand is covered by the “platinum 10 min” emergency criterion. The results include 7 Level III hospitals, 17 Level II hospitals, and 22 Level I hospitals. In theory, the model ensures that each demand can be covered by at least one trauma center; however, a few remote demands cannot be covered in 10 min because of the limited medical resources in Shenzhen.

Figure 7 shows the administrative regions of the chosen low-level trauma centers. The figure clearly shows that more trauma cases occur in Baoan and Longgang, where the number of low-level trauma centers should be 8 and 13, respectively. In Yantian and Dapeng, where fewer trauma cases occur, the number of low-level trauma centers should be 2 and 3, respectively.

In addition,

Table 4. Information about low-level trauma centers.

Name	Level	Average Time (min)	District
Nanshan District People’s Hospital	Level III	6.72	Nanshan
Second People’s Hospital	Level III	7.34	Futian
Longgang Centre Hospital	Level III	6.26	Longgang
Dapeng District People’s Hospital	Level III	6.09	Dapeng
Xinming Hospital	Level III	6.82	Guangming
Xinan Hospital	Level III	7.54	Baoan
Baohe Hospital	Level III	7.60	Longgang
Shekou People’s Hospital	Level II	6.15	Nanshan
Liuhua Hospital	Level II	5.36	Luohu
Baoan District Xixiang People’s Hospital	Level II	6.26	Baoan
Longgang District Second People’s Hospital	Level II	5.80	Longgang
Longgang District Henggang People’s Hospital	Level II	5.19	Longgang
Longhua District People’s Hospital	Level II	5.32	Longhua
Pingshan District Pishan People’s Hospital	Level II	5.45	Pingshan
Baoan District Fuyong People’s Hospital	Level II	6.22	Baoan
Longgang District Pinghu People’s Hospital	Level II	5.88	Longgang
Longhua District Guanlan People’s Hospital	Level II	5.69	Longhua
Baoan District Shajin People’s Hospital	Level II	6.18	Baoan
Guangming District People’s Hospital	Level II	5.99	Guangming
Baoan District Songgangn People’s Hospital	Level II	6.54	Baoan
Guangming District Centre Hospital	Level II	6.79	Guangming
Yantian District People’s Hospital	Level II	5.00	Yantian
Ranai Hospital	Level I	6.36	Futian
Dapeng District Nanao People’s Hospital	Level I	4.39	Dapeng
Peace Hospital	Level I	6.10	Futian
Friendship Hospital	Level I	7.00	Futian
Futian District Second People’s Hospital	Level I	7.03	Futian
Nanshan District Xili Hospital	Level I	7.09	Nanshan
Renkang Hospital	Level I	7.46	Luohu
Yantian District Yangang Hospital	Level I	5.71	Yantian
Longgang District Shawan People’s Hospital	Level I	6.24	Longgang
Lily Hospital	Level I	6.70	Longgang
Dapeng District Kuiyong People’s Hospital	Level I	4.87	Dapeng
Baoxing Hospital	Level I	6.70	Longgang
Longan Hospital	Level I	7.04	Longhua
Xuexiang Hospital	Level I	6.58	Longgang
Longji Hospital	Level I	6.98	Longhua
Baoan District Shiyan People’s Hospital	Level I	6.16	Baoan
Shenlian Hospital	Level I	5.05	Longgang
Jianfeng Hospital	Level I	5.45	Pingshan
Dragon City Hospital	Level I	5.47	Longgang
Guangsheng Hospital	Level I	7.18	Baoan
Renan Hospital	Level I	6.62	Longgang
Pingshan District Maternal and Child Health Hospital	Level I	4.37	Pingshan
Fuya Hospital	Level I	6.28	Baoan
Longgang District Pingdi People’s Hospital	Level I	4.02	Longgang

lists information about the selected low-level trauma centers. The average time indicates the average response time between trauma demands and the trauma centers. A majority of the chosen low-level trauma centers have a quick response time, ensuring that basic emergency services can be provided in time.

3.2. High-Level Trauma Centers

The multi-objective model described in Section 2.3.1 was used to solve the location problem of high-level trauma centers, in which the trauma demands were based on demand data, the facility data included 46 low-level trauma centers, the specified travel time threshold was 30 min and the number of high-level trauma centers was eight (because of the rule that specifies one trauma center is required per 2 million people [4]). The ACO algorithm described in Section 2.3.2 was used to solve the model (implemented in MATLAB) using the following parameters: iternum = 100, m = 50, α = 1, β = 5, ρ = 0.1, and Q = 1.5. Note that a standard deviation was used for each parameter value before the calculation to reduce problems stemming from the use of different units.

Figure 8 shows the convergence curves of z (the comprehensive coverage capability of high-level trauma centers). The value of z increased gradually from 1 to 82 iterations. Additionally, an optimal solution of z = 8.359 was obtained after 83 iterations. Then, the solution displayed a continuously stable state until the end of the simulation was reached. Figure 9 shows the optimal solution for high-level trauma centers, which includes four Level III hospitals and four Level II hospitals.

Information about the identified high-level trauma center locations is listed in

Table 5. Information about high-level trauma centers.

Name	Level	Average Time (min)	Coverage Ratio	District
Nanshan District People’s Hospital	Level III	19.45	30.43%	Nanshan
Second People’s Hospital	Level III	19.22	52.17%	Futian
Longgang Centre Hospital	Level III	18.19	26.09%	Longgang
Dapeng District People’s Hospital	Level III	21.49	19.57%	Dapeng
Liuhua Hospital	Level II	20.53	39.13%	Luohu
Baoan District Xixiang People’s Hospital	Level II	21.75	23.91%	Baoan
Baoan District Songgangn People’s Hospital	Level II	20.12	17.39%	Baoan
Longhua District Guanlan People’s Hospital	Level II	22.24	43.48%	Longhua

, in which the average time indicates the average response time between the high-level trauma centers and low-level trauma centers, and the coverage ratio indicates the ratio of the number of low-level trauma centers that each high-level trauma center covers to the sum of low-level trauma centers. The average response time of each high-level trauma center is approximately 20 min; thus, high-level services are provided to serious trauma patients within the “Golden Hour”. Meanwhile, the average coverage ratio is approximately one third, indicating that each high-level trauma center has a wide coverage area.

High-level trauma centers not only offer emergency services to trauma patients nearby but also provide comprehensive services for serious trauma patients who cannot be treated by low-level trauma centers. The allocation of trauma centers for emergency services should increase the efficiency of trauma patient management. Each low-level trauma center treats the nearest high-level trauma center as a higher-grade medical institution. Figure 10 and

Table 6. List of emergency service allocations of trauma centers.

High-Level Trauma Center	Low-Level Trauma Center
Nanshan District People’s Hospital	Shekou People’s Hospital
	Nanshan District Xili Hospital
Second People’s Hospital	Longgang District Second People’s Hospital
	Ranai Hospital
	Peace Hospital
	Friendship Hospital
	Futian District Second People’s Hospital
	Renkang Hospital
	Lily Hospital
Longgang Centre Hospital	Baohe Hospital
	Longgang District Henggang People’s Hospital
	Pingshan District Pishan People’s Hospital
	Jianfeng Hospital
	Dragon City Hospital
	Renan Hospital
	Pingshan District Maternal and Child Health Hospital
	Longgang District Pingdi People’s Hospital
	Baoxing Hospital
Dapeng District People’s Hospital	Yantian District People’s Hospital
	Dapeng District Nanao People’s Hospital
	Yantian District Yangang Hospital
	Dapeng District Kuiyong People’s Hospital
Liuhua Hospital	Longgang District Shawan People’s Hospital
Baoan District Xixiang People’s Hospital	Xinan Hospital
Baoan District Songgangn People’s Hospital	Xinming Hospital
Longhua District Guanlan People’s Hospital	Baoan District Fuyong People’s Hospital
	Baoan District Shajin People’s Hospital
	Guangming District People’s Hospital
	Guangming District Centre Hospital
	Guangsheng Hospital
	Fuya Hospital
	Longgang District Pinghu People’s Hospital
	Longhua District Guanlan People’s Hospital
	Longan Hospital
	Xuexiang Hospital
	Longji Hospital
	Baoan District Shiyan People’s Hospital
	Shenlian Hospital

show the emergency services allocation among both high- and low-level trauma centers, as well as a detailed list of the hospitals.

3.3. Discussion

A trauma emergency system that includes 46 low-level trauma centers and 8 high-level trauma centers was constructed using the hierarchical location-allocation model and an ACO algorithm. The chosen trauma centers exhibit a spatially uniform distribution, and the locations of all the trauma centers cover 54.04% of the city’s area within 10 min, 85.71% of the city’s area within 20 min, and 93.19% of the city’s area within 30 min (Figure 11). Most areas in Futian, Luohu, Longhua, Nanshan and Longgang, which have high trauma demands, can be covered within the specified emergency response time, while only 6.81% of the entire area of the city is out of the emergency range within the specified time. This occurs because of terrain limitations and the limited medical resources in those areas. In conclusion, the trauma emergency system can cover a large service area and provide a quick service time, ensuring that trauma emergencies are treated in a timely manner and that the rate of death and disability due to trauma is reduced. This result shows the feasibility of the proposed hierarchical location-allocation model.

4. Conclusions

The main contribution of this study was to propose a hierarchical model combining a Location Set Covering Model and a multi-objective model that considers the factors of response, coverage, treatment and cost of optimum location candidates for two levels of trauma centers in Shenzhen, China. The previous location models of trauma centers were mostly single-objective models that focused on only a single aspect, such as coverage area maximization. By contrast, the model in this study is more comprehensive because it combines geographical and medical knowledge to provide a more reasonable result. Simultaneously, the location-allocation model based on large trauma demands would generally result in a large number of polynomials in the process of solving the location problem. However, a heuristic method was used to reduce the calculation time and determine the spatial solution. Ultimately, the location-allocation result for trauma centers based on this model and method showed that a large service area can be covered and a quick service time can be provided, which will play a significant role in improving the trauma emergency system in Shenzhen. Thus, the hierarchical model and heuristic method provide a new framework for optimizing the locations of trauma centers.

However, because of data limitations, the trauma demands in this study were based only on trauma data from 2013, and trauma demands have increased in recent years. Therefore, the location-allocation result for the trauma centers shown here may not provide emergency service for each demand in a timely manner (such as in the 6.81% of blind spots). Trauma centers would need to be added and the trauma emergency system would need to be improved to meet the new trauma demands. Meanwhile, a further analysis of additional trauma demands may identify new location objectives, making the model even more reasonable and reliable.