Response Analysis of Twin Tunnels Under Horizontal and Vertical Earthquake Loading
-
摘要: 基于有效的土-结相互作用有限元数值模拟方法，利用有限元软件ABAQUS对水平及竖向地震共同作用下双线盾构隧道的地震响应进行分析研究。地震动输入选取近场地震Loma、ChiChi、Mammoth和WoLong的基岩水平及竖向加速度时程记录。结果表明，不同近场地震记录对隧道结构的作用不同，隧道的地震反应与场地性质及地震动的频谱特性密切相关。对比隧道在水平及竖向地震动共同作用下的响应与单向水平地震动作用下的响应，发现隧道的最大地震附加内力及其分布均发生较大的变化，在隧道结构抗震设计中需引起重视。另外，分析中还考虑了在双向地震动共同作用下，隧道间距、土-结接触面的摩擦系数、土-结相对刚度、输入的地震记录和竖向地震动相对强度对隧道地震响应的影响等，研究结果对隧道工程的抗震设计具有一定的参考价值。Abstract: Based on an effective Finite Element procedure for soil-structure interaction, the seismic response of twin tunnels under combined horizontal and vertical earthquake loading were analyzed by using the ABAQUS package. The horizontal and vertical acceleration time history records at bedrock from the near-field earthquakes of Loma, ChiChi, Mammoth and WoLong were selected as the input ground motions. The results showed that the seismic responses of the tunnels were different as the input earthquake records changed, and the seismic responses of the tunnels were related to the field properties and the vibration frequency characteristics of the input ground motions. By comparing the seismic responses under bi-directional and unidirectional earthquake loadings, it was found that both the distributions and the magnitudes of the seismic additional internal forces were changed greatly, which should be taken into account in the seismic design of tunnels. In addition, parametric studies were carried out on the effects of the space between two tunnels, friction coefficient of contact surface between soil and tunnels, relative stiffness of the soil-structure system, input earthquake motion and relative intensity of vertical earthquake motion. The results may provide a good reference for the seismic design of tunnels.
-
引言
地震预警是利用电磁波传播速度远大于地震动传播速度以及破坏性地震波（S波或面波）的传播速度小于破坏性较小的P波速度的原理，在地震发生后尽可能短的时间内，确定地震基本参数或直接估计预警目标区的地震动和破坏程度，在破坏性地震波到来前，对可能破坏区发出地震警报。地震预警系统可以快速地进行地震定位并根据震中距不同，提供几秒、十几秒、甚至几十秒的时间，将预警信息在破坏性地震波到达之前发送至用户，以帮助用户采取应急措施，达到减灾的目的（马强，2008；殷海涛等，2012；张红才等，2013）。目前，日本、墨西哥、罗马和中国台湾地区已经建立了地震预警系统并提供预警服务，美国、瑞士、意大利和中国大陆则已建立了试验性质的地震预警系统（Allen等，2009；赵兵，2011）。
天津是中国北方最大的港口城市和环渤海经济中心，是国家防震减灾的重点设防城市，具有发生强烈地震的背景和条件，历史上曾多次遭遇特别重大的地震灾害。1679年9月2日三河-平谷8级地震和1976年7月28日唐山7.8级地震都使天津市遭受了严重的破坏。其中，唐山7.8级地震使宁河和汉沽的地震烈度达到Ⅸ度，塘沽及市区大部分地区地震烈度达到Ⅷ度，从而使天津成为中国唯一遭受Ⅷ度地震灾害的特大型城市。因此，建设天津简易烈度计地震预警试验区，能够建立完善有效的地震预警机制，进一步提高天津地区的地震灾害防御能力，尽可能地避免未来强震造成的人员伤亡和经济损失。
1. 试验区建设概况
天津简易烈度计地震预警试验区的建设目标是在天津地区新建80个简易烈度计观测点，与天津行政区内具备实时传输能力的测震台站和强震动台站共同组建天津地震预警观测网络，其网络规模如图 1所示。观测网络东西孔径约127km，南北孔径约160km，平均台站间距约10km。结合天津地震观测台网现有布局特征，选择观测条件较好的中小学校内平房或二层楼房的一层作为简易烈度仪观测场址。仪器固定在观测室不高于30cm的承重墙上，采用锚固方式安装，如图 2所示。数据传输采用移动4G无线VPDN方式，VPDN为虚拟专用拨号网络，是建立在VPN基础上且隶属与于VPN的1种业务。根据天津市远程网络传输需求，烈度计端路由设备首先通过4G无线网络接入附近移动运营商的无线网，再通过移动运营商VPDN专用网接入地震行业网，实现节点端与地震系统内网的互联互通。移动4G无线VPDN组网流程如图 3所示。
2. 简易烈度计性能
综合考虑仪器的安装方式、性能指标以及京津冀地震烈度速报与预警示范区的整体规划，天津简易烈度计地震预警试验区的台站观测仪器选用GL-P2B型和MI3000型三分向型简易烈度计，均具备16位以上分辨率，水平向测量范围-2—2g，垂直向测量范围-1—3g，动态范围大于80dB。传感器采用微机械加速度传感器（MEMS），MEMS加速度传感器为微电子技术与微机械工程结合发展的1种惯性器件，与传统力加速度传感器相比，具有体积小、成本低等特点。将MEMS传感器应用于地震仪器中，可以较低的成本提高台网密度，并可实现烈度速报与地震预警功能（Chung等，2011；Wu，2013；D'Alessandro等，2014）。2种仪器的性能指标见表 1，仪器的幅频特性曲线见图 4。
表 1 简易烈度计性能指标Table 1. Performance index of simple intensity meter型号 数量 采集延迟 分辨率 动态范围 加速度平坦频率 线性误差 GL-P2B 40 ＜1s 28位 ＞80dB 0.1—20Hz ＜1% MI3000 40 ＜1s 16位 ＞80dB 0.1—20Hz ＜1% 
图 4 简易烈度计幅频特性曲线Figure 4. Amplitude frequency characteristics of simple intensity meter3. 地震预警数据处理系统
首先，简易烈度计台站、强震台站和测震台站的数据实时汇集到地震数据汇集系统，再通过行业内网将数据实时传输至地震预警数据处理系统，数据传输流程如图 5所示。当地震发生时，地震预警数据处理系统利用多个台站收到的地震初至波信号，使用STA/LTA算法检查震相触发，使用AIC算法精确拾取震相到时，使用“着未着”算法快速分析多个台站震相到时数据，确定地震发生的时刻和位置，使用P波最大振幅和台站震中距估算地震震级，从而获得地震基本参数，并向地震波尚未到达的地区发布准确的地震预警信号。
图 5 地震预警试验区数据传输流程Figure 5. Data transfer process of simple intensity meter in early warning experimental region2018年8月16日13时50分20.97秒，天津市蓟州区发生2.2级地震，距离震中较近的简易烈度计台站、强震实时台站和测震台站均监测到了P波信号。地震预警数据处理系统采用256字节的数据打包模式进行数据传输，地震发生后8s发出第一报，发震时刻为2018年8月16日13时50分20.7秒，震中地点为天津蓟州区，震级为2.0级（转换震级），震中地点距天津市中心约94km，横波预计到达时间为12s，触发台站39个（图 6）。由于此次地震的震级较小，尚需经历更大的地震对天津简易烈度计地震预警试验区的效能进行检验。
4. 简易烈度计数据记录
简易烈度计观测的数据存储在地震数据处理中心服务器上，采用定长512字节MiniSeed格式，对于任一时段的数据，可利用JOPENS6.0数据处理系统转换为SEED、ASCII等多种数据文件。在数据处理过程中，采用python语言编写相关的计算程序，主要使用了python中的Numpy、Matplotlib、Basemap、Obspy等扩展包，有效提高了系统的开发效率和数据处理的工作效率。
在2018年8月16日天津蓟州2.2级地震中，位于天津蓟州附近的简易烈度计台站A0055、A0056、A0058、A0063均获取了观测数据。根据《仪器地震烈度计算（征求意见稿）》1的要求，对0.1—10Hz频率范围内的观测数据进行带通滤波处理，然后三分向合成计算最大峰值加速度。计算得到距离震中4km的A0058台站的最大峰值加速度为7.4gal；距离震中10km的A0055台站的最大峰值加速度为6.46gal；距离震中10km的A0063台站的最大峰值加速度为6.78gal；距离震中15km的A0056台站的最大峰值加速度为6.27gal。
1 中国地震局，2015.仪器地震烈度计算暂行规程.
由于简易烈度计台站没有同址的测震或强震台站，故选择台站距离相邻、方位相差较小的强震台站BAJ与简易烈度计台站A0056进行波形数据对比。BAJ台距离震中12km，三分向合成的最大峰值加速度为6.89gal，使用的地震计为SLJ-100型加速度计，测量的频带范围为0.1—30Hz，采样率为100Hz。简易烈度计台站A0056使用GL-P2B型MEMES加速度计，测量频带范围为0.1—30Hz，采样率为100Hz。观测台站分布如图 7所示。
图 7 蓟州2.2级地震震中附近的地震观测台站Figure 7. Location of seismic observation stations near the epicenter of the Jizhou MS 2.2 earthquake由于观测频带不同，采用4阶巴特沃斯（Buttorworth）带通滤波器对强震台BAJ与简易烈度计台A0056的0.1—20Hz观测数据进行滤波，结果如图 8所示。强震台BAJ三分向记录的波形峰值加速度（PGA）分别为3.99gal、3.88gal和1.43gal，简易烈度计台A0056三分向记录的波形峰值加速度（PGA）分别为1.23gal、1.28gal和0.78gal。2个台站的峰值加速度和波形数据时程均不相同，这是由于台站位置不同造成的，其中强震台BAJ距离震中12km，简易烈度计台A0056距离震中15km，二者直线距离相差3km，且相对震中的方位有差异，但二者NS向与EW向波形存在一定的相似性。
图 8 强震台站BAJ（a）与烈度计台站A0056（b）原始记录波形对比Figure 8. Comparison of original records of the strong ground motion station BAJ (a) to that of the simple intensity meter station A0056 (b)观测仪器的记录噪声水平在一定程度上能够直接反应仪器的记录能力，尤其是记录中小地震事件的能力。为此，截取了简易烈度计台站A0056无明显干扰的10分钟的数据记录，利用Welch算法计算加速度噪声功率谱（Peterson，1993），结果如图 9所示。由图可以看出，简易烈度计台A0056的加速度噪声功率谱处于高噪声模型（New High Noise Model）和低噪声模型（New Low Noise Model）之上，且远高于高噪声模型。由此可知，相对于测震台站对地脉动噪声的完整记录（在NHNM和NLNM之间），简易烈度计台站记录的噪声数据完全是仪器的自噪声，天然噪声则淹没在仪器的自噪声中无法分辨。由于简易烈度计的仪器自噪声已经超出地脉动水平，无法观测记录到有效地脉动，且仪器通常布设在人员比较密集的建筑物内，因此，简易烈度计台站的记录数据信噪比较低，地震数据初至震相会淹没在背景噪声中。
图 9 简易烈度计台站的噪声功率谱密度Figure 9. Power spectra density based on seismic noise data of simple intensity meter各类地震观测仪器的观测特点不同，测震台站主要服务于地球物理学的研究，旨在了解地球内部的构造和地震活动性，为此提供各种地震参数，包括地震基本参数、地震活动性、地震波传播、地震趋势等，地震观测的场地条件要求为基岩、深井且背景噪声应尽量小。强震观测主要为地震工程学研究和结构抗震提供资料，旨在测定地面和建筑物在地震作用下的运动过程，包括地震影响场、地面及结构的地震动特征，主要针对的是近场强震，仪器频带范围为-2—2g，通常要求背景噪声小于1cm/s2，城区或乡村均可布设观测仪器。而简易烈度计是以MEMS传感器为主的非专业强震仪，主要针对的是近场强震，要求背景噪声小于2cm/s2，记录质量较低，一般布设在自由地表或2层以下低矮结构建筑物内1。简易烈度计具有成本低、安装要求不高、维护简单等特点，适于高密度布设，以提高区域地震预警和烈度速报能力。
1 马强，2017.地震预警与烈度速报系统的试验示范.浙江省地震局内部报告.
5. 讨论与结论
本文对天津简易烈度计地震预警试验区的建设、台站分布情况、仪器性能指标、数据处理及数据记录情况进行了介绍，探讨了简易烈度计在预警台网中的应用，获得了以下初步认识：
（1）通过建设天津简易烈度计地震预警试验区，加密了天津强震台网的密度，初步实现了天津及邻区的地震预警功能。
（2）由测震台站、强震动实时传输台站及简易烈度计台站组成的天津地震预警观测网络，拓宽了各台网的数据应用范围，不仅可应用于地震应急快速响应，还可应用于地震工程和科学研究。
（3）简易烈度计台站造价低、建设周期短，易于大面积布设或规模化应用，可用较少的投资实现更好的应用（姚会琴等，2012），是现有测震台站和强震台站的有益补充。
（4）简易烈度计记录的噪声数据完全是仪器的自噪声，高于高噪声模型，但对于震中附近满足一定信噪比条件的台站，可获得较准确的初至震相到时信息，对于地震预警信息的准确产出具有意义。
（5）今后还需进一步积累地震数据，结合天津及邻区的地震分布特点优化地震预警参数模型，以期提升该地区的防震减灾能力以及公共服务水平。
-
图 1 有限元模型示意图及各土层动剪切模量比和阻尼比随剪应变变化图
Figure 1. Schematic diagram of the finite element model and the relation curves of dynamic shear module ratio and damp ration shear strain for soils
图 2 近场地震记录的水平地震动、竖向地震动加速度时程和水平加速度反应谱
Figure 2. Acceleration time histories of the horizontal and vertical ground motions and horizontal acceleration response spectra of the near-field seismic records
图 3 Loma地震记录双向地震动作用下不同隧道间距的附加内力
Figure 3. Additional internal forces of tunnel linings with different tunnel spacings under the loading of Loma seismic record
图 4 不同的土-结接触面摩擦系数下的隧道最大附加内力变化曲线（Loma地震记录作用下）
Figure 4. The variation curves of maximum additional internal forces of the tunnel with different friction coefficients of the contact surface between the soil and tunnel(Under the loading of Loma seismic record)
图 5 隧道最大附加内力与隧道刚度的关系（Loma地震记录作用下）
Figure 5. The relationship between the maximum additional internal forces of the lining and the lining stiffness(Under the loading of Loma seismic record)
图 6 不同近场地震记录作用下隧道的附加弯矩
Figure 6. The additional bending moments of the lining under the loadings of different near-field seismic records
图 7 不同近场地震记录作用下隧道的附加轴力
Figure 7. The additional axial forces of the lining under the loading of different near-field seismic records
图 8 土层加速度放大系数和标准加速度反应谱分析
Figure 8. The analysis of soil acceleration amplification coefficient and standard acceleration response spectrum
图 9 不同近场地震记录的水平地震动作用下，土层自由场各深度处的水平加速度放大系数
Figure 9. The horizontal acceleration amplification coefficients of the free field under the loadings of the horizontal ground motions of different near-field seismic records
图 10 隧道在Wolong地震记录的不同竖向地震动相对强度作用下的附加内力
Figure 10. Additional internal forces of the lining with different V/H ratios of the WoLong seismic inputs
图 11 Wolong地震记录作用下不同的竖向地震动相对强度下隧道的最大附加内力和土层自由场最大平均剪应变
Figure 11. Variation trends of maximum additional internal forces of the lining and maximum average shear strain of the free field under loading of Wolong seismic record
图 12 不同近场地震记录下隧道在单向和双向地震荷载作用下最大附加内力
Figure 12. The maximum additional internal forces of the lining under unidirectional and bidirectional loads resulted by different near-field earthquake records
表 1 土层的相应参数
Table 1. Corresponding parameters of soil layers
土层类别 厚度/m 初始剪切模量/MPa 密度/×103 kg·m-3 A 2B r0 /×10-4 泊松比 黏性土 18 50.67 1.94 1.12 0.82 9.00 0.35 砂质黏性土 4 100.84 1.98 1.00 0.76 7.38 0.32 中粗砂 5 149.73 1.98 1.00 0.70 9.67 0.30 砂层 13 300.84 1.98 1.10 0.74 9.50 0.30 砾质黏性土 16 361.57 1.96 1.20 0.74 10.03 0.30 
下载: 导出CSV
表 2 隧道管片的相应参数
Table 2. Corresponding parameters of the tunnel segments
杨氏模量/GPa 泊松比 密度/kg·m-3 外径/m 厚度/m 35.5 0.2 2600 6.2 0.35
下载: 导出CSV
表 3 近场地震记录加速度幅值
Table 3. Acceleration amplitudes of the near-field seismic records
近场地震记录 水平向加速度峰值/g 竖直向加速度峰值/g 竖直/水平 Arias intensity/m·s-1 Loma 0.2 0.106 0.53 0.32 ChiChi 0.2 0.132 0.66 0.88 Mammoth 0.2 0.056 0.28 0.09 WoLong 0.2 0.200 1.00 0.54
下载: 导出CSV
-
陈国兴, 陈苏, 杜修力等, 2016. 城市地下结构抗震研究进展. 防灾减灾工程学报, 36(1): 1-23. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201601002.htm 崔光耀, 王明年, 林国进等, 2011. 汶川地震区典型公路隧道衬砌震害类型统计分析. 中国地质灾害与防治学报, 22(1): 122-127. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGDH201101022.htm 杜修力, 李洋, 许成顺等, 2018a. 1995年日本阪神地震大开地铁车站震害原因及成灾机理分析研究进展. 岩土工程学报, 40(2): 223-236. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201802003.htm 杜修力, 许紫刚, 袁雪纯等, 2018b. 地震动峰值位移和峰值速度对地下结构地震反应的影响. 震灾防御技术, 13(2): 293-303. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20180205&journal_id=zzfyjs 黄宏伟, 徐凌, 严佳梁等, 2006. 盾构隧道横向刚度有效率研究. 岩土工程学报, 28(1): 11-18. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200601001.htm 刘如山, 胡少卿, 石宏彬, 2007. 地下结构抗震计算中拟静力法的地震荷载施加方法研究. 岩土工程学报, 29(2): 237-242. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200702016.htm 楼梦麟, 陈清军, 1999. 侧向边界对桩基地震反应影响的研究. 上海: 同济大学出版社. 苗雨, 孙甜粲, 李威, 2014. 基于多尺度几何有限元模型分析的地铁盾构隧道地震反应研究. 地震工程与工程振动, 34(S1): 915-918. https://www.cnki.com.cn/Article/CJFDTOTAL-DGGC2014S1145.htm 赵丁凤, 阮滨, 陈国兴等, 2017. 基于Davidenkov骨架曲线模型的修正不规则加卸载准则与等效剪应变算法及其验证. 岩土工程学报, 39(5): 888-895. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201705018.htm 赵武胜, 何先志, 陈卫忠等, 2012. 盾构隧道地震响应分析方法及工程应用. 岩土力学, 33(8): 2415-2421. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201208029.htm 中华人民共和国住房和城乡建设部, 2014. GB 50909-2014城市轨道交通结构抗震设计规范. 北京: 中国标准出版社. 朱育才, 2012. 近断层脉冲型地震动作用下隧道抗减震研究. 重庆: 重庆交通大学. Abate G., Massimino M. R., 2017. Numerical modelling of the seismic response of a tunnel-soil-aboveground building system in Catania (Italy). Bulletin of Earthquake Engineering, 15(1): 469-491. Balendra T., Thambiratnam D. P., Koh C. G., et al., 1984. Dynamic response of twin circular tunnels due to incident SH-waves. Earthquake Engineering and Structural Dynamics, 12(2): 181-201. doi: 10.1002/eqe.4290120204 Cilingir U., Madabhushi S P G., 2011. A model study on the effects of input motion on the seismic behaviour of tunnels. Soil Dynamics and Earthquake Engineering, 31(3): 452-462. http://www.sciencedirect.com/science/article/pii/S026772611000240X Kouretzis G. P., Sloan S. W., Carter J. P., 2013. Effect of interface friction on tunnel liner internal forces due to seismic S- and P-wave propagation. Soil Dynamics and Earthquake Engineering, 46: 41-51. http://www.sciencedirect.com/science/article/pii/S0267726112002692 Law H K., Lam I P., 2001. Application of periodic boundary for large pile group. Journal of Geotechnical and Geoenvironmental Engineering, 127(10): 889-892. http://www.researchgate.net/publication/245293547_Application_of_Periodic_Boundary_for_Large_Pile_Group Lee K. M., Ge X. M., 2001. The equivalence of a jointed shield-driven tunnel lining to a continuous ring structure. Canadian Geotechnical Journal, 38(3): 461-483. doi: 10.1139/t00-107 Lu J. C., Elgamal A., Yan L. J., et al., 2011. Large-scale numerical modeling in geotechnical earthquake engineering. International Journal of Geomechanics, 11(6): 490-503. http://www.zhangqiaokeyan.com/academic-journal-foreign_other_thesis/020416988579.html Penzien J., Wu C. L., 1998. Stresses in linings of bored tunnels. Earthquake Engineering Structural Dynamics, 27(3): 283-300. doi: 10.1002/(SICI)1096-9845(199803)27:3<283::AID-EQE732>3.0.CO;2-T/pdf Sedarat H., Kozak A., Hashash Y. M. A., et al., 2009. Contact interface in seismic analysis of circular tunnels. Tunnelling and Underground Space Technology, 24(4): 482-490. http://www.sciencedirect.com/science/article/pii/S0886779808001132 Sharma S., Judd W. R., 1991. Underground opening damage from earthquakes. Engineering Geology, 30(3-4): 263-276. http://www.onacademic.com/detail/journal_1000034559003210_4d4d.html Wang J. N., 1993. Seismic design of tunnels: a state-of-the-art approach, monograph, monograph 7, parsons. New York: Brinckerhoff, Quade and Douglas Inc. 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. http://www.sciencedirect.com/science/article/pii/S0886779801000475 Zhang L., Liu H. B., 2017. Seismic response of clay-pile-raft-superstructure systems subjected to far-field ground motions. Soil Dynamics and Earthquake Engineering, 101: 209-224. http://www.sciencedirect.com/science/article/pii/S0267726117302956 -
下载:
点击查看大图
计量
- 文章访问数: 151
- HTML全文浏览量: 68
- PDF下载量: 5
- 被引次数: 0
