CN106969749A - A kind of detection method of deformation of cross section of subway tunnel - Google Patents
A kind of detection method of deformation of cross section of subway tunnel Download PDFInfo
- Publication number
- CN106969749A CN106969749A CN201710286331.XA CN201710286331A CN106969749A CN 106969749 A CN106969749 A CN 106969749A CN 201710286331 A CN201710286331 A CN 201710286331A CN 106969749 A CN106969749 A CN 106969749A
- Authority
- CN
- China
- Prior art keywords
- section
- cross
- tunnel
- point
- deformation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C7/00—Tracing profiles
- G01C7/06—Tracing profiles of cavities, e.g. tunnels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
The present invention relates to a kind of detection method of deformation of cross section of subway tunnel, methods described comprises the following steps:Using tunnel track center line as origin, the rectangular coordinate system of subway tunnel section is set up;DATA REASONING is carried out to tunnel cross-section according to goniometer, laser range finder and laser scanner, coordinate of all measurement points of coordinate offset angle α and subway tunnel section under rectangular coordinate system is obtained, the wherein coordinate of ith measurement point is (xi,yi);Basic measurement data is carried out curve fitting by ellipse fitting algorithm, the form parameter of the tunnel cross-section of coincidence loss requirement is obtained;Judge whether the form parameter of tunnel cross-section and the deviation of history form parameter exceed deviation threshold, if then showing that obvious deformation occurs for subway tunnel section, obvious deformation does not occur for subway tunnel section if otherwise showing.Compared with prior art, the present invention has the advantages that measuring accuracy is high, measurement is comprehensive and is easy to implement.
Description
Technical field
The present invention relates to technical field of civil engineering, more particularly, to a kind of detection method of deformation of cross section of subway tunnel.
Background technology
At present, urban rail transit in China is quickly grown, and subway has been opened in existing more than ten of city, and subway turns into citizen
The important way of trip, ensures that subway circulation safety responsibility is great.
Safety problem during shield tunnel construction and metro operation is mainly manifested in the following aspects:Crouched under tunnel soft
Soil layer softens under Long-term Vibration load action causes sedimentation;Tunnel adjacent building construction activities cause tunnel deformation;Interval tunnel
Sleeping soil layer soil erosion damaging property linear deformation under road;When tunnel is by different soil property situations area, soil body physics
Matter difference causes tunnel structure to produce differential settlement under action of long-term load;Shield Tunnel in Soft Soil differential settlement causes pipe
Piece seam is opened, subgrade of metro and section of jurisdiction are peeled off and ftractureed, and portion percolating water occurs, emits mud and local failure, influences tunnel
Normal use and metro operation.
The monitoring technology to the constructing tunnel stage domestic at present comparative maturity, but the monitoring weight during runing tunnel
Visual range degree is far from enough.In fact, operation stage is because time span is big, influence factor is complicated, disaster social influence big, tunnel
Health monitoring should more be paid attention to.Tunnel health monitoring includes tunnel structure corrosion monitoring, structural deformation monitoring, structural internal force measure
With ambient conditions monitoring, wherein especially structural deformation monitoring is extremely important, its Contents for Monitoring be mainly tunnel vertical sedimentation,
Horizontal displacement and the convergent deformation of section.
For the detection of the convergent deformation of tunnel cross-section, current existing several method includes artificial section survey, three-dimensional
Laser scanning measurement and the measurement of Tunnel testing dolly, not only precision is low but also spends the time long for artificial detection, 3 D laser scanning
Although mensuration has reached high-acruracy survey, but the requirement to measuring instrument is high, is not easy to realize and cost is too high, pass through tunnel
It is to detect crack and infiltration situation inside tunnel cross-section by image recognition that the detection dolly method that measures in road is mostly at present,
At present temporarily without relatively good method, the detection parameter obtained using Tunnel testing dolly realizes the deformation to tunnel cross-section
The judgement of situation.
The content of the invention
The purpose of the present invention is to provide a kind of detection method of deformation of cross section of subway tunnel regarding to the issue above.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of detection method of deformation of cross section of subway tunnel, methods described comprises the following steps:
1) using tunnel track center line as origin, the rectangular coordinate system of subway tunnel section is set up;
2) DATA REASONING is carried out to tunnel cross-section according to goniometer, laser range finder and laser scanner, obtains coordinate inclined
All measurement points of angle α and subway tunnel section are moved in step 1) coordinate under the rectangular coordinate system set up, wherein ith measurement
The coordinate of point is (xi,yi);
3) by ellipse fitting algorithm to step 2) obtained basic measurement data carries out curve fitting, and obtains coincidence loss
It is required that tunnel cross-section form parameter;
4) judgment step 3) whether the deviation of the obtained form parameter of tunnel cross-section and history form parameter exceed deviation
Threshold value, if then showing that obvious deformation occurs for subway tunnel section, if otherwise showing, obvious deformation does not occur for subway tunnel section.
The step 3) be specially:
31) according to all measurement points in step 1) coordinate under the rectangular coordinate system set up, carry out basic fitting and obtain tunnel
The initial elliptic equation of road section;
32) all measurement points are asked in step 31) orthogonal reference point on obtained initial elliptic equation;
33) according to step 32) obtained orthogonal reference point, the elliptic equation after being optimized by iterative;
34) fitting parameter of the elliptic equation after optimization is asked for, judges whether coincidence loss is required fitting parameter, if then
Using the elliptic equation after optimization as tunnel cross-section form parameter, if otherwise return to step 33).
The step 32) be specially:The point using measurement point as primary iteration, by Newton iteration method to orthogonal correlated condition
Solve the coordinate for obtaining orthogonal reference point, the orthogonal correlated condition is specially:
f2(x', y')=b2x'(yi-y')-a2y'(xi-x')2=0
Wherein, a and b are the basic parameter of initial elliptic equation, and (x', y') is the coordinate of orthogonal reference point, (xi,yi)
For the coordinate of measurement point.
The step 33) be specially:
331) coordinate X' of the orthogonal reference point under natural system of coordinates is obtained by coordinate transform;
332) using step 331) orthogonal reference point under obtained natural system of coordinates, according to minimizing geometric distance condition,
The point coordinates X on the elliptic equation after optimization is asked for by iteration;
333) according to step 332) point coordinates on elliptic equation after obtained optimization, carry out ellipse fitting and optimized
Elliptic equation afterwards.
The coordinate transform is specially:
X'=R-1x'+Xc
Wherein, x' is coordinate of the orthogonal reference point under rectangular coordinate system, XcSat naturally for the origin of rectangular coordinate system
Coordinate under mark system.
The minimizing geometric distance condition is specially:
MinE=min { (X-X')T(X-X')}。
The fitting parameter includes fitting precision σ and error of fitting Δ ε.
The step 3) also include rejecting the abnormity point in curve fitting process, be specially:
351) according to the curvilinear equation obtained after curve matching, residual error of each measurement point relative to curvilinear equation is calculated;
352) by step 351) in the corresponding residual error of obtained all measurement points constitute residual vector, pass through sample median
The method of inspection finds the exceptional value in residual vector;
353) by step 352) in the corresponding measurement point of obtained exceptional value reject, re-start curve matching, accorded with
Close the form parameter of the tunnel cross-section of error requirements.
All measurement points of the subway tunnel section include laser range finder and scan obtained high-acruracy survey point and swash
The low precision measure point that photoscanner is obtained.
It is specially with time interval t, the time interval t between the high-acruracy survey point and low precision measure point:
Wherein, c is the light velocity, and x is the distance between laser range finder and laser scanner, and V is deformation of cross section of subway tunnel
Detection means pace.
Compared with prior art, the invention has the advantages that:
(1), will be using inclinometer, laser range finder and laser scanner to tunnel cross-section measurement by ellipse fitting algorithm
Data carry out calculating integration, by the shape of tunnel cross-section by digitized representations, so as to by the shape amount of tunnel cross-section
Change, the size of deformation is determined by being compared with historical data, the real-time monitoring to tunnel cross-section deformation situation is realized.
(2) ellipse fitting algorithm proposed by the present invention, on the basis of original ellipse fitting, by asking for orthogonal reference point
The fitting result being initially obtained is optimized, so as to obtain optimal ellipse fitting result, that is, actual conditions are best suited
Tunnel cross-section form parameter.
(3) during orthogonal reference point is asked for, the elliptic equation after being optimized by iterative, by repeatedly changing
In generation, the precision for the elliptic equation asked for is further increased, so as to improve the precision judged finally for tunnel cross-section deformation.
(4) during being optimized to elliptic equation, the fitting precision and error of fitting of elliptic equation have also been asked for
,, can be with when the error of fitting result is excessive so as to there is the judge of quantization to final fitting result as fitting parameter
Reduce the error of fitting result by increasing iterations, improve fitting precision.
(5) after curve matching, also the larger point of error in measurement point is rejected by recalculating, and rejected
After re-start fitting, further improve fitting precision.
(6) when recording measurement point, it is contemplated that caused due to the distance between laser range finder and laser scanner
Measurement point by introducing time interval is carried out low precision measure point and high-acruracy survey point not in conplane situation
It is unified, it is ensured that the point of laser range finder and laser scanner measurement is the measurement point under same plane.
Brief description of the drawings
Fig. 1 is flow chart of the method for the present invention.
Embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.The present embodiment is with technical solution of the present invention
Premised on implemented, give detailed embodiment and specific operating process, but protection scope of the present invention is not limited to
Following embodiments.
As shown in figure 1, the present embodiment proposes a kind of detection method of deformation of cross section of subway tunnel, comprise the following steps:
1) using tunnel track center line as origin, the rectangular coordinate system of subway tunnel section is set up;
2) DATA REASONING is carried out to tunnel cross-section according to goniometer, laser range finder and laser scanner, obtains coordinate inclined
All measurement points of angle α and subway tunnel section are moved in step 1) coordinate under the rectangular coordinate system set up, wherein ith measurement
The coordinate of point is (xi,yi), all measurement points of subway tunnel section include laser range finder and scan obtained high-acruracy survey point
The low precision measure point obtained with laser scanner.
With time interval t between high-acruracy survey point and low precision measure point, time interval t is specially:
Wherein, c is the light velocity, and x is the distance between laser range finder and laser scanner, and V is deformation of cross section of subway tunnel
Detection means pace;
3) by ellipse fitting algorithm to step 2) obtained basic measurement data carries out curve fitting, and obtains coincidence loss
It is required that tunnel cross-section form parameter:
31) according to all measurement points in step 1) coordinate under the rectangular coordinate system set up, carry out basic fitting and obtain tunnel
The initial elliptic equation of road section;
32) all measurement points are asked in step 31) orthogonal reference point on obtained initial elliptic equation, be specially:With
Measurement point carries out solving the seat for obtaining orthogonal reference point by Newton iteration method as primary iteration point to orthogonal correlated condition
Mark, orthogonal correlated condition is specially:
f2(x', y')=b2x'(yi-y')-a2y'(xi-x')2=0
Wherein, a and b are the basic parameter of initial elliptic equation, and (x', y') is the coordinate of orthogonal reference point, (xi,yi)
For the coordinate of measurement point;
33) according to step 32) obtained orthogonal reference point, the elliptic equation after being optimized by iterative:
331) coordinate X' of the orthogonal reference point under natural system of coordinates is obtained by coordinate transform, coordinate transform is specially:
X'=R-1x'+Xc
Wherein, x' is coordinate of the orthogonal reference point under rectangular coordinate system, XcSat naturally for the origin of rectangular coordinate system
Coordinate under mark system;
332) using step 331) orthogonal reference point under obtained natural system of coordinates, according to minimizing geometric distance condition,
The point coordinates X on the elliptic equation after optimization is asked for by iteration, wherein, minimizing geometric distance condition is specially:
MinE=min { (X-X')T(X-X')};
333) according to step 332) point coordinates on elliptic equation after obtained optimization, carry out ellipse fitting and optimized
Elliptic equation afterwards;
34) fitting parameter (including fitting precision σ and error of fitting Δ ε) of the elliptic equation after optimization is asked for, judges to intend
Close parameter whether coincidence loss requirement, if then using the elliptic equation after optimization as tunnel cross-section form parameter, if otherwise
Return to step 33);
Step 3) also include rejecting the abnormity point in curve fitting process, be specially:
351) according to the curvilinear equation obtained after curve matching, residual error of each measurement point relative to curvilinear equation is calculated;
352) by step 351) in the corresponding residual error of obtained all measurement points constitute residual vector, pass through sample median
The method of inspection finds the exceptional value in residual vector;
353) by step 352) in the corresponding measurement point of obtained exceptional value reject, re-start curve matching, accorded with
Close the form parameter of the tunnel cross-section of error requirements;
4) judgment step 3) whether the deviation of the obtained form parameter of tunnel cross-section and history form parameter exceed deviation
Threshold value, if then showing that obvious deformation occurs for subway tunnel section, if otherwise showing, obvious deformation does not occur for subway tunnel section.
According to above-mentioned steps, deformation inspection is carried out to tunnel cross-section with the railcar equipped with laser range finder and laser scanner
Exemplified by survey, detailed process is as follows:
If the positional information done in the tunnel that laser range finder and laser scanner are obtained can be calculated as coordinate value,
It is that the relative position information with laser range finder and laser scanner is put in tunnel by these values, and laser range finder and laser
The relative position (the namely position of tunnel track center line) of scanner and detection car is constant, is then put in tunnel relative to tunnel
The relative position of track centre can convert and obtain, and with reference to the data of goniometer, finally calculate and obtain coordinate value.This coordinate value
It is in a rectangular coordinate system using tunnel track center line as origin.By ellipse fitting algorithm, tunnel can be obtained and broken
The SECTION EQUATION in face.The data of vertically and horizontally convergence monitoring point are tried to achieve by equation.
Ellipse fitting algorithm
At present, the conventional method of ellipse fitting is the minimum value that minimum seeks object function by least square method, so that
Determine the elliptical shape parameter at section of jurisdiction interface.The method of ellipse fitting uses geometric distance method.
The definition of geometric distance is point to oval beeline.The most short point of distance to fitting sampling point Xi on ellipse claims
For point X ' i nearest reference point, orthogonal reference point is also.Then just can be with table for giving the geometric distance fitting problems of n point
Geometric distance minimization problem is shown as, is expressed as follows with mathematical expression:
MinE=min { (X-X')T(X-X') } wherein X and X' be expressed as give n sample point and its it is corresponding just
Hand over the column vector of related point coordinates, i.e. X=[X1,X2,…,Xn], X '=[X1‘,X‘2,…,X‘n].Geometric distance is intended as can be seen here
Close ellipse to first have to calculate sample point X orthogonal reference points in initial fitted ellipse, acquisition optimization is then being solved to above formula
Elliptic parameter.
Reference axis xoy is introduced, the origin of the coordinate system is located at (Xc,Yc), and around point selection α angles, then the coordinate introduced
It is that variation relation between xoy and former coordinate system XOY is as follows:
X=R (X-Xc)
X=R-1x+Xc
Wherein,
Under xoy coordinate systems, elliptic equation can be expressed asKnow by orthogonal correlation
The oval tangent line of point is vertical with the line of sample point with orthogonal reference point, therefore can obtain two following equations:
f2(x, y)=b2x(yi-y)-a2y(xi-x)2=0
2 equations of the above are also referred to as orthogonal correlated condition.
Orthogonal correlated condition can be solved with Newton iteration method:
Set up Newton iteration:
xk+1=xk+Δx
The initial value x of iterationoSampled point x can be usedi, through overtesting, enough accuracy is just can obtain by 2-3 iteration
Orthogonal reference point.But when to point xiIn elliptical center, or oval two axles it is equal when, above-mentioned matrix can be singular matrix,
Now, the orthogonal reference point on ellipse is not unique.
Solve above formula, it is therefore desirable toThat is 2J (X-X ')=0, wherein J is Jacobian matrix:By
Gauss regular equation, can be obtained,
J|kΔ α=(X-X ') |k
αk+1=αk+Δα
Above iteration only needs to just can obtain for 2-4 times the oval parameter of enough accuracy.
The rejecting of abnormity point
When electronic equipment, cable or other objects for being installed in laser range finder and laser scanner scans to tunnel inner wall
When, it may appear that measurement abnormity point.When there is abnormity point in sampling point, fitting algorithm all can be by than large effect.Therefore such as
What removes identification and rejecting abnormalities point is critically important.In tunnel fitting, abnormity point takes the sampling point on tunnel duct piece.
The method being first fitted according to prior elliptical, obtains a set of optimal coefficient, and calculates the residual error each put with this coefficient,
Residual error a little constitute a residual vector.Because exceptional value here is usual unnecessary one, so using sample middle position
The number method of inspection, finds exceptional value from residual vector, and then the rough error point corresponding to generation exceptional value is rejected.
Wherein sample median method of inspection principle is as follows:
If x1,x2,…,xnIt is X independent same distribution sample, and x(1),x(2),…,x(n)For its order statistic, sample is used
This median med { XiEstimation central overall position, X(m)med{XiAnd med { Xi}X1It is then both sides extremum and center up and down
The difference of position, by X(i)-med{XiTake median as yardstick again, obtain the test statistics on test of outlier.
Or
Wherein
TmAnd Tm *It is the upper and lower side test of outlier statistic of sample median method of inspection, DT respectivelymIt is sample median
The bilateral test of outlier statistic of method of inspection.Wherein TmAnd Tm *It is respectively used to examine X(m)And X(1)Whether exceptional value.
The advantage of the method is available with ellipse fitting algorithm and finds a set of fitting coefficient with robustness, Ran Houli
It can be found out after exceptional value, rejecting abnormalities value with sample median, recycle ellipse fitting algorithm to be fitted, obtain optimal plan
Close effect.
Fitting algorithm precision
The precision of cross section correct can be represented according to the standard deviation in error theory analysis method:
Wherein viRefer to residual vector, n counts for measurement.
Fitting algorithm error
Oval form error computational methods are as follows:
It is respectively that X, Y-axis set up coordinate system with oval long and short axle using the center of least square ellipse as origin.Measurement
Point pt(xt, yt) deviate least square ellipse deviation delta riFor:
Wherein:ai=tan-1(yi/xi)
Elliptical shape error is:
Δ ε=max { Δ ri, i=1,2, L, n } and-min { Δ ri, i=1,2, L, n }
All the sensors in railcar due to installation site be not in approximately the same plane, so need by when
Between, mounting distance parameter all synchronizes processing, the available data of all the sensors (its as shown in table 1 to obtained data
Middle mounting coordinate is three-dimensional coordinate).
Each sensing data of table 1
Sensor type | Time parameter | Installation site | Data |
Laser scanner 1 | 50Hz | (0,0,0) | 190 ° of totally 760 dot profile data |
Laser scanner 2 | 50Hz | (0,0,0) | 190 ° of totally 760 dot profile data |
Laser range finder 1 | 50Hz | (a1,b1,c1) | Single-point distance value |
Laser range finder 2 | 50Hz | (a2,b2,c2) | Single-point distance value |
Laser range finder 3 | 50Hz | (a3,b3,c3) | Single-point distance value |
Inclinator 1 | 50Hz | With the level of laser range finder 1 | The setting angle of laser range finder 1 |
Inclinator 2 | 50Hz | With the level of laser range finder 3 | The setting angle of laser range finder 3 |
Inclinator 3 | 50Hz | With railcar surface level | Orbit angle is poor |
Measuring wheel | 50kHz | (a4,b4,c4) | Railcar operating range |
Positioning shooting machine | 4Hz | (a4,b5,c4) | Section circumferential weld number, circumferential weld position |
As can be seen from the table, laser range finder, measuring wheel, positioning shooting machine installation site it is not same in laser scanning
One point, so the data that they are collected in the same time can not be used for the profile value for calculating current section, it is necessary to pass through
Time and mounting distance obtain the value that all the sensors of same position synchronization match after being calculated, with these values
Profiled outline value can accurately be calculated.
Assuming that railcar is at the uniform velocity being advanced, laser range finder, measuring wheel and laser scanner can pass through following public affairs
Formula carries out data syn-chronization:
Speed is obtained by measuring wheel:V=L/T
L is operating range, and T is running time, and a speed was calculated at interval of one second.
The time interval of three laser range finders and measuring wheel and laser scanner data is:T=cx/V
Go to find the data at all the sensors same time point by this interval time.
Claims (10)
1. a kind of detection method of deformation of cross section of subway tunnel, it is characterised in that methods described comprises the following steps:
1) using tunnel track center line as origin, the rectangular coordinate system of subway tunnel section is set up;
2) DATA REASONING is carried out to tunnel cross-section according to goniometer, laser range finder and laser scanner, obtains coordinate offset angle α
With all measurement points of subway tunnel section in step 1) coordinate under the rectangular coordinate system set up, wherein ith measurement point
Coordinate is (xi,yi);
3) by ellipse fitting algorithm to step 2) obtained basic measurement data carries out curve fitting, and obtains coincidence loss requirement
Tunnel cross-section form parameter;
4) judgment step 3) whether the deviation of the obtained form parameter of tunnel cross-section and history form parameter exceed deviation threshold,
If then showing that obvious deformation occurs for subway tunnel section, if otherwise showing, obvious deformation does not occur for subway tunnel section.
2. the detection method of deformation of cross section of subway tunnel according to claim 1, it is characterised in that the step 3) it is specific
For:
31) according to all measurement points in step 1) coordinate under the rectangular coordinate system set up, carry out basic fitting and obtain tunnel breaking
The initial elliptic equation in face;
32) all measurement points are asked in step 31) orthogonal reference point on obtained initial elliptic equation;
33) according to step 32) obtained orthogonal reference point, the elliptic equation after being optimized by iterative;
34) fitting parameter of the elliptic equation after optimization is asked for, judges whether coincidence loss is required fitting parameter, if then will be excellent
Elliptic equation after change as tunnel cross-section form parameter, if otherwise return to step 33).
3. the detection method of deformation of cross section of subway tunnel according to claim 2, it is characterised in that the step 32) tool
Body is:The point using measurement point as primary iteration, to orthogonal correlated condition solve obtaining orthogonal correlation by Newton iteration method
The coordinate of point, the orthogonal correlated condition is specially:
f2(x', y')=b2x'(yi-y')-a2y'(xi-x')2=0
Wherein, a and b are the basic parameter of initial elliptic equation, and (x', y') is the coordinate of orthogonal reference point, (xi,yi) it is to survey
Measure the coordinate of point.
4. the detection method of deformation of cross section of subway tunnel according to claim 2, it is characterised in that the step 33) tool
Body is:
331) coordinate X' of the orthogonal reference point under natural system of coordinates is obtained by coordinate transform;
332) using step 331) orthogonal reference point under obtained natural system of coordinates, according to minimizing geometric distance condition, pass through
Iteration asks for the point coordinates X on the elliptic equation after optimization;
333) according to step 332) point coordinates on elliptic equation after obtained optimization, carry out after ellipse fitting optimized
Elliptic equation.
5. the detection method of deformation of cross section of subway tunnel according to claim 4, it is characterised in that the coordinate transform tool
Body is:
X'=R-1x'+Xc
Wherein, x' is coordinate of the orthogonal reference point under rectangular coordinate system, XcFor rectangular coordinate system origin under natural system of coordinates
Coordinate.
6. the detection method of deformation of cross section of subway tunnel according to claim 4, it is characterised in that the minimum geometry away from
It is specially from condition:
MinE=min { (X-X')T(X-X')}。
7. the detection method of deformation of cross section of subway tunnel according to claim 2, it is characterised in that the fitting parameter bag
Include fitting precision σ and error of fitting Δ ε.
8. the detection method of deformation of cross section of subway tunnel according to claim 1, it is characterised in that the step 3) also wrap
The abnormity point rejected in curve fitting process is included, is specially:
351) according to the curvilinear equation obtained after curve matching, residual error of each measurement point relative to curvilinear equation is calculated;
352) by step 351) in the corresponding residual error of obtained all measurement points constitute residual vector, examined by sample median
Method finds the exceptional value in residual vector;
353) by step 352) in the corresponding measurement point of obtained exceptional value reject, re-start curve matching, obtain meeting mistake
The form parameter for the tunnel cross-section that difference is required.
9. the detection method of deformation of cross section of subway tunnel according to claim 1, it is characterised in that the subway tunnel breaks
All measurement points in face include laser range finder and scan the low precision survey that obtained high-acruracy survey point and laser scanner are obtained
Amount point.
10. the detection method of deformation of cross section of subway tunnel according to claim 9, it is characterised in that the high accuracy is surveyed
It is specially with time interval t, the time interval t between amount point and low precision measure point:
Wherein, c is the light velocity, and x is the distance between laser range finder and laser scanner, and V is the inspection of deformation of cross section of subway tunnel
Survey the pace of device.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710286331.XA CN106969749B (en) | 2017-04-27 | 2017-04-27 | A kind of detection method of deformation of cross section of subway tunnel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710286331.XA CN106969749B (en) | 2017-04-27 | 2017-04-27 | A kind of detection method of deformation of cross section of subway tunnel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106969749A true CN106969749A (en) | 2017-07-21 |
CN106969749B CN106969749B (en) | 2019-12-03 |
Family
ID=59333427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710286331.XA Active CN106969749B (en) | 2017-04-27 | 2017-04-27 | A kind of detection method of deformation of cross section of subway tunnel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106969749B (en) |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107957241A (en) * | 2017-12-29 | 2018-04-24 | 石家庄铁道大学 | Subway tunnel section center of circle determining device and method |
CN108036765A (en) * | 2017-11-22 | 2018-05-15 | 山东大学 | A kind of tunnel clearance convergence monitoring measurement system and method based on three base-line methods |
CN108896014A (en) * | 2018-05-11 | 2018-11-27 | 同济大学 | A kind of tunnel cross-section measurement device and its measurement method |
CN108917638A (en) * | 2018-09-25 | 2018-11-30 | 浙江科技学院 | The mounted measuring apparatus for vehicle of subway tunnel 3 d deformation monitoring based on datum tool |
CN109029350A (en) * | 2018-08-02 | 2018-12-18 | 南京航空航天大学 | A kind of tunnel axis extracts and section Convergence analysis method and device |
CN109186480A (en) * | 2018-09-19 | 2019-01-11 | 成都理工大学 | Tunnel surrounding scanning and observation system based on double-shielded TBM technique |
CN110044289A (en) * | 2019-04-17 | 2019-07-23 | 上海同岩土木工程科技股份有限公司 | The multibreak face convergent deformation video monitoring device in tunnel and method based on autozoom |
CN110160463A (en) * | 2019-05-17 | 2019-08-23 | 中国电建集团西北勘测设计研究院有限公司 | A kind of subway tunnel out-of-roundness detection method based on static laser scanning |
CN110309525A (en) * | 2019-03-22 | 2019-10-08 | 北京北科安地科技发展有限公司 | A kind of side slope geometric distortion and destroy trend calculation method |
CN110455211A (en) * | 2019-08-19 | 2019-11-15 | 云南航天工程物探检测股份有限公司 | A kind of automatic monitoring measurement method based on laser section ranging |
CN110588707A (en) * | 2019-09-12 | 2019-12-20 | 江苏航天大为科技股份有限公司 | Method for detecting rail transit tunnel clearance |
CN110986815A (en) * | 2020-03-05 | 2020-04-10 | 浙江交工集团股份有限公司 | Tunnel construction monitoring and measuring method based on three-dimensional laser point cloud |
CN111024037A (en) * | 2019-09-25 | 2020-04-17 | 北京工业大学 | Circular shield tunnel section fitting method based on three-dimensional laser scanning technology |
CN111412851A (en) * | 2020-04-13 | 2020-07-14 | 成都大亦科技有限公司 | Method for measuring deformation based on laser |
CN112097669A (en) * | 2020-11-17 | 2020-12-18 | 南京派光智慧感知信息技术有限公司 | Method for monitoring deformation of structure in tunnel based on laser ranging |
CN112862790A (en) * | 2021-02-18 | 2021-05-28 | 中国矿业大学(北京) | Subway tunnel crack positioning device and method based on linear array camera |
CN113074694A (en) * | 2021-03-31 | 2021-07-06 | 中煤科工集团重庆研究院有限公司 | Automatic monitoring device for tunnel section deformation |
CN113932762A (en) * | 2021-10-19 | 2022-01-14 | 广东电网有限责任公司 | Cable deformation measuring method and device and computer storage medium |
CN114485539A (en) * | 2022-01-04 | 2022-05-13 | 中铁第四勘察设计院集团有限公司 | Tunnel section deformation measuring method and device and storage medium |
CN117169898A (en) * | 2023-09-07 | 2023-12-05 | 南京安透可智能***有限公司 | Oval fitting-based underground pipe network branch pipe and size detection method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246007A (en) * | 2007-02-12 | 2008-08-20 | 上海地铁运营有限公司 | Tunnel convergence monitoring method |
CN101685012A (en) * | 2008-09-27 | 2010-03-31 | 上海市机械施工有限公司 | Measuring and calculating method of subway tunnel section |
CN103900528A (en) * | 2012-12-28 | 2014-07-02 | 朱志洁 | Detection method for three-dimensional space section roundness of large component |
CN104792274A (en) * | 2015-04-02 | 2015-07-22 | 同济大学 | Method for measuring circular tunnel convergence deformation |
-
2017
- 2017-04-27 CN CN201710286331.XA patent/CN106969749B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101246007A (en) * | 2007-02-12 | 2008-08-20 | 上海地铁运营有限公司 | Tunnel convergence monitoring method |
CN101685012A (en) * | 2008-09-27 | 2010-03-31 | 上海市机械施工有限公司 | Measuring and calculating method of subway tunnel section |
CN103900528A (en) * | 2012-12-28 | 2014-07-02 | 朱志洁 | Detection method for three-dimensional space section roundness of large component |
CN104792274A (en) * | 2015-04-02 | 2015-07-22 | 同济大学 | Method for measuring circular tunnel convergence deformation |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108036765A (en) * | 2017-11-22 | 2018-05-15 | 山东大学 | A kind of tunnel clearance convergence monitoring measurement system and method based on three base-line methods |
CN107957241B (en) * | 2017-12-29 | 2019-10-25 | 石家庄铁道大学 | Subway tunnel section center of circle determining device and method |
CN107957241A (en) * | 2017-12-29 | 2018-04-24 | 石家庄铁道大学 | Subway tunnel section center of circle determining device and method |
CN108896014A (en) * | 2018-05-11 | 2018-11-27 | 同济大学 | A kind of tunnel cross-section measurement device and its measurement method |
CN109029350A (en) * | 2018-08-02 | 2018-12-18 | 南京航空航天大学 | A kind of tunnel axis extracts and section Convergence analysis method and device |
CN109029350B (en) * | 2018-08-02 | 2023-05-23 | 南京航空航天大学 | Tunnel axis extraction and section convergence analysis method and device |
CN109186480A (en) * | 2018-09-19 | 2019-01-11 | 成都理工大学 | Tunnel surrounding scanning and observation system based on double-shielded TBM technique |
CN112525093A (en) * | 2018-09-19 | 2021-03-19 | 成都理工大学 | System for building tunnel three-dimensional model |
CN112525093B (en) * | 2018-09-19 | 2022-07-22 | 成都理工大学 | System for establishing tunnel three-dimensional model based on double-shield TBM (Tunnel boring machine) process |
CN108917638A (en) * | 2018-09-25 | 2018-11-30 | 浙江科技学院 | The mounted measuring apparatus for vehicle of subway tunnel 3 d deformation monitoring based on datum tool |
CN110309525A (en) * | 2019-03-22 | 2019-10-08 | 北京北科安地科技发展有限公司 | A kind of side slope geometric distortion and destroy trend calculation method |
CN110044289A (en) * | 2019-04-17 | 2019-07-23 | 上海同岩土木工程科技股份有限公司 | The multibreak face convergent deformation video monitoring device in tunnel and method based on autozoom |
CN110160463A (en) * | 2019-05-17 | 2019-08-23 | 中国电建集团西北勘测设计研究院有限公司 | A kind of subway tunnel out-of-roundness detection method based on static laser scanning |
CN110455211A (en) * | 2019-08-19 | 2019-11-15 | 云南航天工程物探检测股份有限公司 | A kind of automatic monitoring measurement method based on laser section ranging |
CN110588707A (en) * | 2019-09-12 | 2019-12-20 | 江苏航天大为科技股份有限公司 | Method for detecting rail transit tunnel clearance |
CN111024037A (en) * | 2019-09-25 | 2020-04-17 | 北京工业大学 | Circular shield tunnel section fitting method based on three-dimensional laser scanning technology |
CN110986815A (en) * | 2020-03-05 | 2020-04-10 | 浙江交工集团股份有限公司 | Tunnel construction monitoring and measuring method based on three-dimensional laser point cloud |
CN111412851A (en) * | 2020-04-13 | 2020-07-14 | 成都大亦科技有限公司 | Method for measuring deformation based on laser |
CN112097669A (en) * | 2020-11-17 | 2020-12-18 | 南京派光智慧感知信息技术有限公司 | Method for monitoring deformation of structure in tunnel based on laser ranging |
CN112862790A (en) * | 2021-02-18 | 2021-05-28 | 中国矿业大学(北京) | Subway tunnel crack positioning device and method based on linear array camera |
CN112862790B (en) * | 2021-02-18 | 2023-08-22 | 中国矿业大学(北京) | Subway tunnel crack positioning device and method based on linear array camera |
CN113074694A (en) * | 2021-03-31 | 2021-07-06 | 中煤科工集团重庆研究院有限公司 | Automatic monitoring device for tunnel section deformation |
CN113932762A (en) * | 2021-10-19 | 2022-01-14 | 广东电网有限责任公司 | Cable deformation measuring method and device and computer storage medium |
CN114485539A (en) * | 2022-01-04 | 2022-05-13 | 中铁第四勘察设计院集团有限公司 | Tunnel section deformation measuring method and device and storage medium |
CN114485539B (en) * | 2022-01-04 | 2024-05-03 | 中铁第四勘察设计院集团有限公司 | Tunnel section deformation measuring method, device and storage medium |
CN117169898A (en) * | 2023-09-07 | 2023-12-05 | 南京安透可智能***有限公司 | Oval fitting-based underground pipe network branch pipe and size detection method |
Also Published As
Publication number | Publication date |
---|---|
CN106969749B (en) | 2019-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106969749A (en) | A kind of detection method of deformation of cross section of subway tunnel | |
CN106401643B (en) | Tunnel Overbreak & Underbreak detection method based on laser point cloud | |
CN111307031B (en) | Buried pipeline safety state monitoring and early warning method | |
CN108489403B (en) | Rapid and fine evaluation method for joint attitude of surface mine slope rock mass based on three-dimensional laser scanning | |
CN109708615A (en) | A kind of subway tunnel limit dynamic testing method based on laser scanning | |
WO2018090348A1 (en) | Automatic collision avoidance method for adjacent wells | |
CN108731640A (en) | Metro Clearance detection method based on point cloud data and detecting system | |
CN109902425A (en) | The tunnel cross-section extracting method of ground formula point cloud | |
CN106091951A (en) | A kind of municipal rail train wheel rim parameter on-line detecting system and method | |
CN104808248A (en) | Attitude of stratum detecting method based on electric imaging logging full hole image | |
CN109409719B (en) | Method for evaluating vulnerability of long-distance oil and gas pipeline under landslide action | |
CN103675920A (en) | Nondestructive test method for depth and horizontal position of hidden crack of road base | |
CN106080661A (en) | Train guide rail profile undulatory wear measuring method | |
CN105835902A (en) | Method for detecting diameter of wheel based on laser displacement sensors | |
CN111811466B (en) | Method for measuring parameters and maintenance of existing railway ballast railway line | |
CN111895962A (en) | Method for monitoring deformation of operation subway station | |
CN115112044A (en) | Wheel set size measurement method based on light spot cloud data of multi-line structure | |
CN107552994B (en) | A kind of weld penetration, molten wide and the high detection method of weldering | |
CN116128834A (en) | Rail transit platform limit rapid detection method based on three-dimensional laser scanning | |
Du et al. | Dislocation detection of shield tunnel based on dense cross-sectional point clouds | |
CN106767421A (en) | Motor-car vehicle body critical size detecting system solution based on multi-vision visual | |
CN106918812A (en) | Underwater foundation detection means and method based on supersonic sounding and Three Dimensional Reconfiguration | |
CN107139968A (en) | Wheel out of round degree detection means and method based on laser displacement sensor | |
CN105318860A (en) | Building sedimentation observation method and structure | |
CN111455787B (en) | Pavement detection system based on pavement three-dimensional digitization |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A detection method of section deformation of subway tunnel Effective date of registration: 20220310 Granted publication date: 20191203 Pledgee: Agricultural Bank of China Limited Shanghai Wujiaochang sub branch Pledgor: SHANGHAI TONGYAN CIVIL ENGINEERING TECHNOLOGY CO.,LTD. Registration number: Y2022310000048 |