4. 地球模型与重力场：地球形状与重力模型、地球自转与哥里奥利力、垂线偏差与重力异常

做惯导这么多年，我越来越觉得一个道理：导航精度，很大程度上取决于你对地球的理解有多深。你想想看，我们所有的加速度计和陀螺仪，都是在“地球”这个参考系里工作的。如果连地球长什么样、怎么转、重力怎么分布都搞不清楚，那后面的计算全是空中楼阁。

这一节，我们就来啃啃这块硬骨头。别怕，我会用我踩过的坑，帮你把路铺平。

4.1 地球形状：不是球，也不是椭球

小时候我们都学过“地球是圆的”。但做惯导的人要是真把它当球，那可就麻烦了。实际上，地球更像个梨形——赤道鼓起来，两极扁下去。

为什么？因为地球自转。赤道上的物质离心力大，被甩出去了。这个“扁”的程度，专业术语叫扁率。

我个人习惯用 WGS-84 模型。这是目前全球卫星导航系统（GPS）用的标准。它把地球描述成一个旋转椭球体。说白了，就是一个椭圆绕短轴转一圈得到的形状。

嗯，这里要注意：大地水准面和参考椭球面不是一回事。大地水准面是重力等势面，它更接近真实海平面。而参考椭球面是我们为了计算方便，人为拟合出来的光滑曲面。两者之间的差距，就是后面要讲的“垂线偏差”和“重力异常”。

4.2 重力模型：g 不是常数

很多初学者以为重力加速度 g 就是 9.8 m/s²。我刚开始做项目时也这么想，结果在高原地区测试，导航误差大得离谱。

实际上，g 随纬度和高度变化。公式长这样：

g(φ, h) = g₀(φ) * (1 - 2h/R₀)

其中：

• g₀(φ) 是海平面上的重力，随纬度变化
• h 是海拔高度
• R₀ 是地球平均半径

更精确的 正常重力公式（索米里安公式）是：

g₀(φ) = 9.780327 * (1 + 0.0053024 * sin²φ - 0.0000058 * sin²2φ)

你看，在赤道（φ=0°）和两极（φ=90°），g 能差出 0.05 m/s² 左右。这个量级，对于高精度惯导来说，绝对不能忽略。

4.3 地球自转与哥里奥利力

地球在转，我们也在转。这个“转”对惯导的影响，主要体现在两个地方：

1. 地球自转角速度：ωₑ = 7.292115 × 10⁻⁵ rad/s
2. 哥里奥利力：运动物体在旋转参考系中受到的“假想力”

哥里奥利力的公式很简单：

F_c = -2m (ω × v)

其中 ω 是地球自转角速度矢量，v 是物体相对地球的速度。

这个力有什么用？说白了，它会让北半球运动的物体向右偏，南半球向左偏。导弹、炮弹、甚至洋流都受它影响。

你想想看，如果飞机以 900 km/h 的速度飞行，哥里奥利加速度能达到 0.03 m/s²。这个量级，足以让惯导系统在几分钟内偏离航线几百米。

4.4 垂线偏差与重力异常

这两个概念，是惯导领域里最容易让人头疼的。我尽量说清楚。

垂线偏差：真实重力方向（铅垂线方向）与参考椭球法线方向之间的夹角。说白了，就是“你以为的垂直”和“真正的垂直”之间的偏差。

重力异常：真实重力值与正常重力值之间的差值。

为什么会这样？因为地球内部质量分布不均匀。有的地方密度大（比如铁矿区），重力就大；有的地方密度小（比如海洋盆地），重力就小。

这两个量对惯导的影响，我列个表：

误差源	典型量级	对惯导的影响
垂线偏差	1″ ~ 30″（角秒）	引起水平姿态误差，进而导致位置误差
重力异常	10 ~ 100 mGal	引起垂直通道发散，影响高度精度

1 mGal = 10⁻⁵ m/s²。你看，这个量级虽然小，但对于高精度惯导（比如潜艇惯导），是必须补偿的。

我记得有一次做船载惯导测试，在某个海域，垂线偏差达到了 15″。如果不补偿，惯导系统给出的水平位置误差，在航行 1 小时后会超过 500 米。后来我们用了当地的重力图，才把精度拉回来。

4.5 知识体系总览

这一节内容比较多，我画了张图，帮你理清思路：

这张图把这一节的核心逻辑串起来了。你看，从地球形状出发，引出重力模型；从地球自转出发，引出哥里奥利力；最后，两者交汇到垂线偏差和重力异常——这是惯导系统必须面对的现实误差源。

做惯导，说白了就是跟这些“不完美”打交道。地球不是完美的球，重力不是完美的常数，自转也不是完美的匀速。但正是这些不完美，才让我们的工作有了价值。

---