警戒海域反潛直升機(jī)優(yōu)化調(diào)度研究

曾 斌, 姚 路, 李厚樸

(1. 海軍工程大學(xué)管理工程與裝備經(jīng)濟(jì)系, 湖北 武漢 430033;2. 海軍工程大學(xué)導(dǎo)航工程系, 湖北 武漢 430033)

0 引 言

反潛戰(zhàn)是海軍作戰(zhàn)的主要樣式之一,是同對(duì)方潛艇兵力作斗爭(zhēng)的海上作戰(zhàn)。由于海洋環(huán)境復(fù)雜多變,潛艇具有非常強(qiáng)的隱蔽性,給反潛作戰(zhàn)帶來(lái)很大難度。隨著無(wú)人水下航行器的發(fā)展,反潛的地位和作用更為重要,其中艦載反潛直升機(jī)由于其具有飛行速度快、搜潛范圍廣、隱蔽性好以及攻擊效率高等優(yōu)點(diǎn),是海上艦艇編隊(duì)對(duì)潛作戰(zhàn)的主要兵力[1]。

從作戰(zhàn)目標(biāo)上看反潛分為護(hù)航反潛和區(qū)域反潛,本文主要針對(duì)區(qū)域反潛,區(qū)域反潛又細(xì)分為:① 反潛封鎖區(qū),使用陣地防潛器材和機(jī)動(dòng)反潛兵力以阻止對(duì)方潛艇通過(guò)的海區(qū);② 海洋反潛區(qū),將海洋劃分成若干個(gè)反潛區(qū),各區(qū)編組一定的機(jī)動(dòng)反潛兵力,分區(qū)負(fù)責(zé)反潛搜索和跟蹤;③ 反潛警戒線,在沿海和島嶼附近的海岸聲納組網(wǎng),以便提供反潛預(yù)警,并輔助反潛兵力進(jìn)行反潛搜索。

近年來(lái),反潛模式也由平臺(tái)中心戰(zhàn)向網(wǎng)絡(luò)中心戰(zhàn)轉(zhuǎn)型[2],以此縮短反潛觀察-調(diào)整-決策-行動(dòng)閉環(huán)的時(shí)間。通過(guò)各類傳感器等預(yù)警偵察手段,獲得準(zhǔn)確的戰(zhàn)場(chǎng)態(tài)勢(shì);利用反潛機(jī)的快速機(jī)動(dòng)能力,實(shí)施召喚式快速搜索[3]。

反潛警戒海域柵格化為若干監(jiān)測(cè)區(qū),在水下聲納網(wǎng)絡(luò)[4-5]、巡邏反潛、偵察衛(wèi)星、無(wú)線電偵聽(tīng)配合下,為重點(diǎn)海域(海上武器試驗(yàn)區(qū)、重要航道、海上保障基地、艦艇編隊(duì)等)提供反潛預(yù)警;當(dāng)發(fā)現(xiàn)水下可疑目標(biāo)后向反潛指揮部門(mén)發(fā)出警報(bào),警報(bào)中包含該目標(biāo)隨處區(qū)域以及反潛需求;指揮部門(mén)接到警報(bào)后,再調(diào)度我方反潛直升機(jī)到達(dá)該監(jiān)測(cè)區(qū)實(shí)施搜潛、跟蹤、攻擊或驅(qū)逐,反潛機(jī)完成任務(wù)后返回待召。

隨著無(wú)人水下航行器(unmanned underwater vehicle,UUV)的發(fā)展,反潛作戰(zhàn)面臨的挑戰(zhàn)更加艱巨。藍(lán)方往往頻繁地從多處對(duì)紅方警戒海域進(jìn)行水下侵入,再加上水下環(huán)境復(fù)雜[6],水下預(yù)警系統(tǒng)也存在誤報(bào)的情況,而反潛直升機(jī)資源數(shù)量有限,傳統(tǒng)的就近派遣策略無(wú)法滿足這種多次數(shù)多分散地點(diǎn)的反潛需求。因此，優(yōu)化的調(diào)度策略更顯重要。

直升機(jī)反潛是當(dāng)前軍事領(lǐng)域研究熱點(diǎn)問(wèn)題之一。例如,文獻(xiàn)[7]利用貝葉斯網(wǎng)絡(luò)推理潛艇意圖。文獻(xiàn)[8]提出了適合于雙機(jī)搜潛的“前堵后追”法。文獻(xiàn)[9]提出了一種直升機(jī)協(xié)同搜潛的陣型最優(yōu)化方法。文獻(xiàn)[10]研究了聯(lián)合反潛作戰(zhàn)中反潛機(jī)的搜索模式問(wèn)題。文獻(xiàn)[11]將搜索海域柵格化,給出了艦載直升機(jī)吊放聲納區(qū)域反潛搜索算法。文獻(xiàn)[12]研究了反潛機(jī)應(yīng)召搜索時(shí)聲納浮標(biāo)陣的評(píng)估模型。文獻(xiàn)[13]從協(xié)同反潛角度研究了對(duì)潛艇的探測(cè)模型。以上文獻(xiàn)主要研究反潛直升機(jī)起飛后的搜潛策略,但反映出反潛作戰(zhàn)的背景特征,如潛艇位置的隨機(jī)性和搜索海域柵格化等思路,對(duì)本文具有較大啟發(fā)。

文獻(xiàn)[14]以提高覆蓋率為目標(biāo),建立了多艦聯(lián)合搜潛的概率方程,并提出基于概率距離函數(shù)的優(yōu)化算法來(lái)求解。文獻(xiàn)[15]研究了瀕海戰(zhàn)斗艦的反潛規(guī)劃問(wèn)題。文獻(xiàn)[16]研究了應(yīng)召反潛中,水面艦艇和反潛機(jī)協(xié)同反潛中拖曳聲納和反潛直升機(jī)的部署問(wèn)題。文獻(xiàn)[17]利用無(wú)人艇執(zhí)行搜潛任務(wù),提出了改進(jìn)的遺傳算法來(lái)提高累積探測(cè)概率。文獻(xiàn)[18]研究了反潛直升機(jī)的優(yōu)化攻擊地點(diǎn)。文獻(xiàn)[19]設(shè)計(jì)了反潛機(jī)的探測(cè)距離估算模型。以上文獻(xiàn)在反潛監(jiān)測(cè)概率和優(yōu)化部署方面對(duì)本文具有一定啟發(fā)。

