Earthquake prevention and disaster reduction thematic map drawing method and system based on disaster-bearing body vulnerability degree analysis and prediction technology

Abstract

The invention discloses an earthquake prevention and disaster reduction thematic map drawing method and system based on a disaster-bearing body vulnerability degree analysis and prediction technology, and relates to the technical field of earthquake fortification drawing. Comprising the steps of calculating an elliptical region range influenced by an earthquake by taking an epicenter position as a central point; calculating an earthquake influence field, and drawing the earthquake influence field on a map; analyzing administrative regions influenced by earthquakes and earthquake intensity conditions covered by each administrative region by utilizing spatial analysis of a GIS (Geographic Information System); analyzing loss conditions of the disaster-bearing body under different intensities to form a vulnerable matrix; predicting the vulnerability degree of each affected administrative region disaster-bearing body; a program is written, and the geographic coordinates of the influence field are drawn on canvas; calculating and acquiring a range of a corresponding basic geographic information base map in the map component; and drawing map symbols of various earthquake-proof and disaster-reduction thematic maps to form various earthquake-proof and disaster-reduction thematic maps, and performing drawing output according to set output parameters. According to the invention, earthquake-proof and disaster-reduction mapping can be completed in a convenient manner, and the efficiency of earthquake-proof and disaster-reduction mapping is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of earthquake-resistant fortification mapping, and in particular to a method and system for drawing earthquake-resistant and disaster-reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies. Background Art

[0002] Earthquake disasters are frequent. Their suddenness and destructiveness have caused huge losses of life and property to the people and seriously hindered economic development and social stability.

[0003] The model is used to calculate the impact range of the earthquake field and the earthquake intensity in different areas. GIS spatial analysis technology is used to obtain the overall situation of earthquake damage in each affected administrative area. An algorithm based on five damage levels is constructed, combined with regional historical earthquake data. After a series of analyses and evaluations, the vulnerability matrix of the disaster-prone body and the degree of damage in the current earthquake are obtained. Based on this technology, the mapping and SVG technologies are fully utilized to draw the symbols and analysis and prediction results related to seismic fortification onto the required map sheets, generate thematic maps such as population distribution maps, traffic distribution maps, school distribution maps, and hazard source distribution maps in the earthquake zone, and can output high-definition mapping results to provide auxiliary decision-making support for managers.

[0004] In the field of earthquake-resistant fortification, mapping is mostly done manually. GIS software, such as ArcGIS, can be used to assist mapping, but the degree of automation is still not high. Moreover, assessment and mapping can only be completed after an earthquake occurs through traditional means such as on-site assessment. There is no technical means for prediction. In order to assist decision makers in pre-perceiving the possible consequences of an earthquake, it is particularly important to provide a mapping method based on historical data analysis and prediction technology.

[0005] Therefore, it is an urgent problem for technical personnel in this field to propose a method and system for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies to solve the difficulties existing in the prior art. Summary of the invention

[0006] In view of this, the present invention provides a method and system for drawing earthquake prevention and disaster reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, which can improve the timeliness of preliminary assessment when an earthquake occurs and the efficiency of earthquake prevention and disaster reduction mapping.

[0007] In order to achieve the above object, the present invention adopts the following technical solution:

[0008] A method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, comprising:

[0009] S101, taking the epicenter as the center point, calculating the range of the elliptical area affected by the earthquake;

[0010] S102, calculate the earthquake impact field with each 0.1 intensity as a unit, and draw it on the map using a program; use GIS spatial analysis to analyze the administrative areas affected by the earthquake and the earthquake intensity covered by each administrative area;

[0011] S103. Based on earthquake zoning and historical earthquake data, analyze the loss of disaster-prone bodies under different intensities and form a vulnerability matrix;

[0012] S104. Based on the earthquake intensity and vulnerability matrix of the administrative area, predict the vulnerability of the disaster-bearing bodies in each affected administrative area;

[0013] S105, write a program to draw the geographic coordinates of the impact field on the canvas according to the scale of the map with the help of SVG technology; according to the map sheet, use the program to calculate and obtain the range of the corresponding basic geographic information base map in the map component;

[0014] S106. According to the relevant standards for earthquake prevention and disaster reduction, use the program to draw mapping components, evaluate the impact of earthquakes on various topics based on the vulnerability of disaster-bearing bodies in administrative areas and GIS spatial analysis, draw various types of earthquake prevention and disaster reduction thematic map symbols, form various types of earthquake prevention and disaster reduction thematic maps, and perform mapping output according to the set output parameters.

[0015] In the above method, optionally, in S101, calculating the range of the elliptical area affected by the earthquake specifically includes:

[0016] Taking the epicenter as the center point, calculate the length of the affected longitude direction, i.e. the major semi-axis;

[0017] The formula for calculating the semi-major axis is:

[0018] R ai =e (i-A-B*M) / C -R0

[0019] Among them, R a i is the semi-major axis value of the earthquake influence field when intensity = i; i is intensity; M is magnitude; A, B, C are constant parameters for calculating the semi-major axis of the earthquake influence field, and R0 is the correction parameter for the semi-major axis of the earthquake influence field;

[0020] Taking the epicenter as the center point, calculate the length of the affected latitude direction, i.e. the minor semi-axis;

[0021] The formula for calculating the minor semi-axis is:

[0022]

