Geological disaster hidden danger identification method and system based on insar technology

Abstract

The application discloses a geological disaster hidden danger identification method and system based on InSAR technology, and is applied to the technical field of intelligent monitoring, and the method comprises the following steps: acquiring SBAS deformation amounts of multiple monitoring points in a target region, and dividing the target region into multiple monitoring grids according to the topography of the target region; acquiring a monitoring point closest to a geometric center point of a monitoring grid as a reference monitoring point, and taking the SBAS deformation amount of the reference monitoring point as a reference deformation amount of the monitoring grid; and judging a geological disaster hidden danger risk of the monitoring grid according to the difference between the SBAS deformation amounts of all monitoring points in the monitoring grid and the reference deformation amount. According to the method, the error generated in the InSAR technology is synchronously removed through relative data analysis, the actual deformation condition in a certain region can be effectively analyzed to find out the geological disaster hidden danger, the safety of the geological disaster risk body is improved, and the risk body can be extracted and early warned in time.

Description

Technical Field

[0001] This invention relates to intelligent monitoring technology, specifically to a method and system for identifying geological hazard risks based on InSAR technology. Background Technology

[0002] SBAS (Small Baseline Subset) data is a surface deformation product obtained by performing InSAR interferometry, phase unwrapping, and time-series inversion on multiple SAR images of the same area. It belongs to the post-processing output of the InSAR technology chain. Although InSAR technology can capture surface deformation at the millimeter level, atmospheric delay and DEM errors can easily introduce centimeter-level spurious deformations in the identification of geological hazards. Furthermore, the long revisit period and the failure of phase unwrapping in large gradient deformation areas often lead to missed detections.

[0003] In the prior art, Chinese patent application number CN202510168657.7 discloses a time-series InSAR method suitable for multi-dimensional deformation monitoring of landslides. First, it acquires the deformation observations of the landslide along the line of sight. Then, it inputs these observations into a monitoring model containing a fractal composite structural deformation constraint model, thereby obtaining the deformation time series of the landslide body in the vertical, east-west, downslope, and perpendicular-downslope directions. This method has significant advantages, as it can simultaneously acquire the deformation time series characteristics of the landslide along four dimensions: vertical, east-west, downslope, and perpendicular-downslope. The constructed model covers linear, acceleration, periodic, high-frequency, and thermal expansion deformations, fully considering the complexity of landslide deformation and exhibiting strong adaptability.

[0004] As can be seen from the existing technology, the existing technology focuses on the time series effect of the data generated by InSAR. The time series effect is affected by the atmospheric delay and DEM error mentioned above, which will produce huge errors and affect the final analysis results. Summary of the Invention

[0005] In order to at least overcome the above-mentioned shortcomings in the prior art, the purpose of this application is to provide a method and system for identifying geological hazard risks based on InSAR technology.

[0006] In a first aspect, embodiments of this application provide a method for identifying geological hazard risks based on InSAR technology, including:

[0007] Obtain the SBAS deformation of multiple monitoring points in the target area, and divide the target area into multiple monitoring grids according to the terrain of the target area;

[0008] The monitoring point closest to the geometric center of the monitoring grid is obtained as the reference monitoring point, and the SBAS deformation of the reference monitoring point is used as the reference deformation of the monitoring grid.

[0009] The geological hazard risk of the monitoring grid is determined based on the difference between the SBAS deformation of all monitoring points within the monitoring grid and the baseline deformation.

[0010] In one possible implementation, dividing the target area into multiple monitoring grids based on the terrain of the target area includes:

[0011] Obtain the direction of movement of the target area during the disaster, and construct multiple vertical grid lines along the direction of movement;

[0012] Multiple horizontal grid lines are constructed in directions orthogonal to the vertical grid lines;

[0013] The area enclosed by the vertical grid lines, the horizontal grid lines, and the boundary line of the target area is defined as the monitoring grid.

[0014] In one possible implementation, assessing the geological hazard risk of the monitoring grid includes:

[0015] Calculate the physical distance from all monitoring points within the monitoring grid to the reference monitoring point, and calculate the absolute value of the difference between each SBAS deformation and the reference deformation as the deformation difference value;

[0016] The ratio of the deformation difference value to the corresponding physical distance is calculated as the risk difference value of the corresponding monitoring point;

[0017] The average risk difference value of all monitoring points within the monitoring grid is calculated as the risk warning value of the monitoring grid.

[0018] When the risk warning value exceeds the preset warning value, it is determined that there is a risk of geological disaster in the monitoring grid.

[0019] In one possible implementation, the risk warning value is calculated according to the following formula:

[0020]

[0021] In the formula, R is the risk warning value of the monitoring grid, R i Let x be the risk difference value for the i-th monitoring point. i Let y be the horizontal coordinate value of the i-th monitoring point. i Let D be the vertical coordinate value of the i-th monitoring point. i Let x0 be the SBAS deformation of the i-th monitoring point, y0 be the horizontal coordinate of the reference monitoring point, y0 be the vertical coordinate of the reference monitoring point, D0 be the reference deformation of the reference monitoring point, and n be the number of monitoring points in the monitoring grid.

[0022] In one possible implementation, the acquisition of the benchmark monitoring points includes:

[0023] Obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point;

[0024] The monitoring point with the smallest physical distance is taken as the baseline monitoring point.

[0025] Secondly, this application also provides a geological hazard identification system based on InSAR technology, including:

[0026] The acquisition unit is configured to acquire the SBAS deformation of multiple monitoring points in the target area, and divide the target area into multiple monitoring grids according to the terrain of the target area;

[0027] The reference unit is configured to obtain the monitoring point closest to the geometric center point of the monitoring grid as the reference monitoring point, and use the SBAS deformation of the reference monitoring point as the reference deformation of the monitoring grid.

[0028] The judgment unit is configured to determine the geological hazard risk of the monitoring grid based on the difference between the SBAS deformation of all monitoring points within the monitoring grid and the reference deformation.

[0029] In one possible implementation, the acquisition unit is further configured as follows:

[0030] Obtain the direction of movement of the target area during the disaster, and construct multiple vertical grid lines along the direction of movement;

[0031] Multiple horizontal grid lines are constructed in directions orthogonal to the vertical grid lines;

[0032] The area enclosed by the vertical grid lines, the horizontal grid lines, and the boundary line of the target area is defined as the monitoring grid.

[0033] In one possible implementation, the determination unit is further configured as follows:

[0034] Calculate the physical distance from all monitoring points within the monitoring grid to the reference monitoring point, and calculate the absolute value of the difference between each SBAS deformation and the reference deformation as the deformation difference value;

[0035] The ratio of the deformation difference value to the corresponding physical distance is calculated as the risk difference value of the corresponding monitoring point;

[0036] The average risk difference value of all monitoring points within the monitoring grid is calculated as the risk warning value of the monitoring grid.

[0037] When the risk warning value exceeds the preset warning value, it is determined that there is a risk of geological disaster in the monitoring grid.

[0038] In one possible implementation, the determination unit is further configured as follows:

[0039] The risk warning value is calculated according to the following formula:

[0040]

[0041] In the formula, R is the risk warning value of the monitoring grid, R i Let x be the risk difference value for the i-th monitoring point. i Let y be the horizontal coordinate value of the i-th monitoring point. i Let D be the vertical coordinate value of the i-th monitoring point. i Let x0 be the SBAS deformation of the i-th monitoring point, y0 be the horizontal coordinate of the reference monitoring point, y0 be the vertical coordinate of the reference monitoring point, D0 be the reference deformation of the reference monitoring point, and n be the number of monitoring points in the monitoring grid.

[0042] In one possible implementation, the reference unit is further configured as follows:

[0043] Obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point;

[0044] The monitoring point with the smallest physical distance is taken as the baseline monitoring point.

[0045] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0046] This invention relates to a geological hazard identification method and system based on InSAR technology. By using relative data analysis, the errors generated by InSAR technology are eliminated synchronously. This method can effectively analyze the actual deformation in a certain area to identify geological hazard risks, improve the safety of geological hazard risk bodies, and enable timely extraction and early warning of risk bodies. Attached Figure Description

[0047] The accompanying drawings, which are included to provide a further understanding of embodiments of the invention and form part of this application, do not constitute a limitation thereof. In the drawings:

[0048] Figure 1 This is a schematic diagram of the method steps in an embodiment of this application;

[0049] Figure 2 This is a schematic diagram of an embodiment of this application. Detailed Implementation

[0050] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. It should be understood that the accompanying drawings in this application are for illustrative and descriptive purposes only and are not intended to limit the scope of protection of this application. Furthermore, it should be understood that the schematic drawings are not drawn to scale. The flowcharts used in this application illustrate operations implemented according to some embodiments of this application. It should be understood that the operations in the flowcharts may not be implemented in sequence, and steps without logical contextual relationships may be reversed or implemented simultaneously. In addition, those skilled in the art, guided by the content of this application, may add one or more other operations to the flowcharts, or remove one or more operations from the flowcharts.

[0051] Furthermore, the described embodiments are merely some, not all, of the embodiments of this application. The components of the embodiments of this application described and illustrated herein can typically be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely to illustrate selected embodiments of the application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application.

[0052] Please refer to the figure for a flowchart illustrating the geological hazard identification method based on InSAR technology provided in this embodiment of the invention. Further, the geological hazard identification method based on InSAR technology may specifically include the contents described in steps S1-S3.

[0053] S1: Obtain the SBAS deformation of multiple monitoring points in the target area, and divide the target area into multiple monitoring grids according to the terrain of the target area;

[0054] S2: Obtain the monitoring point closest to the geometric center of the monitoring grid as the reference monitoring point, and use the SBAS deformation of the reference monitoring point as the reference deformation of the monitoring grid;

[0055] S3: Determine the geological hazard risk of the monitoring grid based on the difference between the SBAS deformation of all monitoring points within the monitoring grid and the baseline deformation.

[0056] In implementing this embodiment, it is necessary to first obtain the SBAS deformation of multiple monitoring points in the target area using InSAR technology. This can be achieved through techniques such as joint decomposition of ascending and descending orbit observations, which are mature existing technologies. This embodiment does not impose many limitations and can refer to the standard T / CAGHP 013—2018 Guidelines for InSAR Monitoring Technology of Geological Hazards for monitoring. Furthermore, the target area described in this embodiment is generally a complete geological hazard risk body, such as a complete landslide. To reduce errors introduced by InSAR technology, this embodiment divides the target area into multiple monitoring grids and evaluates each grid individually. For InSAR technology, atmospheric delay and DEM error generally remain consistent within a monitoring grid during satellite imaging. Therefore, independently evaluating each monitoring grid can effectively reduce evaluation errors.

[0057] In this embodiment, it is necessary to first obtain benchmark monitoring points. Each benchmark monitoring point is located at the position closest to the geometric center of the monitoring grid in terms of physical distance. Based on the benchmark deformation generated by the benchmark monitoring points and the deformation of other monitoring points within the grid, the overall deformation within the monitoring grid can be analyzed, thereby determining the geological hazard risk of the monitoring grid. This invention uses relative data analysis to simultaneously eliminate errors generated by InSAR technology, which can effectively analyze the actual deformation in a certain area to identify geological hazard risks, improve the safety of geological hazard risk bodies, and enable timely extraction and early warning of risk bodies.

[0058] In one possible implementation, dividing the target area into multiple monitoring grids based on the terrain of the target area includes:

[0059] Obtain the direction of movement of the target area during the disaster, and construct multiple vertical grid lines along the direction of movement;

[0060] Multiple horizontal grid lines are constructed in directions orthogonal to the vertical grid lines;

[0061] The area enclosed by the vertical grid lines, the horizontal grid lines, and the boundary line of the target area is defined as the monitoring grid.

