基于3DDelaunay剖分算法的重力建模与分析

地质体的重力建模是正确解释和应用重力资料的关键问题之一.针对非规则形状变密度的三度体,本文提出了基于3D Delaunay剖分算法的重力建模方法.采用3D Delaunay剖分算法将三维目标地质体分解为若干变密度四面体体元,推导了基于四面体体元的重力正演公式,建立了剩余密度值与重力异常值的线性方程组;以变密度的长方体和倾斜台阶组合体为例,比较分析了常规块体算法和3D Delaunay剖分算法应用于重力正演的有效性,并采用共轭梯度法加密度约束条件对非规则形状变密度的倾斜台阶组合体进行了密度反演.计算结果验证了本文方法的正确性和有效性.基于3D Delaunay剖分算法的重力建模可应用于存在褶皱、断层、裂缝等复杂地质体的重力正反演计算.     

基于泊松小波径向基函数融合多源数据的局部重力场建模

融合多源数据的高精度、高分辨率的局部重力场建模是物理大地测量学的前沿和热点问题.本文研究了基于径向基函数融合多源数据的局部重力场建模方法,利用Monte-Carlo方差分量估计实现了不同类型的观测数据的合理定权,引入了最小标准差法确定基函数的适宜网络,分析了地形因素对于基函数网络确定及局部重力场建模精度的影响.以泊松小波基函数为构造基函数,结合残差地形模型,融合实测的陆地重力异常、船载重力异常及航空重力扰动数据构建了局部区域陆海统一的似大地水准面模型.研究结果表明:引入残差地形模型平滑了地形质量引入的高频扰动信号,简化了基函数的网络设计;并提高了重力似大地水准面的精度,平原地区其精度提高了4 mm,地形起伏较大的山区其精度提高了约5 cm.总体而言,基于"三步法"构建的局部重力似大地水准面在荷兰、比利时及德国相关区域,其精度分别达到1.12 cm、2.80 cm以及2.92 cm.     

卫星重力梯度数据的归算与格网化处理

给出引力梯度数据归算的二阶径向改正公式,计算结果表明选用高精度参考重力场模型可以有效控制归算误差的影响。比较加权平均法、Shepard曲面拟合法和最小二乘配置法用于卫星重力梯度数据格网化处理的精度和适用性,采用不同噪声背景的径向引力梯度数据进行格网化处理,计算结果表明：最小二乘配置法相比加权平均法和Shepard曲面拟合法具有明显优势,其格网化精度可满足高精度重力场恢复的需要,建议在实际计算中采用最小二乘配置法进行格网化处理。     

地形横向密度扰动对于区域大地水准面建模的影响

分别利用Helmert和KTH法模拟试算了不同地形起伏区域内横向密度扰动对于大地水准面的影响。结果表明,在地形起伏较小的平原或丘陵地带,横向密度扰动的影响一般不会超过cm级,在精度允许的范围内可以忽略其影响;在地形起伏较大的山区,横向密度扰动的影响可达cm甚至dm级,对于cm级大地水准面精化而言,需要考虑其影响。此外,基于Helmert与KTH法解算的结果在地形起伏剧烈的山区差异较大,由于Helmert法计算地形效应时忽略了外区地形质量的影响,其解算结果可能存在偏差。总体而言,在地形起伏剧烈的山区建模时,需利用多种计算方法基于实测数据分别进行试算,通过计算的变密度大地水准面与实测的GPS水准数据进行比较分析,得到适合该区域的大地水准面建模方法。     

OpenMP并行算法在卫星重力场模型反演中的应用

利用卫星重力数据反演地球重力场需要解决重力场模型的高效计算问题。分析了最小二乘直接法求解重力场模型涉及的密集型计算任务,基于OpenMP实现了卫星重力场模型直接求解的并行算法。利用30天、5秒采样间隔的沿轨扰动位T和径向扰动重力梯度Tr数据,分别反演了60阶次的卫星重力场模型,计算结果表明,OpenMP并行算法能够有效提高直接法求解卫星重力场模型的计算效率,并具有很好的稳定性。     

�������������������ص�FFT�����㷨�Ƚ�

???????????????????????????????Possion???????????????????????FFT????????????????????????????????λ???EGM2008??????????????????????????????????????????о???????????????????????????￡???Possion?????Ч????á?     

基于多尺度边缘的重力异常分离方法

利用多尺度边缘分析理论,对重力异常的小波多重分解法作了改进,提出了基于多尺度边缘的重力异常分离方法。实验表明,改进方法除保留了多重分解法的优点外,还能够同时实现异常的纵向和横向分离,而且可以根据每个场源对应的多尺度边缘的分布情况确定小波分解的阶数,有利于单一异常特征的提取。利用该方法对琉球俯冲带西南端区域的船测重力资料进行了处理分析,结果与该区域的地质构造具有很好的对应关系。     

���������ݶȻָ�������������ʱ���о�

???????????????????????Ч???????????????????????????????????? ???GOCE??????????????????????????????????????????????????????????????????????????????С????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????С?????????????????????????????????????????     

利用改进FLIC算法进行重力矢量场可视化

针对固定滤波长度公式得到的LIC图像显示效果不理想的缺点,利用纹理流线上所有点亮度值的统计方差来衡量可视化精度,提出了一种滤波长度随矢量变化的改进FLIC算法。用颜色信息表示重力异常大小或者真实垂线偏差大小,纹理方向表示重力方向,基于改进FLIC算法实现了重力矢量场的可视化,结果表明,重力矢量场的可视化具有表达多维信息和矢量方向细节、图形直观易读等优势。     

中国大地测量研究进展(2019-2023)EI北大核心CSCD

2023年7月11日至2023年7月20日,第28届国际大地测量学和地球物理学联合会(The International Union of Geodesy and Geophysics,IUGG)大会在德国柏林举行。按照IUGG传统,国际大地测量协会中国委员会(CNC-IAG)组织国内十余家单位编写了“中国大地测量国家报告(2019—2023)”,分别总结了2019至2023年4年期间的中国大地测量各分支学科研究进展。本文主要归纳和总结了中国大地测量学科近几年的整体进展,侧重各领域代表性研究进展,主要内容包括基准框架、综合PNT与弹性PNT、重力场与垂直基准、GNSS精密产品、多源传感器组合导航和海洋大地测量6个研究方向。此外,结合国际大地测量及相关交叉学科的发展趋势,对我国大地测量学科未来发展提出了几点建议。