A method and system for selecting the location of a disaster relief supply depot based on earthquakes of magnitude 7 or above

Abstract

The present invention discloses a method and system for the site selection of disaster relief material reserve warehouses based on earthquakes of magnitude 7 or above, relating to the technical field of the site selection of disaster relief material reserve warehouses, including: determining an intensity area where the earthquake magnitude is 7 or above and the earthquake intensity is 8 or above according to the earthquake literature materials and earthquake intensity attenuation models of the research area; using remote sensing building recognition technology and eight-connected mode to process the intensity area to obtain the disaster-affected area; determining the farthest safe flight distance when helicopters transport disaster relief materials and determining the candidate areas for the site selection of disaster relief material reserve warehouses within the research area; screening the candidate areas for the site selection according to the constraint relationship between the farthest safe flight distance and the disaster-affected area, as well as natural disasters, terrain indices, traffic conditions, urban planning, etc. of the research area, so as to determine the optimal site selection area. The address area selected by the present invention can achieve the purpose of practical feasibility.

Description

Technical Field

[0001] The present invention relates to the technical field of the location selection of disaster relief material reserve warehouses, and particularly to a method and system for the location selection of disaster relief material reserve warehouses based on earthquakes above magnitude 7. Background Art

[0002] After an earthquake above magnitude 7 occurs, the roads leading to the disaster-stricken areas will be damaged to varying degrees, and even the disaster relief material reserve warehouses in the city will become ineffective within a short period of time, which results in the inability to deliver disaster relief materials to the disaster-stricken areas through conventional means. After an earthquake above magnitude 7 occurs, disaster relief material reserve warehouses located in appropriate positions can obviously provide effective livelihood material guarantees for the victims in the disaster-stricken areas in a timely and sufficient manner. Therefore, the location selection of large-scale disaster relief material reserve warehouses is very crucial for coping with earthquakes above magnitude 7.

[0003] Currently, the methods for the location selection of disaster relief material reserve warehouses mainly are: based on methods such as the maximum area coverage method, the minimum time method, and the P-median method, taking factors such as the construction cost, transportation cost, and procurement cost of disaster relief material reserve warehouses as constraint conditions, and using genetic algorithms and improved genetic algorithms to construct location selection method models. However, the above-mentioned methods for the location selection of disaster relief material reserve warehouses are all applied to scenarios where the roads are not damaged and the disaster relief costs are highly emphasized, and do not meet the actual requirements. Summary of the Invention

[0004] The purpose of the present invention is to provide a method and system for the location selection of disaster relief material reserve warehouses based on earthquakes above magnitude 7, so as to achieve the purpose of practical feasibility.

[0005] To achieve the above purpose, the present invention provides the following solutions:

[0006] A method for the location selection of disaster relief material reserve warehouses based on earthquakes above magnitude 7, comprising:

[0007] Determining the earthquake risk areas with earthquake magnitudes above magnitude 7 in the research area according to the earthquake literature materials of the research area;

[0008] Determining the intensity areas with earthquake intensities above magnitude 8 in the research area according to the earthquake risk areas and the earthquake intensity attenuation model;

[0009] Identifying the buildings in the intensity areas by using remote sensing building identification technology to generate a building distribution map, and then processing the building distribution map in an eight-connected manner to obtain the disaster-stricken areas;

[0010] Determining the farthest safe flight distance when a helicopter transports disaster relief materials and determining the location selection candidate areas of disaster relief material reserve warehouses in the research area;

[0011] Screen the candidate site selection areas according to the constraints to determine the optimal site selection area; the constraints at least include the constraint relationship between the farthest safe flight distance and the disaster area, as well as natural disasters, terrain index, traffic conditions and urban planning in the research area.

[0012] Optionally, according to the seismic literature data of the research area, determine the seismic risk areas with a seismic magnitude of more than 7 in the research area, specifically including:

[0013] Determine the candidate seismic risk areas in the research area according to the prediction results and statistical laws of seismic scientific research in the research area;

[0014] Determine the maximum seismic magnitude corresponding to each candidate seismic risk area;

[0015] Compare the maximum seismic magnitude with a set threshold, and determine the candidate seismic risk area corresponding to the maximum seismic magnitude greater than or equal to the set threshold as the seismic risk area with a seismic magnitude of more than 7; where the set threshold is 7.

[0016] Optionally, the seismic risk area is an elliptical area; according to the seismic risk area and the seismic intensity attenuation model, determine the intensity areas with a seismic intensity of more than 8 in the research area, specifically including:

[0017] Determine the elliptical isoseismal line when the seismic intensity is 8 according to the seismic risk area and the seismic intensity attenuation model;

[0018] After adding the elliptical isoseismal line to the seismic risk area and then making an elliptical enclosure, obtain the intensity areas with a seismic intensity of more than 8.

[0019] Optionally, the expression of the seismic intensity attenuation model is:

[0020] I = 6.8053 + 1.2972Ms - 4.7603log(R a + 22);

[0021] I = 5.3315 + 1.2013Ms - 4.1917log(R b + 10);

[0022] where I is the seismic intensity, I = 8; Ms is the seismic magnitude; R a and R b are the major semi-axis length and minor semi-axis length of the elliptical isoseismal line when the seismic intensity is 8, respectively.

