Tunnel displacement measurement method and tunnel displacement measurement system

By employing pseudo-targets within the tunnel and using point cloud data to limit displacement measurement regions, the method addresses high installation costs and processing challenges, achieving rapid and accurate tunnel displacement measurement.

JP7829889B2Active Publication Date: 2026-03-16KAJIMA CORP +1
View PDF 6 Cites 0 Cited by

Patent Information

Authority / Receiving Office
JP · JP
Patent Type
Patents
Current Assignee / Owner
Filing Date
2024-06-18
Publication Date
2026-03-16

AI Technical Summary

Technical Problem

Existing tunnel displacement measurement methods using 3D laser scanners face challenges with high installation costs and difficulty in high-speed processing due to the need for numerous targets and processing of countless points.

Method used

A method utilizing point cloud data from a 3D laser scanner to set pseudo-targets within the tunnel, limiting displacement measurement to specific regions, and using these pseudo-targets for high-speed and accurate displacement measurement by comparing position information over time.

Benefits of technology

Enables rapid, cost-effective, and precise measurement of tunnel internal displacement without the need for dedicated targets, improving processing speed and accuracy.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure 0007829889000001
    Figure 0007829889000001
  • Figure 0007829889000002
    Figure 0007829889000002
  • Figure 0007829889000003
    Figure 0007829889000003
Patent Text Reader

Abstract

To provide a method, a device and a system capable of measuring the displacement of the inner space of a tunnel at high speed.SOLUTION: 1) In the initial measurement step, the unevenness feature of the target 100a to 100e is acquired. In addition, the point cloud of the inner surface of the tunnel including the target 100a to 100e is obtained as the first point cloud, and the unevenness feature is extracted from the first point cloud to obtain the first position information of the target 100a to 100e. In the displacement measurement step (2) performed after a predetermined time elapses, a search region including a position specified based on the first positional information and having a larger area than the targets 100a to 100e is set, point cloud of the search region is acquired, and the targets 100a to 100e are extracted from the point cloud of the search region based on the unevenness feature information to acquire the second positional information. Then, the displacement of the inner space of the tunnel is measured based on the difference between the first position information and the second position information.SELECTED DRAWING: Figure 17
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to a technique for measuring the settlement of the tunnel ceiling and the displacement of the internal cavity shape during tunnel excavation work.

Background Art

[0002] In tunnel excavation work using the NATM (New Austrian Tunneling Method), blasting is carried out with explosives for each tunnel face. After removing the generated slip, jacking, support work, primary lining work, and installation of rock bolts are performed. Such construction is repeated at a pitch of about 1.5 m for tunneling. In tunnel excavation work, displacement measurement (A measurement) of the ground and support members is regularly performed to ensure the safety and quality of the work. In A measurement, the internal cavity displacement measurement is carried out for the purpose of grasping the behavior of the surrounding ground and the deformation mode of the support, confirming the safety of the construction and the validity of the support, and considering the timing of placing the lining. Conventionally, in A measurement, a control section is set at predetermined intervals, and the coordinates of targets provided at the top and side walls of the tunnel in the control section are measured by a distance measuring instrument such as a total station (see, for example, Patent Document 1).

[0003] A measurement aims to measure the displacement of the tunnel wall surface in order to predict the displacement of the surrounding wall surface based on the displacement at the target installation location. A technique for measuring the tunnel wall surface shape by a 3D laser scanner is known. For example, Patent Document 2 discloses a method for measuring the displacement of the internal cavity of a tunnel in which a plurality of targets are installed on the tunnel wall surface, the shape of the tunnel wall surface including the targets is measured by a 3D laser scanner, and the tunnel cross-sectional shape data at the same location is compared based on the coordinates of the targets to track the displacement of the tunnel (see Patent Document 2). However, when installing a plurality of targets on the tunnel wall surface, there are problems such as complicated installation work and high cost of the targets.

[0004] Therefore, a technique is known that allows for the measurement of tunnel internal displacement without using a target by utilizing scan data from a 3D laser scanner to determine internal displacement at countless points (see Patent Document 3). According to this, since there is no need to attach a target to the measurement cross-section, the measurement of internal displacement can be performed simply and efficiently. However, the tunnel internal displacement measurement method in Patent Document 3 has the problem that high-speed processing is difficult because it determines internal displacement at countless points measured by a 3D laser scanner. [Prior art documents] [Patent Documents]

[0005] [Patent Document 1] Japanese Patent Publication No. 2004-170164 [Patent Document 2] Japanese Patent Publication No. 2014-002027 [Patent Document 3] Japanese Patent Publication No. 2012-058167 [Overview of the project] [Problems that the invention aims to solve]

[0006] In view of the above circumstances, the present invention provides a method for measuring the displacement inside a tunnel at high speed. ,oh The purpose is to provide a calling system. [Means for solving the problem]

[0007] As a result of diligent research, the inventors have found that, rather than using a 3D laser scanner to perform surface measurement, applying point cloud data obtained from a 3D laser scanner to a point measurement technique enables simple and high-speed displacement measurement.

[0008] In other words, the tunnel displacement measurement method of the present invention is a method for measuring the displacement of the tunnel interior over time by setting a target using point cloud data of the inner surface of the tunnel acquired by a 3D laser scanner, and consists of 1) an initial measurement step and 2) a displacement measurement step, the 1) initial measurement step being to acquire surface feature data, which is the surface feature of the target, After obtaining the surface feature data, 1) A displacement measurement step performed after a predetermined time has elapsed includes: acquiring point cloud data of the inner surface of the tunnel including the target as first point cloud data, extracting surface feature data from the first point cloud data to obtain first position information of the target; setting a search area larger in area than the target including the position identified based on the first position information, acquiring point cloud data of the search area, extracting the target from the point cloud data of the search area based on surface feature data to obtain second position information, and measuring the displacement of the tunnel interior based on the difference between the first position information and the second position information. The target is a physical target component installed on the wall surface inside the tunnel. It is characterized by the following: Another method 1) The initial measurement step includes acquiring point cloud data of the tunnel interior as first point cloud data, cutting out a predetermined area of ​​the tunnel interior as a target from the first point cloud data, extracting surface feature data which is the surface feature of the target, and registering the first position information and surface feature data of the target, and is performed after a predetermined time has elapsed. 2) The displacement measurement step includes setting a search area that is larger in area than the target and includes the position identified based on the first position information, acquiring point cloud data of the search area, extracting the target from the point cloud data of the search area based on the surface feature data to acquire second position information, and measuring the displacement of the tunnel interior based on the difference between the first position information and the second position information. That's fine. .

