INSAR-based Structural Safety Risk Assessment and Early Warning Safety Monitoring System

Abstract

The present invention discloses a safety risk assessment and early warning safety monitoring system for structures based on INSAR, which relates to the technical field of risk assessment; elevation data of the structure and elevation data of the terrain are obtained through INSAR technology; an initial structure model is established according to the elevation data of the structure, the initial structure model is divided into several regions, and corresponding risk values for each region are generated according to the elevation data of the terrain; a safety value is set, and a stable structure model is generated according to the risk value and the safety value; a monitoring period is set, and a real-time structure model is regularly obtained according to the monitoring period; it is detected whether there are contact points between the real-time structure model and the stable structure model, if so, an early warning information and a risk detection period are generated; a periodic structure model is regularly generated according to the risk detection period; the elevation difference of the periodic structure model is obtained, and an elevation change rate is generated according to the elevation difference; a risk information and a new risk detection period are generated according to the elevation change rate.

Description

Technical Field

[0001] The present invention relates to the technical field, specifically a safety risk assessment and early warning safety monitoring system for structures based on INSAR. Background Art

[0002] In recent years, due to factors such as natural disasters, human disturbances, and aging, the building safety problems of structures have become increasingly prominent globally. Traditional detection methods, such as visual inspections and structural monitoring, often require a large amount of manpower, time, and cost, and have limitations in early identification of potential hazards.

[0003] As a non-contact, large-scale, and real-time monitoring method, remote sensing technology has shown great potential in the field of structure safety assessment in recent years. Among them, the Interferometric Synthetic Aperture Radar (INSAR) technology has become an effective tool for measuring the deformation changes of structures due to its high resolution, insensitivity to weather conditions and lighting, etc. The monitoring system based on INSAR technology can accurately capture the tiny displacements of structures, which helps to effectively identify potential safety hazards.

[0004] However, the method of only judging whether there are safety hazards in a structure by monitoring whether the structure deforms is too single. How to incorporate the terrain around the structure into the influencing factors of the structure's safety hazards, so as to more comprehensively and accurately evaluate whether there are safety risks in the structure, is an issue that we need to consider currently. Summary of the Invention

[0005] In order to solve the above problems, the purpose of the present invention is to provide a safety risk assessment and early warning safety monitoring system for structures based on INSAR.

[0006] The purpose of the present invention can be achieved by the following technical solutions: A safety risk assessment and early warning safety monitoring system for structures based on INSAR, including a control center, and the control center is communicatively connected to an elevation data acquisition module, an elevation data processing and analysis module, a structure risk detection module, a structure risk assessment module, and an early warning module;

[0007] The elevation data acquisition module is used to obtain the elevation data of the structure and the elevation data of the terrain around the structure through INSAR technology;

[0008] The elevation data processing and analysis module is used to establish an initial structure model based on the elevation data of the structure, divide the initial structure model into several regions, generate corresponding risk values for each region according to the elevation data of the terrain; set a safety value, and generate a stable structure model according to the risk values and the safety value of each region;

[0009] The structure risk detection module is used to set a monitoring period and regularly obtain a real-time structure model according to the monitoring period; detect whether there are contact points between the real-time structure model and the stable structure model. If there are, corresponding warning information and a risk detection period are generated. If not, no operation is required;

[0010] The structure risk assessment module is used to regularly generate a periodic structure model according to the risk detection period; obtain the elevation difference of the periodic structure model and generate an elevation change rate according to the elevation difference; generate corresponding risk information according to the elevation change rate, and generate a new risk detection period according to the elevation change rate. The new risk detection period overwrites the original risk detection period. Under the new risk detection period, a new elevation difference and elevation change rate are obtained, and then corresponding latest risk information is generated;

[0011] The warning module is used to process warning information and risk information.

[0012] Furthermore, the process of obtaining the elevation data of the structure and the elevation data of the terrain around the structure by INSAR technology includes:

[0013] Set a number of collection points on the structure and a number of collection points on the surrounding terrain;

[0014] The elevation data is the distance from the collection point along the plumb line direction to the absolute base surface;

[0015] Obtain the elevation data of the structure and the elevation data of the terrain around the structure by a synthetic aperture radar according to INSAR technology.

[0016] Furthermore, the process of establishing an initial structure model according to the elevation data of the structure and generating corresponding risk values for each region according to the elevation data of the terrain includes:

[0017] Establish a three-dimensional coordinate system, map the elevation data of the structure into the three-dimensional coordinate system to generate a number of object height coordinates; generate an initial structure model according to the number of object height coordinates;

[0018] Set a stress threshold;

[0019] Conduct stress analysis on the initial structure model according to the finite element analysis method to generate a stress distribution diagram. Divide the initial structure model into several regions according to the stress distribution diagram. Obtain the regions in the stress distribution diagram where the stress is less than the stress threshold, mark the regions as vulnerable regions, and then divide the regions outside the vulnerable regions into several stable regions;

[0020] Map the elevation data of the terrain into the three-dimensional coordinate system to generate a number of ground height coordinates, and generate a terrain model according to the number of ground height coordinates;

