地球重力场、重力、重力梯度在三维直角坐标系中的表达式.pdfVIP

第29卷第5期 武汉大学学报·信息科学版 V01．29No．5 2004年5月 GeomaticsandInformationScienceofWuhan Mav2004· University 文章编号：1671—8860(2004)05—0377—03 文献标识码：A 地球重力场、重力、重力梯度在 三维直角坐标系中的表达式 陈俊勇1 (1 国家测绘局办公室，北京市百万庄，100830) 摘要：现代卫星重力测量主要利用星栽GPS接收机、加速度计、星栽测距仪等来确定重力卫星的轨道，削弱 非保守力的干扰，由此根据卫星的位置、速度及其变率来确定地球重力场。而上述GPS等星栽仪器所提供的 数据，包括卫星轨道坐标及其速率、扰动加速度、星间距离及其变率，都是以三维直角坐标(z．y，z)的形式表 示的．因此．地球重力场、重力和重力梯度在三维直角坐标系中的表达式在卫星重力解算中具有实际意义。 关键词：地球重力场；重力；重力梯度；三维直角坐标系 中图法分类号：P223．0；P223．6 近年来，卫星重力测量的实践证明了，它在求 定地球重力场，特别是中、长波长方面的重要作 1 地球重力场在球坐标系和三维直 用。这些新一代卫星重力测量项目与以往不同的 角坐标系中的表达式 重要特点，就是利用星载GPS接收机、K波段测 距仪等来确定重力卫星的轨道及其位置、速度和 众所周知，地球重力场在球坐标系中的表示 它们的变率，通过星载加速度计数据来排除非保 式U(a，9，，．)为‘1|： 守力对卫星轨道的扰动，综合利用这些数据就可 以推算地球重力场。这些星载仪器所提供的数据 7,7C=o＼7，罱 (1) 包括卫星轨道坐标及其速率，通常都是以(地心) +s…sin(rnA))P。。(sin妒) 三维直角坐标的形式表示的。因此，导出地球重 式中，A、饥r是点的球坐标；GM是地心引力常数； 力场、重力和重力梯度在三维直角坐标系中的表 R是地球平均半径；P…(sin9)是完全规格化的” 达式，在卫星重力解算中具有实际意义。本文所 阶m次勒让德函数；f…和R．。是地球重力场的完 推导的公式都是相应于ITRF定义下的地心地固全规格化的球谐函数的系数。相应于式(1)，地球 三维直角坐标系的．以下简称为三维直角坐标系。 重力场在三维直角坐标系中的表达式为∞1： …A1 ” (H—m)72 1／2． U(x，Y，z)=咝∑Rn∑‰。∑6。m．tz…2‘(z2+Y2+Z2)k--n ， n：0 m=0 女=0 ，+2 ～ 一． (一1)72 {气。∑(一1)7fl。it：1x,,r_Ziy加+‰。∑(一1)I (2) t。0 -20 9．“l‘，1／．m--(2i+1)Y2”1}