Geological disaster risk comprehensive assessment method, device, equipment and medium

Abstract

The invention relates to the technical field of geological disaster assessment, in particular to a geological disaster risk comprehensive assessment method, device and equipment and a medium, and the method comprises the steps: obtaining the landform data and rock-soil geological data of a to-be-assessed region, and the type and range of a proposed project; based on the landform data, performing landform unit division on the to-be-evaluated area to obtain DEM grids, the DEM grids being elevation data stored in a grid form, and each grid unit corresponding to an elevation value; the rock and soil geological data is matched with the DEM grid, and a lithology-elevation comprehensive map is generated; enabling the proposed project type and the proposed project range to correspond to the lithology-elevation comprehensive graph, and generating a risk value; and determining an engineering protection measure based on the risk value. The geological risk assessment method has the effect of improving the accuracy of geological risk assessment.

Description

Technical Field

[0001] The present application relates to the technical field of geological hazard assessment, and in particular to a method, device, equipment and medium for comprehensive geological hazard risk assessment. Background Art

[0002] Geological disasters refer to geological actions or phenomena caused by natural or human factors, which cause losses to human life and property and damage to the environment. The distribution and changes of geological disasters in time and space are both subject to the natural environment and related to human activities. They are often the result of the interaction between humans and nature.

[0003] Existing geological hazard assessment methods all use models based on historical data and empirical formulas, which cannot fully reflect the complexity and uncertainty of geological hazards. At the same time, the selection and setting of model parameters also have certain subjectivity and uncertainty, which in turn affects the accuracy of the assessment results. Summary of the Invention

[0004] In order to improve the accuracy of geological risk assessment, the present application provides a method, device, equipment and medium for comprehensive assessment of geological hazard risks.

[0005] In the first aspect, the present application provides a method for comprehensive assessment of geological hazard risk, which adopts the following technical solutions: A comprehensive geological hazard risk assessment method, comprising: Obtain topographic and geomorphological data, geotechnical data, and the type and scope of the proposed project in the area to be assessed; Based on the topographic data, the area to be assessed is divided into geomorphic units to obtain a DEM grid, wherein the DEM grid is elevation data stored in a grid format, and each grid cell corresponds to an elevation value; Matching the rock and soil geological data with the DEM grid to generate a lithology-elevation composite map; Matching the proposed project type and the proposed project scope with the lithology-elevation comprehensive map to generate a risk value; Engineering protection measures are determined based on the risk value.

[0006] By adopting the above technical solutions, since topographic and geomorphic data reflect the surface morphological characteristics, and rock and soil geological data reflect the geological structure and geotechnical properties of the region, the type and scope of the proposed project provide specific application scenarios for risk assessment. Therefore, integrating multi-source data such as topographic and geomorphic data, rock and soil geological data, and the type and scope of the proposed project can more comprehensively and accurately assess the risk of geological hazards and avoid the limitations of single data assessment. Based on the topographic and geomorphic data, geomorphic units are divided to obtain DEM grids. Continuous topographic and geomorphic data can be integrated into regular grid units, which is convenient for subsequent data processing. The DEM grid can also provide a more intuitive understanding of the terrain undulations and topographic changes in the area to be assessed. Matching the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map is to combine geotechnical properties with topographic elevation information, thereby comprehensively reflecting the distribution of geological characteristics of the area to be assessed at different elevations, which helps to identify areas with potential geological disaster risks. Specifically, the geological characteristics of the area to be assessed are refined, making hazard assessments more accurate and enabling the identification of potential small-scale geological hazard risk areas. Because different project types and scopes have different impacts on the geological environment, mapping the proposed project type and scope with the lithology-elevation map allows for the assessment of the potential risks of the proposed project under different geological conditions. The magnitude of the risk value reflects the likelihood and severity of a geological hazard. Determining engineering protection measures based on the risk value allows for targeted development of engineering protection measures based on the assessed risk value, improving protection effectiveness, reducing protection costs, and more effectively ensuring the safety of the proposed project.

[0007] In another possible implementation, dividing the area to be assessed into geomorphic units based on the topographic data includes: Determine the slope data of the area to be assessed from the topographic data; Dividing the slope data into multiple slope levels; The ratio of the area corresponding to each slope grade to the total area of ​​the DEM grid is calculated to obtain the spatial distribution ratio of each geomorphic unit; The spatial distribution ratio and the multiple slope levels are used as a basis for DEM grid classification to perform landform unit division.

[0008] By adopting the above technical solution, since slope data is one of the key characteristic parameters of topography, it has an important impact on the occurrence of geological disasters, surface runoff, soil erosion and other processes. Therefore, accurately determining the slope data of the area to be evaluated from the topography and geomorphology data is the basis for the subsequent division of geomorphic units and geological disaster risk assessment. Since different slope levels correspond to different geomorphic characteristics and geological disaster risk levels, the slope data is divided into multiple slope levels, and areas with different slope ranges can be classified to facilitate subsequent analysis. Among them, the division of slope levels can be determined according to the actual situation and relevant standards of the area to be evaluated. By calculating the ratio of the area corresponding to each slope level to the total area of ​​the DEM grid, the spatial distribution ratio of each geomorphic unit is obtained, and the spatial distribution of different slope levels in the area to be evaluated can be understood, that is, the area ratio of the area occupied by different slope levels. The spatial distribution ratio can provide more comprehensive information for subsequent geomorphic unit division and geological hazard risk assessment. Geomorphic unit division is an important part of the comprehensive assessment of geological hazard risk. By dividing areas with similar terrain features into the same geomorphic unit, the geological hazard risk characteristics of different geomorphic units can be better analyzed.

[0009] In another possible implementation, matching the rock and soil geological data with the DEM grid to generate a lithology-elevation integrated map includes: Determining lithology data based on the rock and soil geological data, wherein the lithology data includes lithology type, shear strength, and weathering degree; Superimposing the lithology type with the DEM grid to obtain a superimposed DEM grid, wherein each grid cell in the superimposed DEM grid corresponds to a lithology data; Calculating a stability index of each cell based on the shear strength, the degree of weathering, and the elevation weight of each grid cell; The DEM grid is rendered in a graded manner according to the stability index to generate a lithology-elevation integrated map.

[0010] Through the above technical solution, lithology type reflects the material composition and basic properties of the rock mass; shear strength measures the rock mass's ability to resist shear failure and has a significant impact on slope stability and other factors; and weathering degree reflects the physical and chemical changes in the rock mass under the influence of the natural environment. Rock mass with different degrees of weathering has significantly different engineering properties. Therefore, it is necessary to clearly define lithology data, including lithology type, shear strength, and weathering degree. The lithology type is overlaid with the DEM grid, so that each grid cell in the DEM grid corresponds to a lithology data point. Because the distribution of different lithologies may vary with terrain elevation, combining lithology data with terrain elevation information allows for a visual display of topographic and elevation distribution patterns. The stability index reflects the stability of the rock mass under current terrain and lithology conditions. By considering shear strength, weathering degree, and elevation, a more comprehensive assessment of the stability of the area under assessment can be achieved. Rendering can display areas with different stability indices using different colors or symbols, visually demonstrating the stability distribution of the area under assessment. This helps personnel quickly understand the stability status of the area under assessment.