[0009] In this way, instead of using all of the acquired point cloud data for displacement measurement, limiting the comparison to the target range enables rapid processing.

[0010] 1) The initial measurement step and 2) the displacement measurement step each comprise the following steps.

[0011] 1) Initial measurement step 1-A) Scanner coordinate acquisition step to obtain the position coordinates of a 3D laser scanner. 1-B) A first point cloud data acquisition step in which point cloud data of the sprayed concrete surface inside the tunnel is acquired as the first point cloud data. 1-C) A pseudo-target setting step in which at least three regions are extracted from the first point cloud data and set as pseudo-targets. 1-D) A first surface feature data extraction step in which surface features based on the plane of the pseudo-target are extracted as first surface feature data. 1-E) A first registration step in which the position coordinates of a pseudo-target, which is first position information, and first surface feature data are registered in association with the position coordinates of a 3D laser scanner.

[0012] 2) Displacement measurement step 2-A) A second point cloud data acquisition step in which point cloud data of the sprayed concrete surface is acquired as second point cloud data. 2-B) A search area setting step in which the search area is set from the second point cloud data based on the position coordinates of the pseudo-target. 2-C) A second surface feature data extraction step in which surface features based on the plane of the pseudo-target are extracted from the search area as second surface feature data. 2-D) A pseudo-target identification step in which the position coordinates of a pseudo-target, which is second positional information, are identified based on the similarity between the surface features of the first surface feature data and the partial surface features of the second surface feature data. 2-E) A second registration step in which the position coordinates of the identified pseudo-target are registered in association with the position coordinates of the 3D laser scanner. 2-F) A displacement measurement step that measures the displacement of the tunnel interior based on the changes in the position coordinates of a pseudo-target set in the initial measurement step and the pseudo-target identified in the displacement measurement step.

[0013] By using a pseudo target, it is possible to measure the displacement of the internal shape of a tunnel at low cost, simply, and with high precision without using an existing target (prism). Generally, for crown settlement measurement and internal displacement measurement, targets are provided every 10 - 30 m, and measurements are carried out at a frequency of twice a day in many cases and at least once a week in few cases. When installing a new pseudo target, for the location where the new pseudo target is installed, 1) perform the initial measurement step, and for the locations where the initial measurement has already been completed, 2) perform the displacement measurement step. By repeating the measurements in this way, the shape displacement of the tunnel interior can be accurately measured.

[0014] 1 - A) For the scanner coordinate acquisition step, a wide range of methods can be used to acquire the position coordinates of the 3D laser scanner in the tunnel interior. For example, the position coordinates of the 3D laser scanner can be acquired by using a target for the laser scanner attached to a point with known position coordinates, or when the 3D laser scanner has the function as a total station, the absolute coordinates or tunnel coordinates can be acquired using the function as a total station. The position coordinates to be acquired can be either absolute coordinates (global coordinates) or tunnel coordinates (local coordinates).

[0015] Rather than regarding all the point cloud data acquired in the 1 - B) first point cloud data acquisition step as innumerable targets, in the 1 - C) pseudo target setting step, by cutting out three or more regions from the point cloud data and setting them as pseudo targets, the locations where displacement is to be measured can be limited, and the processing speed of displacement measurement can be improved. In order to improve the processing speed while maintaining the measurement accuracy, the number of regions to be cut out is preferably 3 - 5 for each measurement cross - section. Since the surface of the sprayed concrete in the tunnel inner space has a fine uneven shape, it is possible to cut out a unique area from the point cloud data. Here, the "surface of the sprayed concrete in the tunnel inner space" refers to the part of the tunnel wall surface or the support and protection work after blasting where concrete spraying has been performed. The "surface of the sprayed concrete" does not have a strict meaning and includes parts where the concrete spraying is not sufficient. Therefore, the rock bolts fixed after the concrete spraying are also included. As the three-dimensional laser scanner, a LiDAR (Light Detection And Ranging) scanner is preferably used. Each area regarded as a pseudo-target is an area cut out in a predetermined shape such as a polygon such as a square, a circle, an elliptical shape, etc., and is cut out in a substantially square shape with a vertical and horizontal dimension of 150 mm each, for example.

[0016] 1-E) The position coordinates of the pseudo-target registered in the first registration step are preferably the first point cloud data included in the pseudo-target, but may also be the centroid, normal vector, and plane data described later.

[0017] 1 -C) The pseudo-target setting step preferably consists of a first cutting step of cutting out an area, a first centroid calculation step of calculating the centroid of the pseudo-target based on the point cloud data included in the cut-out area, a first normal vector calculation step of calculating the normal vector of the pseudo-target based on the point cloud data and the centroid, and a first flattening step of making the pseudo-target into plane data based on the centroid and the normal vector. The normal vector can be calculated as a line segment passing through the centroid coordinates from, for example, the eigenvector of the covariance matrix of the point cloud data using a known method related to the processing of the point cloud data. By calculating the centroid and the normal vector and performing the flattening process instead of directly using the point cloud data as a comparison target, it becomes easy to extract feature data, and it is possible to identify and measure the displacement of the pseudo-target quickly and with high accuracy.

[0018] 1-D) The first surface feature data extraction step preferably involves calculating cross-sectional data from the relative positions of each point in the normal direction for at least one cross-section in the planar data. Calculating cross-sectional data facilitates comparison of surface feature data. Preferably, there are multiple cross-sectional data to be calculated, and more preferably, there are multiple cross-sectional data with different orientations.

[0019] 2 -B) Preferably, the search area setting step includes: a second extraction step of extracting a search area from second point cloud data based on the position coordinates of a pseudo-target set in the initial measurement step; a second centroid calculation step of calculating the centroid of the search area based on the second point cloud data included in the extracted search area; a second normal calculation step of calculating the normal of the search area based on the second point cloud data and the centroid; and a second planarization step of converting the search area into planar data based on the centroid and the normal. Including such steps makes it easier to identify the pseudo-target.

[0020] 2 -C) The second surface feature data extraction step preferably involves calculating cross-sectional data from the relative positions of each point in the normal direction for at least one cross-section in the planar data. This makes it easier to identify pseudo-targets by comparing the first surface feature data with the second surface feature data.

