Method and device for generating central axis of drainage groove of barrier lake

Abstract

The invention provides a method and device for generating a central axis of a drainage groove of a barrier lake, and the method comprises the steps: determining a main water flow path point of a barrier body from a plurality of grid units corresponding to the barrier body based on the topographic data of the barrier body corresponding to the barrier lake; wherein the grid units are obtained by dividing a region corresponding to the weir plug body; determining a stability index value of each grid unit according to monitoring data of a plurality of monitoring points on the weir body; and according to the main water flow path point of the barrier body and the stability index value of each grid unit, generating the central axis of the drainage groove of the barrier lake in the geologic model corresponding to the barrier body. Therefore, the central axis of the drainage groove of the barrier lake can be automatically and efficiently generated based on the main water flow path point and the stability index value of the barrier body, and powerful support is provided for emergency management or hydraulic engineering design of barrier lake disasters.

Description

Technical Field

[0001] The present application relates to the field of computer technology, and in particular to a method and device for generating a central axis of a drainage channel of a barrier lake. Background Art

[0002] A barrier lake is often formed by a river being blocked by natural disasters such as landslides and debris flows. The rapid increase in its water storage capacity may cause secondary disasters such as dam breaches and flooding, posing a huge threat to the safety of life and property in the downstream area. In the process of emergency rescue, how to quickly and safely reduce the water level of the barrier lake is the key to controlling the danger.

[0003] In the relevant technology, in the emergency rescue plan for the landslide lake, the method of digging drainage channels is usually adopted. However, this process often relies on manual experience or simple terrain analysis to determine the location, depth and width of the drainage channels. It lacks accurate decision-making support and is prone to problems such as unreasonable paths and excessive excavation. On the one hand, it increases the difficulty of construction. On the other hand, excavation in the unstable area of ​​the landslide body will increase the risk of secondary disasters such as dam breach. Summary of the invention

[0004] The present application provides a method and device for generating the central axis of a drainage channel of a barrier lake, so as to solve one of the technical problems in the related art at least to a certain extent. The technical solution of the present application is as follows:

[0005] According to one aspect of the present application, a method for generating a central axis of a drainage channel of a barrier lake is provided, comprising:

[0006] Based on the topographic data of the barrier body corresponding to the barrier lake, determining the main water flow path point of the barrier body from a plurality of grid units corresponding to the barrier body; wherein the grid units are obtained by dividing the area corresponding to the barrier body;

[0007] Determining the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body;

[0008] According to the main water flow path points of the barrier body and the stability index values ​​of each grid unit, the central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body.

[0009] According to another aspect of the present application, another device for generating a central axis of a drainage channel of a barrier lake is provided, comprising:

[0010] A first determination module is used to determine the main water flow path point of the damming body from a plurality of grid cells corresponding to the damming body based on the topographic data of the damming body corresponding to the damming lake; wherein the grid cells are obtained by dividing the area corresponding to the damming body;

[0011] A second determination module is used to determine the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body;

[0012] A generation module is used to generate the central axis of the drainage channel of the barrier lake in the geological model corresponding to the barrier body according to the main water flow path points of the barrier body and the stability index value of each grid unit.

[0013] According to another aspect of the present application, an electronic device is provided, characterized in that it includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and when the processor executes the program, the method for generating the central axis of the drainage channel of the barrier lake proposed in the above aspect of the present application is implemented.

[0014] According to another aspect of the present application, a non-transitory computer-readable storage medium of computer instructions is provided, wherein the computer instructions are used to enable the computer to execute the method for generating the central axis of the drainage channel of the barrier lake proposed in the above-mentioned aspect of the present application.

[0015] According to another aspect of the present application, a computer program product is provided, including a computer program, which, when executed by a processor, implements the method for generating the central axis of a drainage channel of a barrier lake proposed in the above aspect of the present application.

[0016] The technical solution provided by the embodiments of the present application brings at least the following beneficial effects:

[0017] Based on the topographic data of the dammed body corresponding to the dammed lake, the main water flow path point of the dammed body is determined from multiple grid units corresponding to the dammed body; wherein the grid unit is obtained by dividing the area corresponding to the dammed body; according to the monitoring data of multiple monitoring points on the dammed body, the stability index value of each grid unit is determined; according to the main water flow path point of the dammed body and the stability index value of each grid unit, the central axis of the drainage channel of the dammed lake is generated in the geological model corresponding to the dammed body. In this way, it is possible to automatically, efficiently and accurately generate the central axis of the drainage channel of the dammed lake based on the effectively determined main water flow path point and stability index value of the dammed body, providing strong support for emergency management of dammed lake disasters or water conservancy project design.

[0018] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0020] Figure 1A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in an embodiment of the present application;

[0021] Figure 2 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application;

[0022] Figure 3 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application;

[0023] Figure 4 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application;

[0024] Figure 5 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application;

[0025] Figure 6 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application;

[0026] Figure 7 A schematic diagram of a grid unit corresponding to the dam body provided in this application;

[0027] Figure 8 A schematic structural diagram of a device for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application;

[0028] Fig. 9 A block diagram of an exemplary electronic device suitable for implementing embodiments of the present invention is shown. DETAILED DESCRIPTION

[0029] Embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0030] The following describes the method, device, equipment and medium for generating the central axis of the drainage channel of the barrier lake in an embodiment of the present application with reference to the accompanying drawings.

[0031] Figure 1 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in one embodiment of the present application.

[0032] The embodiment of the present application takes the method for generating the central axis of the drainage trough of a landslide lake as an example, in which the method is configured in a device for generating the central axis of the drainage trough of a landslide lake. The device for generating the central axis of the drainage trough of a landslide lake can be applied to any electronic device so that the electronic device can perform the function of generating the central axis of the drainage trough of a landslide lake.

[0033] Among them, the electronic device can be any device with computing capabilities, such as a personal computer (PC), an industrial computer, a host computer, a mobile terminal, a server, etc. The mobile terminal can be, for example, a mobile phone, a tablet computer, a personal digital assistant, a wearable device, and other hardware devices with various operating systems, touch screens and / or display screens.

[0034] like Figure 1 As shown, the method for generating the central axis of the drainage channel of the barrier lake may include the following steps:

[0035] Step 101, based on the topographic data of the damming body corresponding to the damming lake, determine the main water flow path points of the damming body from a plurality of grid cells corresponding to the damming body.

[0036] The terrain data may include but are not limited to terrain elevation data. It should be noted that the terrain data may be data in three-dimensional space.

[0037] The grid unit may be obtained by dividing the area corresponding to the dam body.

[0038] As an example, a digital elevation model (DEM) can be constructed based on the terrain elevation data in the terrain data of the dam body corresponding to the landslide lake; the terrain elevation data of the dam body is read, and the area corresponding to the dam body in the digital elevation model is divided into grid units according to the required resolution (for example, 12.5 meters) to obtain multiple grid units.

[0039] Among them, the main water flow path point is used to indicate the main flow path of the water flow.

[0040] In the embodiment of the present application, the main water flow path points of the damming body can be determined from a plurality of grid cells corresponding to the damming body based on the topographic data of the damming body corresponding to the damming lake.

[0041] Step 102: determining the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body.

[0042] Among them, the monitoring data may include but are not limited to the horizontal displacement, vertical settlement, moisture content, pore water pressure, etc. of the corresponding monitoring point, and this application does not impose any restrictions on this.

[0043] Among them, the horizontal displacement can be used to measure the horizontal movement degree of the landslide body at the monitoring point; the vertical settlement can be used to measure the vertical settlement degree of the landslide body at the monitoring point.

[0044] The stability index value can be used to measure the stability of the dam body of the corresponding grid unit.

[0045] As a possible implementation method, monitoring equipment may be deployed at each monitoring point, so that any monitoring point may be monitored by the corresponding monitoring equipment to obtain monitoring data of the monitoring point.

[0046] It is important to understand that moisture content can be used to indicate the proportion of water content inside the soil, which affects the strength and stability of the soil, while pore water pressure is used to indicate the pressure of water in the pores, which affects the effective stress of the soil.

[0047] In the embodiment of the present application, the stability index value of each grid unit can be determined based on the monitoring data of multiple monitoring points on the landslide body, that is, the stability index value of the landslide body area corresponding to each grid unit can be determined.

[0048] Step 103, based on the main water flow path points of the barrier body and the stability index values ​​of each grid unit, a central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body.

[0049] In order to accurately obtain the geological model corresponding to the dammed body, optionally, a digital elevation model can be constructed based on the terrain elevation data in the terrain data of the dammed body corresponding to the dammed lake, and the digital elevation model can be determined as the geological model corresponding to the dammed body.

[0050] In the embodiment of the present application, the central axis of the drainage channel of the barrier lake can be determined and generated in the geological model corresponding to the barrier body according to the main water flow path points of the barrier body and the stability index value of each grid unit.

[0051] The method for generating the central axis of the drainage channel of the barrier lake in the embodiment of the present application determines the main water flow path point of the barrier body from multiple grid units corresponding to the barrier body based on the terrain data of the barrier body corresponding to the barrier lake; wherein the grid unit is obtained by dividing the area corresponding to the barrier body; according to the monitoring data of multiple monitoring points on the barrier body, the stability index value of each grid unit is determined; according to the main water flow path point of the barrier body and the stability index value of each grid unit, the central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body. In this way, the central axis of the drainage channel of the barrier lake can be automatically and efficiently generated based on the main water flow path point and stability index value of the barrier body, providing strong support for emergency management of barrier lake disasters or water conservancy project design.

[0052] In order to clearly illustrate how, in the above embodiments of the present application, the main water flow path points of the dam body are determined from multiple grid units corresponding to the dam body based on the topographic data of the dam body corresponding to the dam lake, the present application also proposes a method for generating the central axis of the drainage channel of the dam lake.

[0053] Figure 2 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application.

[0054] like Figure 2 As shown, the method for generating the central axis of the drainage channel of the barrier lake may include the following steps:

[0055] Step 201, for any grid cell, determine the elevation information of the grid cell according to the topographic data of the dam body.

[0056] It should be noted that the explanation of the grid unit in step 101 is also applicable to this embodiment and will not be repeated here.

[0057] Optionally, the elevation information of the grid cell may include terrain elevation data of various positions (such as a center position, a grid cell vertex position, a boundary position, etc.) in the damming body area corresponding to the grid cell.

[0058] In the embodiment of the present application, for any grid unit, the elevation information of the grid unit can be determined according to the topographic data of the dam body.

[0059] Step 202: Determine the target elevation corresponding to the grid cell according to the elevation information of the grid cell.

[0060] As an example, when the elevation information of the grid cell includes terrain elevation data of each position in the damming body area corresponding to the grid cell, the terrain elevation data of the center position may be determined as the target elevation corresponding to the grid cell.

[0061] As another example, when the elevation information of a grid cell includes terrain elevation data of each position in the landslide dam area corresponding to the grid cell, the average value of the terrain elevation data of each position corresponding to the grid cell can be determined as the target elevation corresponding to the grid cell.

[0062] Step 203, based on the target elevation corresponding to each grid cell, determine the main water flow path point of the damming body from the multiple grid cells corresponding to the damming body.

[0063] As a possible implementation, Figure 3 As shown, the following steps can be used to determine the main flow path point of the dam body:

[0064] Step 2031, for any grid cell, according to the target elevation corresponding to the grid cell and the target elevations corresponding to a plurality of grid cells adjacent to the grid cell, determine the grid cell to which the water flow of the grid cell is directed.

[0065] Specifically, for any grid cell, the grid cell and a plurality of grid cells adjacent to the grid cell, the grid cell with the smallest target elevation, can be determined as the grid cell to which the water flow of the grid cell is directed.

[0066] It should be noted that the grid cell to which the grid cell water flow is directed may be the grid cell itself, or may be a grid cell among a plurality of grid cells adjacent to the grid cell.

[0067] As an example, assuming that for grid cell 1, the grid cells adjacent to the grid cell include grid cell 2, grid cell 3, and grid cell 4, the grid cell with the smallest target elevation can be determined from grid cell 1, grid cell 2, grid cell 3, and grid cell 4, and the grid cell with the smallest target elevation can be determined as the grid cell to which the water flow of grid cell 1 is directed. For example, when the grid cell with the smallest target elevation is grid cell 1, the grid cell to which the water flow of grid cell 1 is directed is the grid cell itself. For another example, when the grid cell with the smallest target elevation is grid cell 3, the grid cell 3 to which the water flow of grid cell 1 is directed is a grid cell among the multiple grid cells adjacent to grid cell 1.

[0068] Step 2032: for any grid unit, count the number of first grid units whose water flows point to the grid unit among the multiple grid units of the dam body, and determine the accumulated flow value of the grid unit.

[0069] It should be noted that the flow accumulation value of a grid cell is used to indicate the number of upstream grid cells through which water flows. The larger the flow accumulation value of a grid cell is, the larger the catchment area of ​​the grid cell is.

[0070] In an embodiment of the present application, for any grid unit, a grid unit in which the water flow direction is directed toward the grid unit is determined from multiple grid units corresponding to the dam body, and the grid unit in which the water flow direction is directed toward the grid unit is determined as the first grid unit; then, the number of first grid units in which the water flow direction is directed toward the grid unit is counted to obtain the accumulated flow value of the grid unit.

[0071] As an example, assuming that the number of grid cells corresponding to the dam body is 500 grid cells, wherein the 500 grid cells include grid cell a, the number of the first grid cells whose water flows are directed toward grid cell a among the 500 grid cells is counted. For example, the number of the first grid cells whose water flows are directed toward grid cell a is 30, then 30 is determined as the accumulated flow value of grid cell a.

[0072] Step 2033: According to the accumulated flow values ​​of each grid unit, a second grid unit is selected from the multiple grid units of the dam body, and the second grid unit is determined as a main water flow path point.

[0073] It should be noted that the present application does not limit the number of second grid units.

[0074] That is to say, the main water flow path point can be one or more, and this application does not impose any restrictions on this.

[0075] In order to accurately select the second grid cell from multiple grid cells, as a possible implementation method, for any grid cell, if the flow accumulation value of the grid cell is greater than a set threshold, the grid cell can be determined to be the second grid cell. The set threshold can be pre-set, such as 25, 20, etc., and this application does not limit this.

[0076] Therefore, based on the target elevation of each grid cell, the water flow direction of each grid cell is determined, and then, according to the water flow direction of each grid cell, the cumulative flow value of each grid cell is determined. Finally, according to the cumulative flow value of each grid cell, the main water flow path point corresponding to the dam body is effectively determined.

[0077] Step 204: Determine the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body.

[0078] Step 205, based on the main water flow path points of the barrier body and the stability index values ​​of each grid unit, the central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body.

[0079] It should be noted that the execution process of step 204 to step 205 can refer to the execution process of any embodiment of the present application and will not be repeated here.

[0080] The method for generating the central axis of the drainage channel of the barrier lake in the embodiment of the present application determines the elevation information of the grid unit according to the terrain data of the barrier body for any grid unit; determines the target elevation corresponding to the grid unit according to the elevation information of the grid unit; and determines the main water flow path point of the barrier body from the multiple grid units corresponding to the barrier body based on the target elevation corresponding to each grid unit. Thus, based on the target elevation of each grid unit, the main water flow path point of the barrier body is determined from the multiple grid units corresponding to the barrier body, thereby improving the accuracy and effectiveness of the central axis of the drainage channel of the barrier lake.

[0081] Based on the above embodiments of the present application, when the number of main water flow path points is multiple, in order to clearly explain how to determine and generate the central axis of the drainage trough of the impounded lake based on the main water flow path points of the impounded body and the stability index value of each grid unit, the present application also proposes a method for generating the central axis of the drainage trough of the impounded lake.

[0082] Figure 4 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application.

[0083] like Figure 4 As shown, the method for generating the central axis of the drainage channel of the barrier lake may include the following steps:

[0084] Step 401, for any grid cell, determine the elevation information of the grid cell according to the topographic data of the dam body.

[0085] Step 402, determining the target elevation corresponding to the grid cell according to the elevation information of the grid cell.

[0086] Step 403, based on the target elevation corresponding to each grid cell, determine the main water flow path point of the damming body from the multiple grid cells corresponding to the damming body.

[0087] Step 404: Determine the stability index value of each grid unit based on the monitoring data of multiple monitoring points on the dam body.

[0088] It should be noted that the execution process of steps 401 to 404 can refer to the execution process of any embodiment of the present application and will not be repeated here.

[0089] Step 405: for any grid cell, determine the target weight corresponding to the grid cell according to the flow accumulation value and the stability index value of the grid cell.

[0090] As a possible implementation method, for any grid cell, the stability index value of the grid cell is reverse normalized to obtain the first sub-weight corresponding to the grid cell; the inverse of the cumulative flow value of the grid cell is calculated to obtain the second sub-weight corresponding to the grid cell; the first sub-weight and the second sub-weight corresponding to the grid cell are weighted summed to obtain the target weight corresponding to the grid cell.

[0091] As an example, for the ith grid cell among the plurality of grid cells, the stability index value of the ith grid cell is S i , the cumulative flow value is A i , the stability index value S of the i-th grid unit can be calculated according to the following formula: i Perform reverse normalization to obtain the first sub-weight ω corresponding to the i-th grid cell stability :

[0092]

[0093] Among them, S max and S min Respectively represent the maximum and minimum values ​​of the stability index;

[0094] The flow accumulation value A of the i-th grid unit can be calculated according to the following formula: i Find the reciprocal and get the second sub-weight ω corresponding to the i-th grid cell terrain :

[0095]

[0096] The target weight ω can be obtained by weighted summing the first sub-weight and the second sub-weight corresponding to the i-th grid cell according to the following formula: m :

[0097] ω m =ω1*ω stability +ω2*ω terrain : (3)

[0098] Among them, the present application does not limit the values ​​of ω1 and ω2. For example, the value of ω1 is 0.6 and the value of ω2 is 0.4, or the value of ω1 is 0.5 and the value of ω2 is 0.5.

[0099] Therefore, relying on the flow accumulation value and stability index value of the grid cell, the target weight corresponding to the grid cell can be effectively obtained. That is to say, when determining the target weight of the grid cell, the catchment area and soil stability of the grid cell are combined, and the actual situation of the dam area corresponding to the grid cell is comprehensively considered, so that the determination of the target weight of the grid cell is more accurate.

[0100] Step 406, determining the starting point of the central axis of the drainage channel of the barrier lake from multiple main water flow path points, and determining the end point of the central axis of the drainage channel.

[0101] As an example, the main water flow path point in the area where the landslide lake is located upstream of the landslide body among multiple main water flow path points can be determined as the starting point of the central axis of the drainage channel of the landslide lake, and any area in the low-lying area downstream of the landslide body can be determined as the end point of the central axis of the drainage channel.

[0102] Step 407: Based on the target weight corresponding to each grid unit, a central axis of the drainage channel from the starting point to the end point is generated in the geological model corresponding to the dam body.

[0103] As a possible implementation method, a path search algorithm can be used to generate the central axis of the drainage channel from the starting point to the end point in the geological model corresponding to the landslide body based on the target weight corresponding to each grid unit.

