|
1996年9月10日，联合国大会通过了《全面禁止核试验条约》（全面禁试条约）。国际监测系统建立了监测技术手段，重点监测地震波、水下声、次声和放射性核素，以确保有效执行《条约》。在各类监测台站中，根据水下深海声定位测距（sofar）信道良好的传输能力，水声站可以通过声压传感器有效地拾取水下sofar通道中的微弱声压信号，或者通过地震检波器将陡峭岛屿上的声波转换成地震信号，可以用较少的监测设备实现对广大海域的监测，因此，它已成为监测小岛屿和沿海大陆地下核爆炸和水下核爆炸等突发事件最直接、最有效的技术手段之一。同时，这种对世界海洋的连续和实时监测也可以为研究地震、火山爆发、鲸鱼叫声、南极冰山破裂和其他事件提供可靠的数据支持。IMS水声监测网络由11个水声台站组成，包括6个基于水听器阵列接收的水听器台站和5个基于地震拾取接收的T相台站。
由于潜艇隐蔽性强、自给自足性强、续航能力强、作战半径大，作为一种新型武器在现代战争中发挥着越来越重要的作用，基于潜艇的侦查和反侦察研究也越来越受到重视。潜艇探测一直是各国海军高度关注的领域，也是水声工作者的重要研究方向。随着潜艇减振降噪技术研究的深入和综合治理的逐步实施，潜艇机械噪声、螺旋桨噪声、水动力噪声等主要噪声源的辐射源水平大大降低，安静的潜艇和无声的速度等新概念不断涌现。虽然各种新的探测技术和实现方法，如航磁潜艇探测、涡流场探测技术和红外探测技术也被应用于潜艇探测，但传统的声学探测方法仍然是潜艇探测的主要方式。IMS水声站可全天值班，实现连续监测。利用水声站对潜艇进行实时、连续的监测具有重要的军事意义。 为了评估IMS水听器站对潜艇目标和爆炸事件的监测能力，本文将从声纳方程出发，分析IMS水听器站附近的噪声水平，调查潜艇辐射噪声源水平，利用基于抛物线方程的声场计算方法预测传播损失，并估计IMS位于尾流岛。HA11水听器站对潜艇和爆炸具有不同的探测能力。本文首先介绍了水听器台站的基本情况，然后从被动声纳方程、声场模拟方法和HA11附近环境噪声水平三个相关层次介绍了IMS水声台站监测能力的分析方法。根据仿真和实测数据，给出了IMS水听器站对爆炸声源和潜艇活动的监测范围。 一、IMS水声台站 IMS的11个水声台站已经全部投入使用，包含6个位于水下的H-相(水听器)台站和5个位于海岸或岛屿的T-相(地震)台站，分布在9个不同的国家，如图1所示，台站位置的详细信息在表1中给出。水听器台站主要分布在南半球，唯一在北半球的水听器台站为位于威客岛的HA11。
图1 IMS水声台站位置 表1 IMS水声台站具体位置
除HA01台站外，所有水听器台站皆位于相对较小的海岛上，由两组电缆构成，每组电缆具有3个水听器，且传感器之间彼此相距近乎2km，以获得较大的阵列孔径与监测性能。三元水听器组被置于海岛两侧以避免海岛地形阻挡效应的影响。信号采样率为250Hz，使用24bit量化。在8∼100Hz频带内，水听器灵敏度起伏不超过3dB。以HA11为例，该台站及水听器位置如图2所示。北侧水听器深度约730m，海深约1430m；南侧水听器深度为740m，海深约1180m。图2中的地形信息来源于GEBCO。
图2 威客岛台站HA11及水听器阵位置 二、IMS水声台站监测能力 分析方法 ⒈被动声呐方程 声呐方程从能量角度综合了声呐参数与声呐性能的联系，它是声呐设计和作战使用的依据，在水声工程中有十分重要的应用。被动探测的声呐方程可以表示为 DT＝SL−TL−NL＋DI ⑴ 式⑴中，SL为处理带宽内的噪声源辐射声源级，也就是水下目标的辐射噪声级，其单位为dB re 1µPa2，对声呐系统效能进行理论分析时，未知目标的SL可按照目标类型设置为一些典型值；NL为处理带宽内的海洋环境噪声级，单位为dB re 1µPa2；DI为接收指向性因子，对于单个无指向性阵元的声呐系统来说，可以认为DI＝0；DT为检测阈，是对于预定置信级下，接收机输入端所需要的接收带宽内信噪比；TL是传播损失，它与海洋环境直接有关，在对声呐系统效能进行理论分析时，TL可以使用水声传播模型计算得到。定义优质因数： FOM＝SL−NL＋DI−DT ⑵ 通过分析FOM和TL的关系，即可对声呐系统效能进行理论分析，其中最为关键的就是声呐探测距离。 ⒉基于抛物方程方法的声场计算方法 处理地形和声速剖面等环境参数随距离变化的海洋波导环境中的声传播计算方法有多种，比如适用于二维情况的水平射线方法、绝热简正波方法和N×2D抛物方程方法；解决三维传播问题的绝热简正波-射线方法以及绝热简正波-抛物方程方法。抛物方程方法由于其距离上的递进算法，适合解决水平变化情况下的声场计算问题。 本文采用基于分裂步进Pade近似抛物方程方法，计算距离水平变化的声场。涉及三维计算时，采用N×2D模式进行计算。时域声压p(r,t)满足波动方程：
第一部分表示向外传播的发散声波，第二部分表示向内传播的声波。一般假设向外传播的能量占主要地位，反向散射的声波能量较小可以忽略不计，故可以得到一阶常微分方程的形式：
可进一步进行数值计算，分别对每一项进行求解后再求和即可得到声场的声压值。 利用抛物方程方法计算以HA11台站北侧水听器阵为中心的声场传播损失分布(水层内全部深度上的平均)如图3所示，传播损失显示范围为70∼170dB。图3表明，除了部分被岛屿，如夏威夷岛，遮挡的区域，声能量基本覆盖了北太平洋的绝大部分海域。第二岛链以外的西太平洋区域，传播损失一般小于120dB；第一岛链和第二岛链之间，传播损失在120∼140dB之间；南海东北部，接近吕宋海峡处，传播损失在130∼150dB之间。计算得到的传播损失表明了可被探测到事件的最低声能量限度。
图3 以HA11(威客岛)北侧水听器阵为中心的传播损失分布(1月份，全海深平均的传播损失) IMS的6个水听器台站的监测能力同时显示在图4中，图中使用了每组台站的两组水听器。在大洋中绝大部分区域，声能量不足140dB的事件也可以被水听器台站探测到。由于岛屿等的遮挡，水听器台站对大陆附近的海域没有监测能力，对这些海域的监测可依赖于地震台站等其他监测手段。
图4 IMS水听器台站的探测能力 ⒉背景噪声分析 水听器台站处的环境噪声级对以声呐方程为基础的声呐系统探测能力分析是非常重要的参数。Lawrence分析了水声台站HA01、HA03、HA08的数据，得到了3个台站一年内的平均噪声谱级。不同台站不同水听器组附近的环境噪声差异明显，HA01台站在澳大利亚西海岸近岸，噪声相对较高，频率10Hz处的噪声谱级为75∼88dB。Harris等使用HA11台站数据分析鲸鱼的分布时，通过台站实测数据，得到HA11台站附近的噪声级。该数据比Wenz统计的噪声谱级高10∼20dB。Wenz噪声谱级总结20Hz附近的噪声谱级范围大约在50∼100dB，具体的数值与地理位置、气象条件等有关。 本文根据HA11台站几个小时的实测数据分析噪声信号功率谱，台站北侧的水听器N1和南侧的水听器S1的功率谱密度如图5所示，其他水听器的结果也类似，在10∼100Hz范围内，噪声谱级约70∼80dB。在后续的计算分析中取HA11台站水听器附近的噪声功率谱谱级为75dB。
图5 HA11台站南北两侧噪声信号的噪声谱级 作为世界级的渔场以及世界最为繁忙的航线区域之一，我国邻近海域渔船、运输船等十分众多，导致我国周边海域环境噪声级远大于CTBT水声台站所处海洋环境噪声。以2014年中国南海三亚附近非休渔期和东中国海台湾东北附近休渔期测量的海洋环境噪声为例，其噪声谱级分别如图6和图7所示。相比较于我国周边海域中未被航船污染且休渔期的环境噪声，HA11台站附近噪声谱级低10∼25dB左右，更是远低于非休渔期我国海域环境噪声。CTBT台站大多位于深海大洋处的小岛，远离大陆，航船较少，且深海环境较为稳定，具有较小的环境噪声级，为水声台站全球大洋监测提供了有效保障。考虑到我国周边海域环境噪声较高，CTBT台站对我国周边海域探测能力也较差。
图6 中国南海三亚附近海域环境噪声测量值
图7 中国台湾东北海域休渔期下的环境噪声测量值 三、监测能力分析 ⒈对潜艇的监测能力分析 对声呐系统的探测性能进行分析本质上就是比较传播损失TL与优质因子FOM的关系，以此得到声呐系统的探测距离或者探测概率。潜艇目标一般可根据声源级进行分类，如Miasniko给出了部分潜艇的源极数据，如图8所示，对常规柴油动力潜艇，在10∼100Hz频段，潜艇源级功率谱密度约为110∼135dB。随着潜艇技术的进步，新建潜艇的辐射噪声功率谱密度逐渐降低。有关文献中根据频率范围5∼200Hz内最大功率谱密度对潜艇进行分类，定义噪声(Noisy)潜艇的源级为140dB，安静型(Quite)潜艇的源级为120dB，极安静型(Veryquite)潜艇的源级为100dB。以上数值是潜艇在“极安静工况(Ultraquietoperation)”下的结果，此时，螺旋桨和电力系统引起的噪声是主要的噪声源。当潜艇提高航速时，源级会相应地增加，潜艇巡航时，速度为8kn，辐射噪声源级比“安静工况”下高约5∼10dB，在航速15kn时，会高15∼20dB。 设潜艇在10∼100Hz频段内的SL为100∼140dB，当NL＝75dB时，可计算得到无指向性水听器相应的FOM为25∼65dB。
图8 常规潜艇(柴电动力)的源级 使用基于抛物方程方法的水声传播模型计算水听器台站HA11附近的传播损失。海底地形选用如图2所示GEBCO提供的数据，声速剖面选用全球海洋地图集(WOA)提供的月平均数据。根据互易原理，以水听器深度为声源深度，在每个方位角上计算距离-深度的二维传播损失。声场分布的部分结果如图9所示，图9(a)∼图9(c)分别对应HA11台站北侧水听器在315°方位角的传播损失结果，图9(e)∼图9(g)分别为225°方位角的传播损失，频率依次是10Hz、50Hz和100Hz。从图9中可以看出，随着频率的增加，声场的空间分布起伏越发剧烈且海底地形对声场有一定的遮蔽作用。
图9 方位角315°和225°以HA11北侧水听器位置作为声源位置的传播损失二维分布(使用1月水文数据，频率依次为10Hz、50Hz、100Hz) 图10为HA11北侧水听器附近的声场在不同深度上的二维分布(10∼100Hz宽带平均的结果)。深度变化导致的声场变化明显，从而导致潜艇在不同深度上探测距离有明显变化。图11为不同月份的结果对比，由于水听器处于深海声道轴附近，深海的声速剖面比较稳定，表层的声速变化对声场的影响较小，所以不同月份引起的探测能力变化很小，可忽略不计。
图10 以HA11北侧水听器位置为中心的声场分布(10∼100Hz宽带结果，依次为不同深度上的水平切面)
图11 不同月份以HA11北侧水听器位置为中心的声场分布(100m深度，10∼100Hz宽带结果) 从图10可以看出，对于FOM＝60dB的情况，台站的探测距离也只局限在水听器附近非常小的距离范围内。深海近距离的声场传播损失可以近似用球面衰减(20lgr)来计算。当优质因子FOM＝60dB时，球面衰减的距离仅仅为1km，图10所示的距离范围(半径100km)对于潜艇探测能力估计来说过大，不利于进行准确的图像分析。将图10放大到半径10km的范围，如图12(a)、图12(b)分别为潜艇位于100m和300m深度时的传播损失分布。从图12中可以看出，传播损失存在不同方向上的差异，在某一方向上，传播损失也不一定是单调递减的。
图12 以HA11北侧水听器位置为中心的声场分布(10∼100Hz宽带结果) 表2和表3给出了HA11台站北侧和南侧水听器对不同源级的潜艇在不同下潜深度(100m、300m和500m)的探测距离估计，计算时使用10∼100Hz的宽带传播损失。表中的探测距离是所有方向上的平均结果，每个方向上的探测距离为大于指定FOM值的传播损失对应的最大距离。从表中可以看出，当FOM≤50dB时，对任意深度的潜艇目标均没有探测能力，此时潜艇的源级为125dB(10∼100Hz频带内的谱级)，已经高于安静型潜艇在低速时的数值。利用潜艇噪声限值模型，可计算得到当低频连续谱谱级为125dB时，其线谱的最大值约为131dB，此时对于500m深度的目标探测距离仅400m。当FOM进一步提高时，才有了一定的探测能力，例如FOM＝60dB时，对300m深度的目标探测距离为0.9km。 表2 不同潜艇源级下的探测距离估计(HA11,North)
表3 不同潜艇源级下的探测距离估计(HA11,South)
图13给出了0°方位角上，HA11台站北侧和南侧的声传播二维分布结果，图中的等值线形状随着距离的增加会发生改变。从图中可以看出不同深度上的探测能力不同，在声道轴以上探测距离随着潜艇深度的增加而增加，这是由于水听器位于深海声道轴上，能量主要分布在声道轴附近，越接近声道轴，能量越强，可探测距离就越大。不同方位角上的探测能力不尽一致(如图12所示)，表2和表3为全部方位角的平均，这种平均结果会导致在某些FOM值上的探测距离随目标深度的规律发生改变，如表3中FOM＝65dB时的结果。
图13 以HA11南北站位水听器位置为中心方位角0°处的声场分布(10∼100Hz宽带平均) ⒉对水下爆炸事件的监测能力分析 水下突发事件经常伴随水下爆炸，水下爆炸产生的声信号一般包含冲击波及多个气泡脉动成分，如图14所示。Cole在文献中对水下爆炸现象进行了分析归纳，系统总结和研究了水下冲击波形成及传播理论、冲击波测量、界面效应等，形成了一套较为完整的水下爆炸理论，建立了一定范围内爆炸流场中冲击波压力峰值、比冲量及能量密度的计算公式，此理论被大量试验所验证，并有很高的精度。
图14 水下爆炸产生的冲击波及气泡脉动 图15为1kg三硝基甲苯(TNT)当量的爆炸声源的源级谱级均值(考虑了冲击波和气泡脉动)，分别在100m和300m深度爆炸的结果。在100Hz频点上，谱级约为203dB。爆炸产生的冲击波和气泡脉动之间的时延导致频谱上出现一定的干涉结构，频率小于100Hz时，不同爆炸深度的声弹声源级谱级有明显差异。在估计水听器台站对水下爆炸的监测能力时，认为1kgTNT当量的声源谱级为205dB，考虑冲击波和气泡脉动的总能量，爆炸当量每扩大10倍，源级增加6.4dB，则得到不同爆炸当量的源级谱级，如表4所示。
图15 1kgTNT当量爆炸声源级均值 表4 爆炸源级与当量的关系 表4中的数据是以爆炸当量为kg量级的实测数据外推得到，对于大当量的水下爆炸来说，具有参考意义。但对于实际情况，水下爆炸产生的水下声源级与多种因素有关，如爆炸深度、爆炸类型。水下突发事件伴随的爆炸，一般相当于100kgTNT当量量级，如2017年11月15日失事的阿根廷潜艇圣胡安号相关的爆炸。对于水下核爆的情况，其和化学爆炸也有区别。与化学爆炸相比，核爆释放的能量更多转换为冲击波。化学爆炸产生的气泡中主要是爆炸产物，可以认为是同质(Homogeneous)的；而核爆产生的气泡包含蒸汽，是不同质(Non-homogeneous)的。在气泡振荡过程中，气泡-水界面上的泰勒不稳定性(Taylorinstability)使水进入气泡并冷却气泡，导致更多的能量损失，使得气泡的振荡次数减少。化学爆炸中，首次气泡的40%能量保留在二次气泡中，而核爆中首次气泡的能量只有8%保留在二次气泡中。这导致核爆的气泡脉动次数比化学爆炸的明显要少，如图16所示。
图16 气泡能量损失及脉动次数 水下爆炸产生的冲击波和气泡脉动是远距离传播的主要形式，不同爆炸类型及爆炸深度都会影响实际的源级。当核爆发生的深度较浅时，产生的气泡不完整，大量能量传递到大气中，也会降低其产生的水下声场强度，从而影响水声台站对此类事件的监测能力。 根据表4中的数据，对比前面计算得到的传播损失图像，对于1kgTNT当量的爆炸，HA11台站的监测能力覆盖了北太平洋的大部分区域。对于更大当量的水下爆炸，监测能力可部分覆盖中国南海东北部及中国东海东部部分海区。
8 T/ q' q8 a" J2 u四、结论 + V- F6 ^- G9 E" G K" V1 e3 G
9 p8 E; m5 p& i2 g( W( r! b1 B
IMS水声台站利用深海信道实现了大范围的核爆监测。由于潜艇的辐射噪声源级很低，现有IMS水听器台站对潜艇的探测能力非常弱。只有当潜艇本身源级较高或者因航速提高而源级增大时，才可能有1km左右的探测能力。为了提高对潜艇的被动水声探测能力，有必要利用接收系统的阵处理增益或者时间增益，以提高声呐系统的优质因子FOM。另外，当潜艇必须通过位于深海声道轴附近的接收水听器时，根据本文的仿真结果，应该选择较浅的工作深度，并选择噪声级较小的工况(航速)。而对于水下爆炸引起的突发事件等，由于水下爆炸的声源级较高，水中的传播损失低，IMS水听器台站对水下爆炸事件的监测能力可覆盖绝大部分大洋区域。 6 F8 p$ P& m8 b
| 
J R* i) l# W6 P: [/ c! u$ C