[0021] 2 -D) The pseudo-target identification step preferably involves comparing the cross-sectional data calculated in 1-D) the first surface feature data extraction step with the cross-sectional data calculated in 2-C) the second surface feature data extraction step to identify a pseudo-target based on the similarity of the relative positions of each point in the normal direction.

[0022] The Step 2, the cutting step, may also involve cutting out an area with an area approximately four times or less the size of the pseudo-target, centered on the position coordinates of the pseudo-target set in step 1) the initial measurement step. In the second extraction step, the search area is extracted from the second point cloud data based on the position coordinates of the pseudo-target registered in the first registration step. However, since the pseudo-target may move or change shape after the initial measurement, it is necessary to extract the area with a predetermined margin. Therefore, by setting the length and width of the pseudo-target to approximately twice its original size and extracting an area with an area of ​​approximately four times or less its original size, it is possible to prevent the pseudo-target from becoming unidentifiable due to its movement or change in shape.

[0023] 2 -F) Preferably, the change in the position coordinates of the pseudo-target in the displacement measurement step is the change in the position coordinates of the first point cloud data and the second point cloud data included in the pseudo-target. A change in position coordinates refers to a change in the three-dimensional position coordinates of each point. By comparing only the point cloud data included in the pseudo-target, rather than all of the acquired point cloud data, high-speed and high-precision displacement measurement can be achieved. Alternatively, the displacement of the tunnel's internal shape may be measured based on changes in the center of gravity, normals, and planar data, or the displacement may be measured based on changes in the first and second surface feature data.

[0024] doubt The target's position coordinates can be absolute coordinates or tunnel coordinates.

[0025] doubtThe pseudo-targets are preferably set in at least one of the following areas in tunnel excavation work using the NATM (New Austrian Tunneling Method): areas where A, B, and C lines can be measured, and areas where D line can be measured, based on the location or shape of the tunnel cavity, support structures, or rock bolts. By setting pseudo-targets in such areas, A measurement can be easily performed. The setting may be done automatically based on the location or shape of the tunnel cavity, support structures, or rock bolts, or it may be set manually. In addition, even outside of these areas, for example, areas with significant irregularities on the tunnel cavity surface may be detected and set automatically or by visual inspection by a worker.

[0026] The target member has protrusions and the shape of the tunnel cavity has irregular characteristics. The target is preferably the anchoring plate of a rock bolt installed on the inner surface of the tunnel. The target is also preferably the steel support structure installed on the inner surface of the tunnel. This eliminates the need to install a dedicated target for measurement. Furthermore, if the target has protrusions at multiple locations in the vertical direction, the characteristic uneven shape of the protrusions on the target's surface allows for accurate determination of horizontal and vertical displacements as changes in position coordinates. This makes it possible to measure the displacement of the tunnel's internal shape inexpensively, simply, and with high accuracy. In other words, the target The components are installed on the walls inside the tunnel and have a textured surface on the portion exposed to the inner surface. Any metal object will do. 。

[0027] Three-dimensional laser scanners are installed along each of the walls on both sides in the tunnel width direction. Each three-dimensional laser scanner acquires point cloud data from the wall on which it is installed and the wall on the opposite side. Topographic feature data, which represents the surface characteristics of the target, can be obtained from the imaging results obtained by imaging the target with an imaging device.

[0028] Tunnel displacement measurement of the present invention system This device uses point cloud data of the tunnel's inner surface acquired by a 3D laser scanner to set a target and measure the displacement of the tunnel's interior space over time. A system comprising: and a physical target member installed as a target on the wall surface inside the tunnel. And, The device is It consists of an initial measurement unit and a displacement measurement unit, and the initial measurement unit acquires surface feature data, which is the surface feature of the target. After obtaining the surface feature data,The system is characterized by acquiring point cloud data of the tunnel's inner surface, including the target, as the first point cloud data, extracting surface feature data from the first point cloud data to obtain the first position information of the target, setting a search area larger in area than the target and including the position identified based on the first position information, acquiring point cloud data of the search area, extracting the target from the point cloud data of the search area based on surface feature data to obtain the second position information, and measuring the displacement of the tunnel's interior based on the difference between the first and second position information. ru.

[0029] As another device The initial measurement unit acquires point cloud data of the tunnel's inner surface as first point cloud data, extracts a predetermined area of ​​the tunnel's inner surface as a target from the first point cloud data, extracts surface feature data which is the surface feature of the target, and registers the first position information and surface feature data of the target. The displacement measurement unit sets a search area that includes the position identified based on the first position information and has a larger area than the target, acquires point cloud data of the search area, extracts the target from the point cloud data of the search area based on the surface feature data to acquire second position information, and measures the displacement of the tunnel's interior based on the difference between the first position information and the second position information. It is also acceptable to do so. . The initial measurement unit and the displacement measurement unit each include the following means.

[0030] 1) Initial measurement section 1-a) A means for acquiring scanner coordinates to obtain the position coordinates of a 3D laser scanner. 1-b) A first point cloud data acquisition means for acquiring point cloud data of the sprayed concrete surface inside the tunnel as first point cloud data. 1-c) A pseudo-target setting means for extracting at least three regions from the first point cloud data and setting them as pseudo-targets. 1-d) A first surface feature data extraction means for extracting surface features based on the plane of a pseudo-target as first surface feature data. 1-e) A first registration means for registering the position coordinates of a pseudo-target, which is first position information, and first surface feature data, in association with the position coordinates of a 3D laser scanner. 2) Displacement measurement unit 2-a) A second point cloud data acquisition means for acquiring point cloud data of the sprayed concrete surface as a second point cloud data. 2-b) Search area setting means for setting a search area from second point cloud data based on the position coordinates of a pseudo-target. 2-c) A second surface feature data extraction means for extracting surface features based on the plane of a pseudo-target from the search area as second surface feature data. 2-d) A pseudo-target identification means that identifies the position coordinates of a pseudo-target, which is second position information, based on the similarity between the surface features of the first surface feature data and the partial surface features in the second surface feature data. 2-e) A second registration means for registering the position coordinates of an identified pseudo-target in association with the position coordinates of a 3D laser scanner. 2-f) Displacement measuring means for measuring the displacement of the space inside a tunnel based on the changes in the position coordinates of a pseudo-target set in the initial measurement means and the pseudo-target identified in the displacement measuring means.

