Slope sliding surface determination method, device and equipment and storage medium

Abstract

The invention discloses a method, a device and equipment for determining a slope sliding surface and a storage medium, and relates to the technical field of electrical digital data processing. Respectively drawing a projection vector line of each coordinate point on the monitoring section based on a geological drawing; respectively acquiring a vertical line of each projection vector line and acquiring all intersection points among all the vertical lines; deleting outliers in all the intersection points, and defining the reserved intersection points as effective intersection points; encircling all the effective intersection points into a surface through a triangulation algorithm; obtaining the circle center of the surface; calculating the distance from the circle center to the coordinate point of the deformation abrupt change inclinometer hole, and calculating the radius of the sliding surface according to the Euclidean distance; and drawing a slope sliding surface through the circle center and the radius. According to the method, all the vertical lines and the intersection points of all the vertical lines are comprehensively considered, and the outlier intersection points with too far relative distances are judged and deleted, so that the drawn sliding surface is closer to or coincides with the reality.

Description

Technical Field

[0001] The present application relates to the technical field of electronic digital data processing, and particularly relates to a method, device, equipment and storage medium for determining a landslide sliding surface. Background Art

[0002] Landslide prevention and control refers to the prevention and control measures taken when it is impossible to avoid landslide areas or unstable slope areas during engineering construction. Landslide stability analysis is the basic work of landslide prevention and control and the premise for subsequent prevention and control work. According to the characteristics of slope rock and soil mass, corresponding measures are taken to prevent slope instability to ensure the safety of production and personnel. Landslide stability analysis is the basic work of landslide prevention and control and the premise for subsequent prevention and control work. The significance of studying landslide stability lies in that it can not only provide a scientific theoretical basis for engineering construction, but also play an important guiding role in the early warning and prediction of the development trend of landslides.

[0003] Most of China's hydropower projects are concentrated in the southwestern region, where the geological conditions are complex. During the engineering construction and water storage process, landslide geological disasters often occur on the reservoir bank slopes, and the safety risk problem is prominent, bringing great potential safety hazards to the construction, operation and maintenance management of hydropower projects. Most traditional methods for searching the sliding surface of landslide geological disasters use the limit equilibrium method or the finite element method for analysis, but these methods are very sensitive to material parameters and solution methods, resulting in a certain error between the searched landslide sliding surface and the actual result. Summary of the Invention

[0004] The main purpose of the present application is to provide a method, device, equipment and storage medium for determining a landslide sliding surface to solve the problem that there is a certain error between the searched landslide sliding surface and the actual result in the prior art.

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

[0006] A method for determining a landslide sliding surface, the determination method is based on a plurality of inclinometer holes reserved on a monitoring section, and the determination method includes:

[0007] Step S1, respectively obtain the monitoring data of each inclinometer hole based on a plurality of natural moments, and respectively obtain the displacement vector of each monitoring data by the borehole inclinometry method;

[0008] Step S2, mark the inclinometer holes with displacement vectors exceeding a preset displacement threshold as deformed and mutated inclinometer holes, and respectively obtain the coordinate points of each deformed and mutated inclinometer hole;

[0009] Step S3, respectively output each coordinate point to the geological drawing of the monitoring section;

[0010] Step S4, based on the geological drawing, project vector lines are respectively drawn at each coordinate point on the monitoring section;

[0011] Step S5, perpendicular lines of each projection vector line are respectively obtained, and all intersection points between all perpendicular lines are obtained;

[0012] Step S6, based on a preset strategy, outlier points among all intersection points are deleted, and the remaining intersection points are defined as valid intersection points;

[0013] Step S7, all valid intersection points are enclosed into a surface through a triangulation algorithm;

[0014] Step S8, the centroid of the surface is obtained, and the centroid is the center of the circle of the landslide surface;

[0015] Step S9, the distances from the center of the circle to the coordinate points of the deformation mutation inclinometer holes are calculated respectively, and the radius of the sliding surface is calculated according to the Euclidean distance;

[0016] Step S10, the landslide surface is drawn through the center of the circle and the radius;

[0017] As a further improvement of the present application, step S6, based on a preset strategy, outlier points among all intersection points are deleted, and the remaining intersection points are defined as valid intersection points, including:

[0018] Step S61, define the intersection point dataset U of all intersection points as U = (P1, P2, …, P k , …, P m ), where m is the number of all intersection points;

[0019] Step S62, divide the intersection point dataset U horizontally and vertically. Based on the horizontal division, the intersection point abscissa dataset Ux = (P 1x , P 2x , …, P kx , …, P mx ) is obtained. Based on the vertical division, the intersection point ordinate dataset Uy = (P 1y , P 2y , …, P ky , …, P my ) is obtained;

[0020] Step S63, calculate the expectation μ x and standard deviation σ x of the intersection point abscissa dataset Ux, and the expectation μ y and standard deviation σ y of the intersection point ordinate dataset;

[0021] Step S64, based on the expectation μ x , standard deviation σ x , expectation μy , standard deviation σ y Determine whether the standardized variable of each intersection point is an outlier respectively;

