阅读说明：本技术 一种高烈度地震区滑坡震后稳定性评估方法及系统 (Method and system for evaluating stability after landslide in high-intensity seismic area ) 是由 王珣 刘勇 伏坤 高柏松 陈伟志 徐鑫 袁焦 裴起帆 潘兆马 杨科 杨学锋 姚 于 2021-09-28 设计创作，主要内容包括：本发明公开了一种高烈度地震区滑坡震后稳定性评估方法及系统,属于抗滑桩长期服役性能监测技术领域。包括：在滑坡的不同区域布设若干个加速度传感器；所述加速度传感器采集地震波的加速度数据,并将所述加速度数据发送至数据采集及传输系统,所述数据采集及传输系统再将触发阈值后的加速度数据发送至监测中心；所述监测中心获取到触发阈值后的加速度数据后,对所述加速度数据进行自适应卡尔曼滤波降噪处理；所述监测中心再对降噪处理后的所述加速度数据进行解算,确定坡体损伤位置,并进一步推导计算得出坡体稳定性系数降低比例,得到损伤识别结果。具有操作简便、成本低、可靠性高等优点,在高烈度频发地震区有着广阔的应用前景。(The invention discloses a method and a system for evaluating stability after landslide in a high-intensity seismic area and belongs to the technical field of long-term service performance monitoring of slide-resistant piles. The method comprises the following steps: arranging a plurality of acceleration sensors in different areas of a landslide; the acceleration sensor collects acceleration data of seismic waves and sends the acceleration data to the data collection and transmission system, and the data collection and transmission system sends the acceleration data after triggering a threshold value to the monitoring center; after the monitoring center acquires the acceleration data after the triggering threshold value, carrying out self-adaptive Kalman filtering noise reduction processing on the acceleration data; and the monitoring center then calculates the acceleration data after noise reduction processing to determine the damage position of the slope body, and further deduces and calculates to obtain the reduction ratio of the stability coefficient of the slope body, so as to obtain a damage identification result. The method has the advantages of simple and convenient operation, low cost, high reliability and the like, and has wide application prospect in high-intensity frequent earthquake regions.)

1. A method for evaluating stability after landslide in a high-intensity seismic area is characterized by comprising the following steps:

s1: arranging a plurality of acceleration sensors in different areas of a landslide;

s2: the acceleration sensor collects acceleration data of seismic waves and sends the acceleration data to the data collection and transmission system, and the data collection and transmission system sends the acceleration data after triggering a threshold value to the monitoring center;

s3: after acquiring the acceleration data after the trigger threshold, the monitoring center performs adaptive Kalman filtering noise reduction processing on the acceleration data after the trigger threshold;

s4: and the monitoring center then calculates the acceleration data after noise reduction to obtain the natural vibration frequency f of each monitoring point position, determines the damage position of the slope body according to the reduction degree of the natural vibration frequency, and further deduces and calculates to obtain the reduction proportion of the stability coefficient of the slope body to obtain a damage identification result.

2. The method for evaluating the stability after landslide in the high intensity seismic region according to claim 1, wherein the step S1 comprises: and comprehensively determining the layout position of the acceleration sensor according to the landslide scale, the geological conditions and the potential disaster mode.

3. The method for evaluating the stability after landslide in a high intensity seismic region according to claim 1, wherein the step S2 further comprises capturing the arrival time of the seismic P wave by a long-short time averaging method by the acceleration sensor, namely: STA/LTA, wherein:

in the formula, STA_iFor short-time average of the acceleration data at time i, LTA_iIs the long-time average value at time i, CF (i) is the characteristic function value of the acceleration data at time i, N_staNumber of recorded points included in time window as short-time average value, N_ltaThe number of recorded points included in the long time average time window, and x (i) the original signal amplitude at time i.

4. The method as claimed in claim 3, wherein the threshold is set according to environmental noise, the value obtained by the STA/LTA recursive calculation is compared with the threshold, and the time when the first time is greater than the threshold is determined as the time when the seismic P-wave arrives at the location.

5. The method for evaluating the stability after landslide in a high intensity seismic region according to claim 1, wherein the step S3 comprises: the theoretical output acceleration of the acceleration sensor is selected as a one-dimensional state quantity to be estimated, the system state is considered to be stable except for system noise, and a signal noise reduction model of the acceleration sensor is as follows:

X_k+1＝X_k+w_k

