大同-阳高地震重力异常f

华昌才, 郑金涵, 刘端法

华昌才, 郑金涵, 刘端法. 1992: 大同-阳高地震重力异常f. 地震学报, 14(3): 369-372.
引用本文: 华昌才, 郑金涵, 刘端法. 1992: 大同-阳高地震重力异常f. 地震学报, 14(3): 369-372.

大同-阳高地震重力异常f

  • 摘要: 1989年10月18日,在大同-阳高之间发生了中强震群:一个6.1级及六个5级地震,震中区约在东经11350',北纬3953'.这次地震在极震区引起了Ⅶ度的破坏,受灾较重的地区为册田(Ⅷ度).在震前的1982年,我们曾在距震中约10km 的西坪地震台实施过流动重力测量;地震发生之后,在国家地震局科技司的支持下立即进行复测.这两次测量,使我们可以比较地震区重力场是否发生过变化,并为震后的震情监视提供背景资料.
  • 黄土是一种第四纪沉积物,广泛分布于世界各个国家,其面积约占陆地总面积的9.3%.我国是黄土占地面积最大的国家,约占国土面积的6%.甘肃引洮供水一期工程中10#—15#隧洞为黄土隧洞,场地黄土分为两层,上部为上更新统风积马兰黄土(eolQ23),下部为冲洪积粉质壤土(al-plQ13),两者的工程地质性状基本相似,分析时统称为Q3黄土.黄土材料在干燥状态下具有较高的强度,黄土隧洞也具有较好的稳定性,然而黄土遇水后强度降低,从而导致隧洞结构破坏而发生坍塌等事故.研究人员对黄土围岩地下结构在地震作用下的破坏情况进行了研究分析. 例如:郑颖人等(2010)为了得到地震作用下无衬砌隧洞的力学规律,利用有限差分程序,同时考虑隧洞的拉破坏和剪破坏,得出了隧洞顶部首先出现局部拉破坏,然后出现侧边整体破坏的破坏机制;陈立伟等(2007)将黄土视为理想弹塑性材料,采用Druke-Prager破坏准则,分析黄土暗穴在地震作用下的稳定性; 邓国华等(2007)基于已提出的Q3黄土初始切线模量与含水量的关系,对一典型的黄土隧洞进行增湿变形分析,得到了典型黄土隧洞在围岩增湿情况下围岩的变形、地面沉降和衬砌受力等结果.

    本文以前人研究为基础,结合黄土抗剪强度指标和初始弹性模量的变化,模拟了增湿条件下黄土隧洞在水平地震作用下,自振圆频率、质量阻尼系数、刚度阻尼系数和隧洞顶部地震加速度时程的变化规律,以及衬砌的薄弱位置及其主应力和最大位移变化等规律.

    隧洞横断面形状为类马蹄形,宽×高为5.5 m×6 m,衬砌厚度为0.5 m.为了减少计算时间,左右边界及底部边界计算范围分别取隧洞宽和高的5倍,上部边界至隧洞顶部为30 m,隧洞上部边界为自由边界,左右边界为竖向约束,底部为固定绞支座(高峰,任侠,2001程选生,郑颖人,2011).由于沿隧洞长度方向尺寸远比隧洞横向宽度大,故按照平面应变问题计算.计算模型如图 1所示.

    图  1  黄土隧洞计算模型
    Figure  1.  The calculation model of the loess tunnel

    引洮一期总干渠黄土主要分布在内官香泉盆地南部的黄土梁峁地貌区,10#—15#隧洞均有分布.该土层分为两层,上层为厚度5—25 m的马兰黄土(eolQ23),下层为粉质壤土(al-plQ13),两者的工程地质性状相似,在分析时统称为Q3黄土.该黄土液限为28%—31.6%,塑限为16.9%—22.1%,黏聚力为18.7—31.1 kPa,内摩擦角为22.7°—28.6°,压缩模量为5.56—14.93 MPa,弹性模量为36.6—142.3 MPa.

    工程中能使黄土含水量增加的原因主要为降雨和地下水位的上升,黄土增湿后的一个显著特点就是其抗剪强度指标发生巨大变化.随着含水量的增加,其黏聚力骤减,但内摩擦角的变化并不是很显著.Drucker-Prager屈服准则能准确描述岩石、土壤等材料,在有限元分析中可得到较为精确的结果.黄土的塑性行为被假定为理想弹塑性. 因此本文假设黄土为理想弹塑性材料,利用Drucker-Prager屈服准则,采用有效应力的变化来描述黄土强度的变化. 增湿黄土的有效黏聚力与含水量的关系为(张茂花等,2006)

    式中:c′为有效黏聚力,单位为kPa;w为含水量.

    弹性模量是黄土隧洞抗震分析最主要的参数之一,其变化与黄土损伤有很大的关系. 基于已提出的Q3黄土初始弹性模量与含水量的关系(林斌,2005),可以初步判断Q3黄土在不同含水量情况下的初始弹性模量为

    式中:E为Q3黄土初始弹性模量,单位为MPa;w为黄土初始含水量.

    根据已有理论公式,结合引洮隧洞黄土的特有性质,取黄土塑限范围内含水量分别为7.42%、9.36%、11.95%、13.25%和15.10%(对超出塑限含水量的黄土围岩隧洞作了相同的分析,其数值结果不收敛,说明黄土围岩含水量超过其塑限后的时程反应并不稳定而发生破坏)的黄土围岩隧洞进行地震时程分析.据陈存礼等(2006)研究结果,可以不考虑含水量对黄土阻尼比的影响.计算中假定阻尼比和泊松比均为常数,隧洞衬砌为C30混凝土,容量为2.5 kN/m3,弹性模量为30 GPa,泊松比为0.20.黄土围岩部分计算参数张茂花等(2006)结果(表 1),其中含水量9.36%和13.25%为其它含水量参数的线性内插值.

    表  1  黄土隧洞围岩计算参数
    Table  1.  The calculation parameters of tunnel in loess surrounding rock
    wγ/(kN·m-3)E/MPac/kPaφμξ
    7.42%1.6413057.824.10.30.12
    9.36%1.6111839.525.20.30.12
    11.95%1.5910031.026.00.30.12
    13.25%1.569720.226.50.30.12
    15.10%1.549317.227.10.30.12
    注: w为含水量(下同), γ为容量, E为初始弹性模量, c为黏聚力, φ为内摩擦角, μ为泊松比, ξ为阻尼比.
    下载: 导出CSV 
    | 显示表格

    在进行地震加速度时程分析时,应考虑场地特征和近、远震情况,合理选用输入模型的地震加速度.地下工程的抗震问题实际上是非常复杂的,应考虑基岩深度以及基岩面上的实际地震波等情况. 为简单起见,本文选取能代表Ⅱ,Ⅲ类场地的El-Centro地震波NS向量,其峰值加速度出现在2.12 s,为3.418 m/s2. 本文按照烈度为Ⅷ度的地震(峰值加速度为1.96 m/s2)进行计算,如图 2所示.对其原始数据波进行调整,使其峰值加速度为1.96 m/s2. 参考翁祝梅和毛丽华(2006)的研究结果,地下结构抗震验算只考虑水平向地震波即可.

    图  2  El-Centro地震波加速度示意图
    Figure  2.  Schematic diagram of El-Centro seismic wave acceleration

    地震时程分析过程为非线性求解,在水平地震作用下,黄土隧洞在(tt)时刻的动力平衡方程为

    式中:at+Δt,ått和ätt分别为黄土隧洞在(tt)时刻的位移向量、速度向量和加速度向量;MCK分别为黄土隧洞的质量矩阵、阻尼矩阵和刚度矩阵;(äg)tt为外部输入的黄土隧洞的地震加速度时程.采用Newmark隐式积分法,假设

    将式(4)、(5)带入式(3),简化后得

    由式(6)解出att后,代入式(4)解出ätt,最后代入式(5)解出ått

    黄土隧洞隔离体的总阻尼矩阵采用瑞雷阻尼矩阵,即

    式中,η为质量阻尼系数,λ为刚度阻尼系数.ηλ可分别由下式求解:

    式中:ξiξj分别为第i阶和第j阶对应的阻尼比. 为分析简单,假定阻尼比ξi=ξj=ξωiωj分别为由模态分析得到的第i阶和第j阶的自振圆频率.

    首先对黄土隧洞模型隔离体进行模态分析,分析其前8阶自振圆频率.增湿后各含水量所对应的圆频率如表 2所示. 随着含水量的增加,同一阶的圆频率呈减小的趋势,且减小的趋势在含水量w≥11.95%时变慢.自振周期与自振圆频率的关系为

    表  2  黄土隧洞围岩含水量与圆频率关系表
    Table  2.  The relationship between water content and circular frequency of tunnel in loess surrounding rock
    wω1ω2ω3ω4ω5ω6ω7ω8
    7.42%0.689981.55692.07022.46392.70833.31173.44903.4992
    9.36%0.663581.49741.99082.36952.60483.18513.31723.3648
    11.95%0.614901.38761.84462.19552.41392.95143.07413.1172
    13.25%0.611421.37981.83412.18312.40032.93483.05683.0994
    15.10%0.602591.35991.80762.15162.36572.89253.01283.0545
    注: ω1ω8分别为第1—8阶圆频率.
    下载: 导出CSV 
    | 显示表格

    可以看出,随着含水量的增加,自振圆频率ω的减小,自振周期T相应增大.

    表 2可以看出,ω1ω2变化最大,而ω6以后的圆频率变化则变小,故本文选取ω1ω6. 根据式(8)和式(9)求出增湿后各含水量下的质量阻尼系数η和刚度阻尼系数λ(表 3),然后根据式(7)计算瑞雷阻尼矩阵.

    表  3  黄土隧洞围岩含水量与阻尼系数关系表
    Table  3.  The relationship between water content and damping coefficients of tunnel in loess surrounding rock
    wηλ
    7.42%0.1370.060
    9.36%0.1320.062
    11.95%0.1220.067
    13.25%0.1210.068
    15.10%0.1200.069
    注: η为质量阻尼系数, λ为刚度阻尼系数.
    下载: 导出CSV 
    | 显示表格

    表 3可以看出,随着含水量的增加,质量阻尼系数减小,刚度阻尼系数变大,且均在含水量w≥11.95%后变化减缓.

    输入水平向El-Centro地震波加速度后,得到不同含水量情况下黄土隧洞顶部的地震峰值加速度,如表 4所示.可以看出:在Ⅷ度地震下,各含水量情况下黄土隧洞顶部峰值加速度发生时间均为2.22 s,相比输入加速度发生时间出现了滞后现象;在含水量为7.42%和9.36%情况下,地震加速度出现了放大现象; 而后3组含水量情况下,峰值加速度均比输入加速度小.根据表 4数据可以拟合出一条关于含水量w2多项式的加速度-含水量曲线,如图 3所示. 由图中曲线得出临界含水量为10.12%,即在此含水量下黄土隧洞顶部地震加速度不会出现增大或减小的现象.

    表  4  黄土隧洞顶部含水量与地震峰值加速度关系表
    Table  4.  The relationship between water content and seismic peak acceleration of the loess tunnel roof
    wt/sA/(m·s-2)
    7.42%2.222.029
    9.36%2.221.978
    11.95%2.221.943
    13.25%2.221.940
    15.10%2.221.933
    注: A为地震峰值加速度.
    下载: 导出CSV 
    | 显示表格
    图  3  黄土隧洞顶部地震峰值加速度随含水量的变化
    Figure  3.  Variation of seismic peak acceleration with water content of the loess tunnel roof

    根据所得的各含水量所对应的黄土隧洞顶部峰值加速度拟合曲线,可得到拟合方程为

    式中,w为黄土围岩含水量,A为地震峰值加速度.

    高波等(2009)研究发现,隧洞的拱肩处出现几乎贯通整条隧洞的纵向裂缝;由田瑞瑞(2012)研究结果可知,有衬砌隧洞的塑性区域出现在两拱肩和拱脚位置,隧洞的最大、最小主应力及最大位移均出现在拱肩和拱脚处,该结果与实际工程发生地震破坏时的情况和数值模拟结果相似.由此可以看出,隧洞的薄弱位置在隧洞左右拱肩及左右拱脚位置.表 5给出了各含水量情况下,左、右拱肩和拱脚处的主应力及最大位移.

    表  5  黄土隧洞衬砌主应力及最大位移
    Table  5.  The principal stress and maximum displacement of the loess tunnel lining
    w衬砌位置σmax/MPaσmin/MPas/cm
    7.42%左拱肩3.02-2.642.72
    右拱肩2.64-3.022.70
    左拱脚2.60-2.972.69
    右拱脚2.96-2.612.67
    9.36%左拱肩2.97-2.632.95
    右拱肩2.62-3.002.92
    左拱脚2.58-2.942.91
    右拱脚2.92-2.582.89
    11.95%左拱肩2.75-2.553.52
    右拱肩2.52-2.823.50
    左拱脚2.46-2.573.48
    右拱脚2.56-2.473.46
    13.25%左拱肩2.09-1.883.90
    右拱肩1.80-2.163.87
    左拱脚1.73-2.103.85
    右拱脚2.02-1.833.83
    15.10%左拱肩1.92-1.674.12
    右拱肩1.58-1.994.10
    左拱脚1.51-1.844.08
    右拱脚1.84-1.624.06
    注:σmax为最大主应力,σmin为最小主应力,s为最大位移.
    下载: 导出CSV 
    | 显示表格

    在水平地震加速度作用下,随着含水量的增加,左、右拱肩处主应力绝对值减小,但各位置处最大位移值增大.当黄土围岩含水量增湿到塑性范围时,黄土隧洞在地震作用下表现为塑性破坏而发生数值模拟不收敛.取左、右拱肩处的位置为黄土隧洞分析点,发现当左拱肩出现最大主应力时,右拱肩出现最小主应力,且其绝对值大小基本相等,发生时间也相同.

    图 4给出了含水量w=9.36%时隧洞左、右拱肩的主应力图. 图中纵轴上部为正向主应力,下部为负向主应力. 可以看出,在同一含水量情况下,随着时间的递增,左拱肩首先出现正向的主应力峰值,同时右拱肩出现负向的主应力峰值,左右拱肩主应力关于σ=0对称相等.由此可以得出,黄土隧洞衬砌材料的破坏原因是由于地震作用下,最大、最小主应力交替出现而使材料发生疲劳损伤所造成的.

    图  4  左拱肩(a)和右拱肩(b)主应力图
    Figure  4.  The principal stress of left- (a) and right-spandrel (b)

    本文采用地震动力时程分析法,基于初始弹性模量和抗剪强度指标的变化,模拟了黄土围岩增湿后隧洞一系列的动力特性.数值模拟结果显示:

    1)随着含水量的增加,同一阶的自振圆频率呈减小的趋势,且减小的趋势在含水量w≥11.95%时变慢,相应的自振周期随着含水量的增加而增大;

    2)随着含水量的增加,质量阻尼系数η减小,刚度阻尼系数λ变大;

    3)随着含水量的增加,相对于输入地震加速度,数值模拟分析得到的隧洞顶部加速度先大于输入地震加速度,随后又小于输入地震加速度;

    4)衬砌拱肩和拱脚位置为抗震薄弱位置,与汶川地震后公路隧道震害结果相类似,且随着含水量的增加,拱肩和拱脚处主应力绝对值呈减小而最大位移呈增大趋势.

    基于上述分析结果,本文提出了黄土围岩增湿后隧洞顶部加速度临界含水量概念,并求得引洮总干渠黄土隧洞顶部加速度临界含水量为10.12%;以含水量w=9.36%时左、右拱肩主应力图为例,分析了衬砌破坏的原因是由于水平地震作用下衬砌拱肩最大、最小主应力交替出现而使材料发生疲劳损伤造成的.因此,加强黄土隧洞拱肩和拱脚位置的支护对于其抗震性能的提升是必不可少的.此外,还应避免在隧洞施工及运营过程中,黄土围岩因外界因素使其含水量增大而导致其抗震性能降低.

计量
  • 文章访问数:  926
  • HTML全文浏览量:  9
  • PDF下载量:  128
  • 被引次数: 0
出版历程
  • 发布日期:  2011-09-01

目录

/

返回文章
返回