1、2019中国城市轨道交通中国城市轨道交通BIM技术发展大会技术发展大会 基于基于BIMBIM的轨道交通地质的轨道交通地质 环境模型及风险评估环境模型及风险评估 1 汇报人： 武志毅 导 师： 高 燕 副教授 中山大学地球科学与工程学 院 2019.11.29 中山大学地球科学与工程学院中山大学地球科学与工程学院 现有教师及专职研究人员百余人，其中，中科院院士2名，教 育部长江特聘教授、国家杰出青年基金、国家高层次人才特殊 支持计划（“万人计划”）入选者、国家中组部青年千人计划 人才、珠江学者特聘教授、珠江科技新星创新人才等15人次， 中山大学百人计划引进人才近20人。拥有广东省重点实验室2 个2、工程中心1个。 拥有博士后流动站、理学工学博士点、硕士点。主要研究方向 大地构造、地球化学、地球物理、资源矿产、海洋地质灾害与 重大工程建设、地下工程、边坡、泥石流等地质灾害。 2 张培震院士张培震院士高锐院士高锐院士 王岳军院长王岳军院长 学院创建于1924年，是全 国建系较早、历史悠久的 综合性大学地质学系之一。 目录 研究意义及技术路线 基于BIM轨道交通地质环境载体 模型 动态模糊可变评价模型 总结 3 研究研究意义及技术路线意义及技术路线 4 研究意义及技术路线 研究意义 技术路线 轨道交通发展现状轨道交通发展现状 3883.6 240.8 98.5 502 246.1 57.9 3、3.9 0 500 1000 1500 2000 2500 3000 3500 4000 4500 地铁地铁 轻轨轻轨 单轨单轨 市域快轨市域快轨 现代有轨电车现代有轨电车 磁浮交通磁浮交通 APM 里程数（单位：里程数（单位：km） 77.20 4.80 2.00 10 4.901.100.10 地铁地铁轻轨轻轨 单轨单轨市域快轨市域快轨 现代有轨电车现代有轨电车磁浮交通磁浮交通 APM 轨道轨道交通各制式的里程数交通各制式的里程数条形图条形图轨道轨道交通各制式占比图交通各制式占比图 研究意义 5 轨道交通发展现状轨道交通发展现状 835 999 1430 1713 2064 2746 314、73 3618 4153 5033 19.64 43.14 19.79 20.49 33.04 15.55 14.02 14.79 21.19 0 5 10 15 20 25 30 35 40 45 0 1000 2000 3000 4000 5000 6000 2008200920102011201220132014201520162017 增长率（增长率（%） 轨道交通运营里程（轨道交通运营里程（km） 年份（年）年份（年） 轨道交通运营里程（轨道交通运营里程（km）增长率（增长率（%） 2008-2017年我国城市轨道交通运营里程发展趋势图年我国城市轨道交通运营里程发展趋势图 6 研究意5、义 轨道交通面临地质风险轨道交通面临地质风险 西安地裂缝分布图（杨新安 2015） 广州市九号线沿线区域地质 构造图（乐璐 2013） 7 研究意义 广州市十四号线一期工程：岩溶、基岩风化不均及风化深 槽、孤石、软土、风化岩及残积土；广州市轨道交通八号 线北延段：断裂、岩溶、基岩风化不均及风化深槽、软土。 轨道交通面临地质风险轨道交通面临地质风险 广州市十四号线一期工程基 岩地质图（杨珊 2016） 广州市八号线北延段线路基 岩地质图（乐璐 2013） 8 研究意义 由于轨道交通的快速发展，地质环境的复杂性，地质灾害 事故频发。 南京 2007北京 2007广州 2013佛山 2018 经过事6、故缘由的分析：发现事故多与工程所处的地质环境 有关。 轨道交通地质灾害事故轨道交通地质灾害事故 渗水坍塌地质结构复杂风化深槽工地突发透水 9 研究意义 轨道交通隧道工程不良地质风险评价研究现状 年份年份作者作者研究对象研究对象研究成果研究成果 2011李明 隧道和地下工 程 基于隧道和地下工程运营阶段的结构变异、衬砌开裂等因素建立了 健康动态诊断系统，实现了隧道和地下工程评价模型的可变性 2012张萌地铁 针对地铁建设期间的土性、岩性等10个风险因素，建立了基于模糊 综合评价法的风险评价模型 2016王鸣涛地下工程 发现针对上软下硬的复合地层的围岩等级划分的欠缺，指出了对于 未来围岩等级划分的7、方向 2017赵文娟地下工程建立了典型上软下硬地层隧道围岩稳定性的定性定量化标准 2018 Shi S. S. Xie X. K. Bu L. 山岭隧道 为评估喀斯特地区隧道施工中突水灾害，建立了包括地下水位、岩 溶层厚度、平均温度、地下水循环特征、断层构造等11个指标的突 水评价模型 2018 Yang X. L. Zhang S. 公路隧道 为有效评价和控制突水的发生，建立了包括隧道压头、围岩渗透系 数、围岩破碎程度、节理面相对角度和隧道截面尺寸5个指标的隧 道突水风险评估模型 2018 Huang X. Li S. C. Xu Z. H. 公路隧道 权重分析表明，围岩等级，支护测量，地层8、岩性，不利地质，施工 方法，爆破技术，先进地质预报和地下水位是主要的控制因素 10 研究意义 整体技术路线 11 动态模糊可变评价动态模糊可变评价 模型模型 施工工施工工法选择法选择 系统系统 基于基于BIM轨道交通轨道交通 地质环境载体模型地质环境载体模型 风险等级风险等级 低低 成本优化分析成本优化分析 模型模型 最最 优优 决决 策策 法法 多目标优化理论多目标优化理论 设计方案合理设计方案合理 地质环境安全地质环境安全 建议修改设计方案，给出处建议修改设计方案，给出处 理不良地质方案成本对比报理不良地质方案成本对比报 告告 地质环境危险，给出风险预地质环境危险，给出风险预 警报告警报告9、 否否 是是 风险分析评价风险分析评价 基于基于BIM轨道交通地质轨道交通地质 环境风险决策模型环境风险决策模型 地质数据导出地质数据导出 基于基于BIMBIM轨道交通地质环境轨道交通地质环境 载体模型载体模型 研究现状 技术路线 地质环境载体模型搭建 基于BIM轨道交通地质环境载体模型 12 基于基于BIM的轨道交通地质载体模型的轨道交通地质载体模型 年份年份作者作者研究成果研究成果模型模型 2014 中国电 建集团 在AU大师会上，展示了他们利用Civil 3D软件二次开发和 其他BIM软件建立的水电站坝体全生命周期的BIM模型 2014林孝城 采用建筑设计专业广泛使用的 Revit Ar10、chitecture 核心 建模软件，对“金山勘察综合试验厂房项目”岩土工程勘 察成果进行三维可视化的初步实践 2016刘 续 在Revit软件中建立了大量的岩土族构建，如：锚索、立柱 桩、挂网喷砼等，以武汉亚洲医院基坑项目为依托，进行 了BIM的完整实施，探讨了模型建立方式与施工应用方法 2017张国锋 以某工程的地质环境三维建模为例介绍了数据采集和预处 理、表面模型构建、实体模型的构建的基本流程 2017饶嘉谊 通过 BIM 技术建立三维地质模型，模拟土层与岩层的分布， 同时根据已有设计条件建立桩基础模型。二者结合，实现 基于 BIM 的三维地质模型与桩长校核应用 研究现状 13 Revi11、t elements Model element s View- specific elements Datum element s Hosts Model components Walls Floors Roofs Ceiling s Stairs Window s Doors Furnit ure Grids Levels Refrenc e plances Annotatio n elements Details Text notes Tags Symbols Dimension s Detail lines Filled regions 2D detail component s 建模12、工作量大且繁琐，效率 较低； 对存在透镜体、尖灭的岩土 地层建模较为困难；断层、 褶皱等复杂地质构造； 缺乏与岩土工程勘察领域相 关的计算与模型输出接口。 解决方法解决方法 建立适合岩土工程勘察的岩土层族库； 利用Revit API二次开发解决与岩土工程勘察成果的接口问题。（Revit API二 次开发的开发工具：Microsoft Visual Studio（简称VS）、Microsoft Visual Studio Tools for Applications，简称VSTA） Revit Revit 数据库对象分数据库对象分 类类 Revit Revit 的不足的不足 研究现状 14 在V13、isual Studio2015中新建类库工程 添加对Revit API接口装配文件及System.Windows.Forms 的引用 添加Revit API的命名 空间引用 添加Windows窗体控件的命 名空间引用 在程序为命令类添加控制命令的事务模式和更新模式 新建类，从IExternalCommand接口派生，重建Execute() 方法，进行编译 是否编译 成功 是 生成.dll文件，通过Add-In Manager加载装配 文件 运行程序查看结果 功能是否满 足要求 是 程序开发结束 否 否 Revit二次开发流程（高原二次开发流程（高原 2017） 建立表面基准三角网建立表面基准14、三角网 Kriging插值Delaunay三角化算法 生成三维地层生成三维地层 三棱柱自适应族自适应族节点匹配 地质环境载体模型地质环境载体模型搭建搭建技术路线技术路线 钻孔数据预处理钻孔数据预处理 钻孔数据收集钻孔族建模区域划分 断层面生成断层面生成 断层面为平面 最小二乘法拟合 断层面为曲面 Kriging插值+Delaunay 三角化算法 技术路线 15 全国重要地质钻孔数据库服务平台（全国重要地质钻孔数据库服务平台（1： 300） 序列序列钻孔编号钻孔编号XY真实真实X真实真实Y 1DHS04380.350.3587.587.5 2DS317913.25250812.5 3DS319815、1.14.852751212.5 4DHS04390.056.2512.51562.5 5DS31721.28.73002175 6DS31962.28.55502125 7DS31873.48.98502225 8DS31683.410.158502537.5 9DS32153.80.295050 10DS32295.650.751412.5187.5 11DS32325.52.51375625 12DS31935.952.51487.5625 13DS31866.85.1517001287.5 14DS32074.96.112251525 15DS31734.67.511501875 1616、DS31636.558.751637.52187.5 17DS31568.710.221752550 18DS31759.659.852412.52462.5 19DS31391111.1527502787.5 20DS316412.810.532002625 21DHS043413.210.4533002612.5 22DS318813.37.233251800 23DS317614.3735751750 24DS320511.654.952912.51237.5 25DS323011.054.12762.51025 26DS322810.7426751000 27DS320013.83.717、53450937.5 28DS321812.17.530251875 29DS325413.70.53425125 30DS320687.220001800 钻孔钻孔数据收集数据收集 钻孔数据预处理 16 现编号现编号DS3197原编号原编号图幅号图幅号f49g022085（海珠广场）（海珠广场） X坐标坐标Y坐标坐标孔口标高孔口标高孔深（孔深（m）稳定水位（稳定水位（m）钻孔类型钻孔类型开孔日期开孔日期终孔日期终孔日期 3.2437.2 资料来源资料来源工程名称工程名称中大祈乐苑教师住宅楼中大祈乐苑教师住宅楼 地层时代地层时代 地层代号地层代号 层号层号层底标高（层底标高（m） 层底深度层底18、深度 （m） 层厚（层厚（m）柱状图柱状图 取样与分析取样与分析 描述描述备注备注 组组段层段层深度（深度（m）测试测试 QPQP110.24 3.00 3.00 杂填土：杂色，由砖、砼碎块及杂物等组 成，松散 2-0.36 3.60 0.60 淤泥：灰黑色，含少量粉细砂及腐木，具 腐臭味，流塑； E2E2h 3-0.96 4.20 0.60 细砂：浅灰白色，饱和，颗粒较均匀，松 散； 4-3.26 6.50 2.30 粉质黏土：暗紫红色，含少量中细砂，硬 塑，为原岩风化残积土； K1K1bh 5-10.56 13.80 7.30 粉砂岩：暗紫红色，强风化，岩芯呈半岩 半土状，残留岩块，手可折19、断，遇水软化； 6-18.16 21.40 7.60 粉砂岩：暗紫红色，裂隙发育，中等风化， 岩芯较完整，呈块状、短柱状、柱状； 7-20.16 23.40 2.00 粉砂岩：暗紫红色，微风化，岩芯完整， 呈柱状、长柱状； 8-24.86 28.10 4.70 粉砂岩：暗紫红色，裂隙极发育，中等风 化，岩芯较完整，呈块状、短柱状、柱状； 9-26.16 29.40 1.30 粉砂岩：暗紫红色，强风化，岩芯呈半岩 半土状，残留岩块，手可折断，遇水软化； 10-27.36 30.60 1.20 粉砂岩：暗紫红色，裂隙较发育，中等风 化；岩芯较完整，呈块状、短柱状、柱状； 11-33.96 37.220、0 6.60 粉砂岩：暗紫红色，微风化，岩芯完整， 呈柱状、长柱状； 广州城市钻孔柱状图（以钻孔广州城市钻孔柱状图（以钻孔2 2为例）为例） 钻孔数据预处理 17 地层参数设置界面地层参数设置界面 钻孔族模型示意图（钻孔族模型示意图（1：200） 钻孔钻孔族族建模建模 区域区域划分划分 区域划分区域划分 钻孔数据预处理 18 选取随机点数据选取随机点数据 计算均值点，拟合计算均值点，拟合 拟合模型选取拟合模型选取 得出偏基台值得出偏基台值 Math.NET 输入原始输入原始数据数据 指数模型指数模型 输出插值后数据输出插值后数据 KrgingKrging插值法内插插值法内插 Krging插值法21、内插 19 KrgingKrging插值法插值法内插程序实现内插程序实现 Krging插值法内插 20 Delaunay三角形插值法流程图三角形插值法流程图 DelaunayDelaunay三角插值法生成地层三角插值法生成地层 Delaunay三角形插值三角形插值法法程序实现程序实现 Delaunay三角插值 法 21 地层生成地层生成 地层模型（地层模型（1：200） 地层生成 22 构造空间直线构造空间直线确定确定断层面方程的拟合点断层面方程的拟合点 断层断层面生成面生成 断层面生成 23 断层面两边区域的地层延伸到断层面断层面两边区域的地层延伸到断层面 平面：平面：基于最基于最 小二乘法22、拟合小二乘法拟合 曲曲面：面：KrigingKriging插值插值 法法+Delaunay+Delaunay插值法插值法 动态模糊可变评价模型动态模糊可变评价模型 建模技术路线 动态模糊可变评价模型的搭建 动态模糊可变评价模型 24 轨道交通轨道交通 不良地质不良地质 风险因素风险因素 集集 权权 重重 确确 定定 技术路线 25 轨道轨道交通不良地质风险辨识图交通不良地质风险辨识图 轨道交通不良地质风险源轨道交通不良地质风险源 围围 岩岩 等等 级级 复复 合合 地地 层层 地地 下下 水水 有有 害害 气气 体体 特特 殊殊 土土 地地 震震 断断 层层 褶褶 皱皱 地地 裂裂 缝缝 岩岩23、 爆爆 岩岩 溶溶 岩岩 体体 基基 本本 质质 量量 指指 标标 地地 层层 黏黏 聚聚 力力 地地 层层 内内 摩摩 擦擦 角角 隧隧 道道 直直 径径 隧隧 道道 埋埋 深深 渗渗 透透 系系 数数 P H 水水 位位 水水 压压 含含 量量 压压 力力 储储 气气 层层 厚厚 度度 湿湿 陷陷 量量 湿湿 陷陷 性性 土土 厚厚 度度 自自 由由 膨膨 胀胀 率率 液液 化化 指指 数数 溶溶 陷陷 系系 数数 总总 溶溶 陷陷 量量 盐盐 胀胀 系系 数数 硫硫 酸酸 钠钠 含含 量量 总总 盐盐 胀胀 量量 土土 中中 总总 盐盐 量量 设设 防防 烈烈 度度 偏偏 压压 角角 度24、度 断断 层层 面面 倾倾 角角 隧隧 道道 距距 断断 层层 距距 离离 隧隧 道道 与与 褶褶 皱皱 相相 对对 位位 置置 地地 裂裂 缝缝 发发 育育 程程 度度 地地 裂裂 缝缝 发发 展展 速速 率率 地地 裂裂 缝缝 埋埋 深深 隧隧 道道 距距 地地 裂裂 缝缝 距距 离离 相相 交交 角角 度度 原原 始始 地地 应应 力力 所所 处处 位位 置置 溶溶 洞洞 距距 隧隧 道道 地地 面面 压压 力力 水水 头头 最最 小小 净净 距距 溶溶 洞洞 长长 轴轴 岩岩 层层 倾倾 角角 岩岩 性性 （ 碳碳 酸酸 钙钙 含含 量）量） 隔隔 水水 板板 岩岩 体体 厚厚 度度 25、轨道交通不良地质风险因素集 26 轨道交通不良地质风险因素集（以断层，岩溶为例）轨道交通不良地质风险因素集（以断层，岩溶为例） 风险等级风险等级（微风险）（微风险）（弱风险）（弱风险） （中风险）（中风险） （强风险）（强风险）参考文献参考文献 断 层 偏压角度()0101020203030王鲁南， 2012； 李文华， 2013 断层面倾角()7590607530602122/3102/3 岩 溶 所处位置 两者中心的连线和垂 向夹角45 两者中心 的连线和 垂向夹角 45 两者中心的 连线和水平 向夹角 45 薛亚东， 2017； 李继伟， 2010 溶洞距隧道地面压 力水头(m) 60 26、最小净距(S/D)10.7510.50.7500.5 溶洞长轴（m）07710101414 岩层倾角()0108090 1025或 6580 2565 岩性(CaCO3含量%)05550507575 隔水板岩体厚度(m) 00.5或 d30 10305100.55 轨道交通不良地质风险因素集 27 权重确定权重确定 权重确定方法权重确定方法优点优点缺点缺点 专家咨询法 不至于出现属性权重与属 性实际重要程度相悖的情 况 较大的主观经验性,专家的 选择将直接影响到结果的准 确性 层次分析法系统性，实用性，简洁性 其标度问题一直未得到良好 的解决 主成分分析法 以实测资料为基础，比较 客观 缺乏对27、不同诊断指标在物理 力学关系上的重要性的考虑 熵值法 客观性 ：相对那些主观赋 值法精度较高；适应性 缺乏各指标之间的横向比较， 权重依赖于样本 人工神经网络 可以学习和自适应不确定 的系统，能同时处理定量 和定性知识 需要较多的实测数据 权重确定方法优劣对比表权重确定方法优劣对比表 权重确定 28 BPBP神经网络赋权法神经网络赋权法 BP网络结构图网络结构图 读取训练数据 对数据归一化 构造神经网络 设置训练参数 对检测样品仿真 误差是否 在允许范 围内 是 否 matlab计算权重步骤计算权重步骤 权重确定 29 乘积标度赋权法乘积标度赋权法 对n个指标的重要性进行排序； 对指标的进行两28、两比较 对该指标层的所有指标进行归一化处理，得到指标层的 总权重。 重要性重要性标度值（标度值（A:B）权重值（权重值（A:B） 指标A、B“相同”1:10.5:0.5 指标A比指标B“稍微大”1.354:10.575:0.425 指标A重要性是指标B重要性的n 个“稍微大” 几种标度法对“稍微大”的标 度值分别为3、1.286、1.5和 1.277，一般两者不应大于1.5， 则(1.5+1.277+1.286)/3=1.354 权重确定 30 BP神经网络赋权 法所得的权重 复合复合地层的组合地层的组合权重权重： 乘积标度赋权所 得的权重 以复合地层为例 =0.63 =0.37 组合赋权组合29、赋权 权重确定 31 轨道交通轨道交通 不良地质不良地质 风险因素风险因素 集集 权权 重重 确确 定定 动态模糊可变评价模型技术路线 32 样本特征值函数向量样本特征值函数向量 针对围岩等级的辨识因子为BQ， 式中Rc为岩石单轴饱和抗压强度，其随时间变化规律（以花岗 石为例）为： 复合地层随时间变化的辨识因子包括黏聚力和內摩擦角（以高 庙子膨润土为例），其随时间的变化规律为： C=907.7exp(- 0.03t)+664.4exp(0.005t) =4.4exp(-0.135t)+11.2exp(- 0.001t) 动态模糊可变评价模型 33 二级指标隶属度矩阵 相对差异函数 范围值区间c30、，d吸引域区间a，b确定M位值 acMdb 二级指标层隶属度矩阵求取二级指标层隶属度矩阵求取 动态模糊可变评价模型 34 一级指标层、目标层隶属度矩阵、级别特征化一级指标层、目标层隶属度矩阵、级别特征化 级别级别特征值化：特征值化： 风险等级风险等级(微风险微风险)(弱风险弱风险)(中风险中风险)(强风险强风险) 特征值H0.8510.60.850.350.600.35 轨道交通隧道工程不良地质风险等级表轨道交通隧道工程不良地质风险等级表 动态模糊可变评价模型 35 区段区段参数取值参数取值H风险等级风险等级评价等级评价等级 1 =1，p=1（0.172 0.251 0.309 0.268）031、.547 =1，p=2（0.164 0.242 0.311 0.283）0.537 =2，p=1（0.161 0.364 0.220 0.255）0.574 =2，p=2（0.151 0.337 0.246 0.266）0.559 2 =1，p=1（0.175 0.291 0.376 0.158）0.592 =1，p=2（0.167 0.275 0.380 0.178）0.578 =2，p=1（0.171 0.344 0.401 0.084）0.626 =2，p=2（0.159 0.317 0.421 0.103）0.607 3 =1，p=1（0.175 0.291 0.332 0.202）032、.579 =1，p=2（0.170 0.278 0.334 0.218）0.568 =2，p=1（0.152 0.368 0.316 0.164）0.598 =2，p=2（0.162 0.354 0.309 0.175）0.596 4 =1，p=1（0.230 0.301 0.327 0.142）0.628 =1，p=2（0.178 0.286 0.357 0.179）0.586 =2，p=1（0.219 0.376 0.309 0.096）0.655 =2，p=2（0.179 0.323 0.384 0.114）0.616 广州某拟建广州某拟建轨道轨道交通各区段地质风险等级评价表交通各区段地33、质风险等级评价表 模型检验模型检验 模型检验 36 区段区段一级指标一级指标风险等级风险等级 1 围岩等级（0.000 0.000 0.429 0.571）0.304 复合地层（0.360 0.374 0.155 0.111）0.724 地下水（0.225 0.513 0.238 0.025）0.714 地震（0.000 0.000 0.500 0.500）0.325 2 围岩等级（0.000 0.167 0.667 0.166）0.467 复合地层（0.221 0.497 0.144 0.138）0.674 地下水（0.375 0.326 0.268 0.031）0.744 地震（0.00034、 0.000 0.500 0.500）0.325 3 围岩等级（0.000 0.000 0.500 0.500）0.325 复合地层（0.402 0.353 0.131 0.114）0.740 地下水（0.213 0.446 0.292 0.048）0.684 地震（0.000 0.000 0.500 0.500）0.325 断层（0.139 0.631 0.230 0.000）0.706 岩溶（0.179 0.303 0.444 0.074）0.621 4 围岩等级（0.000 0.167 0.667 0.166）0.467 复合地层（0.386 0.361 0.120 0.133）0.7235、8 地下水（0.225 0.475 0.269 0.031）0.703 地震（0.000 0.000 0.500 0.500）0.325 断层（0.052 0.514 0.415 0.019）0.625 岩溶（0.369 0.264 0.121 0.246）0.660 广州某拟建广州某拟建轨道交通轨道交通各各区段区段一级指标风险一级指标风险等级等级评价评价表表 模型检验 37 区区 段段 一级一级 指标指标 二级指标二级指标 风险风险 等级等级 4岩溶 溶洞所处位置（0.000 0.000 0.000 1.000） 0.175 溶洞距隧道地面压 力水头(m) （0.787 0.213 0.0036、0 0.000） 0.941 最小净距(S/D)（1.000 0.000 0.000 0.000）1 溶洞长轴（m）（0.250 0.625 0.125 0.000） 0.763 岩层倾角()（0.000 0.000 0.440 0.560） 0.307 岩性(CaCO3含量%) （0.214 0.643 0.143 0.000） 0.748 隔水板岩体厚度(m) （0.065 0.565 0.370 0.000） 0.650 广州某拟建广州某拟建轨道轨道交通区段交通区段4 4部分部分二级指标二级指标风险风险等级等级评价评价表表 模型检验 38 总结总结 总结 39 总结 40 u基于Revit软件进行二次开发，有效地解决了Revit软件中没 有三维地质环境载体建模模块的问题，解决了Revit软件中 断层平面、曲面建模问题。 u基于可变模糊集理论，建立了针对轨道交通隧道工程设计、 施工，运维阶段的模糊可变评价模型，加入轨道交通隧道工 程运维期的风险因素集，即可实现轨道交通隧道工程全生命 周期的风险等级评价。 u将时间因素集加入模糊可变评价模型，实现了隧道工程某些 地质风险的动态评估。 THANKS 谢谢聆听 41