鞏 軍,胡 濤,趙 俊,安德紅

(1.海軍工程大學(xué) 管理工程與裝備經(jīng)濟(jì)系， 武漢 430033；2.中國人民解放軍92108部隊(duì)， 浙江 舟山 316000)

1 引言

現(xiàn)代艦船具有大型化、復(fù)雜化、精密化的特點(diǎn)，系統(tǒng)間關(guān)聯(lián)程度更加緊密，設(shè)備技術(shù)狀態(tài)更加復(fù)雜。艦船建造過程中，隨著項(xiàng)目的推進(jìn)，各類風(fēng)險(xiǎn)在不同任務(wù)、流程與工序間多向傳導(dǎo)與耦合，極易擴(kuò)散放大并產(chǎn)生連鎖反應(yīng)。然而，質(zhì)量監(jiān)督部門在風(fēng)險(xiǎn)分析時(shí)缺乏科學(xué)有效的風(fēng)險(xiǎn)傳導(dǎo)分析與評價(jià)方法，風(fēng)險(xiǎn)定量評價(jià)和定級方法粒度較粗，難以對風(fēng)險(xiǎn)進(jìn)行精確計(jì)算與比較，難以有效辨識(shí)和挖掘各類風(fēng)險(xiǎn)的關(guān)聯(lián)關(guān)系與傳導(dǎo)機(jī)理，難以在上游或風(fēng)險(xiǎn)事件發(fā)生結(jié)點(diǎn)及時(shí)發(fā)現(xiàn)、有效釋放與化解各類風(fēng)險(xiǎn)。

風(fēng)險(xiǎn)傳導(dǎo)的概念最早出現(xiàn)在金融領(lǐng)域，并得到了廣泛的推廣與應(yīng)用。作為一條重要分支，國內(nèi)外諸多專家學(xué)者對項(xiàng)目風(fēng)險(xiǎn)傳導(dǎo)理論及其應(yīng)用進(jìn)行了系統(tǒng)的研究。Zegordi等[1]針對供應(yīng)鏈風(fēng)險(xiǎn)給出了基于Petri 網(wǎng)的中斷風(fēng)險(xiǎn)傳導(dǎo)評價(jià)方法；Ghadge[2]基于系統(tǒng)思維對項(xiàng)目風(fēng)險(xiǎn)傳導(dǎo)路徑與過程進(jìn)行了仿真研究；楊敏等[3]綜合運(yùn)用因果貝葉斯網(wǎng)絡(luò)、云模型、影響圖理論、蒙特卡洛等方法給出了一整套項(xiàng)目風(fēng)險(xiǎn)傳導(dǎo)建模、分析與仿真方法。在裝備風(fēng)險(xiǎn)傳導(dǎo)評估領(lǐng)域，孫贇和王瑛等[4,5]基于圖形評審技術(shù)和機(jī)會(huì)理論構(gòu)建了不確定性隨機(jī)多傳導(dǎo)參量網(wǎng)絡(luò)模型和評審技術(shù)模型，具有同時(shí)處理復(fù)雜裝備風(fēng)險(xiǎn)因素為隨機(jī)變量和不確定性變量的優(yōu)勢；李超等在此基礎(chǔ)上[6]運(yùn)用系統(tǒng)耦合理論、極大熵方法對GERT模型的參數(shù)和變量進(jìn)行了優(yōu)化；白焱等[7]提出了基于模糊云模型的風(fēng)險(xiǎn)傳導(dǎo)評估步驟，將風(fēng)險(xiǎn)因素云模型綜合計(jì)算后的風(fēng)險(xiǎn)量與發(fā)生器產(chǎn)生的二維正態(tài)云標(biāo)尺進(jìn)行比較排序并最終完成裝備研制風(fēng)險(xiǎn)的精細(xì)化評價(jià)；劉艷瓊[8]以武器研制項(xiàng)目為背景提出了全面風(fēng)險(xiǎn)圖景的改進(jìn)影響圖構(gòu)建方法、分析框架及其對應(yīng)評估模型與算法，給出了一套武器準(zhǔn)備研制階段性能風(fēng)險(xiǎn)、進(jìn)度風(fēng)險(xiǎn)和費(fèi)用風(fēng)險(xiǎn)的識(shí)別與評估方法。

上述關(guān)于風(fēng)險(xiǎn)傳導(dǎo)相關(guān)研究的不足之處主要體現(xiàn)在：① 缺乏一套適用于大型復(fù)雜武器裝備生產(chǎn)過程或多階段任務(wù)風(fēng)險(xiǎn)識(shí)別、傳導(dǎo)、評價(jià)與定級的理論與方法；② 專家基于經(jīng)驗(yàn)和歷史數(shù)據(jù)給出的風(fēng)險(xiǎn)評價(jià)信息融合過程難以有效解決用精確隸屬度表示模糊事物的不徹底問題，導(dǎo)致風(fēng)險(xiǎn)評價(jià)精度和準(zhǔn)確度受到較大影響；③ 傳統(tǒng)風(fēng)險(xiǎn)定級法硬性分級粒度較粗，忽視了評價(jià)邊界的模糊性和不確定性，對于集結(jié)后風(fēng)險(xiǎn)信息難以給出定性判斷。