[0023] Optionally, the specific method for determining the farthest safe flight distance when the helicopter transports disaster relief supplies includes:

[0024] Determine the farthest safe flight distance of the helicopter when transporting disaster relief supplies according to the formula D max =(TSR - FLR) / 2;

[0025] Among them, the constraint conditions for the safe flight of the helicopter are:

[0026] v = MCS;

[0027] q = MP;

[0028] In the formula, D max is the farthest safe flight distance, TSR is the total standard range, FLR is the fuel limit range, MCS is the maximum cruise speed, in km / h, MP is the maximum load, v is the speed of the helicopter during safe flight, and q is the load of the helicopter during safe flight.

[0029] Optionally, determine the candidate areas for the location of the disaster relief supply depot in the research area, specifically including:

[0030] Determine the first area and / or the second area in the research area as the candidate areas for the location of the disaster relief supply depot; the first area is the area determined with the highway entrance as the center and a radius of 5 km; the second area is the area determined with the railway freight station as the center and a radius of 5 km.

[0031] Optionally, the natural disaster includes at least flood; the constraint conditions also include the per capita GDP statistical results; screen the candidate areas for the location according to the constraint conditions to determine the best location area, specifically including:

[0032] Calculate multiple first target candidate areas based on the farthest safe flight distance and the disaster-stricken area, and determine the first target candidate areas that meet the candidate areas for the location as the second target candidate areas;

[0033] Screen the second target candidate areas according to the flood areas and the intensity areas above magnitude 8 in the research area, and determine the second target candidate areas that meet the non-flood areas and the non-intensity areas above magnitude 8 in the research area as the third target candidate areas;

[0034] Screen the third target candidate areas according to the per capita GDP statistical results of the research area in the past two years, and determine the third target candidate areas that meet the per capita GDP statistical results of the research area in the past two years as the fourth target candidate areas;

[0035] Screen the fourth target candidate areas according to the terrain index of the research area, and determine the fourth target candidate areas that meet the terrain index constraints of the research area as the fifth target candidate areas;

[0036] Screen the fifth target candidate area according to the traffic conditions of the research area, and determine the fifth target candidate area that meets the traffic condition constraints of the research area as the sixth target candidate area;

[0037] Screen the sixth target candidate area according to the urban planning of the research area, and determine the sixth target candidate area that meets the urban planning constraints of the research area as the optimal site selection area.

[0038] A disaster relief material reserve site selection system based on earthquakes above magnitude 7, comprising:

[0039] An earthquake risk area determination module, configured to determine an earthquake risk area with a magnitude of 7 or above within the research area according to earthquake literature materials of the research area;

[0040] An intensity area determination module, configured to determine an intensity area with an intensity of 8 or above within the research area according to the earthquake risk area and the earthquake intensity attenuation model;

[0041] A disaster-affected area generation module, configured to use remote sensing building identification technology to identify buildings within the intensity area to generate a building distribution map, and then process the building distribution map in an eight-connected manner to obtain a disaster-affected area;

[0042] A farthest safe flight distance and site selection candidate area determination module, configured to determine the farthest safe flight distance when a helicopter transports disaster relief materials and determine a site selection candidate area for a disaster relief material reserve within the research area;

[0043] An optimal site selection area determination module, configured to screen the site selection candidate area according to constraint conditions to determine an optimal site selection area; the constraint conditions at least include the constraint relationship between the farthest safe flight distance and the disaster-affected area, and natural disasters, terrain index, traffic conditions, and urban planning of the research area.

[0044] Optionally, the earthquake risk area determination module specifically includes:

[0045] An earthquake risk candidate area unit, configured to determine an earthquake risk candidate area within the research area according to the prediction results and statistical laws of earthquake scientific research in the research area;

[0046] A maximum earthquake magnitude determination unit, configured to determine the maximum earthquake magnitude corresponding to each earthquake risk candidate area;

[0047] An earthquake risk area determination unit is used to compare the maximum earthquake magnitude with a set threshold, and determine the earthquake risk candidate area corresponding to the maximum earthquake magnitude greater than or equal to the set threshold as an earthquake risk area with an earthquake magnitude above level 7; wherein, the set threshold is 7.

[0048] Optionally, the natural disaster at least includes flood; the constraint condition further includes the per capita GDP statistical result; the optimal site selection area determination module specifically includes:

[0049] A second target candidate area determination unit is used to calculate a plurality of first target candidate areas according to the farthest safe flight distance and the disaster area, and determine the first target candidate areas that meet the site selection candidate area as the second target candidate areas;

[0050] A third target candidate area determination unit is used to screen the second target candidate areas according to the flood areas and intensity areas above level 8 within the research area, and determine the second target candidate areas that meet the non-flood areas and non-intensity areas above level 8 within the research area as the third target candidate areas;

[0051] A fourth target candidate area determination unit is used to screen the third target candidate areas according to the per capita GDP statistical result of the research area in the past two years, and determine the third target candidate areas that meet the per capita GDP statistical result of the research area in the past two years as the fourth target candidate areas;

[0052] A fifth target candidate area determination unit is used to screen the fourth target candidate areas according to the terrain index of the research area, and determine the fourth target candidate areas that meet the terrain index constraint of the research area as the fifth target candidate areas;

