基于復雜網(wǎng)絡(luò)理論的裝備保障網(wǎng)絡(luò)模型研究

張 勇， 楊宏偉， 楊學強， 黃 俊

（裝甲兵工程學院 技術(shù)保障工程系，北京 100072）

裝備保障網(wǎng)絡(luò)是以裝備保障設(shè)施為依托，將各種保障資源按照一定的要求和原則合理部署而最終形成的一個網(wǎng)絡(luò)化布局的保障體系，是各級保障部門完成保障任務的物質(zhì)基礎(chǔ)［1］.新的歷史條件下，隨著綜合保障、一體化保障模式和模塊化編組形式的出現(xiàn)以及信息化技術(shù)的推動，可以看作裝備保障網(wǎng)絡(luò)節(jié)點的裝備保障實體數(shù)量眾多、級別不同、分布各異，如何準確刻畫和描述其組織形式并發(fā)現(xiàn)其中蘊含的特點規(guī)律，是新時期認識網(wǎng)絡(luò)化裝備保障體系復雜性的基礎(chǔ)性工作.

本文從整體角度出發(fā)，運用復雜網(wǎng)絡(luò)相關(guān)理論［2－3］來構(gòu)建裝備保障網(wǎng)絡(luò)的組織結(jié)構(gòu)，將其表示為具有復雜拓撲結(jié)構(gòu)和動力學行為、由大量的節(jié)點通過相互之間的作用關(guān)系連接而成的網(wǎng)絡(luò)圖，為信息化條件下研究網(wǎng)絡(luò)化裝備保障系統(tǒng)的結(jié)構(gòu)和行為提供一種行之有效的方法和途徑.本文主要從兩個方面進行研究：一是有別于傳統(tǒng)軍隊編制體制條件下具有嚴格層次關(guān)系的樹狀網(wǎng)絡(luò)和文獻［4］提出的具有“越級指揮”特征的確定性網(wǎng)絡(luò)模型，構(gòu)建與現(xiàn)實更加符合的裝備保障網(wǎng)絡(luò)演化模型；二是對網(wǎng)絡(luò)結(jié)構(gòu)參數(shù)進行靈敏度分析，對實際裝備保障網(wǎng)絡(luò)的拓撲結(jié)構(gòu)和網(wǎng)絡(luò)特性提供更加精確和細致的描述.

1 裝備保障網(wǎng)絡(luò)模型

1.1 概念模型

按照復雜網(wǎng)絡(luò)和裝備保障體系建設(shè)的相關(guān)觀點［5－7］，構(gòu)建裝備保障網(wǎng)絡(luò)拓撲模型過程中有以下幾點假設(shè)和前提：

a.節(jié)點指在一定的時間、空間范圍內(nèi)，具有發(fā)出或接收裝備保障資源功能的保障實體.本文中定義為某一戰(zhàn)區(qū)范圍內(nèi)所有的維修機構(gòu)、器材保障機構(gòu)等.

b.邊代表節(jié)點（對象）之間的相互作用、相互關(guān)聯(lián)［8］.原則上，裝備保障網(wǎng)絡(luò)各節(jié)點之間可通過公路、鐵路、水路、航線及管道等實體進行支援、協(xié)同及資源交換等保障行動，一般認為2個節(jié)點之間最多可能有1條邊，由于實體之間的聯(lián)系不考慮邊的方向性問題（雙向連接），所有節(jié)點之間的連接為無向邊.

另外，為了區(qū)別于公交網(wǎng)絡(luò)中站點停靠、換乘方式［9］，本文考慮各保障實體之間的業(yè)務關(guān)系，即部隊之間嚴格的層次結(jié)構(gòu)、上下級支援關(guān)系及同級協(xié)同關(guān)系等，設(shè)定以下不同的連接方式：

a.樹狀結(jié)構(gòu).嚴格的層級關(guān)系，即傳統(tǒng)的裝備保障體系基本按照“軍區(qū)－軍－師－團”的力量配置，呈樹狀延伸，橫向之間聯(lián)系少，表現(xiàn)為長的平均路徑和低的聚合度，產(chǎn)生的結(jié)果是信息傳遞層級多，溝通連接復雜，易被分割破壞，整體保障效率低下.好處是結(jié)構(gòu)簡單，易于管理控制，如圖1（a）和1（b）所示.

