水声学研究中，声场的正确预报及海底声学参数的准确反演都需要预知海底浅层沉积层的等效分层结构^[1-7]。由于沉积层等效分层结构获取的困难性，在最新的一些基于模式选择^[8]和跨维度贝叶斯算法^[9-11]的地声参数反演技术中，沉积层层数及参数也作为未知变量进行估计。这些技术由于计算量大、难以有效用于实际工程应用中，另外，跨维度贝叶斯技术得到的沉积层结构极其复杂而难为常用的声场预报模型所利用。通过分析浅层剖面仪的走航测量结果来获取实验区域浅层沉积层的等效分层结构是目前最常用的方法^[12-13]。目前浅剖数据处理大都通过人工观察的方式进行，无法将沉积层的等效结构与GPS位置信息有效结合起来，故而无法直接用于水下三维声场准确建模和预测的应用中, 也无法与其他测绘系统协同合作。本文利用混合联合信号处理方法，从浅层剖面结果中自动提取沉积层的等效垂直分层结构，为声场预报和地声参数反演提供可靠的沉积层先验信息。

1 沉积层的影响及浅层剖面仪

作为水下声传播的下边界，浅层沉积层的结构及声学参数对声传播有着显著的影响。图 1给出了3种不同沉积层结构类型下简正波传播模型^[14]预测的声场分布图。仿真计算中，声源频率为300 Hz，声源深度为10 m，海水深度为65 m，水中声速为1 500 m/s。

图 1给出了利用不同沉积层分层的海底地声模型计算得到的声场二维分布图，其中图 1(a)为半无限海底模型，其声学参数纵波声速、密度和吸收系数分别为1 550 m/s、1.8 g/cm³和0.15 dB/波长。图 1(b)为一层10 m的沉积层覆盖在图 1(a)的半无限空间的计算结果，其中沉积层的声学参数纵波声速、密度和吸收系数分别为1 475 m/s、1.6 g/cm³和0.06 dB/波长。图 1(c)的环境在图 1(b)的沉积层下面加入了一厚度为5 m的沉积层，新加入的沉积层声学参数纵波声速、密度和吸收系数为1 590 m/s、1.8 g/cm³和0.8 dB/波长，海底半无限空间的声学参数不变。

从图 1可以看出，海底沉积层结构及声学参数的不同，导致仿真的水中声场空间结构存在明显差别。在地声参数反演中，沉积层模型经常假设为一固定模型，如果模型与实际情况不符，则由匹配场反演技术得到的声学参数可能无法为声呐性能评估等应用提供可靠的环境信息。

在实际海洋环境中，表层沉积层会受到河流冲积、沉积化过程等影响，其垂直结构及空间分布十分复杂，为获得可靠的声场预报和地声参数反演结果，通过一定的信号处理方法获得沉积层等效结构的先验信息是十分必要的。

浅层剖面仪是在海水深度回声探测仪的基础上发展而来的，后来在港口建设、海底管道铺设及航道疏浚中等海岸工程中得到广泛应用^[15-16]。根据发射信号频率、声源级及沉积层底质的不同，浅层剖面仪的探测深度范围一般为十几米至几十米不等。图 2给出了海上浅层剖面仪拖曳工作示意图。一般情况下，浅层剖面仪的拖体部分随测量船移动，在设置了一定的信号频率、脉冲宽度、增益等参数后，拖体上的发射换能器会向水下垂直发射声波信号。通过累积接收到的海底反射回波的时间和强度，进而获得可用来描绘浅层沉积层的剖面图。如何从剖面中自动提取可以表征沉积层等效分层结构的信息是本文研究的重要目的。

2016年，开展了利用浅层剖面仪进行沉积层结构获取的海上综合实验。实验所用的浅层剖面仪为Edgetech 3100p，测量船按照4~6 kn的速度沿设定航线航行，浅层剖面仪拖体拖曳深度大致维持在水下5 m，发射信号脉冲形式为2~12 kHz的线性调频信号。

