人工岛建设相比于内陆工程建设,具有施工环境复杂、不确定因素多及技术难度高等特点。考虑到此类工程规模庞大、安全指标要求高且管理困难,必须对工程进行系统的风险分析,识别出各工作阶段的风险因素,并作出综合性评价。
国内对于人工岛建设风险的研究取得了较多的成果。宋君妍等[1]通过对人工岛开发建设过程的风险分析与对策研究,提出了有效的风险预防及控制措施。冯民学等[2]对洋口港区人工岛进行了雷击风险计算和风险区等级划分。王彦林等[3]对港珠澳大桥东、西人工岛施工图设计进行了风险评估,详细讨论了安全风险的防范建议。以上大多注重于整体工程的综合评价,对于各分项工程风险识别方法的思考并不明确。WBS-RBS是一种广泛应用于风险分析的方法,已在较多实际工程的风险评价中得到了应用。陈国华等[4]运用WBS-RBS方法对跨海桥梁施工作业和风险层次进行了分解,确定出了对应工程的风险等级和主要影响因素。贾俊峰等[5]通过构建WBS-RBS结构提高了土建工程施工风险评估的准确性。 笔者结合港珠澳大桥人工岛工程实例,运用WBS-RBS风险辨识方法和AHP方法计算并分析工程建设过程中的风险因素及权重,并按照风险度理论及风险标准对各阶段风险水平进行综合评价,详细分析了各分项工程中的风险因素,规避了其他方法分析不全面的弊端,使该工程的风险识别更加清晰和直观,从而提高风险评价结果的可靠性。
港珠澳大桥东西人工岛是连接东西桥梁和海底隧道的重要结合点,更是集办公、监控、隧道救援养护等多功能为一体的综合项目,为港珠澳大桥工程全线路段的重点配套项目。东、西人工岛均采用椭圆形布设,东人工岛岛长625 m,宽225 m,总面积为10.316万m2;西人工岛岛长625 m,宽185 m,总面积为9.796万m2。
由于人工岛施工范围内有800万吨海床淤泥,施工团队采用了“钢筒围岛”方案:在陆地上预先制造120个直径22.5 m、高度55 m、重量达550 t的巨型圆形钢筒,通过船只将其直接固定在海床上,然后在钢筒合围的中间区域进行填土造岛,其施工工法和八锤联动液压振动锤为世界首创[6]。人工岛分别通过岛桥结合部和岛隧结合部实现与桥、隧的连接。
港珠澳大桥人工岛建设地处外海,技术新、经验少,不确定因素较多,容易引起风险的因素也很多,同一风险因素在不同的建设阶段中也会不同。基于各风险因素的时效性和动态特性,对工程分别进行项目结构分解和风险分解,并构造WBS-RBS矩阵对人工岛建设过程中可能存在的风险进行识别。
项目分解(WBS)[7-8]是从全局出发,将一个整体项目分解成若干个相互独立、易于描述的作业单元,可以深入了解和牢牢把握项目实施的具体细节。项目分解的方法比较灵活,可以从实施过程、空间位置、功能或要素等方面入手进行分解。考虑到港珠澳大桥人工岛建设项目的特点,按实施过程进行WBS分解,如图1所示。
图1 港珠澳大桥人工岛建设WBS分解结果
Figure 1 Result of WBS decomposition of Hong Kong-Zhuhai-Macao Bridge
风险分解(RBS)[7-8]是将工程实施过程中可能遇到的风险逐层分解,从而得到不同层次的子风险。在参照相关工程的基础上,充分考虑各方面致险因素,并结合生态文明建设思想,将风险分解为自然环境风险、海洋工程风险、岛内施工风险、环境污染风险和其他风险共5大类,并逐层进一步分解,如图2所示。
图2 港珠澳大桥人工岛建设RBS分解结果
Figure 2 Result of RBS decomposition of Hong Kong-Zhuhai-Macao Bridge
为充分考虑港珠澳大桥人工岛建设工作的实际需要,并使分解结果便于计算,可建立WBS-RBS风险识别矩阵,对各级子工作中所存在的风险因素进行有效识别和清晰反馈,进而全方位、多角度地反映该建设过程中的风险状况。
基于以上理论,构造WBS-RBS风险识别矩阵A。矩阵A的列分别对应各2级子工作,行分别对应各2级子风险。对于确实存在的风险因素,矩阵对应元素设为1。对于比较模糊、可能存在的风险项,也认定其存在,然后根据后期的专家打分确定其权重和重要性。而对于子工作中几乎不可能发生的或对整个工程影响甚小的风险项,则可进行忽略,即将对应元素设为0。
A=
层次分析法(AHP)[9]适用于定性问题的定量分析,可以将复杂问题逐层分解为更加系统、直观的递阶层次结构。结合WBS-RBS判断矩阵,采用层次分析法,确定各分项元素相对于目标的重要程度。具体实施过程如下:
(1)构建各风险元素相对于准则层的判断矩阵:
B=(bij)n×n,
(1)
式中:bij表示对于上一层对象而言,风险因素Ri相对于Rj的重要程度。
bij取值可根据美国运筹学家T.L.Saaty提出的1~9标度法则来确定[7]。bij具体取值及意义见表1。若i=j,则bij=1;反之,bij=1/bji。
表1 bij取值及意义
Table 1 Value and significance of bij
bij取值意义1i与j一样重要3i比j稍微重要5i比j明显重要7i比j强烈重要9i比j极其重要2、4、6、8上述中间状态的对应取值
(2)根据判断矩阵B,确定各准则下风险因素的相对权重。首先求出判断矩阵B的最大特征值λmax,再求出最大特征值对应的特征向量U=(u1,u2,…,un),并对该特征向量进行归一化处理得到向量V,V=(v1,v2,… vn)。归一化处理后向量中的值即为其对应元素的相对权重[10]。向量V的计算如式(2)所示:
(2)
(3)利用一致性指标CR对矩阵进行一致性检验。当CR小于0.1时,认为判断矩阵的一致性满足要求[11]。
(3)
(4)
式中:CI为相容性指标;λmax为判断矩阵的最大特征值;n为判断矩阵的阶数;RI为随机性指标,可参照美国运筹学家T.L.Saaty给定的矩阵阶数与随机性指标参考值[12],见表2。
表2 随机性指标RI取值
Table 2 Value of the mean random consistency index RI
阶数12345678RI000.580.91.121.241.321.41
风险度[11]是一种定量化地描述风险概率以及损失程度的分析方法。此方法认为风险度D可以根据该风险发生的可能性程度L和所造成的损失程度C来确定,其计算公式为:
D=L×C,
(5)
式中:L、C的取值为1~5,其具体取值所对应的含义见表3;D的取值可反映风险的等级,这里采用5级风险标准来划分风险等级[13],具体取值及描述见表4。
表3 L和C取值及含义
Table 3 Value and description of L and C
取值L值含义C值含义1几乎不会发生几乎不引起损失2发生可能性较小引起较小损失3可能会发生引起一般损失4发生可能性较大引起较大损失5极易发生引起巨大损失
表4 风险等级D取值及语言描述
Table 4 Value of risk rank D and description
D取值风险等级风险水平1~31级低风险3~62级较低风险6~103级中等风险10~164级较高风险16~255级高风险
鉴于专家打分法是构建风险判断矩阵及风险度计算的重要环节和依据,在具体评估时应选择具有丰富的外海造岛理论知识和工作经验的专家进行打分,以实际调研资料和现有成果为依托,充分论证各风险因素对工程目标的影响程度。
港珠澳大桥人工岛建设各层因素权重和风险度计算具体步骤如下:
(1)以矩阵A的风险识别结果为依据,运用专家打分法对各2级子工作中的风险因素进行两两比较,构建判断矩阵,并进行一致性检验,进而确定出各风险因素相对于2级子工作的权重。
(2)根据表3中L、C的取值及相应含义,咨询专家并查阅相关资料,对各参数进行赋值,进而计算各2级子风险相对于各2级子工作的风险度。
(3)由各2级子风险的权重及对应的风险度进行加权运算,得出各2级子工作的风险度,公式如下:
(6)
式中:Dij表示对应工作的风险度;λijl表示各风险因素在所属工作下的相对权重;Dijl表示各风险因素在所属工作下的风险度;n为各子工作中所包含的风险因素个数。
(4)根据步骤(3)中各2级子工作的权重及风险度,再次进行加权运算,求出各1级子工作的风险度,以此类推,最终求出该工程的总风险度。
限于文章篇幅,仅以“钢圆筒振沉”(W12)这个2级子工作为例,来说明此方法的应用。由矩阵A可知,W12中所具有的风险因素有R11、R21、R22、R23、R33、R42、R44,对此7项风险因素进行整理并打分,构建判断矩阵。计算该判断矩阵的最大特征值和对应的特征向量,由式(3)和(4)计算出一致性指标为0.051,说明该矩阵一致性满足要求。对该判断矩阵最大特征值所对应的特征向量进行归一化处理,得各2级子风险相对于2级子工作Wij的权重分别为0.050、0.044、0.204、0.434、0.166、0.069、0.030。采用专家打分法并由式(5)计算出各子风险的风险度分别为8、9、3、4、12、6、8,则2级子工作W12的风险度为5.635。以此类推,可求出所有2级子工作的权重,见表5。
表5 各级工作风险度及权重计算结果
Table 5 Calculation result of risk degree and weight in each work
1级子工作1级子工作权重1级子工作风险度2级子工作2级子工作权重2级子工作风险度岸壁构造W10.5396.143陆域填充W20.1637.013配套设施建设W30.2979.344基槽开挖W110.0978.280钢圆筒振沉W120.5435.635挤密沙桩打设W130.2555.809抛石斜护坡建造W140.1037.634地基处理W210.2806.373沙料回填W220.14011.549止水设计W230.3396.210排水设计W240.2396.241挡浪墙建设W310.0546.223护岸结构建设W320.1079.472岛桥结合部建设W330.2139.490岛隧结合部建设W340.50510.160工作场地建设W350.1196.924
由表5可以看出,在岸壁构造工作中,基槽开挖W11的风险度最高;陆域填充工作中,沙料回填W22的风险度最高;配套设施建设工作中,岛隧结合部建设W34的风险度最高。
根据表5中各2级子工作的权重及风险度,通过加权运算,可计算出各1级子工作的风险度,再结合各1级子工作所占权重,最终得出港珠澳大桥人工岛建设项目的总体风险度D=7.236。由评价结果和表5可知,该项目总体风险水平为3级中等风险,属可接受风险,总体上处于可控状态。其中,2级子工作沙料回填W22和岛隧结合部建设W34的风险水平为4级较高风险,需要给予重点关注。
港珠澳大桥人工岛建设工程系统庞大、技术难度高。通过构造WBS-RBS矩阵,将该工程进行项目结构分解和风险分解,并对项目建设过程中各风险因素进行识别和重要性比较,运用AHP方法对各风险因素进行定量计算,确定其权重和各级子工作的风险度,进而计算出该工程总体风险度和风险等级。评价结果和工程调研结果基本一致,可为其他同类工程的风险识别、定量分析和风险防范措施的制定提供借鉴。
此外,工程建设过程中各风险因素的相对重要性是随工程的推进不断变化的,可在风险分析中引入时间因素,使风险评价结果更加科学可靠。
[1] 宋君妍, 许朝旭, 桑大钊. 油田人工岛开发建设风险控制措施分析[J]. 安全、健康和环境, 2013, 13(1):45-49.
[2] 冯民学, 焦雪, 钟颖颖, 等. 江苏洋口港区雷电特性及风险区划[J]. 气象科学, 2012, 32(1):86-93.
[3] 王彦林, 吴泽生. 港珠澳大桥东、西人工岛施工图设计安全风险评估[J]. 施工技术, 2013, 42(11):64-68.
[4] 陈国华, 吴武生, 徐三元, 等. 基于WBS-RBS与AHP的跨海桥梁工程施工HSE风险评价[J]. 中国安全科学学报, 2013, 23(9):51.
[5] 贾俊峰, 梁青槐. WBS-RBS与AHP方法在土建工程施工安全风险评估中的应用[J]. 中国安全科学学报, 2005, 15(7):101.
[6] 王彦林, 闫禹. 港珠澳大桥外海人工岛快速成岛技术[J]. 施工技术, 2012, 41(8):47-51.
[7] 田林钢, 马成功, 王绪. 基于AHP-熵权法的水利工程业主风险模糊评价[J]. 人民黄河, 2017, 39(12):117-122.
[8] LI Z K, WANG T, GE W, et al. Risk analysis of earth-rock dam breach based on dynamic bayesian network[J]. Water, 2019, 11(11): 2305.
[9] 葛巍, 李宗坤, 王文姣, 等. 基于WBS-RBS和AHP的土石坝施工期风险评估[J]. 人民黄河, 2013, 35(6):121-123.
[10] 项贻强, 吴强强, 张婷婷. 基于AHP-FCE模型的桥梁设计风险评估研究[J]. 土木工程学报, 2010,43(S2):275-280.
[11] 李宗坤, 张亚东, 宋浩静, 等. 基于施工进度计划的建设工程施工期风险分析[J]. 水力发电学报, 2015, 34(6):204-212.
[12] 刘铮, 吴小萍, 杨立国. 层次分析法及消去与选择转换法在铁路环评中的应用[J]. 郑州大学学报(工学版), 2007, 28(3):56-60.
[13] 李宗坤, 吴赛, 李定斌, 等. 南水北调中线工程潮河段征迁安置风险分析[J]. 郑州大学学报(工学版), 2016, 37(4):57-61.