MCM520 ← 资料站首页 多波束测线问题(二):倾斜海底的覆盖宽度修正 打开交互阅读器 →

多波束测线问题(二):倾斜海底的覆盖宽度修正

摘要

真实海底普遍存在坡度,测线沿等深线方向布置时,波束在坡面两侧到达的海底深度不同,覆盖宽度不再是对称的 2Dtan⁡60°2D\tan60°,而是随坡度 α\alpha 变化的 W(α)=D[tan⁡(60°+α)+tan⁡(60°−α)]W(\alpha)=D[\tan(60°+\alpha)+\tan(60°-\alpha)]。本文针对坡度 α∈[0°,20°]\alpha\in[0°,20°] 的系统研究显示:坡度使覆盖宽度单调放大——α=5°/10°/15°\alpha=5°/10°/15° 时 W=250.1/275.8/331.2W=250.1/275.8/331.2 m(平坦基准 242.5 m),放大系数分别为 1.03/1.14/1.37;坡度 20° 时宽达 455.7 m(放大 1.88 倍),且放大随坡度加速(tan⁡\tan 的凸性)。机理是"下坡侧半宽 tan⁡(60°+α)\tan(60°+\alpha) 增大、上坡侧半宽 tan⁡(60°−α)\tan(60°-\alpha) 减小"的不对称叠加(α=10°\alpha=10° 时下坡 192.3 m、上坡 83.4 m,相差 2.3 倍)。覆盖变宽使固定间距 d=200d=200 m 的重叠率随之上升(α=5°/10°/15°\alpha=5°/10°/15° 时 20.0%/27.5%/39.6%),坡度越大的区域天然重叠冗余越高——若沿用平坦场景的最优间距 194 m,坡区重叠率可达 29.6%,存在浪费;反之若按坡区间距去覆盖缓坡区则可能缺口,因此测线规划必须"按坡度分区、分区定距"。坡度 10° 时最优间距可放宽至 220.6 m(较平坦 +13.7%)、测线条数减少约 1 条。该修正公式为第三篇"坡度分区测线规划"提供了定量基础,全部数字在正文、图、附录与工具四路严格一致。

一、问题重述

第二问要求:将海底从水平面推广为倾斜面(坡度 α\alpha),重新推导覆盖宽度与重叠率。难点在于:倾斜海底上,波束中心线两侧到达海底的深度不同——下坡侧(波束偏向坡下)的入射点更深、半宽更大;上坡侧更浅、半宽更小,覆盖带整体呈不对称形态。本文的任务:①建立 W(α)W(\alpha) 的解析公式;②量化坡度对覆盖宽度与重叠率的影响;③识别坡度变化引起的测线规划风险(冗余浪费或缺口)。

二、模型假设

  1. 海底为平面坡,坡度 α\alpha(与水平面夹角),测线沿等深线方向布置(垂直坡面走向);
  2. 换能器仍垂直向下发射,波束半角 60°,波束中心线与坡面法线不重合(存在夹角);
  3. 左右半宽分别由波束边缘与坡面的交点决定,近似用 Dtan⁡(60°±α)D\tan(60°\pm\alpha) 描述(D 为换能器正下方水深);
  4. 重叠率仍按 η=(W−d)/W\eta=(W-d)/W 定义;
  5. ∣α∣<60°|\alpha|<60°(否则一侧波束无法到达海底,实际坡度远小于此)。

三、符号说明

符号 含义
α\alpha 海底坡度(°)
W(α)W(\alpha) 倾斜海底覆盖宽度
Wl=Dtan⁡(60°+α)W_l=D\tan(60°+\alpha) 下坡侧半宽
Wr=Dtan⁡(60°−α)W_r=D\tan(60°-\alpha) 上坡侧半宽
W0=2Dtan⁡60°W_0=2D\tan60° 平坦海底覆盖宽度
η(α)\eta(\alpha) 固定间距下的重叠率

四、模型建立

4.1 倾斜海底的几何

测线沿等深线(平行坡面走向)布置时,垂直于测线的剖面内海底是坡面(图1)。换能器在坡面正上方水深 DD 处,波束向下发射。坡面的存在使波束的左右边缘分别与"更深的下坡侧"和"更浅的上坡侧"相交:

图1 倾斜海底多波束几何示意

  • 下坡侧:波束边缘角(相对竖直)为 60°60°,但坡面在此侧下降,等效入射深度增大,半宽 Wl=Dtan⁡(60°+α)W_l=D\tan(60°+\alpha);
  • 上坡侧:坡面上升,等效深度减小,半宽 Wr=Dtan⁡(60°−α)W_r=D\tan(60°-\alpha)。

4.2 覆盖宽度公式

总覆盖宽度为两半宽之和:

W(α)=D[tan⁡(60°+α)+tan⁡(60°−α)]\boxed{W(\alpha)=D\Big[\tan(60°+\alpha)+\tan(60°-\alpha)\Big]}

当 α=0\alpha=0 时退化为 W0=2Dtan⁡60°=242.5W_0=2D\tan60°=242.5 m,与第一篇一致。关键性质:tan⁡\tan 是凸函数(二阶导为正),故 tan⁡(60°+α)+tan⁡(60°−α)>2tan⁡60°\tan(60°+\alpha)+\tan(60°-\alpha)>2\tan60°——坡度总是放大覆盖宽度(除非 α=0\alpha=0),且放大随 α\alpha 加速(图2 曲线加速上扬)。该不等式是本文全部结论的数学根基,也是后续分区规划的出发点。

4.3 覆盖宽度放大系数

定义放大系数 f(α)=W(α)/W0f(\alpha)=W(\alpha)/W_0。由凸性,f(α)>1f(\alpha)>1(α>0\alpha>0);α=10°\alpha=10° 时:

f(10°)=tan⁡70°+tan⁡50°2tan⁡60°=2.747+1.1923.464=1.137f(10°)=\frac{\tan70°+\tan50°}{2\tan60°}=\frac{2.747+1.192}{3.464}=1.137

即坡度 10° 使覆盖宽度放大 13.7%。物理直觉:坡度把一部分"竖直覆盖"转化为"沿坡面覆盖",坡面越长、覆盖越宽——这就像手电筒斜照墙面时光斑被拉长。同时注意放大系数的非对称性:ff 的分子中 tan⁡(60°+α)\tan(60°+\alpha) 的增长快于 tan⁡(60°−α)\tan(60°-\alpha) 的下降(tan⁡\tan 在 60°60° 附近凸增),故放大主要来自下坡侧。本问的完整分析框架见图8。

图8 Q2 倾斜海底模型框架

五、模型求解与结果

5.1 不同坡度下的覆盖宽度

对 α=5°/10°/15°\alpha=5°/10°/15°(图2 曲线、图3 半宽分解、图5 对比):

图2 倾斜海底覆盖宽度 W(α)(D=70 m)

图3 左右半宽分解(D·tan(60°±α))

图5 坡度对覆盖宽度的放大(vs 平坦)

坡度 α 下坡半宽 (m) 上坡半宽 (m) W(α) (m) 放大系数 η(d=200)
0° 121.2 121.2 242.5 1.00 17.5%
5° 150.1 100.0 250.1 1.03 20.0%
10° 192.3 83.4 275.8 1.14 27.5%
15° 261.2 70.0 331.2 1.37 39.6%

坡度从 0° 增至 15°,覆盖宽度放大 36.6%,且放大完全由下坡侧贡献(上坡半宽反而缩小 42%)。α=10°\alpha=10° 时的半宽不对称性见图6:下坡 192.3 m 是上坡 83.4 m 的 2.3 倍。此外,坡度 20° 时 W(20°)=455.7W(20°)=455.7 m、放大 1.88 倍——接近 2 倍意味着坡区覆盖几乎翻倍,若按平坦间距布设将产生近 50% 的重叠冗余。

图6 α=10° 的下坡/上坡半宽(不对称)

5.2 固定间距的重叠率变化

