CN108319566B - 基于向上延拓的航空重力点对点向下延拓解析方法 - Google Patents
基于向上延拓的航空重力点对点向下延拓解析方法 Download PDFInfo
- Publication number
- CN108319566B CN108319566B CN201810053787.6A CN201810053787A CN108319566B CN 108319566 B CN108319566 B CN 108319566B CN 201810053787 A CN201810053787 A CN 201810053787A CN 108319566 B CN108319566 B CN 108319566B
- Authority
- CN
- China
- Prior art keywords
- point
- continuation
- gravity
- downward
- upward
- 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
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于向上延拓的航空重力点对点向下延拓解析方法,其主要技术特点是包括:依据泰勒级数展开思想,建立航空重力向上延拓与向下延拓之间的点对点解析延拓模型;利用Poisson积分方程进行航空重力稳向上延拓,基于点对点解析延拓模型进行航空重力向下延拓反的解算。本发明设计合理,实现了航空重力向下延拓反问题的稳定解算,有效地解决了向下延拓解算固有的不适定性问题,同时利用超高阶位模型EGM2008建立的模拟标准场数据对最小二乘向下延拓解析模型解算结果的合理性和有效性进行了数值验证,试验证明本发明实用易行,具有较高的应用价值。
Description
技术领域
本发明属于海洋重力技术领域,尤其是一种基于向上延拓的航空重力点对点向下延拓解析方法。
背景技术
航空重力位场向上和向下延拓是地球物理和大地测量数据分析应用中最重要的技术环节之一。航空重力位场向下延拓的目的是将数据观测面外推到更接近于位场源体,通过改善观测数据的信噪比,增强地球内部浅层质量分布异常的映射强度,来凸显地球重力场的局部变化特征,以提高位场数据解释推断的可靠性。但航空重力向下延拓在数学上属于典型的不适定反问题,其延拓算子对高频噪声具有明显的放大作用,很小的观测噪声也会引起延拓解算结果严重偏离真实的问题。
发明内容
本发明的目的在于克服现有技术的不足,提出一种基于向上延拓的航空重力点对点向下延拓解析方法,解决航空重力向下延拓的不适定问题。
本发明解决其技术问题是采取以下技术方案实现的:
一种基于向上延拓的航空重力点对点向下延拓解析方法,包括以下步骤:
步骤1、依据泰勒级数展开思想,建立航空重力向上延拓与向下延拓之间的点对点解析延拓模型;
步骤2、利用Poisson积分方程进行航空重力稳向上延拓,基于点对点解析延拓模型进行航空重力向下延拓反的解算;
所述步骤2采用如下Poisson积分方程进行航空重力稳向上延拓:
式中,rp=R+hp,rq=R+hq,R为地球椭球平均半径;Δgpq代表与计算点Q相对应的飞行高度面上的重力异常;Δgp为飞行高度面流动点的重力异常;为计算点与流动点之间的空间距离,ψ为计算点与流动点之间的球面角距。
所述步骤1的实现方法为:
设海拔高度为ho的待求地球表面O点的重力异常为Δgo,已知对应飞行高度面hp上的空中重力异常为Δgp,则由泰勒级数展开模型可得到Δgo和Δgp之间的关系为:
设位于飞行高度面上方、海拔高度为hq处的空中重力异常为Δgq,令Δhqp=hq-hp,则根据重力场解析延拓理论,将Δgq和Δgp之间的关系表示为:
式中,δΔgpq代表重力异常Δgp到Δgq的向上延拓改正数;
取Δhqp=Δhpo,将上述两公式相加并处理得到:
其中,R2n(Δhpo)为向下延拓改正数的余项,设:
R2n(Δhpo)=0
则得到如下点对点解析延拓模型:
Δgo=2Δgp-Δgq
式中,航空重力异常Δgp是已知的,Δgq可依据Poisson积分方程由Δgp向上延拓计算得到。
本发明的优点和积极效果是:
1、本发明依据泰勒级数展开模型,建立航空重力向上延拓与向下延拓之间的点对点解析延拓模型;利用Poisson积分方程进行航空重力稳向上延拓,基于点对点解析延拓模型实现航空重力向下延拓反问题的稳定解算,有效地解决了向下延拓解算固有的不适定性问题。
2、本发明利用超高阶位模型EGM2008建立的模拟标准场数据对最小二乘向下延拓解析模型解算结果的合理性和有效性进行了数值验证,试验证明本发明实用易行,具有较高的应用价值。
具体实施方式
以下通过实施例对本发明做进一步详述。
一种基于向上延拓的航空重力点对点向下延拓解析方法,包括以下步骤:
步骤1、依据泰勒级数展开思想,建立航空重力向上延拓与向下延拓之间的点对点解析延拓模型。
在本步骤中,假设海拔高度为ho的待求地球表面O点的重力异常为Δgo,已知对应飞行高度面hp(P点)上的空中重力异常为Δgp,则由泰勒级数展开模型可知,Δgo和Δgp之间的关系可表示为:
假设位于飞行高度面上方、海拔高度为hq(Q点)处的空中重力异常为Δgq,令Δhqp=hq-hp,那么,根据重力场解析延拓理论,同样可以将Δgq和Δgp之间的关系表示为:
式中,δΔgpq代表重力异常Δgp到Δgq的向上延拓改正数,其它符号意义同前。取Δhqp=Δhpo,那么将式(1)和(2)两式相加并略加整理得:
这里称R2n(Δhpo)为向下延拓改正数的余项,由于垂向梯度随着阶次的增大而急剧减小,故可忽略该余项的影响,即令:
R2n(Δhpo)=0 (5)
那么就得到如下近似的点对点向下解析延拓计算式(点对点解析延拓模型):
Δgo=2Δgp-Δgq (6)
式中,航空重力异常Δgp是已知的,Δgq可依据Poisson积分方程由Δgp向上延拓计算得到。
步骤2、利用Poisson积分方程进行航空重力稳向上延拓,基于点对点解析延拓模型进行航空重力向下延拓反的求解。
本步骤采用的Poisson积分方程为:
式中,rp=R+hp,rq=R+hq,R为地球椭球平均半径;Δgpq代表与计算点Q相对应的飞行高度面上的重力异常;Δgp为飞行高度面流动点的重力异常;为计算点与流动点之间的空间距离,ψ为计算点与流动点之间的球面角距。使用式(7)进行向上延拓计算可避开常见的积分奇异性问题。
下面采用超高阶位模型作为标准场开展数值计算检验及分析比较研究。以美国本土西部一个3°×3°区块(37°N~40°N;λ:248°E~251°E)作为试验区,选用EGM2008位模型模拟产生不同高度面2′×2′网格重力异常及其垂向偏导数的“真值”。选择美国本土作为试验区,是考虑到EGM2008位模型在美国地区具有更高的逼近度。该区块属于地形变化比较剧烈的大山区,试验效果具有一定的代表性。由EGM2008位模型计算空中网格重力异常的公式为:
这里取R=6371km,hi=ikm(i=0,1,…,10),使用EGM2008模型(以360阶次作为参考场)分别计算11个对应于ri=R+hi球面上的2′×2′网格重力异常“真值”Δgi(tru),同时计算5km高度面(r5=R+h5)上的垂向偏导数。表1列出了4个不同高度面上的EGM2008位模型(361~2160阶次)残差重力异常的统计结果,表2列出了5km高度面垂向偏导数计算值的统计结果。本发明采用的试验方案具体设计为:以5km高度面作为航空重力测量的观测面,即假设Δg5(tru)是已知的观测量,依次采用不同的向下延拓模型,由Δg5(tru)分别计算4km、3km、2km、1km和0km高度面上的重力异常将计算值与相对应的“真值”Δgi(tru)(i=0,…,4)作比较,从而获得对应延拓模型解算精度的评估参数。
表1 不同高度面位模型残差重力异常统计结果(单位:10-5m﹒s-2)
高度(km) | 最小值 | 最大值 | 平均值 | 均方根 |
0 | -107.09 | 189.19 | -0.44 | 34.28 |
1 | -88.22 | 162.05 | -0.39 | 30.26 |
3 | -68.34 | 120.33 | -0.33 | 24.07 |
5 | -55.22 | 90.66 | -0.28 | 19.56 |
表2 5km高度面垂向偏导数统计结果
偏导数 | 最小值 | 最大值 | 平均值 | 均方根 |
一阶(mGal/km) | -3.45985 | 2.09489 | 0.00589 | 0.94238 |
二阶(mGal/km<sup>2</sup>) | -0.28745 | 0.49396 | 0.00049 | 0.09966 |
三阶(mGal/km<sup>3</sup>) | -0.08566 | 0.04196 | -0.00022 | 0.01449 |
四阶(mGal/km<sup>4</sup>) | -0.01119 | 0.01833 | 0.00006 | 0.00287 |
为了检验点对点延拓模型解式(6)的计算效果,首先使用向上延拓积分式(7),由Δg5(tru)依据公式(8)分别计算6km、7km、8km、9km和10km高度面上的重力异常求计算值与相对应“真值”Δgi(tru)(i=6,…,10)的互差,可得到向上延拓计算模型的精度评价,具体结果见表3。其中,积分半径统一取为ψ0=30′。为了减小积分边缘效应对评估结果的影响,计算区域边缘30′范围内的数据不参加对比分析,这个区域的数据也不再参加下一步的向下延拓计算检验。
表3 向上延拓模型计算精度检核(单位:10-5m˙s-2)
由表3结果看出,当观测数据不存在噪声干扰时,向上延拓模型计算结果可以达到很高的精度水平,即使是1km的延拓高度差,计算误差也不超过0.3mGal。分别将前面计算得到的6个高度面的向上延拓重力异常和飞行高度面上的观测重力异常Δg5(tru)代入向下延拓计算式(6),可求得相对应延拓高度差下的向下延拓重力异常:
表4 点对点向下延拓模型计算精度检核(单位:10-5m˙s-2)
高度差(km) | 最小值 | 最大值 | 平均值 | 均方根 |
1 | -0.34 | 0.28 | -0.01 | 0.10 |
2 | -1.86 | 1.01 | -0.01 | 0.36 |
3 | -4.45 | 2.50 | -0.01 | 0.88 |
4 | -8.19 | 4.61 | -0.01 | 1.61 |
5 | -13.21 | 7.40 | -0.02 | 2.59 |
表4中延拓高度差为ikm对应于计算面的高度为(5-i)km。由表4结果看出,虽然向上延拓模型的计算误差随计算高度的增大而减小,但点对点向下延拓模型的计算精度随延拓高度差的增大而降低,说明向下延拓模型的截断误差对延拓计算结果有较大影响。当延拓高度差小于3km时,该模型的计算精度优于1mGal。
为了进一步考察数据观测噪声对向下延拓计算结果的影响,这里在模拟观测量Δg5(tru)中分别加入±1mGal、±3mGal和±5mGal的白噪声,对应生成三组带噪声的观测量,并重复前面的计算和比对检核过程。加入观测噪声后,与表3相对应的向上延拓模型计算精度检核结果(这里只列出互差均方根值(rms),下同)如表5所示,与表4相对应的向下延拓模型计算精度检核结果如表6所示。
表5 噪声影响下的向上延拓误差均方根值(单位:10-5m˙s-2)
高度(km) | 6 | 7 | 8 | 9 | 10 |
噪声=1mGal | 0.26 | 0.21 | 0.17 | 0.14 | 0.12 |
噪声=3mGal | 0.70 | 0.56 | 0.45 | 0.37 | 0.31 |
噪声=5mGal | 1.18 | 0.95 | 0.78 | 0.65 | 0.54 |
表6 噪声影响下的点对点向下延拓误差均方根值(单位:10-5m˙s-2)
高度差(km) | 1 | 2 | 3 | 4 | 5 |
噪声=1mGal | 0.25 | 0.41 | 0.89 | 1.61 | 2.58 |
噪声=3mGal | 0.70 | 0.67 | 1.00 | 1.67 | 2.62 |
噪声=5mGal | 1.16 | 0.97 | 1.10 | 1.66 | 2.58 |
对比表3和表5数值结果可以看出,数据观测噪声对向上延拓计算结果的影响很小,几乎不改变计算参数的变化形态,说明向上延拓模型作为低通滤波器的作用相当明显。进一步对比表4和表6结果可以看出,得益于向上延拓模型抑制高频噪声的优良特性,点对点向下延拓模型计算结果也几乎不受数据观测噪声的影响,在延拓高度差小于3km的高度内,仍可获得1mGal左右的计算精度,证明本发明给出的技术方案实用易行,具有较高的应用价值。
需要强调的是,本发明所述的实施例是说明性的,而不是限定性的,因此本发明包括并不限于具体实施方式中所述的实施例,凡是由本领域技术人员根据本发明的技术方案得出的其他实施方式,同样属于本发明保护的范围。
Claims (2)
2.根据权利要求1所述基于向上延拓的航空重力点对点向下延拓解析方法,其特征在于:所述步骤1的实现方法为:
设海拔高度为ho的待求地球表面O点的重力异常为Δgo,已知对应飞行高度面hp上的空中重力异常为Δgp,则由泰勒级数展开模型可得到Δgo和Δgp之间的关系为:
设位于飞行高度面上方、海拔高度为hq处的空中重力异常为Δgq,令Δhqp=hq-hp,则根据重力场解析延拓理论,将Δgq和Δgp之间的关系表示为:
式中,δΔgpq代表重力异常Δgp到Δgq的向上延拓改正数;
取Δhqp=Δhpo,将上述两公式相加并处理得到:
其中,R2n(Δhpo)为向下延拓改正数的余项,设:
R2n(Δhpo)=0
则得到如下点对点解析延拓模型:
Δgo=2Δgp-Δgq
式中,航空重力异常Δgp是已知的,Δgq可依据Poisson积分方程由Δgp向上延拓计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810053787.6A CN108319566B (zh) | 2018-01-19 | 2018-01-19 | 基于向上延拓的航空重力点对点向下延拓解析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810053787.6A CN108319566B (zh) | 2018-01-19 | 2018-01-19 | 基于向上延拓的航空重力点对点向下延拓解析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108319566A CN108319566A (zh) | 2018-07-24 |
CN108319566B true CN108319566B (zh) | 2021-03-16 |
Family
ID=62895145
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810053787.6A Active CN108319566B (zh) | 2018-01-19 | 2018-01-19 | 基于向上延拓的航空重力点对点向下延拓解析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108319566B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109856691B (zh) * | 2019-03-14 | 2020-08-14 | 中国人民解放军61540部队 | 一种基于梯度法的航空重力矢量向下延拓方法及系统 |
CN110133749B (zh) * | 2019-05-30 | 2021-06-22 | 中国地质科学院 | 一种地质资源勘探中重磁数据处理方法及其系统 |
CN112949049B (zh) * | 2021-02-08 | 2021-11-30 | 中国人民解放军92859部队 | 一种利用带限思想计算重力异常低阶径向导数的方法 |
CN112836378B (zh) * | 2021-02-08 | 2023-01-10 | 中国人民解放军92859部队 | 基于Poisson理论计算外部重力异常垂直梯度中央区效应的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258128A (zh) * | 2013-05-08 | 2013-08-21 | 浙江大学 | 一种地磁场空间延拓算法的评估方法 |
CN103399350A (zh) * | 2013-07-29 | 2013-11-20 | 中国人民解放军国防科学技术大学 | 一种基于积分迭代算法的航空重力向下延拓方法 |
CN107291659A (zh) * | 2017-05-16 | 2017-10-24 | 哈尔滨工程大学 | 平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6993433B2 (en) * | 1999-04-02 | 2006-01-31 | Conocophillips Company | Modeling gravity and tensor gravity data using poisson's equation for airborne, surface and borehole applications |
WO2008060900A2 (en) * | 2006-11-10 | 2008-05-22 | Cqgt, Llc | Systems and methods for dynamic display of orders |
-
2018
- 2018-01-19 CN CN201810053787.6A patent/CN108319566B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258128A (zh) * | 2013-05-08 | 2013-08-21 | 浙江大学 | 一种地磁场空间延拓算法的评估方法 |
CN103399350A (zh) * | 2013-07-29 | 2013-11-20 | 中国人民解放军国防科学技术大学 | 一种基于积分迭代算法的航空重力向下延拓方法 |
CN107291659A (zh) * | 2017-05-16 | 2017-10-24 | 哈尔滨工程大学 | 平面地磁异常场一步向上延拓平面模量梯度场的递归余弦变换法 |
Non-Patent Citations (3)
Title |
---|
"利用改进泰勒级数法对磁异常场进行向下延拓";翟国君 等;《海洋测绘》;20120715;第32卷(第4期);第17-21页 * |
"顾及地形效应的地面重力向上延拓模型分析与检验";刘敏 等;《武汉大学学报•信息科学版》;20180105;第43卷(第1期);第113-119页 * |
Boyang Zhou." Far-zone contributions of airborne gravity anomalies" upward/downward continuation".《Geodesy and Geodynamics》.2016,第7卷(第6期), * |
Also Published As
Publication number | Publication date |
---|---|
CN108319566A (zh) | 2018-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108415879B (zh) | 基于向上延拓的航空重力最小二乘向下延拓解析方法 | |
CN108319566B (zh) | 基于向上延拓的航空重力点对点向下延拓解析方法 | |
CN103399350B (zh) | 一种基于积分迭代算法的航空重力向下延拓方法 | |
Zeng et al. | An adaptive iterative method for downward continuation of potential-field data from a horizontal plane | |
Hu et al. | Three-dimensional surface displacements from InSAR and GPS measurements with variance component estimation | |
CN110389391B (zh) | 一种基于空间域的重磁位场解析延拓方法 | |
CN109214050B (zh) | 一种极区垂线偏差无奇异性详细计算模型及其建模方法 | |
CN108680966A (zh) | 海洋可控源电磁勘探噪声降噪效果评估方法 | |
CN111257956A (zh) | 一种基于Matlab的区域似大地水准面精化方法 | |
CN111580174B (zh) | 一种基于帕德近似的重磁数据向下延拓方法 | |
CN113686329A (zh) | 一种基于地磁数据的垂直高度位场测量方法 | |
CN110989002A (zh) | 地空时域电磁系统浅部低阻异常体高精度数据解释方法 | |
JP2017138124A (ja) | 重力偏差データの表層密度値推定方法 | |
CN107942374A (zh) | 绕射波场提取方法和装置 | |
CN110161582B (zh) | 空中与地面数据结合的重力换算方法及系统 | |
CN111427096A (zh) | 全张量重力梯度仪数据质量评价和滤波处理方法 | |
CN113985490A (zh) | 利用地形和地壳密度数据进行地表重力仿真的方法及装置 | |
CN113960635A (zh) | 一种顾及日变化的对流层延迟改正方法 | |
CN108492204B (zh) | 一种高精度空地井不同空间重磁数据变换方法 | |
Pahlevi et al. | Determination of gravimetric geoid model in Sulawesi–Indonesia | |
CN112949049B (zh) | 一种利用带限思想计算重力异常低阶径向导数的方法 | |
CN117128921A (zh) | 一种陆地区域高程异常的确定方法 | |
CN108957571B (zh) | 一种航空重力数据插值、扩边和下延一体化方法 | |
CN113761457B (zh) | 基于测量的重力异常数据提取局部重力异常的方法 | |
CN106970426B (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 |