鑒于此，本文提出一種基于直覺云和影響圖的艦船建造風(fēng)險(xiǎn)評價(jià)方法。創(chuàng)新之處主要有：① 提出了基于WBS-RBS的風(fēng)險(xiǎn)識(shí)別方法，給出了雙語言變量決策信息轉(zhuǎn)化為直覺正態(tài)云模型定量信息的步驟；② 通過修正風(fēng)險(xiǎn)結(jié)點(diǎn)和價(jià)值結(jié)點(diǎn)，提出了影響圖結(jié)點(diǎn)間耦合關(guān)系的計(jì)算規(guī)則，從拓?fù)鋵?、函?shù)層和數(shù)值層等3個(gè)層面給出了全景式風(fēng)險(xiǎn)傳導(dǎo)影響圖的構(gòu)建及其求解方法；③ 基于蒙特拉洛模擬思想，通過將直覺正態(tài)云風(fēng)險(xiǎn)影響值信息生成大量的云滴樣本的方式將云滴樣本均值作為最優(yōu)估計(jì)值，從而完成風(fēng)險(xiǎn)影響值排序，同時(shí)，基于貼近度的思想給出了一種快速、簡便、準(zhǔn)確的風(fēng)險(xiǎn)定級方法。

2 預(yù)備知識(shí)

2.1 艦船建造風(fēng)險(xiǎn)的涵義

艦船建造風(fēng)險(xiǎn)是指在艦船建造過程中由于技術(shù)、管理、資源、環(huán)境等原因產(chǎn)生各類風(fēng)險(xiǎn)事件并可能造成艦船建造質(zhì)量不過關(guān)、進(jìn)度拖期、費(fèi)用增加等損失，從而對項(xiàng)目目標(biāo)產(chǎn)生影響。用函數(shù)形式可以表示為：

R=F(P,C)

(1)

式(1)中：R為風(fēng)險(xiǎn)影響值；P為風(fēng)險(xiǎn)的發(fā)生概率；C為風(fēng)險(xiǎn)發(fā)生產(chǎn)生某種后果后可能造成的風(fēng)險(xiǎn)損失，風(fēng)險(xiǎn)損失包括質(zhì)量損失CQ、進(jìn)度損失CS和費(fèi)用損失CP。

2.2 直覺正態(tài)云的運(yùn)算法則

李德毅院士于1995年首次提出云模型。在云模型中，正態(tài)云最具有普適性，其隨機(jī)隸屬度在對應(yīng)論域上的分布具有正態(tài)分布的特征。但正態(tài)云模型難以處理某些概念的非隸屬度和猶豫特性，楊惡惡等[9]在正態(tài)云的基礎(chǔ)上提出了直覺正態(tài)云模型，實(shí)現(xiàn)了定性概念模糊性、隨機(jī)性和猶豫度的高效融合。

定義1[9]設(shè)U={x}是一個(gè)用確定數(shù)值表示的定量論域，C是論域U上的定性概念，在U上與C對應(yīng)的直覺正態(tài)云為Y=(〈Ex,ρ,ν〉,En,He)，3個(gè)數(shù)字特征Ex、En和He分別為正態(tài)云模型的期望、熵和超熵，用以定量化反映定性概念的整體特性，ρ和ν為隸屬限和非隸屬限，分別表示當(dāng)x=Ex時(shí)對應(yīng)隸屬度可能取值的上限和下限。

給定正實(shí)數(shù)λ與論域U中的2朵直覺正態(tài)云Y1=(〈Ex1,ρ1,ν1〉,En1,He1)和Y2=(〈Ex2,ρ2,ν2〉,En2,He2)，根據(jù)云模型和直覺模糊數(shù)的計(jì)算規(guī)則，并依據(jù)文獻(xiàn)[9-11]，提出直覺正態(tài)云的集結(jié)運(yùn)算法則，如表1所示。

表1 直覺正態(tài)云的集結(jié)運(yùn)算法則

2.3 風(fēng)險(xiǎn)事件的影響圖表示方法

影響圖(influence diagram，ID)是Howard 和Matheson于1984年在決策樹理論基礎(chǔ)上提出的一種用以表達(dá)不確定信息和求解復(fù)雜決策問題的圖模型。ID由結(jié)點(diǎn)集和有向弧構(gòu)成，結(jié)點(diǎn)集包括決策結(jié)點(diǎn)、機(jī)會(huì)結(jié)點(diǎn)和價(jià)值結(jié)點(diǎn)，能夠?qū)⒋鎯?chǔ)的風(fēng)險(xiǎn)信息通過有向弧進(jìn)行傳導(dǎo)與集結(jié)。因此，ID可用于風(fēng)險(xiǎn)耦合類計(jì)算和風(fēng)險(xiǎn)影響值綜合分析。

