A method for monitoring horizontal time-series displacement of super high-rise buildings based on heterogeneous atmospheric correction

Abstract

The present invention relates to a method and system for monitoring horizontal time-series displacement of super high-rise buildings based on inhomogeneous atmosphere correction, and belongs to the field of super high-rise building health monitoring. The method is based on the relative position relationship between the geometric features of the super high-rise building to be monitored and the deployment of the GB-SAR instrument, calculates the change of atmospheric refractive index under the influence of inhomogeneous meteorological factors, and constructs a mathematical expression of horizontal displacement after optimization and correction of inhomogeneous atmospheric parameters; based on fitting methods such as the least square method or the quadratic curve method, the monitoring error correction number caused by the energy attenuation of the radar wave with the increase of the transmission distance is obtained, and then a non-homogeneous atmospheric parameter correction optimization model taking into account the energy attenuation of the radar wave is constructed. The model takes into account the GB-SAR radar wave energy attenuation and the disturbance of inhomogeneous meteorological factors, so it can reduce the influence of inhomogeneous meteorological factors and radar wave energy attenuation on the monitoring accuracy of the horizontal time-series displacement of GB-SAR super high-rise buildings, thereby improving the monitoring accuracy of the horizontal time-series displacement of super high-rise buildings.

Description

Technical Field

[0001] The present invention relates to the technical field of super high-rise building health monitoring, and in particular to a super high-rise building horizontal time-series displacement monitoring method based on inhomogeneous atmosphere correction. Background Art

[0002] After super high-rise buildings are put into use, as the use time increases, due to the influence of complex coupling factors such as abnormal loads, material aging, and component defects, the structural performance such as bearing capacity and deformation capacity is prone to decrease. In order to grasp the health status of super high-rise buildings in real time, it is necessary to conduct high-frequency dynamic monitoring of their structural status. The horizontal time-series displacement of super high-rise buildings can be used as an important basis for structural damage identification and structural status performance assessment. Ground-Based Synthetic Aperture Radar (GB-SAR) can achieve the purpose of dynamic monitoring of the horizontal displacement of super high-rise buildings by transmitting and receiving radar waves. It has the advantages of easy installation and operation, high sampling frequency (1000HZ), high spatial resolution (0.01mm), non-contact, all-weather operation, and large amount of information acquisition. However, since the parameters such as temperature, humidity, and pressure in the atmospheric troposphere are constantly changing in time and space, the propagation path of GB-SAR radar waves is very easy to be bent, resulting in a large amount of atmospheric disturbance noise in the monitoring signal, which greatly restricts the improvement of GB-SAR measurement accuracy. Therefore, in order to obtain high-precision horizontal time-series displacement of super-high-rise buildings, it is necessary to correct the disturbance errors caused by meteorological factors during the monitoring process.

[0003] At present, the main methods for atmospheric disturbance error correction are as follows:

[0004] (1) The manual control point correction method first obtains the distribution characteristics of meteorological factors in the field of view based on the displacement change information of stable target points or deployed reflectors in the study area, and then combines the positional relationship between the monitoring points and the control points to achieve the purpose of reducing the impact of meteorological factors on GB-SAR monitoring accuracy. This method does not require information such as temperature, humidity, and air pressure, and is simple to operate, but it is only applicable to single, long-term static observations. When there are no obvious control points or the control points are difficult to deploy, this method is difficult to implement. Therefore, the manual control point method is not suitable for regular monitoring of horizontal time-series displacements of super-high-rise buildings.

[0005] (2) Persistent Scatterer (PS) correction method: By selecting PS points with high phase stability and high coherence in the scene, a permanent scatterer network is constructed. The atmospheric correction model is constructed with the temporal phase value of the PS point as a reference to correct the influence of atmospheric errors on the remaining points, thereby obtaining accurate target deformation values. The effectiveness of this method is affected by the accuracy of the PS point selection and requires a large amount of image data within the research scope. Due to the limitations of imaging angle and maximum monitoring range, it is difficult to extract effective PS points in the GB-SAR super-high-rise building horizontal temporal displacement monitoring process. Therefore, the application of the PS correction method has great limitations.

[0006] (3) Atmospheric parameter correction method: First, the atmospheric refractive index is calculated based on the atmospheric temperature, humidity and air pressure data in the field of view using the corresponding atmospheric parameter correction model. The refractive index is then used to compensate for the influence of the atmosphere, thereby improving the monitoring accuracy. Experiments have found that this method can improve the monitoring accuracy by about 50%, and is currently a relatively effective method for correcting atmospheric disturbance errors.

[0007] However, the traditional atmospheric parameter correction method still has the following problems: on the one hand, it is necessary to assume that the atmosphere is uniform along the propagation path of the GB-SAR radar wave, that is, only the temporal changes of meteorological factors in the field of view are considered, and the spatial changes are ignored. However, given the height attributes of super-high-rise buildings, the spatial changes of meteorological factors in the observation scene are also very large, and the impact on monitoring accuracy cannot be ignored; on the other hand, since super-high-rise building monitoring mostly uses the high points of the structure as target points, the transmission distance of GB-SAR radar waves is relatively long. During the propagation process, the radar wave energy will attenuate as the propagation distance increases, thereby increasing the GB-SAR super-high-rise building horizontal time series displacement monitoring error. Summary of the invention

