文章快速检索    
  地震地磁观测与研究  2023, Vol. 44 Issue (4): 68-76  DOI: 10.3969/j.issn.1003-3246.2023.04.011
0

引用本文  

单晓琳, 徐长银, 李姗姗, 等. 库尔勒地震台痕量氢气浓度观测数据小波分析[J]. 地震地磁观测与研究, 2023, 44(4): 68-76. DOI: 10.3969/j.issn.1003-3246.2023.04.011.
SHAN Xiaolin, XU Changyin, LI Shanshan, et al. Wavelet analysis of trace hydrogen concentration data from Korla Seismic Station[J]. Seismological and Geomagnetic Observation and Research, 2023, 44(4): 68-76. DOI: 10.3969/j.issn.1003-3246.2023.04.011.

作者简介

单晓琳(1989-), 女, 主要从事地球物理观测工作。E-mail: 1090558023@qq.com

文章历史

本文收到日期:2023-02-27
库尔勒地震台痕量氢气浓度观测数据小波分析
单晓琳   徐长银   李姗姗   张剑锋     
中国新疆维吾尔自治区 841000 库尔勒地震监测中心站
摘要:以2015—2022年库尔勒地震台ATG-6118H型痕量氢分析仪观测数据为研究对象,对新疆2016年12月8日呼图壁MS 6.2、2017年8月9日精河MS 6.6两次地震利用db4小波进行多尺度小波分析,以探究震前痕量氢气浓度的高频短临前兆异常特征。此外,对新疆2018年10月16日精河MS 5.4、2019年11月12日库车MS 5.6、2020年1月16日库车MS 5.7地震进行震例验证。结果表明,距库尔勒痕量氢气观测点400 km范围内,在强震发生前3—5个月痕量氢气浓度观测数据高频细节分析出现异常,异常特征表现为痕量氢气浓度变化大于2倍均方差。
关键词小波分析    痕量氢    资料分析    
Wavelet analysis of trace hydrogen concentration data from Korla Seismic Station
SHAN Xiaolin   XU Changyin   LI Shanshan   ZHANG Jianfeng     
Korla Seismic Station, Xinjiang Uygur Autonomous Region 841000, China
Abstract: Based the observation data of ATG-6118H trace hydrogen analyzer at Korla Seismic Station from 2015 to 2022, db4 wavelets were used to analyze the Hutubi MS 6.2 earthquake on December 8, 2016 and Jinghe MS 6.6 earthquake on August 9, 2017 in Xinjiang to investigate the characteristics of short precursory anomalies of trace hydrogen gas concentration before the earthquakes. In addition, seismic validations were carried out for the earthquakes on October 16, 2018 in Jinghe (MS 5.4), November 12, 2019 in Kuqa (MS 5.6), and January 16, 2020 in Kuqa (MS 5.7). The results show that within a 400 km range of the Korla trace hydrogen observation point, anomalies appeared in the high-frequency detail analysis of trace hydrogen gas concentration observation data 3-5 months before the occurrence of strong earthquakes. These anomalies manifested as the change in trace hydrogen gas concentration greater than twice the standard deviation.
Key words: wavelet analysis    trace hydrogen    data analysis    
0 引言

目前,国内用于地震前兆资料数据处理的软件和方法较多,但前兆数据的干扰因素也较多,数据处理中,对于因外界环境影响、人为干扰、不明原因的脉冲干扰、地震同震记录以及仪器零漂等引起的趋势性变化而造成的数据波动,软件没有相应处理功能,并且,按照规定可不予处理。这些干扰信息却能影响数据分析结果,使得震前的一些短临趋势性变化不易识别。由于小波分析对数字化地震前兆资料干扰的识别与消除,以及对不同频率信息的识别功能较强(冉启文等,2002)。近年来,小波变换理论在地震资料分析中得到广泛应用,宋治平等(2003),将小波理论用于分析地震前兆观测资料,并取得一些研究结果(耿杰等,2009),徐长银等,2017)。为加强新疆天山中段监视区的震情跟踪工作,2013年11月经现场实地勘选,在库尔勒市霍拉山观测山洞架设了杭州超距科技有限公司研制的ATG-6118H型痕量氢在线自动分析仪,开展断层氢气定点连续观测。本文在前人研究的基础上,对2015—2022年库尔勒地震台痕量氢观测资料进行研究,采用Mapsis软件,运用小波分析对痕量氢浓度进行多尺度分析,重点分析距观测点400 km范围内2次中强地震前痕量氢浓度观测数据中的短期异常。

1 断层气痕量氢观测概况 1.1 观测区地质构造背景概况