[0031] The tunnel displacement measurement system of the present invention comprises a tunnel displacement measurement device and A physical target component is installed as a target on the wall surface inside the tunnel, It holds. The target member has protrusions and has an uneven shape in relation to the interior of the tunnel. Alternatively, the target member is a metal fitting installed on the wall surface of the tunnel interior, with an uneven shape on the portion exposed to the inner surface. If the tunnel is constructed using the NATM method, the target is the anchoring plate of the rock bolt or the steel support structure. If the tunnel is constructed using the shield tunneling method, the target is the gripping hardware or joint hardware installed on the inner surface of the tunnel segment piece. [Effects of the Invention]

[0032] The present invention: Method for measuring tunnel displacement ,oh According to the Yobi system, it has the effect of being able to measure the displacement inside the tunnel at high speed. [Brief explanation of the drawing]

[0033] [Figure 1] Functional block diagram of the tunnel displacement measuring device [Figure 2] Schematic diagram of the tunnel displacement measurement device. [Figure 3] Schematic flowchart of tunnel displacement measurement method [Figure 4] Diagram (1) illustrating the initial measurement step and displacement measurement step. [Figure 5] Diagram (2) illustrating the initial measurement step and displacement measurement step. [Figure 6] Initial measurement flow chart [Figure 7] Image diagram of the placement of simulated targets [Figure 8] Pseudo-target setting flowchart [Figure 9] Diagram illustrating the steps for setting a simulated target. [Figure 10] Image diagram of simulated target and feature data [Figure 11] Feature data extraction diagram [Figure 12] Displacement measurement flowchart [Figure 13] Search Area Setting Flowchart [Figure 14] Image diagram of the search area setting [Figure 15] Diagram illustrating the steps for identifying a pseudo-target (1) [Figure 16] Diagram illustrating the pseudo-target identification step (2) [Figure 17] Schematic diagram of the tunnel displacement measurement system. [Figure 18] Image of a rock bolt plate [Figure 19] Target image [Figure 20] Schematic diagram showing the relationship between the number of protrusions and measurement accuracy. [Figure 21] Image of a protrusion [Figure 22] Schematic diagram of the tunnel displacement measurement system. [Figure 23] Conceptual diagram of a protrusion on a steel support structure. [Modes for carrying out the invention]

[0034] Hereinafter, an example of an embodiment of the present invention will be described in detail with reference to the drawings. It should be noted that the scope of the present invention is not limited to the following embodiments or illustrated examples, and numerous modifications and variations are possible. [Examples]

[0035] Figure 1 shows a functional block diagram of the tunnel displacement measuring device. Figure 2 shows an image diagram of the configuration of the tunnel displacement measuring device. As shown in Figure 1, the tunnel displacement measuring device 1 consists of an initial measurement unit 2 and a displacement measuring unit 3. The initial measurement unit 2 includes a scanner coordinate acquisition means 21, a first point cloud data acquisition means 22, a pseudo-target setting means 23, a first unevenness feature data extraction means 24, and a first registration means 26. The displacement measuring unit 3 includes a second point cloud data acquisition means 32, a search area setting means 33, a second unevenness feature data extraction means 34, a pseudo-target identification means 35, a second registration means 36, and a displacement measuring means 37. Furthermore, as shown in Figure 2, the tunnel displacement measuring device 1 consists of a LiDAR scanner 4, a computer 5, and scanner targets (11a to 11c), with the LiDAR scanner 4 and the computer 5 connected by wire or wireless. The scanner targets (11a to 11c) are attached to points in the interior of the tunnel 7 whose position coordinates are known. The scanner targets (11a to 11c) can be installed, for example, on the wall surface 72 of the tunnel 7.

[0036] The scanner coordinate acquisition means 21 acquires the position coordinates of the 3D laser scanner. A LiDAR scanner 4 is used as the 3D laser scanner. In this embodiment, the scanner coordinate acquisition means 21 acquires the absolute coordinates of the LiDAR scanner 4 using scanner targets (11a to 11c). The first point cloud data acquisition means 22 acquires point cloud data of the sprayed concrete surface (inner surface of the tunnel) inside the tunnel using the 3D laser scanner. The pseudo-target setting means 23 extracts at least three predetermined regions from the first point cloud data and sets them as pseudo-targets (targets). The first surface feature data extraction means 24 extracts first surface feature data from the pseudo-targets. The first registration means 26 registers the position coordinates (first position information) of the pseudo-targets and the first surface feature data in association with the position coordinates of the 3D laser scanner.

[0037] The second point cloud data acquisition means 32 acquires point cloud data of the sprayed concrete surface inside the tunnel using a 3D laser scanner. The search area setting means 33 sets the search area of ​​a pseudo-target (target) from the second point cloud data. The search area includes a location identified based on the position coordinates (first position information) of the pseudo-target and is a larger area than the pseudo-target. The second surface feature data extraction means 34 extracts second surface feature data from the point cloud data of the search area. The pseudo-target identification means 35 extracts a pseudo-target based on the similarity between the first surface feature data and the second surface feature data, identifies its position coordinates (second position information), and acquires it. The second registration means 36 registers the identified position coordinates of the pseudo-target in association with the position coordinates of the 3D laser scanner. The displacement measuring means 37 measures the displacement of the tunnel's internal shape based on the change in the position coordinates of a pseudo-target set in the initial measurement unit 2 and the pseudo-target identified in the displacement measuring unit 3 (the difference between the first position information and the second position information).

[0038] In this embodiment, the LiDAR scanner 4 includes a first point cloud data acquisition means 22 and a second point cloud data acquisition means 32, and the computer 5 includes a pseudo-target setting means 23, a first surface feature data extraction means 24, a first registration means 26, a second point cloud data acquisition means 32, a search area setting means 33, a second surface feature data extraction means 34, a pseudo-target identification means 35, a second registration means 36, and a displacement measuring means 37.