2.3.1風(fēng)險(xiǎn)結(jié)點(diǎn)及其表述方法

風(fēng)險(xiǎn)結(jié)點(diǎn)能夠存儲(chǔ)風(fēng)險(xiǎn)事件的概率信息、風(fēng)險(xiǎn)損失信息和綜合風(fēng)險(xiǎn)影響值信息，用橢圓形表達(dá)。根據(jù)風(fēng)險(xiǎn)損失的類型，可分為質(zhì)量、進(jìn)度和費(fèi)用等3類風(fēng)險(xiǎn)結(jié)點(diǎn)。

在風(fēng)險(xiǎn)傳導(dǎo)影響圖模型構(gòu)建時(shí)，假設(shè)風(fēng)險(xiǎn)結(jié)點(diǎn)狀態(tài)與信息僅受到父風(fēng)險(xiǎn)結(jié)點(diǎn)的影響。當(dāng)出現(xiàn)多個(gè)父風(fēng)險(xiǎn)結(jié)點(diǎn)時(shí)，存在父結(jié)點(diǎn)相互耦合后將風(fēng)險(xiǎn)向子節(jié)點(diǎn)傳遞的現(xiàn)象。這種耦合關(guān)系通常包括強(qiáng)耦合、弱耦合和純耦合。一般情況下，進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)之間屬于純耦合關(guān)系，總的風(fēng)險(xiǎn)影響值為其父風(fēng)險(xiǎn)結(jié)點(diǎn)影響值的加權(quán)求和。對于質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)，需要通過專家經(jīng)驗(yàn)對不同風(fēng)險(xiǎn)結(jié)點(diǎn)之間的耦合關(guān)系進(jìn)行評價(jià)后再將風(fēng)險(xiǎn)影響值進(jìn)行綜合集成。

圖1 風(fēng)險(xiǎn)結(jié)點(diǎn)強(qiáng)耦合關(guān)系示意圖

2.3.2價(jià)值結(jié)點(diǎn)及其表述方法

在風(fēng)險(xiǎn)影響值綜合集成過程中，針對FQ、FS和FP對應(yīng)引入3類價(jià)值結(jié)點(diǎn)——質(zhì)量價(jià)值結(jié)點(diǎn)VQ、進(jìn)度價(jià)值結(jié)點(diǎn)VS和費(fèi)用價(jià)值結(jié)點(diǎn)VP，使其實(shí)現(xiàn)風(fēng)險(xiǎn)影響值數(shù)據(jù)對接，再將其進(jìn)一步合成綜合價(jià)值結(jié)點(diǎn)V，如圖2所示。

圖2 風(fēng)險(xiǎn)結(jié)點(diǎn)合成示意圖

(2)

(3)

(4)

(5)

(6)

(7)

3 基于WBS-RBS的艦船建造風(fēng)險(xiǎn)識(shí)別方法

艦船建造過程中存在大量的風(fēng)險(xiǎn)傳導(dǎo)、擴(kuò)散與反饋現(xiàn)象，風(fēng)險(xiǎn)識(shí)別是風(fēng)險(xiǎn)傳導(dǎo)分析的前提，其核心工作：一是如何確定最佳風(fēng)險(xiǎn)辨識(shí)單元；二是如何構(gòu)建風(fēng)險(xiǎn)識(shí)別矩陣，從而為清晰刻畫最佳風(fēng)險(xiǎn)辨識(shí)單元之間的傳導(dǎo)規(guī)律提供前提。