库尔勒痕量氢观测点位于南天山地区的北轮台断裂附近(图 1),该区地处天山构造带和塔里木盆地的交汇部位,天山山体与两侧盆地之间存在由山体向盆地和坳陷逆冲的逆断裂(邓起东等, 2001, 杨晓平等,2006),北轮台断裂(亦被称为霍拉山山前断裂)即为天山与两侧盆地的分割断裂,同时也是库车逆断裂-褶皱带的根部控制断裂,其向南的逆冲缩短速率为(4.7±2.5 mm /a)(杨少敏等,2008),第四纪以来该条断裂带活动性明显。由于天山与塔里木盆地剧烈的挤压碰撞作用,导致该区在地质构造(邓起东等,2001)、地震活动(曲延军等, 2010, 王筱荣等,2010)、地壳介质结构特征(Kumar et al,2005Li et al,2023刘玉虎等, 2011, 刘洁等, 2007, 李金等,2021)等方面均显示出较为独特的特征,是天山地区强震前兆观测的绝佳位置。库尔勒地震台痕量氢观测点位于库尔勒市北霍拉山,在库尔勒地震台形变前兆综合观测“U”形山洞内,仪器架设在距东洞口约150 m的跨断层观测室内,之间连接4道船舱门。

图 1 库尔勒地震台痕量氢浓度观测点与地震震中分布 Fig.1 Trace hydrogen concentration and epicenter distribution at Korla Seismic Station

库尔勒地震台断层氢气观测孔由人工开挖,孔深1.7 m,裸孔直径1.0 m,开挖至基岩。观测孔内放置1根直径250 mm、长135 cm的PVC管,管底50 cm处打有气孔(孔径10 mm),连接集气花管。PVC管上部通过4个变径管头(长35 cm)相连,连接处用PVC胶粘接,直径最终为25 mm,而后用橡皮管塞密封。将氢气观测软管插入橡皮管塞中,为利于气体通过,在集气管周围50 cm铺设一层透气性好的砾石层,再将沙石回填,孔表层用水泥沙浆密封(向阳等,2018)(图 2)。

图 2 断层气集气装置剖面示意图 Fig.2 Schematic diagram of fault gas collecting device
1.2 氢气浓度影响因素分析

氢气是最轻的气体,在常温常压下,无色、透明、无味、难溶于水、密度最小、极易燃烧。震例总结显示,中强地震前后水井中溶解氢气浓度均快速增大,主要异常特征为浓度高值异常,且异常持续时间较短,变化幅度较大(范雪芳等, 2012, 向阳等,2018)。但氢气浓度异常变化并非都是地震前兆异常,也受到一些干扰因素的影响,气象因素是主要干扰因素之一(常秋君等,1993)。如气温变化可以引起近地表岩石热弹性应变周期变化,进而使岩石裂隙中的气体浓度发生周期性变化(闫玮等,2019);气压变化可以改变岩石孔隙中的压力梯度及断层中气体的迁移速率,使得断层氢气浓度发生变化(李杰等,2004闫玮等,2019)。

自2013年11月库尔勒地震台开始痕量氢观测以来,氢气浓度数据连续且稳定,呈现较规律的年变趋势,即夏低冬高的变化趋势。选取2016年8、12月痕量氢浓度原始数据及气压、温度进行分析(图 3)。

图 3 库尔勒地震台断层氢气浓度、气压、温度整点值 (a)2016年8月;(b)2016年12月 Fig.3 Fault hydrogen concentration, pressure and temperature at Korla Seismic Station

图 3可见,库尔勒地震台断层氢气浓度无明显日变化规律,曲线存在单峰、双峰等多种变化形态。相比而言,氢气浓度日变形态与气压日变形态间具有一定同步性,但与温度日变形态不一致。这表明气压的日变形态对库尔勒地震台断层氢气浓度存在一定影响,温度的日变形态则影响不大。

虽然库尔勒地震台断层氢气浓度日变形态受气压影响较大,受温度影响较小,但氢气浓度与气压、温度的年变特征间均存在一定关系。图 4为2019—2022年库尔勒地震台断层氢气浓度预处理数据。由图 4可见,每年1—3月、11—12月氢气浓度测值相对较高且变化平稳,浓度约为1.2×10-6;4—8月氢气浓度下降,且下降幅度较大,测值约为0.6×10-6,其中,6—8月测值降到最低,变化范围0.1×10-6—0.2×10-6

图 4 库尔勒地震台断层氢气浓度年变化 Fig.4 Annual variation of hydrogen concentration in faults at the Korla Seismic Station
1.3 氢气浓度与气压、温度间的相关性

