Validity Verification and Correction Method of Horizontal and Vertical Spectrum Ratio Method of Microtremor in Site Effect Study
-
摘要: 近年来，地脉动水平/竖向谱比法广泛应用于获取场地卓越频率工作中，但在指示地震动放大系数方面是偏小的。利用KiK-net台网403个台站共21万余组数据，通过对比分析地脉动水平/竖向谱比MHVR和地震动水平/竖向谱比EHVR，发现二者的谱形和峰值频率基本一致，但幅值存在差异。根据MHVR谱形将场地分为6类，建立将MHVR修正为EHVR的经验方法。当MHVR有可解释的峰值时，实测EHVR和修正后MHVR预测值一致性提高。地震动峰值加速度在40 cm/s2以下时，EHVR受场地非线性的影响较小，利用MHVR及其修正方法可有效估计EHVR。Abstract: In recent years, the horizontal/vertical spectrum ratio method of microtremor has been widely used to obtain the site dominant frequency, but it is usually too small to indicate the amplification coefficient of ground motion. Based on more than 210,000 sets of data from 403 stations of KiK-net, the horizontal/vertical spectrum ratio of microtremor (MHVR) and horizontal/vertical spectrum ratio of earthquake ground motion (EHVR) are compared and analyzed. It is found that the spectrum shape and peak frequency of the two are consistent, but the amplitude is different. According to the MHVR spectrum, the site was divided into six categories, and an empirical method was established to modify MHVR into EHVR. When MHVR has an interpretable peak value, the consistency between the measured EHVR and the revised MHVR predicted value is significantly improved. When the peak acceleration of ground motion is less than 40 cm/s2, EHVR is less affected by site nonlinearity. Therefore, MHVR and its correction method could be used to estimate EHVR effectively.
-
Key words:
- Strong ground motion /
- Site effect /
- Microtremor /
- Horizontal/Vertical spectrum ratio
-
引言
2017年5月11日，新疆维吾尔自治区喀什地区塔什库尔干县发生MS 5.5地震，地震造成8人死亡、31人受伤，造成房屋及设施破坏，直接经济损失共20.05亿元（侯建盛等，2017）。
此次地震灾区主要涉及喀什地区塔什库尔干县科克亚尔柯尔克孜民族乡、塔合曼乡、提孜那普乡、塔什库尔干乡、班迪尔乡、巴扎达什牧林场（行政隶属班迪尔乡）、瓦恰乡、达布达尔乡等9个乡镇。灾区面积3288km2，受灾人口26486人，9285户，由于房屋毁坏和较大程度破坏造成失去住所人数共计16194人，4753户。
1. 地震构造背景及基本情况
1.1 地震烈度
此次地震震中位于新疆喀什地区塔什库尔干县塔什库尔干乡，宏观震中位于塔什库尔干乡库孜滚村，为Ⅷ度异常点。通过对灾区9个乡（镇、场）的69个调查点展开实地调查，得到的烈度图等震线长轴呈北北西走向分布（图 1）。Ⅶ度区面积227km2，长轴28km，短轴8km，涉及塔什库尔干镇（含县城）和塔什库尔干乡；Ⅵ度区面积3061km2，长轴100km，短轴43km，涉及科克亚尔柯尔克孜民族乡、塔合曼乡、提孜那普乡、塔什库尔干乡、班迪尔乡、巴扎达什牧林场（行政隶属班迪尔乡）、瓦恰乡、达布达尔乡等8个乡（场）；Ⅵ度区及以上总面积为3288km2。
1 新疆地震局，2017.新疆塔什库尔干5.5级地震灾害损失评估报告.
1.2 发震构造
震区位于帕米尔高原塔什库尔干断陷谷地，该谷地是由青藏高原西北帕米尔构造结内部塔什库尔干拉张系晚新生代以来的拉张作用形成的盆地，其南北狭长，东西分布海拔为4000—5000m的高山。
震区内塔什库尔干断裂成型于华力西时期，有长期的演化发育史。大部分在喜马拉雅期重新复活，该断裂控制着塔什库尔干盆地的形成与演化，此次塔什库尔干MS 5.5地震就发生在塔什库尔干断裂带上（图 2）。
1 新疆地震局，2017.新疆塔什库尔干5.5级地震灾害损失评估报告.
从本次地震的新构造运动分区背景来看，新近纪以来震区所处的西昆仑隆起区隆起幅度大致在2—7km，该隆起区第四纪以来的隆起幅度和速率分别为1200—1700m和10—13mm/a。
1.3 震区场地条件
震区位于帕米尔东北—西昆仑区段，西昆仑山体呈北西—南东走向，平均海拔5000—6000m，主要山峰偏于西部。主峰公格尔山海拔7649m，慕士塔格山为7509m，山体宽厚高大，南北不对称，北坡长而陡峭，与海拔1000m多的塔里木盆地相邻，高差4000m。帕米尔高原实际上并非平坦的高原面，由几组山脉和山脉之间宽阔的谷地和盆地构成。
此次震中位于塔什库尔干谷地内，发育有塔什库尔干河，震区附近谷地与两侧高差700—1200m，谷地总体走向近南北，宽数千米，呈狭长状，谷地两岸冰碛物堆积及冲洪积堆积发育，村庄沿河流阶地及山前冲洪积扇分布（刘军等，2014），该地貌单元内场地类别为Ⅰ类，如图 3所示。
2. 震害特征分析
2.1 地震灾害特征分析
本次地震涉及影响范围内的房屋结构类型主要包括简易房（土石木结构）、砖混结构、砖木结构及少量的框架结构。简易房按照承重墙体可分成两类，其中一类主要分布在山前洪积扇倾斜平原地带，多为当地居民就地取材而建，屋顶结构为先搭建房梁后在其上搁置短木条作为椽子，在椽子上铺设草席后覆盖房泥，部分老旧房屋屋顶房泥较厚，墙体由卵石、粉土砌筑而成，粘结强度极差，加之施工质量和场地条件的影响，造成一定数量的毁坏和大面积破坏，是导致本次地震造成人员伤亡的主要原因；另一类主要分布在塔什库尔干河两岸阶地上，此类房屋多为土坯砌筑而成。由于该地区经济落后，交通极为不便，建造成本高，在县城存在大量2000年左右建设的土石木房屋。在地震中，第一类房屋大面积倒塌，房屋倒塌基本为整体性倾覆，第二类土坯房倒塌相对略少，多为局部倒塌，2种不具抗震性能的房屋破坏面积较严重，计算时均列入土石木结构房屋，该类房屋在乡镇Ⅶ度区毁坏达39.1%，在县城的Ⅶ度区毁坏达24.3%。震区各类结构房屋面积如表 1所示。
表 1 震区各类结构房屋面积（单位：m2）Table 1. Total areas of various kinds of structures in the earthquake area (unit: m2)行政区 土木结构面积 砖木结构面积 砖混结构面积 框架结构面积 总面积 县城 9000 114700 171135 325700 620535 乡镇 814770 46954 47636 0 909360 乡镇中的砖木结构房屋为近年新建居住用房，设有构造措施，抗震能力好；老旧砖木房屋未经抗震设防，砌筑工艺和质量较差。震区老旧砖木房屋严重破坏现象为房屋承重墙体大面积剪切裂缝或外闪，局部屋顶塌落；中等破坏现象主要为墙体斜向或竖向开裂，宽度约1mm，但延伸长度较长，由墙体顶部延伸至底部。Ⅶ度区严重破坏以上达24.3%，但无整体倒塌房屋，该类结构房屋未造成人员伤亡。
砖混结构房屋主要是乡（镇）公用房屋，或县城的居住用房。2010年后建设的砖混结构房屋抗震能力较好，地震后出现一定数量墙体细微开裂。2000年前建设砖混结构办公楼设防烈度低，在地震中造成一定数量严重破坏。震区典型房屋灾害如图 4所示。
框架结构多为2010年后新建办公用房，抗震能力好，未产生结构性破坏，但出现大面积填充墙开裂，修复量大。
通过对震区69个调查点进行均匀抽样调查（孙景江等，2011），最后得到本次地震震区房屋破坏比，如表 2所示。
表 2 震区各类结构房屋的破坏情况Table 2. Statistial results of building damages of various structures in the earthquake area行政区 单位 毁坏 严重破坏 中等破坏 轻微破坏 破坏合计 不具备修复价值 县城 m2 26156 54751 74645 341811 497363 99568 间 1308 2738 3732 17091 24869 4979 户 327 685 933 4273 6218 1245 乡镇 m2 80958 151041 194444 260836 687279 280610 间 4048 7552 9722 13042 34364 14031 户 1012 1888 2431 3261 8592 3508 2.2 基础设施与生命线工程的破坏
在市政设施方面，地震造成42km供排水管道、34km供暖管道和4座供热站受损，造成供水厂、污水厂氧化池及构筑物受损；交通系统方面，共86km道路损毁，350m隧道严重损坏，9座桥涵、135处涵洞和20km道路防护损坏。塔什库尔干县村庄道路局部塌陷情况如图 5所示。
在水利系统方面，地震造成114km水渠受损，对震区农作物灌溉造成一定影响。在达布达尔乡，草场水渠破坏造成库什吾尼可尔村、恰特尔塔什村、土拉村及库什吾尼可尔村等30余户、约0.12km2棉作地灌溉受影响，对震区居民的收入造成一定影响，恢复时间需要1个月左右。塔什库尔干县达布达尔乡阿特加依里村草场水渠地基失稳，在地震作用下完全破坏（图 6）。
图 6 布达尔乡阿特加依里村草场水渠地基破坏情况Figure 6. Damage of water channel at Atejiayili village of Dabudaer town在电力系统方面，地震造成63座（110kV、220kV）塔基局部受损，1座110kV变电所和11座35kV变电所受损。电力系统破坏造成塔什库尔干乡布依阿勒村、加隆且特村和吐尔得库勒村等近60户居民用电中断，经过5天的抢修，断电区域已经恢复供电。
在通讯系统方面，地震造成移动、联通、电信公司7个核心机房受损（图 7），3座铁塔损坏，5km光缆倒伏。通讯系统破坏造成塔什库尔县城辖区内部分居民通信不稳定、少数居民固定电话不能呼入进户，经过2天的紧急抢修，通讯基本恢复正常。
图 7 塔什库尔干县通信系统中心机房横向裂缝Figure 7. Transverse cracks in communication room of Taxkorgan county3. 安居富民工程减灾效益分析
在近年来新疆发生的历次破坏性地震中，震区建设的安居富民工程和抗震安居房（张勇，2005）在减少人员伤亡和经济损失中发挥了显著作用（谭明等，2014）。在此次地震中，塔什库尔干县绝大部分倒塌的房屋为老旧的简易房，造成人员伤亡的房屋均为土石木房屋，宏观震中附近的安居富民房屋均完好。新疆大规模实施安居富民工程后，建造的农居符合设计规范要求，无一受到毁坏或者严重破坏，抗震性能得到检验（唐丽华等，2016）。
震后通过对灾区进行抽样调查，统计了安居房及安居工程改造的土木、砖木及砖混结构房屋所占比例，并根据灾区安居房未改造前的土木、砖木及砖混结构房屋数据，结合新疆地区安居富民房震害矩阵，对塔什库尔干5.5级地震的减灾效益进行了计算和分析。在地震灾害损失评估中，将抗震安居房和安居富民房面积替换为改造前的简易房屋面积，计算抗震房减灾效益（刘军等，2016），具体数据见表 3。
表 3 塔什库尔干县震区安居富民房减灾效益对比Table 3. Statistical results of reducing damage with anti-seismic living room project类别 未进行安居工程改造损失 实际损失 减少损失 受伤人数 68 31 37 死亡人数 34 8 26 受灾人数 53438 26486 26952 房屋直接经济损失/亿元 68.8 20.05 38.75 需紧急安置人数 36783 16194 20589 恢复重建费用/亿元 88.3 29.34 42.96 4. 结语
（1）此次地震属于浅源中强地震，震源深度8km，地面振动强。极震区位于地震断裂上方，灾害破坏较集中，对震中附近的库孜滚村造成了毁灭性破坏，与同级别地震相比较灾情较重。
（2）塔什库尔干盆地是1个冰碛堆积盆地，其地下沉积物质具有强烈的不均一性，这种分选性极差的场地地基条件对地震动有一定放大效应，因此造成县城城区内的砖混结构房屋出现了不同程度的破坏，多数框架结构房屋填充墙大面积开裂。
（3）灾区位于帕米尔高原，自然条件恶劣，资源匮乏，经济落后，建设成本高，自建房屋质量差，抗震能力低，也是本次地震震级不大、震害较重的重要原因。
-
图 3 地脉动H/V谱比曲线稳定所需时长
Figure 3. The stabilization time of microtremor H/V different frequency points
图 4 地震记录S波段、P波段、尾波段、全波段和地脉动H/V谱比较
Figure 4. Comparison of H/V spectra of S waves, P waves, coda, full waves and microtremors
图 6 综合考虑谱形差异和卓越频率差异的经验修正系数
Figure 6. Empirical correction coefficient based on spectral shape differences and predominant frequency differences
图 7 MHVR、EHVR与pHVR比较
Figure 7. Comparison of microtremor H/V, measured earthquake H/V and predicted earthquake H/V
-
陈棋福, 刘澜波, 王伟君等, 2008. 利用地脉动探测北京城区的地震动场地响应. 科学通报, 53（18）: 2229—2235. 郭明珠, 谢礼立, 高尔根等, 1999. 利用地脉动进行场地反应分析研究综述. 世界地震工程, 15（3）: 14—19Guo M. Z. , Xie L. L. , Gao E. G. , et al. , 1999. Review for analysis of site response by microtremos. World Information on Earthquake Engineering, 15(3): 14—19. (in Chinese) 郭明珠, 宋泽清, 2000. 论地脉动场地动力特性分析中的Nakamura方法. 世界地震工程, 16（2）: 88—92Guo M. Z. , Song Z. Q. , 2000. On Nakamura method for dynamic analysis of ground characteristics with microtremor. World Information on Earthquake Engineering, 16(2): 88—92. (in Chinese) 李红光, 冷崴, 2019. H/V谱比方法在强震记录应用中几个影响因素的初步研究. 地震工程学报, 41（1）: 169—176Li H. G. , Leng W. , 2019. Preliminary study on the influencing factors of the horizontal-to-vertical spectral ratio applied in strong motion. China Earthquake Engineering Journal, 41(1): 169—176. (in Chinese) 李红玉, 于湘伟, 章文波, 2017. 场地效应研究的传递台站线性反演法——以芦山地震为例. 地球物理学报, 60（1）: 198—211Li H. Y. , Yu X. W. , Zhang W. B. , 2017. Site effects estimation by transfer-station generalized inversion method——A case study of the Lushan earthquake. Chinese Journal of Geophysics, 60(1): 198—211. (in Chinese) 李小军, 2001. 对近年大震震害现象与工程地震问题研究的思考. 国际地震动态, （8）: 26—32Li X. J. , 2001. Comments on phenomena of great earthquake damages and engineering seismological researches. Recent Developments in World Seismology, (8): 26—32. (in Chinese) 李小军, 李娜, 王巨科等, 2021. 场地地震动水平/竖向谱比与地表/基底谱比差异及修正水平/竖向谱比法研究. 震灾防御技术, 16（1）: 81—90Li X. J. , Li N. , Wang J. K. , et al. , 2021. Difference between horizontal-to-vertical spectral ratio and surface-to-bedrock spectral ratio of strong-motion and modified horizontal-to-vertical spectral ratio method. Technology for Earthquake Disaster Prevention, 16(1): 81—90. (in Chinese) 林国良, 张潜, 崔建文等, 2019. 利用地脉动HVSR研究2014年鲁甸6.5级地震场地效应. 地震研究, 42（4）: 531—537Lin G. L. , Zhang Q. , Cui J. W. , et al. , 2019. Determining the site effects of the 2014 Ludian MS6.5 earthquake using HVSR microtremor method. Seismological Research, 42(4): 531—537. (in Chinese) 卢滔, 周正华, 周雍年等, 2006. 关于Nakamura方法有效性的讨论. 地震工程与工程振动, 26（1）: 43—48Lu T. , Zhou Z. H. , Zhou Y. N. , et al. , 2006. Discussion on validation of Nakamura's technique. Earthquake Engineering and Engineering Vibration, 26(1): 43—48. (in Chinese) 罗桂纯, 2015. 基于强震动观测的场地效应非线性特性研究. 北京: 中国地震局地球物理研究所.Luo G. C. , 2015. Study on nonlinearity of site effect using strong-motion recordings. Beijing: Institute of Geophysics, China Earthquake Administration. (in Chinese) 马强, 金星, 李山有等, 2013. 用于地震预警的P波震相到时自动拾取. 地球物理学报, 56（7）: 2313—2321Ma Q. , Jin X. , Li S. Y. , et al. , 2013. Automatic P-arrival detection for earthquake early warning. Chinese Journal of Geophysics, 56(7): 2313—2321. (in Chinese) 欧阳行艳, 章文波, 2008. 利用强震记录进行场地反应分析研究综述. 世界地震工程, 24（3）: 118—126Ouyang X. Y. , Zhang W. B. , 2008. Review on analysis of site effect by strong motion data. World Earthquake Engineering, 24(3): 118—126. (in Chinese) 荣棉水, 李小军, 王振明等, 2016. HVSR方法用于地震作用下场地效应分析的适用性研究. 地球物理学报, 59（8）: 2878—2891Rong M. S. , Li X. J. , Wang Z. M. , et al. , 2016. Applicability of HVSR in analysis of site-effects caused by earthquakes. Chinese Journal of Geophysics, 59(8): 2878—2891. (in Chinese) 师黎静, 陈盛扬, 2020. 基于地脉动单点谱比的场地特征参数测定方法适用性研究. 振动与冲击, 39（11）: 138—145Shi L. J. , Chen S. Y. , 2020. The applicability of site characteristic parameters measurement based on micro-tremor's H/V spectra. Journal of Vibration and Shock, 39(11): 138—145. (in Chinese) 陶夏新, 刘曾武, 郭明珠等, 2001. 工程场地条件评定中的地脉动研究. 地震工程与工程振动, 21（4）: 18—23Tao X. X. , Liu Z. W. , Guo M. Z. , et al. , 2001. A review of microtremor study in engineering site rating. Earthquake Engineering and Engineering Vibration, 21(4): 18—23. (in Chinese) 姚鑫鑫, 任叶飞, 温瑞智等, 2019. 强震动记录H/V谱比法计算处理的若干关键环节. 震灾防御技术, 14（4）: 719—730Yao X. X. , Ren Y. F. , Wen R. Z. , et al. , 2019. Some technical notes on the data processing of the spectral ratio based on the strong-motion records. Technology for Earthquake Disaster Prevention, 14(4): 719—730. (in Chinese) 章文波, 周雍年, 谢礼立, 2001. 场地放大效应的估计. 地震工程与工程振动, 21（4）: 1—9Zhang W. B. , Zhou Y. N. , Xie L. L. , 2001. Estimation of absolute site amplification. Earthquake Engineering and Engineering Vibration, 21(4): 1—9. (in Chinese) Bard P. Y. , 2008. Foreword: the H/V technique: capabilities and limitations based on the results of the SESAME project. Bulletin of Earthquake Engineering, 6(1): 1—2. doi: 10.1007/s10518-008-9059-4 Bignardi S. , Yezzi A. J. , Fiussello S. , et al. , 2018. OpenHVSR-Processing toolkit: enhanced HVSR processing of distributed microtremor measurements and spatial variation of their informative content. Computers & Geosciences, 120: 10—20. Bonnefoy-Claudet S., Cotton F., Bard P. Y., 2006. The nature of noise wavefield and its applications for site effects studies: a literature review. Earth-Science Reviews, 79(3—4): 205—227. Cultrera G. , De Rubeis V. , Theodoulidis N. , et al. , 2014. Statistical correlation of earthquake and ambient noise spectral ratios. Bulletin of Earthquake Engineering, 12(4): 1493—1514. SESAME European project, 2004. Guidelines for the implementation of the H/V spectral ratio technique on ambient vibrations measurements, processing and interpretation: report of European Commission Research General Directorate. Brussels: European Commission. Field E. H. , Jacob K. H. , 1995. A comparison and test of various site-response estimation techniques, including three that are not reference-site dependent. Bulletin of the seismological society of America, 85(4): 1127—1143. Haile M. , Seo K. , Kurita K. , et al. , 1997. Study of site effects in kobe area using microtremors. Journal of Physics of the Earth, 45(2): 121—133. doi: 10.4294/jpe1952.45.121 Ibs-von Seht M. , Wohlenberg J. , 1999. Microtremor measurements used to map thickness of soft sediments. Bulletin of the Seismological Society of America, 89(1): 250—259. doi: 10.1785/BSSA0890010250 Kawase H. , Nagashima F. , Nakano K. , et al. , 2019. Direct evaluation of S-wave amplification factors from microtremor H/V ratios: double empirical corrections to “Nakamura” method. Soil Dynamics and Earthquake Engineering, 126: 105067. doi: 10.1016/j.soildyn.2018.01.049 Konno K. , Ohmachi T. , 1998. Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of microtremor. Bulletin of the Seismological Society of America, 88(1): 228—241. doi: 10.1785/BSSA0880010228 Lermo J. , Chávez-García F. J. , 1993. Site effect evaluation using spectral ratios with only one station. Bulletin of the Seismological Society of America, 83(5): 1574—1594. doi: 10.1785/BSSA0830051574 Lermo J., Chávez-García F. J., 1994 a. Are microtremors useful in site response evaluation? Bulletin of the Seismological Society of America, 84(5): 1350—1364. Lermo J. , Chávez-García F. J. , 1994 b. Site effect evaluation at Mexico city: dominant period and relative amplification from strong motion and microtremor records. Soil Dynamics and Earthquake Engineering, 13(6): 413—423. doi: 10.1016/0267-7261(94)90012-4 Mucciarelli M., Gallipoli M. R., 2001. A critical review of 10 years of microtremor HVSR technique. Bollettino di Geofisica Teorica ed Applicata, 42(3—4): 255—266. Nakamura Y., 1989. A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Railway Technical Research Institute, Quarterly Reports, 30(1): 25—33. Nakamura Y., 2019. What is the nakamura method? Seismological Research Letters, 90(4): 1437—1443. Satoh T. , Kawase H. , Matsushima S. , 2001. Differences between site characteristics obtained from microtremors, S-waves, P-waves, and codas. Bulletin of the Seismological Society of America, 91(2): 313—334. doi: 10.1785/0119990149 Seed H. B. , Idriss I. M. , 1969. Influence of soil conditions on ground motions during earthquakes. Journal of the Soil Mechanics and Foundations Division, 95(1): 99—137. doi: 10.1061/JSFEAQ.0001260 Upadhayay B. , Mori S. , 2013. Comparison of H/V ratios and predominant frequencies inferred from microtremors and earthquake motions in ehime prefecture, Japan. International Journal of Landslide and Environment, 1(1): 117—118. Xu R. , Wang L. M. , 2021. The horizontal-to-vertical spectral ratio and its applications. EURASIP Journal on Advances in Signal Processing, 2021: 75. doi: 10.1186/s13634-021-00765-z -
下载:
下载:
点击查看大图
计量
- 文章访问数: 164
- HTML全文浏览量: 21
- PDF下载量: 24
- 被引次数: 0
