扩展功能
文章信息
- 付金, 徐筱, 陈秉欣, 张泽辉, 张翾
- FU Jin, XU Xiao, CHEN Bing-xin, ZHANG Ze-hui, ZHANG Xuan
- 明挖公路隧道与邻接地铁协同施工风险耦合综合评价
- Comprehensive Evaluation of Risk Coupling of Collaborative Construction of Open Excavation Highway Tunnel and Adjacent Subway
- 公路交通科技, 2022, 39(8): 159-165
- Journal of Highway and Transportation Research and Denelopment, 2022, 39(8): 159-165
- 10.3969/j.issn.1002-0268.2022.08.021
-
文章历史
- 收稿日期: 2022-04-05
2. 交通运输部公路科学研究院, 北京 100088
2. Research Institute of Highway, Ministry of Transport, Beijing 100088, China
随着基础设施建设的不断推进,邻接工程并行建设与协同施工的情形不断涌现[1]。对于面向全域的复杂环境下邻接工程协同施工体系,施工安全管理需求突出,开展全域协同施工风险精细管控意义重大[2-3]。
对于协同施工体系,施工事故往往伴随着多种风险因素共同作用,在其施工风险综合评价及精细管控中需考虑风险间的耦合作用[4]。目前,风险评估中风险耦合作用渐被关注,结合实践应用,形成了初步的评价方法。其中,王乾坤等[5]采用C-OWA算子和相互作用矩阵法构建了地铁深基坑施工风险耦合评价方法; 文艳芳等[6]基于PSR-IAHP模型构建了地铁隧道施工坍塌风险耦合体系; 沈阳等[7]采用指标折减法考虑了超高层施工风险因素耦合。然而,对面向全域的复杂环境下明挖公路隧道与邻接地铁协同施工风险耦合研究还鲜有报道[8],将风险耦合引入协同施工风险研究中,以期探明面向全域的复杂环境下协同施工风险因素之间的耦合作用并量化耦合强度。
城市公路隧道施工中事故的发生往往是周围环境、施工、设计、现场管理等多种致灾因素共同作用的结果,传统的分析方法多以单一风险源为出发点,评估单一风险源诱发事故的概率,并未从事故的机理出发,讨论多种致灾因素共同致灾的耦合机制[9-12]。明挖公路隧道与地铁工程协同施工过程中安全事故致灾因素众多,相互之间的耦合关系也较为复杂。通过耦合路径和耦合程度的度量实现对明挖公路隧道与邻接地铁协同施工风险事件产生机理的阐释,同时为施工风险防控提供对策与建议,从源头上切除或削弱风险耦合事故的产生条件,避免强耦合效应的产生。
本研究针对明挖公路隧道与邻接地铁协同施工风险耦合综合评价,构建了风险耦合评价指标体系,明确了风险源及其耦合信息分析框架,建立了明挖公路隧道与邻接地铁协同施工风险耦合综合评价模型。结合深圳公常路明挖公路隧道与邻接地铁6号线支线协同施工项目,进行了实例分析,度量了风险耦合作用程度及各风险因素之间的耦合效应,以期为面向全域的明挖公路隧道与邻接地铁协同施工项目的风险评价与安全管理提供支撑,同时为后续近接同期、协同施工系统风险耦合分析提供参考。
1 风险耦合评价指标体系影响明挖公路隧道与邻接地铁工程协同施工的风险事件主要包括施工变形超限、支护体系破坏(支撑失稳、基坑崩塌、围护结构位移过大、坑底隆起)、护壁基底渗漏突涌破坏(土体渗透破坏、降水引起周围地面沉降、基底突涌)、施工工期延误等,风险耦合研究首先要通过风险识别与因子分析,建立风险指标体系。
采用WBS工作分解法[13],分析明挖下穿公路隧道与邻接地铁工程协同施工安全事故原因。如明挖公路隧道基坑开挖卸荷过程中,深基坑工程各施工程序、工艺错误(开挖高差太大,开挖土方的空间尺寸、暴露时间等时空效应考虑不足,各工种、工序之间协调和保护不周等),支护结构设计不合理(土体强度指标选择不合理、围护支撑体系强度不够等),邻接基坑开挖、降水的影响,周围环境监测不合理,施工单位施工质量差,施工管理混乱等,多种因素导致安全风险事件发生,即事故发生通常是一系列的失误、缺陷共同作用的结果。对于面向全域明挖公路隧道与邻接地铁工程协同施工风险耦合的复杂性,通过采用理论与应用相结合等研究方法,基于风险耦合致灾机理,明确风险源及其耦合信息的获取方式。运用风险耦合模型,获取风险状态的准确估计[14]。本研究将明挖公路隧道与邻接地铁工程协同施工安全风险主要因素分为同步施工因素、开挖因素、围护结构因素、降水止水因素和项目管理因素这5个方面,如表 1所示。
| 1级指标 | 2级指标 |
| 同步施工因素F1 | 近接全域施工资料完善度不足,特殊地段针对性设计不足,与邻接基坑近接程度较近,近接工程施工组织顺序不当,近接工程施工方法不当,近接控制支护措施不足,监测点布置不合理,监测误警/漏警 |
| 开挖因素F2 | 土层信息不完整,基坑周围环境信息复杂,基坑深度较深,开挖方式不合理,开挖工序不合理,机械设备故障,土方回填不合理 |
| 围护结构因素F3 | 基坑支护方案选择不当,内支撑间距不合理,围护结构入土深度不够,围护结构刚度/强度不够,围护结构施工不合理,围护结构施工质量不高,支护不及时,围护结构失效,坡顶堆载过大,坑内土体加固不到位,基底回弹隆起 |
| 降水止水因素F4 | 坑外降水不力,地下水或雨水浸泡,水压过大,积水使土体力学性能降低,地下水位下降引起地表沉降,防排水措施有效性不足,地下管线破裂,围护体渗漏 |
| 项目管理因素F5 | 近接工程协同管理水平不高,管理模式混乱,施工现场管理水平不高,人员数量及水平不高,国家或地方政策改变,技术装备不齐全,信息反馈不及时,作业标准程度不高,规章制度不健全 |
2 明挖公路隧道与邻接地铁协同施工风险耦合模型
耦合度模型可对风险因素间耦合作用的强弱进行定量评价并较为直观地展现各风险因素间的耦合作用状态[15]。明挖公路隧道与邻接地铁协同施工风险耦合模型的基本原理:在明确风险耦合致灾机理(耦合信息)及评价指标体系基础上,确定风险指标的权重和特征值,构建功效函数及耦合度函数,计算系统耦合度,获取风险耦合状态的准确估计,得到明挖公路隧道与邻接地铁协同施工安全评价体系。
2.1 确定指标权重及指标值基于层次分析法(AHP)来计算指标权重,即将定性问题与定量数值相结合,采用1-9比例标度法构建明挖公路隧道协同施工安全风险耦合评价指标的判断矩阵并计算指标权重ωi。
采用DEMATEL-逆向云模型相结合的方法测定风险特征的定量数值[16],得到云特征值Ex,En,He。其中,Ex为期望值,En为熵,He为超熵,主要计算过程为:
先计算样本均值和方差:
|
(1) |
|
(2) |
式中xi为样本数据;n为样本数量。
再得到期望值、熵及超熵:
|
(3) |
|
(4) |
|
(5) |
基于DEMATEL确立各风险因素之间因果关系和影响程度,构建评估矩阵。通过对评估矩阵的计算,评价明挖公路隧道协同施工安全风险耦合状态。
2.2 确定风险因素功效耦合度模型能够直观、定量地表达出协同施工安全风险耦合评价体系内风险因素的强弱关系,其中通过功效函数来量化不同风险因素对导致事故(风险体系)的影响程度(综合作用水平)。功效函数模型为:
|
(6) |
式中,Uij为明挖公路隧道与邻接地铁协同施工风险耦合体系的功效系数(第i类因素的第j个因子),即风险因子对施工安全的影响程度,其中正功效表示关注因素和安全风险正相关,其值越大对施工安全的影响越大;Exij为前述期望值;Aij和Bij为风险评估矩阵序参量的上、下限值。
基于得到的功效系数Uij,采用线性加权和法,可得各风险因素的功效Ui:
|
(7) |
以此量化各风险因素对施工安全风险耦合系统的综合作用水平。
2.3 耦合度模型构建将容量耦合模型引入明挖公路隧道与邻接地铁协同施工风险耦合系统[17-18],以耦合度来度量风险耦合(因素相互作用)的程度。基于前述得到的风险因素的功效Ui,构建明挖公路隧道与邻接地铁协同施工系统的耦合度模型:
|
(8) |
式中,Cm为整个系统风险耦合的耦合度值;m为明挖公路隧道与邻接地铁工程协同施工系统风险耦合涉及的风险因素数。
依据耦合模型中耦合状态的划分原则,结合明挖公路隧道施工安全风险评价需求,将协同施工安全风险耦合评价等级划分为较低水平、中等水平和高等水平3个级别[10]。由式(9)可知,耦合度的值Cm∈[0, 1]。当Cm=0时,表示风险耦合作用关系最弱;Cm=1时,表示风险耦合作用关系最强。具体来说,Cm∈[0, 0.3]时风险耦合综合水平较低;Cm∈(0.3, 0.7]时风险耦合水平中等;Cm∈(0.7, 1]时风险耦合水平高。
3 实例分析本研究以深圳市公常路明挖下穿公路隧道与邻接地铁6号线支线协同施工项目为例,该协同施工项目中公常路下穿隧道全长3.56 km,其中地下道路长2.6 km(封闭段2.22 km,敞开段0.425 km)。6号线支线项目包括中山大学站、武汉大学站2座地铁车站主体(开挖深度分别为20 m和18.5 m左右)及出入口和区间隧道,施工难度较大。本公路隧道距离中山大学站最小水平净距3 m,距离武汉大学站5 m,且2项工程同期施工。地铁武大站C出入口横下穿本公路隧道,地铁中大站出入口上跨本隧道,地铁盾构隧道与本项目隧道最小水平净距5 m。
该项目位于深圳市西北部的光明新区的公常路左幅,近东西向布置,为建(构)筑物密集地段。项目实施过程中就施工组织,尤其是相邻深基坑开挖等交叉施工,与地铁方进行了多次协调,达成了合理的解决方案,保障了施工安全。
3.1 权重计算结合项目现场施工实际情况,构建风险因素判断矩阵,如1级指标5类风险因素间的风险判断矩阵为:
|
采用MCE软件中的AHP模块计算判断矩阵及权重数值,明挖公路隧道与邻接地铁协同施工风险指标体系权重计算结果如表 2所示。
| 1级指标 | 2级指标 | 1级指标权重 | 2级指标权重 |
| 同步施工因素风险 | 近接全域地铁施工资料完善度不足 | 0.557 7 | 0.016 6 |
| 特殊地段针对性设计不足 | 0.311 7 | ||
| 与临近基坑近接程度较近 | 0.034 3 | ||
| 近接工程施工组织顺序不当 | 0.048 3 | ||
| 近接工程施工方法不当 | 0.047 3 | ||
| 近接控制支护措施不足 | 0.352 7 | ||
| 监测点布置不合理 | 0.077 | ||
| 监测误警/漏警 | 0.112 | ||
| 开挖因素风险 | 土层信息不完整 | 0.087 8 | 0.054 2 |
| 基坑周围环境信息复杂 | 0.087 7 | ||
| 基坑深度较深 | 0.076 2 | ||
| 开挖方式不合理 | 0.201 7 | ||
| 开挖工序不合理 | 0.201 7 | ||
| 机械设备故障 | 0.020 1 | ||
| 土方回填不合理 | 0.358 4 | ||
| 围护结构因素风险 | 基坑支护方案选择不当 | 0.183 8 | 0.027 |
| 内支撑间距不合理 | 0.095 9 | ||
| 围护结构入土深度不够 | 0.062 4 | ||
| 围护结构刚度/强度不够 | 0.077 9 | ||
| 围护结构施工不合理 | 0.023 1 | ||
| 围护结构施工质量不高 | 0.032 3 | ||
| 支护不及时 | 0.212 4 | ||
| 围护结构失效 | 0.308 2 | ||
| 坡顶堆载过大 | 0.014 8 | ||
| 坑内土体加固不到位 | 0.019 1 | ||
| 基底回弹隆起 | 0.127 | ||
| 降水止水因素风险 | 坑外降水不力 | 0.141 2 | 0.032 7 |
| 地下水或雨水浸泡 | 0.106 5 | ||
| 水压过大 | 0.023 1 | ||
| 积水使土体力学性能降低 | 0.231 9 | ||
| 地下水位下降引起地表沉降 | 0.158 5 | ||
| 防排水措施有效性不足 | 0.071 3 | ||
| 地下管线破裂 | 0.047 9 | ||
| 围护体渗漏 | 0.328 | ||
| 管理因素风险 | 近接工程协同管理水平不高 | 0.029 5 | 0.309 |
| 管理模式混乱 | 0.219 9 | ||
| 施工现场管理水平不高 | 0.156 8 | ||
| 人员数量及水平不高 | 0.073 6 | ||
| 国家或地方政策改变 | 0.0195 | ||
| 技术装备不齐全 | 0.0357 | ||
| 信息反馈不及时 | 0.0517 | ||
| 作业标准程度不高 | 0.1075 | ||
| 规章制度不健全 | 0.0262 |
3.2 风险指标值
根据上述建立的项目施工安全风险评价体系,基于DEMATEL-逆向云模型,对各风险因素的直接影响程度赋值,标定该协同项目施工安全风险指标值。数据来源于相关决策试验,通过确定风险指标的得分,利用云发生器计算指标的风险值。基于式(1)~式(5),得到的云模型结果如表 3所示,这里仅示出1级指标标定结果,其中Aij=1,Bij=0。
| 1级指标 | 权重值 | 云模型 |
| 同步施工因素风险F1 | 0.557 7 | (0.928,0.056,0.054) |
| 开挖因素风险F2 | 0.087 8 | (0.485,0.217,0.212) |
| 围护结构因素风险F3 | 0.183 8 | (0.48,0.486,0.485) |
| 降水止水因素风险F4 | 0.141 2 | (0.48,0.361,0.359) |
| 管理因素风险F5 | 0.029 5 | (0.23,0.496,0.495) |
3.3 耦合度计算
结合指标权重及指标值,基于式(7)给出风险因素的功效。结合指标功效,基于式(8),可得到该明挖公路隧道与邻接地铁工程协同施工不同风险耦合状态下的耦合度,这里仅从1级指标角度进行说明,具体结果如表 4所示。C (F1, F2)为F1(同步施工因素)和F2(开挖因素)之间的风险耦合(因素相互作用)的程度,其他同理。
| 耦合情况 | 风险耦合度 | ||||
| 双因素 | C (F1, F2) = 0.677 | C (F1, F3) = 0.674 | C (F1, F4) = 0.663 | C (F1, F5) = 0.466 | C (F2, F3) = 0.487 |
| C (F2, F4) = 0.479 | C (F2, F5) = 0.337 | C (F3, F4) = 0.477 | C (F3, F5) = 0.335 | C (F4, F5) = 0.330 | |
| 三因素 | C (F1, F2, F3) = 0.604 | C (F1, F2, F4) = 0.598 | C (F1, F2, F5) = 0.473 | C (F1, F3, F4) = 0.595 | C (F1, F3, F5) = 0.471 |
| C (F1, F4, F5) = 0.466 | C (F2, F3, F4) = 0.480 | C (F2, F3, F5) = 0.379 | C (F2, F4, F5) = 0.375 | C (F3, F4, F5) = 0.374 | |
| 四因素 | C (F1, F2, F3, F4) = 0.566 | C (F1, F2, F3, F5)= 0.474 | C (F1, F3, F4, F5)= 0.469 | C (F2, F3, F4, F5)= 0.399 | — |
| 五因素 | C (F1, F2, F3, F4, F5)= 0.473 | — | — | — | — |
3.4 风险耦合分析
基于表 4得到的风险耦合状态量化值,来综合评价明挖公路隧道与邻接地铁工程协同施工风险耦合状态:(1)风险耦合体系中双因素耦合均为中等强度耦合,F1(同步施工因素)和F2(开挖因素)之间的耦合度最大,为0.677,接近高强度耦合。(2)在双因素、三因素、四因素耦合中,均表现出当有F1(同步施工因素)参与的耦合时,耦合度均较大,这表明同步施工因素对明挖公路隧道与邻接地铁工程协同施工安全风险影响最为显著。施工中开挖、围护、降止水等技术方案的确定及施工现场管理均应注重考虑同步施工因素的影响,降低同步施工风险因素的耦合作用程度,保障施工安全。(3)F1-F2耦合度最大表明,对于明挖公路隧道全域邻接复杂交叉协同施工项目,邻接深基坑同步开挖风险耦合程度最为显著。其中跟进、先行等相互间施工关系、施工时机、施工区段的间距、基坑开挖方式和工序、土方回填等因素相互作用(影响)明显,相互间风险耦合显著。施工中应重点关注相邻深基坑同步施工及基坑开挖二者风险耦合作用,加强沟通协调,与地铁方对接施工组织情况,切实考虑与地铁交叉施工的影响,预判存在的问题。对深基坑开挖选取适宜的开挖工法,即综合考虑2类因素,寻求合理的解决方案。(4)基于耦合度模型的评价结果对因素序参量的依赖程度较高,得到合理的风险因素序参量(因素功效)取值较为关键,后续研究可结合相近工程对应风险因素导致的风险事件发生的概率统计分析来给出因素功效。
4 结论(1) 通过对施工过程中施工变形超限、支护体系破坏等主要风险事件特征及机理解析,从同步施工、开挖、围护结构、降水止水和项目管理5个方面建立了明挖公路隧道与邻接地铁协同施工风险耦合评价指标体系,构建了耦合信息分析框架。引入耦合度模型,基于各因素权重、风险值及功效,构建了明挖公路隧道与邻接地铁协同施工风险耦合综合评价模型,度量了风险耦合(因素相互作用)的程度。
(2) 采用理论与应用相结合的方法,结合深圳公常路明挖公路隧道与邻接地铁6号线支线协同施工项目,进行了实例研究。基于得到的不同风险耦合状态下的耦合度量化值,综合评价了明挖公路隧道与邻接地铁工程协同施工风险耦合状态,结合项目现场情况验证了模型的可行性。
(3) 明挖公路隧道与邻接地铁协同施工风险耦合体系中,同步施工因素参与耦合时,耦合度均较大,即同步施工因素对施工安全影响较大。同步施工因素和开挖因素二者耦合度最大,即邻接深基坑同步开挖风险耦合程度最为显著。施工中开挖、围护、降止水等技术方案的确定及施工现场管理均应注重考虑同步施工因素的影响,同时对于相邻深基坑同步开挖,应选取合理的施工组织与开挖工法,降低风险耦合作用,保障施工安全。
| [1] |
薛晔, 刘耀龙, 张涛涛. 耦合灾害风险的形成机理研究[J]. 自然灾害学报, 2013(2): 44-50. XUE Ye, LIU Yao-long, ZHANG Tao-tao. Research on Formation Mechanism of Coupled Disaster Risk[J]. Journal of Natural Disasters, 2013(2): 44-50. |
| [2] |
胡兴俊. 建设项目施工安全风险耦合机理研究[D]. 上海: 上海工程技术大学, 2016. HU Xing-jun. Study on Coupling Mechanism of Safety Risk in Construction Projects[D]. Shanghai: Shanghai University of Engineering Science, 2016. |
| [3] |
王勇, 赖芨宇, 陈秋兰, 等. 我国地铁施工事故统计分析与研究[J]. 工程管理学报, 2018, 32(4): 70-74. WANG Yong, LAI Ji-yu, CHEN Qiu-lan, et al. A Statistical Analysis of Metro Construction Accidents in China[J]. Journal of Engineering Management, 2018, 32(4): 70-74. |
| [4] |
陈绍清, 熊思斯, 何朝远, 等. 地铁深基坑坍塌事故安全风险分析[J]. 安全与环境学报, 2020, 20(1): 52-58. CHEN Shao-qing, XIONG Si-si, HE Chao-yuan, et al. Safety Risk Analysis for the Deep Foundation Pit Collapse Accidents in the Subway Construction System[J]. Journal of Safety and Environment, 2020, 20(1): 52-58. |
| [5] |
王乾坤, 亢显卫, 朱科. 地铁深基坑施工风险耦合评价方法[J]. 东北大学学报(自然科学版), 2021, 42(8): 1152-1158. WANG Qian-kun, KANG Xian-wei, ZHU Ke. Coupling Evaluation Method of the Construction Risk for Subway Deep Foundation Pit[J]. Journal of Northeastern University (Natural Science Edition), 2021, 42(8): 1152-1158. |
| [6] |
文艳芳, 陈敬配. 地铁隧道施工坍塌风险耦合机理研究[J]. 地下空间与工程学报, 2021, 17(3): 943-952. WEN Yan-fang, CHEN Jing-pei. Research on Coupling Mechanism of Collapse Risk in Metro Tunnel Construction[J]. Chinese Journal of Underground Space and Engineering, 2021, 17(3): 943-952. |
| [7] |
沈阳, 徐磊, 郑冠雨, 等. 考虑风险因素耦合的超高层施工预警方法研究[J]. 郑州大学学报(工学版), 2021, 42(4): 98-104. SHEN Yang, XU Lei, ZHENG Guan-yu, et al. Research on Early-warning Method for Super High-rise Building Construction Considering Coupling Effect of Risk Factors[J]. Journal of Zhengzhou University (Engineering Science Edition), 2021, 42(4): 98-104. |
| [8] |
黄芮. 深基坑开挖对近邻隧道的影响及其风险传递研究[D]. 北京: 北京交通大学, 2019. HUANG Rui. Research on Effect of Deep Foundation Pit Excavation on Adjacent Tunnel and Corresponding Risk Transfer[D]. Beijing: Beijing Jiaotong University, 2019. |
| [9] |
徐涛. 多因素耦合作用下的水下隧道盾构施工安全风险控制研究[D]. 重庆: 重庆交通大学, 2016. XU Tao. Study on Safety Risk Control of Shield Construction in Underwater Tunnel with Multiple Factor Coupling[D]. Chongqing: Chongqing Jiaotong University, 2016. |
| [10] |
马伶伶. 城市轨道交通运营风险耦合研究[D]. 北京: 北京交通大学, 2014. MA Ling-ling. Study on Operational Risk Coupling of Urban Rail Transit[D]. Beijing: Beijing Jiaotong University, 2014. |
| [11] |
侯公羽, 刘伟, 李乐, 等. 多风险因素耦合的地铁施工安全系统脆弱性分析[J]. 土木工程学报, 2022, 55(2): 111-119. HOU Gong-yu, LIU Wei, LI Le, et al. Vulnerability Analysis of the Subway Construction Safety System with Coupled Multiple Risk Factors[J]. China Civil Engineering Journal, 2022, 55(2): 111-119. |
| [12] |
PAN H, GOU J, WAN Z, et al. Research on Coupling Degree Model of Safety Risk System for Tunnel Construction in Subway Shield Zone[J].
Mathematical Problems in Engineering, 2019(5): 1-19.
|
| [13] |
王志超. 基于故障树的西安地铁深基坑工程施工安全风险识别[D]. 西安: 西安工业大学, 2017. WANG Zhi-chao. Safety Risk Identification of Xi'an Metro Deep Foundation Pit Construction Based on Fault Tree[D]. Xi'an: Xi'an Technological University, 2017. |
| [14] |
刘升福. 基于战略转换的战略风险研究[J]. 管理科学, 2004(2): 32-38. LIU Sheng-fu. Research on the Strategic Risk Based on Strategic Transition[J]. Journal of Management Science, 2004(2): 32-38. |
| [15] |
李新春, 韩继磊, 乔万贯. 煤矿事故风险因子耦合机理研究[J]. 煤矿安全, 2015, 46(8): 240-242. LI Xin-chun, HAN Ji-lei, QIAO Wan-guan. Study on Coupling Mechanism of Mine Accident Risk Factors[J]. Safety in Coal Mines, 2015, 46(8): 240-242. |
| [16] |
孙宏才. 网络层次分析法与决策科学[M]. 北京: 国防工业出版社, 2011. SUN Hong-cai. Network Analytic Hierarchy Process and Decision Science[M]. Beijing: National Defense Industry Press, 2011. |
| [17] |
李杰, 滕佳颖, 王富. 城市多突发事件耦合交通预警指标研究[J]. 武汉工程大学学报, 2012, 34(3): 60-63. LI Jie, TENG Jia-ying, WANG Fu. Survey of City Transportation Early Warning Indicator Based on Multi-emergency Coupling[J]. Journal of Wuhan Institute of Technology, 2012, 34(3): 60-63. |
| [18] |
陈梦捷. 地铁隧道盾构区间施工风险耦合模型研究[D]. 成都: 西南石油大学, 2016. CHEN Meng-jie. Study on Risk Coupling Model of Shield Section Construction of Subway Tunnel[D]. Chendu: Southwest Petroleum University, 2016. |
2022, Vol. 39

