TWI475244B - Slope displacement detection method - Google Patents

Description

Slope displacement detection method

The invention relates to a method for detecting the displacement of a slope, in particular to a method for detecting the displacement of a slope suitable for monitoring the surface change of the slope, the displacement and the change of the groundwater level, the change of the sliding surface of the slope, and the like.

The situation of earth collapse has been heard, especially in the face of a gradual change in the climate system, which has caused amazing rainfall, whether it is the occurrence of earth collapse or earth-rock flow. Among them, for a slope with a potential collapse risk, it is located at a location where there is no use value, and usually collapses or not does not cause too much loss. However, if the collapse occurs in the catchment area, there will be a large amount of earth and stone deposits, resulting in a decrease in the life of the reservoir and a decrease in the amount of water stored. In particular, if such a slope is located above the house, the danger is the loss of people's lives and property.

Taiwan’s metropolitan area is densely populated, and excessive use of land is generally regarded as one of the culprit of land disasters. In the construction of residential areas on hillsides, in addition to the strict control of public sector units, it is very important for builders to understand geological conditions. Among them, the design of the construction method and the overall consideration of disaster prevention and mitigation facilities are all important factors to preserve the safety of future residents. At present, the construction of a large number of residential buildings in Taiwan is a high-risk design. Earth's collapse or earth-rock flow caused by sudden showers during the typhoon season often endangers the lives and property of residents.

However, in the conventional measurement and monitoring methods, it is often only possible to measure and monitor the unique parameters, such as using a GPS receiver for displacement monitoring, or using ground resistance measurement for analysis. Moreover, conventional physical measurement monitoring points have the following disadvantages such as costly, vulnerable, and under-represented monitoring bottlenecks, and conventionally, blind spots with difficulty in grasping accuracy and manpower are required. It can be seen that how to achieve an accurate and detailed, fast measurement, low cost, manpower saving, and more integrated measurement methods to cross-check the slope detection method is really an urgent industry. need.

The main object of the present invention is to provide a method for detecting the displacement of a slope, which can integrate a plurality of measurement methods to perform cross-match analysis, and for surface changes, displacements, groundwater level changes, slope sliding surface changes, etc. of dangerous slope terrain. Monitoring is performed to provide more accurate and detailed measurement results, and the invention has the advantages of rapid measurement, low cost, and labor saving.

In order to achieve the above object, the present invention is a slope displacement detecting method, comprising the following steps: step (A) respectively setting a complex conjugate scanning sphere, and a plurality of GPS receivers at different positions on one side of the slope, to a 3D mine The scanner scans the slope and covers the complex conjugate scan ball, while the 3D laser scanner generates a first 3D numerical terrain data; the complex GPS receiver calculates its coordinates and height, and stores it as A first GPS data. Then, after step (B) is separated by a specific time, the 3D laser scanner scans the slope again and covers the complex conjugate scan ball, and the 3D laser scanner generates a second 3D numerical terrain data; The GPS receiver again calculates its coordinates and height again and stores it as a second GPS data.

Then, step (C) nests and analyzes the first 3D numerical terrain data and the second 3D numerical terrain data, thereby obtaining at least one terrain variation region; and comparing and analyzing the first GPS data and the second GPS data to obtain at least one standard a height variation point; step (D) comparing at least one terrain variation region with at least one standard height variation point; and step (E), when the alignment is matched, routing the plurality of electrodes to perform at least one terrain variation region for resistance measurement Approent Resistivity of at least one terrain variation region.