另外,國(guó)外軍事研究人員,特別是美國(guó)海軍也針對(duì)反潛直升機(jī)的作戰(zhàn)使用展開(kāi)了相關(guān)研究。文獻(xiàn)[20]設(shè)計(jì)了一個(gè)上下文驅(qū)動(dòng)的決策支持工具來(lái)解決反潛調(diào)度和規(guī)劃問(wèn)題,提出了一個(gè)隱馬爾可夫模型來(lái)對(duì)反潛資源和搜索規(guī)劃建模,并提出了二階段算法(進(jìn)化算法和拍賣算法)來(lái)求解。文獻(xiàn)[21-22]利用零和博弈對(duì)反潛作戰(zhàn)建模,利用線性規(guī)劃算法計(jì)算反潛資源(反潛機(jī)或艦艇)的優(yōu)化部署問(wèn)題。文獻(xiàn)[23]利用非零和網(wǎng)絡(luò)阻斷博弈對(duì)區(qū)域反潛建模,并采用強(qiáng)Stackelberg均衡策略解決了反潛資源的分配問(wèn)題。文獻(xiàn)[24-25]以提高覆蓋率為目標(biāo),研究了反潛規(guī)劃問(wèn)題,并運(yùn)用解析法建立了反潛責(zé)任區(qū)防潛兵力的配置模型。文獻(xiàn)[26]研究了帶吊放聲納的反潛直升機(jī)搜索模式和吊放頻率。文獻(xiàn)[27]進(jìn)一步研究了反潛直升機(jī)搜索策略以及油料和載荷的優(yōu)化配置問(wèn)題。但現(xiàn)有文獻(xiàn)缺少?gòu)乃戮W(wǎng)絡(luò)中心戰(zhàn)角度研究多監(jiān)測(cè)區(qū)-多反潛機(jī)的調(diào)度問(wèn)題。

本文的研究目標(biāo)是:把警戒海域劃分為多個(gè)監(jiān)測(cè)區(qū),通過(guò)歷史數(shù)據(jù)或經(jīng)驗(yàn)預(yù)估監(jiān)測(cè)區(qū)的報(bào)警率和期望搜潛時(shí)間,在此基礎(chǔ)上,計(jì)算得到優(yōu)化的多反潛直升機(jī)部署地點(diǎn)以及分派策略(優(yōu)先順序),從而達(dá)到最大覆蓋警戒海域的目的。

多監(jiān)測(cè)區(qū)-多反潛機(jī)調(diào)度中存在的兩個(gè)問(wèn)題增加了優(yōu)化過(guò)程的復(fù)雜性:① 需要綜合考慮不同類型反潛直升機(jī)的調(diào)度,不同類型反潛直升機(jī)各有自己的優(yōu)缺點(diǎn),如何在反潛作戰(zhàn)中相互配合,揚(yáng)其所長(zhǎng),避其所短,是調(diào)度模型的一個(gè)重點(diǎn)研究目標(biāo);② 各反潛直升機(jī)的負(fù)載平衡問(wèn)題,如果調(diào)度(部署和分派策略)不當(dāng),將導(dǎo)致某一架或幾架直升機(jī)工作負(fù)載過(guò)大,這會(huì)加大飛行員的疲勞度,同時(shí)也會(huì)增加飛機(jī)的故障率,從而降低整個(gè)反潛系統(tǒng)的可用率。

因此,本文的主要?jiǎng)?chuàng)新點(diǎn)和貢獻(xiàn)是從水下網(wǎng)絡(luò)中心戰(zhàn)的角度出發(fā),為了達(dá)到最大覆蓋警戒海域的研究目標(biāo),特別針對(duì)以上兩個(gè)問(wèn)題,提出反潛直升機(jī)部署和調(diào)度模型及求解算法。在技術(shù)上借鑒了超立方體排隊(duì)模型來(lái)研究反潛調(diào)度問(wèn)題,該模型主要用于分析應(yīng)急服務(wù)系統(tǒng)的最大覆蓋問(wèn)題,在緊急醫(yī)療設(shè)施[28]和急救車部署[29-30]方面取得了較好的性能。但相關(guān)研究沒(méi)有考慮不確定情況下,多優(yōu)先級(jí)任務(wù)同時(shí)出現(xiàn)時(shí)不同服務(wù)單元的配置問(wèn)題,這也是本文的研究難點(diǎn)。

本文分兩個(gè)階段進(jìn)行優(yōu)化：第1階段建立超立方體反潛調(diào)度排隊(duì)模型和求解算法,以輸入的報(bào)警率作為事件到達(dá)率,期望響應(yīng)時(shí)間為平均服務(wù)時(shí)間,輸出校正系數(shù)和直升機(jī)繁忙概率等統(tǒng)計(jì)結(jié)果；第2階段建立反潛調(diào)度規(guī)劃模型,以第1階段輸出作為規(guī)劃模型的輸入?yún)?shù),從而得到反潛直升機(jī)的部署及分派策略。

1 排隊(duì)模型及求解算法

超立方體排隊(duì)模型是對(duì)常規(guī)排隊(duì)模型的擴(kuò)展,適用于應(yīng)急反應(yīng)系統(tǒng)的概率選址,反映了服務(wù)器與客戶需求之間更為復(fù)雜的分配關(guān)系,對(duì)多服務(wù)器排隊(duì)系統(tǒng)進(jìn)行了擴(kuò)展,增加了服務(wù)器狀態(tài)空間的描述,能夠更為精確地刻畫(huà)排隊(duì)系統(tǒng)的時(shí)空復(fù)雜性[31]。某一時(shí)間段內(nèi)服務(wù)器可能處于忙(1)或閑(0)兩種狀態(tài),系統(tǒng)中多個(gè)服務(wù)器的狀態(tài)組合構(gòu)成了系統(tǒng)狀態(tài)空間。例如狀態(tài){1,1,0}表達(dá)了一個(gè)三服務(wù)器系統(tǒng),其中1號(hào)服務(wù)器空閑,2號(hào)和3號(hào)服務(wù)器繁忙,從中可以看出一個(gè)三服務(wù)器系統(tǒng)的狀態(tài)空間可以由立方體的8個(gè)頂點(diǎn)表示,如果系統(tǒng)中服務(wù)器數(shù)量多于3個(gè),則稱之為超立方體。如果給定排隊(duì)系統(tǒng)參數(shù),例如服務(wù)器數(shù)量、分管區(qū)域數(shù)量(客戶數(shù)量)、客戶請(qǐng)求到達(dá)率以及服務(wù)器的平均服務(wù)時(shí)間等,利用超立方體排隊(duì)模型可以計(jì)算穩(wěn)態(tài)狀態(tài)下系統(tǒng)的一些關(guān)鍵性能指標(biāo),包括服務(wù)器的繁忙概率、工作強(qiáng)度以及服務(wù)器與客戶之間的分配概率等。

反潛直升機(jī)數(shù)量有限,在警戒海域執(zhí)行大規(guī)模搜潛任務(wù)時(shí)不能處于隨時(shí)可用狀態(tài),所以計(jì)算其繁忙概率對(duì)于應(yīng)召搜潛調(diào)度而言非常重要。另外直升機(jī)執(zhí)行反潛任務(wù)時(shí)存在協(xié)作關(guān)系,因此本文還需要計(jì)算校正系數(shù)來(lái)補(bǔ)償排隊(duì)系統(tǒng)中服務(wù)器無(wú)關(guān)性的假設(shè)。

1.1 參數(shù)描述

由于排隊(duì)模型中引用參數(shù)較多,在這里首先統(tǒng)一進(jìn)行描述。

nA(nB)為系統(tǒng)中A類型(B類型)反潛機(jī)數(shù)量，nA+nB=n;

LA(LB)為A類型(B類型)反潛機(jī)候選部署地點(diǎn)集合,LA+LB=L,包括多個(gè)地點(diǎn)l∈L;

Z為警戒海域內(nèi)監(jiān)測(cè)網(wǎng)格區(qū)集合,當(dāng)預(yù)警系統(tǒng)發(fā)現(xiàn)可疑目標(biāo)后,向反潛調(diào)度系統(tǒng)發(fā)出包括目標(biāo)大致所在區(qū)域(監(jiān)測(cè)區(qū)編號(hào))的警報(bào),由調(diào)度系統(tǒng)通知反潛機(jī)前往該區(qū)域?qū)嵤┧褲?包括多個(gè)監(jiān)測(cè)區(qū)z∈Z;

T為報(bào)警類別集合,本文把警報(bào)類別分為4類:① 只能分配給A型機(jī)處理的警報(bào)TA;② 如果A型機(jī)都忙時(shí),可以分配給B型機(jī)處理的警報(bào)TAB,即優(yōu)先要求A類型機(jī);③ 如果B型機(jī)都忙時(shí),可以分配給A型機(jī)處理的警報(bào)TBA;④ 只能分配給B型機(jī)處理的警報(bào)TB,TA、TB、TAB、TBA?T;

λt為類型為t的警報(bào)發(fā)生率,t∈T;

λtz為在第z監(jiān)測(cè)區(qū)發(fā)出的類型為t的警報(bào)發(fā)生率,λt=∑z∈Zλtz,z∈Z,t∈T;

λ為系統(tǒng)級(jí)所有警報(bào)發(fā)生率,λ=∑z∈Z,t∈Tλtz;

τAlzt(τBlzt)為l地點(diǎn)的A型機(jī)(B型機(jī))執(zhí)行z區(qū)t類別警報(bào)任務(wù)的期望響應(yīng)時(shí)間(從收到警報(bào)到抵達(dá)目標(biāo)區(qū)域z所花費(fèi)的時(shí)間);

τA(τB)為系統(tǒng)級(jí)A型機(jī)(B型機(jī))的平均響應(yīng)時(shí)間;

θt為t類別警報(bào)對(duì)響應(yīng)時(shí)間指標(biāo)的要求,本文以規(guī)定響應(yīng)時(shí)間內(nèi)得到處理的t類別警報(bào)數(shù)量比例來(lái)作為衡量指標(biāo);

ρA(ρB)為A型機(jī)(B型機(jī))的工作強(qiáng)度;

pAl(pBl)為部署在l地點(diǎn)的A型機(jī)(B型機(jī))的繁忙概率;

pA(pB)為系統(tǒng)級(jí)A型機(jī)(B型機(jī))的繁忙概率,pA=(∑l∈LApAl)/nA;

aAztk(aBztk)為A型機(jī)(B型機(jī))的排序調(diào)度表,表示當(dāng)z區(qū)域發(fā)出t類別警報(bào)時(shí),作為第k順位分配的反潛機(jī)地點(diǎn),如果t∈TA∪TAB,k=1,2， …,nA,如果t∈TBA,則k=nB+1,nB+2,…,nB+nA;

dAlztk(dBlztk)為A型機(jī)(B型機(jī))的分配概率,表示當(dāng)z區(qū)域發(fā)出t類別警報(bào)時(shí),作為第k順位的l地點(diǎn)A型機(jī)(B型機(jī))被分派到處理該警報(bào)的概率;

PAk(PBk)為由于k個(gè)A型機(jī)(B型機(jī))的繁忙,造成警報(bào)無(wú)法處理的概率,即損失率;

RAlzt(RBlzt)為表示當(dāng)z區(qū)域發(fā)出t類別警報(bào)時(shí),l地點(diǎn)的A型機(jī)(B型機(jī))能夠在規(guī)定響應(yīng)時(shí)間(與警報(bào)類別有關(guān))內(nèi)到達(dá)z區(qū)域執(zhí)行搜潛任務(wù)的概率;

QA(nA,ρA,k)或QB(nB,ρB,k)為超立方體模型的校正系數(shù),可看作輸入?yún)?shù)為反潛機(jī)數(shù)量nA(nB),反潛機(jī)工作強(qiáng)度ρA(ρB)以及派遣順位k的函數(shù)。

本文假設(shè)z區(qū)域t類型警報(bào)的到達(dá)率服從參數(shù)為λtz的泊松分布,如果某警報(bào)需要多架反潛機(jī)協(xié)同反潛,可以發(fā)出多個(gè)虛擬警報(bào)來(lái)處理。

1.2 排隊(duì)模型

為了避免重復(fù),本文中主要描述A型機(jī)的數(shù)學(xué)模型,B型機(jī)的模型與A型機(jī)相同。

1.2.1 校正系數(shù)的計(jì)算模型

當(dāng)出現(xiàn)警報(bào)時(shí),在前面連續(xù)j架反潛機(jī)都繁忙,第k+1順位的反潛機(jī)空閑的概率是在排隊(duì)論服務(wù)器無(wú)關(guān)性假設(shè)下計(jì)算得到的,在反潛系統(tǒng)中,反潛機(jī)的調(diào)度往往存在協(xié)作關(guān)系,因此需要在模型中加入無(wú)關(guān)性校正系數(shù)QA(nA,ρA,k)來(lái)修正[32]。設(shè)PAnA表示所有的nA架A型反潛機(jī)繁忙的概率,PA0表示沒(méi)有A型機(jī)繁忙,即所有A型機(jī)空閑的概率,設(shè)k=0,1,…,nA-1,則有

QA(nA,ρA,k)=

(1)

1.2.2 分配概率的計(jì)算模型

在排隊(duì)模型中,如果要計(jì)算服務(wù)器的繁忙概率,必須先計(jì)算服務(wù)器-客戶對(duì)的分配概率。在本文中,服務(wù)器表示為反潛機(jī),客戶可表示為監(jiān)測(cè)區(qū),所以分配概率既是dAlztk或dBlztk的值。分配概率與排序調(diào)度表aztk有關(guān),如果a1A2=3,表示1號(hào)區(qū)域發(fā)出只能由A型機(jī)反潛的警報(bào)時(shí),位于3號(hào)地點(diǎn)的反潛機(jī)將作為第2順位被分配去處理。另外分配概率與警報(bào)類別有關(guān),因?yàn)閷?duì)于TBA類別的警報(bào),必須要B型機(jī)都繁忙時(shí)才輪到分配A型機(jī)去處理,所以對(duì)于處理TA∪TAB類別的警報(bào)和處理TBA類別的警報(bào),其分配概率的計(jì)算是不同的。特別是當(dāng)t∈TB時(shí),因?yàn)樵擃惥瘓?bào)只能由B型機(jī)處理,所以這時(shí)dAlztk=0。綜上所述,dAlztk的計(jì)算公式分以下兩種情況:

(1)如果t∈TA∪TAB,z∈Z,k=1,2,…,nA,l∈LA,則

(2)

(2)如果t∈TBA,z∈Z,k=nB+1,nB+2,…,nB+nA,l∈LA,則

(3)

1.2.3 繁忙概率的計(jì)算模型

地點(diǎn)l的A型機(jī)繁忙概率pAl計(jì)算公式為

(4)

式中,rAl表示地點(diǎn)l的A型機(jī)分配執(zhí)行搜索任務(wù)的頻率,計(jì)算公式為

(5)

式(5)中第二項(xiàng)表示地點(diǎn)l的A型機(jī)分派給TBA類別警報(bào)任務(wù)的頻率,同第1.2.2節(jié)分派概率的計(jì)算,加入PBnB項(xiàng)的原因是,在t∈TBA時(shí),所有的B型機(jī)都繁忙時(shí),才只好找A型機(jī)調(diào)度。而且這時(shí)A型機(jī)的分派優(yōu)先級(jí)順序都排在B型機(jī)后面,所以優(yōu)先序號(hào)j從nB+1開(kāi)始。

1.2.4 平均響應(yīng)時(shí)間的計(jì)算模型

平均響應(yīng)時(shí)間的計(jì)算公式如下:

(6)

式中,平均響應(yīng)時(shí)間是單架機(jī)響應(yīng)時(shí)間τAlzt與歸一化的分派概率dAlztk/(1-PAnA)的乘積,乘積結(jié)果再用警報(bào)到達(dá)頻率(λtz/λ)加權(quán)。

1.3 排隊(duì)模型的求解算法

由于穩(wěn)態(tài)時(shí)第1.2節(jié)的排隊(duì)模型無(wú)法求得精確解,只能采用近似算法迭代計(jì)算。近似算法步驟如下。

算法輸入?yún)?shù):①λtz,警報(bào)發(fā)生率;②τAlzt,l地點(diǎn)的A型機(jī)響應(yīng)時(shí)間;③aAztk,可由第2節(jié)規(guī)劃模型求解算法輸入。

初始化:初始化平均響應(yīng)時(shí)間為

(7)

按埃爾朗損失率公式初始化所有A型機(jī)繁忙的概率,即

(8)

式中,A型機(jī)初始空閑概率PA0,0為

(9)

設(shè)Lz為距離監(jiān)測(cè)區(qū)z最近的反潛機(jī)地點(diǎn)集合,則反潛機(jī)的初始繁忙概率為

(10)

迭代計(jì)數(shù)器h設(shè)為0。

步驟 1利用上一次迭代計(jì)算的平均服務(wù)時(shí)間τA,h-1(如果為第1次迭代,則利用初始化的τA,0)更新?lián)p失率,計(jì)算公式為

(11)

式中，A型機(jī)空閑概率為

(12)

步驟 2設(shè)k=0,1,…,nA-1,按式(1)計(jì)算本次迭代的校正系數(shù)QA,h(nA,ρA,h,k),其中工作強(qiáng)度ρA,h=λτA,h-1/nA。

步驟 3按式(2)和式(3)計(jì)算本次迭代h的分配概率dAlztk,h,需要的參數(shù)中繁忙概率使用上一迭代的計(jì)算值,校正系數(shù)取步驟2的計(jì)算結(jié)果。

步驟 4按式(6)計(jì)算本次迭代的平均響應(yīng)時(shí)間τA,h,需要的參數(shù)中分派概率使用步驟4的計(jì)算結(jié)果,損失率使用步驟1的計(jì)算結(jié)果。

步驟 5按式(4)和式(5)計(jì)算當(dāng)前h迭代時(shí)的反潛機(jī)繁忙概率pAl,h,需要的參數(shù)中優(yōu)先級(jí)排序的繁忙概率使用上一迭代的計(jì)算值,校正系數(shù)取步驟2的計(jì)算結(jié)果。

步驟 6若對(duì)于所有地點(diǎn)的A類反潛機(jī),其前后兩次迭代的計(jì)算結(jié)果大于某一閾值ε,即|pAl,h-pAl,h-1|>ε,其中l(wèi)∈LA,則跳轉(zhuǎn)步驟1繼續(xù)迭代,否則返回繁忙概率和校正系數(shù)。

2 規(guī)劃模型

反潛機(jī)調(diào)度規(guī)劃模型的目標(biāo)函數(shù)是最大化規(guī)定時(shí)間內(nèi)完成高優(yōu)先級(jí)反潛任務(wù)的數(shù)量比例,并且在此同時(shí)實(shí)現(xiàn)對(duì)低優(yōu)先級(jí)警報(bào)的一定(給定閾值內(nèi))覆蓋率,這樣折中實(shí)現(xiàn)對(duì)所有警報(bào)的及時(shí)處理。

2.1 排隊(duì)模型相關(guān)參數(shù)的定義

排隊(duì)模型近似求解算法的輸出用作規(guī)劃模型的輸入?yún)?shù)。規(guī)劃模型目標(biāo)函數(shù)和約束條件中用到的兩個(gè)系數(shù)與超立方體排隊(duì)模型輸出相關(guān)。一個(gè)系數(shù)捕捉系統(tǒng)分派策略滿足時(shí)限要求的性能,本文稱之為及時(shí)率;另一個(gè)系數(shù)捕捉系統(tǒng)的工作負(fù)載。

及時(shí)率系數(shù)cAlztk或cBlztk表示當(dāng)z區(qū)域發(fā)出t類別警報(bào)時(shí),作為第k順位的l地點(diǎn)的A型機(jī)(B型機(jī))能夠在規(guī)定響應(yīng)時(shí)間(與警報(bào)類別有關(guān))內(nèi)到達(dá)z區(qū)域的數(shù)量比例,為反映監(jiān)測(cè)區(qū)的優(yōu)先級(jí),本文使用監(jiān)測(cè)區(qū)z的警報(bào)率(λtz/λt)對(duì)及時(shí)率系數(shù)加權(quán)。假設(shè)同一類型反潛機(jī)的繁忙概率相同時(shí),及時(shí)率的計(jì)算類似于分派概率,其計(jì)算公式為

(1)如果t∈TA∪TAB,z∈Z,k=1,2,…,nA,l∈LA,則

cAlztk=(λtz/λt)RAlzt[QA(nA,ρA,k)(1-pA)(pA)k-1]

(13)

(2)如果t∈TBA,z∈Z,k=nB+1,nB+2,…,nB+nA,l∈LA,則

cAlztk=(λtz/λt)RAlztPBnB[QA(nA,ρA,k)(1-pA)(pA)k-1]

(14)

式中,校正系數(shù)與A型機(jī)繁忙概率pA都為超立方體排隊(duì)模型的計(jì)算結(jié)果。及時(shí)率系數(shù)計(jì)算過(guò)程中考慮了兩個(gè)不確定因素的影響:① 在校正系數(shù)補(bǔ)償下反潛機(jī)可用性的不確定性,比地點(diǎn)l更優(yōu)先的k個(gè)反潛機(jī)繁忙時(shí)至少有一個(gè)反潛機(jī)空閑的概率為(1-pA)(pA)k-1,再乘以校正系數(shù)QA,表示有k-1架反潛機(jī)優(yōu)先于l調(diào)度,但l仍然被調(diào)度的概率;②RAlzt表示飛行和搜潛時(shí)間的不確定性。

同樣地,負(fù)載系數(shù)gAlztk或gBlztk,分別表示A型機(jī)或B型機(jī)的工作負(fù)載,A型機(jī)負(fù)載系數(shù)計(jì)算公式為

(1)如果t∈TA∪TAB,z∈Z,k=1,2,…,nA,l∈LA,則

gAlztk=(λtz/λt)τAlzt[QA(nA,ρA,k)(1-pA)(pA)k-1]

(15)

(2)如果t∈TBA,z∈Z,k=nB+1,nB+2,…,nB+nA,l∈LA,則

gAlztk=(λtz/λt)τAlztPBnB[QA(nA,ρA,k)(1-pA)(pA)k-1]

(16)

式中,τAlzt表示l地點(diǎn)的A型機(jī)執(zhí)行z區(qū)t類別警報(bào)任務(wù)的期望響應(yīng)(服務(wù))時(shí)間。B型機(jī)計(jì)算公式類似,負(fù)載系數(shù)用于平衡同一類型反潛機(jī)的工作負(fù)載。

2.2 規(guī)劃模型的定義

首先定義反潛規(guī)劃模型的決策變量,為了完成反潛調(diào)度的兩個(gè)主要目標(biāo),反潛規(guī)劃模型相應(yīng)也包含兩套決策變量。布爾變量hAl和hBl決定A型機(jī)和B型機(jī)的部署地點(diǎn),如果hAl=1,表示有A型機(jī)部署在l地點(diǎn)(l∈LA),否則hAl為0。hBl與hAl相似,如果hBl=1,表示有B型機(jī)部署在l地點(diǎn)(l∈LB),否則hBl為0。

第二套決策變量為mAlztk和mBlztk也為布爾變量,表示A型機(jī)和B型機(jī)分派至監(jiān)測(cè)區(qū)的排序調(diào)度表,與aAztk或aBztk可以相互轉(zhuǎn)化,一方面用于表示對(duì)于不同類別警報(bào)的分派順序,例如TA類別警報(bào)任務(wù)只能由A型機(jī)執(zhí)行,TAB類別警報(bào)任務(wù)先由A型機(jī)執(zhí)行,如果沒(méi)有A型機(jī)空閑時(shí)才交由B型機(jī)執(zhí)行,TB和TBA也是類似。另外,也可以表示當(dāng)發(fā)生某警報(bào)時(shí),同一類型反潛機(jī)內(nèi)部的分派順序,如果mAlztk=1,表示當(dāng)z區(qū)域發(fā)出t類別警報(bào)時(shí),地點(diǎn)l的A型機(jī)作為第k順位分配,如果t∈TA∪TAB,k=1,2,…,nA,如果t∈TBA,則k=nB+1,nB+2， …,nB+nA;其他情況時(shí)mAlztk=0。mBlztk的設(shè)計(jì)與mAlztk相似。

以上決策變量的設(shè)計(jì)把反潛調(diào)度規(guī)劃轉(zhuǎn)化為二元線性規(guī)劃模型,目標(biāo)函數(shù)為

(17)

該目標(biāo)函數(shù)的目的是最大化高優(yōu)先級(jí)警報(bào)TA的期望覆蓋率,在反潛調(diào)度規(guī)劃模型中,只有A型機(jī)可以執(zhí)行TA警報(bào)任務(wù)。

約束關(guān)系如下:

(18)

nB+1,nB+2,…,nB+nA

(19)

(20)