b.全連通結(jié)構(gòu).參照層次結(jié)構(gòu)特征，近似于完全連通網(wǎng)絡(luò)，上級保障實體對所屬范圍內(nèi)的保障實體實行完全支援連接，同級之間存在協(xié)同關(guān)系，其好處是系統(tǒng)魯棒性強，抗毀性突出，信息傳遞層級少，體系效率較高.缺點是結(jié)構(gòu)復雜，不利于管理控制，如圖1（d）所示.

c.混合結(jié)構(gòu).介于上述兩種結(jié)構(gòu)之間，具體表現(xiàn)為，上下級別保障實體之間根據(jù)隸屬關(guān)系的不同存在不同概率的支援保障關(guān)系，同一級別保障實體根據(jù)地域關(guān)系的不同存在不同的協(xié)同保障概率，不同級別之間存在一定的“越級支援”的保障模式（見圖1（c））.該方式一定程度反映了信息化條件下指揮控制的靈活性，又不失軍隊層次結(jié)構(gòu)特征，能較為真實地反映實際裝備保障網(wǎng)絡(luò)的連接和形成情況.

圖1 裝備保障網(wǎng)絡(luò)示意圖Fig.1 Equipment support network

1.2 生成機制建模

通過上述分析，借鑒樹狀網(wǎng)絡(luò)結(jié)構(gòu)和文獻［4］所論述的確定性模型的思想，并對其適當進行改進，確定裝備保障網(wǎng)絡(luò)生成機制如下：

a.從迭代步數(shù)s＝0開始，即從一個節(jié)點出發(fā)；

b.當s＝1時，在第一個節(jié)點下方的兩邊對稱地分布n個節(jié)點，且使這n個節(jié)點與第一個節(jié)點相連；

c.當s＝2時，按照相同的方法，在新生長出來的n個節(jié)點的下方分別再對稱生長出來n個節(jié)點，以一定概率p實現(xiàn)不同等級（同級或上下級）節(jié)點之間的連接.

依照這種方法不斷生長，直到生成預期的網(wǎng)絡(luò)模型，如圖2所示.

圖2 生成機制演示圖Fig.2 Mechanism of evolving model

p代表網(wǎng)絡(luò)隨機連接概率，文中指裝備保障網(wǎng)絡(luò)中沒有直接邊相連的節(jié)點對（vi，vj）生成連接邊的概率.在此，結(jié)合無標度網(wǎng)絡(luò)的生成機制和裝備保障網(wǎng)絡(luò)的現(xiàn)實背景，從節(jié)點和邊這2個方面給出連接規(guī)則影響因子.

a.節(jié)點層面影響因子pn.節(jié)點層面，包括結(jié)構(gòu)屬性pns和物理屬性pnp，公式為

式中，ki，kj分別表示節(jié)點vi，vj的 度；si，sj分別表示節(jié)點vi，vj的強度（負荷）；［k］，［s］為相對應的網(wǎng)絡(luò)平均值.

b.邊層面影響因子pe.邊層面主要表征該節(jié)點對之間的距離，與節(jié)點類似，包括結(jié)構(gòu)屬性pes和物理屬性pep，公式為

式中，li，lj分別表示節(jié)點對（vi，vj）之 間的結(jié)構(gòu)距離；di，dj分別表示節(jié)點對（vi，vj）之間的現(xiàn)實距離；［l］，［d］為相對應的網(wǎng)絡(luò)平均值.

定義隨機連接概率

式中，α為隨機連接因子.

特別地，當α＝0時，本文所提出的模型就演化為樹狀網(wǎng)絡(luò)結(jié)構(gòu)模型.另外，當α增大到一定程度時，裝備保障網(wǎng)絡(luò)也可以形成全連通網(wǎng)絡(luò).

綜上，通過引入隨機連接概率p可以實現(xiàn)裝備保障網(wǎng)絡(luò)在樹狀層次結(jié)構(gòu)的基礎(chǔ)上根據(jù)地域、隸屬關(guān)系及資源擁有量等屬性進行擇優(yōu)連接，從而構(gòu)造出與現(xiàn)實情況更加符合的裝備保障演化網(wǎng)絡(luò).該模型很好地刻畫了真實網(wǎng)絡(luò)的“擇優(yōu)連接”和“就地就近”的保障原則，即無論從結(jié)構(gòu)層面還是物理屬性層面，裝備保障網(wǎng)絡(luò)更加傾向于網(wǎng)絡(luò)中的“富者”（度、強度大）、距離近（結(jié)構(gòu)距離、實際距離）的節(jié)點對相連接，在式（3）中設(shè)定連接概率p與節(jié)點重要度pn成正比，而與節(jié)點距離pe成反比關(guān)系，為了修正兩者之間的關(guān)系，引入修正系數(shù)α，也稱隨機連接因子，可表征裝備保障網(wǎng)絡(luò)中不同節(jié)點對邊之間的隨機連接程度.

