Structure safety risk assessment and early warning safety monitoring system based on INSAR

Abstract

The invention discloses a structure safety risk assessment and early warning safety monitoring system based on INSAR, and relates to the technical field of risk assessment. The elevation data of the structure and the elevation data of the terrain are obtained through the INSAR technology; establishing an initial structure model according to the elevation data of the structure, dividing the initial structure model into a plurality of areas, and generating a risk value corresponding to each area according to the elevation data of the terrain; setting a safety value, and generating a stable structure model according to the risk value and the safety value; setting a monitoring period, and regularly obtaining a real-time structure model according to the monitoring period; detecting whether a contact point exists between the real-time structure model and the stable structure model or not, and if yes, generating early warning information and a risk detection period; periodically generating a periodic structure model according to the risk detection period; acquiring an elevation difference value of the periodic structure model, and generating an elevation change rate according to the elevation difference value; and generating risk information and a new risk detection period according to the elevation change rate.

Description

Technical Field

[0001] The invention relates to the technical field, and in particular to a structure safety risk assessment and early warning safety monitoring system based on INSAR. Background Art

[0002] In recent years, due to natural disasters, human disturbances, aging and other factors, the construction safety issues of structures around the world have become increasingly prominent. Traditional detection methods, such as visual inspection and structural monitoring, often require a lot of manpower, time and cost, and have limitations in identifying hidden dangers at an early stage.

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

[0004] However, the method of judging whether a structure has safety hazards only by monitoring whether the structure has deformed is too simplistic. How to incorporate the terrain around the structure into the factors affecting the safety hazards of the structure, so as to more comprehensively and accurately assess whether the structure has safety risks, is an issue we need to consider at present. Summary of the invention

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

[0006] The object of the present invention can be achieved by the following technical scheme: a structure safety risk assessment and early warning safety monitoring system based on INSAR, comprising a control center, wherein the control center is communicatively 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 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 according to the elevation data of the structure, divide the initial structure model into several areas, generate corresponding risk values ​​for each area according to the elevation data of the terrain; set a safety value, and generate a stable structure model according to the risk value and safety value of each area;

[0009] The structure risk detection module is used to set a monitoring cycle and regularly obtain a real-time structure model according to the monitoring cycle; detect whether there is a contact point between the real-time structure model and the stable structure model, and if so, generate corresponding warning information and risk detection cycle; if not, no operation is required;

[0010] The structure risk assessment module is used to regularly generate a periodic structure model according to a risk detection cycle; 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 cycle according to the elevation change rate, overwrite the original risk detection cycle with the new risk detection cycle, obtain new elevation difference and elevation change rate under the new risk detection cycle, and then generate corresponding latest risk information;

[0011] The early warning module is used to process early 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 means of INSAR technology includes:

[0013] Set up 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 acquisition point to the absolute base surface along the plumb line direction;

[0015] The elevation data of the structure and the elevation data of the terrain around the structure are obtained through synthetic aperture radar based on INSAR technology.

[0016] Furthermore, an initial structure model is established according to the elevation data of the structure, the initial structure model is divided into several areas, and the process of generating corresponding risk values ​​of each area 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, and generate a number of object height coordinates; generate an initial structure model according to the number of object height coordinates;

[0018] Setting stress thresholds;

[0019] Performing stress analysis on the initial structure model according to the finite element analysis method to generate a stress distribution map, dividing the initial structure model into a number of regions according to the stress distribution map, obtaining a region in the stress distribution map where the stress is less than a stress threshold, marking the region as a fragile region, and then dividing the region outside the fragile region into a number of stable regions;

[0020] Mapping the terrain elevation data into a three-dimensional coordinate system, generating a number of ground height coordinates, and generating a terrain model based on the number of ground height coordinates;

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

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

[0023] Obtain the coordinates of the hazard line closest to the initial structure model, and obtain the distance ε from the hazard line to the initial structure model according to the coordinates;

[0024] Set the safety distance δ;

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

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

[0027] Acquire several coordinates corresponding to the ground-hazard line, and generate corresponding mapping coordinates according to the horizontal coordinates and vertical coordinates corresponding to the several coordinates; generate the unit radiation range corresponding to the mapping coordinate with the mapping coordinate as the center and the radiation diameter value as the radius; generate the radiation range of the ground-hazard line according to the unit radiation range corresponding to each mapping coordinate;