[0039] Figure 3 shows a schematic flowchart of the tunnel displacement measurement method. As shown in Figure 3, the tunnel displacement measurement method in this embodiment first performs an initial measurement (Step S01: Initial Measurement Step), and then performs a displacement measurement (Step S02: Displacement Measurement Step). Displacement measurements are performed at a frequency of once a week to twice a day. Specifically, this is determined based on the distance between the measurement position and the tunnel face, and the displacement velocity of the internal displacement. Displacement measurements are performed repeatedly, and once all displacement measurements are completed (Step S03), the displacement measurement is terminated. Figures 4 and 5 illustrate the initial measurement step and the displacement measurement step. As shown in Figure 4 or Figure 5, when the initial measurement step (step S01) is performed for range R2, if the initial measurement step (step S01) has already been completed for range R1 as shown in Figure 5, the initial measurement step (step S01) is performed for range R2, and the displacement measurement step (step S02) is performed for range R1. After a predetermined time has elapsed, as the tunnel excavation work progresses, if a pseudo-target is set up for a new range R3, the initial measurement step (step S01) is performed for range R3, and the displacement measurement step (step S02) is performed for the ranges (R1, R2) where the initial measurement has already been completed. In this way, displacement measurements can be repeatedly performed while increasing the number of pseudo-targets.

[0040] (Regarding the initial measurement step) Figure 6 shows the initial measurement flow diagram. As shown in Figure 6, first, the position coordinates of the 3D laser scanner are obtained (step S11: scanner coordinate acquisition step). As mentioned above, in this embodiment, the absolute coordinates of the LiDAR scanner 4 are obtained using scanner targets (11a~11c). A 3D laser scanner is used to acquire point cloud data of the sprayed concrete surface inside the tunnel (Step S12: First point cloud data acquisition step). In this embodiment, a LiDAR scanner 4 is used to acquire point cloud data of the tunnel face 71 and wall surface 72 in front of the LiDAR scanner 4. Concrete is sprayed onto the wall surface 72. The range of point cloud data acquisition here is not limited to the range set as a pseudo-target; for example, as shown in Figure 5, if it is possible to acquire point cloud data for range R, point cloud data is acquired for the entire range R.

[0041] At least three regions are extracted from the point cloud data and set as pseudo-targets (Step S13: Pseudo-target setting step). In this embodiment, as shown in Figure 2, five regions are extracted to include the measurement cross-section 70 and set as pseudo-targets (6a~6e). Figure 7 shows an image diagram of the arrangement of pseudo-targets. Regarding the placement of the pseudo-targets, as shown in Figure 7, pseudo-targets 6a and 6e are placed at positions where the D measurement line can be measured. In addition, pseudo-targets 6b and 6c are placed at positions where the A measurement line can be measured, pseudo-targets 6c and 6d are placed at positions where the B measurement line can be measured, and pseudo-targets 6b and 6d are placed at positions where the C measurement line can be measured.

[0042] Here, we will explain the method for setting a pseudo-target. Figure 8 shows a flowchart for setting a pseudo-target. Figure 9 is an explanatory diagram of the pseudo-target setting steps, where (1) shows region extraction, (2) shows centroid calculation, (3) shows normal vector calculation, and (4) shows the image of planarization. As shown in Figure 8, in the pseudo-target setting step (step S13), first, a region is extracted from the first point cloud data (step S131: first extraction step). Figure 9(1) shows the point cloud 8 in the extracted region. Next, as shown in Figure 9(2), the centroid G of the pseudo-target is calculated based on the point cloud data included in the extracted region (step S132: first centroid calculation step). Based on the point cloud data and centroid, the normal vector of the pseudo-target is calculated (step S133: first normal vector calculation step). The normal vector can be calculated using known methods for processing point cloud data, for example, as a line segment passing through the centroid coordinates from the eigenvectors of the covariance matrix of the point cloud data. Figure 9(3) illustrates the normal direction N. Finally, the pseudo-target is converted into planar data based on the centroid and normal (step S134: first planarization step). Figure 9(4) illustrates the planarized plane P.

[0043] Next, as shown in Figure 6, first surface feature data is extracted from the pseudo-target (Step S14: First surface feature data extraction step). The first surface feature data extraction step calculates cross-sectional data from the relative positions of each point with respect to the normal direction N for at least one cross-section in the planar data. Figure 10 is an illustrative diagram of a pseudo-target and surface feature data, where (1) shows the state in which a pseudo-target is set and (2) shows the state in which surface feature data is extracted. Figure 11 is an illustrative diagram of the extraction of surface feature data. As shown in Figure 10(1), after a pseudo-target 6 is set for the point cloud 80 acquired with respect to the wall surface 72, multiple cross-sectional data are calculated as shown in Figure 11. Specifically, in this embodiment, cross-sectional data for the vertical cross-section (C1~C3) and cross-sectional data for the horizontal cross-section (C4~C6) of the pseudo-target 6 are calculated. The vertical direction corresponds to, for example, the circumferential direction of the tunnel cross-section perpendicular to the tunnel axis direction, and the horizontal direction corresponds to, for example, the tunnel axis direction. In addition, unlike this embodiment, more cross-sectional data may be calculated, or, for example, cross-sectional data in an oblique direction may be calculated. The oblique direction is a direction inclined with respect to the vertical and horizontal directions.

[0044] The position coordinates of the pseudo-target and the first surface feature data are registered in association with the position coordinates of the 3D laser scanner (Step S15: First registration step). In this embodiment, the absolute coordinates of the point cloud 8 in the pseudo-target (6a~6e), as well as the centroid, normal, and plane data are registered as the position coordinates of the pseudo-target based on the absolute coordinates of the LiDAR scanner 4, but it is also acceptable to register only the absolute coordinates of the point cloud 8 in the pseudo-target (6a~6e).

[0045] (Regarding the displacement measurement step) Figure 12 shows the displacement measurement flowchart. As shown in Figure 12, a 3D laser scanner is used to acquire point cloud data of the sprayed concrete surface inside the tunnel (Step S21: Second point cloud data acquisition step). The method for acquiring the point cloud data is the same as in the first point cloud data acquisition step (Step S12).

[0046] Next, the search area for the pseudo-target is set from the second point cloud data (Step S22: Search Area Setting Step). Figure 13 shows the search area setting flowchart. Figure 14 is an image diagram of the setting of the search area 60, and Figures 15 and 16 are explanatory diagrams of the pseudo-target identification step. As shown in Figures 13 and 14(1), first, based on the position coordinates 60a of the pseudo-target 6 set in the initial measurement step (step S01), the search area 60 is extracted from the second point cloud data 81 (step S221: second extraction step). Next, second point cloud data included in the extracted search area 60 is acquired, and the centroid of the search area 60 is calculated based on this point cloud data (Step S222: Second centroid calculation step). Based on the second point cloud data and the centroid, the normal vector of the search area 60 is calculated (Step S223: Second normal vector calculation step). Based on the centroid and normal vector, the search area 60 is converted into planar data (Step S224: Second planarization step). This process is the same as the pseudo-target setting step shown in Figure 9. Figure 14(1) shows the state in which the search area 60 has been converted into planar data.

