Next Article in Journal
Improved Fault Diagnosis in Hydraulic Systems with Gated Convolutional Autoencoder and Partially Simulated Data
Next Article in Special Issue
Railway Overhead Contact System Point Cloud Classification
Previous Article in Journal
Machine-Learning Classification of Soil Bulk Density in Salt Marsh Environments
Previous Article in Special Issue
Three-Dimensional Linear Restoration of a Tunnel Based on Measured Track and Uncontrolled Mobile Laser Scanning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Method for Tunnel Displacements Calculation Based on Mobile Tunnel Monitoring System

Beijing Advanced Innovation Center for Imaging Theory and Technology, Key Laboratory of 3D Information Acquisition and Application, Ministry of Education, College of Resource Environment and Tourism, Academy for Multidisciplinary Studies, Capital Normal University, Beijing 100048, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(13), 4407; https://doi.org/10.3390/s21134407
Submission received: 17 March 2021 / Revised: 25 May 2021 / Accepted: 23 June 2021 / Published: 27 June 2021
(This article belongs to the Special Issue Selected Papers from The Sixth National LiDAR Conference)

Abstract

:
Efficient, high-precision, and automatic measurement of tunnel structural changes is the key to ensuring the safe operation of subways. Conventional manual, static, and discrete measurements cannot meet the requirements of rapid and full-section detection in subway construction and operation. Mobile laser scanning technology is the primary method for tunnel detection. Herein, we propose a method to calculate shield tunnel displacements of a full cross-section tunnel. The point cloud data, obtained via a mobile tunnel deformation detection system, were fitted, projected, and interpolated to generate an orthophoto image. Combined with the cumulative characteristics of the tunnel gray gradient, the longitudinal ring seam of the tunnel was identified, while the Canny algorithm and Hough line detection algorithm identified the transverse seam. The symmetrical vertical foot method and cross-section superposition analysis were used to calculate the circumferential and radial displacements, respectively. The proposed displacement calculation method achieves automatic recognition of a ring seam, reduces human–computer interaction, and is fast, intelligent, and accurate. Furthermore, the description of the tunnel deformation location and deformation amount is more quantitative and specific. These results confirm the significance of shield tunnel displacement monitoring based on mobile monitoring systems in tunnel disease monitoring.

1. Introduction

In 2020, a total of 44 cities in China will put into operation 7715.31 km of urban rail transit lines. New lines and sections with a length of 1349.7 km have been put into operation in 26 cities, which presents a new record. With the rapid development of urban subways and high-speed railways, regular tunnel disease monitoring is essential to ensure tunnel safety.
The development of surveying and mapping technology is constantly modifying the operation methods of conventional measurement, such as the rapid development of three-dimensional (3D) laser scanning technology and vehicle mobile measurement systems, which play increasingly important roles in many fields with advantages such as high precision, speed mass, and non-contact measurement [1]. Mobile 3D laser scanning technology has been used in street view scanning [2], urban road extraction [3], road storage, building modeling, infrastructure management, and mobility evaluation [4,5,6], and has become an important and popular technology for modern surveying and mapping. Owing to characteristics such as high efficiency, good quality, and little influence from the surrounding environment, it is feasible to introduce mobile measurements into tunnels.
At present, the primary need for tunnel deformation monitoring is to obtain the vertical and horizontal displacements, cross-sectional and convergence deformations, cracks, and tunnel leakages to determine whether the tunnel has deformations, and subsequently, the deformation size and period. There are several reasons for tunnel displacement, including segment selection, shield machine propulsion speed, shield posture, and bolt hole diameter [7,8].
The mainstream tunnel deformation monitoring methods include total station and caliper measurement, with the measurement accuracy of the total station being a maximum of 1 mm. Based on the 3D coordinate non-contact measurement technology of free and automatic total stations, the corresponding mathematical adjustment model is established, and the average error of each 3D coordinate measurement point may reach a maximum of 0.8 mm [9,10,11]. However, owing to the large investment of surveyors and the complex environment within the tunnel, instrument operation is difficult and inefficient, making it difficult to achieve full section detection [12]. Targets for station-type 3D laser scanning need to be set, and an image mosaic is needed between adjacent stations. The scanning time for a single station is about 10 min, the data redundancy is large, and the working efficiency and post-processing is relatively low [13,14,15,16,17]. To overcome these disadvantages, mobile measurement equipment can be employed. Such equipment can quickly obtain high-quality 3D point clouds of the track and surrounding environment, mileage of the tunnel segment, accurate rail profile, and more information. It has flexible production operation coupled with high speed and image resolution, reliable data results, and comprehensive information [18]. The measurement speed of the mobile laser measurement system used in this study can reach a maximum of 5.4 km/h, which is 5–10 times higher than that of traditional measurement methods.
The current representative tunnel mobile measurement equipment is the GRP 5000 tunnel scanning measurement system developed by Amberg and Leica, known as SiTrack: One Mobile track scanning system [19,20]. Mobile tunnel inspection vehicles are also becoming popular in China, for example, the railway mobile measurement system (RMMS-01) developed by Wuhan University and the tunnel mobile measurement software and hardware integration system (TMMS) developed by Tongji University [21]. The TMMS for tunnel movement measurement developed by Tongji University and Wuhan University can fit the tunnel cross-section; however, it only supports the convergence diameter measurement in the clearance convergence measurement and cannot calculate the displacement. GRP 5000, SiTrack: One, and RMMS-01 are equipped with GNSS (Global Navigation Satellite System) positioning and inertial navigation systems, which can provide absolute coordinate measurement, but the hardware cost is high. The above four kinds of equipment are the mainstream mobile measurement systems in the market. However, some construction companies have customized the mobile measurement monitoring vehicles according to their own needs and develop control software only for their own use, for example, the tunnel mobile measurement system developed by the Shanghai Geotechnical Investigation and Design Engineering Consulting (Group) Co., Ltd. (Shanghai, China).
Tunnel deformation has always been a major problem affecting tunnel safety. Many scholars have studied tunnel deformation, such as the influence of ground building pile foundation on underground tunnel [22], installing support in tunnel soil layer to reduce tunnel deformation [23], or using convergence-confinement method to analyze tunnel deformation [24]. Some have also established mathematical models for the surrounding soil to analyze the settlement of the tunnel and the surface [25,26]. At this stage, many scholars have studied tunnel convergence measurements; however, our focus is on extracting the section center by tunnel section fitting and disease detection based on tunnel orthophoto images [27,28,29,30]. For example, the calculation of the tunnel convergence diameter based on point cloud data and the calculation of the shield tunnel convergence diameter by denoising, slicing, and ellipse fitting of the point cloud data. The calculation accuracy is approximately 2 mm [31,32]. Several studies have also explored the calculation of tunnel clearance, chord length, guide height, and pull-out value of the catenary based on mobile laser data and achieved a calculation accuracy of 3 mm [33]. An ultrasonic distance measuring instrument was installed in the same section of the shield tunnel to calculate the chord length of the tunnel, leading to a calculation accuracy of 2 mm [34]. According to the surrounding geological environment, some studies have established tunnel and surface building deformation prediction models, whose prediction accuracy can reach 75% [35,36]. The generation algorithm and quality evaluation of tunnel orthophoto images based on mobile laser systems have also been studied. This method uses orthophoto images for color classification and visualization of the deformation. Although this method intuitively reflects tunnel deformation, it cannot quantitatively reflect the deformation, which has certain limitations [37]. In engineering practice, manual inspection and caliper measurement are often used to check the displacement, hence the efficiency is low. At present, research on the calculation method and accuracy of shield tunnel displacements is still in the nascent stages. In this study, a mobile laser scanning system was used for data acquisition and an orthophoto image was used to achieve the automatic recognition of the circumferential and radial seam and the circumferential and radial displacements were calculated. Due to the limitation of hardware, the calculation method introduced in this paper is only applicable to the calculation of railway tunnels constructed by the shield method.
Next, we introduce the software and hardware of the mobile monitoring system. In Section 3, the generation method of orthophotos and the recognition method of circumferential gap are introduced as well as the method of calculation of radial and circumferential displacements. In Section 4, the precision of the calculation results is evaluated. In Section 5 and Section 6, future research prospects and suggestions, and the final research results of this paper are discussed.

