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

1.本发明涉及大地测量技术领域，尤其涉及一种基于张量不变理论的卫星引力梯度观测值粗差探测的单次差方法。


背景技术：

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


技术实现要素：

4.为解决上述技术问题，本发明的目的是提供一种基于张量不变理论的卫星引力梯度粗差探测方法，该方法利用张量不变量的幅值不变特性和放大效应，解决目前卫星引力梯度观测值的粗差探测效能偏低而导致的卫星观测值质量差乃至无法使用的问题。
5.本发明的目的通过以下的技术方案来实现：
6.一种基于张量不变理论的卫星引力梯度粗差探测方法，包括：
7.步骤a整体探测卫星引力梯度观测张量粗差；
8.步骤b独立探测卫星引力梯度观测张量各分量粗差。
9.与现有技术相比，本发明的一个或多个实施例可以具有如下优点：
10.本发明公开的卫星引力梯度观测值粗差探测的单次差方法，以张量不变理论为基础，将卫星引力梯度观测值的张量不变特性以及第二不变量和第三不变量的放大效应应用于卫星引力梯度观测值的粗差探测中，实现了卫星引力梯度观测张量整体粗差的可靠性探测；通过引入先验的重力场模型，可独立实现引力梯度张量主对角线和非对角线六个分量的粗差探测。
附图说明
11.图1是基于张量不变理论的卫星引力梯度粗差探测方法流程图；
12.图2是卫星引力梯度观测张量粗差的整体探测流程图；
13.图3是卫星引力梯度观测张量各分量粗差的独立探测流程图。
具体实施方式
14.为使本发明的目的、技术方案和优点更加清楚，下面将结合实施例及附图对本发明作进一步详细的描述。
15.如图1所示，为基于张量不变理论的卫星引力梯度粗差探测方法，包括：
16.步骤10整体探测卫星引力梯度观测张量粗差；
17.步骤20独立探测卫星引力梯度观测张量各分量粗差。
18.如图2所示，上述步骤10对卫星引力梯度观测张量粗差的整体探测具体包括以下步骤：
19.步骤101建立卫星引力梯度观测值的张量不变系统{i1,i2,i3}；
20.卫星引力梯度观测值张量不变系统{i1,i2,i3}的表达式为：
21.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]
局部指北坐标下卫星引力梯度张量先验模型值转换至卫星梯度仪坐标系，需要经过以下几个转换过程：局部指北坐标系
→
地固坐标系
→
惯性坐标系
→
梯度仪坐标系，转换公式为：
[0051][0052]
式中：和分别为坐标系转换前后的卫星引力梯度张量，r为坐标转换矩阵。
[0053]
将观测历元k局部指北坐标下卫星引力梯度张量的先验模型值以及该时刻的坐标转换矩阵rk代入(11)，即可获得观测历元k卫星梯度仪坐标系下卫星引力梯度张量先验模型值
[0054]
所述步骤202具体包括：
[0055]
步骤2021卫星引力梯度观测张量v
11
分量粗差的独立探测；
[0056]
将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度张量先验模型值以及卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0057][0058][0059]
对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0060][0061][0062]
设定粗差探测窗口m和显著因子α，利用如下准则进行卫星引力梯度观测张量v
11
分量粗差的探测：
[0063][0064]
[0065]
对于不满足准则条件的引力梯度张量v
11
分量即认为存在粗差，由此，通过以上过程实现了对引力梯度观测张量v
11
分量的粗差探测。
[0066]
步骤2022卫星引力梯度观测张量v
22
分量粗差的独立探测；
[0067]
将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度张量先验模型值以及卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0068][0069][0070]
对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0071][0072][0073]
设定粗差探测窗口m和显著因子α，利用如下准则进行卫星引力梯度观测张量v
22
分量粗差的探测：
[0074][0075][0076]
对于不满足准则条件的引力梯度张量v
22
分量即认为存在粗差，由此，通过以上过程实现了对引力梯度观测张量v
22
分量的粗差探测。
[0077]
步骤2023卫星引力梯度观测张量v
33
分量粗差的独立探测；
[0078]
将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度张量先验模型值以及卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0079][0080][0081]
对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0082]
[0083][0084]
设定粗差探测窗口m和显著因子α，利用如下准则进行卫星引力梯度观测张量v
33
分量粗差的探测：
[0085][0086][0087]
对于不满足准则条件的引力梯度张量v
33
分量即认为存在粗差，由此，通过以上过程实现了对引力梯度观测张量v
33
分量的粗差探测。
[0088]
步骤2024卫星引力梯度观测张量v
12
分量粗差的独立探测；
[0089]
将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度张量先验模型值以及卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0090][0091][0092]
对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0093][0094][0095]
设定粗差探测窗口m和显著因子α，利用如下准则进行卫星引力梯度观测张量v
12
分量粗差的探测：
[0096][0097][0098]
对于不满足准则条件的引力梯度张量v
12
分量即认为存在粗差，由此，通过以上过程实现了对引力梯度观测张量v
12
分量的粗差探测。
[0099]
步骤2025卫星引力梯度观测张量v
13
分量粗差的独立探测；
[0100]
将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度张量先验模型值以及卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0101]
[0102][0103]
对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0104][0105][0106]
设定粗差探测窗口m和显著因子α，利用如下准则进行卫星引力梯度观测张量v
13
分量粗差的探测：
[0107][0108][0109]
对于不满足准则条件的引力梯度张量v
13
分量即认为存在粗差，由此，通过以上过程实现了对引力梯度观测张量v
13
分量的粗差探测。
[0110]
步骤2026卫星引力梯度观测张量v
23
分量粗差的独立探测；
[0111]
将卫星梯度仪坐标系下不同观测历元k的卫星引力梯度张量先验模型值以及卫星引力梯度观测值分量代入张量不变系统的计算公式(2)和(3)中，得到不同观测历元k的卫星引力梯度第二不变量和卫星引力梯度第三不变量：
[0112][0113][0114]
对相邻历元的张量不变量作差，计算得到卫星引力梯度第二不变量和卫星引力梯度第三不变量的单次差：
[0115][0116][0117]
设定粗差探测窗口m和显著因子α，利用如下准则进行卫星引力梯度观测张量v
23
分量粗差的探测：
[0118][0119][0120]
对于不满足准则条件的引力梯度张量v
23
分量即认为存在粗差，由此，通过以上过程实现了对引力梯度观测张量v
23
分量的粗差探测。
[0121]
虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本
发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。