断层氢气浓度受多种因素的影响,如断层开合、采样深度、采样孔周围地温、湿度、气象三要素、观测环境等(闫玮等,2022)。由于条件所限,库尔勒地震台断层氢气观测点没有独立的气象要素监测仪器,因此使用氢气观测仪器显示的温度,该数据是仪器所在山洞环境的温度。测点位于地下山洞观测室内,受干扰影响较小,多年温度恒定,为16.0—16.5 ℃,气压为887—918 hPa,氢气浓度、温度、气压均呈现较好的季节变化特征。因测点位于亚欧大陆内部,其温带大陆性气候特征较显著,干旱少雨,且蒸发量较大,降雨影响较小,故忽略不计。因此,重点分析库尔勒地震台断层氢气浓度与气压、气温间的关系。

2015年1月1日至2022年12月31日库尔勒地震台痕量氢浓度日均值、温度、气压的变化如图 5(a)5(b)6(a)6(b)所示。利用一元线性回归分析分别得到痕量氢浓度与温度、气压间的相关系数,与气压间的相关系数为0.642 8,呈正相关[图 5(c)]。

图 5 库尔勒断层氢气浓度(a)、气压(b)及二者间的相关性(c) Fig.5 Hydrogen concentration (a), pressure (b) and their correlation (c) at the Korla Fault
图 6 库尔勒断层氢气浓度(a)、温度(b)及二者间的相关性(c) Fig.6 Hydrogen concentration (a), temperature (B) and their correlation (c) along the Korla Fault
2 小波分析理论 2.1 小波变换理论

小波变换是一种对于信号的时间—尺度(时间—频率)分析方法,适合分析非平稳信号和提取信号的局部特征,对不同频率的信息识别功能较强,对趋势性变化和短临异常能有效进行分离。其主要特点是通过变换能够充分突出某些方面的特征,可以较好地抑制作用于地壳表面和观测系统的多种干扰因素,有助于提取异常信号,是地震前兆数据处理的一种有效方法(郑治真等,2001万永革等,2003)。对于数字信号可近似表示为

$ f(x) \cong A_j f(x)=a_k^j f(x)+d_k^j f(x) $ (1)

式中,$a_k^j f(x), d_k^j f(x)$分别为信号f(x)在分辨率为j时的近似部分、细节部分;正交展开系数$a_k^j, d_k^j$分别为离散近似部分、离散细节部分。利用小波分析方法,可以从近似部分确定趋势变化信息,从细节部分识别短临异常(万永革等,2003)。

2.2 小波基函数选取

小波基函数即小波分析的基,不是唯一存在的,小波基函数的选取没有明确标准,应用小波变换对同一信号进行分析时,所选用的小波基函数不同,可使分析结果差异较大。在选择合适的小波基进行信号处理时需考虑正则性、消失矩及处理信号与小波基间的相似性等因素。一般情况下,应选择与异常波形相似且震荡较平缓、连续、光滑的小波基函数(李杰等,2004)。Daubechies小波是离散正交小波,可表示为dbN,其中,N为小波的阶。dbN不具对称性(即非线性相位),没有明确的解析表达式(郑治真等,2001)。运用小波分析法对痕量氢气浓度进行多尺度分析时,通过尺度因子的改变,将信号分解为低频、高频信息,能够将数据的趋势性变化与细节变化分开,可有效地抑制噪声等干扰因素,提取震前相关的前兆异常信息(刘水莲等,2010)。研究发现,利用db4小波提取断层氢气浓度异常数据效果较好。因此,采用db4小波对库尔勒地震台痕量氢气浓度资料进行多尺度小波分析,主要以2016年12月8日新疆呼图壁MS 6.2、2017年8月9日新疆精河MS 6.6两次强震为例,地震发震前短临异常信息特征。

3 痕量氢气浓度资料分析处理 3.1 小波变换

采用Mapsis软件,选取2016年1月1日至2017年12月31日库尔勒地震台痕量氢气浓度日均值原始数据进行预处理和小波分析,以2016年12月8日呼图壁MS 6.2、2017年8月9日精河MS 6.6两次地震为例,选取db4小波对其分解,分解尺度为5,对低频、高频信号进行分离,分别得到趋势部分(低频)和细节部分(高频)的图像(图 78)。分解后低频信息主要反映趋势性变化过程,细节部分主要反映高频突跳信息(江崇昆等, 2013, 杨绍富等,2011)。由图 78可见,随着分解层数的增加,数字信号中的干扰等因素逐渐被分离,数据曲线反映的局部细节信息逐渐减少,曲线也随之越来越清晰、光滑。

