基于4种技术方法的地下流体资料预报效能探讨

2023-03-06胡小静付虹张翔李琼高文斐

中国地震 2023年4期

胡小静付虹张翔李琼高文斐

云南省地震局,昆明 650224

0 引言

地下流体普遍存在于地壳介质中,被国内外学者认为是最重要的短临前兆观测手段之一(张炜等,1988; 刘耀炜,2006;Wang et al,2009;),在长期的地震预报实践中,所记录到的各类流体异常在大震短临预测过程中发挥了重要的作用(刘耀炜等,2002; 吴忠良等,2007; 晏锐等,2018; 田雷等,2022)。众多研究结果表明,地下流体资料对孕震区和外围区介质变化的响应非常灵敏,地震的孕育和发生过程始终伴随有地下物质运移、能量传输和条件改变,从而引起不同含水层之间的交互渗入等,进而导致地下水物理、化学参量和逸出气体浓度的异常变化(车用太等,1997;Montgomery et al,2003;Zhou et al,2021; 李营等,2022;Chen et al,2022)。

研究结果显示,2014年鲁甸6.5级地震发生前,昭通地区地下流体观测记录到了大量的地下流体异常,对滇东北地区乃至川滇交界东部地区的震情趋势判定发挥着重要的作用(刘耀炜等,2015)。因此,本文对昭通地区地下流体资料开展震情跟踪技术方法预报效能计算,尽可能从定量化的角度来评估观测资料在震情预测方面发挥的作用,同时对这些异常产生的孕震机理进行讨论。

1 资料概况

昭通地区目前正在观测的全部地下流体资料包括19个观测点的47个观测项,其中模拟观测17项,包括水氡观测5项、二氧化碳观测5项、离子观测4项、电导率观测2项、pH观测1项,数字化观测30项,均为水温水位观测。梳理全部观测资料后,巧家老店水温、盐津水温、永善莲峰水温及水位、鲁甸水温及水位由于观测时间较短,无法参与计算,其余41个测项全部参与了阈值法、速率变化法、趋势转折法和破年变法4种方法的预报效能计算。参与本次预报效能检验的震例为2000年以来昭通地区及周边200km范围内的5级以上历史地震,具体震例分布情况见图1,测点基本信息见表1。

图1 昭通地区地下流体观测点及震例分布

2 计算方法

分别采用阈值法、速率变化法、趋势转折法和破年变法4种计算方法对昭通地区地下流体观测资料进行异常提取,并计算其R值进行预报效能检验。

(1)阈值法,该方法直接从原始观测时序数据中识别出异常信息,也称为原始曲线法。首先求取整个时序数据的标准差,计算最佳阈值(标准差倍数n); 然后查询超出该阈值的异常数据,得到观测数据中的高(低)值异常及其时间点。

(2)速率变化法,又称差分法,是一种压抑长周期、突出短周期变化的线性滤波。针对原始观测数据的一阶差分序列,使用平稳序列的均方差作为异常控制线,或使用给定的阈值作为异常控制线来判别时间序列中的异常点,突出那些突跳性或离散度较大的异常。当有一组数据(x1,x2,…,xn)时,可计算相邻2个数据之间的差值(Δxi),并依此按下式计算差分值(绝对值)的平均值

(1)

(3)趋势转折法,该方法中的观测数据趋势项基于线性拟合方法提取,利用最小二乘法对观测时序数据进行线性拟合,得到观测数据的最佳转折时间点,根据线性拟合速率值出现变化或变化量超出给定阈值来判定是否为趋势转折异常。具体如下:给定数据点(xi,yi)(i=0,1,…,m),假定在Kk(xk,yk)(k=1,2,…,m-1)点处出现趋势转折,利用最小二乘法求得Kk点前后数据线性变化的斜率kk1和kk2,相应的残差为ek1和ek2,残差和(ek1+ek2)最小值所对应的点即为最佳数据转折点。

(4)破年变法,年变形态或幅度异于多年平均变化的异常,称之为破年变异常,一般采用矩平分析法识别该类异常。具体为:用多年观测数据的均值序列(五日、旬、月)之间的相关程度来求观测点值的正常年动态曲线,再以原来不同年时间段的均值序列和正常年动态曲线的余差序列提取异常指标,用其剩余标准差的2倍或3倍作为阈值控制线。

(5)R值评分,本次预报效能评估是同时进行时间和空间预测的方法,其R值评分计算公式为(许绍燮,1989)

(2)

当R>0时,表示该方法预测效能高于随机预测,具有预测意义; 当R≥R0时,表明该方法预测效能通过置信度97.5%的显著性检验,具有显著预测意义。

根据历史震例的选取原则,空间网格数的确定以震中为圆心,5级地震半径为200km的范围,6级地震半径为300km的范围,7级地震半径为300km的范围。

3 预报效能计算结果

针对昭通地区现有的43项地下流体观测数据,分别采用4种震情跟踪技术方法提取异常,并计算其预报效能(具体结果见表1),不同方法具体评估结果如下。

表1