k=nA+1,k=nA+2,…,nA+nB

(21)

(22)

(23)

(24)

(25)

(26)

(27)

(28)

(29)

(30)

mAlztk≤hAl,?l∈LA,z∈Z,t∈TA∪TAB,

k=1,2,…,nA

(31)

mAlztk≤hAl,?l∈LA,z∈Z,t∈TBA,k=nB+1,

nB+2,…,nB+nA

(32)

mBlztk≤hBl,?l∈LB,z∈Z,t∈TB∪TBA,

k=1,2，…,nB

(33)

mmlztk≤hBl,?l∈LB,z∈Z,t∈TAB,k=nA+1,

nA+2,…,nA+nB

(34)

(35)

(36)

(37)

(38)

式(18)表示對(duì)于警報(bào)t∈TA∪TAB,在排序調(diào)度表中同一順位只能選擇一架A型機(jī)執(zhí)行。式(19)表示對(duì)于警報(bào)t∈TBA,首先分派B型機(jī)執(zhí)行,如果B型機(jī)都繁忙才分派給A型機(jī),而且在排序調(diào)度表中同一優(yōu)先級(jí)順位只能選擇一架A型機(jī)執(zhí)行。式(20)和式(21)與式(18)和式(19)表示的意義相同,但針對(duì)的是B型機(jī)。式(22)～式(25)表示某反潛機(jī)l不能同時(shí)是某監(jiān)測(cè)區(qū)z的多個(gè)順位分派機(jī),比如不能既是z的第1順位分派機(jī)又是z的第2順位分派機(jī)。式(26)～式(28)保證低優(yōu)先級(jí)的警報(bào)(TAB,TBA,TB)的期望覆蓋率大于某個(gè)閾值θt。式(29)和式(30)為反潛機(jī)的數(shù)量約束。式(31)～式(34)表示只有已部署到指定地點(diǎn)的反潛機(jī)才能被調(diào)度。式(35)和式(36)約束同一型號(hào)反潛機(jī)之間負(fù)載平衡,即負(fù)載差不高于某一閾值δA(δB)。式(37)和式(38)約束反潛機(jī)的分派監(jiān)測(cè)區(qū)相鄰,即如果z的鄰居z′分派給某反潛機(jī)l分管,則z才能被該反潛機(jī)分管。

本文中,由式(18)～式(34)構(gòu)成的模型為基本調(diào)度模型,基本模型增加了式(35)和式(36)后為負(fù)載平衡模型,基本模型增加了式(37)和式(38)后增加了分管區(qū)域相鄰的約束。所有約束構(gòu)成了加強(qiáng)模型。另外假設(shè)反潛機(jī)在我方警戒海域可以就近補(bǔ)充油料,所以約束關(guān)系中沒(méi)有包括油料限制。

3 調(diào)度算法及應(yīng)用場(chǎng)景

由于排隊(duì)模型求解算法中反潛機(jī)的繁忙概率與分派策略aAztk(aBztk)有關(guān),而aAztk(aBztk)可以通過(guò)mAlztk(mBlztk)轉(zhuǎn)換得到,規(guī)劃模型又需要排隊(duì)模型中的校正系數(shù)和繁忙概率,所以本文設(shè)計(jì)了迭代算法求解規(guī)劃模型,其計(jì)算步驟如下。

步驟 1初始化;

步驟 1.1利用“分配最近空閑反潛機(jī)”策略初始化初始排序調(diào)度表aAztk(aBztk);

步驟 1.2調(diào)用第1.3節(jié)排隊(duì)模型求解算法得到繁忙概率和校正系數(shù);

步驟 1.3設(shè)置初始不平衡容忍度δ0=1;

步驟 1.4迭代次數(shù)i設(shè)為1;

步驟 2利用上次迭代(或初始階段)生成的繁忙概率和校正系數(shù)求解規(guī)劃模型,得到新的排序調(diào)度表aAztk,i;

步驟 3利用新的排序調(diào)度表aAztk,i,調(diào)用排隊(duì)模型求解算法,得到新的繁忙概率pA,i和校正系數(shù)Qi,設(shè)置δA,i=1/2(maxpA,i-minpA,i);

步驟 4算法包括2個(gè)迭代停止條件。第1個(gè)停止條件為式(33)和式(34)無(wú)法滿足,這時(shí)再減少不平衡容忍度δ會(huì)導(dǎo)致無(wú)法得到可行解;第2個(gè)停止條件為繁忙概率小于某個(gè)閾值。如果算法停止條件都沒(méi)有達(dá)到,增加迭代次數(shù)i并跳轉(zhuǎn)至步驟2;

步驟 5返回最新迭代計(jì)算的排序調(diào)度表和反潛機(jī)部署地點(diǎn)。

具體過(guò)程如圖1所示。

圖1 柵格化反潛警戒區(qū)示意圖Fig.1 Schematic diagram of grid antisubmarine warning area

本算法的應(yīng)用時(shí)機(jī)為需要對(duì)某海域布防警戒前,具體應(yīng)用步驟描述如下。

步驟 1類似圖1,把需要布防海域柵格化為若干監(jiān)測(cè)區(qū)z,同時(shí)篩選出適合部署反潛直升機(jī)的候選地點(diǎn)l;

步驟 2通過(guò)歷史數(shù)據(jù)、地點(diǎn)重要程度或敵方威脅程度預(yù)估各個(gè)監(jiān)測(cè)區(qū)的報(bào)警率λtz,并根據(jù)l和z之間距離和反潛機(jī)速度估算出期望響應(yīng)時(shí)間τAlzt(或τBlzt),如果可能也可以估計(jì)期望響應(yīng)時(shí)間的不確定性參數(shù)RAlzt或RBlzt(缺省為1),θt根據(jù)相關(guān)條例推導(dǎo)指定;

步驟 3步驟1和步驟2的結(jié)果作為調(diào)度算法輸入?yún)?shù)進(jìn)行迭代計(jì)算。由于調(diào)度算法的步驟1需要計(jì)算第2.2節(jié)的規(guī)劃模型(二元線性規(guī)劃),而該規(guī)劃模型中用到的及時(shí)率和工作負(fù)載2個(gè)系數(shù)需要用到排隊(duì)模型輸出的繁忙概率和校正系數(shù),所以需要嵌入調(diào)用第1.3節(jié)的排隊(duì)模型求解算法;

步驟 4當(dāng)調(diào)度算法退出時(shí),生成兩個(gè)決策變量的計(jì)算結(jié)果,反潛機(jī)的部署地點(diǎn)(hAl和hBl)以及排序調(diào)度表(mAlztk和mBlztk)。

指揮人員可以以本算法結(jié)果作為決策支持建議,在推薦地點(diǎn)部署反潛機(jī)搭載平臺(tái)(驅(qū)逐艦或護(hù)衛(wèi)艦)。如果某監(jiān)測(cè)區(qū)出現(xiàn)威脅目標(biāo)時(shí),能夠以排序調(diào)度表作為調(diào)度策略基礎(chǔ)方案,按優(yōu)先級(jí)次序(如果高優(yōu)先級(jí)直升機(jī)已執(zhí)行任務(wù),則依次派遣低優(yōu)先級(jí)直升機(jī))派遣對(duì)應(yīng)反潛機(jī)執(zhí)行對(duì)應(yīng)任務(wù)。

4 仿真實(shí)驗(yàn)

本文仿真背景案例采用了類似圖1的柵格化警戒海域背景,整個(gè)海域分為85個(gè)監(jiān)測(cè)區(qū),其中32個(gè)區(qū)域可以作為反潛直升機(jī)搭載平臺(tái)的候選部署地點(diǎn),假設(shè)有10架A類機(jī)和4架B類機(jī)可供調(diào)度。從安全角度考慮,本文不討論警戒海域和反潛機(jī)的具體參數(shù)。

在參數(shù)設(shè)置方面,期望響應(yīng)時(shí)間τAlzt以小時(shí)為單位,按反潛機(jī)的部署地點(diǎn)l和報(bào)警區(qū)域z之間的距離除以反潛機(jī)的速度計(jì)算。τBlzt的計(jì)算與A類機(jī)類似,但考慮B類機(jī)可能采取雙機(jī)協(xié)同,且綜合性能較弱,所以其計(jì)算得到的響應(yīng)時(shí)間乘以2作為懲罰因子。警報(bào)類別的概率分布為:TA為75%,TAB為20%,TBA為5%,仿真中不生成TB警報(bào)。系統(tǒng)級(jí)所有警報(bào)發(fā)生率λ=3 h,每一個(gè)監(jiān)測(cè)區(qū)域的警報(bào)發(fā)生率相等。

當(dāng)前算法在聯(lián)想筆記本上即可運(yùn)行,CPU為四核Intel-i5,基礎(chǔ)頻率1.6 GHZ,內(nèi)存8 G,算法在VS2012下用C++實(shí)現(xiàn),規(guī)劃求解部分調(diào)用IBM ILOG CPLEX 12.9。容忍度設(shè)置為0.01時(shí),平均運(yùn)行時(shí)間統(tǒng)計(jì)如下:基本模型耗時(shí)約0.6 min,負(fù)載平衡模型耗時(shí)約6.2 min,加強(qiáng)模型耗時(shí)約11.3 min,3個(gè)規(guī)劃模型中調(diào)用排隊(duì)算法更新繁忙概率和校正系數(shù)的時(shí)間分別為0.1 min、0.65 min和0.6 min。

圖2和圖3為部署了全部10架A型機(jī)和4架B型機(jī)(B型機(jī)可以為雙機(jī)協(xié)同)時(shí),在迭代調(diào)度算法中調(diào)用排隊(duì)模型求解算法得到的繁忙概率pAl(pBl),其中,l表示10架A型機(jī)和4架B型機(jī)的部署地點(diǎn)。圖2中A型機(jī)的繁忙概率最低為0.052,最高為0.105,平均為0.075;圖3中B型機(jī)的繁忙概率最低為0.225,最高為0.357,平均為0.310,所有機(jī)型的繁忙概率不超過(guò)其平均值的5%。

圖2 不同部署地點(diǎn)的A型機(jī)繁忙概率Fig.2 Busy probability of type A aircraft in different deployment locations

圖3 不同部署地點(diǎn)的B型機(jī)繁忙概率Fig.3 Busy probability of type B aircraft in different deployment locations

圖4和圖5為排隊(duì)模型求解算法得到的不同派遣順位對(duì)應(yīng)的校正系數(shù)QA(nA,ρA,k)和QB(nB,ρB,k),其中nA=10,nB=4,ρA和ρB的計(jì)算見(jiàn)排隊(duì)模型求解算法的步驟2。

圖4 A型機(jī)不同派遣順位對(duì)應(yīng)的校正系數(shù)Fig.4 Correction coefficient of different dispatch order for type A aircraft

圖5 B型機(jī)不同派遣順位對(duì)應(yīng)的校正系數(shù)Fig.5 Correction coefficient of different dispatch order for type B aircraft

從圖4和圖5中可以看出QA(nA,ρA,0)=QB(nB,ρB,0)=1,這是因?yàn)楫?dāng)所有反潛機(jī)空閑時(shí),能夠派遣反潛機(jī)的概率為1。

圖6統(tǒng)計(jì)了TA類別警報(bào)發(fā)出后,分派A型機(jī)處理時(shí)得到的目標(biāo)函數(shù)值,即TA警報(bào)的覆蓋率。在這個(gè)仿真實(shí)驗(yàn)中,B型機(jī)數(shù)量固定為4,即nB=4,改變A型機(jī)數(shù)量nA,從中分析當(dāng)增加或減少A型機(jī)數(shù)量時(shí),調(diào)度算法的性能變化情況,除此以外,圖6中還對(duì)警報(bào)發(fā)生率不同時(shí),隨A型機(jī)數(shù)量變化的算法性能進(jìn)行了敏感性分析。

圖6 不同警報(bào)發(fā)生率下的目標(biāo)函數(shù)值敏感性分析Fig.6 Sensitivity analysis of objective function value under different alarm rates

圖6的結(jié)果可以用于幫助指揮人員判別需要多少數(shù)量的A型機(jī)才能達(dá)到指定的監(jiān)測(cè)覆蓋率指標(biāo)。例如,當(dāng)警報(bào)發(fā)生率λ介于2.5和3.5之間時(shí),如果需要在警戒海域達(dá)到90%以上的警報(bào)覆蓋率,則需要部署6架A型機(jī)即可滿足要求,但如果想在同樣的警報(bào)發(fā)生率情況下達(dá)到95%的覆蓋率,A型機(jī)數(shù)量至少需要達(dá)到9架。而且從圖6中還可以發(fā)現(xiàn)隨著A型機(jī)數(shù)量增加,目標(biāo)函數(shù)值的增加逐漸減緩。例如當(dāng)λ=3.0時(shí),A型機(jī)數(shù)量從2架增加到5架時(shí),覆蓋性能呈現(xiàn)一個(gè)較大的增長(zhǎng),從53%增加到了90%,而當(dāng)A型機(jī)數(shù)量從8架增加到10架時(shí),性能增長(zhǎng)幅度較小,從94%增加到95%。

