Slope three-dimensional geologic model establishment and stability calculation method based on airborne geophysical prospecting

Abstract

The invention discloses a side slope three-dimensional geologic model establishment and stability calculation method based on airborne geophysical prospecting, and belongs to the technical field of side slope landslide early warning, and the method comprises the following steps: a, obtaining side slope surface elevation and underground resistivity spatial data; b, rock-soil body characteristics are obtained, and resistivity values corresponding to the rock-soil body characteristics are extracted for calibration; c, extracting surface elevation and aeronautical geophysical prospecting resistivity spatial point data through spatial analysis three-dimensional software, extracting corresponding interface points according to a calibrated range, constructing a side slope three-dimensional structural surface from points to surfaces, and constructing a side slope three-dimensional geologic model in combination with the ground surface; and d, establishing a slope three-dimensional grid stability calculation model, and calculating a slope stability coefficient. According to the method, the slope stability coefficient can be accurately calculated, so that the three-dimensional stability of the slope can be rapidly and accurately evaluated.

Description

Technical Field

[0001] The invention relates to the technical field of slope landslide early warning, and in particular to a method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration. Background Art

[0002] With the rapid development of drones, aerial geophysical exploration can quickly replace traditional drilling and ground geophysical exploration to achieve the exploration of soft structural surface information such as bedrock surfaces and groundwater levels over a large area, providing real and reliable underground data for the construction of three-dimensional geological models of slopes.

[0003] How to efficiently use the aerial geophysical data in a wide area to accurately construct a wide-area three-dimensional geological model of the slope, and then select an appropriate stability calculation method to evaluate the three-dimensional stability of the slope is of great significance for landslide early warning.

[0004] A Chinese patent document with publication number CN119514744A and publication date of February 25, 2025 discloses a machine learning prediction method for high and steep slope landslides in mines based on drone technology and multi-source data fusion, including the following steps: The first step is data collection: using drones to perform terrain simulation flights to collect multi-source data of the slope area, including high-resolution digital surface models taken by drones; building a three-dimensional geological model based on historical mine research reports and data, using mine geological planing maps and geological drilling data; cataloging the specific location and size of large landslide points based on historical mine landslide data and field surveys; The second step is data fusion: after removing non-topographic features, convert it into digital elevation model data, obtain the surface refinement model of the open-pit mine in the digital elevation model, use the surface refinement model as a constraint to trim the surface of the three-dimensional geological model, and analyze the elevation, curvature and slope factors; apply the non-continuous factors affecting lithology and faults to the three-dimensional geological model for analysis; fuse the analyzed factors to construct complete slope data; The third step is model training: select an appropriate machine learning algorithm for training, input the complete slope data fused in the second step into the training model, perform training and verification, and thus generate a prediction model for slope landslide susceptibility; Step 4: Model optimization: Compare the results of different machine learning models, select the two models with the highest accuracy, and use soft voting to propose a high-precision coupled machine learning model to improve prediction accuracy. Step 5. Result generation and application: The optimized model is used to predict the landslide susceptibility of new slope data. The output of the model includes the probability of landslide occurrence, the scope of impact and the time when it may occur. A landslide susceptibility assessment report is generated and the prediction results are visualized as maps or charts for reference by engineers and decision makers.

[0005] The patent document discloses a machine learning prediction method for landslides on steep mine slopes based on drone technology and multi-source data fusion, which realizes dynamic risk management and ensures the safety and sustainable development of mine production. However, since the slope stability coefficient cannot be accurately calculated, the three-dimensional stability of the slope cannot be quickly and accurately evaluated. Summary of the invention

[0006] In order to overcome the defects of the above-mentioned prior art, the present invention provides a method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration. The present invention can accurately calculate the slope stability coefficient, thereby quickly and accurately evaluating the three-dimensional stability of the slope.

[0007] The present invention is achieved through the following technical solutions: The three-dimensional geological model of the slope and the stability calculation method based on airborne geophysical exploration include the following steps: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0008] In a, the surface elevation of the slope is obtained by drone aerial photography.

