A multi-parameter environmental coupling monitoring system for loess slopes

Abstract

The present invention relates to the field of slope monitoring technology, and in particular to a loess slope multi-parameter environmental coupling monitoring system, comprising: a data acquisition unit; a state determination unit for determining a slope state according to a slope regeneration index and a slope angle turbulence coefficient; a region division unit for determining a slope region division method according to the slope state; a region selection unit for performing monitoring analysis on each slope region to determine whether to use the slope region as a monitoring region; a monitoring determination unit for determining a monitoring point setting method of a monitoring region according to setting conditions; a first setting unit for determining an associated state of an abnormal sub-region according to an abnormal variation coefficient and an associated abnormal region reference value, and determining a processing method according to the associated state of the abnormal sub-region; and a second setting unit for determining an optimization method according to an abnormal aggregation region type. The present invention can improve the accuracy of loess slope data monitoring.

Description

Technical Field

[0001] The invention relates to the technical field of slope monitoring, and in particular to a loess slope multi-parameter environmental coupling monitoring system. Background Art

[0002] Due to the loose soil structure and developed vertical joints, as well as the influence of long-term weathering and human activities, loess slopes have poor stability and are often threatened by geological disasters such as landslides and collapses. However, due to the complex terrain of loess slopes, monitoring methods are limited and it is impossible to accurately warn of landslide risks, which seriously threatens the safety of life and property of people in the surrounding areas. Therefore, how to improve the accuracy of loess slope monitoring and timely discover and warn of geological disaster risks is a technical problem that technical personnel in this field need to solve urgently.

[0003] Chinese patent publication number CN113295212A discloses a landslide overall stability safety early warning system that uses the synergistic effect of multiple monitoring points, including: historical information collection module, data collection module, environmental data collection module, central control module, data processing module, data analysis module, index determination module, weight determination module, comprehensive evaluation module and early warning module; in the evaluation and monitoring and early warning of slope stability, the deformation, rainfall and geological influence of the displacement monitoring points in the main sliding area of ​​the slope are comprehensively considered, and multiple qualitative and quantitative indicators that can reflect the overall stability of the slope are flexibly selected, and the weights are determined respectively, so that the overall stability of the slope can be comprehensively divided into stages and monitored and warned. It can be seen that the above technical solution has the following problems: the monitoring point setting method is single, and the monitoring point position is not adaptively adjusted according to the actual state of the slope, resulting in poor data monitoring effect. Summary of the invention

[0004] To this end, the present invention provides a loess slope multi-parameter environmental coupling monitoring system to overcome the problem in the prior art that the monitoring point positions are not adaptively adjusted according to the actual state of the slope, resulting in poor data monitoring effect.

[0005] To achieve the above object, the present invention provides a loess slope multi-parameter environmental coupling monitoring system, comprising:

[0006] A data collection unit, used to collect slope data;

[0007] A state determination unit connected to the data acquisition unit and used to determine the state of the slope according to the slope regeneration index and the slope angle turbulence coefficient;

[0008] A region division unit, which is connected to the data acquisition unit and the state determination unit respectively, and is used to determine a slope region division method according to the slope state, wherein the region division method is to perform associated division according to the crack influence coefficient and the distance reference value or to perform uniform division according to the slope creep coefficient;

[0009] A region selection unit, which is connected to the region division unit, is used to perform monitoring analysis on each slope region in descending order of collapse thresholds to determine whether to use the slope region as a monitoring region;

[0010] A monitoring determination unit connected to the area selection unit is used to determine a monitoring point setting method of the monitoring area according to the setting conditions, wherein the monitoring point setting method is to determine a processing method according to the association state of the abnormal sub-area or to determine an optimization method according to the type of abnormal aggregation area;

[0011] A first setting unit, which is connected to the monitoring and determining unit, is used to determine the association state of the abnormal sub-region according to the abnormal variation coefficient and the reference value of the associated abnormal region, and determine a processing method according to the association state of the abnormal sub-region, and the processing method is to set monitoring points according to the associated abnormal region, or to set monitoring points at intervals between adjacent abnormal sub-regions;

[0012] The second setting unit is connected to the monitoring and judgment unit, and is used to determine the type of abnormal aggregation area according to the convexity and concavity ratio, and determine the optimization method according to the type of abnormal aggregation area. The optimization method is to determine the number of monitoring points in the abnormal aggregation area according to the similarity of abnormal sub-areas in the first type of aggregation area or according to the abnormal dislocation coefficient of the second type of aggregation area.

[0013] Furthermore, the state determination unit determines the slope state according to the slope regeneration index and the slope angle turbulence coefficient, and the slope state includes:

[0014] A first slope state in which a slope regeneration index is greater than or equal to a preset slope regeneration index or a slope angle turbulence coefficient is greater than or equal to a preset slope angle turbulence coefficient;

[0015] A second slope state in which the slope regeneration index is less than a preset slope regeneration index and the slope angle turbulence coefficient is less than a preset slope angle turbulence coefficient.

[0016] Furthermore, the area division unit responds to the slope state to divide the slope into areas, wherein:

[0017] The area division unit responds to the first slope state, and the area division method is to perform associated division according to the crack influence coefficient and the distance reference value. In the associated division, the minimum rectangular area that can contain each crack in a single associated crack combination is recorded as a slope area;

[0018] The area division unit responds to the second slope state. The area division method is uniform division according to the slope creep coefficient. In the uniform division, the slope is divided into several rectangular areas with the same area and shape. Each rectangular area is recorded as a slope area. The number of slope areas is positively correlated with the slope creep coefficient.

[0019] Furthermore, the slope creep coefficient is determined according to the number of associated slopes, wherein:

[0020] If the number of associated slopes is greater than or equal to the preset number of associated slopes, the slope creep coefficient is determined according to the associated slope radiation value and the sub-influence threshold;

[0021] If the number of associated slopes is less than the preset number of associated slopes, the slope creep coefficient is determined according to the sub-influence threshold.

[0022] Furthermore, when the region selection unit performs monitoring and analysis on a single slope region, it determines whether to use the slope region as a monitoring region according to the radiation coefficient and the crack assessment value, wherein:

[0023] The slope area with a radiation coefficient less than the preset radiation coefficient or a crack assessment value greater than or equal to the preset crack assessment value is taken as the monitoring area.

[0024] Furthermore, the monitoring determination unit responds to the setting conditions to determine the setting mode of the monitoring points in the monitoring area, wherein:

[0025] The setting condition for the response of the monitoring and judgment unit is that the crack dislocation coefficient is greater than or equal to the preset crack dislocation coefficient or the abnormal distribution threshold is greater than or equal to the preset abnormal distribution threshold, and the monitoring point setting method is to determine the processing method according to the correlation state of the abnormal sub-area;

[0026] The setting condition for the response of the monitoring and judgment unit is that the crack dislocation coefficient is less than the preset crack dislocation coefficient and the abnormal distribution threshold is less than the preset abnormal distribution threshold. The monitoring point setting method is to determine the optimization method according to the type of abnormal aggregation area.

[0027] Furthermore, the first setting unit responds to the association status of the abnormal sub-area to determine the processing method, wherein:

[0028] The first setting unit responds to the first preset association state, and the processing method is to set monitoring points according to the associated abnormal area. When setting the monitoring points according to the associated abnormal area, the associated abnormal area is determined according to the abnormal distance, and the monitoring point is set at the central position corresponding to each associated abnormal area;

[0029] The first setting unit responds to the second preset association state, and the processing method is to set monitoring points at intervals between adjacent abnormal sub-regions. When the monitoring points are set at intervals between adjacent abnormal sub-regions, a monitoring point is set at the center position of each abnormal sub-region, and a monitoring point is set at a corresponding adjacent position of each abnormal sub-region;

[0030] The adjacent position corresponding to a single abnormal sub-region is the midpoint of a line connecting the center position of the abnormal sub-region and the center position of an adjacent abnormal sub-region corresponding to the abnormal sub-region.

[0031] Furthermore, the first setting unit determines the association state of the abnormal sub-region according to the abnormal variation coefficient and the associated abnormal region reference value, and the association state includes:

[0032] A first preset association state in which the abnormal variation coefficient is less than a preset abnormal variation coefficient or the associated abnormal region reference value is greater than or equal to a preset associated abnormal region reference value;

[0033] A second preset association state in which the abnormal variation coefficient is greater than or equal to a preset abnormal variation coefficient and the associated abnormal region reference value is less than a preset associated abnormal region reference value.

[0034] Furthermore, the second setting unit determines the type of abnormal aggregation area according to the convexity-concavity ratio, and the type of abnormal aggregation area includes:

[0035] A type of abnormal clustering area whose convexity-concavity ratio is less than a preset convexity-concavity ratio, and a type of abnormal clustering area whose convexity-concavity ratio is greater than or equal to the preset convexity-concavity ratio.

[0036] Furthermore, the second setting unit determines the optimization method according to the type of abnormal aggregation area, wherein:

[0037] For a type of abnormal clustering area, the optimization method is to determine the number of monitoring points in the abnormal clustering area according to the similarity of the abnormal sub-areas in the type of clustering area;

[0038] For the second type of abnormal aggregation area, the optimization method is to determine the number of monitoring points in the abnormal aggregation area according to the abnormal dislocation coefficient of the second type of aggregation area;

[0039] There is a negative correlation between the number of monitoring points and the similarity of the abnormal sub-regions in the first type of clustering area, and there is a positive correlation between the number of monitoring points and the abnormal dislocation coefficient of the second type of clustering area.

[0040] Compared with the prior art, the beneficial effect of the present invention lies in that, in the technical scheme of the present invention, the slope state is determined according to the slope regeneration index and the slope angle turbulence coefficient, and the risk degree of geological disasters occurring in the loess slope is effectively reflected by the slope regeneration index and the slope angle turbulence coefficient, and then different slope area division methods are adaptively selected according to the slope state, so that the selection of the slope area division method is more in line with the actual application scenario, which is conducive to improving the pertinence and efficiency of monitoring, reducing unnecessary waste of monitoring resources while ensuring the monitoring effect, and thus improving the accuracy of monitoring.

[0041] Furthermore, when monitoring and analyzing a single slope area in the present invention, the radiation coefficient and crack assessment value are used to effectively reflect the degree of influence of the monitoring area on the monitoring effect of the slope area and the crack situation, and then determine whether to use the slope area as a monitoring area based on the radiation coefficient and the crack assessment value, so that the determination of the monitoring area is more in line with the actual application scenario, while ensuring that high-risk areas are fully monitored and covered, the monitoring cost is reduced, thereby improving the monitoring efficiency of loess slopes.

[0042] Furthermore, the present invention determines the setting method of monitoring points in the monitoring area according to the setting conditions, and effectively reflects the dislocation of cracks in the monitoring area and the distribution of abnormal sub-areas through the setting conditions, and then adaptively selects different monitoring point setting methods according to the setting conditions, so that the selection of monitoring point setting methods is more in line with the actual application scenario, avoiding the problem of poor monitoring accuracy caused by unreasonable monitoring point setting, ensuring the reasonable configuration of monitoring points, reducing the setting of unnecessary monitoring points, and thus improving the monitoring accuracy while reducing the monitoring cost.

