[发明专利]一种工程进度中不确定风险的分析方法

权利要求书

1.一种工程进度中不确定风险的分析方法，其特征在于，包括步骤：

步骤一、获取项目所在地近五年内至少1个自然年的降雨气象资料，包括每日或半日降雨量情况及雨量大小；

步骤二、对影响工期的每日或半日降雨量大小进行定义：当降雨日内统计的12小时内雨量大于5毫米或24小时内雨量大于10毫米时，工期影响1天；

步骤三、蒙特卡罗法模拟所需数据的生成：假设第i月第k日12小时内的降雨量为m_1k，24小时内的降雨量为m_2k，当日是否影响工期记为X_k，当m_1k＞5或m_2k＞10时，X_k记为1，否则记为0，则有该月的降雨数据生成一组数X₁，X₂，...，X_n，n为该月的天数；

步骤四、蒙特卡罗法模拟及降雨影响计算：基于降雨数据生成的一组数X₁，X₂，...，X_n，n为该月的天数，采用蒙特卡罗法进行模拟m次，生成的模拟结果为X′₁，X'₂，...，X'_m，则有影响工期的降雨发生概率降雨影响天数为np_i；

步骤五、工程进度不确定风险分析及风险表达式推求：

由于降雨存在不确定性，对工期的影响亦存在一定不确定性，对不确定性引起的风险通过概率的大小进行判定，计算项目实际工期超过计划工期的概率；

若计划工期为T_p，绝对工期为T_a，降雨影响工期为T_m，则项目不能按计划完工的风险可表达为：

P_r＝P(T_m+T_aT_p)    (1)

由概率学可知，任意一天施工是否受降雨影响服从参数为p的(0-1)分布，即X_k～(0-1)，其中p为降雨影响概率；

设包含降雨影响工期在内的工期为T，降雨影响工期T_m～B(T，p)；借助T内各月的模拟试验确定式(1)的P_r，设第i月天数为T_i，则第i月降雨影响工期∑X_ik～B(T_i，p_i)，由中心极限定理，∑X_ik近似服从正态分布，即其中：

μ_i＝T_ip_i      (2)

式中:μ_i—第i月降雨影响天数的期望值；

—第i月降雨影响天数的方差；

若项目的工期有s个月，则降雨影响工期T_m＝∑(X_1k+X_2k+...+X_sk)，且由正态分布的可加性，得：

其中：

μ_s＝μ_s1+μ_s2+...+μ_ss

＝T₁p₁+T₂p₂+...+T_sp_s

≈T₀(p₁+p₂+...+p_s)       (5)

式中:μ_s、μ_s1、μ_s2、μ_ss—s个月内、第1月、第2月、第s月降雨影响天数的期望值；

—s个月内、第1月、第2月、第s月降雨影响天数的方差；

p₁、p₂…p_s—第1月、第2月…第s月的降雨影响概率；

T₁、T₂…T_s—第1月、第2月…第s月天数，可取30天；

T₀—各月天数均值，可取30天；

则由式(1)，可得：

P_r＝1-P(T_m≤T_p-T_a)

＝1-Φ(((T_p-T_a)-μ_s)/σ_s)      (7)；

步骤六、项目计划工期内降雨影响天数期望值、降雨影响天数方差计算：

利用步骤五中公式(5)、公式(6)计算得到项目计划工期内降雨影响天数期望值μ_s、降雨影响天数方差

步骤七、工程不能按期完工风险率计算：

将计划工期T_p、绝对工期T_a、降雨影响天数期望值μ_s、降雨影响天数方差代入公式(7)，查正态分布表便可得到工程不能按照计划完工的概率。