利用航空重力梯度数据恢复区域重力场

发布时间：2017-09-14 14:30

  本文关键词：利用航空重力梯度数据恢复区域重力场

  更多相关文章： 地球重力场 航空重力梯度测量 矩谐分析 向下延拓 数值模拟

【摘要】：地球重力场作为地球的基本物理场之一，制约着地球系统内部物质的运动和变化。同时，地球系统内部的所有物理事件均在地球重力场中有所反映。因此高精度的地球重力场信息的获取对于固体地球物理、海洋学等相关科学领域的研究有举足轻重的作用。随着空间技术的发展，以飞机作为观测载体探测地球重力场的航空重力测量技术得到了不断发展，进而地球重力场的观测精度也随着航空重力测量系统的发展不断提高。这使得获取更高分辨率更高精度的地球重力位模型成为可能。 从理论上讲，航空重力梯度测量（AGG:Airborne Gravity Gradiometry）所能获得的重力观测值的精度和分辨率比航空标量或矢量重力测量所得结果更高，在地球物理勘探等领域具有重要应用价值。因此，当前对航空重力梯度测量数据处理理论与方法的研究更具有现实意义。且对于区域重力场逼近而言，矩谐分析（RHA: Rectangular HarmonicAnalysis）是一种理想的实用方法。利用AGG数据进行矩谐分析，有效控制边缘效应，通过适当的正则化方法克服向下延拓等病态问题。 在此背景下，本文研究基于航空重力梯度测量数据恢复区域重力场的理论和方法，编写了AGG数据处理程序，开发了数值模拟软件，不仅为AGG实测数据处理和相关应用研究奠定基础，，也为发展我国自主的航空重力梯度测量系统以及利用AGG数据研制高精度的地球重力位模型积累经验。 本文的研究内容及创新点主要包括以下几个方面： 1.系统总结和归纳地球重力场与相关学科领域的关系，对地球重力场逼近理论的发展进行多层面的概述，概括介绍了航空重力梯度测量系统的发展历程，对AGG数据处理及恢复地球重力场的国内外研究进展进行总结和评述。 2.详细研究航空重力梯度测量所涉及基本理论，包括探究过程中用到的坐标系统、重力梯度张量在不同坐标间的转换关系、AGG测量误差及消弱措施等。此外，对本文数据处理中未涉及的数据预处理进行了概括性的阐述。 3.对航空重力梯度数据恢复区域重力场中所遇到的关键性技术进行了详细研究和分析。提出了航空重力梯度解算的矩谐分析法，导出了求解扰动位矩谐系数的表达式，利用重力梯度张量解算矩谐系数，构建重力位模型。重点研究了Tikhonov正则化、岭估计、TSVD正则化方法的基本原理和数学模型，并结合适当的正则化参数选取方法解决了航空重力梯度观测值在向下延拓过程中法方程的病态问题，为软件编制提供算法依据。且通过扩展计算区域的范围方法降低了边缘效应，数值计算表明：当扩展参数取值在[50,60]km之间，可显著提高矩谐分析的解算精度；正则化参数选取方法中，GCV获得的最优参数相对稳定。 4.因无实测航空重力梯度数据，为验证本文提出的恢复区域重力场模型的合理性，只能借助模拟数据来进行试算工作。本文构建了模拟航空重力梯度测量数据的数学模型，进行AGG观测数据仿真；利用仿真值结合AGG数学模型进行位系数计算，并计算模型扰动重力；最后由EGM2008地球重力位模型结合超高阶地球重力场参数计算软件，进行扰动重力真值计算。数值计算表明：基于TSVD的延拓值与真实值的差值在±20mGal以内，大部分计算点处的差值不超过±5mGal，TSVD正则化方法能够实现重力场信号的稳定向下延拓，本文提出的航空重力梯度数据恢复区域重力场解算方案可行。 5.独立编制了一套利用航空重力梯度数据恢复区域重力场的软件系统，包含有航空重力梯度测量基本理论（参考系统、坐标转换）、数据模拟、AGG矩谐分析模型、正则化方法解算及重力场参数赋值等多个模块。
【关键词】：地球重力场 航空重力梯度测量 矩谐分析 向下延拓 数值模拟
【学位授予单位】：兰州交通大学
【学位级别】：硕士
【学位授予年份】：2014
【分类号】：P223
【目录】：

• 摘要4-6
• Abstract6-11
• 1 绪论11-18
• 1.1 引言11-12
• 1.2 航空重力探测技术12-16
• 1.2.1 航空重力梯度测量系统的研究进展12-14
• 1.2.2 航空重力梯度测量数据处理的研究进展14-15
• 1.2.3 恢复区域重力场的研究进展15-16
• 1.3 本文研究目标与研究内容16-18
• 2 航空重力梯度测量基本理论与数据预处理18-30
• 2.1 引言18
• 2.2 航空重力梯度测量的基本理论18-21
• 2.2.1 航空重力梯度测量系统18-19
• 2.2.2 重力梯度张量19-21
• 2.3 坐标系统与坐标转换21-24
• 2.3.1 常用坐标系统21-23
• 2.3.2 重力梯度张量实用转换公式23-24
• 2.4 数据预处理方法24-27
• 2.4.1 飞机速度和加速的确定24-25
• 2.4.2 低通滤波25-26
• 2.4.3 归算至平均飞行高度面26-27
• 2.5 航空重力梯度测量误差分析27-28
• 2.5.1 主要误差源27-28
• 2.5.2 消弱措施28
• 2.6 本章小结28-30
• 3 利用航空重力梯度数据恢复区域地球重力场的理论和方法30-42
• 3.1 引言30
• 3.2 基于矩谐分析的航空重力梯度解算30-34
• 3.2.1 矩谐分析的基本理论30-31
• 3.2.2 重力梯度张量的矩谐展开31-33
• 3.2.3 最小二乘求解33
• 3.2.4 边缘效应33-34
• 3.3 病态性分析34
• 3.4 正则化解法及参数选取34-40
• 3.4.1 正则化法34-37
• 3.4.2 正则化参数选取37-40
• 3.5 数值计算与分析40-41
• 3.6 本章小结41-42
• 4 航空重力梯度观测数据的模拟与试算结果42-52
• 4.1 引言42
• 4.2 航空重力梯度观测数据的模拟42-44
• 4.2.1 扰动重力梯度张量数学模型的建立42-43
• 4.2.2 Legendre 函数及其导数的计算43-44
• 4.3 模拟试算结果与分析44-51
• 4.3.1 模拟观测数据44-46
• 4.3.2 恢复区域重力场信息46-51
• 4.4 本章小结51-52
• 5 结论与展望52-55
• 5.1 结论52-53
• 5.2 存在的问题与展望53-55
• 致谢55-56
• 参考文献56-59
• 攻读硕士学位期间的主要工作和成果59

【相似文献】

中国期刊全文数据库 前10条

1 贾明;杨功流;;基于超导的全张量重力梯度敏感头设计[J];中国惯性技术学报;2012年02期

2 张传定,吴晓平,陆仲连;全张量重力梯度数据的谱表示方法[J];测绘学报;2000年04期

3 吴星;;卫星重力梯度数据处理理论与方法[J];测绘学报;2010年02期

4 曹书锦;朱自强;鲁光银;曾思红;郭文斌;;基于单元细分H-自适应有限元全张量重力梯度正演[J];地球物理学进展;2010年03期

5 孙鹏飞;吴燕冈;;重力梯度法中台阶倾角的确定与实际应用[J];金属矿山;2011年11期

6 钱东;刘繁明;李艳;张迎发;;导航用重力梯度基准图构建方法的比较研究[J];测绘学报;2011年06期

7 孙鹏飞;吴燕冈;杨春成;韩兆红;范美宁;;重力梯度法的角点位置最优化选取研究(英文)[J];Applied Geophysics;2011年04期

8 吴星;王凯;冯炜;汪涛;;基于非全张量卫星重力梯度数据的张量不变量法[J];地球物理学报;2011年04期

9 王东明;利用地球重力位模型计算重力和重力梯度[J];地球物理学报;1999年S1期

10 魏兴;荆心;刘繁明;钱东;;利用重力梯度数据的高精度地形估计[J];东北林业大学学报;2013年12期

中国重要会议论文全文数据库 前10条

1 郑伟;许厚泽;钟敏;员美娟;周旭华;彭碧波;;我国将来卫星重力梯度计划的实施[A];中国地球物理·2009[C];2009年

2 刘利军;;聚焦反演应用于重力梯度独立分量[A];中国地球物理学会第二十七届年会论文集[C];2011年

3 吴星;张传定;;卫星重力梯度数据的模拟研究[A];中国地球物理学会第二十四届年会论文集[C];2008年

4 王虎彪;王勇;许大欣;戴全发;;基于重力梯度的水下辅助导航仿真[A];中国地球物理学会第二十四届年会论文集[C];2008年

5 李红蕾;方剑;;代数重构算法在重力梯度张量数据反演中的应用[A];中国地球物理2013——第二十六专题论文集[C];2013年

6 方明;赵德军;;重力梯度张量的延拓[A];第二十一届海洋测绘综合性学术研讨会论文集[C];2009年

7 王萌;郑元满;姚长利;;应用重力梯度张量方法定位场源体中心及推断其走向[A];中国地球物理2013——第十八专题论文集[C];2013年

8 吴怿昊;罗志才;钟波;吴云龙;;利用地球重力场模型校准卫星重力梯度数据[A];中国地球物理学会第二十七届年会论文集[C];2011年

9 吴星;李建伟;张传定;;卫星重力梯度分量复数组合的向上和向下延拓[A];中国地球物理·2009[C];2009年

10 许闯;罗志才;林旭;李振海;;GOCE卫星重力梯度观测数据的海潮影响[A];中国地球物理学会第二十七届年会论文集[C];2011年

中国博士学位论文全文数据库 前2条

1 吴星;卫星重力梯度数据处理理论与方法[D];解放军信息工程大学;2009年

2 刘凤鸣;重力梯度辅助导航定位技术研究[D];哈尔滨工程大学;2010年

中国硕士学位论文全文数据库 前10条

1 李德禹;水下重力梯度导航关键技术研究[D];华中科技大学;2009年

2 孙鹏飞;重力梯度法的应用及改进[D];吉林大学;2008年

3 聂慧梓;重力梯度张量三维聚焦反演方法研究[D];中国地质大学（北京）;2014年

4 陈少华;基于预条件共轭梯度法的重力梯度张量反演研究[D];中南大学;2012年

5 陈涛;三维重力梯度快速正反演研究[D];中国地质大学（北京）;2014年

6 李伟;利用航空重力梯度数据恢复区域重力场[D];兰州交通大学;2014年

7 夏朝龙;顾及地形、大气影响的GOCE卫星SGG数据到参考基准的归算研究[D];长安大学;2013年

8 王灿;重力梯度张量解析信号的欧拉反褶积[D];中南大学;2014年

9 郭文斌;重力梯度张量的拟BP神经网络反演[D];中南大学;2011年

10 孙亚力;基于GOCE卫星研究大气空间分布特性[D];长安大学;2012年

本文编号：850575