[0047] Second surface feature data is extracted from the search area 60 (Step S23: Second surface feature data extraction step). In Figure 14(2), the second surface feature data is the vertical cross-section of the search area 60 (C 10 ~C 30 Cross-sectional data relating to ) and the cross-sectional data in the lateral direction (C 40 ~C 60 This image shows the calculated cross-sectional data for ). As shown in Figure 15(1), the pseudo-target 6 is extracted and its position coordinates are identified based on the similarity between the first surface feature data of the pseudo-target 6 and the second surface feature data of the search area 60 (Step S24: Pseudo-target identification step). In this embodiment, as shown in Figure 15(2), cross section C 20 Section C2 and section C 50The similarity to section C5 is determined to be high. Then, as shown in Figure 16(1), a pseudo-target 6 is extracted from the position coordinates of the region with high similarity (A1, A2), as shown in Figure 16(2), and its position coordinates are identified and obtained.

[0048] The position coordinates of the identified pseudo-target are registered in association with the position coordinates of the 3D laser scanner (Step S25: Second registration step). The method of registering the coordinates is the same as in the first registration step (Step S15). In the initial measurement step, the displacement of the tunnel's internal shape is measured based on the changes in the position coordinates of the simulated target set in the initial measurement step and the simulated target identified in the displacement measurement step (Step S26: Displacement Measurement Step). In this embodiment, the point cloud data 8 of the simulated target 6 shown in Figure 10(2) and the point cloud data 8a of the simulated target 6 shown in Figure 16(2) are compared to measure the displacement of the tunnel's internal shape. [Examples]

[0049] In Example 1, point cloud data (pseudo-target) is used as the target, but the target may be a physical target component installed in the tunnel. However, to eliminate the effort of installing the target, it is preferable that the target is not dedicated to measurement but is used in conjunction with other components normally installed in the tunnel.

[0050] One example is the tunnel displacement measurement system 10 shown in Figure 17, which includes a tunnel displacement measurement device 1a and targets 100a to 100e installed on the wall surface 72 of the tunnel 7. Targets 100a to 100e are positioned in the same locations as the pseudo-targets 6a to 6e in each range R1, R2, R3, ... (see Figure 5) in the tunnel axial direction.

[0051] Furthermore, LiDAR scanners 4 are installed along each of the walls 72 on both sides in the tunnel width direction, and each LiDAR scanner 4 is configured to acquire point cloud data of the wall 72 on the opposite side from the wall 72 on which the LiDAR scanner 4 is installed. The tunnel width direction is perpendicular to the tunnel axis direction in the plane. In Figure 17, the tunnel axis direction is indicated by symbol A, and the tunnel width direction is indicated by symbol B.

[0052] Figure 18 shows a target 100a. In this example, a plate for anchoring a rock bolt 200 is used as the target 100a. The rock bolt 200 is driven from the wall surface 72 of the tunnel 7 into the ground behind it, and the plate of the rock bolt 200, i.e., the target 100a, is placed on the sprayed concrete of the wall surface 72. The end of the rock bolt 200 protrudes from the outer surface of the target 100a, and a nut 210 is tightened onto the protruding portion.

[0053] The target 100a is a plate-shaped metal object with projections 110 on its outer surface, giving it an uneven surface. The projections 110 are continuous horizontal ridges and are arranged in multiple parallel rows vertically within the target 100a. More specifically, there are two rows of projections 110 on the upper part of the target 100a and two rows on the lower part, with through holes (not shown) provided between the two upper rows of projections 110 and the two lower rows of projections 110 for passing the ends of the lock bolts 200. The cross-section of the projections 110 is rectangular.

[0054] Other targets 100b to 100e have similar shapes, and only their orientation during installation differs depending on the position of targets 100b to 100e. For example, target 100c at the top of tunnel 7 shown in Figure 17 is positioned with its plate surface horizontal, and the protrusions 110 are arranged in multiple rows with spacing in the tunnel width direction B.

[0055] The functions of the tunnel displacement measuring device 1a are substantially the same as those of the tunnel displacement measuring device 1 described in Figure 1. However, the initial measurement unit 2 replaces the pseudo-target setting means 23 and the first surface feature data extraction means 24 with a surface feature data acquisition means that accepts input of imaging results obtained in advance from an imaging device (not shown) such as a 3D laser scanner that has imaged targets 100a to 100e, and acquires surface feature data, which is the surface feature of targets 100a to 100e, from the imaging results, and a position information acquisition means that extracts surface feature data from point cloud data (first point cloud data) of the inner surface of the tunnel including targets 100a to 100e and acquires the position coordinates (first position information) of targets 100a to 100e. The first registration means 26 registers the position information and surface feature data of targets 100a to 100e.

[0056] The displacement measurement unit 3 is equipped with a target identification means that replaces the pseudo-target identification means 35 with a target identification means that extracts targets 100a to 100e from the point cloud data of the search area 60 and identifies and acquires their position coordinates (second position information). The search area 60 includes the positions identified based on the position coordinates (first position information) of targets 100a to 100e acquired by the position information acquisition means, and is a larger area than targets 100a to 100e. The second registration means 36 and the displacement measurement means 37 perform the same processing on the area of ​​targets 100a to 100e instead of the pseudo-targets 6a to 6e.

[0057] The initial measurement flow in this embodiment is basically the same as that described in Figure 6. However, in this embodiment, the computer 5 receives input of imaging results obtained in advance by a separate imaging device (not shown) that has imaged targets 100a to 100e, and obtains surface feature data of targets 100a to 100e from these imaging results. Then, instead of performing the pseudo-target setting step (step S13) and the first surface feature data extraction step (step S14) in Figure 6, surface feature data is extracted from the point cloud data of the tunnel interior including targets 100a to 100e obtained in the first point cloud data acquisition step (step S12) to obtain the position coordinates of targets 100a to 100e. Subsequently, in the first registration step (step S15), the position coordinates and surface feature data of targets 100a to 100e are registered.