[0053] A sixth target candidate area determination unit is used to screen the fifth target candidate areas according to the traffic conditions of the research area, and determine the fifth target candidate areas that meet the traffic condition constraint of the research area as the sixth target candidate areas;

[0054] An optimal site selection area determination unit is used to screen the sixth target candidate areas according to the urban planning of the research area, and determine the sixth target candidate areas that meet the urban planning constraint of the research area as the optimal site selection area.

[0055] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0056] First, the present invention determines the intensity areas with an earthquake intensity above level 8 through the earthquake literature materials of the research area, which is more conducive to quickly realizing the selection of the disaster area.

[0057] Second, the present invention identifies the buildings in the disaster area through the remote sensing building identification technology, and sets tags in the eight-connectivity manner, effectively ensuring the accuracy of the selection of the disaster area.

[0058] Third, usually, the transportation of disaster relief roads is restricted, while the present invention solves the farthest safe flight distance according to the helicopter model, effectively improving the safety and efficiency of the transportation of disaster relief materials.

[0059] Fourth, the present invention restricts the terrain index, disaster situation, municipal planning, etc. of the candidate site selection area, effectively ensuring the practical feasibility of the site selection. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required to be used in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0061] Figure 1 It is a flowchart of a method for selecting a site for a disaster relief material storage depot based on an earthquake of magnitude 7 or above provided by an embodiment of the present invention;

[0062] Figure 2 It is a result map of the selection of intensity regions of magnitude 8 or above provided by an embodiment of the present invention;

[0063] Figure 3 It is a result map of the marked disaster area provided by an embodiment of the present invention;

[0064] Figure 4 It is a result map of the candidate area after excluding natural disaster areas provided by an embodiment of the present invention;

[0065] Figure 5 It is a result map of the selected reference area according to the urban planning index provided by an embodiment of the present invention;

[0066] Figure 6 It is a structural diagram of a system for selecting a site for a disaster relief material storage depot based on an earthquake of magnitude 7 or above provided by an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0067] The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, rather than all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0068] At present, the site selection method of the disaster relief material reserve adopted combines less with the prediction results and statistical laws of earthquake scientific research, and does not consider problems such as the inaccessibility of roads after an earthquake and the failure of nearby disaster relief material reserves. In view of this, the present invention provides a site selection method and system for disaster relief material reserves based on earthquakes above magnitude 7. First, obtain the disaster-stricken areas with an earthquake intensity above magnitude 8 and above magnitude 7 through the earthquake literature materials in the research area; secondly, determine the farthest safe flight distance through helicopter equipment; then take the areas with a radius of 5 km centered on the highway entrances and railway freight stations in the research area as candidate site selection areas; finally, conduct site selection according to the candidate site selection areas, disaster-stricken areas, and in combination with conditions such as the farthest safe flight distance, disasters, helicopter takeoff and landing, and terrain index, to ensure that the site selection of large-scale disaster relief material reserves can meet the actual needs, and after an earthquake above magnitude 7 occurs in the area, the disaster-stricken areas can also obtain the assistance of livelihood materials in a short time.

[0069] To make the above objects, features, and advantages of the present invention more obvious and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0070] Embodiment 2

[0071] This embodiment provides a site selection method for disaster relief material reserves based on earthquakes above magnitude 7, as Figure 1 shown, including the following steps:

[0072] Step 100: According to the earthquake literature materials in the research area, determine the earthquake risk area in the research area with an earthquake magnitude above magnitude 7.

[0073] This step is specifically:

[0074] First, according to the prediction results and statistical laws of earthquake scientific research in the research area, determine the earthquake risk candidate areas in the research area.

[0075] Secondly, determine the maximum earthquake magnitude corresponding to each earthquake risk candidate area;

[0076] Compare the maximum earthquake magnitude with a set threshold, and determine the earthquake risk candidate area corresponding to the maximum earthquake magnitude greater than or equal to the set threshold as the earthquake risk area with an earthquake magnitude above magnitude 7.

[0077] Among them, the set threshold is 7; the earthquake risk area is an elliptical area.

[0078] Step 200: According to the earthquake risk area and the earthquake intensity attenuation model, determine the intensity area in the research area with an earthquake intensity above magnitude 8.

[0079] This step is specifically:

[0080] Step A: Determine the elliptical isoseismal line when the seismic intensity is 8 degrees according to the seismic risk area and the seismic intensity attenuation model.

[0081] One example is: Determine the extended boundary when the seismic intensity is 8 degrees according to the boundary position of the seismic risk area.

[0082] Determine the elliptical isoseismal line when the seismic intensity is 8 degrees according to the extended boundary and the seismic intensity attenuation model when the seismic intensity is 8 degrees.

[0083] The expression of the seismic intensity attenuation model is:

[0084] I = 6.8053 + 1.2972Ms - 4.7603log(R a + 22);

[0085] I = 5.3315 + 1.2013Ms - 4.1917log(R b + 10);

[0086] where I is the seismic intensity, I = 8; Ms is the magnitude of the earthquake; R a and R b are the major semi - axis length and minor semi - axis length of the elliptical isoseismal line when the seismic intensity is 8 degrees, respectively.

[0087] Step B: After adding the elliptical isoseismal line to the seismic risk area and then performing an elliptical enclosure, obtain the intensity area where the seismic intensity is above 8 degrees.