[0011] In another possible implementation, the generating of the risk value by matching the proposed project type and the proposed project scope with the lithology-elevation comprehensive map includes: Superimposing the proposed project scope with the lithology-elevation comprehensive map to obtain a lithology-elevation comprehensive map of the project impact area; Dividing the proposed project into regions to obtain a plurality of proposed project sub-regions; Obtain the classification parameters and project attributes of the proposed project in the proposed project sub-area and calculate the environmental impact coefficient; Calculate based on the infrastructure coverage rate and the type of the proposed project to obtain a project impact coefficient; Obtain the stability index of the proposed project based on the lithology-elevation comprehensive map of the project impact area; The risk value is calculated based on the stability index of the proposed project, the project impact coefficient and the environmental impact coefficient.

[0012] By employing the above technical solution, the proposed project scope is overlaid with the lithology-elevation comprehensive map to clarify the specific lithology and elevation distribution of the area where the proposed project is located. The project impact coefficient and environmental impact coefficient are determined based on the proposed project type to quantify the impact of different project types on the surrounding environment. This facilitates a more comprehensive assessment of the risks posed by the proposed project. The proposed project stability index, derived from the lithology-elevation comprehensive map of the project's impact area, is used to assess the stability of the proposed project under specific lithology and elevation conditions. The stability index comprehensively considers the impact of lithology and elevation on project stability, providing a direct reflection of the project's safety status under natural environmental conditions. Furthermore, the risk value, through the proposed project stability index, project impact coefficient, and environmental impact coefficient, takes into account multiple factors, including the stability of the proposed project itself, the impact of the proposed project construction on the surrounding environment, and the impact of the natural environment on the proposed project construction. By calculating the risk value, the risk level of the proposed project can be more accurately assessed, avoiding the limitations of single-factor assessments.

[0013] In another possible implementation, determining engineering protection measures based on the risk value includes: Classifying the risk value into risk levels to obtain at least one risk level; formulating a protective measures mapping table based on the at least one risk level and the engineering protective measures; Extracting risk values ​​for subunits of the proposed project from a comprehensive lithology-elevation map of the project's affected area; Corresponding engineering protection measures are determined in the protection measures mapping table based on the risk value of the proposed engineering subunit.

[0014] By employing the above technical solution, risk values ​​are classified into risk levels to provide a clearer understanding of the range of different risk levels, facilitating the subsequent development of targeted protective measures. Different risk levels represent varying degrees of potential harm. This classification allows for a more intuitive assessment of the severity of the risk, thus providing a basis for developing protective measures. Risk levels are determined using specific standards or methods. Risk levels can be determined by reference to industry standards, historical data, or expert experience. A protective measure mapping table, based on at least one risk level and project protective measures, establishes a correspondence between risk levels and specific protective measures. Once the risk level is determined, the corresponding protective measures can be quickly and accurately identified, improving the efficiency of protective measure development. Extracting risk values ​​for proposed project subunits from a comprehensive lithology-elevation map of the project's impact area enables a more detailed risk assessment of the proposed project. Proposed projects typically consist of multiple subunits, each with varying lithology and elevation conditions, resulting in varying risk values. Extracting subunit risk values ​​provides a more accurate understanding of the risk profile of each subunit, providing a basis for developing targeted protective measures. By identifying the corresponding project protection measures in the protection measures mapping table based on the risk values ​​of the proposed project subunits, a specific protection plan can be developed for each proposed project subunit. By matching the risk values ​​of the proposed project subunits with the protection measures mapping table, appropriate protection measures can be quickly identified, ensuring their effectiveness and relevance. This fully considers the actual conditions of different proposed project subunits and improves the relevance of protection measures.

[0015] In another possible implementation, the step of rendering the DEM grid in a hierarchical manner according to the stability index to generate a lithology-elevation integrated map includes: Determine the lithology code corresponding to each sub-grid unit in the DEM grid by using the lithology type and the DEM grid, wherein the lithology code is the number of each lithology, and the sub-grid unit is a subset of the DEM grid, that is, a plurality of the sub-grid units constitute the DEM grid; Classifying the stability index into grades to obtain a stability grade; Mapping the elevation value to the R channel, mapping the lithology code to the G channel, and mapping the stability grade to the B channel; Merge the R channel, G channel, and B channel to generate an RGB image; The RGB image is determined as a lithology-elevation comprehensive map.

[0016] By adopting the above technical solution, the lithology type is associated with the DEM grid data, and a unique lithology code is assigned to each sub-grid unit. The stability index is divided into levels to intuitively display the stability status of the proposed project, with different levels representing different levels of risk. The R, G, and B channel data are merged to generate an RGB image, which can realize centralized information display and comprehensive analysis, avoiding the tediousness of processing data separately and improving information processing efficiency. The generated RGB image can intuitively display information such as elevation, lithology, and stability, facilitating spatial analysis and visualization. It can intuitively display the lithology distribution, elevation changes, and stability status of the proposed project, providing an important reference for geological hazard assessment.

[0017] In another possible implementation, formulating a protective measure mapping table based on the at least one risk level and the engineering protective measure includes: Based on the project attributes, determine at least one engineering protection measure; determining at least one risk value interval based on the risk value and the at least one risk level; The at least one engineering protection measure is recorded as a row, and the at least one risk value interval is recorded as a column to generate a protection measure mapping table.

[0018] By adopting the above technical solution, different types of projects have different characteristics and requirements. Determining protective measures based on project attributes ensures targeted and effective measures. The risk value quantifies the degree of risk, while the risk level is the standard for classifying and categorizing risk values. By analyzing the risk value and risk level, different risk ranges can be determined to better differentiate risks. A protective measures mapping table is generated, with project protective measures as rows and risk ranges as columns. This mapping table clearly displays the protective measures corresponding to different risk ranges, making it easier for engineering personnel to query during actual operations.

[0019] In a second aspect, the present application provides a method and device for comprehensive geological hazard assessment, which adopts the following technical solution: A method and device for comprehensive assessment of geological hazard risk, comprising: Data acquisition module, used to obtain topographic data, geomorphological data, rock and soil geological data of the area to be assessed, as well as the type and scope of the proposed project; A unit division module is used to divide the area to be evaluated into landform units based on the topographic data to obtain a DEM grid, wherein the DEM grid is elevation data stored in a grid format, and each grid cell corresponds to an elevation value; A comprehensive map generation module is used to match the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map; a risk value determination module, configured to correspond the type of the proposed project, the scope of the proposed project, and the lithology-elevation comprehensive map to generate a risk value; A measure determination module is used to determine engineering protection measures based on the risk value.

