李明超,段志龍

(天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072)

耦合風(fēng)險(xiǎn)因素的大型地下廠房施工進(jìn)度仿真分析

李明超,段志龍

(天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室,天津 300072)

針對(duì)大型地下廠房施工影響因素眾多、關(guān)系復(fù)雜的問(wèn)題,從影響每個(gè)工序時(shí)間的風(fēng)險(xiǎn)事件出發(fā),量化風(fēng)險(xiǎn)事件,建立耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真模型,突破了地下廠房施工進(jìn)度仿真?zhèn)鹘y(tǒng)模型中以工程整體為對(duì)象進(jìn)行風(fēng)險(xiǎn)識(shí)別的不足,能夠從工序底層考慮風(fēng)險(xiǎn)因素對(duì)施工工期的影響。進(jìn)一步結(jié)合蒙特卡羅模擬方法,實(shí)現(xiàn)了耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真計(jì)算與優(yōu)化控制分析。工程實(shí)例分析表明,該方法不僅可計(jì)算考慮各種不確定風(fēng)險(xiǎn)因素下的大型地下廠房施工工期,而且能夠深入分析工程進(jìn)度對(duì)關(guān)鍵風(fēng)險(xiǎn)因素的敏感性,可為工程風(fēng)險(xiǎn)的防范和施工進(jìn)度的控制分析提供依據(jù)。

水電站;地下廠房;施工進(jìn)度仿真;風(fēng)險(xiǎn)控制;蒙特卡羅方法;敏感性分析

隨著我國(guó)水利水電建設(shè)事業(yè)的飛速發(fā)展和施工技術(shù)的不斷改進(jìn),一批大型的水電站已經(jīng)建成或正在建設(shè),同時(shí)越來(lái)越多地將水電站廠房布置在地下,其規(guī)模也愈來(lái)愈大,如三峽、溪洛渡等水電站。地下廠房洞室群的施工是一項(xiàng)影響因素眾多、關(guān)系復(fù)雜的系統(tǒng)工程,無(wú)論在開(kāi)挖、支護(hù)、襯砌和施工組織等方面都與地面露天作業(yè)有很大的差別,受到不良地質(zhì)構(gòu)造、施工技術(shù)水平、施工機(jī)械配置等眾多因素影響與限制;各單項(xiàng)洞室的施工又是許多個(gè)作業(yè)組成的反復(fù)循環(huán)的過(guò)程,同時(shí)每個(gè)作業(yè)的工作時(shí)間大多是隨機(jī)的,每個(gè)循環(huán)都可以產(chǎn)生隨機(jī)排隊(duì)擁堵現(xiàn)象;地下洞室群縱橫交錯(cuò),布置密集,高差大,施工通道少,在整個(gè)地下洞室群系統(tǒng)的施工過(guò)程中,不但各工序配合與相互干擾錯(cuò)綜復(fù)雜,而且在安排各個(gè)洞室施工先后順序及隧洞施工順序時(shí),需要顧及交通運(yùn)輸?shù)葐?wèn)題的影響[1]。這些不確定性的風(fēng)險(xiǎn)因素往往使得大型地下廠房的施工進(jìn)度難以滿足計(jì)劃安排,導(dǎo)致實(shí)際工期的延長(zhǎng)而造成巨大的社會(huì)和經(jīng)濟(jì)損失。