[0009] In a, the underground resistivity spatial data is obtained by semi-aerial electromagnetic method.

[0010] In b, the rock and soil characteristics include lithological information of the rock and soil strata, thickness of the overlying soil layer and moisture content of the soil.

[0011] In the above d, the anti-slip force is calculated by formula 2; Formula 2; in, For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; For the The internal friction angle of the structural surface of a three-dimensional geological grid, unit: °; For the The cohesion of the structural surface of a three-dimensional geological grid, in kPa; is the unit length of the three-dimensional geological grid, in meters; It is the unit width of the three-dimensional geological grid, in meters.

[0012] In the above d, the sliding force is calculated by formula 3; Formula 3; in, For the The sliding force caused by the gravity of the three-dimensional geological grid and other external forces, unit kN, For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN.

[0013] The deadweight of the geological grid is calculated by equation 4; Formula 4; in, For the The self-weight of a three-dimensional geological grid, in kN; For the The soil mass of a three-dimensional geological grid, in kN / m³; For the The center height of a three-dimensional geological grid, in meters; is the unit width of the three-dimensional geological grid, in meters; is the unit length of the three-dimensional geological grid, in meters; For the The terrain slope of a 3D geological grid, in degrees.

[0014] The horizontal additional force is calculated by formula 5; Formula 5; in, For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; is the horizontal additional force of the three-dimensional geological model of the slope, in kN; It is the angle between the inclination of the 3D geological grid and the horizontal additional force of the 3D geological model of the slope, in degrees.

[0015] The beneficial effects of the present invention are mainly manifested in the following aspects: 1. Compared with the prior art, the present invention can accurately calculate the slope stability coefficient, thereby quickly and accurately evaluating the three-dimensional stability of the slope.

[0016] 2. The present invention combines the refined surface elevation data and aerial geophysical data, and based on the 3D modeling software, quickly constructs a real 3D geological model of the slope in a wide area, which can intuitively and accurately reflect the structural characteristics of the slope, such as the base-cover interface, groundwater level and soil moisture content, and provide potential sliding weak structural surfaces for slope stability calculation, which is beneficial to the regional 3D stability evaluation of slopes.

[0017] 3. The present invention takes into account the real geometric model and potential sliding direction of the three-dimensional geological grid body, assumes that each grid body has only one sliding direction, the lateral interface is in a stress equilibrium state, and the sliding force is parallel to the interface layer. The stability coefficient of each grid body is calculated through static equilibrium conditions, thereby obtaining the stability results of the slope in a wide area. It can be applied to the three-dimensional geological model obtained by aerial geophysical exploration, and can quickly and accurately evaluate the three-dimensional stability of the slope.

[0018] 4. The present invention, based on the entity's three-dimensional geological model, makes the calculation object more realistic and the calculation result more accurate and reliable without affecting the calculation time and complexity.

[0019] 5. The present invention improves the authenticity of the three-dimensional geological model of the slope, thereby helping to improve the accuracy of the calculation results of the three-dimensional stability of the slope. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments: Figure 1 It is a schematic diagram of the process of the present invention; Figure 2 It is a three-dimensional spatial schematic diagram of the slope based on aerial geophysical exploration of the present invention; Figure 3 It is a schematic diagram of a three-dimensional geological model of a slope calibrated based on aerial geophysical prospecting according to the present invention; Figure 4 It is a schematic diagram of the force analysis of the three-dimensional geological grid of the slope of the present invention. DETAILED DESCRIPTION

[0021] Example 1 The three-dimensional geological model of the slope and the stability calculation method based on airborne geophysical exploration include the following steps: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0022] This embodiment is the most basic implementation method. Compared with the prior art, it can accurately calculate the slope stability coefficient, thereby quickly and accurately evaluating the three-dimensional stability of the slope.