風(fēng)險(xiǎn)存在于各層級任務(wù)與工序中，并隨任務(wù)推進(jìn)和工序流程傳導(dǎo)，任務(wù)與工序是風(fēng)險(xiǎn)存在并進(jìn)行傳導(dǎo)的基本載體。最佳風(fēng)險(xiǎn)辨識(shí)單元是風(fēng)險(xiǎn)傳導(dǎo)過程中的最小任務(wù)單元，為提高風(fēng)險(xiǎn)評價(jià)的可操行和精準(zhǔn)度，需要首先確定風(fēng)險(xiǎn)分析的最小任務(wù)單元層級。WBS-RBS法[13]是一種將WBS和RBS相結(jié)合的風(fēng)險(xiǎn)識(shí)別方法，通過WBS和RBS將項(xiàng)目和風(fēng)險(xiǎn)逐層依次分解，最終將分解之后的任務(wù)單元與風(fēng)險(xiǎn)單元進(jìn)行對應(yīng)，從而能夠全面、系統(tǒng)地識(shí)別風(fēng)險(xiǎn)。艦船建造過程通常采用WBS將整個(gè)過程由高至低逐層分解細(xì)化并得到最小任務(wù)單元或任務(wù)清單[14]，通常分解為艦船、建造過程的五大節(jié)點(diǎn)、五大節(jié)點(diǎn)中的關(guān)鍵工序、關(guān)鍵工序中的關(guān)鍵任務(wù)、可執(zhí)行的任務(wù)包等5個(gè)層級。在風(fēng)險(xiǎn)分析時(shí)，如果將最小任務(wù)單元定位于第5層——可執(zhí)行的任務(wù)包，不僅會(huì)使風(fēng)險(xiǎn)分解得過于零散，增加風(fēng)險(xiǎn)之間傳導(dǎo)與耦合關(guān)系定性分析的難度，還將進(jìn)一步使風(fēng)險(xiǎn)評價(jià)模型更加復(fù)雜；如果將最小任務(wù)單元定位于第3層——五大節(jié)點(diǎn)中的關(guān)鍵工序，會(huì)造成風(fēng)險(xiǎn)分析的粒度不夠精細(xì)，難以準(zhǔn)確分析風(fēng)險(xiǎn)之間的傳導(dǎo)與相互影響關(guān)系。綜合分析后，將最小任務(wù)單元定位于第4層——關(guān)鍵工序中的關(guān)鍵任務(wù)，即最佳風(fēng)險(xiǎn)辨識(shí)單元為第4層中的關(guān)鍵任務(wù)。

確定最佳風(fēng)險(xiǎn)辨識(shí)單元后，需要構(gòu)建風(fēng)險(xiǎn)識(shí)別矩陣，以刻畫風(fēng)險(xiǎn)辨識(shí)單元之間的傳導(dǎo)關(guān)系。為便于風(fēng)險(xiǎn)分析并建立局部與整體之間的連接，這里提出建造模塊工序的概念。建造模塊工序是指艦船建造WBS分解結(jié)構(gòu)第3層中相互獨(dú)立、互不包含的關(guān)鍵工序及其所屬的第4層各項(xiàng)關(guān)鍵任務(wù)。建造模塊工序最大的特點(diǎn)是獨(dú)立性和時(shí)序性，一個(gè)模塊工序具有相對獨(dú)立的檢驗(yàn)驗(yàn)收節(jié)點(diǎn)和相對確定的資源需求，同時(shí)也是各項(xiàng)目共同存在的通用關(guān)鍵工序，各模塊工序間存在時(shí)序或流程關(guān)系。在艦船建造WBS分解的第4層中，各項(xiàng)關(guān)鍵任務(wù)在項(xiàng)目內(nèi)相互影響與傳導(dǎo)，這種影響與傳導(dǎo)是典型的因果關(guān)系，具有明顯的時(shí)序特征。因此，以建造模塊工序?yàn)閱卧?，基于RBS方法逐一識(shí)別模塊工序內(nèi)存在的各層次風(fēng)險(xiǎn)，依據(jù)各模塊工序之間的流程關(guān)系形成風(fēng)險(xiǎn)識(shí)別矩陣。風(fēng)險(xiǎn)評價(jià)時(shí)，通過分析模塊工序中所包含關(guān)鍵任務(wù)之間的相互關(guān)系，可以確定項(xiàng)目內(nèi)風(fēng)險(xiǎn)傳導(dǎo)規(guī)律和路徑。因此，可以通過進(jìn)度管理優(yōu)化模型得到縱向范圍內(nèi)各模塊工序及其所屬風(fēng)險(xiǎn)辨識(shí)單元之間的時(shí)序關(guān)系，從而形成風(fēng)險(xiǎn)識(shí)別矩陣的列向量。

根據(jù)上述思路，艦船建造風(fēng)險(xiǎn)識(shí)別的步驟如下。

步驟1根據(jù)艦船建造風(fēng)險(xiǎn)識(shí)別的目標(biāo)要求，確定風(fēng)險(xiǎn)識(shí)別的范圍與對象，確定模塊工序；

步驟2以項(xiàng)目的建造模塊工序?yàn)榉治鰡卧?，將模塊內(nèi)所有的最小任務(wù)單元可能存在的風(fēng)險(xiǎn)逐步分解細(xì)化、歸類整理，形成風(fēng)險(xiǎn)分解結(jié)構(gòu)；

圖3 WBS-RBS風(fēng)險(xiǎn)識(shí)別矩陣元素框圖

4 艦船建造風(fēng)險(xiǎn)傳導(dǎo)評價(jià)方法

構(gòu)建風(fēng)險(xiǎn)識(shí)別矩陣后，需要厘清風(fēng)險(xiǎn)傳導(dǎo)關(guān)系與規(guī)律。各模塊工序間存在時(shí)序或流程關(guān)系，具有相對獨(dú)立的檢驗(yàn)驗(yàn)收節(jié)點(diǎn)。因此，風(fēng)險(xiǎn)評價(jià)應(yīng)以模塊工序?yàn)檠芯繂卧?，基于風(fēng)險(xiǎn)識(shí)別矩陣分析單元內(nèi)和單元之間風(fēng)險(xiǎn)的耦合與傳導(dǎo)關(guān)系。

4.1 問題描述

現(xiàn)要根據(jù)上述條件求得：① 模塊工序β內(nèi)各風(fēng)險(xiǎn)事件發(fā)生后質(zhì)量風(fēng)險(xiǎn)、進(jìn)度風(fēng)險(xiǎn)與費(fèi)用風(fēng)險(xiǎn)的影響值；② 整個(gè)建造項(xiàng)目中質(zhì)量風(fēng)險(xiǎn)、進(jìn)度風(fēng)險(xiǎn)與費(fèi)用風(fēng)險(xiǎn)的影響值；③ 總風(fēng)險(xiǎn)影響值。

4.2 風(fēng)險(xiǎn)評價(jià)語言信息轉(zhuǎn)化方法

4.2.1基于黃金分割法完成s→(Ex,En,He)

黃金分割法可將專家的定性語言精準(zhǔn)便捷地轉(zhuǎn)化為云模型的定量信息，其基本原理為[15]：

將區(qū)間[Xmin,Xmax]均分為2個(gè)部分，n′為S中的語言標(biāo)度數(shù)量，在區(qū)間上生成與語言標(biāo)度逐一對應(yīng)的n′朵云，設(shè)Y0(Ex0,En0,He0)為處于中間位置的云，則其左右相鄰的云分別為Y-1(Ex-1,En-1,He-1)，Y+1(Ex+1,En+1,He+1)；Y-2(Ex-2,En-2,He-2)，Y+2(Ex+2,En+2,He+2)，…；Y-(n-1)/2(Ex-(n-1)/2,En-(n-1)/2,He-(n-1)/2)，Y+(n-1)/2(Ex+(n-1)/2,En+(n-1)/2,He+(n-1)/2)。由于黃金分割法的思想源于對線段的分割，設(shè)線段的兩端分別為前一朵云與最后一朵云的Ex，后一朵云的Ex為靠近中間云線段的0.382倍對應(yīng)值，前后兩朵云En與He之比為0.618，據(jù)此可生成云的數(shù)字特征。

3.2.2完成h→(ρ,ν)

在H={h1,h2,…h(huán)k′,…,hn′}中，n′為H中的語言標(biāo)度數(shù)量，現(xiàn)要將語言隸屬度變量hk′轉(zhuǎn)化為一個(gè)隸屬度區(qū)間(ρk′,νk′)。由于H中的語言標(biāo)度是均勻分布的，其對應(yīng)隸屬度值域可分為n′個(gè)子區(qū)間，故有：[16]

基于上述2個(gè)環(huán)節(jié)，即完成風(fēng)險(xiǎn)評價(jià)基礎(chǔ)數(shù)據(jù)采集與轉(zhuǎn)化。

4.3 風(fēng)險(xiǎn)傳導(dǎo)評價(jià)方法

為了能夠全面掌握風(fēng)險(xiǎn)事件之間的相互關(guān)系及其在項(xiàng)目中的傳導(dǎo)規(guī)律，需要構(gòu)建全景式風(fēng)險(xiǎn)傳導(dǎo)影響圖，從拓?fù)鋵印⒑瘮?shù)層和數(shù)值層細(xì)致研究風(fēng)險(xiǎn)事件的具體特征并進(jìn)行定量分析。

4.3.1全景式影響圖拓?fù)鋵訕?gòu)建

全景式影響圖拓?fù)鋵拥臉?gòu)建分為以下3個(gè)步驟。

2) 以模塊工序?yàn)閱卧獦?gòu)建風(fēng)險(xiǎn)傳導(dǎo)局部影響圖。針對模塊工序β內(nèi)所有風(fēng)險(xiǎn)事件對應(yīng)的3類風(fēng)險(xiǎn)結(jié)點(diǎn)，依據(jù)圖1構(gòu)建q個(gè)模塊工序的局部影響圖。

3) 構(gòu)建全景式影響圖。根據(jù)逐層整合的順序，依據(jù)模塊工序的時(shí)序整合局部影響圖，得到FQ、FS和FP，并對應(yīng)至VQ、VS和VP，最后綜合價(jià)值結(jié)點(diǎn)V。

4.3.2全景式影響圖函數(shù)層和數(shù)值層確定

2) 確定風(fēng)險(xiǎn)結(jié)點(diǎn)間耦合關(guān)系。進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)之間屬于純耦合關(guān)系，其耦合系數(shù)為0，對于質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)，需要通過專家經(jīng)驗(yàn)對不同風(fēng)險(xiǎn)結(jié)點(diǎn)之間的耦合關(guān)系進(jìn)行評價(jià)后再將風(fēng)險(xiǎn)影響值進(jìn)行綜合集成。

4.3.3風(fēng)險(xiǎn)傳導(dǎo)分析步驟

在風(fēng)險(xiǎn)傳導(dǎo)過程中，每一風(fēng)險(xiǎn)事件發(fā)生后均可能造成質(zhì)量、進(jìn)度和費(fèi)用損失與影響，本文假設(shè)風(fēng)險(xiǎn)事件發(fā)生后產(chǎn)生的質(zhì)量損失、進(jìn)度損失和費(fèi)用損失彼此獨(dú)立。風(fēng)險(xiǎn)傳導(dǎo)評價(jià)時(shí)分別計(jì)算每一風(fēng)險(xiǎn)事件產(chǎn)生的質(zhì)量、進(jìn)度和費(fèi)用損失影響值后再進(jìn)行綜合。

1) 計(jì)算不同模塊工序中各風(fēng)險(xiǎn)事件對應(yīng)風(fēng)險(xiǎn)結(jié)點(diǎn)的風(fēng)險(xiǎn)影響值，其計(jì)算公式如下：

(8)

2) 計(jì)算各進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)的風(fēng)險(xiǎn)影響值及其貢獻(xiàn)度。

步驟1：由式(9)計(jì)算模塊工序β(β=1,2,…,q)進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)的風(fēng)險(xiǎn)影響值。

(9)

由于直覺正態(tài)云模型沒有“÷”運(yùn)算，運(yùn)用式(6)和式(7)進(jìn)行直覺正態(tài)云信息的“÷”運(yùn)算前需要對數(shù)據(jù)進(jìn)行轉(zhuǎn)化處理。具體過程參照4.1節(jié)“風(fēng)險(xiǎn)影響值比較方法”，將式中分子、分母的直覺正態(tài)云數(shù)據(jù)借助計(jì)算機(jī)程序，通過云滴生成算法獲取足夠的云滴樣本，將樣本均值作為分子、分母的最優(yōu)估計(jì)值后再進(jìn)行“÷”運(yùn)算。

步驟3：由式(10)、式(11)計(jì)算整個(gè)項(xiàng)目的進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)影響值及其貢獻(xiàn)度。

(10)

(11)

3) 計(jì)算各質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)的風(fēng)險(xiǎn)影響值及其貢獻(xiàn)度。

步驟1：獲取模塊工序1中各二級風(fēng)險(xiǎn)事件對應(yīng)的風(fēng)險(xiǎn)結(jié)點(diǎn)耦合關(guān)系數(shù)據(jù)。

(12)

(13)

需要注意的是，由于緊前模塊工序β-1(β=2,3,…,q)的質(zhì)量風(fēng)險(xiǎn)將整體傳導(dǎo)至緊后模塊工序β，并與其質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)發(fā)生耦合關(guān)系。從模塊工序2開始計(jì)算風(fēng)險(xiǎn)影響值時(shí)需要專家對緊前模塊工序β-1的風(fēng)險(xiǎn)結(jié)點(diǎn)與緊后模塊工序中各風(fēng)險(xiǎn)結(jié)點(diǎn)的耦合關(guān)系進(jìn)行分析。

4) 計(jì)算總風(fēng)險(xiǎn)影響值。

將FQ、FS和FP對應(yīng)到價(jià)值結(jié)點(diǎn)VQ、VS和VP，使其實(shí)現(xiàn)風(fēng)險(xiǎn)影響值數(shù)據(jù)對接，并將其進(jìn)一步合成綜合價(jià)值結(jié)點(diǎn)V。由式(14)可求得總風(fēng)險(xiǎn)影響值，并將其存儲(chǔ)于V中。

R=ψQRQ⊕ψSRS⊕ψPRP

(14)

5) 對各級風(fēng)險(xiǎn)進(jìn)行比較定級，并提出相應(yīng)對策。

5 基于蒙特卡洛模擬的風(fēng)險(xiǎn)比較與定級方法

5.1 風(fēng)險(xiǎn)影響值比較方法

(14)

根據(jù)定義2和式(14)，可以比較2朵直覺正態(tài)云所代表風(fēng)險(xiǎn)事件的影響值大小。

5.2 風(fēng)險(xiǎn)定級方法

假設(shè)Rd=(〈Exd,ρd,νd〉,End,Hed)是待定級風(fēng)險(xiǎn)事件的風(fēng)險(xiǎn)影響值，將風(fēng)險(xiǎn)由低到高分為t′個(gè)等級，分別為f1、f2、…、ft′，則風(fēng)險(xiǎn)定級步驟為：

