Model Test Study on Deformation Transfer between Overlying Soil Layer and Tunnel under Reverse Fault Action
-
摘要: 基于自行设计的错动模型试验装置和数字图像相关技术开展1∶80模型试验，研究60°倾角逆断层错动作用导致上覆土层剪切破裂的过程。依托数字图像相关技术非接触式全场测量的优势，分析上覆土层与隧道相互作用对上覆土层剪切破裂扩展、上覆土层变形和地表变形的影响，总结60°倾角逆断层错动作用下上覆土层与隧道之间的变形传递形式。研究结果表明，与自由场试验结果相比，由于隧道与上覆土层变形不同步，剪切破裂遇到隧道时会产生分叉，即隧道能够偏移上覆土层剪切破裂路径；在逆断层作用下，由于上覆土层与隧道力学性能存在差异，二者不能同步变形，为适应剪切区上覆土层的大变形，隧道周边土体会出现脱空，不利于隧道抗震；与自由场试验结果相比，隧道在影响破裂路径的同时能够将剪切变形扩散到更宽的区域，并在地表产生更大范围的陡坎。Abstract: Based on the self-designed dislocation model test device and DIC technology, the 1∶80 model test was carried out to study the process of shear fracture of overlying soil caused by the action of a 60° reverse fault. Relying on the advantages of non-contact whole-field measurement of DIC technology, the influence of the interaction between the overlying soil layer and the tunnel on the shear fracture propagation of the overlying soil, the deformation of the overlying soil and the surface deformation is emphatically analyzed, and the deformation transfer mechanism between the overlying soil layer and the tunnel is summarized. The test results show that, compared with the free field test results, because the rigidity of the tunnel is far greater than that of the overburden, the shear fracture and the tunnel encounter will produce bifurcation, that is, the tunnel can offset the shear fracture path of the overburden; Under the action of reverse fault, the overlying soil layer and the tunnel cannot deform synchronously due to the difference of mechanical properties. To adapt to the large deformation of the overlying soil layer in the shear zone, the soil around the tunnel will appear void, which is unfavorable to the seismic resistance of the tunnel. Compared with the free field test results, the tunnel can spread the shear deformation to a wider area while affecting the fracture path, and produce a wider range of scarps on the surface.
-
Key words:
- Tunnel /
- Fault zone /
- Fault dislocation /
- Shear deformation /
- Interaction /
- Digital image correlation (DIC)
-
引言
地壳速度模型和地震研究有着密不可分的联系，如何得到合适的地壳速度模型是地震学中的经典问题。准确的速度模型对于地震定位、震源机制、矩张量反演（赵韬等，2016）等研究至关重要，同时也可以约束地壳介质的组分。区域地震台网的日常工作，就是为了获得可靠的地震基本参数（发震时刻、震中位置、震源深度和震级），同时也为地震预测和地震学研究工作提供基础资料。而地壳速度模型在很大程度上决定了地震的定位精度，合理的地壳速度模型有助于对地震测定精度进行可靠估计（朱元清等, 1997, 2002；吕作勇等，2014；罗佳宏等，2017）。长期以来，地震工作者用各种方法深入研究地壳速度模型，以进一步揭示地震活动和地球内部的奥秘。
20世纪80至90年代，多位研究者利用地震模拟记录资料及爆破资料在陕西地区开展了地壳速度模型相关研究工作（姜家兰，1985；丁韫玉等，1987）。目前陕西地区采用的地壳速度模型（1985模型）就是姜家兰等（1985）利用陕西省7次爆破和邻近地区72个地震模拟记录资料计算得到的。近年来，随着数字地震学的发展，一些学者利用数字地震资料对陕西一些局部地区的速度模型进行了研究，如：利用联合走时与波形反演技术研究陕、甘、宁、青地壳三维速度结构（张元生等，1998；郭瑛霞等，2017），渭河断陷盆地及邻近地区地震活动的深部背景（王卫东等，2004），玛沁—兰州—靖边地震测深剖面速度结构的初步研究（李松林等，2002），渭河盆地及邻区地壳深部结构特征研究（任隽等，2012），秦岭造山带与邻域华北克拉通和扬子克拉通的壳、幔精细速度结构与深层过程（滕吉文等，2014）等。尽管利用这些方法，获得了陕西部分区域比较精确的地壳速度结构，但是，这些结果很难较好地描述陕西地区整体速度结构特征。在实际应用中以彩色图形和等值线方式给出的二维剖面和三维速度结构更适用于科学研究。然而，目前测震台网地震分析工作中，定位程序还难以直接使用这类研究成果。
陕西1985地壳速度模型和现今的研究成果存在一定的差异，如：1985模型的康拉德界面为16km，滕吉文等（2014）及任隽等（2012）认为陕西地区的康拉德界面在20km以上。一直以来，采用1985模型进行地震速报和编目工作，与中国地震台网中心及8个邻省台网的地震定位结果差异较大，给地震监测预报以及应急救援等工作带来困惑。因此，急需建立更符合陕西区域地质特征的地壳速度模型。
本研究利用2009年1月—2014年10月间陕西测震台网的1560个地震，53843条震相数据，通过区域速度拟合、速度稳定性分析、折合走时分析，结合前人爆破测线、重力测深和层析成像等结果得到了陕西地区地壳速度模型中各层速度和厚度的平均值和扰动范围。在此基础上，采用Hyposat定位程序，在扰动范围内对满足地震精度条件的地震资料进行“试错”，批量测定其参数。然后，参照中国地震局正式地震编目结果，选取结果最优、总体残差最小的一组参数作为陕西地区地壳速度模型。得到模型后，过滤掉参与速度模型计算的地震资料，另外选取2015年以后的地震资料，运用批量定位残差比较、震源深度比较、实测爆破结果对比等方法对得到的地壳速度模型进行验证，最终确定陕西地区地壳速度模型（2015模型）。
1. 计算方法
拟合方法的总体思路（尹欣欣等，2017）为：根据陕西测震台网各台实际记录的走时和震中距资料，采用线性和滑动窗2种方法进行拟合。拟合前对数据进行预处理，滑动窗在整个资料包含的震中距范围内，以150km为固定窗长，并按50km的间隔进行滑动线性拟合，最后根据每个窗所获得的结果进行加权平均求出最后的速度结果。单个窗内样本数少于平均样本数的80%，则不参与加权平均计算。
加权平均的做法是：首先计算所有单个窗内获得结果的总平均值；然后将每个窗内的结果与平均值之比的倒数作权重系数，并将权重系数乘以相应窗对应的值后再平均。物理意义为偏离平均值越远的值，权重就越小。
本文线性拟合是求2组变量T、△之间的线性函数关系：
$$ T = A + {\rm{\Delta }}/{V_N} $$ (1) 式中，T是震相走时；△是震中距；VN是拟合速度；A是常数。
利用最小二乘方法，可求得A和VN，使得每个样本观测点与拟合线性方程之间的偏差最小。
根据线性拟合的速度结果，可以用折合走时的方法来调整初始模型并确定模型的试错范围。即给定地壳速度模型参数的扰动范围，确保真实地壳速度模型参数值包含在可变化的范围内。
计算折合走TZ时可用下式表示：
$$ {T_Z} = {T_L} - \Delta /V $$ (2) 式中，TL是实际观测走时；△是震中距；V是波速。用折合走时可对波速进行调整，使其更接近理论值。当V与实际地层速度变化一致时，观测值和理论值吻合，图像结果将汇聚成水平条带状。基于这个规律，我们可以将实际观测走时与理论计算走时叠加到一张图上分析。
通过以上工作，获得了陕西地区地壳速度模型中各层P波速度、厚度的平均值和扰动范围。在扰动范围内，用Hyposat定位程序并采用“穷举法”对地震资料进行“试错”，分2步进行：①取各层的速度变化步长为0.1km/s，层厚变化步长为1km；②缩小其扰动范围，取各层速度变化步长为0.01km/s，层厚的变化步长为0.5km，批量测定这些地震事件的空间位置。参照正式地震编目结果，选取结果最优、总体残差最小的一组参数作为区域地壳速度模型。然后固定P波速度及各层厚度，使用地震资料，将波速比在1.68—1.75、间隔0.01取值进行批量定位，取平均定位残差最小的波速比为最优波速比。
地球内部随着密度变大，温度升高，S波速度比P波速度增加慢，波速比会随着深度变大，求得合适的分层波速比，就可以换算得到分层的S波速度。取得上述地震模型后，采用Hyposat定位程序，分2步计算得到各层的波速比：
第1步：将已经求得的地壳速度模型第1、2层波速比固定为最优波速比，“试错”莫霍面波速比，将综合残差最小且震中差也最小所对应的波速比，确定为莫霍面波速比；
第2步：固定第1层和莫霍面波速比，“试错”第2层波速比，不断统计定位残差和检验震中位置，选取残差最小且满足震中误差更小的模型为下一轮的初始值。不断重复以上过程，通过调整各层波速比，最后综合选取总体残差和震中差均为最小的参数，作为最终陕西地区地壳速度模型。
2. 资料及地震事件筛选
经过“十五”测震台网及“灾后恢复重建”建设，陕西省目前管理的国家和区域测震台站达52个，其中包括4个国家数字地震台，48个区域数字地震台；汇集周边20个台站（其中四川3个，湖北4个，山西3个，河南2个，宁夏2个，甘肃3个，内蒙古3个）及榆林市地方地震台7个，即同时接收79个台站的数据。陕西省内59个台站较均匀地分布于全省境内，平均台间距约60km，北部地区台站间距较大，约80km，东部及南部老震区，为加强地震监测，台站分布较为密集，台站间距约为50km。
考虑到汶川余震序列过多，若采用全部选取的资料会导致权重过大，对最终结果造成不利影响。为此，我们对资料进行了筛选，选取震相清楚、定位结果比较可靠的地震用于速度结构研究，剔除了位置相近、震相数量较少或震相不清楚的地震，同时考虑了地震完整性（王平等，2014）。最终，共选取2组地震资料进行研究：第1组选自2009年1月—2014年10月，陕西台网记录台站数N≥10的1560个地震事件（其中Pg震相21170个，Pn震相8258个，Pb震相106个，Sg震相21038个，Sn震相3377个）（图 1（a））。第2组资料从第1组资料中选择陕西台网记录台站数不少于同期台网台站总数（包括本省和接入的临省台站）80%的52次地震事件（其中Pg震相1495个，Pn震相903个，Pb震相106个，Sg震相1334个，Sn震相335个）（图 1（b））。第1组数据主要用于确定初始模型，包括各层速度、厚度和扰动范围以及数据稳定性分析；第2组数据用于最终地壳速度模型的确定。
3. 初始模型及最优模型的确定
参考前人对陕西地区地壳速度模型的研究成果，并依据地震、爆破及塌陷的震相速度拟合曲线与折合走时曲线等结果，对上地壳速度VPg，下地壳速度VPb、上地幔速度VPn、上地壳厚度H1和下地壳厚度H2等5个参数进行合理调整，建立初始模型。调整原则为：删除与理论走时曲线相差4s的震相数据；在折合走时图中，调整厚度使实际震相数据位于理论走时曲线中间；调整层速度、震源深度使观测数据与理论走时曲线相一致。
依据震相特征，对现有编目资料中未标注的Pb震相进行识别。考虑到震源深度的影响，天然地震拟合选取震中距大于50km的台站。根据陕西省测震台网跨度、分布及考虑地球扁率等的影响，震中距上限选取550km。
基于第1、2组数据及抽取的18次非天然地震，进行Pg、Sg、Pb、Pn和Sn震相速度拟合，得到了5种震相对应的VPg、VSg、VPb、VPn和VSn拟合曲线（图 2）。从拟合结果可知，第1组数据（图 2（a））、第2组数据（图 2（b））和非天然地震（图 2（c））的Pg波速VPg分别为6.06km/s、6.06km/s和6.10km/s，平均为6.07km/s；Pn波速VPn分别为7.82km/s、7.99km/s和8.09km/s，平均为7.97km/s；Sg波速VSg分别为3.56km/s和3.58km/s，平均为3.57km/s；Sn波速VSn分别为4.47km/s和4.46km/s，平均为4.47km/s；Pb波速VPb同为6.55km/s；P波与S波的波速比VPg/VSg≈1.7。
对第1组数据进行速度稳定性分析，分别考察VPg、VSg、VPn和VSn在下面3种情况下的稳定性（表 1）：①随着震中距增大时的变化情况；②震中距以150km为窗口，50km为步长，速度值的变化情况；③震中距以200km为窗口，50km为步长，速度值的变化情况。结果表明：拟合的趋势性是一致的，但具体的速度范围存在一定差别；速度的平均值比较接近，线性拟合的速度具有一定的稳定性。
表 1 3种情况下速度稳定性Table 1. Velocity stability in three casesVPg/km·s－1 VSg/km·s－1 VPn/km·s－1 VSn/km·s－1 随区域增大 速度范围 6.07—6.08 3.53—3.56 7.05—7.80 3.98—4.46 平均速度 6.08 3.55 7.43 4.22 不同区域范围
（150km）速度范围 6.03—6.12 3.53—3.61 7.05—8.18 3.98—4.63 平均速度 6.08 3.57 7.62 4.30 不同区域范围
（200km）速度范围 5.94—6.10 3.54—3.60 7.05—8.14 4.12—4.62 平均速度 6.02 3.57 7.60 4.37 基于陕西1985模型（表 2），使用公式（2）分别计算第1、2组地震事件数据和抽取的18次非天然地震的P波折合走时曲线（图 3（a）—（c）为陕西模型结果，图 3（d）—（f）为速度模型调整后的结果），调整后的模型为：VPg＝6.03km/s，VPb＝6.68km/s，VPn＝7.95km/s，H1＝19km，H2＝21km。从图（3）可知，Pg、Pb和Pn震相理论折合走时和实测折合走时大致平行，但实测折合走时与理论折合走时水平位置有一定偏离。折合走时曲线显示VPg、VPb、VPn、H1和H2的实际值与理论值不吻合。不同震源深度速度拟合对比表明，VPg实际走时与理论走时不符，需要增大VPg值；VPb实际走时也与理论走时不符，需要增大VPb值，并调整地壳各层厚度。因此，原陕西模型与陕西区域地质特征存在一定的差异。根据陕西模型获得的速度拟合曲线、折合走时曲线以及不同震源深度速度拟合对比结果，得到了陕西地区速度初始模型（表 3）。通过不同深度速度理论与实际走时的对比可知，初始模型还需进一步修正。
表 2 陕西1985地壳速度模型Table 2. Crustal velocity model of Shaanxi Province in 1985VP/km·s－1 VS/km·s－1 H/km 5.90 3.44 16 6.50 3.70 24 7.95 4.55 综合天然地震和非天然地震数据速度拟合曲线及折合走时曲线，得到陕西地区地壳速度模型各层的平均值及其扰动范围（表 3）。
表 3 初始模型Table 3. Initial model地壳分层 VP/km·s－1 扰动范围/km·s－1 H/km 扰动范围/km 上地壳 6.09 5.8—6.1 24 16—24 下地壳 6.67 6.4—6.8 15 上地幔 7.99 7.8—8.1 38—42 基于初始模型，对第2组数据利用Hyposat定位方法进行批处理计算，获得的最优模型做为陕西模型的P波结构；通过多轮“试错”获得各层的波速比，从而得到S波速度结构。综合选取总体残差和震中差均为最小的参数，作为最终的陕西地区地壳速度模型（表 4），简称2015模型。结果表明2015模型和1985模型有较大差异（图 4）。
表 4 2015地壳速度模型Table 4. Crustal velocity model in 2015地壳分层 P波速度/km·s－1 S波速度/km·s－1 厚度/km 上地壳 6.09 3.58 24 下地壳 6.67 3.83 15 莫霍面 7.99 4.46 
图 4 陕西地区2015模型与1985模型对比图Figure 4. Contrast of the 2015 model to the 1985 model in the Shaanxi region针对第2组数据，使用Hyposat方法分别对1985模型和2015模型进行定位，得到了2种模型的定位走时残差（图 5（a））。从图中可以看出，1985模型走时残差分布范围为0.182—1.659s，均值为0.494s； 2015模型走时残差分布范围为0.129—0.605s，均值为0.306s。对比结果显示，2015模型定位走时残差均值降低了0.188s。
图 5 走时残差及震中差对比图Figure 5. Comparison of travel time residual and epicentral difference采用2015模型，使用Hyposat方法对52个地震定位，其定位震中结果与1985模型结果对比，分别得到了2015模型和1985模型（Hyposat方法）定位结果与编目结果的震中差。图 5（b）中可以看出，2015模型与编目定位震中差范围在0.18—10.14km，均值为1.92km；1985模型与编目震中差范围在0.28—11.76km，均值为2.54km，新模型比旧模型平均减少0.62km。由于采用同一种定位方法，差异来自速度模型，说明2015模型确定的震中位置更稳定。
4. 实测爆破检验
2015年1月8日01:40陕西蒲城发生ML 1.6爆破。随后，陕西省地震局进行了核实，确定了爆破地点。此次爆破在500km范围内有55个台站记录到较清晰震相，共记录Pg到时37个，Pn到时32个。采用同样震相和定位方法、不同速度模型的定位结果见表 5。
表 5 不同速度模型的实测爆破定位结果Table 5. Positioning results of measured blasting with different velocity models发震时刻 纬度/°N 经度/°E 深度/km 与实测震中差/km 残差/s 说明 01:40:14.523 34.81 109.39 实测结果 01:40:14.090 34.81 109.37 1 1.94 0.497 1985模型 01:40:15.400 34.81 109.39 1 1.14 0.321 2015模型 检验结果表明2015模型震中位置更接近实测爆破位置，且残差更小，更符合陕西地区的地质构造特征。
5. 结论
通过对陕西地区地壳速度模型的研究及确定，得到以下结论：
（1）滕吉文等（2014）及任隽等（2012）认为陕西地区的康拉德界面在20km以上，本文基于陕西地震台网走时资料建立的新速度模型的康拉德界面深度为24km，表明目前陕西测震台网使用的1985模型康拉德界面偏浅。
（2）以前人研究成果为基础，依据地震、爆破及塌陷的震相速度拟合曲线与折合走时曲线等结果，确定初始模型及扰动范围。利用Hyposat定位法确定的P波最优模型为VPg＝6.09km/s，VPn＝6.67km/s，VPn＝7.99km/s，H1＝24km，H2＝15km；采用“试错”法确定的波速比分别为上地壳1.70，下地壳1.74，上地幔1.79。
（3）利用新模型，采用Hyposat定位软件对地震进行重定位，结果表明，2015模型比1985模型定位残差平均低近0.2s，震中差平均低0.62km。利用实测爆破对模型进行检验，结果表明2015模型更接近实测爆破结果，且残差更小，进一步表明了新模型更符合陕西地区的地质构造特征。
致谢: 感谢中国地震局监测预报司、中国地震台网中心、上海市地震局、陕西省地震局和“全国区域一维速度模型建设及推广使用”项目组全体成员，对本文研究工作给予的支持和帮助。 -
图 8 上覆土层位移云图 （错动位移30 mm）
Figure 8. Nephogram of overlying soil layer displacement （Dislocation 30 mm）
图 9 上覆土层剪应变云图（错动位移30 mm）
Figure 9. Shearing strain nephogram of overlying soil layer （Dislocation 30 mm）
图 10 模型试验与数值模拟结果对比验证（错动位移30 mm）
Figure 10. Comparison and verification of model test and numerical simulation （Dislocation 30 mm）
表 1 ISO标准砂基本物理参数
Table 1. Physical mechanical parameters of ISO standard sand
相对密度$ {G}_{{\rm{s}}} $ 最大孔隙比$ {e}_{{\rm{max}}} $ 最小孔隙比$ {e}_{{\rm{min}}} $ 样品的累计粒度分布百分数达到50%时所对应的粒径$ {d}_{50} $ 不均匀系数$ {C}_{{\rm{u}}} $ 曲率系数$ {C}_{{\rm{c}}} $ 2.643 0.848 0.519 0.21 1.542 1.004 
下载: 导出CSV
表 2 上覆土层剪切破裂带扩展关键参数
Table 2. Key parameters of shear fracture zone expansion of overlying soil layer
工况 有无隧道 水平传播距离/mm 地表梯度 地表扩展角/（°） 自由场 无 180.0 −0.35 19 隧道 有 220.0 −0.75、−1.18 37、50
下载: 导出CSV
表 3 地表变形关键参数
Table 3. Key parameters of surface deformation
工况 影响区范围/$ \mathrm{m}\mathrm{m} $ 地表位移曲线拐点位置/$ \mathrm{m}\mathrm{m} $ 地表位移曲线拐点位置倾角
/（°）上盘边界 下盘边界 自由场 400 400 −150 11 133 −267 隧道 467 467 −233、−83 4、9 133 −333
下载: 导出CSV
-
安韶, 陶连金, 边金等, 2020. 跨活动断裂带城市浅埋地铁隧道结构两阶段设计方法研究. 中南大学学报（自然科学版）, 51（9）: 2558—2570An S. , Tao L. J. , Bian J. , et al. , 2020. Study on two-level design method of urban shallow subway tunnel structure crossing active fault. Journal of Central South University (Science and Technology), 51(9): 2558—2570. (in Chinese) 程斌, 李得睿, 2022. 基于退相关DIC的疲劳裂纹全局动态测量方法. 力学学报, 54（4）: 1040—1050Cheng B. , Li D. R. , 2022. Full-field dynamic measurement method for fatigue cracks based on decorrelation DIC. Chinese Journal of Theoretical and Applied Mechanics, 54(4): 1040—1050. (in Chinese) 蒋建平, 章杨松, 罗国煜等, 2002. 南京地铁地基下稳定性因素分析及对策. 地下空间, 22（1）: 42—44Jiang J. P. , Zhang Y. S. , Luo G. Y. , 2022. Analysis on unstable factors for ground of Nanjing metro and their countermeasures. Underground Space, 22(1): 42—44. (in Chinese) 李晓博, 张亮, 2020. 乌鲁木齐地铁1号线穿越断裂带的设计与施工. 都市快轨交通, 33（1）: 70—76Li X. B. , Zhang L. , 2020. Design and construction of Urumqi metro line 1 crossing fault zone. Urban Rapid Rail Transit, 33(1): 70—76. (in Chinese) 刘学增, 林亮伦, 2011.75°倾角逆断层黏滑错动对公路隧道影响的模型试验研究. 岩石力学与工程学报, 30（12）: 2523—2530Liu X. Z. , Lin L. L. , 2011. Research on model experiment of effect of thrust fault with 75° dip angle stick-slip dislocation on highway tunnel. Chinese Journal of Rock Mechanics and Engineering, 30(12): 2523—2530. (in Chinese) 石吉森, 凌道盛, 徐泽龙等, 2018. 倾斜场地中逆断层错动对上覆土体影响的模型试验研究. 工程力学, 35（7）: 194—207Shi J. S. , Ling D. S. , Xu Z. L. , et al. , 2018. Model testing study on the influence of reverse faulting on overlaying soil under an inclined ground. Engineering Mechanics, 35(7): 194—207. (in Chinese) 徐锡伟, 赵伯明, 马胜利等, 2011. 活动断层地震灾害预测方法与应用. 北京: 科学出版社. Chang Y. Y. , Lee C. J. , Huang W. C. , et al. , 2015. Evolution of the surface deformation profile and subsurface distortion zone during reverse faulting through overburden sand. Engineering Geology, 184: 52—70. doi: 10.1016/j.enggeo.2014.10.023 Chen Z. Y. , Jia P. , 2019. Three-dimensional analysis of effects of ground loss on static and seismic response of shafts. Tunnelling and Underground Space Technology, 92: 103067. doi: 10.1016/j.tust.2019.103067 Lee J. W. , Hamada M. , 2005. An experimental study on earthquake fault rupture propagation through a sandy soil deposit. Structural Engineering / Earthquake Engineering, 22(1): 1 s—13 s. doi: 10.2208/jsceseee.22.1s Lin M. L., Chung C. F., Jeng F. S., 2006. Deformation of overburden soil induced by thrust fault slip. Engineering Geology, 88(1—2): 70—89. Lin M. L., Chung C. F., Jeng F. S., et al., 2007. The deformation of overburden soil induced by thrust faulting and its impact on underground tunnels. Engineering Geology, 92(3—4): 110—132. Pan B. , 2018. Digital image correlation for surface deformation measurement: historical developments, recent advances and future goals. Measurement Science and Technology, 29(8): 082001. doi: 10.1088/1361-6501/aac55b Rechenmacher A. L. , Finno R. J. , 2004. Digital image correlation to evaluate shear banding in dilative sands. Geotechnical Testing Journal, 27(1): 13—22. Sabagh M. , Ghalandarzadeh A. , 2020. Centrifugal modeling of continuous shallow tunnels at active normal faults intersection. Transportation Geotechnics, 22: 100325. doi: 10.1016/j.trgeo.2020.100325 Take W. A. , 2015. Thirty-sixth Canadian geotechnical colloquium: advances in visualization of geotechnical processes through digital image correlation. Canadian Geotechnical Journal, 52(9): 1199—1220. doi: 10.1139/cgj-2014-0080 Wang W. L. , Wang T. T. , Su J. J. , et al. , 2001. Assessment of damage in mountain tunnels due to the Taiwan Chi-Chi earthquake. Tunnelling and Underground Space Technology, 16(3): 133—150. doi: 10.1016/S0886-7798(01)00047-5 -
下载:
点击查看大图
计量
- 文章访问数: 141
- HTML全文浏览量: 26
- PDF下载量: 20
- 被引次数: 0
