金厂沟梁金矿1#与35#竖井七中段的贯通测量研究

内蒙古科技大学本科生毕业设计说明书题目：金厂沟梁金矿1#与35#竖井七中段的贯通测量研究学生姓名：学号：专业：测绘工程班级：09-1指导教师：金厂沟梁金矿1#与35#竖井七中段的贯通测量研究摘要贯通测量是目前矿山测量中重要的工作之一，贯通测量的精度要求相对于井下其它测量应该相对精确，对于主要运输巷道的贯通，它的意义不仅在于技术上的精度要求，更重要是其直接关系到矿山的生产效益。本论文主要完成了金厂沟梁矿区两井间主要运输巷道贯通的方案设计。首先通过对矿区地面控制网的检核，进行了1#与35#竖井的近井点重测与检核、水准点的往返测量、矿井的联系测量等工作，然后布设井下经纬仪导线测量开切点坐标位置及方位角，分析其在水平方向和竖直方向上对贯通点K的影响，并进行误差预计的计算工作，选择合理的设计方案。最后根据合理的贯通方案对待贯通巷道进行施测，其主要测量工作内容包括：地面控制测量，联系测量，井下导线测量等。关键词：贯通测量，误差预计，平面测量，高程测量AbstractThroughmeasurementisoneoftheimportantworkoftheminesurveying，Throughthemeasurementaccuracyrequirementsrelativetootherundergroundmeasurementshouldberelativelyaccurate，Forthroughthemainhaulageroadway,itssignificanceliesnotonlyinthetechnicalrequirementsofprecision,moreimportantisthatitisdirectlyrelatedtothemineproductionefficiency.ThispapermainlydescribesthedesignschemeofJinchanggouliangMineWelltwomaintransportroadway,first,Checkingthegroundcontrolnetworkofminingarea,thenearwellpoint1#and35#shaftweightmeasurementandverification,benchmarkreturnmeasurement,contactmeasurementworkofmine,ThenthelayoutofUndergroundTheodolitetraverseatangentialpointpositioncoordinatesandazimuthangle,andanalyzesitsinfluenceonthevalueofKinthehorizontaldirectionandtheverticaldirection,andthecalculationerrorisexpected,areasonablechoiceofdesignplan.Finally,accordingtotheschemeweremeasuredthroughpropertreatmentthroughthetunnel,themainworkincludes:themeasurementofgroundcontrolsurvey,contactmeasurement,undergroundtraversesurveyetc..Keywords:throughsurvey,errorestimates,;levelmeasurementheightmeasurement目录TOC\o"1-3"\h\u12362摘要 ②用水准测量或经纬仪三角高程测量连测两端巷道中的已知高程控制点，求出高程闭合差，它也实际上反映了贯通高程测量的精度。5.4.2贯通后巷道中腰线的调整(1)中线的调整巷道贯通后，实际的中线偏差如果在规定的容许范围之内，对于次要的巷道我们只需要将最后几架棚子加以修改即可[17]。但是由于本次贯通属于主要的运输巷道，我们必须将与相遇点一定距离的两端中心线A和B相连接，用新标定的中心线代替原来两端的中线和如图5.8所示图5.8运输巷道贯通后中线的调整(2)腰线的调整如果贯通偏差在实际中很小时，我们可以依据距离和实测高差计算出最后一段巷道的坡度，重新标定出新的腰线。在平巷中，如果高差在贯通时较大，可相应的延长调整巷道的的距离[18]。如图5.9所示。实际测得贯通高程的偏差为60mm，我们可以从贯通相遇点向两端各后退30m，仍然和该处的原有的腰线点相连接。就可以得到被重新调整后的腰线，其坡度就从原设计的4%变为了3%，这样就可以有效的减缓巷道的坡度，使主要贯通巷道适合运输矿石。图5.9运输巷道贯通后腰线的调整第六章实测精度的分析在贯通测量中，尤其是针对本矿山这种主要运输巷道的大型贯通测量工程，都必须进行两次、三次甚至更多次的测量，这样的积累测量才能得到大量的实测资料，使我们根据这些实测资料进行精度的分析，以便评定实测成果的精度是否满足贯通的要求。虽然根据两次、三次的成果估算出来的定向和导线中误差数据很少，估算的方法也不是非常的严格，然而，我们这样独立测量的实际偏差也在一定程度上反映了实测成果的质量问题，以此来作为我们衡量实测误差的参考标准[19]。6.1终边方位角和终点坐标井下测量导线采用全站仪，1#竖井起始点至35#竖井A点导线全长887.001m，共10个测站，要求测角中误差，导线相对精度为1/12000，导线均独立进行测量了三次，其最终边的方位角和点位中误差如表6.1所示表6.1井下终边方位角和终点坐标值次数123平均最终边方位角+20-12-8最终点坐标X（m）672.419672.380672.383672.394（m）+0.025-0.014-0.011Y（m）355.862355.890355.864355.872（m）-0.010+0.018-0.008（m）0.0270.0230.014导线相对闭合差1/328521/385651/63357导线的实际测角中误差式（6.1）导线终点位置中误差式（6.2）井上两井间的导线连测采用全站仪测量，从1#竖井近井点至35#竖井近井点连测，导线全长967.545m，共7测站，要求测角的中误差为，导线相对精度为1/12000，导线均独立进行测量了三次，其最终边的方位角和点位中误差如表6.2所示。表6.2井上导线终边方位角和终点坐标次数123平均最终边方位角+9-3-6最终点坐标X（m）671.058671.032671.036671.042（m）+0.016-0.010-0.006Y（m）330.027330.047330.031330.035（m）-0.008+0.012-0.004（m）0.0180.0160.007导线相对闭合差1/537521/604721/138221导线的实际测角中误差式（6.3）导线终点位置中误差式（6.4）6.2井下起始边方位角的中误差本次工程属于一井独立定向工作，两个矿井的定向测量独立完成，因此应该分别对57#竖井和15#竖井的井下起始边的方位角中误差分别计算。1#竖井定向对起始边方位角中误差的影响:表6.3井下起始边方位角中误差计算次数起始边方向计算1+214412-152253-636平均702由定向测量误差引起的井下起始边方位角的中误差式（6.5）35#竖井定向对起始边方位角中误差的影响表6.4井下起始边方位角中误差计算次数起始边方向计算1+183242-101003-864平均488式（6.6）地面导线，定向及井下导线引起的总的点位误差式（6.7）通过井下导线的三次重复测量，起始边和其实方向与原测量资料数据对比，起始边的方向闭合差为，坐标闭合差，，，通过此次结果显示，在贯通方案中误差预计的工作很准确，但是，同时也反映出，无论是误差预计工作，还是实测精度分析的工作，其只能对贯通巷道的误差做一个范围内的估算，并不能分析出确切的数值。第七章总结与展望在矿山测量中贯通测量占据了主导的地位，而贯通测量中，贯通的误差预计的工作是必不可少的。只有在贯通方案中，将误差预计工作做好，才能保证矿山巷道的正确贯通。此次设计论文，我主要是针对了金厂沟梁金矿的1#与35#竖井七中段的主要运输巷道的贯通做了详细的误差预计工作。本次矿山的贯通目的在于将矿井的提升量最大化，而此次贯通工程又属于主要的运输巷道，贯通的距离也比较远，因此我们将井上的已知控制点及井下的起始边方位角都做了精确的复查工作，确保此次贯通的顺利完成。此次贯通的设计的方案及误差预计的工作都非常的准确，在确定了原有的已知数据正确的情况下，我们依据矿山的测量规范，对井下的巷道进行了经纬仪导线的布设，并在施工阶段，对贯通巷道的中腰线标定也完全按着规范要求进行测量，因此，在贯通结束后，我们队实际施测的精度进行了评定工作，也从客观上反映了此次贯通工程误差预计工作的正确性，完成了矿山的生产建设需求。矿山贯通测量工程在我国的发展历程对比国外是非常缓慢的，因此，我们这一代的测量知识分子，应该对此进行深入的研究与开拓。我国是一个发展中的大国，对于矿山资源的需求也是非常迫切的，我们只有做好矿山测量的工作，才能确保矿山生产的顺利进行，才能加快矿山的发展速度。参考文献[1]汪云甲;郭达志;邓喀中;卞正富;吴立新;张书毕;杜培军.我国矿山测量学科的发展与创新[J],测绘通报,2005，2：1-6[2]王恕，隋振义.金厂沟梁金矿志[M]，内蒙古民族出版社，1998.11,1-30[3]付新启.测量学[M]，北京理工大学出版社，2006[4]罗志清.测量学[M]，云南大学出版社，2010[5]郭玉社.矿井测量与矿图[M]，北京：化学工业出版社，2007.7，123[6]张国良.矿山测量学[M]，江苏：中国矿业大学出版社，2006.[7]朱红侠.矿山测量[M]，重庆大学，2010.1[8]李天文.GPS原理及应用[M]，北京：科学出版社，2003[9]崔希民.测量学教程[M]，北京：煤炭工业出版社，2009[10]姚德言.金属矿山两中段一次定向实践[J]，采矿技术.2012，5：88-90[11]苏玉明.矿山测量中高程导入方法的研究与探讨[J]，科技风，2012，7：21-21[12]崔学敏[1]赵雪松[1]姚继华[1]王永[2].巷道中线标定的方法及实践[J]，科技信息，2010，2：69-70[13]王观宇.贯通巷道开切点位置计算公式的推导及其应用[J]，矿山测量，1992，4：10-22[14]周小安.工程测量[M]，西南交通大学出版社，2007[15]唐其武，冯明伟.煤矿专业基础知识读本[M]，重庆大学出版社，2012,27-28[16]焦亨余.工程测量[M]，重庆大学出版社,2010[17]王丽叶.贯通后实际偏差的测定和中腰线的调整[J]，科技风，2012，1：15[18]张国良.矿山测量学[M]，中国矿业大学出版社发行科,2006.8.1[19]陈步尚，陈国山.矿山测量技术[M]，北京：冶金工业出版社，2009.8[20]郝福恒.生产矿井测量[M]，内蒙古：内蒙古大学大学出版社，2008附录A井上下导线误差预计图附录B1#与35#竖井巷道贯通图附录CAnewmethodforminingdeformationmonitoringwithGPS-RTKGAOJing-xiang1,2,LIUChao2,WANGJian1,2,LIZeng-ke2,MENGXiang-chao21.KeyLaboratoryforLandEnvironmentandDisasterMonitoringofStateBureauofSurveyingandMapping,ChinaUniversityofMiningandTechnology,Xuzhou221116,China;2.SchoolofEnvironmentandSpatialInformatics,ChinaUniversityofMiningandTechnology,Xuzhou221116,ChinaReceived19June2011;accepted10November2011Abstract:Basedonrangingintersectiontheory,anewmethodwhichissimpleandeasytooperatewasproposedfordatacollectionintheminesurfacedeformationmonitoringwithGPS-RTKcenteringrodmeasurements.Itcanfullyeliminatetheinevitableshakingerrorandtheverticaldeflection,andtosomeextentweakenthemultipatheffectontheestimatesofcoordinatesinarelativelyshortperiodoftime,usinghigh-frequencyobservations.Theresultsshowthatthree-dimensionalcoordinateswithaheightaccuracybetterthan1cm,horizontalaccuracybetterthan2−4cmcanbeachievedthroughonly15−30scontinuousobservationby20Hzhigh-frequencyandeffectivelyimprovethemeasurementaccuracyandefficiencyofRTK,fullysatisfyingthehigh-speedandhigh-precisiondataacquisitioninminesurfacesubsidencedeformationmonitoring.Keywords:GPS;roverpole;multipathdelay;deformationmonitoringInordertodevelop,utilizeandprotectthecoalresourcesrationally,andpromoteandensurethedevelopmentofthecoalindustry,theChinesegovernmenthasenactedandimprovedsomerulesandregulations.Akeyproblemoftherulesisthathowtoconfirmtherangeofprotectingcoalcolumnbydifferentcoalminingmethods,burialdepth,seamthickness,roofandfloorlithologiccharacters,geologicstructure[1].Typically,aseriesofobservationstationsarrangedatregularintervals(15−25m),measuringrelativedisplacementsofthestationsbythetraditionaltechniques(leveling,totalstation,etc),wereestablishedtosummarizethechangelawofparametersofsurfacemovementforproductionandconstructionofmine[2].However,thesurfacedeformationisexpensivetocarryoutorcannotbemeasuredbytraditionaltechniquesinsomecoalmineareas,especiallyinthemountain.Hence,itisnecessarytoemployeffectivetechniquestoobtainthedeformationdata.Recentyears,interferometricsyntheticapertureradar(InSAR)[3−4],three-dimensional(3D)laserscanner[5−6],digitalclose-rangephotogrammetry[7−8]havereceivedincreasinglygreatattentionsandhaveachievedfirstfruits.However,atpresenttheyhaveboundednesstobeusedfordeformationmonitoringofminesurfaceandrelatedtheoriesshouldbecontinuouslyperfectedinfutureresearch.RecentadvantagesinGPSreceivers,antennasanddataprocessingsoftwares,havemadeGPSRTKtechniqueaveryeffectivetoolfordeformationmonitoringincoalminingareas[9],butthedata-acquiredqualitystillneedstobeimproved.Withthemainerrorsourcesincludingroverpoledeflectionofthevertical,un-modeledsystematicerror(grosserror,multipatherror,etc)andtheheighttransformationerror,planeaccuracycanonlyreach5−10cm,theheightaccuracycanonly3−5cm,whichrestrictthefurtherapplicationofthistechnology.Usingtripodinsteadofroverpoleandwithlongerobservationtotakefulladvantageofthecyclicalnatureofmultipathcaneffectivelyeliminateorreducethesystemerror,butforalargenumberofmininggrounddeformationobservationstations,themethodmeansahugeamountoftasks.Inthepresentwork,asimplehigh-precisionGPS-RTKmeasurementwasproposedwiththeroverpole.Thatis,thepositionofthestationcanbeobtainedbyrangingintersectionwiththecoordinateseriescollectedbycontinuouslyshakingGPSantenna.Therelevanttheorywasdeduced,theaccuracyfactorindifferentdirectionswasgiven,andthefeasibilityofincreasingsamplingfrequencyandreducingtheobservationtimetoweakentheeffectofmultipathwasanalyzedinthismodeTwoexperimentsweredonetoevaluatetheprecisionandaccuracyofthemodels.2TheoreticalfoundationUsually,roverpoleisusedinsteadoftripodtoimproveoperationalefficiency,butthedeflectionandtheshakingerrorswillappearinevitably.Whilethelengthoftheroverpoleis2mandthedeflectionangleoftheverticalis3°,thethehorizontalandtheverticaldeflectionscanreach104.7and2.8mm,respectively,whicharemainrestrictivefactorsofGPSRTKtechniqueapplicabletominingsubsidencemonitoringandshouldbecorrected.Withoutconsiderationofobservationerrors,thecoordinateseriescollectedbycontinuouslyshakingGPSantennawereonasphereasdiscretepoints(Fig.1).Therefore,thepositionsofthestationscanbeobtainedbyrangingintersectionwiththelengthoftheroverpoleasthesidelength.Fig.1SchematicdiagramofrangingintersectionAtepochi,observationequationis(1)where[X(i)Y(i)Z(i)]iscoordinateoftheGPSantennaatepochi;[abc]iscoordinateoftheunkownpoint;ristherealdistancefromtheantanaphasecentertotheunkownpoint;δmultipath(i)isthemultipatherrors;ε(i)istheobservationnoiseatepochi.Let(xisthevectorrepresentationof[abc]).Thesecond-orderTaylorseriesexpansionfortheleftpartofequation(1)aboutx0(whichis[a0b0c0])is(2)where,whichisHessematrix,canbetreatedasaresidualfortheinsignificantimpactofthequadratictermsontheparameterestimationwithalongdistanceof2×104kmbetweenthesatelliteandthestationinaGPSsurveying.Inthiswork,itisnecessarytoevaluatetheinfluenceofthequadratictermsastheroverpoleisonly2mlong.TEUNISSEN[10]presentedaborderestimateofthesizeofthequadraticterms:(3)whereλminandλmaxaretheminimumandmaximumeigenvaluesoftheHessematrix,respectively,andthemaximumpossibleeigenvaluecouldbe1/r(risthelengthoftheroverpole),themaximumpossibleeigenvaluecouldbe0,thentheinfluencevaluesofthequadratictermsis[11](4)whereδxistheobservationerrorsofthereceiver.Figure2givestherelationshipbetweenthemaximumnonlinearerrorandthereceiverpositioningerrors.Withtheinfluenceofmultipatheffect,underthemostadverseconditions,thepositioningerrorisabout5cmandthemaximumnonlinearerrorisabout0.6mmforL1frequency.So,onlythefirst-orderTaylorseriesexpansionisadequateandthenonlinearerrorcanbeignored.Equation(1)canbeexpressedas[12]:Fig.2Estimationofthemaximumnonlinearerror(5)where[l(i)m(i)n(i)]isunitvectorcomponentinthethreeorthogonaldirections(N,E,U)ofthedistancebetweentheGPSantennaandstationpointatepochi;[δaδbδc]Tisthecorrectedvalueforthecoordinatesoftheunknownpoint;ρ0(i)istheapproximationofgeometricdistancebetweentheantennatothesatellite.N-errorequationsareobtainedlikeequationbelow,whentheobservationepochisnV=A·X−L（6）whereV=[V(1)V(2)···V(n)]T;；Accordingtotheprincipleoftheleastsquaremethod,theco-factormatrixofunknownparameteris（7）ResemblingthetheoryofGPSpointpositioning,introducetheDOPindifferentdirection:Afterthefieldexperiment,wecanfindthatthecenteringrodof2mcansucceedtowobbleintheanglebetween70°and90°withground.InordertosolvethechangeofDOPindifferentobservationtime,thecontinuouslyuniformdistributionfunctionisused.AsshowninFig.3,50groupsdataweregatheredandeverygrouphad300epochs.TheDOPdecreaseswhentheobservationtimeincreases.BecausetheswayrangeofGPSreceiverissmall,theDOPVisfarlessthantheDOPP.Whiletheepochis110,theDOPVdecreasesto0.1andtheDOPHdecreasesto0.8.Theresultindicatesthatthenewmethodcanweakentheinfluenceofthereceiverpositionerrortounknownpointcoordinateparameterestimationandimprovetheheightmeasurementprecisioneffectively.ThemainpositionerrorofGPSreceiverismultipath.MultipathisaphenomenonwherebyasignalisreflectedordiffractedfromnearbyobstaclesandFig.3DOPvaluesofdifferentepochsarrivesatareceiver’santennaviatwoormoredifferentpaths(Fig.4).Multipathdelayduetosinglereflectedsignalcomponentcanbedescribedas[13]:(8)where[14] (9)AccordingtoEq.(9),thephaseexcursioncausedbyreceivermultipathisdeterminedbytherelativepositionofreceiver,satelliteanddiffractingobjectinspace.TinyGPSreceiverchangewillbringtolargedegradationofmultipathrelativity[15].Inordertodecreasetheobservationtimeandimprovethemeasurementefficiency,thesamplingfrequencyshouldbeincreasedandtheswayrateshouldbeaccelerated. Fig.4MultipathphenomenonofGPSsignal3ExperimentsToexaminetheperformanceoftheprocedure,anextensivestudywasconductedonadatasetcollectedfromdifferentbaselinesinTheNationalScienceParkofCUMT(CNUSP)onOct.31st,2010(DOY304GPSweek1608)from2:00−4:45UTC(10:00−12:45localtime).Theexperiencewasperformedbythetripod,traditionalcenteringrodandthenewmethodofthispaperunderfourdifferentshelteredconditions.Thedatawereobtainedatthesamplingrateof1Hzand20Hz;themeasuretimeis180s.TheenvironmentoffourpointsisshowninFig.5.Figure6showsthecalculationheightintheobservationtimeof150,60,30and15sbydifferentFig.5Testsindifferentenvironments:(a)Point1;(b)Point2;(c)Point3;(d)Point4Fig.6Heightsbasedondifferentperiods(1Hz):(a)150s;(b)60s;(c)30s;(d)15smethods.Thesamplingrateis1Hz.Whiletheobservationtimeisdecreasing,thebiasofthethreemethodsincreasesinvariousdegreesandthenewmethodofthispaperisinfluencedlargely.Onthecontrary,thenewmethodgetsthemostexcellentresultthattheprecisionisbetterthan1cmintheobservationtimeof150sand60s.Figure7showsthedifferentpointFig.7Horizontalcoordinatesofdifferentmethods(150s,1Hz)Fig.8Heightsbasedondifferentperiods(20Hz):(a)15s;(b)30sFig.9Convergenceratesofcoordinates(No.1point-2ndperiod):(a)1Hz;(b)20Hzplanecoordinatesobservedbyvariousmethodsat1Hzandin150s.Thetripodmethodcangetthebestoutcomeandtheplanbiasofnewmethodisbetween2cmand3cm,whichisbetterthantraditionalcenteringrodmethod.Thedatasamplingfrequencyischangedindifferentobservationlengthtotesttheabovemethods.Figure8showsthecontrastiveheightresultofdifferentmethodsat20Hz.Theobservationtimeis15sand30s.Becausethetripodandtraditionalcenteringrodpositiondoesnotvaryandmultipathperiodisalmostfixed,thetwomethodsareabletoimprovetheprecisiontinily.Thenewmethodcangettheheightwithanaccuracyof1cm,whichimprovestheobservationprecision.Becausethenewmethodcanweakentheinfluenceofthereceiverpositionerrortounknownpointcoordinateparameterestimation,itcangetbetteroutcomethanthetripodmethodandtraditionalcenteringrod.Futheranalysisisstudyingthefeasibilityofnewmethodbyincreasingthesamplingfrequencyandreducingtheobservationtime.Thethree-dimensionalcoordinateofthefirstpointobservedindifferentsamplingfrequencyisshowninFig.9.Theobservetimeis150s.Theresultindicatesthatconstringencyrateat20Hzsamplingfrequencyisfasterthan1Hzsamplingfrequency,especiallyinN-directionandinE-direction.Figure10showstheplanecoordinatesofdifferentmethodsat20Hz.Theobservationtimeis15sand30s.Thenewmethodcangettheplanmetricpositionwithanaccuracybetween3cmand4cm,whichimprovestheobservationprecisionlargelythantraditionalcenteringrod.Fig.10Horizontalcoordinatesbasedondifferentperiods(20Hz):(a)15s;(b)30s4Conclusions1)Inordertoimprovetheaccuracyandreliabilityofminesurfacesubsidencemonitoring,GPSRTKtechniquewasintroducedandthemechanismofthemainsystematicerrorswasanalyzedsubsequently,basedondifferenttheoriesandmethodologies.2)Thenewmethodcanfullyeliminatetheinevitableshakingerrorandtheverticaldeflection,andtosomeextentweakenthemultipatheffect,whichimprovestheprecisionofestimates.3)Theexperienceprovesthatagoodresultcanbeobtainedafterincreasingthesamplingfrequencyanddecreasingtheobservationtime.4)Three-dimensionalcoordinateswithheightaccuracybetterthan1cm,horizontalaccuracybetterthan2−4cmcanbeachievedthroughonly15−30scontinuousobservationby20Hzhigh-frequencyandeffectivelyimprovethemeasurementaccuracyandefficiencyofRTKinminingareaobservationstationofsurfacemovement.5)Themethodshouldbestudiedfuturetoremovethemutationerrorofcenteringrodanddiffractionandlongperiodmultipath.References[1]TheProfessionalStandardsCompilationGroupofPeople’sRepublicofChina.Regulationsforsettingcoalpillarandminingpressedcoalunderbuildings,waterbodies,railways,andmainroadways[S].Beijing:ChinaCoalIndustryPublishingHouse,2000.(inChinese)[2]HEGuo-qing,YANGLun,LINGGeng-di.Miningsubsidencesciences[M].Xuzhou:ChinaUniversityofMiningandTechnologyPress,1991.(inChinese)[3]LUNDGRENP,USAIS,SANSOSTIE,LANARIR,TESAUROM,FORNAROG,BERARDINOP.ModelingsurfacedeformationobservedwithsyntheticapertureradarinterferometryatCampiFlegreicaldera[J].JournalofGeophysicalResearch,2001,106(9):355−366.[4]PERSIKIZ,JURAD.ERSSARinterferometryforlandsubsidencedetectionincoalminingareas[J].EarthObservationQuarterly,2001,63:25−29.[5]LOVELLJL,JUPPDLB,CULVENORDS,COOPSNC.Usingairborneandground-basedranginglidartomeasurecanopystructureinaustralianforests[J].CanadianJournalofRemoteSensing,2003,29(5):607−622.[6]ABELLANA,JABOYEDOFFM,OPPIKOFERT,VILAPLANAJM.Detectionofmillimetricdeformationusingaterrestriallaserscanner:experimentandapplicationtoarockfallevent[J].NaturalHazardsandEarthSystemSciences,2009,9:365−372.[7]RENWei-zhong,KOUXing-jian,LINGHao-mei.Applicationofdigitalcose-rangephotogrammetryindeformationmeasurementofmodeltest[J].ChineseJournalofRockMechanicsandEngineering,2004,23(3):436−440.(inChinese)[8]YANGHua-chao,DENGKa-zhong,GUOGuang-li.Monitoringtechniquefordeformationmeasurementofsimilarmaterialmodelwithdigitalclose-rangephotogrammetry[J].JournalofChinaCoalSociety,2006,31(3):292−295.(inChinese)[9]ZHAOXing-wang,YUXue-xiang.StudyonapplicationofGPS-RTKinsurfacemovement[J].CoalTechnology,2007,26(9):104−106.(inChinese)[10]TEUNISSENPJG.Nonlinearinversionofgeodeticandgeophysicaldata:diagnosingnonlinearity[C]//BRUNNERFK,RIZOSC.DevelopmentsinFour-DimensionalGeodesy.Berlin:Springer-Verlag,1987:241−264.[11]YANGGuang.IntegratedGPSandpseudolitepositioningfordeformationmonitoringsystem[D].Nanjing:HohaiUniversity,2004.(inChinese)[12]LIUJi-yu.GPSsatellitenavigationprincipleandmethod[M].Beijing:SciencePress,2008.(inChinese)[13]YUANLin-guo,HUANGDing-fa,DINGXiao-li,XIONGYong-lian,ZHONGPing,LICheng-gang.OntheinfluenceofsignalmultipatheffectsinGPScarrierphasesurveying[J].ActaGeodaeticaetCartographicaSinica,2004,33(3):210−215.(inChinese)[14]HUANGDing-fa,DINGXiao-li,CHENYong-qi,GONGTao,XIONGYong-liang.WaveletfiltersbasedseparationofGPSmultipatheffectsandengineeringstructurevibrations[J].ActaGeodaeticaatCartographicaSinica,2001,30(1):36−41.(inChinese)[15]DAIWu-jiao,DINGXiao-li,ZHUJian-jun.StudyonmultipatheffectinstructurehealthmonitoringusingGPS[J].JournalofGeodesyandGeodynamics,2008,28(1):65−71.(inChinese)对GPS-RTK矿区变形监测的一种新方法1.国家测绘局土地环境与灾害监测重点实验室，中国矿业大学、徐州221116,中国2.环境与测绘学院，中国矿业大学，徐州221116，中国收到2011年6月19日，发表于2011年11月10日；摘要基于测距交会理论，针对矿山地表变形监测时数据采集提出了一种利用GPS-RTK中心杆技术的简单易操作的新方法。它可以完全消除不可避免的晃动误差和垂直偏转，通过使用高频率的观察方法能在一定程度上削弱多路径效应对较短的时间坐标的估计。结果表明，通过仅仅15−30秒20赫兹的高频的连续观测有效提高RTK测量精度和效率，即可实现用垂直精度高于1厘米水平精度高于2−4厘米的三维坐标，充分满足高速高精度数据采集在煤矿地面沉降变形监测。关键词：GPS；车杆；多径延迟；变形监测；1引言为了合理发展，利用和保护煤炭资源，促进和保障煤炭工业的发展，中国政府制定和完善的一些规则和条例。该规则的一个关键问题是，如何通过确定不同的采煤方法、埋藏深度、煤层厚度、顶底板岩性、地质结构，保护煤柱的范围[1]。通常情况下，定期安排一系列观测站（15−25米），由传统的技术测量站的相对位移（水平，总台，等），总结地表移动参数生产和建设矿井变化规律[2]。然而，在一些煤矿，尤其是在山区，依靠传统的测量技术进行地表的变形测量。因此，有必要采用有效的技术获取变形的数据。近年来，合成孔径雷达干涉测量，三维（3D）激光扫描，数字近景摄影测量，越来越受到人们的重视并且已经取得了初步的成果。然而，目前它们在矿山地表及相关理论方面还具有一定的局限性，应该在今后的研究中进行完善。GPS接收机，天线和数据处理软件的新优势，使得GPS-RTK技术在煤矿的地表变形监测中成为非常有效的测量工具，但是，获得数据的质量方面还有待提高。主要误差来源包括垂直车杆挠度，非模型系统误差（人为误差，多路径误差，等）和高程转换误差，由于平面精度只能达到5−10厘米，高程精度只能3−5厘米，制约了该技术的进一步应用。使用三脚架来代替需要长时间观测的车级可以充分利用多径的周期性并有效地消除或减少系统误差，但对矿区地面大量的变形观测站而言，这种测量方法的任务量太大。在目前的工作中，提出了车杆这样一个简单的高精度GPS-RTK测量方法。即，就是说，站的位置可以通过GPS不断摇晃的天线，在坐标系下用测距交会的方法获得。经过相关理论的推导，给出了在不同方向上的精度的因素，并且通过这一模式还能得出提高采样频率和减少观测时间减弱多路径衰落的影响这一结论。这两个实验是进行评估模型的精度和准确度的。2理论基础通常，用车杆代替三脚架来提高操作的效率，但挠度和振动误差将不可避免。然而车杆的长度是2米，垂直偏转角度为3°，其水平和垂直偏差分别可达104.7毫米和2.8毫米，是GPSRTK技术适用于矿山地表沉陷监测的主要制约因素，应予以纠正。不考虑观测误差的情况下，通过连续振动的GPS天线收集到的一系列坐标作为一个球体上的离散点（图1）。因此，测站的位置可以通过具有一定长度的车杆作为边长，利用测距交会获得。图1测距交会示意图在时段i内，观测方程是：（1）[x（i）Y（i）Z（i）]是时段i的GPS天线坐标；[a,b,c]是未知点坐标；r是从天线相位中心到未知点的实际距离；δmultipath(i)是多路径误差；ε（i）是在时段i内观测噪声令（x是[a,b,c]的向量表示）。将方程（1）左边部分按（为[a0b0c0]）的二阶泰勒级数展开约为是（2）如果用表示黑塞矩阵，可以视为在距离长达2×104公里的卫星和空间站之间的GPS测量的二次项参数的估计影响不大的残差。在这项工作中，当车杆只有2米长时评价的二次项的影响9是必要的。伊尼森[10]提出了二次项大小的边界估计。（3）如果λmin和maxλ分别是的黑塞矩阵的最小和最大特征值，且最大特征值可能是1/R（R是罗孚杆的长度），最大可能特征值是0，那么二次项的影响值是[11]（4）δX是接收机的观测误差图2给出了最大非线性误差和接收机的定位误差之间的关系。由于多径效应的影响，在最恶劣的条件下，定位误差约为5厘米，最大非线性误差约为0.6毫米的L1频率。所以，仅展开一阶泰勒级数，非线性误差可以忽略不计。方程（1）可以表示为[12]图2最大非线性误差估计（5）[l（i）m（i）n（i）]是GPS天线和站点之间距离在时刻的i在三个正交方向（N，E，U）的单位矢量分量；[δaδbδc]T是未知点坐标的修正值；ρ0（i）是天线与卫星之间的近似几何距离。当观测历元是N，得到如下n阶误差方程V=A·X−L（6）当V=[V(1)V(2)···V(n)]T;；根据最小二乘原理原理，未知参数的有限因子矩阵是（7）类似于GPS单点定位理论，在不同的方向引入DOP：;;;通过试验，我们可以发现，2米的中心杆可以与地面形成70°至90°的夹角。为了解决DOP在不同观测时间的变化的问题，采用了连续均匀分布函数的方法。如图3所示，收集50组数据并且每一组有300个时段的DOP。观察时间的增加时DOP减少。因为GPS接收机的摆动幅度小，所以DOPV远小于DOPP。当时间段长为110时，DOPV降低到0.1，DOPH减小到0.8。结果表明，新方法可以减弱接收器的位置误差对未知点坐标参数估计的影响，有效提高高程测量精度。GPS接收机的主要位置误差多路径误差。多路径是指一个信号从附近的障碍反射或衍射到达接收器的天线通过两个或两个以上的不同的路径现象（图4）。图3不同时段的DOP值由于单个反射信号造成的多径时间延迟可描述为[13]：（8）在[14](9)根据式（9），接收机多径所引起的相位偏移由接收机相对位置、卫星和衍射空间物体决定。微型GPS接收机的变化将大大的降低多径相关性[15]。为了减少观测时间和提高测量效率，应增加采样频率和加速摆率。图4GPS信号的多路径现象3实验为了检查的程序的性能，2010年10月31日（doy304GPS周1608）从下午2:00−4:45UTC（10−当地时间12时45分），在中国矿业大学的国家科技园进行了一个收集不同基线。进行的实验是由三脚架，传统对中杆和本文的新方法，在四不同的庇护条件。获得的数据是通过1赫兹到20赫兹的采样率；测量时间为180s，四中环境如图5所示。图5在不同环境中的测试：（a）1点，（b）2点，（c）3点，（d）4点图6不同时段的高程（1赫兹）:(a)150s;(b)60s;(c)30s;(d)15s图6显示了用不同的方法观测150，60，30和15s的时间高程的计算结果。采样率为1Hz。当观测时间减少，第三种方法的偏差不同程度的增加，并且对本文叙述的新方法的影响很大。相反，如果让新的方法得到的最优结果观测精度必须高于1厘米，观测时间为150秒和60秒。图7显示了在1赫兹和150s的观测条件得到不同的点面坐标。三脚架的方法可以得到最好的结果和新方法的预计偏差在2厘米和3之间厘米，这是优于传统的对中杆的方法。图7不同方法横坐标（150s，1赫兹）通过测试以上方法，不同的改变观测时间的长度同样改变数据采样频率。图8显示为在20赫兹时不同方法的高程对比结果。观测时间为15秒和30秒。由于三脚架和传统的对中杆的位置不变化并且多径时间几乎是固定的，两个方法可以微微地提高精度。新的方法可以在高程上精确1cm，从而提高了观测精度。因为新的方法能够减弱接收机的位置误差对未知点坐标参数估计的影响，它可以比三脚架方法与传统的对中杆得到更好的结果。