張 昕,苑德春,潘 璠
(1.軍事交通學院 聯(lián)合投送系,天津300161;2.軍事交通學院 軍事交通運輸研究所,天津300161)
戰(zhàn)略投送能力是衡量一個國家未來作戰(zhàn)能力的重要指標,而戰(zhàn)略空運則是國家戰(zhàn)略投送的重要組成部分,是奪取軍事制高點的戰(zhàn)略籌碼。一方面,隨著國家綜合實力的增強以及戰(zhàn)略格局的宏觀運籌,發(fā)生外敵對我軍實施大規(guī)模陸上入侵的可能性越來越小,但在與國家核心利益密切相關(guān)的熱點方向,發(fā)生局部戰(zhàn)爭的可能性增大。另一方面,隨著國家利益拓展,能源與資源的穩(wěn)定供給對國家經(jīng)濟的長遠發(fā)展影響重大;頻發(fā)的重大自然災害對人民生命財產(chǎn)構(gòu)成嚴重威脅;社會轉(zhuǎn)型期的矛盾突顯,突發(fā)公共事件對社會穩(wěn)定構(gòu)成嚴重沖擊。面對傳統(tǒng)安全威脅與非傳統(tǒng)安全威脅,建設(shè)強大的航空戰(zhàn)略投送力量十分迫切。
一直以來,國內(nèi)對航空戰(zhàn)略投送兵力分配的研究停留在點對點式的直達航線,投送兵力的規(guī)模想定也偏小。但在實際的投送過程中,能夠稱之為戰(zhàn)略投送的,必須是遠距離、大規(guī)模的多種機型混合運輸。因此,本文主要研究帶延誤的多個機場航線網(wǎng)絡中多種機型的兵力分配模型。
在戰(zhàn)略空運過程中,假設(shè)以軍、民航大型機場為空運平臺,現(xiàn)有開通航線為投送網(wǎng)絡,飛行計劃以戰(zhàn)略空運任務為第一優(yōu)先,同時考慮空域流量的影響和人員在中轉(zhuǎn)機場整備集結(jié)再次登機的時間延誤,所有分配的兵力任務都試圖選擇最短航路到達目的地。
給定帶延誤的多機型航空戰(zhàn)略投送兵力分配模型參數(shù)(見表1)。

表1 帶延誤的多機型航空戰(zhàn)略投送兵力分配模型參數(shù)

中轉(zhuǎn)航段a和b在連接機場網(wǎng)絡中連接機場ω 的航段k上,當且僅當航段a和b同樣在該航段上,由此滿足以下條件:

航段投送量滿足以下條件:

航段投送量與航路投送量滿足以下關(guān)系[2]:

考慮中轉(zhuǎn)機場的時間延誤,中轉(zhuǎn)流量與航段流量需滿足以下關(guān)系:

空域流量阻抗則可以修正為航段阻抗與中轉(zhuǎn)延誤阻抗2 個部分之和,即

式(7)表明,機場的空域阻抗是由連接該機場的各條航段的流量阻抗及各個中轉(zhuǎn)方向的時間延誤阻抗累計而成。

也就是說,機型m在航段a上的投送量對機型l在航段b上的流量阻抗影響一般不等于機型l在航段b上的投送量對機型m在航段a上的流量阻抗的影響。
同理,機型m在中轉(zhuǎn)機場c→d上的投送量對機型l在中轉(zhuǎn)機場g→h上的延誤影響一般不等于機型l在中轉(zhuǎn)機場g→h上的投送量對機型m在中轉(zhuǎn)機場c→d上的延誤影響。需要指出,這里將時間延誤t和空域阻抗d分開處理,并認為前者是航段流量x的函數(shù)t(x),后者是中轉(zhuǎn)流量y的函數(shù)d(y)。由式(3)、(5)、(6)可見,x和y之間具有內(nèi)在關(guān)系(由航段中轉(zhuǎn)之間天生的拓撲關(guān)系決定),t和y之間以及d和x之間有間接關(guān)系,但是這些關(guān)系的正確性不需要用額外約束條件加以強制,而是經(jīng)由下文中分配算法的特點自動地加以保證。
給定一個存在多種機型飛行的機場航線有向網(wǎng)絡G(N,A)[4],有

可變換為