2. Mobile Tunnel Monitoring System (MTDS)

The system mainly includes a scanner, mobile monitoring car body, and acquisition and processing software. The scanner employs three brands of scanners, namely, Faro, Leica, and Z + F, and supports the auxiliary mileage encoder. The equipment was driven by a motor having five gears with adjustable speed. The mobile tunnel monitoring system relies primarily on the scanners and mobile measuring vehicles for scanning. Further, a tunnel laser scanning acquisition and processing software was developed with a net native interface for the Z + F, Faro, and Leica scanners. By controlling the start and stop of the scanners and coordinating through the mobile monitoring equipment, laser scanning data can be quickly collected in the tunnel. This paper introduces some efficient algorithms, such as robust estimation least-squares, robust noise elimination, real-time signal processing algorithms, and surface computational geometry algorithms. Furthermore, we used a central processing unit (CPU) multithread parallel computing algorithm for data solution and point cloud rendering to ensure the overall efficiency and high precision of the software and achieve the integration of the tunnel point cloud data acquisition and post-processing. The system is illustrated in Figure 1.
Commonly used scanners in mobile monitoring systems are the Faro focus 3D 120, the Leica scanning P30/P40, and the Z + F PROFILER 9012. Their angle measurement accuracy are 19″, 8″, and 0.02°, respectively. Their ranging accuracy and scanning speed are ±2, ±1.2, and ±1 mm, and 97,600, 1,000,000, and 10,160,000 pixel/s, respectively.

3. Methods

In this study we introduce an orthophoto generation algorithm, using orthophoto through peak detection algorithm, Canny edge detection, and Hough line detection algorithm to identify the capping block and longitudinal ring seam. Thereafter a cross-section superposition method is used to calculate the circular displacement using the local symmetrical perpendicular foot method to calculate the radial displacement. The calculation flowchart is shown in Figure 2.

3.1. Orthophoto Generation

3.1.1. Point Cloud Fitting Circle

After denoising, the point cloud is sliced according to the mileage, with the forward direction of the measuring equipment as y. Thereafter, each slice can be approximated as a circle and regarded as a plane composed of X and Z coordinates. Denoising was followed by fitting the point cloud of each slice via the least-squares method, and the fitting results provided the center coordinates and the radius of the circle. The point cloud was projected according to the fitting results.
The definition equation of the circle is used to establish the model
d i s t ( P , A ) + d i s t ( P , B ) = D I S T
where P is a point on the circle and A and B are the foci. Three points on the circle were randomly selected to build the circle model, the distance between each point and the two foci as well as the difference between each point and the DIST was calculated. This point is in line with the model when the difference is less than a certain threshold. The model with the largest number of points is the best-circle model. According to the points that meet the conditions, the general equation of the circle is as follows,
x i 2 + y i 2 + A x + B y + C = 0 ,
Coefficient fitting was carried out according to the obtained coincidence points. For the circle fitting of the least-squares method, the optimal objective function of the error squared is:
E = i = 0 n ( x i 2 + y i 2 + Ax i + By i + C ) 2 ,
According to the principle of the least-squares method, the parameters A, B, and C should be minimized. Finally, after the fitting the values best-fit center coordinates ( X 0 , Y 0 ) and radius R were obtained. The results are illustrated in Figure 3.
X 0 = A 2 , Y 0 = B 2 , r = 1 2 A 2 + B 2 4 C ,

3.1.2. Point Cloud Projection According to the Tunnel Arc Length

After fitting the circle, each tunnel slice corresponded to a set of fitting circle parameters. Considering that the cloud coordinate of a point in the slice is ( X , Z ) , the center coordinate of the circle is ( X R , Z R ), and the radius of the fitting circle is R. In the next part of this paper, the tunnel slice is called tunnel section. The angle θ formed by the section points and the center coordinate of the circle uses the arctangent function as follows:
θ = arctan ( Y Y R X X R )
Then according to the arc length formula of the circle:
l = θ R
Therefore, each point in the slice cloud calculates the arc length l one by one, and each slice corresponds to the mileage m, thus the 3D coordinates are projected into plane coordinates (m, l). A point cloud expansion diagram is shown in Figure 4.