1.3 算法實現(xiàn)

具體算法實現(xiàn)如下：

a.輸入迭代步數(shù)s和子節(jié)點數(shù)目n，生成對應的樹狀結(jié)構(gòu)網(wǎng)絡(luò)Gt；

b.根據(jù)不同的規(guī)則對網(wǎng)絡(luò)節(jié)點和邊相應賦權(quán)，如依據(jù)統(tǒng)計數(shù)據(jù)或隨機賦權(quán)，形成加權(quán)初始裝備保障網(wǎng)絡(luò)Gw；

c.計算Gt中所有沒有邊連接的節(jié)點對（vi，vj）集合，根據(jù)式（1）～（3）計算連接概率p，并依此隨機生成連接邊，直至生成最終加權(quán)裝備保障網(wǎng)絡(luò)G；

文中圖3（見下頁）即是作者根據(jù)典型的部隊編制體制，分戰(zhàn)區(qū)（ZQ）－集團軍（J）－師（S）－團（T）這4級保障實體節(jié)點，上下級別之間采取3∶3編制，即戰(zhàn)區(qū)內(nèi)轄3個集團軍，1個集團軍又轄3個編制師，依次類推.取s＝3，n＝3，α＝0.2時，基于 Matlab 7.0和Netdraw等相關(guān)分析軟件生成的裝備保障網(wǎng)絡(luò)拓撲圖.

2 模型特征參數(shù)分析

2.1 度分布P（k）

度分布表示網(wǎng)絡(luò)的節(jié)點度恰好為k的概率.網(wǎng)絡(luò)中節(jié)點i的度ki定義為與該節(jié)點連接的其它節(jié)點的數(shù)目，也可視為與這一節(jié)點連接的邊的數(shù)量.節(jié)點的度分布反映了網(wǎng)絡(luò)的整體性質(zhì)，根據(jù)復雜網(wǎng)絡(luò)無標度特性［3］：少部分節(jié)點的度遠遠高于絕大部分的節(jié)點的度；度分布概率在雙對數(shù)坐標下顯示出了較明顯的冪律特征.

調(diào)節(jié)文中模型的各個影響參數(shù)，得到裝備保障網(wǎng)絡(luò)模型度分布變化圖，如圖4所示.其中，圖4（a）分析了迭代步數(shù)s不同條件下裝備保障網(wǎng)絡(luò)雙對數(shù)度分布曲線，不難發(fā)現(xiàn)，在子節(jié)點數(shù)目n＝3和隨機連接因子α＝0.02的前提下，裝備保障網(wǎng)絡(luò)P（k）曲線呈現(xiàn)出較為明顯的無標度特征，P（k）～k－γ，γ≈1.5.與之對應，如果改變子節(jié)點數(shù)目，當n比較小的時候（文中為n≤5），裝備保障網(wǎng)絡(luò)尚具有一定的無標度特性，隨著子節(jié)點數(shù)目的進一步增加，無標度特性逐漸減弱，在這種情況下，整體裝備保障網(wǎng)絡(luò)傾向于隨機連接，而不是擇優(yōu)連接.另外，這一特性，在圖4（c）中表現(xiàn)得更加明顯，即隨著隨機連接因子α的增大，裝備保障網(wǎng)絡(luò)無標度特性逐步消失而隨機網(wǎng)絡(luò)特性逐步凸顯，這也可以從圖4（d）得到更為直觀清晰的認識.

值得一提的是，無尺度網(wǎng)絡(luò)的典型特征是面對隨機攻擊時具有很高的承受能力（魯棒性），而對關(guān)鍵節(jié)點進行攻擊則表現(xiàn)出極大的脆弱性，與之對應的全連通網(wǎng)絡(luò)具備較高的抗毀性.

2.2 平均聚集系數(shù)C