[0023] Among them, R biis the value of the short semi-axis of the earthquake influence field when the intensity is i; i is the intensity; M is the magnitude; A2, B2, C2 calculate the constant parameters of the long semi-axis of the earthquake influence field, R 20 Correction parameters of the minor semi-axis of the earthquake influence field.

[0024] In the above method, optionally, in S102, the earthquake impact field with each 0.1 intensity as a unit is calculated and drawn on a map using a program, specifically:

[0025] Write a program to loop with 0.1 intensity as a step, the minimum intensity is 6, and the maximum intensity is less than 11. In each loop, substitute the current intensity value i into the major and minor semi-axis calculation formula, and calculate the earthquake influence field R of the current loop with the constant parameters of the current earthquake magnitude and the location of the epicenter. ai , R bi value;

[0026] If the current intensity is calculated, R ai or R bi If it is not a positive number, the current maximum earthquake intensity is the intensity value of the previous cycle, and the cycle ends;

[0027] If R ai and R bi is a positive number, the program uses the map API to convert ai is the major semi-axis, R bi The elliptical earthquake influence field line R is the minor semi-axis i Draw on the map, i+0.1; repeat until the cycle ends, forming several elliptical earthquake influence field lines on the map.

[0028] In the above method, optionally, in S103, the vulnerability matrix of all disaster-sustaining bodies in the area under different intensities is calculated, and the formula is:

[0029]

[0030] Among them, k is the disaster-bearing body; R ki is the vulnerability matrix of disaster-bearing body k under intensity i; i is intensity;

[0031] The calculation formula of the vulnerability matrix of the disaster-bearing body k under intensity i is:

[0032]

[0033] The vulnerability matrix of the disaster-bearing body k under all intensities is:

[0034]

[0035] Each element of the vulnerability matrix of the disaster-bearing body k under intensity i is in fractional form, which is the ratio of quantity to area based on 100. The numerator is the relative quantity of the current damage level, and the denominator is the relative area of ​​the current damage level. The calculation formula is:

[0036]

[0037] m is different damage levels, there are five levels in total, corresponding to a, b, c, d, and e;

[0038] N ai ~N ei is the relative number of different damage levels of disaster-bearing body k under intensity i;

[0039] A ai ~A ei is the relative area of ​​disaster-bearing body k at different damage levels under intensity i;

[0040] C i is the constant parameter for calculating the seismic performance index of the disaster-bearing body k under intensity i, and takes values ​​of 2 to 1.65 for intensities 6 to 11 respectively;

[0041] E i is the constant parameter for calculating the seismic performance index of the disaster-bearing body k under intensity i, and the values ​​for intensities 6 to 11 are 0.007 to 0.002 respectively;

[0042] F s is the resistance of the disaster-bearing body k, is the comprehensive weight index of the disaster-bearing body k;

[0043] N mi is the relative number of disaster-bearing bodies at m damage level under intensity i;

[0044] N mi0 is the relative number of disaster-bearing bodies at damage level m calculated based on historical earthquake data;

[0045] A mi is the relative area of ​​the disaster-bearing body at m damage level under intensity i;

[0046] A mi0 is the relative area of ​​the disaster-bearing body at damage level m calculated based on historical earthquake data.

[0047] The above method is optional. The weighted sum of the vulnerability matrices of all hazard-prone bodies in a single administrative area at different intensities is as follows:

[0048]

[0049] Among them, F jis the weighted sum of the vulnerability matrices of all hazard-prone objects in administrative area j at different intensities; k is the hazard-prone object in the administrative area.

[0050] In the above method, optionally, in S104, the comprehensive vulnerability analysis prediction result is used to calculate the weighted sum of the damage ratios under each affected intensity for each affected administrative area, and the formula is:

[0051]

[0052] k is the disaster-prone body within the administrative area; A ki is the degree of damage of the disaster-bearing body k under intensity i, and is the weighted value of the disaster-bearing body k under different damage levels; p i is the area of ​​the administrative region affected by the earthquake impact field of intensity i; P is the total area of ​​the administrative region;

[0053] The area p of the administrative region affected by the earthquake impact field of intensity i i The calculation formula is:

[0054] P i =π*R ai *R bi .

[0055] In the above method, optionally, in S105, according to the map sheet, a program is used to calculate and obtain the range of the corresponding basic geographic information base map in the map component, and the specific method is:

[0056] The elliptical area of ​​the impact field is obtained with the epicenter as the center point;

[0057] Determine the coordinates of the four corners of the map display range based on the geographic coordinate range of the outermost circle of the impact field;

[0058] Determine the tile level of the vector map according to the printing accuracy required for mapping;

[0059] The program controls the map to automatically zoom to the map level that meets the conditions, uses canvas conversion technology to convert the required geographic base map into a high-definition picture for mapping, and uses SVG technology to draw it to the bottom layer of mapping;

[0060] Restore the map to the range displayed before the program automatically took control.

[0061] The above method is optional. In S106, the impact of earthquake on each topic is evaluated, and map symbols of various earthquake prevention and disaster reduction thematic maps are drawn to form various earthquake prevention and disaster reduction thematic maps, which specifically includes:

[0062] According to the vulnerability analysis and prediction results, combined with the population and economic kilometer grid data of the affected administrative areas and basic geographic information data, the casualties and direct economic losses are calculated, and the affected hazard sources and protection targets are analyzed;