[0023] Example 2 The three-dimensional geological model of the slope and the stability calculation method based on airborne geophysical exploration include the following steps: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0024] In a, the surface elevation of the slope is obtained by drone aerial photography.

[0025] In a, the underground resistivity spatial data is obtained by semi-aerial electromagnetic method.

[0026] In b, the rock and soil characteristics include lithological information of the rock and soil strata, thickness of the overlying soil layer and moisture content of the soil.

[0027] This embodiment is a preferred implementation method. It combines the refined surface elevation data and aerial geophysical data, and quickly constructs a real three-dimensional geological model of the slope in a wide area based on the three-dimensional modeling software. It can intuitively and accurately reflect the structural characteristics of the slope, such as the base-cover interface, groundwater level and soil moisture content in the slope, and provide potential sliding weak structural surfaces for slope stability calculation, which is beneficial to the regional three-dimensional stability evaluation of slopes.

[0028] Example 3 The three-dimensional geological model of the slope and the stability calculation method based on airborne geophysical exploration include the following steps: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Extract the surface elevation and aerial geophysical resistivity spatial point data through spatial analysis 3D software, extract the corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from the point to the surface, and construct the 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0029] In a, the surface elevation of the slope is obtained by drone aerial photography.

[0030] In a, the underground resistivity spatial data is obtained by semi-aerial electromagnetic method.

[0031] In b, the rock and soil characteristics include lithological information of the rock and soil strata, thickness of the overlying soil layer and moisture content of the soil.

[0032] In the above d, the anti-slip force is calculated by formula 2; Formula 2; in, For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; For the The internal friction angle of the structural surface of a three-dimensional geological grid, unit: °; For the The cohesion of the structural surface of a three-dimensional geological grid, in kPa; is the unit length of the three-dimensional geological grid, in meters; It is the unit width of the three-dimensional geological grid, in meters.

[0033] This embodiment is another preferred implementation method, which takes into account the real geometric model of the three-dimensional geological grid body and the potential sliding direction. It is assumed that each grid body has only one sliding direction, the lateral interface is in a stress equilibrium state, and the sliding force is parallel to the interface layer. The stability coefficient of each grid body is calculated through static equilibrium conditions, thereby obtaining the stability results of the slope in a wide area. It can be applied to the three-dimensional geological model obtained by aerial geophysical exploration, and can quickly and accurately evaluate the three-dimensional stability of the slope.

[0034] Example 4 The three-dimensional geological model of the slope and the stability calculation method based on airborne geophysical exploration include the following steps: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0035] In a, the surface elevation of the slope is obtained by drone aerial photography.

[0036] In a, the underground resistivity spatial data is obtained by semi-aerial electromagnetic method.

[0037] In b, the rock and soil characteristics include lithological information of the rock and soil strata, thickness of the overlying soil layer and moisture content of the soil.

[0038] In the above d, the anti-slip force is calculated by formula 2; Formula 2; in, For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; For the The internal friction angle of the structural surface of a three-dimensional geological grid, unit: °; For the The cohesion of the structural surface of a three-dimensional geological grid, in kPa; is the unit length of the three-dimensional geological grid, in meters; It is the unit width of the three-dimensional geological grid, in meters.

[0039] In the above d, the sliding force is calculated by formula 3; Formula 3; in, For the The sliding force caused by the gravity of the three-dimensional geological grid and other external forces, unit kN, For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN.

[0040] The deadweight of the geological grid is calculated by equation 4; Formula 4; in, For the The self-weight of a three-dimensional geological grid, in kN; For the The soil mass of a three-dimensional geological grid, in kN / m³; For the The center height of a three-dimensional geological grid, in meters; is the unit width of the three-dimensional geological grid, in meters; is the unit length of the three-dimensional geological grid, in meters; For the The terrain slope of a 3D geological grid, in degrees.

[0041] This embodiment is another preferred implementation method. Based on the entity's three-dimensional geological model, it makes the calculation object more realistic and the calculation result more accurate and reliable without affecting the calculation time and complexity.