3.1.3. Building Grid Index

In the last part, we have projected the coordinates of the point cloud to the plane. By traversing these points, we can find the maximum value X m a x and minimum value X m i n of the abscissa and the maximum value Y m a x and minimum value Y m i n of the ordinate. The plane grid is established according to the X m a x , X m i n and Y m a x , and Y m i n . The grid width is set to 0.01 m. Thus, the values of the grid horizontal total number a, grid vertical total number b, and grid normalization parameter n (normalization of point cloud intensity) are as follows:
a = X m a x X m i n 0.01 b = Y m a x Y m i n 0.01 n = 1 I m a x I m i n
The coordinates of the point cloud projected by the arc length in the upper section are recorded as ( X i , Y i ) . At the same time, X i X m i n 0.01 and ( Y i Y m i n ) 0.01 are the x- and y-direction index value of the points, respectively. Through the establishment of the grid, every point in the point cloud is placed into the grid, and each point matches a corresponding grid index value.

3.1.4. Grid to Image

An image matrix was established. The data type in the matrix was 64-bit single channel data, and the number of rows and columns were consistent with the grid index. The purpose of this step was to assign gray value to the image matrix. There are three scenarios in the projection.
  • There is only one point in the grid index
    If there is a point in the grid index corresponding to the row and column number of the matrix, and the intensity value of the point is I, then the gray value of the point g = 255 × I n is stored in the row and column numbers corresponding to the matrix.
  • There are two or more points in the grid
    If there are two or more points in the grid, let the grid center coordinate be ( X 0 , Y 0 ) and the point coordinate be ( X i , Y i ) . The Euclidean distance between the grid center and point was calculated as follows:
    d = ( X 0 X i ) 2 + ( Y 0 Y i ) 2
    If there is a point such that d < 10 6 , then the gray value g = 255 × I n corresponding to that point is stored in the row and column numbers corresponding to the image matrix.
  • No point in grid corresponding to matrix
    If there are no points in the grid corresponding to the row and column numbers of the matrix, making d < 10 6 , then the inverse distance weighted average interpolation (IDW) method was used. The reason of using IDW algorithm is to assign gray value to the grid (that is, image pixels) to prevent the existence of black spots in the image. The points in the point cloud data are regarded as discrete points in the plane, assuming that there are N discrete points ( X 1 , Y 1 ) , ( X 2 , Y 2 ) , …, ( X n , Y n ) . The grid center is regarded as the predicted interpolation point Z ( X , Y ) . The Euclidean distance formula was used as the distance function formula:
    h i = ( X X i ) 2 + ( Y Y i ) 2
    where ( X , Y ) is the interpolation point coordinate, and ( X i , Y i ) is the coordinate of each discrete point. When the power parameter is determined, the weight function is used to calculate the weight of each discrete point.
    w i = h i p i = 1 n h i p
    where p is the power parameter, and n is the number of discrete points; in this study, p = 3. The function for calculating the interpolation point is as follows:
    Z ( X , Y ) = i = 1 n w i × Z ( X i , Y i )
    where Z ( X , Y ) is the intensity value of the interpolation point, and the gray value g of the point is stored in the row and column number corresponding to the matrix. The calculation method for g is as follows:
    g = 255 × Z ( X , Y ) n
The gray orthophoto image in the tunnel can be obtained by outputting the image matrix. The orthophoto effect of the tunnel image is shown in Figure 5.

3.2. Longitudinal Ring Seam Identification

3.2.1. Longitudinal Ring Seam Detection

In this study, the peak detection algorithm was used to detect the longitudinal ring seam, and the gray value of the orthophoto image was searched longitudinally. The tunnel seam is dark in the image (the gray value is small), and the seam between the rings is perpendicular to the X-axis. Using this feature, the automatic recognition of the longitudinal seam can be realized. First, the gray value difference between columns i + 1 and i − 1 column is calculated row by row according to the image matrix.
Δ G i = G i + 1 G i 1 ,
After comparing i + 1 and i − 1 columns, the occurrence times of Δ G i < 0 were counted as t i , and the t value of each column of the matrix was recorded and treated as an array. Many peaks will appear in this array. According to the principle that the gray value at the ring seam is small and the longitudinal circumferential seam is perpendicular to the X axis, the column corresponding to the peak value is the position of the longitudinal ring seam. To avoid the influence of cracks and other ancillary tunnel facilities, the data of 0.75 m on either side of the first determined ring joint were not processed, thus eliminating identification errors. The above steps were repeated, and after several iterations, all the ring joints were identified; however, there is a possibility of gross errors. Therefore, we used the root mean square error (RMSE) to eliminate gross errors. Let x 1 , x 2 , …, x n be the column index of the calculated longitudinal ring seams, and md be the arithmetic average of the interval column numbers between the longitudinal ring seams.
m d = x 1 + x 2 + + x n n
Z = i = 1 n ( x i m d ) 2 n ,
where Z is the root mean square of x 1 , x 2 , …, x n . If | x i m d | > 2 Z is satisfied, then if I is not the corresponding column number of the longitudinal ring seam, it will be eliminated. After the circular seam is identified, the image pixel position is inversely indexed to the point cloud section data such that the position of the ring seam in the point cloud can be determined. The results of the longitudinal ring seam identification are shown in Figure 6. Through annular seam detection, the orthophoto image can be stored in different rings, which is convenient for later transverse seam detection and tunnel displacement calculation.

3.2.2. Transverse Seam Inspection of Capping Block