[0008] In order to solve or at least alleviate the above problems, the present invention proposes a method for monitoring the horizontal time-series displacement of super high-rise buildings based on inhomogeneous atmosphere correction, so as to reduce the influence of inhomogeneous meteorological factors and radar wave energy attenuation on the accuracy of GB-SAR monitoring of the horizontal time-series displacement of super high-rise buildings, and improve the accuracy of monitoring the horizontal time-series displacement of super high-rise buildings.

[0009] To achieve the above object, the present invention provides the following solutions:

[0010] A method for monitoring horizontal time-series displacement of super high-rise buildings based on inhomogeneous atmosphere correction, comprising:

[0011] Deploy the GB-SAR instrument on one side of the super high-rise building to be monitored, perform projection calculation based on the geometric characteristics of the super high-rise building and the relative position relationship of the GB-SAR instrument deployment, and construct a mathematical expression for the horizontal displacement of the super high-rise building after optimization and correction of inhomogeneous atmospheric parameters;

[0012] Based on the fitting method, the monitoring error correction factor caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases is obtained;

[0013] Based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number, a non-homogeneous atmospheric parameter correction optimization model taking into account the radar wave energy attenuation is constructed;

[0014] Obtain the radial displacement of super high-rise buildings monitored by the GB-SAR instrument and the radial monitoring distance from the GB-SAR instrument to the target monitoring point;

[0015] Calculate the atmospheric refractive index at the target height according to the target height of the target monitoring point;

[0016] Calculate the radial displacement and radial monitoring distance after optimization correction of inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance;

[0017] The radial displacement and radial monitoring distance corrected by the optimization of the heterogeneous atmospheric parameters are substituted into the heterogeneous atmospheric parameter correction optimization model to calculate the horizontal time-series displacement of the super high-rise building.

[0018] Optionally, the projection calculation is performed based on the relative position relationship between the geometric characteristics of the super high-rise building and the deployment of the GB-SAR instrument to obtain a mathematical expression for the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters, including:

[0019] Projection calculation is performed based on the relative position relationship between the geometric characteristics of super high-rise buildings and the deployment of GB-SAR instruments, and a mathematical expression for the horizontal displacement of super high-rise buildings after optimization and correction of inhomogeneous atmospheric parameters is constructed. Wherein, Δs' is the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; S is the horizontal distance between the deployment position of the GB-SAR instrument and the super high-rise building to be monitored; Δd' is the radial displacement after optimization and correction of the inhomogeneous atmospheric parameters; and D' is the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters.

[0020] Optionally, constructing a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number specifically includes:

[0021] Based on the mathematical expression of horizontal displacement of super high-rise buildings and the monitoring error correction number v(r), construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation Where ΔS is the horizontal time-series displacement of the super high-rise building.

[0022] Optionally, calculating the atmospheric refractive index at the target height according to the target height of the target monitoring point specifically includes:

[0023] According to the target height h of the target monitoring point, the atmospheric refractive index N(h) at the target height is calculated using the formula N(h)=N0exp(-h / h0); wherein N0 is the near-surface bottom value of the refractive index, and h0 is the near-surface elevation value;

[0024] According to the atmospheric refractive index N(h) at the target height, the formula N(h)=(n h -1)×10 6 Calculate the atmospheric refractive index n at the target height h .

[0025] Optionally, the calculating the radial displacement and radial monitoring distance after optimization correction of inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance specifically includes:

[0026] According to the atmospheric refractive index n at the target height h and the radial displacement Δd, using the formula Δd'=Δd / n h Calculate the radial displacement Δd' after optimization correction of inhomogeneous atmospheric parameters;

[0027] According to the atmospheric refractive index n at the target height h , using the formula Calculate the radial monitoring distance D' after optimization correction of inhomogeneous atmospheric parameters; wherein Δh is the height change value of each refraction unit divided within the radial monitoring distance; n0 is the atmospheric refractive index of the near-surface refraction unit; β0 is the radar wave incident angle of the near-surface refraction unit; n i is the atmospheric refractive index of the i-th refractive unit.

[0028] A super high-rise building horizontal time-series displacement monitoring system based on inhomogeneous atmosphere correction, comprising:

[0029] The module for constructing the mathematical expression of horizontal displacement is used to deploy the GB-SAR instrument on one side of the super-high-rise building to be monitored, perform projection calculation based on the geometric characteristics of the super-high-rise building and the relative position relationship of the GB-SAR instrument deployment, and construct the mathematical expression of the horizontal displacement of the super-high-rise building after optimization and correction of the inhomogeneous atmospheric parameters;

[0030] The correction number acquisition module is used to obtain the monitoring error correction number caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases based on the fitting method;

[0031] A non-homogeneous atmospheric parameter correction optimization model construction module, used to construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number;

[0032] The monitoring data acquisition module is used to obtain the radial displacement of the super high-rise building monitored by the GB-SAR instrument and the radial monitoring distance from the GB-SAR instrument to the target monitoring point;

[0033] A refractive index calculation module, used for calculating the atmospheric refractive index at the target height according to the target height of the target monitoring point;

[0034] An atmospheric refractive index corrected data calculation module, used for calculating the radial displacement and radial monitoring distance corrected by inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance;

