CN102607506B  Free stationing transformation monitoring method of highfill airport side slope unit set total station  Google Patents
Free stationing transformation monitoring method of highfill airport side slope unit set total station Download PDFInfo
 Publication number
 CN102607506B CN102607506B CN 201210051569 CN201210051569A CN102607506B CN 102607506 B CN102607506 B CN 102607506B CN 201210051569 CN201210051569 CN 201210051569 CN 201210051569 A CN201210051569 A CN 201210051569A CN 102607506 B CN102607506 B CN 102607506B
 Authority
 CN
 China
 Prior art keywords
 centerdot
 survey station
 coordinate
 coordinate system
 monitoring
 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.)
 Expired  Fee Related
Links
 230000001131 transforming Effects 0.000 title abstract 2
 238000005259 measurement Methods 0.000 claims abstract description 42
 238000006243 chemical reaction Methods 0.000 claims description 14
 239000011159 matrix material Substances 0.000 claims description 10
 238000004364 calculation method Methods 0.000 claims description 9
 230000000875 corresponding Effects 0.000 claims description 8
 230000000694 effects Effects 0.000 claims description 4
 238000003672 processing method Methods 0.000 claims description 3
 235000020127 ayran Nutrition 0.000 claims description 2
 230000001276 controlling effect Effects 0.000 claims description 2
 229910052757 nitrogen Inorganic materials 0.000 claims description 2
 229910052698 phosphorus Inorganic materials 0.000 claims description 2
 230000014616 translation Effects 0.000 claims description 2
 238000000034 method Methods 0.000 description 14
 238000006073 displacement reaction Methods 0.000 description 9
 238000010276 construction Methods 0.000 description 8
 238000005516 engineering process Methods 0.000 description 5
 230000003068 static Effects 0.000 description 3
 208000000682 Megaloblastic Anemia Diseases 0.000 description 2
 238000010586 diagram Methods 0.000 description 2
 239000002689 soil Substances 0.000 description 2
 230000015572 biosynthetic process Effects 0.000 description 1
 230000002950 deficient Effects 0.000 description 1
 238000005755 formation reaction Methods 0.000 description 1
 230000002706 hydrostatic Effects 0.000 description 1
 230000001360 synchronised Effects 0.000 description 1
