利用扰动重力数据反演海底地形的高斯曲面函数估计方法

利用2017年1月—2020年12月的SARAL卫星测高数据反演得到南海海域1′×1′格网扰动重力数据,通过与船载重力数据比较,精度达到5.5 mGal。提出了利用重力数据及高斯曲面函数求解区域特征参数进而估计海底地形模型的方法,利用SARAL卫星反演获得的重力数据在南海海域开展了计算试验。结果表明,利用ETOPO-1先验模型,在10×10格网一组划分条件下估计得到的1′×1′海底地形精度相对于先验模型提高了约10 m;利用DTU18先验模型,在9×9格网一组划分条件下估计得到的1′×1′海底地形精度相对于先验模型提高了约9 m。计算结果表明,利用卫星测高反演得到的重力数据,通过高斯曲面函数解算获得的5个特征参数,可以在一定程度上代表相应区域海底地形的曲面特征,进而可在不依赖外部实测数据条件下对先验海底模型进行精化。对于曲面估计而言,格网划分越小,曲面函数就越能反映区域变化特征。因此,对于未来卫星测高技术发展而言,更高分辨率的重力场探测技术有望继续提升海底地形细节的反演能力。     

航空和地面重力测量数据的一体化处理方法

在三维经典扰动位协方差模型基础上,推导获得局部重力异常协方差模型、局部扰动重力协方差模型,以及扰动重力和重力异常的互协方差模型,这3个模型在忽略小项后是一致的,即航空扰动重力数据与地面重力异常数据的联合处理采用一个模型即可,大大简化了处理流程.以最小二乘配置理论及本文推导的协方差模型为基础,构建了航空扰动重力测线数据与地面重力异常数据一体化处理的具体方法,利用该方法可推估地面待求点的重力异常数据,且处理过程自动完成了航空扰动重力数据的向下延拓,同时也大大削弱了航空数据中存在的系统偏差影响.在某试验区开展了一体化处理试验,结果表明,采用一体化处理得到的重力异常数据比单独使用5'分辨率航空重力数据推估获得的数据精度提高2.4 mGal,比单独使用2.5'分辨率地面重力数据推估获得的数据精度提高1.2 mGal.     

航空重力测量的分辨率和精度分析

研究了航空重力测量分辨率与飞行高度的关系,表明对于300 km/h的飞行速度和2.5 km的飞行高度在山区和平地可恢复的最小波长分辨率分别为9 km和14 km,在此高度该频段重力异常的衰减率约为50%.探讨了低通滤波器截止频率对航空重力测量沿线分辨率和精度的影响,对于大同航空重力测量,滤波尺度为150,200,250 s时,沿线的半波长分辨率分别为7.5,10和12.5 km,相应的精度分别为7.5,6.2和5.5 mGal.     

基于D-P算法和最优路径的海山地理实体边界划定方法

海底地理实体的划定与命名是当前国际海洋权益领域的热点研究问题之一。然而,由于海底地理实体边界的量化界定技术缺乏,导致其界线确定不可避免地存在人为性,为此,本文提出了一种基于D-P算法和最优路径的海底地理实体边界划定方法。将高分辨率水深模型转换为二维水深矩阵,从横向与纵向进行剖面分析,采用极值点简化和D-P算法进行二次简化的方法,通过坡度变化定位和基部高程定位等基本条件综合判断,从而实现了单体海山的山体与其基部的自动划分;在此基础上,采用路径寻优的方法实现了连体型海山的自动分割。该方法在南海海底地理实体划定中进行了验证,取得了良好的应用效果。     

方位角对斜路径对流层延迟的影响研究

针对4种不同气候条件的地区，分析不同高度角时方位角对斜路径对流层延迟的影响，进而提出一种对流层延迟水平梯度的精化方法。实验表明，该方法可更真实地反映对流层延迟在水平方向的非对称性，有效提高VLBI解算精度，基线可重复度提高约4%，天顶对流层延迟估计精度提高约15%。
     

GRACE时变重力场各向异性高斯组合滤波方法

受测量误差等因素影响，直接使用GRACE时变重力场模型的地表质量变化反演结果呈现严重条带噪声，必须采用滤波消除。本文对不同滤波方法进行了试验分析，以信噪比最大为准则，确定了不同滤波方法的最优滤波参数，并在此基础上提出了一种各向异性组合滤波方法。该方法根据时变重力场模型球谐系数误差特性，结合各向异性高斯滤波和均方根滤波特点，对精度较高的低次项系数采用较大权重以保留更多有效信号，而对精度较差的高次项系数采用较小权重以压制噪声。不同于传统的两步法组合滤波，该方法仅需进行一步滤波处理。试验结果表明，本文提出的各向异性组合滤波方法计算步骤简单，能够有效消除条带噪声；与单一滤波和传统两步法组合滤波方法相比，提高了反演结果信噪比，保留了更多真实信号。     

IGDR数据用于高度计定标先期计算的可行性

针对拟开展的Jason-3卫星高度计定标试验,分析了利用IGDR数据先期计算定标结果的可行性。介绍了基于GPS浮标的直接绝对定标方法,设计了IGDR数据用于定标计算的可行性评估方法,利用经过全球4大定标站海域的数据比较了IGDR、GDR数据计算定标结果的差异。结果表明:通过两种数据计算的定标结果之间没有明显的系统性偏差,其差异的不确定性约为1cm;相对于传统定标站GDR定标结果 3cm量级的标准差,该不确定性对于定标结果的一致性影响不大,在定标试验中可使用IGDR先期进行定标计算。     

重力异常阶方差模型的构建及在扰动场元频谱特征计算中的应用

利用最新的全球引力位模型-EGM2008对经典的重力异常阶方差模型进行了分析比较,分析表明,经典的阶方差模型由于限于当时的观测条件,已经不能准确地描述扰动场元在各个频段的频谱分布。在Moritz阶方差模型基础上,利用EGM2008位模型获得的2160阶阶方差重新构建了新的分段重力异常阶方差模型-TSD模型,该模型与EGM2008位模型计算的阶方差比较其标准差和均值分别为0.25mgal2 、0.0 。利用TSD模型计算了不同频段内大地水准面高、重力异常、扰动重力、垂线偏差四个重力场扰动场元的频谱特征,计算结果表明：扰动场元频谱分布较之传统分析结果有较大的变化,其中重力异常、扰动重力及垂线偏差在中、低频部分的能量有明显的增加而高频及甚高频部分的比重有明显的减少。     

GRACE RL06与RL05时变重力场模型数据初步比较分析

简要介绍CSR、GFZ和JPL机构的GRACE RL06时变重力场模型数据,并对比分析RL06和RL05数据的解算模型。从全球陆地水储量变化反演结果、时变重力场模型阶方差和C20项时间序列3个方面,对2004-01~2014-11期间RL06和RL05时变重力场模型数据进行对比分析。结果表明,GRACE RL06时变重力场模型数据质量的和精度较RL05有明显提高,其全球陆地水储量变化反演结果去条带噪声效果更好、信噪比更高;在高阶项部分,RL06模型数据的阶方差小于RL05;RL06模型数据的C20项时间序列幅值变化小于RL05,与SLR所得C20项数据也更接近。相同条件下,采用CSR RL06模型阶方差最小,利用RL06模型所得全球陆地水储量变化反演结果信噪比值最大。