[0063] Based on this information and data, and in accordance with the relevant specifications for earthquake disaster reduction mapping, a program is used to draw mapping components, including map names, scales, projections, and compasses;

[0064] The program draws various earthquake disaster reduction map symbols, including faults, fault zones, and historical earthquake examples;

[0065] Generate thematic maps such as population distribution map, economic distribution map, traffic distribution map, school distribution map, hospital distribution map, etc. in the earthquake zone.

[0066] A system for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, which implements any of the above-mentioned methods for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, comprising:

[0067] Earthquake impact field generation module: Based on earthquake information, calculate the coverage of different intensities and use the program to generate the earthquake impact field;

[0068] Affected administrative district analysis module: Based on the earthquake impact field, GIS spatial analysis is used to analyze the administrative districts affected by the earthquake and the earthquake intensity covered by each administrative district;

[0069] Single disaster-bearing body vulnerability matrix calculation module: calculates the vulnerability matrix of various types of disaster-bearing bodies under different intensities and different damage levels;

[0070] Single administrative district vulnerability assessment module: assesses the weighted sum of the vulnerability of all hazard-prone objects within the administrative district under each affected intensity;

[0071] Vulnerability comprehensive analysis and prediction module: Based on the vulnerability assessment results of each administrative district, combined with the population, economic kilometer grid data and basic geographic information data of the affected administrative districts, the impact of the earthquake is analyzed and predicted;

[0072] Map making module: Based on the analysis and prediction results of the impact of this earthquake, combined with relevant specifications for earthquake prevention and disaster reduction mapping, draw mapping components and various map symbols to form various earthquake prevention and disaster reduction thematic maps.

[0073] In the above system, optionally, in the vulnerability comprehensive analysis and prediction module: the program uses a loop method to loop according to the prefecture-level administrative districts, and in each loop, the data query capability of GIS is used to obtain the township-level administrative districts included in the current prefecture-level administrative district, and a second loop is performed on this;

[0074] In the secondary cycle, combined with the analysis results of the affected administrative area analysis module, it is determined whether the secondary cycle object is the administrative area affected by the earthquake. If so, the weighted sum of the vulnerability matrices of all disaster-bearing bodies at different intensities F j , combined with the current administrative area population, economic kilometer grid data and basic geographic information data, analyze and predict the impact of the earthquake on the current administrative area, including casualties, economic losses, affected gas stations, gas stations, hazardous chemical industrial product production units, hospitals, schools, stations and other hazardous sources and protection targets; if not, enter the next secondary cycle;

[0075] After completing the two-level cycle, we get a list of the impacts of the earthquake on the prefecture-level administrative regions and specific data on the impacts. Based on these data and using historical earthquake data, the program makes a comprehensive assessment and generates the results of the impact of the earthquake, including a statistical table of casualties and economic losses in the affected prefecture-level administrative regions, and identifies the affected hazards and protection targets on the map.

[0076] It can be seen from the above technical solution that, compared with the prior art, the present invention provides a method and system for drawing earthquake disaster prevention and reduction thematic maps based on disaster-bearing body vulnerability analysis and prediction technology, which has the following beneficial effects: the present invention can solve the problems of existing mapping methods that are highly dependent on evaluation results and excessive manual drawing in the mapping process while meeting the relevant specifications for earthquake disaster prevention and reduction mapping, thereby improving the timeliness, drawing efficiency and intelligence level of earthquake disaster prevention and reduction thematic maps. BRIEF DESCRIPTION OF THE DRAWINGS

[0077] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0078] Figure 1 A flowchart of a method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies provided by the present invention;

[0079] Figure 2 A flowchart of the intensity field analysis program provided by the present invention;

[0080] Figure 3 A flowchart of the analysis and prediction process provided by the present invention;

[0081] Figure 4 A flowchart of a program for obtaining a basic geographic information base map provided by the present invention;

[0082] Figure 5 A structural block diagram of a system for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies provided by the present invention;

[0083] Figure 6 The semi-long and semi-short axis values ​​of various intensities when an earthquake with a magnitude of 8.0 occurred in XX City provided by the present invention;

[0084] Figure 7 The earthquake impact field provided by the present invention when an earthquake with a magnitude of 8.0 occurs in the city of XX;

[0085] Figure 8 The vulnerability matrix of different disaster-bearing bodies under the intensity range when an earthquake with a magnitude of 8.0 occurs in XX City provided by the present invention;

[0086] Fig. 9 The thematic base map data when an earthquake with a magnitude of 8.0 occurred in XX city provided by the present invention;

[0087] Fig.10 , Fig.11 , Fig.12 The present invention provides various earthquake prevention and disaster reduction thematic maps when an earthquake with a magnitude of 8.0 occurs in XX City. DETAILED DESCRIPTION

[0088] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0089] In this application, relational terms such as first and second, etc. are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. The terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of more restrictions, the elements defined by the sentence "comprise one..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0090] Reference Figure 1 As shown, the present invention discloses a method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, comprising:

[0091] S101. Taking the epicenter as the center point, calculate the range of the elliptical area affected by the earthquake, i.e. the length of the major semi-axis and the length of the minor semi-axis;

[0092] S102, calculate the earthquake impact field with each 0.1 intensity as a unit, and draw it on the map using a program; use GIS spatial analysis to analyze the administrative areas affected by the earthquake and the earthquake intensity covered by each administrative area;

[0093] S103. Based on earthquake zoning and historical earthquake data, analyze the loss of disaster-prone bodies under different intensities and form a vulnerability matrix;