The first 3D numerical terrain data may include a first coordinate system and a first point cloud digital model, and the second 3D numerical terrain data may include a second coordinate system and a second point cloud digital model. The nested analysis of the first 3D numerical topographic data and the second 3D numerical topographical data in step (C) may include the following steps: (C1) comparing the first 3D numerical topographic data with the second 3D numerical topographical data in the complex conjugate sweep The position of the relative relationship of the aiming ball changes, and the point where the relative relationship position is not changed is taken as a target conversion conjugate point; (C2) the common point when the coordinate conversion conjugate point is used as the coordinate system for conversion, the second The second coordinate system of the 3D numerical terrain data is converted into the first coordinate system of the first 3D numerical terrain data; the first point cloud digital model of the first 3D numerical terrain data, and the second point cloud digit of the second 3D numerical terrain data The model is transformed into a first terrain numerical model (DEM) and a second topographic numerical model (DEM) using the irregular triangulation rule (TIN); (C3) the second topographic numerical model (DEM) is imported into the first coordinate system. And calculating the distance between the vertices of all triangular meshes in the first topographic numerical model and the vertical axis direction of the vertices of the vertices and the coordinates of the coordinates of the second topographic numerical model; and (C4) analyzing the equidistances, wherein the greater than the average distance The area formed by the quantity is at least one terrain variation area. In addition, the step (C3) further comprises using an irregular triangulation method to interpolate the equidistant distance into a numerical model of the variation amount, as the analysis and presentation of the difference between the two models.

In addition, the first GPS data in the step (A) may further include a first reference coordinate of one of the satellite receiving stations, and the second GPS data in the step (B) may further include a second reference coordinate of one of the satellite receiving stations; The comparison between the first GPS data and the second GPS data in step (C) is based on the comparison between the first reference coordinate and the second reference coordinate. The purpose of using the satellite receiving station as the reference coordinate is to measure the stable and high-precision point as the source of the coordinates of the whole area, so as to avoid the doubt that the GPS position of the site may change.

Furthermore, the first GPS data and the second GPS data in the step (C) can be compared by a multi-station-multi-time adjustment method; wherein the first GPS data and the second GPS data are processed first. After that, the equation matrix is adjusted as a whole by the merging method. The invention utilizes the multi-station-multi-time adjustment method to compare and analyze the main factors, and as the GPS observation time increases, and the sampling interval is shortened, the observation data is gradually increased. Therefore, in order to solve such a large amount of observation data and to take into consideration the rigor of the adjustment mode, the present invention facilitates the calculation of data by the multi-station-multi-time adjustment method.

Furthermore, in the step (E), the plurality of electrodes may include a pair of current poles and a pair of potential poles, and the ground resistance measurement system is to continuously flow current or a low frequency alternating direct current into at least one terrain by the pair of current poles. In the subsurface of the variability region, an artificial electric field is established, and the difference in conductivity between the interlayer media is used, and the potential difference is measured through the pair of potential poles, thereby obtaining the apparent resistivity of the formation, thereby estimating the underground formation. Conductivity distribution.

Moreover, in the step (E), a first ground resistance measurement and a second ground resistance measurement may be respectively performed, and the first plurality of electrodes in the first ground resistance measurement are arranged in a Wenner Array. The second electrode is measured in a Schlumberger Array. Wherein, the Wenner Array method is that the pair of current poles are aligned with the pair of potential poles, and the center point thereof is symmetrically arranged, wherein the distance of the pair of current poles is the distance of the pair of potential poles. Three times. However, the advantage of the Winner alignment method is that the measured value is relatively stable, and the resolution of the vertical change of the measured data is high. In addition, the Schlumberger Array is a line in which the pair of current poles are aligned with the pair of potential poles, and the center point thereof is symmetrically arranged, and the pair of current poles and the pair of potential poles are gradually positioned. Expand outward. The main advantage of the Rambigi alignment method is that the accuracy of the instrument is not high, and it is easy to calculate the apparent resistivity.

In conclusion, the present invention can simultaneously perform cross-measurement and analysis comparison using a Wenner Array electrode alignment method and a Schlumberger Array electrode alignment method. Among them, the Wenner Array electrode alignment method is most suitable for slope geological monitoring and can obtain the best vertical direction analysis, and then the Schlumberger Array electrode arrangement method is better. The deep data, the two methods combined and reversed to obtain the electrical profile with high shallow and deep resolution, and then monitor the underground geological conditions of the slope. However, the invention is not limited thereto, and may include a Pole-Pole Array, a Dipole-Dipole Array, or other equivalent alignment methods. Preferably, in step (A), the 3D laser scanner can further capture the slope and generate a first image; and in step (B), the 3D laser scanner further captures the slope and produces a second image. And comparing the at least one terrain variation region and the at least one elevation height variation point in the step (D), and comparing the correspondence between the first image and the second image. Accordingly, the present invention not only performs laser scanning by a 3D laser scanner, but also performs further comparison through shooting of a live photo, thereby avoiding the occurrence of misjudgment.

In addition, the authoring step (C) can perform the nesting analysis by a main controller, and the step (E) is performed by the grounding resistance measurement by a geoelectric surveying instrument, and is also analyzed by the main controller. However, the main controller of the present invention may be a computer workstation, a professional analytical computing host, or other equivalent device.

Please refer to FIG. 1 and FIG. 2 simultaneously. FIG. 1 is a schematic diagram of a measurement system for a slope displacement detecting method according to a preferred embodiment of the present invention. 2 is a main flow chart of a preferred embodiment of the slope displacement detecting method of the present invention. A slope displacement detecting method of the embodiment includes the following steps: First, step (A) first sets three conjugate scanning balls b1, b2, b3, and four GPS receivers G1, G2, G3, G4. At a different position on the side slope SP, and scanning the slope SP with a 3D laser scanner 2, the scanning range covers the plurality of conjugate scanning balls b1, b2, b3. In this embodiment, only four GPS receivers G1, G2, G3, and G4 are disposed, but if the measurement is more accurate, more GPS receivers can be disposed.

However, the scanning principle of the 3D laser scanner used in this embodiment is a time-of-flight method, that is, the scanning head is obtained by calculating the round-trip flight time of the laser hitting the point to be measured and reflected back to the sensor. The distance to the point to be measured is measured, and the coordinate position of the point to be measured is calculated accordingly. Accordingly, the 3D laser scanner 2 generates a first 3D numerical terrain data; and the four GPS receivers G1, G2, G3, and G4 calculate the coordinates and height of the location, and store them as a first A GPS data. In addition, the 3D laser scanner 2 uses the built-in CCD to capture the slope SP and produces a first image.

Moreover, the first GPS data further includes a first reference coordinate of one of the satellite receiving stations GS. The satellite receiving station used in the embodiment is a satellite receiving station North Port Station (PKGN), and the satellite receiving station is used as a reference coordinate. The purpose of this is to induce stable and high-precision points (satellite receiving station GS) as the source of the coordinates of the whole region, so as to avoid the doubt that the GPS position of the site may change.

Next, after step (B) is separated by a specific time, in the present embodiment, the specific time is measured in units of one month. The slope SP is again scanned by the 3D laser scanner 2 and also covers three conjugate scan balls b1, b2, b3, and the 3D laser scanner 2 generates a second 3D numerical topographical data. In addition, the four GPS receivers G1, G2, G3, and G4 calculate their coordinates and heights again, and store them as a second GPS data. Among them, three conjugate scanning balls b1, b2, b3, and four GPS receivers G1, G2, G3, and G4 are disposed at the same position. Similarly, the second GPS data also includes a coordinate measured by the Ministry of the Interior satellite receiving station - North Port Station (PKGN) as a second reference coordinate, and the 3D laser scanner 2 also shoots the slope SP and generates one Second image.

Please refer to FIG. 3 together. FIG. 3 is a flow chart of step (C) of a preferred embodiment of the slope displacement detecting method of the present invention. Furthermore, step (C) nests the first 3D numerical terrain data and the second 3D numerical topographic data to obtain a terrain variation area CA1. The first 3D numerical terrain data includes a first coordinate system and a first point cloud digital model, and the second 3D numerical terrain data includes a second coordinate system and a second point cloud digital model. However, the nested analysis of the first 3D numerical topographic data and the second 3D numerical topographical data in step (C) comprises the following steps: (C1) comparing the first 3D numerical topographic data with the third 3D numerical topographical data in three conjugates The positions of the relative positions of the scanning balls b1, b2, and b3 are changed, and the point where the relative relationship position is not changed is taken as a target conversion conjugate point.