[0043] Furthermore, in the present invention, the association status of the abnormal sub-region is determined according to the abnormal variation coefficient and the reference value of the associated abnormal region, and the degree of association of the abnormal sub-region is effectively reflected by the abnormal variation coefficient and the reference value of the associated abnormal region, and then different processing methods are adaptively selected according to the association status of the abnormal sub-region, so that the selection of processing method is more in line with the actual application scenario, avoiding the problem of too few monitoring points to meet the monitoring needs when the abnormal sub-region is far away, thereby improving the accuracy of loess slope monitoring.

[0044] Furthermore, in the present invention, the type of abnormal aggregation area is determined according to the convexity and concavity ratio, and the type of abnormal sub-area in the abnormal aggregation area is effectively reflected through the concavity ratio, and then different optimization methods are adaptively selected according to the type of abnormal aggregation area, so that the selection of optimization method is more in line with the actual application scenario, avoiding the problem that the monitoring points set when the abnormal aggregation area is more dangerous cannot meet the monitoring needs, thereby improving the accuracy of the monitoring effect. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 It is a unit connection diagram of the loess slope multi-parameter environmental coupling monitoring system of the present invention;

[0046] Figure 2 This is a flow chart of the present invention for determining a slope area division method according to the slope state;

[0047] Figure 3 It is a flow chart of the method for setting monitoring points in a monitoring area according to setting conditions of the present invention;

[0048] Figure 4 The flowchart of the present invention is to determine the processing method according to the association status of the abnormal sub-region. DETAILED DESCRIPTION

[0049] In order to make the objects and advantages of the present invention more clearly understood, the present invention is further described below in conjunction with embodiments; it should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0050] The preferred embodiments of the present invention are described below with reference to the accompanying drawings. It should be understood by those skilled in the art that these embodiments are only used to explain the technical principles of the present invention and are not intended to limit the protection scope of the present invention.

[0051] It should be noted that, in the description of the present invention, terms such as "up", "down", "left", "right", "inside" and "outside" indicating directions or positional relationships are based on the directions or positional relationships shown in the drawings. This is merely for the convenience of description and does not indicate or imply that the device or element must have a specific orientation, be constructed and operated in a specific orientation. Therefore, it cannot be understood as a limitation on the present invention.

[0052] In addition, it should be noted that in the description of the present invention, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two components. For those skilled in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0053] See also Figures 1 to 4 As shown, the present invention provides a loess slope multi-parameter environmental coupling monitoring system, comprising:

[0054] A data collection unit, used to collect slope data;

[0055] A state determination unit connected to the data acquisition unit and used to determine the state of the slope according to the slope regeneration index and the slope angle turbulence coefficient;

[0056] A region division unit, which is connected to the data acquisition unit and the state determination unit respectively, and is used to determine a slope region division method according to the slope state, wherein the region division method is to perform associated division according to the crack influence coefficient and the distance reference value or to perform uniform division according to the slope creep coefficient;

[0057] A region selection unit, which is connected to the region division unit, is used to perform monitoring analysis on each slope region in descending order of collapse thresholds to determine whether to use the slope region as a monitoring region;

[0058] A monitoring determination unit connected to the area selection unit is used to determine a monitoring point setting method of the monitoring area according to the setting conditions, wherein the monitoring point setting method is to determine a processing method according to the association state of the abnormal sub-area or to determine an optimization method according to the type of abnormal aggregation area;

[0059] A first setting unit, which is connected to the monitoring and determining unit, is used to determine the association state of the abnormal sub-region according to the abnormal variation coefficient and the reference value of the associated abnormal region, and determine a processing method according to the association state of the abnormal sub-region, and the processing method is to set monitoring points according to the associated abnormal region, or to set monitoring points at intervals between adjacent abnormal sub-regions;

[0060] The second setting unit is connected to the monitoring and judgment unit, and is used to determine the type of abnormal aggregation area according to the convexity and concavity ratio, and determine the optimization method according to the type of abnormal aggregation area. The optimization method is to determine the number of monitoring points in the abnormal aggregation area according to the similarity of abnormal sub-areas in the first type of aggregation area or according to the abnormal dislocation coefficient of the second type of aggregation area.

[0061] The application scenario of the present invention is loess slope monitoring. The slope data includes but is not limited to crack location, crack morphology and slope height. The target loess slope in the present invention is the loess slope that needs to be monitored. In the present invention, several historical records are correspondingly set. Any historical record records the coupling threshold, slope regeneration index, slope angle turbulence coefficient, radiation coefficient, crack assessment value, abnormal spacing, abnormal variation coefficient and associated abnormal area reference value in at least one historical process of target loess slope monitoring, and each historical record corresponds to a qualified mark, which records whether the accuracy of the target loess slope monitoring meets user requirements. The qualified mark can be recorded manually.

[0062] After setting the monitoring points, the present invention monitors the slope parameters at each monitoring point. The slope parameters include displacement, displacement rate, stress, rainfall and groundwater level, and performs geological disaster warning when the coupling threshold is greater than the preset coupling threshold. The calculation formula of the coupling threshold δ is: ,in, is the monitoring mean corresponding to the i-th slope parameter. The monitoring mean is the average value of the values ​​corresponding to a single monitoring parameter at each monitoring point. is the weight coefficient corresponding to the i-th slope parameter, i=1,2,…,n, n is the number of slope parameters. It should be noted that the number, category and weight coefficient of slope parameters can be determined by the user according to actual needs. Here, a setting method for the number, category and weight coefficient of slope parameters is provided. The number of slope parameters is 3, and the categories of slope parameters are displacement, stress and rainfall. The weight coefficients corresponding to displacement, stress and rainfall are 0.4, 0.4 and 0.2 respectively. Then the coupling threshold , , as well as They are the monitored mean values ​​corresponding to displacement, stress and rainfall respectively; the value of the preset coupling threshold can be determined by the user according to the actual application scenario. The greater the user's demand for reducing the impact of slope disasters, the larger the value of the preset coupling threshold. A value of a preset coupling threshold is provided, and the minimum value of the coupling threshold corresponding to the historical record that can meet the user's needs is recorded as the preset coupling threshold.

