题目展示
单波束测探是用声波在水中的传播特性来测量水体深度的技术。声波在均匀介质中作匀速直线传播,在不同界面上产生反射,利用这一原理,从测量船换能器垂直向海底发射声波信号,并记录从声波发射到信号接收的传播时间,通过声波在海水中的传播速度和传播时间计算出海水的深度,其工作原理如图1所示。由于单波束测深过程中采取单点连续的测量方法,因此,其测深数据分布的特点是,沿航迹的数据十分密集,而在测线间没有数据。
多波束测深条带的覆盖宽度
W
W
W随换能器开角
θ
\theta
θ和水深
D
D
D 的变化而变化。若测线相互平行且海底地形平坦,则相邻条带之间的重叠率定义为
η
=
1
−
d
w
\eta=1-\frac{d}{w}
η=1−wd,其中
d
d
d为相邻两条测线的间距,
W
W
W为条带的覆盖宽度(图4)。若
η
<
0
\eta<0
η<0,则表示漏测。为保证测量的便利性和数据的完整性,相邻条带之间应有10%~20%的重叠率。
但真实海底地形起伏变化大,若采用海区平均水深设计测线间隔,虽然条带之间的平均重叠率可以满足要求,但在水深较浅处会出现漏测的情况(图5),影响测量质量;若采用海区最浅处水深设计测线间隔,虽然最浅处的重叠率可以满足要求,但在水深较深处会出现重叠过多的情况(图6),数据冗余量大,影响测量效率。
问题 1与测线方向垂直的平面和海底坡面的交线构成一条与水平面夹角为
α
\alpha
α的斜线(图7),称
α
\alpha
α为坡度。请建立多波束测深的覆盖宽度及相邻条带之间重叠率的数学模型。
若多波束换能器的开角为
12
0
。
120^{。}
120。,坡度为
1.
5
。
1.5^{。}
1.5。,海域中心点处的海水深度为70m,利用上述模型计算表1中所列位置的指标值,将结果以表1的格式放在正文中,同时保存到result.xlsx文件中。
测线距离中心点处的距离/m | -800 | -600 | -400 | -200 | 0 | 200 | 400 | 600 | 800 |
---|---|---|---|---|---|---|---|---|---|
海水深度/m | |||||||||
覆盖宽度/m | |||||||||
与前一条测线的重叠率/% | ----- |
问题 2 请考虑一个矩形待测海域(图8),测线方向与海底坡面的法向在水平面上投影的夹角为
β
\beta
β,请建立多波束测深覆盖宽度的数学模型。
若多波束换能器的开角为
12
0
。
120^{。}
120。,坡度为
1.
5
。
1.5^{。}
1.5。,海域中心点处的海水深度为120m,利用上述模型计算表 2 中所列位置多波束测深的覆盖宽度,将结果以表 2 的格式放在正文中,同时保存到result2.xlsx文件中。
问题 3考虑一个南北长2海里、东西宽4海里的矩形海域内,海域中心点处的海水深度为110m,西浅东深,坡度为
1.
5
。
1.5^{。}
1.5。,多波束换能器的开角为
12
0
。
120^{。}
120。。请设计一组测量长度最短、可完全覆盖整个待测海域的测线,且相邻条带之间的重叠率满足10%~20%的要求。
问题 4海水深度数据(附件.xlsx)是若干年前某海域(南北长 5 海里、东西宽 4 海里)单波束测量的测深数据,现希望利用这组数据为多波束测量船的测量布线提供帮助。在设计测线时,有如下要求:(1)沿测线扫描形成的条带尽可能地覆盖整个待测海域;(2)相邻条带之间的重叠率尽量控制在20%以下;(3)测线的总长度尽可能短。在设计出具体的测线后,请计算如下指标:(1)测线的总长度;(2)漏测海区占总待测海域面积的百分比;(3)在重叠区域中,重叠率超过20%部分的总长度。
注 在附件中,横、纵坐标的单位是海里,海水深度的单位是米。1海里=1852米。
附件 海水深度数据
摘要
本文针对多波束测线问题进行研究,通过遗传算法,遍历搜索法,基于Matlab软件,建立单目标规划模型,旨在得出测量海域的最佳测线方案。
针对问题一,本文通过寻找平面内的相似关系,求解不同位置的海水深度
D
D
D。根据计算所得海水深度,结合正弦定理将覆盖宽度
L
L
L分为两部分进行求解:
L
=
D
s
i
n
β
(
1
c
o
s
(
α
+
β
)
+
1
c
o
s
(
α
−
β
)
)
L=Dsin\beta(\frac{1}{cos(\alpha+\beta)}+\frac{1}{cos(\alpha-\beta)})
L=Dsinβ(cos(α+β)1+cos(α−β)1)
结合重叠率定义计算各点处的值,具体计算结果见文中表5-1。
针对问题二,本文对测线与坡面法向量的投影夹角进行分区间讨论,将问题一看作问题二的一种特解。对该海域进行局部分析,利用相似关系对相邻测线的覆盖宽度进行转化求解,建立空间直角坐标系确立覆盖宽度
w
w
w的模型为:
w
=
(
x
3
−
x
1
)
2
+
(
y
3
−
y
1
)
2
+
(
z
3
−
z
1
)
2
w
1
s
i
n
α
D
0
w=\frac{\sqrt{(x_{3}-x_{1})^{2}+(y_{3}-y_{1})^{2}+(z_{3}-z_{1})^{2}}w_{1}sin\alpha}{D_{0}}
w=D0(x3−x1)2+(y3−y1)2+(z3−z1)2w1sinα
&emsop;其中
w
1
w_{1}
w1为中心覆盖宽度,
α
\alpha
α为坡度。结合提供数据,计算结果见文中表5-2。
针对问题三,本文将该问看作对问题二航行角度的研究。当测线方向为南北方向时,以总测线长度为目标,具体约束条件见文中公式5-23,建立单目标规划模型:
m
i
n
L
=
3704
n
minL=3704n
minL=3704n
在
M
a
t
l
a
b
Matlab
Matlab中通过遗传算法生成测线条数进行迭代求解,得出最优解为68海里。依据问题二的计算结果,将测线方向分为东西测线、南北测线和倾斜测斜分类进行研究。当测线方向为东西方向时,覆盖宽度的平均重叠率为37.5%,不满足范围。当测线为倾斜方向时,不存在一夹角区间使得测线全程均在重叠率范围内,舍去。
针对问题四,本文将该问看作为问题二的实际应用。通过遍历搜索法筛选海脊的坐标数据,绘制散点图。通过建立一次拟合模型,将海脊拟合为 作为测线基准方向,其残差范数为6.772,
R
2
=
0.9165
R^{2}=0.9165
R2=0.9165,该模型具有较好的精度。将两侧的山坡通过拟合的方法调整为两个平面进行计算。以总测线长度为目标函数,以重叠率和海域面积覆盖率为限制条件,建立目标规划模型:
m
i
n
L
=
∑
i
=
1
m
i
n
n
l
0
x
i
7544.7
+
∑
i
=
1
m
i
n
m
l
0
15089.4
−
x
i
7544.7
minL=\sum_{i=1}^{minn}l_{0}\frac{x_{i}}{7544.7}+\sum_{i=1}^{minm}l_{0}\frac{15089.4-x_{i}}{7544.7}
minL=i=1∑minnl07544.7xi+i=1∑minml07544.715089.4−xi
通过
M
a
t
l
a
b
Matlab
Matlab计算,得出总测线长度为351.90海里,漏测海区占总待测海域面积的百分比为0.0068%。
一、问题背景与重述
1.1问题背景
近年来,声波的功能逐渐被挖掘,它通过电声转换和信息处理发展成了一项方便准确的测距技术,且被广泛应用于海底深度的探测实践中。但单波束测探在实际探测中却显露出成本高且效率低的弊端。多波束测深系统基于单个声波束测距技术发展起来,大大提高了海底深度测探行动的效率。但由于漏测和重复测量以及海底环境的不稳定性,优化设计测线方案成了一个重大问题。国外研究机构以多波束测深系统的更高精度、更大覆盖范围、更强水体探测等为目标优化系统结构,且各方面都在缓慢进步。为提高测量效率,建立一种多波束测探系统的方案是很有必要的。
本文针对多波束测探系统的结构进行研究分析,建立数学模型,在一定限制条件下,改变测线间隔、测线方向夹角度数等,基于前人的研究,对深海大深度声场的声传播特性研究,获取深海声场的声传播规律已有研究[1],从而找到最优多波束测量方案。这对未来优化多波束测探系统,以提高海底测量行动效率和准确率具有重要意义。优化多波束测探,使得漏测区域尽可能小,从而提高准确率;使得重叠区域面积尽可能小,从而提高探测效率,且较于单波束测探技术,大大减少了时间和成本。
1.2问题重述
本文通过建立数学模型,解决以下问题:
1.根据题意建立关于多波束测深的覆盖宽度及相邻条带之间重叠率的数学模型,并求出限定条件下所列位置的指标值;
2.根据题意建立关于多波束测深覆盖宽度的数学模型,并求出限定条件下所列位置多波束测深的覆盖宽度数据;
3.在给定矩形水域和相关条件下,根据要求设计一组测量长度最短、可完全覆盖整个待测海域的测线;
4.基于已有的海水深度数据,设计多波束测量布线,给出最优条件下测线最短总长度、最小漏测海区占总待测海域面积的百分比以及重叠率超过20%部分的最小总长度。
二、问题分析
2.1对于问题一的分析
问题一要求给出多波束测深的覆盖宽度及相邻条带之间重叠率的数学模型,并给出所列位置的指标值。本文对一般情况入手对覆盖宽度和重叠率进行分析,以相邻条带间的限制重叠率为限制条件,结合重叠率定义式和覆盖宽度计算式,建立一般情况下的目标函数。在确立目标函数后,结合问题一给出的参数,给出距离中心点不同位置的海水深度、覆盖宽度以及重叠率。
2.2对于问题二的分析
问题二要求建立多波束测深覆盖宽度的数学模型,本文对该种情况下的三维立体试图进行分析。通过对不同投影面进行研究,寻找海底坡度、测线方向与海底坡面的法向在水平面上投影的夹角以及不同测量点对应海水深度与覆盖宽度之间的关系,建立关系模型。在确立模型后将问题二所给参数代入,计算不同点处的覆盖宽度。
2.3对于问题三的分析
问题三要求设计一条测量长度最短且可以完全覆盖整个待测海域的测线。本文在确保重叠率满足变化范围的前提下,以总测线最短为目标,通过调整测量船的偏转角度确立最适合的测线。选取垂直于海域边界的测线和倾斜的测线,分别对其计算,计算在满足重叠率前提下的最短测线距离。通过比较,得出测量长度最短、可完全覆盖整个待测海域的最佳测线方案。
2.4对于问题四的分析
问题四要求依据附件中的海水深度数据给出总长度最短方案下的测线长度、漏测海区占比以及重叠率超出部分占总长的百分比。本文结合问题二计算结果,比较不同的测线方向与海底坡面法向夹角的值,结合覆盖率不超过20%、扫描条带尽可能覆盖整个海域的限制条件,建立模型对其进行分析。依据模型的分析结果,选择合适的方向设计测量方案,同时计算所要求解的数据。
三、模型假设
结合本题的实际,为确保模型求解的准确性和合理性,本文排除一些因素的干扰,提出以下几点假设:
1.假设海水可以视为均匀介质;
2.假设声波在传播过程中不会被其他物质吸收;
3.假设测量船沿测线方向做匀速直线运动,不受外界因素干扰。
四、符号说明
为了便于问题的求解,本文给出以下符号说明:
符号 | 说明 |
---|---|
Δ L i \Delta L_{i} ΔLi | 测量船第 i i i个点与第 i − 1 i-1 i−1个点间的覆盖距离之差 |
α \alpha α | 海底坡面与水平面构成的夹角 |
β \beta β | 多波束换能器的对称角 |
L 1 L_{1} L1 | 测线左侧的覆盖宽度 |
l i l_{i} li | 第 i i i条测线的测线长度 |
D D D | 海水深度 |
ε \varepsilon ε | 偏转角 |
w w w | 覆盖宽度 |
L 2 L_{2} L2 | 测线右侧的覆盖宽度 |
Δ L \Delta L ΔL | 覆盖宽度的重叠部分 |
η \eta η | 重叠部分所占前一次的比率 |
D i D_{i} Di | 测量船在第 i i i个点的海水深度 |
d i d_{i} di | 第 i i i条测线与第 i + 1 i+1 i+1条测线间的间距 |
w i w_{i} wi | 每条测线在观测面的覆盖宽度在水平面的投影长度 |
五、模型的建立与求解
经过以上的分析和准备,本文将逐步建立以下数学模型,进一步阐述模型的实际建立过程。
5.1问题一的建立与求解
本文先对只有一艘测量船的情况进行分析,结合多波束换能器的开角以及海底坡度,建立海水深度的线性回归模型,进而可以得出不同位置的覆盖宽度。在单个测量船的基础上,分别计算相邻条带之间的重叠率,得出一般状态下的模型。结合题目参数,计算问题 1 中距中心点不同距离的海水深度,覆盖宽度和重叠率。
对单个测量船的分析
对于单个测量船,在海底为倾斜坡度的情况下,其信号的覆盖宽度为在斜坡上的一段倾斜线段。为可以从一般情况推广至特殊情况,本文先以初始测量点为坐标原点,以测量船行进方向垂直的左侧为
x
x
x轴正方向,以海域中心点所在直线为
y
y
y轴,建立平面直角坐标系。
图中, β \beta β代表多波束换能器的对称角, α \alpha α代表海底坡面与水平面构成的夹角,即坡度,虚线代表坡面的一条水平面, D D D代表海水深度。 A A A点为原始测量船所在直线与海底坡面的交点,记为(0, − D -D −D),结合海底坡面角度 α \alpha α,建立线性回归模型: D = ∣ y ∣ = ∣ x t a n α − D ∣ D=|y|=|xtan\alpha-D| D=∣y∣=∣xtanα−D∣ 由模型可得测量船在不同位置时的海水深度。依据求出的海水深度,将覆盖宽度从测线处分为左右两部分,结合正弦定理计算两边的线段长度,进而得出总体的覆盖宽度。 L 1 s i n β = D c o s ( α + β ) \frac{L_{1}}{sin\beta}=\frac{D}{cos(\alpha+\beta)} sinβL1=cos(α+β)D 其中, L 1 L_{1} L1为测线左侧的覆盖宽度。对于海水宽度所对应的角度,本文过 B B B点分别做水平线和铅垂线,依据内错角相等和同位角相等,得出 γ = π / 2 − β − θ \gamma=\pi/2-\beta-\theta γ=π/2−β−θ。结合诱导公式,得出计算时的正弦定理表达式。同理可得: L 2 s i n β = D c o s ( α − β ) \frac{L_{2}}{sin\beta}=\frac{D}{cos(\alpha-\beta)} sinβL2=cos(α−β)D 其中, L 2 L_{2} L2为测线右侧的覆盖宽度,其余参数含义均与上式一致。结合上述计算过程,得出覆盖宽度 L = L 1 + L 2 L=L_{1}+L_{2} L=L1+L2,代入参数即为: L = D s i n β ( 1 c o s ( α + β ) + 1 c o s ( α − β ) ) L=Dsin\beta(\frac{1}{cos(\alpha+\beta)}+\frac{1}{cos(\alpha-\beta)}) L=Dsinβ(cos(α+β)1+cos(α−β)1) 由上述模型可得出在不同位置时的覆盖宽度,进而实现不同位置海水深度和覆盖宽度的准确计算,建立一般情况下的模型式。
5.1.2确立相邻区域重叠率的分析
在得出单个测量船海水深度和覆盖宽度的情况下,本文进一步对相邻两测量船的覆盖宽度重叠率进行分析。
上图为海底存在坡面时的覆盖距离重叠了示意图。其中,
η
\eta
η为重叠部分所占前一次的比率,
w
w
w为覆盖距离。题中给出当测线相互平行且海底地形平坦时,规定相邻条带之间的重叠率定义为:
η
=
1
−
d
w
\eta=1-\frac{d}{w}
η=1−wd 但对于该种情况,上式显然不适合作为重叠率计算公式。结合上图覆盖宽度、测线间距和重叠率之间的关系,将重叠率定义为重叠部分面积与上一次覆盖面积之比:
η
=
Δ
L
w
\eta=\frac{\Delta L}{w}
η=wΔL 其中,
Δ
L
\Delta L
ΔL代表覆盖宽度的重叠部分。结合各点处的测量数据,确立最终适合海底为坡面的覆盖率计算式,将上式完善为:
η
i
=
Δ
L
i
D
i
s
i
n
β
(
1
c
o
s
(
α
+
β
)
+
1
c
o
s
(
α
−
β
)
)
\eta_{i}=\frac{\Delta L_{i}}{D_{i}sin\beta(\frac{1}{cos(\alpha+\beta)}+\frac{1}{cos(\alpha-\beta)})}
ηi=Disinβ(cos(α+β)1+cos(α−β)1)ΔLi 其中,
D
i
D_{i}
Di代表测量船在第
i
i
i个点的海水深度,
Δ
L
i
\Delta L_{i}
ΔLi代表测量船第
i
i
i个点与第
i
−
1
i-1
i−1个点间的覆盖距离之差。在计算
Δ
L
i
\Delta L_{i}
ΔLi的数值时,本文通过确立交点的方式进行运算。假设左交点
C
(
x
1
i
,
y
1
i
)
C(x_{1i},y_{1i})
C(x1i,y1i),右交点
D
(
x
2
i
,
y
2
i
)
D(x_{2i},y_{2i})
D(x2i,y2i),假设测量点为
(
x
i
,
0
)
(x_{i},0)
(xi,0),海底坡面角
α
\alpha
α,多波束换能器对称角
β
\beta
β。
根据上述分析,分别建立线性回归方程并联立,得出
C
C
C点和
D
D
D点的坐标如下。
{
y
=
1
t
a
n
(
α
+
β
)
(
x
−
x
i
)
y
=
x
t
a
n
α
−
D
⇒
C
(
x
1
i
,
y
1
i
)
\begin{cases} y=\frac{1}{tan(\alpha+\beta)}(x-x_{i})\\ y=xtan\alpha-D\\ \end{cases} \Rightarrow C(x_{1i},y_{1i})
{y=tan(α+β)1(x−xi)y=xtanα−D⇒C(x1i,y1i)
{
y
=
1
t
a
n
β
(
x
−
x
i
)
y
=
x
t
a
n
α
−
D
⇒
C
(
x
2
i
,
y
2
i
)
\begin{cases} y=\frac{1}{tan\beta}(x-x_{i})\\ y=xtan\alpha-D\\ \end{cases} \Rightarrow C(x_{2i},y_{2i})
{y=tanβ1(x−xi)y=xtanα−D⇒C(x2i,y2i) 随后依据两点间距离公式,得出重叠宽度距离为:
Δ
L
i
=
(
x
2
i
−
x
1
i
)
2
+
(
y
2
i
−
y
1
i
)
2
\Delta L_{i}=\sqrt{(x_{2i}-x_{1i})^{2}+(y_{2i}-y_{1i})^{2}}
ΔLi=(x2i−x1i)2+(y2i−y1i)2 依据上述计算公式,可以得出一般情况下的覆盖宽度重叠率计算公式。
5.1.3对问题一参数的计算
依据上文建立的海水深度、覆盖宽度和与前一条测线的重叠率计算模型,结合问题一中的数据,即 α = 1. 5 。 , β = 6 0 。 \alpha=1.5^{。},\beta=60^{。} α=1.5。,β=60。,中心点的海水深度为 70 m 70m 70m,代入模型得出计算结果为:
测线距离中心点处的距离/m | -800 | -600 | -400 | -200 | 0 | 200 | 400 | 600 | 800 |
---|---|---|---|---|---|---|---|---|---|
海水深度/m | 90.95 | 85.71 | 80.47 | 75.24 | 70.00 | 64.76 | 59.53 | 54.29 | 49.05 |
覆盖宽度/m | 315.81 | 297.63 | 279.44 | 261.26 | 243.07 | 224.88 | 206.70 | 188.51 | 170.33 |
与前一条测线的重叠率/% | —— | 33.64 | 29.59 | 25.00 | 19.78 | 13.78 | 6.81 | 1.39 | -11.17 |
由表可得,从深水区至浅水区,覆盖宽度随海水深度的减小而减小,两相邻测线间的重叠率同时也随着海水深度的减小而减小,最终会缩减为0,即相邻覆盖宽度不重合。
上图为海底面上各测量位的端点,从图中也可明显看出随着海水深度的下降,覆盖宽度同时在减小。从整体上来看,无论是左交点还是右交点,都符合先疏后密的规律。
5.2问题二的建立与求解
5.2.1确立中心点覆盖宽度
假设海域中心点处的海水深度为
D
D
D,为便于计算研究,本文以正北方向为
y
y
y轴正方向,以正东方向为
x
x
x轴正方向,以海底中心点与海面测量船连线所在直线为
z
z
z轴,建立空间直角坐标系。
图中,
β
\beta
β代表测线方向与海底坡面的法向在水平面上投影的夹角,
α
\alpha
α代表坡度,
D
0
D_{0}
D0代表该测量点处的海水深度。本文对中心点处的覆盖宽度进行计算,当测量船位于海域中心点时,计算最外侧波束所在直线方程,得出与海底所在平面两个交点
A
(
x
1
,
y
1
,
z
1
)
B
(
x
2
,
y
2
,
z
2
)
A(x_{1},y_{1},z_{1}) B(x_{2},y_{2},z_{2})
A(x1,y1,z1)B(x2,y2,z2)。在计算两个交点时,同样采用联立方程组的方法:
{
x
=
y
t
a
n
β
=
z
t
a
n
θ
2
z
=
y
t
a
n
α
−
D
0
⇒
A
(
x
1
,
y
1
,
z
1
)
\begin{cases} x=ytan\beta=ztan\frac{\theta}{2}\\ z=ytan\alpha-D_{0}\\ \end{cases} \Rightarrow A(x_{1},y_{1},z_{1})
{x=ytanβ=ztan2θz=ytanα−D0⇒A(x1,y1,z1)
{
x
=
−
y
t
a
n
β
=
z
t
a
n
θ
2
z
=
y
t
a
n
α
−
D
0
⇒
B
(
x
2
,
y
2
,
z
2
)
\begin{cases} x=-ytan\beta=ztan\frac{\theta}{2}\\ z=ytan\alpha-D_{0}\\ \end{cases} \Rightarrow B(x_{2},y_{2},z_{2})
{x=−ytanβ=ztan2θz=ytanα−D0⇒B(x2,y2,z2) 上述计算中,仅考虑了
β
\beta
β在
(
0
,
β
/
2
)
(0,\beta/2)
(0,β/2)区间内的情况,但实际中
β
\beta
β的值应该为
[
0
,
π
]
[0,\pi]
[0,π]内的任意一个值。对此,本文对其他角度进行分析,得出不同情况下最外侧波束方程的方程表达式,以确立交点坐标。
{
x
t
a
n
β
=
y
±
1
=
z
t
a
n
θ
2
,
0
<
β
<
π
/
2
x
1
=
y
±
cot
β
=
z
t
a
n
θ
2
,
π
/
2
<
β
<
π
x
−
t
a
n
β
=
y
±
1
=
z
t
a
n
θ
2
,
π
<
β
<
3
π
/
2
x
−
1
=
y
±
cot
β
=
z
t
a
n
θ
2
,
3
π
/
2
<
β
<
2
π
\begin{cases} \frac{x}{tan\beta}=\frac{y}{\pm1}=\frac{z}{tan\frac{\theta}{2}},0<\beta<\pi/2\\ \frac{x}{1}=\frac{y}{\pm\cot\beta}=\frac{z}{tan\frac{\theta}{2}},\pi/2<\beta<\pi\\ \frac{x}{-tan\beta}=\frac{y}{\pm1}=\frac{z}{tan\frac{\theta}{2}},\pi<\beta<3\pi/2\\ \frac{x}{-1}=\frac{y}{\pm\cot\beta}=\frac{z}{tan\frac{\theta}{2}},3\pi/2<\beta<2\pi\\ \end{cases}
⎩
⎨
⎧tanβx=±1y=tan2θz,0<β<π/21x=±cotβy=tan2θz,π/2<β<π−tanβx=±1y=tan2θz,π<β<3π/2−1x=±cotβy=tan2θz,3π/2<β<2π 其中,当
β
=
0
。
,
18
0
。
\beta=0^{。},180^{。}
β=0。,180。时,代入问题一的模型进行求解,当
β
=
9
0
。
,
27
0
。
\beta=90^{。},270^{。}
β=90。,270。时,此时的覆盖宽度不随测线变化而变化,挑选中心处的覆盖宽度进行数学分析即可得出结果。对于其他非特殊情况,利用平面两点间距离公式,得出覆盖宽度
w
1
w_{1}
w1:
w
1
=
(
x
2
−
x
1
)
2
+
(
y
2
−
1
)
2
+
(
z
2
−
z
1
)
2
w_{1}=\sqrt{(x_{2}-x_{1})^{2}+(y_{2}-_{1})^{2}+(z_{2}-z_{1})^{2}}
w1=(x2−x1)2+(y2−1)2+(z2−z1)2 由此可得中心点处的覆盖宽度,为下文研究不同测量点的覆盖宽度做好基础。
5.2.2对两点间覆盖宽度关系的研究
在确立中心点覆盖宽度后,采用观测截面的方法对测线方向取相邻位置进行研究。
图中位置1表示中心点处的覆盖距离,位置2表示沿着测线方向选取的对照组,
β
\beta
β代表测线方向与海底坡面的法向在水平面上投影的夹角。将位置1和位置2的局部立体图单独划出进行分析,以得出一般状态下的模型表达式。
若不考虑海底坡面,则沿着测线方向的测量点与相应覆盖宽度可看作为一个平面三角形,同时将该三角形平移至位置2,构成一个三棱柱,其横向棱长为距离
d
d
d。绘制立体图如图5-6所示。
在此基础上,以
F
F
F点为基准点将海底坡面平面加入,与三棱柱的四个面形成交线,其中在
J
H
G
JHG
JHG平面内的交线为位置2的覆盖宽度。只需求出位置2的覆盖宽度表达式,即可得出一般情况下的多波束测深覆盖宽度模型。
5.2.3确立多波束测深覆盖宽度模型
为便于计算,本文以图5-4所建立的坐标系为基础,采用降维分析的方法求解图5-6中的各点坐标。通过研究,本文认为对于不同的 β \beta β值,在进行公式计算时略有不同。当由已知条件可得 J ( d s i n β , d c o s β , 0 ) J(dsin\beta,dcos\beta,0) J(dsinβ,dcosβ,0),通过对立体图分析,得出: w 1 w = I F J K = J G J K = D 0 K M = D 0 F K s i n α \frac{w_{1}}{w}=\frac{IF}{JK}=\frac{JG}{JK}=\frac{D_{0}}{KM}=\frac{D_{0}}{FKsin\alpha} ww1=JKIF=JKJG=KMD0=FKsinαD0 对于该相似关系在图中显示为用红色线所围成的平面 M K I N MKIN MKIN。因此,在得出中心点覆盖宽度的前提下,仅需求出 K K K点坐标即可得出覆盖宽度的模型。对于 K K K点坐标,采用线面方程联立的方法进行求解,得出直线 I F IF IF上一点为: ( 1 , t a n β , c o t θ 2 1 + t a n 2 β ) \Big(1,tan\beta,cot\frac{\theta}{2}\sqrt{1+tan^{2}\beta}\Big) (1,tanβ,cot2θ1+tan2β) 其中,若 β \beta β的范围处于 ( π / 2 , π ) ( π , 3 π / 2 ) ( 3 π / 2 , 2 π ) (\pi/2,\pi)(\pi,3\pi/2)(3\pi/2,2\pi) (π/2,π)(π,3π/2)(3π/2,2π)内时,依次将 β \beta β值变为 β − π / 2 , β − π , β − 3 π / 2 \beta-\pi/2,\beta-\pi,\beta-3\pi/2 β−π/2,β−π,β−3π/2求解,得出直线 J G JG JG上一点的坐标为: ( 1 + d s i n β , t a n β + d t a n β , c o t θ 2 1 + t a n 2 β ) \Big(1+dsin\beta,tan\beta+dtan\beta,cot\frac{\theta}{2}\sqrt{1+tan^{2}\beta}\Big) (1+dsinβ,tanβ+dtanβ,cot2θ1+tan2β) 依据 J J J点和 G G G点,得出直线 J G JG JG的表达式,联立平面: { l J G y = x t a n α − D ⇒ k ( x 3 , y 3 , z 3 ) \begin{cases} l_{JG}\\ y=xtan\alpha-D\\ \end{cases} \Rightarrow k(x_{3},y_{3},z_{3}) {lJGy=xtanα−D⇒k(x3,y3,z3) 得出覆盖宽度的模型表达式为: w = s i n α ( x 2 − x 1 ) 2 + ( y 2 − 1 ) 2 + ( z 2 − z 1 ) 2 ( x 3 − x 1 ) 2 + ( y 3 − 1 ) 2 + ( z 3 − z 1 ) 2 D 0 w=\frac{sin\alpha\sqrt{(x_{2}-x_{1})^{2}+(y_{2}-_{1})^{2}+(z_{2}-z_{1})^{2}}\sqrt{(x_{3}-x_{1})^{2}+(y_{3}-_{1})^{2}+(z_{3}-z_{1})^{2}}}{D_{0}} w=D0sinα(x2−x1)2+(y2−1)2+(z2−z1)2(x3−x1)2+(y3−1)2+(z3−z1)2
5.2.4对不同位置覆盖宽度的计算
依据上文建立的海水深度、覆盖宽度和与前一条测线的重叠率计算模型,结合问题一中的数据,即
α
=
1.
5
。
,
β
=
12
0
。
\alpha=1.5^{。},\beta=120^{。}
α=1.5。,β=120。,中心点的海深度为120m,代入模型得出计算结果。
由上表可得当
β
=
9
0
。
,
27
0
。
\beta=90^{。},270^{。}
β=90。,270。时各点处的覆盖宽度均相等,当
β
∈
(
0
,
π
/
2
)
∪
(
3
π
/
2
,
2
π
)
\beta\in{(0,\pi/2)\cup(3\pi/2,2\pi)}
β∈(0,π/2)∪(3π/2,2π)时,覆盖宽度随距离的增大而增大。当
β
=
0
。
\beta=0^{。}
β=0。时,覆盖距离的增长速度最大,当距离中心
2.1
2.1
2.1海里时,有最大覆盖宽度
770.33
770.33
770.33m。当
β
∈
(
π
/
2
,
π
)
∪
(
π
,
3
π
/
2
)
\beta\in{(\pi/2,\pi)\cup(\pi,3\pi/2)}
β∈(π/2,π)∪(π,3π/2)时,覆盖宽度随距离的增大而减小。
问题三的建立与求解
本文在求解问题三时,将测线分为垂直边和与边成夹角情况进行讨论。规定测线方向与海底坡面的法向在水平面上投影的夹角为 β \beta β。当 β ∈ ( 0 , π / 2 ) ∪ ( π / 2 , 2 π ) \beta\in{(0,\pi/2)\cup(\pi/2,2\pi)} β∈(0,π/2)∪(π/2,2π)在范围 内时,表示与边成一夹角;当 β = 0 。 , 9 0 。 , 18 0 。 \beta=0^{。},90^{。},180^{。} β=0。,90。,180。时,表示与海域边界垂直。分别对这五种情况进行讨论,结合覆盖率范围进行目标规划,得出最优解。
5.3.1对东西方向测线的研究
当
β
=
0
。
,
9
0
。
,
18
0
。
\beta=0^{。},90^{。},180^{。}
β=0。,90。,180。时,可以将
0
。
0^{。}
0。和
18
0
。
180^{。}
180。视为同一种情况,即测量船沿东西方向的测线航行。若测量船的测线为东西方向的一系列平行线时,对于每一段测线所扫过的覆盖面积均相等,仅需对相邻两段测线进行分析,结合重叠率范围即可。
结合问题二计算结果,即该种情况下随着距离浅海区越近,覆盖宽度越小的结论,本文选择在近海区相邻测线覆盖宽度重叠率为
0
0
0%作为最短测线间距。若近海处覆盖重叠率小于
0
0
0%,则无法确保完全覆盖整个待测海域;若近海处覆盖重叠率大于
0
0
0%,则会造成测量长度增长,故选择在近海处覆盖宽度重叠率为
0
0
0%。
由上图可得,此时海域中心点处的海水深度为 110 110 110米,得出深水区的海水深度为 220 220 220米。依据近海区的覆盖宽度重叠率为 0 0 0%这一前置条件,确立相邻测线距离为: d = D 0 t a n θ 2 d=D_{0}tan\frac{\theta}{2} d=D0tan2θ 其中, D 0 D_{0} D0为中心处海水深度, d d d表示相邻测线间距,计算可得 d = 110 3 d=110\sqrt{3} d=1103m。在此截面上,可以将海底看做为一个地形平坦的模型,依据题中重叠率定义式: η = 1 − d w \eta=1-\frac{d}{w} η=1−wd 计算得出此时 w = 440 3 w=440\sqrt{3} w=4403m,重叠率 η = 75 \eta=75 η=75%,但考虑到浅海区的重叠率 η = 0 \eta=0 η=0%,本文在此采用平均值的方法,计算得出平均值 η ˉ = 37.5 \bar{\eta}=37.5 ηˉ=37.5%,仍大于 20 20 20%,故舍测线为东西方向时候的方案。
5.3.2对南北方向测线的研究
当测线沿着南北方向,即 β = 9 0 。 \beta=90^{。} β=90。时,对于同一测线上的所有观测点,其到海底坡面的距离均保持不变,故只需对重叠率和测线数目进行研究即可。本文选取东西截面上的剖面图进行分析,假设测量总长度为 L L L,有 i i i条测线,每条测线在海底基准面覆盖宽度的投影长度为 w i w_{i} wi,第i条测线与第i+1条测线间的间距为 d i d_{i} di,建立目标规划模型: m i n L = 3704 n minL=3704n minL=3704n 将限制条件与上述目标规划模型进行整合,其限制条件为:
s
.
t
.
{
∑
i
=
1
n
w
i
≥
3704
i
≥
1
,
n
≥
1
x
1
=
220
3
≈
381.05
w
i
=
∣
J
i
1
−
J
i
2
∣
η
i
=
1
−
d
i
w
i
+
1
0.1
≤
η
i
≤
0.2
x
i
+
1
=
d
i
+
x
i
s.t.\begin{cases} \sum\limits_{i=1}^{n}w_{i}\ge3704\\ i\ge1,n\ge1\\ x_{1}=220\sqrt{3}\approx381.05\\ w_{i}=|J_{i1}-J_{i2}|\\ \eta_{i}=1-\frac{d_{i}}{w_{i+1}}\\ 0.1\le\eta_{i}\le0.2\\ x_{i+1}=d_{i}+x_{i}\\ \end{cases}
s.t.⎩
⎨
⎧i=1∑nwi≥3704i≥1,n≥1x1=2203≈381.05wi=∣Ji1−Ji2∣ηi=1−wi+1di0.1≤ηi≤0.2xi+1=di+xi 对于该截面,每一个测线点的坐标为
(
x
i
,
220
)
(x_{i},220)
(xi,220),坡面方程为
y
=
x
t
a
n
1.
5
。
y=xtan1.5^{。}
y=xtan1.5。。在上述约束中,
J
i
1
,
J
i
2
J_{i1},J_{i2}
Ji1,Ji2分别为第
i
i
i个测线点的左边界探测线交点和右交点横坐标。
结合上述截面坐标图,对于第 i i i个点的 J i 1 , J i 2 J_{i1},J_{i2} Ji1,Ji2计算表达式为: { y i = 3 3 ( x − x i ) + 220 y = x t a n 1. 5 。 ⇒ J i 1 \begin{cases} y_{i}=\frac{\sqrt{3}}{3}(x-x_{i})+220\\ y=xtan1.5^{。}\\ \end{cases} \Rightarrow{J_{i1}} {yi=33(x−xi)+220y=xtan1.5。⇒Ji1 { y i = − 3 3 ( x − x i ) + 220 y = x t a n 1. 5 。 ⇒ J i 2 \begin{cases} y_{i}=-\frac{\sqrt{3}}{3}(x-x_{i})+220\\ y=xtan1.5^{。}\\ \end{cases} \Rightarrow{J_{i2}} {yi=−33(x−xi)+220y=xtan1.5。⇒Ji2 在确立其计算式后,本文在 M a t l a b Matlab Matlab中采用遗传算法进行计算。通过采用自然数编码的遗传算法对模型进行求解,以种群适应度值的大小来判别个体的优势,运用选择、交叉及变异运算对种群进化改良,产生最优个体后,得到最优解,按照如下流程求解:
通过上述流程,得出n的最优解为34,即总测线长度为68海里。
5.3.3与海域边呈夹角时的研究
当测线方向与坡面法向的夹角
β
\beta
β处于 范围内时,对于同一条测线,测量船在不同的位置具有不同的重叠率。由问题二的计算结果可得当
β
\beta
β为
18
0
。
180^{。}
180。时,覆盖宽度在同一距离上的最小,即意味着需要更多的航线实现未知海域的测量。同时,当夹角不为
9
0
。
90^{。}
90。时,对于单条测线在不同位置存在不同的重叠率,如测线为东西方向时,其平均重叠率为
37.5
37.5
37.5%,仍大于重叠率限制范围。本文首先对
β
\beta
β的适合区间进行研究。
结合问题三中海底地形为西深东浅,本文限定东边浅海区的重叠率为
10
10
10%,西边深海区的重叠率为
20
20
20%进行研究,确保测线全程重叠率均满足要求,避免出现重叠过多的情况。且当满足这一条件时,即可得出
β
\beta
β的上限角度和下限角度,两者关于
9
0
。
90^{。}
90。角度对称。为简化计算过程,本文选取上限角度进行研究。
假设测线经过海域中点,要使得从深海处到浅海处的重叠率从
20
20
20%变为
10
10
10%,可以先过海域中点做一条
β
=
9
0
。
\beta=90^{。}
β=90。的线,将该线作为基准线。假设偏转角为
ε
\varepsilon
ε,过中点在水平面的投影点再做一条直线
y
=
t
a
n
(
ε
+
π
/
2
)
(
x
−
1
)
+
2
y=tan(\varepsilon+\pi/2)(x-1)+2
y=tan(ε+π/2)(x−1)+2。在此,本文将重叠率视为直线变化率,即:在基准线上距离中心点相等的两个点作垂直线,分别与偏转直线相交于两点,随后将水平面上的线段投影至坡面所形成的两段线段长度之比为
1
:
2
1:2
1:2。选取
x
x
x的值为
10
10
10和
−
10
-10
−10,则在坡面的投影点为
t
1
(
10
,
2
,
2
t
a
n
1.
5
。
)
,
t
2
(
10
,
2
−
9
c
o
t
ε
,
(
2
−
9
ε
)
t
a
n
1.
5
。
)
t_{1}\big(10,2,2tan1.5^{。}\big),t_{2}\big(10,2-9cot\varepsilon,(2-9\varepsilon)tan1.5^{。}\big)
t1(10,2,2tan1.5。),t2(10,2−9cotε,(2−9ε)tan1.5。),同理得出另外两点,依据距离公式:
2
(
9
c
o
t
ε
)
2
+
(
9
c
o
t
ε
t
a
n
1.
5
。
)
2
=
(
11
c
o
t
ε
)
2
+
(
11
c
o
t
ε
t
a
n
1.
5
。
)
2
2\sqrt{(9cot\varepsilon)^{2}+(9cot\varepsilon tan1.5^{。})^{2}}=\sqrt{(11cot\varepsilon)^{2}+(11cot\varepsilon tan1.5^{。})^{2}}
2(9cotε)2+(9cotεtan1.5。)2=(11cotε)2+(11cotεtan1.5。)2 通过计算,得出
ε
=
9
0
。
\varepsilon=90^{。}
ε=90。,与上文所规定的偏转范围相悖,故不存在这样这一个偏转角,使得测线全程的重叠率满足范围,故舍去该种情况。
5.3.4问题三小结
当测线方向与坡面法向的夹角
β
\beta
β在
β
∈
(
0
,
π
/
2
)
∪
(
π
/
2
,
2
π
)
\beta\in{(0,\pi/2)\cup(\pi/2,2\pi)}
β∈(0,π/2)∪(π/2,2π)范围内时,本文得出以下定理:对于该种情况下的夹角,不存在一个角度,满足测线重叠率在各点处均在变化范围内。
当测线方向为东西方向时,平均重叠率为
37.5
37.5
37.5%,超过变化范围,该种方案不成立。
当测线方向为南北方向时,存在总测线长度最优解为
68
68
68海里。
问题四的建立与求解
对海底地形的分析
本文先对附件中的数据进行数据预处理,将预处理后的数据导入
M
a
t
l
a
b
Matlab
Matlab。以横向坐标为
x
x
x值,以纵向坐标为
y
y
y值,以海水深度为
z
z
z值,绘制等高线地形图如下:
对测线方程的确立
为得出海脊最高点的坐标数据,本文采用遍历搜索法,按照如下步骤进行分析:
对遍历结果进行降维,即舍去每一个点的
z
z
z坐标值,得出在水平面上的坐标。将降维后的点坐标绘制散点图为:
通过分析,海沟最低点的分布特征在总体上呈现上升的趋势,且测线均为直线。故本文对以上散点图进行一次拟合,确立拟合模型为 y = k x + b y=kx+b y=kx+b:
经计算,拟合模型的表达式为:
y
=
1.364
x
−
0.5567
,
R
2
=
0.9165
y=1.364x-0.5567,R^{2}=0.9165
y=1.364x−0.5567,R2=0.9165,具有较好的拟合度。为进一步确立该拟合的模型准确性,本文通过残差范数进一步检验。通过计算,该模型的残差范数为
6.772
6.772
6.772,当残差范数在区
[
2
,
7
]
[2,7]
[2,7]内认为模型较为合理,故改模型的精度较高。
考虑到此时的单位为海里,为确保单位一致,将上式修改为
y
=
1.364
x
−
1031
y=1.364x-1031
y=1.364x−1031,故选取测量船的测线分别按照平行于该方程和垂直于该方程的方向进行测量,通过修改截距确立不同的方程。
5.4.3建立目标规划模型
为满足题目要求,本文以重叠率为约束条件,以测线总长度和海域内海底覆盖面积尽可能大为目标函数,建立目标规划模型:
m
i
n
L
=
∑
i
=
1
n
l
i
minL=\sum\limits_{i=1}^{n}l_{i}
minL=i=1∑nli
其中,
L
L
L为测线长度,
n
n
n为测线个数,
l
i
l_{i}
li为第
i
i
i条测线的测线长度。为便于代数计算,本文以拟合方程所在水平面为基准面,分别对两侧的不规则坡面进行平面拟合,计算两平面与基准面所成夹角。
上图为对两侧坡面进行拟合后的立体图,经计算,得出
α
1
=
0.64
4
。
,
α
2
=
1.540
5
。
\alpha_{1}=0.644^{。},\alpha_{2}=1.5405^{。}
α1=0.644。,α2=1.5405。。本文先对平行于拟合方程的测线方向进行研究。若测线均平行于拟合方程方向,则该种情况下的漏测海区占总待测海域面积的百分比为
0
0
0%。
对于该种情况,为使得测线总长度尽可能短,此时的相邻覆盖宽度重叠率均保持在
0
0
0%。若小于
0
0
0%则会造成海域出现漏测的情况;若大于
0
0
0%会造成测线长度边长,故选择重叠率均为
0
0
0%。此时,重叠率超过
20
20
20%部分的总长度为
0
0
0。综合上述情况,本文认为该种情况即为最佳方案,对于本题,仅需求解最短总测线长度。
对于该种情况下的总测线长度,本文选取垂直于 在水平面投影直线的平面作为观测面进行分析。同时结合该种情况下重叠率为
0
0
0,则说明由坡度变化所引起的重叠率变化在此时已被理想化,即可以不考虑。故可以简化计算观测面如下图:
在上图中,为计算两测线间的距离,本文重新建立平面直角坐标系如上图所示。依据两拟合平面端点投影点到 l l l在水平面投影线的距离,计算出 0.64 4 。 0.644^{。} 0.644。坡度在观测面的距离为 7544.7 m 7544.7m 7544.7m, 坡度在观测面的距离为 7332.5 m 7332.5m 7332.5m。接着以A点为基准探测点,向两边寻找合适的探测点并记录左边个数 n n n和右边个数 m m m,并确定其横坐标。 m i n ( n + m − 1 ) min(n+m-1) min(n+m−1)
{
m
,
m
≥
1
∑
i
=
1
n
w
i
≥
7544.7
∑
i
=
1
n
w
i
≥
7332.5
x
1
=
7544.7
w
i
=
∣
J
i
1
−
J
i
2
∣
当在左边时,
x
i
+
1
=
x
i
−
2
(
x
i
−
J
i
1
)
当在右边时,
x
i
+
1
=
x
i
+
2
(
x
i
−
J
i
2
)
d
i
=
∣
x
i
+
1
−
x
i
∣
\begin{cases} m,m\ge1\\ \sum\limits_{i=1}^{n}w_{i}\ge7544.7\\ \sum\limits_{i=1}^{n}w_{i}\ge7332.5\\ x_{1}=7544.7\\ w_{i}=|J_{i1}-J_{i2}|\\ 当在左边时,x_{i+1}=x_{i}-2(x_{i}-J_{i1})\\ 当在右边时,x_{i+1}=x_{i}+2(x_{i}-J_{i2})\\ d_{i}=|x_{i+1}-x_{i}|\\ \end{cases}
⎩
⎨
⎧m,m≥1i=1∑nwi≥7544.7i=1∑nwi≥7332.5x1=7544.7wi=∣Ji1−Ji2∣当在左边时,xi+1=xi−2(xi−Ji1)当在右边时,xi+1=xi+2(xi−Ji2)di=∣xi+1−xi∣ 其中,
w
i
w_{i}
wi为每条测线在观测面的覆盖宽度在水平面的投影长度,
d
i
d_{i}
di为第
i
i
i条测线与第
i
+
1
i+1
i+1条测线间的间距。
J
i
1
,
J
i
2
J_{i1},J_{i2}
Ji1,Ji2分别为第
i
i
i个测线点的左边界探测线交点和右交点横坐标。其计算式为(以在右边时的测线为例):
{
y
i
=
3
3
(
x
−
x
i
)
y
=
−
0.0153
x
+
50.8
⇒
J
i
1
\begin{cases} y_{i}=\frac{\sqrt{3}}{3}(x-x_{i})\\ y=-0.0153x+50.8\\ \end{cases} \Rightarrow{J_{i1}}
{yi=33(x−xi)y=−0.0153x+50.8⇒Ji1
{
y
i
=
−
3
3
(
x
−
x
i
)
y
=
−
0.0153
x
+
50.8
⇒
J
i
2
\begin{cases} y_{i}=-\frac{\sqrt{3}}{3}(x-x_{i})\\ y=-0.0153x+50.8\\ \end{cases} \Rightarrow{J_{i2}}
{yi=−33(x−xi)y=−0.0153x+50.8⇒Ji2 通过计算,得出
n
=
57
,
m
=
36
n=57,m=36
n=57,m=36,此时共有
92
92
92条测线,结合计算过程中的
x
i
x_{i}
xi,有:
l
1
=
∑
i
=
1
58
w
i
=
∑
i
=
1
58
l
0
x
i
7544.7
l_{1}=\sum\limits_{i=1}^{58}w_{i}=\sum\limits_{i=1}^{58}l_{0}\frac{x_{i}}{7544.7}
l1=i=1∑58wi=i=1∑58l07544.7xi
l
2
=
∑
i
=
1
37
w
i
=
∑
i
=
1
37
l
0
15089.4
−
x
i
7544.7
l_{2}=\sum\limits_{i=1}^{37}w_{i}=\sum\limits_{i=1}^{37}l_{0}\frac{15089.4-x_{i}}{7544.7}
l2=i=1∑37wi=i=1∑37l07544.715089.4−xi 经计算,得出总海里数为
351.90
351.90
351.90海里。下图为最终测线方向:
为确保漏测海区占总待测海域面积的百分比,本文对观测截面上的边界处进行分析,选取两侧最外围的测线投影点作为研究对象。假设观测截面上外围点分别为
g
1
g_{1}
g1和
g
2
g_{2}
g2。
对于点
g
1
g_{1}
g1,其
J
i
1
=
58.19684
J_{i1}=58.19684
Ji1=58.19684,与边界
x
=
0
x=0
x=0相比,存在
58.19684
58.19684
58.19684的距离没有测量到。对于点
g
2
g_{2}
g2,其
J
i
1
=
15042.45998
J_{i1}=15042.45998
Ji1=15042.45998,与边界
x
=
14877.2
x=14877.2
x=14877.2相比,已超出测量范围,故对于该侧已全部测量到。对于左侧的边界测线,该线距离边界的距离为
58.19684
m
58.19684 m
58.19684m,测线长度为
160
m
160 m
160m,所围成的漏测面积为
4655.7472
4655.7472
4655.7472平方米。又该海区的总面积为
20
20
20平方海里,更新漏测海区占总待测海域面积的百分比为
0.0068
0.0068
0.0068%。