基于gps的海湾大桥水平井增注线拟合研究.docxVIP

基于gps的海湾大桥水平井增注线拟合研究 青岛海湾大桥位于青岛主城区300国道杨家社入口处，横跨胶州湾海域，西至黄岛红石崖。全长约35.4公里，其中海上段周长26.75公里。在大桥区域内利用VRS技术建设了连续运行GPS参考站,并设有带一等水准的首级GPS控制网,红岛、青岛、黄岛共分布15个测点,如图1所示,为大桥建设提供高精度测量基准。 当前,水准测量仍然是建立高程控制基准的主要方法,但在胶州湾海域利用传统水准测量传递水准十分困难,且费用昂贵,所以考虑利用GPS水准方法。GPS水准方法的关键在于似大地水准面的拟合。目前,国内外拟合似大地水准主要采用的几何方法有:绘制等值线图、解析内插、曲面拟合。本文针对青岛海湾大桥测点的分布情况,利用曲面拟合似大地水准面,研究青岛海湾大桥GPS水准代替局部高精度水准的问题,并推广到其他类似桥梁建设和工程测量中。 1 似效果测量辅助 随着GPS定位技术的应用和发展,平面控制测量能获得相当高的精度,而高程精度还有待提高。GPS直接测量的高程与工程建设所需高程,二者坐标系统不一致,转换时会受到区域性似大地水准面精度等因素影响,于是越来越多的人开始关注综合 利用GPS测量和水准测量资料提高似大地水准面精度的方法。 GPS高程为相对于特定参考椭球面的大地高H,测区联测点的水准为以似大地水准面为基准的正常高h,似大地水准面上点量测到参考椭球面的距离被称为高程异常δ。利用测区联测点先算出测点δ,再与测点坐标建立函数模型,模拟该地区似大地水准面,然后内插即可算出测区内任一点δ,最后利用式(1)进行高程转换,将大地高转换为正常高。 h=H-δ(1) 青岛海湾大桥整个测区呈三角片状分布,高程异常无剧烈变化,似大地水准面变化平缓,适宜曲面拟合。 2 移动拟合物结构 目前普遍采用的曲面拟合法主要有最小二乘法拟合、移动最小二乘法拟合、多曲面函数拟合。其中多曲面函数拟合不适合青岛海湾大桥这样地势平坦且测点偏少的地区,所以本文主要采用前两种方法。 2.1 调查区参数模型分析 传统最小二乘法拟合一般取多项式: p(xI,yI)=a0+a1xI+a2yI+a3xIyI+ a4x2IΙ2+a5y2IΙ2+… (2) 主要有三参数(前三项)、四参数(前四项)、六参数(前六项)以及更高次方多项式模型。一般高次多项式在数据点较少的情况下,往往会引起误差,且内插结果没很大改进。此测区首级控制网联测点偏少,利用统计检验法针对本测区具体情况将前三种参数模型比较,找出最佳多项式。 以六参数与四参数比较为例分析统计检验法。表1中:Ω1、Ω2分别表示六参数和四参数残差平方和;R表示残差改正项,且R1=Ω2-Ω1,R与Ω相互独立,可构造F分布的统计量。 F=R/2Ω/(n?6)F=R/2Ω/(n-6)(3) 式(3)中:分子自由度为六参数比四参数多出参数; 分母自由度是六参数平差函数多余观测值,本例中n=15。对F分布做右尾检验,给定显著水平α=0.05,用α和分子、分母自由度查F分布表,得右尾分位值F′。表1中方案1为六参数与四参数比较的各项指标,F′F表示六参数比四参数多引入的参数对拟合起显著作用。四参数与三参数的比较见表1中方案2 ,可见F′F表示四参数比三参数多引入的参数对拟合效果无显著提高,所以四参数与三参数拟合效果相当。六参数与三参数的比较见表1中方案3。 综合分析本测区,选择二阶六参数多项式。 2.2 不同约束域权函数的拟合精度 相对于传统最小二乘法,它引入紧支概念,认为x处节点值只受x附近影响域内节点影响,并在影响域上定义权函数。对于移动最小二乘法拟合GPS水准,根据前面对最小二乘法多项式的讨论,本测区采用二阶六参数多项式作为基函数,所以关键就是权函数的确定。 不同权函数对计算结果有较大影响,其选择有以下规则:测区中选择某点为中心点,半径为rI的圆形影响域,权函数在此域外为0;拟合点到中心点距离为dI,权函数随dI增加而单调递减;权函数在域内满足常态性且随rI趋近0,其趋近于Diracdelta函数。权函数一般使用样条函数,基本通式为式(4)。 其中:t—权函数阶数,S=dIrI,akS=dΙrΙ,ak—待求系数。利用边界条件计算出t=5,6,7时权函数表达式。其中:t=m0+m1+1(m00,m10),m0、m1为S分别等于0、1两种情况下权函数为1的偏导阶数。 本测区只有15个带有GPS数据的水准点,分成几个影响域时拟合点数太少,所以取整个区域作为一个影响域。任取除2,8,14外12个点作为拟合点,10号作为中心点,其他点到此点最大距离为半径r。从表2中看出,本测区在t=5,m0=3,m1=1条件下的权函数使移动最小二乘法拟合精度最好。 分别利用最小二乘法和移动最小二乘法拟合全部测点,精度