[0094] S104. Based on the earthquake intensity situation and vulnerability matrix of the administrative area, predict the vulnerability of the disaster-bearing bodies in each affected administrative area;

[0095] S105, write a program to draw the geographic coordinates of the impact field on the canvas according to the scale of the map with the help of SVG technology; according to the map sheet, use the program to calculate and obtain the range of the corresponding basic geographic information base map in the map component;

[0096] S106. According to the relevant standards for earthquake prevention and disaster reduction, use the program to draw mapping components, evaluate the impact of earthquakes on various topics (population, economy, etc.) based on the vulnerability of the disaster-bearing bodies in the administrative area and GIS spatial analysis, draw various types of earthquake prevention and disaster reduction thematic map symbols, form various types of earthquake prevention and disaster reduction thematic maps, and perform mapping output according to the set output parameters.

[0097] Furthermore, in S101, the range of the elliptical area affected by the earthquake is calculated, specifically including:

[0098] Taking the epicenter as the center point, calculate the length of the affected longitude direction, i.e. the major semi-axis;

[0099] The calculation formula for the major semi-axis is:

[0100] R ai =e (i-A-B*M) / C -R0

[0101] Among them, R ai is the semi-long axis value of the earthquake influence field when intensity = i; i is intensity; M is magnitude; A, B, C are constant parameters for calculating the semi-long axis of the earthquake influence field, and R0 is the correction parameter for the semi-long axis of the earthquake influence field; different regions have different values. According to relevant standards and specifications, there are two division methods, the Chinese four-division method and the seven-division method. The parameters of each division method are different for each region, as shown in Table 1:

[0102] Table 1 Parameters of each region in the four-division method of China

[0103] Partition parameters A B C <![CDATA[R0]]> Eastern earthquake zone 5.7123 1.3626 -4.2903 25 Moderate to strong earthquake zone 5.8410 1.0710 -3.6570 15 Xinjiang 5.6018 1.4347 -4.4899 25 Qinghai-Tibet Region 6.458 1.2746 -4.4709 25

[0104] Taking the epicenter as the center point, calculate the length of the affected latitude direction, i.e. the minor semi-axis;

[0105] The formula for calculating the minor semi-axis is:

[0106]

[0107] Among them, R bi is the value of the short semi-axis of the earthquake influence field when the intensity is i; i is the intensity; M is the magnitude; A2, B2, C2 calculate the constant parameters of the long semi-axis of the earthquake influence field, R 20 Correction parameters of the minor semi-axis of the earthquake influence field.

[0108] Further, see Figure 2 As shown, in S102, the earthquake impact field with each 0.1 intensity as a unit is calculated and drawn on the map using a program, specifically:

[0109] Write a program to loop with 0.1 intensity as a step, the minimum intensity is 6, and the maximum intensity is less than 11. In each loop, substitute the current intensity value i into the major and minor semi-axis calculation formula, and calculate the earthquake influence field R of the current loop with the constant parameters of the current earthquake magnitude and the location of the epicenter. ai , R bi value;

[0110] If the current intensity is calculated, R ai or R bi If it is not a positive number, the current maximum earthquake intensity is the intensity value of the previous cycle, which means that the current maximum earthquake intensity is the intensity value of the previous cycle, that is, i-0.1, and the cycle ends;

[0111] If R ai and R bi is a positive number, the program uses the map API to convert ai is the major semi-axis, R bi The elliptical earthquake influence field line R is the minor semi-axis i Draw on the map, i+0.1; repeat until the cycle ends, forming several elliptical earthquake influence field lines on the map.

[0112] Further, see Figure 3 As shown, in S103, the vulnerability matrix of all disaster-bearing bodies in the area under different intensities is calculated, and the formula is:

[0113]

[0114] Among them, k is the disaster-bearing body, including multi-story masonry houses, multi-story reinforced concrete houses, high-rise buildings, factory buildings, and other buildings (including urban villages); Rki is the vulnerability matrix of the disaster-bearing body k under intensity i; i is the intensity, ranging from 6 to 11, with a step of 0.1;

[0115] The calculation formula of the vulnerability matrix of the disaster-bearing body k under intensity i is:

[0116]

[0117] The vulnerability matrix of the disaster-bearing body k under all intensities is:

[0118]

[0119] Each element of the vulnerability matrix of the disaster-bearing body k under intensity i is in fractional form, which is the ratio of quantity to area based on 100. The numerator is the relative quantity of the current damage level, and the denominator is the relative area of ​​the current damage level. The calculation formula is:

[0120]

[0121] m is different levels of damage, including basically intact, slightly damaged, moderately damaged, severely damaged, and destroyed, corresponding to a, b, c, d, and e respectively;

[0122] N ai ~N ei is the relative number of different damage levels of disaster-bearing body k under intensity i;

[0123] A ai ~A ei is the relative area of ​​disaster-bearing body k at different damage levels under intensity i;

[0124] C i is the constant parameter for calculating the seismic performance index of the disaster-bearing body k under intensity i, and takes values ​​of 2 to 1.65 for intensities 6 to 11 respectively;

[0125] E i is the constant parameter for calculating the seismic performance index of the disaster-bearing body k under intensity i, and the values ​​for intensities 6 to 11 are 0.007 to 0.002 respectively;

