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

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

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

根据观测面和延拓面测量数据的Poisson积分平面近似关系,结合快速傅立叶变换算法,将向下延拓转换到频率域进行计算,并采用迭代Tikhonov正则化方法,克服计算的不稳定性,提高计算结果的精度,实现了航空重力测量数据的向下延拓。最后采用模拟航空重力测量数据验证了该算法的有效性,取得了较好的延拓结果。     

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

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

航空重力向下延拓的最小二乘配置Tikhonov正则化法

基于最小二乘配置法向下延拓航空重力的过程中,由于协方差矩阵严重病态,影响延拓结果的稳定性和精度。针对这一问题,提出了航空重力向下延拓的最小二乘配置Tikhonov正则化法。基于全球协方差函数模型建立航空重力数据与地面重力数据的协方差关系,引入基于广义交叉验证法,选择正则化参数的Tikhonov正则化法改善协方差矩阵的病态性,抑制观测噪声对延拓结果的放大影响。基于EGM2008重力场模型,设计了山区、丘陵和海域3种不同地形区域的航空重力数据向下延拓的仿真实验,实验结果验证了该方法的有效性。     

航空重力测量数据向下延拓基于信噪比的正则化方法的研究

本文研究了正则化方法在航空重力测量数据向下延拓问题中的应用。首先对这种不适定问题的线性模型,分析了设计阵的复共线性结构与其对参数估计危害之间的关系,利用参数LS估计的信噪比提取了各个参数是否受到复共线性严重危害的信息,从而在一定程度上揭示了设计阵复共线性结构的特征。然后提出了基于信噪比的正则化方法（SNR）,以信噪比为依据构造正则化矩阵,以极小化均方误差为目标选取正则化参数。本文构造正则化矩阵无需利用附加物理或先验信息,这对于在缺乏此类信息的情况下运用正则化方法提供了新的手段。最后进行了数值试验,结果表明,本文提出的新方法（SNR）比普通的正则化方法（OR）在滤噪和保真方面表现更佳。     

基于矩谐分析的航空重力向下延拓

提出基于矩谐分析的航空重力向下延拓方法,以重力扰动作为基本观测值,给出基于矩谐分析的向下延拓模型和算法。利用EGM2008重力位模型设计模拟数值试验,对比研究矩谐分析、直接法和基于广义岭估计的逆泊松积分法,分别采用这3种方法将飞行高度处含高斯白噪声的2.5′×2.5′重力扰动向下延拓至大地水准面,与真实值作外部检验。数值比较结果表明:矩谐分析在延拓精度、稳定性和边界效应等方面都要优于直接法和基于广义岭估计的逆泊松积分法,能取得良好的向下延拓效果。     

航空重力数据向下延拓的波数域迭代Tikhonov正则化方法

航空重力数据向下延拓是重力场数据联合处理的重要步骤。研究了Tikhonov正则化方法以及迭代Tikhonov正则化方法在航空重力数据向下延拓中的应用,指出Tikhonov正则化方法在利用GCV法或L曲线法选取正则化参数时存在的不足以及空间域迭代法迭代不收敛的问题。波数域Tikhonov迭代正则化法的引入,有效解决了上述问题。算例结果表明,波数域迭代法迭代过程收敛,且计算精度高、速度快,值得广泛应用于航空重力数据的向下延拓。     

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

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

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

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

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

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

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

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

航空重力向下延拓的误差影响分析

本文从分析航空重力向下延拓过程中偶然误差和系统误差的变化特性入手,进而提出处理办法。首先,利用试验说明移去恢复法局限性,同时表明需处理系统误差和偶然误差的必要性。然后,采用理论推演和数值模拟计算分别估计了系统误差和偶然误差影响,试验结果发现:系统误差影响和偶然误差影响均与数据格网间隔、向下延拓高度呈线性关系,当格网化间隔较小和延拓高度较高时系统误差影响和偶然误差影响较大。最后,提出使用半参数模型和正则化算法的两步法估计系统误差和减弱偶然误差影响,试验结果说明两步法处理向下延拓各类误差影响优于仅用半参数模型或正则化算法的结果,在试验数据的偶然误差标准差为2×10~(-5) m/s~2、恒值系统误差3×10~(-5) m/s~2和变值系统误差标准差约1.3×10~(-5) m/s~2时,以及向下延拓高度6.3 km和格网间隔6′的条件下,两步法向下延拓结果的精度可达2.3×10~(-5) m/s~2。     

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

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

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

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

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

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