贯通测量误差预计新公式.docVIP

导线测量对贯通影响的误差预计新公式 杜宁　王莉　张显云 摘要　鉴于大型贯通测量中对误差预计精度的要求以及计算机在误差预计中的广泛应用，本文给出了一组适用于导线测量对贯通误差预计的电算化严密公式。关键词　导线测量　贯通　误差预计 1　前　言　　目前，在进行贯通误差预计时，导线在贯通面水平重要方向，即垂直于巷道中线方向——横向（X方向）上的误差预计公式为： 式中　　mβ——导线的测角中误差，以秒为单位；　　ΣR2y——各导线点至贯通面的垂直距离的平方和；　　ρ——取206265＂；　　——导线边的相对中误差；　　Σdx2——各导线边在贯通面上投影长度平方的总和。　　而上述公式存在以下几个问题：　　a　手工从贯通测量设计图上量取ΣR2y、Σdx2，不方便，速度慢，且受作图与量取误差的限制；　　b　不利于用计算机进行贯通误差预计以及贯通相遇点最佳位置的选取。　　为此，笔者推导以下电算化误差预计严密公式。2　公式推导2.1　导线在贯通点K处的误差计算式　　如图1所示，K为贯通点，X方向为贯通面水平重要方向，巷道（或隧道）在未贯通前，贯通面两端的导线为支导线（导线Ⅰ、导线Ⅱ）。 图　1 　　a　支导线I在贯通点K处的方差及协方差计算式　　把K点看作为支导线I的终点，则有： 　　(1) 　　而导线任意边i的方位角是所测角度的函数，即 ? (2) 上列两式中　　XB——支导线I起算点B的已知X坐标；　　YB——支导线I起算点B的已知Y坐标；　　αAB——支导线I起算边AB的已知坐标方位角；　　βj——支导线I各转折角；　　　　li——支导线I各边的边长，如ln1表示导线点n1至贯通点K的边长；　　αi——支导线I各边的坐标方位角，如αn1表示导线点n1至贯通点K的坐标方位角；　　n1——支导线I的总点数，起算点B编为1号，贯通点K前最后一个导线点编为n1号。　　考虑到式（2），对式（1）进行全微分，得： 　　(3) 　　运用协方差传播定律，并顾及到观测量之间相互独立，方差阵为如下一对角阵： 得： 　　(4) 　　b　同理，得支导线Ⅱ在贯通点K处的方差及协方差计算式 　　(5) 　　c　导线在贯通点K处的方差及协方差计算式　　因，由支导线Ⅰ、Ⅱ分别推算到K点时，在X与Y方向上的偏差为： 　　将上式写成矩阵形式，为：　　F=AX式中 而，　　所以，运用协方差传播定律，得支导线Ⅰ、Ⅱ对K点综合影响的方差及协方差为： 　　(6) 上式中　　n——导线Ⅰ、Ⅱ的导线点数之总和，即n=n1+n2。　　d　导线在贯通点K处的误差计算式　　以σ2xK、σ2yK、σxyK、σ2l1、…、σ2ln、σ2β1、…、σ2βn的估值m2xK、m2yK、mxyK、m2l1、…、m2ln、m2β1、…、m2βn代入式（6），得： 　　(7) 式中　　n——导线的总点数，n=n1+n2；　　αi——导线各边的坐标方位角；　　xK、yK——贯通点K的X、Y坐标；　　xi、yi——导线各点的X、Y坐标；　　mli——导线各边的测边中误差；　　mβi——导线各转折角的测角中误差；　　ρ——取206265＂。如测边采用相对中误差，并顾及m2β1=…=m2βn=m2β，则由上式可得： 　　(8) 式中　　dxi——导线各边在X方向上的坐标增量；　　dyi——导线各边在Y方向上的坐标增量；　　——导线边的边长相对中误差；　　其它同式（7）。2.2　导线在贯通面横向（X方向）上的位差　　如图2所示，K点为贯通点，X—O—Y坐标系为测量坐标系，X—K—Y坐标系为K点处假定坐标系，X轴垂直于巷道中线方向，Y轴为巷道中线方向。 图　2 　　因在给定方位角φ方向上的位差计算式为：　　m2φ=m2o(Qxxcos2φ+Qyysin2φ+Qxysin2φ）　　所以，由导线引起的K点在X方向上的位差计算式如下： 　　(9) 式中　　γ——X轴在测量坐标系X—O—Y中的坐标方位角，γ=α中线±90°。2.3　误差预计公式　　取2倍中误差作为极限误差。　　a　导线测量一次时，误差预计公式为： 　　(10) 　　b　导线复测n次时，误差预计公式为： 　　(11) 3　实　例　　图3所示为某矿一已贯通巷道，贯通导线总长约5.0km（其中地面布设为长约1.7km的二级导线，井下为长约3.3km的15＂基本控制导线），平均边长200m。水平角观测采用T2经纬仪，测距采用REDmini（5+5ppm）型测距仪，导线进行复测。取mβ地面=±10＂、mβ井下=±15＂、 图　3 　　a　原设计采用老公式进行误差预计，得MX预=±0.130m。　　b　现采用本文公式，运用Quick BASIC语言编程计算，得　　MX预