[0021] Divide the terrain model into several elevation units according to the ground elevation coordinates;

[0022] Obtain the slope points of each elevation unit, generate a terrain linear model based on the slope points, and then obtain the ground danger line through the terrain linear model;

[0023] Obtain the coordinate closest to the initial structure model on the ground danger line, and obtain the distance ε from this ground danger line to the initial structure model according to the coordinate;

[0024] Set a safety distance δ;

[0025] If the distance from the ground danger line to the initial structure model is less than the safety distance, mark this ground danger line as an obvious danger line; otherwise, mark this ground danger line as a hidden danger line;

[0026] Set the radiation radius value and the radiation main coefficient;

[0027] Obtain a number of coordinates corresponding to the ground danger line, and generate corresponding mapping coordinates according to the abscissa and ordinate corresponding to each of the number of coordinates; use the mapping coordinate as the center and the radiation radius value as the radius to generate the unit radiation range corresponding to this mapping coordinate; generate the radiation range of this ground danger line according to the unit radiation ranges corresponding to each mapping coordinate;

[0028] Assume that there is a coordinate A(x, y, z) on the ground danger line, then the mapping coordinate of coordinate A is (x, y);

[0029] The radiation range includes the main radiation range and the secondary radiation range;

[0030] Use the mapping coordinate as the center and the value after multiplying the radiation radius value by the radiation main coefficient as the radius to generate the corresponding main radiation range, and mark the range other than the main radiation range within the radiation range as the secondary radiation range;

[0031] Obtain the radiation range of the obvious danger line, and obtain the radiation range of the hidden danger line;

[0032] Generate the risk value corresponding to each area of the initial structure model according to the respective radiation ranges of the obvious danger line and the hidden danger line;

[0033] Set the vulnerability coefficient θ1 and the stability coefficient θ2, set the obvious danger coefficient σ1 and the hidden danger coefficient σ2, set the main radiation coefficient τ1 and the secondary radiation coefficient τ2;

[0034] Count the ground danger lines corresponding to each area within the radiation range, and then obtain the risk value Fx corresponding to each area; among them, the calculation formula for the risk value corresponding to the area is: where num 地危is the number of ground danger lines corresponding to the area. When the area is a vulnerable area, θ is θ1; when the area is a stable area, θ is θ2. When the ground danger line is an obvious danger line, σ is σ1; when the ground danger line is a hidden danger line, σ is σ2. When the area is within the main radiation range of the ground danger line, τ is τ1; when the area is within the secondary radiation range of the ground danger line, τ is τ2. When the area is within both the main and secondary radiation ranges, τ is τ1 + τ2, ε i is the distance from the ground danger line to the initial structure model.

[0035] Further, the process of obtaining the slope points of each ground height unit, generating a terrain linear model based on the slope points, and then obtaining the ground danger line through the terrain linear model is as follows:

[0036] Obtain the coordinates existing in the ground height unit, obtain the average value of the abscissa, the average value of the ordinate, and the average value of the vertical coordinate according to the coordinates, and generate slope points based on the average value of the abscissa, the average value of the ordinate, and the average value of the vertical coordinate; the slope point = (average value of the abscissa, average value of the ordinate, average value of the vertical coordinate);

[0037] Assume that there are coordinates 1(x1, y1, z1), coordinates 2(x2, y2, z2), coordinates 3(x3, y3, z3), ……, coordinates n (x n , y n , z n ), then according to

[0038] And according to the average value of the abscissa,

[0039] Step a1: Set an effective terrain range;

[0040] Step a2: When there are other slope points in the effective terrain range of the slope point, the slope point and the other slope points generate a terrain unit line;

[0041] Step a3: Obtain the corresponding terrain unit line of the other slope points according to the principle of obtaining the terrain unit line in step a2;

[0042] Step a4: Repeat step a3 until no terrain unit line can be obtained; generate a terrain line according to the terrain unit line, and then generate a terrain linear model according to the terrain line;

[0043] Set a segment length threshold, a slope threshold α1, and a slope threshold α2, where α1 < α2;

[0044] Obtain the coordinates of the two ends of the terrain line, and generate a terrain vector according to the coordinates;

[0045] Obtain the modulus length of the terrain vector. If the modulus length is greater than the segment length threshold, mark this terrain line as a dangerous line; otherwise, no operation is required.

[0046] Obtain the cosine values corresponding to adjacent terrain unit lines in the dangerous line and the cosine values corresponding to the terrain unit lines at both ends.

[0047] Obtain the cosine values corresponding to adjacent terrain unit lines in the dangerous line to generate a unit cosine set; obtain the number of cosine values in the unit cosine set that are within [α1, α2]. If the proportion of this number is greater than or equal to 50%, mark this dangerous line as a severely dangerous line.

[0048] If the cosine values corresponding to the terrain unit lines at both ends are greater than or equal to α2, mark this dangerous line as a severely dangerous line.

[0049] Furthermore, set a safety value. The process of generating a stable structure model based on the risk values and safety values of each region includes:

[0050] Generate corresponding safety values according to the stress analysis results of each region.