1.
Canny image binarization
Each shield tunnel ring consists of six segments, one capping block (KP), two adjacent blocks (CP, BP), and three standard blocks (A1p, A2p, A3p). According to the design data of the tunnel, the segment and installation angles of the capping block can be determined. The block diagram of the tunnel segment is shown in Figure 7a, and the position of the tunnel capping block is shown in Figure 7b. According to the corresponding angle of the capping block in the tunnel, the point cloud data can be preprocessed, and the orthophoto image of the capping block can be generated to prevent the interference of other tunnel ancillary facilities.
In this study, the Canny operator was used for image binarization. At first, the Canny operator used the Gaussian filtering method to smoothen the pixels in the image neighborhood, and thereafter the pixels in different positions in the neighborhood were assigned different weights.
While smoothing the image, most of the overall gray distribution characteristics of the image are retained. The 3D Gaussian distribution function is expressed as follows:
G ( x , y ) = 1 2 π σ 2 e x 2 + y 2 2 σ 2 ,
where σ represents the standard deviation, and the squares of S x and S y represent the distance between other pixels and the center pixel in the neighborhood, respectively. Subsequently, the gradient operator (Sobel) was used to calculate the gradient amplitude and direction:
G ( x , y ) = S x 2 + S y 2
R ( x , y ) = arctan ( S x S y )
where S x and S y are the convolutions of the gray level and grayscale in the X and Y direction, of the original image, respectively, G ( x , y ) is the gradient amplitude, and R ( x , y ) is the gradient direction. The edge intensity and direction at each point were estimated by calculating the gradient amplitude and direction. According to the gradient direction, the gradient amplitude was non-maximally suppressed, and finally, the double threshold method was used to detect and remove certain noise and false edges. The edge detection result of the Canny operator is shown in Figure 8a.
2.
Hough line detection
Hough line detection transforms the line in the image space to a point in the parameter space and solves the detection problem through statistical characteristics. Specifically, if the pixels in an image form a straight line, then the curves corresponding to the pixel coordinate values (x, y) in the parameter space must intersect at a point, such that we only transform all the pixel coordinate values in the image into curves in the parameter space, and then detect the intersection of the curves in the parameter space to determine the straight line. The expression for the Hough parameter space is as follows:
r = x i cos θ + y i sin θ ,
Theoretically, a point corresponds to innumerable lines in the parameter space; however, in practical applications, the number of lines (that is, a limited number of directions) must be limited for calculation. Therefore, a straight line can be drawn by discretizing θ, calculating r according to the coordinates (x, y) of the point, counting the occurrence times of (R, θ), and selecting the two (R, θ) values that were repeated the most. In this study, Hough line detection was used to identify the two transverse ring seams of the capping block. According to the segment design data, the ring seam position information of the other four segments can be inferred. The results of the Hough line detection are shown in Figure 8b.

3.3. Tunnel Displacement Calculation

3.3.1. Tunnel Radial Displacement Calculation

The radial displacement of shield tunnel refers to the phenomenon that the inner walls of adjacent segments in the same ring are not on a curved surface after the tunnel segments are assembled. In this study, the local symmetrical vertical foot method is used to solve the tunnel radial displacement. According to the tunnel transverse joint detection results, the real point cloud coordinates of the capping block circumferential joints are obtained via reverse retrieval of the point cloud. The middle position of the longitudinal circumferential seam can better reflect the radial dislocation of the segment. Thus, the section in the middle of the transverse joint (the section in the middle of the tunnel segment) is selected to calculate the radial displacements. The joint positions of other segments are calculated from the segment design data. Five adjacent sections in the middle of the circumferential seam are selected for circle fitting using the least square method, and five center coordinates ( X 0 , Y 0 ) are obtained. Each section can be regarded as a circle in the plane Cartesian coordinate system. Taking the capping block joint as an example, we set the joint position point coordinate as A ( X A , Y A ) , make a straight line from the center of the circle to the capping block joint position and perform the same operation for the five sections. Then, the functional expression of the line is as follows:
x X A X 0 X A = y Y 0 Y A Y 0
Furthermore, we select the nearest point on the left side of point A as B and the point on the right side of point A as C. To prevent the influence of noise points, we do not require point B and point C. Therefore, we need to set a threshold. The distance difference between points B and C and joint position point A ( X A , Y A ) should be less than 10 cm. Points B and C lead the perpendicular foot to the straight line x X A X 0 X A = y Y 0 Y A Y 0 (the straight line from the center of the circle to the position of the joint obtained in the previous step), and the perpendicular foot point is marked as B and C , respectively. The distance between the two perpendicular feet is the displacement value between the two segments, as shown in Figure 9b. The same calculation method was used for the five adjacent sections and an average value was obtained. The same process is followed for the other transverse joints, for example, BP–A1P. The radial displacement diagram is shown in Figure 9a.

3.3.2. Circumferential Displacement Calculation

The circumferential displacement of the shield tunnel is the phenomenon of displacement between two adjacent rings in the forward direction of the tunnel after the tunnel segments are assembled. Section superposition method is used to calculate circumferential displacement. According to the longitudinal ring detection results, the position information (odometer value) of the circumferential seam is obtained. In the point cloud, the nearest 10 adjacent sections on the left and right sides of the longitudinal ring seam are retrieved. After eliminating noise points, the first section on the left side of the circumferential seam is selected for circle fitting to obtain the section center ( X 0 , Y 0 ) and radius R. Then, a section is found on the corresponding symmetrical position on the right side of the circumferential seam. The two sections are superposed according to the center of the fitting circle. The superposition result is shown in Figure 10. The distance from each point on the left and right section to the fitting center, respectively, is calculated. The distance is calculated as:
d = ( X 0 X i ) 2 + ( Y 0 Y i ) 2
The distance from each point to the center of the left section is subtracted from the distance from the corresponding point of the right section to the center of the circle. If the difference between the distance d and the fitting radius R is more than 5 mm, and the arc length of the continuous displacement with a difference exceeding 5 mm is more than 0.5 m, the starting and ending angles of the displacement, the arc length of the displacement, and the arithmetic mean value of each displacement segment are counted. In the previous step, 10 sections were taken from the left and right sections, and a total of 10 groups of data were calculated in the same way. The average value of 10 sets of data is output as the result of circumferential displacement. The schematic cross sectional superposition diagram is shown in Figure 10.

4. Method Validation

4.1. Data Sources