步驟1：根據(jù)云的生成方法，將定性等級語言f1、f2、…、ft′轉(zhuǎn)化為正態(tài)云模型的定量信息，fl→Yl，其中Yl=(Exl,Enl,Hel)，l=1,2,…,t′。

步驟2：將正態(tài)云信息Yl看作專家對評價(jià)結(jié)果信心水平最高時(shí)的直覺正態(tài)云信息Yl′=(〈Exl,ρk′,νk′〉,Enl,Hel)，根據(jù)h→(ρ,ν)轉(zhuǎn)化方法求解ρk′,νk′。

6 算例分析

船臺(tái)是造船廠最重要的生產(chǎn)設(shè)施，船臺(tái)施工是組織難度最大、投入勞動(dòng)及各項(xiàng)資源最多的階段。由于大量的工序交叉進(jìn)行，各類風(fēng)險(xiǎn)隨著項(xiàng)目的推進(jìn)多向傳導(dǎo)，極易擴(kuò)散放大并產(chǎn)生連鎖反應(yīng)。因此，針對即將進(jìn)入船臺(tái)施工階段的某型艦船，軍事代表室亟需對存在的質(zhì)量、進(jìn)度和費(fèi)用風(fēng)險(xiǎn)進(jìn)行分析和評估。

6.1 船臺(tái)施工階段風(fēng)險(xiǎn)識(shí)別

1) 根據(jù)最佳風(fēng)險(xiǎn)辨識(shí)單元的識(shí)別原則與方法，將船臺(tái)施工階段劃分為總段形成模塊、全船形成模塊、船臺(tái)舾裝模塊和軸系施工模塊。

2) 軍事代表室前期在綜合考慮艦船建造的實(shí)際工作和專家訪談的基礎(chǔ)上，根據(jù)風(fēng)險(xiǎn)引發(fā)的原因進(jìn)行分類，將風(fēng)險(xiǎn)事件一級指標(biāo)分為技術(shù)風(fēng)險(xiǎn)(A1)、管理風(fēng)險(xiǎn)(A2)、資源風(fēng)險(xiǎn)(A3)和環(huán)境風(fēng)險(xiǎn)(A4)，對應(yīng)二級指標(biāo)分別為：技術(shù)狀態(tài)變更風(fēng)險(xiǎn)(A11)、新技術(shù)應(yīng)用風(fēng)險(xiǎn)(A12)、技術(shù)指標(biāo)風(fēng)險(xiǎn)(A13)、質(zhì)量問題處理風(fēng)險(xiǎn)(A21)、檢驗(yàn)驗(yàn)收風(fēng)險(xiǎn)(A22)、現(xiàn)場管理風(fēng)險(xiǎn)(A23)、人員配置風(fēng)險(xiǎn)(A31)、設(shè)備材料風(fēng)險(xiǎn)(A32)、現(xiàn)場條件風(fēng)險(xiǎn)(A41)。

3) 以模塊工序?yàn)榉治鰡卧瑢⒖赡艽嬖诘娘L(fēng)險(xiǎn)逐步歸類整理，形成風(fēng)險(xiǎn)分解結(jié)構(gòu)；同時(shí)，將WBS分解后的模塊工序和RBS的下層風(fēng)險(xiǎn)分別作為行和列構(gòu)建風(fēng)險(xiǎn)識(shí)別矩陣，有關(guān)元素如表2所示。

表2 WBS-RBS風(fēng)險(xiǎn)識(shí)別矩陣有關(guān)元素

6.2 船臺(tái)施工階段風(fēng)險(xiǎn)數(shù)據(jù)采集與轉(zhuǎn)化

現(xiàn)請3位專家e={e1,e2,e3}對各模塊工序中各風(fēng)險(xiǎn)事件發(fā)生概率及其造成的質(zhì)量、進(jìn)度和費(fèi)用損失等風(fēng)險(xiǎn)數(shù)據(jù)進(jìn)行評價(jià)。專家權(quán)重向量為q={0.3,0.4,0.3}。以總段形成模塊中風(fēng)險(xiǎn)A21發(fā)生概率轉(zhuǎn)化過程為例進(jìn)行說明。3位專家給出風(fēng)險(xiǎn)A21的發(fā)生概率雙語言評價(jià)信息如下：

假設(shè)本文對風(fēng)險(xiǎn)事件概率等級和后果等級均采用5標(biāo)度法進(jìn)行描述，給定區(qū)間為[0,1]，He0=0.01。根據(jù)4.2.1節(jié)內(nèi)容，基于黃金分割法計(jì)算云數(shù)字特征，得到語言評價(jià)標(biāo)度對應(yīng)語義信息的云模型見表3。

表3 語言評價(jià)標(biāo)度與語義信息云模型對應(yīng)表