[0063] Specifically, the state determination unit determines the slope state according to the slope regeneration index and the slope angle turbulence coefficient, and the slope state includes:

[0064] A first slope state in which a slope regeneration index is greater than or equal to a preset slope regeneration index or a slope angle turbulence coefficient is greater than or equal to a preset slope angle turbulence coefficient;

[0065] A second slope state in which the slope regeneration index is less than a preset slope regeneration index and the slope angle turbulence coefficient is less than a preset slope angle turbulence coefficient.

[0066] Wherein, slope regeneration index = number of sliding surfaces + overlapping area of ​​sliding surfaces, the number of sliding surfaces is the total amount of sliding surfaces in the target loess slope, the overlapping area of ​​sliding surfaces is the average value of the reference areas corresponding to each sliding surface, for a single sliding surface, the sliding surface is recorded as the target sliding surface, and the other sliding surfaces except the target sliding surface are recorded as reference sliding surfaces, the reference area corresponding to the target sliding surface is the maximum value of the overlapping areas of the target sliding surface and each reference sliding surface, the sliding surface is the interface between two adjacent soil layers in the target loess slope, and the soil layers include but are not limited to silt, silty clay and sandy soil, and the user can obtain the number, location and area of ​​sliding surfaces through geological survey, which is easy for technicians in this field to understand and will not be described in detail;

[0067] The slope angle turbulence coefficient is the standard deviation of the slope angles corresponding to each slope surface. The slope surface refers to the area on the surface of the target loess slope where the angle with the horizontal plane does not change. The slope angle corresponding to a single slope surface is the angle between the slope surface and the horizontal plane.

[0068] The values ​​of the preset slope regeneration index and the preset slope angle turbulence coefficient can be determined by the user according to the actual application scenario. The greater the user's demand for associated division based on the crack influence coefficient and the distance reference value, the more values ​​of a preset slope regeneration index and a preset slope angle turbulence coefficient are provided, and the historical records that are evenly divided according to the slope creep coefficient are detected. The average value of the slope regeneration index corresponding to the historical records that can meet the user's needs is recorded as the preset slope regeneration index, and the average value of the slope angle turbulence coefficient corresponding to the historical records that can meet the user's needs is recorded as the preset slope angle turbulence coefficient.

[0069] Specifically, the area division unit responds to the slope state in a slope area division manner, wherein:

[0070] The area division unit responds to the first slope state, and the area division method is to perform associated division according to the crack influence coefficient and the distance reference value. In the associated division, the minimum rectangular area that can contain each crack in a single associated crack combination is recorded as a slope area;

[0071] The area division unit responds to the second slope state. The area division method is uniform division according to the slope creep coefficient. In the uniform division, the slope is divided into several rectangular areas with the same area and shape. Each rectangular area is recorded as a slope area. The number of slope areas is positively correlated with the slope creep coefficient.

[0072] The associated crack combination is determined according to the crack influence coefficient and the distance reference value, wherein the associated analysis is performed on each crack, and when the combined analysis is performed on a single crack, the crack is recorded as a target crack, and other cracks other than the target crack that are not recorded in the associated crack combination are recorded as reference cracks, and a combination of reference cracks whose crack influence coefficient with the target crack is less than a preset crack influence coefficient and whose distance reference value is less than a preset distance reference value and the target crack is recorded as an associated crack combination, and the associated analysis is continued for the cracks that are not recorded in the associated crack combination until all cracks are recorded in the associated crack combination;

[0073] The method for confirming the crack influence coefficient is that, for any two cracks, the larger value of the sub-influence coefficients corresponding to the two cracks is recorded as A1, and the smaller value is recorded as A2. The crack influence coefficient = (A1-A2) / A1,

[0074] The sub-influence coefficient is confirmed as follows:

[0075] If the impact angle is less than the preset impact angle, the sub-influence coefficient = crack threshold - impact angle;

[0076] If the influence angle is greater than or equal to the preset influence angle, the sub-influence coefficient is positively correlated with the crack threshold;

[0077] The method for confirming the influence angle and crack threshold is that, for a single crack, the crack is recorded as the target crack. If the target crack is a straight line, the influence angle is the minimum angle between the target crack and the baseline. If the target crack is a curved line, the influence angle is the minimum angle between the line connecting the two end points of the target crack and the baseline. The baseline corresponding to the target crack is the straight line at the intersection of the slope corresponding to the target crack and the adjacent slope. The crack threshold corresponding to the target crack = crack length + maximum crack width + maximum crack depth, the maximum crack width is the maximum width of the target crack, and the maximum crack depth is the average of the shortest distance from the height of the surface on both sides of the target crack to the lowest point of the target crack. The distance reference value is the shortest distance between the center positions corresponding to the two cracks, and the center position corresponding to a single crack is the center of the circumscribed circle of the crack.

[0078] The values ​​of the preset influence angle, the preset crack influence coefficient and the preset distance reference value can be determined by the user according to the actual application scenario. The larger the value of the preset influence angle, the greater the user's demand for determining the sub-influence coefficient based on the crack threshold and the influence angle. A value of the preset influence angle is provided, and the preset influence angle is 70°. The greater the user's demand for the degree of crack correlation within the associated crack combination, the smaller the values ​​of the preset crack influence coefficient and the preset distance reference value. A value of the preset crack influence coefficient and the preset distance reference value is provided, and the preset crack influence coefficient is 30%, and the preset distance reference value is 100m.

