CN111365076A - Shield tunnel linear early warning method - Google Patents

Shield tunnel linear early warning method Download PDF

Info

Publication number
CN111365076A
CN111365076A CN202010189479.3A CN202010189479A CN111365076A CN 111365076 A CN111365076 A CN 111365076A CN 202010189479 A CN202010189479 A CN 202010189479A CN 111365076 A CN111365076 A CN 111365076A
Authority
CN
China
Prior art keywords
ring
segment
current
actual
offset
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
Application number
CN202010189479.3A
Other languages
Chinese (zh)
Other versions
CN111365076B (en
Inventor
安国勇
仇峰涛
徐宏
杨永强
刘丹
安欢
王勇
韩兴博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Railway First Engineering Group Co Ltd
Intelligent Technology Branch of China Railway First Engineering Group Co Ltd
Urban Rail Transit Engineering Co Ltd of China Railway First Engineering Group Co Ltd
Original Assignee
China Railway First Engineering Group Co Ltd
Intelligent Technology Branch of China Railway First Engineering Group Co Ltd
Urban Rail Transit Engineering Co Ltd of China Railway First Engineering Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Railway First Engineering Group Co Ltd, Intelligent Technology Branch of China Railway First Engineering Group Co Ltd, Urban Rail Transit Engineering Co Ltd of China Railway First Engineering Group Co Ltd filed Critical China Railway First Engineering Group Co Ltd
Priority to CN202010189479.3A priority Critical patent/CN111365076B/en
Publication of CN111365076A publication Critical patent/CN111365076A/en
Application granted granted Critical
Publication of CN111365076B publication Critical patent/CN111365076B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18—Special adaptations of signalling or alarm devices
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D11/00—Lining tunnels, galleries or other underground cavities, e.g. large underground chambers; Linings therefor; Making such linings in situ, e.g. by assembling
    • E21D11/04—Lining with building materials
    • E21D11/08—Lining with building materials with preformed concrete slabs
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Abstract

The invention discloses a shield tunnel line shape early warning method, which comprises the following steps: firstly, calculating the actual accumulated rotation angle of the current pipe sheet ring I; secondly, calculating the actual rambling offset of the current tube sheet ring I; thirdly, calculating the residual error of the base counting offset deviation; fourthly, calculating the standard deviation of the residual error of the offset of each ring of the pipe sheet ring; fifthly, calculating a residual coefficient of each ring of the tube sheet ring; sixthly, calculating a residual coefficient oscillation value; and seventhly, early warning of the shield tunnel line shape. The method comprises the steps of calculating the included angle between the axes of two adjacent ring pipe pieces, accumulating the included angle to obtain the actual accumulated rotation angle of the current pipe piece ring, calculating the actual rampart offset of each ring pipe piece, solving the residual standard deviation and the residual coefficient of the rampart offset deviation, obtaining the residual coefficient oscillation value and giving out the linear early warning signal of the shield tunnel, wherein the data acquisition does not depend on a total station, the measurement accumulated error is not generated, the early warning is timely and efficient, the precision is high, and the accident of serious deviation of the splicing route of the shield pipe pieces is reduced.

Description

Shield tunnel linear early warning method
Technical Field
The invention belongs to the technical field of shield tunnel line shape early warning, and particularly relates to a shield tunnel line shape early warning method.
Background
In the shield construction propelling process, the actual assembling line shape and the original design line shape of the shield segment can deviate under the influences of excavation face soil pressure, surrounding rock pressure, stratum, shield machine operator level and the like. This problem is also known as the snake-like problem of shield tunneling. When the deviation is small, a certain deviation rectifying measure is timely adopted to enable the shield tunneling to meet the original design line shape requirement. However, when the deviation is large, the shield deviation rectification has no success possibility. Therefore, the shield machine cannot reach the receiving well according to the route, the workload is increased, the project is delayed, and the construction risk is increased. Therefore, the deviation between the actual splicing position of the shield segment and the original design position is warned in time, and shield construction personnel are reminded to find out the construction deviation in time, so that the deviation correction is implemented as early as possible, and the method has important significance for ensuring the smooth construction of the tunnel.
Most of the existing shield construction propulsion is guided and early-warned by a shield machine laser guide system, the working principle of the shield machine laser guide system is that the direction of the tunneling machine is calculated according to the angle relation between guide laser emitted by a total station and a target unit, and the deviation of the shield machine and the design axis is converted by combining the actual application coordinate. Before starting, the design line calculation data is imported into a guide system, and the coordinate precision of a base point and the like are the core and the key of tunnel construction control. Aiming at the problem, the conventional solution at present is to manually measure and record the actual axis coordinate data and compare the actual data with the original design coordinates. When the deviation value exceeds a certain range, the method reports to a superior administrative department for processing, the method has subjectivity and timely ductility, the axis deviation state cannot be timely and efficiently processed, meanwhile, the data measurement depends on the precision of the total station, once the precision of a control point of the total station is not enough or an error occurs in one measurement, errors are introduced to all subsequent data, the input data of a laser guide system of the shield machine is wrong, and the guide of the laser guide system of the shield machine is wrong and cannot be known by self.
Disclosure of Invention
The invention aims to solve the technical problem that the prior art is not enough, and provides a shield tunnel linear early warning method, which is characterized in that the included angle between the axes of two adjacent ring pipe segments is calculated and accumulated to obtain the actual accumulated rotation angle of the current pipe segment ring, then the actual rampart offset of each ring pipe segment is calculated, the residual standard difference and the residual coefficient of the rampart offset deviation are obtained, the residual coefficient oscillation value is obtained, and a shield tunnel linear early warning signal is given, the data acquisition does not depend on a total station, the measurement accumulated error is not generated, the early warning is timely and efficient, the precision is high, the serious shield segment splicing route deviation accident is reduced, and the method is convenient to popularize and use.
In order to solve the technical problems, the invention adopts the technical scheme that: a shield tunnel linear early warning method is characterized by comprising the following steps:
step one, according to a formula
Figure BDA0002415344810000021
Calculating the actual cumulative rotational angle β of the current segment ring IIWherein, βI' is the axis of the current segment ring IActual angle of rotation between line and ring axis of the segment immediately preceding it, βi' is the actual rotation angle between the axis of the first I ring segment of the current segment ring I and the axis of the first I +1 ring segment of the current segment ring I, I and I are integers,
Figure BDA0002415344810000022
for the actual half wedge amount of the current segment ring I,
Figure BDA0002415344810000023
the actual half wedge-shaped amount of the ring segment I in front of the current segment ring I is shown, and R is the segment radius;
step two, according to a formula HzI＝∑LIβICalculating the actual rampart offset H of the current tube sheet ring IzIWherein L isIThe central ring width of the current tube sheet ring I;
step three, calculating the residual error of the base gauge offset deviation, wherein the process is as follows:
step 301, according to formula yz＝HzI-HszICalculating the deviation y between the actual rampart offset of the current tube sheet ring I and the design rampart offset of the current tube sheet ring IzWherein H isszICalculating the offset for the design of the current tube sheet ring I;
step 302, performing linear regression on the segment ring number and the segment base gauge offset deviation of each ring to obtain a linear regression function Y between the segment ring number and the segment base gauge offset deviationz＝δ0+δ1X + epsilon, wherein X is the independent variable of the number of the segment ring, delta0Is a constant number, δ1Is a regression coefficient, ε is the residual error of the rampart offset, YzCalculating a standoff distance deviation dependent variable for the segment base;
according to a linear regression function Y between a segment ring number and a segment rampart offset deviationz＝δ0+δ1X + epsilon obtains the residual error of the rampart offset deviation of each ring of the tube sheet ring;
step four, according to the formula
Figure BDA0002415344810000031
Each ring of the computing pipe sheet ringStandard deviation sigma of residual error of rampart offset deviationXWherein, epsilonXIs the residual error of the rampart offset corresponding to the segment ring number, x is the accumulated ring number,
Figure BDA0002415344810000032
the residual mean value of the base gauge offset deviation under the accumulated ring number is obtained;
step five, according to the formula
Figure BDA0002415344810000033
Calculating residual error coefficient s of each ring of tube sheet ringX；
Step six, according to the formula
Figure BDA0002415344810000034
Calculating residual coefficient oscillation value ZD, wherein XqIs the initial segment ring number, sqIs the residual coefficient, X, corresponding to the ring number of the starting segmentJTo stop the ring number of the segment, sJIs the residual coefficient, s, corresponding to the ring number of the cut-off segmentmaxFor the maximum residual coefficient under the accumulated ring number, XmaxThe segment ring number corresponding to the maximum residual error coefficient under the accumulated ring number;
seventhly, shield tunnel line shape early warning: according to the formula
Figure BDA0002415344810000035
Calculating a shield segment splicing linear fitting error characteristic value P;
when the characteristic value P of the splicing linear fitting error of the shield segments is less than 50, the shield tunnel line shape is normal, and no early warning signal is sent out;
when the characteristic value of the shield segment splicing linear fitting error is more than or equal to 50 and less than 67, the splicing quality of the shield segments is controlled, and the shield tunnel linear early warning is yellow early warning;
when the characteristic value 67 of the splicing linear fitting error of the shield segments is not more than P and less than 100, indicating that the splicing quality of the shield segments deviates from a design route, and at the moment, the linear early warning of the shield tunnel is an orange early warning;
when the characteristic value P of the shield segment splicing linear fitting error is more than or equal to 100, the shield segment splicing quality is unqualified, and at the moment, the shield tunnel linear early warning is red early warning and needs to be shut down for on-site investigation.
The linear early warning method for the shield tunnel is characterized by comprising the following steps: the actual half-wedge amount in the first step comprises an actual horizontal half-wedge amount and an actual vertical half-wedge amount, the actual horizontal half-wedge amount comprises an actual horizontal left half-wedge amount and an actual horizontal right half-wedge amount, and the actual vertical half-wedge amount comprises an actual vertical top half-wedge amount and an actual vertical bottom half-wedge amount; the actual accumulated rotation angle of the current tube sheet ring I comprises an actual horizontal accumulated rotation angle of the current tube sheet ring I and an actual vertical accumulated rotation angle of the current tube sheet ring I;
the actual rampart offset of the current tube sheet ring I in the step two comprises an actual horizontal rampart offset of the current tube sheet ring I and an actual vertical rampart offset of the current tube sheet ring I, and the design rampart offset of the current tube sheet ring I comprises a design horizontal rampart offset of the current tube sheet ring I and a design vertical rampart offset of the current tube sheet ring I;
in the third step, the residual error of the base gauge support distance deviation of each ring of the tube sheet ring comprises the residual error of the horizontal base gauge support distance deviation of each ring of the tube sheet ring and the residual error of the vertical base gauge support distance deviation of each ring of the tube sheet ring;
the standard deviation of the residual error of the offset of each ring of the tube sheet ring in the fourth step comprises the standard deviation of the residual error of the offset of each ring of the horizontal base of the tube sheet ring and the standard deviation of the residual error of the offset of each ring of the vertical base of the tube sheet ring;
the residual error coefficients of each ring of the tube sheet ring in the step five comprise a horizontal residual error coefficient of each ring of the tube sheet ring and a vertical residual error coefficient of each ring of the tube sheet ring;
the residual error coefficient oscillation value in the sixth step comprises a horizontal residual error coefficient oscillation value and a vertical residual error coefficient oscillation value;
and seventhly, performing shield tunnel line shape early warning on shield segment assembly horizontal error early warning and shield segment assembly vertical error early warning.
The linear early warning method for the shield tunnel is characterized by comprising the following steps: in the step one, the current pipe sheet ring is a 16-point pipe sheet ring;
when the current pipe sheet ring is a standard ring, according to a formula
Figure BDA0002415344810000041
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000042
And actual vertical bottom half wedge amount
Figure BDA0002415344810000043
When the current pipe sheet ring is a left-turning ring, according to a formula
Figure BDA0002415344810000044
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000045
And actual vertical bottom half wedge amount
Figure BDA0002415344810000046
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000047
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000051
When the current pipe sheet ring is a right-turning ring, according to a formula
Figure BDA0002415344810000052
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000053
And actual vertical bottom half wedge amount
Figure BDA0002415344810000054
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000055
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000056
The linear early warning method for the shield tunnel is characterized by comprising the following steps: in the step one, the current pipe sheet ring is a 10-point pipe sheet ring;
when the current pipe sheet ring is a standard ring, according to a formula
Figure BDA0002415344810000057
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000058
And actual vertical bottom half wedge amount
Figure BDA0002415344810000059
When the current pipe sheet ring is a left-turning ring, according to a formula
Figure BDA00024153448100000510
Calculating the half wedge amount on the right side of the actual level
Figure BDA00024153448100000511
And actual vertical bottom half wedge amount
Figure BDA00024153448100000512
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA00024153448100000513
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA00024153448100000514
When the current pipe sheet ring is a right-turning ring, according to a formula
Figure BDA0002415344810000061
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000062
And actual vertical bottom half wedge amount
Figure BDA0002415344810000063
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000064
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000065
The linear early warning method for the shield tunnel is characterized by comprising the following steps: in the step one, the current pipe sheet ring is a universal pipe sheet ring;
according to the formula
Figure BDA0002415344810000066
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000067
And actual vertical bottom half wedge amount
Figure BDA0002415344810000068
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000069
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA00024153448100000610
Compared with the prior art, the invention has the following advantages:
1. the method calculates the included angle between the axes of two adjacent ring duct pieces and accumulates the included angle to obtain the actual accumulated corner of the current duct piece ring, the actual accumulated corner of each ring duct piece is obtained by accumulating the axis corners of the front 100 ring duct pieces of the current ring duct piece, the calculated actual accumulated corner can be compared with the designed accumulated corner, whether the assembling trend of the duct pieces of the tunnel duct is correct or not is reflected, the observation reliability is high, for the duct pieces with less than hundreds of rings, the hundreds of rings are completely supplemented according to the negative ring supplementing principle, the gaps of the supplemented duct pieces are all 0, the effect that the actual tunnel duct pieces are completely matched with the designed tunnel duct pieces is achieved, and the method is convenient to popularize and use.
2. The method calculates the actual rammeter supporting distance of each ring of pipe pieces, the actual rammeter supporting distance of each ring of pipe pieces is obtained by the accumulation of the actual accumulated half wedge-shaped amount of each ring of pipe pieces, the residual error of the actual rammeter supporting distance is obtained through linear regression, the approximation effect is good, the residual error standard deviation of the rammeter supporting distance deviation is obtained, the aggregation degree of data is observed, the maximum residual error coefficient under the accumulated ring number and the pipe piece ring number corresponding to the maximum residual error coefficient under the accumulated ring number are searched by drawing a residual error coefficient distribution diagram, the area of the composition is obtained by a three-point method through the residual error coefficient corresponding to the initial pipe piece ring number, the maximum residual error coefficient under the accumulated ring number and the residual error coefficient corresponding to the cut-off pipe piece ring number, the residual error coefficient oscillation value is obtained, the reliability and stability are realized.
3. The method has simple steps, utilizes the residual error coefficient oscillation value to calculate the linear early warning signal of the shield tunnel, has no dependence on a total station for data acquisition, does not generate measurement accumulated errors, has timely and efficient early warning and high precision, reduces the accident of serious deviation of the splicing route of the shield segments, and is convenient for popularization and use.
In summary, the invention calculates the included angle between the axes of two adjacent ring duct pieces and accumulates the included angle to obtain the actual accumulated rotation angle of the current duct piece ring, then calculates the actual stacking support distance of each ring duct piece, calculates the residual standard deviation and the residual coefficient of the stacking support distance deviation, obtains the residual coefficient oscillation value and gives the shield tunnel linear early warning signal, the data acquisition does not depend on a total station, the measurement accumulated error is not generated, the early warning is timely and efficient, the precision is high, the accident of serious deviation of the shield duct piece assembling route is reduced, and the method is convenient to popularize and use.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
FIG. 1 is a block diagram of the process flow of the present invention.
Detailed Description
As shown in fig. 1, the shield tunnel line shape early warning method of the present invention includes the following steps:
step one, according to a formula
Figure BDA0002415344810000071
Calculating the actual cumulative rotational angle β of the current segment ring IIWherein, βI' is the actual angle of rotation between the axis of the current segment ring I and the axis of the segment ring immediately preceding it, βi' is the actual rotation between the axis of the first I ring segment of the current segment ring I and the axis of the first I +1 ring segment of the current segment ring IThe angles, I and I are integers,
Figure BDA0002415344810000081
for the actual half wedge amount of the current segment ring I,
Figure BDA0002415344810000082
the actual half wedge-shaped amount of the ring segment I in front of the current segment ring I is shown, and R is the segment radius;
it should be noted that, the actual accumulated corner of the current segment ring is obtained by calculating the included angle between the axes of two adjacent two ring segments and accumulating the included angle, the actual accumulated corner of each ring segment is obtained by accumulating the axis corners of the previous 100 ring segments of the current ring segment, the calculated actual accumulated corner can be compared with the design accumulated corner, whether the assembling trend of the segments of the reaction tunnel is correct or not is reflected, the observation reliability is high, for the segments which are not enough for hundreds of rings, the hundreds of rings are supplemented according to the negative ring supplementing principle, the gaps of the segments which are supplemented are all 0, and the effect that the actual tunnel segment is completely matched with the design tunnel segment is achieved.
Step two, according to a formula HzI＝∑LIβICalculating the actual rampart offset H of the current tube sheet ring IzIWherein L isIThe central ring width of the current tube sheet ring I;
step three, calculating the residual error of the base gauge offset deviation, wherein the process is as follows:
step 301, according to formula yz＝HzI-HszICalculating the deviation y between the actual rampart offset of the current tube sheet ring I and the design rampart offset of the current tube sheet ring IzWherein H isszICalculating the offset for the design of the current tube sheet ring I;
step 302, performing linear regression on the segment ring number and the segment base gauge offset deviation of each ring to obtain a linear regression function Y between the segment ring number and the segment base gauge offset deviationz＝δ0+δ1X + epsilon, wherein X is the independent variable of the number of the segment ring, delta0Is a constant number, δ1Is a regression coefficient, ε is the residual error of the rampart offset, YzGauge stand offset variation for duct piece stackingAn amount;
according to a linear regression function Y between a segment ring number and a segment rampart offset deviationz＝δ0+δ1X + epsilon obtains the residual error of the rampart offset deviation of each ring of the tube sheet ring;
step four, according to the formula
Figure BDA0002415344810000083
Calculating standard deviation sigma of residual error of each ring of pipe sheet ring rampart offsetXWherein, epsilonXIs the residual error of the rampart offset corresponding to the segment ring number, x is the accumulated ring number,
Figure BDA0002415344810000084
the residual mean value of the base gauge offset deviation under the accumulated ring number is obtained;
step five, according to the formula
Figure BDA0002415344810000085
Calculating residual error coefficient s of each ring of tube sheet ringX；
The method includes the steps of calculating the actual rampart supporting distance of each ring of pipe pieces, obtaining the residual error of the actual rampart supporting distance through linear regression, obtaining a good approximation effect, obtaining the residual error standard deviation of the rampart supporting distance deviation, observing the aggregation degree of data, finding the pipe piece ring number corresponding to the maximum residual error coefficient under the accumulated ring number and the maximum residual error coefficient under the accumulated ring number by drawing a residual error coefficient distribution diagram, obtaining the area of the pipe piece by a three-point method through the residual error coefficient corresponding to the initial pipe piece ring number, the maximum residual error coefficient under the accumulated ring number and the residual error coefficient corresponding to the cut-off pipe piece ring number, obtaining the residual error coefficient oscillation value, and obtaining the reliable and stable use effect.
Step six, according to the formula
Figure BDA0002415344810000091
Calculating residual coefficient oscillation value ZD, wherein XqIs the initial segment ring number, sqCorresponding to the ring number of the starting segmentResidual coefficient, XJTo stop the ring number of the segment, sJIs the residual coefficient, s, corresponding to the ring number of the cut-off segmentmaxFor the maximum residual coefficient under the accumulated ring number, XmaxThe segment ring number corresponding to the maximum residual error coefficient under the accumulated ring number;
seventhly, shield tunnel line shape early warning: according to the formula
Figure BDA0002415344810000092
Calculating a shield segment splicing linear fitting error characteristic value P;
when the characteristic value P of the splicing linear fitting error of the shield segments is less than 50, the shield tunnel line shape is normal, and no early warning signal is sent out;
when the characteristic value of the shield segment splicing linear fitting error is more than or equal to 50 and less than 67, the splicing quality of the shield segments is controlled, and the shield tunnel linear early warning is yellow early warning;
when the characteristic value 67 of the splicing linear fitting error of the shield segments is not more than P and less than 100, indicating that the splicing quality of the shield segments deviates from a design route, and at the moment, the linear early warning of the shield tunnel is an orange early warning;
when the characteristic value P of the shield segment splicing linear fitting error is more than or equal to 100, the shield segment splicing quality is unqualified, and at the moment, the shield tunnel linear early warning is red early warning and needs to be shut down for on-site investigation.
The method has simple steps, calculates the linear early warning signal of the shield tunnel by utilizing the residual coefficient oscillation value, obtains data without depending on a total station, does not generate measurement accumulated error, has timely and efficient early warning and high precision, and reduces the accident of serious deviation of the splicing route of the shield segments.
In this embodiment, the actual half-wedge amount in the first step includes an actual horizontal half-wedge amount and an actual vertical half-wedge amount, the actual horizontal half-wedge amount includes an actual horizontal left half-wedge amount and an actual horizontal right half-wedge amount, and the actual vertical half-wedge amount includes an actual vertical top half-wedge amount and an actual vertical bottom half-wedge amount; the actual accumulated rotation angle of the current tube sheet ring I comprises an actual horizontal accumulated rotation angle of the current tube sheet ring I and an actual vertical accumulated rotation angle of the current tube sheet ring I, and the calculation processes of the actual horizontal accumulated rotation angle of the current tube sheet ring I and the actual vertical accumulated rotation angle of the current tube sheet ring I are the same;
in the second step, the actual rampart offset of the current tube sheet ring I comprises an actual horizontal rampart offset of the current tube sheet ring I and an actual vertical rampart offset of the current tube sheet ring I, the calculation processes of the actual horizontal rampart offset of the current tube sheet ring I and the actual vertical rampart offset of the current tube sheet ring I are the same, and the design rampart offset of the current tube sheet ring I comprises a design horizontal rampart offset of the current tube sheet ring I and a design vertical rampart offset of the current tube sheet ring I;
in the third step, the residual error of the rampart gauge support distance deviation of each ring of the tube sheet ring comprises the residual error of the horizontal rampart gauge support distance deviation of each ring of the tube sheet ring and the residual error of the vertical rampart gauge support distance deviation of each ring of the tube sheet ring, and the calculation processes of the residual error of the horizontal rampart gauge support distance deviation of each ring of the tube sheet ring and the residual error of the vertical rampart gauge support distance deviation of each ring of the tube sheet ring are the same;
the standard deviation of the residual error of the offset deviation of each ring of the tube sheet ring in the fourth step comprises the standard deviation of the residual error of the offset deviation of each ring of the horizontal base of the tube sheet ring and the standard deviation of the residual error of the offset deviation of each ring of the vertical base of the tube sheet ring, and the calculation processes of the standard deviation of the residual error of the offset deviation of each ring of the horizontal base of the tube sheet ring and the standard deviation of the residual error of the offset deviation of each ring of the vertical base of the tube sheet ring are the same;
the residual error coefficients of each ring of the tube sheet ring in the step five comprise a horizontal residual error coefficient of each ring of the tube sheet ring and a vertical residual error coefficient of each ring of the tube sheet ring, and the calculation processes of the horizontal residual error coefficient of each ring of the tube sheet ring and the vertical residual error coefficient of each ring of the tube sheet ring are the same;
the residual error coefficient oscillation value in the sixth step comprises a horizontal residual error coefficient oscillation value and a vertical residual error coefficient oscillation value, and the calculation processes of the horizontal residual error coefficient oscillation value and the vertical residual error coefficient oscillation value are the same;
and seventhly, performing shield tunnel linear early warning by using the shield segment splicing horizontal error early warning and the shield segment splicing vertical error early warning, wherein the shield segment splicing horizontal error early warning and the shield segment splicing vertical error early warning have the same process.
In this embodiment, in the first step, the current tube sheet ring is a 16-point tube sheet ring;
when the current pipe sheet ring is a standard ring, according to a formula
Figure BDA0002415344810000111
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000112
And actual vertical bottom half wedge amount
Figure BDA0002415344810000113
When the current pipe sheet ring is a left-turning ring, according to a formula
Figure BDA0002415344810000114
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000115
And actual vertical bottom half wedge amount
Figure BDA0002415344810000116
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000117
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000118
When the current pipe sheet ring is a right-turning ring, according to a formula
Figure BDA0002415344810000119
Calculating realityHorizontal right half wedge amount
Figure BDA00024153448100001110
And actual vertical bottom half wedge amount
Figure BDA00024153448100001111
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA00024153448100001112
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA00024153448100001113
In this embodiment, in the first step, the current tube sheet ring is a 10-point tube sheet ring;
when the current pipe sheet ring is a standard ring, according to a formula
Figure BDA00024153448100001114
Calculating the half wedge amount on the right side of the actual level
Figure BDA00024153448100001115
And actual vertical bottom half wedge amount
Figure BDA00024153448100001116
When the current pipe sheet ring is a left-turning ring, according to a formula
Figure BDA00024153448100001117
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000121
And actual vertical bottom half wedge amount
Figure BDA0002415344810000122
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000123
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000124
When the current pipe sheet ring is a right-turning ring, according to a formula
Figure BDA0002415344810000125
Calculating the half wedge amount on the right side of the actual level
Figure BDA0002415344810000126
And actual vertical bottom half wedge amount
Figure BDA0002415344810000127
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA0002415344810000128
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000129
In this embodiment, in the first step, the current pipe piece ring is a general pipe piece ring;
according to the formula
Figure BDA00024153448100001210
Calculating the half wedge amount on the right side of the actual level
Figure BDA00024153448100001211
And actual vertical bottom half wedge amount
Figure BDA00024153448100001212
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure BDA00024153448100001213
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure BDA0002415344810000131
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and all simple modifications, changes and equivalent structural changes made to the above embodiment according to the technical spirit of the present invention still fall within the protection scope of the technical solution of the present invention.

Claims (5)

1. A shield tunnel linear early warning method is characterized by comprising the following steps:
step one, according to a formula
Figure FDA0002415344800000011
Calculating the actual cumulative rotational angle β of the current segment ring IIWherein, β'IIs the actual angle of rotation between the axis of the current segment ring I and the axis of the segment ring immediately preceding it, β'iIs the actual corner between the axis of the front I ring segment of the current segment ring I and the axis of the front I +1 ring segment of the current segment ring I, wherein I and I are integers,
Figure FDA0002415344800000012
for the actual half wedge amount of the current segment ring I,
Figure FDA0002415344800000013
the actual half wedge-shaped amount of the ring segment I in front of the current segment ring I is shown, and R is the segment radius;
step two, according to a formula HzI＝∑LIβICalculating the actual rampart offset H of the current tube sheet ring IzIWherein L isIThe central ring width of the current tube sheet ring I;
step three, calculating the residual error of the base gauge offset deviation, wherein the process is as follows:
step 301, according to formula yz＝HzI-HszICalculating the deviation y between the actual rampart offset of the current tube sheet ring I and the design rampart offset of the current tube sheet ring IzWherein H isszICalculating the offset for the design of the current tube sheet ring I;
step 302, performing linear regression on the segment ring number and the segment base gauge offset deviation of each ring to obtain a linear regression function Y between the segment ring number and the segment base gauge offset deviationz＝δ0+δ1X + epsilon, wherein X is the independent variable of the number of the segment ring, delta0Is a constant number, δ1Is a regression coefficient, ε is the residual error of the rampart offset, YzCalculating a standoff distance deviation dependent variable for the segment base;
according to a linear regression function Y between a segment ring number and a segment rampart offset deviationz＝δ0+δ1X + epsilon obtains the residual error of the rampart offset deviation of each ring of the tube sheet ring;
step four, according to the formula
Figure FDA0002415344800000014
Calculating standard deviation sigma of residual error of each ring of pipe sheet ring rampart offsetXWherein, epsilonXIs the residual error of the rampart offset corresponding to the segment ring number, x is the accumulated ring number,
Figure FDA0002415344800000015
the residual mean value of the base gauge offset deviation under the accumulated ring number is obtained;
step five, according to the formula
Figure FDA0002415344800000016
Calculating residual error coefficient s of each ring of tube sheet ringX；
Step six, according to the formula
Figure FDA0002415344800000021
Calculating residual coefficient oscillation value ZD, wherein XqIs the initial segment ring number, sqIs the residual coefficient, X, corresponding to the ring number of the starting segmentJTo stop the ring number of the segment, sJIs the residual coefficient, s, corresponding to the ring number of the cut-off segmentmaxFor the maximum residual coefficient under the accumulated ring number, XmaxThe segment ring number corresponding to the maximum residual error coefficient under the accumulated ring number;
seventhly, shield tunnel line shape early warning: according to the formula
Figure FDA0002415344800000022
Calculating a shield segment splicing linear fitting error characteristic value P;
when the characteristic value P of the splicing linear fitting error of the shield segments is less than 50, the shield tunnel line shape is normal, and no early warning signal is sent out;
when the characteristic value of the shield segment splicing linear fitting error is more than or equal to 50 and less than 67, the splicing quality of the shield segments is controlled, and the shield tunnel linear early warning is yellow early warning;
when the characteristic value 67 of the splicing linear fitting error of the shield segments is not more than P and less than 100, indicating that the splicing quality of the shield segments deviates from a design route, and at the moment, the linear early warning of the shield tunnel is an orange early warning;
when the characteristic value P of the shield segment splicing linear fitting error is more than or equal to 100, the shield segment splicing quality is unqualified, and at the moment, the shield tunnel linear early warning is red early warning and needs to be shut down for on-site investigation.
2. The shield tunnel line-shaped early warning method according to claim 1, characterized in that: the actual half-wedge amount in the first step comprises an actual horizontal half-wedge amount and an actual vertical half-wedge amount, the actual horizontal half-wedge amount comprises an actual horizontal left half-wedge amount and an actual horizontal right half-wedge amount, and the actual vertical half-wedge amount comprises an actual vertical top half-wedge amount and an actual vertical bottom half-wedge amount; the actual accumulated rotation angle of the current tube sheet ring I comprises an actual horizontal accumulated rotation angle of the current tube sheet ring I and an actual vertical accumulated rotation angle of the current tube sheet ring I;
the actual rampart offset of the current tube sheet ring I in the step two comprises an actual horizontal rampart offset of the current tube sheet ring I and an actual vertical rampart offset of the current tube sheet ring I, and the design rampart offset of the current tube sheet ring I comprises a design horizontal rampart offset of the current tube sheet ring I and a design vertical rampart offset of the current tube sheet ring I;
in the third step, the residual error of the base gauge support distance deviation of each ring of the tube sheet ring comprises the residual error of the horizontal base gauge support distance deviation of each ring of the tube sheet ring and the residual error of the vertical base gauge support distance deviation of each ring of the tube sheet ring;
the standard deviation of the residual error of the offset of each ring of the tube sheet ring in the fourth step comprises the standard deviation of the residual error of the offset of each ring of the horizontal base of the tube sheet ring and the standard deviation of the residual error of the offset of each ring of the vertical base of the tube sheet ring;
the residual error coefficients of each ring of the tube sheet ring in the step five comprise a horizontal residual error coefficient of each ring of the tube sheet ring and a vertical residual error coefficient of each ring of the tube sheet ring;
the residual error coefficient oscillation value in the sixth step comprises a horizontal residual error coefficient oscillation value and a vertical residual error coefficient oscillation value;
and seventhly, performing shield tunnel line shape early warning on shield segment assembly horizontal error early warning and shield segment assembly vertical error early warning.
3. The shield tunnel line-shaped early warning method according to claim 2, characterized in that: in the step one, the current pipe sheet ring is a 16-point pipe sheet ring;
when the current pipe sheet ring is a standard ring, according to a formula
Figure FDA0002415344800000031
Calculating the half wedge amount on the right side of the actual level
Figure FDA0002415344800000032
And actual vertical bottom half wedge amount
Figure FDA0002415344800000033
When the current pipe sheet ring is a left-turning ring, according to a formula
Figure FDA0002415344800000034
Calculating the half wedge amount on the right side of the actual level
Figure FDA0002415344800000035
And actual vertical bottom half wedge amount
Figure FDA0002415344800000036
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure FDA0002415344800000037
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure FDA0002415344800000038
When the current pipe sheet ring is a right-turning ring, according to a formula
Figure FDA0002415344800000039
Calculating the half wedge amount on the right side of the actual level
Figure FDA0002415344800000041
And actual vertical bottom half wedge amount
Figure FDA0002415344800000042
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure FDA0002415344800000043
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure FDA0002415344800000044
4. The shield tunnel line-shaped early warning method according to claim 2, characterized in that: in the step one, the current pipe sheet ring is a 10-point pipe sheet ring;
when the current pipe sheet ring is a standard ring, according to a formula
Figure FDA0002415344800000045
Calculating the half wedge amount on the right side of the actual level
Figure FDA0002415344800000046
And actual vertical bottom half wedge amount
Figure FDA0002415344800000047
When the current pipe sheet ring is a left-turning ring, according to a formula
Figure FDA0002415344800000048
Calculating the half wedge amount on the right side of the actual level
Figure FDA0002415344800000049
And actual vertical bottom half wedge amount
Figure FDA00024153448000000410
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure FDA00024153448000000411
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure FDA00024153448000000412
When the current pipe sheet ring is a right-turning ring, according to a formula
Figure FDA00024153448000000413
Calculating the half wedge amount on the right side of the actual level
Figure FDA00024153448000000414
And actual vertical bottom half wedge amount
Figure FDA00024153448000000415
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure FDA0002415344800000051
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure FDA0002415344800000052
5. The shield tunnel line-shaped early warning method according to claim 2, characterized in that: in the step one, the current pipe sheet ring is a universal pipe sheet ring;
according to the formula
Figure FDA0002415344800000053
Calculating the half wedge amount on the right side of the actual level
Figure FDA0002415344800000054
And actual vertical bottom half wedge amount
Figure FDA0002415344800000055
Wherein, delta is the maximum wedge-shaped amount of the current segment ring design, thetayThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the leftmost point and the circle center of the segment is
Figure FDA0002415344800000056
n is the number of the current segment ring assembling points, thetadThe included angle between the connecting line of the current segment assembling point position and the circle center and the connecting line of the segment highest point and the circle center is
Figure FDA0002415344800000057
CN202010189479.3A 2020-03-18 2020-03-18 Shield tunnel linear early warning method Active CN111365076B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010189479.3A CN111365076B (en) 2020-03-18 2020-03-18 Shield tunnel linear early warning method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010189479.3A CN111365076B (en) 2020-03-18 2020-03-18 Shield tunnel linear early warning method

Publications (2)

Publication Number Publication Date
CN111365076A true CN111365076A (en) 2020-07-03
CN111365076B CN111365076B (en) 2021-05-25

Family

ID=71204656

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010189479.3A Active CN111365076B (en) 2020-03-18 2020-03-18 Shield tunnel linear early warning method

Country Status (1)

Country Link
CN (1) CN111365076B (en)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006009417A (en) * 2004-06-25 2006-01-12 Okumura Corp Tunnel lining correcting method
CN101446201A (en) * 2008-12-19 2009-06-03 上海隧道工程股份有限公司 Method for rectifying and selecting tunnel segment
CN102495961A (en) * 2011-12-05 2012-06-13 山东电力研究院 Simplified statistical control method of measurement process
CN105909270A (en) * 2016-04-19 2016-08-31 上海建科工程咨询有限公司 Shield machine axis control system
CN106089242A (en) * 2016-07-18 2016-11-09 广州地铁设计研究院有限公司 A kind of duct pieces of shield tunnel universality typesetting type selecting and the method for the most assembled type selecting
CN106767402A (en) * 2016-11-30 2017-05-31 华中科技大学 A kind of shield tunnel apparent mass detection method and system
CN106934159A (en) * 2017-03-14 2017-07-07 中交公局厦门工程有限公司 A kind of assembled point bit digitizing type selecting assembling method of duct pieces of shield tunnel
CN108397201A (en) * 2018-01-31 2018-08-14 中铁二十二局集团第工程有限公司 The safe construction method of rock tunnel(ling) machine
CN108868807A (en) * 2018-09-07 2018-11-23 上海隧道工程有限公司 The intelligent control method of shield driving correction
CN109164470A (en) * 2018-10-31 2019-01-08 西安思彼斯信息科技有限公司 The method and computation of real-time high-precision deformation monitoring in a kind of boring construction
CN109630132A (en) * 2018-11-21 2019-04-16 中铁十二局集团有限公司 Round jacked piles with small section outline excavation controls construction method
CN110135047A (en) * 2019-05-10 2019-08-16 上海隧道工程有限公司 Universal segment type selecting method for correcting error

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006009417A (en) * 2004-06-25 2006-01-12 Okumura Corp Tunnel lining correcting method
CN101446201A (en) * 2008-12-19 2009-06-03 上海隧道工程股份有限公司 Method for rectifying and selecting tunnel segment
CN102495961A (en) * 2011-12-05 2012-06-13 山东电力研究院 Simplified statistical control method of measurement process
CN105909270A (en) * 2016-04-19 2016-08-31 上海建科工程咨询有限公司 Shield machine axis control system
CN106089242A (en) * 2016-07-18 2016-11-09 广州地铁设计研究院有限公司 A kind of duct pieces of shield tunnel universality typesetting type selecting and the method for the most assembled type selecting
CN106767402A (en) * 2016-11-30 2017-05-31 华中科技大学 A kind of shield tunnel apparent mass detection method and system
CN106934159A (en) * 2017-03-14 2017-07-07 中交公局厦门工程有限公司 A kind of assembled point bit digitizing type selecting assembling method of duct pieces of shield tunnel
CN108397201A (en) * 2018-01-31 2018-08-14 中铁二十二局集团第工程有限公司 The safe construction method of rock tunnel(ling) machine
CN108868807A (en) * 2018-09-07 2018-11-23 上海隧道工程有限公司 The intelligent control method of shield driving correction
CN109164470A (en) * 2018-10-31 2019-01-08 西安思彼斯信息科技有限公司 The method and computation of real-time high-precision deformation monitoring in a kind of boring construction
CN109630132A (en) * 2018-11-21 2019-04-16 中铁十二局集团有限公司 Round jacked piles with small section outline excavation controls construction method
CN110135047A (en) * 2019-05-10 2019-08-16 上海隧道工程有限公司 Universal segment type selecting method for correcting error

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
丁凯等: "《量测》", 31 May 2014, 黄河水利出版社 *
朱旭: "地铁盾构施工中的若干测量手段及方法", 《中小企业管理与科技(下旬刊)》 *

Also Published As

Publication number Publication date
CN111365076B (en) 2021-05-25

Similar Documents

Publication Publication Date Title
CN105974398B (en) A kind of preceding leveling method of radar laser target calibration
US20230112991A1 (en) Method of high-precision 3d reconstruction of existing railway track lines based on uav multi-view images
CN102779345A (en) Point cloud precise registering method based on gravity center Euclidean distance
CN104596420B (en) The accurate measurement method of laser tracker measuring basis prism square center position
CN110044326A (en) Mountainous area highway application Trigonometric Leveling
WO2022262571A1 (en) System for automated measurement of levelness of end surface of tunnel ring
CN107167119B (en) Data processing method for projection deformation
CN110411375B (en) Three-dimensional imaging method based on passive millimeter wave/terahertz imaging technology
CN107179533A (en) A kind of airborne LiDAR systematic errors Self-checking method of multi-parameter
CN101413785A (en) Error compensation method of positioning system based on double-rotating laser plane transmitter network
CN113295049A (en) Carrier rocket aiming method and device
CN110441760B (en) Wide-range seabed topographic map expansion composition method based on prior topographic map
CN111365076B (en) Shield tunnel linear early warning method
CN109211183B (en) Corner intersection measuring method for long and large tunnel deep shaft
CN106066177A (en) Shield automatic guiding system measures manual review's method
CN112857329B (en) Existing railway turnout center measuring method and system, storage medium and electronic equipment
CN111721260B (en) High-precision light beam method settlement measurement method based on i-angle error correction of level gauge
CN109655048A (en) Total station Fast Fixed-point line-putting method
CN110440769B (en) Method for measuring longitudinal and transverse offsets of positioning line
CN104949692B (en) Multibeam sounding system installs the computational methods of corrected value
CN108362493A (en) A kind of numerically-controlled machine tool linear axis angular errors rapid detection method
CN111854715A (en) Pipe ring flatness measuring method based on total station
CN101209461A (en) Steel tube punching machine three-roller centring device center alignment degree on-line measurement method
CN111174790B (en) Method for forming topographic profile tracking path
CN108637037A (en) A kind of method of steel cold straightener verification straightening roll levelness

Legal Events

Date Code Title Description
PB01 Publication
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