[0020] By adopting the above technical solution, since topographic and geomorphic data reflect surface morphological characteristics, geotechnical data reflect the regional geological structure and geotechnical properties, and the type and scope of the proposed project provide specific application scenarios for hazard assessment. Therefore, the data acquisition module acquires topographic and geomorphic data, geotechnical data, and the type and scope of the proposed project. This multi-source data integration enables a more comprehensive and accurate assessment of geological hazard risk, avoiding the limitations of single-data assessment. The unit division module divides the topographic and geomorphic data into geomorphic units, generating a DEM grid. This integrates continuous topographic and geomorphic data into regular grid cells, facilitating subsequent data processing. The DEM grid also provides a more intuitive understanding of the topographic undulations and topographic changes in the area being assessed. The comprehensive map generation module matches the geotechnical data with the DEM grid to generate a lithology-elevation comprehensive map. This map combines geotechnical properties with topographic elevation information to comprehensively reflect the distribution of geological characteristics at different elevations in the area being assessed, helping to identify areas of potential geological hazard risk. Specifically, the geologic characteristics of the area to be assessed are refined, making risk assessment more accurate and enabling the identification of potential small-scale geological hazard risk areas. Because different project types and scopes have varying impacts on the geological environment, the risk value determination module maps the proposed project type and scope to the lithology-elevation integrated map, enabling assessment of the potential risks of the proposed project under different geological conditions. The magnitude of the risk value reflects the likelihood and severity of a geological hazard. The measures determination module determines project protection measures based on the risk value. This allows for targeted engineering protection measures to be formulated based on the assessed risk value, improving protection effectiveness, reducing protection costs, and more effectively ensuring the safety of the proposed project.

[0021] In another possible implementation, when dividing the area to be assessed into landform units based on the topographic data, the unit division module is specifically configured to: Determine the slope data of the area to be assessed from the topographic data; Dividing the slope data into multiple slope levels; The ratio of the area corresponding to each slope grade to the total area of ​​the DEM grid is calculated to obtain the spatial distribution ratio of each geomorphic unit; The spatial distribution ratio and the multiple slope levels are used as a basis for DEM grid classification to perform landform unit division.

[0022] In another possible implementation, when matching the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map, the comprehensive map generation module is specifically used to: Determining lithology data based on the rock and soil geological data, wherein the lithology data includes lithology type, shear strength, and weathering degree; Superimposing the lithology type with the DEM grid to obtain a superimposed DEM grid, wherein each grid cell in the superimposed DEM grid corresponds to a lithology data; Calculating a stability index of each cell based on the shear strength, the degree of weathering, and the elevation weight of each grid cell; The DEM grid is rendered in a graded manner according to the stability index to generate a lithology-elevation integrated map.

[0023] In another possible implementation, when the risk value determination module generates the risk value by matching the proposed project type and the proposed project scope with the lithology-elevation integrated map, it is specifically configured to: Superimposing the proposed project scope with the lithology-elevation comprehensive map to obtain a lithology-elevation comprehensive map of the project impact area; Superimposing the proposed project scope with the lithology-elevation comprehensive map to obtain a lithology-elevation comprehensive map of the project impact area; Dividing the proposed project into regions to obtain a plurality of proposed project sub-regions; Obtain the classification parameters and project attributes of the proposed project in the proposed project sub-area and calculate the environmental impact coefficient; Calculate based on the infrastructure coverage rate and the type of the proposed project to obtain a project impact coefficient; Obtain the stability index of the proposed project based on the lithology-elevation comprehensive map of the project impact area; The risk value is calculated based on the stability index of the proposed project, the project impact coefficient and the environmental impact coefficient.

[0024] In another possible implementation, when determining the engineering protection measures based on the risk value, the measure determination module is specifically configured to: Classifying the risk value into risk levels to obtain at least one risk level; formulating a protective measures mapping table based on the at least one risk level and the engineering protective measures; Extracting risk values ​​for subunits of the proposed project from a comprehensive lithology-elevation map of the project's affected area; Corresponding engineering protection measures are determined in the protection measures mapping table based on the risk value of the proposed engineering subunit.

[0025] In another possible implementation, when the comprehensive map generation module performs hierarchical rendering on the DEM grid according to the stability index to generate the lithology-elevation comprehensive map, it is specifically used to: Determine the lithology code corresponding to each sub-grid unit in the DEM grid by using the lithology type and the DEM grid, wherein the lithology code is the number of each lithology, and the sub-grid unit is a subset of the DEM grid, that is, a plurality of the sub-grid units constitute the DEM grid; Classifying the stability index into grades to obtain a stability grade; Mapping the elevation value to the R channel, mapping the lithology code to the G channel, and mapping the stability grade to the B channel; Merge the R channel, G channel, and B channel to generate an RGB image; The RGB image is determined as a lithology-elevation comprehensive map.

[0026] In another possible implementation, when formulating the protective measure mapping table based on the at least one risk level and the engineering protective measure, the measure determination module is specifically configured to: Based on the project attributes, determine at least one engineering protection measure; determining at least one risk value interval based on the risk value and the at least one risk level; The at least one engineering protection measure is recorded as a row, and the at least one risk value interval is recorded as a column to generate a protection measure mapping table.

[0027] In a third aspect, the present application provides an electronic device, which adopts the following technical solution: An electronic device, comprising: at least one processor; Memory; At least one application, wherein at least one application is stored in a memory and configured to be executed by at least one processor, and at least one is configured to: execute a comprehensive geological hazard risk assessment method shown in any possible implementation of the first aspect.

[0028] In a fourth aspect, the present application provides a computer-readable storage medium, which adopts the following technical solution: A computer-readable storage medium, when the computer program is executed in a computer, causes the computer to execute a comprehensive geological hazard risk assessment method as described in any one of the first aspects. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 It is a flow chart of a comprehensive geological hazard risk assessment method in an embodiment of the present application.

[0030] Figure 2 It is a flow chart of a comprehensive geological disaster risk assessment device in an embodiment of the present application.

[0031] Figure 3 This is a flow chart of an electronic device in an embodiment of the present application. DETAILED DESCRIPTION

[0032] The following is combined with Figure 1-3 This application is described in further detail.

[0033] After reading this specification, those skilled in the art may make non-creative modifications to this embodiment as needed, but such modifications are protected by patent law as long as they fall within the scope of the claims of this application.

[0034] To make the purpose, technical solutions, and advantages of the embodiments of this application more clear, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.

[0035] In this document, the term "and / or" simply describes a relationship between related objects, indicating that three possible relationships exist. For example, "A and / or B" can represent: A exists alone, A and B exist simultaneously, or B exists alone. Furthermore, the character " / " in this document, unless otherwise specified, generally indicates an "or" relationship between the related objects.

[0036] The embodiments of the present application are described in further detail below with reference to the accompanying drawings.