in the formula, g₀Is the acceleration of gravity; k is the accelerometer voltage scale factor, X_kThe state quantity of the system at the moment k is obtained; w is a_kIs the system noise; v. of_k+1To measure noise; z_k+1The measured value is the system measured value at the moment k + 1.

6. The method for evaluating stability after landslide in a high intensity seismic region according to claim 1, wherein the step S4 further comprises:

s41: the self-oscillation frequency f of each point position in the slope body is obtained through a transfer function method, the transfer function is a complex variable function taking the frequency omega as an independent variable, and the function structure is as follows:

H_a(ω,z_j)＝R(ω)+I(ω)i

in the formula, H_a(ω,z_j) The transfer function is adopted, R (omega) is the real frequency characteristic of the transfer function, I (omega) is the imaginary frequency characteristic of the transfer function, and I is an imaginary number unit;

s42: carrying out statistical analysis on the natural vibration frequency f of each monitored part after the earthquake, and judging the damage degree of each part of the slope body according to the descending proportion of the natural vibration frequency f;

s43: analyzing the stability factor F of the landslide before earthquake_sStability coefficient F of weak stable stage after earthquake_s,dThe stability coefficient is calculated by the following equation, at this time, regardless of the frictional force action:

according to the intensity reduction method, c and phi simultaneously divide the reduction coefficient w:

c_d＝cw

tanφ_d＝wtanφ

wherein, the viscous equivalent stiffness is introduced, the change of the natural vibration frequency f can lead to the reduction of the cohesive force c, and the reduction relation is as follows:

wherein T' and T are respectively anti-slip force and downward-slip force, W is the self weight of the slope body, theta is the break angle of the slip body, c and phi are respectively cohesive force and internal friction angle, W is the reduction coefficient, k is the rigidity, m is the mass, L is the length of the slip surface, phi is the length of the slip surface_dAnd phi₀The internal friction angles of the slope before and after an earthquake are respectively;

defining the reduction coefficient R of the stability index of the landslide after earthquake as follows:

R＝(1-λ)²

in the formula, R is the stability index reduction coefficient of the landslide after the earthquake, and lambda is the natural vibration frequency F reduction rate at the monitoring point position of the landslide after the earthquake, and F_sAnd F_s,dStability factor, f, of landslides before and after an earthquake₀And f_dThe self-vibration frequencies of the landslide mass before and after the earthquake are respectively.

7. The method for evaluating stability after landslide in a high intensity seismic region according to claim 6, wherein said step S41 further comprises the expression of said transfer function being:

in the formula, G_XX(ω) is the self-power spectrum of the acceleration response at the foot of the landslide; g_XY(ω，z_j) A cross power spectrum of the acceleration response of the slope foot of the landslide and the acceleration response of each point position in the slope body; z is a radical of_jResponding to the acceleration of each point position in the slope body; h_a(ω,z_j) Each peak value of (a) corresponds to ω_kApproximately equal to the natural frequency f of each order of the measuring point.

8. The evaluation system of the stability evaluation method after landslide in the high intensity seismic region according to the claims 1-7, comprising a plurality of acceleration sensors, a data acquisition and transmission system and a monitoring center;

an acceleration sensor: the acceleration sensor is in communication connection with the data acquisition and transmission system and is used for acquiring acceleration data of seismic waves and transmitting the acceleration data to the data acquisition and transmission system;

data acquisition and transmission system: the data acquisition and transmission system is in communication connection with the monitoring center and is used for transmitting the acceleration data after triggering the threshold value to the monitoring center;

a monitoring center: the system is used for carrying out self-adaptive Kalman filtering denoising processing on the received acceleration data after the triggering threshold, resolving the denoising processed acceleration data, obtaining the natural vibration frequency f of each monitoring point, determining the damage position of the slope body according to the reduction degree of the natural vibration frequency, and further deducing and calculating to obtain the reduction ratio of the stability coefficient of the slope body.

9. The high intensity seismic area landslide post-shock stability assessment system according to claim 8, wherein said acceleration sensors are mounted in the landslide principal axis direction and at the landslide foot.

Technical Field

The invention relates to the technical field of monitoring of long-term service performance of slide-resistant piles, in particular to a method and a system for evaluating stability after landslide in a high-intensity seismic area.

Background

At present, the landslide stability analysis under high-intensity seismic load is mostly carried out by adopting a numerical simulation mode, such as document 1: patent "a method for predicting three-dimensional slope stability under earthquake load" (CN103135128B, zhou xiao ping), document 2: patent "a method for predicting landslide occurrence time of three-dimensional seismic slope based on sliding displacement analysis" (CN104182646B, zhou xiao ping), document 3: patent "a novel layer slope earthquake damage imbalance calculation statistics and method based on energy demodulation frequency spectrum" (CN108445536A, lugge).