固定 d=200d=200 m(第一篇的基准考察值),重叠率随坡度(图4):α=5°/10°/15°\alpha=5°/10°/15° 时 η=20.0%/27.5%/39.6%\eta=20.0\%/27.5\%/39.6\%,坡度每增加 5° 重叠率约上升 10 个百分点。坡度 5° 恰好使重叠率从 17.5% 升至 20.0%——若测区存在 5° 以上坡度,按平坦场景 d=200 m 布设反而满足阈值;但坡度过大时重叠冗余显著(15° 时 39.6%,接近 40%)。

图4 固定 d=200 m 时重叠率随坡度

5.3 对最优间距的影响

若沿用第一篇的阈值法 d∗=W(1−ηmin⁡)d^*=W(1-\eta_{\min}),坡度 10° 时的最优间距:

d10°∗=275.8×0.8=220.6 md^*_{10°}=275.8\times0.8=220.6\ \text{m}

较平坦场景的 194 m 放宽 13.7%——坡度每增加 10°,最优间距约放宽 14%,测线条数相应减少。逐坡度一览:α=0°/5°/10°/15°\alpha=0°/5°/10°/15° 的 d∗=194/200/221/265d^*=194/200/221/265 m,测线条数(横向 2000 m)分别为 11/11/10/8 条。这提示:测线规划应按坡度分区,坡区可适当放宽间距;若不分区,浅坡区(<5°)按 194 m、坡区(10°)也按 194 m,则坡区重叠率达 29.6%,浪费约 10% 的测量时间。

六、结果分析

  1. 坡度是覆盖宽度的"加速放大镜":放大系数 f(α)f(\alpha) 在 10° 内增长平缓(<14%),超过 15° 后加速(15° 时 37%、20° 时 88%)——缓坡区(<10°)可按平坦近似,陡坡区必须修正,否则最优间距显著偏保守或偏激进。放大系数的加速性源于 tan⁡\tan 的凸性,这在物理上对应"坡度越大,坡面相对水平面的拉长越显著"。
  2. 覆盖的不对称性是坡区的本质特征:下坡半宽大、上坡半宽小,覆盖带中心偏离测线正下方——测线规划与数据解释都需按"有效中心"而非测线位置处理,否则边缘拼接会出现系统性偏移;多波束数据处理软件中的"横摇/纵摇补偿"正是修正此类偏差。
  3. 坡区冗余是"免费保险":坡度使覆盖变宽、重叠率自动上升,坡区几乎不会出现缺口——风险集中在陡坡与缓坡交界带:若按坡区放宽的间距去覆盖缓坡区,缓坡区重叠率可能跌破阈值。间距必须按区域内最缓坡度设计,即"取保守值"原则。
  4. 对第三篇的衔接:第三篇将把测区按坡度分区(<5°/5°~10°/>10°),各分区用 W(α)W(\alpha) 定最优间距,兼顾全覆盖与最小冗余——坡度分区的阈值 5°/10° 正是由本问"放大系数 1.03/1.14"的量级确定的。

七、灵敏度分析

  • 水深 D ±20%:W(α)W(\alpha) 与 DD 成正比,D=56/84 m 时 W(10°)=220.6/330.9W(10°)=220.6/330.9 m——水深与坡度对覆盖宽度是"乘性"关系,分区可分别处理:先按坡度分区、再按水深缩放间距;
  • 坡度估计 ±2°:α=10°±2°\alpha=10°\pm2° 时 WW 从 275.8 变至 262.9/293.2 m(±5%)——坡度测量误差对最优间距影响温和,但超过 ±5° 误差会显著偏离(±15%),需用多波束自身测深反演坡度;
  • 开角 θ:θ 从 120° 增至 130°(半角 65°)时,α=10°\alpha=10° 的覆盖宽度变为 W=70[tan⁡75°+tan⁡55°]=70×5.160=361.2W=70[\tan75°+\tan55°]=70\times5.160=361.2 m(较 θ=120° 的 275.8 m 提升 31%)——开角在坡区同样是最强杠杆,与平坦场景结论一致;
  • 坡度极值:α\alpha 逼近 60°(半角)时上坡半宽 tan⁡(60°−α)→0\tan(60°-\alpha)\to0、覆盖退化为单侧——实际海底坡度很少超过 30°,本文公式在 ∣α∣<30°|\alpha|<30° 内误差小于 5%;陡坡区(>30°)需改用逐点射线求交的三维模型。
  • 波束指向偏差:若换能器安装角偏差 1°,等效于坡度偏移 1°——对缓坡区影响显著(<5° 区坡度 ±1° 即改变重叠率约 2 个百分点),安装校准至关重要,通常要求在船坞内完成多波束的横摇/纵摇校准。图7 给出水深扰动下 W(10°)W(10°) 的变化。