[0088] One example is: Add the major semi - axis and minor semi - axis lengths of the elliptical isoseismal line with a seismic intensity of 8 degrees above, below, left, and right of the seismic risk area, and then perform an elliptical enclosure to obtain the Figure 2 intensity area where the seismic intensity is above 8 degrees as shown. The intensity area where the seismic intensity is above 8 degrees is an elliptical area.

[0089] Step 300: Use the remote - sensing building identification technology to identify the buildings in the intensity area to generate a building distribution map, and then process the building distribution map in an eight - connected manner to obtain the disaster - affected area. The disaster - affected area is also an elliptical area.

[0090] One example is: Select the IGBP global vegetation classification system, take urban and built - up areas (classification code 13) as the main targets for building area extraction, fill the pixels, design a buffer around this area, then discriminate the isolated pixels around the city through visual interpretation, and include ground object types such as villages and factories in the scope of the disaster - bearing body, and finally generate a building distribution map.

[0091] To eliminate the influence of isolated buildings, erosion and dilation operations are performed on the building distribution map, and the eight-connected method is used as the judgment basis. The connected regions are set as disaster-bearing bodies, that is, the disaster-affected area shown in Figure 3 is obtained.

[0092] Step 400: Determine the farthest safe flight distance when a helicopter transports disaster relief supplies and determine the candidate areas for the location of the disaster relief supply depot within the research area.

[0093] One example is: Given the model of the disaster relief supply transport helicopter, considering the requirements for the safe flight of the helicopter, calculate the farthest safe flight distance when transporting disaster relief supplies.

[0094] The calculation formula for the farthest safe flight distance is as follows;

[0095] D max =(TSR - FLR) / 2;

[0096] The constraint condition for the safe flight of the helicopter is: v = MCS;

[0097] q = MP;

[0098] In the formula, D max is the farthest safe flight distance, TSR is the total standard range, FLR is the fuel limit range, MCS is the maximum cruising speed, in km / h, MP is the maximum load, v is the speed of the helicopter during safe flight, and q is the load of the helicopter during safe flight.

[0099] One example is: Determine the first area and / or the second area within the research area as the candidate areas for the location of the disaster relief supply depot; the first area is the area determined with the highway entrance as the center and a radius of 5 km; the second area is the area determined with the railway freight station as the center and a radius of 5 km.

[0100] Step 500: Screen the candidate areas for the location according to the constraint conditions to determine the optimal location area; the constraint conditions at least include the constraint relationship between the farthest safe flight distance and the disaster-affected area, as well as the natural disasters, terrain index, traffic conditions and urban planning of the research area.

[0101] Among them, the natural disasters at least include floods; the constraint conditions also include the statistical results of per capita GDP; this step specifically includes:

[0102] Step a: Calculate multiple first target candidate areas based on the furthest safe flight distance and the disaster area, and check whether the obtained first target candidate areas meet the requirements, that is, whether they are within the candidate site selection area; if they do not meet the requirements, abandon the first target candidate area; if they meet the requirements, save the first target candidate area and determine the first target candidate areas that meet the candidate site selection area as second target candidate areas.

[0103] Further, the formula for calculating the first target candidate area based on the furthest safe flight distance and the disaster area (the disaster area is also an elliptical area) is as follows:

[0104] a = 6378137m;

[0105] b = 6356752m;

[0106]

[0107]

[0108]

[0109]

[0110] y = K(L - L 0 );

[0111] where L is the longitude, B is the latitude, L 0 is the standard longitude (customized), B 0 is the standard latitude (customized), (x, y) are the coordinates after the conversion of (L, B), a is the semi-major axis of the ellipsoid, b is the semi-minor axis of the ellipsoid, e is the first eccentricity, e' is the second eccentricity, is the radius of curvature of the prime vertical, K is the radius of the parallel circle at the longitude (L 0 , B 0 ).

[0112] Constraint conditions:

[0113]

[0114] In the formula, n represents the number of disaster areas; (x i , y i ) are the coordinates of the disaster area, and (x, y) are the coordinates of the first target candidate area.

[0115] Step b: Considering the flood areas and intensity areas above magnitude 8 within the study area, count whether the second target candidate area meets the non-flood areas and non-intensity areas above magnitude 8 within the study area; if it does not meet the requirements, abandon the second target candidate area; if it meets the requirements, save the second target candidate area, and determine the second target candidate area that meets the non-flood areas and non-intensity areas above magnitude 8 within the study area as the third target candidate area, specifically as Figure 4 shown.

[0116] Further, in step b, the flood area is the area where the flood statistics are ≥ 0.3 times / a.

[0117] Step c: Considering the per capita GDP statistical results of the study area within two years, count whether the third target candidate area meets the requirements; if it does not meet the requirements, abandon the third target candidate area; if it meets the requirements, save the third target candidate area, and determine the third target candidate area that meets the per capita GDP statistical results of the study area within two years as the fourth target candidate area.

[0118] Further, in step c, the per capita GDP statistical result here is not the last one in the city ranking.

[0119] Step d: Considering the topographic index of the study area, check whether the fourth target candidate area meets the requirements; if it does not meet the requirements, abandon the fourth target candidate area; if it meets the requirements, save the fourth target candidate area, and determine the fourth target candidate area that meets the topographic index constraint conditions of the study area as the fifth target candidate area.