[0058] Since targets 100a to 100e have protrusions 110, in the surface feature data, the surface features caused by the protrusions 110 appear in the vertical cross-sectional data (see symbols C1 to C3 in Figure 11), and the height of the protrusions and the depth of the recesses in the surface features appear in the horizontal cross-sectional data (see symbols C4 to C6 in Figure 11). It is desirable that the cross-section used to calculate the cross-sectional data is a position that avoids the ends of the lock bolts 200 and the nuts 210.

[0059] The displacement measurement flow in this embodiment is basically the same as the flow in Embodiment 1 described in Figure 12. However, in this embodiment, instead of identifying the position coordinates of pseudo-targets 6a to 6e in the pseudo-target identification step (step S24) in Figure 12, the position coordinates (second position information) of the region of targets 100a to 100e are identified. The target identification step extracts regions with similar surface features to those of targets 100a to 100e from the search region 60 based on the similarity between the surface feature data of targets 100a to 100e and the surface feature data of the search region 60, and identifies and obtains their position coordinates. The target identification method may be the same template matching method as in Embodiment 1, but is not limited to this; any method that searches for surface features due to the protrusions 110 of targets 100a to 100e within the search region 60 is acceptable, and various known matching methods can be applied.

[0060] In this embodiment, the placement of protrusions 110 on the outer surfaces of targets 100a to 100e allows for accurate positioning of these targets 100a to 100e. That is, as shown in Figure 19(a), if the surface of target 100a is smooth, there are no characteristic features in the shape of the tunnel interior including target 100a. However, as in this embodiment, when protrusions 110 are formed on target 100a, characteristic irregularities are created in the shape of the tunnel interior including target 100a, as shown by the solid line portion in Figure 19(b).

[0061] In this embodiment, depending on the location of these uneven features, even when horizontal or vertical displacement occurs in the tunnel width direction B of the target 100a, the horizontal position along the tunnel width direction B and the vertical position of the target 100a after displacement can be accurately determined. In particular, if the surface of the target 100a is smooth as shown in Figure 19(a), even if vertical displacement occurs in the target 100a as shown by arrow a (the shape of the tunnel interior is shown by a dashed line), the shape of the tunnel interior hardly changes, making it difficult to determine the vertical displacement of the target 100a. In contrast, in this embodiment, as shown in Figure 19(b), when vertical displacement occurs in the target 100a, the position of the protrusion 110 clearly changes, thereby making it easy to determine the vertical displacement of the target 100a.

[0062] In the second registration step (step S25) and the displacement measurement step (step S26), processing is performed on the region of targets 100a to 100e instead of the pseudo-targets 6a to 6e. In the displacement measurement step (step S26), the vertical displacement of targets 100a to 100e, i.e., the vertical displacement of the tunnel interior, and the horizontal displacement of targets 100a to 100e in the tunnel width direction B, i.e., the horizontal displacement of the tunnel interior in the tunnel width direction B, can be determined by the change in the position coordinates of targets 100a to 100e (the difference between the first position information and the second position information).

[0063] Thus, in this embodiment, by using the plates of the rock bolt 200 as targets 100a to 100e, there is no need to set up a dedicated target for measurement in the tunnel 7, and the measurement work is simplified. Furthermore, while some dedicated targets for measurement use prisms or mirrors, there is a risk of damage when the tunnel cutting edge is blasted, but there is no such risk with targets 100a to 100e (plates of the rock bolt 200) in this embodiment. Therefore, it is possible to measure the displacement immediately after blasting.

[0064] Furthermore, since target 100a has protrusions 110 on its outer surface, the horizontal and vertical displacements of the area of ​​target 100a can be measured accurately due to the uneven surface characteristics of the protrusions 110. This is also true for the other targets 100b to 100e.

[0065] Furthermore, since the projection 110 of target 100a is a ridge provided in multiple vertical stages, the vertical displacement of the area of ​​target 100a can be measured with greater accuracy. The same applies to targets 100b, 100d, and 100e. For target 100c at the top of tunnel 7, the projection 110 is arranged in multiple rows with spacing in the tunnel width direction B, allowing for accurate measurement of the horizontal displacement in the tunnel width direction B.

[0066] The shapes of targets 100a to 100e are not limited to those shown above. For example, in the example shown in Figure 18, targets 100a to 100e are provided with four protrusions 110, but the number of protrusions 110 is not particularly limited. However, in terms of accuracy, it is desirable to provide three or more protrusions 110. Figure 20 briefly shows the results of the inventor's study on the relationship between the number of protrusions 110 and measurement accuracy, showing that accuracy improves when there are three or more protrusions 110.

[0067] Furthermore, while the projections 110 of target 100a are continuous ridges in the horizontal direction, the shape of the projections 110 is not limited to this. For example, as shown in target 100a' in Figure 21, which is a modified example of target 100a in Figure 18, the projections 110' may be provided discretely with intervals between them in the vertical, horizontal, and vertical directions. However, having the projections 110 continuous in the horizontal direction makes it easier to set up a large number of vertical cross-sectional data (see symbols C1 to C3 in Figure 11) where the unevenness characteristics caused by the projections 110 appear, thus contributing to improved measurement accuracy.

[0068] Furthermore, the protrusion height of the projection 110 can also be freely set, taking into consideration measurement accuracy and ease of handling. For example, increasing the protrusion height of the projection 110 improves measurement accuracy because the convex part becomes clearly visible in the cross-sectional data, but the projection 110 is more easily damaged during transportation and installation, and the space required for storage also increases. It is also possible to change the pattern of the projection 110 for each individual target, making it possible to identify each target by its pattern.

[0069] In addition, as a physical target member installed inside the tunnel, steel support structures 73 (metal fittings) installed inside the tunnel may be used as targets 100a to 100e, as shown in the tunnel displacement measurement system 10a in Figure 22. In this case, as described above, a plate having a projection 110 may be attached to the outer surface of the steel support structure 73, and the range including the projection 110 of the steel support structure 73 may be designated as targets 100a to 100e. Alternatively, as shown in target 100a in Figure 23, the steel support structure 73 itself may have the projection 110 on its outer surface. In this case as well, target 100a is the range including the projection 110 of the steel support structure 73.

[0070] In this embodiment, five targets 100a to 100e are provided for each measurement cross-section 70, but some may be omitted. For example, targets 100b to 100d can be omitted, and measurements can be performed at these locations using the aforementioned pseudo-targets. [Examples]

