2. 风场基础

无人机在真实大气中飞行时,所遭遇的风并不是一个恒定不变的矢量,而是随高度、时间与空间位置不断变化的复杂场。要建立可用的风场估计与抗风控制方法,首先必须理解大气边界层、湍流、阵风、风切变,以及平均风与脉动风的基本概念,并掌握几种典型的工程化风场模型。本节从物理结构出发,逐步过渡到可用于仿真与控制器设计的数学描述。

一、大气边界层:风场的第一层结构

地球表面与自由大气之间,存在一层受地面摩擦显著影响的大气区域,称为大气边界层(Atmospheric Boundary Layer, ABL)。其厚度通常在白天对流条件下可达 1~2 km,夜间稳定层结时可能只有几十到几百米。无人机绝大多数任务高度在 0~500 m 范围内,因此几乎始终运行在边界层内部。

边界层内风速随高度增加而增大,风向也会随高度发生偏转(埃克曼螺线)。按热力稳定度,边界层可分为对流边界层、中性边界层和稳定边界层。中性条件下,风速廓线最接近对数律;对流条件下,热浮力产生强湍流混合;稳定条件下,湍流被抑制,风切变和低空急流更突出。对无人机而言,这意味着同一地点不同时刻的风场模型参数可能完全不同。

二、平均风与脉动风:风场的分解视角

工程上通常把风速分解为两部分:

  • 平均风:在给定时间窗口(如 10 min 或 1 min)内的统计平均值,决定无人机的位置漂移与配平迎角。
  • 脉动风:瞬时风速与平均风之差,反映湍流引起的随机扰动,决定姿态抖动、控制带宽与能量消耗。

用数学表达即 u(t) = U + u'(t),其中 U 为平均风,u'(t) 为零均值的脉动分量。平均风是控制器前馈与轨迹修正的主要依据,脉动风则通过反馈与鲁棒性设计来抑制。二者在频域上的分离,是风场估计滤波器设计的核心出发点。

平均风决定“飞到哪里”,脉动风决定“飞得稳不稳”。抗风控制必须同时处理这两个尺度的问题。

三、湍流、阵风与风切变

湍流是边界层内不规则、多尺度的涡旋运动,具有随机性、耗散性和宽频带特征。湍流强度常用 TI = σ_u / U 表示,其中 σ_u 为脉动风速标准差。低空湍流强度可达 0.1~0.3,对小型无人机影响显著。

阵风是短时间内风速的突然变化,通常用阵风因子 G = u_max / U 描述。阵风具有明显的方向性和持续时间,常见模型有离散阵风(1-cos 型)和连续阵风(Dryden、von Kármán 谱)。阵风是导致无人机瞬时姿态偏差和位置超调的主要外因。

风切变是风速或风向在空间上的急剧变化,包括垂直风切变和水平风切变。低空风切变对起降阶段威胁最大,可能造成升力突变和操纵失效。垂直风切变常用幂律或对数律描述:

幂律模型:U(z) = U_ref * (z / z_ref)^α
对数律模型:U(z) = (u* / κ) * ln(z / z0)

其中:
α   —— 风切变指数,开阔地约 0.14,城市约 0.3
u*  —— 摩擦速度
κ   —— 卡门常数,约 0.4
z0  —— 地表粗糙度长度

四、典型风场模型与参数对照

在仿真与控制器设计中,常用的风场模型包括平均风廓线模型、Dryden 湍流模型、von Kármán 湍流模型、离散阵风模型和风切变模型。它们分别描述不同尺度的风现象,通常组合使用。

模型 描述对象 关键参数 典型用途
幂律/对数律 平均风垂直廓线 α、z0、U_ref 航迹规划、前馈补偿
Dryden 谱 连续湍流脉动 σ_u、L_u、L_w 频域控制器设计
von Kármán 谱 连续湍流脉动 σ、L、频率 高保真仿真
1-cos 离散阵风 单次阵风事件 幅值、周期 时域响应测试
风切变模型 空间风速梯度 切变率、方向 起降阶段鲁棒性验证

五、案例:小型多旋翼在 50 m 高度的风场组合

假设某多旋翼在开阔地 50 m 高度悬停,地面 10 m 处平均风速 4 m/s,风切变指数 α=0.14,湍流强度 TI=0.15。则 50 m 处平均风约为:

U(50) = 4 * (50 / 10)^0.14 ≈ 4 * 1.258 ≈ 5.03 m/s
脉动标准差 σ_u ≈ 0.15 * 5.03 ≈ 0.75 m/s

若叠加幅值 3 m/s、周期 4 s 的 1-cos 阵风,则瞬时风速峰值可达约 8 m/s。此时控制器若仅按平均风配平,姿态角会出现明显波动,位置误差可能超过 1 m。合理的做法是:用平均风做前馈倾角补偿,用 Dryden 谱设计陷波与带宽匹配,用阵风模型验证瞬态超调。

在风场估计中,平均风与脉动风的分离窗口不宜过短,否则会把阵风误判为平均风;也不宜过长,否则会丢失风切变信息。工程上常取 1~10 s 滑动窗口,并结合 IMU 与 GPS 数据做互补滤波。

六、常见误区与复盘

  • 误区一:把平均风当作恒定风。实际平均风随时间缓慢变化,长时间任务需在线更新。
  • 误区二:忽略风向脉动。风向变化会改变无人机相对气流角,影响升力与侧滑。
  • 误区三:只用一种湍流模型。Dryden 与 von Kármán 在不同频段表现不同,应交叉验证。
  • 误区四:风切变只在起降阶段考虑。山地、城市峡谷中巡航阶段同样可能遭遇强切变。

复盘时应重点检查:风场模型参数是否与实际气象条件匹配、估计窗口是否合理、控制器是否在阵风与切变组合下仍保持稳定。只有把平均风、脉动风、阵风和风切变统一到同一套风场描述框架中,后续的风场估计与抗风控制才有可靠的基础。

大气边界层风场结构示意图 地面 / 地表粗糙度 z₀ 大气边界层 (ABL) 自由大气(无地面摩擦影响) 高度 z 平均风廓线 U(z) 脉动风 u'(t) 无人机 风切变 ∂U/∂z 阵风 G = u_max/U 1-cos / Dryden 湍流 TI = σ_u/U 多尺度涡旋 u(t) = U + u'(t) | U(z) = U_ref·(z/z_ref)^α | TI = σ_u/U | G = u_max/U