[0051] Set a deformation coefficient ω1 and a deformation coefficient ω2, where ω1 < ω2.

[0052] Generate corresponding deformation values Bx for each region according to the risk values and safety values of each region; the calculation formula for the deformation value corresponding to the region is: Bx = ω1 × safety value + ω2 × risk value.

[0053] Generate a stable structure model based on the deformation values and the initial structure model.

[0054] Furthermore, set a monitoring period. The process of regularly obtaining a real-time structure model according to the monitoring period includes:

[0055] Set a monitoring period.

[0056] According to the monitoring period, regularly obtain the elevation data of the structure through a synthetic aperture radar based on the INSAR technology, and generate a real-time structure model according to the process of obtaining the initial structure model through the elevation data.

[0057] Furthermore, detect whether there are contact points between the real-time structure model and the stable structure model. If so, the process of generating corresponding warning information and a risk detection period includes:

[0058] Generate an x - y cross-section according to the horizontal axis and the vertical axis in the three-dimensional coordinate system.

[0059] Obtain the z maximum value and the z minimum value in the real-time structure model.

[0060] Obtain the z maximum value and the z minimum value in the stable structure model.

[0061] Compare the maximum z value of the real-time structure model with the maximum z value of the stable structure model, and mark the larger of the two maximum z values as the scanned z outer value; compare the minimum z value of the real-time structure model with the minimum z value of the stable structure model, and mark the smaller of the two minimum z values as the scanned z inner value; generate a scanned z interval based on the scanned z outer value and the scanned z inner value; the scanned z interval = [scanned z inner value, scanned z outer value];

[0062] Scan the real-time structure model and the stable structure model according to the scanned z interval through the x-y cross-section. During the scanning process, obtain the coordinates of the real-time structure model and the stable structure model on the x-y cross-section respectively, and mark the coordinates as contour coordinates;

[0063] Obtain the area of the real-time structure model on the x-y cross-section according to the contour coordinates of the real-time structure model, and mark the area as the contour area;

[0064] According to the process of obtaining the contour area of the real-time structure model, obtain the corresponding contour area of the stable structure model;

[0065] If there is a coincidence between the contour coordinates of the real-time structure model and the contour coordinates of the stable structure model, it is determined that there is a contact point between the real-time structure model and the stable structure model;

[0066] If the contour area of the real-time structure model ≥ the contour area of the stable structure model, it is determined that there is a contact point between the real-time structure model and the stable structure model;

[0067] Otherwise, it is determined that there is no contact point between the real-time structure model and the stable structure model;

[0068] If there is a contact point between the real-time structure model and the stable structure model, generate a warning message and a risk detection period;

[0069] If there is no contact point between the real-time structure model and the stable structure model, no operation is required.

[0070] Furthermore, the process of regularly generating a periodic structure model according to the risk detection period, obtaining the elevation difference of the periodic structure model, and generating an elevation change rate based on the elevation difference includes:

[0071] Regularly obtain the elevation data of the structure according to the risk detection period through synthetic aperture radar based on the INSAR technology, map the elevation data to a three-dimensional coordinate system, generate a number of object height coordinates, and generate a periodic structure model according to the number of object height coordinates;

[0072] Generate a number of monitoring coordinates based on the heights of the objects, and number the monitoring coordinates, where the numbers can be Dot1, Dot2, ……, Dot n , where n is a natural number greater than 0;

[0073] Obtain the elevation data corresponding to each of the monitoring coordinates in the periodic structure model, and mark the elevation data as the initial elevation value; the Among them, the initial elevation value R is the initial elevation value corresponding to the Rth periodic structure model, and the elevation data R-i is the elevation data corresponding to the monitoring coordinate Dot i in the Rth periodic structure model;

[0074] Mark the first generated periodic structure model as the standard structure model;

[0075] Obtain the elevation difference of the periodic structure model according to the initial elevation value of the standard structure model;

[0076] The formula for obtaining the elevation difference is:

[0077] Generate an elevation change rate according to the elevation difference, the

[0078] Furthermore, the process of generating corresponding risk information according to the elevation change rate and generating a new risk detection period includes:

[0079] Obtain the maximum value of (elevation data R-i - elevation data 1-i ) / elevation data 1-i in the elevation change rate, and mark the maximum value as the elevation substitution value;

[0080] Set a change threshold φ1 and a change threshold φ2, where φ1 < φ2;

[0081] When the elevation substitution value < φ1, generate first-level risk information according to the elevation change rate, and generate a new risk detection period, and the duration of the new risk detection period is

[0082] When φ1 ≤ elevation substitution value < φ2, generate second-level risk information according to the elevation change rate, and generate a new risk detection period, and the duration of the new risk detection period is

[0083] When the elevation substitution value ≥ φ2, then generate third-level risk information according to the elevation change rate, and generate a new risk detection period, and the duration of the new risk detection period is

[0084] Further, the process of processing warning information and risk information includes:

[0085] When a warning information is received, then send a warning signal to the user, and inform the user that there is a problem with the stability of the structure, and let the user make corresponding handling;