[0022] Step S65, delete the intersection points determined to be outliers, and define the remaining intersection points as the effective intersection points.

[0023] As a further improvement of the present application, in step S64, based on the expected value μ x , standard deviation σ x , expected value μ y , standard deviation σ y Determine whether the standardized variable of each intersection point is an outlier respectively, including:

[0024] Step S641, when and , then determine that P k is an effective intersection point;

[0025] Step S642, when or , then determine that P k is an outlier.

[0026] As a further improvement of the present application, in step S8, obtain the centroid of the surface, and the centroid is the center of the circle of the slope sliding surface, including:

[0027] Step S81, define the center coordinates of the slope sliding surface as O(x o , y o );

[0028] Step S82, calculate the coordinate values of the center coordinates through formula (1):

[0029]

[0030] Among them, ∫∫xdA is the double integral of the product of the abscissa of all effective intersection points and the area element dA, ∫∫ydA is the double integral of the product of the ordinate of all effective intersection points and the area element dA, and ∫∫dA is the area of the surface.

[0031] As a further improvement of the present application, in step S9, calculate the distances from the center of the circle to the coordinate points of the deformation mutation inclinometer holes respectively, and calculate the sliding surface radius according to the Euclidean distance, including:

[0032] Step S91, locate the coordinate points of all deformation mutation inclinometer holes as the deformation data set C = (C1, C2,..., C n ), where n is the number of coordinate points of all deformation mutation inclinometer holes;

[0033] Step S92, calculate the center coordinates O(xo , y o ), the distances to each coordinate point in the deformed dataset

[0034] Step S93, arrange all the distances in ascending order to form a distance dataset D = (D1, D2, …, D n );

[0035] Step S94, sort the order of each coordinate point in the deformed dataset C according to the order of the corresponding distances in the distance dataset D to obtain the sorted deformed dataset C′ = (C1′, C2′, …, C n ′);

[0036] Step S95, based on the sorted deformed dataset C′ = (C1′, C2′, …, C n ′), calculate the radius of the slope sliding surface according to Equation (2) or Equation (3):

[0037] When the number of elements in D is odd:

[0038]

[0039] where x ( ′ n+1 ), y / 2 ′ ( ), n+1 ) / 2 are the abscissa and ordinate of C ( ′ n+1 ), / 2 respectively;

[0040] When the number of elements in D is even:

[0041]

[0042] where x′ n / 2 , y′ n / 2 are the abscissa and ordinate of C n ′ / 2 respectively;

[0043] As a further improvement of the present application, in step S10, draw the slope sliding surface through the center of the circle and the radius. After that, it includes:

[0044] Step S11, output the slope sliding surface on the geological drawing and send it to an external visualization terminal

[0045] To achieve the above object, the present application also provides the following technical solutions:

[0046] A device for determining the sliding surface of a slope body, which is applied to the above-mentioned determination method, and is characterized in that the device includes:

[0047] A monitoring data and displacement vector acquisition module, configured to respectively acquire the monitoring data of each inclinometer borehole based on a plurality of natural moments, and respectively acquire the displacement vector of each monitoring data by the borehole inclinometry method;

[0048] A deformed abrupt change inclinometer borehole coordinate point acquisition module, configured to mark the inclinometer borehole with a displacement vector exceeding a preset displacement threshold as a deformed abrupt change inclinometer borehole, and respectively acquire the coordinate points of each deformed abrupt change inclinometer borehole;

[0049] A coordinate point output module, configured to respectively output each coordinate point to the geological drawing of the monitoring section;

[0050] A projection vector line drawing module, configured to respectively draw the projection vector lines of each coordinate point on the monitoring section based on the geological drawing;

[0051] A perpendicular line intersection point acquisition module, configured to respectively acquire the perpendicular lines of each projection vector line and acquire all the intersection points between all the perpendicular lines;

[0052] An effective intersection point screening module, configured to delete the outlier points among all the intersection points based on a preset strategy, and define the remaining intersection points as effective intersection points;

[0053] An effective intersection point triangulation module, configured to enclose all the effective intersection points into a surface by a triangulation algorithm;

[0054] A sliding surface center acquisition module, configured to acquire the centroid of the surface, and the centroid is the center of the slope body sliding surface;

[0055] A sliding surface radius acquisition module, configured to calculate the distances from the center to the coordinate points of the deformed abrupt change inclinometer boreholes respectively, and calculate the sliding surface radius according to the Euclidean distance;

[0056] A slope body sliding surface drawing module, configured to draw the slope body sliding surface through the center and the radius.

[0057] To achieve the above object, the present application also provides the following technical solutions:

[0058] An electronic device, including a processor and a memory coupled to the processor, where the memory stores program instructions executable by the processor; when the processor executes the program instructions stored in the memory, the above-mentioned method for determining the slope body sliding surface is implemented.

[0059] To achieve the above object, the present application also provides the following technical solutions:

[0060] A storage medium stores program instructions, and when the program instructions are executed by a processor, the method for determining the landslide sliding surface as described above can be implemented.