Abstract
The invention relates to a free stationing transformation monitoring method of a highfill airport side slope unit set total station, and is used for effectively solving the accurate monitoring problems that a highfill airport side slope is long and narrow, is a strip shape and is linear extension, and the graph intensity is poor. The monitoring method comprises the following steps of: dividing the involved range into a first measuring region, a second measuring region and a third measuring region from large to small, and setting a control point, a monitoring point and a turning point in each measure region, wherein the position of a survey station is freely selected; establishing a control net of the whole monitoring region, carrying out substation measurement on the control points and the monitoring points in the measuring regions, establishing a posture and oriented relation between adjacent survey stations, meanwhile, realizing larger expansion in an instrument measurement range; and taking the acquired horizon angle, a vertical angle and slant distance as observed values, taking a survey station position parameter and a monitoring point coordinate as an adjustment parameter so as to carry out parameter adjustment, and obtaining the precision result of the survey station position parameter and the monitoring point coordinate. The monitoring method provided by the invention has the advantages of high measuring accuracy, flexibility and convenience, convenience for carrying and good portability.
Description
Technical field
The present invention relates to the Aerodrome Construction engineering, particularly the Free Station deformation monitoring method of a kind of high embankment airport side slope separate unit total powerstation.
Background technology
Because aviation development, air resource is more and more nervous, certainly will cause airport addressing to a certain extent " the empty ground that determines ", because State owned land is tightly controlled policy, certainly will cause airport construction " to be climbed mountains and gone to sea ", airport construction has to be chosen on the geology unfavorable foundation thus, such as the soft soil foundation of southeastern coast, westsouthwest Mountainous Region high fill foundation etc., is commonly referred to high embankment airport on the airport that this class ground is built.The southeastern coast soft soil foundation in order to satisfy the designing requirement of flood frequency flood control level, need to carry out embankment usually; The mountain area high fill foundation, because mistake, cheuch are in length and breadth between alpine terrain undulation, river, mound, must cut the mountain and fill out gully, carry out on the spot a cubic meter balance, form filling body, because the formation of filling body has destroyed original equilibrium state on natural ground, filling body itself also deforms owing to many factors simultaneously, the form of expression is: settlement of foundation, slope displacement etc.
Because the anisotropy of distortion genesis mechanism and the uncertainty that Deformation Theory is calculated, Stability Checking Calculation is analyzed, deformation monitoring is as the only resource of the substantial amount that obtains to deform, safety evaluatio for engineering seems particularly important, it also is important means informationbased, scientific construction, can verify on the one hand the rationality of ground Treatment Design scheme, can be on the other hand the foundation that next step construction provides science.
Existing deformation monitoring method is divided into horizontal displacement monitoring, perpendicular displacement monitoring and threedimensional monitoring by the observation purpose:
Horizontal displacement monitoring has following several method: method of tension wire alignment, and collimation line method, method of laser alignment, just, the reversed pendulum method, forward intersection and precise traverse method etc.The perpendicular displacement monitoring mainly contains geometric leveling method and hydrostatic levelling method (communicating pipe method).More than various monitoring methods be with the horizontal shift of deformation point and respectively testing of perpendicular displacement.Along with the development of surveying instrument and measuring technique, this problem has obtained solution basically.Developed at present the measuring system of the observation deformation point horizontal shift of energy realtime continuous and perpendicular displacement, because these systematic surveys is the threeD displacement value of deformation point, so be called " threeD displacement monitoring system ", can be divided into GPS method, Distance Intersection method, polar coordinates method of difference and forward intersection by its principle and observation procedure.Although static Relative Difference method is feasible, signal is affected seriously by high slope GPS (comprising " the many antennas of a machine " technology), and instrument is expensive simultaneously, and is excessive for unitem cost input, lays under the multiple spot condition in net, and observed efficiency is very low.Distance Intersection method and polar coordinates method of difference all need special measuring system to cooperate, and the whole system expense is higher.In the deformation monitoring of high embankment airport side slope, use most often forward intersection, the error source of the method has: angle error, the variation of intersection angle and graphic structure, base length, external condition etc.The method labour intensity is large in general, the measurements and calculations process is complicated, computational accuracy is undesirable.The reason of searching to the bottom is to have occurred the artificial error in many places in observation process, such as the error of centralization of instrument and reflecting prism, Setting error, instrument high measurement error, prism high measurement error etc.
Said method exists: observation condition institute is restricted greatly, the monitoring system Meteorological is higher, the accidental error source is many, the Measurement and Data Processing process complicated, and its Monitoring Result does not have realtime usually.
Simultaneously, the Free Station with Total Station electronic Thacheometer method of mentioning in the pertinent literature is only measured two reference mark, and there is following problem in composition control net not:
(1), the distribution method of control net lacks the assessment of the strength of figure factor;
(2), lack the robust estimation of control net in its initial observation adjustment result, control the antirough error capabilities of net;
(3), do not carry out deciding power in the Monitoring Data adjustment process, effect is bad, and measuring accuracy is poor, does not satisfy the actual needs of airport construction, and therefore, the improvement and bring new ideas on the monitoring method is imperative.
Summary of the invention
For abovementioned situation, for overcoming the defective of prior art, the present invention's purpose just provides the Free Station deformation monitoring method of a kind of high embankment airport side slope separate unit total powerstation, it is poor effectively to solve measuring accuracy, can not satisfy high embankment airport side slope and be that long and narrow, ribbon, straight line extend, the poor monitoring problem of strength of figure.
The technical scheme that the present invention solves is may further comprise the steps:
1, high embankment airport side slope subregion
The scope that high embankment airport side slope is comprised is divided into from large to small first, second, third and establishes measured rectangular area, station, and the length of side of each rectangular area is 200300m; Establish reference mark, monitoring point and turning point in each measured zone, the position of survey station can freely be selected, but is located at the centre of measured zone as far as possible, so that measuring error is regular even distribution;
2, measure
The resulting coordinate system of control survey is called control net coordinate system, and every station formed coordinate system of surveying instrument three axles is called the survey station coordinate system, measures in two steps:
(1) control survey: before deformation monitoring begins, should set up the control net of whole monitored area, the establishment of control net mode is GPS free net or classical triangulateration network, adopt rectangular space coordinate, shorter for the common length of side of deformation monitoring, do not carry out projection, with the error effect of avoiding distortion of projection to be produced;
(2) deformation monitoring: Free Station is measured, and substation is carried out in the reference mark in the measured zone and monitoring point measure, and method is, and the one, single station is measured, and instrument is motionless, measures all targets, the zone that suitable monitoring range is less; The 2nd, turn the station and measure (be commonly referred to and leapfrog), finishing the oneshot measurement task needs the repeatedly position of mobile instrument, can avoid the impact of external environment, improves sighting condition; By the orientation point more than 3 (being the common point of adjacent two station duplicate measurementss) is measured, set up attitude, directional relation between the adjacent survey station, realize simultaneously the larger expansion of apparatus measures scope, also can avoid measuring accuracy to increase and fast reducing with distance, be fit to measured zone in a big way;
3, data are processed
The horizontal angle, vertical angle and the oblique distance that measure to obtain as observed reading, are carried out parameter adjustment with survey station location parameter and monitoring point coordinate as adjustment parameter, obtain the precise results of survey station location parameter and monitoring point coordinate.
Monitoring method precision of the present invention is high, the measuring system that consists of is convenient, flexible, volume is little, lightweight, be convenient to move, have good portability, adapted to fully that slope monitoring place, airport is long and narrow, construction disturbs frequently, needs flexible, the motordriven needs of establishing the station.
Description of drawings
Fig. 1 is Free Station deformation monitoring schematic diagram of the present invention.
Fig. 2 is that common point is measured sign picture between adjacent 2 survey station points of the present invention.
Fig. 3 is polar coordinate measurement sign picture of the present invention.
Fig. 4 is the transformational relation figure between survey station of the present invention and control net.
Fig. 5 is certain embodiment monitoring point, airport distribution plan of the present invention.
Fig. 6 is choppedoff head control net net type figure of the present invention.
Fig. 7 is encryption control net net type figure of the present invention.
Embodiment
Below in conjunction with drawings and Examples the specific embodiment of the present invention is elaborated.
The present invention is realized by following steps in implementation:
1, high embankment airport side slope subregion
As shown in Figure 1, be separate unit Free Station with Total Station electronic Thacheometer deformation monitoring method schematic diagram, its intermediate cam represents the reference mark, circle represents the monitoring point, square frame represents turning point, the instrument icon represents the survey station position, and the scope that high embankment airport side slope is comprised is divided into from large to small first, second, third and establishes measured rectangular area, station, and the length of side of each rectangular area is 200300m; If monitoring range is large again, measure the measured zone of separating and continue to increase, each survey station is measured all reference mark, monitoring point and the turning point in its scope, and the position of survey station can freely be selected, but as far as possible in the centre of measured zone, can make like this measuring error be regular even distribution;
The total station survey precision is relevant with angle measurement and range error, the error of angle measurement and range finding is less on the positional accuracy impact in the time of closely, but along with distance increases, the range error impact obviously increases, the Free Station deformation monitoring method is divided into the zone that less, the suitable total powerstation of several scopes is monitored with larger monitored area, thereby effectively guaranteed measuring accuracy, fully adapted to high embankment airport side slope and be that long and narrow, ribbon, straight line extend, the poor situation of strength of figure;
2, measure
The resulting coordinate system of control survey is called control net coordinate system, and every station formed coordinate system of surveying instrument three axles is called the survey station coordinate system, measures in two steps:
(1) control survey: before deformation monitoring begins, set up first the control net of whole monitored area, the establishment of control net mode is GPS free net or classical triangulateration network, adopt rectangular space coordinate, shorter because of the length of side for deformation monitoring, do not carry out projection, with the error effect of avoiding distortion of projection to be produced;
(2) deformation monitoring, Free Station are measured, and substation is carried out in the reference mark in the measured zone and monitoring point measure, and method is, when measured zone is less, adopt single station to measure, and instrument is motionless, measures all targets, and measuring radius is in the 300m; When measured zone in a big way, employing turns the station and measures (be commonly referred to and leapfrog), finish the oneshot measurement task and need the repeatedly position of mobile instrument, can avoid the impact of external environment, improve sighting condition, by the orientation point more than 3 (being the common point of adjacent two station duplicate measurementss) is measured, set up attitude, directional relation between the adjacent survey station, realize simultaneously the larger expansion of apparatus measures scope, also can avoid measuring accuracy to increase and fast reducing with distance;
3, data are processed
The horizontal angle, vertical angle and the oblique distance that measure to obtain as observed reading, are carried out parameter adjustment with survey station location parameter and monitoring point coordinate as adjustment parameter, obtain the precise results of survey station location parameter and monitoring point coordinate; Said process is similar to the bundle adjustment in photogrammetric, each survey station parameter is similar to the inside and outside parameter in photogrammetric, the coordinate of monitoring point then is the unknown point coordinate in photogrammetric, the data processing method that the many survey stations data unification that therefore is referred to as to resolve based on bundle adjustment is resolved;
A, net adjustment
The net adjustment is after observation finishes, eliminate incongruent data, assess measuring accuracy, obtain the important means of coordinate, traditional net adjustment can be divided into the adjustment of condition equation, parameter adjustment two large classes, parameter adjustment (reaching the parameter adjustment with constraint condition) is convenient to computer aided calculation, specifically (as shown in Figure 2):
Two instruments are observed 1 point simultaneously, and 6 observed readings are arranged, if to the simultaneously observation of n point, 6n observed reading just arranged; There is 3n unknown number in n unknown point, two instrument relative orientations have other 7 unknown numbers (3 rotation parameters, 3 translation parameterss, 1 scale factor), if resolve the relative position of instrument, should be so that following condition establishment: 6n＞7+3n, be n＞2.5, for this reason during actual measurement, require adjacent survey station that common point more than 3 is arranged, in Practical Project, highprecision monitoring net can be laid 5～15 common points usually, increase excess observation component, and network point distribution has good geometry, to improve reliability, to reduce the impact of measuring error;
Process by excess observation being carried out least square during adjustment, try to achieve best instrument position and the best coordinates of spatial attitude and spatial point, so that the quadratic sum of observed reading correction is minimum, its error equation is nonlinear, need to carry out repeatedly iteration and just can reach final requirement;
B, data processing method
(1) coordinate Calculation (as shown in Figure 3) under the survey station coordinate system
Survey station i comprises horizontal angle H to the observed reading of arbitrfary point P
_{Ip}, vertical angle V
_{Ip}With oblique distance D
_{Ip}, then put the coordinate of P under survey station i and be denoted as (X
_{IP}, Y
_{IP}, Z
_{IP}), computing formula is as follows:
(2) the conversion parameter summary between survey station 1 coordinate system and all the other survey stations calculates
Can calculate the conversion parameter between two adjacent survey station i and i+1, i.e. translation parameters according to boolean Sha seven parameter models (realizing by 3 basic rotations of coordinate system, 3 translations and 1 yardstick convergentdivergent)
Rotation parameter
With scale factor k:
(X
_{i}, Y
_{i}, Z
_{i}) be the coordinate of survey station i, (X
_{I+1}, Y
_{I+1}, Z
_{I+1}) be survey station i+1 coordinate,
Be the angle of rotating around X, Y, Z axis successively, N
_{Ii+1}It is rotation parameter
Corresponding rotation matrix;
The transformational relation that can be got between survey station 1 and survey station i by (2) formula is:
Be denoted as:
Wherein:
Scale factor k mainly is because two coordinate systems adopt different length standards to cause, perhaps the testee factor such as expand with heat and contract with cold causes, if the length standard of two coordinate systems is identical, scale factor k is fixed as 1, otherwise can be calculated by following formula the summary value of k:
Wherein: (x
_{i}, y
_{i}, z
_{i}) be the coordinate under the survey station coordinate system, (x "
_{i}, y "
_{i}, z "
_{i}) for controlling the coordinate under the net coordinate system;
(3) transformational relation calculates between each survey station coordinate system and control net coordinate system
By (3) formula, can with all measurement point coordinate conversion under survey station 1 coordinate system, to measure at least 3 not reference mark of conllinear during monitoring, obtain the summary conversion parameter between control net coordinate system and survey station 1 coordinate system, i.e. translation parameters
And rotation parameter
Then have:
Wherein: (X
_{c}, Y
_{c}, Z
_{c}) coordinate under the expression control net coordinate system, (X
_{1}, Y
_{1}, Z
_{1}) be the coordinate under survey station 1 coordinate system, N
_{C1}It is rotation parameter
Corresponding rotation matrix is namely controlled the net coordinate and is tied to survey station 1;
By the transformational relation between (3) and (5) Shi Kede control net coordinate system and At any points i:
Wherein: (X
_{i}, Y
_{i}, Z
_{i}) coordinate under the expression survey station i coordinate system,
Be the translation parameters of survey station 1 to survey station i, N
_{1i}It is rotation parameter
Corresponding rotation matrix, namely survey station 1 to survey station i;
(4) accurately resolve position relationship and monitoring point coordinate between survey station
For survey station i, can be calculated the summary value of the conversion parameter between control net coordinate system and survey station coordinate system by (6) formula
In the actual measurement process, the reference field of control net coordinate system and the reference field of survey station coordinate system are geoid surface, so the conversion parameter between two coordinate systems will only have the rotation angle around Z axis
As shown in Figure 4, O wherein
_{C}X
_{C}Y
_{C}Z
_{C}Be control net coordinate system, O
_{i}X
_{i}Y
_{i}Z
_{i}Be survey station i coordinate system, O '
_{C}X '
_{C}Y '
_{C}Z '
_{C}Coordinate system is by coordinate system O
_{C}X
_{C}Y
_{C}Z
_{C}Move to O
_{i}, plane O
_{C}X
_{C}Y
_{C}, O '
_{C}X '
_{C}Y '
_{C}And O
_{i}X
_{i}Y
_{i}All parallel with surface level, O then
_{i}X
_{i}Y
_{i}Z
_{i}By coordinate system O '
_{C}X '
_{C}Y '
_{C}Z '
_{C}Rotation obtains the rotation angle around Z axis
The measured point is divided into reference mark and monitoring point, and the coordinate of monitoring point j under control net coordinate system is denoted as (X
_{Cj}, Y
_{Cj}, Z
_{Cj}), for whole monitoring net, the unknown number that find the solution comprises the conversion parameter between survey station and control net coordinate system
And the unknown number coordinate (X of monitoring point
_{Ck}, Y
_{Ck}, Z
_{Ck}); Survey station i comprises horizontal angle Hip, vertical angle Vip and oblique distance Dip to the observed reading of measurement point P, and following relational expression is arranged:
Then can get following error equation by parameter adjustment:
Order:
Wherein: X
_{Ck}
_{0},
Y
_{Ck}
_{0},
Z
_{Ck}
_{0},
Be X
_{Ck},
Y
_{Ck},
Z
_{Ck},
Corresponding approximate value, and the unit of angle adopts second, and the unit of length adopts millimeter, can avoid like this error equation coefficient difference larger, and following relational expression is arranged:
Wherein: H
_{Ip}
_{0}Be H
_{Ip}Corresponding approximate value, then:
If observation station be the reference mark then
${a}_{3}^{ip}={a}_{4}^{ip}={b}_{4}^{ip}={b}_{5}^{ip}={b}_{6}^{ip}={c}_{4}^{ip}={c}_{5}^{ip}={c}_{6}^{ip}=0;$
For all survey stations, can set up error equation V=AX+L, observed reading horizontal angle Hip, vertical angle Vip and oblique distance Dip to decide the power expression formula as follows:
Wherein: m
_{β}Be the nominal accuracy of instrument angle measurement, unit is second; m
_{D}Be distance accuracy, unit is millimeter; A, b are respectively fixed error and the proportional error coefficient of total powerstation nominal;
By least square method V
^{T}PV=is minimum, forms equation:
$\left\{\begin{array}{c}\mathrm{NX}+W=0\\ X={N}^{1}W\end{array}\right.,$ Wherein
$\left\{\begin{array}{c}N={A}^{T}\mathrm{PA}\\ W={A}^{T}\mathrm{PL}\end{array}\right.,$ Thus, adjustment is resolved and can be obtained the survey station unknown number
With monitoring point coordinate (X
_{Ck}, Y
_{Ck}, Z
_{Ck}) exact value;
(5) accuracy assessment
If n is the error equation number, t is the unknown number number, m
_{i}It is the precision valuation of i parameter; Q
_{Ii}Be the capable i column data of i on the inverse of weight matrix Q diagonal line; By adjustment result residual error V, but error in the unit of account power
Wherein
Unknown number inverse of weight matrix Q=(A
^{T}PA)
^{1}, parameters precision is estimated
C, accuracy comparison test (being combined in the application in the specific embodiment, the reliability of sufficient proof the inventive method and accuracy)
(1) subregion is layouted
Certain airport is positioned at China the Northeast, belongs to Airport in Mountain Region, 22 meters of average height of embankments, 35 meters of maximum height of embankments, side slope is than 1: 0.75, more than 800 meters of whole length of slopes, as shown in Figure 5, monitoring section is laid 10 at reference mark altogether, is distributed in around the whole deformed slope, and period is K01K10; Lay altogether 42 by the designing requirement monitoring point, minute side slope top, centre and bottom 3 rows, 14 row are laid, and period is B01B42, and row are spaced apart 50 meters;
(2) control survey
Because the GPS technology has bearing accuracy height, good in economic efficiency, the easy and simple to handle and characteristics such as degree of freedom is large of arranging net, military affairs at home and abroad, the department such as civilian have obtained using very widely, especially in geodetic surveying and deformation monitoring field, publish from Chinese Geology Publishing House, Liu Ping write " newly organized mapping and measurement data quick checking disposal route checks in quick checking technical application handbook: the monitoring point horizontal direction precision that the 6hGPS Monitoring Data is resolved in the Geheyan Reservoir dam monitoring is ± 0.5mm, and the elevation directional precision is ± 1.0mm; The monitoring point horizontal direction precision that 1～2hGPS Monitoring Data is resolved is ± 1.0mm that the elevation directional precision is ± 1.5mm;
According to abovementioned situation, adopt the GPS method, divided for three steps carried out: the high precision Fixed Initial Point is measured, measurement is netted in choppedoff head control and encrypt the measurement of control net, in the GPS processing procedure, be divided into highprecision GPS measurement, conventional GPS measurement, the subject matter of its solution is: the foundation of control net, the monitoring of control net absolute deformation, the monitoring of control net relative deformation;
High precision Fixed Initial Point and choppedoff head control net are measured:
Choppedoff head is controlled net net type as shown in Figure 6, and it is as follows that high precision Fixed Initial Point and choppedoff head control net are measured the employing scheme:
(a) according to actual conditions, form this project choppedoff head control net by K001, K005, K006, K010, selection K010 is reference point, K001, K005, K006 are point to be located;
(b) translocation BJFS (BeiJing, China), DAEJ (Taejon, Korea), the international tracking station of 3 IGS of SUWN (Suwon, Korea) adopt the final sophisticated product of IGS, the final sophisticated product of CODE;
(c) adopt Bernese 5.0 softwares to carry out data and process, calculation method is two poor phase place static solutions, short Baselines pattern;
(d) achievement adopts geocentric coordinate X, Y, Z;
(e) all the other observation technologies adopt the requirement of National GPS specifications of surveys C level network technology;
Encrypting the control net measures:
Detail network net type is encrypted the measurement of control net and is adopted conventional GPS to measure static Relative Difference method as shown in Figure 7, and scheme is as follows:
(a) resolve software and adopt Leica LGO 7.0;
(b) grid connect mode adopts Bian Lianshi;
(c) requirement of National GPS specifications of surveys D level network technology is adopted in all the other observations;
Flow process is resolved in control survey:
(a) the fixing calculating coordinate K010 of BJFS tracking station reference point coordinate;
(b) fixedly each point coordinate is netted in K010 calculating coordinate K001, K005, the control of K006 choppedoff head;
(c) adopt K010, K001, K005, K006 geocentric coordinate directly to retrain adjustment and resolve each point coordinate in the detail network;
(d) geocentric coordinate is transformed into the airport coordinate system;
Final calculation result is as shown in table 1:
Table 1 airport coordinate system calculation result
Call the roll  Type  X/m  Y/m  Z/m  RMS/mm 
K001  Adjusted value  2409.1665  3642.0733  45.3931  0.05 
K002  Adjusted value  2412.8010  3515.4494  36.2717  Fixed value 
K003  Adjusted value  2408.1534  3315.6843  32.8461  0.06 
K004  Adjusted value  2392.0331  3173.7755  35.9157  0.06 
K005  Adjusted value  2255.2726  3112.0514  33.8521  Fixed value 
K006  Adjusted value  2292.3103  2910.9994  46.2114  Fixed value 
K007  Adjusted value  2518.8474  3084.7063  47.3788  0.06 
K008  Adjusted value  2523.6450  3189.1383  48.4965  0.04 
K009  Adjusted value  2508.7601  3534.8179  49.5865  0.04 
K010  Adjusted value  2517.5188  3735.8074  59.5505  Fixed value 
(3) deformation monitoring
The deformation monitoring scheme adopts Leica company's T CA1800 type total powerstation (as shown in Figure 5), lay altogether 4 survey stations, the scope that different big or small square frames surround is respectively first, second, the 3rd, the 4th establishes measured rectangular area, station, the length of side is 200～300 meters, if monitoring range is large again, measuring the zone of separating continues to increase, each is established the station and measures all interior reference mark of its scope, monitoring point and interim turning point (can be laid interim turning point when the interval, monitoring point is larger, strengthen network structure), at least duplicate measurements two row points between adjacent two stations, i.e. at least 6 common points, the position of survey station can freely be selected, but as far as possible in the centre of measured zone, can make like this measuring error be regular even distribution, realtime Accurate Measurement temperature and air pressure before every station is measured, and correct in instrument;
(4) resolve
First measurement result is as shown in table 2, and Using such method can obtain the result of each issue observation;
The first monitoring point coordinate of table 2 and precision (fragment)
Simultaneously, in order to verify the correctness of the inventive method, also the result of total powerstation and GPS measurement contrasted, measurement result such as table 3, comparing result is as shown in table 4:
Table 3GPS and total station survey result
Table 4GPS and total station survey result are relatively poor
Period  dX/mm  dY/mm  dH/mm  Plane deviation/mm  The point biased poor/mm 
K02  1.3  2.0  1.9  2.4  3.0 
K04  2.7  0.7  0.3  2.8  2.8 
K08  0.0  2.2  0.0  2.2  2.2 
K09  2.8  0.7  1.9  2.9  3.5 
Mean value  1.7  1.4  1.0  2.6  2.9 
Can find out from above result of calculation, use the Free Station method and be 2.9mm with the mean value of the biased difference of GPS method measurement point, the mean value of planimetric coordinates deviation is 2.6mm, comprehensively it seems, GPS method and Free Station deformation monitoring method measured result grid deviation in this paper are in 3.0mm, can satisfy the accuracy requirement of deformation measurement secondary in " building deformation measurement standard " (JGJ 82007) fully, namely satisfy the monitoring requirement of high embankment airport side slope.
Clearly show by abovementioned situation, the present invention is directed to the deformation monitoring problem of high embankment airport side slope, on the basis of summing up classic method, Free Station deformation monitoring method based on the separate unit total powerstation has been proposed, the measuring system that the method forms, volume is little, lightweight, be convenient to move, have good portability, adapted to fully that slope monitoring place, airport is long and narrow, construction is disturbed frequent, need flexible, motordriven characteristics of establishing the station, also carried out ratio of precision pair with the GPS method, proved absolutely the correctness of the inventive method, the inventive method has the following advantages:
(1) laying at reference mark requires lowly, and single survey station does not require the measurement number at reference mark, as long as whole survey stations can observe 4 reference mark altogether.
(2) Free Station, alignment error of instrument, amount high level error have been avoided, guarantee measurement of higher degree precision, the survey station position freely arranges, and namely instrument does not need to be centered on some fixing reference mark, has avoided the error of centralization of instrument, do not need to measure the antenna height of instrument in the measuring process, height can be any, is as the criterion with suitable observation, avoided height of instrument to measure error.
(3) threedimensional measurement, make plane and the synchronization gain of measurement of higher degree data, improve observed efficiency, traditional airport Monitoring of Slope Deformation method is processed respectively by plane and elevation, because plane, dissimilar instruments is adopted in the measurement of higher degree, different working specifications and computing method, so that measuring process is more loaded down with trivial details, inefficiency, sometimes unavoidably measuring accuracy is caused damage, comparatively speaking, threedimensional measurement then utilizes a kind of surveying instrument to obtain simultaneously threedimensional coordinate, measurement data is synchronous, be the innovation on the slope monitoring method of high embankment airport, extend to highway, railway, use on the correlation engineering High Slope Monitoring such as dam, huge economic and social benefit is arranged.
Claims (1)
1. the Free Station deformation monitoring method of a high embankment airport side slope separate unit total powerstation is characterized in that, may further comprise the steps:
One, high embankment airport side slope subregion
The scope that high embankment airport side slope is comprised is divided into from large to small first, second, third and establishes measured rectangular area, station, and each regional length of side is 200300m; If monitoring range is large again, measure the measured zone of separating and continue to increase, each survey station is measured all reference mark, monitoring point and the turning point in its scope, and the position of survey station can freely be selected, in the centre of measured zone, make measuring error be regular even distribution as far as possible;
Two, measure
The resulting coordinate system of control survey is called control net coordinate system, and every station formed coordinate system of surveying instrument three axles is called the survey station coordinate system, measures in two steps:
(1) control survey: before deformation monitoring begins, set up first the control net of whole monitored area, the establishment of control net mode is GPS free net or classical triangulateration network, adopt rectangular space coordinate, shorter because of the length of side for deformation monitoring, do not carry out projection, with the error effect of avoiding distortion of projection to be produced;
(2) deformation monitoring, Free Station are measured, and substation is carried out in the reference mark in the measured zone and monitoring point measure, and method is, when measured zone is less, adopt single station to measure, and instrument is motionless, measures all targets, and measuring radius is in the 300m; When measured zone in a big way, employing turns the station and measures, finish the oneshot measurement task and need the repeatedly position of mobile instrument, can avoid the impact of external environment, improve sighting condition, by the orientation point more than 3 is measured, set up attitude, directional relation between the adjacent survey station, realize simultaneously the larger expansion of apparatus measures scope, avoid measuring accuracy to increase and fast reducing with distance;
Three, data are processed
The horizontal angle, vertical angle and the oblique distance that measure to obtain as observed reading, are carried out parameter adjustment with survey station location parameter and monitoring point coordinate as adjustment parameter, obtain the precise results of survey station location parameter and monitoring point coordinate;
A, net adjustment
The net adjustment is after observation finishes, eliminates incongruent data, assesses measuring accuracy, obtains coordinate, specifically:
Two instruments are observed 1 point simultaneously, and 6 observed readings are arranged, if to the simultaneously observation of n point, 6n observed reading just arranged; There is 3n unknown number in n unknown point, and two instrument relative orientations have other 7 unknown numbers, 3 rotation parameters, 3 translation parameterss, 1 scale factor, according to 6n＞7+3n, n＞2.5, resolve the relative position of instrument, adjacent survey station has the common point more than 3, and highprecision monitoring net is laid 5～15 common points, increases excess observation component, and network point distribution has geometry, to improve reliability, to reduce the impact of measuring error;
By excess observation is carried out Least Square in Processing, try to achieve the coordinate of instrument position and spatial attitude and spatial point during adjustment, so that the quadratic sum of observed reading correction is minimum, its error equation is nonlinear, through repeatedly iteration realization;
B, data processing method
(1) coordinate Calculation under the survey station coordinate system
Survey station i comprises horizontal angle H to the observed reading of arbitrfary point P
_{Ip}, vertical angle V
_{Ip}With oblique distance D
_{Ip}, then put the coordinate of P under survey station i and be denoted as (X
_{IP}, Y
_{IP}, Z
_{IP}), computing formula is as follows:
(2) the conversion parameter summary between survey station 1 coordinate system and all the other survey stations calculates
According to boolean Sha seven parameter models, by 3 basic rotations of coordinate system, 3 translations and 1 yardstick convergentdivergent, calculate the conversion parameter between two adjacent survey station i and i+1, comprise translation parameters
Rotation parameter
With scale factor k:
(X
_{i}, Y
_{i}, Z
_{i}) be the coordinate of survey station i, (X
_{I+1}, Y
_{I+1}, Z
_{I+1}) be survey station i+1 coordinate,
Be the angle of rotating around X, Y, Z axis successively, N
_{Ii+1}It is rotation parameter
Corresponding rotation matrix;
Transformational relation between survey station 1 and survey station i is:
Be denoted as:
Wherein:
Scale factor k mainly is because two coordinate systems adopt different length standards to cause, and perhaps the testee factor such as expand with heat and contract with cold causes, if the length standard of two coordinate systems is identical, scale factor k is fixed as 1, is calculated the summary value of k by following formula:
Wherein: (x
_{i}, y
_{i}, z
_{i}) be the coordinate under the survey station coordinate system, (x "
_{i}, y "
_{i}, z "
_{i}) for controlling the coordinate under the net coordinate system;
(3) transformational relation calculates between each survey station coordinate system and control net coordinate system
By
$({X}_{1},{Y}_{1},{Z}_{1})={N}_{1i}\·({X}_{i},{Y}_{i},{Z}_{i})+({X}_{0}^{1i},{Y}_{0}^{1i},{Z}_{0}^{1i})$ Formula under survey station 1 coordinate system, will measure at least 3 not reference mark of conllinear with all measurement point coordinate conversion during monitoring, obtain the summary conversion parameter between control net coordinate system and survey station 1 coordinate system, translation parameters
And rotation parameter
Then have:
Wherein: (X
_{c}, Y
_{c}, Z
_{c}) coordinate under the expression control net coordinate system, (X
_{1}, Y
_{1}, Z
_{1}) be the coordinate under survey station 1 coordinate system, N
_{C1}It is rotation parameter
Corresponding rotation matrix is namely controlled the net coordinate and is tied to survey station 1;
By
$({X}_{1},{Y}_{1},{Z}_{1})={N}_{1i}\·({X}_{i},{Y}_{i},{Z}_{i})+({X}_{0}^{1i},{Y}_{0}^{1i},{Z}_{0}^{1i})$ With
$({X}_{c},{Y}_{c},{Z}_{c})={N}_{c1}\·({X}_{1},{Y}_{1},{Z}_{1})+({X}_{0}^{c1},{Y}_{0}^{c1},{Z}_{0}^{c1})$ Formula must be controlled the transformational relation between net coordinate system and At any points i:
Wherein: (X
_{i}, Y
_{i}, Z
_{i}) coordinate under the expression survey station i coordinate system,
Be the translation parameters of survey station 1 to survey station i, N
_{1i}It is rotation parameter
Corresponding rotation matrix, namely survey station 1 to survey station i;
(4) accurately resolve position relationship and monitoring point coordinate between survey station
For survey station i, by
$({X}_{c},{Y}_{c},{Z}_{c})={N}_{c1}\·{N}_{1i}\·({X}_{i},{Y}_{i},{Z}_{i})+{N}_{c1}\·({X}_{0}^{1i},{Y}_{0}^{1i},{Z}_{0}^{1i})+({X}_{0}^{c1},{Y}_{0}^{c1},{Z}_{0}^{c1})$ Formula calculates the summary value of the conversion parameter between control net coordinate system and survey station coordinate system
In measurement, the reference field of control net coordinate system and the reference field of survey station coordinate system are geoid surface, so the conversion parameter between two coordinate systems will only have the rotation angle around Z axis
O wherein
_{C}X
_{C}Y
_{C}Z
_{C}Be control net coordinate system, O
_{i}X
_{i}Y
_{i}Z
_{i}Be survey station i coordinate system, O '
_{C}X '
_{C}Y '
_{C}Z '
_{C}Coordinate system is by coordinate system O
_{C}X
_{C}Y
_{C}Z
_{C}Move to O
_{i}, plane O
_{C}X
_{C}Y
_{C}, O '
_{C}X '
_{C}Y '
_{C}And O
_{i}X
_{i}Y
_{i}All parallel with surface level, O then
_{i}X
_{i}Y
_{i}Z
_{i}By coordinate system O '
_{C}X '
_{C}Y '
_{C}Z '
_{C}Rotation obtains the rotation angle around Z axis
The measured point is divided into reference mark and monitoring point, and the coordinate of monitoring point j under control net coordinate system is denoted as (X
_{Cj}, Y
_{Cj}, Z
_{Cj}), for whole monitoring net, the unknown number that find the solution comprises the conversion parameter between survey station and control net coordinate system
And the unknown number coordinate (X of monitoring point
_{Ck}, Y
_{Ck}, Z
_{Ck}); Survey station i comprises horizontal angle Hip, vertical angle Vip and oblique distance Dip to the observed reading of measurement point P, and following relational expression is arranged:
Then can get following error equation by parameter adjustment:
Order:
Wherein: X
_{Ck}
_{0},
Y
_{Ck}
_{0},
Z
_{Ck}
_{0},
Be X
_{Ck},
Y
_{Ck},
Z
_{Ck},
Corresponding approximate value, and the unit of angle adopts second, and the unit of length adopts millimeter, can avoid like this error equation coefficient difference larger, and following relational expression is arranged:
Wherein: H
_{Ip}
_{0}Be H
_{Ip}Corresponding approximate value, then:
If observation station be the reference mark then
${a}_{3}^{ip}={a}_{4}^{ip}={b}_{4}^{ip}={b}_{5}^{ip}={b}_{6}^{ip}={c}_{4}^{ip}={c}_{5}^{ip}={c}_{6}^{ip}=0;$
For all survey stations, can set up error equation V=AX+L, observed reading horizontal angle Hip, vertical angle Vip and oblique distance Dip to decide the power expression formula as follows:
Wherein: m
_{β}Be the nominal accuracy of instrument angle measurement, unit is second; m
_{D}Be distance accuracy, unit is millimeter; A, b are respectively fixed error and the proportional error coefficient of total powerstation nominal;
By least square method V
^{T}PV=is minimum, forms equation:
$\left\{\begin{array}{c}\mathrm{NX}+W=0\\ X={N}^{1}W\end{array}\right.,$ Wherein
$\left\{\begin{array}{c}N={A}^{T}\mathrm{PA}\\ W={A}^{T}\mathrm{PL}\end{array}\right.,$ Thus, adjustment is resolved and can be obtained the survey station unknown number
With monitoring point coordinate (X
_{Ck}, Y
_{Ck}, Z
_{Ck}) exact value;
(5) accuracy assessment
If n is the error equation number, t is the unknown number number, m
_{i}It is the precision valuation of i parameter; Q
_{Ii}Be the capable i column data of i on the inverse of weight matrix Q diagonal line; By adjustment result residual error V, but error in the unit of account power
Wherein
Unknown number inverse of weight matrix Q=(A
^{T}PA)
^{1}, parameters precision is estimated
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN 201210051569 CN102607506B (en)  20120301  20120301  Free stationing transformation monitoring method of highfill airport side slope unit set total station 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN 201210051569 CN102607506B (en)  20120301  20120301  Free stationing transformation monitoring method of highfill airport side slope unit set total station 
Publications (2)
Publication Number  Publication Date 

CN102607506A CN102607506A (en)  20120725 
CN102607506B true CN102607506B (en)  20131030 
Family
ID=46525111
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN 201210051569 Expired  Fee Related CN102607506B (en)  20120301  20120301  Free stationing transformation monitoring method of highfill airport side slope unit set total station 
Country Status (1)
Country  Link 

CN (1)  CN102607506B (en) 
Families Citing this family (14)
Publication number  Priority date  Publication date  Assignee  Title 

CN103471519B (en) *  20130927  20160907  国家电网公司  A kind of method applied without the deformation of prism photoelectricity total station survey power transmission and transformation pylon 
CN104729489B (en) *  20150420  20170301  中铁二局工程有限公司  One kind how public turning point lateral intersection traverse survey method 
CN105115473B (en) *  20150914  20170926  四川大学  The measuring method of underground chamber cross section 3 D deformation 
CN106599483B (en) *  20161216  20190827  鞍钢集团矿业有限公司  Processing method based on robot measurement openpit slope monitoring panel data 
CN108195341B (en) *  20171120  20201127  同济大学  Automatic deformation monitoring method for instrument station transfer without need of relearning 
CN108981659A (en) *  20180727  20181211  中国十七冶集团有限公司  A kind of measurement method of fullautomatic GPS displacement observation equipment 
CN110044342B (en) *  20190419  20210115  武汉地震计量检定与测量工程研究院有限公司  Threedimensional intersection measuring method 
CN110788857B (en) *  20191018  20210525  北京建筑大学  Automatic point searching and aiming method of measuring robot and measuring robot 
CN110806194A (en) *  20191204  20200218  中国十七冶集团有限公司  Steel truss deformation monitoring and measuring method 
CN110906902A (en) *  20191205  20200324  国电大渡河流域水电开发有限公司  Twostep correction method for slope deformation monitoring data 
CN111473737A (en) *  20200413  20200731  中国三峡建设管理有限公司  Method for measuring and setting precision control network for freely setting underground engineering of hydropower station 
CN112504214B (en) *  20201216  20210824  北京建工集团有限责任公司  Monitoring device and monitoring method for deformation point of steel member 
CN112906198B (en) *  20210121  20211022  西南交通大学  Method for processing nonlinear least square adjustment data of arbitrary plane measurement control network 
CN114001704B (en) *  20211229  20220329  四川中水成勘院测绘工程有限责任公司  Engineering construction deformation monitoring method based on free station setting reduction algorithm 
Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

JP3028198B2 (en) *  19960328  20000404  昭和アルミニウム株式会社  Fixture cover device 
CN101629799A (en) *  20090818  20100120  中国科学院武汉岩土力学研究所  Nonintervisibility high and steep side slope deformation monitoring method and device thereof 
CN102002955A (en) *  20101109  20110406  河海大学  Slope deformation monitoring and arranging method 
CN102252646A (en) *  20110415  20111123  中国水利水电科学研究院  Dam and side slope threedimensional continuous deformation monitoring system 

2012
 20120301 CN CN 201210051569 patent/CN102607506B/en not_active Expired  Fee Related
Patent Citations (4)
Publication number  Priority date  Publication date  Assignee  Title 

JP3028198B2 (en) *  19960328  20000404  昭和アルミニウム株式会社  Fixture cover device 
CN101629799A (en) *  20090818  20100120  中国科学院武汉岩土力学研究所  Nonintervisibility high and steep side slope deformation monitoring method and device thereof 
CN102002955A (en) *  20101109  20110406  河海大学  Slope deformation monitoring and arranging method 
CN102252646A (en) *  20110415  20111123  中国水利水电科学研究院  Dam and side slope threedimensional continuous deformation monitoring system 
NonPatent Citations (3)
Title 

JP特许第3028198号B2 2000.02.04 
全站仪坐标差分法在高边坡变形监测中的应用;陈子进等;《重庆建筑大学学报》;20050630;第27卷(第3期);第130134页 * 
陈子进等.全站仪坐标差分法在高边坡变形监测中的应用.《重庆建筑大学学报》.2005,第27卷(第3期),第130134页. 
Also Published As
Publication number  Publication date 

CN102607506A (en)  20120725 
Similar Documents
Publication  Publication Date  Title 

CN102607506B (en)  Free stationing transformation monitoring method of highfill airport side slope unit set total station  
CN102998690B (en)  Attitude angle direct resolving method based on global position system (GPS) carrier wave doubledifference equation  
CN103499340A (en)  Measurement device and measurement method for vertical greatheight difference height transmission  
CN101493324A (en)  Region quasigeoid refining method based on CQG2000  
CN106959456A (en)  A kind of GNSS SURVEYING CONTROL NETWORKs Accuracy Estimation  
CN104864858A (en)  Laying method and device of engineering control network  
CN103983219A (en)  Insuit measuring method of large size planeness  
CN104808225A (en)  Measurement method, correction method and measurement device of singlepoint satellite positioning errors  
CN108151698A (en)  A kind of antenna rotation center Calibration Method for intersecting method based on axis  
Muller et al.  Integrated velocity field from ground and satellite geodetic techniques: application to Arenal volcano  
CN106597504A (en)  Measurement system and method for building construction  
CN110261876A (en)  The unrelated GNSS of high precision position monitors virtual reference method  
CN108168504B (en)  Antenna rotation center calibration method based on fitting sphere center method  
Shouny et al.  Evaluating the performance of using PPKGPS technique in producing topographic contour map  
CN104330078B (en)  Combined measuring method based on threepoint resection model  
CN105388494A (en)  Laser ranging positioning method for RTK receiver  
Chrzanowski  Modern surveying techniques for mining and civil engineering  
CN102121825A (en)  Method for measuring rivercrossing levelling altitude difference by using GPS (global positioning system)  
Badejo et al.  Replacing orthometric heights with ellipsoidal heights in engineering surveys  
Bagherbandi et al.  Physical and Geometric Effects on the Classical Geodetic Observations in SmallScale Control Networks  
CN103884319A (en)  Total stationfree tower footing crosssection diagram measuring method and equipment  
Shrivastava et al.  Geoid Models for Indian Territory  
Nikonov et al.  Geodetic base for the construction of groundbased facilities in a tropical climate  
GarcíaAsenjo et al.  Establishment of a multipurpose 3D geodetic reference frame for deformation monitoring in Cortes de Pallás (Spain)  
Khalil  Alternative Solutions for RTKGPS Applications in Building and Road Constructions 
Legal Events
Date  Code  Title  Description 

C06  Publication  
PB01  Publication  
C10  Entry into substantive examination  
SE01  Entry into force of request for substantive examination  
C14  Grant of patent or utility model  
GR01  Patent grant  
C17  Cessation of patent right  
CF01  Termination of patent right due to nonpayment of annual fee 
Granted publication date: 20131030 Termination date: 20140301 