[0086] When a first-level risk information is received, then send the elevation change rate to the user, and inform the user that the structure has a slight deformation and there is a slight safety risk problem, and let the user make corresponding handling;

[0087] When a second-level risk information is received, then send the elevation change rate to the user, and inform the user that the structure has a medium deformation and there is a general safety risk problem, and let the user make corresponding handling;

[0088] When a third-level risk information is received, then send the elevation change rate to the user, and inform the user that the structure has a serious deformation and there is a serious safety risk problem, and let the user make corresponding handling.

[0089] Compared with the prior art, the beneficial effects of the present invention are:

[0090] 1. Obtain the elevation data of the terrain, generate corresponding risk values for each area according to the elevation data of the terrain, and set safety values. Generate a stable structure model according to the risk values and safety values of each area. By taking the terrain around the structure as a safety risk factor affecting the structure, a model for judging whether the structure is stable is generated, and whether the structure is stable is judged more comprehensively and accurately according to the model;

[0091] 2. By obtaining the elevation change rate of the structure, that is, the change rate of the elevation data of the structure, the deformation of the structure is visualized, and a new risk detection period is generated according to the elevation change rate, so as to detect the change of the structure more timely. BRIEF DESCRIPTION OF THE DRAWINGS

[0092] Figure 1 is the schematic diagram of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0093] As Figure 1 shown, the structure safety risk assessment and warning safety monitoring system based on INSAR includes a control center, and the control center is connected with an elevation data acquisition module, an elevation data processing and analysis module, a structure risk detection module, a structure risk assessment module and a warning module;

[0094] The elevation data acquisition module is used to obtain the elevation data of the structure and the elevation data of the terrain around the structure through the INSAR technology;

[0095] The elevation data processing and analysis module is used to establish an initial structure model based on the elevation data of the structure, divide the initial structure model into several regions, generate corresponding risk values for each region according to the elevation data of the terrain; set a safety value, and generate a stable structure model according to the risk values and safety values of each region;

[0096] The structure risk detection module is used to set a monitoring period and regularly obtain a real-time structure model according to the monitoring period; detect whether there are contact points between the real-time structure model and the stable structure model. If there are, generate corresponding warning information and a risk detection period. If not, no operation is required;

[0097] The structure risk assessment module is used to regularly generate a periodic structure model according to the risk detection period; obtain the elevation difference of the periodic structure model, and generate an elevation change rate according to the elevation difference; generate corresponding risk information according to the elevation change rate, and generate a new risk detection period according to the elevation change rate, overwrite the original risk detection period with the new risk detection period, obtain a new elevation difference and elevation change rate under the new risk detection period, and then generate corresponding latest risk information;

[0098] The warning module is used to process warning information and risk information;

[0099] It should be further noted that in the specific implementation process, the process of obtaining the elevation data of the structure and the elevation data of the terrain around the structure through the INSAR technology includes:

[0100] Set several collection points on the structure and several collection points on the surrounding terrain;

[0101] The elevation data is the distance from the collection point along the plumb line direction to the absolute base surface;

[0102] Obtain the elevation data of the structure and the elevation data of the terrain around the structure through a synthetic aperture radar according to the INSAR technology;

[0103] It should be further noted that in the specific implementation process, the process of establishing an initial structure model based on the elevation data of the structure, dividing the initial structure model into several regions, and generating corresponding risk values for each region according to the elevation data of the terrain includes:

[0104] Establish a three-dimensional coordinate system, map the elevation data of the structure into the three-dimensional coordinate system to generate several object height coordinates; generate an initial structure model according to the several object height coordinates;

[0105] Set a stress threshold;

[0106] Perform stress analysis on the initial structure model according to the finite element analysis method to generate a stress distribution map. Divide the initial structure model into several regions according to the stress distribution map, obtain the regions in the stress distribution map where the stress is less than the stress threshold, mark the regions as vulnerable regions, and then divide the regions outside the vulnerable regions into several stable regions;

[0107] Map the elevation data of the terrain into a three-dimensional coordinate system to generate several ground elevation coordinates, and generate a terrain model according to the several ground elevation coordinates;

[0108] Divide the terrain model into several ground elevation units according to the ground elevation coordinates;

[0109] Obtain the slope points of each ground elevation unit, generate a terrain linear model according to the slope points, and then obtain the ground danger line through the terrain linear model;

[0110] Obtain the coordinates in the ground danger line that are closest to the initial structure model, and obtain the distance ε from this ground danger line to the initial structure model according to the coordinates;

[0111] Set a safety distance δ;

[0112] If the distance from the ground danger line to the initial structure model is less than the safety distance, mark this ground danger line as an obvious danger line; otherwise, mark this ground danger line as a hidden danger line;

[0113] Set a radiation radius value and a radiation main coefficient;

[0114] Obtain several coordinates corresponding to the ground danger line, and generate corresponding mapped coordinates according to the abscissas and ordinates of the several coordinates; Use the mapped coordinates as the center and the radiation radius value as the radius to generate the unit radiation range corresponding to this mapped coordinate; Generate the radiation range of this ground danger line according to the unit radiation ranges corresponding to each mapped coordinate;