However, there is no method for long-term online monitoring and evaluation of landslide stability in a high-intensity seismic area, where the high-intensity seismic area is an area where faults are active, a large earthquake may occur, the intensity is 7 degrees or more, and the requirement for earthquake prevention is high. Dynamic equation (K-omega) before damage by undamped n-degree-of-freedom system²M) phi is 0, the rigidity of the system is reduced if the anti-slide pile is locally damaged after an earthquake, and the rigidity of the system is in direct proportion to the natural vibration frequency of the system.

Disclosure of Invention

The invention aims to overcome the defect that no method for long-term online monitoring and evaluation of landslide stability in a high-intensity seismic area exists in the prior art, and provides a system for evaluating the stability after landslide in the high-intensity seismic area.

In order to achieve the above purpose, the invention provides the following technical scheme:

a method for evaluating stability after landslide in a high-intensity seismic area comprises the following steps:

s1: arranging a plurality of acceleration sensors in different areas of a landslide;

s2: the acceleration sensor collects acceleration data of seismic waves and sends the acceleration data to the data collection and transmission system, and the data collection and transmission system sends the acceleration data after triggering a threshold value to the monitoring center;

s3: after acquiring the acceleration data after the trigger threshold, the monitoring center performs adaptive Kalman filtering noise reduction processing on the acceleration data after the trigger threshold;

s4: and the monitoring center then calculates the acceleration data after noise reduction to obtain the natural vibration frequency f of each monitoring point position, determines the damage position of the slope body according to the reduction degree of the natural vibration frequency, and further deduces and calculates to obtain the reduction proportion of the stability coefficient of the slope body to obtain a damage identification result.

By adopting the technical scheme, the dynamic response of each point position in the landslide body is monitored, the self-vibration frequency of the point position is obtained, the damage state of the landslide body after earthquake is judged, and the stability of the landslide body after earthquake is evaluated.

As a preferable aspect of the present invention, step S1 includes: the method comprises the steps of comprehensively determining the arrangement positions of acceleration sensors according to the landslide scale, geological conditions and potential disaster modes, wherein the acceleration sensors are arranged in a slope body above a bedrock surface at equal intervals, a plurality of deep holes are generally designed along the direction of a main shaft, the acceleration sensors are arranged, and 1 acceleration sensor is arranged at a stable point near a slope toe.

As a preferable embodiment of the present invention, the step S2 further includes capturing, by the acceleration sensor, an arrival time of the seismic P-wave by a long-short time averaging method, that is: STA/LTA, wherein:

in the formula, STA_iFor short-time average of the acceleration data at time i, LTA_iIs the long-time average value at time i, CF (i) is the characteristic function value of the acceleration data at time i, N_staNumber of recorded points included in time window as short-time average value, N_ltaThe number of recorded points included in the long time average time window, and x (i) the original signal amplitude at time i.

In a preferred embodiment of the present invention, the Threshold (THR) is set to 4 to 10, and the time when the value obtained by the STA/LTA recursive calculation is first greater than the threshold is compared with the threshold requirement (THR), and it is determined that the seismic P-wave has arrived at the location.

As a preferable embodiment of the present invention, the step S3 includes: the theoretical output acceleration of the acceleration sensor is selected as a one-dimensional state quantity to be estimated, the system state is considered to be kept stable except for system noise, the direct measurement is the output voltage of the accelerometer containing measurement noise, and the conversion of the theoretical value of the acceleration and the voltage measurement value can be realized through the scale factor of the accelerometer, so that the signal noise reduction model of the acceleration sensor is as follows:

X_k+1＝X_k+w_k

in the formula, g₀Is the acceleration of gravity; k is the accelerometer voltage scale factor, X_kThe state quantity of the system at the moment k is obtained; w is a_kIs the system noise; v. of_k+1To measure noise; z_k+1The measured value is the system measured value at the moment k + 1.

As a preferable embodiment of the present invention, the step S4 further includes:

s41: the self-oscillation frequency f of each point position in the slope body is obtained through a transfer function method, the transfer function is a complex variable function taking the frequency omega as an independent variable, and the function structure is as follows:

H_a(ω,z_j)＝R(ω)+I(ω)i