[0062] When implementing the embodiments of this application, please refer to Figure 2 This paper provides an example of a specific method for dividing a monitoring grid, which requires first obtaining the direction of movement of the target area during a disaster. For landslides, this is generally the direction of the sliding surface. Figure 2In the image, the landslide location is marked with a green line, and the entire landslide area constitutes the target region. The possible sliding direction of the landslide is then determined. In this example, the possible sliding direction is from the top of the slope to the toe of the slope near the river, so a horizontal grid is constructed based on this. Simultaneously, a vertical grid orthogonal to the horizontal grid lines is constructed, dividing the area into multiple monitoring grids. From... Figure 2 The data shows monitoring points within different monitoring grids, represented by points ranging from blue to red. The closer the color is to red, the greater the detected deformation. It should be understood that the detected deformation is a cumulative deformation over a time period, not a real-time figure, because the satellite needs to fly over the area to generate data. It should also be understood that the specific spacing of the horizontal and vertical grid lines is generally based on experience, such as a width range of 50-100m, to ensure that atmospheric influences within this range are as consistent as possible. Figure 2 The grid width shown is 70m.

[0063] In one possible implementation, assessing the geological hazard risk of the monitoring grid includes:

[0064] Calculate the physical distance from all monitoring points within the monitoring grid to the reference monitoring point, and calculate the absolute value of the difference between each SBAS deformation and the reference deformation as the deformation difference value;

[0065] The ratio of the deformation difference value to the corresponding physical distance is calculated as the risk difference value of the corresponding monitoring point;

[0066] The average risk difference value of all monitoring points within the monitoring grid is calculated as the risk warning value of the monitoring grid.

[0067] When the risk warning value exceeds the preset warning value, it is determined that there is a risk of geological disaster in the monitoring grid.

[0068] In the implementation of this application embodiment, a relative deformation calculation method is adopted in judging the geological disaster risk of the monitoring grid. For the same monitoring grid, the error value of SBAS deformation is relatively similar. Therefore, it is necessary to first calculate the physical distance from different monitoring points to the reference monitoring point and the absolute difference between the SBAS deformation of each monitoring point and the reference deformation. This absolute difference can effectively eliminate the error within the monitoring grid. Therefore, the absolute difference can better express the deformation difference between different monitoring points and near the center point of the monitoring grid.

[0069] In this embodiment, since the deformation of a monitoring grid is continuous, calculating the ratio of the aforementioned deformation difference to the specific physical distance can express the deformation of the landslide per unit distance, which in turn reflects the stability of the area. Simultaneously, this ratio can effectively reflect factors affecting the deformation differences at different points within the monitoring grid, such as slope steepness and vegetation cover. It is a comprehensive indicator that effectively expresses the stability within the monitoring grid. It should be understood that the satellite image generated by InSAR technology is a two-dimensional graphic. Therefore, the physical distance described in this embodiment is a physical distance within a two-dimensional graphic, calculated using horizontal and vertical coordinates.

[0070] In the implementation of this application embodiment, the risk warning value within the monitoring grid can be obtained by calculating the average of the risk difference values ​​of all monitoring points within the monitoring grid. Then, by combining this risk warning value with the statistically obtained alert value, the risk situation within the monitoring grid can be determined. The alert value needs to be selected based on the period of the aforementioned SBAS deformation, with a preferred periodic alert value of 8.3 mm / km in annual units.

[0071] In one possible implementation, the risk warning value is calculated according to the following formula:

[0072]

[0073] In the formula, R is the risk warning value of the monitoring grid, R i Let x be the risk difference value for the i-th monitoring point. i Let y be the horizontal coordinate value of the i-th monitoring point. i Let D be the vertical coordinate value of the i-th monitoring point. i Let x0 be the SBAS deformation of the i-th monitoring point, y0 be the horizontal coordinate of the reference monitoring point, y0 be the vertical coordinate of the reference monitoring point, D0 be the reference deformation of the reference monitoring point, and n be the number of monitoring points in the monitoring grid.

[0074] In the implementation of this application embodiment, a specific risk warning value calculation scheme is provided, wherein the horizontal coordinate value is generally selected from the coordinate value along the latitude direction, and the vertical coordinate value is generally selected from the coordinate value along the longitude direction.

[0075] In one possible implementation, the acquisition of the benchmark monitoring points includes:

[0076] Obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point;

[0077] The monitoring point with the smallest physical distance is taken as the baseline monitoring point.

[0078] Based on the same inventive concept, this application also provides a geological hazard hazard identification system based on InSAR technology, including:

[0079] The acquisition unit is configured to acquire the SBAS deformation of multiple monitoring points in the target area, and divide the target area into multiple monitoring grids according to the terrain of the target area;

[0080] The reference unit is configured to obtain the monitoring point closest to the geometric center point of the monitoring grid as the reference monitoring point, and use the SBAS deformation of the reference monitoring point as the reference deformation of the monitoring grid.

[0081] The judgment unit is configured to determine the geological hazard risk of the monitoring grid based on the difference between the SBAS deformation of all monitoring points within the monitoring grid and the reference deformation.

[0082] In one possible implementation, the acquisition unit is further configured as follows:

[0083] Obtain the direction of movement of the target area during the disaster, and construct multiple vertical grid lines along the direction of movement;

[0084] Multiple horizontal grid lines are constructed in directions orthogonal to the vertical grid lines;

[0085] The area enclosed by the vertical grid lines, the horizontal grid lines, and the boundary line of the target area is defined as the monitoring grid.

[0086] In one possible implementation, the determination unit is further configured as follows:

[0087] Calculate the physical distance from all monitoring points within the monitoring grid to the reference monitoring point, and calculate the absolute value of the difference between each SBAS deformation and the reference deformation as the deformation difference value;

[0088] The ratio of the deformation difference value to the corresponding physical distance is calculated as the risk difference value of the corresponding monitoring point;

[0089] The average risk difference value of all monitoring points within the monitoring grid is calculated as the risk warning value of the monitoring grid.

[0090] When the risk warning value exceeds the preset warning value, it is determined that there is a risk of geological disaster in the monitoring grid.

[0091] In one possible implementation, the determination unit is further configured as follows:

[0092] The risk warning value is calculated according to the following formula:

[0093]

[0094] In the formula, R is the risk warning value of the monitoring grid, R i Let x be the risk difference value for the i-th monitoring point. i Let y be the horizontal coordinate value of the i-th monitoring point. i Let D be the vertical coordinate value of the i-th monitoring point. i Let x0 be the SBAS deformation of the i-th monitoring point, y0 be the horizontal coordinate of the reference monitoring point, y0 be the vertical coordinate of the reference monitoring point, D0 be the reference deformation of the reference monitoring point, and n be the number of monitoring points in the monitoring grid.

[0095] In one possible implementation, the reference unit is further configured as follows:

[0096] Obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point;

[0097] The monitoring point with the smallest physical distance is taken as the baseline monitoring point.

[0098] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0099] In the several embodiments provided in this application, it should be understood that the disclosed apparatus and methods can be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative. For instance, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. In addition, the mutual coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices or units, or may be electrical, mechanical or other forms of connection.

[0100] The units described as separate components may or may not be physically separate. As will be apparent to those skilled in the art, the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

[0101] Furthermore, the functional units in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit.

[0102] If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or grid device, etc.) to execute all or part of the steps of the methods described in the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0103] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for identifying geological hazard risks based on InSAR technology, characterized in that, include: Obtain the SBAS deformation of multiple monitoring points in the target area, and divide the target area into multiple monitoring grids according to the terrain of the target area; The monitoring point closest to the geometric center of the monitoring grid is obtained as the reference monitoring point, and the SBAS deformation of the reference monitoring point is used as the reference deformation of the monitoring grid. The geological hazard risk of the monitoring grid is determined based on the difference between the SBAS deformation of all monitoring points within the monitoring grid and the benchmark deformation. The assessment of geological hazard risks in the monitoring grid includes: Calculate the physical distance from all monitoring points within the monitoring grid to the reference monitoring point, and calculate the absolute value of the difference between each SBAS deformation and the reference deformation as the deformation difference value; The ratio of the deformation difference value to the corresponding physical distance is calculated as the risk difference value of the corresponding monitoring point; The average risk difference value of all monitoring points within the monitoring grid is calculated as the risk warning value of the monitoring grid. When the risk warning value exceeds the preset warning value, it is determined that there is a risk of geological disaster in the monitoring grid.

