常春光,赵 腾
(沈阳建筑大学管理学院,辽宁 沈阳110168)
随着当今建造技术的发展,传统建筑产品复杂繁琐的建筑工艺已经不能满足人们对高品质建筑产品的需求,装配式建筑利用加工好的建筑配件直接装配完工,具有建造速度快、产品质量高、节能环保等优点[1],因此迅速在世界各地得以推广。对于建筑行业来说,安全事故经济损失数额庞大,不仅会消耗大量的资源,同时还会影响装配式建筑产业未来的发展,因此,装配式建筑安全施工对保障建筑安全、引领装配式建筑产业积极健康发展具有重要意义。目前,国内外专家学者对装配式建筑的安全风险评估研究较多,吴凤平等[2]深入分析了混凝土装配式住宅施工危险源及存在的安全问题;黄桂林等[3]利用可拓学理论评价施工过程中的安全风险状况;李高锋等[4]对装配式住宅施工安全影响因素及相应管理应对措施进行了研究;李颖等[5]在施工危险源分析的基础上,对预制装配式建筑施工的安全和质量问题进行了分析和评估;王志强等[6]运用事故树(Fault Tree Analysis,FTA)-集对分析法(Set Pair Analysis,SPA)-灰色聚类法评估装配式建筑施工过程中的危险因子。上述学者虽然采用了不同方法对装配式建筑施工安全管理问题进行了大量研究,但是大多局限于风险的识别以及评估分析,而对于经济层面及安全投入与风险控制问题之间优化及权衡等方面的研究较少,在保证装配式建筑工程项目处于合理安全等级的前提下,建筑企业如何在风险识别的基础上进一步经济合理地将一定数额的安全投入进行最优化分配已成为目前亟待解决的问题。
鉴于此,本研究从事故致因理论出发,对装配式建筑施工安全事故风险进行识别和分析,同时结合Cobb-Douglas函数构造数学模型,将一定合理数额的安全投入通过最优化的分配方式分配于各类别安全投入之中,从而使装配式建筑施工安全事故发生概率降为最低,最后借助工程算例进行分析,证明该模型的适用性,从而为建设企业得到最优事故预防效果以及最优安全决策提供依据。
装配式建筑施工,潜在性风险因素较多,同时,大量的吊装以及高空作业增加了施工过程中的风险因素[7]。根据国家住房和城乡建设部印发的相关文件,并结合装配式建筑施工的特点,梳理出装配式建筑施工过程中5类高频率事故,依次为高处坠落事故、物体打击事故、坍塌事故、起重伤害以及机械伤害。通过选取典型的装配式建筑作为分析样本,并依据装配式建筑施工流程,反演分析得出相关事故的关键影响因素[8],如表1所示。
表1 装配式建筑施工安全关键影响因素
安全投入,是为满足参与施工生产的人员、机械、设施、环境等具备相应的安全保障而作出的投入,具体包括预防性安全投入和损失性安全投入两种,且两者之间为负相关[9-10]。依据国家《高危行业企业安全生产费用财务管理暂行办法》,预防性安全投入可分为日常安全管理、劳动用品、安全技术措施、建筑施工卫生措施、宣传教育五大类型投入。损失性安全投入即事故损失,由建筑安全事故的发生而导致建筑企业所付出的代价。
“ALARP”(As Low As Reasonably Practicable)原则,即“最低合理可行”原则。在建筑系统风险评价过程中,首先要具备科学合理的风险参照系,ALARP原则可提供项目风险的判定标准,由此依据该原则确定项目建设的风险可接受水平,已被国外建筑企业普遍采用[11-12]。该原则通过不可容忍线和可忽略线将项目风险划分为风险严重区、ALARP区和可忽略区,如图1所示。
图1 风险评价的ALARP原则
(1)如果评估所得风险度落入风险严重区域,此时该装配式建筑项目的风险是不可以接受的,必须加大安全投入,健全管控体系以降低风险。
(2)如果评估所得风险度落入风险ALARP区域,此时该装配式建筑项目的风险在可接受范围之内,需要对该项目进行进一步的成本风险分析,尽可能使项目风险度向可忽略区域靠近[13]。
(3)如果评估所得风险度落入风险可忽略区域,此时该装配式建筑项目的风险在可接受范围内,从经济角度出发,无需加大安全投入。
事故树分析(Fault Tree Analysis,FTA)是一种逻辑演绎的系统评价方法,主要通过逻辑推理、正反演绎等方式,来辨识和评价系统的危险性[14-15]。该方法可同时进行定性与定量分析,在用FTA法分析问题时将顶上事件定义为T,基本事件定义为X,在具体运用时分析及解决问题的程序如图2所示。
图2 FTA法分析及解决问题的程序
定量分析即依据逻辑演绎所得基本事件的发生概率,通过计算转化得出顶上事件概率值[16]。定量计算时常采用最小割集法,过程如下:
设某事故树有K个最小割集:E1,E2,…,Er,…,Ek,则有
(1)
顶上事件发生的概率为
(2)
化简得出顶上事件概率为
(3)
式(3)中:r、s、k为最小割集的序号,r三、装配式建筑安全投入优化模型的建立
1.事故损失模型的建立
柯布-道格拉斯(Cobb-Douglas)生产函数是在一定技术条件下投入与产出关系的数学模型。它是通过建立指数关系衡量系统生产的经济模型,在经济学中使用最为广泛。依据参考文献[17],将C-D生产函数引入建筑业,测算损失性安全投入(事故损失),并建立计量模型
(4)
式中:Y为事故损失;ce为安全教育类投入;ct为安全技术措施类投入;ch为建筑施工卫生措施类投入;cl为劳保用品类投入;cr为日常安全管理类投入;ae,at,ah,al,ar为各类安全投入弹性系数;A为施工技术发展系数。
根据装配式建筑特点及合理决策要求,并结合我国建筑企业实际情况,该函数应满足以下条件:
①Y≥0,针对装配式建筑而言,事故损失额不可能一直为0。
按照事故树定量计算的相关定义,结合式(3)Pi为各基本事件Xi发生的概率,可得装配式建筑施工安全事故树顶上事件T发生的概率为
(5)
因为装配式建筑安全施工事故树中各基本事件的发生概率与安全投入呈负相关对应,且随着安全投入的增加,顶上事件发生概率的降低率随之减小,且当二者之一趋于0时,另一项为无穷大。这种关系特点与生产函数曲线“龚伯兹曲线”(当0<αi<1,bi>1时)的趋势相一致。因此,可借助龚伯兹曲线描述基本事件发生概率与安全投入之间的函数关系,如图3所示。
图3 龚伯兹曲线走势
利用龚伯兹曲线拟合基本事件发生概率Pi与其对应的各项关键风险因素预防所需的具体安全投入ci之间的关系,可得
(6)
式中:i=1,2,…,n,Ki,ai,bi均为参数,且有0<αi<1,bi>1,参数值依赖以往装配式建筑安全事故统计数据,利用曲线拟合加以确定。
将式(6)代入式(5),可得
(7)
借助目标规划的思想,寻求最优化的安全投入分配方式,从而使项目风险度最低,即事故树顶上事件发生概率最小,所以该模型目标函数为minP。
限制条件推导:
为保证装配式建筑项目施工安全,依照ALARP原则,预防性安全投入总额需限定在安全区间[M,N]内,故根据此约束条件有
(8)
从建筑企业自身利益出发,事故损失需尽可能小,即对于损失性安全投入需设定上限值,此处假设上限值为R,则有
(9)
此外,根据《施工企业安全生产管理规范GB50656—2011》,并结合我国装配式建筑安全事故发生特点,要求企业在实际决策过程中分配额需有侧重,即各类别安全总投入之间存在一个排序,记为
ce≤ch≤cr≤cl≤ct
(10)
结合实际情况,建筑施工安全技术措施经费由地方政府规定,即相应的安全技术措施经费需设定取值下限,故有
ct≥T,cl≥L,cr≥S
(11)
式中:T,L,S为各地方政府规定的相应安全技术措施费取值下限。
综上所述,以minP为目标函数,建立非线性规划数学模型
在实际运行中,借助运筹学中常用计算最优解的Lingo软件求解模型,得各项预防基本事件发生的安全投入ci的值,此时可以获得事故树顶上事件的最小发生概率,同时将各ci值代入式(6),以获取每个基本事件的发生概率。然后按类别将所求得的安全投入ci分别归入5类安全投入中,即可求得最优安全投入分配方案。
北京市某建筑企业计划完成6 083万元的装配式住宅建设项目,该企业近10年各项安全投入和事故损失Y的数据如表2所示,利用该模型对该项目的安全投入结构进行优化。
表2 该企业安全投入数据 万元
根据式(4),将各项安全投入与事故损失取对数,即该函数关系可变式为
lnY=lnA+lncl+lnct+lncr+lnch+lnce
(12)
将该建筑企业近10年安全投入与事故损失取对数后,通过数据统计分析软件SPSS24.0进行多元线性回归可得
Y=49.367ce-0.062ct-0.624ch-0.173cl-0.091cr-0.047
(13)
依据装配式建筑施工安全事故关键风险因素的划分,各类别安全投入额与预防基本事件发生的安全投入ci之间有以下关系
将此关系代入式(13),该函数模型可以转化为事故损失额与各项预防基本事件发生的安全投入ci之间的函数关系
(14)
经检验,上述模型同时满足3项约束条件,即该事故损失模型成立。
根据相关文献及数据记载,同时结合表1,运用事故树分析法(FTA)建立装配式建筑施工安全事故树,如图4所示。
注:X1为风险意识淡薄;X2为安全知识欠缺;X3为偶然性失误;X4为脚手架搭设不合格;X5为机械设备的检修与维护不到位;X6为构件出厂质量不达标;X7为吊装设备操作失误或机具故障;X8为装配区域施工杂乱;X9为天气、照明情况恶劣;X10为安全防护设施已损坏;X11为防火、防电、防水措施不严格;X12为未配备有效防护措施;X13为现场安全监理风险;X14为日常安全检查疏漏;X15为管理组织风险。图4 装配式建筑施工安全事故树
根据ALARP原则,预先确定合理可行的项目安全投入区间,由于我国装配式建筑事故统计体系尚未完善,缺乏相关方面文献及数据记录资料,可行性区域的标准安全投入区间[M,N]难于确定,参照发达国家风险度划分标准,并结合本国国情和本企业实际情况,此处作出合理假设:其预防性安全投入的下限取近10年的最小值,即M=153.53,上限取近10年的最大值,即N=230.41。同时,根据海因里希法则(Heinrich′s Law)[18],安全事故损失总额是事故直接损失的5倍,由此该建筑企业在过去10年中的事故损失最小值为35.80,设想通过该模型对安全投入进行优化后,事故损失上限为35.80万元,即R=35.80,满足风险度落在可行性区域。
依据《北京市人民政府贯彻落实国务院关于进一步加强安全生产工作决定的若干意见》,要求以建筑安装工程造价为计取依据,房屋建筑工程安全生产费用提取标准为2%,其中,安全技术措施保障资金为40%,劳动防护用品保障资金为30%,日常安全管理专项资金为30%。由此有
该装配式建筑安全投入优化分配模型为
根据模型目标函数,该装配式建筑施工项目安全事故发生的最小概率为4.77%,计算如下:
在使装配式建筑施工安全事故发生概率最低这个目标函数下求得15项的基础事件的安全投入,并将上述15项预防基本事件发生的安全投入归到安全投入五大类别中,得出最终计算结果,即安全教育投入11.64万元、安全技术措施投入66.87万元、工业卫生措施投入2.49万元、劳动保护用品投入43.60万元以及日常安全管理投入38.93万元。
表3 安全投入计算结果及类别划分
将该装配式建筑施工安全投入优化模型计算出的最优安全投入与企业原安全投入相对比(见表4),差额率在10%左右为微小变动,在实际工程安全投入过程中虽可略作调整,但工业卫生措施投入74.13%、日常安全管理投入57.17%以及安全教育投入58.50%,差额率较大。从计算结果看,安全教育投入、日常安全管理投入及工业卫生措施投入在项目建设的风险控制方面发挥了不可或缺的作用,但企业对其重要性认识不足,导致企业对于此方面的投入较少,远远达不到实际需求,说明管理者对于安全投入的决策缺乏科学性和高效性,需重新调整分配方式,以达到降低事故发生概率的最佳效果。
表4 最佳安全投入与企业原投入计算比较
(1)在应用本模型进行实例分析时,需根据ALARP原则确定能够使风险度落在ALARP区域内的合理安全投入区间,但由于我国建筑行业对此尚未提供确切标准,所以作合理假设,缺少确定合理安全总投入区间[M,N]的推理和判断。
(2)FTA法在构建事故树时,基本事件的推理分析依赖于进行龚伯兹曲线拟合的数据的多少,由于我国建筑数据统计系统尚不完善,所以本研究中的事故树仍不够全面,精细度有待提高,基本因素事件还有待深入推理和挖掘。
(3)模型利用FTA法所分析的基本事件多是直接因素,在保障项目安全建设过程中,还存在安全科研投入和安全保险投入等,但通过FTA法不易获得。