in the formula, H_a(ω,z_j) The transfer function is adopted, R (omega) is the real frequency characteristic of the transfer function, I (omega) is the imaginary frequency characteristic of the transfer function, and I is an imaginary number unit;

the principle is that the acceleration response of the stable point of the landslide foot and the acceleration response of each point in the landslide body are subjected to cross power spectrum analysis, the acceleration response of the stable point of the landslide foot is subjected to self power spectrum analysis, and a transfer function, namely an absolute transfer function, at the measured point is obtained;

s42: carrying out statistical analysis on the natural vibration frequency f of each monitored part after the earthquake, and judging the damage degree of each part of the slope body according to the descending proportion of the natural vibration frequency f;

s43: analyzing the stability factor F of the landslide before earthquake_sStability coefficient F of weak stable stage after earthquake_s,dThe stability coefficient is calculated by the following equation, at this time, regardless of the frictional force action:

according to the intensity reduction method, c and phi simultaneously divide the reduction coefficient w:

c_d＝cw

tanφ_d＝wtanφ

wherein, the viscous equivalent stiffness is introduced, the change of the natural vibration frequency f can lead to the reduction of the cohesive force c, and the reduction relation is as follows:

wherein T' and T are respectively anti-slip force and downward-slip force, W is the self weight of the slope body, theta is the break angle of the slip body, c and phi are respectively cohesive force and internal friction angle, W is the reduction coefficient, k is the rigidity, m is the mass, L is the length of the slip surface, phi is the length of the slip surface_dAnd phi₀The internal friction angles of the slope before and after an earthquake are respectively;

defining the reduction coefficient R of the stability index of the landslide after earthquake as follows:

R＝(1-λ)²

in the formula, R is the stability index reduction coefficient of the landslide after the earthquake, and lambda is the natural vibration frequency F reduction rate at the monitoring point position of the landslide after the earthquake, and F_sAnd F_s,dStability factor, f, of landslides before and after an earthquake₀And f_dThe self-vibration frequencies of the landslide mass before and after the earthquake are respectively.

As a preferable embodiment of the present invention, the step S41 further includes that the transfer function expression is:

in the formula, G_XX(ω) is the self-power spectrum of the acceleration response at the foot of the landslide; g_XY(ω，z_j) A cross power spectrum of the acceleration response of the slope foot of the landslide and the acceleration response of each point position in the slope body; z is a radical of_jResponding to the acceleration of each point position in the slope body; wherein H_a(ω,z_j) Each peak value of (a) corresponds to ω_kApproximately equal to the natural frequency f of each order of the measuring point.

The evaluation system of the stability evaluation method after landslide in the high-intensity seismic area comprises a plurality of acceleration sensors, a data acquisition and transmission system and a monitoring center;

an acceleration sensor: the acceleration sensor is in communication connection with the data acquisition and transmission system and is used for acquiring acceleration data of seismic waves and transmitting the acceleration data to the data acquisition and transmission system;

data acquisition and transmission system: the data acquisition and transmission system is in communication connection with the monitoring center and is used for transmitting the acceleration data after triggering the threshold value to the monitoring center;

a monitoring center: the system is used for carrying out self-adaptive Kalman filtering denoising processing on the received acceleration data after the triggering threshold, resolving the denoising processed acceleration data, obtaining the natural vibration frequency f of each monitoring point, determining the damage position of the slope body according to the reduction degree of the natural vibration frequency, and further deducing and calculating to obtain the reduction ratio of the stability coefficient of the slope body.

By adopting the technical scheme, the method has the advantages of simplicity and convenience in operation, low cost, high reliability and the like, and has a wide application prospect in high-intensity frequent earthquake regions.

As a preferable aspect of the present invention, the acceleration sensor is installed in a direction of a main axis of the landslide and at a foot of the landslide.

Compared with the prior art, the invention has the beneficial effects that: the method has the advantages of simplicity and convenience in operation, low cost, high reliability and the like, and has wide application prospect in high-intensity frequent earthquake regions.

Drawings

Fig. 1 is a flowchart of a method for evaluating stability after landslide in a high-intensity seismic area according to embodiment 1 of the present invention;

FIG. 2 is a cross-sectional view of the arrangement of acceleration sensors in the method for evaluating the stability after landslide in a high-intensity seismic region according to embodiment 1 of the present invention;

FIG. 3 is a diagram illustrating the change of the natural frequency of the landslide body after loading seismic waves of different seismic magnitudes in the high-intensity seismic region landslide post-earthquake stability evaluation method according to embodiment 1 of the present invention;

