航空重力测量数据向下延拓中空间协方差函数特性研究

本文通过对航空重力测量数据的分析,建立起具有空间相关特性的空间协方差函数模型。利用所建立的模型,将航空重力测量试验数据采用最小二乘配置的方法向下延拓到地面,与地面已知重力数据进行比较。由比较结果,分析了空间协方差函数模型对航空重力测量数据向下延拓结果的影响,并表明所建立的空间协方差函数模型,满足协方差函数特性,且在局部区域使用能够达到满意的精度要求。     

向下延拓航空重力数据的Tikhonov双参数正则化法

为了避免正则化参数对向下延拓过程可靠成分的修正影响,提出了Tikhonov双参数正则化法。引进截断参数,将法矩阵的奇异值分为相对较大的奇异值（可靠部分）和相对较小的奇异值（不可靠部分）;引进正则化参数,只对法矩阵的小奇异值进行修正,以抑制高频误差对向下延拓解的影响。采用改进的广义交互确认法（GCV）确定截断参数和正则化参数。基于EGM2008重力场模型仿真了一组航空重力数据,验证了该方法对航空重力数据向下延拓过程的有效性。     

重力异常向上延拓严密改化模型及向下延拓应用

重力异常向上延拓全球积分模型在航空重力测量数据质量评估和向下延拓迭代计算等领域具有广泛的应用。为了消除积分核函数奇异性影响,需要对该模型进行基于积分恒等式的移去-恢复转换及全球积分域的分区改化处理。在此过程中,传统改化处理方法往往忽略了全球积分过渡到局域积分引起的积分恒等式偏差影响,从而导致不必要的计算模型误差,最终影响向上延拓计算结果的可靠性,甚至影响向下延拓迭代解算结果的稳定性。针对此问题,本文开展了重力异常向上延拓积分模型改化及向下延拓应用分析研究,依据实测数据保障条件和积分恒等式适用条件要求,导出了重力异常向上延拓积分模型的分步改化公式,提出了补偿传统改化模型缺陷的修正公式,并将最终的严密改化模型应用于重力异常向下延拓迭代解算。使用超高阶地球位模型EGM2008作为标准位场开展数值计算检验,分别对重力异常向上延拓分步改化模型的计算精度及在向下延拓迭代解算中的应用效果进行了检核评估,验证了采用严密改化模型的必要性和有效性。     

虚拟压缩恢复法在向下延拓问题中的应用

阐述了基于虚拟压缩恢复法的重力场向下延拓,构造了点质量处于内部球内部、内部球外部和在内部球内外各有一个点质量三种点质量异常模型,试验结果支持基于虚拟压缩恢复法的向下延拓。     

重力场的局部虚拟向下延拓以及利用EGM96模型的模拟实验检验

阐述了基于虚拟压缩恢复法导出的虚拟向下延拓法,该法从理论上解决了重力场的向下延拓问题。为了检验理论的正确性,选取EGM96模型进行了模拟实验检验,实验结果支持局部虚拟向下延拓法。同时,对观测值加进随机误差的情形也进行了模拟实验检验,结果表明,局部虚拟向下延拓几乎不损失精度。     

航空重力数据向下延拓的倒锥法(上)——方法的提出和存在性

基于离散 Possion积分提出了一种灵活直接的向下延拓新方法 ,证明了方法的存在性 ,利用模拟数据进行的数值实验取得了很好的效果。     

验证引力位虚拟压缩恢复法以及虚拟向下延拓

采用15′×15′网格离散化手段对均质旋转对称椭球体产生的引力位场进行了数值模拟计算,在相对精度优于10-4的水平上,验证了引力位虚拟压缩恢复法以及重力场虚拟向下延拓法的可靠性和有效性。     

重力异常向下延拓的虚拟压缩恢复法与Bjerhammar虚拟球法的比较

边值问题一直以来都是物理大地测量中的重点和难点问题。随着卫星重力测量和航空重力测量应用的日益广泛,研究重力数据精密稳健的向下延拓方法在重力场研究中越来越重要。本文分别对基于等效虚拟法的虚拟胝缩恢复法（虚拟向下延拓法）和Bjerhammar虚拟球泫进行了简单介绍,并采用模拟数据对两种方法的延拓结果进行了对比。     

基于多尺度边缘约束的重力场信号的向下延拓

随着航空重力测量技术和卫星重力技术的发展,越来越多的空间重力资料可用于地球重力场的确定,但是由于空间重力观测值不在地球表面或大地水准面上,为了联合地面重力资料通常需要对空间重力数据进行向下延拓。由于重力场对高度敏感,向下延拓过程是不适定的,观测     

重力与磁力测量数据向下延拓中最优正则化参数确定方法研究

向下延拓是重、磁测量数据处理的关键步骤之一,然而,向下延拓是一个典型的不适定问题,需要采用正则化方法实现有效延拓,因此,正则化参数的确定是重、磁测量数据向下延拓正则化方法研究中最重要内容。本文根据观测面和延拓面测量数据的Poisson积分平面近似关系,结合快速傅里叶变换算法,将其转换到频率域进行计算,提高了计算速度,为了克服计算的不稳定性并进一步提高计算结果的精度,引入Landweber正则化迭代法,在此基础上采用L曲线法研究了最优正则化参数的确定,最后采用模型磁测数据验证了所确定的正则化参数的有效性,并取得了较好的延拓结果。     

顾及地形效应的重力向下延拓模型分析与检验

向下延拓是航空重力测量数据实际应用中必不可少的技术环节。向下延拓属于不适定反问题,其解算过程具有较大的不确定性,故该问题一直是大地测量领域国内外学者的研究热点。本文深入分析研究了当前国内外最具代表性的3种向下延拓计算模型的技术特点和适用条件,提出了应用超高阶位模型、局部地形改正和移去—恢复技术顾及地形效应,以及位场延拓结果球面化曲面的工程化方法,重点探讨了计算模型的稳定性及数据观测误差对延拓计算结果的影响。通过理论分析、数值仿真和实测数据计算等手段,定量评估了不同向下延拓模型的解算精度及其可靠性。其主要结论是:传统逆Poisson积分模型解严重受制于输入数据观测噪声的干扰,在现有作业条件下,该模型至多只能用于1km以下高度的延拓解算;频谱截断积分和位模型加地改两种延拓新模型具有良好的计算稳定性,完全适用于2′分辨率和5km飞行高度条件下的航空重力测量数据向下延拓解算,其延拓计算精度可达2×10~(-5) m/s~2,可满足各方面实际应用需求。     

航空重力向下延拓病态问题的求解

提出将广义岭估计用于求解航空重力向下延拓病态问题,研究了求解逆Poisson积分问题的3种正则化方法：Tikhonov正则化、岭估计和广义岭估计。利用EGM2008地球位模型设计模拟数值实验,将飞行高度处含白噪声的2.5′×2.5′重力扰动向下延拓至大地水准面上,与参考值作外部检验,全面检验、比较了各向下延拓方法的可靠性、精度和稳定性,数值结果表明基于多个最优正则化参数的广义岭估计在延拓精度、稳定性和抗差性等方面要显著优于基于单个正则化参数的Tikhonov法和岭估计。     

航空重力向下延拓分组修正的正则化解法

在航空重力向下延拓中,针对病态性对解算结果各个部分影响的不同进行分组修正,提出了分组修正的正则化解法。利用信噪比评估病态性影响,得到将参数分组的方式;依分组修正思想构造正则化矩阵;通过极小化均方误差选取正则化参数。基于EGM2008重力场模型仿真一组航空重力数据,验证了该方法对航空重力数据向下延拓过程的有效性,并与另外3种方法作比较。结果表明,新方法具有更高的精度。     

航空重力测量数据向下延拓的正则化算法及其谱分解

基于Poisson积分方程,提出了以地面重力观测值作为控制并顾及外区影响的向下解析延拓数学模型,推导了向下解析延拓的谱分解式,在频域内分析了造成向下延拓结果不稳定的原因,进而给出了向下延拓的正则化算法,并讨论了向下延拓中的地形影响.通过对我国首次航空重力测量试验数据的处理表明,提出的方法可获得稳定、精确的向下延拓结果.     

基于广义补偿最小二乘法的航空重力向下延拓

航空重力向下延拓是病态问题,而广义补偿最小二乘法可以很好地克服病态性。研究了基于广义补偿最小二乘法的逆Poisson积分的航空重力向下延拓模型,并设计实验方案,将EGM2008地球位模型计算的重力异常作为仿真实验的数据,分别用最小二乘、Tikhonov正则化、广义补偿最小二乘3种方法求解,对其精度及仿真效果进行比较。结果表明,广义补偿最小二乘方法精度高,对仿真效果有显著提高。     

航空重力测量向下延拓破损点的搜寻与修复

航空重力测量数据向下解析延拓通常是一个病态方程的求解问题,从空中延拓到地面上的重力异常数据因受病态干扰质量不等。把其中质量较差的部分挑选出来,然后再进行有针对性的处理,这种做法的好处是尽量减少对于质量较好数据的破坏,以保持其原有的完好性。提出了两种破损点修复的方法:一是部分岭估计,可在不增加任何新观测的情况下实施;二是部分贝叶斯估计,利用已有的地面数据对破损点做有节制的修复。两种方法均具有较强的可操作性。     

联合使用位模型和地形信息的陆区航空重力向下延拓方法

为了规避传统逆Poisson积分向下延拓解算过程的不适定性问题,借鉴导航定位中的"差分"概念,利用超高阶位模型直接计算海域航空重力测量向下延拓改正数的方法。本文在此基础上提出联合使用重力位模型和地形高数据,计算陆部航空重力向下延拓总改正数的改进方案,以飞行高度面与地面对应点的位模型差分信息表征总改正数的中长波分量,以相对应的局部地形改正差分修正量表征总改正数的中高频成分,从而实现航空重力数据点对点向地面的全频段延拓。在地形变化不同区域,联合使用EGM2008位模型、地面实测重力和高分辨率高程数据进行了实际数值计算和精度评估,验证了该方法的有效性。     

基于改进高斯-牛顿法的位场向下延拓

本文针对位场向下延拓的不适定问题,在分析最优化算法中高斯-牛顿法基本原理及滤波函数滤波特性的基础之上,提出基于正则参数指数递增计算方法和残差最小步长准则的改进高斯-牛顿法。基于理论重力模型和航磁实测数据的对比实验表明,改进后得到的自适应迭代法具有相对较高的位场向下延拓精度和很好的收敛性。