2 多处理方法联合提取沉积层结构

浅层剖面仪用于海上测量时会与海水深度仪和GPS记录同时进行。现有浅层剖面结果一般由专业人士通过手动调节增益进行沉积层结构信息的观察提取，无法与水深测量结果及位置信息相结合，因而很难描述沉积层的空间分布规律，也难以直接用于地声模型的建模中。本文通过联合图像统计增强、多尺度线条滤波、小波分解及重构和图像二值化的处理方式，从浅层剖面仪的输出图像中自动提取浅层沉积层的等效结构信息。

2.1 图像统计增强技术

由于环境噪声、海底界面混响、接收灵敏度等多种影响，浅层剖面仪的输出的浅剖图像中的结构会通常会比较模糊。这种模糊对于具有多分层结构海底的浅剖图像会更加明显，由于界面深度越深，其反射的能量越少，这导致接收到图像信息的能量在界面处与临近深度的对比度很低。

图 3(a)为浅层剖面仪的原始输出数据图像，图中可以看出在深度为55 m处有一较强的黑色线条，在63 m处有一较为不明显的灰色线条。这2个线条分别为海水与沉积层的交界面和沉积层与基底层的交界面。这些水平线条的识别和提取是获得沉积层垂直分层结构的重要前提，提高线条与图像其他位置上的对比度是需要进行预处理。为获得更加清晰的剖面图像，首先利用统计理论的图像均衡技术^[17]以增强界面的对比度:

$ P = \min \left| {{c_1}\left( {T\left( k \right)} \right)-{c_0}\left( k \right)} \right| $

(1)

式中：c₀是原始图像的累积统计直方图；c₁是在所有强度下图像的累积统计直方图；T代表对图像灰度的一种线性处理。

图 3(b)为经过统计均衡后的浅剖图像，相比于图 3(a)的原始图像，统计均衡后的图像在界面处有更强的能量聚积并且线条连续。这些连续的线条为后续线条提取提供了更好的图像输入，利于沉积层结构的准确提取。

2.2 多尺度线条滤波器

沉积层的分层结构大致呈水平分布的连续线条。由于海底体积混响和界面阻抗差的距离非均匀性，浅剖图像中的水平线条的宽度随深度和距离在发生着一定的变化，而且存在间断点，这为线条识别和提取的信号处理增加了一定的难度。

为应对浅剖图像中表征沉积层结构线条变化的特点，利用一种可以在一定尺度范围内检测线性结构的多尺度线性滤波器^[18]。该滤波器可以在给定尺度范围内输出对处理图像L中线性结构的最大响应：

$ M = \mathop {\max }\limits_{{s_{\min }} \le s \le {s_{\max }}} M\left( s \right) $

(2)

式中s_min和s_max是设定的最小和最大的滤波器尺度。二维多尺度线条滤波器通过计算Hessian矩阵在尺度s下的本征值来计算局部结构的响应：