[0120] Further, in step d, the formula for calculating the topographic index is as follows;

[0121] Calculation formula for elevation mean:

[0122]

[0123] Calculation formula for elevation variance:

[0124]

[0125] Calculation formula for elevation standard deviation:

[0126]

[0127] Among them, Xi is the elevation of each pixel, and N is the number of pixels in the area.

[0128] The implementation process of the constraint method for the topographic index is as follows:

[0129] Taking a 5×5 window as a unit, the slope of the units contained in the candidate area needs to have units with a slope less than 25%, and the smaller the slope, the more preferred; the average elevation of the units contained in the candidate area is in the top 50% of all units in the study area; the standard deviation of the elevation of the units contained in the candidate area is in the bottom 15% of all units in the study area.

[0130] Step e: Screen the fifth target candidate area according to the traffic conditions of the research area, and determine the fifth target candidate area that meets the traffic condition constraints of the research area as the sixth target candidate area.

[0131] For example: Considering the factor of convenient transportation, calculate the road distances from the fifth target candidate area to the coordinates of each selected area, and sort them from small to large according to the sum of each road distance, and select the first m fifth target candidate areas corresponding to the sum of each road distance, such as the top 10, and save the fifth target candidate areas that meet the requirements, that is, the sixth target candidate area is obtained.

[0132] Further, in step d, the calculation process of the coordinates of the selected area is as follows:

[0133] Take the road transportation hub, the train freight station (or highway intersection) of railway transportation, and the airport closest to the fifth target candidate area as the selected area.

[0134] During the selection process of the selected area, it is necessary to exclude the natural disaster hazard areas involved above. The selection of the transportation hub point needs to be the intersection of two expressways and have at least 2 directions and be determined as a non-seismic hazard area direction.

[0135] First, take the area with the fifth target candidate area as the center of the highway intersection as an example, and calculate the road distances from the fifth target candidate area to the coordinates of each selected area, as shown below:

[0136] a = 6378137m;

[0137] b = 6356752m;

[0138]

[0139]

[0140]

[0141]

[0142] y = K(L - L 0 )

[0143]

[0144] Z = min∑j＝J ∑ f＝F ∑ g＝G ∑ p＝P (h j，f d j，f +h j，g d j，g +h j，p d j，p );

[0145] Decision variables:

[0146]

[0147]

[0148]

[0149] Constraint conditions:

[0150] ∑ f∈F h j，f = 1, j ∈ J;

[0151] ∑ g∈G h j，g = 1, j ∈ J;

[0152] ∑ p∈P h j，p = 1, j ∈ J;

[0153] Symbols for the meanings of variables used in the study:

[0154] d is the sum of the distances of each point of the road network geographical information;

[0155] Q is the sum of the quantities of the road network geographical information;

[0156] q is the road network information point;

[0157] Z is the objective function, that is, the sum of the shortest distances;

[0158] j is the fifth target candidate area, j ∈ J, and J is the set of the fifth target candidate areas;

[0159] h j，f 、h j，g 、h j，p are binary decision variables

[0160] f is a transportation hub, f ∈ F, and F is the set of transportation hubs

[0161] g is an airport, g ∈ G, and G is the set of airports

[0162] p is a railway freight station, p ∈ P, and P is the set of railway freight stations

[0163] d j，f 、dj，g , d j，p is the distance from each selected region to the fifth target candidate region.

[0164] Step f: Screen the sixth target candidate region according to the urban planning of the research region, and determine the sixth target candidate region that meets the urban planning constraints of the research region as the optimal site selection region, specifically as Figure 5 shown.

[0165] For example: Based on the saved candidate regions, estimate the reference candidate regions according to the urban planning indicators; obtain the ideal candidate regions considering the maximum distance from the fire source.

[0166] The urban planning indicators are:

[0167] Taking Baoshan City as an example, according to the DEM slope data of Baoshan City, combined with the requirements of the "Code for Vertical Planning of Urban Land Use", factors restricting the construction of disaster relief material reserve warehouses such as airport clearance restrictions, hydrological conditions, distribution of mineral resources, and helicopter takeoff and landing are incorporated to obtain the regional construction land situation.

[0168] It should be noted that during the process of determining the optimal site selection region, the result output at a certain step may be an empty set, which may cause the inability to screen out the optimal site selection region. In this case, first, reduce the seismic intensity of the intensity reduction region. If the optimal site selection region still cannot be screened out, the remote sensing identification building technology needs to be replaced, and other identification technologies are used to obtain the human activity region, that is, update the disaster-stricken region. If the optimal site selection region still cannot be screened out, the research region needs to be replaced.

[0169] A method for selecting a location for a disaster relief supply depot based on earthquakes of magnitude 7 or above provided by an embodiment of the present invention determines an intensity area with an earthquake intensity of 8 or above according to the prediction of earthquake danger areas; identifies the buildings within the intensity area of 8 or above using remote sensing building identification technology and uses them as disaster-stricken areas; determines the farthest safe distance when a helicopter transports disaster relief supplies according to the helicopter equipment for post-earthquake disaster relief supply transportation. Set the area within the research area with highway entrances and railway freight stations as the center and a radius of 5 km as the candidate location area. Under the condition of the farthest safe distance when a helicopter transports disaster relief supplies, find the candidate areas that meet the distance requirements for reaching all disaster-stricken areas; based on the historical flood statistics data of the research area, divide the candidate areas that meet the frequency requirements; based on the statistical data of the per capita GDP of each county in the city where the candidate area is located, obtain the ranking of the county where the candidate area is located in the city, and obtain the candidate areas that meet the ranking requirements; under the condition of considering the convenience of helicopter takeoff and landing, check and find the candidate areas that meet the terrain index requirements. If the candidate area meets the above requirements, output the land transportation distance of each candidate area to the traffic convenience index. Under the condition of considering the urban planning index, bring the top 20% into the planning map for judgment. The embodiment of the present invention can obtain the location of the disaster relief supply depot based on earthquakes of magnitude 7 or above, ensuring that the location selection work can cope with the transportation and rescue of disaster relief supplies for earthquakes of magnitude 7 or above within the region.

[0170] Embodiment 2

[0171] As Figure 6 shown, a location selection system for a disaster relief supply depot based on earthquakes of magnitude 7 or above provided by an embodiment of the present invention includes:

[0172] An earthquake risk area determination module 1, configured to determine an earthquake risk area with a magnitude of 7 or above within the research area according to the earthquake literature materials of the research area.

[0173] An intensity area determination module 2, configured to determine an intensity area with an earthquake intensity of 8 or above within the research area according to the earthquake risk area and the earthquake intensity attenuation model.

[0174] A disaster-stricken area generation module 3, configured to identify the buildings within the intensity area using remote sensing building identification technology to generate a building distribution map, and then process the building distribution map in an eight-connected manner to obtain a disaster-stricken area.

[0175] A farthest safe flight distance and candidate location area determination module 4, configured to determine the farthest safe flight distance when a helicopter transports disaster relief supplies and determine the candidate location area of the disaster relief supply depot within the research area.

[0176] The optimal site selection area determination module 5 is used to screen the candidate site selection areas according to the constraint conditions to determine the optimal site selection area; the constraint conditions at least include the constraint relationship between the farthest safe flight distance and the disaster area, as well as the natural disasters, terrain index, traffic conditions and urban planning of the research area.

[0177] Among them, the earthquake risk area determination module 1 specifically includes:

[0178] The earthquake risk candidate area unit is used to determine the earthquake risk candidate areas within the research area according to the prediction results and statistical laws of earthquake scientific research in the research area.

[0179] The maximum earthquake magnitude determination unit is used to determine the maximum earthquake magnitude corresponding to each earthquake risk candidate area.

[0180] The earthquake risk area determination unit is used to compare the maximum earthquake magnitude with a set threshold, and determine the earthquake risk candidate area corresponding to the maximum earthquake magnitude greater than or equal to the set threshold as an earthquake risk area with an earthquake magnitude above 7; where the set threshold is 7.

[0181] The natural disasters at least include floods; the constraint conditions also include the per capita GDP statistical results; the optimal site selection area determination module 5 specifically includes:

[0182] The second target candidate area determination unit is used to calculate a plurality of first target candidate areas according to the farthest safe flight distance and the disaster area, and determine the first target candidate areas that meet the candidate site selection areas as the second target candidate areas.

[0183] The third target candidate area determination unit is used to screen the second target candidate areas according to the flood areas and intensity areas above 8 degrees within the research area, and determine the second target candidate areas that meet the non-flood areas and non-intensity areas above 8 degrees within the research area as the third target candidate areas.

[0184] The fourth target candidate area determination unit is used to screen the third target candidate areas according to the per capita GDP statistical results of the research area in the past two years, and determine the third target candidate areas that meet the per capita GDP statistical results of the research area in the past two years as the fourth target candidate areas.

[0185] The fifth target candidate area determination unit is used to screen the fourth target candidate areas according to the terrain index of the research area, and determine the fourth target candidate areas that meet the terrain index constraints of the research area as the fifth target candidate areas.

[0186] The sixth target candidate area determination unit is configured to screen the fifth target candidate area according to the traffic conditions of the research area, and determine the fifth target candidate area that meets the traffic condition constraints of the research area as the sixth target candidate area.

[0187] The optimal site selection area determination unit is configured to screen the sixth target candidate area according to the urban planning of the research area, and determine the sixth target candidate area that meets the urban planning constraints of the research area as the optimal site selection area.

[0188] The present invention determines the earthquake risk area according to literature materials, sets the maximum earthquake magnitude and the destructive intensity level (Level VIII) that may occur in the risk area; determines the intensity level expansion boundary according to the boundary position of the determined earthquake risk area; based on the earthquake intensity attenuation model of the research area, calculates the axis distance of the earthquake Level VIII intensity area, and the intensity area is elliptical; according to the intensity area, uses remote sensing recognition technology to find the buildings in the intensity area; gives the disaster area label by the eight-connectivity method; gives the type of the rescue material transportation helicopter; considering the safety flight requirements of the helicopter, calculates the farthest safe flight distance of the helicopter; gives the areas within 5 km centered on the highway entrances and railway freight stations in the research area as the candidate areas for site selection; based on the disaster area label and the farthest safe flight distance; considering the flood statistics results of the research area and the Level VIII intensity area, counts whether the candidate areas meet the requirements; considering the county per capita GDP statistics results within two years, counts whether the candidate areas meet the requirements; considering the terrain index, checks and calculates whether the candidate areas meet the requirements; considering the traffic convenience factor, calculates the road distances from the candidate areas to the coordinates of each selected area; sorts according to the sum of each; based on the saved candidate areas, estimates the reference candidate areas according to the urban planning indicators; considering the case of maximizing the distance from the fire source, finally obtains the ideal candidate areas. The present invention, based on the disaster area and the farthest safe distance of the helicopter, takes the disaster situation, economic situation, and terrain index as constraint conditions, traffic convenience as the preferred condition, and urban planning factors as reference conditions, to ensure the effectiveness of the site selection standard for the rescue material reserve depot and the actual declared site selection.