2. The geological hazard identification method based on InSAR technology according to claim 1, characterized in that, The target area is divided into multiple monitoring grids based on its terrain, including: Obtain the direction of movement of the target area during the disaster, and construct multiple vertical grid lines along the direction of movement; Multiple horizontal grid lines are constructed in directions orthogonal to the vertical grid lines; The area enclosed by the vertical grid lines, the horizontal grid lines, and the boundary line of the target area is defined as the monitoring grid.

3. The geological hazard identification method based on InSAR technology according to claim 1, characterized in that, The risk warning value is calculated according to the following formula: ； In the formula, R is the risk warning value of the monitoring grid, R i Let x be the risk difference value for the i-th monitoring point. i Let y be the horizontal coordinate value of the i-th monitoring point. i Let D be the vertical coordinate value of the i-th monitoring point. i Let x0 be the SBAS deformation of the i-th monitoring point, y0 be the horizontal coordinate of the reference monitoring point, y0 be the vertical coordinate of the reference monitoring point, D0 be the reference deformation of the reference monitoring point, and n be the number of monitoring points in the monitoring grid.

4. The geological hazard identification method based on InSAR technology according to claim 1, characterized in that, The acquisition of the benchmark monitoring points includes: Obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point; The monitoring point with the smallest physical distance is taken as the baseline monitoring point.

5. A geological hazard identification system based on InSAR technology, characterized in that, include: The acquisition unit is configured to acquire the SBAS deformation of multiple monitoring points in the target area, and divide the target area into multiple monitoring grids according to the terrain of the target area; The reference unit is configured to obtain the monitoring point closest to the geometric center point of the monitoring grid as the reference monitoring point, and use the SBAS deformation of the reference monitoring point as the reference deformation of the monitoring grid. The judgment unit is configured to judge the geological hazard risk of the monitoring grid based on the difference between the SBAS deformation of all monitoring points in the monitoring grid and the reference deformation. The judgment unit is further configured to: Calculate the physical distance from all monitoring points within the monitoring grid to the reference monitoring point, and calculate the absolute value of the difference between each SBAS deformation and the reference deformation as the deformation difference value; The ratio of the deformation difference value to the corresponding physical distance is calculated as the risk difference value of the corresponding monitoring point; The average risk difference value of all monitoring points within the monitoring grid is calculated as the risk warning value of the monitoring grid. When the risk warning value exceeds the preset warning value, it is determined that there is a risk of geological disaster in the monitoring grid.

6. The geological hazard identification system based on InSAR technology according to claim 5, characterized in that, The acquisition unit is further configured to: Obtain the direction of movement of the target area during the disaster, and construct multiple vertical grid lines along the direction of movement; Multiple horizontal grid lines are constructed in directions orthogonal to the vertical grid lines; The area enclosed by the vertical grid lines, the horizontal grid lines, and the boundary line of the target area is defined as the monitoring grid.

7. The geological hazard identification system based on InSAR technology according to claim 5, characterized in that, The judgment unit is further configured to: The risk warning value is calculated according to the following formula: ； In the formula, R is the risk warning value of the monitoring grid, R i Let x be the risk difference value for the i-th monitoring point. i Let y be the horizontal coordinate value of the i-th monitoring point. i Let D be the vertical coordinate value of the i-th monitoring point. i Let x0 be the SBAS deformation of the i-th monitoring point, y0 be the horizontal coordinate of the reference monitoring point, y0 be the vertical coordinate of the reference monitoring point, D0 be the reference deformation of the reference monitoring point, and n be the number of monitoring points in the monitoring grid.

8. The geological hazard identification system based on InSAR technology according to claim 5, characterized in that, The reference unit is also configured to: Obtain the coordinates of the geometric center point of the monitoring grid, and calculate the physical distance from all monitoring points in the monitoring grid to the geometric center point; The monitoring point with the smallest physical distance is taken as the baseline monitoring point.

Images