根據(jù)4.2.2節(jié)h→(ρ,ν)轉(zhuǎn)化方法，語言隸屬度標(biāo)度描述如表4所示。

表4 語言隸屬度標(biāo)度描述

6.3 船臺(tái)施工階段風(fēng)險(xiǎn)傳導(dǎo)評價(jià)

6.3.1全景式影響圖構(gòu)建

圖4 風(fēng)險(xiǎn)傳導(dǎo)全景式影響框圖

在此基礎(chǔ)上，由專家給出不同風(fēng)險(xiǎn)事件所對應(yīng)的質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)間耦合關(guān)系。

6.3.2計(jì)算各進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)的風(fēng)險(xiǎn)影響值

各進(jìn)度和費(fèi)用風(fēng)險(xiǎn)結(jié)點(diǎn)的風(fēng)險(xiǎn)影響值，相關(guān)數(shù)據(jù)計(jì)算如下：

RS=(〈0.491,0.542,0.958〉,0.020,0.006)

RP=(〈0.339,0.439,0.920〉,0.013,0.005)

6.3.3計(jì)算各質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)的相關(guān)數(shù)據(jù)

各質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)的相關(guān)數(shù)據(jù)計(jì)算如下：

總風(fēng)險(xiǎn)影響值為：

R=(〈0.488,0.530,0.936〉,0.012,0.006)

6.3.4風(fēng)險(xiǎn)比較與定級

表5 基于蒙特卡洛模擬的風(fēng)險(xiǎn)等級表

通過分析，可以得到以下結(jié)論：

1) 由表6可知，船臺(tái)施工整個(gè)階段總體風(fēng)險(xiǎn)等級為3級，其中，質(zhì)量風(fēng)險(xiǎn)和進(jìn)度風(fēng)險(xiǎn)等級均為3級，費(fèi)用風(fēng)險(xiǎn)等級為2級，總體風(fēng)險(xiǎn)可控，質(zhì)量風(fēng)險(xiǎn)最高；模塊工序2和工序3產(chǎn)生的質(zhì)量風(fēng)險(xiǎn)較高，等級均為4級，模塊工序2產(chǎn)生的進(jìn)度風(fēng)險(xiǎn)較高，等級為4級。因此，軍事代表室應(yīng)重點(diǎn)關(guān)注模塊工序2、工序3產(chǎn)生的質(zhì)量風(fēng)險(xiǎn)和模塊工序2產(chǎn)生的進(jìn)度風(fēng)險(xiǎn)。

表6 風(fēng)險(xiǎn)結(jié)點(diǎn)影響值相關(guān)數(shù)據(jù)匯總表

船臺(tái)施工階段整個(gè)過程中質(zhì)量風(fēng)險(xiǎn)較高的原因主要在于部分結(jié)點(diǎn)之間存在強(qiáng)耦合關(guān)系，緊前模塊工序β-1(β=2,3,4)的質(zhì)量風(fēng)險(xiǎn)將整體傳導(dǎo)至緊后模塊工序β，并與其質(zhì)量風(fēng)險(xiǎn)結(jié)點(diǎn)發(fā)生耦合關(guān)系，造成模塊1～4的質(zhì)量風(fēng)險(xiǎn)等級較高，各結(jié)點(diǎn)質(zhì)量風(fēng)險(xiǎn)的貢獻(xiàn)度也呈現(xiàn)出不斷增長的趨勢。

表7 質(zhì)量風(fēng)險(xiǎn)影響值

表8 計(jì)分統(tǒng)計(jì)

表9 20次模擬統(tǒng)計(jì)

圖5 資源風(fēng)險(xiǎn)等級云仿真圖

7 結(jié)論

提出了基于WBS-RBS的風(fēng)險(xiǎn)識(shí)別步驟、風(fēng)險(xiǎn)識(shí)別矩陣的構(gòu)建方法；為了厘清風(fēng)險(xiǎn)傳導(dǎo)關(guān)系與規(guī)律，將雙語言變量決策信息轉(zhuǎn)化為直覺正態(tài)云模型定量信息，基于全景式影響圖搭建了風(fēng)險(xiǎn)傳導(dǎo)評價(jià)框架；給出了基于蒙特卡洛模擬的風(fēng)險(xiǎn)比較與定級方法，形成了一套完備的艦船建造過程風(fēng)險(xiǎn)識(shí)別、傳導(dǎo)、評價(jià)與定級的方法；通過算例分析，驗(yàn)證了方法的有效性與實(shí)用性。該方法可為軍事代表室有效開展合同履行監(jiān)督風(fēng)險(xiǎn)分析工作提供指導(dǎo)，并可用于大型復(fù)雜武器裝備生產(chǎn)過程或多階段任務(wù)風(fēng)險(xiǎn)評價(jià)領(lǐng)域相關(guān)工作。