Next, the step (C2) converts the second coordinate system of the second 3D numerical terrain data into the first coordinate system of the first 3D numerical terrain data by using the coordinate conversion conjugate point as a common point of the coordinate system for conversion. In this embodiment, the first coordinate system in the first 3D numerical terrain is first selected as the coordinate system of the reference station, and the second coordinate system in the second 3D numerical terrain is converted to the first coordinate system, and the conjugate sweep The aiming ball is used as a common point in the conversion of the two-coordinate system, and is used as a conversion parameter. Among them, the coordinate conversion is described in detail below.

Please refer to FIG. 4 at the same time. FIG. 4 is a schematic diagram of a spatial coordinate system between the 3D laser scanner 2 and the conjugate scanning balls b1, b2, and b3 of the present invention. Where S is the position of any 3D laser scanner 2, point P is the position of the conjugate scan balls b1, b2, b3, and point O is the selected 3D laser scanner reference station position. ρ is the slant distance between S and P, α is the vertical angle between S and P, and θ is the horizontal angle. The mathematical formula for converting the second coordinate system to the first coordinate system can be written as the following equation (1), which uses three or more conjugate scan balls to convert the known coordinates (Lichti, 2002).

=[ x _p y _p z _p ] ^T : the coordinate vector of the conjugate scan sphere in the first coordinate system;

=[ X _p Y _p Z _p ] ^T : coordinate vector of the conjugate scan sphere in the second coordinate system;

=[ X _s Y _s Z _s ] ^T : The second coordinate system is converted to the coordinate vector of the first coordinate system;

M is a rotation matrix that rotates (ω, φ, κ) angles around the X, Y, and Z axes, respectively;

In addition, the first point cloud digital model in the first 3D numerical terrain data and the second point cloud digit model in the second 3D numerical terrain data are respectively converted into a first terrain numerical model by using the irregular triangulation rule (TIN) ( DEM), and a second topographic numerical model (DEM), as shown in Figures 5A and 5B, respectively. 5A is a point cloud model diagram, and FIG. 5B is a topographic numerical model diagram.

Then, step (C3) and second terrain numerical model (DEM) are merged into the first coordinate system, and all triangular mesh vertex coordinates and vertex coordinates of the first topographic numerical model are calculated and the vertical axis direction and the second The distance ΔZ _{n of the} coordinates of the intersection of the topographic numerical model is also the single point variation. Further, the step (C4) analyzes the equidistance ΔZ _n , wherein the region constituted by the average distance amount is the terrain variation region CA1.

Further details of the operation are as follows: In the numerical terrain model formed by the irregular triangulation, the plane equation of any triangular surface is assumed to be AX + BY + CZ + D + = 0. However, the three vertices constituting the triangular face are I( X ₁ , Y ₁ , Z ₁ ), J( X ₂ , Y ₂ , Z ₂ ), K( X ₃ , Y ₃ , Z ₃ ), respectively.

D =- AX ₁ - BY ₁ - CZ ₁ .

Any Z in the triangle face =- X - Y - .....................(2);

Another assumption is that the space point P ( a ₁ , b ₁ , c ₁ ) and the point Q ( a ₂ , b ₂ , c ₂ ) form a space straight line equation.

If the spatial line L is parallel to the Z axis, then a ₁ = a ₂ and b ₁ = b ₂

The distance from any point D ( X _d , Y _d , Z _d ) on the spatial straight line L to the intersection point E ( X _e , Y _e , Z _e ) between the spatial straight line and the above triangular surface is

S =| Z _e - Z _d |, Z _e =- X _d - Y _d - ........................(3)

S in equation (3) is the amount of difference between the two scan models, as shown in FIG. among them:

Z _e - Z _d >0, indicating that the surface changes convex, Z _e - Z _d =0, indicating no change in the surface, Z _e - Z _d <0, indicating that the surface changes concave.

In addition, in the embodiment, the step (C3) further comprises interpolating the equidistant distance Δ Z _n into a variation quantity numerical model by using the irregular triangulation method, as a representation of the two-model post-stacking difference quantity analysis and representation.

In addition, in step (C), the first GPS data and the second GPS data are compared by a multi-station-multi-time adjustment method; wherein, the first GPS data and the second GPS data are first analyzed. After processing, the equation matrix is adjusted as a whole by the merging method.

Accordingly, the present embodiment utilizes a multi-station-multi-time adjustment method to analyze the main factors, and as the GPS observation time increases, the GPS observation layout increases, and the sampling interval is shortened, the observation data is gradually increased. . Therefore, in order to solve such a large amount of observation data, and to take into account the rigor of the adjustment mode, it is convenient to use the multi-station-multi-time adjustment method for data calculation.

Next, after the step (C), the step (D) is performed, that is, whether the terrain variation region CA1 and the coordinate height variation point G4 are matched. That is to say, in the present embodiment, as shown in FIG. 1, the coordinate height variation point G4 falls within the terrain variation region CA1, and thus is completely in conformity. Certainly, the present embodiment simultaneously refers to the correspondence between the first image and the second image to perform further comparison analysis through the photograph taken on the spot, and the occurrence of the false positive situation can be avoided.

Finally, in step (E), when the comparison is met, four electrodes E1, E2, E3, and E4 are disposed in the topographical variation region CA1 and the ground resistance measurement is performed by a geoelectric surveying device 4 to obtain a terrain variation region CA1. The apparent resistivity of the formation (Apparent Resistivity). In this embodiment, the four electrodes E1, E2, E3, and E4 include a pair of current electrodes E1, E4, and a pair of potential electrodes E2, E3. Wherein, the ground resistance measurement is performed by the pair of current poles E1, E4 of the power module 41 to continuously flow or a low frequency alternating direct current to the ground of the terrain variation area CA1 to establish an artificial electric field. At the same time, there is a difference in conductivity between the inter-layer media. The geoelectric surveying instrument 4 measures the potential difference for the pair of potential poles E2, E3, thereby obtaining the apparent resistivity of the formation, and thereby estimating the conductivity distribution of the subterranean formation.

It is particularly worth mentioning that in the conventional resistance monitoring, only one type of electrode arrangement is used, and the obtained data may not be complete enough, which often leads to erroneous judgment. However, the present embodiment uses different electrode arrangement methods to completely reflect the underground structure distribution, comprehensively combine the advantages of various electrode alignment methods, and increase the integrity of the observed values, and the multi-variability of the slope formation monitoring of the present invention is also domestically The first case. In this embodiment, step (E) performs a first ground resistance measurement and a second ground resistance measurement, wherein the four electrodes E1, E2, E3, and E4 in the first ground resistance measurement are in Wennai. The arrangement of the Wenner Array is arranged. In the second ground resistance measurement, the plurality of electrodes E1, E2, E3, and E4 are arranged in a Schlumberger Array.

Please refer to FIG. 7 and FIG. 8 simultaneously. FIG. 7 is a schematic diagram of a Winner arrangement electrode arrangement method according to a preferred embodiment of the present invention, and FIG. 8 is a Schlumberg arrangement electrode arrangement according to a preferred embodiment of the present invention. Schematic diagram of the law. As shown in the figure, the Wenner Array (Wenner Array) is arranged in a line with the pair of potential poles E1, E4 and the pair of potential poles E2, E3, and the center point is symmetrically arranged, wherein the pair of currents The distance between the poles E1 and E4 is three times the distance between the pair of potential poles E2 and E3, and E1E2 = E2E3 = E3E4 = a. Therefore, the depth that can be detected is about the distance of E2E3. When the E1E4 and E2E3 spacing is increased, the shallow to deep formation information can be obtained successively. However, the advantage of the Winner alignment method is that the measured value is relatively stable, and the resolution of the vertical change of the measured data is high. Most of the resistivity of the measured outcrop adopts the electrode arrangement method, but the disadvantage is that it is restricted by the terrain and is measured. It is quite time consuming, and the deeper data contains a higher noise ratio.

