地球潮汐负荷影响与形变监测计算系统(Geodetic Computation for Earth Tide, Load Effects and Deformation Monitoring) ETideLoad4.0，一种地球物理大地测量监测科学计算的大型Windows程序包，主要由大地测量全要素各类潮汐影响计算，非潮汐大地测量时序处理分析，高分负荷形变场及时变重力场逼近计算，CORS/InSAR协同监测与地面稳定性估计，以及大地测量数据编辑计算与可视化五大子系统有机构成。 采用一致的地球物理模型和统一的数值标准，构造相容的大地测量与地球动力学算法，精确计算地面及固体地球外部各种大地测量的各类潮汐与非潮汐形变效应；统一多源异构大地测量时空监测基准和参考历元，构建其相互之间几何物理大地测量约束，精化多种监测量之间时空动力学关系，推进空天地海多种监测技术协同。为几何物理大地测量时空监测基准构建维持，多源异构地球监测数据深度融合，固体地球形变效应监测计算，地表水文环境、地质灾害与地面稳定性变化监测，提供一套科学实用的大地测量地球动力学算法系统。 |
（1）严格采用科学统一的标准与解析相容的算法，精密计算地面及固体地球外部空间点（泛指海洋、低空和卫星等不与地球固连的空间点）全要素大地测量观测量和参数的固体潮、海潮及大气潮负荷效应，非潮汐极移效应，永久潮汐影响与地球质心运动，实现地面全要素大地测量各类潮汐影响全球预报。
（2）计算大气、海平面、土壤水、江河湖库水和冰川冰盖雪山等地表环境负荷非潮汐变化导致的全球或区域非潮汐负荷形变场及时变重力场；综合多种大地测量监测量时序，基于固体地球形变理论，同化地表环境负荷观测数据，监测区域陆地水时空变化，精化高分负荷形变场及时变重力场时序。
（3）构建区域统一、长期稳定、高抗差性能的几何物理时空监测基准，开展CORS网、多源异构InSAR及多种大地测量协同监测与解析深度融合的科学计算；按客观自然规律，由大地测量几何物理场时序，构造地面稳定性降低定量辨识准则，实现地面稳定性变化空间无缝、时间持续的定量跟踪监测。
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ETideLoad采用自定义格式的五种形式大地测量时间序列数据，包括地面监测量时间序列、监测网站点记录时间序列、监测网观测量记录时间序列、监测量（向量）数值格网时间序列和球谐系数模型时间序列。
（1）地面监测量时间序列。一个地面监测量时序文件可存储某一站点、某条基线或路线多种属性监测量时间序列数据，各种属性的采样历元时刻完全相同。如CORS坐标解时间序列、固体潮站观测与分析成果时间序列、GNSS基线解时间序列等。
（2）监测网站点记录时间序列。一个监测网站点记录时序文件可存储大地监测网多个地面站点一种类型监测量的时间序列数据。如CORS站网坐标解（一维）、水准网点高程、固体潮站网、验潮站网时间序列，以及时序InSAR监测量等。
（3）监测网观测量记录时间序列。一个大地测量监测网观测量记录时序文件可存储CORS网（GNSS网）一维基线解、水准控制网水准路线高差观测量或重力控制网重力段差观测量时间序列数据。
（4）数字模型（向量）格网时间序列。由一组带时间（历元）信息的某一种类型监测量（向量）格网文件构成，每个格网头文件第7个属性约定为该格网数据的采样历元时刻。如陆地水等效水高、海平面变化格网时序，各种地面负荷形变场格网时序，时变重力场格网时间序列等。
（5）球谐系数模型时间序列。由一组带时间信息的球谐系数（Stokes系数）模型文件构成，用于表示全球负荷变化时间序列、全球负荷形变场时间序列和全球重力场变化时间序列。每个文件头占住一行，至少包含地心引力常数GM（×10¹⁴m³/s²）、地球椭球长半轴a（m）和采样历元时刻（ETideLoad约定的长整数格式）3个属性。GM,a称为球谐系数模型的尺度参数。
因技术需要，ETideLoad约定了大地测量监测量概念。
☆ ETideLoad4.0自定义
大地测量监测量定义为当前历元时刻的大地测量观测量或参数，与其一段时间内观测量或参数的平均值或某一参考历元时刻大地测量观测量或参数之差。大地测量监测量通常用大地测量观测量或参数变化（特指时间差分）表示。如地面重力变化、地倾斜向量变化分别表示地面重力监测量和地倾斜监测量。
☆ 监测量类型及单位约定
（1）高程异常/大地水准面变化（mm，毫米）、地面重力/扰动重力变化（μGal，微伽）、地倾斜/垂线偏差变化（mas，0.001ʺ/毫角秒）；
（2）地面水平位移（东向/北向，mm）、地面径向位移（大地高变化mm）、地面正（常）高变化（mm）；
（3）扰动重力梯度变化（10μE）、水平重力梯度变化（北向/东向10μE）；
（4）卫星轨道重力位摄动（0.1m²/s²）、重力摄动及其三维分量变化（μGal）、扰动重力梯度摄动对角线三分量变化（10μE）；
（5）陆地水负荷等效水高变化（cm）、海平面变化（cm）、海洋潮高（cm）和大气压变化（hPa）。
☆ 监测量向量方向约定
（1）地倾斜/垂线偏差变化SW。第一分量指向南，第二分量指向西，与地面重力/扰动重力方向构成右手直角坐标系，即自然坐标系；
（2）站点水平位移EN。第一分量指向东，第二分量指向北，与径向位移（大地高变化）方向构成右手直角坐标系，即东北天坐标系；
（3）水平重力梯度变化NE。第一分量指向北，第二分量指向东，与扰动重力梯度方向（铅垂线方向）构成右手直角坐标系；
（4）分潮调和常数。第一分量为同相幅值（余弦分量），第二分量为异相幅值（正弦分量）。
 功能：按输入时序文件中地点和时刻，计算地面各种大地测量参数的固体潮影响。 可供选择的大地测量参数类型包括：高程异常/大地水准面（mm）、扰动重力（μGal）、地倾斜北南分量（ms）、地倾斜东西分量（ms），地面重力（μGal）、地面东方向（mm）、地面北方向（mm）、地面径向（大地高，mm）或地面正常高（mm）。 输入：待计算的时间序列文件。 输出：计算结果时间序列文件（如图2），在原时间序列文件记录的基础上增加一列或若干列固体潮影响时间序列值，保留两位有效数字。 图2~3为重庆万州某地面站点3天的固体潮影响曲线。图中显示，正常高固体潮影响（最大幅值约300mm）与大地高固体潮影响、大地水准面固体潮影响（最大幅值约600mm）异相（即符号相反）。
输入：待计算的时间序列文件。程序要求头文件中计算点的高程为正（常）高。 参数设置：输入球谐系数最大计算阶数。程序自动选择海潮负荷球谐系数模型最大阶数和输入最大阶数中的最小值作为计算阶数。 （1）计算速度取决于海潮球谐系数模型阶数、分潮数量和时序长度。 （2）图2为重庆万州某地5天海潮负荷影响。与固体潮影响不同，正常高与大地高负荷影响同相，正常高海潮负荷影响的幅值约为大地高的1.75倍。 （3）图3、4为浙江温州沿海某站点5天的海潮负荷影响曲线。比较图2、图3可以看出，浙江温州沿海地区海潮负荷影响是重庆万州内陆的5倍以上。
|