FIG. 4 is a graph illustrating the degradation of the stability of a landslide body under different dynamic peak acceleration seismic conditions according to the method for evaluating the stability after landslide in a high-intensity seismic region in the embodiment 1 of the present invention;

fig. 5 is a structural diagram of a system for evaluating stability after landslide in a high-intensity seismic area according to this embodiment 2.

Detailed Description

The present invention will be described in further detail with reference to test examples and specific embodiments. It should be understood that the scope of the above-described subject matter is not limited to the following examples, and any techniques implemented based on the disclosure of the present invention are within the scope of the present invention.

Example 1

A method for evaluating stability after landslide in a high-intensity seismic region is disclosed as shown in FIG. 1, and comprises the following steps:

s1: arranging a plurality of acceleration sensors in different areas of a landslide;

s2: the acceleration sensor collects acceleration data of seismic waves and sends the acceleration data to the data collection and transmission system, and the data collection and transmission system sends the acceleration data after triggering a threshold value to the monitoring center;

s3: after acquiring the acceleration data after the trigger threshold, the monitoring center performs adaptive Kalman filtering noise reduction processing on the acceleration data after the trigger threshold;

s4: and the monitoring center then calculates the acceleration data after noise reduction to obtain the natural vibration frequency f of each monitoring point position, determines the damage position of the slope body according to the reduction degree of the natural vibration frequency, and further deduces and calculates to obtain the reduction proportion of the stability coefficient of the slope body to obtain a damage identification result.

Step S1 includes: the method comprises the steps of comprehensively determining the arrangement positions of acceleration sensors according to the landslide scale size, geological conditions and potential disaster modes, arranging the acceleration sensors in a slope body above a bedrock surface at equal intervals, designing a plurality of deep holes along the main shaft direction to arrange the acceleration sensors, arranging 1 acceleration sensor at a stable point near a slope toe, arranging 9 acceleration sensors J1-J9 in the slope body of a certain frequent earthquake area as shown in FIG. 2, and assuming that 7 times of earthquakes with different dynamic peak accelerations of 0.1 g-0.7 g occur successively in the landslide in recent years.

The method for determining whether the interior is damaged or not and the stability reduction condition after landslide and earthquake comprises the following specific steps:

the step S2 further includes that the acceleration sensor captures the arrival time of the seismic P-wave by a long-short time averaging method, that is: STA/LTA, wherein:

in the formula, STA_iFor short-time average of the acceleration data at time i, LTA_iIs the long-time average value at time i, CF (i) is the characteristic function value of the acceleration data at time i, N_staNumber of recorded points included in time window as short-time average value, N_ltaThe number of recorded points included in the long time average time window, and x (i) the original signal amplitude at time i, where the short window takes 0.2S and the long window takes 10S.

And setting the Threshold (THR) according to the environmental noise, comparing the value obtained by the STA/LTA recursive calculation with the Threshold (THR) when the Threshold (THR) is set to be 8, and judging that the earthquake P wave reaches the time when the time is greater than the threshold for the first time.

The step S3 includes: the theoretical output acceleration of the acceleration sensor is selected as a one-dimensional state quantity to be estimated, the system state is considered to be kept stable except for system noise, the direct measurement is the output voltage of the accelerometer containing measurement noise, and the conversion of the theoretical value of the acceleration and the voltage measurement value can be realized through the scale factor of the accelerometer, so that the signal noise reduction model of the acceleration sensor is as follows:

X_k+1＝X_k+w_k

in the formula, g₀Is the acceleration of gravity; k is the accelerometer voltage scale factor, X_kThe state quantity of the system at the moment k is obtained; w is a_kIs the system noise; v. of_k+1To measure noise; z_k+1The measured value is the system measured value at the moment k + 1.

The step S4 further includes:

s41: obtaining the natural vibration frequency f of each point position in the slope body by a transfer function method₀As shown in fig. 3, taking the monitoring point of J4 as an example, the natural frequency f of the slip mass decreases to a certain extent as the magnitude of the vibration increases, which indicates that the internal damage of the slip mass occurs, the transfer function is a complex function with the frequency ω as an independent variable, and the function structure is:

H_a(ω,z_j)＝R(ω)+I(ω)i

in the formula, H_a(ω,z_j) The transfer function is adopted, R (omega) is the real frequency characteristic of the transfer function, I (omega) is the imaginary frequency characteristic of the transfer function, and I is an imaginary number unit;