图 7 2016—2017年痕量氢气浓度数据多尺度小波分析低频曲线 (a)原始图像;(b)、(c)、(d)、(e)、(f)经尺度1、2、3、4、5分解后趋势部分的图像 Fig.7 Multi-scale wavelet analysis of low frequency curve for trace hydrogen concentration data from 2016 to 2017
图 8 2016—2017年痕量氢气浓度数据多尺度小波分析高频曲线 (a)原始图像;(b)、(c)、(d)、(e)、(f)经尺度1、2、3、4、5分解后细节部分的图像 Fig.8 Multi-scale wavelet analysis of high frequency curve for trace hydrogen concentration data from 2016 to 2017

图 78可见,随着分解层数的增加,低频部分噪声、干扰等因素逐渐消失,曲线的年变特征逐渐明显,曲线亦逐渐光滑。但高频信号被剥离的同时,也剥离了构造作用下引起的地震短临前兆信息,因此,以下重点探讨小波变换对库尔勒痕量氢观测资料高频信号的短临异常,从细节部分识别的信号幅值变化识别地震短临异常特征。

3.2 高频短期异常变化分析

图 8(b)(f)可见,随着分离层数的增加,分解出的细节信号周期性越明显,序列变化越不平稳。当尺度为5时,细节信号部分的周期相位和幅值发生明显变化,这可能与该频率范围内正常信号叠加了相同频率的其他信号有关,而叠加的信号可能包含在构造应力作用下,断层氢气在活动断裂带的逸出速率发生变化的地震前兆异常。因此,选取尺度为5的曲线进行高频短期异常变化分析(图 9)。由图 9可见,从2016年9月和2017年5月开始,曲线分别出现了显著的超正常幅度的变化,主要表现为信号幅值、年周期相位、年周期幅值均超线。2016年12月8日呼图壁MS 6.2、2017年8月9日精河MS 6.6两次强震前,痕量氢气浓度变化超出了2倍均方差。笔者认为这2次明显的超线可能是2次6级地震前兆短期异常,并且地震一般发生在高频异常后的3—5个月。

图 9 2016—2017年尺度为5的小波分析高频部分曲线 (a)原始图像;(b)尺度为5的细节信号 Fig.9 Wavelet analysis of the high frequency part curve with scale 5 in 2016-2017
3.3 震例检验

选取震中距库尔勒霍拉山观测点380 km的2018年10月16日精河MS 5.4,距观测点260 km的2019年11月12日库车MS 5.6、2020年1月16日库车MS 5.7地震,选用2018年1月1日至2020年12月31日痕量氢气浓度日均值进行异常分析,利用db4小波对其分解,分解尺度为5(图 10)。由图 10可见,2018年4月中旬、2019年8月中旬、2019年10月中旬细节信号的部分周期相位、幅值发生显著变化,痕量氢气浓度变化超出2倍均方差。笔者认为是地震前兆,并且异常出现后3—5个月相继发生精河、库车地震。研究认为,对于震中距观测点400 km范围内的中强地震,依据不同的震级、震中距,痕量氢气浓度在强震发生前3—5个月小波高频细节分析出现异常,异常特征表现为痕量氢气浓度变化超出了2倍均方差。

图 10 2018—2020年尺度为5的小波分析高频部分曲线 (a)原始图像;(b)尺度为5的细节信号 Fig.10 The wavelet analysis of the high frequency partial curve with scale 5 in 2018-2020
4 结论与讨论

(1)应用小波变换分析了库尔勒地震台2015—2022年痕量氢气浓度观测资料,利用db4小波对痕量氢气浓度进行多尺度分离,对比分析低频、高频信息,发现痕量氢气浓度高频短期异常变化特征明显。认为小波分析可以作为识别与消除地震前兆资料中的干扰因素、提取趋势性异常与短期异常的有力手段

(2)结合所选震例,在距观测点400 km范围内,对痕量氢气浓度高频细节的分析显示较明显的异常特征,地震一般在异常出现前3—5个月左右发生。

(3)库尔勒地震台痕量氢气浓度年变特征较明显,影响因素主要是气象因素(气压、温度),主要特征是夏低冬高的变化趋势。

(4)库尔勒地震台痕量氢气浓度与气压间的相关系数为0.642 8,呈正相关,既随着气压的增大,痕量氢气浓度会随之逐渐上升;痕量氢气浓度与温度间的相关系数为–0.463 9,呈负相关,既随着温度的上升,痕量氢气浓度随之逐渐下降。

(5)库尔勒地震台断层氢气浓度变化表现为夏低冬高。笔者认为造成这种现象的原因可能是观测点周围土壤表层局部潮湿所致,即当气温在0℃以下时,潮湿土壤中的水分冻结,形成冻土层,当冻土层达到一定厚度时可形成良好的封闭系统。

