工程网络计划蒙特卡洛仿真及进度风险分析.docxVIP

工程网络计划蒙特卡洛仿真及进度风险分析 工程网络计划蒙特卡洛仿真及进度风险分析 第21卷　第4期 文章编号:1006-9348(2021)04-0143-05 计　算　机　仿　真 2021年4月 工程网络计划蒙特卡洛仿真及进度风险分析 王仁超,欧阳斌,褚春超 (天津大学建筑工程学院,天津300072) 摘要:运用网络计划蒙特卡洛仿真(MCS),从宏观的项目层和微观的工序层分别进行进度风险分析。在项目层,仿真计算获取总工期的统计分布,据此计算某一计划工期的项目完工概率或完工风险,工期;在工序层,通过计算工序关键概率(ACP)和关键指标(ACI),掌握关键线路、关键工序分布情况,外,用VisualC++开发了一个网络计划MCS软件,关键词:网络计划;蒙特卡洛仿真;进度风险;项目总工期;;中图分类号:TP391.9　　文献标识码:A ,CPM)和计划评审技术(PERT)。CPM适用于工程中各工序间逻辑关系和持续时间都确定的情况。然而,实际工程项目中通常包含诸多不确定性因素,进而导致工程进度计划存在一定的风险,因而有必要对进度风险进行分析。PERT为各工序持续时间不确定而逻辑关系确定的情况,提供了简捷的处理方法,即:基于“三时估计法”计算各工序持续时间的期望值D及方差σ,即D=(a+4m+b)/6,σ=(b-a)/36(其中a、最可能、最乐观的工序估计持续时间);m、b分别为最悲观、然后可按CPM处理,并计算方差、完工概率等。然而,这些简化处理基于以下简化假定 。而这些可分别借助各工序ACP和 ACI指标分析来实现。因此,本文提出进度风险分析应兼顾 宏观的项目层与微观的工序层,即关注项目整体进度风险的同时也应重视对风险源工序的辨识,并开发相应网络计划 MSC软件予以实现。 2　网络计划MCS基本原理 2.1　网络计划时间参数的计算 网络计划按图形表现形式可分为单代号网络(AoE)和双代号网络(AoN)。由于单代号网络具有建模直观而简捷、网络图易于绘制、便于计算机编程实现等优点,本文研究采用这种形式。假定项目共有n项工序,为建立一个具有单一始点和终点、无回环的项目网络模型,统一添加一个虚拟的开工工序和完工工序。网络中各节点(工序)时间参数及工期的计算公式如下[1]: 工序最早开始时间ES: ESi=max{ESh+Dh}(i=1,…,n+1 :①各工序持续时间服从β分 布,且相互独立;②网络计划图中只有一条线路(期望关键线路)占支配地位,其它线路成为关键线路的概率可忽略不计;③网络计划中线路的持续时间服从正态分布。而事实上,没有任何的经验证明工序持续时间的最适合形式是β分布 实际中,由于工序持续时间的不确定性,几乎所有的线路理论上都有可能成为关键线路[4],导致实际网络中经常出现两条或多条关键、准关键线路的情况[3];线路的持续时间服从正态分布的假定以线路间相互独立以及中心极限定理适用为前提,与工程实际存在差距。 工序最早完成时间EF: EFi=ESi+Di 针对上述问题,通常在网络计划中引入MCS。然而,纵观国内网络计划MCS及其进度风险分析的相关研究[3][5][6],往往只关注于整个项目完工概率或风险值。这些分析对于宏观决策固然有着重要的作用,但网络计划的具体实施与管理是以工序为生命线的,倘若进度风险分析仅局限于项目层显然是片面的、缺乏实践指导意义的。对于计划的具体实施者和管理者,不仅需要了解网络计划内部的关键线路和关键工序的分布情况,还需甄别出其中的高风险工序作为进度风 网络计划的计算工期Tc: Tc=EFn+1 工序最迟完成时间LF: LFn+1=EFn+1 LFi=min{LFj-Dj}(i=n,…,0 工序最迟开始时间LS: LSi=LFi-Di 收稿日期:2021-04-22 工序总时差TF: TFi=LFi-EFi=LSi-ESi(6) F(T)(也称为完工概率曲线)。据此可求(8)式近似表达式: Pr=1-P(T≤Tp)=1-F3(Tp) 工序自由时差FF: FFi=min{ESj}-EFi 于是,由完工概率曲线可进一步获得进度风险曲线。利用完工概率曲线和进度风险曲线,根据总工期T与完工概率或进度风险率的对应关系,则可进行如下进度风险分析: 1)输入指定某一总工期,查询相应完工概率和进度风 上述各式中:下标i(i=0,…,n+1)表示各工序的ID号,其中0、n+1分别为虚拟开工、完工工序