[0037] The embodiment of the present application provides a method for comprehensive assessment of geological hazard risk, which is executed by an electronic device, which can be a server or a terminal device, wherein the server can be an independent physical server, a server cluster or a distributed system composed of multiple physical servers, or a cloud server providing cloud computing services. The terminal device can be a smart phone, a tablet computer, a laptop computer, a desktop computer, etc., but is not limited thereto. The terminal device and the server can be directly or indirectly connected via wired or wireless communication, which is not limited in the embodiment of the present application. Figure 1As shown, the method includes: step S101, step S102, step S103, step S104, step S105, wherein, Step S101: Obtain the topographic data, geomorphological data, rock and soil geological data, and the type and scope of the proposed project in the area to be assessed. Specifically, topographic data can be obtained through geospatial data clouds, satellite imagery, LiDAR, and drone aerial surveys. Geotechnical geodetic data can be obtained by accessing geodetic databases or through on-site surveys. The type and scope of the proposed project are determined by the construction company or design institute.

[0038] Step S102: Based on the topographic data, the area to be assessed is divided into geomorphic units to obtain a DEM grid.

[0039] Among them, the DEM grid is elevation data stored in raster form, and each grid cell corresponds to an elevation value.

[0040] In the embodiment of the present application, the topographic and geomorphic data are clipped according to the evaluation area, that is, the obtained topographic and geomorphic data are data of the area to be evaluated.

[0041] When demarcating geomorphic units, it is necessary to calculate the slope of each grid cell, determine the slope direction of each grid cell, calculate the surface curvature, and then analyze the rate of change of elevation with distance. The geomorphic units are converted into vector polygons, and then the rasterization tools of GIS software are used to convert the vector polygons into raster data.

[0042] For example, obtain regional DEM data and crop it to 500 meters on either side of the bridge centerline. Slopes are: 0°-10° (flat), 10°-30° (gentle slope), >30° (steep slope), and directions are: north, south, east, and west. Five landform units are identified: flat, gentle slope, steep slope, ridge, and valley. Generate a 5-meter resolution DEM grid, with each grid cell corresponding to a landform unit type (e.g., 1 = flat, 2 = gentle slope).

[0043] Step S103: Match the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map.

[0044] In this embodiment of the present application, the acquired geotechnical data needs to be rasterized. That is, the borehole lithologic data is interpolated into a raster format to form a lithologic grid with the same resolution as the DEM grid. The borehole lithologic data can be sandstone or shale. Using GIS overlay analysis, the lithologic grid is spatially aligned with the DEM grid, and each grid cell is assigned both a lithologic type and an elevation value to generate a comprehensive map.

[0045] For example, consider geotechnical data for the distribution of lithology along a tunnel, with sandstone at 0-20 meters and shale at 20-50 meters, at a resolution of 10 meters, and covering 200 meters on either side of the tunnel centerline. Interpolate the geotechnical data into a continuous lithology grid with a resolution of 10 meters. Ensure that the spatial extent of the lithology grid matches that of the DEM grid. Generate a 10-meter-resolution lithology-elevation composite map, displaying the lithology type and elevation value for each grid cell.

[0046] Step S104: Match the proposed project type and the proposed project scope with the lithology-elevation comprehensive map to generate a risk value.

[0047] In the embodiment of the present application, since different types of projects have different sensitivities to geological conditions, in order to avoid construction problems caused by geological problems, it is necessary to correspond the project type with the lithology-elevation comprehensive map and calculate the risk value, where the risk value represents the risk of construction in the area to be evaluated, and the higher the risk value, the greater the construction risk.

[0048] Step S105: Determine engineering protection measures based on the risk value.

[0049] In the embodiment of the present application, by quantifying the risk value and matching the protection measures, different risk levels correspond to different protection intensities, thereby reducing the risk of engineering accidents.

[0050] For example, select measures based on the risk ranking table:

[0051] When the risk value is 17.35, the risk level is high, and the corresponding protective measures are underground continuous wall + grouting reinforcement.

[0052] In another possible implementation method, the landform unit division of the assessment area based on topographic data includes: Determine the slope data of the area to be assessed from the topographic data; Divide into multiple slope levels according to slope data; The ratio of the area corresponding to each slope grade to the total area of ​​the DEM grid is calculated to obtain the spatial distribution ratio of each geomorphic unit; The spatial distribution ratio and multiple slope levels are used as the basis for DEM raster classification to divide the geomorphic units.

[0053] By adopting the above technical solution, since slope directly affects ecological processes such as soil erosion, landslide risk, and vegetation distribution, it is a basic parameter for landform classification. Therefore, in order to make landform classification more accurate, the slope data of the area to be assessed is determined from the topographic and geomorphological data, and the continuous terrain elevation data is converted into a quantitative indicator that can be analyzed. Since the slope data is directly used, when the slope data of the area to be assessed continuously increases or decreases, the landform type cannot be clearly determined. Therefore, the slope data is classified into levels to determine the landform type of the area to be assessed. By calculating the spatial distribution ratio, the dominant landform in the area to be assessed can be further clarified.

[0054] When calculating the slope, the formula slope = arctan (Δz / Δx) × 180 / π is used; Among them, Δz is the elevation and Δx is the DEM grid resolution.

[0055] For example, the slope is graded according to the Soil Erosion Classification Standard, namely:

[0056] Map each slope range to a corresponding level, i.e. 0-5 is level 1, 6-15 is level 2, and so on.

[0057] The area to be assessed is 10 km × 10 km, the DEM grid resolution is 30 m, and the total number of grids is: 10,000 × 10,000 / (30 × 30) ≈ 111,111 grids. One grid has a slope of 22°, corresponding to a slope of 3. The number of grids with a slope of 1 is 66,667, and the area occupied by grids with a slope of 1 is 60 km. 2 , the area occupied by the grid with slope grade 2 is 30km 2 , the area occupied by the grid with slope grade 1 is 10km 2 Since the area of ​​the area to be evaluated is 10km×10km, the spatial distribution ratio of slope grade 1 is 60 / 100=60%, the spatial distribution ratio of slope grade 2 is 30 / 100=30%, and the spatial distribution ratio of slope grade 3 is 10 / 100=10%. Each grid cell in the DEM grid is assigned a slope grade value, that is, the grid in the plain area is assigned a value of 1, the grid in the hilly area is assigned a value of 2, and the grid in the low mountain area is assigned a value of 3. After the assignment, the plain accounts for 60% of the DEM grid, the hilly area accounts for 30%, and the low mountain area accounts for 30%. It can be obtained that the landform of the area to be evaluated is dominated by plains.

[0058] Another possible implementation method is to match the geotechnical data with the DEM grid to generate a lithology-elevation integrated map including: Determine lithology data based on rock and soil geological data, including lithology type, shear strength, and weathering degree; The lithology type is superimposed on the DEM grid to obtain a superimposed DEM grid, in which each grid cell corresponds to a lithology data; Calculate the stability index of each cell based on the shear strength, weathering degree, and elevation weight of each grid cell; The DEM grid is rendered in grades according to the stability index to generate a lithology-elevation integrated map.

[0059] In the examples of this application, lithologic data is fundamental to geotechnical engineering and geological hazard assessment. Specifically, determining lithologic type is the basis for subsequent shear strength and weathering analysis. Shear strength reflects the rock's ability to resist shear failure and is directly controlled by lithologic type, making it a core metric for stability analysis. Weathering reflects the degree of rock degradation due to environmental factors and directly affects its mechanical properties. Therefore, lithologic data can more accurately assess the risk level of the area being assessed.

[0060] Furthermore, lithology represents the fundamental geological characteristics of a region, while shear strength reflects the differences in stability of different lithologies under elevation changes. Weathering degree can alter topography, which in turn affects the elevation weight. It should be noted that the degree of weathering can range from unweathered, weakly weathered, moderately weathered, strongly weathered, and fully weathered. To facilitate the calculation of the stability index, the degree of weathering is assigned a value. This assignment can be accomplished through mapping or based on historical data. For example, unweathered is 1, weakly weathered is 0.9, moderately weathered is 0.8, strongly weathered is 0.7, and fully weathered is 0.6. Elevation weights can be determined based on historical data: for elevations ≤ 300m, the weight is 1; for elevations 300m < ≤ 800m, the weight is 1.5; and for elevations > 800m, the weight is 2.

[0061] For example, if the area to be assessed is 10 km × 10 km and the DEM grid resolution is 30 m, the lithology data is converted to a raster format with a 30 m resolution. Using GIS tools, the lithology types are overlaid on the DEM grid pixels one by one, so that each grid cell in the DEM grid corresponds to a lithology type. For example, the lithology type corresponding to grid cell (1,1) is granite. Shear strength = 100 kPa, weathering degree = 0.8, and elevation = 500 m.

[0062] The stability index is calculated based on the formula = =100 / (1.5×0.8)=83.33; Rendering the DEM grid in different colors according to the stability index can visually display the rock and soil geological conditions in the lithology-elevation integrated map.

[0063] In another possible implementation method, the proposed project type and scope are mapped to the lithology-elevation integrated map to generate risk values ​​including: Superimpose the proposed project scope with the lithology-elevation comprehensive map to obtain a lithology-elevation comprehensive map of the project impact area; Dividing the proposed project into regions to obtain a plurality of proposed project sub-regions; Obtain the classification parameters and project attributes of the proposed project in the proposed project sub-area and calculate the environmental impact coefficient; Calculating based on the proposed project sub-area to obtain the infrastructure coverage rate of the proposed project; Calculate based on the infrastructure coverage rate and the type of the proposed project to obtain a project impact coefficient; Obtain the stability index of the proposed project based on the lithology-elevation comprehensive map of the project impact area; The risk value is calculated based on the stability index, engineering impact coefficient and environmental impact coefficient of the proposed project.

[0064] By employing the above technical solution and leveraging the spatial analysis capabilities of Geographic Information System (GIS) software, the proposed project scope is overlaid with a comprehensive lithology-elevation map. This overlay extracts lithology-elevation information within the project's impact area, creating a comprehensive lithology-elevation map of the project's impact area. Because different types of projects have varying degrees of impact on the surrounding geology, the project impact coefficient is determined based on the proposed project type, referring to relevant engineering specifications, empirical data, or expert evaluations. Different environmental conditions also affect different types of proposed projects to varying degrees. Therefore, the environmental impact coefficient is determined based on the proposed project type, taking into account the environmental characteristics of the project area, such as topography, climate, and ecological environment.

[0065] Obtain the classification parameters and project attributes of the proposed project in the proposed project sub-area and calculate the environmental impact coefficient according to the following formula: Calculate the adaptability index λ of the proposed project:

[0066] For example, a proposed expressway project in a mountainous area encompasses three sub-areas (A, B, and C). The weight vectors for the project's influencing factors are terrain slope (x1), annual rainfall (x2), and geotechnical stability (x3). The project type classification parameter C = {1.2, 1.5, 1.8}, and the project attribute P = {1.0, 1.3, 1.6}. For sub-area A, the road coverage ratio Beff = 0.6, and the pipeline coverage ratio Beff = 0.4. For sub-area B, the road coverage ratio Beff = 0.8, and the pipeline coverage ratio Beff = 0.5. For sub-area C, the road coverage ratio Beff = 0.3, and the pipeline coverage ratio Beff = 0.2.

[0067] Assume that the matrix of factors affecting the proposed project is:

[0068] Among them, the environmental impact coefficient of the proposed project is the average value of the environmental impact coefficients of the sub-areas of the proposed project.

[0069] Furthermore, the project impact coefficient is calculated based on the infrastructure coverage rate and the type of the proposed project. The specific calculation steps are as follows: Determine the infrastructure coverage Beff:

[0070] Calculation of engineering impact coefficient: G=β×P×(1-Beff)+γ×Zeff; Among them, β and γ are weight factors, that is, β+γ=1; P is the engineering property; Zeff is the geological environment sensitivity factor.

[0071] Specifically, calculate the infrastructure coverage Beff: Sub-area A: Similarly, sub-area B=0.65 and sub-area C=0.25 are calculated.

[0072] Calculate the engineering impact coefficient: Assume β=0.6, γ=0.4, and geological sensitivity factor Zeff=0.8: G A =0.6×1.0×(1−0.5)+0.4×0.8=0.3+0.32=0.62, calculate G in the same way B =0.56, G C =0.74.

[0073] The engineering impact coefficient of the proposed project is the average of the engineering impact coefficients of the proposed project sub-areas.

[0074] Suppose the proposed project is a tunnel with a stability index of 83.33 and a project impact coefficient of 0.6. Since tunnel construction may impact the surrounding groundwater and ecosystem, the environmental impact coefficient is 0.7. Since the stability index is a calculated value and not a standardized value, for uniformity, it is converted to a value between 0 and 1: 83.33 / 100 = 0.8333. Since the stability index, project impact coefficient, and environmental impact coefficient have different impacts on the proposed project, to improve accuracy, each is calculated with its own weight. For example, if the weight of the stability index is 0.3, the weight of the project impact coefficient is 0.3, and the weight of the environmental impact coefficient is 0.4, then the risk value is 0.8333 × 0.3 + 0.6 × 0.3 + 0.7 × 0.4 = 0.7099.

[0075] In another possible implementation method, determining engineering protection measures based on risk values ​​includes: Divide the risk value into risk levels to obtain at least one risk level; Develop a mapping table of protective measures based on at least one risk level and engineering protective measures; Extract risk values ​​for proposed project subunits from the lithology-elevation composite map of the project impact area; Based on the risk value of the proposed project subunit, the corresponding engineering protection measures are determined in the protection measures mapping table.