聚集系數(shù)是專門用來衡量網(wǎng)絡(luò)節(jié)點集聚程度的一個重要參數(shù)，屬局部特征量［10］.一般地，聚集系數(shù)越大，表示個體間的聯(lián)系愈緊密.節(jié)點i的聚集系數(shù)Ci定義為節(jié)點i的鄰接點之間實際存在的邊數(shù)與所有可能的邊數(shù)的比值.一般地，假設(shè)網(wǎng)絡(luò)中1個節(jié)點i有ki條邊將它與其它節(jié)點連接，這ki個節(jié)點之間最多可能有ki（ki－1）／2條邊.而這ki個節(jié)點之間實際存在的邊數(shù)Ei和總的可能邊數(shù)ki（ki－1）／2之比就定義為節(jié)點的聚類系數(shù)Ci，即

網(wǎng)絡(luò)的平均聚集系數(shù)C指所有節(jié)點的聚集系數(shù)Ci的均值，反映了網(wǎng)絡(luò)整體的內(nèi)聚性.

很明顯，0≤C≤1.C＝0當且僅當所有的節(jié)點為孤立節(jié)點；C＝1當且僅當網(wǎng)絡(luò)是全局耦合的.

圖3 裝備保障網(wǎng)絡(luò)拓撲結(jié)構(gòu)圖Fig.3 Topology of equipment support network

圖4 裝備保障網(wǎng)絡(luò)的度分布圖Fig.4 Degree distribution of equipment support network

圖5 裝備保障網(wǎng)絡(luò)的聚集系數(shù)分布Fig.5 Clustering coefficient of equipment support network

如圖5所示，裝備保障網(wǎng)絡(luò)的聚集系數(shù)隨著迭代步數(shù)、子節(jié)點數(shù)目和隨機連接因子的增大而增大，意即隨著網(wǎng)絡(luò)中元素（節(jié)點和邊）的增加，節(jié)點之間的聯(lián)系更加趨于緊密，這符合典型的小世界特征［2］.特別地，在迭代到第3步時，聚集系數(shù)曲線出現(xiàn)了較為明顯的拐點（見圖5（a）），而當子節(jié)點數(shù)目小于4之前，聚集系數(shù)表現(xiàn)出較強的增長趨勢，隨后趨于平緩（見圖5（b））.另外，根據(jù)本文中所賦予的權(quán)值，當隨機連接因子小于0.2時（見圖5（c）），聚集系數(shù)增長迅速，之后呈現(xiàn)波動式增長.經(jīng)測算，當α≈13，C≈1，即此時形成全連通圖.

2.3 平均路徑長度L

平均路徑長度指所有節(jié)點之間的最短路徑的均值.一般定義節(jié)點i和j之間的距離dij為連接兩者所要經(jīng)歷的邊的最小數(shù)目，網(wǎng)絡(luò)的平均路徑長度是所有節(jié)點對之間距離的平均值.從路長的分布可以看出節(jié)點間的遠近程度，其中任意兩點間的最大距離定義為網(wǎng)絡(luò)的直徑D.

D描述的是網(wǎng)絡(luò)節(jié)點間的分離程度，也就是說網(wǎng)絡(luò)有多小，文中裝備保障網(wǎng)絡(luò)的直徑為3.

網(wǎng)絡(luò)的平均路徑長度L定義為任意2個節(jié)點之間的距離的平均值，即

式中，N為網(wǎng)絡(luò)節(jié)點數(shù).

平均路徑長度可以衡量裝備保障網(wǎng)絡(luò)中各節(jié)點之間的傳輸特性和連通性，一般當L比較小時，信息、物質(zhì)或能量交換的速度較快.如對器材保障而言，降低網(wǎng)絡(luò)的網(wǎng)絡(luò)直徑，即降低器材的最大周轉(zhuǎn)次數(shù)；降低網(wǎng)絡(luò)的平均最短距離，即降低整個網(wǎng)絡(luò)中器材的平均周轉(zhuǎn)次數(shù).通過計算得出，文中圖3所示裝備保障網(wǎng)絡(luò)的平均路徑長度為1.895，即在該網(wǎng)絡(luò)中任意2個保障節(jié)點之間傳輸物資的話，平均經(jīng)過不到2次周轉(zhuǎn)就可到達.圖6即為裝備保障網(wǎng)絡(luò)在迭代步數(shù)、子節(jié)點數(shù)目和隨機連接因子變化時平均路徑長度變化情況.

圖6 裝備保障網(wǎng)絡(luò)平均最短距離Fig.6 Aver path length of equipment support network