[0035] The horizontal time-series displacement calculation module is used to substitute the radial displacement and radial monitoring distance corrected by the non-homogeneous atmospheric parameter optimization into the non-homogeneous atmospheric parameter correction optimization model to calculate the horizontal time-series displacement of the super high-rise building.

[0036] Optionally, the horizontal displacement mathematical expression construction module specifically includes:

[0037] The horizontal displacement mathematical expression acquisition unit is used to perform projection calculation based on the relative position relationship between the geometric characteristics of super high-rise buildings and the deployment of GB-SAR instruments, and to construct the mathematical expression of the horizontal displacement of super high-rise buildings after optimization and correction of inhomogeneous atmospheric parameters. Wherein, Δs' is the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; S is the horizontal distance between the deployment position of the GB-SAR instrument and the super high-rise building to be monitored; Δd' is the radial displacement after optimization and correction of the inhomogeneous atmospheric parameters; and D' is the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters.

[0038] Optionally, the heterogeneous atmospheric parameter correction optimization model construction module specifically includes:

[0039] A non-homogeneous atmospheric parameter correction optimization model construction unit is used to construct a super high-rise building horizontal displacement mathematical expression based on the above and the monitoring error correction number v(r), construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation Where ΔS is the horizontal time-series displacement of the super high-rise building.

[0040] Optionally, the refractive index calculation module specifically includes:

[0041] The atmospheric refractive index calculation unit is used to calculate the atmospheric refractive index N(h) at the target height h of the target monitoring point using the formula N(h)=N0exp(-h / h0); wherein N0 is the near-surface bottom value of the refractive index, and h0 is the near-surface elevation value;

[0042] The refractive index calculation unit is used to calculate the atmospheric refractive index N(h) at the target height using the formula N(h)=(n h -1)×10 6 Calculate the atmospheric refractive index n at the target height h .

[0043] Optionally, the atmospheric refractive index corrected data calculation module specifically includes:

[0044] A radial displacement calculation unit is used to calculate the atmospheric refractive index n at the target height. h and the radial displacement Δd, using the formula Δd'=Δd / n h Calculate the radial displacement Δd' after optimization correction of inhomogeneous atmospheric parameters;

[0045] The radial monitoring distance calculation unit is used to calculate the atmospheric refractive index n at the target height. h , using the formula Calculate the radial monitoring distance D' after optimization correction of inhomogeneous atmospheric parameters; wherein Δh is the height change value of each refraction unit divided within the radial monitoring distance; n0 is the atmospheric refractive index of the near-surface refraction unit; β0 is the radar wave incident angle of the near-surface refraction unit; n i is the atmospheric refractive index of the i-th refractive unit.

[0046] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0047] The present invention provides a method and system for monitoring the horizontal time-series displacement of a super-high-rise building based on inhomogeneous atmosphere correction. The method comprises: deploying a GB-SAR instrument on one side of a super-high-rise building to be monitored, performing projection calculation based on the relative position relationship between the geometric characteristics of the super-high-rise building and the deployment of the GB-SAR instrument, and constructing a mathematical expression for the horizontal displacement of the super-high-rise building after optimization and correction of inhomogeneous atmosphere parameters; obtaining a monitoring error correction number caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases based on a fitting method; constructing a linear time-series displacement correction number based on the mathematical expression for the horizontal displacement of the super-high-rise building and the monitoring error correction number, which takes into account the energy attenuation of radar waves. The invention provides an optimization model for correction of non-homogeneous atmospheric parameters with quantitative attenuation; obtains the radial displacement of super-high-rise buildings monitored by GB-SAR instruments and the radial monitoring distance from the GB-SAR instrument to the target monitoring point; calculates the atmospheric refractive index at the target height according to the target height of the target monitoring point; calculates the radial displacement and radial monitoring distance after optimization correction of non-homogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance; substitutes the radial displacement and radial monitoring distance after optimization correction of non-homogeneous atmospheric parameters into the optimization model for correction of non-homogeneous atmospheric parameters, and calculates the horizontal time-series displacement of super-high-rise buildings. The method of the invention can reduce the influence of non-homogeneous meteorological factors and radar wave energy attenuation on the monitoring accuracy of horizontal time-series displacement of super-high-rise buildings by GB-SAR, and improve the monitoring accuracy of horizontal time-series displacement of super-high-rise buildings. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0049] Figure 1 This is a flow chart of a method for monitoring horizontal time-series displacement of a super high-rise building based on inhomogeneous atmosphere correction according to the present invention;

[0050] Figure 2 A schematic diagram of radial displacement correction for heterogeneous meteorological factors provided by an embodiment of the present invention;

[0051] Figure 3 A schematic diagram of horizontal displacement projection calculation provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0052] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0053] The purpose of the present invention is to provide a method for monitoring the horizontal time-series displacement of super high-rise buildings based on heterogeneous atmosphere correction, so as to reduce the influence of heterogeneous meteorological factors and radar wave energy attenuation on the accuracy of GB-SAR super high-rise building horizontal time-series displacement monitoring, and improve the accuracy of super high-rise building horizontal time-series displacement monitoring.