the principle is that the acceleration response of the stable point of the landslide foot and the acceleration response of each point in the landslide body are subjected to cross power spectrum analysis, the acceleration response of the stable point of the landslide foot is subjected to self power spectrum analysis, and a transfer function, namely an absolute transfer function, at the measured point is obtained;

s42: carrying out statistical analysis on the natural vibration frequency f of each monitored part after the earthquake, and judging the damage degree of each part of the slope body according to the descending proportion of the natural vibration frequency f;

s43: analyzing the stability factor F of the landslide before earthquake_sStability coefficient F of weak stable stage after earthquake_s,dThe stability coefficient is calculated by the following equation, at this time, regardless of the frictional force action:

according to the intensity reduction method, c and phi simultaneously divide the reduction coefficient w:

c_d＝cw

tanφ_d＝wtanφ

wherein, the viscous equivalent stiffness is introduced, the change of the natural vibration frequency f can lead to the reduction of the cohesive force c, and the reduction relation is as follows:

wherein T' and T are respectively anti-slip force and downward-slip force, W is the self weight of the slope body, theta is the break angle of the slip body, c and phi are respectively cohesive force and internal friction angle, W is the reduction coefficient, k is the rigidity, m is the mass, L is the length of the slip surface, phi is the length of the slip surface_dAnd phi₀The internal friction angles of the slope before and after an earthquake are respectively;

defining the reduction coefficient R of the stability index of the landslide after earthquake as follows:

R＝(1-λ)²

in the formula, R is the stability index reduction coefficient of the landslide after the earthquake, and lambda is the natural vibration frequency F reduction rate at the monitoring point position of the landslide after the earthquake, and F_sAnd F_s,dStability factor, f, of landslides before and after an earthquake₀And f_dThe self-vibration frequencies of the landslide mass before and after the earthquake are respectively.

Using the reduction factor index R ═ (1-lambda)²And calculating the stability coefficient of the landslide after the seismic waves of different levels are loaded, and assuming that the stability coefficient of the landslide before the seismic waves is 1.1, and the slope body is in a basic stable state. According to the J4 point monitoring result and the calculation result of the landslide stability reduction coefficient after earthquake as shown in FIG. 4, it can be seen that the reduction coefficient is continuously increased along with the increase of the earthquake grade, namely the landslide stability coefficient is continuously reduced, the final reduction coefficient reaches 85.6%, and the landslide stability coefficient is reduced after the earthquake in the area where the monitoring point is located for several times in different gradesTo 1.1 × 85.6% ═ 0.95. And further calculating the reduction coefficient of each monitored point position, substituting a stability number calculation formula into the calculation to finally calculate the stability coefficient to be 1.03, and judging that the landslide body is in an under-stable state.

The step S41 further includes that the transfer function expression is:

in the formula, G_XX(ω) is the self-power spectrum of the acceleration response at the foot of the landslide; g_XY(ω，z_j) A cross power spectrum of the acceleration response of the slope foot of the landslide and the acceleration response of each point position in the slope body; z is a radical of_jResponding to the acceleration of each point position in the slope body; wherein H_a(ω,z_j) Each peak value of (a) corresponds to ω_kApproximately equal to the natural frequency f of each order of the measuring point.

And finally, carrying out next treatment scheme design aiming at the ground according to the damage identification result of the landslide body and combining with field investigation.

Example 2

The evaluation system of the stability evaluation method after landslide in the high-intensity seismic area comprises a plurality of acceleration sensors, a data acquisition and transmission system and a monitoring center, wherein the acceleration sensors are connected with the monitoring center;

an acceleration sensor: the acceleration sensor is in communication connection with the data acquisition and transmission system and is used for acquiring acceleration data of seismic waves and transmitting the acceleration data to the data acquisition and transmission system;

data acquisition and transmission system: the data acquisition and transmission system is in communication connection with the monitoring center and is used for transmitting the acceleration data after triggering the threshold value to the monitoring center;

a monitoring center: the system is used for carrying out self-adaptive Kalman filtering denoising processing on the received acceleration data after the triggering threshold, resolving the denoising processed acceleration data, obtaining the natural vibration frequency f of each monitoring point, determining the damage position of the slope body according to the reduction degree of the natural vibration frequency, and further deducing and calculating to obtain the reduction ratio of the stability coefficient of the slope body.

The acceleration sensor is installed in the direction of a main shaft of the landslide and at the position of a foot of the landslide.

The method has the advantages of simplicity and convenience in operation, low cost, high reliability and the like, and has wide application prospect in high-intensity frequent earthquake regions.

The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention.