随着人类文明的发展，海洋的战略地位日益突显，各国积极开发利用海洋资源和空间，采用各种手段对海洋的水下环境进行探测。 声波是目前在海洋中唯一能够进行远距离传播的能量形式，因此，在众多的海洋水下探测设备中，声波探测是其主要使用的技术手段。
目前，水下声探测设备种类多，数量大，已广泛用于海洋的水下探测和调查研究[1-8],例如海洋的工程地质勘探、海底地形地貌测量等。 但水下声探测也是外国调查船和水下文物盗捞船进行非法调查和水下文物探测、定位的主要技术手段，这些海洋水下声探测活动可以通过探测信号的接收声特征进行区分，因此，开展海洋水下声探测信号的采集研究对非法调查的监视取证和水下文物保护具有重要意义。
在海洋声学中，水声信号处理的研究驱动主要来自军事需求，重点关注水下目标的辐射噪声特性、水下目标的声回波特性、水下声场的信息获取与处理等[9-14],另外，一些海洋生物的发声特性也因其军事和海洋生态效应逐渐被关注[15-18],对于水下声探测设备发射声信号的分类、被动采集、分析和识别方面的研究则极少，一直以来，人们主要关注水下声探测设备在海洋调查和海洋工程勘探中的使用[1-8],即强调发射声信号对海洋环境的主动探测过程。 目前，关于海洋水下声探测信号的分析研究方面，有基于希尔伯特黄变换的 C-BOOM 浅剖信号分析，通过集合经验模态分解，实现对未知频带的 C-BOOM 浅剖信号的滤波[19]; 还有对多波束测深仪与侧扫声纳信号回波检测技术的分析，从理论上给出了两种仪器的海底回波信号检测方法[20]. 本文则立足实际海洋工程作业过程中对海洋水下声探测信号的被动测量，在未知作业设备目标信息的条件下，实现对海上测量数据的分析、提取，结合对已知的水下声探测设备信号声特征的分类，来实现对海洋水下声探测活动的初步识别。
水下声探测设备的生产厂家或研制人员均会给出设备一些主要的技术参数，包括声学特性参数。 但海洋是一个随机时变、空变的复杂声信道[21],水下声探测设备的发射声信号经海洋这个水声信道调制后会产生畸变和信息损失，并受到海洋中其他声信号的干扰，导致接收的水下声探测信号无法直接区分、识别。 因此，本文依据设备给出的声信号的频带范围对主要的海洋水下声探测信号进行分类，给出各类水下声探测信号的主要声特征； 采集了 4 种代表性的海洋水下声探测信号，通过对 4 种实测的海洋水下声探测信号数据的分析、研究，给出了对应类型的海洋水下声探测信号的常规分析方法； 根据已知的水下声探测信号的声特征实现对分析目标信号的识别。
1 海洋水下声探测信号分类
外国调查船、水下文物盗捞船所使用的主要水下声探测装备中，会根据应用类型不同选择不同的声源和声信号，而根据声源及声信号所使用的频段大致可分为两大类---宽带中低频声源和窄带（ 单频或双频为主） 高频声源。
浅地层剖面测量、海底地震测量、水声场测量等主要使用宽带中低频声源进行探测，其使用的中低频声源有电火花声源、枪震源（ 包括气枪、水枪、蒸汽枪和枪阵等） 、炸药震源、剖面仪和 BOOMER 震源等[1-6]. 其中浅地层剖面测量使用的频率相对较高，且信号类型确知，目前各型浅地层剖面仪的声源频率主要在几百 Hz 到十几 kHz; 枪震源和炸药震源的主能量频率最低，一般在几十 Hz 到几百 Hz,为宽带脉冲信号； 电火花声源和 BOOMER 震源的主能量区则一般在几十 Hz 到几 kHz,也是宽带脉冲信号。
多波束海底地形测量、侧扫声纳测量及测深侧扫声纳测量等地形地貌测量设备主要使用窄带（ 以单频或双频为主） 高频声源，信号类型确知，为高频短脉冲信号，其换能器探头工作频段主要在几十kHz 到几百 kHz[6-8].
2 海洋水下声探测信号分析
水下声探测设备的声源特性不同，其声探测信号的分析方法也不同。 对于宽带中低频声探测信号，易受船舶辐射噪声、风雨噪声等海洋环境背景噪声的干扰，需从时域、频域、时频域或空域等采用多种方法进行分析和识别。 对于窄带（ 单频或双频为主）高频声探测信号，其声信号频率高，信号形式稳定，抗干扰能力较强，在功率谱分析的基础上，采用传统的频域滤波方法就可分析、识别。 本研究将分别以两类声探测信号中的典型信号为例，探讨两类信号的分析、识别和提取方法。 下面首先对声探测信号数据的来源与海上实验测量情况进行简要介绍。
2. 1 海上实验简介
利用地质部门进行海上地质工程勘探的机会，项目组人员对地质工程勘探过程中发射的海洋水下声探测信号进行了采集。 采集的声探测信号类型有浅地层剖面仪信号、电火花声源信号、双频侧扫声纳信号和单频测深信号。 实验地点在福建省莆田市秀屿区外的兴化湾海域，测量时水深约 6 m,天气晴朗，海况 2级左右。 测量的过程中收发同船，测量仪器的连接及船上的收发位置如图1、2 所示：
图 1 中，水下声探测信号发声后经海底海面反射，由 RESON-TC4014 水听器接收，转化为电信号送入滤波调理放大器（ SR640） ,最后送入数据采集存储系统（ NI6122 高频水声信号采集卡） . 考虑到SR640 带宽不够，在采集窄带高频信号时将其设置为直通。
图 2 中，船舶发动机位于各水下声探测信号和接收水听器之间。 电火花声源距接收水听器约 16m,船舶发动机距接收水听器约 5 m. 测深仪距接收水听器约 5 m,侧扫声纳距接收水听器约 8 m,浅地层剖面仪距接收水听器约 8 m. 整个实验测量过程中均有船舶发动机噪声干扰。
2. 2 宽带中低频声探测信号分析
出海采集了宽带中低频声探测信号中的浅地层剖面仪信号和电火花声源信号。 浅地层剖面仪以强指向性指向海底，因此采集到的浅剖信号波形主要是海底海面反射信号； 电火花声源无指向性，但实验海域水深仅6 m,收发距离16 m,因此采集到的电火花声源信号波形应是直达波和海底海面反射信号的叠加。
虽然有干扰，但从图 3 中可以看到似乎存在一个周期性发射的信号。 从图 4 中可以看出浅地层剖面仪信号主要在几 kHz 到十几 kHz 范围，将图 4 横坐标放大后可以确定浅地层剖面仪信号的主能量区约在 2 ~12 kHz. 为此将含噪原始信号进行 1. 2 ~14. 0 kHz 的带通滤波，带通滤波后的时域波形如图5、6 所示。
从图 5 中可以看出，经带通滤波后干扰噪声大大减少。 图 6 为浅地层剖面仪信号经带通滤波后的单个时域波形样本，从图中可以看出明显的信号疏密变化，这是典型的线性扫频调制特征。 因此，分析结果无论是从频域（ 图 4） 还是从时域（ 图 6） 上来判断，均与已知的浅地层剖面仪信号特征是一致的。电火花声源原始测量信号的时域波形和时频分析分别如图 7、8 所示。
从图 7 的时域图中可以看到干扰，但也存在一个周期性发射的信号。 从图 8 的时频分析中可以知道电火花声源信号是宽频带信号，但采集的含噪信号中，在 5 kHz 以下和 16 kHz 附近存在强干扰。 由文献[22]可知，电火花声源信号的主能量区一般在几百 Hz 到几 kHz,因此，对于 5 kHz 以下的强干扰无法通过频域滤波的方法滤除。 由图 7 可知，时域的幅度滤波方法也不可行。 针对这类特殊的信号提取问题，需从空域角度分析，采用小波尺度相关滤波及其改进算法[23]进行处理，有效的提取了电火花声源信号。 尺度相关滤波后的信号时域波形及时、频分析如图 9、10 所示。从图中可以看出，经小波尺度相关滤波后，5kHz 以下和 16 kHz 附近的强干扰得到了有效抑制，仅在信号间隔期内存在少量干扰。
2. 3 窄带高频声探测信号分析
出海采集的窄带高频声探测信号中包含单频的测深信号和双频的侧扫声纳信号。 测深仪和侧扫声纳均具有强指向性，因此，接收信号不是直达波，而是海底或者海底海面的多次反射信号，相比初始样本信号会有一定的畸变。 接收信号的时域波形和功率谱分析分别如图 11、12 所示。
由于低频的背景干扰幅度较大，无法从图 11中辨识信号的存在。 从图 12 的功率谱分析中可以看出，在 135 kHz 和 205 kHz 附近存在两个强频谱峰值。 为此，对图 11 的时域信号分两次进行带通滤波，通频带分别为： 110 ~150 kHz 和 195 ~215 kHz,得到的时域波形如图 13、14 所示。
从时域上看，图 13、14 均存在周期性的发射信号。 已知所用的侧扫声纳为 KLEIN SYSTEM 3000,通过查阅该设备的出厂参数可知其信号频率为 132、445 kHz. 本次实验数据采集设备的最大采样频率为 500 kHz,且滤波调理放大器 SR640 的带宽不够，在采集窄带高频信号时设为直通，所以对于 445kHz 的信号会出现欠采样，在此不予讨论。 因此，图13 带通滤波后得到的应是 132 kHz 的侧扫声纳信号。 已知所用的测深仪为 HY1600,查阅该设备的出厂参数可知其信号频率为208 kHz,因此，图14 带通滤波后得到的应是 208 kHz 的测深仪信号。
3 结语
本研究根据水下声探测信号的主要声特征，对主要海洋水下声探测信号进行了分类，并分为宽带中、低频信号和窄带高频信号，其中宽带中、低频信号又可分为信号类型确定的中、低频信号和宽带脉冲中、低频信号。 在福建兴化湾海域采集了浅地层剖面仪信号、电火花声源信号、侧扫声纳信号和单频测深信号 4 种典型的海洋水下声探测活动信号样本，并进行了分析、提取，得到如下结论：
（ 1） 浅地层剖面仪信号是信号类型确定的中、低频信号的典型代表，其信号波形为线性扫频调制脉冲。 针对这类信号的分析、处理和识别，需先进行初步的时、频分析，掌握初步的频域特性，然后针对性的进行带通滤波，提取目标信号，最后结合已知的海洋水下声探测信号的声特征，进行信号识别。
（ 2） 电火花声源信号是宽带脉冲中、低频信号的典型代表，该类信号没有确知的信号形式，以宽带的冲击波为主，针对这类信号的分析、处理和识别，可以从空域角度分析，采用小波尺度相关滤波的方法进行处理，提取和识别目标信号。
（ 3） 侧扫声纳信号、单频测深仪信号是窄带高频信号的典型代表。 这类信号频率很高，以单频或双频为主，信号类型（ 波形） 固定。 因此，在对信号进行初步的时、频分析后，再对信号进行带通滤波，最后结合已知的海洋水下声探测信号的声特征，提取和识别目标信号。
本文数据来源于实际海洋工程勘探作业过程中的被动观测，通过数据分析给出了代表性的海洋水下声探测信号的分析、提取和识别方法，可为海上辨识外国调查船和水下文物盗捞船的非法水下声探测活动提供参考和帮助。 随着技术的发展，对于本文可能没有涉及到的海洋水下声探测信号，还需根据实测数据情况进行分析。
参考文献：
[1] 裴彦良，王揆洋，李官保，等。 海洋工程地震勘探震源及其应用研究[J]. 石油仪器，2007,21（ 2） : 20-23.
[2] 万芃，吴衡，王劲松，等。 海洋高分辨反射地震勘探震源的技术特征[J]. 地质装备，2010,11（ 3） : 21-28.
[3] 左公宁。 浅地层反射法勘探震源[J]. 电工电能新技术，1993（ 4） : 17-21.
[4] 金翔龙。 海洋地球物理研究与海底探测声学技术的发展[J]. 地球物理学进展，2007,22（ 4） : 1 243-1 249.
[5] 李一保，张玉芬，刘玉兰，等。 浅地层剖面仪在海洋工程中的应用[J]. 工程地球物理学报，2007,4（ 1） : 4-8.
[6] 吴自银，郑玉龙，初凤友，等。 海底浅表层信息声探测技术研究现状及发展[J]. 地球科学进展，2005,20（ 11） : 1 210-1 217.
[7] 刘忠臣，周兴华，陈义兰，等。 浅水多波束系统及其最新技术发展[J]. 海洋测绘，2005,25（ 6） : 67-70.
[8] 王朝英。 测深侧扫声呐仿真研究[D]. 哈尔滨： 哈尔滨工程大学，2007: 1-9.
[9] 王之程，陈宗岐，于沨，等。 舰船噪声测量与分析[M]. 北京： 国防工业出版社，2004: 201-212.
[10] 李启虎。 中国水声信号处理 30 年[J]. 中国科学： 物理学 力学 天文学 （ 海洋声学专刊） ,2013,43（ s1） : 32-41.
[11] 丁玉薇。 被动声纳目标识别技术的现状与发展[J]. 声学技术，2004,23（ 4） : 253-257.
