CN112965127B - 一种基于重力异常计算外部扰动重力径向分量的方法 - Google Patents
一种基于重力异常计算外部扰动重力径向分量的方法 Download PDFInfo
- Publication number
- CN112965127B CN112965127B CN202110180191.4A CN202110180191A CN112965127B CN 112965127 B CN112965127 B CN 112965127B CN 202110180191 A CN202110180191 A CN 202110180191A CN 112965127 B CN112965127 B CN 112965127B
- Authority
- CN
- China
- Prior art keywords
- gravity
- radial component
- anomaly
- point
- kernel function
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种基于重力异常计算外部扰动重力径向分量的方法,其主要技术特点是:利用移去‑恢复技术,从重力异常观测值中移去参考重力异常,得到残差重力异常;从积分式核函数中移去与参考场对应阶次的核函数球谐表达式,得到截断核函数,使之与残差重力异常的频谱匹配;基于截断核函数和残差重力异常的局域积分,得到残差外部扰动重力径向分量;利用全球重力场位模型高阶信息进行远区效应补偿;恢复参考外部扰动重力径向分量;利用积分恒等式转换,补偿由全球积分过渡到局域积分引起的模型误差,得到计算点处高精度的外部扰动重力径向分量。本发明解决了基于重力异常计算外部扰动重力径向分量的不确定问题,提高外部扰动重力径向分量的计算精度,可广泛用于物理大地测量领域。
Description
技术领域
本发明属于大地测量学与测量工程技术领域,尤其是一种基于重力异常计算外部扰动重力径向分量的方法。
背景技术
外部扰动重力径向分量是地球重力场逼近建模研究内容的重要组成部分,也是解算大地测量边值问题的主要应用目标之一,在航空航天器飞行轨迹精密计算和空间科学技术研究中具有重要的应用价值。
在基于重力异常计算外部扰动重力径向分量的积分式中,当计算点趋近于观测数据格网点时积分核函数会发生奇异,导致利用积分式计算超低空外部扰动重力径向分量时,会出现由核函数奇异引起的奇异性问题,无法保证计算结果的有效性。同时,在实际应用中受观测数据覆盖范围限制,无法做到全球覆盖,实际计算过程中外部扰动重力径向分量的全球积分式需要进行适用观测数据保障条件的改化,以保证计算结果的可靠性。目前迫切需要解决基于重力异常计算外部扰动重力径向分量的不确定问题,以提高外部扰动重力径向分量的计算精度。
发明内容
本发明的目的在于克服现有技术的不足,提出一种基于重力异常计算外部扰动重力径向分量的方法,解决基于重力异常计算外部扰动重力径向分量的不确定问题,提高外部扰动重力径向分量的计算精度。
本发明解决其技术问题是采取以下技术方案实现的:
一种基于重力异常计算外部扰动重力径向分量的方法,包括以下步骤:
步骤1、建立基于重力异常计算外部扰动重力径向分量δgr的积分式,恢复积分核函数的零阶和一阶项影响,得到完整核函数Ftr(r,ψ);从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换恢复外部空间计算点在球面上的投影点处的重力异常的全球影响;
步骤2、顾及实测数据局域保障条件,引入全球重力场位模型,利用移去-恢复技术,从重力异常观测值中移去参考重力异常,得到残差重力异常;从积分式核函数中移去与参考场对应阶次的核函数球谐表达式,得到截断核函数,使之与残差重力异常的频谱匹配;基于截断核函数和残差重力异常的局域积分,得到残差外部扰动重力径向分量;利用全球重力场位模型高阶信息进行远区效应补偿;恢复参考外部扰动重力径向分量;
步骤3、利用积分恒等式转换,补偿由全球积分过渡到局域积分引起的模型误差,得到计算点处高精度的外部扰动重力径向分量。
而且,所述步骤1的具体实现方法为:
建立基于重力异常计算外部扰动重力径向分量δgr的积分式:
式中,Δg为球面上流动点处的已知观测重力异常;R为地球椭球平均半径;r为计算点地心向径;为计算点的纬度和经度;为流动点的纬度和经度;σ为单位球面;dσ为单位球面的面积元;ψ为计算点至流动点之间的球面角距;是计算点至流动点之间的空间距离;Fr(r,ψ)为积分核函数;
恢复积分核函数的零阶和一阶项影响,得到完整核函数Ftr(r,ψ):
采用移去-恢复技术,从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换,恢复外部空间计算点在球面上的投影点处的重力异常的影响,以削弱积分核函数奇异性的影响,得到外部扰动重力径向分量计算公式:
式中ΔgRp为外部空间计算点P在球面上的投影点PR处的重力异常。
而且,所述步骤2处理后得到的外部扰动重力径向分量计算公式为:
式中;δΔg为移去了N阶位模型参考场的流动点残差重力异常;δΔgRp为移去了N阶位模型参考场的外部空间计算点P在球面上的投影点PR处的残差重力异常;为扰动重力径向分量截断核函数;为扰动重力径向分量远区效应计算值,δgrefr为由N阶参考场位模型计算得到的参考扰动重力径向分量;
其中,Pn(cosψ)为n阶勒让德函数;
式中,GM为地球引力常数;L为用于补偿远区效应的高阶次重力位模型的最高阶数;N代表由位模型定义的参考场最高阶数;Tn代表地球扰动位n阶拉普拉斯面球谐函数;Qn(δgr)为扰动重力径向分量积分核截断系数;为完全规格化缔合勒让德函数;和为完全规格化地球位系数。
而且,所述步骤3处理后得到的外部扰动重力径向分量计算公式为:
本发明的优点和积极效果是:
本发明设计合理,其针对基于重力异常计算外部扰动重力径向分量的不确定问题,采用移去-恢复运算和积分恒等式转换技术,同时顾及实测数据保障条件,建立基于重力异常计算外部扰动重力径向分量的无奇异局域算法,解决了基于重力异常计算外部扰动重力径向分量的不确定问题,提高外部扰动重力径向分量的计算精度,可广泛用于物理大地测量领域。
具体实施方式
本发明的设计思想是:(1)联合采用移去-恢复运算和积分恒等转换技术,恢复积分核函数的零阶和一阶项影响;从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换,恢复外部空间计算点在球面上的投影点处的重力异常的全球影响,以消除积分核函数奇异性的影响;(2)顾及实测数据局域保障条件,引入全球重力场位模型,利用移去-恢复技术,从重力异常观测值中移去参考重力异常,以得到残差重力异常;从积分式核函数中移去与参考场对应阶次的核函数球谐表达式,以得到截断核函数,使之与残差重力异常的频谱匹配;基于截断核函数和残差重力异常的局域积分得到残差外部重力异常;利用全球重力场位模型高阶信息进行远区效应补偿,以削弱远区截断误差的影响;恢复参考外部扰动重力径向分量。(3)利用积分恒等式转换,补偿由全球积分过渡到局域积分引起的模型误差,得到计算点处高精度的外部扰动重力径向分量。
基于上述设计思想,本发明提出一种基于重力异常计算外部扰动重力径向分量的方法,包括以下步骤:
步骤1、建立基于重力异常计算外部扰动重力径向分量δgr的积分式,恢复积分核函数的零阶和一阶项影响,得到完整核函数Ftr(r,ψ);从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换,恢复外部空间计算点在球面上的投影点处的重力异常的全球影响,以消除积分核函数奇异性的影响。本步骤的具体实现方法为:
建立基于重力异常计算外部扰动重力径向分量δgr的积分式:
式中,Δg为球面上流动点处的已知观测重力异常;R为地球椭球平均半径;r为计算点地心向径;为计算点的纬度和经度;为流动点的纬度和经度;σ为单位球面;dσ为单位球面的面积元;ψ为计算点至流动点之间的球面角距;是计算点至流动点之间的空间距离;Fr(r,ψ)为积分核函数。
恢复积分核函数的零阶和一阶项影响,得到完整核函数Ftr(r,ψ):
在式(3)中,当计算点趋近于数据点时,即当r→R和ψ→0时,会出现分母项l→0,积分核函数Ftr(r,ψ)发生奇异,这说明利用式(3)计算地球外部超低空扰动重力径向分量时,会出现由核函数奇异性引起的不确定性问题,无法保证计算结果的可靠性和有效性。为了消除式(1)的奇异性,确保外部扰动重力径向分量的整体连续性,采用移去-恢复技术,从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换,恢复外部空间计算点在球面上的投影点处的重力异常的影响,以削弱积分核函数奇异性的影响;可得到如下公式:
式中ΔgRp为外部空间计算点P在球面上的投影点PR处的重力异常。
将公式(4)可化为:
经上述移去恢复变换后,使用式(5)计算外部扰动重力径向分量,不仅可以避免奇异积分的影响,同时可确保积分计算值从地球外部到球边界面的连续性。
步骤2、由于计算外部扰动重力径向分量要求全球积分,但在实际应用中受观测数据覆盖范围限制,无法做到全球覆盖,因此,实际计算过程中外部扰动重力径向分量的全球积分式需要进行适用观测数据保障条件的局域化处理,以保证计算结果的可靠性。本发明顾及实测数据局域保障条件,引入全球重力场位模型,利用移去-恢复技术,从重力异常观测值中移去参考重力异常,以得到残差重力异常;从积分式核函数中移去与参考场对应阶次的核函数球谐表达式,以得到截断核函数,使之与残差重力异常的频谱匹配;基于截断核函数和残差重力异常的局域积分得到残差外部扰动重力径向分量;利用全球重力场位模型高阶信息进行远区效应补偿,以削弱远区截断误差的影响;恢复参考外部扰动重力径向分量。具体实现方法为:
将式(5)化为:
式中;δΔg为移去了N阶位模型参考场的流动点残差重力异常;δΔgRp为移去了N阶位模型参考场的流外部空间计算点P在球面上的投影点PR处的残差重力异常;为扰动重力径向分量截断核函数;为扰动重力径向分量远区效应计算值,δgrefr为由N阶参考场位模型计算得到的参考扰动重力径向分量;其它符号意义同前。
式中,Pn(cosψ)为n阶勒让德(Legendre)函数,其它符号意义同前。
式中,GM为地球引力常数;L为用于补偿远区效应的高阶次重力位模型的最高阶数;N代表由位模型定义的参考场最高阶数;Tn代表地球扰动位n阶拉普拉斯(Laplace)面球谐函数;Qn(δgr)为扰动重力径向分量积分核截断系数;为完全规格化缔合勒让德函数;和为完全规格化地球位系数。
步骤3、在既消除积分奇异性又顾及实测数据局域保障条件的式(6)中,由于远区效应只代表积分项δΔgq在远区的补偿,并未顾及另一积分项δΔgRp在远区的影响。对此利用积分恒等式转换,补偿由全球积分过渡到局域积分引起的模型误差,得到计算点处高精度的外部扰动重力径向分量。本步骤的具体实现方法为:
将式(6)转换为:
下面采用超高阶位模型EGM2008作为数值计算检验的参考标准场,用于模拟产生地球表面1′×1′网格重力异常观测量“真值”及地球外部不同高度面上的1′×1′网格扰动重力径向分量理论“真值”。为了体现检验结果的代表性,这里特意选取重力异常场变化比较剧烈的马里亚纳海沟作为试验区,具体覆盖范围为:6°×6°(10°N~16°N;λ:142°E~148°E)。选取r=R+h,R=6371km,使用标准场模型EGM2008分别计算对应于9个高度面上的1′×1′网格扰动重力径向分量理论“真值”δgtri(i=1,2,…,9),每个高度面对应360×360=129600个网格点数据,9个高度分别取为:hi=0km,0.1km,0.3km,1km,3km,5km,10km,30km,50km。表1列出了其中的4个高度面上的扰动重力径向分量理论“真值”和球面上的重力异常观测量“真值”的统计结果。
表1由EGM2008模型计算得到的重力异常和扰动重力统计结果(单位:mGal)
为了对比分析本发明算法的计算效果,采用球面上的1′×1′网格重力异常“真值”Δgt作为观测量,基于本发明算法对前面选定的试验区对应于9个高度面上的1′×1′网格扰动重力径向分量进行计算分析;其中,传统算法是指直接使用式(1)作为基础计算模型,并对全球积分域作了分区处理,但在实施近区计算时,扣除掉计算点所在的1′×1′数据块,以避免出现奇异积分问题。将计算值分别与相对应的理论“真值”Δgtri作比较,可获得改精度评估信息,具体比对结果列于表2。这里积分半径统一取为ψ0=2°,为了减小积分边缘效应对评估结果的影响,表2只列出中心区2°×2°方块内的比对结果。
表2由本发明计算得到的9个高度面扰动重力径向分量与“真值”的比较(单位:mGal)
对比表2和表1结果可以看出,传统算法在超低空高度段的误差量值甚至超过了径向分量自身大小,说明传统算法在超低空高度段是失效的,只有在5km以上计算高度才是可用的;本发明算法从理论上消除了积分奇异性和数值不连续性影响,并在超低空高度段取得了比传统模型好得多的计算精度,在所有9个高度面,计算值与比对基准“真值”的最大互差均不超过1.2mGal,均方根值不超过0.4mGal,验证了本发明算法的必要性和有效性,表明新方法实用易行,具有较高的应用价值。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (4)
1.一种基于重力异常计算外部扰动重力径向分量的方法,其特征在于:包括以下步骤:
步骤1、建立基于重力异常计算外部扰动重力径向分量δgr的积分式,恢复积分核函数的零阶和一阶项影响,得到完整核函数Ftr(r,ψ);从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换恢复外部空间计算点在球面上的投影点处的重力异常的全球影响;
步骤2、顾及实测数据局域保障条件,引入全球重力场位模型,利用移去-恢复技术,从重力异常观测值中移去参考重力异常,得到残差重力异常;从积分式核函数中移去与参考场对应阶次的核函数球谐表达式,得到截断核函数,使之与残差重力异常的频谱匹配;基于截断核函数和残差重力异常的局域积分,得到残差外部扰动重力径向分量;利用全球重力场位模型高阶信息进行远区效应补偿;恢复参考外部扰动重力径向分量;
步骤3、利用积分恒等式转换,补偿由全球积分过渡到局域积分引起的模型误差,得到计算点处高精度的外部扰动重力径向分量。
2.根据权利要求1所述的一种基于重力异常计算外部扰动重力径向分量的方法,其特征在于:所述步骤1的具体实现方法为:
建立基于重力异常计算外部扰动重力径向分量δgr的积分式:
式中,Δg为球面上流动点处的已知观测重力异常;R为地球椭球平均半径;r为计算点地心向径;σ为单位球面;dσ为单位球面的面积元;ψ为计算点至流动点之间的球面角距;是计算点至流动点之间的空间距离;Fr(r,ψ)为积分核函数;
恢复积分核函数的零阶和一阶项影响,得到完整核函数Ftr(r,ψ):
采用移去-恢复技术,从流动点重力异常中移去外部空间计算点在球面上的投影点处的重力异常,通过积分恒等变换,恢复外部空间计算点在球面上的投影点处的重力异常的影响,以削弱积分核函数奇异性的影响,得到外部扰动重力径向分量计算公式:
式中ΔgRp为外部空间计算点P在球面上的投影点PR处的重力异常。
3.根据权利要求2所述的一种基于重力异常计算外部扰动重力径向分量的方法,其特征在于:所述步骤2处理后得到的外部扰动重力径向分量计算公式为:
式中;δΔg为移去了N阶位模型参考场的流动点残差重力异常;δΔgRp为移去了N阶位模型参考场的外部空间计算点P在球面上的投影点PR处的残差重力异常;为扰动重力径向分量截断核函数;为扰动重力径向分量远区效应计算值,δgrefr为由N阶参考场位模型计算得到的参考扰动重力径向分量;
其中,Pn(cosψ)为n阶勒让德函数;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110180191.4A CN112965127B (zh) | 2021-02-08 | 2021-02-08 | 一种基于重力异常计算外部扰动重力径向分量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110180191.4A CN112965127B (zh) | 2021-02-08 | 2021-02-08 | 一种基于重力异常计算外部扰动重力径向分量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112965127A CN112965127A (zh) | 2021-06-15 |
CN112965127B true CN112965127B (zh) | 2022-05-31 |
Family
ID=76284623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110180191.4A Active CN112965127B (zh) | 2021-02-08 | 2021-02-08 | 一种基于重力异常计算外部扰动重力径向分量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112965127B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115236759B (zh) * | 2022-02-28 | 2023-09-05 | 中国人民解放军战略支援部队信息工程大学 | 一种确定地球重力场的六边形网格剖分方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7113868B2 (en) * | 2004-08-11 | 2006-09-26 | Bell Geospace, Inc. | Method and system for processing geophysical survey data |
GB201008993D0 (en) * | 2010-05-28 | 2010-07-14 | Arkex Ltd | Processing geophysical data |
CN108415879B (zh) * | 2018-01-19 | 2021-04-06 | 中国人民解放军92859部队 | 基于向上延拓的航空重力最小二乘向下延拓解析方法 |
CN108267792B (zh) * | 2018-04-13 | 2019-07-12 | 武汉大学 | 全球重力场模型反演方法 |
CN108873091B (zh) * | 2018-06-29 | 2019-10-25 | 中国人民解放军61540部队 | 卫星重力梯度全张量恢复地球重力场的确定方法及系统 |
CN109856691B (zh) * | 2019-03-14 | 2020-08-14 | 中国人民解放军61540部队 | 一种基于梯度法的航空重力矢量向下延拓方法及系统 |
-
2021
- 2021-02-08 CN CN202110180191.4A patent/CN112965127B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN112965127A (zh) | 2021-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110487301B (zh) | 一种雷达辅助机载捷联惯性导航系统初始对准方法 | |
CN110631574B (zh) | 一种惯性/里程计/rtk多信息融合方法 | |
CN109425339B (zh) | 一种基于惯性技术的考虑杆臂效应的舰船升沉误差补偿方法 | |
CN109507706B (zh) | 一种gps信号丢失的预测定位方法 | |
CN116182795B (zh) | 普速铁路纵断面精密测量方法 | |
CN112965127B (zh) | 一种基于重力异常计算外部扰动重力径向分量的方法 | |
CN111722295B (zh) | 一种水下捷联式重力测量数据处理方法 | |
CN112965124B (zh) | 一种顾及局域保障条件计算外部重力异常垂直梯度的方法 | |
CN110631573B (zh) | 一种惯性/里程计/全站仪多信息融合方法 | |
CN110058324B (zh) | 利用重力场模型的捷联式重力仪水平分量误差修正方法 | |
CN113900069A (zh) | 一种基于干涉成像高度计的垂线偏差计算方法及其系统 | |
CN113008229A (zh) | 一种基于低成本车载传感器的分布式自主组合导航方法 | |
CN112949049B (zh) | 一种利用带限思想计算重力异常低阶径向导数的方法 | |
CN112965125B (zh) | 一种基于重力异常计算外部扰动重力东向分量的方法 | |
CN112965123B (zh) | 一种基于重力异常计算外部扰动重力北向分量的方法 | |
CN112965128B (zh) | 一种无奇异性顾及局域保障条件计算外部重力异常的方法 | |
CN115790645A (zh) | 一种车载组合导航系统轮速计误差在线估计和补偿方法 | |
CN112987118B (zh) | 一种利用带限思想计算重力异常高阶径向导数的方法 | |
CN112818285A (zh) | 一种计算外部扰动重力北向分量中央区效应的方法 | |
CN113124863B (zh) | 一种基于kld采样的混合粒子联邦滤波数据处理方法 | |
CN113985491B (zh) | 一种基于多源数据的重力场模型精化方法及系统 | |
CN112965126B (zh) | 一种计算外部扰动重力东向分量中央区效应的方法 | |
CN117948986B (zh) | 一种极区因子图构建方法和极区组合导航方法 | |
CN114543799B (zh) | 一种抗差联邦卡尔曼滤波方法、设备与系统 | |
CN111854728B (zh) | 一种基于广义相对熵的容错滤波方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |