从零到一:为Bernese5.2高精度解算定制专属的潮汐改正文件
1. 为什么需要定制潮汐改正文件第一次用Bernese5.2处理非IGS测站数据时我发现解算结果总是比预期差2-3毫米。排查了半天才发现原来软件自带的EXAMPLE.BLQ和EXAMPLE.ATL文件里只包含少数几个示范测站的潮汐参数。这就好比用别人的眼镜看世界——度数不对画面自然模糊。潮汐改正主要包含两大类型海洋潮汐加载Ocean Tide Loading海水质量变化导致的地壳形变常用FES系列模型计算大气潮汐加载Atmospheric Tide Loading大气压力变化引起的地表位移通过S1/S2潮汐分量建模实测对比显示对于沿海站点忽略潮汐改正可能导致高程方向出现8-15mm的偏差。即便是内陆站点影响也可能达到3-5mm。这个量级对高精度GNSS解算来说已经足够让结果从优秀变成勉强可用。2. 海洋潮汐文件.BLQ定制指南2.1 模型选择FES2004还是FES2014b在Chalmers大学提供的在线工具里我发现模型下拉菜单有十多个选项。经过实测对比模型版本适用场景计算点位数限制典型耗时FES2004全球通用单次≤50个点3-5分钟FES2014b近海区域精度更高单次≤20个点5-8分钟EOT20侧重欧洲区域无明确限制即时返回建议优先选择FES2014b它的空间分辨率从0.5°提升到0.125°特别适合东亚这种潮汐复杂的海域。不过要注意计算超过20个站点时需要分批提交。2.2 分步操作指南访问计算工具打开http://holt.oso.chalmers.se/loading注意需使用学术网络参数设置模型选择FES2014b输出格式BLQ默认参考框架保持ITRF2014输入坐标# 格式站名 经度(度) 纬度(度) 高程(m) BJFS 116.1954 39.6086 87.5 WUHN 114.3574 30.5315 25.3提交计算填写有效邮箱后点击Submit通常5分钟内会收到含附件的邮件收到邮件后用文本编辑器打开附件你会看到类似这样的内容BJFS FES2014b 11 .00000 .00000 .00000 .00000 .00000 .00000 .00000 .00000 .00000 .00061 .00102 .00048 .00022 .00005 .00001 .00000 .00000 .00000 ...将这段内容追加到EXAMPLE.BLQ文件末尾即可。记得每个站点的数据块要以站名和模型标识开头。3. 大气潮汐文件.ATL生成技巧3.1 通过Bernese GUI生成在Service → Coordinate Tools菜单里选择Extract atmospheric tidal loading cofficients功能时有几个关键参数需要注意潮汐格网文件建议使用最新的s1_s2_def_ce.dat软件安装包内提供输入坐标文件确保是包含测站ECEF坐标的CRD文件质心修正采用默认的RAY_PONTE模型即可我遇到过的一个典型报错是Grid file not found这通常是因为格网文件路径包含中文或空格文件权限设置问题Linux/Mac下需chmod 755文件编码格式非ASCII3.2 手动编辑ATL文件对于需要特殊处理的站点可以直接编辑ATL文件。其结构如下WUHN 1.234e-3 5.678e-4 ... [12个系数]每个站点占一行包含站名必须与CRD文件完全一致12个S1/S2潮汐系数顺序不能错单位米曾经有个项目需要处理北极站点发现自动生成的系数异常。后来手动参考了IERS Technical Note 36中的修正公式才得到合理结果。4. 实战问题排查手册4.1 常见错误代码解析错误提示可能原因解决方案BLQ_READ_ERR文件格式错误检查站名是否含空格ATL_MISSING_STA坐标文件与ATL不匹配确认CRD中站点名称一致TIDE_MODEL_MISMATCH模型版本冲突统一使用FES2014b4.2 精度验证方法建议按这个流程验证改正效果先不做潮汐改正解算一次添加BLQ文件后解算再添加ATL文件解算对比三次解算的NEU坐标差异以我处理的某沿海CORS站为例无改正时高程方向周期性波动达±12mm仅加BLQ后波动减至±5mm两者都用后稳定在±2mm内5. 高级应用场景5.1 动态潮汐改正对于长时间跨度解算可以考虑按季节生成多组BLQ文件在BPE脚本中添加条件判断IF ${YEAR} 2010 THEN SET V_BLQINF FES2004 ELSE SET V_BLQINF FES2014b ENDIF5.2 混合模型策略遇到某些近海站点时我发现深海区域用FES2014b浅海区域用EOT20陆地部分用ATLAS 通过脚本自动合并不同模型的BLQ数据最终平面精度提升了约30%最后提醒一点每次更新潮汐文件后建议清空WORK目录下的缓存文件否则可能读到旧数据。这个坑我踩过三次每次都要浪费半天时间排查。