[0054] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0055] Figure 1 This is a flow chart of a method for monitoring horizontal time-series displacement of super high-rise buildings based on heterogeneous atmosphere correction according to the present invention. Figure 1 The present invention provides a method for monitoring horizontal time-series displacement of super high-rise buildings based on heterogeneous atmosphere correction, which specifically includes:

[0056] Step 1: Deploy the GB-SAR instrument on one side of the super high-rise building to be monitored, perform projection calculation based on the geometric characteristics of the super high-rise building and the relative position relationship of the GB-SAR instrument deployment, and construct a mathematical expression for the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters.

[0057] Figure 2 Schematic diagram of radial displacement correction for heterogeneous meteorological factors provided by an embodiment of the present invention. Figure 2 As shown, the present invention monitors the horizontal time-series displacement of a super high-rise building by deploying a GB-SAR instrument on one side of the super high-rise building. Assume that the height of the super high-rise building is H; point Q is the point with the strongest reflection in a certain distance resolution unit, and it is used as the target monitoring point. The radial monitoring distance from the GB-SAR instrument to point Q is D, and the actual radial monitoring distance (i.e., the radial monitoring distance corrected by the optimization of inhomogeneous atmospheric parameters) is D'. If the meteorological factors on the propagation path of the radar wave are uniformly distributed, and the refractive index is n, the relationship between D and D' is:

[0058] D'=D / n(1)

[0059] When considering the disturbance of heterogeneous meteorological factors, the method of the present invention assumes that the meteorological factors within the height range of Δh (Δh=D / m, m is a positive integer greater than or equal to 1) are uniformly distributed, but the meteorological factors change continuously in units of height Δh to form a refraction unit. Figure 2As shown, the refractive index values ​​of different refractive units are represented by n0, n1, ... n m ; d0,d1,…d m is the GB-SAR radial monitoring distance in different refraction units; β0,β1,…β m are the incident angle and refraction angle of the radar wave at different refraction units. For ease of analysis, Figure 2 The medium refraction unit is an amplification of the actual situation. When Δh→0, the refractive index is close to the real change state. Based on the continuous and random characteristics of the atmospheric medium, the atmospheric refractive index N(h) decreases with increasing altitude. The specific change process can be expressed as:

[0060] N(h)=N0 exp(-h / h0) (2)

[0061] Where N0 is the near-surface bottom value of the refractive index, which can be calculated by monitoring the surface temperature in degrees Celsius, atmospheric pressure and water vapor pressure; h0 is the near-surface elevation value, and h is the target height.

[0062] The relationship between the refractive index N and the refractive index n is as follows: (3), so the atmospheric refractive index n of each refractive unit can be obtained by the atmospheric refractive index: m (m=1,2,3…):

[0063] N=(n-1)×10 6 (3)

[0064] According to the law of refraction, the radial displacement from the GB-SAR instrument to point Q in each refraction unit is first obtained as:

[0065]

[0066] The radial displacement from the GB-SAR instrument to point Q in each refraction unit after optimization and correction of inhomogeneous atmospheric parameters is:

[0067]

[0068] where d i is the radial displacement within the i-th refractive unit; d i ' is the radial displacement corrected by the optimization of inhomogeneous atmospheric parameters in the i-th refraction unit; n i is the atmospheric refractive index of the i-th refractive unit.

[0069] Therefore, the radial monitoring distance from the GB-SAR instrument to point Q after optimization and correction of inhomogeneous atmospheric parameters is:

[0070]

[0071] In order to obtain the GB-SAR horizontal displacement corrected by the inhomogeneous atmospheric refractivity, it is necessary to perform projection calculation based on the relative position relationship between the geometric characteristics of the super-high-rise buildings to be monitored and the deployment of the GB-SAR instrument. Figure 3 Schematic diagram of horizontal displacement projection calculation provided by an embodiment of the present invention. Figure 3 As shown, point G is the deployment position of the GB-SAR instrument. It is assumed that point Q is the point with the strongest reflection in a certain distance resolution unit of the super high-rise building (i.e., the monitoring target point), and the super high-rise building is horizontally offset to point Q1 under the action of dynamic loading such as wind load. The monitoring distances from the GB-SAR instrument to points Q and Q1 are defined as D and D1, respectively. Considering that the vertical displacement of the super high-rise building structure is very small, it is assumed that point Q1 is still within the original distance resolution. In order to obtain the projection theoretical calculation model of the horizontal displacement Δs of the super high-rise building and the radial displacement Δd monitored by GB-SAR, a parallel line GQ1 is made through point Q. At this time there are:

[0072]

[0073] Where S is the horizontal distance between the GB-SAR instrument deployment location and the monitored super-high-rise building.

[0074] And from Figure 3 From the geometric relationship, we can know that:

[0075]

[0076] Substituting formula (8) into formula (7) yields:

[0077] (S-Δs) 2 +D 2 -S 2 =(D-Δd) 2 (9)

[0078] Simplifying formula (9) yields:

[0079]

[0080] The Δd and D monitored by GB-SAR are corrected for atmospheric refractivity and converted to Δd' and D' respectively, where Δd' = Δd / n m , D' can be obtained according to formula (6). Therefore, considering the heterogeneous atmospheric disturbance, the mathematical expression of the horizontal displacement Δs' of the super high-rise building after optimization and correction of the heterogeneous atmospheric parameters is:

[0081]

[0082] Wherein, Δs' is the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; S is the horizontal distance between the deployment position of the GB-SAR instrument and the super high-rise building to be monitored; Δd' is the radial displacement after optimization and correction of the inhomogeneous atmospheric parameters; and D' is the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters.

[0083] Step 2: Based on the fitting method, the monitoring error correction factor caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases is obtained.

[0084] Considering that the energy of GB-SAR radar waves will decay as the propagation distance increases, thus affecting the GB-SAR super-high-rise building horizontal time-series displacement monitoring accuracy, the present invention obtains the correction number v(r) (called the monitoring error correction number) of the horizontal time-series displacement monitoring error caused by the attenuation of GB-SAR radar wave energy as the transmission distance increases based on the fitting method, and constructs a non-homogeneous atmospheric parameter correction optimization model taking into account the radar wave energy attenuation based on this. The fitting method used to obtain the monitoring error correction number v(r) caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases can be the least square method or the quadratic curve method. Of course, when obtaining the radar wave transmission distance correction number v(r), in order to avoid the low accuracy of the correction number v(r) based on the least square method or the quadratic curve method, a variety of fitting methods such as exponential, linear, logarithmic, polynomial, etc. can also be used instead.

[0085] Step 3: Based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number, a non-homogeneous atmospheric parameter correction optimization model taking into account the radar wave energy attenuation is constructed.

[0086] Based on the mathematical expression (11) of the horizontal displacement of the super high-rise building and the monitoring error correction number v(r), the present invention constructs a non-homogeneous atmospheric parameter correction optimization model taking into account the radar wave energy attenuation as follows:

[0087]

[0088] Where ΔS is the horizontal time-series displacement of the super high-rise building.

[0089] In addition to the influence of heterogeneous atmospheric disturbances, the energy attenuation of GB-SAR radar waves with increasing transmission distance will also affect the GB-SAR horizontal displacement monitoring results. Therefore, formula (12) adds the correction factor v(r) for the influence of GB-SAR radar wave energy attenuation with increasing transmission distance on horizontal displacement monitoring on the basis of formula (11), thereby constructing a heterogeneous atmospheric parameter correction optimization model that takes into account radar wave energy attenuation.

[0090] Step 4: Obtain the radial displacement of the super-high-rise building monitored by the GB-SAR instrument and the radial monitoring distance from the GB-SAR instrument to the target monitoring point.

[0091] In the process of horizontal time-series displacement monitoring of super high-rise buildings, the radial displacement Δd of the super high-rise building monitored by the GB-SAR instrument at a certain moment and the radial monitoring distance D from the GB-SAR instrument to the target monitoring point Q are obtained.

[0092] Step 5: Calculate the atmospheric refractive index at the target height according to the target height of the target monitoring point.

[0093] The calculation of the atmospheric refractive index at the target height can refer to formulas (2) and (3). The specific calculation process is as follows.

[0094] According to the target height h of the target monitoring point, the atmospheric refractive index N(h) at the target height is calculated using formula (13):

[0095] N(h)=N0 exp(-h / h0) (13)

[0096] Where N0 is the refractive index near the bottom of the surface, and h0 is the near-surface elevation.

[0097] Specifically,

[0098]