目前工程項(xiàng)目施工進(jìn)度計(jì)劃與分析主要采用關(guān)鍵線路法(critical path method,CPM)和計(jì)劃評(píng)審技術(shù)(programevaluationandreviewtechnique, PERT)[2]。在CPM中,假定各子工序的施工時(shí)間為一確定值,而實(shí)際上子工序受許多的不確定性風(fēng)險(xiǎn)事件影響。雖然PERT給出了3種可能的時(shí)間估計(jì),即樂(lè)觀的時(shí)間估計(jì)、悲觀的時(shí)間估計(jì)和最可能的時(shí)間估計(jì),彌補(bǔ)了CPM的不足,但還是比較粗糙。Littlefield等[3]認(rèn)為,工序的持續(xù)時(shí)間也可服從其他隨機(jī)分布,如三角分布、Gamma分布等;Hulett[4]通過(guò)分析認(rèn)為三角分布更能反映工序時(shí)間的不確定性。王卓甫等[5]基于影響施工進(jìn)度的不確定性風(fēng)險(xiǎn)事件的分布情況,采用蒙特卡羅方法進(jìn)行仿真計(jì)算;Wang等[6]運(yùn)用蒙特卡羅方法和敏感性分析對(duì)工程全生命期進(jìn)行了管理分析;余建星等[7]提出了在工程風(fēng)險(xiǎn)分析中確定風(fēng)險(xiǎn)評(píng)級(jí)指標(biāo)及標(biāo)準(zhǔn)的原則和方法;Kulk等[8]采用邏輯回歸分析對(duì)風(fēng)險(xiǎn)因素進(jìn)行量化;Liu等[9]在EPC工程中使用了模糊理論來(lái)量化工程中的風(fēng)險(xiǎn)事件,考慮了事件發(fā)生與影響之間的維度問(wèn)題;杜修力等[10]將網(wǎng)絡(luò)分析法應(yīng)用到地下工程風(fēng)險(xiǎn)評(píng)估中,對(duì)各風(fēng)險(xiǎn)因素進(jìn)行分析和運(yùn)算,為后期風(fēng)險(xiǎn)的管理和控制提供了科學(xué)依據(jù);劉東海等[11]通過(guò)柔性循環(huán)網(wǎng)絡(luò)模型,突破了傳統(tǒng)PERT中只考慮工序持續(xù)時(shí)間不確定性的限制,全面考慮了工序的持續(xù)時(shí)間柔性、順序柔性和操作柔性;鐘登華等[12]通過(guò)對(duì)比削峰填谷模型、最小方差模型和遺傳算法模型對(duì)實(shí)際工程的優(yōu)化結(jié)果,得出遺傳算法模型在處理工期一定、資源均衡的資源進(jìn)度計(jì)劃問(wèn)題上,相對(duì)最優(yōu);這些研究成果從影響工程施工進(jìn)度的不同因素角度來(lái)模擬分析工程進(jìn)度安排和優(yōu)化方案,取得了相應(yīng)的成果;但其中對(duì)風(fēng)險(xiǎn)的識(shí)別主要是以工程整體為對(duì)象,較少?gòu)娘L(fēng)險(xiǎn)驅(qū)動(dòng)影響工序時(shí)間的角度出發(fā)[13-14]進(jìn)行施工進(jìn)度方案分析。本文從影響底層工序工期的角度進(jìn)行風(fēng)險(xiǎn)識(shí)別,建立耦合風(fēng)險(xiǎn)的進(jìn)度仿真模型,確定風(fēng)險(xiǎn)的發(fā)生概率及影響值,以風(fēng)險(xiǎn)驅(qū)動(dòng)理論進(jìn)行仿真,并對(duì)工程施工進(jìn)度進(jìn)行優(yōu)化分析。

1 仿真建模

1.1 地下廠房施工過(guò)程風(fēng)險(xiǎn)因素分解