参考文献
常秋君, 何跟巧, 郭玉英, 等. 断层气定点观测主要影响因素的研究[J]. 西北地震学报, 1993, 15(4): 82-89.
邓起东, 冯先岳, 张培震, 等. 天山活动构造[M]. 北京: 地震出版社, 2000: 42-60.
范雪芳, 黄春玲, 刘国俊, 等. 山西夏县痕量氢观测资料的初步分析[J]. 山西地震, 2012(3): 7-12.
耿杰, 陈安方, 潘双进. 用小波变换方法提取地下流体观测异常信息[J]. 地震研究, 2009, 32(1): 12-17.
江崇昆, 杨绍富, 徐长银, 等. 利用小波方法分析新43泉气氡异常特征[J]. 内陆地震, 2013, 27(3): 257-262.
李杰, 李明海, 刘希强, 等. 利用小波变换方法分析跨断层形变异常[J]. 大地测量与地球动力学, 2004, 24(3): 73-78.
李金, 高原, 王琼. 天山构造带上地壳介质各向异性分区特征[J]. 中国科学: 地球科学, 2021, 51(4): 582-597.
刘洁, 刘启元, 郭飚, 等. 中国境内天山上地幔小尺度对流与造山作用[J]. 中国科学D辑: 地球科学, 2007, 37(6): 728-735.
刘水莲, 刘礼诚, 黄跃进, 等. 小波分析在永安地震台前兆观测资料处理中的应用[J]. 华北地震科学, 2010, 28(3): 53-57.
刘玉虎, 刘兴旺, 郑建京, 等. 天山南北地块构造演化与地幔对流耦合动力机制[J]. 地球物理学进展, 2011, 26(5): 1544-1556.
曲延军, 王琼, 聂晓红. 新疆南部构造区带与地震活动状态研究[J]. 内陆地震, 2010, 24(4): 298-306.
冉启文, 谭立英. 小波分析与分数傅里叶变换及应用[M]. 北京: 国防工业出版社, 2002: 1-124.
宋治平, 武安绪, 王梅, 等. 小波分析方法在形变数字化资料处理中的应用[J]. 大地测量与地球动力学, 2003, 23(4): 21-270.
万永革, 齐福荣, 孟晓春, 等. 中国大陆及华北地区地震资料的小波分析[J]. 大地测量与地球动力学, 2003, 23(4): 28-33.
王筱荣, 王琼. 南天山地震带地震平静与中强震关系再研究[J]. 西北地震学报, 2010, 32(3): 279-285.
向阳, 孙小龙, 高小其, 等. 新疆库尔勒断层氢气浓度的影响因素及其地震预测的潜在效能评价[J]. 中国地震, 2018, 34(1): 48-59.
徐长银, 杨绍富, 孙海军, 等. 2016年12月8日呼图壁MS 6.2地震库尔勒体应变异常分析[J]. 内陆地震, 2017, 31(3): 284-289.
闫玮, 李晓东, 阿里木江·麦麦提依明, 等. 新疆伽师55井水温梯度实验及结果应用[J]. 内陆地震, 2019, 33(4): 337-344.
闫玮, 朱成英, 李新勇, 等. 新疆伽师MS 6.4地震前后柯坪断裂土壤气地球化学特征变化[J]. 内陆地震, 2022, 36(4): 362-368.
杨绍富, 徐长银. 小波分析在库尔勒地震台气氡资料处理中的应用[J]. 内陆地震, 2011, 25(4): 329-333.
杨少敏, 李杰, 王琪. GPS研究天山现今变形与断层活动[J]. 中国科学D辑: 地球科学, 2008, 38(7): 872-880.
杨晓平, 冉勇康, 宋方敏, 等. 西南天山柯坪逆冲推覆构造带的地壳缩短分析[J]. 地震地质, 2006, 28(2): 194-204.
郑治真, 沈萍, 杨选辉, 等. 小波变换及其MATLAB工具的应用[M]. 北京: 地震出版社, 2001: 294-301.
Kumar P, Yuan X, Kind R, et al. The lithosphere-asthenosphere boundary in the Tien Shan-Karakoram region from S receiver functions: evidence for continental subduction [J]. Geophysical Research Letters, 2005, 32, L07305, doi: 10.1029/2004GL022291, 1-4.
Li J, Jiang H K, Wang Q. Crustal attenuation structure of the Tianshan tectonic belt and its spatiotemporal variations[J]. Front. Earth Sci, 2023, 11: 1094151. DOI:10.3389/feart.2023.1094151