论文部分内容阅读
大地测量学的基本科学任务之一是精确确定地球重力场及其时变。高精度高分辨率的重力场模型不仅是地球内部物质分布、变化的反映,也是人们认知、利用和改造地球物理的基础。径向基函数是一种局部重力场建模方法,在空间域和频率域都有较好的局部化特征,可以以多种观测数据作为输入数据,广泛地应用于地球重力场建模领域。现阶段径向基函数应用过程中仍有很多问题未解决,如径向基函数甄选、模型参数求解、实测数据合理化使用、求逆过程中的不适定性以及扰动引力赋值精度和效率不能兼顾等。本文从径向基函数模型基础理论出发,针对基函数的建模流程以及基函数在局部重力场建模中的应用展开讨论,主要工作与结论如下:1.针对径向基函数利用离散数据建模时设计矩阵的不适定问题,在建模过程中引入Tikhonov正则化方法,消除了离散数据分布不均导致的设计矩阵的病态性,使得直接利用离散数据建模成为可能。实验结果表明,在实验区内当建模数据无污染时,正则化方法和格网化方法精度相当;在3mGal的噪声污染条件下,正则化方法较格网化方法精度提升27.9%。说明正则化方法可以提高模型的稳定性和抗干扰性。2.提出了一种可用于局部(似)大地水准面拟合的位置固定、埋深自由的半约束点质量核模型。利用空间距离最近点最优原则建立了解算模型参数的迭代算法,并顾及周围点影响提出固定加权点方法和自适应法。实验结果表明相比于CoKriging方法,实验区内的精度分别提升13.7%~72.6%。相比于格网化基函数模型,本文所提算法有系统偏差小、无奇异性、所需数据少以及计算效率高等优点。实验表明,利用本文所提方法进行(似)大地水准面拟合是可行的,自适应加权半约束点质量模型效能最优。自适应方法在保持精度的前提下,实验区内相比于简单半约束模型计算效率提升了 1.42~18.83倍。3.讨论了径向基函数多尺度建模的可行性。以传统经验模型参数为例给出了多尺度点质量核基函数的建模方法和各层模型参数,利用数值实验分析了多尺度模型的计算效能。实验证明多尺度模型精度较单尺度模型高,在空中扰动引力赋值应用时,三层模型可以兼顾计算效率和精度,三层组合模型在3Km高度层精度能达到0.5mGal。4.提出点质量核径向基函数大范围解算快速算法。引入窗口移动控制法将模型的设计矩阵转化为稀疏矩阵,并结合稀疏矩阵的压缩存储方法极大地降低了数据存储量。最后利用MKL数学库中的稀疏矩阵计算库提升了方程求逆的效率,扩大了模型的解算范围。实验表明在同一台计算设备上模型解算半径扩展了约4.2倍。在基础型计算机上可完成半径27°(分辨率5’)的基函数模型解算。5.针对极区球坐标格网变形问题,提出基于换极法的极区点质量核径向基函数建模方法。在基函数建模过程中引入球坐标换极法,将两极地区数据网络和基函数网络重新分布,解决了两极地区传统基函数不适用问题,增强了两极地区点质量核基函数模型的稳定性,同时方便了模型程序的统一设计。6.建立了全球“拼接”点质量核基函数模型。为适应远距离飞行器的扰动引力赋值需求,结合局部多尺度点质量核模型的特性,采用“分割-拼接”法建立全球范围的多尺度模型。该模型顾及每一层点质量基函数的边界效应和计算设备效能设定拼接准则和数据存储形式。实验证明,该模型可在全球范围内实现扰动引力赋值服务,3Km高度层精度可控制在0.5mGal。