[0061] In this application, the monitoring data of each inclinometer hole at several natural moments are obtained respectively, and the displacement vectors of each monitoring data are obtained respectively by the borehole inclinometry method; the inclinometer holes with displacement vectors exceeding the preset displacement threshold are marked as deformed and mutated inclinometer holes, and the coordinate points of each deformed and mutated inclinometer hole are obtained respectively; each coordinate point is output to the geological drawing of the monitoring section respectively; the projection vector lines of each coordinate point are drawn on the monitoring section respectively based on the geological drawing; the perpendicular lines of each projection vector line are obtained respectively and all the intersection points between all the perpendicular lines are obtained; the outlier points among all the intersection points are deleted based on the preset strategy, and the remaining intersection points are defined as valid intersection points; all the valid intersection points are enclosed into a surface by the triangulation algorithm; the centroid of the surface is obtained, and the centroid is the center of the landslide sliding surface; the distances from the center to the coordinate points of the deformed and mutated inclinometer holes are calculated respectively, and the sliding surface radius is calculated according to the Euclidean distance; the landslide sliding surface is drawn through the center and the radius. This application comprehensively considers all the perpendicular lines and the intersection points of all the perpendicular lines, determines and deletes the outlier intersection points with relatively large distances, so that when drawing the sliding surface, every mutation point can be taken into account as much as possible and the mutation points with greater influence can be excluded. Compared with the traditional strength reduction method and limit equilibrium method, the sliding surface determined by this application is closer to or coincides with the actual situation. BRIEF DESCRIPTION OF THE DRAWINGS

[0062] Figure 1 It is a schematic flowchart of the steps of an embodiment of the method for determining the landslide sliding surface of this application;

[0063] Figure 2 It is a schematic diagram of the functional modules of an embodiment of the device for determining the landslide sliding surface of this application;

[0064] Figure 3 It is a schematic structural diagram of an embodiment of the electronic device of this application;

[0065] Figure 4 It is a schematic structural diagram of an embodiment of the storage medium of this application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0066] Next, the technical solutions in the embodiments of this application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of this application. Obviously, the described embodiments are only a part of the embodiments of this application, rather than all the embodiments. Based on the embodiments in this application, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of this application.

[0067] The terms "first", "second", and "third" in this application are only used for descriptive purposes and should not be construed as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, features defined with "first", "second", and "third" may explicitly or implicitly include at least one such feature. In the description of this application, the meaning of "a plurality" is at least two, such as two, three, etc., unless otherwise specifically and clearly defined. All directional indications (such as up, down, left, right, front, back...) in the embodiments of this application are only used to explain the relative positional relationship, movement conditions, etc. between components in a specific posture (as shown in the drawings). If the specific posture changes, then the directional indications will also change accordingly. In addition, the terms "comprise" and "have" and any variations thereof are intended to cover non-exclusive inclusion. For example, a process, method, system, product, or device that includes a series of steps or units is not limited to the listed steps or units, but optionally also includes steps or units not listed, or optionally also includes other steps or units inherent to these processes, methods, products, or devices.

[0068] Reference to "embodiments" in this context means that a particular feature, structure, or characteristic described in connection with the embodiments can be included in at least one embodiment of the application. The phrase appears in various places in the specification and does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment mutually exclusive with other embodiments. Those skilled in the art will explicitly and implicitly understand that the embodiments described herein can be combined with other embodiments.

[0069] As Figure 1 shown, this embodiment provides an embodiment of a method for determining the landslide surface. In this embodiment, the determination method is based on a number of inclinometer holes reserved on the monitoring section.

[0070] Preferably, an inclinometer hole is a method for measuring the deformation of the ground soil layer, determining the foundation deformation and the horizontal shape of the stratum, and is usually used for monitoring building structures, subways, roads, bridges, and other infrastructure projects, and can also be applied in fields such as groundwater hydrogeological exploration. The principle of an inclinometer borehole is to monitor and record the deformation of the stratum by setting measuring instruments such as inclinometers and deformeters in the borehole, so as to analyze the deformation characteristics of the soil mass and further evaluate the foundation stability of the project.

[0071] Preferably, inclinometer holes are usually realized by a borehole inclinometer:

[0072] 1) Install a borehole inclinometer.

[0073] ① Select a suitable borehole diameter to ensure that the borehole diameter can accommodate the inclinometer and its cable.

[0074] ②Install the inclinometer in the borehole pipe and ensure that the inclinometer can work properly.

[0075] ③Insert the cable and connect the data collector.

[0076] 2) Calibrate the borehole inclinometer.

[0077] ①Install the level and adjust the horizontal position of the borehole inclinometer.

[0078] ②Perform zero calibration so that the output of the borehole inclinometer is zero when it is in the horizontal position.

[0079] ③Calibrate the angle, adjust according to the reading difference, and ensure that the measured angle data is accurate.

[0080] 3) Conduct the measurement.

[0081] ①Turn on the data collector and start recording the angle data of the borehole inclinometer.

[0082] ②Record at different depth intervals and ensure that the depth intervals are consistent.

[0083] ③After the measurement is completed, turn off the data collector, export the data and process it.

[0084] 4) Data processing.

[0085] ①Import the data into professional data processing software.

[0086] ②Conduct data analysis and generate a measurement report.

[0087] ③Preliminarily exclude the data that does not meet the requirements and retain the accurate and reliable data.

[0088] Specifically, the determination method includes the following steps:

[0089] Step S1, respectively obtain the monitoring data of each inclinometer borehole based on several natural moments, and respectively obtain the displacement vector of each monitoring data by the borehole inclinometry method.

[0090] Preferably, the borehole inclinometry method can measure the vertical displacement and horizontal displacement of the soil mass, and then convert the vertical and horizontal displacement lengths into vectors.

[0091] Step S2, mark the inclinometer boreholes with displacement vectors exceeding the preset displacement threshold as deformed and mutated inclinometer boreholes, and respectively obtain the coordinate points of each deformed and mutated inclinometer borehole.

[0092] Preferably, the preset displacement threshold can be set to 10 mm, that is, the length of the historical displacement vector does not exceed 10 mm.

[0093] Step S3, respectively output each coordinate point to the geological drawing of the monitoring section.

[0094] Preferably, the geological drawing is a slope profile drawing generated during preliminary exploration.

[0095] Step S4: Based on the geological drawing, draw the projection vector lines of each coordinate point on the monitoring section respectively.

[0096] Step S5: Obtain the perpendicular lines of each projection vector line respectively and obtain all the intersection points between all the perpendicular lines.

[0097] Step S6: Based on a preset strategy, delete the outlier points among all the intersection points, and define the remaining intersection points as valid intersection points.

[0098] Step S7: Enclose all the valid intersection points into a surface through the triangulation algorithm.

[0099] Step S8: Obtain the centroid of the surface, and the centroid is the center of the circle of the slope sliding surface.

[0100] It can be understood that the centroid is the center of a shape.

[0101] Step S9: Calculate the distances from the center of the circle to the coordinate points of the deformation mutation inclinometer holes respectively, and calculate the radius of the sliding surface according to the Euclidean distance.

[0102] Step S10: Draw the slope sliding surface through the center of the circle and the radius.

[0103] Furthermore, in step S6, based on a preset strategy, delete the outlier points among all the intersection points, and define the remaining intersection points as valid intersection points, including:

[0104] Step S61: Define the intersection point data set U of all the intersection points as U=(P1, P2, …, P k , …, P m ), where m is the number of all the intersection points.

[0105] Step S62: Divide the intersection point data set U horizontally and vertically. Based on the horizontal division, obtain the intersection point abscissa data set Ux=(P 1x , P 2x , …, P kx , …, P mx ), and based on the vertical division, obtain the intersection point ordinate data set Uy=(P 1y , P 2y , …, P ky , …, P my ).

[0106] Step S63: Calculate the expectation μ x and the standard deviation σ x of the intersection point abscissa data set Ux, and the expectation μ y, standard deviation σ y .

[0107] Step S64, based on the expectation μ x , standard deviation σ x , expectation μ y , standard deviation σ y Judge whether the standardized variable of each intersection point is an outlier respectively.

[0108] Step S65, delete the intersection points judged as outliers, and define the remaining intersection points as valid intersection points.

[0109] Furthermore, in step S64, based on the expectation μ x , standard deviation σ x , expectation μ y , standard deviation σ y Judge whether the standardized variable of each intersection point is an outlier respectively, including:

[0110] Step S641, when and , then determine that P k is a valid intersection point.

[0111] Step S642, when or , then determine that P k is an outlier.

[0112] Furthermore, in step S8, obtain the centroid of the surface, and the centroid is the center of the circle of the slope sliding surface, including:

[0113] Step S81, define the center coordinates of the slope sliding surface as O(x o , y o ).

[0114] Step S82, calculate the coordinate values of the center coordinates through formula (1):

[0115]

[0116] Among them, ∫∫xdA is the double integral of the product of the abscissa of all valid intersection points and the area element dA, ∫∫ydA is the double integral of the product of the ordinate of all valid intersection points and the area element dA, and ∫∫dA is the area of the surface.

[0117] Furthermore, in step S9, calculate the distances from the center of the circle to the coordinate points of the deformation mutation inclinometer holes respectively, and calculate the sliding surface radius according to the Euclidean distance, including:

[0118] Step S91, locate the coordinate points of all deformation mutation inclinometer holes as the deformation data set C = (C1, C2,..., C n), where n is the number of coordinate points of all deformed and mutated inclinometer holes.

[0119] Step S92, calculate the distances from the center coordinates O(x o , y o ) to each coordinate point in the deformation dataset.

[0120] Step S93, arrange all the distances in ascending order to form a distance dataset D = (D1, D2, …, D n ).

[0121] Step S94, sort the order of each coordinate point in the deformation dataset C according to the order of the corresponding distances in the distance dataset D to obtain the sorted deformation dataset C′ = (C1′, C2′, …, C n ′).

[0122] Step S95, based on the sorted deformation dataset C′ = (C1′, C2′, …, C n ′), calculate the radius of the slope sliding surface according to Equation (2) or Equation (3):

[0123] When the number of elements in D is odd:

[0124]

[0125] where x ( ′ n+1 ), / 2 y ( ′ n+1 ), / 2 are the abscissa and ordinate of C ( ′ n+1 ), / 2 in sequence.

[0126] When the number of elements in D is even:

[0127]

[0128] where x′ n / 2 , y′ n / 2 are the abscissa and ordinate of C n ′ / 2 in sequence.

[0129] Furthermore, in step S10, draw the slope sliding surface through the center and the radius. After that, it includes:

[0130] Step S11: Output the slope sliding surface on the geological drawing respectively and send it to an external visualization terminal. In this embodiment, by respectively obtaining the monitoring data of each inclinometer borehole based on several natural moments, and respectively obtaining the displacement vectors of each monitoring data through the borehole inclinometry method; marking the inclinometer boreholes with displacement vectors exceeding the preset displacement threshold as deformed and mutated inclinometer boreholes, and respectively obtaining the coordinate points of each deformed and mutated inclinometer borehole; respectively outputting each coordinate point to the geological drawing of the monitoring section; respectively drawing the projection vector lines of each coordinate point on the monitoring section; respectively obtaining the perpendicular lines of each projection vector line and obtaining all the intersection points between all the perpendicular lines; deleting the outlier points among all the intersection points based on a preset strategy, and defining the remaining intersection points as valid intersection points; enclosing all the valid intersection points into a surface through the triangulation algorithm; obtaining the centroid of the surface, and the centroid is the center of the circle of the slope sliding surface; calculating the distances from the center of the circle to the coordinate points of the deformed and mutated inclinometer boreholes respectively, and calculating the sliding surface radius according to the Euclidean distance; drawing the slope sliding surface through the center of the circle and the radius. This embodiment comprehensively considers all the perpendicular lines and the intersection points of all the perpendicular lines, determines and deletes the outlier intersection points with relatively large distances, so that when drawing the sliding surface, every mutation point can be taken into account as much as possible and the mutation points with greater influence can be excluded. Compared with the traditional strength reduction method and limit equilibrium method, the sliding surface determined in this embodiment is closer to or coincides with the actual situation.

[0131] As Figure 2 shown, this embodiment provides an embodiment of the device for determining the slope sliding surface. In this embodiment, the determining device is applied to the determining method in the above-mentioned embodiment.

[0132] Specifically, the determining device includes a monitoring data and displacement vector obtaining module 1, a deformed and mutated inclinometer borehole coordinate point obtaining module 2, a coordinate point output module 3, a projection vector line drawing module 4, a perpendicular line intersection point obtaining module 5, a valid intersection point screening module 6, a valid intersection point triangulation module 7, a slope sliding surface center obtaining module 8, a slope sliding surface radius obtaining module 9, and a slope sliding surface drawing module 10 that are electrically connected in sequence.

[0133] Among them, the monitoring data and displacement vector acquisition module 1 is used to respectively acquire the monitoring data of each inclinometer borehole based on several natural moments, and respectively acquire the displacement vectors of each monitoring data through the borehole inclinometry method; the deformed and mutated inclinometer borehole coordinate point acquisition module 2 is used to mark the inclinometer boreholes with displacement vectors exceeding the preset displacement threshold as deformed and mutated inclinometer boreholes, and respectively acquire the coordinate points of each deformed and mutated inclinometer borehole; the coordinate point output module 3 is used to respectively output each coordinate point to the geological drawing of the monitoring section; the projection vector line drawing module 4 is used to respectively draw the projection vector lines of each coordinate point on the monitoring section based on the geological drawing; the perpendicular intersection point acquisition module 5 is used to respectively acquire the perpendiculars of each projection vector line and acquire all the intersection points between all the perpendiculars; the valid intersection point screening module 6 is used to delete the outliers among all the intersection points based on a preset strategy, and define the remaining intersection points as valid intersection points; the valid intersection point triangulation module 7 is used to enclose all the valid intersection points into a surface through a triangulation algorithm; the sliding surface center acquisition module 8 is used to acquire the centroid of the surface, and the centroid is the center of the slope sliding surface; the sliding surface radius acquisition module 9 is used to calculate the distances from the center to the coordinate points of the deformed and mutated inclinometer boreholes respectively, and calculate the sliding surface radius according to the Euclidean distance; the slope sliding surface drawing module 10 is used to draw the slope sliding surface through the center and the radius.

[0134] Further, the valid intersection point screening module 6 specifically includes a first valid intersection point screening sub-module, a second valid intersection point screening sub-module, a third valid intersection point screening sub-module, a fourth valid intersection point screening sub-module, and a fifth valid intersection point screening sub-module that are electrically connected in sequence; the first valid intersection point screening sub-module is electrically connected to the perpendicular intersection point acquisition module 5, and the fifth valid intersection point screening sub-module is electrically connected to the valid intersection point triangulation module 7.

[0135] Among them, the first valid intersection point screening sub-module is used to define the intersection point data set U of all the intersection points as U=(P1, P2, …, P k , …, P m ), where m is the number of all the intersection points; the second valid intersection point screening sub-module is used to divide the intersection point data set U horizontally and vertically, and based on the horizontal division, obtain the intersection point abscissa data set Ux=(P 1x , P 2x , …, P kx , …, P mx ), and based on the vertical division, obtain the intersection point ordinate data set Uy=(P 1y , P 2y , …, P ky , …, P my ); the third valid intersection point screening sub-module is used to calculate the expectation μ x of the intersection point abscissa data set Ux, the standard deviation σ x , and the expectation μ of the intersection point ordinate data sety , standard deviation σ y ; The fourth effective intersection point screening sub-module is used to based on the expected value μ x , standard deviation σ x , expected value μ y , standard deviation σ y respectively determine whether the standardized variable of each intersection point is an outlier; The fifth effective intersection point screening sub-module is used to delete the intersection points determined as outliers and define the remaining intersection points as effective intersection points.

[0136] Furthermore, the fourth effective intersection point screening sub-module specifically includes a first effective intersection point screening unit and a second effective intersection point screening unit that are electrically connected in sequence; The first effective intersection point screening unit is electrically connected to the third effective intersection point screening sub-module, and the second effective intersection point screening unit is electrically connected to the fifth effective intersection point screening sub-module.

[0137] Among them, the first effective intersection point screening unit is used when and , then determine that P k is an effective intersection point; The second effective intersection point screening unit is used when or , then determine that P k is an outlier.

[0138] Furthermore, the sliding surface center obtaining module 8 specifically includes a first sliding surface center obtaining sub-module and a second sliding surface center obtaining sub-module that are electrically connected in sequence. The first sliding surface center obtaining sub-module is electrically connected to the effective intersection point triangulation module 7, and the second sliding surface center obtaining sub-module is electrically connected to the sliding surface radius obtaining module 9.

[0139] Among them, the first sliding surface center obtaining sub-module is used to define the center coordinates of the slope sliding surface as O(x o , y o ).

[0140] The second sliding surface center obtaining sub-module is used to calculate the coordinate value of the center coordinates through Equation (1):

[0141]

[0142] Among them, ∫∫xdA is the double integral of the product of the abscissa of all effective intersection points and the area element dA, ∫∫ydA is the double integral of the product of the ordinate of all effective intersection points and the area element dA, and ∫∫dA is the area of the surface.

[0143] Further, the sliding surface radius acquisition module 9 specifically includes a first sliding surface radius acquisition sub-module, a second sliding surface radius acquisition sub-module, a third sliding surface radius acquisition sub-module, a fourth sliding surface radius acquisition sub-module, and a fifth sliding surface radius acquisition sub-module that are electrically connected in sequence; the first sliding surface radius acquisition sub-module is electrically connected to the second sliding surface center acquisition sub-module, and the fifth sliding surface radius acquisition sub-module is electrically connected to the slope body sliding surface drawing module 10.

[0144] Among them, the first sliding surface radius acquisition sub-module is used to locate the coordinate points of all deformed mutation inclinometer holes as the deformation data set C = (C1, C2, …, C n ), where n is the number of coordinate points of all deformed mutation inclinometer holes; the second sliding surface radius acquisition sub-module is used to calculate the distances from the center coordinates O(x o , y o ) to each coordinate point in the deformation data set; the third sliding surface radius acquisition sub-module is used to arrange all the distances in ascending order to form the distance data set D = (D1, D2, …, D n ); the fourth sliding surface radius acquisition sub-module is used to sort the order of each coordinate point in the deformation data set C according to the order of the corresponding distances in the distance data set D to obtain the sorted deformation data set C′ = (C1′, C2′, …, C n ′); the fifth sliding surface radius acquisition sub-module is used to calculate the radius of the slope body sliding surface based on the sorted deformation data set C′ = (C1′, C2′, …, C n ′) according to Equation (2) or Equation (3):

[0145] When the number of elements in D is odd:

[0146]

[0147] Among them, x ( ′ n+1 ) / 2 , y ( ′ n+1 ) / 2 are the abscissa and ordinate of C ( ′ n+1 ) / 2 in sequence.

[0148] When the number of elements in D is even:

[0149]

[0150] Among them, x′ n / 2 , y′ n / 2 are the abscissa and ordinate of C n ′ / 2 in sequence.

[0151] Furthermore, the determination device further includes a slope sliding surface sending module electrically connected to the slope sliding surface drawing module 10.

[0152] The slope sliding surface sending module is configured to respectively output the slope sliding surface on a geological drawing and send it to an external visualization terminal.

[0153] It should be noted that this embodiment is a functional module embodiment based on the above method embodiment. For the preferred, extended, limited, exemplified, and principle description parts of this embodiment, refer to the above embodiment, and this embodiment will not be elaborated herein.

[0154] In this embodiment, by respectively obtaining the monitoring data of each inclinometer borehole based on a number of natural moments, and respectively obtaining the displacement vectors of each monitoring data through the borehole inclinometry method; marking the inclinometer boreholes with displacement vectors exceeding the preset displacement threshold as deformed mutation inclinometer boreholes, and respectively obtaining the coordinate points of each deformed mutation inclinometer borehole; respectively outputting each coordinate point to the geological drawing of the monitoring section; respectively drawing the projection vector lines of each coordinate point on the monitoring section; respectively obtaining the perpendicular lines of each projection vector line and obtaining all the intersection points between all the perpendicular lines; deleting the outlier points among all the intersection points based on a preset strategy, and defining the remaining intersection points as valid intersection points; enclosing all the valid intersection points into a surface through a triangulation algorithm; obtaining the centroid of the surface, and the centroid is the center of the circle of the slope sliding surface; calculating the distances from the center of the circle to the coordinate points of the deformed mutation inclinometer boreholes respectively, and calculating the sliding surface radius according to the Euclidean distance; drawing the slope sliding surface through the center of the circle and the radius. This embodiment comprehensively considers all the perpendicular lines and the intersection points of all the perpendicular lines, determines and deletes the outlier intersection points with relatively large distances, so that when drawing the sliding surface, it can take into account every mutation point as much as possible and exclude the mutation points with greater influence. Compared with the traditional strength reduction method and limit equilibrium method, the sliding surface determined in this embodiment is closer to or coincides with the actual situation.

[0155] Figure 3 It is a schematic structural diagram of an electronic device according to an embodiment of the present application. As Figure 3 shown, the electronic device 11 includes a processor 111 and a memory 112 coupled to the processor 111.

[0156] The memory 112 stores program instructions for implementing the method for determining the slope sliding surface according to any of the above embodiments.

[0157] The processor 111 is configured to execute the program instructions stored in the memory 112 to determine the slope sliding surface.

[0158] Among them, the processor 111 can also be referred to as a CPU (Central Processing Unit). The processor 111 may be an integrated circuit chip with signal processing capabilities. The processor 111 can also be a general-purpose processor, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components. The general-purpose processor can be a microprocessor or the processor can also be any conventional processor, etc.

[0159] Furthermore, Figure 4 is a schematic structural diagram of a storage medium according to an embodiment of the present application. Refer to Figure 4 , the storage medium 12 of the embodiment of the present application stores program instructions 121 that can implement all of the above methods. Among them, the program instructions 121 can be stored in the above storage medium in the form of a software product, including several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) or a processor to execute all or part of the steps of the methods described in various embodiments of the present application. The foregoing storage medium includes: various media that can store program codes such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs, or terminal devices such as computers, servers, mobile phones, and tablets.

[0160] In several embodiments provided by the present application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of units is only a logical function division, and there can be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed mutual coupling or direct coupling or communication connection can be through some interfaces, and the indirect coupling or communication connection of devices or units can be in an electrical, mechanical, or other form.

[0161] In addition, in each embodiment of the present application, each functional unit may be integrated into one processing unit, may exist independently as individual units physically, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of a software functional unit. The above is only the implementation manner of the present application, and does not limit the patent scope of the present application. Any equivalent structure or equivalent process transformation made by using the content of the specification and drawings of the present application, or directly or indirectly applied in other related technical fields, shall be included in the patent protection scope of the present application by the same token.

[0162] The specific implementation manners of the present application have been described in detail above, but they are only examples, and the present application is not limited to the specific implementation manners described above. For those skilled in the art, any equivalent modification or substitution of the present application is also within the scope of the present application. Therefore, all equivalent transformations, modifications, improvements, etc. made without departing from the spirit and principles of the present application should be covered by the scope of the present application.

Claims

1. A method for determining a slope sliding surface, the method being based on a number of inclinometer holes reserved on a monitoring section, characterized in that: The determination method comprises: Step S1, respectively obtaining monitoring data of each inclinometer hole based on a number of natural moments, and respectively obtaining a displacement vector of each monitoring data by a borehole inclinometer method; Step S2, marking the inclinometer holes whose displacement vectors exceed a preset displacement threshold as deformation mutation inclinometer holes, and obtaining the coordinate points of each deformation mutation inclinometer hole respectively; Step S3, outputting each coordinate point to the geological drawing of the monitoring section respectively; Step S4, drawing the projection vector line of each coordinate point on the monitoring section based on the geological drawing; Step S5, respectively obtaining the perpendicular line of each projection vector line and obtaining all the intersection points between all the perpendicular lines; Step S6, deleting outlier points from all intersection points based on a preset strategy, and defining the remaining intersection points as valid intersection points; Step S7, all valid intersection points are enclosed into a surface through a triangulation algorithm; Step S8, obtaining the centroid of the surface, which is the center of the slope sliding surface; Step S9, calculating the distances from the center of the circle to the coordinate points of the deformation mutation inclination hole, and calculating the sliding surface radius according to the Euclidean distance; Step S10: draw the slope sliding surface through the center of the circle and the radius.

