CN108846195A - 一种极区扰动重力无奇异性详细计算模型及其建模方法 - Google Patents
一种极区扰动重力无奇异性详细计算模型及其建模方法 Download PDFInfo
- Publication number
- CN108846195A CN108846195A CN201810591809.4A CN201810591809A CN108846195A CN 108846195 A CN108846195 A CN 108846195A CN 201810591809 A CN201810591809 A CN 201810591809A CN 108846195 A CN108846195 A CN 108846195A
- Authority
- CN
- China
- Prior art keywords
- formula
- singularity
- disturbing gravity
- gravity
- disturbing
- 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.)
- Pending
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/06—Power analysis or power optimisation
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
本发明提供一种极区扰动重力无θ奇异性的详细计算模型,即从Legendre函数的基本性质出发,研究得出Legendre函数对θ的一阶导数和一类多项式的无奇异性计算公式,再将其代入扰动重力的传统计算模型中,同时充分考虑了m等于0、1以及其它量时的情况,最终建立了极区扰动重力无θ奇异性的详细计算模型。该发明使得扰动重力计算模型更接近于研究对象的物理特性,并且公式形式简单,稳定性和实用性较好,彻底解决了两极地区扰动重力计算公式存在的θ奇异性难题,为飞行器飞越两极及其附近地区时的扰动重力计算奠定了理论基础,也可利用该发明提供的模型计算全球任何区域的扰动重力结果。
Description
技术领域:
本发明公开了一种极区扰动重力无奇异性详细计算模型及其建模方法,属于大地测量领域
背景技术
在地球物理问题的求解中,为了使数学模型更接近于研究对象的物理特性,坐标形式的选择往往取决于所研究对象边界面的几何形状。因此,在扰动重力的表示中,更倾向于采用地心球坐标(r,θ,λ),这使得扰动重力的计算变得更加简单实用。
扰动重力三分量的地心球坐标表达式为
其中,δρ、δλ分别表示扰动重力的径向、纬向和经向分量;fM表示地球重力常数;ρ表示计算点的地心向径;θ和λ分别表示计算点的余纬和地心经度;R表示地球平均半径;表示完全正常化地球扰动重力位系数;n和m分别表示球谐系数的阶和次;表示完全正常化缔合Legendre函数。
公式(2)中的的递推计算公式为
式中是kroneker符号,当m=0时,当m≠0时,
将(4)式代入(2)式,就可以得到扰动重力的传统计算模型。
从公式(2)-(4)可以看出,当飞行器飞越两极及其附近地区时,由于分母上的sinθ=0或者接近于0,会使得扰动重力的2个分量的计算出现无穷大。然而,地球两极及其附近地区作为地球表面上的普通区域,扰动重力的计算值应该与其它地球表面一样是有限值,不同区域,量级不会太大的差距。因此,这就是极区扰动重力计算存在的θ奇异性问题。随着重力测量技术的发展,从理论上彻底解决球坐标下极区扰动重力计算存在的θ奇异性问题就显得越来越重要。
为了解决极区扰动重力计算存在的θ奇异性问题,传统方法主要是直角坐标表示法。
直角坐标表示法的基本原理是指将扰动重力直接表示成地心直角坐标的形式,因而分母上不会含有等于0或接近于0的项,避免了扰动重力计算存在的θ奇异性问题。但也由此引入了其他的问题:公式形式和参数比较复杂,对递推计算的稳定性和精度也有较大影响。故本专利提出一种极区扰动重力无θ奇异性的详细计算模型。
发明内容
本发明提供一种极区扰动重力无θ奇异性的详细计算模型,即从Legendre函数的基本性质出发,研究得出Legendre函数对θ的一阶导数和一类多项式的无奇异性计算公式,再将其代入扰动重力的传统计算模型中,同时充分考虑了m等于0、1以及其它量时的情况,最终建立了极区扰动重力无θ奇异性的详细计算模型。该发明使得扰动重力计算模型更接近于研究对象的物理特性,并且公式形式简单,稳定性和实用性较好,彻底解决了两极地区扰动重力计算公式存在的θ奇异性难题,为飞行器飞越两极及其附近地区时的扰动重力计算奠定了理论基础,也可利用该发明提供的模型计算全球任何区域的扰动重力结果。
参阅图7,本发明提出的一种极区扰动重力无奇异性详细计算模型的建立过程主要步骤如下:
第一步:计算得到Legendre函数对θ的一阶导数的无奇异性结果。
用表示伴随Legendre函数,其中Pn(x)是n阶Legendre多项式,将x用cosθ来代替,可得
(5)式两端对θ求导,得
因为,Legendre多项式Pn(x)满足下列关系
结合(6)式和(7)式,可得
对(8)式中的Legendre函数进行规格化处理,得到Legendre函数对θ的一阶导数的无奇异性计算公式为:
第二步:计算得到一类多项式的无奇异性结果。
Legendre函数满足下列伴随Legendre方程
由此可得
在(8)式两端继续对θ求导,得
将(8)式代入(13)式,可以得到Legendre函数二阶导数的计算公式
将(8)式和(14)式代入(12)式,可得
根据Legendre函数的标准前向行推法,有
因此,
将(17)式和(18)式代入(15)式,可得
对(20)式中的Legendre函数进行规格化处理,得到多项式无奇异性计算公式为:
第三步:将上述公式代入扰动重力各分量的传统计算模型中,并充分考虑m等于0、1以及其它量时的情况,建立极区扰动重力无θ奇异性详细计算模型。
将(9)、(10)、(20)式代入(1)-(3)式,可得到扰动重力各分量的无θ奇异性计算公式为:
公式(21)-(23)即为扰动重力在极区的无θ奇异性详细计算模型。
在对本发明新建模型有效性进行测试基础上,本发明与现有技术相比体现出如下优点:
a.彻底解决了极区扰动重力计算公式存在的θ奇异性难题;
b.新建立的扰动重力各分量无θ奇异性详细计算模型,公式形式简单,稳定性和实用性较好;
c.新建立的扰动重力无θ奇异性详细计算模型,为飞行器飞越两极及其附近地区时的扰动重力计算奠定了理论基础,也可利用该发明提供的模型计算全球任何区域的扰动重力结果。
附图说明:
图1是实施例中南极地区扰动重力分量的等值线图;
图2是实施例中南极地区扰动重力δλ分量的等值线图;
图3是实施例中南极地区扰动重力δρ分量的等值线图;
图4是实施例中北极地区扰动重力分量的等值线图;
图5是实施例中北极地区扰动重力δλ分量的等值线图;
图6是实施例中北极地区扰动重力δρ分量的等值线图;
图7是本发明提出的极区扰动重力无θ奇异性详细计算模型建立基本流程。
具体实施方式:
假设飞行器飞行在距离地面300km的被动段轨道上,利用360阶EGM2008地球重力场模型,计算了南、北极圈(66.5°)以内30′×30′格网大小的扰动重力格网点值。两极地区扰动重力各分量的统计结果如表1、2所示:
表1南极地区扰动重力三分量统计结果(mGal)
表2北极地区扰动重力三分量统计结果(mGal)
从表1、2可以看出,扰动重力分量计算结果的标准差最小,因此其结果的变化幅度较小,δλ分量次之,δρ分量变化幅度最大,这是由于地球重力场对扰动重力δρ分量影响最大;从各分量的平均值来看,δλ分量沿格林尼治子午线的对称性最好。
为了更加直观地分析扰动重力各分量的变化趋势,本文绘出了南、北极圈以内扰动重力三分量的等值线图。
图1~6反映了扰动重力各分量的变化趋势,可以看出,扰动重力各分量基本上沿格林尼治子午线对称分布,其中δλ分量的对称性最好,这与表1、2的分析结果一致;北极地区比南极地区的扰动重力各分量的对称性要好,这是因为北极地区陆地面积很少,地形变化比较平缓的原因造成的;δρ分量比δλ分量变化更加剧烈,这与前面表1和表2的分析结果一致;各分量在陆地部分比海洋部分的变化更加剧烈,这是因为该地区陆地上地势起伏较大,地形情况复杂,而该地区海洋地形变化较为平缓的原因造成的。
图1~3反映出扰动重力各分量的变化比较均匀,这是因为南极地区大部分范围都是陆地,而图4~6反映出扰动重力各分量的变化比较复杂,从北极圈到北极点变化程度由剧烈到平缓,这是因为北极地区经历了从陆地到海洋的变化过程。
Claims (2)
1.一种极区扰动重力无奇异性详细计算模型,其特征在于,所述模型为:
式中,δρ、δλ分别表示扰动重力的径向、纬向和经向分量;fM表示地球重力常数;ρ表示计算点的地心向径;θ和λ分别表示计算点的余纬和地心经度;R表示地球平均半径;表示完全正常化地球扰动重力位系数;n和m分别表示球谐系数的阶和次;表示完全正常化缔合Legendre函数。
2.一种如权利要求1所述的极区扰动重力无奇异性详细计算模型的建模方法,其特征在于,包括如下步骤:
第一步:计算得到Legendre函数对θ的一阶导数的无奇异性结果;
用表示伴随Legendre函数,其中Pn(x)是n阶Legendre多项式,将x用cosθ来代替,可得
(5)式两端对θ求导,得
因为,Legendre多项式Pn(x)满足下列关系
结合(6)式和(7)式,可得
对(8)式中的Legendre函数进行规格化处理,得到Legendre函数对θ的一阶导数的无奇异性计算公式为:
第二步:计算得到一类多项式的无奇异性结果;
Legendre函数满足下列伴随Legendre方程
由此可得
在(8)式两端继续对θ求导,得
将(8)式代入(13)式,可以得到Legendre函数二阶导数的计算公式
将(8)式和(14)式代入(12)式,可得
根据Legendre函数的标准前向行推法,有
因此,
将(17)式和(18)式代入(15)式,可得
对(20)式中的Legendre函数进行规格化处理,得到多项式无奇异性计算公式为:
第三步:将上述公式代入扰动重力各分量的传统计算模型中,并充分考虑m等于0、1以及其它量时的情况,建立极区扰动重力无θ奇异性详细计算模型;
将(9)、(10)、(20)式代入(1)-(3)式,可得到扰动重力各分量的无θ奇异性计算公式为:
公式(21)-(23)即为扰动重力在极区的无θ奇异性详细计算模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810591809.4A CN108846195A (zh) | 2018-06-11 | 2018-06-11 | 一种极区扰动重力无奇异性详细计算模型及其建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810591809.4A CN108846195A (zh) | 2018-06-11 | 2018-06-11 | 一种极区扰动重力无奇异性详细计算模型及其建模方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108846195A true CN108846195A (zh) | 2018-11-20 |
Family
ID=64210681
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810591809.4A Pending CN108846195A (zh) | 2018-06-11 | 2018-06-11 | 一种极区扰动重力无奇异性详细计算模型及其建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108846195A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113885101A (zh) * | 2021-09-28 | 2022-01-04 | 中国船舶重工集团公司第七0七研究所 | 一种基于椭球模型构建重力梯度基准图方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100036614A1 (en) * | 2005-12-12 | 2010-02-11 | Hannes Georges Zuecher | Locating oil or gas passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background noise, including optional differentiation of oil, locatinggas and water |
CN104834320A (zh) * | 2015-03-23 | 2015-08-12 | 中国人民解放军信息工程大学 | 一种空间分层扰动引力场网格模型快速构建方法 |
-
2018
- 2018-06-11 CN CN201810591809.4A patent/CN108846195A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100036614A1 (en) * | 2005-12-12 | 2010-02-11 | Hannes Georges Zuecher | Locating oil or gas passively by observing a porous oil and gas saturated system giving off its characteristic resonance response to ambient background noise, including optional differentiation of oil, locatinggas and water |
CN104834320A (zh) * | 2015-03-23 | 2015-08-12 | 中国人民解放军信息工程大学 | 一种空间分层扰动引力场网格模型快速构建方法 |
Non-Patent Citations (1)
Title |
---|
刘晓刚等: "弹道扰动引力无奇异性计算模型的建立", 《地球物理学进展》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113885101A (zh) * | 2021-09-28 | 2022-01-04 | 中国船舶重工集团公司第七0七研究所 | 一种基于椭球模型构建重力梯度基准图方法 |
CN113885101B (zh) * | 2021-09-28 | 2023-12-12 | 中国船舶重工集团公司第七0七研究所 | 一种基于椭球模型构建重力梯度基准图方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | A multiresolution hierarchical classification algorithm for filtering airborne LiDAR data | |
Slivinski et al. | A hybrid particle–ensemble Kalman filter for Lagrangian data assimilation | |
CN103675790A (zh) | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 | |
Redpath et al. | Accuracy assessment for mapping glacier flow velocity and detecting flow dynamics from ASTER satellite imagery: Tasman Glacier, New Zealand | |
Maclachlan et al. | Quantitative geomorphological analysis of drumlins in the Peterborough drumlin field, Ontario, Canada | |
CN107844757A (zh) | 一种利用数字高程模型提取流域内河道宽度的方法 | |
CN109214050A (zh) | 一种极区垂线偏差无奇异性详细计算模型及其建模方法 | |
Nayak et al. | An extended DBC approach by using maximum Euclidian distance for fractal dimension of color images | |
CN107561592B (zh) | 一种密度为多项式的任意多面体的重力场正演方法 | |
Patel et al. | Terrain characterization using SRTM data | |
CN103617336A (zh) | 一种航空噪声等值线图的绘制方法 | |
Fuss et al. | DEM Fusion using a modified k-means clustering algorithm | |
Favalli et al. | Digital elevation model construction from structured topographic data: the DEST algorithm | |
Uddin et al. | Two dimensional hydrodynamic modelling of Northern Bay of Bengal coastal waters | |
CN105573963A (zh) | 一种电离层水平不均匀结构重构方法 | |
CN108846195A (zh) | 一种极区扰动重力无奇异性详细计算模型及其建模方法 | |
Penna et al. | Statistical characterization of unworked and water-worked gravel-bed roughness structures | |
CN105678716B (zh) | 一种地基sar大气干扰相位校正方法及装置 | |
Ben-Artzi et al. | Hyperbolic conservation laws on the sphere. A geometry-compatible finite volume scheme | |
US20190078944A1 (en) | Method for estimating geothermal gradient and device for estimating geothermal gradient | |
Debella-Gilo | Bare-earth extraction and dtm generation from photogrammetric point clouds with a partial use of an existing lower resolution DTM | |
Iovino et al. | Evaluation of the CMCC global eddying ocean model for the Ocean Model Intercomparison Project (OMIP2) | |
Li et al. | Extraction and analysis of elevation changes in Antarctic ice sheet from CryoSat-2 and Sentinel-3 radar altimeters | |
Elshambaky | Enhancing the predictability of least-squares collocation through the integration with least-squares-support vector machine | |
JIANG et al. | A study on the correlation of ridge line and slope with peak ground velocity amplification factor |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181120 |