Shield excavation face deformation monitoring method and system for undercrossing airport runway

Abstract

The invention provides a shield excavation face deformation monitoring method and system for undercrossing an airport runway, and the method comprises the steps: generating a pressure load response characteristic based on the pressure fluctuation time sequence data of a cutter head hydraulic cylinder in a shield tunneling machine in combination with dynamic load distribution data and static load distribution data; the deformation track of the shield excavation face and the displacement change of a bottom soil body undercrossing the airport runway are obtained, point cloud registration processing is carried out, registered data are obtained, and two-dimensional deformation monitoring data of the shield excavation face and the bottom soil body in the horizontal direction and the vertical direction are extracted; based on the pressure load response characteristics, the two-dimensional deformation monitoring data and the surface settlement monitoring data of the undercrossing airport runway, obtaining a deformation trend prediction result of the shield excavation face; according to the method, high-precision real-time monitoring and trend prediction of collaborative deformation of the shield excavation face and the runway structure in the underpass airport runway scene are achieved, and runway operation safety and shield construction accurate control are guaranteed.

Claims

1. A method for monitoring deformation of shield excavation surface for underpass of airport runway, characterized in that: include: Based on the pressure fluctuation time series data of the cutterhead hydraulic cylinder in the shield machine, combined with the dynamic load distribution data of the airport runway surface during the aircraft take-off and landing stage and the static load distribution data during the shutdown stage, the pressure load response characteristics associated with the shield machine's propulsion direction are generated; Obtaining the deformation trajectory of the shield excavation surface and the displacement change of the bottom soil under the airport runway, performing point cloud registration processing on the deformation trajectory and the displacement change to obtain registered data; Extracting two-dimensional deformation monitoring data of the shield excavation surface and the bottom soil in the horizontal direction and the vertical direction from the registered data; The pressure load response characteristics, the two-dimensional deformation monitoring data and the surface settlement monitoring data of the underpass airport runway are correlated and analyzed to obtain the deformation trend prediction result of the shield excavation face.

2. The method according to claim 1, characterized in that Obtain the deformation trajectory of the shield excavation surface and the displacement change of the bottom soil under the airport runway, perform point cloud registration processing on the deformation trajectory and the displacement change, and obtain the registered data, including: Generate deformation trajectory data based on the deformation trajectory of the shield excavation surface scanned by the laser scanning device in a preset annular scanning path; Based on the bottom soil data corresponding to a plurality of different times scanned by the laser scanning device in a preset period, the coordinate offset of the bottom soil in the preset scanning period is identified, and the displacement change data is generated based on the coordinate offset; According to the advancement distance of the shield machine, the scanning frequency of the laser scanning device is adjusted to generate synchronous deformation trajectory data and synchronous displacement change data; Taking the axis of the shield machine as a reference, a three-dimensional coordinate system is constructed, and the synchronous deformation trajectory data and the synchronous displacement change data are mapped to the three-dimensional coordinate system to generate a deformation trajectory point cloud and a displacement change point cloud; On the inner wall surface of the shield machine cutout ring, high reflectivity marking points are set to dynamically compensate the deformation trajectory point cloud and the displacement change point cloud to generate a compensated deformation trajectory point cloud and a compensated displacement change point cloud; The position offset of the compensated deformation trajectory point cloud in the horizontal propulsion direction of the shield machine and the height change of the compensated displacement change point cloud in the vertical direction are calculated to generate registered data.

