一种基于张量不变理论的卫星引力梯度粗差探测方法技术

本发明专利技术公开了一种基于张量不变理论的卫星引力梯度粗差探测方法，该方法包括整体探测卫星引力梯度观测张量粗差；独立探测卫星引力梯度观测张量各分量粗差。本发明专利技术公开的卫星引力梯度观测值粗差探测的单次差方法，以张量不变理论为基础，将卫星引力梯度观测值的张量不变特性以及第二不变量和第三不变量的放大效应应用于卫星引力梯度观测值的粗差探测中，实现了卫星引力梯度观测张量整体粗差的可靠性探测；通过引入先验的重力场模型，可独立实现引力梯度张量主对角线和非对角线六个分量的粗差探测。粗差探测。粗差探测。

【技术实现步骤摘要】
一种基于张量不变理论的卫星引力梯度粗差探测方法


[0001]本专利技术涉及大地测量
，尤其涉及一种基于张量不变理论的卫星引力梯度观测值粗差探测的单次差方法。

技术介绍

[0002]卫星引力梯度数据对于确定地球重力场中短波精细结构至关重要。空间物理环境的变化导致卫星引力梯度观测值中不可避免的含有一定的粗差。观测值粗差的存在影响观测值的数学统计特性，带来最小二乘平差等理论方法的不适用性和解算过程的病态，使得卫星引力梯度观测数据的可利用性变差，制约了地球重力场中短波频谱精度的提升。为此，卫星引力梯度观测值需要进行粗差探测和剔除。
[0003]卫星引力梯度观测值的粗差探测方法多基于卫星引力梯度观测张量某一分量，利用各种数理统计方法开展。该处理过程未考虑卫星引力梯度观测张量各分量存在的内在关系特性，无法对卫星引力梯度张量进行整体探测，且探测过程较为繁琐，各种数理统计方法往往对某一类粗差较为敏感，对各类粗差的兼容性探测能力较差，从而导致粗差探测效率和探测精度欠佳，间接带来卫星数据的使用性能差，影响了地球重力场中短波精细结构的精确测定。

技术实现思路

