Characteristics of seismic coda attenuation in Longtan reservoir of Guangxi region
-
摘要:
按照高信噪比和计算要求,从 2006—2013 年广西龙滩水库地震台网的10个子台记录到的龙滩库区近场数字地震波形中挑选出 3 822 条三分向记录,经滤波和消除环境噪声后,采用Sato单次散射模型,计算了地震波传播路径上的尾波Q ( f )值,对Q ( f )与频率f之间的关系进行了拟合,并在此基础上分析讨论了库区尾波衰减的时空特征. 结果显示,龙滩库区总体呈低Q0值和依赖性指数η较高的特点,这两个参数值均低于浙江珊溪水库,与其地震活动程度高于后者的事实相一致。尾波Q值的空间分布显示,周围小震密集分布的台站和历史地震活跃地区的台站的Q0值低于相对于周围小震分布稀疏的台站的Q0值。库区尾波衰减参数与尾波采样体之间的关系显示,采样深度越深,尾波参数Q0值越低,因此认为该区域可能存在深部高衰减层。衰减参数随时间的变化过程表明,蓄水后地震尾波衰减参数仍主要反映了库区原有地壳内部介质的非均匀性和地震波吸收特性,库区地壳介质特性无显著变化,但局部敏感点有所调整。
Abstract:This paper firstly selected 3 822 three-component records from 10 stations of Long-tan reservoir seismic network from 2006 to 2013 according to the high signal-to-noise ratio and calculation requirements. Based on the Sato single scattering model, this paper calculated the coda Q( f ) after filtering the data and eliminating ambient noises, and then obtained the relationship between coda Q( f ) and frequency f. Finally we anayzed the spatio-temporal characteristics of seismic coda attenuation in the Longtan reservoir area. The results show that the Longtan reservoir area is characterized by low Q0 and high dependence index η, and its Q0 value and dependence index η are lower than that of Shanxi Reservoir, which is consistent with the fact that the its seismic activity is higher than the latter. The spatial distribution of the coda Q value shows that Q0 values of the stations with surrounding small earthquakes densely distributed and the stations in the areas with active historical earthquakes are lower than those of the stations with sparse distribution of surrounding small earthquakes. Further analyses on the relationship between attenuation parameters and sampling depth of the reservoir area suggest that there exist deep high attenuation layers beneath the studied area. By analyzing the change process of attenuation parameters with time, it is considered that after reservoir impoundment the coda attenuation parameters still mainly reflect the non-uniformity of the original crustal internal medium of the reservoir area, and the properties of the crustal medium have not been changed significantly, but some local sensitive points have been adjusted.
-
Keywords:
- Longtan reservoir /
- reservoir earthquakes /
- coda /
- attenuation parameter
-
-
图 1 龙滩水库地质构造以及蓄水后地震震中和台站分布图
F1:八腊断裂;F2:天峨—金谷断裂带;F3:望谟—逻西断裂带;F4:乐业—甘田断裂;F5:罗甸—望谟断裂;F6:里落—白塘—浪全断裂带
Figure 1. Geological structure of Longtan reservoir as well as location of epicenters after impoundment and stations distribution
F1:Bala fault;F2:Tiane-Jingu fault zone;F3:Wangmo-Luoxi fault zone;F4:Leye-Gantian fault zone;F5:Luodian-Wangmo fault zone;F6:Liluo-Baitang-Langquan fault zone
图 2 龙滩库区尾波衰减参数计算示例
(a) EW向原始波形记录,t0表示发震时刻,t1表示尾波截断点时刻(t1−t0=60 s);(b) 频率f=11×(1±1/3) Hz时的波形滤波示例;(c) 频率f为4—18 Hz时拟合波形与实际波形的相关系数;(d) Q( f )与频率f的函数关系
Figure 2. Demonstrations of calculation process of coda wave attenuation parameters in Longtan reservoir area
(a) E-W component original waveform record,where t0 is the origin time,t1 denotes cutoff point of coda waves used in the calculation (t0−t1=60 s);(b) Waveform filtered demonstration with f =11×(1±1/3) Hz;(c) Correlation coefficient between fitted waveforms and actual ones when f is 4−18 Hz;(d) Functional relationship between Q( f ) and f
表 1 龙滩水库地震台网各台站尾波衰减参数统计
Table 1 Statistics on coda attenuation parameters of the stations of Longtan reservoir seismic network
台站 台基岩性 Q0 衰减参数 记录数 δQ0 η δη 天峨 灰岩 50.38 9.96 0.85 0.07 270 仁顶 灰岩 49.00 7.55 0.88 0.06 198 坡皇洞 灰岩 53.46 11.05 0.83 0.08 492 甲龙 砂岩 52.83 10.82 0.82 0.07 430 八广 砂岩 55.07 14.95 0.85 0.10 628 里纳 砂岩 52.30 12.09 0.85 0.08 517 凉风坳 砂岩 48.92 10.09 0.86 0.07 354 小苗坡 砂岩 45.67 7.76 0.90 0.07 443 罗甸 灰岩 52.57 12.04 0.86 0.09 212 坪上寨 砂泥岩 48.63 9.73 0.83 0.06 278 -
郭培兰,姚宏,袁媛. 2006. 龙滩水库地震危险性分析[J]. 高原地震,18(4):17–23. doi: 10.3969/j.issn.1005-586X.2006.04.003 Guo P L,Yao H,Yuan Y. 2006. Analysis on potential seismic risk in Longtan reservoir[J]. Earthquake Research in Plateau,18(4):17–23 (in Chinese).
华卫,赵翠萍,陈章立,郑斯华. 2009. 龙滩水库地区P波、S波和尾波衰减[J]. 地震学报,31(6):620–628. doi: 10.3321/j.issn:0253-3782.2009.06.003 Hua W,Zhao C P,Chen Z L,Zheng S H. 2009. Attenuation of P,S and coda waves in Longtan reservoir region[J]. Acta Seismologica Sinica,31(6):620–628 (in Chinese).
梁国经,李仲巍,郑双凤,刘俊清,刘达峰,綦伟. 2011. 龙岗火山区尾波Q值的初步研究[J]. 地震地质,33(1):114–122. doi: 10.3969/j.issn.0253-4967.2011.01.011 Liang G J,Li Z W,Zheng S F,Liu J Q,Liu D F,Qi W. 2011. Primary study on coda Q value in Longgang volcano area[J]. Seismology and Geology,33(1):114–122 (in Chinese).
刘希强,石玉燕,曲均浩,周彦文,李红,李铂,张坤. 2009. 品质因子的尾波测定方法讨论[J]. 中国地震,25(1):11–23. doi: 10.3969/j.issn.1001-4683.2009.01.002 Liu X Q,Shi Y Y,Qu J H,Zhou Y W,Li H,Li B,Zhang K. 2009. Discussion on coda wave measurement methods for quality factor[J]. Earthquake Research in China,25(1):11–23 (in Chinese).
潘建雄. 1989. 红水河龙滩水库诱发地震地质条件的探讨[J]. 地震地质,11(4):91–99. Pan J X. 1989. The geological environment for induced earthquake in Longtan reservoir of Hongshui river[J]. Seismology and Geology,11(4):91–99 (in Chinese).
史水平,于学庆,龙政强,蒙陆斌,黄树生. 2010. 龙滩库区蓄水后地震性质分析[J]. 地震地磁观测与研究,31(3):40–45. doi: 10.3969/j.issn.1003-3246.2010.03.007 Shi S P,Yu X Q,Long Z Q,Meng L B,Huang S S. 2010. Analysis on the nature of earthquake in Longtan area after the reservoir storage water[J]. Seismological and Geomagnetic Observation and Research,31(3):40–45 (in Chinese).
史水平,周斌,黄树生,陆丽娟,蒙陆斌,张帆. 2015. 2013年广西平果震群尾波衰减特征研究[J]. 中国地震,31(4):710–720. doi: 10.3969/j.issn.1001-4683.2015.04.011 Shi S P,Zhou B,Huang S S,Lu L J,Meng L B,Zhang F. 2015. Research on characteristics of seismic wave attenuation for the earthquake swarm occurring in Pingguo county,Guangxi Zhuang Autonomous Region in 2013[J]. Earthquake Research in China,31(4):710–720 (in Chinese).
阎春恒,周斌,郭培兰,孙学军,文翔. 2016. 结合波形互相关技术的龙滩水库地震双差定位[J]. 地震研究,39(3):427–435. doi: 10.3969/j.issn.1000-0666.2016.03.010 Yan C H,Zhou B,Guo P L,Sun X J,Wen X. 2016. Double-difference relocation of earthquakes in Longtan reservoir combined with waveform cross-correlation technique[J]. Journal of Seismological Research,39(3):427–435 (in Chinese).
姚宏,陈鑫,黄树生,龙政强. 2008. 龙滩水电工程数字遥测地震台网监测能力检验[J]. 地震地磁观测与研究,29(4):62–66. doi: 10.3969/j.issn.1003-3246.2008.04.012 Yao H,Chen X,Huang S S,Long Z Q. 2008. The inspection of monitoring capability of the digital telemetry seismograph network in Longtan hydropower project[J]. Seismological and Geomagnetic Observation and Research,29(4):62–66 (in Chinese).
曾健,陈天长,韩渭宾,范军,郑斯华,长谷川昭,堀内茂木,松本聪. 1997. 成都遥测台网区域尾波衰减特性及Qc值研究[J]. 地震学报,19(3):246–253. Zeng J,Chen T C,Han W B,Fan J,Zheng S H,Hasegawa A,Horiuchi S,Matsumoto S. 1997. Attenuation of coda waves and Qc value beneath the Chengdu telemetered seismic network[J]. Acta Seismologica Sinica,19(3):246–253 (in Chinese).
张锦玲,朱新运,马起杨,金春华. 2018. 应用Sato模型探究宁夏区域尾波Q值特征[J]. 地震,38(4):74–83. doi: 10.3969/j.issn.1000-3274.2018.04.007 Zhang J L,Zhu X Y,Ma Q Y,Jin C H. 2018. Characteristics of Q value in Ningxia region based on Sato model[J]. Earthquake,38(4):74–83 (in Chinese).
张天中,高龙生,张卫平. 1990. 滇西试验场区的Q值及其随时间窗的变化[J]. 地震学报,12(1):12–21. Zhang T Z,Gao L S,Zhang W P. 1990. The Q-factor and its time window dependence in the area of western-Yunnan seismological experiment site[J]. Acta Seismologica Sinica,12(1):12–21 (in Chinese).
周斌,孙峰,阎春恒,薛世峰,史水平. 2014. 龙滩水库诱发地震三维孔隙弹性有限元数值模拟[J]. 地球物理学报,57(9):2846–2868. doi: 10.6038/cjg20140911 Zhou B,Sun F,Yan C H,Xue S F,Shi S P. 2014. 3D-poreelastic finite element numerical simulation of Longtan reservoir-induced seismicity[J]. Chinese Journal of Geophysics,57(9):2846–2868 (in Chinese).
周连庆,赵翠萍,修济刚,陈章立,郑斯华. 2008. 利用天然地震研究地壳Q值的方法和进展[J]. 国际地震动态,(2):1–11. doi: 10.3969/j.issn.0253-4975.2008.02.001 Zhou L Q,Zhao C P,Xiu J G,Chen Z L,Zheng S H. 2008. Methods and developments of research on crustal Q value by using earthquakes[J]. Recent Developments in World Seismology,(2):1–11 (in Chinese).
朱新运. 2011. 衰减、场地响应等地震波传播相关信息综合研究[D]. 北京: 中国地震局地球物理研究所: 66−92. Zhu X Y. 2011. Comprehensive Study on Seismic Wave Propagation Related Information Such as Attenuation and Site Response[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 66−92 (in Chinese).
朱新运,杨钢宇,张震峰. 2005. 基于Sato模型的近震S波尾波Q值求解及分析软件研制[J]. 地震地磁观测与研究,26(3):63–70. doi: 10.3969/j.issn.1003-3246.2005.03.010 Zhu X Y,Yang G Y,Zhang Z F. 2005. Developing on the data processing software of Q value of local equathquake[J]. Seismological and Geomagnetic Observation and Research,26(3):63–70 (in Chinese).
朱新运,刘杰,张帆. 2006. 基于Aki模型的近震S波尾波Q值求解及分析软件研制[J]. 地震研究,29(1):76–80. doi: 10.3969/j.issn.1000-0666.2006.01.015 Zhu X Y,Liu J,Zhang F. 2006. Development of Q-value calculating and processing software using S-wave coda in local earthquakes based on Aki model[J]. Journal of Seismological Research,29(1):76–80 (in Chinese).
朱新运,于俊谊,张帆. 2013. 浙江珊溪水库库区地震波衰减特征研究[J]. 地震学报,35(2):199–208. doi: 10.3969/j.issn.0253-3782.2013.02.006 Zhu X Y,Yu J Y,Zhang F. 2013. Characteristics of seismic wave attenuation in Shanxi reservoir area of Zhejiang Province[J]. Acta Seismologica Sinica,35(2):199–208 (in Chinese).
Aki K. 1969. Analysis of the seismic coda of local earthquakes as scattered waves[J]. J Geophys Res,74(2):615–631. doi: 10.1029/JB074i002p00615
Aki K. 1980. Attenuation of shear-waves in the lithosphere for frequencies from 0.05 to 25 Hz[J]. Phys Earth Planet Inter,21(1):50–60. doi: 10.1016/0031-9201(80)90019-9
Aki K,Chouet B. 1975. Origin of coda waves:Source,attenuation,and scattering effects[J]. J Geophys Res,80(23):3322–3342. doi: 10.1029/JB080i023p03322
Canas J A,Mitchell B J. 1978. Lateral variation of surface-wave anelastic attenuation across the Pacific[J]. Bull Seismol Soc Am,68(6):1637–1650.
Carpenter P J,Sanford A R. 1985. Apparent Q for upper crustal rocks of the central Rio grande rift[J]. J Geophys Res,90(B10):8661–8674. doi: 10.1029/JB090iB10p08661
Castro R R,Rebollar C J,Inzunza L,Orozco L,Sánchez J,Gálvez O,Farfán F J,Méndez I. 1997. Direct body-wave Q estimates in northern Baja California,Mexico[J]. Phys Earth Planet Inter,103(1/2):33–38.
Drouet S,Souriau A,Cotton F. 2005. Attenuation,seismic moments,and site effects for weak-motion events:Application to the Pyrenees[J]. Bull Seismol Soc Am,95(5):1731–1748. doi: 10.1785/0120040105
Horasan G,Boztepe-Güney A. 2004. S-wave attenuation in the Sea of Marmara,Turkey[J]. Phys Earth Planet Inter,142(3/4):215–224.
Hwang H J,Mitchell B J. 1987. Shear velocities,Q β,and the frequency dependence of Q β in stable and tectonically active regions from surface wave observations[J]. Geophys J Int,90(3):575–613. doi: 10.1111/j.1365-246X.1987.tb00743.x
Jacobson R S,Jr Shor G G,Bée M. 1984. A comparison of velocity and attenuation between the Nicobar and Bengal Deep Sea Fans[J]. J Geophys Res,89(B7):6181–6196. doi: 10.1029/JB089iB07p06181
Lee W B,Solomon S C. 1979. Simultaneous inversion of surface-wave phase velocity and attenuation:Rayleigh and Love waves over continental and oceanic paths[J]. Bull Seismol Soc Am,69(1):65–95.
Nava F A,García-Arthur R,Castro R R,Suárez C,Márquez B,Núñez-Cornú F,Saavedra G,Toscano R. 1999. S wave attenuation in the coastal region of Jalisco-Colima,México[J]. Phys Earth Planet Inter,115(3/4):247–257.
Pulli J J. 1984. Attenuation of coda waves in New England[J]. Bull Seismol Soc Am,74(4):1149–1166.
Sanders C O,Nixon L D. 1995. S wave attenuation structure in Long Valley Caldera,California,from three-component S-to-P amplitude ratio data[J]. J Geophy Res,100(B7):12395–12404. doi: 10.1029/95JB00962
Sato H. 1977. Energy propagation including scattering effects single isotropic scattering approximation[J]. J Phys Earth,25(1):27–41. doi: 10.4294/jpe1952.25.27
Wong V,Rebollar C J,Munguía L. 2001. Attenuation of coda waves at the Tres Virgenes volcanic area,Baja California Sur,Mexico[J]. Bull Seismol Soc Am,91(4):683–693. doi: 10.1785/0120000025
Zucca J J,Hutchings L J,Kasameyer P W. 1994. Seismic velocity and attenuation structure of the Geysers geothermal field,California[J]. Geothermics,23(2):111–126. doi: 10.1016/0375-6505(94)90033-7
-
期刊类型引用(1)
1. 都燊,俞言祥,肖亮. 利用密集台阵资料研究宾川盆地浅层介质尾波Q值. 地震学报. 2023(06): 1011-1024 . 本站查看
其他类型引用(2)