As for the Schlumberger Array, the pair of current poles E1, E4 are aligned with the pair of potential poles E2, E3, and the center point is symmetrically arranged, and the pair of current poles E1, E4 The position of the pair of potential electrodes E2, E3 gradually expands outward. However, when the half spread is gradually increased, the resistivity change of the formation from shallow to deep can be obtained. The main advantage of the Schrambji arrangement is that the accuracy of the instrument is not high, and the calculation of the apparent resistivity is easy, but its disadvantage is that the detection is more labor-intensive, and each time the potential is moved, the local heterogeneity of the shallow surface will be underground. Mixed signal causes poor quality of data and leads to misjudgment.

In this embodiment, all of the above measurements are performed by a main controller 3 after all the measurement data is acquired, and then the nested analysis or the comparison analysis is performed. However, the main controller 3 of this embodiment may be a general desktop computer, a notebook computer, a computer workstation, a professional analysis computing host, or other equivalent device.

The above-mentioned embodiments are merely examples for convenience of description, and the scope of the claims is intended to be limited to the above embodiments.

2. . . 3D laser scanner

3. . . main controller

4. . . Geoelectric surveyor

41. . . Power module

B1, b2, b3. . . Conjugate scan ball

SP. . . Slope

G1, G2, G3, G4. . . GPS receiver

E1, E2, E3, E4. . . electrode

GS. . . Satellite receiving station

CA1. . . Terrain variation region

G4. . . Coordinate height variation point

P. . . Conjugate scan ball position

S. . . The position of any 3D laser scanner

ρ. . . Slope between S and P

α. . . Vertical angle between S and P

θ. . . Horizontal angle

O. . . Selected 3D laser scanner reference station position

1 is a schematic diagram of a measurement system in accordance with a preferred embodiment of the present invention.

2 is a main flow chart of a preferred embodiment of the present invention.

3 is a flow chart of step (C) of a preferred embodiment of the present invention.

4 is a schematic diagram of a spatial coordinate system between a 3D laser scanner and a conjugate scanning sphere of the present invention.

Figure 5A is a diagram of a point cloud model.

Figure 5B is a topographical numerical model diagram.

Fig. 6 is a schematic diagram showing the analysis of the difference in the amount of the secondary scanning surface.

Fig. 7 is a schematic view showing the arrangement of electrodes of the Winner arrangement method according to a preferred embodiment of the present invention.

Fig. 8 is a schematic view showing the arrangement of electrodes of the Schramb arrangement method according to a preferred embodiment of the present invention.

Claims (8)