[0042] Example 5 The three-dimensional geological model of the slope and the stability calculation method based on airborne geophysical exploration include the following steps: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0043] In a, the surface elevation of the slope is obtained by drone aerial photography.

[0044] In a, the underground resistivity spatial data is obtained by semi-aerial electromagnetic method.

[0045] In b, the rock and soil characteristics include lithological information of the rock and soil strata, thickness of the overlying soil layer and moisture content of the soil.

[0046] In the above d, the anti-slip force is calculated by formula 2; Formula 2; in, For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; For the The internal friction angle of the structural surface of a three-dimensional geological grid, unit: °; For the The cohesion of the structural surface of a three-dimensional geological grid, in kPa; is the unit length of the three-dimensional geological grid, in meters; It is the unit width of the three-dimensional geological grid, in meters.

[0047] In the above d, the sliding force is calculated by formula 3; Formula 3; in, For the The sliding force caused by the gravity of the three-dimensional geological grid and other external forces, unit kN, For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN.

[0048] The deadweight of the geological grid is calculated by equation 4; Formula 4; in, For the The self-weight of a three-dimensional geological grid, in kN; For the The soil mass of a three-dimensional geological grid, in kN / m³; For the The center height of a three-dimensional geological grid, in meters; is the unit width of the three-dimensional geological grid, in meters; is the unit length of the three-dimensional geological grid, in meters; For the The terrain slope of a 3D geological grid, in degrees.

[0049] The horizontal additional force is calculated by formula 5; Formula 5; in, For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; is the horizontal additional force of the three-dimensional geological model of the slope, in kN; It is the angle between the inclination of the 3D geological grid and the horizontal additional force of the 3D geological model of the slope, in degrees.

[0050] This embodiment is the best implementation method, which helps to improve the accuracy of the three-dimensional stability calculation results of the slope by improving the authenticity of the three-dimensional geological model of the slope.

[0051] The present invention is described below with specific examples: See also Figure 1-Figure 4 Taking the hilly slope with an area of ​​about 4km² in the loess area as an example, the lithology from new to old is the light yellow loess-like sub-clay of the Holocene alluvial deposits, the Upper Pleistocene Malan loess, the mudstone of the lower part of the Neogene, sandy conglomerate and sandy mudstone.

[0052] The calculation method of the present invention is used for calculation, and the specific steps are as follows: S1. Obtain the spatial data of slope surface elevation and underground resistivity, where the geophysical data retains four columns of information: horizontal coordinate, vertical coordinate, elevation and resistivity value; S2. Select the controlling section in the area, use natural outcrops, boreholes and ground high-density resistivity method to obtain the characteristics of the rock and soil, extract the resistivity value corresponding to the characteristics of the rock and soil for calibration; For the calibration of resistivity value, the relationship between resistivity and saturation is obtained based on indoor tests, and the saturation is converted into the bulk density of a single three-dimensional geological grid. , obtain the spatial soil bulk density distribution, and then substitute it into each grid body for calculation; For sandstone or gravel soil layers with less clay content, The model is used to describe the relationship between the resistivity and saturation of the grid soil. For silty clay or clay layer containing clay particles, the Model is used to describe the relationship between grid soil resistivity and saturation; Model, namely Archie model: Formula 6; in, For the The saturation of a three-dimensional geological grid, unit: % For the Soil resistivity of a three-dimensional geological grid, in units ; is the water resistivity, unit ; is the soil porosity, unit: % is the soil cementation coefficient, dimensionless; is the soil property coefficient, dimensionless; is the soil saturation index, dimensionless; Model, namely the Waksman-Smits model: Formula 7; in, is the cation equivalent conductivity in the soil, in units of ; is the cation exchange capacity per unit pore volume in the soil, in units , reflecting the conductivity of clay; S3. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface. The 3D modeling software uses ArcGIS, Res2Dinv, Suffer or Python software; ArcGIS software is geographic information system software; Res2Dinv software is two-dimensional resistivity inversion software; Suffer software is scientific drawing software; Python is a programming language; S4, establishing a three-dimensional grid stability calculation model for the slope, and calculating the slope stability coefficient by formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity and other external forces of a three-dimensional geological grid, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