Data from Tianjin and Nanchang Metro systems were used in this verification. The Tianjin metro tunnel is a circular shield tunnel with a ring width and inner diameter of 1.5 and 5.5 m, respectively. The tunnel was put into operation for operation safety monitoring. The Nanchang metro tunnel is also a circular shield tunnel with ring width and inner diameter of 1.5 and 5.4 m, respectively. The tunnel has been put into operation. The tunnel in the test section was built many years ago. Due to the displacement of the shield tunnel segments caused by geological changes, there are circumferential and radial displacements. The circumferential displacement is calculated using the section superposition method, and the radial dislocation is calculated using the local symmetric foot method. The mobile monitoring vehicle travels at a constant speed in the tunnel. Simultaneously, the acquisition software is used to control the start and stop of the scanner. The same vehicle speed and the same type of scanner were used in the same tunnel section for repeated scanning. Leica P16, Faro Focus 3D 120 and Z + F 9012 scanners (as shown in the Figure 11) were used for scanning. Leica P16, Faro, and Z + F 9012 operation photos are shown in Figure 11a–c, respectively.

4.2. Circumferential Displacement Round-Trip Precision Verification

In this study, Nanchang metro data were used for circumferential displacement precision verification. Circumferential displacement is the displacement of two adjacent rings in the mileage direction and the calculation result represents the average value of the displacement of two adjacent rings. Circumferential displacement is calculated by section superposition method. The length of the tunnel used for the experiment was 50 m. The mobile tunnel measurement system with Leica P16 was used to carry out the round-trip measurement at a speed of 0.5 m/s, and the round-trip comparison of 35 groups of circumferential displacement was carried out. The maximum absolute deviation of the round-trip measurement was 2.9 mm, and the arithmetic mean deviation was 1.4 mm. The results are shown in Table 1.

4.3. Radial Displacement Round-Trip Precision Verification

The Tianjin test data were used to verify the repeatability precision of radial displacement. Radial displacement is the phenomenon that the inner walls of adjacent segments in the same ring are not on the same arc surface and stagger with each other. The shield tunnel wall consists of six segments, one top block (KP), two adjacent blocks (CP, BP) and three standard blocks (A1p, A2p, A3p). The radial displacement is calculated by the vertical foot method. The test length was approximately 200 m. Because there is cable shielding in the Tianjin tunnel, only KPBP, BP–A1P, and A1P–A2P were used to verify the repeatability precision of radial displacement. The A1P–A2P, BP–A1P, and KP–BP displacements are shown in Figure 12, Figure 13 and Figure 14, respectively. According to the comparison of round-trip data, the KP–BP round-trip maximum deviation was 3.0 mm with an average deviation of 1.2 mm, the BP–A1P round-trip maximum deviation was 3.0 mm with an average deviation of 1.4 mm, while the A1P–A2P round-trip maximum deviation was 3.0 mm with an average deviation of 1.3 mm. The round-trip measurement precision was within 3.0 mm.

4.4. Radial Displacement Comparative Precision Verification

In this study, Focus 3D 120 and Z + F 9012 scanners were used to scan the 50 m shield tunnel of Tianjin Metro. The purpose of this experiment was to use different brands of scanners to scan in the same tunnel and compare the difference. Using this data, 21 groups of radial dislocation and 35 groups of circumferential displacement were calculated. The circumferential displacement was verified by KP–BP and BP–A1P. The difference of circumferential displacement is shown in Figure 15, and that of radial displacement is shown in Figure 16 and Figure 17. According to the comparison of data from different brands of scanners, the maximum circumferential displacement difference was 3.0 mm with an average difference of 1.96 mm, the KP–BP maximum difference was 3.0 mm with an average difference of 1.12 mm, and the BP–A1P maximum difference was 3.0 mm with an average difference of 1.6 mm. The contrast difference precision was within 3.0 mm.

4.5. Radial Displacement Absolute Accuracy Verification

Tianjin data were used to verify the absolute accuracy of the radial displacement algorithm. Because the BP–A1P ring seam was close to the ground and easy to measure, the BP–A1P ring seam was used to verify the absolute accuracy. In this study, 10 groups of BP–A1P radial displacement were measured using vernier calipers and compared with the radial displacement calculated by the local symmetrical vertical foot method. The comparison results are shown in Table 2. The maximum radial displacement difference was 3.0 mm with an average difference of 1.98 mm.
This study used three methods to verify the precision of the algorithm, which are round-trip repeatability precision verification, comparison verification of different brands of scanners in the same tunnel section, and absolute accuracy verification of vernier calipers measurement value compared with the measurement value of this method. These results showed that the calculation accuracy of the circumferential and radial displacements achieved good results. Through the test comparison, the average value of repeatability difference of circumferential displacement was 1.4 mm, the average value of repeatability difference of three groups of radial displacement was 1.3 mm, the difference of circumferential displacement contrast using different scanners was 1.96 mm, and the average difference of radial displacement contrast between the two groups was 1.36 mm. The absolute accuracy of radial displacement verified by vernier calipers was within 3 mm. The maximum difference of all validation was less than 3 mm, which meets the Metro deformation monitoring precision requirements.

5. Discussion