水電站地下廠房系統(tǒng)布置較為集中,主廠房、主變室和尾調(diào)室等規(guī)模很大;同時(shí),施工支洞在空間布置上縱橫交錯(cuò),形成多處交叉路口。地下洞室群施工十分復(fù)雜,作業(yè)面多,多個(gè)工序平行作業(yè),且機(jī)械化程度高,施工中產(chǎn)生的有害物質(zhì)不僅有爆破后的氣體,還有裝載機(jī)和汽車(chē)排放的氣體、煙氣,混凝土作業(yè)中產(chǎn)生的粉塵等,需要及時(shí)通過(guò)通風(fēng)系統(tǒng)排出地下廠房;對(duì)施工運(yùn)輸問(wèn)題若不能合理地規(guī)劃或組織安排不當(dāng),極易造成車(chē)輛運(yùn)行不暢,在交通要道或瓶頸部位往往形成車(chē)輛堵塞而影響施工的正常進(jìn)行,嚴(yán)重影響工期。由于地下廠房施工的復(fù)雜性,從每一個(gè)工序出發(fā),找到影響其工期的風(fēng)險(xiǎn)因素,能夠更為客觀地反映其實(shí)際施工過(guò)程。為了能夠識(shí)別地下廠房的施工風(fēng)險(xiǎn),本文采用分解分析法[15],將復(fù)雜的大系統(tǒng)分解成若干便于識(shí)別的簡(jiǎn)單小系統(tǒng),即將地下廠房的每個(gè)施工工序按照施工順序再分解成更小的更容易識(shí)別的工序,同時(shí)結(jié)合施工經(jīng)驗(yàn),從而識(shí)別出施工中存在的各種主要風(fēng)險(xiǎn)。圖1為經(jīng)分析總結(jié)后得出的影響地下廠房施工進(jìn)度的主要風(fēng)險(xiǎn)因素,各主要風(fēng)險(xiǎn)因素的具體說(shuō)明見(jiàn)表1,表中的風(fēng)險(xiǎn)因素可適用于同類(lèi)工程。

圖1 地下廠房施工風(fēng)險(xiǎn)因素分解結(jié)構(gòu)示意圖

表1 主要風(fēng)險(xiǎn)因素

1.2 地下廠房施工進(jìn)度仿真模型

大型地下廠房開(kāi)挖施工是一個(gè)相當(dāng)復(fù)雜的過(guò)程,作業(yè)面多,鉆孔、爆破、出渣、運(yùn)輸、噴錨支護(hù)、二次襯砌等多個(gè)工序平行作業(yè),在進(jìn)行耦合風(fēng)險(xiǎn)的施工進(jìn)度模擬計(jì)算分析時(shí),要詳細(xì)分析其開(kāi)挖的過(guò)程,建立如圖2所示的開(kāi)挖工序網(wǎng)絡(luò)層。該網(wǎng)絡(luò)層能夠從底層關(guān)聯(lián)工序時(shí)間的不確定因素,考慮風(fēng)險(xiǎn)因素對(duì)施工工期的影響,是建立仿真模型的基礎(chǔ)。

圖2 地下廠房開(kāi)挖工序網(wǎng)絡(luò)層

式中:Rk為工序k所對(duì)應(yīng)的風(fēng)險(xiǎn)因子;Pi、Ii分別為與工序k相關(guān)的第i個(gè)風(fēng)險(xiǎn)事件的發(fā)生概率和影響值;n為影響工序k的風(fēng)險(xiǎn)事件的總數(shù)。

2.2 考慮風(fēng)險(xiǎn)的蒙特卡羅模擬分析

蒙特卡羅模擬方法彌補(bǔ)了傳統(tǒng)方法受人為因素影響大的不足,以概率統(tǒng)計(jì)理論為基礎(chǔ),以隨機(jī)抽樣為主要手段,利用隨機(jī)數(shù)產(chǎn)生抽樣樣本進(jìn)行統(tǒng)計(jì)分析,以求得統(tǒng)計(jì)值作為待解問(wèn)題的數(shù)值解[16]。耦合風(fēng)險(xiǎn)因素的蒙特卡羅模擬步驟如下:

a.在每一次迭代過(guò)程中,假定第i個(gè)風(fēng)險(xiǎn)事件發(fā)生的概率為Pi,則在[0,1]之間生成一個(gè)隨機(jī)數(shù)ri,如果ri＜Pi,則認(rèn)為該風(fēng)險(xiǎn)發(fā)生,同時(shí)認(rèn)為風(fēng)險(xiǎn)事件對(duì)工序的影響值Ii具有不確定性,由其三角分布確定。