A slope displacement detecting method comprises the following steps: (A) respectively setting a complex conjugate scanning sphere and a plurality of GPS receivers at different positions on one side of the slope, and scanning the slope with a 3D laser scanner And covering the complex conjugate scan ball, the 3D laser scanner generates a first 3D numerical terrain data; the complex GPS receiver calculates its coordinates and height respectively, and stores it as a first GPS data; B) after a certain time interval, the 3D laser scanner scans the slope again and covers the complex conjugate scan ball, and the 3D laser scanner generates a second 3D numerical terrain data; the complex GPS The receiver again calculates its coordinates and height, and stores it as a second GPS data; (C) nesting and analyzing the first 3D numerical terrain data and the second 3D numerical topographic data to obtain at least one terrain variation region; Aligning the first GPS data with the second GPS data to obtain at least one standard height variation point; (D) comparing the at least one terrain variation region with the at least one standard height variation point; and (E) when the at least a terrain variation region and the at least one landmark height variation When the points are coincident, the plurality of electrodes are disposed to perform resistance measurement on the at least one terrain variation region to obtain a formation apparent resistivity of the at least one terrain variation region; wherein the first 3D numerical topographic data includes a first coordinate system, a first point cloud digital model, the second 3D numerical terrain data The method includes a second coordinate system and a second point cloud digital model. The step (C) of nesting and analyzing the first 3D numerical topographic data and the second 3D numerical topographical data includes the following steps: (C1) comparing The relative position relationship between the first 3D numerical topographic data and the complex conjugate scanning sphere in the second 3D numerical topographical data is changed, and the point where the relative relationship position is not changed is taken as a target conversion conjugate point; (C2) Converting the second coordinate system of the second 3D numerical terrain data into the first coordinate system of the first 3D numerical terrain data by using the coordinate conversion conjugate point as a common point when converting the coordinate system; The first point cloud digital model of the 3D numerical terrain data and the second point cloud digital model of the second 3D numerical terrain data are respectively converted into a first terrain numerical model by using an irregular triangular network rule (TIN), and a first (2) a topographic numerical model; (C3) the second topographical numerical model is merged into the first coordinate system, and all triangular mesh vertex coordinates in the first topographic numerical model and the vertical axis direction of the vertex coordinate extension are calculated a distance between the coordinates of the intersection of the second topographical numerical model; and (C4) analyzing the equidistance, wherein the region formed by the greater than the average distance is the at least one terrain variation region; wherein, in step (E), a first ground resistance is respectively performed Measuring, and a second ground resistance measurement, wherein the plurality of electrodes are arranged in a Wenner Array in the first ground resistance measurement, wherein the plurality of electrodes are in the second ground resistance measurement The Schlumberger Array is arranged. The method for detecting a slope displacement according to claim 1, wherein the step (C3) further comprises interpolating the equidistant distance into a variable quantity numerical model by using an irregular triangulation method. The method for detecting a slope displacement according to claim 1, wherein the first GPS data in the step (A) further comprises a first reference coordinate of one of the satellite receiving stations, and the step (B) The second GPS data further includes a second reference coordinate of the satellite receiving station; wherein, in the step (C), the first GPS data and the second GPS data are compared by the first reference coordinate and the second The base coordinates are the comparison benchmarks. The method for detecting a slope displacement according to claim 1, wherein the first GPS data and the second GPS data in the step (C) are compared by a multi-station-multi-time adjustment method. The analysis comprises: first performing data processing of the first GPS data and the second GPS data, and then adjusting the equation matrix as a whole by a combination method. The method for detecting a slope displacement according to claim 1, wherein the plurality of electrodes in the step (E) includes a pair of current poles and a pair of potential poles, wherein the ground resistance is measured by the pair of current poles A continuous current or a low frequency alternating direct current is passed into the ground of the at least one terrain variation region, and the potential difference is measured through the pair of potential poles. The method for detecting a slope displacement according to the first aspect of the invention, wherein the pair of current poles is arranged in a line with the pair of potential poles, and the center point thereof is symmetrically arranged, wherein The distance between the pair of current poles is three times the distance of the pair of potential poles; the Schramb arrangement method is to align the pair of current poles with the pair of potential poles, and take the center point thereof symmetrically. Current pole and the pair of potential poles The position gradually expands outward. The method for detecting a slope displacement according to claim 1, wherein the 3D laser scanner in the step (A) further captures the slope and generates a first image; in the step (B) The 3D laser scanner further captures the slope and generates a second image; and in the step (D), comparing the at least one terrain variation region with at least one standard height variation point, and comparing the first The image, and the correspondence of the second image. The method for detecting a slope displacement according to claim 1, wherein the step (C) is performed by a main controller, and the step (E) is performed by a geoelectric surveying device. Resistance measurement.