CN103093101A - 基于重力梯度误差模型原理的卫星重力反演方法 - Google Patents
基于重力梯度误差模型原理的卫星重力反演方法 Download PDFInfo
- Publication number
- CN103093101A CN103093101A CN2013100241728A CN201310024172A CN103093101A CN 103093101 A CN103093101 A CN 103093101A CN 2013100241728 A CN2013100241728 A CN 2013100241728A CN 201310024172 A CN201310024172 A CN 201310024172A CN 103093101 A CN103093101 A CN 103093101A
- Authority
- CN
- China
- Prior art keywords
- theta
- lambda
- gravity gradient
- sigma
- overbar
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明涉及一种地球重力场精密测量的方法,特别是一种基于重力梯度误差模型原理的卫星重力反演方法;通过分析卫星重力梯度垂向张量误差、水平张量误差和相关张量误差对累计大地水准面精度的联合影响建立新型重力梯度误差模型,进而精确和快速反演地球重力场;该方法卫星重力梯度反演精度高,地球重力场解算速度快,易于重力梯度卫星系统误差分析,卫星观测方程物理含义明确,计算机性能要求低;重力梯度误差模型法是反演高精度和高空间分辨率地球重力场的有效方法。
Description
技术领域
本发明涉及卫星重力梯度学、大地测量学、地球物理学、航空航天等交叉技术领域,特别是涉及一种通过分析卫星重力梯度垂向张量误差、水平张量误差和相关张量误差对累计大地水准面精度的联合影响建立新型重力梯度误差模型,进而精确和快速反演地球重力场的方法。
背景技术
地球重力场及其时变反映地球表层及内部物质的空间分布、运动和变化,同时决定着大地水准面的起伏和变化。因此,确定地球重力场的精细结构及其时变不仅是大地测量学、地球动力学、海洋学、冰川学、空间科学、国防建设等的需求,同时也将为寻求资源、保护环境和预测灾害提供重要的信息资源。
不同于GRACE(Gravity Recovery and Climate Experiment)双星高精度感测地球中长波重力场,欧空局(ESA)提出了专用于地球中短波重力场精密探测的GOCE(Gravity Field and Steady-State Ocean Circulation Explorer)卫星重力梯度(SGG)计划。如图1所示,GOCE卫星已于2009年3月17日成功发射升空,采用近圆(轨道离心率0.001)、极地(轨道倾角96.5°)和太阳同步轨道,经过3年的飞行计划,轨道高度由250km降为240km。GOCE卫星采用卫星跟踪卫星高低模式和卫星重力梯度模式的结合(SST-HL/SGG),除基于高轨道的GPS和GLONASS卫星对低轨道的GOCE进行精密跟踪定位(定轨精度1cm),同时利用定位于卫星质心处的星载重力梯度仪(测量精度3×10-12/s2)高精度测量卫星轨道高度处引力位的二阶导数。GOCE采用了非保守力补偿技术(Drag-free),首先利用重力梯度仪测量由非保守力(大气阻力、太阳光压、地球辐射压、轨道高度和姿态控制力等)引起的卫星质心的线性加速度与卫星平台的角加速度;最后,结合卫星平台姿态测量数据,通过无阻尼离子微推进器补偿卫星受到的非保守力。由于卫星重力梯度测量数据中的非保守力效应得到了有效扣除,因此进一步提高了地球重力场反演的精度和空间分辨率。自20世纪初匈牙利物理学家设计出第一台重力梯度仪(扭称)以来,重力梯度仪经历了从单轴旋转到三轴定向,从室温到低温(低于4.2k),从扭力、静电悬浮、超导到冷原子干涉的发展过程,测量精度日益提高。由于地球重力场信号随卫星轨道高度的增加而急剧衰减(Re/r)l,基于分析卫星轨道运动仅适合于确定地球中长波重力场,而卫星重力梯度是直接测定地球引力位的二次微分,其结果将球谐系数放大了l2倍,因此可有效抑制地球引力位随高度的衰减效应,进而高精度感测地球中高频重力场信号。欧空局独立研制的GOCE卫星原计划于2004年6月发射,由于星载三维静电悬浮重力梯度仪未能达到预期精度指标3×10-12/s2(单个加速度计的分辨率超过10-13m/s2,较GRACE卫星加速度计分辨率高约3个数量级)以及重力梯度卫星整体系统研制的困难性,因此距成功发射为止已推迟至少6次之多。
基于GOCE卫星重力梯度测量计划预计于2014年前结束,以及为了进一步提高地球重力场中短波信号的探测精度,目前国际众多科研机构正积极推动GOCE Follow-On卫星重力梯度测量计划的成功实施。我国相关研究机构紧跟国际卫星重力梯度测量的热点和动态,正积极投身于卫星重力梯度测量计划的需求论证和载荷预研之中。目前,国内外众多学者在基于卫星重力梯度技术反演地球重力场的理论和方法等方面已开展了广泛研究。不同于前人的研究,本发明通过分析卫星重力梯度垂向张量误差、水平张量误差和相关张量误差对累计大地水准面精度的联合影响建立新型重力梯度误差模型,进而精确和快速反演地球重力场。本发明不仅可为我国下一代地球重力场模型的精确建立提供理论基础和技术保证,同时对卫星重力梯度反演技术的发展方向具有一定借鉴意义。
在现有技术中已有通过分析一维垂向分量Vzz和三维全张量Vij重力梯度对累计大地水准面精度的影响,建立卫星重力梯度解析误差模型,进而精确和快速估计GOCE地球重力场精度的方法,该方法假设全球均匀分布的重力梯度观测值有N0个,并假设N0个观测值的误差满足正态分布随机特性,大量数据的平均可有效降低噪声,因此地球引力位系数的方差正比于1/N0。同时该方法未考虑垂向重力梯度和水平重力梯度之间所存在的较强相关性,未考虑基于相关性重力梯度误差反演累计大地水准面精度,因此其反演精度仍未达到预期的要求。
发明内容
本发明的目的是:通过分析卫星重力梯度垂向张量误差、水平张量误差和相关张量误差对累计大地水准面精度的联合影响,建立新型重力梯度误差模型,进而精确和快速反演地球重力场。
为达到上述目的,本发明采用了如下技术方案:
一种基于重力梯度误差模型的卫星重力反演方法,包括如下步骤:
步骤1,通过重力梯度卫星的星载重力梯度仪采集卫星重力梯度误差数据δTxyz;
步骤2,建立重力梯度误差模型,具体包括:
步骤2.1,在地固系中,将地球扰动位T(r,θ,λ)按球谐函数展开,分别对重力梯度卫星轨道位置矢量r的三个分量x,y,z进行二阶求导,其中θ和λ分别表示重力梯度卫星的地心余纬度和地心经度,将地球扰动位T(r,θ,λ)分别对r,θ,λ进行一阶求导和二阶求导;
步骤2.2,通过地球扰动位T(r,θ,λ)分别对x,y,z的二阶导数和分别对r,θ,λ一阶导数、二阶导数以及Legendre函数及其一阶导数和二阶导数得到待求规格化地球引力位系数和与一维垂向重力梯度Tzz的关系,将此关系表示为垂向重力梯度公式;
步骤2.3,利用所得到的垂向重力梯度公式以及地球引力位系数精度与累计大地水准面精度的关系获得基于垂向重力梯度误差数据δTzz反演累计大地水准面精度的误差模型;
步骤2.4,基于球谐函数的正交性,得到待求规格化地球引力位系数和与一维水平重力梯度Txx和Tyy的关系,将此关系表示为水平重力梯度公式,利用所得到的水平重力梯度公式以及地球引力位系数精度与累计大地水准面精度的关系分别获得基于水平重力梯度误差数据δTxx和δTyy反演累计大地水准面精度的误差模型;
步骤2.5,根据垂向重力梯度Tzz与水平重力梯度Txx和Tyy三者之间的非相互独立性,联合垂向重力梯度误差数据δTzz反演累计大地水准面精度的误差模型与水平重力梯度误差数据δTxx和δTyy反演累计大地水准面精度的误差模型,得到基于相关性重力梯度误差数据δTz-x-y反演累计大地水准面精度的误差模型,忽略非对角张量对地球重力场精度的影响,得到基于卫星重力梯度全张量误差数据δTxyz反演累计大地水准面精度的误差模型,以此作为重力梯度误差模型;
步骤3,基于所述重力梯度误差模型反演地球重力场精度,具体包括:
步骤3.1,在地球表面绘制网格;其次,按照重力梯度卫星的轨道在地球表面的轨迹点位置依次加入卫星重力梯度误差数据δTxyz,将分布于地球表面的卫星重力梯度误差数据δTxyz平均归算于划分的网格点δTxyz(φ,λ)处;
步骤3.2,将δTxyz(φ,λ)按球谐函数展开,得到δTxyz(φ,λ)按球谐函数展开的系数与卫星重力梯度误差数据在各阶处方差之间的关系,将卫星重力梯度误差数据在各阶处方差代入所述重力梯度误差模型进而反演地球重力场。
本发明是基于新型重力梯度误差模型法有利于快速反演高精度和高空间分辨率地球重力场的特点而设计的,优点为:
1)卫星重力梯度反演精度高;
2)地球重力场解算速度快;
3)易于重力梯度卫星系统误差分析;
4)卫星观测方程物理含义明确;
5)计算机性能要求低。
附图说明
图1表示GOCE卫星重力梯度计划。
图2表示基于新型重力梯度误差模型法反演累计大地水准面精度对比。
具体实施方式
以下结合附图,对本发明的具体实施方式作进一步的说明。
基于重力梯度误差模型原理的卫星重力反演方法包含下列步骤:
步骤一:重力梯度卫星数据采集
通过重力梯度卫星的星载重力梯度仪采集卫星重力梯度误差数据δTxyz。
步骤二:重力梯度误差模型建立
在地固系中,地球扰动位T(r,θ,λ)按球谐函数展开的表达式为
其中,GM表示地球质量M和万有引力常数G之积,Re表示地球的平均半径,L表示球函数展开的最大阶数;表示卫星的地心半径,x,y,z分别表示卫星轨道位置矢量r的三个分量,θ和λ分别表示卫星的地心余纬度和地心经度;表示规格化的Legendre函数,l表示阶数,m表示次数;和表示待求的规格化地球引力位系数。
T(r,θ,λ)分别对x,y,z的二阶导数表示为
其中,地球扰动位二阶导数是对称张量,同时在真空情况下满足Laplace方程表现为无迹性,Txx+Tyy+Tzz=0,因此在9个卫星重力梯度分量中有5个是独立的。全张量重力梯度的9个分量表示为
其中,地球扰动位T(r,θ,λ)分别对r,θ,λ的一阶导数表示为
地球扰动位T(r,θ,λ)分别对r,θ,λ的二阶导数表示为
Legendre函数及其一阶导数和二阶导数表示为
其中,
基于球谐函数的正交性,联合公式(3)和公式(5)可得一维垂向重力梯度公式
其中,Tzz表示一维垂向重力梯度,实际计算时需要离散化数值积分。基于等间隔的Δθ和Δλ在地球表面进行全球经纬网划分,同时将每个网格中的垂向重力梯度值取平均值其中i,j表示网格的经纬度标号。因此,公式(7)可改写为
累计大地水准面精度公式表示为
其中,表示地球引力位系数精度。
联合公式(8)和公式(9),可得基于一维垂向重力梯度误差数据δTzz反演累计大地水准面精度的误差模型
基于球谐函数的正交性,联合公式(3)~(5)可得水平方向重力梯度Txx(yy)公式
联合公式(9)和公式(11),可得分别基于水平重力梯度误差数据δTxx和δTyy反演累计大地水准面精度的误差模型
由于垂向重力梯度Tzz以及水平重力梯度Txx和Tyy非相互独立,而具有较强的相关性,因此,联合公式(10)和公式(12)可得基于相关性重力梯度误差数据δTz-x-y反演累计大地水准面精度的误差模型
如公式(2)所示,在卫星重力梯度的9个张量中,对角张量(垂向分量Tzz和水平分量Txx,Tyy)是主要分量,非对角张量对地球重力场精度的影响相对于对角张量基本可忽略。因此,联合公式(10)、(12)和(13),可得基于卫星重力梯度全张量误差数据δTxyz反演累计大地水准面精度的误差模型
步骤三:卫星重力梯度反演
基于新型卫星重力梯度误差模型法,利用2012年的GOCE-Level-1B卫星重力梯度误差数据δTxyz反演累计大地水准面精度的过程如下
第一,首先以0.3°×0.3°为网格分辨率,在地球表面的经度(0°~360°)和纬度(-90°~90°)范围内绘制网格;其次,按照GOCE卫星轨道在地球表面的轨迹点位置依次加入δTxyz;最后,将分布于地球表面的δTxyz平均归算于划分的网格点δTxyz(φ,λ)处。
第二,将δTxyz(φ,λ)按球谐函数展开为
δTxyz在各阶处的方差表示为
将公式(17)代入公式(14),可有效和快速反演地球重力场精度。
图2表示基于新型重力梯度误差模型法反演累计大地水准面精度对比(公式(14)),虚线表示未加入相关性重力梯度误差δTz-x-y反演累计大地水准面精度,实线表示加入相关性重力梯度误差δTz-x-y(公式(13))反演累计大地水准面精度。通过对图2所示结果的对比研究表明:第一,加入相关性重力梯度误差δTz-x-y反演累计大地水准面精度(实线)较未加入相关性重力梯度误差δTz-x-y反演精度平均提高2~3倍;第二,由于实际上垂向重力梯度Tzz以及水平重力梯度Txx和Tyy并非相互独立,而是具有较强的相关性,因此,在本发明建立的新型重力梯度误差模型中加入相关性重力梯度误差δTz-x-y是进一步提高地球重力场精度的重要因素;第三,新型重力梯度误差模型法是建立下一代高精度、高空间分辨率和高阶次地球重力场模型的有效方法。
Claims (5)
1.一种基于重力梯度误差模型的卫星重力反演方法,其特征在于包括如下步骤:
步骤1,通过重力梯度卫星的星载重力梯度仪采集卫星重力梯度误差数据δTxyz;
步骤2,建立重力梯度误差模型,具体包括:
步骤2.1,在地固系中,将地球扰动位T(r,θ,λ)按球谐函数展开,分别对重力梯度卫星轨道位置矢量r的三个分量x,y,z进行二阶求导,其中θ和λ分别表示重力梯度卫星的地心余纬度和地心经度,将地球扰动位T(r,θ,λ)分别对r,θ,λ进行一阶求导和二阶求导;
步骤2.2,通过地球扰动位T(r,θ,λ)分别对x,y,z的二阶导数和分别对r,θ,λ一阶导数、二阶导数以及Legendre函数及其一阶导数和二阶导数得到待求规格化地球引力位系数和与一维垂向重力梯度Tzz的关系,将此关系表示为垂向重力梯度公式;
步骤2.4,基于球谐函数的正交性,得到待求规格化地球引力位系数和与一维水平重力梯度Txx和Tyy的关系,将此关系表示为水平重力梯度公式,利用所得到的水平重力梯度公式以及地球引力位系数精度与累计大地水准面精度的关系分别获得基于水平重力梯度误差数据δTxx和δTyy反演累计大地水准面精度的误差模型;
步骤2.5,根据垂向重力梯度Tzz与水平重力梯度Txx和Tyy三者之间的非相互独立性,联合垂向重力梯度误差数据δTzz反演累计大地水准面精度的误差模型与水平重力梯度误差数据δTxx和δTyy反演累计大地水准面精度的误差模型,得到基于相关性重力梯度误差数据δTz-x-y反演累计大地水准面精度的误差模型,忽略非对角张量对地球重力场精度的影响,得到基于卫星重力梯度全张量误差数据δTxyz反演累计大地水准面精度的误差模型,以此作为重力梯度误差模型;
步骤3,基于所述重力梯度误差模型反演地球重力场精度,具体包括:
步骤3.1,在地球表面绘制网格;其次,按照重力梯度卫星的轨道在地球表面的轨迹点位置依次加入卫星重力梯度误差数据δTxyz,将分布于地球表面的卫星重力梯度误差数据δTxyz平均归算于划分的网格点δTxyz(φ,λ)处;
2.如权利要求1所述的基于重力梯度误差模型的卫星重力反演方法,其特征在于所述步骤2具体包括:
在地固系中,将地球扰动位T(r,θ,λ)按球谐函数展开表达为
其中,GM表示地球质量M和万有引力常数G之积,Re表示地球的平均半径,L表示球谐函数展开的最大阶数;表示卫星的地心半径,x,y,z分别表示卫星轨道位置矢量r的三个分量,θ和λ分别表示卫星的地心余纬度和地心经度;表示规格化的Legendre函数,l表示阶数,m表示次数;和表示待求的规格化地球引力位系数;
T(r,θ,λ)分别对x,y,z的二阶导数表示为
其中,地球扰动位二阶导数是对称张量,同时在真空情况下满足Laplace方程表现为无迹性,Txx+Tyy+Tzz=0,因此在9个卫星重力梯度分量中有5个是独立的;全张量重力梯度的9个分量表示为
其中,地球扰动位T(r,θ,λ)分别对r,θ,λ的一阶导数表示为
地球扰动位T(r,θ,λ)分别对r,θ,λ的二阶导数表示为
Legendre函数及其一阶导数和二阶导数表示为
其中,
基于球谐函数的正交性,联合公式(3)和公式(5)可得一维垂向重力梯度公式
其中,Tzz表示一维垂向重力梯度,实际计算时需要离散化数值积分;基于等间隔的Δθ和Δλ在地球表面进行全球经纬网划分,同时将每个网格中的垂向重力梯度值取平均值其中i,j表示网格的经纬度标号;因此,公式(7)可改写为
累计大地水准面精度公式表示为
联合公式(8)和公式(9),可得基于一维垂向重力梯度误差数据δTzz反演累计大地水准面精度的误差模型
基于球谐函数的正交性,联合公式(3)~(5)可得水平方向重力梯度Txx(yy)公式
联合公式(9)和公式(11),可得分别基于水平重力梯度误差数据δTxx和δTyy反演累计大地水准面精度的误差模型
由于垂向重力梯度Tzz以及水平重力梯度Txx和Tyy非相互独立,而具有较强的相关性,因此,联合公式(10)和公式(12)可得基于相关性重力梯度误差数据δTz-x-y反演累计大地水准面精度的误差模型
如公式(2)所示,在卫星重力梯度的9个张量中,对角张量(垂向分量Tzz和水平分量Txx,Tyy)是主要分量,非对角张量对地球重力场精度的影响相对于对角张量基本可忽略;因此,联合公式(10)、(12)和(13),可得基于卫星重力梯度全张量误差数据δTxyz反演累计大地水准面精度的误差模型
4.如权利要求3所述的基于重力梯度误差模型的卫星重力反演方法,其特征在于:所确定的网格分辨率为0.3°×0.3°。
5.如权利要求1-4中任意一项所述的基于重力梯度误差模型的卫星重力反演方法,其特征在于:所述重力梯度卫星为GOCE卫星或GOCE Follow-On卫星,优选为GOCE Follow-On卫星。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310024172.8A CN103093101B (zh) | 2013-01-22 | 2013-01-22 | 基于重力梯度误差模型原理的卫星重力反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310024172.8A CN103093101B (zh) | 2013-01-22 | 2013-01-22 | 基于重力梯度误差模型原理的卫星重力反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103093101A true CN103093101A (zh) | 2013-05-08 |
CN103093101B CN103093101B (zh) | 2015-08-26 |
Family
ID=48205660
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310024172.8A Expired - Fee Related CN103093101B (zh) | 2013-01-22 | 2013-01-22 | 基于重力梯度误差模型原理的卫星重力反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103093101B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015042754A1 (zh) * | 2013-09-29 | 2015-04-02 | 清华大学 | 一种低低星星跟踪卫星重力场测量性能解析计算方法 |
CN106529164A (zh) * | 2016-11-03 | 2017-03-22 | 清华大学 | 联合重力卫星获取地下水储量变化值的方法及系统 |
CN108415879A (zh) * | 2018-01-19 | 2018-08-17 | 中国人民解放军92859部队 | 基于向上延拓的航空重力最小二乘向下延拓解析算法 |
CN110967041A (zh) * | 2019-12-18 | 2020-04-07 | 自然资源部国土卫星遥感应用中心 | 基于张量不变理论的卫星引力梯度数据精度的验证方法 |
CN113885101A (zh) * | 2021-09-28 | 2022-01-04 | 中国船舶重工集团公司第七0七研究所 | 一种基于椭球模型构建重力梯度基准图方法 |
CN114089432A (zh) * | 2021-11-10 | 2022-02-25 | 中国地质大学(北京) | 一种利用卫星测高数据反演海洋重力梯度的频率域方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000060379A1 (en) * | 1999-04-02 | 2000-10-12 | Conoco, Inc. | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
CN101975969A (zh) * | 2010-10-19 | 2011-02-16 | 华中科技大学 | 一种基于全张量重力梯度反演的水下目标探测方法 |
CN102262248A (zh) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | 基于双星空间三维插值原理的卫星重力反演方法 |
CN102305949A (zh) * | 2011-06-30 | 2012-01-04 | 中国科学院测量与地球物理研究所 | 利用星间距离插值建立全球重力场模型的方法 |
CN102313905A (zh) * | 2011-07-18 | 2012-01-11 | 中国科学院测量与地球物理研究所 | 基于星间速度插值原理的卫星重力反演方法 |
CN102393535A (zh) * | 2011-07-20 | 2012-03-28 | 中国科学院测量与地球物理研究所 | 基于双星能量插值原理的卫星重力反演方法 |
CN102608668A (zh) * | 2011-12-19 | 2012-07-25 | 华中科技大学 | 重力梯度测量系统及其测量方法 |
-
2013
- 2013-01-22 CN CN201310024172.8A patent/CN103093101B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000060379A1 (en) * | 1999-04-02 | 2000-10-12 | Conoco, Inc. | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
CN101975969A (zh) * | 2010-10-19 | 2011-02-16 | 华中科技大学 | 一种基于全张量重力梯度反演的水下目标探测方法 |
CN102262248A (zh) * | 2011-06-03 | 2011-11-30 | 中国科学院测量与地球物理研究所 | 基于双星空间三维插值原理的卫星重力反演方法 |
CN102305949A (zh) * | 2011-06-30 | 2012-01-04 | 中国科学院测量与地球物理研究所 | 利用星间距离插值建立全球重力场模型的方法 |
CN102313905A (zh) * | 2011-07-18 | 2012-01-11 | 中国科学院测量与地球物理研究所 | 基于星间速度插值原理的卫星重力反演方法 |
CN102393535A (zh) * | 2011-07-20 | 2012-03-28 | 中国科学院测量与地球物理研究所 | 基于双星能量插值原理的卫星重力反演方法 |
CN102608668A (zh) * | 2011-12-19 | 2012-07-25 | 华中科技大学 | 重力梯度测量系统及其测量方法 |
Non-Patent Citations (3)
Title |
---|
L. VÖLGYESI,M. DOBRÓKA,Z. ULTMANN: "Determination of vertical gradients of gravity by series expansion based inversion", 《ACTA GEODAETICA ET GEOPHYSICA HUNGARICA》 * |
王虎彪,王勇,陆洋,鲍李峰: "联合多种测高资料确定西太平洋海域2′×2′重力梯度", 《地球物理学进展》 * |
郇恒飞,吴燕冈,管彦武,杨长保: "利用实测重力垂直梯度反演长白山地区一剖面的深部构造", 《世界地质》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015042754A1 (zh) * | 2013-09-29 | 2015-04-02 | 清华大学 | 一种低低星星跟踪卫星重力场测量性能解析计算方法 |
DE112013007459B4 (de) * | 2013-09-29 | 2017-12-28 | Tsinghua University | Leistungsanalyse-und-berechnungsverfahren für Gravitationsfeldmessung aus Niedrig-Niedrig-Satellite-to-Satellite-Tracking (SST) |
CN106529164A (zh) * | 2016-11-03 | 2017-03-22 | 清华大学 | 联合重力卫星获取地下水储量变化值的方法及系统 |
CN106529164B (zh) * | 2016-11-03 | 2019-01-04 | 清华大学 | 联合重力卫星获取地下水储量变化值的方法及系统 |
CN108415879A (zh) * | 2018-01-19 | 2018-08-17 | 中国人民解放军92859部队 | 基于向上延拓的航空重力最小二乘向下延拓解析算法 |
CN108415879B (zh) * | 2018-01-19 | 2021-04-06 | 中国人民解放军92859部队 | 基于向上延拓的航空重力最小二乘向下延拓解析方法 |
CN110967041A (zh) * | 2019-12-18 | 2020-04-07 | 自然资源部国土卫星遥感应用中心 | 基于张量不变理论的卫星引力梯度数据精度的验证方法 |
CN113885101A (zh) * | 2021-09-28 | 2022-01-04 | 中国船舶重工集团公司第七0七研究所 | 一种基于椭球模型构建重力梯度基准图方法 |
CN113885101B (zh) * | 2021-09-28 | 2023-12-12 | 中国船舶重工集团公司第七0七研究所 | 一种基于椭球模型构建重力梯度基准图方法 |
CN114089432A (zh) * | 2021-11-10 | 2022-02-25 | 中国地质大学(北京) | 一种利用卫星测高数据反演海洋重力梯度的频率域方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103093101B (zh) | 2015-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103076640B (zh) | 利用方差-协方差对角张量原理反演地球重力场的方法 | |
CN103093101B (zh) | 基于重力梯度误差模型原理的卫星重力反演方法 | |
CN102305949B (zh) | 利用星间距离插值建立全球重力场模型的方法 | |
CN102262248B (zh) | 基于双星空间三维插值原理的卫星重力反演方法 | |
Fuchs et al. | Observing coseismic gravity change from the Japan Tohoku‐Oki 2011 earthquake with GOCE gravity gradiometry | |
CN102393535B (zh) | 基于双星能量插值原理的卫星重力反演方法 | |
CN102313905B (zh) | 基于星间速度插值原理的地球重力反演方法 | |
CN103163562B (zh) | 基于滤波原理的卫星重力梯度反演方法 | |
CN103018783B (zh) | 重力卫星编队轨道稳定性优化设计和精密反演地球重力场方法 | |
CN102736118B (zh) | 一种用于全球重力场测量的综合型卫星系统 | |
CN102998713B (zh) | 基于功率谱半解析的卫星重力梯度反演方法 | |
CN103091722B (zh) | 基于载荷误差分析原理的卫星重力反演方法 | |
Pan et al. | GPS imaging of vertical bedrock displacements: quantification of two‐dimensional vertical crustal deformation in China | |
Zheng et al. | Requirements analysis for future satellite gravity mission Improved-GRACE | |
CN103076639B (zh) | 基于残余星间速度原理反演地球重力场的方法 | |
Fuchs et al. | Distributed fault slip model for the 2011 Tohoku‐Oki earthquake from GNSS and GRACE/GOCE satellite gravimetry | |
CN103091723B (zh) | 降低重力卫星质心调整误差对地球重力场精度影响的方法 | |
CN103091721B (zh) | 利用不同轨道倾角卫星联合反演地球重力场的方法 | |
Sun et al. | Progress of research on the earth’s gravity tides and its application in geodynamics in China | |
CN103064128B (zh) | 基于星间距离误差模型的地球重力场恢复方法 | |
Zhou et al. | An improved gravity compensation method for high-precision free-INS based on MEC–BP–AdaBoost | |
Xing et al. | Far‐field coseismic gravity changes related to the 2015 MW7. 8 Nepal (Gorkha) earthquake observed by superconducting gravimeters in mainland China | |
Frédéric et al. | Contribution of GRACE satellite gravimetry in global and regional hydrology, and in ice sheets mass balance | |
El-Ashquer | An improved hybrid local geoid model for Egypt | |
CN105549105B (zh) | 一种短基线相对轨道摄动重力场测量性能的评估方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150826 Termination date: 20160122 |
|
EXPY | Termination of patent right or utility model |