b.根據(jù)式(1)計(jì)算出每個(gè)工序的風(fēng)險(xiǎn)因子。

c.將步驟b中的風(fēng)險(xiǎn)因子加1,再乘以其相應(yīng)工序的工期,即為該次仿真計(jì)算中工序的工期,計(jì)算公式為質(zhì)構(gòu)造”。由于同一類(lèi)型風(fēng)險(xiǎn)事件在不同工程中發(fā)生的概率具有較強(qiáng)的隨機(jī)性,該參數(shù)采用歷史資料統(tǒng)計(jì)和專(zhuān)家打分相結(jié)合的方法來(lái)確定。

b.風(fēng)險(xiǎn)事件的影響值:是指風(fēng)險(xiǎn)發(fā)生時(shí),對(duì)相關(guān)工序工期的影響大小,取值區(qū)間為[-1,1],負(fù)值代表對(duì)工期有正影響,正值對(duì)工期有負(fù)影響。風(fēng)險(xiǎn)事件的影響值采用三角分布表示,該參數(shù)可由表1定義的風(fēng)險(xiǎn)特征值統(tǒng)計(jì)得到,但應(yīng)注意同一種風(fēng)險(xiǎn)事件對(duì)不同工程影響值不同,具體要結(jié)合實(shí)際工程進(jìn)行修正。進(jìn)而引入風(fēng)險(xiǎn)因子,是指總的風(fēng)險(xiǎn)事件對(duì)工序的影響值,計(jì)算公式如下:

采用施工過(guò)程循環(huán)網(wǎng)絡(luò)理論[1],根據(jù)地下廠房實(shí)際施工特點(diǎn),考慮地下廠房施工過(guò)程的主要風(fēng)險(xiǎn)因素,建立其仿真模型。圖3給出了一個(gè)示例模型,該模型將相關(guān)風(fēng)險(xiǎn)因素添加到相應(yīng)的施工工序中,確定影響地下廠房施工進(jìn)度每個(gè)工序風(fēng)險(xiǎn)事件的發(fā)生概率和影響值,即可得到耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真模型。

2 仿真計(jì)算流程

2.1 風(fēng)險(xiǎn)事件定量分析

假定風(fēng)險(xiǎn)事件具有兩個(gè)參數(shù),即風(fēng)險(xiǎn)事件的發(fā)生概率(possibility of occurrence)及影響值(impact value),通過(guò)實(shí)際調(diào)查統(tǒng)計(jì),確定其定量屬性。

a.風(fēng)險(xiǎn)事件的發(fā)生概率:是指在每一次迭代中風(fēng)險(xiǎn)事件發(fā)生的可能性,用百分?jǐn)?shù)表示。實(shí)際工程中,有些風(fēng)險(xiǎn)事件是一定發(fā)生的,如“進(jìn)度計(jì)劃執(zhí)行不佳”;有些風(fēng)險(xiǎn)事件是有可能發(fā)生的,如“不良地

圖3 耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真模型

式中:T′k為考慮風(fēng)險(xiǎn)事件影響后的工序k的工期; Tk為未考慮風(fēng)險(xiǎn)事件影響的工序k的工期。

d.計(jì)算出各個(gè)工序的工期,找出關(guān)鍵路徑,計(jì)算出本次仿真的總工期。

e.重復(fù)上述步驟,得到該工程的多次仿真計(jì)算結(jié)果。

f.確定總工期的分組組距Tg,若仿真計(jì)算結(jié)果中總工期的最大值為T(mén)max,最小值為T(mén)min,l為分組數(shù),則組距Tg的計(jì)算公式為

g.確定區(qū)間邊界,假定T1=Tmin+Tg,T2=T1+Tg,…,Tl=Tl-1+Tg,那么區(qū)間結(jié)果為[Tmin,T1),[T1, T2),…,[Tl-1,Tl)。