[0004]为解决上述技术问题，本专利技术的目的是提供一种基于张量不变理论的卫星引力梯度粗差探测方法，该方法利用张量不变量的幅值不变特性和放大效应，解决目前卫星引力梯度观测值的粗差探测效能偏低而导致的卫星观测值质量差乃至无法使用的问题。
[0005]本专利技术的目的通过以下的技术方案来实现：
[0006]一种基于张量不变理论的卫星引力梯度粗差探测方法，包括：
[0007]步骤A整体探测卫星引力梯度观测张量粗差；
[0008]步骤B独立探测卫星引力梯度观测张量各分量粗差。
[0009]与现有技术相比，本专利技术的一个或多个实施例可以具有如下优点：
[0010]本专利技术公开的卫星引力梯度观测值粗差探测的单次差方法，以张量不变理论为基础，将卫星引力梯度观测值的张量不变特性以及第二不变量和第三不变量的放大效应应用于卫星引力梯度观测值的粗差探测中，实现了卫星引力梯度观测张量整体粗差的可靠性探测；通过引入先验的重力场模型，可独立实现引力梯度张量主对角线和非对角线六个分量的粗差探测。
附图说明
[0011]图1是基于张量不变理论的卫星引力梯度粗差探测方法流程图；
[0012]图2是卫星引力梯度观测张量粗差的整体探测流程图；
[0013]图3是卫星引力梯度观测张量各分量粗差的独立探测流程图。
具体实施方式
[0014]为使本专利技术的目的、技术方案和优点更加清楚，下面将结合实施例及附图对本专利技术作进一步详细的描述。
[0015]如图1所示，为基于张量不变理论的卫星引力梯度粗差探测方法，包括：
[0016]步骤10整体探测卫星引力梯度观测张量粗差；
[0017]步骤20独立探测卫星引力梯度观测张量各分量粗差。
[0018]如图2所示，上述步骤10对卫星引力梯度观测张量粗差的整体探测具体包括以下步骤：
[0019]步骤101建立卫星引力梯度观测值的张量不变系统{I1,I2,I3}；
[0020]卫星引力梯度观测值张量不变系统{I1,I2,I3}的表达式为：
[0021]I1＝V
11
+V
22
+V
33
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0022][0023][0024]式中：I1为卫星引力梯度张量不变系统的第一不变量；I2为卫星引力梯度张量不变系统第二不变量；I3为卫星引力梯度张量不变系统第三不变量；V
11
为卫星梯度仪坐标系下xx方向的卫星引力梯度观测值分量；V
12
为卫星梯度仪坐标系下xy方向的卫星引力梯度观测值分量；V
13
为卫星梯度仪坐标系下xz方向的卫星引力梯度观测值分量；V
22
为卫星梯度仪坐标系下yy方向的卫星引力梯度观测值分量；V
23
为卫星梯度仪坐标系下yz方向的卫星引力梯度观测值分量；V
33
为卫星梯度仪坐标系下zz方向的卫星引力梯度观测值分量。考虑到第一不变量具有迹无关特性，即I1＝V
11
+V
22
+V
33
＝0，本专利技术提出的卫星引力梯度观测值粗差探测的单次差方法仅涉及第二不变量和第三不变量。
[0025]步骤102计算卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差；
[0026]将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0027][0028][0029]对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0030][0031][0032]步骤103卫星引力梯度观测张量整体粗差的探测；
[0033]设定粗差探测窗口M和显著因子α，利用如下准则进行卫星引力梯度观测张量整体粗差的探测：
[0034][0035][0036]对于不满足准则条件的引力梯度张量即认为存在粗差。由此，通过以上过程实现了对引力梯度观测张量的整体粗差探测。
[0037]如图3所示，上述步骤20具体包括：
[0038]步骤201计算卫星梯度仪坐标系下卫星引力梯度张量先验模型值；
[0039]步骤202独立探测卫星引力梯度观测张量各分量粗差。
[0040]所述步骤201具体包括：
[0041]步骤2011局部指北坐标下卫星引力梯度张量先验模型值的计算；
[0042]利用先验重力场模型计算局部指北坐标系下引力梯度张量先验模型值计算公式为：
[0043][0044]式中：GM为地心引力常数，r、θ、λ分别为地心向径、地心余纬和地心经度，R为地球平均半径，n、m为球谐模型展开的阶和次，为引力梯度分量值，λ
ij
、为引力梯度张量系数，其表达式参见表1，i＝1,2,3分别表示i＝x,y,z方向，j＝1,2,3分别表示j＝x,y,z方向；N为先验重力场模型的最高阶数。局部指北坐标系下的引力梯度分量表达式如表1：
[0045]表1
[0046][0047]表1中，为完全规格化缔合勒让德函数，为重力场模型的完全规
格化的引力位球谐系数，和为勒让德函数对地心余纬θ的一、二阶导数。
[0048]将某一观测历元k所对应的地心向径、地心余纬和地心经度以及先验重力场模型系数代入(10)，即可获得观测历元k局部指北坐标下卫星引力梯度张量的先验模型值
[0049]步骤2012卫星梯度仪坐标系下卫星引力梯度张量先验模型值的计算；
[0050]局部指北坐标下卫星引力梯度张量先验模型值转换至卫星梯度仪坐标系，需要经本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于张量不变理论的卫星引力梯度粗差探测方法，其特征在于，所述方法包括以下步骤：步骤A整体探测卫星引力梯度观测张量粗差；步骤B独立探测卫星引力梯度观测张量各分量粗差。2.如权利要求1所述的基于张量不变理论的卫星引力梯度粗差探测方法，其特征在于，所述步骤A具体包括：A1建立卫星引力梯度观测值的张量不变系统{I1,I2,I3}；A2计算卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差；A3卫星引力梯度观测张量整体粗差的探测。3.如权利要求2所述的基于张量不变理论的卫星引力梯度粗差探测方法，其特征在于，所述A1中，卫星引力梯度观测值张量不变系统{I1,I2,I3}的表达式为：I1＝V
11
+V
22
+V
33
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)(1)式中：I1为卫星引力梯度张量不变系统的第一不变量；I2为卫星引力梯度张量不变系统第二不变量；I3为卫星引力梯度张量不变系统第三不变量；V
11
为卫星梯度仪坐标系下xx方向的卫星引力梯度观测值分量；V
12
为卫星梯度仪坐标系下xy方向的卫星引力梯度观测值分量；V
13
为卫星梯度仪坐标系下xz方向的卫星引力梯度观测值分量；V
22
为卫星梯度仪坐标系下yy方向的卫星引力梯度观测值分量；V
23
为卫星梯度仪坐标系下yz方向的卫星引力梯度观测值分量；V
33
为卫星梯度仪坐标系下zz方向的卫星引力梯度观测值分量。4.如权利要求3所述的基于张量不变理论的卫星引力梯度粗差探测方法，其特征在于，所述A2具体包括：将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度观测值分量其中(i,j＝1,2,3)代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：变量和卫星引力梯度第三不变量：对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：三不变量的单次差：其中，(k＝1,2,3,
…
)。5.如权利要求1所述的基于张量不变理论的卫星引力梯度粗差探测方法，其特征在于，所述A3具体包括设定粗差探测窗口M和显著因子α，利用如下准则进行卫星引力梯度观测张量整体粗差
的探测：的探测：对于不满足准则条件的引力梯度张量即认为存在粗差，其中，k＝1,2,3,
…
。6.如权利要求1所述的基于张量不变理论的卫星引力梯度粗差探测方法，其特征在于，所述步骤B...

【专利技术属性】
技术研发人员：朱广彬，常晓涛，刘伟，唐洪钊，谢俊峰，窦显辉，莫凡，陈辉，
申请(专利权)人：自然资源部国土卫星遥感应用中心，
类型：发明
国别省市：