[0053] In this specific example, a three-dimensional grid model of the slope was established, with a terrain elevation accuracy of 0.2m, a resolution of 3.5m for the airborne geophysical points, and a grid width and length of 5m. The interface resistivity values ​​of the overlying loess and loess-like sub-clay and the underlying mudstone, sandy conglomerate and sandy mudstone were obtained by using natural outcrops, boreholes and ground high-density resistivity method, ranging from 200 to 400. , the groundwater level is not obvious; ArcGIS Pro, Res2Dinv, Suffer or Python are used to extract the soil thickness within the example range of 1.6-60m; according to indoor tests, the internal friction angle of the soil layer is 10°, the cohesion is 9kPa, and the horizontal additional force is 0.

Claims

1. A three-dimensional geological model of slopes and a stability calculation method based on aerial geophysical exploration, characterized in that: The following steps are involved: a. Obtaining the spatial data of slope surface elevation and underground resistivity; b. Obtain the characteristics of the rock and soil mass, extract the resistivity value corresponding to the characteristics of the rock and soil mass for calibration; c. Use spatial analysis 3D software to extract surface elevation and aerial geophysical resistivity spatial point data, extract corresponding interface points according to the calibrated range, construct the 3D structural surface of the slope from points to surfaces, and build a 3D geological model of the slope in combination with the ground surface; d. Establish a three-dimensional grid stability calculation model for the slope and calculate the slope stability coefficient using formula 1; Formula 1; in, For the Slope stability coefficient of a three-dimensional geological grid; For the The anti-sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN; For the The sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN.

2. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 1 is characterized in that: In a, the surface elevation of the slope is obtained by drone aerial photography.

3. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 1, characterized in that: In a, the underground resistivity spatial data is obtained by semi-aerial electromagnetic method.

4. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 1 is characterized in that: In b, the rock and soil characteristics include lithological information of the rock and soil strata, thickness of the overlying soil layer and moisture content of the soil.

5. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 1, characterized in that: In the above d, the anti-slip force is calculated by formula 2; Formula 2; in, For the The anti-sliding force caused by the gravity of a three-dimensional geological grid and other external forces, unit kN; For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; For the The internal friction angle of the structural surface of a three-dimensional geological grid, unit: °; For the The cohesion of the structural surface of a three-dimensional geological grid, in kPa; is the unit length of the three-dimensional geological grid, in meters; It is the unit width of the three-dimensional geological grid, in meters.

6. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 1, characterized in that: In the above d, the sliding force is calculated by formula 3; Formula 3; in, For the The sliding force caused by the gravity of the three-dimensional geological grid and other external forces, unit kN, For the The self-weight of a three-dimensional geological grid, in kN; For the The inclination angle of the structural surface of a three-dimensional geological grid, unit: °; For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN.

7. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 5 or 6, characterized in that: The deadweight of the geological grid is calculated by equation 4; Formula 4; in, For the The self-weight of a three-dimensional geological grid, in kN; For the The soil mass of a three-dimensional geological grid, in kN / m³; For the The center height of a three-dimensional geological grid, in meters; is the unit width of the three-dimensional geological grid, in meters; is the unit length of the three-dimensional geological grid, in meters; For the The terrain slope of a 3D geological grid, in degrees.

8. The method for establishing a three-dimensional geological model of a slope and calculating its stability based on aerial geophysical exploration according to claim 5 or 6, characterized in that: The horizontal additional force is calculated by formula 5; Formula 5; in, For the The horizontal additional force of a three-dimensional geological grid body on the slope, unit kN; is the horizontal additional force of the three-dimensional geological model of the slope, in kN; It is the angle between the inclination of the 3D geological grid and the horizontal additional force of the 3D geological model of the slope, in degrees.

Images