In this paper, we proposed a displacement calculation method based on mobile laser scanning, which includes circumferential and radial displacement calculation. Tunnel displacement monitoring plays an important role in current tunnel clearance convergence measurement. At present, the traditional manual inspection method, vernier caliper measurements, total station measurements, and model predictions are mostly used for tunnel displacement monitoring. However, the data acquisition speed of these technologies is low and the information obtained is not comprehensive. Furthermore, they do not provide a quantitative displacement calculation; therefore, data acquisition of full section displacements cannot be achieved. In this paper, the mobile laser scanning method was used to achieve rapid tunnel data acquisition and generate high-definition tunnel orthophotos. The tunnel ring seam recognition algorithm greatly reduces the in-office data processing workload. In the calculation of circumferential displacement, the section superposition method was used to comprehensively evaluate the displacement of two adjacent rings and calculate and output the average displacement. When calculating the radial displacement, the local symmetrical vertical foot method was used to evaluate the displacement between different segments in the same ring. Due to the shielding of tunnel cables, pipelines, and other ancillary facilities, some radial displacement could not be accurately calculated. Because the mobile monitoring vehicle can only run on the track, there is no in-depth study on its applicability. This paper mainly studies the structural deformation calculation of shield tunnel, without in-depth study on the applicability of the platform. However, the calculation method mentioned in this paper is suitable for shield tunnels with various purposes, which is only limited by hardware equipment (such as road tunnel). In future, the platform can be improved to drive on the ground, and the applicability of the algorithm can be explored. Despite its many advantages, the method proposed in this paper also has some limitations, such as the need to obtain tunnel segment design data, image quality dependence, the tunnel point cloud quality and the tunnel design data have a greater dependence, some capping block contours are not sufficiently clear, and the ring seam recognition effect is not ideal. To overcome these limitations and further improve the recognition effect, a deep learning neural network could be used in future work for the capping block and circumferential seam identification.

6. Conclusions

The calculation method proposed in this paper enriches the existing content of tunnel convergence measurement. Combining the orthophoto generation algorithm, peak detection algorithm, image Canny edge detection, and Hough line detection algorithm, the automatic recognition of tunnel ring seams and tunnel capping blocks was achieved. The section superposition and local symmetrical vertical foot methods were used to calculate the circumferential displacement and radial displacement, respectively. The experimental results showed that the repeatability precision of the circumferential and radial displacements was within 3 mm and the results of the difference comparison of circumferential and radial displacement using different scanners and the absolute accuracy verification using vernier calipers were less than 3 mm. The mining method, shield method construction monitoring and operation line track, track bed, concrete structure, engineering structure in construction, tunnel vault sinking, and structure convergence deformation, with an allowable error of ±3 mm, met the accuracy requirements. At present, the calculation of tunnel displacement was always the key and difficult point in tunnel deformation detection. The section superposition method and the local symmetrical vertical foot method proposed in this paper can accurately and comprehensively evaluate the tunnel shape variables, which provides a guarantee for the construction safety and operational safety. The calculation accuracy of the method can achieve good results through comparison. Although the total station method has higher accuracy in tunnel deformation monitoring, this method has lower efficiency and cannot meet the requirement of data acquisition speed. Some studies have explored the convergence diameter calculation method based on the mobile laser scanning method achieving good results; however, the specific segment displacement location cannot be determined. Other studies have established deformation prediction models through the geological environment around the tunnel; however, their prediction accuracy is not ideal. Based on the mobile laser scanning system, the point cloud data in the tunnel can be obtained quickly and conveniently, in contrast with the traditional contact measurement method. The orthophoto image was generated from the point cloud data, and an image processing algorithm was used to automatically identify the position of the circumferential seam and the position of the capping block. Finally, the point cloud data were retrieved in a reverse manner for the calculation of tunnel displacement, which quickly provided large and accurate tunnel displacement data for tunnel monitoring work, greatly reducing the field workload and human–computer interaction, improving the operation efficiency. The current algorithm was heavily dependent on the image quality, the quality of the tunnel point cloud, and the tunnel design data. Furthermore, the recognition effect is unsuitable for certain capping blocks and circumferential seams whose contours are not clear enough. To further improve the recognition effect, we propose the use of deep learning neural networks to recognize the capping blocks and circumferential seams.

Author Contributions