[0104] Specifically, a path search algorithm (such as Dijkstra's algorithm (also known as Dijkstra's algorithm, English name is Dijkstra's Algorithm), A-star algorithm (also known as A* search algorithm, English name is A-star Algorithm)) can be used to start from the starting point in the geological model corresponding to the landslide body, based on the target weight corresponding to each grid unit, automatically search and intelligently generate the optimal path from the starting point to the end point, and use the optimal path as the central axis of the drainage trough.

[0105] The method for generating the central axis of the drainage channel of the barrier lake in the embodiment of the present application determines the target weight corresponding to any grid unit according to the flow accumulation value and stability index value of the grid unit; determines the starting point and end point of the central axis of the drainage channel of the barrier lake from multiple main water flow path points; and generates the central axis of the drainage channel from the starting point to the end point in the geological model corresponding to the barrier body based on the target weight corresponding to each grid unit. Thus, depending on the flow accumulation value and stability index value of the grid unit, the target weight corresponding to the grid unit can be effectively obtained, and then, the automatic generation of the central axis of the drainage channel can be realized based on the target weight corresponding to each grid unit intelligently.

[0106] Based on any of the above embodiments of the present application, when the monitoring data includes the horizontal displacement and vertical settlement of the corresponding monitoring point, in order to clearly explain how to determine the stability index value of each grid unit based on the monitoring data of multiple monitoring points on the landslide body, the present application also proposes a method for generating the central axis of the drainage channel of the landslide lake.

[0107] Figure 5 A schematic flow chart of a method for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application.

[0108] like Figure 5 As shown, based on any of the above embodiments of the present application, the method for generating the central axis of the drainage channel of the barrier lake may further include the following steps:

[0109] Step 501: for any monitoring point, determine the monitoring area of ​​the monitoring point on the landslide body according to the horizontal displacement and vertical settlement corresponding to the monitoring point.

[0110] As an example, for any monitoring point, assume that the horizontal displacement corresponding to the monitoring point is D h , the vertical settlement corresponding to the monitoring point is D v , then the radius R of the monitoring area of ​​the monitoring point on the dam body can be determined according to the following formula:

[0111]

[0112] Wherein, θ is the slope angle corresponding to the monitoring point in the dam body, k is the adjustment factor, and this application does not limit the value of k;

[0113] Finally, a circular area is determined on the landslide body with the location of the monitoring point as the center and R as the radius, and the circular area is determined as the monitoring area of ​​the monitoring point on the landslide body.

[0114] Step 502: Determine the anti-slip force of the monitoring point based on the monitoring data corresponding to the monitoring area.

[0115] Among them, the anti-sliding force can be used to indicate the shear strength of the soil at the corresponding monitoring point (also called shear strength).

[0116] As a possible implementation method, when the monitoring data include water content and pore water pressure, such as Figure 6 As shown, the following steps can be used to determine the anti-slip force of the monitoring point:

[0117] Step 5021: For any grid cell in the monitoring area corresponding to the monitoring point, obtain the slope angle corresponding to the grid cell.

[0118] As an example, for any grid cell in the monitoring area corresponding to the monitoring point, the slope value corresponding to the grid cell can be determined based on the terrain elevation data in the terrain data of the grid cell; and the slope angle corresponding to the grid cell can be determined based on the slope value corresponding to the grid cell. For example, for any grid cell in the monitoring area corresponding to the monitoring point, the elevation change value Δh and the horizontal change distance Δx corresponding to the grid cell can be determined based on the terrain elevation data in the terrain data of the grid cell. Then the slope value corresponding to the grid cell is Δh / Δx. According to the following formula, the slope angle θ corresponding to the grid cell is determined as:

[0119]

[0120] Step 5022: Determine the grid gravity corresponding to the grid cell based on the pore water pressure corresponding to the grid cell.

[0121] The grid gravity can be used to indicate the gravity of the soil in the corresponding grid unit under water-containing conditions.

[0122] As a possible implementation method, the original soil gravity corresponding to the grid cell is determined based on the sliding surface area of ​​the grid cell; wherein the original soil gravity can be used to indicate the gravity of the soil in the corresponding grid cell without water; according to the original soil gravity and pore water pressure corresponding to the grid cell, the grid gravity corresponding to the grid cell is determined.

[0123] As an example, assume that for the i-th grid cell, the sliding surface area of ​​the grid cell is A i , the original soil gravity W of the i-th grid unit can be determined according to the following formula: i :

[0124] W i =A i h i ρg; (6)

[0125] Among them, h i It represents the difference between the elevation corresponding to the grid unit and the elevation of the barrier lake water level, ρ represents the density of the soil, and g represents the acceleration of gravity;

[0126] Finally, based on the original soil gravity W corresponding to the i-th grid cell i and pore water pressure u i , the grid gravity W corresponding to the i-th grid unit can be determined according to the following formula i 'for:

[0127] W i ′=W i -u i A i ; (7)

[0128] It should be noted that the pore water pressure corresponding to the i-th grid unit is the pore water pressure of the monitoring point corresponding to the monitoring area to which the i-th grid unit belongs.

[0129] Step 5023, based on the grid gravity, slope angle and moisture content corresponding to each grid unit, determine the anti-slip force of the monitoring point.

[0130] As a possible implementation method, for any grid cell in the monitoring area corresponding to the monitoring point, first, the cohesion corresponding to the grid cell is determined based on the moisture content corresponding to the grid cell and the sliding surface area of ​​the grid cell; secondly, the friction corresponding to the grid cell is determined based on the grid gravity, slope angle and moisture content corresponding to the grid cell; finally, the anti-slip force of the monitoring point is determined based on the cohesion corresponding to each grid cell and the friction corresponding to each grid cell in the monitoring area corresponding to the monitoring point.

[0131] Still using the above example to illustrate, for the i-th grid cell in the monitoring area corresponding to the monitoring point, first, the following formula can be used based on the moisture content m corresponding to the i-th grid cell: i , and the sliding surface area A of the i-th grid unit i , determine the cohesion F corresponding to the i-th grid unit r,c,i :

[0132] F r,c,i =f1(m i )A i ; (8)

[0133] f1(m i )=c·(1-αm i ); (9)

[0134] Where c represents the initial cohesion of the soil in a dry state, and its value range is generally 5-50 kPa; α is an empirical constant, which represents the influence of water content on cohesion; f1(mi) represents the initial cohesion of the soil under the condition of water content m. i The cohesion of soil under the influence of soil;

[0135] Secondly, the grid gravity W corresponding to the i-th grid unit can be calculated according to the following formula: i ′, slope angle θ i , moisture content m i , determine the friction force F corresponding to the i-th grid unit r ' ,μ,i :

[0136] F r ' ,μ,i =W i ′cosθ i ·f2(m i ); (10)

[0137] f2(m i )=μ·(1-βm i ); (11)

[0138] Among them, μ represents the initial friction coefficient of soil in dry state, and its value range is generally 0.5-0.8; β is an empirical constant, which represents the influence of water content on friction coefficient; f2(mi) represents the influence of water content m on friction coefficient. i Friction coefficient of soil under the influence of soil;

[0139] Finally, the anti-slip force F of the monitoring point can be determined according to the following formula based on the cohesion corresponding to each grid unit in the monitoring area corresponding to the monitoring point and the friction corresponding to each grid unit: r :

[0140]

[0141] Where n is the number of grid cells in the monitoring area corresponding to the monitoring point.

[0142] Therefore, the anti-slip force of the monitoring point can be effectively and accurately determined by combining the cohesion and friction force corresponding to the grid unit in the monitoring area corresponding to the monitoring point.

[0143] Step 503: Obtain a first coefficient corresponding to the monitoring point according to the anti-slip force of the monitoring point.

[0144] Among them, the first coefficient can be used to measure the stability of the corresponding monitoring point.

[0145] As a possible implementation method, the sliding force of the monitoring point is determined based on the grid gravity and slope angle corresponding to each grid unit in the monitoring area corresponding to the monitoring point; and the first coefficient corresponding to the monitoring point is determined according to the sliding force and anti-sliding force of the monitoring point.

[0146] As an example, based on the grid gravity and slope angle corresponding to each grid unit in the monitoring area corresponding to the monitoring point, the sliding force of the monitoring point can be determined according to the following formula:

[0147]

[0148] Among them, F d represents the sliding force at the monitoring point, W i ′ represents the grid gravity corresponding to the i-th grid unit, θ i The slope angle corresponding to the i-th grid cell, n represents the number of grid cells in the monitoring area corresponding to the monitoring point;

[0149] Finally, the sliding force F at the monitoring point can be r and anti-slip force F d , according to the following formula, determine the first coefficient F corresponding to the monitoring point s :

[0150] F s =F r / Fd ; (14)

[0151] Among them, when F s ≥1, it means that the area corresponding to the monitoring point is stable. s When <1, it indicates that the area corresponding to the monitoring point may be unstable.

[0152] Step 504: Determine the stability index value corresponding to each grid unit corresponding to the damming body according to the first coefficient corresponding to each monitoring point.

[0153] As a possible implementation method, for any monitoring point, the first coefficient corresponding to the monitoring point can be determined as the stability index value of each grid cell in the monitoring area of ​​the monitoring point; the inverse distance weighted difference algorithm can be used to determine the stability index value of the grid cell whose stability index value has not been determined among the multiple grid cells corresponding to the landslide body based on the stability index value of each grid cell corresponding to each monitoring point.

[0154] For example, Figure 7 The grid cells shown are grid cells 1 to 16. Assuming that the corresponding stability index values ​​of grid cells 4, 6, 10, and 14 are not determined, the inverse distance weighted difference algorithm is used to perform difference calculation, and the stability index value S4 corresponding to the grid cell 4 is determined according to the following formula:

[0155]

[0156] Wherein, n represents the number of horizontal (or vertical) grid cells in the interpolation area considered when interpolating and filling the stability index value of the grid cells. For example, in this example, n can be 4; S i Represents the stability index value of the i-th grid cell in the horizontal (or vertical) grid cell in the interpolation area; d i Represents the distance between grid cell 4 and the i-th grid cell.

[0157] Thus, the stability index value of each grid unit corresponding to the dam body can be determined.

[0158] The method for generating the central axis of the drainage channel of the barrier lake in the embodiment of the present application determines the monitoring area of ​​the monitoring point on the barrier body for any monitoring point; determines the anti-sliding force of the monitoring point based on the monitoring data corresponding to the monitoring area; wherein the anti-sliding force is used to indicate the shear strength of the soil at the corresponding monitoring point; obtains the first coefficient corresponding to the monitoring point according to the anti-sliding force of the monitoring point; wherein the first coefficient is used to measure the stability of the corresponding monitoring point; determines the stability index value corresponding to each grid unit corresponding to the barrier body according to the first coefficient corresponding to each monitoring point. Thus, by analyzing the anti-sliding force of the area corresponding to each monitoring point, the stability index value corresponding to each grid unit corresponding to the barrier body is effectively determined, thereby improving the accuracy of the stability index value.

[0159] With the above Figures 1 to 6 Corresponding to the method for generating the central axis of the drainage channel of the barrier lake provided in the embodiment, the present application also provides a device for generating the central axis of the drainage channel of the barrier lake. Figures 1 to 6 The embodiment provides a method for generating the central axis of the drainage trough of a barrier lake, so the implementation method of the method for generating the central axis of the drainage trough of a barrier lake is also applicable to the device for generating the central axis of the drainage trough of a barrier lake provided in the embodiment of the present application, and will not be described in detail in the embodiment of the present application.

[0160] Figure 8 A schematic structural diagram of a device for generating a central axis of a drainage channel of a barrier lake provided in another embodiment of the present application.

[0161] like Figure 8 As shown, the device 800 for generating the central axis of the drainage channel of the barrier lake may include: a first determination module 801 , a second determination module 802 and a generation module 803 .

[0162] The first determination module 801 is used to determine the main water flow path points of the dam body from multiple grid cells corresponding to the dam body based on the topographic data of the dam body corresponding to the dammed lake; wherein the grid cells are obtained by dividing the area corresponding to the dam body.

[0163] The second determination module 802 is used to determine the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body.

[0164] The generation module 803 is used to generate the central axis of the drainage channel of the barrier lake in the geological model corresponding to the barrier body according to the main water flow path points of the barrier body and the stability index value of each grid unit.

[0165] In a possible implementation of an embodiment of the present application, the first determination module 801 is used to: for any grid cell, determine the elevation information of the grid cell according to the terrain data of the dam body; determine the target elevation corresponding to the grid cell according to the elevation information of the grid cell; and determine the main water flow path point of the dam body from multiple grid cells corresponding to the dam body based on the target elevations corresponding to each grid cell.

[0166] In a possible implementation of the embodiment of the present application, the first determination module 801 is used to: for any grid cell, determine the grid cell to which the water flow of the grid cell is directed based on the target elevation corresponding to the grid cell and the target elevations corresponding to multiple grid cells adjacent to the grid cell; for any grid cell, count the number of first grid cells in which the water flow is directed to the grid cell among the multiple grid cells of the dam body, and determine the cumulative flow value of the grid cell; based on the cumulative flow value of each grid cell, screen out a second grid cell from the multiple grid cells of the dam body, and determine the second grid cell as the main water flow path point.

[0167] In a possible implementation of an embodiment of the present application, there are multiple main water flow path points; a generation module 803 is used to: for any grid unit, determine the target weight corresponding to the grid unit according to the accumulated flow value and stability index value of the grid unit; determine the starting point of the central axis of the drainage channel of the barrier lake from the multiple main water flow path points, and determine the end point of the central axis of the drainage channel; based on the target weights corresponding to each grid unit, generate the central axis of the drainage channel from the starting point to the end point in the geological model corresponding to the barrier body.

[0168] In a possible implementation of an embodiment of the present application, the monitoring data includes the horizontal displacement and vertical settlement of the corresponding monitoring point; the second determination module 802 is used to: for any monitoring point, determine the monitoring area of ​​the monitoring point on the landslide body according to the horizontal displacement and vertical settlement corresponding to the monitoring point; determine the anti-sliding force of the monitoring point based on the monitoring data corresponding to the monitoring area; wherein the anti-sliding force is used to indicate the shear strength of the soil at the corresponding monitoring point; according to the anti-sliding force of the monitoring point, obtain the first coefficient corresponding to the monitoring point; wherein the first coefficient is used to measure the stability of the corresponding monitoring point; according to the first coefficient corresponding to each monitoring point, determine the stability index value corresponding to each grid unit corresponding to the landslide body.

[0169] In a possible implementation of an embodiment of the present application, the monitoring data includes moisture content and pore water pressure; the second determination module 802 is used to: obtain the slope angle corresponding to the grid cell for any grid cell in the monitoring area corresponding to the monitoring point; determine the grid gravity corresponding to the grid cell based on the pore water pressure corresponding to the grid cell; wherein the grid gravity is used to indicate the gravity of the soil in the corresponding grid cell under water-containing conditions; and determine the anti-sliding force of the monitoring point based on the grid gravity, slope angle and moisture content corresponding to each grid cell.

[0170] In a possible implementation of an embodiment of the present application, the second determination module 802 is used to: determine the sliding force of the monitoring point based on the grid gravity and slope angle corresponding to each grid unit in the monitoring area corresponding to the monitoring point; determine the first coefficient corresponding to the monitoring point according to the sliding force and anti-sliding force of the monitoring point.

[0171] In a possible implementation of an embodiment of the present application, the second determination module 802 is used to: determine the original soil gravity corresponding to the grid cell based on the sliding surface area of ​​the grid cell; wherein the original soil gravity is used to indicate the gravity of the soil in the corresponding grid cell in the absence of water; and determine the grid gravity corresponding to the grid cell based on the original soil gravity and pore water pressure corresponding to the grid cell.

[0172] In a possible implementation of an embodiment of the present application, the second determination module 802 is used to: for any grid cell in the monitoring area corresponding to the monitoring point, determine the cohesion corresponding to the grid cell based on the moisture content corresponding to the grid cell and the sliding surface area of ​​the grid cell; determine the friction corresponding to the grid cell based on the grid gravity, slope angle and moisture content corresponding to the grid cell; determine the anti-slip force of the monitoring point according to the cohesion corresponding to each grid cell and the friction corresponding to each grid cell in the monitoring area corresponding to the monitoring point.

[0173] The device for generating the central axis of the drainage channel of the barrier lake in the embodiment of the present application determines the main water flow path point of the barrier body from multiple grid units corresponding to the barrier body based on the terrain data of the barrier body corresponding to the barrier lake; wherein the grid unit is obtained by dividing the area corresponding to the barrier body; according to the monitoring data of multiple monitoring points on the barrier body, the stability index value of each grid unit is determined; according to the main water flow path point of the barrier body and the stability index value of each grid unit, the central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body. Thus, it is possible to automatically and efficiently generate the central axis of the drainage channel of the barrier lake based on the main water flow path point and stability index value of the barrier body, providing strong support for emergency management of barrier lake disasters or water conservancy project design.

[0174] In order to implement the above-mentioned embodiments, the present application also proposes an electronic device, characterized in that it includes a memory, a processor, and a computer program stored in the memory and run on the processor. When the processor executes the program, it implements the method for generating the central axis of the drainage channel of the barrier lake proposed in any of the aforementioned embodiments of the present application.

[0175] In order to implement the above embodiments, the present application also proposes a non-temporary computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the method for generating the central axis of the drainage channel of a barrier lake as proposed in any of the aforementioned embodiments of the present application is implemented.

[0176] In order to implement the above embodiments, the present application also provides a computer program product. When the instructions in the computer program product are executed by a processor, the method for generating the central axis of the drainage channel of the barrier lake proposed in any of the aforementioned embodiments of the present application is executed.

[0177] like Fig. 9 As shown, the electronic device 12 is in the form of a general purpose computing device. The components of the electronic device 12 may include, but are not limited to: one or more processors or processing units 16, a system memory 28, and a bus 18 that connects various system components (including the system memory 28 and the processing unit 16).

[0178] The bus 18 represents one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, a graphics acceleration port, a processor or a local bus using any of a variety of bus structures. For example, these architectures include but are not limited to Industry Standard Architecture (ISA) bus, Micro Channel Architecture (MAC) bus, Enhanced ISA bus, Video Electronics Standards Association (VESA) local bus and Peripheral Component Interconnection (PCI) bus.

[0179] The electronic device 12 typically includes a variety of computer system readable media. These media can be any available media that can be accessed by the electronic device 12, including volatile and non-volatile media, removable and non-removable media.

[0180] The memory 28 may include computer system readable media in the form of volatile memory, such as random access memory (RAM) 30 and / or cache memory 32. The electronic device 12 may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, the storage system 34 may be used to read and write non-removable, non-volatile magnetic media ( Fig. 9 not shown, usually called a "hard drive"). Although Fig. 9 Not shown in the figure, a disk drive for reading and writing a removable non-volatile disk (e.g., a "floppy disk"), and an optical disk drive for reading and writing a removable non-volatile optical disk (e.g., a compact disc read only memory (CD-ROM), a digital versatile disc read only memory (DVD-ROM), or other optical media) may be provided. In these cases, each drive may be connected to the bus 18 via one or more data medium interfaces. The memory 28 may include at least one program product having a set (e.g., at least one) of program modules that are configured to perform the functions of the various embodiments of the present application.

[0181] A program / utility 40 having a set (at least one) of program modules 42 may be stored, for example, in the memory 28, such program modules 42 including, but not limited to, an operating system, one or more application programs, other program modules, and program data, each of which or some combination may include an implementation of a network environment. The program modules 42 generally perform the functions and / or methods of the embodiments described herein.

[0182] The electronic device 12 may also communicate with one or more external devices 14 (e.g., keyboards, pointing devices, displays 24, etc.), may also communicate with one or more devices that enable a user to interact with the electronic device 12, and / or may communicate with any device that enables the electronic device 12 to communicate with one or more other computing devices (e.g., network cards, modems, etc.). Such communication may be performed via an input / output (I / O) interface 22. Furthermore, the electronic device 12 may also communicate with one or more networks (e.g., a local area network (LAN), a wide area network (WAN), and / or a public network, such as the Internet) via a network adapter 20. As shown, the network adapter 20 communicates with other modules of the electronic device 12 via a bus 18. It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in conjunction with the electronic device 12, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems.

[0183] The processing unit 16 executes various functional applications and data processing by running programs stored in the system memory 28, such as implementing the methods mentioned in the above embodiments.

[0184] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0185] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Therefore, the features defined as "first" and "second" may explicitly or implicitly include at least one of the features. In the description of this application, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0186] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0187] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute the instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or more wires (electronic device), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically, for example, by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways if necessary, and then stored in a computer memory.

