基金项目:[JP2]中国地震局地震科技星火计划(XH15006)、山西省自然基金(201901D111462)和中国地震局震情跟踪定向工作任务(2016020302)联合资助.
通讯作者:李宏伟(1986-),工程师,主要从事地壳形变分析预报工作.E-mail:347719545@qq.com
(1.山西省地震局,山西 太原 030021; 2.中国辐射防护研究院,山西 太原 030006)
(1. Shanxi Earthquake Agency,Taiyuan 030021,Shanxi,China)(2. China Institute for Radiation,Taiyuan 030006,Shanxi,China)
hydrogen in the soil; air temperature; air pressure; dynamic characteristic; influence factors
收集全国15个不同区域、不同地质环境、不同深度测点的H2浓度,以及气温、气压整点值连续观测资料,研究高精度土壤H2浓度日、月动态变化特征,通过计算H2浓度与气温、气压的相关系数,定量分析不同观测孔深气温、气压对H2浓度的影响程度。结果 表明:土壤H2浓度动态变化与气象因素密切相关,观测孔优选深度为4~11 m。
Continuous observation data of hydrogen levels and hourly values of air temperature and air pressure for 15 different regions,geological environments,and depths in China are collected. Daily dynamic characteristic of levels of high-accuracy hydrogen in the soil is studied and influence degree of air temperature and air pressure of different hole depths to hydrogen levels is analyzed quantitatively through the calculation of correlation coefficient of hydrogen levels and different air temperature and atmospheric pressure. The result shows that observation dynamics of levels of hydrogen in the soil and atmosphere factor is closely related and the best depth of observation hole is 4~11 m.
利用地表自由逸出气体以及溶解于水、吸附于土壤中的气体浓度变化来监测预报地震是当前国内外广泛采用的地震预报方法之一(Whiticar et al,1999; Hernández et al,2000)。断层带是深部气体逸出的良好通道,气体在应力或热动力作用下,能沿上下贯通的活断层通道自深部上涌(王永才等,1992; Wakita et al,1980)。粟启初(1992)通过实验指出H2逸出对断层的运动、地震的发生有直接的诱发作用; 邵济安等(2010)研究认为,H2主要来源于地幔,少量地壳来源是由于活动断层产生的局部高温高压趋使中-浅层地热储库中的CH4转化为H2。
在20世纪,有学者对断层气进行了一些研究,但受限于当时的观测仪器精度不够,该研究断断续续。随着国内高精度氢连续观测仪器的研制,断层气研究进一步深入。分析地震监测资料表明,H2浓度是对地震前兆响应的灵敏气体组分,特别是在短临阶段的映震能力方面具有明显的反映(车用太等,2015; 刘耀炜等,2006; 张炜等,1988; 耿杰等,2019)。地震前出现的H2浓度异常特征显著,通常表现为震前几天至几十天的突发性大幅上升。由于H2在土壤及地下水中的背景值较低,因此,在断层带H2浓度观测中,其出现异常的幅度较高,通常表现为数量级的增幅,且H2浓度会随着地震释放能量的增大而增加。车用太等(2015)认为断层带H2浓度观测有可能提高短临预测水平。大量的观测资料表明,断层土壤气观测的影响因素较为单一,主要是气象因素,即气温、气压等,年变动态规律清楚,由此可知,土壤H2浓度的异常变化突出,有较高的映震效能。
断层带土壤H2浓度观测的首次实验研究是在山西开展,后逐步向全国推广。目前,利用痕量氢自动分析仪开展断层带土壤H2浓度观测时间不长,尚处于实践探索阶段,对观测条件和各种干扰因素等研究尚未开始。对于一个新开展的观测项目,必须研究它的动态特征及其影响因素,其研究结果对异常变化的识别,甚至是异常机理解释和预测效能评价都有重要作用。为此,本文收集了全国目前处于运行阶段的痕量H2浓度观测资料,同时,同步观测环境温度(气温)和气压,采用相关分析、回归分析的定量分析方法,研究了气温、气压及观测孔深度对断层带土壤H2浓度的影响。
关于H2浓度观测及其与地震活动的关系研究,已有近40年的历史。国际上普遍认为断层带逸出气中的H2是效果较好的地震前兆观测项目之一,H2的浓度与地震活动有密切关系,异常幅度大,震前异常时间短,短临特征显著,是有可能量化的。车用太等(2015)的研究结果表明,震中距、震级、震源深度不同,H2浓度的异常都有显著差别。从已有的震例来看,在震前10~40 d,H2浓度变化最大。
断层带土壤H2浓度观测点或观测井(泉)的布设地点主要是选择已被确定为未来10年内有可能发生MS≥6.0地震的地区,或近几年来MS≥5.0地震活动频繁的地区,其次还需要参考活动断块与地震活动关系。总体上,断层带土壤H2浓度观测点的布设要求如下:①在断层破碎带上,隐伏断裂带区测点选在断裂带的正上方; ②断层的上盘; ③岩浆岩发育区; ④作为试验观测研究,测点可考虑选在不同性质的断层带和断裂的不同部位,进行对比观测研究; ⑤在有条件的地区,也可打穿断裂带的自流井进行地下水中逸出气H2浓度观测。
国内一般采用气相色谱仪和质谱仪进行H2浓度观测,检出限一般为1×10-6。如纪振红等(2010)系统总结了不同纯度H2浓度的分析方法,主要有气相色谱法和热导池气相色谱法等。但一般水中溶解或土壤中H2浓度的背景值为0.5×10-6,甚至更低,用气相色谱仪等很难检测到H2浓度的背景动态。而杭州电子科技大学和杭州超钜科技有限公司研制的ATG-6118H痕量氢自动分析仪的灵敏度高,检出限为0.05×10-6,稳定性好,可准确连续记录到H2浓度的变化; 该仪器还具有网络实时监控和温度、气压辅助观测功能,这将有效解决采用气相色谱仪等进行观测时仪器灵敏度较低等问题。范雪芳等(2015,2016)对ATG-6118H痕量氢自动分析仪的稳定性及其记录的连续观测数据进行了评估,评估结果表明,分析仪产出的测值稳定,2台仪器平行观测动态一致性很好,能够实现断层带土壤H2浓度的连续观测。因此本文使用ATG-6118H痕量氢自动分析仪进行研究。
2010年山西省首先开始痕量H2浓度观测,目前在全国已布设了30多个测点,数据类型以整点值为主,其次有10 min,20 min和30min值。目前的观测多为断层带土壤气观测,部分为热水井或温泉水逸出气观测,本文仅对断层带土壤H2浓度进行研究。
通过筛选和资料分析处理,剔除无效数据,本文共筛选了15个测点、不同时段连续稳定的H2浓度观测数据,对缺数情况通过内插法进行处理。
为了研究断层带土壤H2浓度的正常动态变化特征,需要按照动态曲线所反映的时间长度,将研究区间划分为多年动态、年动态、月动态、日动态等。这类动态所反映的是在一定时段内某测项的周期变化特征。
常秋君等(1993)对断层气(Rn,CO2)的观测资料进行研究,发现H2测值具有明显的周期性变化,依据H2变化周期的尺度长短,可将其分为若干类型。由于目前H2浓度观测时间比较短,故仅提取H2浓度日、月动态变化进行分析研究(图1,2)。
由图1可看出,H2浓度日动态变化特征清晰,可分为3种类型,即双峰双谷型(类似于正弦型图1a、余弦型图1b)、单峰型(图1c,1d)和高低型(前半天高值、后半天低值的图1e或前半天低值、后半天高值的图1f)。日变幅度为(0.2~2)×10-6,不同测点日变幅度不同。
由于H2浓度的日动态变化可能与气温、气压、固体潮效应有关,但有些观测点的日动态变化,可能是受地温的影响,这需要分析气温与H2浓度的相位是否一致来判别。如果H2浓度日变相位滞后于气温若干小时,根据地温滞后于气温这一判据,就可以推断出H2浓度的变化与地温相关,所以今后的断层气H2观测点应该同步观测集气层位的地温。
使用整点值提取连续一个月的数据资料进行分析,发现H2浓度月动态变化以平稳型为主(图2)。H2浓度日变规律清晰,因此月动态变化表现出高频变化特征,其月变幅度和日变幅度接近。
图1 痕量H2浓度日动态变化类型
Fig.1 Daily dynamic variation type of trace hydrogen
图2 痕量H2浓度月动态变化曲线
Fig.2 Normal monthly dynamic variation curve of trace hydrogen
土壤H2观测点多位于地表浅层(1.25~16 m),观测结果受气温、气压的影响,随着季节变化,H2观测值也可能发生变化; 观测孔深度不同,其受气温、气压影响程度也不同。要使新的测项在短临跟踪和分析预报中发挥作用,必须对该测项的正常动态特征及其影响因素进行研究分析。本文选取了各种类型、不同时段、不同观测孔深观测数据,定量分析气温和气压对H2浓度影响。
计算H2浓度与气温之间的相关系数,结果如表1、图3所示。参与计算的样本数最少为538个,最多为26 268个。从15个测点不同时段、不同数据类型的处理结果可见,在河北张北海流图的H2浓度与气温的相关关系最大,r=0.945 7; [JP2]新疆阿克苏相关系数r=0.656,也比较显著; 在15个测点中,有8个测点呈负相关,占53.3%; 7个测点呈正相关,占46.7%; 相关系数≤-0.200或≥0.200的测点有11个。
在某一显著水平前提下,若相关能够成立,所要求的相关系数大小与样本数有关。当显著水平为0.01,样本数为200时,若相关系数为0.181,则表明相关显著(林记曾,1981)。因此认为气温对H2浓度的影响是明显的。
表1 H2浓度与气温、气压相关关系统计结果
Tab.1 The statistics of correlation results between hydrogen levels air temperature and air pressure
图3 不同观测点H2浓度与气温相关关系
Fig.3 The correlation of hydrogen levels and air temperature at different observation points
采取一元线性回归方法,计算H2浓度与气压的相关关系,计算结果见表1和图4。从表1可见,张北海流图土壤H2浓度与气压之间的相关关系最大,r=0.557 9; 在15个测点中,有11个测点呈负相关,占73.3%; 相关系数≤-0.200或≥0.200的测点有9个。这表明气压对H2浓度的影响也是比较明显的,但其影响的程度小于气温。
河北张北海流图和山西夏县的H2浓度与气温、气压的综合变化曲线如图5所示。从图中可见,在1年内,H2浓度和气温呈正弓型变化,与气压则呈反弓型变化且其变化,趋向都是可对比的。
图4 不同观测点H2浓度与气压相关关系
Fig.4 The correlation of hydrogen levels and air pressure at different observation points
图5 河北张北海流图(a)和山西夏县(b)H2浓度与气温、气压综合变化曲线图
Fig.5 Ocean current chart of Zhangbei,Hebei and composite change chart of hydrogen levels,air temperature and air pressure of Xiaxian,Shanxi
对不同测点H2浓度与气温、气压进行二元回归拟合,拟合曲线如图6所示,相关系数见表1。从图6和表1可见,所有观测点的H2浓度与气温、气压之间一律呈正相关; 有13个测点的相关系数大于0.2; 张北海流图测点相关系数最大,为0.947 32,这表明该测点受气温、气压共同作用,而其他干扰因素很小。由此可知,H2浓度与气温、气压之间的相关关系是显著的。因此,在异常分析时,排除了气温、气压影响,就可以识别异常。
通过一元回归、二元相关分析方法,定量分析不同深度的观测孔气温、气压对H2浓度的影响,观测孔深度见表1,观测孔深度与气温、气压之间关系如图7所示,由图可见:
(1)如果H2浓度与气温为正相关,则两者间呈反变关系,即观测孔深度越大,H2浓度与气温的相关系数越小; 如果H2浓度与气温为负相关,则两者呈正变关系,即观测孔深度越大,H2浓度与气温的相关系数也越大。
(2)如果H2浓度与气压为正相关,则两者呈正变关系,即观测孔深度越大,H2浓度与气压的相关系数越大; 如果H2浓度与气压为负相关,则两者呈反变关系,即观测孔深度越大,H2浓度与气压的相关系数越小。
(3)无论观测孔为任何深度,H2浓度与气温、气压的复相关系数均为正值。
图6 不同测点H2浓度与气温、气压拟合曲线
Fig.6 The fitted curves of hydrogen levels and air temperature and air pressure at different observation points
图7 观测孔深度与气温(a)、气压(b)以及气温和气压(c)相关系数之间关系
Fig.7 The relationship of correlation coefficients between the depth of observation hole and air temperature(a),and air pressure(b), and between the temperature and the air pressure(c)
(4)当观测孔深为4~11 m时,H2浓度与气温、气压之间相关系数的分布较为集中; 若孔深小于4 m,则比较分散。
为了分析断层带土壤H2浓度变化及其影响因素,本文选取不同时段、不同数据类型、不同观测深度、不同观测环境的15个观测点,通过定性和定量分析研究,得出以下结论:
(1)断层带土壤H2浓度具有明显的日变化,涨落现象十分突出,变化规律清晰; 月动态变化则不明显。
(2)气温和气压是土壤H2浓度观测的主要气象影响因素; 就其影响程度而言,气温的影响大于气压。
(3)土壤H2浓度与气温、气压之间的相关关系是比较复杂的。土壤H2浓度与气温或与气压之间,均可呈现正相关或负相关2种情况。但土壤H2浓度与气温、气压之间的复相关系数却为正值。
(4)土壤H2浓度与观测孔深度间的关系比较明显,也比较复杂,4~11 m可视为观测孔的优选深度。
断层带土壤H2观测浓度技术目前处于试验性探索阶段,还有一些问题有待深入研究。但气象因素中对断层气组分最主要的影响是气温,其次是气压,而且气温、气压对土壤H2浓度的影响程度会受到不同测点的地质环境(主要指观测层的土壤、岩石的成分和结构等)、土壤H2测点与断裂带距离等因素的影响。因此,在分析异常与地震的关系时,应将气象因素引起的变化予以识别或排除。目前痕量氢数字化连续观测处于实验阶段,研究成果存在一定的误差,也需做进一步深入分析。