[0099] e=e'-0.000662(t-t')(1+0.001146t')P (15)

[0100] Where, T represents Kelvin temperature, T = 273.15 + t, t is Celsius temperature; t' is atmospheric humidity temperature; P is atmospheric pressure, unit is mmHg; e represents water vapor pressure, unit is mmHg; e' is standard water vapor pressure saturation, The unit is mmHg. These are the data that can be obtained through instrument monitoring at the height h0.

[0101] According to the corresponding relationship between the target height h of the target monitoring point, h0 and Δh, the refraction unit i corresponding to the target height h can be obtained, that is:

[0102] h=h0+i×Δh (16)

[0103] Therefore, the atmospheric refractive index n at the target altitude is h It is the atmospheric refractive index n of the i-th refractive unit corresponding to the target height h i .

[0104] According to the atmospheric refractive index N(h) at the target height, the atmospheric refractive index n at the target height can be calculated using formula (15): h :

[0105] N(h)=(n h-1)×10 6 (17)

[0106] Step 6: Calculate the radial displacement and radial monitoring distance after optimization correction of inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance.

[0107] Based on the above formula (5), the calculation formula of the radial displacement Δd' after optimization correction of heterogeneous atmospheric parameters can be obtained as follows:

[0108] Δd'=Δd / n h (18)

[0109] Based on the above formula (6), the calculation formula of the radial monitoring distance D' after optimization correction of heterogeneous atmospheric parameters can be obtained as follows:

[0110]

[0111] Wherein Δh is the height change value of each refraction unit divided within the radial monitoring distance; n0 is the atmospheric refractive index of the near-surface refraction unit; β0 is the radar wave incident angle of the near-surface refraction unit.

[0112] Step 7: Substitute the radial displacement and radial monitoring distance corrected by the inhomogeneous atmospheric parameter optimization into the inhomogeneous atmospheric parameter correction optimization model to calculate the horizontal time-series displacement of the super high-rise building.

[0113] When performing actual monitoring of the horizontal time-series displacement of a super high-rise building, it is only necessary to substitute the calculated radial displacement Δd' and radial monitoring distance D' after the non-homogeneous atmospheric parameter optimization correction into the non-homogeneous atmospheric parameter correction optimization model (12) to calculate the horizontal time-series displacement monitoring value ΔS of the super high-rise building. In actual monitoring, it is usually necessary to calculate multiple ΔS values ​​within a period of time to form a horizontal time-series displacement monitoring sequence of the super high-rise building to be monitored. According to the horizontal time-series displacement monitoring sequence, the super high-rise building structural damage identification and structural status performance assessment can be carried out, which is an important basis for the health status assessment of the super high-rise building.

[0114] Traditional atmospheric parameter correction methods mostly assume that the atmosphere in the propagation path of GB-SAR radar waves within the observation range is uniformly distributed, that is, only the influence of meteorological factors changing with time in the field of view is considered, and the uneven distribution in space is ignored. Based on the height attribute of super high-rise buildings, the present invention fully considers the characteristics of the inhomogeneous distribution of meteorological factors in the observation scene, and studies the change of atmospheric refractive index under the influence of inhomogeneous meteorological factors based on the relative position relationship between the geometric characteristics of the super high-rise buildings to be monitored and the GB-SAR instrument layout, and obtains the mathematical expression (11) of the horizontal displacement of GB-SAR super high-rise buildings corrected by the inhomogeneous atmospheric refractive index, thereby reducing the influence of inhomogeneous meteorological factors on the accuracy of GB-SAR super high-rise building horizontal time series displacement monitoring.

[0115] For the horizontal time series monitoring of super high-rise buildings, the high points of the structures are often used as the monitoring target points. The transmission distance of GB-SAR radar waves is relatively long. During the propagation process, the radar wave energy will decay as the propagation distance increases, which increases the GB-SAR horizontal time series displacement monitoring error. Therefore, the present invention obtains the monitoring error correction number v(r) caused by the energy decay of the radar wave with the transmission distance based on the least square method or the quadratic curve method, and proposes a non-homogeneous atmospheric parameter correction optimization model (12) that takes into account the radar wave energy decay, thereby reducing the influence of the radar wave energy decay on the monitoring accuracy caused by the long GB-SAR monitoring distance, and improving the credibility of GB-SAR super high-rise building horizontal time series displacement monitoring.

[0116] Based on the method provided by the present invention, the present invention also provides a super high-rise building horizontal time-series displacement monitoring system based on inhomogeneous atmosphere correction, the system comprising:

[0117] The module for constructing the mathematical expression of horizontal displacement is used to deploy the GB-SAR instrument on one side of the super-high-rise building to be monitored, perform projection calculation based on the relative position relationship between the geometric characteristics of the super-high-rise building and the deployment of the GB-SAR instrument, and construct the mathematical expression of the horizontal displacement of the super-high-rise building after optimization and correction of the inhomogeneous atmospheric parameters;

[0118] The correction number acquisition module is used to obtain the monitoring error correction number caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases based on the fitting method;

[0119] A non-homogeneous atmospheric parameter correction optimization model construction module, used to construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number;

[0120] The monitoring data acquisition module is used to obtain the radial displacement of the super high-rise building monitored by the GB-SAR instrument and the radial monitoring distance from the GB-SAR instrument to the target monitoring point;

[0121] A refractive index calculation module, used for calculating the atmospheric refractive index at the target height according to the target height of the target monitoring point;

[0122] An atmospheric refractive index corrected data calculation module, used for calculating the radial displacement and radial monitoring distance corrected by inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance;

[0123] The horizontal time-series displacement calculation module is used to substitute the radial displacement and radial monitoring distance corrected by the non-homogeneous atmospheric parameter optimization into the non-homogeneous atmospheric parameter correction optimization model to calculate the horizontal time-series displacement of the super high-rise building.

[0124] The horizontal displacement mathematical expression construction module specifically includes:

[0125] The horizontal displacement mathematical expression acquisition unit is used to perform projection calculation based on the relative position relationship between the geometric characteristics of super high-rise buildings and the deployment of GB-SAR instruments, and to construct the mathematical expression of the horizontal displacement of super high-rise buildings after optimization and correction of inhomogeneous atmospheric parameters. Wherein, Δs' is the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; S is the horizontal distance between the deployment position of the GB-SAR instrument and the super high-rise building to be monitored; Δd' is the radial displacement after optimization and correction of the inhomogeneous atmospheric parameters; and D' is the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters.

[0126] The heterogeneous atmospheric parameter correction optimization model construction module specifically includes:

[0127] A non-homogeneous atmospheric parameter correction optimization model construction unit is used to construct a super high-rise building horizontal displacement mathematical expression based on the above and the monitoring error correction number v(r), construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation Where ΔS is the horizontal time-series displacement of the super high-rise building.

[0128] The refractive index calculation module specifically includes:

[0129] The atmospheric refractive index calculation unit is used to calculate the atmospheric refractive index N(h) at the target height h of the target monitoring point using the formula N(h)=N0exp(-h / h0); wherein N0 is the near-surface bottom value of the refractive index, and h0 is the near-surface elevation value;

[0130] The refractive index calculation unit is used to calculate the atmospheric refractive index N(h) at the target height using the formula N(h)=(n h -1)×10 6 Calculate the atmospheric refractive index n at the target height h .

[0131] The atmospheric refractive index corrected data calculation module specifically includes:

[0132] A radial displacement calculation unit is used to calculate the atmospheric refractive index n at the target height. h and the radial displacement Δd, using the formula Δd'=Δd / n h Calculate the radial displacement Δd' after optimization correction of inhomogeneous atmospheric parameters;

[0133] The radial monitoring distance calculation unit is used to calculate the atmospheric refractive index n at the target height. h , using the formula Calculate the radial monitoring distance D' after optimization and correction of the inhomogeneous atmospheric parameters; wherein Δh is the height change value of each refraction unit divided within the radial monitoring distance; n0 is the atmospheric refractive index of the near-surface refraction unit; β0 is the radar wave incident angle of the near-surface refraction unit; n i is the atmospheric refractive index of the ith refractive unit.

[0134] The method and system of the present invention are based on the relative position relationship between the geometric characteristics of the super-high-rise building to be monitored and the deployment of the GB-SAR instrument, calculate the change of the atmospheric refractive index under the influence of the inhomogeneous meteorological factors, and obtain the mathematical expression of the horizontal displacement after the optimization and correction of the inhomogeneous atmospheric parameters; based on the fitting method such as the least square method or the quadratic curve method, the monitoring error correction number caused by the energy attenuation of the radar wave with the increase of the transmission distance is obtained, and then the inhomogeneous atmospheric parameter correction optimization model considering the radar wave energy attenuation is constructed. The inhomogeneous atmospheric parameter correction optimization model takes into account the GB-SAR radar wave energy attenuation and the disturbance of the inhomogeneous meteorological factors, so it can reduce the influence of the inhomogeneous meteorological factors and the radar wave energy attenuation on the GB-SAR super-high-rise building horizontal time-series displacement monitoring accuracy, and improve the GB-SAR monitoring atmospheric parameter correction accuracy.

[0135] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0136] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for monitoring horizontal time-series displacement of super high-rise buildings based on inhomogeneous atmosphere correction, characterized in that: include: Deploy the GB-SAR instrument on one side of the super high-rise building to be monitored, perform projection calculation based on the geometric characteristics of the super high-rise building and the relative position relationship of the GB-SAR instrument deployment, and construct a mathematical expression for the horizontal displacement of the super high-rise building after optimization and correction of inhomogeneous atmospheric parameters; The projection calculation is performed based on the relative position relationship between the geometric characteristics of the super high-rise building and the deployment of the GB-SAR instrument, and a mathematical expression for the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters is constructed, including: Projection calculation is performed based on the relative position relationship between the geometric characteristics of super high-rise buildings and the deployment of GB-SAR instruments, and a mathematical expression for the horizontal displacement of super high-rise buildings after optimization and correction of inhomogeneous atmospheric parameters is constructed. Where Δs' is the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; S is the horizontal distance between the deployment position of the GB-SAR instrument and the super high-rise building to be monitored; Δd' is the radial displacement after optimization and correction of the inhomogeneous atmospheric parameters; D' is the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters; Based on the fitting method, the monitoring error correction factor caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases is obtained; Based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number, a non-homogeneous atmospheric parameter correction optimization model taking into account the radar wave energy attenuation is constructed; Obtain the radial displacement of super high-rise buildings monitored by the GB-SAR instrument and the radial monitoring distance from the GB-SAR instrument to the target monitoring point; Calculate the atmospheric refractive index at the target height according to the target height of the target monitoring point; The step of calculating the atmospheric refractive index at the target height according to the target height of the target monitoring point specifically includes: According to the target height h of the target monitoring point, the atmospheric refractive index N(h) at the target height is calculated using the formula N(h)=N0exp(-h / h0); wherein N0 is the near-surface bottom value of the refractive index, and h0 is the near-surface elevation value; According to the atmospheric refractive index N(h) at the target height, the formula N(h)=(n h -1)×10 6 Calculate the atmospheric refractive index n at the target height h ; When considering the disturbance of non-homogeneous meteorological factors, it is assumed that the meteorological factors within the height range of Δh are uniformly distributed, but the meteorological factors change continuously in units of height Δh, forming a refractive unit; when Δh→0, the refractive index is close to the real change state; According to the corresponding relationship between the target height h of the target monitoring point and h0 and Δh, the refractive unit i corresponding to the target height h is known, that is: h = h0 + i × Δh; therefore, the atmospheric refractive index n at the target height h It is the atmospheric refractive index n of the i-th refractive unit corresponding to the target height h i ; Calculate the radial displacement and radial monitoring distance after optimization correction of inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance; The step of calculating the radial displacement and the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance specifically includes: According to the atmospheric refractive index n at the target height h and the radial displacement Δd, using the formula Δd'=Δd / n h Calculate the radial displacement Δd' after optimization correction of inhomogeneous atmospheric parameters; According to the atmospheric refractive index n at the target height h , using the formula Calculate the radial monitoring distance D' after optimization correction of inhomogeneous atmospheric parameters; wherein Δh is the height change value of each refraction unit divided within the radial monitoring distance; n0 is the atmospheric refractive index of the near-surface refraction unit; β0 is the radar wave incident angle of the near-surface refraction unit; n i is the atmospheric refractive index of the i-th refractive unit; The radial displacement and radial monitoring distance corrected by the optimization of the heterogeneous atmospheric parameters are substituted into the heterogeneous atmospheric parameter correction optimization model to calculate the horizontal time-series displacement of the super high-rise building.