同時,航線、航段、中轉(zhuǎn)方向的投送量和阻抗及其關(guān)系需滿足式(4)—(7)。
盡管D-MAMUT 模型結(jié)合了多模式以及延誤的因素,但在本質(zhì)上并沒有改變該模型作為非對稱UE 模型的特點。鑒于精簡對角化算法是求解非對稱UE 模型的有效算法,因此把該算法作為D-MAMUT 模型求解算法的基本框架,在此基礎(chǔ)上得到了考慮多模式和延誤的確定型精簡對角化算法,稱其為D-SAMT 算法。
假設(shè)算法在第n次迭代后,得到了航段流量x(n)和中轉(zhuǎn)流量y(n),這時若機型m在航段Ka上的阻抗函數(shù)除了可變外,其余變量均取值于x(n)中的對應項并保持不變,由此得到第n次迭代后可分離的阻抗函數(shù)形式,即

與此類似,可以得到第n次迭代后可分離的中轉(zhuǎn)延誤函數(shù)形式,即

Step 0:初始化。對?m,基于執(zhí)行1 次AON 網(wǎng)絡加載,從而得到初始航段流量x(1)和中轉(zhuǎn)流量y(1)。
Step 1:更新。令n=1,對?m,利用計算及,得到各個
Step 2:確定迭代方向。對?m,基于及執(zhí)行1 次AON 網(wǎng)絡加載,得到各個附加的航段流量和中轉(zhuǎn)流量
Step 3:確定迭代步長。迭代步長λn由以下一維極值問題決定:

Step 4:確定新迭代點。令x(n+1)=xn+ λn×
假設(shè)航線網(wǎng)絡A 如圖1 所示,由1 號機場向9號機場投送兵力,航線網(wǎng)絡中共有12 條航路,飛行機型為A319 和A340,其兵力投送任務量分別為740 人和2 950 人,各條航路上的阻抗參數(shù)和延誤參數(shù)見表2,各航路的中轉(zhuǎn)參數(shù)見表3。

圖1 航線網(wǎng)絡A 示意圖
航段阻抗函數(shù)的計算如下:

中轉(zhuǎn)延誤參數(shù)的計算如下:

其中


表2 阻抗參數(shù)和延誤參數(shù)

表3 中轉(zhuǎn)參數(shù)
求解中,迭代收斂條件是相鄰2 次解的最大相對偏差σn小于某個預定精度值ε,即

式中ε=0.01。
當ε=0.01 時,算法在迭代44 次后收斂到平衡解(x*,y*)。此時將各個和ycd=標到對應的航段和中轉(zhuǎn)機場附近,分別如圖2、圖3 所示。
平衡狀態(tài)下的航路阻抗和流量見表4。機型A319 在這6 條航路上均有流量存在,并且這些航路阻抗在允許誤差范圍內(nèi)可視為相等(約為116.2);機型A340 在6 條航路上也都有流量存在,航路阻抗也都相等(約為86. 8)。這個結(jié)果吻合D-MAMUT 模型的基本出發(fā)點,即多模式UE 條件。

圖2 機型A319 在航線網(wǎng)絡A 中的兵力分配情況

圖3 機型A340 在航線網(wǎng)絡A 中的兵力分配情況

表4 阻抗與流量
本文研究了帶延誤的多機型航空戰(zhàn)略投送兵力分配模型,通過算例對帶中轉(zhuǎn)延誤的2 種機型的航空網(wǎng)絡進行了兵力分配。但航空戰(zhàn)略投送兵力分配模型非常復雜,除本文考慮的影響因素外,機場的保障能力、運力載客量單位和航線網(wǎng)絡中的民航計劃都將對整個投送任務帶來影響,本文將在后續(xù)繼續(xù)研究。
[1] 陳森發(fā),周晶. 城市交通均衡流分配的多模式模型及算法[J].自動化學報,1993,19(2):190-196.
[2] 任剛,王瑋.交通規(guī)劃中的動態(tài)路網(wǎng)及其模型研究[J].公路交通科技,2002,19(6):108-111.
[3] 楊佩昆,錢林波.交通分配中路段行程實踐研究[J].同濟大學學報,1994,22(1):27-32.
[4] 四兵鋒,高自友.多模式的城市混合交通均衡配流模型及算法[J].公路交通科技,1999,16(1):44-48.