h.統(tǒng)計(jì)落入各個(gè)區(qū)間的頻數(shù),計(jì)算頻率值。

i.根據(jù)步驟h的計(jì)算結(jié)果,畫(huà)出總工期的頻率分布圖和完工概率曲線。

2.3 施工過(guò)程風(fēng)險(xiǎn)因素的敏感性分析

敏感性分析是在項(xiàng)目評(píng)價(jià)和決策中常用的一種不確定分析方法,是用來(lái)測(cè)定一個(gè)或多個(gè)因素的變化對(duì)目標(biāo)的影響程度。地下廠房的施工過(guò)程受多種風(fēng)險(xiǎn)因素影響,但要確定對(duì)施工進(jìn)度影響最大的風(fēng)險(xiǎn)因素,即關(guān)鍵風(fēng)險(xiǎn)事件,需進(jìn)行敏感性分析。敏感性分析方法分為單因素敏感分析和多因素敏感分析,本文采用單因素敏感性分析確定關(guān)鍵風(fēng)險(xiǎn)事件,為后面進(jìn)行有效合理地控制,達(dá)到優(yōu)化施工進(jìn)度的目的打基礎(chǔ)。敏感性分析的步驟如下:

a.針對(duì)項(xiàng)目的主要風(fēng)險(xiǎn)事件,確定其發(fā)生概率及影響值的最大變動(dòng)幅度。

b.在所有風(fēng)險(xiǎn)事件變動(dòng)的幅度內(nèi),確定都可以達(dá)到的一個(gè)變動(dòng)幅度。

c.以步驟b得出的變動(dòng)幅度,逐一改變風(fēng)險(xiǎn)事件的發(fā)生概率和影響值等仿真參數(shù),其他風(fēng)險(xiǎn)事件不變,得到同一風(fēng)險(xiǎn)事件在不同發(fā)生概率下的施工進(jìn)度方案。

d.對(duì)步驟c得到的施工進(jìn)度方案進(jìn)行進(jìn)度仿真,得到施工總工期,用總工期除以變動(dòng)前的工期得到風(fēng)險(xiǎn)事件在變動(dòng)幅度相同的情況下對(duì)總工期的影響值,影響值越大說(shuō)明對(duì)總工期的影響越大。

2.4 影響施工進(jìn)度的風(fēng)險(xiǎn)控制分析

風(fēng)險(xiǎn)控制的方法主要有風(fēng)險(xiǎn)回避、風(fēng)險(xiǎn)降低、風(fēng)險(xiǎn)抵消、風(fēng)險(xiǎn)分離、風(fēng)險(xiǎn)分散、風(fēng)險(xiǎn)轉(zhuǎn)移、風(fēng)險(xiǎn)自留等[17];同時(shí),風(fēng)險(xiǎn)事件主要涉及機(jī)械、人員、管理和外界影響等4個(gè)方面,針對(duì)這4個(gè)方面分別采用相應(yīng)的控制措施,如對(duì)機(jī)械設(shè)備定期進(jìn)行檢修、維護(hù),加強(qiáng)人員培訓(xùn)、提高專(zhuān)業(yè)技能,建立多種管理體系等,可減少風(fēng)險(xiǎn)事件對(duì)總工期的影響。

圖4為耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真總的計(jì)算流程。

圖4 仿真計(jì)算流程

圖5 某水電站地下廠房布置

3 工程實(shí)例

3.1 工程概況