3. The method according to claim 2, characterized in that On the inner wall surface of the shield machine cut ring, high reflectivity marking points are set to dynamically compensate the deformation trajectory point cloud and the displacement change point cloud, and generate the compensated deformation trajectory point cloud and the compensated displacement change point cloud, including: On the inner wall surface of the shield machine cutout ring, a plurality of high reflectivity marking points are arranged, wherein the spatial coordinates of each high reflectivity marking point are pre-calibrated in the three-dimensional coordinate system; Extracting the real-time coordinates of each high-reflectivity marking point from the deformation trajectory point cloud and the displacement change point cloud, and calculating a coordinate offset vector between the real-time coordinates and the calibration space coordinates; Decomposing the coordinate offset vector into a longitudinal offset component and a lateral offset component, so as to perform reverse translation correction on the point cloud coordinates in the deformation trajectory point cloud and the displacement change point cloud, and obtain a corrected deformation trajectory point cloud and a corrected displacement change point cloud; According to the pitch angle data and yaw angle data of the shield machine attitude sensor, rotation compensation is performed on the corrected deformation trajectory point cloud and the corrected displacement change point cloud to obtain a rotation-compensated deformation trajectory point cloud and a rotation-compensated displacement change point cloud; According to the height change of the bottom soil in the rotationally compensated displacement change point cloud, the position offset of the rotationally compensated deformation trajectory point cloud is reversely corrected to generate a compensated deformation trajectory point cloud. According to the position offset of the rotationally compensated deformation trajectory point cloud, the height change of the displacement change point cloud is reversely corrected to generate a compensated displacement change point cloud.

4. The method according to claim 2, characterized in that: Calculating the position offset of the compensated deformation trajectory point cloud in the horizontal advancement direction of the shield machine and the height change of the compensated displacement change point cloud in the vertical direction to generate registered data, including: Along the horizontal advancement direction of the shield machine, the compensated deformation trajectory point cloud is divided into a plurality of horizontal monitoring units, and the average offset of the horizontal axis coordinates of all points in the compensated deformation trajectory point cloud in the horizontal monitoring unit relative to the initial reference horizontal position is counted to generate the position offset in the horizontal advancement direction; Dividing the compensated displacement change point cloud into a plurality of vertical monitoring units along the vertical advancing direction of the shield machine, counting the average height changes of the longitudinal axis coordinates of all points in the compensated displacement change point cloud in the vertical monitoring unit relative to the initial reference vertical position, and generating the height change in the vertical direction; According to the spatial position coordinates of the horizontal monitoring unit and the vertical monitoring unit, the corresponding position offset and height change at the same timestamp are mapped to generate a bidirectional displacement data set; A first ratio of the position offset of each horizontal monitoring unit in the bidirectional displacement data set to the height change of the corresponding vertical monitoring unit is calculated, and when the first ratio exceeds a preset threshold, the position offset is corrected to generate registered data.

5. The method according to claim 1, characterized in that Extracting two-dimensional deformation monitoring data of the shield excavation surface and the bottom soil in the horizontal direction and the vertical direction from the registered data, including: The horizontal advancement direction displacement data set corresponding to the shield excavation surface in the registered data is divided into a plurality of continuous grid units, and the maximum positive offset and the maximum negative offset of the horizontal coordinates of all data points in each grid unit relative to the original horizontal position are counted to generate a horizontal deformation monitoring vector; The vertical displacement data set corresponding to the soil at the bottom of the runway in the registered data is divided into multiple monitoring areas, and the maximum settlement and maximum uplift of the vertical coordinates of all data points in each monitoring area relative to the original height position are counted to generate a vertical deformation monitoring vector; According to the spatial position coordinates of the grid unit and the monitoring area, the horizontal deformation monitoring vector and the vertical deformation monitoring vector at the same timestamp are associated and bound to generate two-dimensional deformation monitoring data to be verified; Based on the mechanical properties of the rigid pavement structure underpassing the airport runway, the positive offset of the horizontal deformation monitoring vector and the settlement of the vertical deformation monitoring vector in the two-dimensional deformation monitoring data are proportionally verified to generate the two-dimensional deformation monitoring data.