[0079] Specifically, the slope creep coefficient is determined according to the number of associated slopes, where:

[0080] If the number of associated slopes is greater than or equal to the preset number of associated slopes, the slope creep coefficient is determined according to the associated slope radiation value and the sub-influence threshold;

[0081] If the number of associated slopes is less than the preset number of associated slopes, the slope creep coefficient is determined according to the sub-influence threshold.

[0082] If the number of associated slopes is greater than or equal to the preset number of associated slopes, the slope creep coefficient = associated slope radiation value + sub-influence threshold;

[0083] If the number of associated slopes is less than the preset number of associated slopes, the slope creep coefficient is positively correlated with the sub-influence threshold;

[0084] The other loess slopes outside the target loess slope are recorded as reference loess slopes, the associated slopes are the reference loess slopes whose shortest distance to the target loess slope is less than the preset shortest distance, and the number of associated slopes is the total number of associated slopes;

[0085] The radiation value of the associated slope is the maximum value of the sub-associated radiation values ​​corresponding to each associated slope. For a single associated slope, the associated slope is recorded as the target associated slope. The sub-associated radiation value corresponding to the target associated slope = (slope regeneration index corresponding to the target associated slope + slope angle turbulence coefficient corresponding to the target associated slope) / the shortest distance from the target associated slope to the target loess slope. The sub-influence threshold = slope regeneration index corresponding to the target loess slope + slope angle turbulence coefficient corresponding to the target loess slope.

[0086] The values ​​of the preset shortest distance and the preset number of associated slopes can be determined by the user according to the actual application scenario. The smaller the value of the preset shortest distance and the larger the value of the preset number of associated slopes, the greater the user's demand for determining the slope creep coefficient based on the sub-influence threshold. The values ​​of the preset shortest distance and the preset number of associated slopes are provided. The preset shortest distance is 1000m, and the preset number of associated slopes is 2.

[0087] Specifically, when the region selection unit performs monitoring and analysis on a single slope region, it determines whether to use the slope region as a monitoring region according to the radiation coefficient and the crack assessment value, wherein:

[0088] The slope area with a radiation coefficient less than the preset radiation coefficient or a crack assessment value greater than or equal to the preset crack assessment value is taken as the monitoring area.

[0089] Wherein, when determining whether to use the slope area as a monitoring area according to the radiation coefficient and the crack assessment value, if the radiation coefficient of a single slope area is greater than or equal to the preset radiation coefficient and the crack assessment value is less than the preset crack assessment value, the slope area will not be used as a monitoring area;

[0090] When monitoring and analyzing each slope area in descending order of collapse threshold, the slope area that has been monitored and analyzed and has been used as the monitoring area is recorded as the reference slope area, and the single slope area that is being monitored and analyzed is recorded as the target slope area; the radiation coefficient corresponding to the target slope area = radiation distance + radiation area, the radiation distance is the minimum value of the shortest distance from the center position corresponding to the target slope area to the center position corresponding to each reference slope area, and the radiation area is the area of ​​the target slope area overlapping with each reference slope area; the crack assessment value is the average value of the sub-influence coefficients corresponding to each crack in the target slope area; the collapse threshold is the average value of the crack thresholds corresponding to each crack in the single slope area;

[0091] The values ​​of the preset radiation coefficient and the preset crack assessment value can be determined by the user according to the actual application scenario. The greater the user's demand for improving the monitoring accuracy, the larger the value of the preset radiation coefficient, and the smaller the value of the preset crack assessment value. A value of the preset radiation coefficient and the preset crack assessment value is provided, and the historical records of setting the slope area as the monitoring area are detected. The average value of the radiation coefficient corresponding to the historical records that can meet the user's needs is recorded as the preset radiation coefficient, and the average value of the crack assessment values ​​corresponding to the historical records that can meet the user's needs is recorded as the preset crack assessment value.

[0092] Specifically, the monitoring determination unit responds to the setting conditions to determine the setting mode of the monitoring points in the monitoring area, wherein:

[0093] The setting condition for the response of the monitoring and judgment unit is that the crack dislocation coefficient is greater than or equal to the preset crack dislocation coefficient or the abnormal distribution threshold is greater than or equal to the preset abnormal distribution threshold, and the monitoring point setting method is to determine the processing method according to the correlation state of the abnormal sub-area;

[0094] The setting condition for the response of the monitoring and judgment unit is that the crack dislocation coefficient is less than the preset crack dislocation coefficient and the abnormal distribution threshold is less than the preset abnormal distribution threshold. The monitoring point setting method is to determine the optimization method according to the type of abnormal aggregation area.

[0095] The setting conditions include a first setting condition and a second setting condition. The first setting condition is that the crack dislocation coefficient is greater than or equal to a preset crack dislocation coefficient or the abnormal distribution threshold is greater than or equal to a preset abnormal distribution threshold. The second setting condition is that the crack dislocation coefficient is less than the preset crack dislocation coefficient and the abnormal distribution threshold is less than the preset abnormal distribution threshold.

[0096] The crack dislocation coefficient is the standard deviation of the dislocation heights corresponding to each crack in a single monitoring area. The dislocation height corresponding to a single crack is the absolute value of the difference between the height of the surface on one side of the crack and the minimum distance to the surface on the other side of the crack.