某水電站地下廠房系統(tǒng)主要由進(jìn)水口、引水隧洞、地下主副廠房、母線洞、主變室、通風(fēng)洞、尾水調(diào)壓室、尾水隧洞等組成,其中主廠房長(zhǎng)396.0 m,上部頂拱跨度31.0m,下部邊墻跨度29.0 m,高77.47 m,規(guī)模較大。各洞室主要位于花崗巖巖體中,受到斷層構(gòu)造的影響。該大型地下廠房的三維布置方案如圖5所示。由于主廠房施工是整個(gè)系統(tǒng)的關(guān)鍵,本文以主廠房為分析對(duì)象,其施工進(jìn)度網(wǎng)絡(luò)仿真模型如圖6所示,其中圓形圖標(biāo)為確定性工序,無(wú)須仿真計(jì)算;方形圖標(biāo)為仿真工序,該工序的持續(xù)時(shí)間將由耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真模型計(jì)算來(lái)確定。

3.2 風(fēng)險(xiǎn)事件定量分析

結(jié)合表1所示的風(fēng)險(xiǎn)特征值,通過(guò)對(duì)歷史資料進(jìn)行統(tǒng)計(jì)分析,并結(jié)合專(zhuān)家打分結(jié)果,得到影響該地下廠房施工進(jìn)度的風(fēng)險(xiǎn)事件的發(fā)生概率和影響值如表2所示。

圖6 主廠房施工進(jìn)度網(wǎng)絡(luò)仿真模型

表2 風(fēng)險(xiǎn)事件參數(shù)

3.3 仿真計(jì)算與分析

基于所建立的主廠房施工進(jìn)度網(wǎng)格仿真模型和相應(yīng)的風(fēng)險(xiǎn)事件及其參數(shù),總工期在1805～2657 d之間,平均值為2144 d,標(biāo)準(zhǔn)差為123 d,圖7為仿真總工期的概率分布圖和完工概率曲線,其概率分布符合β分布。傳統(tǒng)CPM的計(jì)算工期為1783 d,從圖7分析可知,該工期實(shí)際施工難以達(dá)到,所以有必要在開(kāi)工前進(jìn)行風(fēng)險(xiǎn)控制分析,對(duì)各種風(fēng)險(xiǎn)事件進(jìn)一步分析,找到影響施工進(jìn)度的關(guān)鍵風(fēng)險(xiǎn)事件,并提出相應(yīng)的控制措施,為工程決策提供依據(jù)。

對(duì)各種風(fēng)險(xiǎn)事件進(jìn)行敏感性分析,計(jì)算結(jié)果如表3所示,表中的回歸系數(shù)表示風(fēng)險(xiǎn)事件對(duì)總工期的影響大小,即風(fēng)險(xiǎn)事件在相同的變化幅度下引起總工期的變化幅度。從表3可知影響該工程施工進(jìn)度的關(guān)鍵風(fēng)險(xiǎn)事件是“不良地質(zhì)構(gòu)造”,該關(guān)鍵風(fēng)險(xiǎn)事件屬于外界影響,需采用多種技術(shù)手段,對(duì)該地區(qū)的地質(zhì)情況加以詳細(xì)分析,以減少地質(zhì)情況的不確定性對(duì)工期的影響。

圖7 仿真總工期概率分布圖和完工概率曲線

表3 風(fēng)險(xiǎn)事件敏感性分析

實(shí)際施工中可通過(guò)地質(zhì)勘探數(shù)據(jù)和地質(zhì)建模等手段來(lái)提前預(yù)報(bào)施工區(qū)的地質(zhì)情況。結(jié)合之前實(shí)例工程資料,該風(fēng)險(xiǎn)事件的仿真參數(shù)可以控制在一定范圍內(nèi),見(jiàn)表4;經(jīng)過(guò)計(jì)算得到控制前后總工期的完工概率曲線,如圖8所示,總工期的平均值從控制前的2144d縮短到2100d,減少了44d,施工進(jìn)度得到了一定程度的優(yōu)化。

表4 關(guān)鍵風(fēng)險(xiǎn)事件控制前后數(shù)據(jù)

圖8 控制前后總工期的累計(jì)概率曲線

4 結(jié) 語(yǔ)