6. The method according to claim 1, characterized in that The pressure load response characteristics, the two-dimensional deformation monitoring data and the surface settlement monitoring data under the airport runway are correlated and analyzed to obtain the deformation trend prediction results of the shield excavation face, including: Dividing the pressure load response characteristics into a plurality of settlement analysis segments, extracting the peak characteristics and period characteristics of the pressure fluctuation of the cutter head hydraulic cylinder in the settlement analysis segments, and generating a pressure load characteristic vector, wherein the settlement analysis segments correspond to the grid cells in the two-dimensional deformation monitoring data; Divide the surface settlement monitoring data under the airport runway into a plurality of settlement analysis segments, extract the mean and change trend of the surface settlement rate, and generate a surface settlement feature vector, wherein the settlement analysis segment corresponds to the monitoring area in the two-dimensional deformation monitoring data; Establishing a spatiotemporal correspondence between the pressure load characteristic vector, the horizontal deformation monitoring vector and the vertical deformation monitoring vector in the two-dimensional deformation monitoring data, and the surface settlement characteristic vector; Based on the time-space correspondence, construct a joint analysis data set; Generate a linkage status mark according to the joint analysis data set, so as to construct a deformation trend spatial distribution map of the shield excavation face based on the linkage status mark; Based on the distribution characteristics and evolution laws of each area in the deformation trend spatial distribution map, a deformation trend prediction result of the shield excavation face is generated.

7. The method according to claim 6, characterized in that Generate a linkage state mark according to the joint analysis data set, and construct a deformation trend spatial distribution map of the shield excavation face based on the linkage state mark, including: Traversing each analysis unit in the joint analysis data set, calculating a second ratio of a change amplitude of a peak feature in the pressure load feature vector to a change amplitude of a maximum positive offset of the horizontal deformation monitoring vector, and assigning a first type of linkage mark to the corresponding analysis unit if the second ratio is within a first preset interval; Calculating a third ratio of the maximum settlement in the vertical deformation monitoring vector to the mean value in the surface settlement feature vector, and assigning a second type of linkage mark to the corresponding analysis unit if the third ratio is within a second preset interval; Marking the analysis unit assigned with the first type linkage mark and the second type linkage mark as a linkage unit of shield thrust and runway deformation, and mapping the geographical location coordinates of the linkage unit to a preset two-dimensional space grid, wherein the preset two-dimensional space grid is generated based on the coordinates of the three-dimensional coordinate system; According to the density distribution of linkage units in the preset two-dimensional space grid, the linkage units are clustered to obtain a high-voltage load strong correlation area, a deformation hysteresis area, and a runway deformation conduction area; Based on the spatial distribution relationship of the high-pressure load strong correlation area, the deformation hysteresis area and the runway deformation conduction area, the expansion direction and area change of each area are superimposed to generate a deformation trend spatial distribution map of the shield excavation face.

8. A shield excavation face deformation monitoring system for underpass airport runway, characterized in that: include: A generation module is used to generate pressure load response characteristics associated with the propulsion direction of the shield machine based on the pressure fluctuation time series data of the cutter head hydraulic cylinder in the shield machine, combined with the dynamic load distribution data of the airport runway surface during the aircraft take-off and landing phase and the static load distribution data during the shutdown phase; A registration module is used to obtain the deformation trajectory of the shield excavation surface and the displacement change of the bottom soil under the airport runway, and perform point cloud registration processing on the deformation trajectory and the displacement change to obtain registered data; An extraction module, used to extract two-dimensional deformation monitoring data of the shield excavation surface and the bottom soil in the horizontal direction and the vertical direction from the registered data; The analysis module is used to correlate and analyze the pressure load response characteristics, the two-dimensional deformation monitoring data, and the surface settlement monitoring data of the underpass airport runway to obtain the deformation trend prediction result of the shield excavation face.

9. A computing device, characterized in that It comprises a processing component and a storage component; the storage component stores one or more computer instructions; the one or more computer instructions are used to be called and executed by the processing component to implement a shield excavation face deformation monitoring method for crossing an airport runway as described in any one of claims 1 to 7.

10. A computer storage medium, characterized in that: A computer program is stored, and when the computer program is executed by a computer, a shield excavation face deformation monitoring method for underpassing an airport runway as described in any one of claims 1 to 7 is implemented.

Images