地球物理大地测量
大型科学
计算平台

CASM.jpg

高精度重力场逼近与大地水准面计算系统PAGravf4.5
轮播多图

高精度重力场逼近与大地水准面计算系统PAGravf4.5,严格遵循物理大地测量理论与技术要求,通过构造不同高度、各类场元、多种地形影响统一的解析算法体系,提升地球物理重力探勘与重力场数据处理水平;通过构建边值理论框架中重力场积分与径向基函数逼近算法体系,实现由某边界面任一种扰动场元,确定整个外部空间各种场元的计算能力,有效提升重力场全要素建模能力及大地水准面精化水平;研发基于重力场数据和理论方法的特色算法,优化区域高程基准,拓展重力场及高程基准成果的应用服务水平。

用户可根据工程和研究需要,设计个性化的作业流程,灵活组织有关功能模块,完成空天地海重力场数据处理、各种地形影响计算与分析工作,实现复杂多源异质重力场数据的有效综合,确定高精度陆海统一的重力场及大地水准面,精化区域高程基准,开展重力场及高程基准应用服务计算及分析工作。

PAGravf4.5程序包—2022年10月14日更新
PAGravf4.5参考说明书.pdf
39.11MB下载
PAGravf4.5_win64cn.exe
650.0MB下载
PAGravf4.5_win64cn.2
628.02MB下载
两种大地水准面精化简捷流程练习.pdf
640.91KB下载
PAGravf4.5cn help.chm
57.54MB下载

✍请务必将PAGravf4.5_win64cn.exe和PAGravf4.5_win64cn.2一起下载到同一个文件夹后,再安装(Windows10平台)。下载后,若防护软件将.exe隔离,请去隔离区找到文件,恢复信任。

地球潮汐负荷影响与形变监测计算系统ETideLoad4.0
轮播多图

地球潮汐负荷影响与形变监测计算系统ETideLoad4.0,主要由大地测量全要素各类固体地球潮汐影响计算,地面大地测量非潮汐时间序列处理分析,高分负荷形变场及时变重力场逼近计算,CORS/InSAR协同监测与地面稳定性估计,以及大地测量数据编辑计算与可视化五大子系统构成。整个系统在Visual studio 2017 x64集成环境中研发,采用QT C++(界面)、Intel Fortran(核心功能模块)和mathGL C++(数据可视化)代码级混合编程技术,ETideLoad4.0包含50多个win64程序和600多功能模块。

ETideLoad4.0适合大地测量、地球物理、环境灾害、水文气象、卫星动力学、地震与地球动力学等高年级本科生、研究生、科研和工程技术人员。为方便自学、教学与技术培训,绝大多数程序有完整计算样例,全部样例练习完成后(约5个工作日),基本具备独立使用软件系统的能力。

用户可根据工程或研究需要,设计个性化技术流程,灵活组织有关功能模块,完成地面及固体地球外空间各种潮汐/非潮汐影响监测计算,全球或区域形变场/时变重力场/地面稳定性/陆地水/地表动力环境变化监测,以及多源异构大地测量深度融合等工作。

ETideLoad4.0程序包—2022年3月9日更新
ETideLoad4.0参考说明书.pdf
32.86MB下载
ETideLoad4.0_win64cn.exe
650.0MB下载
ETideLoad4.0_win64cn.2
493.14MB下载
ETideLoad4.0cn help.chm
51.87MB下载

✍请务必将ETideLoad4.0_win64cn.exe和ETideLoad4.0_win64cn.2一起下载到同一个文件夹后,再安装。下载后,若防护软件将.exe隔离,请去隔离区找到文件,恢复信任。

ETideLoad典型作业流程
负荷形变场与时变重力场计算

 目的:由区域高分大气压、地表水和海平面变化数据,结合全球大气压、海平面和土壤水变化格网模型,计算高精度区域负荷形变场与重力场变化格网模型。

方法:由区域高分辨率地表环境负荷数据,采用基于全球负荷形变场球谐系数模型的移去恢复法,计算区域非潮汐负荷形变场及时变重力场。

(1)将全球大气压格网转换成地面等效水高格网,由“格网陆地海洋数据分离”,消除全球土壤水和海平面变化格网在海岸线附近的矛盾。

(2)将全球大气压、海平面和土壤水等效水高变化格网模型相加,得到全球等效水高变化格网模型。

(3)利用“全球等效水高格网规格化球谐分析”,由全球等效水高变化格网模型生成全球等效水高规格化球谐系数模型。

(4)整合区域大气压、地表水和局部海平面变化数据,生成区域高分辨率等效水高变化格网模型。

(5)利用“等效水高球谐综合负荷形变场计算”,由全球等效水高规格化球谐系数模型,计算与区域等效水高变化格网规格相同的等效水高模型值。

(6)移去模型值:利用“两格网文件的加权运算”,从区域等效水高变化格网中移去等效水高模型值,得到区域剩余等效水高格网模型。

(7)利用“区域格林函数积分负荷形变场计算”,由区域剩余等效水高格网模型,确定区域剩余负荷影响格网模型。

(8)利用“等效水高球谐综合负荷形变场计算”,计算与区域剩余负荷影响规格相同的全球负荷影响模型值。

(9)恢复模型值:利用“两格网文件的加权运算”,将区域剩余负荷影响格网与全球负荷影响模型值格网相加,得到区域高精度非潮汐负荷影响格网模型,包括各种类型的负荷形变场与重力场变化格网模型。

地球卫星潮汐与非潮汐摄动影响计算

目的:充分利用与综合各种数据资源,精确计算全部类型潮汐与非潮汐地球动力学影响导致的卫星轨道摄动,以期全面提升大地测量卫星的地球重力场、地球变化监测和卫星动力学定轨水平。

(1)利用“卫星轨道摄动潮汐影响计算”,逐一选择潮汐影响类型,分别计算卫星重力位、摄动力的固体潮、海潮大气潮负荷、地球极移和海潮极潮影响。

(2)构造卫星轨道面大地高格网模型,并依据卫星运行时间,构建全球非潮汐海平面、大陆水和大气压(月/半月)变化格网模型时间序列。

(3)利用“全球等效水高格网规格化球谐分析”,由全球非潮汐海平面、大陆水和大气压(月/半月)变化格网模型时间序列,生成相应的非潮汐负荷规格化球谐系数模型时间序列。

(4)利用“球谐系数模型加权运算”,将全球非潮汐海平面、大陆水和大气压球谐系数模型相加,生成全球非潮汐负荷球谐系数模型时间序列。

(5)利用“地球外部空间非潮汐负荷影响计算”,由卫星轨道面大地高格网模型和全球非潮汐负荷球谐系数模型,选择大地坐标格式,分别计算卫星轨道面上摄动位和摄动力三分量非潮汐负荷影响格网时间序列。

(6)利用“格网时序插值任意时刻点值”,由卫星星历文件和卫星轨道面上非潮汐负荷格网时序,按经纬度坐标和时间,内插卫星轨道在实际时刻的非潮汐负荷影响。

CORS网大地高变化负荷形变场精化

目的:基于CORS网站点大地高变化时序,结合大气、地表水及海平面变化等观测资料,通过精化负荷垂直形变、大地水准面、地面重力及地倾斜变化格网时序,实现区域非潮汐负荷形变场及时变重力场的连续监测。

方法:基于全球和区域地表负荷模型的移去恢复法。

(1)计算全球负荷垂直形变场及重力场月变化模型值,以及区域剩余负荷垂直形变场与重力场月变化模型值。

(2)利用“非潮汐时序分析与处理”有关程序,从CORS站天解(或周解)时间中移去线性变化,构造CORS站大地高月变化(非线性变化)时序。

(3)利用“4等效水高球谐综合负荷形变场计算”,从CORS站大地高月变化时间序列中,移去全球负荷影响(大地高月变化、负荷垂直形变)模型值。

(4)继续利用“等效水高球谐综合负荷形变场计算”,选择计算负荷影响类型为“等效水高值”,从区域等效水高变化格网中,移去全球负荷等效水高模型值,生成区域剩余等效水高变化格网模型时间序列。

(5)利用“区域格林函数积分负荷形变场计算”,由区域剩余等效水高变化格网模型时间序列,从CORS站大地高月变化时间序列中,再移去区域剩余负荷影响(大地高月变化)模型值,生成CORS站剩余大地高月变化时序。

(6)利用“站点大地测量参数负荷形变场精化”,由CORS站网各站点剩余大地高月变化时间序列,生成区域剩余等效水高月变化时序。

(7)利用“区域格林函数积分负荷形变场计算”,由区域剩余等效水高月变化格网模型时序,生成区域剩余地面大地高、大地水准面、地面重力与地倾斜月变化精化值格网时序。

(8)将地面大地高、大地水准面、地面重力与地倾斜变化格网时间序列的全球负荷影响模型值、区域剩余负荷影响模型值和区域剩余负荷影响精化值相加,得到区域地面大地高、大地水准面、地面重力与地倾斜月变化格网时间序列。

CORS时序InSAR地面沉降综合监测

目的:通过InSAR监测量野值探测,极浅地表自身变化分离,InSAR监测量CORS网时空监测基准统一,以及InSAR监测量时序分析与时变信号重构等,优化整合CORS网非线性时变信号全天候连续监测能力和时变信号无误差累积优势,与InSAR高精度地面垂直变化空间差异性的高分辨率监测能力,实现统一时空监测基准中高时空分辨率垂直形变(地面沉降)的高精度监测。

技术方案:

(1)依据地面沉降空间高相关特点,由InSAR监测量构造平滑的参考格网,以分离InSAR野值、粗差和突变等非沉降信息信号,分离SAR多路径效应和复杂粗糙地表环境干扰。

(2)依据沉降量与动力源/作用点距离或距离平方近似反比的空间性质,对每期InSAR监测量进行高斯滤波,使其二阶可微,以分离极浅地表层非沉降信号,抑制大气延迟短波扰动及地表自身变化。

(3)利用InSAR监测量与CORS网地面沉降量中长波成分的差异,将InSAR监测量时序统一到CORS网监测基准,以控制InSAR监测量随时间的误差累积,补偿或校正大气延迟与潮汐形变引起的空间变化信息和时变信号,修复空间尺度大于InSAR干涉区的时变信号,将InSAR监测量统一转为正(变化。

20200611(4)依据动力(加速度)是沉降量(位移)对时间的二阶导数,对InSAR沉降量进行时序分析与重构,使其连续可导,以提高沉降量的可靠性和精度水平,实现地面沉降的情势分析与短时推估预报。

说明:

(1)由CORS网地面沉降格网范围覆盖各组InSAR时序监测区,实现多组InSAR监测量时序的空间监测基准统一。

(2)按相同采样时刻,对各组InSAR沉降量进行信号重构,就能实现多组时序InSAR监测成果的深度融合。

CORS网地面稳定性变化定量监测

目的:从CORS网地面高程、重力和地倾斜变化格网时序中,提取地面稳定性变化信号,估计CORS网作用范围内总体地面稳定性,定量跟踪CORS网运行期间,地面稳定性降低的时间与地点、持续作用时间与空间作用范围,揭示地面稳定性时空演变动态。

(1)利用稳定性降低定量辨识准则,由CORS网地面形变场变化格网时序,直接生成地面稳定性变化格网时序。

利用“高程变化格网地面稳定性变化检测”、“重力变化格网地面稳定性变化检测”、“地倾斜变化的地面稳定性变化检测”,分别由地面高程变化、地面重力变化和地倾斜变化格网时间序列,确定3种地面稳定性变化格网时序。

(2)综合CORS网地面稳定性变化与形变场监测成果,定量估计CORS网作用范围内的总体地面稳定性格网模型。

①利用“格网格值水平梯度估计”,选择“采用极坐标表示水平梯度向量”,由CORS网地面沉降速率格网,计算地面沉降速率水平梯度格网模型。

②利用“系列格网格值统计文件构造”,分别统计基于高程、重力和地倾斜变化的地面稳定性变化格网时间序列,获得3种稳定性变化时序的标准差格网模型。

③利用“文件数据统计信息提取”,获得3种稳定性变化时序的标准差格网、地面沉降速率水平梯度格网的平均值和标准差。地面稳定性-附等值线11.png

④灵活运用“两格网文件的加权运算”,将3种稳定性变化时序标准差格网、地面沉降速率水平梯度格值分别减去其平均值,再除以各自标准差,实现3种稳定性变化时序标准差格网和地面沉降速率水平梯度格网的归一化。

⑤利用“两格网文件的加权运算”,以总体稳定性较大值与地质灾害(隐患)点的分布一致性为目标,调整权比关系,将3种稳定性变化时序的标准差和地面沉降速率水平梯度格网加权平均,定量估计CORS网作用范围内总体地面稳定性格网。

⑥以总体地面稳定性格网模型为底图,展绘该地区已有地质灾害(隐患)点位置,评价 CORS站网对整个区域总体地面稳定性的定量估计能力。

总体地面稳定性能刻画地面稳定性空间分布差异。数值越小的地方,表示该地方地面稳定性状态变化越小,可认为其附近地面总体上是稳定的,反之,总体的稳定性较差。

(3)以总体地面稳定性为背景场,优化整合基于地面高程、重力和地倾斜格网的地面稳定性变化,构建CORS网地面稳定性变化格网时序。

①利用“系列格网数据统计信息提取”,分别统计基于地面高程、重力和地倾斜格网的地面稳定性变化格网时序,获得3种时序的平均值和标准差。

②灵活运用“两格网文件的加权运算”,通过将全部格网格值减去平均值,再除以标准差,归一化3种地面稳定性变化时序与总体地面稳定性格网。

目的:统一地面稳定性变化监测基准。

③以地面稳定性降低与已有地质灾害(险情)事件基本一致为优化目标,调整3种地面稳定性变化、总体地面稳定性的权比关系。权值之和等于1。

④加权平均确定CORS站网地面稳定性变化格网时间序列。

⑤地面稳定性变化格网模型时序,反映了地面稳定性变化的当前态势(时间上),以及不同位置上当前态势之间的差异(空间上)。

⑥在整个地面稳定性变化格网模型时间序列中,格值越大,表明该时刻、该地方的地面稳定性降低幅度越大。

⑦无论什么情况下,某时某地(用时序对应时刻的格网所在水平坐标表示)的地面稳定性变化值增大,都表明此时该处的地面稳定性正在加速降低。

⑧在地面稳定性变化格网时序中,格值较小的格网所在地与时序对应时间,能作为CORS网作用范围内整个时间跨度的地面稳定性变化参考基准。

PAGravf典型作业流程
边值问题积分法大地水准面便捷建模流程练习

✍练习目的省略地形影响处理,由地面、海洋和航空重力,以及GNSS水准数据,按边值问题积分法精化大地水准面,快速了解和把握大地水准面精化要领、效果与目标区域重力场及大地水准面特性,客观认识地形影响处理对大地水准面建模的作用。

PAGravf4.5统一处理地面、海洋和航空重力,无须区分,可分别配权。

✍主要数据源:地面、船测和航空重力数据,实测GNSS水准数据。

重力数据文件记录格式:点号/名 经度 纬度…实测大地高…实测重力。

GNSS水准记录格式:点号/点名 经度(度小数) 纬度…实测大地高…实测正(常)高…。

002.jpg

(1)观测数据文件构造。调用[实际测点扰动重力场元计算],由重力数据文件,生成扰动重力观测文件;调用[点值属性加权运算],由GNSS水准数据文件,生成GNSS水准实测高程异常(或大地水准面高)文件。

(2)等位边界面大地高格网构造。计算扰动重力观测点的平均大地高,将360阶模型大地水准面高格网,加上该平均大地高,再减去360阶模型大地水准面高格网的统计平均值,生成等位边界面大地高格网。

此处的等位边界面,用重力点等正高面大地高格网近似表示,也是地面、船测、航空扰动重力的统一归算面。

(3)离散扰动重力归算与格网化。

调用[地球重力场各种场元模型值计算],设置最大计算阶数1440,计算扰动重力点处的模型扰动重力,将重力点处的观测扰动重力与模型扰动重力相减,生成离散残差扰动重力文件,进而调用[异质数据基函数插值格网化],生成残差扰动重力格网。

再次调用[地球重力场各种场元模型值计算],采用完全相等的最大计算阶数1440,由等位边界面大地高格网,生成模型扰动重力格网;进而将残差扰动重力格网与模型扰动重力格网相加(恢复),生成等位边界面上的扰动重力格网。

完成重力点归算与格网化。以下采用参考重力场(前540阶)移去–Hotine积分–恢复法,按边值理论解算地面高程异常或大地水准面高。

(4)调用[地球重力场各种场元模型值计算],最大计算阶数540,输入等位面边界面大地高格网,计算模型扰动重力格网;将扰动重力格网与模型扰动重力格网相减(移去参考重力场),生成等位边界面上残差扰动重力格网。

(5)调用[广义Hotine积分外部高程异常计算],由等位面大地高格网及其面上残差扰动重力格网,输入地面/海面大地高格网时,计算残差地面高程异常;输入360阶模型大地水准面高格网,计算残差大地水准面高格网。

(6)调用[地球重力场各种场元模型值计算],采用相等的最大计算阶数540,输入地面/海面大地高(或模型大地水准面高)格网,计算模型地面高程异常(或大地水准面高)格网;进而将地面高程异常(或大地水准面高)的模型格网与残差格网相加(恢复参考场),生成重力地面高程异常(或大地水准面高)格网。

地面/海面大地高格网(= 数字高程模型 + 模型地面高程异常)。模型地面高程异常、模型大地水准面高用360阶重力场模型计算(略)。要求格网范围大于目标区域范围,以抑制积分边缘效应。

完成边值问题解算,获得目标区域重力大地水准面成果。此后,将GNSS水准与重力大地水准面进行融合。

(7)在步骤(5)(6)中,由GNSS水准实测数据文件,按离散点分别计算GNSS水准点的残差和模型高程异常(大地水准面高),生成GNSS水准实测残差高程异常(或大地水准面高),探测并剔除粗差,计算并移去统计平均值。

统计平均值代表区域高程基准与全球高程基准(重力大地水准面)之间的系统偏差。

GNSS水准点处重力地面高程异常最好不用重力高程异常格网插值,否则由于GNSS水准点与内插点的大地高不等,需增加一项高程异常高差改正。

(8)调用[Possion方程约束GNSS水准融合],生成残差高程异常(大地水准面)格网,再与重力地面高程异常(大地水准面)格网相加,生成实用地面高程异常(大地水准面高)模型格网成果。

完成全部流程计算。PAGravf4.5建议采用残差迭代逼近方案,通过提高步骤(5)残差地面高程异常(大地水准面高)逼近性能,改善积分法重力大地水准面精化水平。

(9)调用[高程异常计算外部扰动场元],由步骤(5)的残差地面高程异常格网,计算等位边界面上残差扰动重力格网;将步骤(4)的残差扰动重力格网与这里的残差扰动重力格网相减,得到剩余残差扰动重力格网;再由剩余残差扰动重力格网,按步骤(5)计算剩余残差地面高程异常。

将地面高程异常的剩余残差值与步骤(5)中残差值相加,就是一次迭代逼近后的残差地面高程异常。一般情况下,1~2次迭代逼近能满足要求。

调节上述流程中的参数设置,使得重力地面高程异常(大地水准面高)与移去平均值后的GNSS水准实测地面高程异常(大地水准面高)尽量接近。

大地水准面径向基函数逼近简捷流程练习

✍练习目的:省略地形影响处理,由地面、海洋和航空重力和GNSS水准数据,按径向基函数逼近法精化区域大地水准面,以便快速了解大地水准面建模要领、效果与目标区域重力场及大地水准面特性,客观认识地形影响处理对大地水准面建模的作用。

省略地形影响处理后,由于不需要额外的延拓归算、格网化和GNSS水准融合流程,大地水准面径向基函数逼近流程非常简单。

✍计算面大地高格网

地面/海面大地高格网(= 数字高程模型 + 模型地面高程异常),用于地面高程异常建模;模型大地水准面高格网,用于大地水准面建模。

模型地面高程异常、模型大地水准面高用360阶重力场模型计算(略)。要求计算面大地高格网范围大于目标区域范围,以抑制边缘效应。

001.jpg

(1)多种异质离散残差观测场元文件构造。

调用[实际测点扰动重力场元计算],由重力数据文件,生成离散扰动重力观测文件;调用[点值属性加权运算];由GNSS水准数据文件,生成GNSS水准实测高程异常(或大地水准面高)文件。

调用[地球重力场各种场元模型值计算],最大计算阶数540,分别计算并移去观测扰动重力和实测GNSS水准的参考重力场模型值,按约定格式要求合并,生成多种异质离散残差观测场元文件。

(2)探测并剔除GNSS水准粗差。

调用[多源异质数据径向基函数重力场建模]程序,选择高程异常为可控观测场元,待估场元类型为高程异常,观测场元贡献率κ=0,输入计算面大地高格网,估计残差地面高程异常(大地水准面高)格网。可累积1次逼近。

从观测量残差点值文件*.chs中分离出GNSS水准残差记录,探测并剔除GNSS水准粗差点,移去GNSS水准残差的统计平均值后,重新按约定格式要求合并,生成新的多种异质离散残差观测场元文件。

统计平均值代表区域高程基准与全球高程基准(重力大地水准面)之间的系统偏差。

(3)残差重力地面高程异常(大地水准面高)SRBF逼近。

调用[多源异质数据径向基函数重力场建模]程序,选择高程异常为可控观测场元,贡献率κ不小于1,待估类型为高程异常,输入计算面大地高格网,估计残差地面高程异常(大地水准面高)格网。可累积1~2次逼近。

(4)恢复参考重力场模型值。调用[地球重力场各种场元模型值计算],最大计算阶数540,由计算面大地高格网,计算模型地面高程异常(大地水准面高)格网,并与残差地面高程异常(大地水准面高)格网相加,得到SRBF逼近法地面高程异常(大地水准面高)格网成果。

完成全部流程计算。计算过程输出的观测量残差点值文件*.chs,可用于粗差探测、质量控制和精度评估。

调节上述流程中的参数设置,使得地面高程异常(大地水准面高)与移去平均值后的GNSS水准实测地面高程异常(大地水准面高)尽量接近。

坐标框架中大地水准面与水准网的统一

 目的:

(1)通过GNSS水准网,将重力大地水准面纳入GNSS定位所在的坐标参考框架中,以满足GNSS代替水准的基本测量条件。

(2)利用GNSS水准实测高程异常高精度短波优势,抑制重力大地水准面短波超短波误差,改善重力大地水准面的局部相对精度。

(3)利用重力大地水准面高精度中长波优势,控制长距离水准传递的误差累积,实现坐标框架中高程基准面、水准控制网成果的统一。

方法:利用GNSS水准残差高程异常,评估大地水准面精度水平,精化高程基准基准零点重力位,计算基准差异。进而以GNSS水准残差高程异常之间的空间关系满足Possion积分方程为约束目标,通过调节积分半径大小来控制空间尺度,分离GNSS水准残差高程异常中满足Possion积分方程的成分,实现GNSS水准与重力大地水准面的解析融合。

(1)由GNSS水准高程异常成果与重力地面高程异常成果,生成离散GNSS水准残差高程异常文件。

①利用地面高程异常径向变化率格网模型,将GNSS水准高程异常,由GNSS水准点归算到重力地面高程异常所在的地面大地高格网模型面上。

②内插GNSS水准点处重力地面高程异常,由归算后的GNSS水准高程异常,计算GNSS水准残差高程异常。

③统计GNSS水准残差高程异常,对GNSS水准残差高程异常配权(初始可用等权),生成离散GNSS水准残差高程异常文件。

(2)调用“地面高程异常成果精度评估”,由GNSS水准残差高程异常,估计并确定重力地面高程异常差误差、实用地面高程异常内部误差、实用地面高程异常差误差曲线与GNSS水准高程异常差误差曲线。

(3)利用“高程基准零点重力位参数计算”,计算区域高程基准与重力大地水准面重力位差W-U,将GNSS水准残差地面高程异常减去基准差异量(W-U )/γ,得到基准统一后的GNSS水准残差地面高程异常文件。

(4)利用“泊松方程约束GNSS水准融合”,计算残差地面高程异常格网,生成剩余GNSS水准残差点值文件。

①统计地面大地高格网,确定格网中心点及平均大地高,调用“5.9过指定点模型等位面构造”,确定过格网中心、大地高为平均大地高的等位面大地高格网。

②调用“泊松方程约束GNSS水准融合”,由GNSS水准残差高程异常和等位面大地高,计算残差地面高程异常格网,生成剩余GNSS水准残差点值文件。

③重复调用“泊松方程约束GNSS水准融合”,将剩余GNSS水准残差作为GNSS水准残差高程异常,计算残差地面高程异常格网的迭代值,直到剩余残差量的标准差和平均值均无明显变化为止。

(5)由剩余GNSS水准残差点值文件,按3倍剩余残差标准差,探测并剔除GNSS水准粗差,分析剩余残差的空间分布性质。

(6)重复作业流程(2)~(5)1~2次(至多3次)后,将最后一次的残差地面高程异常格网迭代值相加,再与重力地面高程异常格网相加,并恢复高程基准差异量(W-U )/γ,生成坐标框架中地面高程异常格网模型。

(7)由剩余GNSS水准残差高程异常点值文件,调用“残差高程异常水准网拟稳平差”,改善水准网长距离高程传递质量。

①将GNSS水准点的正常高,与剩余GNSS水准残差高程异常相加,生成GNSS水准点正常高的新值。

②以全部GNSS水准点为拟稳基准,以流程(3)中最终的权值作为GNSS水准点正常高改正数(基准未知数)的权。

③以GNSS水准点为起算点,全部水准点正常高改正数为未知数,并将其中部分改正数作为基准未知数,由水准网测段文件,构建全部测段的观测方程。

④按附有基准未知数条件的间接平差方法,重新对水准网进行拟稳平差。

⑤平差后,地面高程异常格网纳入GNSS定位所在的坐标参考框架中,并与水准网高程和大地水准面成果实现高度统一。

(8)综合地面高程异常精度评估、泊松方程约束的GNSS水准融合以及水准网拟稳平差结果,评价大地水准面、GNSS水准网与高程基准成果的质量和精度水平。