阈值法计算的预报效能R值结果(图2(a))显示,水富水位、巧家电导率、巧家镁离子、乐居水位、茨院水位共计5个测项R≥R0,通过预测效能置信度97.5%的显著性检验,在日常跟踪分析过程中具有较好的预测意义; 82%的测项具有一定的预测意义; 19%的测项预报效能低于随机预测,在日常跟踪过程中无参考价值。从最佳预测时间来看,预测时段在3个月内的测项占比约61%,3～6个月的测项占比约9%,半年以上的测项占30%,表明利用阈值法提取的异常具有较好的短临预测意义,其次可作为半年以上的年尺度指标。

图2 不同方法预报效能评估结果

速率变化法计算的预报效能R值结果(图2(b))显示,水富水位、水富水温、巧家电导率、巧家镁离子、巧家钙离子、盐津水氡、茨院水位共计7个测项R≥R0,通过预测效能置信 97.5%的显著性检验,在日常跟踪分析过程中具有较好的预测意义; 68%的测项具有一定的预测意义; 16%的测项预报效能低于随机预测,在日常跟踪过程中无参考价值。从最佳预测时间来看,预测时段在3个月以内的测项占比约70%,3～6个月和半年以上的测项分别占12%和18%,表明利用速率变化法提取的异常其短临预测意义更为显著,可作为3个月以内显著的短临指标。

破年变法计算的预报效能R值结果(图2(c))显示,水富水位、水富水温、巧家镁离子、巧家碳酸根、大关电导率、鲁甸二氧化碳、昭阳一中水位、茨院水位、茨院水温共计9个测项R≥R0,通过预测效能置信度97.5%的显著性检验,在日常跟踪分析过程中具有较好的预测意义; 65%的测项具有一定的预测意义; 14%的测项预报效能低于随机预测,在日常跟踪过程中无参考价值。从最佳预测时间来看,预测时段在3个月以内的测项占比56%,具有较好的短临预测意义,6个月及以上的测项占比32%,可作为半年或年度指标。

趋势转折法计算的预报效能R值结果(图2(d))显示,水富水位、巧家钙离子、巧家碳酸根、大关电导率、鲁甸二氧化碳共计5个测项R≥R0,通过预测效能置信度97.5%的显著性检验,在日常跟踪分析过程中具有较好的预测意义; 59%的测项具有一定的预测意义; 约30%的测项预报效能低于随机预测,在日常跟踪过程中无参考价值。从最佳预测时间来看,预测时段在3个月以内的测项占比约70%,表明趋势转折后将进入显著的短临阶段。

综合上述4种方法的结果分析,速率变化法和破年变法提取的异常预报效能更佳; 水富水位、水富水温、巧家电导率、巧家钙离子、巧家镁离子、巧家碳酸氢根、大关电导率、鲁甸二氧化碳、茨院水位、茨院水温、昭阳一中水位等共计11个测项可以作为滇东部北地区年度或日常跟踪的重点测项; 大多数流体异常的最佳预测时段在3个月以内,尤其是速率法提取的异常和趋势转折后会进入较为显著的短临预测阶段; 除此之外,约有22%的测项通过震例评估的预报效能低于随机预测(R<0),对震情跟踪判定的预测意义不大。

4 讨论与结论

4.1 讨论

针对上述提取出的预报效能较好的测点测项,选择最佳的计算方法结果(表2),依据不同手段的物理意义、异常特征以及流体补给来源,并结合已有的研究结果对这些异常产生的孕震机理进行讨论。

表2 昭通地区重点跟踪测项的预报效能信息

从不同观测测项异常的最佳预测时段来看,不同测项的异常具有一定的阶段性特征,其中水温和水化学离子异常大多数出现在震前2个月以内,表现出较为明显的短临特征; 而水位、电导率异常大多数出现在震前半年或一年左右,表现出一定的中期到年度的前兆特征。从提取出的异常形态来看,水温和水化学离子主要表现为测值的速率变化,反映的是震前观测项在短周期范围的高频异常信息; 水位、电导率主要表现为原始测值正常年变超出正常阈值范围或者趋势形态发生改变,主要反映的是观测项在长周期范围的变化过程。整体来看,不同的观测手段呈现出较为显著的时间层次性和形态差异性。

注: (a)水化学Piper图; (b)Na-K-Mg三角图; (c)氢氧同位素组成图。

研究表明,降水等因素造成的地下水动态变化可产生一种作用于地壳岩石的附加流体力,地震前多种流体力的作用过程是环境因子对浅层前兆异常作用的主要力学机制之一,孔隙压力、动水压力以及化学腐蚀等原理,是该异常机理的理论基础。这种浅层诱发前兆异常非常显著的一个特征是异常变化的层次性,即在这种作用力的影响下,会产生不同时间尺度的异常变化(王吉易等,2002)。此外也有研究表明,2014年8月3日鲁甸6.5级地震前,昭通地区的流体异常产生的原因可能是断裂带空隙的开启与闭合,导致观测井含水层承压性发生改变,引起地下水混合作用增强(刘耀炜等,2015),该研究结果与上述机制的理论基础相符。因此,结合昭通地区不同异常所呈现出的较为显著的时间层次性和形态差异性,以及大气降水补给的水化学特征,认为昭通地区的地下流体异常大多可能属于这种地下水诱发的浅层前兆异常。