[0189] The present invention makes up for the gap in the site selection method for large-scale rescue material reserve depots based on earthquakes above magnitude 7, and guarantees the demand for the site selection of rescue material reserve depots. Through the solution of effective constraint conditions and preferred conditions, and then through the analysis of reference conditions, the ideal site selection area for the rescue material reserve depot can be obtained, which can effectively save the manpower, material resources, and time spent on on-site surveys, and provide a method for the currently vacant site selection problem of rescue material reserve depots, helping to quickly promote the site selection work of emergency management personnel for building depots.

[0190] In this specification, the various embodiments are described in a progressive manner. Each embodiment focuses on the differences from other embodiments. For the same or similar parts among the various embodiments, reference can be made to each other. For the system disclosed in the embodiments, since it corresponds to the method disclosed in the embodiments, the description is relatively simple. For the relevant parts, reference can be made to the description in the method section.

[0191] In this article, specific examples are used to elaborate on the principles and implementation manners of the present invention. The descriptions of the above embodiments are only used to help understand the method of the present invention and its core idea. At the same time, for those of ordinary skill in the art, based on the idea of the present invention, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation on the present invention.

Claims

1. A method for selecting the location of a disaster relief supply depot based on earthquakes of magnitude 7 or above, characterized in that, it includes: According to the earthquake literature of the research area, determine the earthquake risk areas in the research area with earthquake magnitudes of 7 or above; According to the earthquake risk areas and the earthquake intensity attenuation model, determine the intensity areas in the research area with earthquake intensities of 8 or above; Use remote sensing building identification technology to identify the buildings in the intensity area to generate a building distribution map, and then use the eight-connected method to process the building distribution map to obtain the disaster area; Determine the farthest safe flight distance when helicopters transport disaster relief supplies and determine the candidate areas for the location of the disaster relief supply depot in the research area; Determine the candidate areas for the location of the disaster relief supply depot in the research area, specifically including: Determine the first area and / or the second area in the research area as the candidate areas for the location of the disaster relief supply depot; the first area is the area determined with the highway entrance as the center and a radius of 5 km; the second area is the area determined with the railway freight station as the center and a radius of 5 km; Screen the candidate areas according to the constraint conditions to determine the optimal location area; the constraint conditions at least include the constraint relationship between the farthest safe flight distance and the disaster area, as well as the natural disasters, terrain index, traffic conditions and urban planning of the research area.

2. A method for selecting the location of a disaster relief supply depot based on earthquakes of magnitude 7 or above according to claim 1, characterized in that, According to the earthquake literature of the research area, determine the earthquake risk areas in the research area with earthquake magnitudes of 7 or above, specifically including: According to the prediction results and statistical laws of earthquake scientific research in the research area, determine the earthquake risk candidate areas in the research area; Determine the maximum earthquake magnitude corresponding to each earthquake risk candidate area; Compare the maximum earthquake magnitude with the set threshold, and determine the earthquake risk candidate area corresponding to the maximum earthquake magnitude greater than or equal to the set threshold as the earthquake risk area with an earthquake magnitude of 7 or above; wherein, the set threshold is 7.

3. A method for selecting the location of a disaster relief supply depot based on earthquakes of magnitude 7 or above according to claim 1, characterized in that, The earthquake risk area is an elliptical area; According to the earthquake risk area and the earthquake intensity attenuation model, determine the intensity areas in the research area with earthquake intensities of 8 or above, specifically including: According to the earthquake risk area and the earthquake intensity attenuation model, determine the elliptical isoseismal line when the earthquake intensity is 8; After adding the elliptical isoseismal line to the earthquake risk area and then performing elliptical delineation, obtain the intensity area with an earthquake intensity of 8 or above.

4. A method for selecting the location of a disaster relief supply depot based on earthquakes of magnitude 7 or above according to claim 3, characterized in that, The expression of the earthquake intensity attenuation model is: ； ； Among them, I is the seismic intensity, I = 8; Ms is the earthquake magnitude; R a and R b are the major semi - axis length and minor semi - axis length of the elliptical isoseismal line when the seismic intensity is 8 degrees respectively.

5. A method for selecting the location of a disaster relief supply depot based on earthquakes of magnitude 7 or above according to claim 1, characterized in that, The determination of the maximum safe flight distance when a helicopter transports disaster relief supplies specifically includes: According to the formula D max = ( TSR - FLR ) / 2 to determine the maximum safe flight distance when a helicopter transports disaster relief supplies; Among them, the constraint conditions for the safe flight of the helicopter are: v = MCS ； q = MP ； In the formula, D max is the farthest safe flight distance, TSR is the total standard flight range, FLR is the fuel quantity limit flight range, MCS is the maximum cruising speed, with the unit of km / h, MP is the maximum load, v is the speed of the helicopter during safe flight, q is the load of the helicopter during safe flight.

