2. 基于海区历史数据的水下滑翔机观测区域筛选

为什么不能直接拿历史数据画图就定区域

海区历史温盐流场数据(再分析产品、多年实测剖面、卫星高度计融合流场)本身只是"场",而滑翔机观测区域筛选要回答的是三个工程问题:哪里梯度大、值得加密观测;哪里流太强、会拖垮航程与续航;哪里必须避让、否则滑翔机可能被流带离任务区或撞上地形。因此本节把历史场数据加工成两类可判读的图层——观测价值层与航行风险层,再叠加得到优先区与避让区。

筛选指标体系与判读阈值

以某目标海区(示例设定为 2°×2° 矩形海域,中心 18°N、114°E,水深 800–4200 m)为例,从历史数据中提取 4 个指标,每个指标归一化到 0–1 后加权:

指标物理含义数据来源示例阈值
表层温度梯度 |∇SST|锋面/涡旋边界,生物活动与混合强多年月均 SST 场>0.08 ℃/km 记为高价值
温跃层深度方差 σ(Zt)层结不稳定,剖面信息量大历史 CTD 剖面统计>15 m 记为高价值
表层流速均值 |U|影响航程与航线保持能力高度计融合流场>0.35 m/s 记为高风险
流速空间切变 |∂U/∂n|强切变区易造成航线偏移与姿态扰动流场梯度计算>0.05 (m/s)/km 记为高风险

综合评分 Score = 0.35·G_T + 0.25·G_Z − 0.25·R_U − 0.15·R_S,其中 G 为价值项、R 为风险项(均归一化)。Score 高且风险项未超阈值的格点进入优先观测区;风险项超阈值且价值项不高的格点进入避让区;其余为常规观测区。

从格点场到区域结论:筛选流程

历史场数据 多年月均 SST 0.25° 网格 历史 CTD 剖面 温跃层统计 融合流场 表层流速 U/V 海底地形 水深/坡度 统一重采样至 0.1°×0.1° 分析格网 时间窗:近 10 年 同月(如 3 月) 指标计算 G_T:|∇SST| 中心差分求梯度 G_Z:σ(Zt) 逐格点剖面统计 R_U:|U| 均值 多年同月平均 R_S:|∂U/∂n| 法向切变 地形坡度 避让硬约束 各指标 min-max 归一化到 0–1 综合评分 Score = 0.35·G_T + 0.25·G_Z − 0.25·R_U − 0.15·R_S 价值项权重 0.60 风险项权重 0.40 硬约束:坡度 > 8° 或水深 < 200 m 直接判为避让 输出:格点评分栅格 区域判读 优先观测区 Score ≥ 0.55 且 R_U、R_S 未超阈 → 加密剖面/锯齿航线 航线避让区 |U| > 0.35 m/s 或切变 > 0.05 (m/s)/km 或坡度/水深硬约束 常规观测区 0.30 ≤ Score < 0.55 → 基础覆盖航线 筛选依据文档 指标图层 + 阈值表 + 分区多边形坐标 + 判读结论 阈值不合理时回退调整权重/阈值

关键计算:以两个格点为例

取分析格网中两个 0.1° 格点做示范(数值为教学示例,非实测):

  • 格点 A(锋面附近):|∇SST| = 0.12 ℃/km,σ(Zt) = 22 m,|U| = 0.18 m/s,切变 = 0.02 (m/s)/km。归一化后 G_T=0.80、G_Z=0.73、R_U=0.26、R_S=0.20,Score = 0.35×0.80 + 0.25×0.73 − 0.25×0.26 − 0.15×0.20 = 0.28 + 0.1825 − 0.065 − 0.03 = 0.3675。风险项均未超阈,但 Score 落在常规区——说明"梯度大"不等于"综合价值高",需结合层结与流场一起看。
  • 格点 B(涡旋西侧):|∇SST| = 0.15 ℃/km,σ(Zt) = 28 m,|U| = 0.22 m/s,切变 = 0.03 (m/s)/km。归一化后 G_T=0.95、G_Z=0.90、R_U=0.32、R_S=0.28,Score = 0.3325 + 0.225 − 0.08 − 0.042 = 0.4355,进入优先观测区候选。
  • 格点 C(西边界流轴):|U| = 0.42 m/s,切变 = 0.07 (m/s)/km,虽 G_T 高,但 R_U、R_S 均超阈,直接判为避让区,不参与优先区竞争。
判读要点:优先区看"价值高且风险可控",避让区看"任一风险项超阈即触发",两者不是同一评分轴上的两端,必须分开判定,否则会把高价值强流区误选为观测区。

区域筛选依据文档(可直接填写)

字段填写内容(示例)
目标海区范围17°N–19°N,113°E–115°E
分析时间窗近 10 年 3 月,0.1° 格网
指标与权重G_T 0.35 / G_Z 0.25 / R_U 0.25 / R_S 0.15
优先区判据Score ≥ 0.55 且 R_U ≤ 0.35 m/s、R_S ≤ 0.05 (m/s)/km
避让区判据|U| > 0.35 m/s 或切变 > 0.05 (m/s)/km 或坡度 > 8° 或水深 < 200 m
优先区多边形P1:18.2–18.6°N,113.4–113.9°E;P2:17.6–18.0°N,114.3–114.7°E
避让区多边形D1:18.7–19.0°N,113.0–113.5°E(西边界流轴)
结论与建议优先区 P1、P2 加密剖面;航线绕行 D1 东侧,保持 ≥0.15° 缓冲

结果验证与常见误判

筛选完成后至少做两项校验:一是阈值敏感性检查——把 Score 阈值在 0.50–0.60 间滑动,若优先区面积剧烈跳变,说明权重或阈值不稳,需回退调整(图中红色虚线即此反馈路径);二是与地形/边界叠加检查——确认优先区不落在坡度突变带或领海边界敏感区。常见误判是把"高梯度"直接等同于"优先区",忽略强流切变对航线保持的影响,导致滑翔机在观测途中被流推离目标格点,实际采样点偏离设计位置。