[0097] The abnormal distribution threshold is the average value of the reference distances corresponding to each abnormal sub-region in a single monitoring area. For a single abnormal sub-region, the abnormal sub-region is recorded as the target abnormal sub-region, and other abnormal sub-regions outside the target abnormal sub-region are recorded as reference abnormal sub-regions. The minimum value of the shortest distances from the target abnormal sub-region to each reference abnormal sub-region is recorded as the reference distance corresponding to the target abnormal sub-region. The abnormal sub-region includes a concave sub-region and a convex sub-region. The concave sub-region is a region that is concave compared with the surrounding area, and the convex sub-region is a region that is convex compared with the surrounding area.

[0098] Specifically, the first setting unit responds to the association status of the abnormal sub-area to determine the processing method, wherein:

[0099] The first setting unit responds to the first preset association state, and the processing method is to set monitoring points according to the associated abnormal area. When setting the monitoring points according to the associated abnormal area, the associated abnormal area is determined according to the abnormal distance, and the monitoring point is set at the central position corresponding to each associated abnormal area;

[0100] The first setting unit responds to the second preset association state, and the processing method is to set monitoring points at intervals between adjacent abnormal sub-regions. When the monitoring points are set at intervals between adjacent abnormal sub-regions, a monitoring point is set at the center position of each abnormal sub-region, and a monitoring point is set at a corresponding adjacent position of each abnormal sub-region;

[0101] The adjacent position corresponding to a single abnormal sub-region is the midpoint of a line connecting the center position of the abnormal sub-region and the center position of an adjacent abnormal sub-region corresponding to the abnormal sub-region.

[0102] Wherein, the associated abnormal area is determined according to the abnormal distance, wherein a combined analysis is performed on each abnormal sub-area in a single monitoring area, when the combined analysis is performed on a single abnormal sub-area, the abnormal sub-area is recorded as a target abnormal sub-area, and other abnormal sub-areas other than the target abnormal sub-area in the monitoring area that are not recorded in the abnormal combination are recorded as reference abnormal sub-areas, and the set of the reference abnormal sub-area whose abnormal distance with the target abnormal sub-area is less than a preset abnormal distance and the target abnormal sub-area is recorded as an abnormal combination, and the combined analysis is continued for the abnormal sub-areas that are not included in the abnormal combination until all abnormal sub-areas are recorded in the abnormal combination, then the combined analysis is stopped, and the minimum rectangle that can contain each abnormal sub-area in the single abnormal combination is recorded as an associated abnormal area;

[0103] The abnormal spacing is the shortest distance between the center points corresponding to any two abnormal sub-areas. The center point corresponding to a single abnormal sub-area is the center of the circumscribed circle corresponding to the abnormal sub-area. The value of the preset abnormal spacing can be determined by the user according to the actual application scenario. The greater the user's demand for improving monitoring accuracy, the smaller the value of the preset abnormal spacing. A value of the preset abnormal spacing is provided, and the historical records of setting monitoring points according to the associated abnormal areas are detected. The average value of the abnormal spacing corresponding to the historical records that can meet the user's needs is recorded as the preset abnormal spacing;

[0104] The center position corresponding to a single associated abnormal area is the center of the circumscribed circle of the associated abnormal area. The method for confirming adjacent abnormal sub-areas is as follows: for a single abnormal sub-area, the abnormal sub-area is recorded as the first target abnormal sub-area, and other abnormal sub-areas outside the first target abnormal sub-area in the monitoring area are recorded as the first reference abnormal sub-area. The adjacent abnormal sub-area corresponding to the first target abnormal sub-area is the first reference abnormal sub-area with the smallest abnormal distance to the first target abnormal sub-area; the center position of the adjacent abnormal sub-area corresponding to a single abnormal sub-area is the center of the circumscribed circle of the adjacent abnormal sub-area corresponding to the abnormal sub-area.

[0105] Specifically, the first setting unit determines the association state of the abnormal sub-region according to the abnormal variation coefficient and the associated abnormal region reference value, and the association state includes:

[0106] A first preset association state in which the abnormal variation coefficient is less than a preset abnormal variation coefficient or the associated abnormal region reference value is greater than or equal to a preset associated abnormal region reference value;

[0107] A second preset association state in which the abnormal variation coefficient is greater than or equal to a preset abnormal variation coefficient and the associated abnormal region reference value is less than a preset associated abnormal region reference value.

[0108] Among them, the abnormal variation coefficient is the average value of the associated abnormal distances corresponding to each associated abnormal area in a single monitoring area. For a single associated abnormal area, the associated abnormal area is recorded as the target associated abnormal area, and each associated abnormal area in the monitoring area except the target associated abnormal area is recorded as the reference associated abnormal area. The associated abnormal distance corresponding to the target associated abnormal area is the minimum value of the shortest distance from the center position corresponding to the target associated abnormal area to the center position corresponding to each reference associated abnormal area; the reference value of the associated abnormal area is the total amount of associated abnormal areas in a single monitoring area;

[0109] The values ​​of the preset abnormality variation coefficient and the preset associated abnormal area reference value can be determined by the user according to the actual application scenario. The larger the value of the preset abnormality variation coefficient and the smaller the value of the preset associated abnormal area reference value, the greater the user's need to set monitoring points according to the associated abnormal area. A value of the preset abnormality variation coefficient and the preset associated abnormal area reference value is provided, and the historical records of setting monitoring points according to the associated abnormal area are detected. The average value of the abnormality variation coefficient corresponding to the historical records that can meet the user's needs is recorded as the preset abnormality variation coefficient, and the average value of the associated abnormal area reference value corresponding to the historical records that can meet the user's needs is recorded as the preset associated abnormal area reference value.

[0110] Specifically, the second setting unit determines the type of abnormal aggregation area according to the convexity-concavity ratio, and the abnormal aggregation area type includes:

[0111] A type of abnormal clustering area whose convexity-concavity ratio is less than a preset convexity-concavity ratio, and a type of abnormal clustering area whose convexity-concavity ratio is greater than or equal to the preset convexity-concavity ratio.

[0112] Among them, the abnormal aggregation area is the minimum rectangle that can contain all the abnormal sub-areas in a single monitoring area. The convexity-concavity ratio is confirmed as follows: for a single abnormal aggregation area, the sum of the areas corresponding to the concave sub-areas in the abnormal aggregation area is recorded as M1, and the sum of the areas corresponding to the convex sub-areas in the monitoring area is recorded as M2. The convexity-concavity ratio = (the smaller value of M1 and M2) / (the larger value of M1 and M2). The value of the preset convexity-concavity ratio can be determined by the user according to the actual application scenario. The larger the value of the preset convexity-concavity ratio, the greater the user's need to determine the abnormal aggregation area as a type of abnormal aggregation area. A preset convexity-concavity ratio value is provided, and the preset convexity-concavity ratio is 50%.

[0113] Specifically, the second setting unit determines the optimization method according to the type of abnormal aggregation area, wherein:

[0114] For a type of abnormal clustering area, the optimization method is to determine the number of monitoring points in the abnormal clustering area according to the similarity of the abnormal sub-areas in the type of clustering area;

[0115] For the second type of abnormal aggregation area, the optimization method is to determine the number of monitoring points in the abnormal aggregation area according to the abnormal dislocation coefficient of the second type of aggregation area;

[0116] There is a negative correlation between the number of monitoring points and the similarity of the abnormal sub-regions in the first type of clustering area, and there is a positive correlation between the number of monitoring points and the abnormal dislocation coefficient of the second type of clustering area.

[0117] Among them, the confirmation method of the similarity of the abnormal sub-regions is as follows: if M1 is greater than M2 in the first type of abnormal aggregation area, the similarity of the abnormal sub-regions = 1 / (the maximum value of the depression coefficients corresponding to each concave sub-region in the first type of abnormal aggregation area - the minimum value of the depression coefficients corresponding to each concave sub-region in the first type of abnormal aggregation area); if M1 is less than M2 in the first type of abnormal aggregation area, the similarity of the abnormal sub-regions = 1 / (the maximum value of the convexity coefficients corresponding to each convex sub-region in the first type of abnormal aggregation area - the minimum value of the convexity coefficients corresponding to each convex sub-region in the first type of abnormal aggregation area); the depression coefficient corresponding to a single concave sub-region = the area of ​​the concave sub-region + the maximum depression depth, the convexity coefficient corresponding to a single convex sub-region = the area of ​​the convex sub-region + the maximum convex height; the maximum depression depth is the maximum depth of the depression of a single concave sub-region, and the maximum convex height is the maximum height of the convexity of a single convex sub-region; the abnormal dislocation coefficient = the standard deviation of the maximum depression depth corresponding to each concave sub-region in the second type of aggregation area + the standard deviation of the maximum convexity height corresponding to each convex sub-region in the second type of aggregation area.

[0118] So far, the technical solutions of the present invention have been described in conjunction with the preferred embodiments shown in the accompanying drawings. However, it is easy for those skilled in the art to understand that the protection scope of the present invention is obviously not limited to these specific embodiments. Without departing from the principle of the present invention, those skilled in the art can make equivalent changes or substitutions to the relevant technical features, and the technical solutions after these changes or substitutions will fall within the protection scope of the present invention.

[0119] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A loess slope multi-parameter environmental coupling monitoring system, characterized in that: include: A data collection unit, used to collect slope data; A state determination unit connected to the data acquisition unit and used to determine the state of the slope according to the slope regeneration index and the slope angle turbulence coefficient; A region division unit, which is connected to the data acquisition unit and the state determination unit respectively, and is used to determine a slope region division method according to the slope state, wherein the region division method is to perform associated division according to the crack influence coefficient and the distance reference value or to perform uniform division according to the slope creep coefficient; A region selection unit, which is connected to the region division unit, is used to perform monitoring analysis on each slope region in descending order of collapse thresholds to determine whether to use the slope region as a monitoring region; A monitoring determination unit connected to the area selection unit is used to determine a monitoring point setting method of the monitoring area according to the setting conditions, wherein the monitoring point setting method is to determine a processing method according to the association state of the abnormal sub-area or to determine an optimization method according to the type of abnormal aggregation area; A first setting unit, which is connected to the monitoring and determining unit, is used to determine the association state of the abnormal sub-region according to the abnormal variation coefficient and the reference value of the associated abnormal region, and determine a processing method according to the association state of the abnormal sub-region, and the processing method is to set monitoring points according to the associated abnormal region, or to set monitoring points at intervals between adjacent abnormal sub-regions; a second setting unit, which is connected to the monitoring and determining unit, and is used to determine the type of abnormal aggregation area according to the convex-concave ratio, and determine an optimization method according to the type of abnormal aggregation area, wherein the optimization method is to determine the number of monitoring points in the abnormal aggregation area according to the similarity of abnormal sub-areas in the first type of aggregation area or according to the abnormal dislocation coefficient of the second type of aggregation area; The area division unit responds to the slope state in a slope area division manner, wherein: The area division unit responds to the first slope state, and the area division method is to perform associated division according to the crack influence coefficient and the distance reference value. In the associated division, the minimum rectangular area that can contain each crack in a single associated crack combination is recorded as a slope area; The area division unit responds to the second slope state. The area division method is uniform division according to the slope creep coefficient. In the uniform division, the slope is divided into a number of rectangular areas with the same area and shape. Each rectangular area is recorded as a slope area. The number of slope areas is positively correlated with the slope creep coefficient. The slope creep coefficient is determined according to the number of associated slopes, where: If the number of associated slopes is greater than or equal to the preset number of associated slopes, the slope creep coefficient is determined according to the associated slope radiation value and the sub-influence threshold; If the number of associated slopes is less than the preset number of associated slopes, the slope creep coefficient is determined according to the sub-influence threshold; When the area selection unit performs monitoring and analysis on a single slope area, it determines whether to use the slope area as a monitoring area according to the radiation coefficient and the crack evaluation value, wherein: The slope area with a radiation coefficient less than the preset radiation coefficient or a crack assessment value greater than or equal to the preset crack assessment value is taken as the monitoring area; The monitoring determination unit responds to the setting conditions to determine the setting mode of the monitoring points in the monitoring area, wherein: The setting condition for the response of the monitoring and judgment unit is that the crack dislocation coefficient is greater than or equal to the preset crack dislocation coefficient or the abnormal distribution threshold is greater than or equal to the preset abnormal distribution threshold, and the monitoring point setting method is to determine the processing method according to the correlation state of the abnormal sub-area; The setting condition of the monitoring and judgment unit response is that the crack dislocation coefficient is less than the preset crack dislocation coefficient and the abnormal distribution threshold is less than the preset abnormal distribution threshold. The monitoring point setting method is to determine the optimization method according to the type of abnormal aggregation area; Slope regeneration index = number of sliding surfaces + overlapping area of ​​sliding surfaces. The number of sliding surfaces is the total number of sliding surfaces in the target loess slope. The overlapping area of ​​sliding surfaces is the average value of the reference areas corresponding to each sliding surface. For a single sliding surface, the sliding surface is recorded as the target sliding surface, and the other sliding surfaces except the target sliding surface are recorded as reference sliding surfaces. The reference area corresponding to the target sliding surface is the maximum value of the overlapping areas of the target sliding surface and each reference sliding surface. The slope angle turbulence coefficient is the standard deviation of the slope angle corresponding to each slope surface; The method for confirming the crack influence coefficient is that, for any two cracks, the larger value of the sub-influence coefficients corresponding to the two cracks is recorded as A1, and the smaller value is recorded as A2. The crack influence coefficient = (A1-A2) / A1, The sub-influence coefficient is confirmed as follows: If the impact angle is less than the preset impact angle, the sub-influence coefficient = crack threshold - impact angle; If the influence angle is greater than or equal to the preset influence angle, the sub-influence coefficient is positively correlated with the crack threshold.

2. The loess slope multi-parameter environmental coupling monitoring system according to claim 1 is characterized in that: The state determination unit determines the slope state according to the slope regeneration index and the slope angle turbulence coefficient. The slope state includes: A first slope state in which a slope regeneration index is greater than or equal to a preset slope regeneration index or a slope angle turbulence coefficient is greater than or equal to a preset slope angle turbulence coefficient; A second slope state in which the slope regeneration index is less than a preset slope regeneration index and the slope angle turbulence coefficient is less than a preset slope angle turbulence coefficient.

3. The loess slope multi-parameter environmental coupling monitoring system according to claim 1 is characterized in that: The first setting unit responds to the association status of the abnormal sub-area to determine the processing method, wherein: The first setting unit responds to the first preset association state, and the processing method is to set monitoring points according to the associated abnormal area. When setting the monitoring points according to the associated abnormal area, the associated abnormal area is determined according to the abnormal distance, and the monitoring point is set at the central position corresponding to each associated abnormal area; The first setting unit responds to the second preset association state, and the processing method is to set monitoring points at intervals between adjacent abnormal sub-regions. When the monitoring points are set at intervals between adjacent abnormal sub-regions, a monitoring point is set at the center position of each abnormal sub-region, and a monitoring point is set at a corresponding adjacent position of each abnormal sub-region; The adjacent position corresponding to a single abnormal sub-region is the midpoint of a line connecting the center position of the abnormal sub-region and the center position of an adjacent abnormal sub-region corresponding to the abnormal sub-region.

4. The loess slope multi-parameter environmental coupling monitoring system according to claim 3 is characterized in that: The first setting unit determines the association state of the abnormal sub-region according to the abnormal variation coefficient and the associated abnormal region reference value, and the association state includes: A first preset association state in which the abnormal variation coefficient is less than a preset abnormal variation coefficient or the associated abnormal region reference value is greater than or equal to a preset associated abnormal region reference value; A second preset association state in which the abnormal variation coefficient is greater than or equal to a preset abnormal variation coefficient and the associated abnormal region reference value is less than a preset associated abnormal region reference value.

5. The loess slope multi-parameter environmental coupling monitoring system according to claim 1 is characterized in that: The second setting unit determines the type of abnormal aggregation area according to the convexity-concavity ratio, and the type of abnormal aggregation area includes: A type of abnormal clustering area whose convexity-concavity ratio is less than a preset convexity-concavity ratio, and a type of abnormal clustering area whose convexity-concavity ratio is greater than or equal to the preset convexity-concavity ratio.

6. The loess slope multi-parameter environmental coupling monitoring system according to claim 5 is characterized in that: The second setting unit determines the optimization method according to the type of abnormal aggregation area, wherein: For a type of abnormal clustering area, the optimization method is to determine the number of monitoring points in the abnormal clustering area according to the similarity of the abnormal sub-areas in the type of clustering area; For the second type of abnormal aggregation area, the optimization method is to determine the number of monitoring points in the abnormal aggregation area according to the abnormal dislocation coefficient of the second type of aggregation area; There is a negative correlation between the number of monitoring points and the similarity of the abnormal sub-regions in the first type of clustering area, and there is a positive correlation between the number of monitoring points and the abnormal dislocation coefficient of the second type of clustering area.

Images