Investigation and Prevention Suggestion of Earthquake Disaster in the Extreme Earthquake Area of MS6.9 Earthquake in Menyuan County, Qinghai Province
-
摘要: 北京时间2022年1月8日1时45分青海省海北藏族自治州门源县(37.77°N,101.26°E)发生6.9级地震,此次地震极震区为距震中4 km的兰新高铁硫磺沟大桥附近(Ⅸ度区,面积约157 km2)。调查结果显示,在极震区大西沟至硫磺沟区域,道路、桥梁及房屋破坏较严重。沿区域内长约22 km的托来山-冷龙岭断裂带附近,道路形成多处裂缝与挤压鼓包,桥梁整体倾斜移位,房屋不同程度破坏。对极震区兰新高铁硫磺沟大桥及周边道路、房屋进行实地调查与震害分析,提出灾后重建及震害防御意见:对于灾区房屋建筑,组织专业技术人员进行详细调查及安全鉴定;建议定期对房屋建筑进行隐患排查及加固;对于兰新高铁硫磺沟大桥,建议原地修建,加设减隔震装置、连梁及柔性限位装置。Abstract: At 1:45 a.m. on Jan 8th, 2022, a magnitude 6.9 earthquake occurred in Menyuan County (37.77°N, 101.26°E), Haibei Tibetan Autonomous Prefecture, Qinghai Province. The epicenter of the earthquake is near the Liuhuanggou Bridge on the Lanxin high-speed railway, which is 4 km away from the epicenter (the area of earthquake intensity 9 is about 157 km2). The survey results show that roads, bridges and houses were seriously damaged from Daxigou to Liuhuanggou in the extreme disaster area. In the vicinity of the 22 km Long Tuolaishan-Lenglongling fault zone in the region, it can be seen that the damage type of several cracks and extrusion bulges were formed on the ground, the bridge tilted and shifted as a whole, and the house was damaged to varying degrees. The field investigations and earthquake disaster analysis were carried out on the Lanxin High-speed Railway Liuhuanggou Bridge and surrounding road and houses in the extreme disaster zone, and relevant suggestion on post-earthquake reconstruction and earthquake disaster prevention were put forward: Conducting detailed house safety appraisal by the professional personnel of housing construction organization in disaster area; Suggesting to conduct hidden danger investigation and reinforcement work on housing construction regularly; For the Liuhuanggou bridge, it is recommended to rebuild in situ with shock-absorbing and isolating devices, limiting and connecting beam device.
-
引言
我国渤海海域的油气资源较为丰富,大部分油气资源集中在浅海大陆架上,占我国石油年产量的1/3(薛永安等,2020),渤海海域布设了大量导管架海洋平台用于开采海上油气资源。渤海浅海海域丰富的油气资源与该地区频繁的地质构造运动密切相关,如郯庐断裂带晚期强烈地质活动形成了大量构造圈闭储藏油气(万桂梅等,2009),但这也导致地震等灾害的发生。历史上,渤海海域共发生有记录的7级以上地震4次,6级以上地震2次,地震海啸1次(王健,2007)。虽然渤海海域纬度较高,但近年来时有台风发生,如2014年的“麦德姆”台风,2019年的“利奇马”台风,2020年的“巴威”台风。导管架平台的设计使用年限一般都会达到15年以上,存在部分超期服役的导管架平台,如果地震和台风同时或先后发生,必然会对平台造成损坏,甚至会造成严重的次生灾害。作为渤海浅海海域海洋平台的主要结构形式,生产中的导管架平台如果发生破坏,当地自然环境必然遭受严重破坏,产生重大安全事故。因此,对地震-台风耦合作用下的导管架平台进行相应的动力响应分析具有重要意义。
国内外学者针对海洋平台的动力响应进行了一系列研究。张光发等(2011)用解析方法分析导管架下水系统三维运动得到的数学模型优化了导管架下水参数的计算;吴家鸣等(2009)对桁架式近海构筑物进行数值分析,结果表明结构物对波浪的响应主要取决于其自身尺寸与波浪环境参数;刘育丰等(2012)对渤海区域海洋平台抗震设防经验进行总结,对不同地震危险性分区给出了设防参数;孙久洋等(2020)对比增量动力法和耐震时程法对强震作用下的导管架进行分析后,认为耐震时程分析法的计算效率更高;董汝博等(2010)基于流固耦合理论对海洋储油罐平台地震反应进行数值模拟,并与简化模型进行对比,结果表明采用流固耦合的方法更加精确;Wu等(2018)通过在海洋平台上加装非线性磁力调制阻尼器进行振动台实验,发现调制阻尼器在宽频段内的抑制作用十分显著且能解决平台在大震下延迟响应等问题;张大勇等(2015)对钻井平台进行冰激动力响应分析,并对冰振安全评估提出了建议;Sharma等(2019)通过在导管架平台上设置调质阻尼器,分析其在地震和冰荷载作用下的动力响应,结果表明设置调质阻尼器能够有效减小平台层的位移;朱本瑞(2014)进行大型导管架平台风洞实验,得到了平台周围风场变化规律;左华楠(2017)分析了导管架平台在恶劣风浪条件下的动力响应,给出了其在不同水深和上浪高度下的响应规律及甲板的应力变化情况。目前,针对地震-台风耦合作用下海洋平台动力响应分析的研究仍有待丰富。
针对导管架平台这类重要工程,在地震、台风等自然灾害作用下的灾害风险评估、隐患排查和安全评估具有一定实际意义。本文针对地震-台风耦合作用下导管架平台的安全评估,开展如下研究:
(1)通过Morison方程及Stokes五阶波理论将波浪力、台风拖曳力施加至结构上,建立运动方程;
(2)基于Abaqus建立地震-台风耦合作用下的导管架平台数值分析模型;
(3)对模型进行模态分析和非线性时程响应计算,对数值模拟后得到的应力、顶层位移角等结果进行分析,并与美国石油学会(American petroleum institute, API)的相关规范和《海上平台状态评定指南》(GD 04—2005)(中国船级社,2005)中所要求的响应参数进行对比,评估导管架平台在地震-台风耦合作用下的安全性能。
1. 地震-台风耦合环境下荷载施加基本理论
1.1 波浪力施加
对于流体中的刚架或桁架等构件直径较波浪波长小的结构,可通过Morison方程求得结构所受的波浪力,波浪力主要包括拖曳力和惯性力。流体的拖曳力与稳态流的速度和特定波浪形式有关,惯性力与波浪的加速度有关。
流体在杆件轴向正交方向的拖曳力为:
$$ {{\boldsymbol{F}}_{\text{D}}}{\text{ = }}\frac{{\text{1}}}{{\text{2}}}\rho {C_{\text{D}}}D\Delta {{\boldsymbol{v}}_{\text{n}}}{\left( {\Delta {{\boldsymbol{v}}_{\text{n}}} \cdot \Delta {{\boldsymbol{v}}_{\text{n}}}} \right)^{\frac{1}{2}}} $$ (1) 流体质点的相对速度沿杆件轴向正交方向分量为:
$$ \Delta {{\boldsymbol{v}}_{\text{n}}} = \Delta {\boldsymbol{v}}{{ - }}\Delta {{\boldsymbol{v}}_{\text{t}}} $$ (2) 流体质点相对速度沿杆件轴向分量为:
$$ \Delta {{\boldsymbol{v}}_{\text{t}}} = \left( {\Delta {\boldsymbol{v}} \cdot {\boldsymbol{t}}} \right){\boldsymbol{t}} $$ (3) 流体质点相对速度为:
$$ \Delta {\boldsymbol{v}} = {{\boldsymbol{v}}_{\text{f}}} - {\alpha _{\text{R}}}{{\boldsymbol{v}}_{\text{p}}} $$ (4) 式中,
$ {{\boldsymbol{v}}_{\text{f}}} $ 为流体质点速度;$ {{\boldsymbol{v}}_{\text{p}}} $ 为导管架结构质点速度;$ {\alpha _{\text{R}}} $ 为结构速度常数;$ {\boldsymbol{t}} $ 为沿着单元轴向的单位向量;$ \rho $ 为流体密度;CD为拖曳系数;D为杆件外径。由于流体惯性引起的单位长度上的惯性力为:
$$ {{\boldsymbol{F}}_{\text{I}}}{\text{ = }}\rho \frac{{{\text{π}} {D^2}}}{4}\left[ {{C_{\text{M}}}{{\boldsymbol{a}}_{{\text{fn}}}} - {C_{\text{A}}}{{\boldsymbol{a}}_{{\text{pn}}}}} \right] $$ (5) 式中,CM为流体惯性系数;CA为附加质量系数;
$ {{\boldsymbol{a}}_{{\text{fn}}}} $ 为流体在杆件轴向正交方向上的加速度;$ {{\boldsymbol{a}}_{{\text{pn}}}} $ 为流体在杆件轴向正交方向上的加速度。1.2 台风荷载的施加
台风与波浪类似,也会对结构产生拖曳力,但流体对结构的拖曳力仅在液面以下产生,而台风对结构的拖曳力在液面以上结构部分产生,且这部分结构上的风速要符合假设的风速廓线,本文采用指数律风速廓线模式。该部分拖曳力以静力的形式施加在结构上。
指数律风速廓线为:
$$ u = u_{\left( t \right)}^0{\left( {\frac{z}{{{z_0}}}} \right)^\alpha } $$ (6) 风作用在结构上的拖曳力为:
$$ F = \left( {\rho /2} \right){u^2}{C_{\text{S}}}A $$ (7) 式中,u为风速;
$ u_{\left( t \right)}^0 $ 为$ {z_0} $ 处根据时间变化的风速;$ z $ 为海面上任一位置距海面的距离,$ z = 0 $ 时为海面;$ {z_0} $ 为某一已知风速位置距海面的距离;α为地面粗糙度指数;ρ为空气密度;CS为形状系数;A为结构的受力面积。1.3 运动方程建立
考虑地震-台风耦合作用下的结构运动方程(何晓宇等,2007):
$$ \left[\boldsymbol M \right]\left\{ {\ddot {\boldsymbol{x}}} \right\} + \left[ \boldsymbol C \right]\left\{ {\dot {\boldsymbol{x}}} \right\} + \left[ \boldsymbol K \right]\left\{ x \right\} = - \left[ \boldsymbol M \right]\left\{ {{{\ddot {\boldsymbol{x}}}_{\text{g}}}} \right\}{\text{ + }}\left\{ {{{\boldsymbol{f}}_{\text{H}}}} \right\} $$ (8) 式中,
$ \left[ \boldsymbol M \right] $ 为结构质量矩阵;$ \left[ \boldsymbol C \right] $ 为结构阻尼矩阵;$ \left[ \boldsymbol K \right] $ 为结构刚度矩阵;$ \left\{ {\ddot {\boldsymbol{x}}} \right\} $ 为结构加速度向量;$ \left\{ {\dot {\boldsymbol{x}}} \right\} $ 为结构速度向量;$ \left\{ \boldsymbol x \right\} $ 为结构位移向量;$ \left\{ {{{\ddot {\boldsymbol{x}}}_{\text{g}}}} \right\} $ 为地震动加速度向量;$ \left\{ {{{\boldsymbol{f}}_{\text{H}}}} \right\} $ 为考虑了台风、波浪的环境荷载向量。环境荷载向量为:
$$ \left\{ {{{\boldsymbol{f}}_{\text{H}}}} \right\} = \left\{ {{{\boldsymbol{F}}_{\text{D}}}} \right\} + \left\{ {{{\boldsymbol{W}}_{\text{D}}}} \right\} + \left\{ {{{\boldsymbol{F}}_{\text{I}}}} \right\} $$ (9) 式中,
$ \left\{ {{{\boldsymbol{F}}_{\text{D}}}} \right\} $ 为流体引起的拖曳力向量;$ \left\{ {{{\boldsymbol{W}}_{\text{D}}}} \right\} $ 为台风引起的拖曳力向量;$ \left\{ {{{\boldsymbol{F}}_{\text{I}}}} \right\} $ 为流体引起的惯性力向量。运动方程建立后,采用ABAQUS中的动力显式分析方法,求解运动方程。
2. 有限元模型建立
2.1 环境及导管架参数
本文选取渤海海域某导管架海洋平台进行分析(图1)。平台采用四腿导管架型式,导管架顶标高5.0 m,底标高−13.4 m,工作点标高6.6 m。主导管采用Φ1350×24钢管,成矩形布置。在标高−4.0 m、−11.9 m之间设竖向斜拉筋。平台甲板共2层,上层平台梁顶标高12.0 m,平台主尺寸为18.5 m×19.0 m;下层平台梁顶标高10.0 m,平台尺寸为7.6 m×7.3 m。
潮位选择该海区不规则半日潮校核高水位3.08 m。对于地震-台风耦合作用的环境,由于波高较大,波面高度和水质点速度采用Stokes五阶波理论,波浪波高选择校核高水位最大可能波高8.0 m,波浪周期为8.6 s,浪向NE,台风风向为强风向NW。根据《海港水文规范》(JTS 145-2-2013)(中华人民共和国交通运输部,2013)的要求,本文流体拖曳力系数取1.2,惯性系数取2.0,地面粗糙度指数根据《建筑结构荷载规范》(GB 50009—2012)(中华人民共和国住房和城乡建设部等,2012)的要求,近海海面取0.12。风荷荷载的形状系数根据American Petroleum Institute(2014)取0.5。
2.2 导管架平台有限元模型
本文采用ABAQUS/AQUA模块进行动力响应分析。AQUA是ABAQUS中专门应用于海洋工程的分析模块,可用于水下结构或部分水下结构上的稳态流、波浪和风荷载计算。波浪和风的参数通过AQUA模块进行录入,拖曳力和惯性力通过子程序进行施加。下部导管架用B31单元进行模拟,上部甲板用S4R单元进行模拟。考虑到甲板上布设有设施和设备,通过质量点来进行模拟。在海床面以下的桩,采用固定约束进行模拟。本模型共包含5314个梁单元,1576个壳单元。导管架平台模型如图2所示。
导管架平台材料选择Q345钢,密度为7 850 kg/m3,弹性模量为206 GPa,泊松比为0.28,屈服应力为345 MPa,考虑材料的塑性变形,塑性阶段的应力-应变关系如图3所示,具体实验数据如表1所示(谢海波等, 2004)。
表 1 Q345钢应力-塑性应变参数Table 1. Stress-plastic strain parameters of Q345应力/MPa 塑性应变 276.00 0 300.48 1.4e-5 320.16 5.5e-5 333.96 1.24e-4 342.24 2.21e-4 345.00 3.45e-4 345.00 0.01338 2.3 模态分析
对导管架平台进行模态分析,前8阶振型的自振频率如表2所示,由表可知前3阶的自振频率较后面振型自振频率小的多。前3阶振型如图4所示,第1振型为y方向弯曲振型,第2振型为x方向弯曲振型,第3振型为扭转振型。
表 2 前8阶自振频率Table 2. The first 8 order natural frequencies振型编号 自振频率/Hz 1 2.8018 2 2.8753 3 2.9002 4 3.3711 5 3.4598 6 6.1433 7 6.1440 8 7.0278 根据《建筑结构荷载规范》(GB 50009—2012),参考高耸结构的基本自振周期经验公式,钢结构取下式经验系数中的较大值:
$$ T{\text{ = }}\left( {{\text{0}}{{.007\sim0}}{\text{.013}}} \right)H $$ (10) 式中,H为结构高度。
本文海洋平台属于钢结构,故取经验系数为0.013,结构高度为25.4 m,按式(10)计算自振周期为0.33 s,即自振频率为3.03 Hz,与模态分析得到的结构自振频率较为接近,说明本文采用的导管架平台建模方法合理可靠。
2.4 动力响应分析
2.4.1 地震动选取
本文选取2008年美国Chino Hills海底地震动CH08,震级Mw5.4,截取包含最大响应的20s地震动时程,并对地震动时程进行滤波、调幅,得到0.10 g、0.15 g、0.20 g和0.40 g共4种地震动强度,其中0.20 g和0.40 g分别为设防地震动强度和罕遇地震动强度,通过基底进行输入。该地震动x、y方向加速度时程如图5所示。对地震动进行频谱分析,得到其加速度傅里叶谱(图6),可见所选地震动的卓越频率与导管架平台的基本频率较为接近,符合最不利设计地震动要求(翟长海等,2005)。
2.4.2 波浪可视化
为使台风作用下的大浪在数值模型中可视化,进而更直观的观察地震-台风耦合作用下导管架平台的响应,本文在海面上建立了一个20 m×50 m的平面,用表面单元进行划分,来模拟波浪的变化。不同时间下波浪的形态如图7所示。
2.4.3 地震-台风耦合工况
近年来,受台风影响的渤海海域风速如表3所示。本文选取2种风速,第1种为近年来渤海海域出现的台风中风速最大的“利奇马”台风,其在山东登陆时风速为23 m/s;第2种为该导管架设计说明书中重现期为50年的3秒钟阵风风速,风速为33.6 m/s。地震动强度选择0.1 g、0.15 g、0.2 g和0.4 g。本文共考虑10种组合工况用于地震-台风耦合作用下的动力响应计算,如表4所示。
表 3 渤海海域风速Table 3. Wind speed of Bohai sea台风名称 风速/m·s−1 布拉万(2012) 20.8 利奇马(2019) 23 巴威(2020) 17.2 表 4 工况表Table 4. Table of working conditions地震动强度/g 风速/m·s−1 23 33.6 0 0.1 工况1 工况2 — 0.15 工况3 工况4 — 0.2 工况5 工况6 — 0.4 工况7 工况8 工况9 0 — 工况10 — 2.4.4 耦合作用下的动力响应
根据《海上平台状态评定指南》(GD 04—2005)结构动力分析评价中的强度要求,参考American Petroleum Institute(2014)相关规定,构件许用应力按下式进行计算:
$$ {F_{\text{t}}} = 0.6{F_{\text{y}}} $$ (11) 式中,Fy为构件屈服应力。
对于Q345钢,按上式计算得到其一般条件下构件许用应力为207 MPa。同时,API规范规定,对于极端环境条件下的构件应力校核,许用应力可相较于一般条件下构件应力校核提升70%,即允许少数构件出现屈服,而没有明显的破坏发生。
10种工况计算得到的最大应力如图8所示,除工况8外,其余9种工况的最大应力位置相同(图10),而工况8最大应力位置在C、D间的桩腿上(图11)。工况8~10的应力云图如图12所示,可以发现在地震-台风耦合作用下的导管架平台最大应力明显高于地震或台风单独作用下的最大应力,此时导管架帽下的立柱应力水平过高,接近屈服。
由于钢管构件上的应力相对较大,本文对80个钢管构件进行统计,以207 MPa为评估指标,得到在20 s的作用时间内不同工况下超过该指标的构件个数,如图9所示。其中,工况2、4、6、7超指标构件集中在导管架帽的4根立柱处,即图11所示E高度处的4根立柱;工况8超指标构件除导管架帽的4根立柱处外,还集中在图11所示C、D高度之间钢管构件处;工况9超指标构件集中在图11所示C、D高度之间;工况10超指标构件集中在导管架帽的2根立柱处。
对比工况1~6 ,当风速为23 m/s时,导管架平台的应力水平较低,不存在超指标构件;当风速提高至33.6 m/s时,导管架平台的应力水平显著提高,超指标构件占比也随着地震动强度增加而提升,0.10 g时占总构件数的6.25%,0.15 g时占总构件数的10%,0.20 g时占总构件数的12.5%。
工况7、8为罕遇地震-台风耦合作用下的极端环境条件。工况7没有构件出现屈服,有5个构件超过评估指标,占总构件数的6.25%;工况8有16个构件应力超过评估指标,占总构件数的20%,其中有1个构件接近材料的屈服强度。相较于地震动强度为0.10 g、0.15 g和0.20 g的情况,在罕遇地震动强度下,风速的增大导致最大应力的增幅相对减弱。
对比工况2、4、6、8发现,在相同的风速条件下,导管架平台应力水平有所上升,但相较同一地震动强度下风速增大,应力水平增幅相对较小。
对比工况8、9、10发现,地震-台风耦合作用下超指标构件个数明显高于地震或台风单独作用下的超指标构件个数。
综上所述,地震-台风耦合作用对导管架平台的安全性能有较大影响,尤其在台风风速较快的情况下,必须考虑耦合作用的影响。
本文参考毛晨曦等(2018)对于通信铁塔的损伤标准,对导管架平台损伤状态进行评估。选取导管架平台顶点相对于基底的位移角(Relatively Displacement Angle, RDA)峰值作为评定导管架平台损伤状态的指标,即RDA为导管架平台顶相对于基底的最大位移与平台总高的比值,按下式进行计算:
$$ {\text{RDA = }}{\Delta _{\max }}/H $$ (12) 将导管架平台顶所能达到的最大顶层位移角,即极限位移角,用
$ {\text{RD}}{{\text{A}}_{{\text{ult}}}} $ 表示;屈服位移角用$ {\text{RD}}{{\text{A}}_{{\text{ye}}}} $ 表示;0.5倍的$ {\text{RD}}{{\text{A}}_{{\text{ye}}}} $ 用$ {\text{RD}}{{\text{A}}_{{\text{ud}}}} $ 表示,采用表5对导管架平台损伤状态进行定义。表 5 导管架平台损伤状态Table 5. Definition of damage state of jacket platform$ {\rm{RDA}} \leqslant {\rm{RD}}{{{A}}_{{\rm{ud}}}} $ $ {\rm{RD}}{{{A}}_{{\rm{ud}}}}{{ < \rm RDA}} \leqslant {\rm{RD}}{\rm{{A}}_{{\rm{ye}}}} $ $ {\rm{RD}}{{\rm{A}}_{{\rm{ye}}}}{{ < \rm RDA}} \leqslant {\rm{RD}}{{\rm{A}}_{{\rm{ult}}}} $ $ {\rm{RDA > RD}}{{\rm{A}}_{{\rm{ult}}}} $ 基本完好 轻微破坏 严重破坏 毁坏 为确定导管架平台的屈服位移角和极限位移角,对平台施加侧向倒三角形荷载,进行静力推覆分析。导管架平台的Pushover曲线如图13所示。
依据导管架平台损伤状态定义,由静力推覆分析得到各损伤状态限值,判断5种工况作用下损伤状态如表6所示。
表 6 各工况下导管架平台损伤状态Table 6. Damage state of jacket platform under each working condition工况 1 2 4 6 8 损伤状态 基本完好 轻微破坏 轻微破坏 轻微破坏 严重破坏 考虑到超指标构件主要集中在高度D、F之间和高度B、D之间,本文将选取的80个钢管构件进一步细分,分为28个上部导管架帽构件和52个下部导管架构件,对工况1、2、4、6和8进行统计,结果如表7所示。结果表明,当导管架平台处于基本完好时,没有构件超过指标限值;当处于轻微破坏时,导管架帽构件超过指标限值的比例小于1/4,导管架构件超过指标限值的比例小于1/20;当处于严重破坏时,导管架帽超过指标限值的比例小于1/4,导管架构件超过指标限值的比例小于1/5(表8)。
表 7 超指标构件数量统计Table 7. Statistics on the number of super-index components工况 导管架帽构件 导管架构件 超指标构件 1 0 0 0 2 5 0 5 4 8 0 8 6 8 2 10 8 8 8 16 表 8 不同损伤状态下超指标构件占比Table 8. Proportion of super-index components under different damage states损伤状态 超指标导管架帽构件比例 超指标导管架构件比例 基本完好 ≤0 ≤0 轻微破坏 ≤1/4 ≤1/20 严重破坏 ≤1/4 ≤1/5 3. 结论
针对近海导管架海洋平台在地震-台风耦合作用下的响应分析结果,可得到以下结论:
(1)地震-台风耦合作用下的导管架平台安全性能明显小于地震或台风单独作用下的安全性能,尤其是在台风风速较快的情况下,平台损伤较为严重。
(2)绝大多数工况下的最大应力均出现在同一构件处,在导管架平台设计时应对结构薄弱环节做适当加强,防止因局部构件的破坏导致整体结构倒塌。
(3)通过研究给出了导管架平台损伤状态的定量描述,为灾害风险评估和现场抢修工作提供了参考依据。
-
表 1 2022年门源MS6.9地震震源参数
Table 1. Source parameters for the 2022 Menyuan MS6.9 earthquakes
来源 纬度 经度 震级/MS 走向/° 倾角/° 滑动角/° 深度/km 中国地震台网中心 37.770°N 101.260°E 6.9 — — — 10.0 USGS 37.815°N 101.278°E 6.6 104 88 15 11.5 GCMT 37.800°N 101.310°E 6.7 104 82 1 14.8 GFZ 37.810°N 101.340°E 6.6 285 82 16 15.0 表 2 2022年门源MS6.9地震强震动记录分析
Table 2. Analysis of strong vibration records for the 2022 Menyuan MS 6.9 earthquakes
台站
名称台站经度 台站纬度 震中距/km 震中
方位角/°PGA/cm·s−2 PGV/cm·s−1 仪器
烈度/°EW NS UD EW NS UD C0028 37.7°N 101.3°E 7.8 324.6 −456.9 445.0 355.3 27.6 23.4 12.0 8.2 C0029 37.6°N 101.2°E 13.7 4.4 210.8 −144.1 132.7 −32.6 16.6 7.6 8.3 C0036 37.8°N 101.1°E 15.1 120.0 −136.0 −134.6 −122.4 −18.7 −30.4 18.7 8.2 C0027 37.6°N 101.4°E 18.0 332.4 −242.6 198.7 −114.3 16.1 9.6 −7.6 7.6 表 3 硫磺沟大桥场地等效剪切波速
Table 3. Equivalent shear wave velocity of Liuhuanggou bridge site
钻孔编号 等效剪切波速/m·s−1 覆盖层厚度/m LHG-3 365 24.0 LHG-7 398 18.0 表 4 极震区不同结构类型房屋震害调查表
Table 4. Investigation on earthquake disaster of buildings with different structures in the extreme areas
调查点 结构类型 砖混结构 轻钢结构 土石结构 砖木结构 数量/间 面积/m2 破坏情况 数量/间 面积/m2 破坏情况 数量/间 面积/m2 破坏情况 数量/间 面积/m2 破坏情况 1 1 20 轻微破坏 — — — — — — — — — 2 — — — 26 1000 16间基本完好
8间轻微破坏
2间中等破坏— — — — — — 3 — — — 46 2 000 42间基本完好
4间轻微破坏— — — — — — 4 — — — 2 30 基本完好 2 200 毁坏 — — — 5 — — — — — — — — — 1 15 毁坏 6 — — — — — — — — — 1 20 毁坏 总计 1 20 — 74 3030 — 2 200 — 2 35 — -
[1] 何文贵, 刘百篪, 袁道阳等, 2000. 冷龙岭活动断裂的滑动速率研究. 西北地震学报, 22(1): 90—97He W. G. , Liu B. C. , Yuan D. Y. , et al. , 2000. Research on slip rates of the Lenglongling active fault zone. Northwestern Seismological Journal, 22(1): 90—97. (in Chinese) [2] 胡朝忠, 杨攀新, 李智敏等, 2016.2016年1月21日青海门源6.4级地震的发震机制探讨. 地球物理学报, 59(5): 1637—1646 doi: 10.6038/cjg20160509Hu C. Z. , Yang P. X. , Li Z. M. , et al. , 2016. Seismogenic mechanism of the 21 January 2016 Menyuan, Qinghai MS6.4 earthquake. Chinese Journal of Geophysics, 59(5): 1637—1646. (in Chinese) doi: 10.6038/cjg20160509 [3] 黄勇, 王君杰, 韩鹏等, 2010. 考虑支座破坏的连续梁桥地震反应分析. 土木工程学报, 43(S2): 217—223Huang Y. , Wang J. J. , Han P. , et al. , 2010. Seismic response analysis of continuous bridges taking account of bearing failure. China Civil Engineering Journal, 43(S2): 217—223. (in Chinese) [4] 蒋一然, 宁杰远, 温景充等, 2022. 高铁地震波场中的多普勒效应及应用. 中国科学: 地球科学, 52(3): 438—449.Jiang Y. R. , Ning J. Y. , Wen J. C. , et al. , 2022. Doppler effect in high-speed rail seismic wavefield and its application. Science China Earth Sciences, 65(3): 414—425. (in Chinese) [5] 李红, 2011. 玉树地震房屋震害的几点启示. 青海师范大学学报(自然科学版), 27(3): 87—89 doi: 10.3969/j.issn.1001-7542.2011.03.023Li H. , 2011. Several views from Yushu earthquake damages on the houses. Journal of Qinghai Normal University (Natural Science), 27(3): 87—89. (in Chinese) doi: 10.3969/j.issn.1001-7542.2011.03.023 [6] 李鑫, 姚生海, 殷翔等, 2021a. 青海玛多7.4级地震极震区震灾调查及分析. 震灾防御技术, 16(3): 429—436Li X. , Yao S. H. , Yin X. , et al. , 2021a. Investigation and analysis of earthquake disaster in the extreme earthquake area of MS7.4 earthquake in Qinghai Province. Technology for Earthquake Disaster Prevention, 16(3): 429—436. (in Chinese) [7] 李鑫, 殷翔, 姚生海等, 2021b. 青海玛多7.4级地震重灾区房屋震灾调查及分析. 地震工程学报, 43(4): 896—902Li X. , Yin X. , Yao S. H. , et al. , 2021b. Investigation and analysis of earthquake disasters of houses in the stricken areas of Maduo M7.4 earthquake in Qinghai Province. China Earthquake Engineering Journal, 43(4): 896—902. (in Chinese) [8] 李振洪, 韩炳权, 刘振江等, 2022. InSAR数据约束下的2016年和2022年青海门源地震震源参数及其滑动分布. 武汉大学学报(信息科学版), 1—15. (2022-01-13). https: //kns. cnki. net/kcms/detail/detail. aspx?dbcode=CAPJ&dbname=CAPJLAST&filename=WHCH20220113000&uniplatform=NZKPT&v=7ezb_b4zijR-Rvx9NvQtho8_rGD3PRIqoo9xZV3eZHzPyNvlFQbRwuh-ybDuW5RV.Li Z. H., Han B. Q., Liu Z. J., et al., 2022. Source parameters and slip distributions of the 2016 and 2022 Menyuan, Qinghai earthquakes constrained by InSAR observations. Geomatics and Information Science of Wuhan University, 1—15. (2022-01-13). https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CAPJ&dbname=CAPJLAST&filename=WHCH20220113000&uniplatform=NZKPT&v=7ezb_b4zijR-Rvx9NvQtho8_rGD3PRIqoo9xZV3eZHzPyNvlFQbRwuh-ybDuW5RV. (in Chinese) [9] 刘佳, 2013. 桥梁球形支座用新型复合滑板材料的结构、性能与机理研究. 北京: 北京化工大学.Liu J., 2013. Strcture, performance and mechanism research of new composite board for bridge spherical bearing. Beijing: Beijing University of Chemical Technology. (in Chinese) [10] 罗春燕, 2011. 玉树地震中房屋结构的震害分析. 山西建筑, 37(26): 53—54 doi: 10.3969/j.issn.1009-6825.2011.26.031Luo C. Y. , 2011. On analysis of earthquake diseases of house structures in Yushu earthquake. Shanxi Architecture, 37(26): 53—54. (in Chinese) doi: 10.3969/j.issn.1009-6825.2011.26.031 [11] 罗全波, 陈学良, 高孟潭等, 2018. 近断层速度脉冲与震源机制的关系浅析. 震灾防御技术, 13(3): 646—661Luo Q. B. , Chen X. L. , Gao M. T. , et al. , 2018. Relationship between Near-fault velocity pulse and focal mechanism. Technology for Earthquake Disaster Prevention, 13(3): 646—661. (in Chinese) [12] 王青桥, 韦晓, 王君杰, 2009. 桥梁桩基震害特点及其破坏机理. 震灾防御技术, 4(2): 167—173 doi: 10.3969/j.issn.1673-5722.2009.02.005Wang Q. Q. , Wei X. , Wang J. J. , 2009. Characteristics and mechanisms of earthquake damage of bridge pile foundation. Technology for Earthquake Disaster Prevention, 4(2): 167—173. (in Chinese) doi: 10.3969/j.issn.1673-5722.2009.02.005 -