[0071] While the tunnel displacement measurement systems 10 and 10a in Example 2 are primarily for mountain tunnels and tunnels constructed using the NATM method, this example will instead describe an example of a displacement measurement system primarily for urban tunnels and tunnels constructed using the shield tunneling method.

[0072] In shield tunneling, the displacement of the interior space of the lining constructed by the shield method is measured. In a shield tunnel, a segment ring, which is a ring-shaped lining, is constructed in the excavated shaft drilled by a shield boring machine by connecting multiple segment pieces in the circumferential direction of the tunnel. Furthermore, multiple segment rings are connected in the axial direction of the tunnel to construct a lining inside the excavated shaft.

[0073] The segment piece has its outer surface (outside) facing the ground of the excavation shaft and its inner surface (inside) facing the interior cavity. The tunnel displacement measurement system in this embodiment measures the displacement of the inner surface of the segment using the tunnel displacement measurement device 1a described above. On the inner surface (inside) of the segment piece, a gripping fitting is exposed and provided, which is gripped by the erector of the shield excavator. The gripping fitting has a recess into which the gripping tool of the erector is inserted, and the gripping fitting itself has an uneven surface. In addition, adjacent segment pieces in the circumferential direction of the tunnel are connected in the circumferential direction of the tunnel via joint fittings provided at the circumferential ends of each segment piece. The joint fittings also have an uneven surface, similar to the gripping fittings.

[0074] In this embodiment, instead of the anchoring plates for rock bolts or steel supports that serve as target members on the inner surface of the tunnel in Embodiment 2, metal fittings that are provided on the inner surface of the segment pieces and have uneven surfaces on the portion exposed to the inner surface of the tunnel are used as target members. For example, gripping fittings and joint fittings are metal fittings that have uneven surfaces. In Embodiment 3 as well, the displacement of the tunnel interior can be measured using the same tunnel displacement measurement method as in Embodiment 2. [Industrial applicability]

[0075] This invention is useful for measuring tunnel top settlement and displacement of the internal shape during tunnel excavation work. [Explanation of Symbols]

[0076] 1,1a Tunnel displacement measuring device 2 Initial measurement unit 3. Displacement Measurement Unit 4 LiDAR scanners 5 Computers 6,6a~6e Simulated targets 7 Tunnel 8,8a,80,81 Point cloud (data) 10,10a Tunnel Displacement Measurement System 11a~11c Scanner Targets 21. Scanner coordinate acquisition means 22 First point cloud data acquisition means 23. Pseudo-target setting means 24. First means for extracting surface feature data 26. First registration method 32 Second point cloud data acquisition means 33 Search area setting means 34. Second means for extracting surface feature data 35. Means for identifying a pseudo-target 36. Second registration method 37 Displacement measuring means 60 Search area 70 Measurement cross section 71. Face of the pipe 72 Wall surface Targets: 100a~100e, 100a' 110,110' protrusion

Claims

1. A method for measuring the displacement of the tunnel interior over time by setting a target using point cloud data of the tunnel interior acquired by a 3D laser scanner, It consists of an initial measurement step and a displacement measurement step. 1) The initial measurement step is: To obtain surface feature data, which is the surface feature of the aforementioned target, The process includes obtaining the aforementioned surface feature data, then obtaining point cloud data of the inner surface of the tunnel including the target as first point cloud data, and extracting the surface feature data from the first point cloud data to obtain first position information of the target, Step 2), which is performed after a predetermined time has elapsed, is: A search area is set that includes a location identified based on the first location information and has a larger area than the target; point cloud data of the search area is acquired; the target is extracted from the point cloud data of the search area based on the surface feature data to acquire second location information. This includes measuring the displacement of the space inside the tunnel based on the difference between the first position information and the second position information, The tunnel displacement measurement method is characterized in that the target is a physical target member installed on the wall surface inside the tunnel.

2. The tunnel displacement measurement method according to Claim 1, characterized in that the target member has protrusions and the shape of the tunnel cavity has the aforementioned uneven features.

3. The tunnel displacement measurement method according to Claim 1, characterized in that the target member is a metal fitting installed on the wall surface inside the tunnel and having the aforementioned uneven surface characteristics in the portion exposed to the inner surface.

4. A three-dimensional laser scanner is provided along each of the walls on both sides in the tunnel width direction, The tunnel displacement measurement method according to any one of claims 1 to 3, characterized in that each three-dimensional laser scanner acquires point cloud data of the wall surface opposite to the wall surface on which the three-dimensional laser scanner is installed.

5. The tunnel displacement measurement method according to any one of claims 1 to 3, characterized in that the surface feature data, which is the surface feature of the target, is obtained from the imaging result obtained by imaging the target with an imaging device.

6. The tunnel displacement measurement method according to claim 2, characterized in that the target has protrusions provided at multiple locations in the vertical direction.

7. A device that uses point cloud data of the tunnel's inner surface acquired by a 3D laser scanner to set a target and measure the displacement of the tunnel's interior space over time, A physical target member is installed as the target on the wall surface inside the tunnel, A tunnel displacement measurement system having, The aforementioned device is It consists of an initial measurement unit and a displacement measurement unit. The initial measurement unit is, To obtain surface feature data, which is the surface feature of the aforementioned target, After acquiring the aforementioned surface feature data, point cloud data of the inner surface of the tunnel including the target is acquired as first point cloud data, and the surface feature data is extracted from the first point cloud data to obtain the first position information of the target. The displacement measurement unit is, A search area is set that includes a location identified based on the first location information and has a larger area than the target; point cloud data of the search area is acquired; the target is extracted from the point cloud data of the search area based on the surface feature data to acquire second location information. A tunnel displacement measurement system characterized by measuring the displacement of the space inside a tunnel based on the difference between the first position information and the second position information.

8. The tunnel displacement measurement system according to claim 7, characterized in that the target member has protrusions and the shape of the tunnel cavity has the aforementioned uneven features.

9. The tunnel displacement measurement system according to claim 7, characterized in that the target member is a metal fitting installed on the wall surface inside the tunnel and having the aforementioned uneven surface characteristics in the portion exposed to the inner surface.

Citation Information

Patent Citations

  • Survey method and program for displaying result thereof

    JP2004170164A

  • System, method and program for measuring displacement information

    JP2006349579A

  • Inner space displacement measurement method and inner space displacement measurement system

    JP2012058167A

  • Method and system for measuring displacement of fluctuation surface

    JP2013238549A

  • Method of measuring displacement of space within tunnel

    JP2014002027A