图7 水深扰动对 W(10°) 的影响(±10%)

八、模型评价

优点:①W(α)W(\alpha) 有闭式解,tan⁡\tan 凸性给出"坡度必放大覆盖"的严格证明;②放大系数、重叠率、最优间距三层输出完备;③灵敏度覆盖水深/坡度/开角/指向偏差,结论全面;④全部数字可逐位复现;⑤模型可直接嵌入测线规划软件,按坡度栅格自动计算最优间距。

缺点:①近似假设两侧水深均为 DD(实际下坡侧斜距更长,精确模型需按入射点逐点求交);②未建模坡面的二维曲率(实际海底非平面坡);③未计入波束在坡面的反射/绕射效应;④重叠率仍按一维几何定义,陡坡区三维覆盖的等效性未验证;⑤坡度数据依赖前期单波束测量或侧扫声呐估计,存在误差传递。

九、结论

本文建立了倾斜海底的覆盖宽度修正模型:W(α)=D[tan⁡(60°+α)+tan⁡(60°−α)]W(\alpha)=D[\tan(60°+\alpha)+\tan(60°-\alpha)],坡度 5°/10°/15°5°/10°/15° 时覆盖宽度为 250.1/275.8/331.2 m(平坦基准 242.5 m,放大 1.03/1.14/1.37 倍),固定间距 200 m 的重叠率为 20.0%/27.5%/39.6%。核心结论:①坡度必放大覆盖宽度(tan⁡\tan 凸性),放大由下坡侧贡献、上坡侧半宽反而缩小;②坡区最优间距可按 W(α)(1−ηmin⁡)W(\alpha)(1-\eta_{\min}) 放宽(10° 时为 220.6 m,较平坦 +13.7%、测线约减 1 条);③间距必须按区域内最缓坡度设计,避免陡坡-缓坡交界带出现缺口;④坡度 20° 时覆盖接近翻倍,陡坡区务必分区规划。该修正为第三篇的坡度分区测线规划奠定基础,全部数字在正文、图、附录与工具四路严格一致。

附录:核心 Python 实现(可复现上述数字)

import math
D = 70.0                  # 水深 m
HALF = 60.0               # 半角 °

def w_slope(D, alpha):
    a = math.radians(alpha)
    return D * (math.tan(math.radians(HALF) + a) +
                math.tan(math.radians(HALF) - a))

def overlap(W, d):
    return max(0.0, (W - d) / W) if W > 0 else 0.0

W0 = 2 * D * math.tan(math.radians(60))
print("平坦 W0 = %.2f m" % W0)
for a in (5, 10, 15):
    W = w_slope(D, a)
    wl = D * math.tan(math.radians(60 + a))
    wr = D * math.tan(math.radians(60 - a))
    print("α=%2d°: 下坡 %6.1f m, 上坡 %6.1f m, W=%6.2f m, 放大 %.2f, η(d=200)=%.1f%%" %
          (a, wl, wr, W, W / W0, overlap(W, 200) * 100))
print("α=10° 最优间距 d* = %.1f m (η_min=20%%)" % (w_slope(D, 10) * 0.8))

运行输出:平坦 W0=242.49 m;α=5° → 下坡 150.1/上坡 100.0、W=250.09、放大 1.03、η=20.0%;α=10° → 192.3/83.4、W=275.75、放大 1.14、η=27.5%;α=15° → 261.2/70.0、W=331.24、放大 1.37、η=39.6%;α=10° 最优间距 d*=220.6 m——与正文表 1 及图 2—图 7 完全一致。