[0115] Assume that there is a coordinate A(x, y, z) in the ground danger line, then the mapped coordinate of coordinate A is (x, y);

[0116] The radiation range includes a main radiation range and a secondary radiation range;

[0117] Generate a corresponding main radiation range with the mapped coordinates as the center and the value after multiplying the radiation radius value by the radiation main coefficient as the radius, and mark the range outside the main radiation range within the radiation range as the secondary radiation range;

[0118] Obtain the radiation range of the obvious danger line and the radiation range of the hidden danger line;

[0119] Generate the risk value corresponding to each region of the initial structure model according to the respective radiation ranges of the obvious danger line and the hidden danger line;

[0120] Set the vulnerability coefficient θ1 and the stability coefficient θ2, set the obvious danger coefficient σ1 and the hidden danger coefficient σ2, set the main radiation coefficient τ1 and the secondary radiation coefficient τ2;

[0121] Count the ground danger lines corresponding to each area within the radiation range, and then obtain the risk value Fx corresponding to each area; among them, the calculation formula for the risk value corresponding to the area is: Among them, num 地危 is the number of ground danger lines corresponding to this area. When this area is a vulnerable area, θ is θ1. When this area is a stable area, θ is θ2. When the ground danger line is an obvious danger line, σ is σ1. When the ground danger line is a hidden danger line, σ is σ2. When this area is within the main radiation range of this ground danger line, τ is τ1. When this area is within the secondary radiation range of this ground danger line, τ is τ2. When this area is within both the main radiation range and the secondary radiation range, τ is τ1 + τ2, and ε i is the distance from this ground danger line to the initial structure model;

[0122] Among them, the process of obtaining the slope points of each ground height unit, generating a terrain linear model based on the slope points, and then obtaining the ground danger lines through the terrain linear model is as follows:

[0123] Obtain the coordinates existing in the ground height unit, obtain the average value of the abscissa, the average value of the ordinate, and the average value of the vertical coordinate according to the coordinates, and generate slope points according to the average value of the abscissa, the average value of the ordinate, and the average value of the vertical coordinate; the slope point = (average value of the abscissa, average value of the ordinate, average value of the vertical coordinate);

[0124] Assume that there are coordinates 1(x1, y1, z1), coordinates 2(x2, y2, z2), coordinates 3(x3, y3, z3), ……, coordinates n (x n , y n , z n ) in the ground height unit, then according to

[0125] And according to the average value of the abscissa,

[0126] Step a1: Set the effective terrain range;

[0127] Step a2: When there are other slope points in the effective terrain range of the slope point, the slope point and the other slope points generate a terrain unit line;

[0128] Step a3: Obtain the terrain unit line corresponding to the other slope points according to the principle of obtaining the terrain unit line in step a2;

[0129] Step a4: Repeat step a3 until no more terrain unit lines can be obtained; generate terrain lines based on the terrain unit lines, and then generate a terrain linear model based on the terrain lines;

[0130] Set a segment length threshold, a slope threshold α1, and a slope threshold α2, where α1 < α2;

[0131] Obtain the coordinates of the two ends of the terrain line, and generate a terrain vector based on the coordinates;

[0132] Obtain the modulus length of the terrain vector. If the modulus length is greater than the segment length threshold, mark the terrain line as a dangerous line; otherwise, no operation is required;

[0133] Obtain the cosine values corresponding to adjacent terrain unit lines in the dangerous line and the cosine values corresponding to the terrain unit lines at the two ends;

[0134] Obtain the cosine values corresponding to adjacent terrain unit lines in the dangerous line to generate a unit cosine set; obtain the number of cosine values in the unit cosine set that are within [α1, α2]. If the proportion of the number is greater than or equal to 50%, mark the dangerous line as a severely dangerous line;

[0135] If the cosine values corresponding to the terrain unit lines at the two ends are greater than or equal to α2, mark the dangerous line as a severely dangerous line;

[0136] It should be further noted that in the specific implementation process, setting a safety value and generating a stable structure model based on the risk value and safety value of each region includes:

[0137] Generate corresponding safety values according to the stress analysis results of each region;

[0138] Set a deformation coefficient ω1 and a deformation coefficient ω2, where ω1 < ω2;

[0139] Generate the corresponding deformation value Bx for each region according to the risk value and safety value of each region; the calculation formula for the deformation value corresponding to the region is: Bx = ω1 × safety value + ω2 × risk value;

[0140] Generate a stable structure model based on the deformation value and the initial structure model;

[0141] It should be further noted that in the specific implementation process, setting a monitoring period and regularly obtaining a real-time structure model according to the monitoring period includes:

[0142] Set a monitoring period;

[0143] Regularly obtain the elevation data of the structure according to the monitoring period through a synthetic aperture radar based on the INSAR technology, and generate a real-time structure model according to the process of obtaining the initial structure model through the elevation data;

[0144] It should be further noted that in the specific implementation process, the process of detecting whether there are contact points between the real-time structure model and the stable structure model, and if so, generating corresponding warning information and a risk detection period includes:

[0145] Generate an x-y cross-section according to the horizontal and vertical axes in the three-dimensional coordinate system;

[0146] Obtain the maximum z value and the minimum z value in the real-time structure model;

[0147] Obtain the maximum z value and the minimum z value in the stable structure model;

[0148] Compare the maximum z value of the real-time structure model with the maximum z value of the stable structure model, and mark the larger maximum z value as the scanned z outer value; compare the minimum z value of the real-time structure model with the minimum z value of the stable structure model, and mark the smaller minimum z value as the scanned z inner value; generate a scanned z interval according to the scanned z outer value and the scanned z inner value; the scanned z interval = [scanned z inner value, scanned z outer value];

[0149] Scan the real-time structure model and the stable structure model according to the scanned z interval through the x-y cross-section. During the scanning process, obtain the coordinates of the real-time structure model and the stable structure model on the x-y cross-section in real time, and mark the coordinates as contour coordinates;

[0150] Obtain the area of the real-time structure model on the x-y cross-section according to the contour coordinates of the real-time structure model, and mark the area as the contour area;

[0151] According to the process of obtaining the contour area of the real-time structure model, obtain the corresponding contour area of the stable structure model;

[0152] If the contour coordinates of the real-time structure model are the same as those of the stable structure model, it is determined that there are contact points between the real-time structure model and the stable structure model;

[0153] If the contour area of the real-time structure model ≥ the contour area of the stable structure model, it is determined that there are contact points between the real-time structure model and the stable structure model;

[0154] Otherwise, it is determined that there are no contact points between the real-time structure model and the stable structure model;

[0155] If there are contact points between the real-time structure model and the stable structure model, generate warning information and a risk detection period;

[0156] If there are no contact points between the real-time structure model and the stable structure model, no operation is required;

[0157] It should be further noted that in the specific implementation process, the process of regularly generating a periodic structure model according to the risk detection period, obtaining the elevation difference of the periodic structure model, and generating an elevation change rate based on the elevation difference includes:

[0158] Regularly obtain the elevation data of the structure according to the INSAR technology through a synthetic aperture radar according to the risk detection period, map the elevation data into a three-dimensional coordinate system, generate a number of object height coordinates, and generate a periodic structure model based on the number of object height coordinates;

[0159] Generate a number of monitoring coordinates according to the number of object height coordinates, and number the number of monitoring coordinates. The numbering can be Dot1, Dot2,..., Dot n , where n is a natural number greater than 0;

[0160] Obtain the elevation data corresponding to each of the number of monitoring coordinates in the periodic structure model, and mark the elevation data as the initial elevation value; the where the initial elevation value R is the initial elevation value corresponding to the Rth periodic structure model, and the elevation data R-i is the elevation data corresponding to the monitoring coordinate Dot i in the Rth periodic structure model;

[0161] Mark the first generated periodic structure model as the standard structure model;

[0162] Obtain the elevation difference of the periodic structure model according to the initial elevation value of the standard structure model;

[0163] The formula for obtaining the elevation difference is:

[0164] Generate an elevation change rate based on the elevation difference, the

[0165] It should be further noted that in the specific implementation process, the process of generating corresponding risk information according to the elevation change rate and generating a new risk detection period according to the elevation change rate includes:

[0166] Obtain the maximum value of (elevation data R-i - elevation data 1-i ) / elevation data 1-i in the elevation change rate, and mark the maximum value as the elevation substitution value;

[0167] Set a change threshold φ1 and a change threshold φ2, where φ1 < φ2;

[0168] When the elevation substitution value < φ1, first-level risk information is generated according to the elevation change rate, and a new risk detection period is generated. The duration of the new risk detection period is

[0169] When φ1 ≤ elevation substitution value < φ2, second-level risk information is generated according to the elevation change rate, and a new risk detection period is generated. The duration of the new risk detection period is

[0170] When the elevation substitution value ≥ φ2, third-level risk information is generated according to the elevation change rate, and a new risk detection period is generated. The duration of the new risk detection period is

[0171] It should be further noted that in the specific implementation process, the process of processing warning information and risk information includes:

[0172] When a warning information is received, a warning signal is sent to the user, and it is informed that there is a problem with the stability of the structure, and the user shall make corresponding handling;

[0173] When a first-level risk information is received, the elevation change rate is sent to the user, and it is informed that the structure has undergone slight deformation and there is a slight safety risk problem, and the user shall make corresponding handling;

[0174] When a second-level risk information is received, the elevation change rate is sent to the user, and it is informed that the structure has undergone medium deformation and there is a general safety risk problem, and the user shall make corresponding handling;

[0175] When a third-level risk information is received, the elevation change rate is sent to the user, and it is informed that the structure has undergone severe deformation and there is a severe safety risk problem, and the user shall make corresponding handling;

[0176] The above embodiments are only used to illustrate the technical method of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those of ordinary skill in the art should understand that the technical method of the present invention can be modified or equivalently replaced without departing from the spirit and scope of the technical method of the present invention.

Claims