本文從影響地下廠房施工工序時(shí)間的風(fēng)險(xiǎn)事件出發(fā),量化風(fēng)險(xiǎn)事件,建立了地下廠房開(kāi)挖網(wǎng)絡(luò)模型,結(jié)合傳統(tǒng)網(wǎng)絡(luò)模型,建立了耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真模型,突破了傳統(tǒng)模型中對(duì)風(fēng)險(xiǎn)分析的不足,能夠從工序底層考慮風(fēng)險(xiǎn)因素對(duì)施工工期的影響。結(jié)合蒙特卡羅模擬方法,提出了耦合風(fēng)險(xiǎn)因素的地下廠房施工進(jìn)度仿真計(jì)算與優(yōu)化控制分析方法。以某水電站地下廠房為實(shí)例開(kāi)展仿真分析,得到了總工期的概率分布圖和完工概率曲線,可確定合理的總工期;通過(guò)敏感性分析確定關(guān)鍵風(fēng)險(xiǎn)事件,并以“不良地質(zhì)構(gòu)造”關(guān)鍵風(fēng)險(xiǎn)為例進(jìn)行控制分析,優(yōu)化了工程進(jìn)度,為工程風(fēng)險(xiǎn)的防范和施工進(jìn)度的控制分析提供了科學(xué)依據(jù)。

[1]鐘登華,劉東海.大型地下洞室群施工系統(tǒng)仿真理論方法與應(yīng)用[M].北京:中國(guó)水利水電出版社,2003.

[2]巴延廷.基于PERT的工程進(jìn)度風(fēng)險(xiǎn)分析研究[D].大連:大連理工大學(xué),2010.

[3]LITTLEFIELD T K,RANDOLPH P H.Reply:an answer to sasieni’s question on PERT times[J].Management Science,1987,33(10):1357-1359.

[4]HULETT D.Practical schedule risk analysis[M].London: MPG Books Group,2009.

[5]王卓甫,陳登星.水利水電施工進(jìn)度計(jì)劃的風(fēng)險(xiǎn)分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,1999,27(4):83-87. (WANGZhuofu,CHENDengxing.Riskanalysisof construction schedule for water resource projects[J]. Journal of Hohai University:Natural Sciences,1999,27 (4):83-87.(in Chinese))

[6]WANG N,CHANG Y,EL-SHEIKH A.Monte Carlo simulation approach to life cycle cost management[J]. Structure and Infrastructure Engineering,2012,8(8):739-746.

[7]余建星,李成.工程風(fēng)險(xiǎn)分析中的風(fēng)險(xiǎn)當(dāng)量及其評(píng)價(jià)標(biāo)準(zhǔn)[J].海洋技術(shù),2004,23(1):48-51.(YU Jianxing,LI Cheng.The risk equivalent and its standard of estimate in engineering risk analysis[J].Ocean Technology,2004,23 (1):48-51.(in Chinese))

[8]KULK G P,PETERS R J,VERHOEF C.Quantifying IT estimation risks[J].Science of Computer Programming, 2009,74(11/12):900-933.

[9]LIU J,WU S,ZHAO X,et al.Application of fuzzy theory in contractor risk assessment under EPC model[C]// International Conference on Information Engineering and ComputerScience.Piscataway,NewJersey:IEEE Publisher,2009:1-5.

[10]杜修力,高云昊,張明聚,等.網(wǎng)絡(luò)分析法及其在地下工程風(fēng)險(xiǎn)分析中的應(yīng)用[J].土木工程學(xué)報(bào),2010,43(增刊2):353-357.(DU Xiuli,GAO Yunhao,ZHANG Mingju,et al.Analytic network process and its application in risk assessment in underground engineering[J].China Civil Engineering Journal,2010,43(Sup2):353-357.(in Chinese))

[11]劉東海,黃培志.隧洞TBM施工過(guò)程柔性循環(huán)網(wǎng)絡(luò)建模與仿真分析[J].土木工程學(xué)報(bào),2010,43(10):115-121.(LIU Donghai,HUANG Peizhi.Flexible cyclic operation network modeling and simulation analysis of TBM construction process[J].China Civil Engineering Journal,2010,43(10):115-121.(in Chinese))