2. The method according to claim 1, characterized in that The method of constructing a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number specifically includes: Based on the mathematical expression of horizontal displacement of super high-rise buildings and the monitoring error correction number v(r), construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation Where ΔS is the horizontal time-series displacement of the super high-rise building.

3. A super high-rise building horizontal time-series displacement monitoring system based on inhomogeneous atmosphere correction, characterized in that: include: The module for constructing the mathematical expression of horizontal displacement is used to deploy the GB-SAR instrument on one side of the super-high-rise building to be monitored, perform projection calculation based on the geometric characteristics of the super-high-rise building and the relative position relationship of the GB-SAR instrument deployment, and construct the mathematical expression of the horizontal displacement of the super-high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; The horizontal displacement mathematical expression building module specifically includes: The horizontal displacement mathematical expression acquisition unit is used to perform projection calculation based on the relative position relationship between the geometric characteristics of super high-rise buildings and the deployment of GB-SAR instruments, and to construct the mathematical expression of the horizontal displacement of super high-rise buildings after optimization and correction of inhomogeneous atmospheric parameters. Where Δs' is the horizontal displacement of the super high-rise building after optimization and correction of the inhomogeneous atmospheric parameters; S is the horizontal distance between the deployment position of the GB-SAR instrument and the super high-rise building to be monitored; Δd' is the radial displacement after optimization and correction of the inhomogeneous atmospheric parameters; D' is the radial monitoring distance after optimization and correction of the inhomogeneous atmospheric parameters; The correction number acquisition module is used to obtain the monitoring error correction number caused by the energy attenuation of GB-SAR radar waves as the transmission distance increases based on the fitting method; A non-homogeneous atmospheric parameter correction optimization model construction module, used to construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation based on the mathematical expression of the horizontal displacement of the super high-rise building and the monitoring error correction number; The monitoring data acquisition module is used to obtain the radial displacement of the super high-rise building monitored by the GB-SAR instrument and the radial monitoring distance from the GB-SAR instrument to the target monitoring point; A refractive index calculation module, used for calculating the atmospheric refractive index at the target height according to the target height of the target monitoring point; The refractive index calculation module specifically includes: The atmospheric refractive index calculation unit is used to calculate the atmospheric refractive index N(h) at the target height h of the target monitoring point using the formula N(h)=N0exp(-h / h0); wherein N0 is the near-surface bottom value of the refractive index, and h0 is the near-surface elevation value; The refractive index calculation unit is used to calculate the atmospheric refractive index N(h) at the target height using the formula N(h)=(n h -1)×10 6 Calculate the atmospheric refractive index n at the target height h ; When considering the disturbance of non-homogeneous meteorological factors, it is assumed that the meteorological factors within the height range of Δh are uniformly distributed, but the meteorological factors change continuously in units of height Δh, forming a refractive unit; when Δh→0, the refractive index is close to the real change state; According to the corresponding relationship between the target height h of the target monitoring point and h0 and Δh, the refractive unit i corresponding to the target height h is known, that is: h = h0 + i × Δh; therefore, the atmospheric refractive index n at the target height h It is the atmospheric refractive index n of the i-th refractive unit corresponding to the target height h i ; An atmospheric refractive index corrected data calculation module, used for calculating the radial displacement and radial monitoring distance corrected by inhomogeneous atmospheric parameters according to the atmospheric refractive index at the target height, the radial displacement and the radial monitoring distance; The atmospheric refractive index corrected data calculation module specifically includes: A radial displacement calculation unit is used to calculate the atmospheric refractive index n at the target height. h and the radial displacement Δd, using the formula Δd'=Δd / n h Calculate the radial displacement Δd' after optimization correction of inhomogeneous atmospheric parameters; The radial monitoring distance calculation unit is used to calculate the atmospheric refractive index n at the target height. h , using the formula Calculate the radial monitoring distance D' after optimization correction of inhomogeneous atmospheric parameters; wherein Δh is the height change value of each refraction unit divided within the radial monitoring distance; n0 is the atmospheric refractive index of the near-surface refraction unit; β0 is the radar wave incident angle of the near-surface refraction unit; n i is the atmospheric refractive index of the i-th refractive unit; The horizontal time-series displacement calculation module is used to substitute the radial displacement and radial monitoring distance corrected by the non-homogeneous atmospheric parameter optimization into the non-homogeneous atmospheric parameter correction optimization model to calculate the horizontal time-series displacement of the super high-rise building.

4. The system according to claim 3, characterized in that The heterogeneous atmospheric parameter correction optimization model construction module specifically includes: A non-homogeneous atmospheric parameter correction optimization model construction unit is used to construct a super high-rise building horizontal displacement mathematical expression based on the above and the monitoring error correction number v(r), construct a non-homogeneous atmospheric parameter correction optimization model taking into account radar wave energy attenuation Where ΔS is the horizontal time-series displacement of the super high-rise building.

Images