下面比較本文提出的優(yōu)化算法與常規(guī)“就近派遣”策略的性能。當(dāng)有警報(bào)發(fā)出后,就近派遣策略按離警報(bào)區(qū)域的距離由近到遠(yuǎn)調(diào)度,很明顯該策略的響應(yīng)時(shí)間最短,因此覆蓋率也相應(yīng)最大。但警報(bào)率較高時(shí)可能導(dǎo)致反潛機(jī)處于繁忙狀態(tài),這時(shí)需要根據(jù)全局態(tài)勢(shì)從長(zhǎng)考慮,合理的做法可能是預(yù)留最近的反潛機(jī)以備更高等級(jí)的警報(bào)。其中就近派遣策略的實(shí)現(xiàn)步驟如下。

步驟 1根據(jù)候選部署點(diǎn)與85個(gè)監(jiān)測(cè)區(qū)之間的平均距離,利用分簇排序方法選擇了A型機(jī)和B型機(jī)的部署點(diǎn);

步驟 2按照監(jiān)測(cè)區(qū)與部署點(diǎn)的距離遠(yuǎn)近生成排序調(diào)度表;

步驟 3在步驟1和步驟2的基礎(chǔ)上,設(shè)置2套決策變量hAl、hBl、mAlztk和mBlztk的相應(yīng)值為1,按式(7)計(jì)算覆蓋率。

圖7顯示了優(yōu)化算法與就近策略的警報(bào)覆蓋率性能指標(biāo),其中B型機(jī)數(shù)量設(shè)置為4,警報(bào)率λ=3,如圖7所示,就近派遣策略也能取得不錯(cuò)的性能,但總體上優(yōu)化算法比就近派遣策略占優(yōu),在反潛機(jī)數(shù)量較小時(shí)(nB=4,nA=2)就近派遣策略與優(yōu)化算法性能相同。進(jìn)一步比較當(dāng)警報(bào)率變化時(shí),從λ=1增加到10,優(yōu)化算法與就近派遣的性能差異。

圖7 就近派遣與優(yōu)化算法的性能對(duì)比(不同A型機(jī)數(shù)量)Fig.7 Performance comparison between the closest dispatch policy and optimal algorithm (different number of type A aircraft)

同樣地就近派遣策略的警報(bào)覆蓋率也較高,圖8中可以看出優(yōu)化算法最大高于就近派遣1.71%,仿真中參數(shù)設(shè)置為nB=4,nA=10。盡管就近派遣策略在警報(bào)覆蓋率上與優(yōu)化算法的差異不大,但優(yōu)化算法的優(yōu)勢(shì)在于構(gòu)建排序調(diào)度表時(shí)避免了局部反潛機(jī)過(guò)忙的影響。例如在某些情況下(指揮人員憑經(jīng)驗(yàn)難以識(shí)別),派遣較遠(yuǎn)的反潛機(jī)與就近派遣相比,可以在整體反潛性能不變的情況下,降低反潛機(jī)的繁忙率,從而延長(zhǎng)整個(gè)反潛機(jī)系統(tǒng)的生命周期。同時(shí),通過(guò)把任務(wù)盡可能平均分配到所有反潛機(jī)上,如圖2和圖3所示,可以提高反潛機(jī)的空閑率,從而較少警報(bào)區(qū)域的等待時(shí)間,最終增加了警報(bào)的期望覆蓋率。而且優(yōu)化算法能夠產(chǎn)生一個(gè)合理的排序調(diào)度表,可以幫助指揮人員判別合適的預(yù)備反潛機(jī)。

圖8 就近派遣與優(yōu)化算法的性能(λ不同)Fig.8 Performance comparison between the closest dispatch policy and optimal algorithm (different λ)

5 結(jié) 論

當(dāng)海上傳感器平臺(tái)監(jiān)測(cè)到水下可疑目標(biāo)并發(fā)出警報(bào)后,如何快速有效地派遣反潛機(jī)至警報(bào)區(qū)域是決定反潛成敗的關(guān)鍵問(wèn)題,如果響應(yīng)過(guò)慢則會(huì)丟失目標(biāo)而貽誤戰(zhàn)機(jī)。為此,本文提出了一個(gè)二元整數(shù)規(guī)劃模型及相應(yīng)求解算法,其主要目標(biāo)包括:① 解決目前兩種主力反潛直升機(jī)的優(yōu)化部署問(wèn)題;② 構(gòu)造排序調(diào)度表以解決反潛機(jī)的分派問(wèn)題。其中,如何利用有限的反潛機(jī)資源警戒大范圍海域是本文的一個(gè)難點(diǎn),這也是應(yīng)召反潛的重點(diǎn)問(wèn)題。為此,建立了一個(gè)排隊(duì)模型計(jì)算直升機(jī)繁忙概率和無(wú)關(guān)性校正系數(shù),以此作為規(guī)劃模型的輸入?yún)?shù),以便更加準(zhǔn)確地反映當(dāng)前態(tài)勢(shì),提高規(guī)劃模型的精確度。

盡管該模型以兩種反潛機(jī)為例進(jìn)行說(shuō)明,但可以利用相似方法擴(kuò)展到多類型反潛機(jī)優(yōu)化調(diào)度上。當(dāng)多處同時(shí)發(fā)現(xiàn)目標(biāo)需要派遣多機(jī)反潛時(shí),也可以考慮應(yīng)用本文研究成果進(jìn)行進(jìn)一步擴(kuò)展。另外,警報(bào)類別或優(yōu)先級(jí)的劃分對(duì)調(diào)度模型的性能影響較大,盡管當(dāng)前有非常多的有關(guān)類別或優(yōu)先級(jí)計(jì)算的模型和算法,但尚缺乏反潛領(lǐng)域應(yīng)用的相關(guān)研究,這也是下一步工作的重點(diǎn)方向。