[0076] By employing the above technical solution, thresholds for different risk levels can be set based on engineering experience, historical data, or industry standards. For example, risk values ​​can be categorized into four levels: low (0-0.3), medium (0.3-0.6), high (0.6-0.8), and extremely high (0.8-1.0). Each project is assigned to the corresponding risk level based on the calculated risk value. For example, a project subunit with a risk value of 0.7099 is classified as high risk. Project protection measures for different risk levels are collected, such as strengthened support, additional drainage facilities, and more stable construction methods. A mapping relationship is established between risk levels and corresponding project protection measures, creating a protection measure mapping table. For example, a high risk level might correspond to measures such as strengthened support and increased monitoring frequency. The proposed project subunit is a grid cell in the lithology-elevation map of the project's impact area. In other words, the lithology-elevation map of the project's impact area is composed of multiple proposed project subunits. Since the risk levels of different construction locations in the same proposed project are different and the required protective measures are also different, in order to determine the protective measures more accurately and avoid hazards caused by construction, each proposed project sub-unit is matched with the project protective measures.

[0077] In another possible implementation, the DEM grid is graded and rendered according to the stability index to generate a lithology-elevation integrated map, including: The lithology code corresponding to each sub-grid unit in the DEM grid is determined by the lithology type and the DEM grid. The lithology code is the number of each lithology, and the sub-grid unit is a subset of the DEM grid, that is, multiple sub-grid units constitute the DEM grid; The stability index is divided into grades to obtain the stability grade; Map the elevation value to the R channel, the lithology code to the G channel, and the stability grade to the B channel; Merge the R channel, G channel and B channel to generate an RGB image; The RGB image is determined as a lithology-elevation integrated map.

[0078] By using the above technical solution, each lithology type is assigned a unique code. The lithology type data is spatially overlaid with the DEM grid data, and the lithology code corresponding to each DEM sub-grid cell is determined by spatial position matching. Assume there are three lithology types: granite (coded 1), limestone (coded 2), and shale (coded 3). At the location of a DEM grid cell, if the spatial overlay shows that the lithology corresponding to that location is limestone, the lithology code for that grid cell is 2.

[0079] Based on the numerical range of the stability index, the stability index is divided into several levels, and each level is assigned a unique identifier. Specifically, the level division method can be determined according to actual needs.

[0080] For example, the stability index ranges from 0 to 100, which is divided into 5 levels, level 1: 0-20 (very unstable); level 2: 21-40 (unstable); level 3: 41-60 (moderately stable); level 4: 61-80 (stable); level 5: 8-100 (very stable).

[0081] The R channel is the red channel, the G channel is the green channel, and the B channel is the blue channel. Specifically, the elevation values ​​are normalized to a range between 0 and 255 and then mapped to the R channel. The lithology codes are directly mapped to the G channel. To ensure clear color distinction between different lithologies, the lithology codes can be processed, such as multiplying them by a coefficient to a range between 0 and 255. The stability level is encoded to a range between 0 and 255 and then mapped to the B channel.

[0082] For example, in the R channel, the elevation value of a grid cell is 500 meters, and the elevation range of the entire DEM area is 0 to 1000 meters. The mapping coefficient is 255. After normalization, the elevation value of the grid cell is mapped to (50 / 1000)×255=127.5, rounded to 128, and mapped to the R channel.

[0083] G channel: Assuming that the lithology code of the grid cell is 2 (limestone), in order to map it to between 0 and 255, it can be multiplied by a coefficient. Assuming the mapping coefficient is 85, then 2×85=170, which is mapped to the G channel.

[0084] B channel: Assume that the stability index of this grid cell is 70, which belongs to level 4 (stable). Its mapping coefficient is 51, and level 4 is encoded as 4×51=204 and mapped to the B channel.

[0085] Combine the values ​​of the R channel, G channel, and B channel into an RGB triplet, perform a synthesis operation on all grid cells, and generate an image matrix containing the RGB values ​​of all grid cells. For example, for the grid cell in the above example, its RGB value is (128, 170, 204). Traverse the entire DEM grid to generate a complete RGB image matrix. Save the generated RGB image matrix as an image file, such as PNG or JPEG format. This image is a lithology-elevation composite map, in which different colors represent different elevation, lithology, and stability combinations, and then intuitively show the workers the degree of danger of the proposed project through different colors.

[0086] It should be noted that the R channel represents elevation, the G channel represents lithology code, and the B channel represents stability level. Therefore, higher elevation values ​​indicate a stronger red component. Different lithology types are represented by varying green intensities. Higher stability indicates a stronger blue component.

[0087] In another possible implementation, formulating a protective measure mapping table based on at least one risk level and engineering protective measures includes: Based on the project attributes, determine at least one engineering protection measure; determining at least one risk value interval based on the risk value and at least one risk level; At least one engineering protection measure is recorded as a row, and at least one risk value interval is recorded as a column to generate a protection measure mapping table.

[0088] By adopting the above technical solutions, due to different project types, different geographical environments, and inconsistent levels of danger, the protective measures that need to be taken are also different. Therefore, corresponding protective measures need to be formulated according to the properties of the project. For example, when the project type is slope protection, for mountain projects, it may be necessary to consider the stability of the slope and adopt measures such as anchor rods and retaining walls. When the project type is flood control measures, for projects located near rivers, flood control measures such as dams and drainage systems need to be considered. When the project type is seismic measures, for projects in earthquake-prone areas, seismic design such as shock absorbers and seismic isolation bearings need to be adopted.

[0089] Based on the risk value, the risk level can be divided. For example, the risk value can be divided into three levels: low risk (0-0.3), medium risk (0.3-0.6), and high risk (0.6-0.8). Create a table with engineering protection measures as rows and risk value ranges as columns. Fill in the corresponding protection measures for each risk value range in the table. The details are as follows:

[0090]

[0091] In a second aspect, the present application provides a method and device for comprehensive geological hazard assessment, which adopts the following technical solution: A method and apparatus 20 for comprehensively assessing geological hazard risk, comprising: Data acquisition module 201, used to obtain topographic data, geomorphological data, rock and soil geological data, and the type and scope of the proposed project in the area to be assessed; The unit division module 202 is used to divide the assessment area into geomorphic units based on the topographic data to obtain a DEM grid. The DEM grid is elevation data stored in a grid format, and each grid cell corresponds to an elevation value. A comprehensive map generation module 203 is used to match the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map; The risk value determination module 204 is used to match the proposed project type and the proposed project scope with the lithology-elevation comprehensive map to generate a risk value; The measure determination module 205 is used to determine engineering protection measures based on the risk value.

[0092] By adopting the above technical solution, since topographic and geomorphic data reflect surface morphological characteristics, geotechnical data reflect the regional geological structure and geotechnical properties, and the type and scope of the proposed project provide specific application scenarios for hazard assessment. Therefore, the data acquisition module acquires topographic and geomorphic data, geotechnical data, and the type and scope of the proposed project. This multi-source data integration enables a more comprehensive and accurate assessment of geological hazard risk, avoiding the limitations of single-data assessment. The unit division module divides the topographic and geomorphic data into geomorphic units, generating a DEM grid. This integrates continuous topographic and geomorphic data into regular grid cells, facilitating subsequent data processing. The DEM grid also provides a more intuitive understanding of the topographic undulations and topographic changes in the area being assessed. The comprehensive map generation module matches the geotechnical data with the DEM grid to generate a lithology-elevation comprehensive map. This map combines geotechnical properties with topographic elevation information to comprehensively reflect the distribution of geological characteristics at different elevations in the area being assessed, helping to identify areas of potential geological hazard risk. Specifically, the geologic characteristics of the area to be assessed are refined, making risk assessment more accurate and enabling the identification of potential small-scale geological hazard risk areas. Because different project types and scopes have varying impacts on the geological environment, the risk value determination module maps the proposed project type and scope to the lithology-elevation integrated map, enabling assessment of the potential risks of the proposed project under different geological conditions. The magnitude of the risk value reflects the likelihood and severity of a geological hazard. The measures determination module determines project protection measures based on the risk value. This allows for targeted engineering protection measures to be formulated based on the assessed risk value, improving protection effectiveness, reducing protection costs, and more effectively ensuring the safety of the proposed project.

[0093] In another possible implementation, when dividing the area to be assessed into landform units based on the topographic data, the unit division module 202 is specifically configured to: Determine the slope data of the area to be assessed from the topographic data; Divide into multiple slope levels according to slope data; The ratio of the area corresponding to each slope grade to the total area of ​​the DEM grid is calculated to obtain the spatial distribution ratio of each geomorphic unit; The spatial distribution ratio and multiple slope levels are used as the basis for DEM raster classification to divide the geomorphic units.

[0094] In another possible implementation, when matching the rock and soil geological data with the DEM grid to generate the lithology-elevation comprehensive map, the comprehensive map generation module 203 is specifically used to: Determine lithology data based on rock and soil geological data, including lithology type, shear strength, and weathering degree; The lithology type is superimposed on the DEM grid to obtain a superimposed DEM grid, in which each grid cell corresponds to a lithology data; Calculate the stability index of each cell based on the shear strength, weathering degree, and elevation weight of each grid cell; The DEM grid is rendered in grades according to the stability index to generate a lithology-elevation integrated map.

[0095] In another possible implementation, when the risk value determination module 204 generates the risk value by mapping the proposed project type and the proposed project scope with the lithology-elevation integrated map, it is specifically configured to: Superimpose the proposed project scope with the lithology-elevation comprehensive map to obtain a lithology-elevation comprehensive map of the project impact area; Determine the engineering impact coefficient and environmental impact coefficient based on the type of project to be built; Obtain the stability index of the proposed project based on the lithology-elevation comprehensive map of the project impact area; The risk value is calculated based on the stability index, engineering impact coefficient and environmental impact coefficient of the proposed project.

[0096] In another possible implementation, when determining engineering protection measures based on the risk value, the measure determination module 205 is specifically configured to: Divide the risk value into risk levels to obtain at least one risk level; Develop a mapping table of protective measures based on at least one risk level and engineering protective measures; Extract risk values ​​for proposed project subunits from the lithology-elevation composite map of the project impact area; Based on the risk value of the proposed project subunit, the corresponding engineering protection measures are determined in the protection measures mapping table.

[0097] In another possible implementation, when the comprehensive map generation module 203 performs hierarchical rendering on the DEM grid according to the stability index to generate the lithology-elevation comprehensive map, it is specifically used to: The lithology code corresponding to each sub-grid unit in the DEM grid is determined by the lithology type and the DEM grid. The lithology code is the number of each lithology, and the sub-grid unit is a subset of the DEM grid, that is, multiple sub-grid units constitute the DEM grid; The stability index is divided into grades to obtain the stability grade; Map the elevation value to the R channel, the lithology code to the G channel, and the stability grade to the B channel; Merge the R channel, G channel and B channel to generate an RGB image; The RGB image is determined as a lithology-elevation integrated map.

[0098] In another possible implementation, when formulating the protective measure mapping table based on at least one risk level and engineering protective measure, the measure determination module 205 is specifically configured to: Based on the project attributes, determine at least one engineering protection measure; determining at least one risk value interval based on the risk value and at least one risk level; At least one engineering protection measure is recorded as a row, and at least one risk value interval is recorded as a column to generate a protection measure mapping table.

[0099] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and units described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0100] An electronic device is provided in an embodiment of the present application, such as Figure 3 As shown, Figure 3 The electronic device 30 shown includes a processor 301 and a memory 303. The processor 301 and the memory 303 are connected, for example, via a bus 302. Optionally, the electronic device 30 may further include a transceiver 304. It should be noted that in actual applications, the number of transceivers 304 is not limited to one, and the structure of the electronic device 30 does not constitute a limitation on the embodiments of the present application.

[0101] Processor 301 may be a CPU (Central Processing Unit), a general-purpose processor, a DSP (Digital Signal Processor), an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array), or other programmable logic device, transistor logic device, hardware component, or any combination thereof. It may implement or execute the various exemplary logic blocks, modules, and circuits described in conjunction with the disclosure of this application. Processor 301 may also be a combination that implements computing functions, such as a combination of one or more microprocessors, or a combination of a DSP and a microprocessor.

[0102] The bus 302 may include a path for transmitting information between the above components. The bus 302 may be a PCI (Peripheral Component Interconnect) bus or an EISA (Extended Industry Standard Architecture) bus. The bus 302 may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 3 Only one thick line is used in the diagram, but it does not mean that there is only one bus or one type of bus.

[0103] The memory 303 may be a ROM (Read Only Memory) or other types of static storage devices that can store static information and instructions, a RAM (Random Access Memory) or other types of dynamic storage devices that can store information and instructions, or an EEPROM (Electrically Erasable Programmable Read Only Memory), a CD-ROM (Compact Disc Read Only Memory) or other optical disk storage, optical disk storage (including compact discs, laser discs, optical discs, digital versatile discs, Blu-ray discs, etc.), a magnetic disk storage medium or other magnetic storage device, or any other medium that can be used to carry or store desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto.

[0104] The memory 303 is used to store application code for executing the solution of the present application, and the execution is controlled by the processor 301. The processor 301 is used to execute the application code stored in the memory 303 to implement the content shown in the above method embodiment.

[0105] Electronic devices include, but are not limited to, mobile terminals such as mobile phones, laptop computers, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), and in-vehicle terminals (e.g., in-vehicle navigation terminals), as well as fixed terminals such as digital TVs and desktop computers. They may also include servers, etc. Figure 3 The electronic device shown is merely an example and should not limit the functions and scope of use of the embodiments of the present application.

