The Preliminary Analysis of Airguns Seismic Data in the MS6.2 Taiwan Strait Earthquake Region
-
摘要: 对2019年在台湾海峡6.2级地震震区布设的N01、NE02测线地震测深剖面的Pg波走时进行层析成像反演,获得测线下方地壳上部二维速度结构。对N01测线单道反射地震测深剖面进行多次波衰减等处理,并与Pg波成像结果进行对比。研究结果表明,采用走时层析成像方法与单道反射地震测深剖面获得沿探测剖面沉积层上地壳基底形态特征等具有较好的对应关系。由于测线穿越多个地质构造单元和多条断裂带,走时层析成像和单道反射地震测深剖面结果综合显示研究区结晶基底面起伏较大,沉积层速度和厚度变化较剧烈,受台湾海峡西部新生代构造活动影响,显示了相应的断裂或不同地质构造单元在上部地壳内的结构特征。Abstract: The Pg wave travel time of N01 and NE02 deep seismic sounding profiles laid in the MS6.2 Taiwan Strait earthquake region in 2019 was inverted by tomography, and the two-dimensional velocity structure of the upper crust below the survey lines was obtained. The N01 single seismic section was processed by multiple attenuation, and compared with the Pg wave imaging results, there is a good correspondence between the morphological characteristics of the basement. As the survey lines pass through multiple geological tectonic units and multiple fault zones, the results of tomography and single seismic section show that the crystalline basement surface in the study area fluctuates greatly, and the velocity and thickness of the sedimentary layer changes violently. Affected by cenozoic tectonic activities in the western Taiwan Strait, it shows the structural characteristics of the corresponding faults or different geological tectonic units in the upper crust.
-
Key words:
- Taiwan strait /
- Deep seismic exploration /
- Basement structure /
- Seismic profile
-
引言
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)。
在电力系统方面,地震造成63座(110kV、220kV)塔基局部受损,1座110kV变电所和11座35kV变电所受损。电力系统破坏造成塔什库尔干乡布依阿勒村、加隆且特村和吐尔得库勒村等近60户居民用电中断,经过5天的抢修,断电区域已经恢复供电。
在通讯系统方面,地震造成移动、联通、电信公司7个核心机房受损(图 7),3座铁塔损坏,5km光缆倒伏。通讯系统破坏造成塔什库尔县城辖区内部分居民通信不稳定、少数居民固定电话不能呼入进户,经过2天的紧急抢修,通讯基本恢复正常。
3. 安居富民工程减灾效益分析
在近年来新疆发生的历次破坏性地震中,震区建设的安居富民工程和抗震安居房(张勇,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)灾区位于帕米尔高原,自然条件恶劣,资源匮乏,经济落后,建设成本高,自建房屋质量差,抗震能力低,也是本次地震震级不大、震害较重的重要原因。
-
[1] 蔡锋, 曹超, 周兴华等, 2013. 中国近海海洋—海底地形地貌. 北京: 海洋出版社. [2] 陈园田, 谢志平, 1996. 台湾海峡的活动断裂与地震活动. 华南地震, 16(1): 57—62.Chen Y. T., Xie Z. P., 1996. Active faults and earthquakes in Taiwan Strait. South China Journal of Seismology, 16(1): 57—62. (in Chinese) [3] 邓起东, 闽伟, 晁洪太等, 2001. 渤海地区新生代构造与地震活动. 见: 卢演俦主编. 新构造与环境. 北京: 地震出版社, 218—233. [4] 丁祥焕, 1999. 福建东南沿海活动断裂与地震. 福州: 福建科学技术出版社. [5] 郭令智, 施央申, 马瑞士, 1983. 西太平洋中、新生代活动大陆边缘和岛弧构造的形成及演化. 地质学报, 57(1): 11—21.Guo L. Z., Shi Y. S., Ma R. S., 1983. On the formation and evolution of the Mesozoic-Cenozoic active continental margin and island arc tectonics of the western Pacific Ocean. Acta Geologica Sinica, 1983, 57(1): 11—21. (in Chinese) [6] 郭令智, 马瑞士, 施央申等, 1998. 论西太平洋活动大陆边缘中—新生代弧后盆地的分类和演化. 成都理工学院学报, 25(2): 134—144.Guo L. Z., Ma R. S., Shi Y. S., et al., 1998. On the generation and evolution of Mesozoic Cenozoic back arc basins of western Pacific active continental margin. Journal of Chengdu University of Technology, 25(2): 134—144. (in Chinese) [7] 黄昭, 王善雄, 2006. 台湾海峡滨海断裂带的构造特征与活动性. 大地测量与地球动力学, 26(3): 16—22.Huang Z., Wang S. X., 2006. Tectonic features and activity of Binhai Fault zone in Taiwan Strait. Journal of Geodesy and Geodynamics, 26(3): 16—22. (in Chinese) [8] 李赶先, 1994. 台湾海峡新生代断裂构造及其形成演化. 热带海洋, 13(2): 17—24.Li G. X., 1994. The formation and evolution of Cenozoic rifts Taiwan Strait. Tropic Oceanology, 13(2): 17—24. (in Chinese) [9] 戚筱俊, 1999. 台湾海峡地形、地质及地震概述. 西部探矿工程, 11(4): 4—6. [10] 单联瑜, 刘连升, 2007. 三维高阶快速步进波前重构旅行时计算精度分析. 石油地球物理勘探, 42(4): 413—417. doi: 10.3321/j.issn:1000-7210.2007.04.010Shan L. Y., Liu L. S., 2007. Analysis on computation precision of wavefront-reconstructing traveltime for 3-D high-order rapid stepping algorithm. Oil Geophysical Prospecting, 42(4): 413—417. (in Chinese) doi: 10.3321/j.issn:1000-7210.2007.04.010 [11] 王本善, 1995. 台湾海峡地区地质构造特征. 南方油气地质, 2(1): 57—64. [12] 王建, 赵明辉, 贺恩远等, 2014. 初至波层析成像的反演参数选取: 以南海中央次海盆三维地震探测数据为例. 热带海洋学报, 33(5): 74—83. doi: 10.3969/j.issn.1009-5470.2014.05.010Wang J., Zhao M. H., He E. Y., et al., 2014. The selection of optimal inversion parameters for first-arrival seismic tomography: an application to 3D seismic data from the central sub-basin of the South China Sea. Journal of Tropical Oceanography, 33(5): 74-83. (in Chinese) doi: 10.3969/j.issn.1009-5470.2014.05.010 [13] 王笋, 丘学林, 方伟华等, 2018. 台湾海峡西南部的海陆联合深地震探测资料特点与处理对策. 热带海洋学报, 37(2): 92—99.Wang S., Qiu X. L., Fang W. H., et al., 2018. Features of the onshore-offshore seismic data in Southwest Taiwan Strait and some countermeasures for data processing. Journal of Tropical Oceanography, 37(2): 92—99. (in Chinese) [14] 熊绍柏, 刘宏兵, 王有学等, 2002. 华南上地壳速度分布与基底、盖层构造研究. 地球物理学报, 45(6): 784—791. doi: 10.3321/j.issn:0001-5733.2002.06.005Xiong S. B., Liu H. B., Wang Y. X., et al., 2002. A study on velocity distribution in upper crust and tectonics of basement and cover in South China. Chinese Journal of Geophysics, 45(6): 784—791. (in Chinese) doi: 10.3321/j.issn:0001-5733.2002.06.005 [15] 闫培, 张艺峰, 李普春等, 2015. 台湾海峡深地震探测: HX9测线试验及初步成果. 中国地震, 31(1): 27—35. doi: 10.3969/j.issn.1001-4683.2015.01.003Yan P., Zhang Y. F., Li P. C., et al., 2015. Deep seismic exploration in the Taiwan Strait: the experiment of HX9 survey line and preliminary results. Earthquake Research in China, 31(1): 27—35. (in Chinese) doi: 10.3969/j.issn.1001-4683.2015.01.003 [16] 杨肖琪, 宋文隆, 陈承惠, 1996. 台湾海峡地质构造特征. 台湾海峡, 15(2): 127—136.Yang X. Q., Song W. L., Chen C. H., 1996. Characteristics of geological structure in Taiwan Strait. Journal of Oceanography in Taiwan Strait, 15(2): 127—136. (in Chinese) [17] 姚道平, 张艺峰, 闫培等, 2016. 台湾海峡大容量气枪震源海陆联测初探. 地震学报, 38(2): 167—178. doi: 10.11939/jass.2016.02.002Yao D. P., Zhang Y. F., Yan P., et al., 2016. A preliminary exploration into onshore-offshore seismic experiment by using large volume air-gun source in Taiwan Strait. Acta Seismologica Sinica, 38(2): 167—178. (in Chinese) doi: 10.11939/jass.2016.02.002 [18] 俞何兴, 陈汝勤, 1996. 台湾海域之沉积盆地. 台北: 台湾编译馆. [19] 俞何兴, 2006. 台湾海域海底地形及其地质意义. 台北: 台湾大学. [20] 詹文欢, 孙宗勋, 唐诚等, 2004. 华南滨海断裂带及其对台湾海峡地震活动的控制作用. 热带海洋学报, 23(4): 19—24. doi: 10.3969/j.issn.1009-5470.2004.04.003Zhan W. H., Sun Z. X., Tang C., et al., 2004. Littoral active fault belt of south China and its control on seismic activity in Taiwan Strait. Journal of Tropical Oceanography, 23(4): 19—24. (in Chinese) doi: 10.3969/j.issn.1009-5470.2004.04.003 [21] 张艺峰, 金星, 郭晓然等, 2017. 台湾海峡地区深部结构探测数据处理及初步结果. 见: 2017中国地球科学联合学术年会论文集. 北京: 中国地球物理学会. [22] Paige C. C., Saunders M. A., 1982. LSQR: An algorithm for sparse linear equations and sparse least squares. ACM Transactions on Mathematical Software, 8(1): 43—71. doi: 10.1145/355984.355989 [23] Rawlinson N., Sambridge M., 2004a. Wave front evolution in strongly heterogeneous layered media using the fast marching method. Geophysical Journal International, 156(3): 631—647. doi: 10.1111/j.1365-246X.2004.02153.x [24] Rawlinson N., Sambridge M., 2004b. Multiple reflection and transmission phases in complex layered media using a multistage fast marching method. Geophysics, 69(5): 1338—1350. doi: 10.1190/1.1801950 [25] Sethian J. A., 1996. Theory, algorithms, and applications of level set methods for propagating interfaces. Acta Numerica, 5: 309—395. doi: 10.1017/S0962492900002671 [26] Sethian J. A., Popovici A. M., 1999. 3-D traveltime computation using the fast marching method. Geophysics, 64(2): 516—523. doi: 10.1190/1.1444558 [27] Zelt C. A., Barton P. J., 1998. Three-dimensional seismic refraction tomography: a comparison of two methods applied to data from the Faeroe Basin. Journal of Geophysical Research: Solid Earth, 103(B4): 7187—7210. doi: 10.1029/97JB03536 -