2. The determination method according to claim 1, characterized in that: Step S6, based on a preset strategy, delete outliers from all intersection points, and define the remaining intersection points as valid intersection points, including: Step S61, define the intersection point data set U = (P1, P2, ..., P k ,…,P m ), where m is the number of all intersection points; Step S62: divide the intersection point data set U into horizontal and vertical directions, and obtain the intersection point abscissa data set Ux=(P 1x ,P 2x ,…,P kx ,…,P mx ), based on the vertical division, the intersection ordinate data set Uy=(P 1y ,P 2y ,…,P ky ,…,P my ); Step S63, calculate the expected μ of the intersection point horizontal coordinate data set Ux x , standard deviation σ x , and the expected μ of the intersection ordinate data set y , standard deviation σ y ; Step S64, based on the expected μ x , standard deviation σ x , expected μ y , standard deviation σ y Determine whether the standardized variable of each intersection point is an outlier; Step S65, deleting the intersection points determined to be outliers, and defining the remaining intersection points as the valid intersection points.

3. The determination method according to claim 2, characterized in that: Step S64, based on the expected μ x , standard deviation σ x , expected μ y , standard deviation σ y Determine whether the standardized variable of each intersection point is an outlier, including: Step S641, when and When P k is a valid intersection point; Step S642, when or When P k For outliers.

4. The determination method according to claim 1, characterized in that: Step S8, obtaining the centroid of the surface, which is the center of the slope sliding surface, includes: Step S81, define the center coordinates of the slope sliding surface as O(x o ,y o ); Step S82, calculating the coordinate value of the center coordinate of the circle by formula (1): Among them, ∫∫xdA is the double integral of the product of the horizontal coordinates of all valid intersection points and the area element dA, ∫∫ydA is the double integral of the product of the vertical coordinates of all valid intersection points and the area element dA, and ∫∫dA is the area of ​​the surface.

5. The determination method according to claim 4, characterized in that: Step S9, calculating the distances from the center of the circle to the coordinate points of the deformation mutation inclination hole, and calculating the sliding surface radius according to the Euclidean distance, including: Step S91: locate the coordinate points of all deformation mutation inclination holes as a deformation data set C = (C1, C2, ..., C n ), where n is the number of coordinate points of all deformation mutation inclinometer holes; Step S92, calculate the center coordinates O(x o ,y o ) to each coordinate point in the deformation data set; Step S93: Arrange all distances in ascending order to form a distance data set D = (D1, D2, ..., D n ); Step S94: sort the order of each coordinate point in the deformed data set C according to the order of the corresponding distances in the distance data set D, and obtain the sorted deformed data set C′=(C1′, C2′, …, C n ′); Step S95: based on the sorted deformed data set C′=(C1′, C2′, …, C n ') Calculate the radius of the slope sliding surface according to formula (2) or formula (3): When the number of elements of D is odd: Among them, x ( ' n+1 ) / 2 ,y ( ' n+1 ) / 2 C ( ' n+1 ) / 2 The horizontal and vertical coordinates of When the number of elements of D is even: Among them, x′ n / 2 , y′ n / 2 C n ' / 2 The horizontal and vertical coordinates of .

6. The determination method according to claim 1, characterized in that: Step S10, drawing the slope sliding surface through the center of the circle and the radius, and then comprising: Step S11, outputting the slope sliding surface on the geological drawing respectively, and sending it to an external visualization terminal.

7. A device for determining a slope sliding surface, the device being applied to the determination method according to any one of claims 1 to 6, characterized in that: The determining device comprises: The monitoring data and displacement vector acquisition module is used to respectively acquire the monitoring data of each inclinometer hole based on a number of natural moments, and respectively acquire the displacement vector of each monitoring data by means of the borehole inclinometer method; A deformation mutation inclinometer hole coordinate point acquisition module is used to mark the inclinometer holes whose displacement vectors exceed a preset displacement threshold as deformation mutation inclinometer holes, and to respectively acquire the coordinate points of each deformation mutation inclinometer hole; A coordinate point output module, used for outputting each coordinate point to the geological drawing of the monitoring section; A projection vector line drawing module, used for drawing the projection vector line of each coordinate point on the monitoring section based on the geological drawing; A vertical line intersection point acquisition module is used to respectively acquire the vertical line of each projection vector line and acquire all the intersection points between all the vertical lines; The effective intersection point screening module is used to delete outliers from all intersection points based on a preset strategy and define the retained intersection points as effective intersection points; The effective intersection triangulation module is used to enclose all effective intersection points into a surface through the triangulation algorithm; A sliding surface center acquisition module is used to acquire the centroid of the surface, which is the center of the slope sliding surface; A sliding surface radius acquisition module is used to calculate the distances from the center of the circle to the coordinate points of the deformation mutation inclination hole, and calculate the sliding surface radius according to the Euclidean distance; The slope sliding surface drawing module is used to draw the slope sliding surface through the center of the circle and the radius.

8. An electronic device, characterized in that: It comprises a processor and a memory coupled to the processor, wherein the memory stores program instructions executable by the processor; when the processor executes the program instructions stored in the memory, the determination method as described in any one of claims 1 to 6 is implemented.

9. A storage medium, characterized in that: The storage medium stores program instructions, and when the program instructions are executed by the processor, the determination method according to any one of claims 1 to 6 can be implemented.

Images