[0028] Assuming that there is a coordinate A (x, y, z) in the geohazard line, the coordinates of coordinate A are (x, y);

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

[0030] The projection coordinates are taken as the center of the circle and the value obtained by multiplying the radiation diameter value by the radiation main coefficient is taken as the radius to generate the corresponding main radiation range, and the range other than the main radiation range within the radiation range is marked as the secondary radiation range;

[0031] Get the radiation range of the obvious danger line and 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 corresponding radiation range of the obvious danger line and the hidden danger line;

[0033] Set the fragility coefficient θ1 and the stability coefficient θ2, set the apparent danger coefficient σ1 and the latent danger coefficient σ2, set the primary radiation coefficient τ1 and the secondary radiation coefficient τ2;

[0034] The corresponding ground hazard lines of each area within the radiation range are counted, and then the corresponding risk value Fx of each area is obtained; the calculation formula of the corresponding risk value of the area is: Among them, num 地危is the number of ground-danger lines corresponding to the area. When the area is a fragile 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 radiation range and the secondary radiation range, τ is τ1+τ2, ε i It is the distance from the danger line to the initial structure model.

[0035] Furthermore, the slope points of each ground height unit are obtained, a terrain linear model is generated according to the slope points, and then the process of obtaining the ground hazard line through the terrain linear model is as follows:

[0036] Obtaining coordinates existing in the ground height unit, obtaining an average value of the horizontal coordinate, an average value of the vertical coordinate and an average value of the vertical coordinate according to the coordinates, and generating a slope point according to the average value of the horizontal coordinate, the average value of the vertical coordinate and the average value of the vertical coordinate; the slope point = (average value of the horizontal coordinate, average value of the vertical coordinate, 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 in the ground height unit n (x n ,y n , z n ), then according to

[0038] And according to the horizontal axis average value,

[0039] Step a1: Set the 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 terrain unit lines;

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

[0042] Step a4: repeat step a3 until no terrain unit lines can be obtained; generate terrain lines according to the terrain unit lines, and then generate terrain linear models according to the terrain lines;

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

[0044] Obtaining coordinates of the beginning and end of the terrain line, and generating a terrain vector according to the coordinates;

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

[0046] Obtain the cosine values ​​of the adjacent terrain unit lines and the cosine values ​​of the terrain unit lines at both ends of the danger line;

[0047] Obtain the cosine values ​​corresponding to the adjacent terrain unit lines in the danger line to generate a unit cosine set; obtain the number of cosine values ​​located in [α1, α2] in the unit cosine set, and if the proportion of the cosine values ​​is greater than or equal to 50%, mark the danger line as a terrain danger line;

[0048] If the corresponding cosine values ​​of the unit terrain lines at the beginning and end are greater than or equal to α2, the danger line is marked as the danger line.

[0049] Furthermore, the process of setting the safety value and generating the stable structure model according to the risk value and safety value of each area includes:

[0050] Generate corresponding safety values ​​based on stress analysis results of each area;

[0051] Set the deformation coefficient ω1 and the deformation coefficient ω2, where ω1<ω2;

[0052] The corresponding deformation value Bx of each area is generated according to the risk value and safety value of each area; wherein the calculation formula of the corresponding deformation value of the area is: Bx = ω1 × safety value + ω2 × risk value;

[0053] A stable structure model is generated according to the deformation value and the initial structure model.

[0054] Furthermore, a monitoring cycle is set, and the process of regularly obtaining the real-time structure model according to the monitoring cycle includes:

[0055] Set monitoring period;

[0056] According to the monitoring period, the elevation data of the structure is regularly acquired by synthetic aperture radar according to the INSAR technology, and the real-time structure model is generated according to the process of acquiring the initial structure model through the elevation data.

[0057] Furthermore, it is detected whether there is a contact point between the real-time structure model and the stable structure model. If so, the process of generating corresponding warning information and risk detection cycle includes:

[0058] Generate an xy cross section according to the horizontal and vertical axes in a three-dimensional coordinate system;

[0059] Obtain the z maximum and z minimum values ​​in the real-time structure model;

[0060] Obtain the maximum and minimum z values ​​in the stable structure model;

[0061] Compare the z maximum value of the real-time structure model with the z maximum value of the stable structure model, and mark the larger z maximum value as the scan z outside value; compare the z minimum value of the real-time structure model with the z minimum value of the stable structure model, and mark the smaller z minimum value as the scan z inside value; generate a scan z interval according to the scan z outside value and the scan z inside value; the scan z interval = [scan z inside value, scan z outside value];

[0062] Scanning the real-time structure model and the stable structure model according to the scanning z interval through the xy cross section, acquiring the coordinates of the real-time structure model and the stable structure model respectively located on the xy cross section in real time during the scanning process, and marking the coordinates as contour coordinates;

[0063] Acquire the area of ​​the real-time structure model at the xy 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, the contour area corresponding to the stable structure model is obtained;

[0065] If the contour coordinates of the real-time structure model are consistent with 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 is greater than or equal to 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 are contact points between the real-time structure model and the stable structure model, early warning information and risk detection cycles are generated;

[0069] If there are no contact points between the live structure model and the stable structure model, no action is required.

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

[0071] According to the risk detection period, the elevation data of the structure is regularly acquired by synthetic aperture radar according to the INSAR technology, the elevation data is mapped into a three-dimensional coordinate system, a plurality of object height coordinates are generated, and a periodic structure model is generated according to the plurality of object height coordinates;

[0072] Generate a plurality of monitoring coordinates according to the plurality of object height coordinates, and number the plurality of monitoring coordinates, which may be Dot1, Dot2, ..., Dot n , where n is a natural number greater than 0;

[0073] Acquire the elevation data corresponding to a plurality of monitoring coordinates in the periodic structure model, and mark the elevation data as the initial elevation value; Among them, the initial value of elevation R is the initial value of the elevation corresponding to the structure model in the Rth period, the elevation data R-i The monitoring coordinate Dot in the structure model of the Rth period i Corresponding elevation data;

[0074] Mark the first periodic structure model generated as a 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 calculation formula for obtaining the elevation difference is:

[0077] The elevation change rate is generated according to the elevation difference.

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

[0079] Obtaining the elevation change rate (elevation data R-i - Elevation data 1-i ) / Elevation data 1-i The maximum value of the elevation is marked as the elevation value;

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

[0081] When the elevation value is <φ1, the first-level risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the length of the new risk detection cycle is twice the length of the original risk detection cycle.

[0082] When φ1≤elevation value<φ2, secondary risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the duration of the new risk detection cycle is twice the duration of the original risk detection cycle.

[0083] When the elevation value is ≥φ2, the third-level risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the length of the new risk detection cycle is twice the length of the original risk detection cycle.

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

[0085] When receiving warning information, it sends a warning signal to the user and informs the user that there is a problem with the stability of the structure, and the user takes corresponding measures;

[0086] When receiving the first-level risk information, the elevation change rate is sent to the user, and the user is informed that the structure has slightly deformed and there is a slight safety risk, and the user is required to handle it accordingly;

[0087] When receiving the secondary risk information, the elevation change rate is sent to the user, and the user is informed that the structure has undergone moderate deformation and there is a general safety risk problem, and the user is required to handle it accordingly;

[0088] When the third-level risk information is received, the elevation change rate will be sent to the user, and the user will be informed that the structure has undergone serious deformation and there is a serious safety risk, and the user will take appropriate measures.

[0089] Compared with the prior art, the present invention has the following beneficial effects:

[0090] 1. Obtain the elevation data of the terrain, generate the corresponding risk value of each area according to the elevation data of the terrain, and set the safety value. Generate a stable structure model according to the risk value and safety value 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. Based on the model, the stability of the structure can be judged more comprehensively and accurately;

[0091] 2. By obtaining the elevation change rate of the structure, that is, the rate of change of the structure's elevation data, the deformation of the structure can be visualized, and a new risk detection cycle is generated based on the elevation change rate, so that changes in the structure can be detected more timely. BRIEF DESCRIPTION OF THE DRAWINGS

[0092] Figure 1 It is a schematic diagram of the present invention. DETAILED DESCRIPTION

[0093] like Figure 1 As shown, the structure safety risk assessment and early 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 an early 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 INSAR technology;

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

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

[0097] The structure risk assessment module is used to regularly generate a periodic structure model according to a risk detection cycle; 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 cycle according to the elevation change rate, overwrite the original risk detection cycle with the new risk detection cycle, obtain new elevation difference and elevation change rate under the new risk detection cycle, and then generate corresponding latest risk information;

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

[0099] It should be further explained 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 INSAR technology includes:

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

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

[0102] Acquire the elevation data of the structure and the elevation data of the terrain around the structure through synthetic aperture radar based on INSAR technology;

[0103] It should be further explained that, in the specific implementation process, the initial structure model is established according to the elevation data of the structure, the initial structure model is divided into several areas, and the process of generating the corresponding risk value of each area 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, and generate a number of object height coordinates; generate an initial structure model according to the number of object height coordinates;

[0105] Setting stress thresholds;

[0106] Performing stress analysis on the initial structure model according to the finite element analysis method to generate a stress distribution map, dividing the initial structure model into a number of regions according to the stress distribution map, obtaining a region in the stress distribution map where the stress is less than a stress threshold, marking the region as a fragile region, and then dividing the region outside the fragile region into a number of stable regions;

[0107] Mapping the terrain elevation data into a three-dimensional coordinate system, generating a number of ground height coordinates, and generating a terrain model based on the number of ground height coordinates;

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

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

[0110] Obtain the coordinates of the hazard line closest to the initial structure model, and obtain the distance ε from the hazard line to the initial structure model according to the coordinates;

[0111] Set the safety distance δ;

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

[0113] Set the radiation path value and radiation main coefficient;

[0114] Acquire several coordinates corresponding to the ground-hazard line, and generate corresponding mapping coordinates according to the horizontal coordinates and vertical coordinates corresponding to the several coordinates; generate the unit radiation range corresponding to the mapping coordinate with the mapping coordinate as the center and the radiation diameter value as the radius; generate the radiation range of the ground-hazard line according to the unit radiation range corresponding to each mapping coordinate;

[0115] Assuming that there is a coordinate A (x, y, z) in the geohazard line, the coordinates of coordinate A are (x, y);

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

[0117] The projection coordinates are taken as the center of the circle and the value obtained by multiplying the radiation diameter value by the radiation main coefficient is taken as the radius to generate the corresponding main radiation range, and the range other than the main radiation range within the radiation range is marked as the secondary radiation range;

[0118] Get 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 area of ​​the initial structure model according to the corresponding radiation range of the obvious danger line and the hidden danger line;

[0120] Set the fragility coefficient θ1 and the stability coefficient θ2, set the apparent danger coefficient σ1 and the latent danger coefficient σ2, set the primary radiation coefficient τ1 and the secondary radiation coefficient τ2;

[0121] The corresponding ground hazard lines of each area within the radiation range are counted, and then the corresponding risk value Fx of each area is obtained; the calculation formula of the corresponding risk value of the area is: Among them, num 地危 is the number of ground-danger lines corresponding to the area. When the area is a fragile 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 radiation range and the secondary radiation range, τ is τ1+τ2, ε i The distance from the danger line to the initial structure model;

[0122] 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 hazard line through the terrain linear model is as follows:

[0123] Obtaining coordinates existing in the ground height unit, obtaining an average value of the horizontal coordinate, an average value of the vertical coordinate and an average value of the vertical coordinate according to the coordinates, and generating a slope point according to the average value of the horizontal coordinate, the average value of the vertical coordinate and the average value of the vertical coordinate; the slope point = (average value of the horizontal coordinate, average value of the vertical coordinate, 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 in the ground height unit n (x n ,y n , z n ), then according to

[0125] And according to the horizontal axis average value,

[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 terrain unit lines;

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

[0129] Step a4: repeat step a3 until no terrain unit lines can be obtained; generate terrain lines according to the terrain unit lines, and then generate terrain linear models according to the terrain lines;

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

[0131] Obtaining coordinates of the beginning and end of the terrain line, and generating a terrain vector according to the coordinates;

[0132] Obtain the modulus of the terrain vector. If the modulus 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 ​​of the adjacent terrain unit lines and the cosine values ​​of the terrain unit lines at both ends of the danger line;

[0134] Obtain the cosine values ​​corresponding to the adjacent terrain unit lines in the danger line to generate a unit cosine set; obtain the number of cosine values ​​located in [α1, α2] in the unit cosine set, and if the proportion of the cosine values ​​is greater than or equal to 50%, mark the danger line as a terrain danger line;

[0135] If the corresponding cosine value of the topographic unit line at the first and last ends is greater than or equal to α2, the danger line is marked as the ground danger line;

[0136] It should be further explained that, in the specific implementation process, the process of setting the safety value and generating a stable structure model according to the risk value and safety value of each area includes:

[0137] Generate corresponding safety values ​​based on stress analysis results of each area;

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

[0139] The corresponding deformation value Bx of each area is generated according to the risk value and safety value of each area; wherein the calculation formula of the corresponding deformation value of the area is: Bx = ω1 × safety value + ω2 × risk value;

[0140] Generate a stable structure model according to the deformation value and the initial structure model;

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

[0142] Set monitoring period;

[0143] acquiring elevation data of the structure by synthetic aperture radar according to the INSAR technology at regular intervals according to the monitoring period, and generating a real-time structure model by the elevation data according to the process of acquiring the initial structure model;

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

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

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

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

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

[0149] Scanning the real-time structure model and the stable structure model according to the scanning z interval through the xy cross section, acquiring the coordinates of the real-time structure model and the stable structure model respectively located on the xy cross section in real time during the scanning process, and marking the coordinates as contour coordinates;

[0150] Acquire the area of ​​the real-time structure model at the xy 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, the contour area corresponding to the stable structure model is obtained;

[0152] If the contour coordinates of the real-time structure model are consistent with 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;

[0153] If the contour area of ​​the real-time structure model is greater than or equal to 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;

[0154] Otherwise, it is determined that there is no contact point 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, early warning information and risk detection cycles are generated;

[0156] If there are no contact points between the live structure model and the stable structure model, no action is required;

[0157] It should be further explained 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 the elevation change rate according to the elevation difference includes:

[0158] According to the risk detection period, the elevation data of the structure is regularly acquired by synthetic aperture radar according to the INSAR technology, the elevation data is mapped into a three-dimensional coordinate system, a plurality of object height coordinates are generated, and a periodic structure model is generated according to the plurality of object height coordinates;

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

[0160] Acquire the elevation data corresponding to a plurality of monitoring coordinates in the periodic structure model, and mark the elevation data as the initial elevation value; The initial value of elevation R is the initial value of the elevation corresponding to the structure model in the Rth period, the elevation data R-i The monitoring coordinate Dot in the structure model of the Rth period i Corresponding elevation data;

[0161] Mark the first periodic structure model generated as a 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 calculation formula for obtaining the elevation difference is:

[0164] The elevation change rate is generated according to the elevation difference.

[0165] It should be further explained 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 cycle according to the elevation change rate includes:

[0166] Obtaining the elevation change rate (elevation data R-i - Elevation data 1-i ) / Elevation data 1-i The maximum value of the elevation is marked as the elevation value;

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

[0168] When the elevation value is <φ1, the first-level risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the length of the new risk detection cycle is twice the length of the original risk detection cycle.

[0169] When φ1≤elevation value<φ2, secondary risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the duration of the new risk detection cycle is twice the duration of the original risk detection cycle.

[0170] When the elevation value is ≥φ2, the third-level risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the length of the new risk detection cycle is twice the length of the original risk detection cycle.

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

[0172] When receiving warning information, it sends a warning signal to the user and informs the user that there is a problem with the stability of the structure, and the user takes corresponding measures;

[0173] When receiving the first-level risk information, the elevation change rate is sent to the user, and the user is informed that the structure has slightly deformed and there is a slight safety risk, and the user is required to handle it accordingly;

[0174] When receiving the secondary risk information, the elevation change rate is sent to the user, and the user is informed that the structure has undergone moderate deformation and there is a general safety risk problem, and the user is required to handle it accordingly;

[0175] When receiving the third-level risk information, the elevation change rate is sent to the user, and the user is informed that the structure has undergone serious deformation and there is a serious safety risk problem, and the user is required to handle it accordingly;

[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 skilled in the art should understand that the technical method of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical method of the present invention.

Claims

1. The INSAR-based structure safety risk assessment and early warning safety monitoring system includes a control center, which is characterized by: The control center is communicatively 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 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 according to the elevation data of the structure, divide the initial structure model into several areas, and generate corresponding risk values ​​for each area according to the elevation data of the terrain; Set safety values ​​and generate stable structure models based on the risk values ​​and safety values ​​of each area; 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 is a contact point between the real-time structure model and the stable structure model. If so, generate corresponding warning information and risk detection cycle. If not, no operation is required. The structure risk assessment module is used to regularly generate a periodic structure model according to a risk detection cycle; Obtain the elevation difference of the periodic structure model, and generate the elevation change rate based on the elevation difference; generate corresponding risk information based on the elevation change rate, and generate a new risk detection period based on the elevation change rate, overwrite the original risk detection period with the new risk detection period, obtain new elevation difference and elevation change rate under the new risk detection period, and then generate the corresponding latest risk information; The early warning module is used to process early warning information and risk information.

2. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 1 is characterized in that: The process of obtaining the elevation data of the structure and the elevation data of the terrain around the structure through INSAR technology includes: Set up 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 acquisition point to the absolute base surface along the plumb line direction; The elevation data of the structure and the elevation data of the terrain around the structure are obtained through synthetic aperture radar based on INSAR technology.

3. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 2 is characterized in that: The process of establishing an initial structure model based on the elevation data of the structure, dividing the initial structure model into several areas, and generating the corresponding risk value of each area based on the elevation data of the terrain includes: Establish a three-dimensional coordinate system, map the elevation data of the structure into the three-dimensional coordinate system, and generate a number of object height coordinates; generate an initial structure model according to the number of object height coordinates; Setting a stress threshold; performing stress analysis on the initial structure model according to the finite element analysis method to generate a stress distribution map, obtaining an area in the stress distribution map where the stress is less than the stress threshold, marking the area as a fragile area, and then dividing the area outside the fragile area into a number of stable areas; Mapping the terrain elevation data into a three-dimensional coordinate system, generating a number of ground height coordinates, and generating a terrain model based on the number of ground height coordinates; The terrain model is divided into several ground height units according to the ground height coordinates; the slope points of each ground height unit are obtained, a terrain linear model is generated according to the slope points, and then the ground hazard line is obtained through the terrain linear model; Obtain the coordinates of the hazard line closest to the initial structure model, and obtain the distance ε from the hazard line to the initial structure model according to the coordinates; Set the safety distance δ; if the distance from the ground danger line to the initial structure model is less than the safety distance, the ground danger line is marked as an obvious danger line; otherwise, the ground danger line is marked as a hidden danger line; Set the radiation path value and radiation main coefficient; The radiation range includes a main radiation range and a secondary radiation range; Obtain the radiation range of the ground hazard line based on the radiation diameter value and the main radiation coefficient; Set the fragility coefficient θ1 and the stability coefficient θ2, set the apparent danger coefficient σ1 and the latent danger coefficient σ2, set the primary radiation coefficient τ1 and the secondary radiation coefficient τ2; The corresponding ground hazard lines of each area within the radiation range are counted, and then the corresponding risk value Fx of each area is obtained; the calculation formula of the corresponding risk value of the area is: Among them, num 地危 is the number of ground-danger lines corresponding to the area. When the area is a fragile 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 radiation range and the secondary radiation range, τ is τ1+τ2, ε i It is the distance from the danger line to the initial structure model.

4. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 3 is characterized in that: 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 hazard line through the terrain linear model is as follows: Acquire the coordinates existing in the ground height unit, acquire the average value of the horizontal coordinate, the average value of the vertical coordinate and the average value of the vertical coordinate according to the coordinates, and generate the slope point according to the average value of the horizontal coordinate, the average value of the vertical coordinate and the average value of the vertical coordinate; Step a1: Set the effective terrain range; 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 terrain unit lines; Step a3: obtaining the terrain unit lines corresponding to the other slope points according to the principle of obtaining the terrain unit lines in step a2; Step a4: repeat step a3 until no terrain unit lines can be obtained; generate terrain lines according to the terrain unit lines, and then generate terrain linear models according to the terrain lines; Set the segment length threshold, slope threshold α1 and slope threshold α2, where α1<α2; Obtaining coordinates of the beginning and end of the terrain line, and generating a terrain vector according to the coordinates; Obtain the modulus of the terrain vector. If the modulus is greater than the segment length threshold, mark the terrain line as a dangerous line. Otherwise, no operation is required. Obtain the cosine values ​​of the adjacent terrain unit lines and the cosine values ​​of the terrain unit lines at both ends of the danger line; Generate a unit cosine set according to the corresponding cosine values ​​of the adjacent terrain unit lines in the danger line; obtain the number of cosine values ​​in the unit cosine set located in [α1, α2], and if the proportion of the said number is greater than or equal to 50%, mark the danger line as a terrain danger line; if the corresponding cosine values ​​of the terrain unit lines at the head and tail ends are greater than or equal to α2, mark the danger line as a terrain danger line.

5. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 4 is characterized in that: The process of setting safety values ​​and generating a stable structure model based on the risk values ​​and safety values ​​of each area includes: Set the deformation coefficient ω1 and the deformation coefficient ω2, where ω1<ω2; The corresponding deformation value Bx of each area is generated according to the risk value and safety value of each area; wherein the calculation formula of the corresponding deformation value of the area is: Bx = ω1 × safety value + ω2 × risk value; A stable structure model is generated according to the deformation value and the initial structure model.

6. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 5 is characterized in that: The process of setting a monitoring cycle and regularly obtaining a real-time structure model according to the monitoring cycle includes: The elevation data of the structure is acquired periodically through synthetic aperture radar according to the monitoring period, and a real-time structure model is generated through the elevation data according to the process of acquiring the initial structure model.

7. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 6 is characterized in that: The process of detecting whether there are contact points between the live structure model and the stable structure model includes: Generate an xy cross section according to the horizontal and vertical axes in a three-dimensional coordinate system; Obtain the z maximum and z minimum values ​​in the real-time structure model; Obtain the maximum and minimum z values ​​in the stable structure model; Generate a scanning z interval according to the z maximum and z minimum values ​​corresponding to the real-time structure model and the stable structure model respectively; scan the real-time structure model and the stable structure model according to the scanning z interval through the xy cross section, obtain the coordinates of the real-time structure model and the stable structure model respectively located on the xy cross section in real time during the scanning process, and mark the coordinates as contour coordinates; According to the contour coordinates of the real-time structure model, the area of ​​the real-time structure model at the xy cross section is obtained, and the area is marked as the contour area; according to the process of obtaining the contour area of ​​the real-time structure model, the contour area corresponding to the stable structure model is obtained; If the contour coordinates of the real-time structure model are consistent with 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; If the contour area of ​​the real-time structure model is greater than or equal to 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; Otherwise, it is determined that there is no contact point between the real-time structure model and the stable structure model.

8. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 7 is characterized in that: The process of regularly generating a periodic structure model according to a 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 acquiring elevation data of the structure according to the risk detection period, generating a number of object height coordinates according to the elevation data, and generating a periodic structure model according to the number of object height coordinates; Generate a plurality of monitoring coordinates according to the plurality of object height coordinates, and number the plurality of monitoring coordinates, which may be Dot1, Dot2, ..., Dot n , where n is a natural number greater than 0; Obtaining elevation data corresponding to a plurality of monitoring coordinates, and marking the elevation data as initial elevation values; The initial value of elevation R is the initial value of the elevation corresponding to the structure model in the Rth period, the elevation data R-i The monitoring coordinate Dot in the structure model of the Rth period i Corresponding elevation data; The calculation formula for obtaining the elevation difference is: The elevation change rate is generated according to the elevation difference.

9. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 8 is characterized in that: The process of generating corresponding risk information according to the elevation change rate and generating a new risk detection cycle according to the elevation change rate includes: Obtaining a maximum value of the elevation change rate, and marking the maximum value as an elevation proxy value; Set the change threshold φ1 and the change threshold φ2, where φ1<φ2; when When the elevation change rate is reached, the first-level risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the duration of the new risk detection cycle is the original risk detection cycle. when When the elevation change rate is reached, the secondary risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the duration of the new risk detection cycle is the original risk detection cycle. when When the elevation change rate is , the third-level risk information is generated according to the elevation change rate, and a new risk detection cycle is generated, and the duration of the new risk detection cycle is the original risk detection cycle.

10. The INSAR-based structure safety risk assessment and early warning safety monitoring system according to claim 9 is characterized in that: The process of handling early warning information and risk information includes: When receiving warning information, it sends a warning signal to the user and informs the user that there is a problem with the stability of the structure, and the user takes corresponding measures; When receiving the first-level risk information, the elevation change rate is sent to the user, and the user is informed that the structure has slightly deformed and there is a slight safety risk, and the user is required to handle it accordingly; When receiving the secondary risk information, the elevation change rate is sent to the user, and the user is informed that the structure has undergone moderate deformation and there is a general safety risk problem, and the user is required to handle it accordingly; When the third-level risk information is received, the elevation change rate will be sent to the user, and the user will be informed that the structure has undergone serious deformation and there is a serious safety risk, and the user will take appropriate measures.

Images