我国陆地垂线偏差精化计算.docVIP

我国陆地垂线偏差精化计算 宁津生1、郭春喜1，2、王斌2、王惠民2 （1．武汉大学，武汉市珞喻路129号， 430079） （2. 国家测绘局大地测量数据处理中心，西安市友谊东路334号 710054） 摘要：我国天文大地网与高精度GPS2000网联合平差时，为满足地面观测数据归算到WGS84椭球的需要，我们采用先进的移去－恢复技术和高阶次地球重力位模型，应用全国重力资料和30″×30″数字地形模型，完成了参加联合平差的48919个大地点相应于WGS84椭球垂线偏差的精化计算，并用115个与天文点重合的GPS点的实测垂线偏差对计算结果进行了外部检验，得到全国垂线偏差子午分量的总体精度为±1.45″，卯酉分量的总体精度为±1.50″，该成果已提供联合平差使用。 关键词： 垂线偏差；WGS84；移去－恢复技术 1 引言 1.1 精化计算的目的 传统大地测量技术建立国家精密平面控制网（又称天文大地网）时，长度、水平方向等观测数据都是以垂线方向为依据完成的，必须对这些观测数据加以改正而归算到相应的参考椭球体上，转化为以参考椭球体法线为依据的归算数据，才能进行天文大地网的平差和大地位置的计算。在进行归算时，必须计算各大地点相对于参考椭球体的垂线偏差分量。 在实施我国天文大地网与高精度GPS2000网联合平差（以下简称联合平差）时，为提高垂线偏差的计算精度，满足天文大地网的观测数据归算到WGS84椭球的要求，需要求解各大地点相应于WGS84椭球的垂线偏差。因此，我们研究了垂线偏差的精化计算方法，并完成了48919个大地点WGS84椭球垂线偏差的计算。 1.2 垂线偏差实测值的应用 垂线偏差实测值由下式定义： （1） 为某天文基准系统的地面纬度、经度，B、L为相应椭球的大地纬度、经度。不管是参考椭球，还是平均椭球，都有其相应系统的实测值。 天文大地网是上世纪80年代以前施测的，其中的天文测量成果当时采用的是FK4星表系统，IAU1964天文常数系统，我国综合时刻1978年系统，天文大地网整体平差前重新联测得到的天文经度基准（经度由上海天文台徐家汇目视等高仪基座中心起算），地极坐标则是我国JYD（1968.0）系统。联合平差时应该采用FK5星表系统，IAU1976天文常数系统和IAU1980章动模型，我国综合时刻1978年系统（重新订正结果），全国天文基本网第二次联测得到的天文经度基准（经度由上海天文台佘山经度原点起算），以及IERS的IRP地极坐标系统。前者统称为旧天文系统，后者统称为新天文系统。 我国有实测天文点2千多个，经纬度测定精度平均为±0.3″,大部分属于旧天文系统。，并进行精度估计。 （5）应用联合平差中确定的天文大地网与GPS网的336个重合点，求出1980西安坐标系与WGS84的转换参数，将经过系统改正后的IAG75椭球垂线偏差值转换为WGS84椭球的垂线偏差。 （6）利用与天文点重合的GPS点的实测垂线偏差，进行外部精度检验。 2 垂线偏差计算数学模型 2.1 重力垂线偏差 采用移去－恢复技术计算： （2） 式中： ，G1项可采用地形改正值代替； ，剩余空间异常，空间重力异常，为用模型计算的重力异常，为采用加过G1项并移去模型重力异常后得到的地面剩余重力异常计算的垂线偏差，为地球重力场模型计算的垂线偏差，为地球平均重力，计算点与流动点间的距离，β1、β2分别为南北向东西向倾角。 2.2 数字积分方法计算 （2）式中： 计算点，流动点 编制计算重力垂线偏差的程序时，采用数字积分形式： （3） 3 计算实施 3.1 重力垂线偏差计算 重力垂线偏差计算时有两种选择，第一种是先完成格网点垂线偏差值的计算，然后内插出大地点上的值；第二种是直接计算全国任意大地点垂线偏差值。本次计算时，采用第二种方法直接完成参加两网联合平差的近48919个大地点垂线偏差值的计算。 在利用数字积分方法计算区域重力垂线偏差时，将整个积分范围定为以计算点为中心的3°×3°区域。在积分计算时，以计算点为中心的5′×5′区域采用30″×30″数据，其他区域采用2.5′×2.5′数据。 利用前面介绍的垂线偏差计算原则与方法，完成了参加联合平差的48919个大地点的重力垂线偏差的计算，并同时完成了作为系统控制的天文点重力垂线偏差的计算。 3.2 系统改正与精度估计 在计算重力垂线偏差时，中除包含数字高程模型、重力资料引起的误差外，还包含着数学模型中椭球系统不同和倾斜项β的影响，因此需要以相应椭球上的实测垂线偏差值作控制对重力垂线偏差进行系统改正以及精度估计。同时，鉴于垂线偏差与地形的变化趋