[0188] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, multiple steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0189] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0190] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0191] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for generating the central axis of a drainage channel of a barrier lake, characterized in that: The method comprises: Based on the topographic data of the barrier body corresponding to the barrier lake, determining the main water flow path point of the barrier body from a plurality of grid units corresponding to the barrier body; wherein the grid units are obtained by dividing the area corresponding to the barrier body; Determining the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body; According to the main water flow path points of the barrier body and the stability index values ​​of each grid unit, the central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body.

2. The method according to claim 1, characterized in that: The determining of the main water flow path points of the damming body from a plurality of grid cells corresponding to the damming body based on the topographic data of the damming body corresponding to the damming lake comprises: For any of the grid cells, determining the elevation information of the grid cell according to the topographic data of the dam body; Determining a target elevation corresponding to the grid cell according to the elevation information of the grid cell; Based on the target elevation corresponding to each of the grid cells, main water flow path points of the damming body are determined from a plurality of grid cells corresponding to the damming body.

3. The method according to claim 2, characterized in that The determining of the main water flow path point of the damming body from a plurality of grid cells corresponding to the damming body based on the target elevation corresponding to each of the grid cells comprises: For any of the grid cells, according to the target elevation corresponding to the grid cell and the target elevations corresponding to a plurality of the grid cells adjacent to the grid cell, determine the grid cell to which the water flow of the grid cell is directed; For any of the grid units, counting the number of first grid units whose water flows point to the grid unit among the multiple grid units of the dam body, and determining a cumulative flow value of the grid unit; According to the accumulated flow values ​​of each of the grid cells, a second grid cell is selected from the plurality of grid cells of the dam body, and the second grid cell is determined as the main water flow path point.