6. A method for selecting a site for a disaster relief supply reserve based on an earthquake of magnitude 7 or above as claimed in claim 1, characterized in that the natural disaster at least includes flood; the constraint condition further includes the statistical result of per capita GDP; The candidate site selection areas are screened according to the constraint conditions to determine the optimal site selection area, which specifically includes: Calculating a plurality of first target candidate areas based on the maximum safe flight distance and the disaster-stricken area, and determining the first target candidate areas that meet the candidate site selection areas as the second target candidate areas; Screening the second target candidate areas according to the flood areas and the intensity areas above magnitude 8 in the research area, and determining the second target candidate areas that meet the non-flood areas and non-intensity areas above magnitude 8 in the research area as the third target candidate areas; Screening the third target candidate areas according to the statistical result of per capita GDP in the research area in the past two years, and determining the third target candidate areas that meet the statistical result of per capita GDP in the research area in the past two years as the fourth target candidate areas; Screening the fourth target candidate areas according to the terrain index of the research area, and determining the fourth target candidate areas that meet the terrain index constraint of the research area as the fifth target candidate areas; Screening the fifth target candidate areas according to the traffic conditions of the research area, and determining the fifth target candidate areas that meet the traffic condition constraint of the research area as the sixth target candidate areas; Screening the sixth target candidate areas according to the urban planning of the research area, and determining the sixth target candidate areas that meet the urban planning constraint of the research area as the optimal site selection area.

7. A system for selecting a site for a disaster relief supply reserve based on an earthquake of magnitude 7 or above, characterized in that it includes: An earthquake risk area determination module, configured to determine an earthquake risk area with a seismic magnitude of 7 or above in the research area according to the earthquake literature materials of the research area; An intensity area determination module, configured to determine an intensity area with an earthquake intensity of 8 or above in the research area according to the earthquake risk area and the earthquake intensity attenuation model; A disaster-stricken area generation module, configured to identify the buildings in the intensity area by using the remote sensing building identification technology to generate a building distribution map, and then process the building distribution map in an eight-connected manner to obtain the disaster-stricken area; A maximum safe flight distance and candidate site selection area determination module, configured to determine the maximum safe flight distance when a helicopter transports disaster relief supplies and determine the candidate site selection area for the disaster relief supply reserve in the research area; Determining the candidate site selection area for the disaster relief supply reserve in the research area specifically includes: Identify the first area and / or the second area within the research area as the candidate areas for the location of the disaster relief supply depot; the first area is the area determined with the highway entrance as the center and a radius of 5 km; the second area is the area determined with the railway freight station as the center and a radius of 5 km; The optimal location area determination module is used to screen the candidate location areas according to the constraint conditions to determine the optimal location area; the constraint conditions at least include the constraint relationship between the farthest safe flight distance and the disaster area, as well as the natural disasters, terrain index, traffic conditions and urban planning of the research area.

8. A method for selecting a location for a disaster relief supply depot based on an earthquake of magnitude 7 or above according to claim 7, characterized in that The earthquake risk area determination module specifically includes: The earthquake risk candidate area unit is used to determine the earthquake risk candidate areas within the research area according to the prediction results and statistical laws of earthquake scientific research in the research area; The maximum earthquake magnitude determination unit is used to determine the maximum earthquake magnitude corresponding to each earthquake risk candidate area; The earthquake risk area determination unit is used to compare the maximum earthquake magnitude with the set threshold, and determine the earthquake risk candidate area corresponding to the maximum earthquake magnitude greater than or equal to the set threshold as the earthquake risk area with an earthquake magnitude of 7 or above; wherein, the set threshold is 7.

9. A method for selecting a location for a disaster relief supply depot based on an earthquake of magnitude 7 or above according to claim 7, characterized in that The natural disasters at least include floods; the constraint conditions also include the per capita GDP statistical results; the optimal location area determination module specifically includes: The second target candidate area determination unit is used to calculate a plurality of first target candidate areas according to the farthest safe flight distance and the disaster area, and determine the first target candidate areas that meet the candidate location areas as the second target candidate areas; The third target candidate area determination unit is used to screen the second target candidate areas according to the flood areas and intensity areas above magnitude 8 within the research area, and determine the second target candidate areas that meet the non-flood areas and non-intensity areas above magnitude 8 within the research area as the third target candidate areas; The fourth target candidate area determination unit is used to screen the third target candidate areas according to the per capita GDP statistical results of the research area in the past two years, and determine the third target candidate areas that meet the per capita GDP statistical results of the research area in the past two years as the fourth target candidate areas; The fifth target candidate area determination unit is used to screen the fourth target candidate areas according to the terrain index of the research area, and determine the fourth target candidate areas that meet the terrain index constraints of the research area as the fifth target candidate areas; The sixth target candidate area determination unit is used to screen the fifth target candidate areas according to the traffic conditions of the research area, and determine the fifth target candidate areas that meet the traffic condition constraints of the research area as the sixth target candidate areas; The optimal site selection area determination unit is configured to screen the sixth target candidate area according to the urban plan of the research area, and determine the sixth target candidate area that meets the urban plan constraints of the research area as the optimal site selection area.

Images