DOI： 10.11817/j.issn.1672-7207.2019.07.019
基于蒙特卡洛和BIM的CRTSⅢ型板式无砟轨道铺设工期仿真
王卫东1, 2,雷晓鸣1, 2,杜香刚3,王继军3,贺小刚4
（1. 中南大学 土木工程学院，湖南 长沙，410075；
2. 中南大学 重载铁路工程结构教育部重点实验室，湖南 长沙，410075；
3. 中国铁道科学研究院 铁道建筑研究所，北京，100081；
4. 中国铁路设计集团有限公司，天津，300142）
摘要:以高速铁路1榀32 m标准简支梁上的CRTSⅢ型轨道板铺设为对象，开展铺板工艺过程及工期的三维仿真模拟研究。根据32 m简支梁上CRTSⅢ型板铺设工艺，将铺设过程分解为16道工序，并在京沈(北京—沈阳)高铁施工现场获取每道工序的持续时间样本各100组，拟合确定各道工序持续时间的概率分布。基于双代号时标网络计划图，应用Revit软件，以蒙特卡洛方法开展工序工期和铺板总工期仿真模拟并建立对施工过程进行动态演示的轨道板铺设BIM模型，确定1榀32 m简支梁上CRTSⅢ型板铺设总工期的概率分布。研究结果表明：应用Revit软件建立具有模拟工期信息的轨道板铺设BIM模型，能形象展示工期仿真结果，实现CRTSⅢ型板铺设进度和施工质量档案的可视化管理。
关键词:CRTSⅢ型板；工程进度；网络计划；蒙特卡洛模拟；建筑信息模型
中图分类号:TU723.3 文献标志码:A 文章编号:1672-7207(2019)07-1655-07
Period-simulation of CRTSⅢ slab track based on Monte Carlo simulation and BIM
WANG Weidong1, 2, LEI Xiaoming1, 2, DU Xianggang3, WANG Jijun3, HE Xiaogang4
(1. School of Civil Engineering, Central South University, Changsha 410075, China；
2. Key Laboratory of Engineering Structures of Heavy Haul Railway of Ministry of Education,Central South University, Changsha 410075, China；
3. Railway Engineering Research Institute, China Academy of Railway Sciences, Beijing 100081, China；
4. China Railway Design Corporation, Tianjin 300142, China)
Abstract: Based on CRTSⅢ slab track on a 32 m standard simple-supported girder of high-speed railway bridge, three-dimensional simulation of the paving process and duration was carried out. According to the paving technology of CRTSⅢ slab track on the 32 m simple-supported girder, the paving process was divided into 16 stages. The duration specimen of process was obtained from the construction site of Beijing—Shenyang high-speed railway, and each has 100 groups, which was used to detemine the duration probability distribution of each process. The process and total durations of the track were simulated with Monte Carlo method based on the double codename network plan diagram, thus the total duration probability distribution of CRTSⅢ slab track on a 32 m simple-supported girder was determined. The BIMmodel for track paving which could simulate the duration information was established by Revit software. The results show that it can display the simulation results of duration vividly and achieve the visual management of CRTSⅢ slab track paving process and construction quality archives.
Key words: CRTSⅢ slab track; project schedule; network plan; Monte Carlo simulation; building information modeling (BIM)
铁路信息化是实现铁路现代化的关键。由于我国铁路结构庞大、联系紧密、各项业务较复杂、建设难度较大，因而采取正确的发展策略显得尤为重要[1]。CRTSⅢ型板铺设工程具有点多线长、单位工程复杂和工艺要求高等特点，其铺设过程不确定因素多，各施工队伍的水平参差不齐，施工管理难度较大[2]，因此，在施工准备阶段要制定科学合理的进度计划，合理安排资源配置，最大限度降低成本。建筑信息模型(BIM)技术在建筑规划、勘察、设计、施工和运营维护全过程中的集成应用能够实现工程建设项目全生命周期数据共享和信息化管理，为项目方案优化和科学决策提供依据，促进建筑业提质增效[3]。目前，铁路信息化及BIM技术应用较广泛，大多处在针对具体专业的开发与应用阶段。由于铁路工程的复杂性、三维协同设计理论与方法未成体系等，BIM技术在铁路项目全生命周期中的价值及优势还有待于进一步提升[4-5]。为此，本文作者将CRTSⅢ型板铺设过程进行BIM建模，有利于进行工期模拟和资源调配，有利于铁路行业轨道板铺设过程BIM技术应用。为合理确定CRTSⅢ型板铺设工期，基于双代号时标网络计划图，以BIM和蒙特卡洛技术开展CRTSⅢ型板铺设工程进度仿真研究。根据CRTSⅢ型板式无砟轨道铺设的施工工艺，明确各道工序；在京沈(北京—沈阳)高速铁路辽宁段32 m简支梁轨道板铺设现场调查，获取每道工序持续时间的100组样本数据，拟合确定轨道板铺设过程中每道工序持续时间的概率分布模型及相应的分布参数；编制双代号时标网络计划图，应用蒙特卡洛(Monte-Carlo)模拟方法对每道工序持续时间进行100万次模拟，得到1榀32 m简支梁上CRTSⅢ型板铺设总工期的概率分布；最后，应用Revit软件对基础构建进行BIM建模，通过搭建网站将前述模拟得到的信息附着在具体模型中，形象展示仿真结果。
图1 研究框架
Fig. 1 Research framework
1 CRTSⅢ型板式无砟轨道铺设工艺WBS分解
本文的整体研究框架如图1所示，其中，WBS工作分解是将项目工作按阶段可交付成果分解成较小的更易于管理的组成部分的过程。
CRTSⅢ型板式无砟轨道的施工工序较多，技术要求高，需要做好详细的施工方案，逐层分解施工流程，确保轨道铺设符合相关规定要求。CRTSⅢ型板铺设过程划分为底座施工和轨道板施工两大部分，在详细调研京沈高速铁路辽宁段32 m简支梁上CRTSⅢ型板(双线共12块板)铺设现场后，本文将32 m简支梁的CRTSⅢ型板式无砟道床铺设过程分解为16道工序，施工流程如图2所示。
2 CRTSⅢ型板式无砟轨道铺设网络计划图
2.1 双代号时标网络计划
施工进度计划图的编制方法分为横道图和网络计划图。横道图是一种简单并且应用较广的传统方法，其最主要的特点是简洁明了，但横道图也存在一些缺点，如各个工序之间的逻辑关系难以表达清楚，不能确定进度计划的关键工作和关键路径，难以适应较大的进度计划系统。鉴于CRTSⅢ型板铺设过程的复杂性，本文采用双代号时标网络计划图进行表示。双代号时标网络计划图兼有网络图与横道图的优点，既能简洁、完整地表达工作间的逻辑关系，又能清楚表示计划的进行过程，便于查询和控制。
图2 CRTSⅢ型板式无砟轨道铺设施工流程
Fig. 2 Laying construction processes of CRTSⅢ slab track
2.2 计划评审技术
在网络计划中，对工作间逻辑关系和工作持续时间的研究经历了从确定到不确定2个阶段。早期在以关键路径法(CPM)为代表的网络计划模型中，工作间逻辑关系和工时都是确定的，工时被假定为常量[6-7]，网络中的参数通过简单的加减计算求得，总工期为关键线路上所有工时之和。而在实际工程中，每道工序的持续时间受诸多因素影响，存在着波动范围，整个项目的关键路径可能随之改变，简单地运用CPM进行工期分析难以反映实际情况。计划评审技术(PERT)是工作间逻辑关系确定而工时不确定的一种网络计划技术，假设各道工序的持续时间服从一定数学分布，近似地用三点估计法算出1个期望值作为工作的持续时间，以此反映活动的“不确定性”[8-10]。本文基于PERT展开，利用蒙特卡洛模拟方法，对网络计划进度进行概率分析，以便为项目管理人员在工程项目进度管理控制方面提供帮助。
2.3 CRTSⅢ型板式无砟轨道铺设工程网络计划图
前面明确了32 m简支梁上CRTSⅢ型板铺设过程中的各道工序以及各工序之间的逻辑关系。根据现场调查，得到各道工序的持续时间范围，得出工作明细表(见表1)，再根据工作明细表绘制双代号时标网络计划图(见图3)。网络图中各项工序的持续时间先采用其时间范围的平均值，之后用蒙特卡洛方法进行模拟。从时标网络图上可直观得到各项工序的持续时间及自由时差，便于决策控制。
表1 工作明细表
Table 1 Work schedules
3 蒙特卡洛模拟
蒙特卡洛模型又称为随机抽样模型，是基于大量实验以确定随机变量数值结果的计算方法，其关键是构造概率过程[11-14]。根据现场收集的样本数据，用MATLAB拟合得到各道工序持续时间的概率分布模型及参数，应用蒙特卡洛方法并基于计划评审技术进行模拟，得到32 m简支梁上CRTSⅢ型板铺设总工期的概率分布，用此分布函数作为工程进度模拟的基准。
图3 CRTSⅢ型板式无砟轨道铺设工程双代号时标网络计划
Fig. 3 Network plan of double code time of CRTSⅢ slab track laying project
3.1 各工序持续时间的样本和概率分布
项目工期由每道工序的持续时间、工序之间的逻辑关系确定，每道工序的持续时间主要由该道工序的工程量确定，可以通过BIM建模软件直接导出各构件的工程量。在工程实施过程中，工期受很多因素的影响，导致其控制存在风险，不能唯一确定。本文根据现场调研结果，取得了在一定资源配置下1个工作面上1榀32 m简支梁上CRTSⅢ型板铺设每道工序持续时间的100个样本值。
根据工程进度安排及施工组织模式，在机械及材料配置满足劳动力需求情况下，1个工作面的施工劳动力初始配置情况如下。
1) 底座板施工人员配备：路基面、梁面、填充面清扫凿毛2人，底座板钢筋制安4人，底座板模板加固4人，底座板模板拆除及倒运3人，底座板混凝土施工6人，混凝土收面2人，混凝土养护2人，其他2人。
2) 无砟轨道轨道板粗铺施工人员配备：龙门吊操作员1人，轨道板吊装人员4人，吊车驾驶员1人，平板车司机1人，中间隔离层铺设及弹性垫层安装2人，轨道板钢筋绑扎4人，轨道板粗铺4人。
3) 精调人员配备：安装精调爪及安装压紧装置3人，精调调整人员4人。
4) 自密实混凝土灌注人员配备：工装安、拆、倒运4人，自密实混凝土灌注5人，斗车司机1人，灌车司机1人，龙门吊司机1人，清扫、养护1人。
以MATLAB软件对每道工序的100组样本数据进行分布拟合，确定了在32 m简支梁上CRTSⅢ型板铺设16道工序持续时间概率分布如表2所示。其概率分布函数主要是正态分布、Beta分布、三角分布和均匀分布[15-16]。
3.2 蒙特卡洛模拟
通过拟合得到各道工序持续时间的概率分布后，基于CRTSⅢ型板式无砟轨道铺设工程双代号时标网络计划图(见图3)，进行100万次轨道板铺设的蒙特卡洛仿真模拟，可以得到1榀32 m简支梁上铺设CRTSⅢ型板的各工序持续时间及总工期的100万个结果。总工期的蒙特卡洛仿真模拟结果如下。
1) 工期概率分布见图4。从图4可见：1榀32 m简支梁上铺设CRTSⅢ型板的总工期概率分布更加符合Beta分布，最短工期为28.56 d，最长工期为35.50 d。
2) 总工期的累计概率分布模拟结果见图5。由图5可知：在33 d内完成工程的概率为99.96%，在32 d内完成工程的概率为90%。
3) 总工期的敏感度分析结果见图6，其中，总工期敏感度指各道工序的自身变化对总工期的影响程序。各工序持续时间对总工期不确定性的影响程度不同，其中自密实混凝土层及底座混凝土层的灌注和养护两道工序对总工期的不确定性影响最大，分别达到38.6%和17.9%。在施工过程中，应该重点控制这2道工序，降低其对总工期的影响。
表2 各工序持续时间概率分布类型及参数
Table 2 Probability distribution types and parameters of each process
图4 工期概率分布
Fig. 4 Probability distribution of construction period
图5 工期累计概率分布
Fig. 5 Cumulative probability distribution of construction period
图6 工期敏感度分析
Fig. 6 Period-sensitive analysis of construction
4 BIM仿真系统
应用Revit软件建立具有模拟工期信息的轨道板铺设BIM模型，以计划工期为时间轴，以双代号时标网络计划图为基础数据，对轨道板施工过程进行动态演示，形象说明工期仿真结果，实现CRTSⅢ型板铺设进度可视化管理。
4.1 无砟道床几何模型
对应CRTSⅢ型板式无砟道床施工过程，建立底座、隔离层及弹性垫层、自密实混凝土层、轨道板共4个主要基础构件的参数化族，再将各族组合成 CRTSⅢ型板式无砟轨道的整体模型[17]，见图7。
图7 CRTSⅢ型板BIM模型
Fig. 7 BIM models of CRTSⅢ slab track
4.2 BIM模型及仿真(施工过程模拟)
将前述过程模拟得到的各工序持续时间、铺板总工期等进度信息添加在每个构件中，然后，以工期为时间轴，以双代号时标网络计划图为基础数据，模拟施工过程进度计划，动态演示底座板施工过程，并实现项目信息动态查询。
该BIM模型可用于实时掌握工程进度，通过与计划进度相比较，及时发现问题，指导现场施工进度安排和人、材、机的安排，减少资源闲置与紧张情况，为工程进度决策提供合理依据，解决了传统施工管理过程中可视化程度低、进度管理与信息查询效率低、信息来源单一等诸多问题。动态展示系统前台采用HTML+CSS进行开发，后台采用PHP+MySQL进行设计与通信[18-19]。系统界面如图8所示。
图8 BIM信息系统界面
Fig. 8 Interface of BIM information system
5 结论
1) 通过对现场调查的100组数据的分析，确定了在32 m简支梁上CRTSⅢ型板铺设16道工序持续时间的概率分布。
2) 以蒙特卡洛模拟方法确定了32 m简支梁上CRTSⅢ型板铺设工期服从正态分布，工期为30~33 d。
3) 各道工序持续时间的随机性对总工期影响程度不一，其中自密实混凝土层和底座混凝土层这2道工序的持续时间对总工期的影响最大。
4) 在不同施工水平下，每道工序持续时间的概率分布及其参数可能存在差异。
参考文献：
[1] 朱建生, 裴坤寿. 21世纪铁路信息化的发展与对策研究[J]. 中国铁道科学, 2000, 21(1): 18-24.
ZHU Jiansheng, PEI Kunshou. Development and Counter- measures of railway informatization in the 21st century[J]. China Railway Science, 2000, 21(1): 18-24.
[2] 王小雪. 基于不确定型网络图的铁路桥梁工程造价模拟方法研究与应用[D]. 长沙: 中南大学土木工程学院, 2013: 11-12.
WANG Xiaoxue. Study and Application of simulation method for railway bridge engineering cost based on uncertain network diagram[D]. Changsha: Central South University.School of Civil Engineering, 2013: 11-12.
[3] 国办发(2017)19号. 国务院办公厅关于促进建筑业持续健康发展的意见[R]. 北京: 国务院办公厅, 2017: 1-5.
State Council (2017) No. 19. Opinions of the general office of the state council on promoting the sustainable and healthy development of the construction industry[R]. Beijing: General Office of the State Council, 2017: 1-5.
[4] 段熙宾. 大型铁路工程 BIM 设计的探索及实现[J]. 铁道标准设计, 2016, 60(12): 124-127.
DUAN Xibin. Exploration and realization of BIM design of large-scale railway engineering[J]. Railway Standard Design, 2016, 60(12): 124-127.
[5] 夏宇, 谭衢霖, 蔡小培, 等. 铁路BIM应用三维线路场景构建研究[J]. 铁路计算机应用, 2018, 27(7): 95-98.
XIA Yu, TAN Qulin, CAI Xiaopei, et al. Railway BIM applied to 3D line scene construction[J]. Railway Computer Application, 2018, 27(7): 95-98.
[6] 谢志强, 韩英杰, 齐永红, 等. 基于关键路径和任务复制的多核调度算法[J]. 国防科技大学学报, 2014, 29(1): 172-177.
XIE Zhiqiang, HAN Yingjie, QI Yonghong, et al. A scheduling algorithm for multi-core based on critical path and task duplication[J]. Journal of National University of Defense Technology, 2014, 29(1): 172-177.
[7] KELLEY J E, WALKER M R. Critical path planning and scheduling[C]//Proceedings of the Eastern Joint Computer Conference. Boston, USA, 1959: 160-173.
[8] CLARK C E. The PERT model for the distribution of an activity time[J]. Operations Research, 1962, 10(3): 405-406.
[9] HAHN E D. Mixture densities for project management activity times: a robust approach to PERT[J]. European Journal of Operational Research, 2008, 188(2): 450-459.
[10] MALCOLM D G, ROSEBOOM J H, CLARK C E, et al. Application of technique for a research and development program evaluation[J]. Operations Research, 1959, 7(5): 646- 669.
[11] 郭琦, 何金静, 胡苗. 基于蒙特卡洛模拟的工程项目网络计划进度风险分析[J]. 项目管理技术, 2013, 11(11): 66-70.
GUO Qi, HE Jinjing, HU Miao. Strategic analysis of network planning progress based on Monte Carlo simulation[J]. Project Management Technology, 2013, 11(11): 66-70.
[12] 邵伟. 蒙特卡洛方法及在一些统计模型中的应用[D]. 济南: 山东大学数学学院, 2012: 1-83.
SHAO Wei. Monte Carlo methods and their applications in some statistical model[D]. Jinan: Shangdong University. School of Mathematics, 2012: 1-83.
[13] 王霞, 张本涛, 马庆. 蒙特卡洛法在投资项目风险分析中的应用[J]. 价值工程, 2011, 30(26): 64-65.
WANG Xia, ZHANG Bentao, MA Qing . Analysis of investment risk based on Monte Carlo method[J]. Value Engineering, 2011, 30(26): 64-65.
[14] 罗丽哲, 胡志根, 刘全, 等. 高坝施工全过程导流风险多维评估[J]. 中南大学学报(自然科学版), 2014, 45(3): 840-846.
LUO Lizhe, HU Zhigen, LIU Quan, et al. Multi-dimensional estimation for risk of diversion system in high-dam construction process[J]. Journal of Central South University(Science and Technology), 2014, 45(3): 840-846.
[15] 徐志军, 张博, 王凯, 等. 基于Beta分布的改进的一次可靠度法研究[J]. 土木工程与管理学报, 2014, 31(3): 17-21.
XU Zhijun, ZHANG Bo, WANG Kai, et al. Study of modified first order reliability method on Beta distribution[J]. Journal of Civil Engineering and Management, 2014, 31(3): 17-21.
[16] 刘银萍. 均匀分布的相关性质以及检验[D]. 扬州: 扬州大学统计学院, 2016: 1-35.
LIU Yinping. Uniform distribution of the correlation properties and testing[D]. Yangzhou: Yangzhou University. Institute of Statistics, 2016: 1-35.
[17] 李彦霖, 王卫东. 基于4D的桥梁施工质量管理信息系统研究与实现[J]. 中南林业科技大学学报, 2011, 31(9): 123-127.
LI Yanlin, WANG Weidong. Study and implementation of bridge construction quality management information system based on 4D[J]. Journal of Central South University of Forestry, 2011, 31(9): 123-127.
[18] 王军. 基于HTML5的网络化认知评测系统的设计与开发[D]. 南京: 东南大学软件学院, 2015: 1-65.
WANG Jun. Design and development of networked cognitive evaluation system based on HTML5[D]. Nanjing: Southeast University. School of Software, 2015: 1-65.
[19] 李红艳, 刘蓉, 陆翼山, 等. 基于相关属性的异构数据融合路由策略[J]. 中南大学学报(自然科学版), 2017, 48(4): 996-1005.
LIU Hongyan, LIU Rong, LU Yishan, et al. Correlation attributes based heterogeneous data fusion routing scheme [J]. Journal of Central South University (Science and Technology), 2017, 48(4): 996-1005.
(编辑 陈灿华)
收稿日期： 2018 -09 -19; 修回日期： 2018 -11 -22
基金项目(Foundation item)：教育部人文社会科学研究项目(13YJA630091)(Project(13YJA630091) supported by Humanities and Social Sciences Research Foundation of Ministry of Education)
通信作者：王卫东，博士，教授，从事道路与铁道工程研究；E-mail: 147745@163.com