4. The method according to claim 3, characterized in that The number of the main water flow path points is multiple; the central axis of the drainage channel of the barrier lake is generated in the geological model corresponding to the barrier body according to the main water flow path points of the barrier body and the stability index value of each grid unit, including: For any of the grid cells, determine the target weight corresponding to the grid cell according to the flow accumulation value and the stability index value of the grid cell; Determine the starting point of the central axis of the drainage channel of the barrier lake from the plurality of main water flow path points, and determine the end point of the central axis of the drainage channel; Based on the target weight corresponding to each of the grid cells, a central axis of the drainage channel from the starting point to the end point is generated in the geological model corresponding to the dam body.

5. The method according to any one of claims 1 to 4, characterized in that: The monitoring data includes the horizontal displacement and vertical settlement of the corresponding monitoring point; and the stability index value of each grid unit is determined according to the monitoring data of multiple monitoring points on the dam body, including: For any of the monitoring points, determining a monitoring area of ​​the monitoring point on the landslide dam according to the horizontal displacement and vertical settlement corresponding to the monitoring point; Determine the anti-sliding force of the monitoring point based on the monitoring data corresponding to the monitoring area; wherein the anti-sliding force is used to indicate the shear strength of the soil at the corresponding monitoring point; According to the anti-slip force of the monitoring point, a first coefficient corresponding to the monitoring point is obtained; wherein the first coefficient is used to measure the stability of the corresponding monitoring point; According to the first coefficient corresponding to each monitoring point, the stability index value corresponding to each grid unit corresponding to the dam body is determined.

6. The method according to claim 5, characterized in that The monitoring data include water content and pore water pressure; The step of determining the anti-slip force of the monitoring point based on the monitoring data corresponding to the monitoring area includes: For any of the grid cells in the monitoring area corresponding to the monitoring point, obtaining the slope angle corresponding to the grid cell; Based on the pore water pressure corresponding to the grid cell, determining the grid gravity corresponding to the grid cell; wherein the grid gravity is used to indicate the gravity of the soil in the corresponding grid cell under the condition of water; The anti-slip force of the monitoring point is determined based on the grid gravity, the slope angle and the moisture content corresponding to each grid unit.

7. The method according to claim 6, characterized in that The step of obtaining a first coefficient corresponding to the monitoring point according to the anti-slip force of the monitoring point includes: Determine the sliding force of the monitoring point based on the grid gravity and the slope angle corresponding to each grid unit in the monitoring area corresponding to the monitoring point; A first coefficient corresponding to the monitoring point is determined according to the sliding force and the anti-sliding force of the monitoring point.

8. The method according to claim 6, characterized in that The determining the grid gravity corresponding to the grid unit based on the pore water pressure corresponding to the grid unit comprises: Based on the sliding surface area of ​​the grid unit, determine the original soil gravity corresponding to the grid unit; wherein the original soil gravity is used to indicate the gravity of the soil in the corresponding grid unit without water; The grid gravity corresponding to the grid unit is determined according to the original soil gravity and pore water pressure corresponding to the grid unit.

9. The method according to claim 8, characterized in that The step of determining the anti-slip force of the monitoring point based on the grid gravity, the slope angle and the moisture content corresponding to each grid unit includes: For any of the grid cells in the monitoring area corresponding to the monitoring point, based on the water content corresponding to the grid cell and the sliding surface area of ​​the grid cell, determining the cohesion corresponding to the grid cell; Determining the friction force corresponding to the grid cell based on the grid gravity, slope angle and water content corresponding to the grid cell; The anti-slip force of the monitoring point is determined according to the cohesion corresponding to each of the grid units and the friction corresponding to each of the grid units in the monitoring area corresponding to the monitoring point.

10. A device for generating the central axis of a drainage channel of a barrier lake, characterized in that: The device comprises: A first determination module is used to determine the main water flow path point of the damming body from a plurality of grid cells corresponding to the damming body based on the topographic data of the damming body corresponding to the damming lake; wherein the grid cells are obtained by dividing the area corresponding to the damming body; A second determination module is used to determine the stability index value of each grid unit according to the monitoring data of multiple monitoring points on the dam body; A generation module is used to generate the central axis of the drainage channel of the barrier lake in the geological model corresponding to the barrier body according to the main water flow path points of the barrier body and the stability index value of each grid unit.

Images