$ \begin{array}{l} M\left( s \right) = \\ \left\{ \begin{array}{l} 0, \;\;\;{\lambda _2} < 0\\ \exp \left( {-\frac{{{R^2}}}{{2{\alpha ^2}}}} \right)\left( {1-\exp \left( {-\frac{{{S^2}}}{{2{\beta ^2}}}} \right)} \right), \;\;\;\;其他 \end{array} \right. \end{array} $

(3)

式中：参数α和β用来调节滤波器对R(椭圆比)和S(信噪比)的敏感度，R和S分别为：

$ R = \frac{{A/\pi }}{{{l^2}}} = \frac{{\left| {{\lambda _1}} \right|}}{{\left| {{\lambda _2}} \right|}} $

(4)

$ S = {\left\| L \right\|_F} = \sqrt {\lambda _1^2 + \lambda _2^2} $

(5)

式中：λ₁和λ₂是Hessian矩阵的本征值，而且满足λ₁>λ₂；参数A和l分别是是用来描述结构特征的最大的椭圆截面积和椭圆长轴距离。S是图像L的F范数，可以看作图像结构信息和图像背景噪声的比值。式(3)所检测的为图像中能量高点的线条结构，如果需要检测图像中能量低点组成的线条结构，只需调整λ₂的正负值。需要指出的是，尺度范围目前是根据数据的采样频率和线条厚度变化设定的，需要根据浅地剖面仪的参数设置进行选择。

图 4为图 3(b)的多尺度滤波后的输出结果，其中滤波器的尺度参数为s_min=2和s_max=16，增加量为2。相比于滤波前的图像，多尺度线条滤波器可有效检测出图像中不同宽度的条纹，并输出界面对比度更高的图像。需要指出的是，滤波后的图像中也出现了许多垂直方向上的线条结构，这部分条纹可能是多尺度滤波器同时检测各方向线条的结果，是水平线条提取的干扰。

2.3 二维小波变换分解与重构

由于垂直干扰的尺度范围与水平线条相近，为了在干扰环境下准确地提取沉积层的垂直分层结构，本文利用二维小波变换的方法^[19-20]滤除图 4中垂直方向的干扰线条并保留和增强表征沉积层分层界面的水平线条。

二维小波变换是一种成熟、有效的图像处理手段，将图像逐层分解成一系列互相独立的、表示不同细节的子图像。利用双正交小波滤波器组的二维小波变换将二维的图像分解为互相独立的方向(水平方向、垂直方向和对角方向)细节子图和平滑的低频图像：

$ {f_0} = {c_A} + {c_H} + {c_V} + {c_D} $

式中，c_A、c_H、c_V和c_D分别表示低频子图像、水平子图像、垂直子图像和对角子图像，如图 5所示。垂直子图像中包含了垂直方向线条的干扰，而水平线条的结构则在水平方向的细节子图像中被保留了下来。

图 5为对图 4的二维小波变换分解的结果，图 5(a)为低频分量子图表示原始图中的低频分量，这一分量中保存了沉积层的主要结构，保证了分层界面在重建时的连续性；图 5(b)为水平分量子图表示原始图中的具有水平分量的部分，这一分量包含了待增强的水平分层界面结构；图 5(c)和(d)分别表示原始图中的垂直和交叉方向的分量子图，这一分量表示沉积层结构中的非水平分层分量与垂直告饶，在对分层结构的重构和提取时是待滤除的干扰。在二维小波变换对多尺度滤波器输出结果进行分解后，对分解子图处理，最后输出利用逆二维小波变换重建图像的流程图如图 6所示。

对水平方向子图像所进行的阈值增强处理如下：

$ {c'_1} = \left\{ \begin{array}{l} 0, \;\;\;\;{c_1} < {T_1}\\ {P_1}, \;\;\;{c_1} \ge {T_1} \end{array} \right. $

(6)

式中：阈值T₁=c₁P₁，P₁为水平细节子图像的最大灰度值，阈值系数0 < c₁ < 1，将水平细节子图像中灰度值小于阈值T₁的点设置为0，灰度值大于阈值T₁的点设置为P₁。

对低频子图像、垂直方向子图像和交叉方向子图像进行的阈值滤波处理如下：

$ {c'_2} = \left\{ \begin{array}{l} 0, \;\;\;\;{c_2} \le {T_2}\\ {c_2}, \;\;\;{c_2} > {T_1} \end{array} \right. $

(7)

式中：阈值T₂=c₂P₂，其中P₂为所处理子图像的最大灰度值，阈值系数0 < c₂ < 1，将所处理子图中灰度值小于阈值T₂的点设置为0，灰度值大于阈值T₂的点保持不变。

在采用阈值滤波和阈值增强对4个子图进行处理后，利用逆二维小波变换对子图进行重建。图 7(a)为逆二维小波变换重建后的结果，相比于图 4，垂直条纹的干扰被完全抑制，图中只保留了表征沉积层垂直分层结构的水平线条，并且线条能量得到增强。

2.4 等效沉积层分层参数提取

为了方便提取等效分层结构参数，需要先将经过上述信号处理的输出图像(图 7(a))进行二值化处理。二值化处理结果如图 7(b)所示，用“1”表示沉积层分层界面(黑色)，其他区域用“0”表示(白色)。因为二值化的图像中水层和海底基底层均用“0”来表示，因而在任意一距离上的沉积层分层界面边界可以通过这一距离上深度方向的循环移位相减来求得。

如图 8所示，二值化后的深度向量可以通过对距离r₀处列向量(深度方向)循环移位相减，得到水平线条区域的边界列向量(参数向量)：

$ {\mathit{\boldsymbol{b}}_{{r_0}}} = \left| {{\mathit{\boldsymbol{a}}_{{r_0}}}-{{\mathit{\boldsymbol{a'}}}_{{r_0}}}} \right| $

(8)

式中：a_r0是距离r₀处深度方向向量；a′_r0是a_r0循环移位后的结果。那么分层沉积层在每一个水平距离上的等效分层数目n可以通过下式求得：

$ n = \frac{N}{2} $

(9)

式中N表示参数向量b_r0中非零值的个数，即沉积层分层界面的边界数。利用沉积层等效分层数可以对某一区域的海底进行统计建模。而分层沉积层中各分层界面的等效深度d_i可以表示为：

$ {d_i} = \frac{{{L_{{\rm{even}}}} + {L_{{\rm{odd}}}}}}{2} $

(10)

式中：L_even表示参数向量b_r0中偶数项位置，即分层界面上边界的位置；L_odd表示向量b_r0中奇数项位置，即分层界面下边界的位置。那么分层沉积层中各分层界面的等效厚度t_j表示为：

$ {t_j} = {L_{{\rm{even}}}}-{L_{{\rm{odd}}}} $

(11)

这一厚度信息通常与沉积层界面混响有关，因而可以进一步寻求分层界面厚度与沉积层介质参数之间的关系，也可以用来计算沉积层分层界面处的反射能量。

3 典型沉积层结构的提取结果

利用上述联合处理方法，对实验中获得的浅剖测线数据进行处理。图 9~11分别给出了3种典型沉积层结构的海底浅剖原始输出图像和等效分层界面提取结果：无明显沉积层分层海底、覆盖单层沉积层的海底和覆盖复杂的多层沉积层的海底。从原始输出图中可以看出，沉积层结构随距离非均匀变换，特别是深层的沉积层结构。而处理结果中，原始图像中可用来表征沉积层界面的特征均能够被准确地提取出来，尤其是在原始图像中不易直观观察到的模糊界面(如图 11所示)，这为海底地声模型的建立提供了更准确的先验信息。

在水下声场预报和地声参数反演中，图 9所示的海底可作为半无限空间海底声学模型，图 10所示的海底可近似作为单层沉积层覆盖在半无限空间的海底声学模型。图 11所示的海底沉积层结构比较复杂，沉积层分层数随距离发生一定的变化，这就需要进行分段地声模型建模才能获得准确的声传播预报和地声参数反演结果。

4 结论

1) 采用图像统计增强对浅剖输出数据进行预处理，然后利用多尺度线条滤波提取到多种尺度下的沉积层分层界面结构。

2) 结合数据特点，对子图处理的二维小波分解与重构方法，抑制了浅剖测量结果中会出现的垂直方向干扰，增强了深层界面的能量。

3) 二值化与移位相减处理自动提取到沉积层等效分层结构参数，包括沉积层等效分层数、沉积层等效深度和沉积层等效厚度。

4) 海上实测数据表明该方法可以在几种典型海底沉积层分层环境下均能准确提取沉积层的等效分层界面。

在下一步研究中，需要结合沉积物沉积过程和海洋环境因素对局部海域沉积层特性做进一步研究。