Conceptualization, Z.Y. and H.S.; methodology, Z.Y.; software, L.D.; validation, Z.Y., H.S.; formal analysis, Z.Y.; investigation, Z.Y.; resources, Z.Y.; data curation, Z.Y.; writing—original draft preparation, Z.Y., H.S.; writing—review and editing, H.S., L.D.; visualization, Z.Y.; supervision, Z.Y.; project administration, Z.Y., R.Z.; funding acquisition, R.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key R&D Program of China (2016YFB1200602-02), National Science Foundation of China (No. 41771482 and No. 41671451).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Authors would like to thank Haizhi Ma for his help and inspiration, which greatly helped them in writing this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yu, C.J.; Sun, R.Y. The Application of Three-Dimensional Laser Scanner and Vehicle Mounted Mobile Measurement System in Engineering Survey. Surv. Mapp. Spat. Geogr. Inf. 2019, 42, 188–190. [Google Scholar]
  2. Holgado-Barco, A.; Riveiro, B.; González-Aguilera, D.; Arias, P. Automatic Inventory of Road Cross-Sections From Mobile Laser Scanning System. Comput. Aid. Civ. Infrastruct. Eng. 2017, 32, 3–17. [Google Scholar] [CrossRef] [Green Version]
  3. Xiao, W.; Vallet, B.; Brédif, M.; Paparoditis, N. Street Environment Change Detection from Mobile Laser Scanning Point Clouds. ISPRS J. Photogramm. 2015, 107, 38–49. [Google Scholar] [CrossRef]
  4. Yang, B.; Dong, Z. A Shape-Based Segmentation Method for Mobile Laser Scanning Point Clouds. ISPRS J. Photogramm. 2013, 81, 19–30. [Google Scholar] [CrossRef]
  5. Zhou, L.; Vosselman, G. Mapping Curbstones in Airborne and Mobile Laser Scanning Data. Int. J. Appl. Earth Obs. Geoinf. 2012, 18, 293–304. [Google Scholar] [CrossRef]
  6. Wang, Z.F.; Yao, L.B.; Zhou, B. Mobile Measurement System Design and Key Technology Research. Bull. Surv. Mapp. Suppl. 2015, 1, 231–234. [Google Scholar]
  7. Chang, X.J. Study on the Control Technology of Segment Relocation in Metro Shield Tunnel Construction. Mod. Urban Rail Transit 2020, 4, 37–40. [Google Scholar]
  8. Ma, J.; Chen, X.W.; Chen, X.H. Analysis of the Influence of Bolt Diameter on the Displacement of Pipe Segments. Pearl River Water Transp. 2020, 7, 91–92. [Google Scholar]
  9. Yuan, H.; Liu, C.; Lu, J.; Deng, C.; Gong, L. Measurement principle and accuracy analysis of three-dimensional free station tunnel non-contact monitoring with total station. Eng. Invest. 2012, 40, 63–67. [Google Scholar]
  10. Li, G.Y.; Fan, B.X. Precision engineering survey technology and its development. Acta Geod. Cartogr. Sin. 2017, 46, 1742–1751. [Google Scholar]
  11. Luo, Y.; Chen, J.; Xi, W.; Zhao, P.; Qiao, X.; Deng, X.; Liu, Q. Analysis of tunnel displacement accuracy with total station. Measurement 2016, 83, 29–37. [Google Scholar] [CrossRef]
  12. Tang, W.Y. The application of 3D laser scanning technology in the re-measurement of subway tunnel axis. Eng. Constr. Des. 2020, 2, 60–62. [Google Scholar]
  13. Tang, K.; Hua, X.H.; Cheng, W.; Xia, Z. Research on the Monitoring Method of Building Deformation Based on Three-Dimensional Laser Scanning. Geogr. Inf. Surv. Mapp. 2013, 38, 54–55. [Google Scholar]
  14. Yin, H.; Feng, Q.H.; Liao, Z.H.; Jiang, L.L.; Cai, J. Disease Tunnel Monitoring Based on 3D Laser Scanning Technology. Undergr. Space Eng. 2014, 10, 895–901. [Google Scholar]
  15. Zhang, M.Z.; Zhang, M.D. Application of 3D Laser Scanning Technology in Deformation Monitoring of Subway Tunnels. City Surv. 2018, 3, 144–147. [Google Scholar]
  16. Cheng, Z.C. Application of 3D Laser Scanning Technology in Deformation Monitoring of Rectangular Tunnel Segments. Construction 2017, 39, 1409–1411. [Google Scholar]
  17. Cheng, Y.J.; Qiu, W.; Lei, J. Automatic extraction of tunnel lining cross-sections from terrestrial laser scanning point clouds. Sensors 2016, 16, 1648. [Google Scholar] [CrossRef]
  18. Gao, H.; Li, K.; Ma, Q.M.; Han, Z.S.; Sun, P.C. Application research of mobile three-dimensional laser measurement system in disease monitoring of subway operating tunnels. Bull. Surv. Mapp. 2019, 8, 101–161. [Google Scholar]
  19. Shen, Z.Y.; Tan, Z. Application Research of grp5000 Tunnel Inspection Car in Shanghai Metro. Shanxi Archit. 2013, 39, 158–159. [Google Scholar]
  20. Dianyou, B. The Application of Three-Dimensional Laser Scanning Technology in the Investigation of Subway Structure Status. Urban Geotech. Invest. Surv. 2019, 4, 126–128. [Google Scholar]
  21. Cui, H.; Ren, X.; Mao, Q.; Hu, Q.; Wang, W. Shield Subway Tunnel Deformation Detection Based on Mobile Laser Scanning. Autom. Constr. 2019, 106, 102889. [Google Scholar] [CrossRef]
  22. Lueprasert, P.; Jongpradist, P.; Jongpradist, P.; Suwansawat, S. Numerical investigation of tunnel deformation due to adjacent loaded pile and pile-soil-tunnel interaction. Tunn. Undergr. Space Technol. 2017, 70, 166–181. [Google Scholar] [CrossRef]
  23. Parsapour, D.; Fahimifar, A. Semi-analytical solution for time-dependent deformations in swelling rocks around circular tunnels. Geosci. J. 2016, 20, 517–528. [Google Scholar] [CrossRef]
  24. Paraskevopoulou, C.; Diederichs, M. Analysis of time-dependent deformation in tunnels using the Convergence-Confinement Method. Tunn. Undergr. Space Eng. 2018, 71, 62–80. [Google Scholar] [CrossRef]
  25. Sainoki, A.; Tabata, S.; Mitri, H.; Fukuda, D.; Kodama, J. Time-dependent tunnel deformations in homogeneous and heterogeneous weak rock formations. Comput. Geotech. 2017, 92, 186–200. [Google Scholar] [CrossRef]
  26. Janda, T.; Šejnoha, M.; Šejnoha, J. Applying Bayesian approach to predict deformations during tunnel construction. Int. J. Numer. Anal. Methods Géoméch. 2018, 42, 1765–1784. [Google Scholar]
  27. Sun, H.; Xu, Z.; Yao, L.; Zhong, R.; Du, L.; Wu, H. Tunnel Monitoring and Measuring System Using Mobile Laser Scanning: Design and Deployment. Remote Sens. 2020, 12, 730. [Google Scholar] [CrossRef] [Green Version]
  28. Du, L.M.; Zhong, R.F.; Sun, H.L.; Zhu, Q.; Zhang, Z. Study of the Integration of the CNU-TS-1 Mobile Tunnel Monitoring System. Sensors 2018, 18, 420. [Google Scholar] [CrossRef] [Green Version]
  29. Arastounia, M. Automated As-Built Model Generation of Subway Tunnels From Mobile Lidar Data. Sensors 2016, 16, 1486. [Google Scholar] [CrossRef] [Green Version]
  30. Du, L.M.; Zhong, R.F.; Sun, H.L.; Liu, Y.J.; Wu, Q. Cross-Section Positioning Based on a Dynamic Mls Tunnel Monitoring System. Photogramm. Rec. 2019, 34, 244–265. [Google Scholar] [CrossRef]
  31. Du, L.; Zhong, R.; Sun, H.; Wu, Q. Automatic Monitoring of Tunnel Deformation Based on High Density Point Clouds Data. ISPRS—International Archives of the Photogrammetry. Remote Sens. Spat. Inf. Sci. 2017, 42, 353. [Google Scholar] [CrossRef] [Green Version]
  32. Sun, J.; Sun, H.L.; Zhong, R.F.; Han, Y.L. Deformation Detection Method of Mine Tunnel Based on Mobile Detection System. Sensors 2020, 20, 5400. [Google Scholar] [CrossRef] [PubMed]
  33. Zhu, Q.; Zhong, R.F.; Sun, H.; Du, L.M.; Zhen, Z. Tunnel Boundary Detection Method of Mobile Laser Detection System. Sci. Surv. Mapp. 2020, 45, 118–125. [Google Scholar]
  34. Xu, D.S.; Zhao, Y.M.; Liu, H.B.; Zhu, H.H. Deformation Monitoring of Metro Tunnel with a New Ultrasonic-Based System. Sensors 2017, 17, 1758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Yertutanol, K.; Akgün, H.; Sopacı, E. Displacement monitoring, displacement verification and stability assessment of the critical sections of the Konak tunnel, İzmir, Turkey. Tunn. Undergr. Space Technol. 2020, 101, 103357. [Google Scholar] [CrossRef]
  36. Jiang, H.; Cheng, J.; Zhang, J.; Zhang, J.; Su, Y.; Zheng, Y. Principle and application of in-situ monitoring system for ground displacement induced by shield tunnelling. Tunn. Undergr. Space Technol. 2021, 112, 103905. [Google Scholar] [CrossRef]
  37. Sun, H.; Liu, S.; Zhong, R.; Du, L. Cross-Section Deformation Analysis and Visualization of Shield Tunnel Based on Mobile Tunnel Monitoring System. Sensors 2020, 20, 1006. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic diagram of the mobile tunnel inspection vehicle.
