Partitioning method for seismic impact zones of tunnels crossing active faults based on energy principle and its applications
-
摘要: 为提高穿越活动断层隧道的抗震性能，从结构的能量响应出发，引入了断层错动因素，优化了地震作用下穿越活动断层隧道的能量响应计算理论，根据输入能量与结构固有能量之间的关系，提出了地震和断层错动作用下隧道能量响应曲线计算方法，建立了穿越活动断层隧道震害影响区分区方法，依托龙洞子隧道F5活动断层段进行应用研究。结果表明：F5活动断层段隧道在地震10 s左右发生失稳破坏，破坏率最大为58%，根据破坏率对穿越断层隧道进行了震害影响分区；通过汶川地震龙洞子隧道F5和龙溪隧道F8活动断层段现场震害资料分析表明，现场震害类型及震害分布情况与震害影响分区较为符合，这说明提出的基于能量原理的穿越活动断层隧道震害影响分区方法可用于工程实际。研究成果可为穿越活动断层隧道抗震设防设计提供参考。Abstract: In order to improve the seismic performance of tunnels crossing active faults, starting from the energy response of structures, the factor of fault misalignment is introduced, the theory for calculating the energy response of tunnels crossing active faults under the action of earthquakes is optimized, and the method for calculating the energy response curve of the tunnels under the action of earthquakes and fault misalignment is put forward according to the relationship between the input energy and the intrinsic energy of the structures. Thus the method for distinguishing the zone of seismic damage impacts in the tunnels crossing active faults is established, and it is applied in the researches on the active fault section of F5 of the Longdongzi Tunnel. The results show that the destabilisation damage of the F5 active fault segment occurs around 10 s after the earthquake, with the maximum damage rate of 58%. The seismic impact zones for the tunnels crossing active faults are partitioned based on the damage rate, Through the analysis of the seismic damage information of the F5 active fault section of Longdongzi Tunnel and the F8 active fault of Longxi Tunnel set in Wenchuan Earthquake, the type and distribution of seismic damages on the site are more consistent with the proposed seismic impact zones. It is shown that the seismic impact zone method for the tunnels crossing active faults based on the energy principle can be used in engineering practice. The research results can provide an reference for the anti-seismic design of tunnels crossing active faults.
-
Keywords:
- tunnel engineering /
- energy principle /
- fault /
- impact zone /
- damage rate
-
0. 引言
北疆供水渠道在进行渠道设计时，不仅仅要考虑经济最优和保证流量最优的问题，还需要考虑渠道的冬季冻胀破坏。梯形渠道占中国输水渠道的比例较大，应用较广。因此，梯形渠道的冻胀问题受到了北方寒区广泛关注。王正中等[1-2]对梯形渠道冻胀模型进行了理论分析，将梯形渠道考虑为在切向冻结力、法向冻胀力的以及衬砌板约束下的简支梁，王正中等[3]、刘旭东等[4]、李爽等[5]用数值模拟的方法模拟梯形断面形式下冻胀的规律，得到了温度场，位移场沿着渠道分布情况。
本文主要研究不同梯形断面情况下，以北疆某供水工程总干渠退水渠段梯形渠道为例，对该渠道用comsol进行数值模拟，对不同的断面渠道进行参数化模拟，得到最佳参数范围，并结合水力最佳断面，得到双优断面形式，为寒区渠道设计提供参考。
1. 力学模型的建立
1.1 基本假设
渠基土冻结时，土体、水和冰之间的相互作用的微观结构及动态过程相当复杂，目前很难准确模拟整个冻胀过程。为了便于分析，对其进行适当的简化，以便抓住影响冻结过程中及冻胀变形的主要特征，主要假设如下[6]：①根据现场及室内试验研究，假设冻土是均匀连续各向同性体；②尽管土的冻胀与其温度、水分、土质密切相关，对具体工程当水分及土壤条件确定时土体最终冻胀主要取决于温度，将水分迁移对冻土的体积影响，以线膨胀系数表达；③根据试验研究假定相变温度在同一种土中和同种外力条件下为常值，即暂取相变温度为0℃[7]；④由于渠道为细长结构，不考虑长方向上土颗粒对温度的影响，选取平面应变问题进行模拟。
1.2 热传导方程
根据以上假设，将温度热传递视为二维瞬态热传递的过程，建立二维热传导控制方程：
(λx∂2T∂x2)+(λy∂2T∂y2)=ρc∂T∂t。 (1) 式中T为温度；λx，λy分别是冻土沿x，y方向的导热系数；ρ为土体的密度；c为土体的比热容；t为时间。
1.3 本构方程
冻土属于冷胀热缩材料，冻土在冻结过程中水冻结成冰，除了原位水冻结体积膨胀，还有从未冻结区向冻结区迁移的水分冻结成冰，本文将原位水以及迁移水冻结成冰的体积膨胀以关于温度T的函数的线膨胀系数表示。因此，土体的本构方程可以表示为
ϵx=σxE−μσyE+α(T−Tref) ，ϵy=σyE−μσxE+α(T−Tref) ，γxy=2(1+μ)Eτxy 。} (2) 式中εx，εy为正应变；γxy为剪应变；σx、σy为正应力；τxy为剪应力；E为弹性模量；μ为泊松比，α为混凝土或者冻土自由冻胀时的热膨胀系数；T为计算温度；Tref为参考温度。
2. 计算模型
2.1 原型概况
本文选择北疆地区某输水渠道退水渠段的梯形渠道（图 1）为研究对象，原型渠道的基本情况如表 1，2所示，并对该模型渠道进行了冻胀模拟分析计算，因原型渠道各部位坡向，土质，水分不同，在模型边界设置时对阴阳破和渠堤分别对温度和土体冻胀率赋值。
表 1 渠道各部位的表面温度以及冻结期Table 1. Surface temperatures and freezing periods of various parts of canal部位 月平均表面温度/(℃) 冻结期、（月-日） 12月 1月 2月 阴坡 -14.92 -18.85 -10.72 11-27—02-27 渠底 -14.56 -16.22 -9.15 11-27—02-26 阳坡 -12.55 -14.75 -10.54 11-27—02-27 表 2 原型渠道基本情况Table 2. Basic information of prototype canal部位 渠床
土质基土干密度/(g·cm-3) 冻深h/cm 冻胀量Δh/cm 冻胀率η/% 阴坡 砂砾石 1.80 171 5.0 2.92 渠底 159 4.4 2.77 阳坡 146 3.7 2.53 2.2 有限元模型
有限元计算模型为：①原型梯形断面渠道冻胀数值模拟；②梯形断面渠道冻胀“参数化”断面渠道数值模拟。
有限元模型如图 2所示，选取各表面温度接近原型渠道，渠道两边设置为热绝缘边界，下边界为固定温度边界取10℃，阴坡边界温度取冻结期平均值-14.8℃，渠底温度取-13.3℃，阳坡温度为-12.6℃。几何模型如图 2所示，渠底长2.0 m，渠坡的横向长度为3.75 m，坡比为1.5，基础底向下取2.5 m，左右边界取0.75 m。模型将衬砌板与冻土作为一个整体进行数值模拟，力学的边界设置时，左右边界设置为辊支撑，底部设置指定位移0。
2.3 参数取值
参考文献[5]中线膨胀系数按照η/Tmin取值，η为冻胀率，Tmin为相应部位月平均表面温度的最小值。由式（2），将T-Tref看成1℃，那么可以将线膨胀系数α考虑跟温度Tmin相关的函数，数值上与冻胀率相等。未冻土与混凝土的导热系数参考文献[3]，冻土的导热系数与土体的未冻水含量相关，冻土融土导热系数比与未冻水含量相关[9]，模型渠道的未冻水含量选取19.4%，因此冻土的导热系数为未冻土导热系数的0.9倍。
3. 计算结果分析
3.1 温度场
图 3可以看出：接近于渠底的温度分布接近于一组平行的直线，0℃等温线以上，由于未冻水冻结导致剧烈相变，该区域导热系数小，因此等温线较密集，温度梯度比较小；0℃等温线以下，导热系数较大，等温线较稀疏。阴坡冻深为162 cm，渠底冻深为146 cm，阳坡的冻深为141 cm。与表 3模型试验较符合，最大相对误差为8.2%。
表 3 材料热力学参数Table 3. Thermodynamic parameters of materials介质 弹性模量E/Pa 泊松比ν 导热系数λ/(W·m-1·℃) 线膨胀系数α/℃-1 混凝土 2.4×1010 0.167 1.580 1.1×10-5 冻土 4.6×107 0.330 1.188[9] 阴坡2.92%、渠底2.77%、阳坡2.53% 未冻土 1.5×107 0.375 1.320 0 3.2 位移场
如图 4所示的是模拟得到的结果与试验结果进行对比，模型试验结果冻胀量边坡较大，渠底较小。模型试验的冻胀量整体比模拟冻胀量大，因为模型试验所用的衬砌板受到边坡的约束较大，且衬砌板厚度小于模拟设定厚度。模拟得到渠底冻胀量比较平缓而试验得到的渠底中间部位冻胀量大而两边小，模拟设定的衬砌板是整体现浇型板，在冻胀力的作用下容易产生整体变形；而模型试验的衬砌板分布较为离散。
3.3 “参数化”位移场模拟结果
根据梯形断面的设计规范要求，应满足：①渠道边坡系数应不小于允许最小边坡系数，应不大于允许最大边坡系数，1≤m≤2。②渠道的宽深比应满足规范给定的要求，0.8 < b/h < 3.5。
根据要求，“参数化”分析模拟了宽深比β一定时，不同边坡比的渠道断面冻胀量情况；和边坡比一定时，不同宽深比的渠道断面冻胀量情况。
由图 5，宽深比一定时，边坡比越小，渠坡的法向冻胀量越大；m=2时，左侧阴坡法向冻胀量最大为5.2 cm，阳坡为1.7 cm。m=1时，阴坡法向冻胀量最大值为3.3 cm。由图 6，边坡比一定时，宽深比对法向冻胀量影响不大。
![]()
图 5 不同坡比法向冻胀量比较Figure 5. Comparison of normal frost-heave quantity under different slopes![]()
图 6 不同宽深比法向冻胀量比较Figure 6. Comparison of normal frost-heave quantity under different breadth depth ratios4. 水力最佳断面
从设计角度考虑水力最佳断面通常是指过水流量一定，所要求的过水断面面积最小，也就是所需要的材料最少，施工且最方便[12]。或者是过水断面面积一定，所通过的流量最大。
考虑最小冻胀量，所以采用边坡比为2的冻胀断面。水力学中考虑流量一定，过水断面面积最小可采用以下公式计算宽深比[13]：
β=2（√1+m2−1)。 (3) 将相关参数代入式（3），可得到考虑水力冻胀最优宽深比为2.472。
5. 结语
本文利用comsol有限元软件，对北疆供水渠道以及渠基土冻胀破坏规律进行有限元分析，与试验数据对比，基本与试验数据相似。对不同坡比以及不同宽深比的渠道进行有限元模拟分析，发现坡比影响冻胀量的大小，坡比越大，渠道冻胀量越小。理论分析受力情况易得坡比越大越能释放变形，符合实际。再结合水力最佳断面的公式，可以得到最佳宽深比。最终的断面形式能有效的防止冻胀，且满足水力断面要求。为北疆地区渠道防冻胀断面设计提供了理论支撑。
-
![]()
图 6 隧道纵向破坏率及震害影响分区图
Figure 6. Longitudinal damage rate and seismic impact zones of tunnel
表 1 计算参数
Table 1 Computational parameters
材料 重度/
(kN·m-3)弹性模量/
GPa泊松比 内摩擦角/
(°)黏聚力/MPa 上下盘Ⅳ级围岩 22 2.5 0.3 35 0.5 破碎带Ⅴ级围岩 20 1.0 0.4 25 0.2 Ⅱ级围岩(基岩) 25 20 0.2 50 1.5 初支 22 23 0.2 — — 二衬 25 28 0.2 — — 
下载: 导出CSV
表 2 监测点能量峰值
Table 2 Peak energy values at monitoring points
(单位：kJ) 能量 拱顶 左拱肩 左边墙 左拱脚 仰拱 结构动能 28.7 20.5 26.1 27.3 18.1 地震输入能 28.2 26.1 20.1 25.1 27.3 错动输入能 7.05 6.52 5.03 6.29 6.82 阻尼能 2.41 2.31 2.04 2.27 2.13 应变能 6.51 5.94 2.21 3.85 6.26
下载: 导出CSV
表 3 监测点破坏率
Table 3 Damage rates at monitoring points
监测点位置 拱顶 左拱肩 左边墙 左拱脚 仰拱 破坏率 56.67% 53.0% 26.33% 50.33% 43.66%
下载: 导出CSV
表 4 震害影响区分区标准
Table 4 Criteria for seismic impact zones
分区 强影响区 中影响区 弱影响区 无影响区 破坏率 ∂≥45% 45%＞∂≥30% 30%＞∂≥10% 10%＞∂
下载: 导出CSV
表 5 震害统计
Table 5 Statistics of seismic damage
震害
类型震害长度/m 震害
比例震害
类型震害长度/m 震害
比例A 155 30.39% E 68 13.33% B 159 31.18% F 12 2.35% C 47 9.22% G 207 40.59% D 69 13.52% — — — 注：G为震害总长度；震害比例为震害长度比统计隧道长度，F5断层附近的统计隧道长度为510 m。
下载: 导出CSV
-
[1] 公路工程抗震规范: JTG B02—2013[S]. 北京: 人民交通出版社, 2014. Specification of Seismic Design for Highway Engineering: JTG B02—2013[S]. Beijing: China Communications Press, 2014. (in Chinese)
[2] 铁路工程抗震设计规范: GB 50111—2006[S]. 北京: 中国计划出版社, 2006. Code for Seismic Design of Railway Engineering: GB 50111—2006[S]. Beijing: China Planning Press, 2006. (in Chinese)
[3] 刘学增, 唐精, 桑运龙, 等. 逆断层黏滑错动对跨断层隧道影响机制的模型试验研究[J]. 隧道建设(中英文), 2020(4): 481-489. LIU Xuezeng, TANG Jing, SANG Yunlong, et al. Model experimental study on influencing mechanism of reverse fault stick-slip dislocation on cross-fault tunnel[J]. Tunnel Construction, 2020(4): 481-489. (in Chinese)
[4] 武世燕. 高烈度地震区铁路隧道震害特征与整治研究[J]. 铁道工程学报, 2023, 40(1): 76-83. WU Shiyan. Research on the damage characteristics and countermeasures of railway tunnel in high-intensity earthquake area[J]. Journal of Railway Engineering Society, 2023, 40(1): 76-83. (in Chinese)
[5] 汪振, 钟紫蓝, 赵密, 等. 正断型断裂模拟及其对山岭隧道影响研究[J]. 岩土工程学报, 2020, 42(10): 1876-1884. WANG Zhen, ZHONG Zilan, ZHAO Mi, et al. Simulation of normal fault rupture and its impact on mountain tunnels[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(10): 1876-1884. (in Chinese)
[6] 唐浪洲, 于丽, 王玉锁, 等. 走滑断层错动量大小对铁路隧道结构安全性影响的数值分析[J]. 现代隧道技术, 2022, 59(1): 214-224. TANG Langzhou, YU Li, WANG Yusuo, et al. Numerical analysis on the effect of strike-slip fault dislocation on the structural safety of railway tunnels[J]. Modern Tunnelling Technology, 2022, 59(1): 214-224. (in Chinese)
[7] 何川, 李林, 张景, 等. 隧道穿越断层破碎带震害机理研究[J]. 岩土工程学报, 2014, 36(3): 427-434. HE Chuan, LI Lin, ZHANG Jing, et al. Seismic damage mechanism of tunnels through fault zones[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(3): 427-434. (in Chinese)
[8] 崔光耀, 王明年, 于丽, 等. 汶川地震公路隧道洞口结构震害分析及震害机理研究[J]. 岩土工程学报, 2013, 35(6): 1084-1091. CUI Guangyao, WANG Mingnian, YU Li, et al. Seismic damage and mechanism of portal structure of highway tunnels in Wenchuan Earthquake[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(6): 1084-1091. (in Chinese)
[9] 田四明, 吴克非, 于丽, 等. 穿越活动断裂带铁路隧道抗震关键技术[J]. 隧道建设(中英文), 2022, 42(8): 1351-1364. TIAN Siming, WU Kefei, YU Li, et al. Key technology of anti-seismic for railway tunnels crossing active fault zone[J]. Tunnel Construction, 2022, 42(8): 1351-1364. (in Chinese)
[10] 周栋, 张成军, 严武建, 等. 2022年青海门源6.9级地震发震构造及甘肃灾区震害特征研究[J]. 世界地震工程, 2023, 39(1): 20-27. ZHOU Dong, ZHANG Chengjun, YAN Wujian, et al. Study on the seismogenic mechanism of the 2022 Menyuan MS6.9 and the seismic damage characteristics in Gansu Province[J]. World Earthquake Engineering, 2023, 39(1): 20-27. (in Chinese)
[11] 尹超. 穿越走滑断层隧道破坏特性及结构抗震性能研究[D]. 成都: 西南交通大学, 2022. YIN Chao. Study on Failure Characteristics and Seismic Performance of Tunnel Crossing Strike-Slip Fault[D]. Chengdu: Southwest Jiaotong University, 2022. (in Chinese)
[12] 张锴. 穿越断层带多破裂面隧道地震动力响应及设防分区研究[J]. 铁道标准设计, 2025, 69(4): 119-128, 138. ZHANG Kai. Research on seismic dynamics response and fortification zoning of tunnels crossing multiple rupture surfaces in fault zones[J]. Railway Standard Design, 2025, 69(4): 119-128, 138. (in Chinese)
[13] 孙风伯. 穿越活动断裂山岭隧道抗位错机理与方法研究[D]. 北京: 北京交通大学, 2018. SUN Fengbo. Study on Dislocation Resistance Mechanism and Method of Mountain Tunnel Crossing Active Fault[D]. Beijing: Beijing Jiaotong University, 2018. (in Chinese)
[14] 崔光耀, 宋博涵, 王明年, 等. 基于能量守恒原理的跨活动断层隧道抗错断设计方法研究[J]. 土木工程学报, 2020, 53(增刊2): 309-314. CUI Guangyao, SONG Bohan, WANG Mingnian, et al. Study on the anti-breaking design method of tunnel across active fault based on the energy balance method[J]. China Civil Engineering Journal, 2020, 53(S2): 309-314. (in Chinese)
[15] 曲宏略, 罗浩, 刘辉, 等. 跨断层隧道震害特性的能量分析方法[J]. 铁道工程学报, 2017, 34(3): 58-62, 91. QU Honglüe, LUO Hao, LIU Hui, et al. An energy analysis approach for seismic characteristics of crossing-fault tunnel[J]. Journal of Railway Engineering Society, 2017, 34(3): 58-62, 91. (in Chinese)
[16] 马至刚, 朱宝龙, 吴红刚, 等. 基于能量分析滑坡-隧道减震优化地震破坏特性试验研究[J]. 岩石力学与工程学报, 2023, 42(4): 879-895. MA Zhigang, ZHU Baolong, WU Honggang, et al. Experimental study on seismic failure characteristics optimization of landslide-tunnel shock absorption based on energy analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2023, 42(4): 879-895. (in Chinese)
[17] ZHANG X, SHEN Y S, ZHOU P F, et al. Partitions of seismic affected zone and characteristics of segmental lining based on energy principle[J]. Journal of Central South University, 2023, 30(1): 227-242.
计量
- 文章访问数: 159
- HTML全文浏览量: 20
- PDF下载量: 38
