Observed travel-time curves by stacking records from Lanzhou small aperture seismic array
-
摘要: 使用兰州小孔径地震台阵记录的近10年地震观测垂直分量波形数据， 采用长、 短时间平均数比值方法(STA/LTA)叠加出适用于青藏高原东北缘地区的观测走时曲线. 结果表明， 兰州小孔径地震台阵独特的地理位置， 基本上能够记录到不同震中距(0—180°)和不同方位的地震事件； 使用不同频率滤波处理之后的地震数据记录叠加出相应的观测走时曲线， 从观测走时曲线中可以识别出不同体波震相(P， PKIKP， PKP， PP， PPP， PcP， ScP， S， SS等)的到时及其观测走时曲线特征. 这对识别地震各种震相， 认识和研究地球内部精细结构等具有非常重要的科学意义.Abstract: Based on the past 10 years vertical observational waveforms recorded by Lanzhou small aperture seismic array, which is located at the northeastern margin of the Tibetan Plateau, the observed travel-time curves are obtained by using STA/LTA technique. The results show that Lanzhou seismic array can record events with different epicentral distances between 0—180° and with different azimuths due to its particular location. By stacking the seismic records filtered with different frequencies, we can get corresponding observed travel-time curves so as to identify different body wave phases, such as P, PKIKP, PKP, PP, PPP, PcP, ScP, S, SS, and to study the characteristics of travel-time curves. It is of scientific significance to identify various seismic phases and further study fine structure of the Earth’s interior.
-
引言
泥火山是在特定地质构造及水文地质环境下产生的一种构造流体地质现象(王道，2000)． 泥火山的出现通常与生油气带有关，某些地区的地表或近地表有天然气或者地下水，而地层是比较松软的泥岩，当地层受到地下压力时，比较松软的物质就会沿着裂缝或断层上升，并在地下水或天然气的挟带下穿透地表的空隙喷出，将携带的泥沙等溅落到地面，就形成泥火山(刘嘉麒，2003)． 由于泥火山喷发时可以将大量有价值的信息带到地表，因此，不少科学家将泥火山称为深度可达12 km的“天赐钻井”(邦华，1990； 樊祺诚等，1998； 栗周熊，2003)．
泥火山喷发是对内部大量气体聚集引起异常高压的释放，其喷发的剧烈程度取决于内部气体的多少或压力的高低． 泥浆和气体聚合到足够的压力才冲出地面，使得绝大多数泥火山为间歇式喷发，在没有地震事件影响时，泥火山有较为稳定的喷发周期，当受到地震等事件影响时喷发周期为不定期(刘若新等，1995； 樊祺诚等，2005； 高小其等，2009； 许建东等，2011)． 台湾屏东泥火山平均每年爆发一次，但有地震发生时，出现了爆发的频率增加、 地点移动、 面积扩大等现象，一年喷发可达4次(Shih，1967； Gieskes et al，1992； 马骏，1999； 叶高华，2003)．
类似泥火山喷发具有较好地震前兆指示意义的研究近年来也显著增加(王道等，1997； Yang et al，2003； 高小其，2009)． 例如1995年5月，克拉玛依市独山子矿区泥火山重新活动(图 1)，5处泥火山口喷出大量稠泥浆、 水和天然气. 喷发时泥浆柱高约20—50 cm，由泥浆形成的泥火山锥的直径达5—6 m，高度大于2 m，泥浆中混有岩石碎块，经沉淀后水样分析，其地下水为高矿化(14.7 g/L)的氯-钠质水． 气体成分主要为甲烷、 氮和二氧化碳等，有的泥火山出口还可见少量油花． 此泥火山自20世纪70年代以来就很少这样活动过，这次喷发一直持续到1995年10月底． 到1996年1月9日，首先发生了沙湾MS5.2地震，而后的3年相继出现了11次MS6.0—6.9、 2次MS7.0—7.9地震(含玛尼地震)的活跃期． 从上面列举的现象可以初步肯定，泥火山喷发和地震的发生都与区域地壳构造应力有关．
![]()
图 1 1995年5—10月新疆独山子泥火山异常喷发照片Figure 1. Photo of Dushanzi mud volcano eruption from May to October of 19951. 艾其沟泥火山观测点简介
为系统捕捉泥火山的地震前兆异常，基于网络技术的视频监控服务，在中国地震局星火科技攻关项目资助下，2011年8月底建成了包含有3个监测点的新疆北天山地区泥火山实时监测系统．
新疆乌苏艾其沟泥火山为3个24小时连续观测点之一(图 2)． 该泥火山是高约8 m、 相距10 m的双泥火山锥，地处准噶尔盆地南缘，在构造上属于北天山山前坳陷的艾其沟背斜及背斜北翼的艾其沟断裂带． 两座泥火山锥中，南侧的泥火山锥喷口直径在1 m左右，北侧的喷口直径近2 m； 观测点设在北侧较大的喷口(44.2°N，84.5°E)． 从该泥火山所喷发的气体、 水质等成分来看，气、 液都来源于地壳深部．
2. MS6.0和MS6.6两次震例现象
新疆乌苏艾其沟泥火山观测点自2011年8月底开始观测以来，监测点周围200 km范围内，于2011年10月16日在新疆精河发生了MS5.0地震(震中距144 km). 两周后的2011年11月1日在新疆尼勒克发生了MS6.0地震(震中距170 km)，2012年6月30日在新疆新源—和静又发生了MS6.6地震(震中距92 km)． 两次MS≥6.0地震前，新疆艾其沟泥火山液面都出现了明显的“背景值—上升—转折—下降—背景值”的宏观异常变化现象. 图 3为乌苏艾其沟泥火山观测点与地震震中分布图．
![]()
图 3 新疆艾其沟泥火山口与三次地震震中分布图Figure 3. Mud volcanic vent at Aiqi valley and distribution of three earthquakes2.1 2011年11月1日新疆尼勒克MS6.0地震
自2011年9月22日起，乌苏艾其沟泥火山观测点的泥浆出现明显外溢现象，这种溢出变化一直持续到9月30日； 进入10月以后，艾其沟泥火山观测点的液面开始逐渐下降，在持续下降过程中，于2011年10月16日在新疆精河发生了MS5.0地震(44.3°N，82.7°E). 15天后于2011年11月1日在新疆伊犁哈萨克自治州尼勒克县、 伊宁县、 巩留县交界处(43.6°N，82.4°E)发生了MS6.0地震，震源深度为28 km(图 3)． 直至11月底，泥火山液面缓慢恢复至背景值(图 4)． 图 5为艾其沟泥火山液面在新疆尼勒克MS6.0地震前后的变化示意图．
![]()
图 4 新疆尼勒克MS6.0地震前后乌苏艾其沟泥火山口监控视频截图Figure 4. Snapshot of Aiqi valley mud volcanic vent before and after Nilka-Gongliu MS6.0 earthquake![]()
图 5 新疆尼勒克2011年MS6.0地震前后乌苏艾其沟泥火山液面变化示意图Figure 5. Dynamic change in liquid level of Aiai valley mud volcano before and after Nilka-Gongliu MS6.0 earthquake in the year 20112.2 2012年6月30日新疆新源—和静MS6.6地震
自2011年12月10日起，艾其沟泥火山观测点的泥浆再次出现明显外溢现象，这种溢出变化一直持续到2012年3月. 进入2012年4月以后，艾其沟泥火山观测点的液面又开始逐渐下降，在持续下降过程中，于2012年6月30日在距离艾其沟泥火山观测点92 km的新疆伊犁哈萨克自治州新源县与巴音郭楞蒙古自治州和静县交界处(43.4°N，84.8°E)发生了MS6.6地震(图 5，6). 目前，该泥火山液面仍未恢复到背景值． 该MS6.6地震前后，艾其沟泥火山口液面又再次出现了上述MS6.0地震前后的类似变化(图 7)．
![]()
图 6 新疆新源-和静MS6.6地震前后乌苏艾其沟泥火山口监控视频截图Figure 6. Snapshot of Aiqi valley mud volcano vent before and after MS6.6 Xinyuan-Hejing earthquake![]()
图 7 新疆新源-和静MS6.6地震前后乌苏艾其沟泥火山液面变化示意图Figure 7. Dynamic change in liquid level of Aiqi valley mud volcano before and after Xinyuan-Hejing MS6.6 earthquake3. 讨论
1)推进泥火山化学项目观测. 地震监测预报系统的发展，离不开科技创新，离不开新技术新方法的开发与应用． 目前，我们对地震前兆的观测基本上都在数米到数百米深，仅少数地震地下流体观测井可达数千米． 某种程度上讲，我们对地球内部的探测技术仍十分有限，只能利用在地表或近地表观测到的物理量，推测震源所在深度的应力应变场． 然而，地表的应力应变场在多大程度上反映深部的应力应变场我们还不得而知(马瑾，马胜利，1995； 李霓等，2000)． 而泥火山喷发时可以将深达12 km的大量有价值的信息带到地表，因此，泥火山的系统观测可以在一定程度上弥补地震前兆观测中深度不够的问题．
考虑到新疆艾其沟泥火山仅仅是新疆北天山地区诸如克拉玛依市独山子泥火山、 沙湾县霍尔果斯泥火山、 乌苏市白杨沟泥火山群、 温泉县天泉泥火山等众多泥火山中的一个点，并且这些泥火山都分布在全国2010—2020年全国地震重点监视防御区内及邻近地区. 从更长远的角度来看，对这类来自地层相对深处的物质增加自动化化学观测项目(如气体总量、 He、 Rn、 CH4等)，并建点成网，这样既可以填补中国大陆地区现今泥火山地球化学系统观测和研究的空白，又对认识泥火山的成因、 拓展具有新疆特色的地震监测和预测新方法等都有重要意义．
2)泥火山活动与地震发生关系的复杂性. 泥火山活动通常情况下具有一定的周期性. 当区域地壳构造应力不断增强时，封闭构造中岩石的孔隙压力逐渐增大，直到压力达到泥火山喷发条件时，泥火山就打破原有周期而重新活动． 同样，这种区域地壳构造应力的不断增强也能导致地震的孕育和发生，因此可以说二者具有同源关系． 地震活动也可以伴随泥火山的活动，大规模的泥火山活动同样可以引发小地震． 但是，我们也同时观测到，新疆艾其沟泥火山在两次MS≥6.0地震前后的变化与已掌握的新疆霍尔果斯泥火山(地质构造上位于北天山山前坳陷带霍尔果斯背斜轴部)震例异常现象有所区别(高小其等，2008)．
我们依据2011年11月1日新疆伊宁县MS6.0地震前后艾其沟泥火山的宏观异常经验，对2012年6月30日MS6.6地震进行了预测. 预测结果表明，预测的地点和震级都正确，但预测时间过期了2.5个月． 该预测结果说明泥火山活动与地震的关系可能相当复杂，目前，很难用某一种简单的模式对其加以总结或科学的解释．
总之，新疆艾其沟泥火山在两次MS≥6.0地震前的类似异常变化，进一步验证了泥火山显著活动期也是泥火山周围一定距离范围内中强以上地震的优势发震时间段. 但是，新疆艾其沟泥火山与霍尔果斯泥火山震例异常现象的不同，显示了泥火山的震兆异常特征可能与其构造条件密切相关．
-
![]()
图 1 兰州小孔径地震台阵示意图.三角形表示台阵各子台的位置
Figure 1. Sketch of Lanzhou small aperture seismic array, where the triangle denotes station and the rectangle denotes the observatory
![]()
图 2 地震震中分布图. 三角形表示台阵的位置， 小圆点表示震中
Figure 2. Epicentral distribution used in this study, where the triangle represents the Lanzhou seismic array and the dots are epicenters
![]()
图 3 不同震中距与地震数目分布图
Figure 3. Histogram of number of events in 0.5° epicentral distance bin for the earthquakes used in this study
![]()
图 4 (a)兰州小孔径台阵10 s低通滤波之后的记录(2010年7月4日日本本州东海岸附近地震，震中距为30.6°，MS6.3). STACK为9个子台的平均叠加结果;(b)STACK作STA/LTA的计算结果;(c)该地震事件0.5°震中距范围内所有地震事件作STA/LTA的叠加平均结果
Figure 4. (a) Records of an earthquake (MS6.3, epicentral distance 30.6°, 2010-07-04) in Coast of Honshu, Japan, on nine stations of Lanzhou seismic array, which were filtered by 10s low-pass filter; (b) The computed STA/LTA for the STACK trace; (c) The average STA/LTA resulting from all traces within a 0.5° distance range
![]()
图 5 垂直分量2 s高通滤波叠加的观测走时曲线(a)和可识别震相IASP91模型理论走时曲线(b)
Figure 5. Vertical stacking of travel-time curves after a 2 s high-pass filter was applied to each trace (a) and travel times of the identifiable seismic phases based on the IASP91 velocity model (b)
![]()
图 6 垂直分量6 s高通滤波叠加的观测走时曲线(a)和可识别震相IASP91模型理论走时曲线(b)
Figure 6. Vertical stacking of travel-time curves after a 6 s high-pass filter was applied to each trace (a) and travel times of the identifiable seismic phases based on the IASP91 velocity model (b)
![]()
图 7 垂直分量10 s低通滤波叠加的观测走时曲线(a)和可识别震相IASP91模型理论走时曲线(b)
Figure 7. Vertical stacking of travel-time curves after a 10 s low-pass filter was applied to each trace (a) and travel times of the identifiable seismic phases based on the IASP91 velocity model (b)
![]()
图 8 垂直分量30 s低通滤波叠加的观测走时曲线(a)和可识别震相IASP91模型理论走时曲线(b)
Figure 8. Vertical stacking of travel-time curves after a 30 s low-pass filter was applied to each trace (a) and travel times of the identifiable seismic phases based on the IASP91 velocity model (b)
![]()
图 9 垂直分量100 s低通滤波叠加的观测走时曲线(a)和可识别震相IASP91模型理论走时曲线(b)
Figure 9. Vertical stacking of travel-time curves after a 100 s low-pass filter was applied to each trace (a) and travel times of the identifiable seismic phases based on the IASP91 velocity model (b)
![]()
图 10 垂直分量5种不同滤波综合叠加的观测走时曲线(a)和可识别震相IASP91模型理论走时曲线(b)
Figure 10. Vertical stacking result after five different filters were applied to each trace (a) and travel times calculated for the IASP91 velocity model showing the identifiable seismic phases (b)
![]()
图 11 兰州小孔径地震台阵极远震记录(2003年 1月22日墨西哥地震， 震中距119.5°， MS7.5). STACK为台阵9个子台叠加记录， 记录中可以识别Pdiff, PKIKP, PP, PPP等震相
Figure 11. The far most earthquake recordings of a MS=7.5 event that occurred in Mexico on January 22, 2003, on Lanzhou seismic array. This event is located 119.5° away. These traces have been low-pass filtered at 10 s and the arrivals of a few phases (Pdiff， PKIKP， PP, PPP) are identified
![]()
图 12 (a) 台湾MS5.7地震(东南方向)和新疆MS6.0地震(西北方向)相同震中距(19.2°)观测波形对比； (b) 相同深度(30 km±3 km)、 不同方位(75°和200°)、 相近震中距(30.5°和31.0°)多次地震叠加的P波观测走时与理论走时对比； (c)相同深度(均为10 km)、 不同方位(136°和250°)、 相同震中距(49.0°)多次地震叠加的P波观测走时与理论走时的对比． 图b， c中垂直黑线为叠加的P波观测走时， 蓝线为理论走时
Figure 12. (a) Comparison of earthquake waveforms with the same epicentral distance (19.2°) between Taiwan MS5.7 (southeast) and Xinjiang MS6.0 (northwest) earthquakes. (b) Comparison of observed travel times with theoretical ones of the multiple earthquakes stacking of P-wave traces. The earthquakes have the same depth (30 km±3 km), the different back azimuth (75° and 200°) and the nearly similar distance (30.5° and 31.0°). (c) Comparison of observed P-wave travel times with theoretical ones of the multiple earthquakes stacking. The earthquakes have the same focal depth (10 km), the different back azimuth (136° and 250°) and the same distance (49.0°). In Figs.2b, c, the black vertical line is the observed stacked travel time for P wave and the blue one is the corresponding theoretical travel time
表 1 不同滤波和STA/LTA参数表
Table 1 Filtering and STA/LTA parameters
滤波 STA/s LTA/s 2 s高通 1 9 6 s高通 2 20 10 s低通 3 30 30 s低通 5 45 100 s低通 10 90 
下载: 导出CSV
-
秦满忠, 张元生, 沈旭章, 魏从信. 2012. 小孔径台阵垂直分量P波接收函数研究及应用[J]. 地震学报, 34 (1): 44-51. Qin M Z, Zhang Y S, Shen X Z, Wei C X. 2012. A study on vertical-component P-wave receiver function of small-aperture seismic array and its application[J]. Acta Seismologica Sinica, 34 (1): 44-51 (in Chinese)
沈旭章, 周蕙兰. 2009. 用PKIKP的前驱震相探测青藏高原东部地幔底部的散射体[J]. 科学通报, 54 (24): 3844-3851. Shen X Z, Zhou H L. 2009. Locating seismic scatterers at the base of the mantle beneath eastern Tibet with PKIKP precursors[J]. Chinese Sci Bull, 54 (24): 3844-3851 (in Chinese)
沈旭章, 梅秀苹, 张淑珍, 秦满忠. 2010. 兰州台阵响应函数及不同方位地震事件FK分析结果[J]. 西北地震学报, 32 (1): 59-64. Shen X Z, Mei X P, Zhang S Z, Qin M Z. 2010. The array response function of Lanzhou seismic array and results of FK analysis for small earthquakes in different azimuths[J]. Northwestern Seismological Journal, 32 (1): 59-64 (in Chinese).
周元泽, 王卓君. 2011. 伊豆-小笠原和汤加地区中地幔多层速度结构差异性研究[J]. 中国科学: D辑, 41 (7): 936-944 Zhou Y Z, Wang Z J. 2011. On the difference of mid-mantle multiple velocity structure beneath Izu-Bonin and Tonga areas[J]. Science in China: Series D, 41 (7): 936-944 (in Chinese)
中国地震局监测预报司. 2007. 地震学与地震观测[M]. 北京: 地震出版社: 172 Department of Earthquake Monitoring and Predition, China Earthquake Administration. 2007. Seismology and Earthquake Observation[M]. Beijing: Seismological Press: 172 (in Chinese)
Astiz L, Earle P, Shearer P. 1996. Global stacking of broadband seismograms[J]. Seismol Res Lett, 67 (4): 8-18.
Earle P S, Shearer P M. 1994. Characterization of global seismograms using an automatic-picking algorithm[J]. Bull Seismol Soc Am, 84 (2): 366-376.
Gaia S, Lapo B, Antonio P. 2006. Global seismic tomography and modern parallel computers[J]. Annals of Geophysics, 49 (4/5): 977-986
Rost S, Thomas C. 2002. Array seismology: Methods and applications[J]. Rev Geophys, 40 (3): 2-1-2-27. doi:10.1029/2000RG000100
Rost S, Thorne M S, Garnero E J. 2006. Imaging global seismic phase arrivals by stacking array processed short-period data[J]. Seismol Res Lett , 77 (6): 697-707
Shearer P M. 1991. Imaging global body wave phases by stacking long-period seismograms[J]. J Geophys Res, 96 (B12): 20353-20364
Shearer P M. 1994. Imaging Earth's seismic response at long periods[J]. Eos Trans AGU, 75 (39): 449-452
Shen X, Ritter J R R. 2010. Small-scale heterogeneities below the Lanzhou CTBTO seismic array, from seismic wavefield fluctuations[J]. J Seismol, 14 (3): 481-493
Walck M C, Clayton R W. 1984. Analysis of upper mantle structure using wave field continuation of P waves[J]. Bull Seismol Soc Am, 74 (5): 1703-1719
-
期刊类型引用(3)
1. 纪春玲，张合，董博，马广庆. 河北红山台宽频带地震记录区域叠加研究. 震灾防御技术. 2023(02): 380-388 . 
百度学术
2. 白永福，高翔，张伟，奚冲霄，王凌. 海拉尔台阵数据功率谱分析. 防灾减灾学报. 2020(02): 49-55 . 百度学术
3. 秦满忠，李顺成，刘旭宙，张淑珍，张元生. 核幔界面衍射波P_(diff)震相最大记录震中距探讨. 地震学报. 2015(02): 249-256 . 本站查看
其他类型引用(0)
计量
- 文章访问数: 627
- HTML全文浏览量: 236
- PDF下载量: 20
- 被引次数: 3
