海上近岸目标威胁评估模型

2019-10-14黄钦龙夏家伟

指挥控制与仿真 2019年5期

黄钦龙,刘忠,夏家伟

(海军工程大学,湖北武汉 430000)

目标威胁评估是指根据当前战场态势信息评估敌方目标的作战能力和威胁程度,进而明确敌方意图的一种手段。作为美国JDL信息融合模型的第三层,威胁评估具有预警感知、威胁计算和威胁排序等功能。合理的目标威胁判断能够为我方指挥决策和作战资源分配提供有效的帮助,缩短决策反应时间,提高威胁目标的拦截成功率[1]。

目前,进行威胁判断的方法主要包括多属性决策理论[2-4]、模糊逻辑法[5-7]、贝叶斯网络法[8-10]以及灰关联法[11-13]。不少学者对其展开研究,并获得可喜的成果。其中,文献[14]提出了一种包含区间数的变权目标威胁评估方法,该方法对于极端条件下的威胁指标引起的系统性误差有良好的调整能力,能够给出合理的威胁评估结果;文献[15]利用指标隶属度和信息熵的方法,求解出各项威胁指标的权重,构建了空间目标威胁评估模型;文献[16]则利用拉格朗日乘子法求得权系数熵最大的指标权重,结合区间数的不同定义形式建立了基于TOPSIS法的威胁判断模型。但现有威胁评估模型多是基于离散数据对目标的能力威胁和态势威胁进行孤立分析,对目标作战意图威胁考虑较少,不能准确反映目标的潜在威胁以及威胁变化趋势。

本文提出了一种基于层次分析法和熵权法相结合的潜在威胁评估方法,该方法结合了层次分析法的主观性分析和熵权法的客观数据判断,通过分析目标的实时连续运动数据,构建了基于威胁机动模型的潜在威胁评估体系。

在港口防御监测中,近岸威胁目标具有隐蔽性强,机动性高的特点,其威胁度波动较大,易造成严重的突发事件,如美国海军驱逐舰“科尔号” 在也门亚丁港遭遇的自杀式爆炸袭击事件[17]。若能充分利用威胁目标的连续运动信息,提前预测其潜在威胁,便可有效提高对威胁目标的拦截成功率。

1 威胁评估指标及其量化处理

1.1 评估指标的选取

在港口防御中,监测区域的海上近岸目标数量多,隐蔽性强,机动性高,威胁程度难以估计,需要依靠雷达和光电等传感器获得的目标状态信息进行判断。对于常见的诸如渔船、货船、汽艇、快艇等海上近岸目标,其威胁影响因素较多,相互关系复杂,且部分信息难以直接获得。为了与实际作战相结合,现将威胁因素分为两方面:属性方面和行为方面。属性方面包括目标的威胁类型,目标武器载荷信息,目标无人化信息等;行为方面包括目标速度信息,目标距离信息,目标航向角信息等。具体如图1。

图1 海上近岸目标威胁评估指标

1.2 评估指标的量化

由1.1节可知，目标的威胁度是由多方面定性和定量信息共同决定的,部分信息难以用数值直接描述,现采用隶属度的方法消除目标信息的量纲和数量级差别,对数据进行规范化处理,将其转换为0～1的无量纲数据。

目标的威胁类型属于定性指标,对于海上近岸威胁目标,可通过敌我识别器或其他手段判别。本文依据常见的目标类型将该指标分为敌船、未知和民船,对应的威胁量化值根据专家评判如表1所示。当目标为敌船时,威胁最为明显,故量化值取最大;当无法有效识别目标类型时,存在较大潜在威胁,量化值取适中大小;对于已识别的民船,仍要有一定的警戒意识,取较小值。

表1 威胁类型量化表

武器载荷也属于定性指标,本文将其分为有明确武器,有潜在武器和未装备武器三种等级,具体信息可根据光电数据利用深度学习识别判断或人在回路参与判断。有明确武器是指可通过人工智能或检测人员直接识别出目标外部武器装备,或者根据数据库信息对比得知;有潜在武器是指目标有隐蔽武器的空间与可能;未装备武器是指已证实目标没有负载武器的空间与性能。其威胁度由高到低,量化值适当依次降低,根据专家相关经验,具体取值如表2所示。

表2 武器载荷信息量化表

无人化指标在这里是指威胁目标的作战平台是否无人,本文将其分为有人系统,未知系统和无人系统。在现代作战中,由于无人系统的非接触和零伤亡,使得其作战成本大大降低,潜在威胁更大。故无人系统,未知系统和有人系统的威胁度依次减小,其量化值取值适当依次降低,根据专家相关经验,具体取值如表3所示。

表3 无人化信息量化表

速度指标的威胁隶属度函数应满足速度越大,威胁度越大的特点,根据实际作战经验,文本采用上升型指数函数的形式,如式(1)。v0为速度威胁阈值,γv为速度威胁衰减系数,αv为最低速度威胁隶属度。根据相关经验,各项取值分别为:v0=5,γv=-0.055,αv=0.2。

(1)

威胁目标离己方保卫目标越近,则其威胁意图越大,故根据作战经验,隶属度函数可选取下降型指数函数形式,如式(2)。d0为距离威胁阈值,γd为距离威胁衰减系数,αd为最小距离威胁隶属度。根据相关经验,各项取值分别为:d0=2,γd=-0.16,αd=0.2。

(2)

威胁目标航向与保卫目标方位越接近,则其威胁度越大。为计算方便,现以威胁目标与保卫目标连线为基准,定义航向偏差角θ,顺时针为正,逆时针为负。θ0为航向偏差角威胁阈值,γθ为角度威胁衰减系数,αθ为最小偏差角威胁隶属度。根据相关经验,各项取值分别为:θ0=90,γθ=-0.03,αθ=0.2。

(3)

根据以上评估指标模型,可计算出威胁隶属度矩阵A

A=(aij)m×n

(4)

式中,m和n表示有m个目标和n个威胁指标,aij为威胁目标i的第j个指标威胁隶属度。

2 基于潜在威胁分析的修正综合评估模型

2.1 层次分析法与熵权法综合求权重

权重的确定方法包括主观赋权法和客观赋权法,主观赋权法偏向于专家的直观经验,客观赋权法更侧重于数据的严谨。现采用主观赋权中的层次分析法和客观赋权中的熵权法相结合求常权重。

层次分析法(AHP)通过将复杂问题的各个因素梳理为相互关联的有序层次,是一种定性与定量相结合的多目标决策分析方法[18]。利用层次分析法求权重具体步骤如下。

1) 构造指标比较矩阵

对于已经构建的指标体系,根据指标层次关系,将相同层级的指标进行两两比较,得到各个指标的重要性量化矩阵B,式中n为指标个数。

(5)

矩阵各个元素bij可根据指标相互重要性取值。

2) 计算相对权重

根据判断矩阵,求出其最大特征根λmax及其对应的特征向量w,方程如下

Bw=λmaxw

(6)

所求特征向量w经归一化处理,即各个指标的相对权重。

3) 一致性检验

以上求得的权重是否合理并符合专家的评判标准,还要通过一致性检验来确定。

熵本来是一个热力学概念,后由C.E.Shannons引入信息论,定义了信息熵,如果某指标信息熵越小,则其能提供的信息量越大,在评价体系中所占的权重应当更大[19]。对于m个目标,n个指标的威胁评价体系,具体权重求解步骤如下。

1) 根据威胁隶属度矩阵A求出标准化矩阵P

由上一节求出的威胁隶属度矩阵A,可求出第j个威胁指标下目标i的威胁占全部m个目标威胁的比重,即矩阵P中各个元素pij为

(7)

2) 求出各个指标的信息熵

(8)

3) 求解各个指标权重

(9)

综上,利用层次分析法求得的权重向量和熵权法求得的权重向量,可得综合权重

(10)

2.2 基于潜在威胁分析的修正综合评估

战场态势时刻不断在变化,对于目标的威胁评估也是一个动态的过程。我们可以根据目标不断变化的运动态势分析其潜在威胁,通过构建基于威胁机动模型,对目标威胁度进行修正评估,如图2所示。

图2 威胁机动模型

威胁目标在执行作战任务时,会进行必要的威胁机动,此时部分威胁指标会连续增大,如航向角不断偏向保卫目标或航速不断增大。即使此时威胁评估指标的隶属度并不大,但根据其作战意图及指标的增长率可预测其后续威胁会不断增大,即当前存在潜在威胁。为更好防范突发威胁事故,应对现有威胁评估系统进行潜在威胁修正。

假设对于传感器实时监测的威胁目标信息,即每隔ΔT监测周期进行抽取分析,可得到T1,T2,T3…Ts时刻m个目标的n项指标的威胁判断信息。如果第i个目标的j项指标威胁隶属度从T1时刻到当前Tk+1时刻连续K个监测周期持续增大,则定义该威胁目标进行了K周期的威胁机动,进行如下综合动态修正威胁评估:

(11)

(12)

式中，aij(TK)为TK时刻第i个目标的j项指标的威胁隶属度,ηj为j项指标的综合权重,rij为基于威胁运动的修正系数,fij(TK+1+k)为根据T1时刻到Tk+1时刻该指标威胁隶属度相对时间的拟合函数,K为监测周期,可根据实际需求取值。

3 威胁评估仿真分析

设某一时间段我方港口监测系统探测到5批近岸海上目标,光电与雷达获取的目标在T1～T4时刻内目标的威胁类型、武器载荷、无人化、速度、距离和航向角等信息如表4所示。

表4 T1～T4时刻目标特征信息

基于潜在威胁分析的海上近岸目标威胁评估的流程如下。

1) 计算目标隶属度矩阵

对属性表进行量化求解,得到目标的隶属度矩阵A(T1),A(T2),A(T3)和A(T4)：

2)结合层次分析法与熵权法综合求权重

对于层次分析法:首先通过专家评议,得到表5的重要性比较表格,然后利用和积法,得到其正规化矩阵B与特征向量w。

表5 指标相对重要性比较

最大特征根为λmax=6.1182,CR=0.01906<0.1满足一致性检验。

对于熵权法:获得目标权重向量为:

故综合常权重:

如果直接利用该常权重计算可得不同时刻威胁度:

各个目标T1～T4时刻的威胁隶属度如图3～6所示。

图3 T1时刻各个目标威胁隶属度

图4 T2时刻各个目标威胁隶属度

图5 T3时刻各个目标威胁隶属度

图6 T4时刻各个目标威胁隶属度

即不同时刻威胁度排序均为D3>D2>D1>D4>D5。

3)基于潜在威胁分析的修正综合评估

根据T1～T4时的威胁隶属度矩阵,现取K=3,k=1,对于符合威胁运动定义的指标,对其采用多项式拟合,求出对应fij(T5),并依据式(12)求得其修正系数rij。本案例中目标D1的速度和角度指标符合威胁运动定义,对其计算后得到r14=1.1692,r16=1.0453,代入修正后得到威胁度

综合修正威胁隶属度如图7所示。

图7 各个目标综合修正威胁隶属度

威胁度排序与修正前所有时刻都大致相似,为D3>D1>D2>D4>D5。主要区别为D1目标的威胁度有所提高,现对于D1目标进行分析:其运动状态为持续加速,并转向保卫目标,具有明显的威胁行为,在各个孤立时刻,其威胁度综合评估都并没有D2目标大,但按照其当前运动趋势,其威胁度会逐步超过目标D2,即D1目标在当前时刻有极大的潜在威胁。此威胁度排序符合实际情况。