[0106] The embodiment of the present application provides a computer-readable storage medium having a computer program stored thereon, which, when run on a computer, enables the computer to execute the corresponding contents in the aforementioned method embodiment. Compared with the related art, in the embodiment of the present application, since the topographic and geomorphic data reflect the surface morphological characteristics, the geotechnical data reflect the geological structure and geotechnical properties of the region, the type of the proposed project and the scope of the proposed project provide a specific application scenario for risk assessment. Therefore, by integrating the topographic and geomorphic data, geotechnical data, and the type of the proposed project and the scope of the proposed project, the risk of geological disasters can be assessed more comprehensively and accurately, avoiding the limitations of single data assessment. Based on the topographic and geomorphic data, the geotechnical data are divided into geomorphic units to obtain a DEM grid, which can integrate the continuous topographic and geomorphic data into regular grid units, which is convenient for subsequent data processing, and the DEM grid can be used to more intuitively understand the terrain undulations and terrain changes in the area to be assessed. Matching geotechnical data with a DEM raster to generate a lithology-elevation integrated map is designed to integrate geotechnical properties with terrain elevation information, comprehensively reflecting the distribution of geological characteristics at different elevations in the area under assessment. This helps identify areas of potential geological hazard risk. This refined representation of the geological characteristics of the area under assessment enables more accurate risk assessment and the identification of potential small-scale geological hazard risk areas. Because different project types and scopes have different impacts on the geological environment, matching the proposed project type and scope with the lithology-elevation integrated map allows for the assessment of potential risks under different geological conditions. The magnitude of the risk value reflects the likelihood and severity of a geological hazard. Determining engineering protective measures based on the risk value allows for targeted development of protective measures based on the assessed risk value, improving protection effectiveness, reducing protection costs, and more effectively ensuring the safety of the proposed project.

[0107] It should be understood that although the steps in the flowcharts of the accompanying drawings are shown in sequence as indicated by the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some of the steps in the flowcharts of the accompanying drawings may include multiple sub-steps or multiple stages, and these sub-steps or stages are not necessarily executed at the same time, but can be executed at different times, and their execution order is not necessarily sequential, but can be executed in turn or alternately with other steps or at least a portion of the sub-steps or stages of other steps.

[0108] The above are only some of the implementation methods of the present application. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present application. These improvements and modifications should also be regarded as the scope of protection of the present application.

Claims

1. A comprehensive assessment method for geological hazard risk, characterized in that: include: Obtain topographic and geomorphological data, geotechnical data, and the type and scope of the proposed project in the area to be assessed; Based on the topographic data, the area to be assessed is divided into geomorphic units to obtain a DEM grid, wherein the DEM grid is elevation data stored in a grid format, and each grid cell corresponds to an elevation value; Matching the rock and soil geological data with the DEM grid to generate a lithology-elevation composite map; Matching the proposed project type and the proposed project scope with the lithology-elevation comprehensive map to generate a risk value; Engineering protection measures are determined based on the risk value.

2. A method for comprehensive assessment of geological hazard risk according to claim 1, characterized in that: The dividing the area to be assessed into geomorphic units based on the topographic data includes: Determine the slope data of the area to be assessed from the topographic data; Dividing the slope data into multiple slope levels; The ratio of the area corresponding to each slope grade to the total area of ​​the DEM grid is calculated to obtain the spatial distribution ratio of each geomorphic unit; The spatial distribution ratio and the multiple slope levels are used as a basis for DEM grid classification to perform landform unit division.

3. A method for comprehensive assessment of geological hazard risk according to claim 1, characterized in that: The matching of the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map includes: Determining lithology data based on the rock and soil geological data, wherein the lithology data includes lithology type, shear strength, and weathering degree; Superimposing the lithology type with the DEM grid to obtain a superimposed DEM grid, wherein each grid cell in the superimposed DEM grid corresponds to a lithology data; Calculating a stability index of each cell based on the shear strength, the degree of weathering, and the elevation weight of each grid cell; The DEM grid is rendered in a graded manner according to the stability index to generate a lithology-elevation integrated map.

4. A method for comprehensive assessment of geological hazard risk according to claim 3, characterized in that: The generating of risk values ​​by matching the proposed project type and the proposed project scope with the lithology-elevation comprehensive map includes: Superimposing the proposed project scope with the lithology-elevation comprehensive map to obtain a lithology-elevation comprehensive map of the project impact area; Dividing the proposed project into regions to obtain a plurality of proposed project sub-regions; Obtain the classification parameters and project attributes of the proposed project in the proposed project sub-area and calculate the environmental impact coefficient; Calculate based on the infrastructure coverage rate and the type of the proposed project to obtain a project impact coefficient; Obtain the stability index of the proposed project based on the lithology-elevation comprehensive map of the project impact area; The risk value is calculated based on the stability index of the proposed project, the project impact coefficient and the environmental impact coefficient.

5. A method for comprehensive assessment of geological hazard risk according to claim 4, characterized in that: Determining engineering protection measures based on the risk value includes: Classifying the risk value into risk levels to obtain at least one risk level; formulating a protective measures mapping table based on the at least one risk level and the engineering protective measures; Extracting risk values ​​for subunits of the proposed project from a comprehensive lithology-elevation map of the project's affected area; Corresponding engineering protection measures are determined in the protection measures mapping table based on the risk value of the proposed engineering subunit.

6. A method for comprehensive assessment of geological hazard risk according to claim 3, characterized in that: The step of rendering the DEM grid in a hierarchical manner according to the stability index to generate a lithology-elevation comprehensive map includes: Determine the lithology code corresponding to each sub-grid unit in the DEM grid by using the lithology type and the DEM grid, wherein the lithology code is the number of each lithology, and the sub-grid unit is a subset of the DEM grid, that is, a plurality of the sub-grid units constitute the DEM grid; Classifying the stability index into grades to obtain a stability grade; Mapping the elevation value to the R channel, mapping the lithology code to the G channel, and mapping the stability grade to the B channel; Merge the R channel, G channel, and B channel to generate an RGB image; The RGB image is determined as a lithology-elevation comprehensive map.

7. A method for comprehensive assessment of geological hazard risk according to claim 5, characterized in that: The formulating a protective measure mapping table based on the at least one risk level and the engineering protective measure includes: Based on the project attributes, determine at least one engineering protection measure; determining at least one risk value interval based on the risk value and the at least one risk level; The at least one engineering protection measure is recorded as a row, and the at least one risk value interval is recorded as a column to generate a protection measure mapping table.

8. A comprehensive geological hazard assessment device, characterized in that: include: Data acquisition module, used to obtain topographic data, geomorphological data, rock and soil geological data of the area to be assessed, as well as the type and scope of the proposed project; A unit division module is used to divide the area to be evaluated into landform units based on the topographic data to obtain a DEM grid, wherein the DEM grid is elevation data stored in a grid format, and each grid cell corresponds to an elevation value; A comprehensive map generation module is used to match the rock and soil geological data with the DEM grid to generate a lithology-elevation comprehensive map; a risk value determination module, configured to correspond the type of the proposed project, the scope of the proposed project, and the lithology-elevation comprehensive map to generate a risk value; A measure determination module is used to determine engineering protection measures based on the risk value.

9. An electronic device, characterized in that: It includes: at least one processor; Memory; At least one application, wherein the at least one application is stored in the memory and is configured to be executed by the at least one processor, and the at least one application is used to execute a comprehensive geological hazard assessment method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed in a computer, the computer is caused to execute a method for comprehensive geological hazard assessment according to any one of claims 1 to 7.

Images