1. An INSAR-based safety risk assessment and early warning safety monitoring system for structures, including a control center, characterized in that, The control center is communicatively connected to an elevation data acquisition module, an elevation data processing and analysis module, a structure risk detection module, a structure risk assessment module, and an early warning module; The elevation data acquisition module is used to obtain the elevation data of the structure and the elevation data of the terrain around the structure through INSAR technology; The elevation data processing and analysis module is used to establish an initial structure model based on the elevation data of the structure, divide the initial structure model into several regions, and generate corresponding risk values for each region according to the elevation data of the terrain; Set a safety value, and generate a stable structure model according to the risk values and the safety value of each region; The structure risk detection module is used to set a monitoring period and regularly obtain a real-time structure model according to the monitoring period; Detect whether there are contact points between the real-time structure model and the stable structure model. If there are, generate corresponding early warning information and a risk detection period. If not, no operation is required; The structure risk assessment module is used to regularly generate a periodic structure model according to the risk detection period; Obtain the elevation difference of the periodic structure model, and generate an elevation change rate according to the elevation difference; generate corresponding risk information according to the elevation change rate, and generate a new risk detection period according to the elevation change rate. Overwrite the original risk detection period with the new risk detection period, obtain a new elevation difference and an elevation change rate under the new risk detection period, and then generate corresponding latest risk information; The early warning module is used to process early warning information and risk information; The process of establishing an initial structure model based on the elevation data of the structure, dividing the initial structure model into several regions, and generating corresponding risk values for each region according to the elevation data of the terrain includes: Establish a three-dimensional coordinate system, map the elevation data of the structure to the three-dimensional coordinate system to generate several object height coordinates; generate an initial structure model according to the several object height coordinates; Set a stress threshold; perform stress analysis on the initial structure model by the finite element analysis method to generate a stress distribution map, obtain the regions in the stress distribution map where the stress is less than the stress threshold, mark the regions as vulnerable regions, and then divide the regions outside the vulnerable regions into several stable regions; Map the elevation data of the terrain to the three-dimensional coordinate system to generate several ground height coordinates, and generate a terrain model according to the several ground height coordinates; Divide the terrain model into several ground height units according to the ground height coordinates; obtain the slope points of each ground height unit, generate a terrain linear model according to the slope points, and then obtain a ground danger line through the terrain linear model; Obtain the coordinates on the ground hazard line that are closest to the initial structure model, and obtain the distance from the ground hazard line to the initial structure model based on the coordinates ; Set a safety distance ; if the distance from the ground hazard line to the initial structure model is less than the safety distance, mark this ground hazard line as an obvious hazard line; otherwise, mark this ground hazard line as a hidden hazard line; Set a radiation radius value and a radiation main coefficient; Obtain several coordinates corresponding to the ground danger line, generate corresponding mapped coordinates according to the abscissa and ordinate corresponding to each of the several coordinates; use the mapped coordinates as the center and the radiation radius value as the radius to generate the unit radiation range corresponding to the mapped coordinates; generate the radiation range of the ground danger line according to the unit radiation ranges corresponding to each mapped coordinate; The radiation range includes a main radiation range and a secondary radiation range; Obtain the radiation range of the ground danger line according to the radiation radius value and the radiation main coefficient; Set the vulnerability coefficient and the stability coefficient , set the obvious danger coefficient and the hidden danger coefficient , set the main radiation coefficient and the secondary radiation coefficient ; Count the corresponding ground danger lines within the radiation range of each area, and then obtain the risk value Fx corresponding to each area. Among them, the calculation formula for the risk value corresponding to the area is: , where is the number of ground danger lines corresponding to the area. When the area is a vulnerable area, is . When the area is a stable area, is . When the ground danger line is an obvious danger line, is . When the ground danger line is a hidden danger line, is . When the area is within the main radiation range of the ground danger line, then is . When the area is within the secondary radiation range of the ground danger line, then is . When the area is within both the main radiation range and the secondary radiation range, then is , is the distance from the ground danger line to the initial structure model; The process of obtaining the slope points of each ground elevation unit, generating a terrain linear model based on the slope points, and then obtaining the dangerous terrain line through the terrain linear model is as follows: Obtain the coordinates existing in the ground elevation unit, obtain the average abscissa, average ordinate, and average vertical coordinate according to the coordinates, and generate slope points based on the average abscissa, average ordinate, and average vertical coordinate; Step a1: Set an effective terrain range; Step a2: When there are other slope points in the effective terrain range of a slope point, the slope point and the other slope points generate a terrain unit line; Step a3: Obtain the corresponding terrain unit lines of the other slope points according to the principle of obtaining the terrain unit lines in step a2; Step a4: Repeat step a3 until no more terrain unit lines can be obtained; generate a terrain line based on the terrain unit lines, and then generate a terrain linear model based on the terrain line; Set the segment length threshold and the slope threshold and the slope threshold , where ; Obtain the coordinates at both ends of the terrain line, and generate a terrain vector according to the coordinates; Obtain the modulus length of the terrain vector. If the modulus length is greater than the segment length threshold, mark the terrain line as a dangerous approaching line; otherwise, no operation is required; Obtain the corresponding cosine values of adjacent terrain unit lines in the dangerous approaching line and the cosine values of the terrain unit lines at both ends; Generate a unit cosine set according to the cosine values corresponding to adjacent terrain unit lines in the danger line; Obtain the number of cosine values in the unit cosine set that are located at If the proportion of the number is greater than or equal to 50%, then mark this danger line as a ground danger line; If the cosine values corresponding to the terrain unit lines at both the head and the tail are greater than or equal to , then mark this danger line as a ground danger line; Set a safety value. The process of generating a stable structure model based on the risk value and safety value of each area includes: Set the deformation coefficient and the deformation coefficient , where ; Generate the corresponding deformation value Bx for each area according to the risk value and safety value of each area; among them, the calculation formula for the deformation value corresponding to the area is: ; Generate a stable structure model based on the deformation value and the initial structure model; Set a monitoring period. The process of regularly obtaining a real-time structure model according to the monitoring period includes: Regularly obtain the elevation data of the structure through a synthetic aperture radar according to the monitoring period, and generate a real-time structure model based on the elevation data according to the process of obtaining the initial structure model; The process of detecting whether there are contact points between the real-time structure model and the stable structure model includes: Generate an x-y cross-section according to the horizontal axis and vertical axis in the three-dimensional coordinate system; Obtain the maximum z value and minimum z value in the real-time structure model; Obtain the maximum z value and minimum z value in the stable structure model; Generate a scanning z interval according to the maximum z value and minimum z value corresponding to the real-time structure model and the stable structure model respectively; scan the real-time structure model and the stable structure model through the x-y cross-section according to the scanning z interval, and obtain the coordinates of the real-time structure model and the stable structure model located on the x-y cross-section in real time during the scanning process, and mark the coordinates as contour coordinates; Obtain the area of the real-time structure model located on the x-y cross-section according to the contour coordinates of the real-time structure model, and mark the area as the contour area; according to the process of obtaining the contour area of the real-time structure model, obtain the corresponding contour area of the stable structure model; If the contour coordinates of the real-time structure model are the same as those of the stable structure model, it is determined that there are contact points between the real-time structure model and the stable structure model; If there is a contour area of the real-time structure model and a contour area of the stable structure model, it is determined that there is a contact point between the real-time structure model and the stable structure model; Otherwise, it is determined that there are no contact points between the real-time structure model and the stable structure model; The process of regularly generating a periodic structure model according to the risk detection period, obtaining the elevation difference of the periodic structure model, and generating an elevation change rate according to the elevation difference includes: Regularly obtain the elevation data of the structure according to the risk detection period, generate a number of object height coordinates based on the elevation data, and generate a periodic structure model based on the number of object height coordinates; Generate a number of monitoring coordinates based on the heights of the several objects, and number the several monitoring coordinates. The numbering can be , , ……, , where n is a natural number greater than 0; Obtain the elevation data corresponding to several monitoring coordinates, and mark the elevation data as the initial elevation value; the , where is the initial elevation value corresponding to the structure model in the R-th cycle, is the monitoring coordinate in the structure model in the R-th cycle corresponding elevation data; The calculation formula for obtaining the elevation difference is: ; Generate an elevation change rate based on the elevation difference, the ; The process of generating corresponding risk information according to the elevation change rate and generating a new risk detection period includes: Obtain the maximum value in the elevation change rate and mark the maximum value as the elevation generation value; Set the change threshold and the change threshold , where ; When occurs, first-level risk information is generated based on the elevation change rate, and a new risk detection period is generated. The duration of the new risk detection period is of the duration of the original risk detection period; When occurs, secondary risk information is generated based on the elevation change rate, and a new risk detection period is generated. The duration of the new risk detection period is of the duration of the original risk detection period; When it is the case, three-level risk information is generated according to the elevation change rate, and a new risk detection period is generated, and the duration of the new risk detection period is of the original risk detection period duration.

2. The INSAR-based safety risk assessment and early warning safety monitoring system for structures according to claim 1, characterized in that The process of obtaining the elevation data of the structure and the elevation data of the terrain around the structure through the INSAR technology includes: Set a number of collection points on the structure and a number of collection points on the surrounding terrain; The elevation data is the distance from the collection point along the vertical line direction to the absolute base surface; Obtain the elevation data of the structure and the elevation data of the terrain around the structure through the synthetic aperture radar according to the INSAR technology.

3. The INSAR-based safety risk assessment and early warning safety monitoring system for structures according to claim 2, characterized in that, The process of processing the warning information and the risk information includes: When receiving the warning information, send a warning signal to the user and inform the user that there is a problem with the stability of the structure, and let the user make corresponding processing; When receiving the first-level risk information, send the elevation change rate to the user and inform the user that the structure has a slight deformation and there is a slight safety risk problem, and let the user make corresponding processing; When receiving the second-level risk information, send the elevation change rate to the user and inform the user that the structure has a medium deformation and there is a general safety risk problem, and let the user make corresponding processing; When receiving the third-level risk information, send the elevation change rate to the user and inform the user that the structure has a serious deformation and there is a serious safety risk problem, and let the user make corresponding processing.

Images