通過圖6可以發(fā)現(xiàn)，裝備保障網(wǎng)絡(luò)平均路徑長度隨著影響因子的增加逐步減小，特別地，當?shù)綌?shù)為3～4，子節(jié)點數(shù)目為3時，L均出現(xiàn)下降的趨勢，這與圖5中的仿真結(jié)果基本是一致的，筆者認為出現(xiàn)這種現(xiàn)象主要是節(jié)點數(shù)目偏少，以及文中節(jié)點權(quán)值效應不明顯所導致.不過整體而言，按照相關(guān)定義［11－12］，裝備保障網(wǎng)絡(luò)仍具有一定的小世界特性.

3 結(jié)束語

基于本文所提出裝備保障網(wǎng)絡(luò)演化模型及相關(guān)分析，得到以下幾點結(jié)論和認識：

a.本文所構(gòu)建模型可以更加真實地描述具有層級結(jié)構(gòu)和隨機連接特性的裝備保障網(wǎng)絡(luò)，引入的3個影響因子可以較為精確地控制生成網(wǎng)絡(luò)規(guī)模和節(jié)點連接的緊密程度，可以有效地說明所構(gòu)建網(wǎng)絡(luò)的層次、跨度及隨機連接程度等具有一定現(xiàn)實意義的衡量指標.

b.通過仿真分析可知，該裝備保障網(wǎng)絡(luò)演化模型反映出一定的無標度和小世界特性，可以作為分析裝備保障網(wǎng)絡(luò)或類似網(wǎng)絡(luò)結(jié)構(gòu)和動力學研究的基礎(chǔ).

c.基于復雜網(wǎng)絡(luò)理論研究裝備保障網(wǎng)絡(luò)具有很強的應用背景，如裝備保障網(wǎng)絡(luò)中戰(zhàn)役支撐點的確定、樞紐點的確定、道路基礎(chǔ)設(shè)施的構(gòu)建、保障效益與連通規(guī)模的關(guān)系、樞紐節(jié)點的堵塞所導致的保障網(wǎng)絡(luò)穩(wěn)定性問題等可以而且有待于進一步分析和研究.

值得說明的是，由于現(xiàn)實數(shù)據(jù)來源的有限性和保密性要求，本文在研究過程中對裝備保障網(wǎng)絡(luò)節(jié)點和邊的權(quán)重均作了一定假設(shè)和限制，這有待于在今后的研究過程中進行驗證和修改.

［1］王文峰.裝備保障網(wǎng)絡(luò)優(yōu)化設(shè)計問題研究［D］.長沙：國防科學技術(shù)大學，2008.

［2］Watts D J，Strogatz S H.Collective dynamics of“small－world”networks［J］.Nature，1998，393（4）：440－442.

［3］Barabási A L，Albert R.Emergence of scaling in random networks［J］.Science，1999，286（5439）：509－512.

［4］Chen Mu，Yu Boming，Xu Peng，et al.A new deterministic complex network model with hierarchical structure［J］.Physica A，2007，385：307－317.

［5］吳金閃，狄增如.從統(tǒng)計物理學看復雜網(wǎng)絡(luò)研究［J］.物理學進展，2004，24（1）：18－23.

［6］汪曉帆，李翔，陳關(guān)榮.復雜網(wǎng)絡(luò)理論及其應用［M］.北京：清華大學出版社，2006.

［7］李勇.物流保障網(wǎng)絡(luò)級聯(lián)失效抗毀性研究［D］.長沙：國防科學技術(shù)大學，2009.

［8］張嗣瀛.復雜網(wǎng)絡(luò)的演化過程，n（n－1）律，自聚律［J］.復雜系統(tǒng)與復雜性科學，2005，2（3）：84－90.

［9］顧前，楊旭華，王萬良，等.基于復雜網(wǎng)絡(luò)的城市公共交通網(wǎng)絡(luò)研究［J］.計算機工程，2008，34（20）：266－268.

［10］Almaas E，Kulkarni R V，Stroud D.Characterizing the structure of small－world networks［J］.Phys Rev Let，2002，88（9）：98－101.

［11］王建偉，榮莉莉.突發(fā)事件的連鎖反應網(wǎng)絡(luò)模型研究［J］.計算機應用研究，2008，25（11）：3288－3291.

［12］劉建國.復雜網(wǎng)絡(luò)模型構(gòu)建及其在知識系統(tǒng)中的應用［D］.大連：大連理工大學，2006.