[0126] F s is the resistance of the disaster-bearing body k, is the comprehensive weight index of the disaster-bearing body k, for example, the resistance of multi-story masonry houses is the comprehensive weight index of the resistance of all multi-story masonry houses within the statistical area, which is affected by factors such as the construction year;

[0127] N mi is the relative number of disaster-bearing bodies at m damage level under intensity i;

[0128] N mi0is the relative number of disaster-bearing bodies at damage level m calculated based on historical earthquake data;

[0129] A mi is the relative area of ​​the disaster-bearing body at m damage level under intensity i;

[0130] A mi0 is the relative area of ​​the disaster-bearing body at damage level m calculated based on historical earthquake data.

[0131] Furthermore, the weighted sum of the vulnerability matrices of all disaster-prone bodies in a single administrative region at different intensities is given by:

[0132]

[0133] Among them, F j is the weighted sum of the vulnerability matrices of all hazard-prone objects in administrative area j at different intensities; k is the hazard-prone object in the administrative area.

[0134] Furthermore, in S104, the comprehensive vulnerability analysis prediction result is used to calculate the weighted sum of the damage ratios under each affected intensity for each affected administrative area, and the formula is:

[0135]

[0136] k is the disaster-prone body within the administrative area; A ki is the degree of damage of the disaster-bearing body k under intensity i, and is the weighted value of the disaster-bearing body k under different damage levels; p i is the area of ​​the administrative region affected by the earthquake impact field of intensity i; P is the total area of ​​the administrative region;

[0137] The area p of the administrative region affected by the earthquake impact field of intensity i i The calculation formula is:

[0138] P i =π*R ai *R bi .

[0139] Further, see Figure 4 As shown, in S105, according to the map sheet, a program is used to calculate and obtain the range of the corresponding basic geographic information base map in the map component, and the specific method is:

[0140] The elliptical area of ​​the impact field is obtained with the epicenter as the center point;

[0141] Determine the coordinates of the four corners of the map display range according to the geographic coordinate range of the outermost circle of the impact field (the earthquake impact field line with the largest range);

[0142] Determine the tile level of the vector map according to the printing accuracy required for cartography to ensure that the precision meets the standards for cartographic output;

[0143] The program controls the map to automatically zoom to the map level that meets the conditions, uses canvas conversion technology to convert the required geographic base map into a high-definition picture for mapping, and uses SVG technology to draw it to the bottom layer of mapping (below the vector layer of the impact field);

[0144] Restore the map to the range displayed before the program automatically took control.

[0145] Furthermore, in S106, the impact of the earthquake on each topic is evaluated, and various earthquake disaster reduction thematic map symbols are drawn to form various earthquake disaster reduction thematic maps, including:

[0146] According to the vulnerability analysis and prediction results, combined with the population and economic kilometer grid data of the affected administrative areas and basic geographic information data, the casualties and direct economic losses are calculated, and the affected hazard sources and protection targets are analyzed;

[0147] Based on this information and data, and in accordance with the relevant specifications for earthquake disaster reduction mapping, a program is used to draw mapping components, including map names, scales, projections, and compasses;

[0148] The program draws various earthquake disaster reduction map symbols, including faults, fault zones, and historical earthquake examples;

[0149] Generate thematic maps such as population distribution map, economic distribution map, traffic distribution map, school distribution map, hospital distribution map, etc. in the earthquake zone.

[0150] and Figure 1 Corresponding to the method shown, the present invention also discloses a system for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of the disaster-bearing body for realizing the method. The structural block diagram is referred to Figure 5 As shown, including:

[0151] Earthquake impact field generation module: used to calculate the coverage area (long and short semi-axis lengths) of different intensities based on earthquake information (epicenter location, magnitude, focal depth), and generate earthquake impact fields using the program;

[0152] Affected administrative district analysis module: used to analyze the administrative districts (township level) affected by the earthquake and the earthquake intensity covered by each administrative district (each 0.1 intensity is a level, calculated one by one) based on the earthquake impact field using GIS spatial analysis;

[0153] Single disaster-bearing body vulnerability matrix calculation module: used to calculate the vulnerability matrix of various types of disaster-bearing bodies under different intensities and different damage levels;

[0154] Single administrative district vulnerability assessment module: used to assess the weighted sum of the vulnerability of all disaster-prone bodies within the administrative district (township level) under each affected intensity;

[0155] Vulnerability comprehensive analysis and prediction module: used to analyze and predict the impact of the earthquake based on the vulnerability assessment results of each administrative area, combined with the population, economic kilometer grid data and basic geographic information data of the affected administrative areas;

[0156] Map drawing module: It is used to draw mapping components and various map symbols according to the analysis and prediction results of the impact of this earthquake, combined with relevant specifications for earthquake prevention and disaster reduction mapping, to form various earthquake prevention and disaster reduction thematic maps.

[0157] The single disaster-bearing body vulnerability matrix calculation module specifically includes:

[0158] The losses caused by earthquake disasters are related to the comprehensive vulnerability of all disaster-bearing bodies within the disaster area. The comprehensive vulnerability is obtained by weighting the vulnerability of a single disaster-bearing body. The vulnerability of a single disaster-bearing body is affected by the building damage level. Different intensities have different damage levels to buildings. Under the same intensity, the damage level is different for different building structures. According to earthquake-related standards and specifications, the damage level is divided into five levels: basically intact, slightly damaged, moderately damaged, severely damaged, and destroyed; the building structure is divided into multi-story masonry houses, multi-story reinforced concrete houses, high-rise buildings, factory buildings, and other (including urban villages) buildings. The vulnerability matrix is ​​calculated by the formula; each element of the vulnerability matrix is ​​in fractional form, which is the ratio of quantity to area based on 100, with the numerator being the relative quantity of the current damage level (based on 100), and the denominator being the relative area of ​​the current damage level (based on 100); the vulnerability matrix of disaster-bearing body k under intensity i is:

[0159]

[0160] The vulnerability matrix of the disaster-bearing body k under all intensities is:

[0161]

[0162] Among them, abcde are the damage levels, basically intact a, slightly damaged b, moderately damaged c, severely damaged d, destroyed e;

[0163] N ai ~N ai is the relative number of different damage levels of disaster-bearing body k under intensity i (relative number based on 100);

[0164] A ai ~A eiIt is the relative area of ​​disaster-bearing body k at different damage levels under intensity i (relative area based on 100).

[0165] Single administrative district vulnerability assessment module: used to assess the weighted sum of the vulnerability of all disaster-prone bodies within the administrative district (township level) under each affected intensity;

[0166] The single administrative district vulnerability assessment module specifically includes:

[0167] Calculate the weighted sum of the vulnerability matrix of all disaster-prone bodies in the administrative area at different intensities. The formula is:

[0168]

[0169] Among them, F j is the weighted sum of the vulnerability matrices of all hazard-bearing bodies in administrative area j at different intensities, and k is the hazard-bearing body in the administrative area.

[0170] R k is the vulnerability matrix of a single hazard-bearing body, and the formula is:

[0171]

[0172] R k That is, the vulnerability matrix of the disaster-bearing body k under all intensities, R ki That is, the vulnerability matrix of disaster-bearing body k under intensity i;

[0173] i is the intensity, and its value range is 6 to 11, with a step of 0.1 until the maximum intensity is reached.

[0174] Vulnerability comprehensive analysis and prediction module: used to analyze and predict the impact of the earthquake based on the vulnerability assessment results of each administrative area, combined with the population, economic kilometer grid data and basic geographic information data of the affected administrative areas;

[0175] Furthermore, in the vulnerability comprehensive analysis and prediction module: the program uses a loop method to loop according to the prefecture-level administrative districts. In each loop, the data query capability of GIS is used to obtain the township-level administrative districts included in the current prefecture-level administrative district, and a second loop is performed on this;

[0176] In the secondary cycle, combined with the analysis results of the affected administrative area analysis module, it is determined whether the secondary cycle object is the administrative area affected by the earthquake. If so, the weighted sum of the vulnerability matrices of all disaster-bearing bodies at different intensities F j, combined with the current administrative area population, economic kilometer grid data and basic geographic information data, analyze and predict the impact of the earthquake on the current administrative area, including casualties, economic losses, affected gas stations, gas stations, hazardous chemical industrial product production units, hospitals, schools, stations and other hazardous sources and protection targets; if not, enter the next secondary cycle;

[0177] After completing the two-level cycle, we get a list of the impacts of the earthquake on the prefecture-level administrative regions and specific data on the impacts. Based on these data and using historical earthquake data, the program makes a comprehensive assessment and generates the results of the impact of the earthquake, including a statistical table of casualties and economic losses in the affected prefecture-level administrative regions, and identifies the affected hazards and protection targets on the map.

[0178] In a specific embodiment, taking the drawing of earthquake disaster prevention and reduction thematic map of XX city as an example:

[0179] When an earthquake occurs in XX City, the main earthquake information is obtained based on the earthquake information pushed by the earthquake early warning network, including the epicenter location, magnitude, and time of occurrence;

[0180] According to the historical data of XX City at that time and the attenuation model corrected by experience, the area affected by the earthquake was calculated and the earthquake impact field was drawn. After the parameters of the attenuation model were corrected using the existing historical earthquake data of XX City, the long and short semi-axis values ​​of the area affected by the earthquake were calculated; Figure 6 As shown in the figure, when an earthquake with a magnitude of 8.0 occurs in XX City, the long and short semi-axis values ​​of each intensity can be calculated by using the attenuation model of eastern China and correcting the parameters with the historical earthquake cases of XX City;

[0181] Reference Figure 7 As shown, the program draws the earthquake impact field on the map based on the algorithm and the long and short semi-axis values ​​of each intensity field line (in units of 0.1 intensity, Figure 7 In order to more clearly present the intensity field on the map, the whole intensity field is displayed); because the highest intensity level calculated by the algorithm based on the comprehensive seismic data of this earthquake is X, the map does not show the XI-level earthquake impact field, and there will be no XI-level intensity affected area for this earthquake. The highest intensity level of each earthquake is calculated from the seismic data, and it is different for each earthquake;

[0182] Reference Figure 8 As shown in the figure, using the algorithm and matrix calculation formula, combined with the relevant information of this earthquake, the vulnerability matrix of different disaster-bearing bodies (including multi-story masonry houses, multi-story reinforced concrete houses, high-rise buildings, factory buildings, and other (including urban villages) buildings) in this earthquake is calculated under the intensity range (VI-X) affected by this earthquake;

[0183] Get the thematic map base map data, refer to Fig. 9 As shown, with the epicenter as the center of the circle and the earthquake impact field analyzed above as the radius of the major and minor axes, the data range of the basic geographic base map is obtained, and the final area delineation is carried out in combination with the map sheet information selected by the user in the mapping, and the data range of the mapping base map calculated by the algorithm;

[0184] Obtain thematic map information, analyze and predict the results of vulnerability, combine the population and economic kilometer grid data of the affected administrative areas with basic geographic information data, calculate the casualties and direct economic losses, and analyze the affected hazard sources and protection targets;

[0185] Draw symbols and thematic map output, according to the relevant specifications of earthquake disaster reduction mapping, use the program to draw mapping components, draw various earthquake disaster reduction map symbols and finally draw the map, refer to Fig.10 , Fig.11 , Fig.12 shown.

[0186] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can refer to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system or system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can refer to the partial description of the method embodiment. The system and system embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without creative work.

[0187] The above description of the disclosed embodiments enables one skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but rather to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, characterized in that: include: S101, taking the epicenter as the center point, calculating the range of the elliptical area affected by the earthquake; S102, calculate the earthquake impact field with each 0.1 intensity as a unit, and draw it on the map using a program; use GIS spatial analysis to analyze the administrative areas affected by the earthquake and the earthquake intensity covered by each administrative area; S103. Based on earthquake zoning and historical earthquake data, analyze the loss of disaster-prone bodies under different intensities and form a vulnerability matrix; S104. Based on the earthquake intensity and vulnerability matrix of the administrative area, predict the vulnerability of the disaster-bearing bodies in each affected administrative area; S105, write a program to draw the geographic coordinates of the impact field on the canvas according to the scale of the map with the help of SVG technology; according to the map sheet, use the program to calculate and obtain the range of the corresponding basic geographic information base map in the map component; S106. According to the relevant standards for earthquake prevention and disaster reduction, use the program to draw mapping components, evaluate the impact of earthquakes on various topics based on the vulnerability of disaster-bearing bodies in administrative areas and GIS spatial analysis, draw various types of earthquake prevention and disaster reduction thematic map symbols, form various types of earthquake prevention and disaster reduction thematic maps, and perform mapping output according to the set output parameters.

2. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 1 is characterized in that: In S101, the range of the elliptical area affected by the earthquake is calculated, including: Taking the epicenter as the center point, calculate the length of the affected longitude direction, i.e. the major semi-axis; The formula for calculating the semi-major axis is: R ai s.e (i-A-B*M) / C -R0 Among them, R ai is the semi-long axis value of the earthquake influence field when intensity = i; i is intensity; M is magnitude; A, B, C are constant parameters for calculating the semi-long axis of the earthquake influence field, and R0 is the correction parameter for the semi-long axis of the earthquake influence field; Taking the epicenter as the center point, calculate the length of the affected latitude direction, i.e. the minor semi-axis; The formula for calculating the minor semi-axis is: Among them, R bi is the value of the short semi-axis of the earthquake influence field when the intensity is i; i is the intensity; M is the magnitude; A2, B2, C2 calculate the constant parameters of the long semi-axis of the earthquake influence field, R 20 Correction parameters of the minor semi-axis of the earthquake influence field.

3. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 1 is characterized in that: In S102, the earthquake impact field with each 0.1 intensity as a unit is calculated and drawn on the map using a program, specifically: Write a program to loop with 0.1 intensity as a step, the minimum intensity is 6, and the maximum intensity is less than 11. In each loop, substitute the current intensity value i into the major and minor semi-axis calculation formula, and calculate the earthquake influence field R of the current loop with the constant parameters of the current earthquake magnitude and the location of the epicenter. ai , R bi value; If the current intensity is calculated, R ai or R bi If it is not a positive number, the current maximum earthquake intensity is the intensity value of the previous cycle, and the cycle ends; If R ai and R bi is a positive number, the program uses the map API to convert ai is the major semi-axis, R bi The elliptical earthquake influence field line R is the minor semi-axis i Draw on the map, i+0.1; repeat until the cycle ends, forming several elliptical earthquake influence field lines on the map.

4. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 1 is characterized in that: In S103, the vulnerability matrix of all disaster-affected bodies in the region under different intensities is calculated, and the formula is: Among them, k is the disaster-bearing body; R ki is the vulnerability matrix of disaster-bearing body k under intensity i; i is intensity; The calculation formula of the vulnerability matrix of the disaster-bearing body k under intensity i is: The vulnerability matrix of the disaster-bearing body k under all intensities is: Each element of the vulnerability matrix of the disaster-bearing body k under intensity i is in fractional form, which is the ratio of quantity to area based on 100. The numerator is the relative quantity of the current damage level, and the denominator is the relative area of ​​the current damage level. The calculation formula is: m is different damage levels, there are five levels in total, corresponding to a, b, c, d, and e; N ai ~N ei is the relative number of different damage levels of disaster-bearing body k under intensity i; A ai ~A ei is the relative area of ​​disaster-bearing body k at different damage levels under intensity i; C i is the constant parameter for calculating the seismic performance index of the disaster-bearing body k under intensity i, and takes values ​​of 2 to 1.65 for intensities 6 to 11 respectively; E i is the constant parameter for calculating the seismic performance index of the disaster-bearing body k under intensity i, and the values ​​for intensities 6 to 11 are 0.007 to 0.002 respectively; F s is the resistance of the disaster-bearing body k, is the comprehensive weight index of the disaster-bearing body k; N mi is the relative number of disaster-bearing bodies at m damage level under intensity i; N mi0 is the relative number of disaster-bearing bodies at damage level m calculated based on historical earthquake data; A mi is the relative area of ​​the disaster-bearing body at m damage level under intensity i; A mi0 is the relative area of ​​the disaster-bearing body at damage level m calculated based on historical earthquake data.

5. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 4 is characterized in that: The weighted sum of the vulnerability matrices of all hazard-bearing bodies in a single administrative region at different intensities is as follows: Among them, F j is the weighted sum of the vulnerability matrices of all hazard-prone objects in administrative area j at different intensities; k is the hazard-prone object in the administrative area.

6. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 1 is characterized in that: In S104, the comprehensive vulnerability analysis prediction results are calculated for each affected administrative area, and the weighted sum of the damage ratios under each affected intensity is calculated using the formula: k is the disaster-prone body within the administrative area; A ki is the degree of damage of the disaster-bearing body k under intensity i, and is the weighted value of the disaster-bearing body k under different damage levels; p i is the area of ​​the administrative region affected by the earthquake impact field of intensity i; P is the total area of ​​the administrative region; The area p of the administrative region affected by the earthquake impact field of intensity i i The calculation formula is: P i =π*R ai *R bi 。 7. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 1 is characterized in that: In S105, according to the map sheet, a program is used to calculate and obtain the range of the corresponding basic geographic information base map in the map component, and the specific method is as follows: The elliptical area of ​​the impact field is obtained with the epicenter as the center point; Determine the coordinates of the four corners of the map display range based on the geographic coordinate range of the outermost circle of the impact field; Determine the tile level of the vector map according to the printing accuracy required for mapping; The program controls the map to automatically zoom to the map level that meets the conditions, uses canvas conversion technology to convert the required geographic base map into a high-definition picture for mapping, and uses SVG technology to draw it to the bottom layer of mapping; Restore the map to the range displayed before the program automatically took control.

8. The method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies according to claim 1 is characterized in that: In S106, the impact of earthquakes on various topics is evaluated, and map symbols of various earthquake disaster prevention and reduction thematic maps are drawn to form various earthquake disaster prevention and reduction thematic maps, including: According to the vulnerability analysis and prediction results, combined with the population and economic kilometer grid data of the affected administrative areas and basic geographic information data, the casualties and direct economic losses are calculated, and the affected hazard sources and protection target information are analyzed; Based on this information and data, and in accordance with the relevant specifications for earthquake disaster reduction mapping, a program is used to draw mapping components, including map names, scales, projections, and compasses; The program draws various earthquake disaster reduction map symbols, including faults, fault zones, and historical earthquake examples; Generate thematic maps such as population distribution map, economic distribution map, traffic distribution map, school distribution map, hospital distribution map, etc. in the earthquake zone.

9. A system for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies, executing a method for drawing earthquake disaster prevention and reduction thematic maps based on the vulnerability analysis and prediction technology of disaster-bearing bodies as claimed in any one of claims 1 to 8, comprising: Earthquake impact field generation module: Based on earthquake information, calculate the coverage of different intensities and use the program to generate the earthquake impact field; Affected administrative district analysis module: Based on the earthquake impact field, GIS spatial analysis is used to analyze the administrative districts affected by the earthquake and the earthquake intensity covered by each administrative district; Single disaster-bearing body vulnerability matrix calculation module: calculates the vulnerability matrix of various types of disaster-bearing bodies under different intensities and different damage levels; Single administrative district vulnerability assessment module: assesses the weighted sum of the vulnerability of all hazard-prone objects within the administrative district under each affected intensity; Vulnerability comprehensive analysis and prediction module: Based on the vulnerability assessment results of each administrative district, combined with the population, economic kilometer grid data and basic geographic information data of the affected administrative districts, the impact of the earthquake is analyzed and predicted; Map making module: Based on the analysis and prediction results of the impact of this earthquake, combined with relevant specifications for earthquake prevention and disaster reduction mapping, draw mapping components and various map symbols to form various earthquake prevention and disaster reduction thematic maps.

10. The earthquake disaster prevention and reduction thematic map drawing system based on the vulnerability analysis and prediction technology of the disaster-bearing body according to claim 9 is characterized in that: In the vulnerability comprehensive analysis and prediction module: the program uses a loop method to loop according to the prefecture-level administrative districts. In each loop, the data query capability of GIS is used to obtain the township-level administrative districts included in the current prefecture-level administrative district, and a second loop is performed on this. In the secondary cycle, combined with the analysis results of the affected administrative area analysis module, it is determined whether the secondary cycle object is the administrative area affected by the earthquake. If so, the weighted sum of the vulnerability matrices of all disaster-bearing bodies at different intensities F j , combined with the current administrative area population, economic kilometer grid data and basic geographic information data, analyze and predict the impact of the earthquake on the current administrative area, including casualties, economic losses, affected gas stations, gas stations, hazardous industrial product production units, hospitals, schools, stations, hazard sources and protection targets; if not, enter the next secondary cycle; After completing the two-level cycle, we get a list of the impacts of the earthquake on the prefecture-level administrative regions and specific data on the impacts. Based on these data and using historical earthquake data, the program makes a comprehensive assessment and generates the results of the impact of the earthquake, including a statistical table of casualties and economic losses in the affected prefecture-level administrative regions, and identifies the affected hazards and protection targets on the map.

Images