[12]鐘登華,李明超,張偉波,等.復(fù)雜工程施工系統(tǒng)資源優(yōu)化模型及其應(yīng)用[J].天津大學(xué)學(xué)報(bào):自然科學(xué)與工程技術(shù)版,2004,37(7):589-594.(ZHONG Denghua,LI Mingchao,ZHANG Weibo,et al.Resource optimization models and application for complex engineering construction system[J].Journal of Tianjin University:Science and Technology,2004,37(7):589-594.(in Chinese))

[13]郭宇,劉爾烈.應(yīng)用蒙特卡羅方法改進(jìn)項(xiàng)目成本風(fēng)險(xiǎn)分析[J].天津大學(xué)學(xué)報(bào):自然科學(xué)與工程技術(shù)版,2002, 35(2):199-202.(GUO Yu,LIU Erlie.Two improvements in risk analysis for project cost using Monte Carlo methods [J].JournalofTianjinUniversity:Scienceand Technology,2002,35(2):199-202.(in Chinese))

[14]金峰.基于風(fēng)險(xiǎn)驅(qū)動(dòng)理論的費(fèi)用風(fēng)險(xiǎn)分析方法[J].工程管理學(xué)報(bào),2012,26(1):75-78.(JIN Feng.Project cost risk analysis based on risk driven theory in international engineeringprojects[J].JournalofEngineering Management,2012,26(1):75-78.(in Chinese))

[15]鄧鵬.青草沙水庫(kù)工程的風(fēng)險(xiǎn)識(shí)別[J].中國(guó)勘察設(shè)計(jì), 2010(9):69-72.(DENG Peng.Risk identification for Qincaosha reservoir works[J].ChineseSurveyand Design,2010(9):69-72.(in Chinese))

[16]李景茹.大型工程施工進(jìn)度分析理論方法與應(yīng)用[D].天津:天津大學(xué),2003.

[17]宋永嘉,張淙皎,田林鋼,等.水利工程風(fēng)險(xiǎn)量化問(wèn)題的探討[J].人民黃河,2004(9):33-34.(SONG Yongjia, ZHANG Congjiao,TIAN Lingang,et al.The discussing on thequestionsaboutriskquantificationofhydraulic engineering[J].Yellow River,2004(9):33-34.in Chineses))

Construction scheduling simulation integrated with risks of large-scale underground powerhouse//

LI Mingchao,DUAN Zhilong
(State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China)

The construction of large-scale underground powerhouse is very complicated and easily affected by risk events. The current methods for analyzing the risks are based on the whole construction,while the risk events directly affect the stages of production.In order to consider the influence of risk events,the integrated risk-scheduling model of large-scale underground powerhouse was established.Further,the process of risk-driven construction scheduling and optimal controlling was introduced by combining the Monte-Carlo method.The results,based on engineering examples,show that that the proposed method can not only calculate the construction schedule of large-scale underground powerhouse with many risk factors,but also achieve in-depth analysis of the sensitivity of key risk factors to the construction schedule.The simulation results can provide effective support to the prevention of project risks and the control of construction schedule.

hydropower station;underground powerhouse;construction schedule simulation;risk control;Monte-Carlo method;sensitivity analysis

TV512

:A

:1006-7647(2014)05-0039-06

10.3880/j.issn.1006-7647.2014.05.008

2013-1118 編輯:熊水斌)

國(guó)家自然科學(xué)基金(51379006);國(guó)家創(chuàng)新研究群體科學(xué)基金(51021004);教育部新世紀(jì)優(yōu)秀人才計(jì)劃(NCET-12-0404)

李明超(1979—),男,湖北仙桃人,副教授,博士,主要從事水利水電工程仿真分析研究。E-mail:lmc@tju.edu.cn