Figure 1. Schematic diagram of the mobile tunnel inspection vehicle.
Sensors 21 04407 g001
Figure 2. Flow chart of the displacement calculation method.
Figure 2. Flow chart of the displacement calculation method.
Sensors 21 04407 g002
Figure 3. Before section fitting (a) and after section fitting (b).
Figure 3. Before section fitting (a) and after section fitting (b).
Sensors 21 04407 g003
Figure 4. Schematic diagram of slice point cloud projection.
Figure 4. Schematic diagram of slice point cloud projection.
Sensors 21 04407 g004
Figure 5. Orthophoto generation effect.
Figure 5. Orthophoto generation effect.
Sensors 21 04407 g005
Figure 6. Recognition results of longitudinal ring seams.
Figure 6. Recognition results of longitudinal ring seams.
Sensors 21 04407 g006
Figure 7. (a) Schematic segmentation diagram. (b) Capping block (KP) position in point cloud.
Figure 7. (a) Schematic segmentation diagram. (b) Capping block (KP) position in point cloud.
Sensors 21 04407 g007
Figure 8. (a) Canny edge detection result. (b) Hough line detection result.
Figure 8. (a) Canny edge detection result. (b) Hough line detection result.
Sensors 21 04407 g008
Figure 9. (a) Schematic radial displacement diagram and site photo of radial dislocation. (b) Local symmetric foot drop method.
Figure 9. (a) Schematic radial displacement diagram and site photo of radial dislocation. (b) Local symmetric foot drop method.
Sensors 21 04407 g009
Figure 10. Sectional superposition results.
Figure 10. Sectional superposition results.
Sensors 21 04407 g010
Figure 11. Operation photos of (a) Leica (b) Faro (c) Z + F.
Figure 11. Operation photos of (a) Leica (b) Faro (c) Z + F.
Sensors 21 04407 g011
Figure 12. Repeatability deviation of A1P–A2P displacement.
Figure 12. Repeatability deviation of A1P–A2P displacement.
Sensors 21 04407 g012
Figure 13. Repeatability deviation of BP–A1P displacement.
Figure 13. Repeatability deviation of BP–A1P displacement.
Sensors 21 04407 g013
Figure 14. Repeatability deviation of KP–BP displacement.
Figure 14. Repeatability deviation of KP–BP displacement.
Sensors 21 04407 g014
Figure 15. Circumferential displacement contrast difference.
Figure 15. Circumferential displacement contrast difference.
Sensors 21 04407 g015
Figure 16. BP–A1P contrast difference.
Figure 16. BP–A1P contrast difference.
Sensors 21 04407 g016
Figure 17. KP–BP contrast difference.
Figure 17. KP–BP contrast difference.
Sensors 21 04407 g017
Table 1. Circumferential displacement round-trip precision comparison.
Table 1. Circumferential displacement round-trip precision comparison.
NumberForward Displacement (mm)Backward Displacement (mm)Absolute Deviation (mm)
17.710.12.4
28.79.40.7
39.712.12.4
47.19.92.8
512.510.71.8
69.19.30.2
710.814.03.2
89.98.61.3
97.210.12.9
107.59.11.6
1110.810.50.3
1215.718.02.3
1310.19.60.5
147.87.80.0
1511.212.00.8
168.69.30.7
177.79.11.4
188.78.00.7
1911.213.82.6
207.79.51.8
2112.111.30.8
2213.211.41.8
2312.112.40.3
248.07.80.2
259.87.42.4
Absolute Mean Deviation (mm)1.4
Table 2. Radial displacement absolute accuracy comparison.
Table 2. Radial displacement absolute accuracy comparison.
NumberMeasured Value of Vernier Calipers (mm)Radial Displacement Measured by Proposed Method (mm)Absolute Deviation (mm)
13.54.51.2
25.33.02.3
32.15.13.0
47.44.92.5
51.53.72.2
63.86.02.2
74.92.22.7
81.91.60.3
93.66.12.5
103.24.10.9
Absolute Mean Deviation (mm)1.98
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yue, Z.; Sun, H.; Zhong, R.; Du, L. Method for Tunnel Displacements Calculation Based on Mobile Tunnel Monitoring System. Sensors 2021, 21, 4407. https://doi.org/10.3390/s21134407

AMA Style

Yue Z, Sun H, Zhong R, Du L. Method for Tunnel Displacements Calculation Based on Mobile Tunnel Monitoring System. Sensors. 2021; 21(13):4407. https://doi.org/10.3390/s21134407

Chicago/Turabian Style

Yue, Zeyu, Haili Sun, Ruofei Zhong, and Liming Du. 2021. "Method for Tunnel Displacements Calculation Based on Mobile Tunnel Monitoring System" Sensors 21, no. 13: 4407. https://doi.org/10.3390/s21134407

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop