1. 引言
自1966年邢台地震后,我国便开展了各项前兆观测工作,地电阻率观测作为最先开展的观测项目之一,大多数台站的地电阻率年变形态呈“冬高夏低”的“正常年变”,台址表层介质的地电阻率相对变化的影响系数为正,部分台站年变呈“冬低夏高”的“反向年变”形态 [1],固原台就是反向年变较为典型的台站之一。同属宁南地区的海原台,其地电阻率是近几年新架设的观测项目,数据质量较为稳定,由于固原和海原地电阻率观测站直线距离约为80公里,相距较近,因此,对两者数据进行相关性对比分析,有助于我们对该地区及邻区的地电阻率观测状况和数据变化信息进行系统性观测分析 [2],从而为宁南地区地电阻率数据在异常信息的提取和研判方面提供一个可信度较高的依据,为以后的地震预报工作提供一个有力的支持。
2. 台站地质构造背景简介
2.1. 固原地电台地质构造背景
固原地电台位于固原市原州区彭堡乡彭堡村,距离固原市区14 km。该台地处宁夏南部山区,属大陆干旱性气候。多年平均降雨量为478 mm,年平均蒸发量为1772 mm左右,年均气温6.2℃,年温差在27℃左右。测区位于清水河压陷型盆地内,清水河盆地展布于香山–桃山、马东山–六盘山山脉与烟筒山–窑山和由“南北古脊梁”构成的残山丘陵、黄土丘陵之间,长约180 km,宽5~13 km,海拔1260~1400 m。基底由古近系寺口子组、清水营组和新近系红柳沟组构成,盆地充填物以第四系冲、洪积相砂砾石–粘质砂土、砂质粘土为主,厚10~400 m,构成清水河盆地的主要含水岩组。盆地内的次级地貌主要为清水河冲积平原,地电观测区在清水河古冲沟塬面上,地势相对平坦。
2.2. 海原地电台地质构造背景
海原地电台测区位于海原盆地西北边缘,海原盆地展布于南、西华山北麓,走向北西西,南界为海原断裂带,北侧可能亦为断裂控制,沉降中心位于山前地带,地质构造上属于全新统冲积物层,海拔1805 m。台站所属地属宁夏中部干旱带,大陆性季风气候明显,特点是春暖迟、夏热短、秋凉早、冬寒长。年均气温7℃,一月均温−6.7℃,七月均温19.7℃。年降水量,多年平均286 mm,最多706 mm,最少325 mm。该场地在地貌上处于海原盆地之内,周围地势平坦开阔,测区为干旱农田,地形高差小于20 m。勘选场地位于海原断裂北侧,距离海原断裂10.5 km。据钻孔资料,盆地第四系厚度超过430 m,表层为上更新统风积黄土,观测环境较好。
3. 宁南地电阻率年变趋势对比分析
固原彭堡地电测区位于清水河压陷型盆地内,存在的主要断裂为海原~六盘山东麓断裂。观测区域南偏西约6 km是海原8.5级地震断裂带,北偏西约5 km是清水河断裂带。海原高台地电台位于海原盆地西北边缘,海原盆地展布于南、西华山北麓,走向北西西,南界为海原断裂带。两台站相距约80 km (如图1所示),均邻近海原断裂带和六盘山断裂带,地理位置特殊,两地地电阻率数据具有很好的对比可观测性 [3],掌握它们联合特性对异常信息的研判能够提供可信度较高的依据,因此以下我们首先对宁南地区地电阻率数据进行整体年变形态进行对比分析。
Figure 1. Distribution of georesistivity stations in southern Ningxia
图1. 宁南地区地电台分布图
3.1. 年变幅度进行统计分析
首先对固原和海原地电阻率逐年的年变幅度进行统计分析,统计公式为
(
为年变幅度,
为地电阻率年变最高值,
为地电阻率年变最低值),发现固原地电阻率的年变幅度整体要高于海原地电阻率的年变幅度,幅度比大于1.5的占到全部统计结果的41.6%,幅度比小于1.0的只有16.6% (如表1所示)。
Table 1. Comparison of annual variation of georesistivity in southern Ningxia (unit: Ω∙m)
表1. 宁南地区地电阻率年变幅度对比表(单位:Ω∙m)
从整体变化趋势来看,固原和海原东西向数据年变幅度趋势变化较为一致,海原地电阻率北南向变化趋势与其它测向相反 [4],进行反向处理后,数据变化趋势与固原地电阻率北南向变化趋势基本吻合(如图2所示);地电阻率东西向两测项数据年变趋势2016年7月之前变化相对一致,自2016年下半年以来变化趋势逐渐分离(如图2所示),其原因可能是固原东西向地电阻率自2016年下半年以来长期受环境因素干扰有关,造成固原东西向地电阻率数据出现显著的下降趋势性变化,且变化幅度较大,致使整体年变幅度较以往有较大幅度变化,呈现破年变变化。
3.2. 宁南地区地电阻率相关性分析
固原地电阻率和海原地电阻率东西向观测数据都是呈现“夏高冬低”的“反向年变”形态 [5],只有海原地电阻率北南向观测数据呈现的是“冬高夏低”的年变形态 [6] [7] [8]。
对于两台站地电阻率数据均呈现正弦波式变化,因此我们对两台站相同测道的数据进行相关性分析,首先对海原地电阻率北南向数据进行反向坐标变换,然后与固原地电阻率北南向数据进行对比分析,由图3可以看出两者数据的变化整体趋势一致,年变态势具有较高的相关性,特别是在一些典型事件如数据的同时下降或上升等破年变现象出现时,可以排除是受测区环境干扰或者设备系统故障等局部因素所造成的数据变化,如图3所示,2017年7、8月间固原和海原数据同时出现的上升和下降破年变变化,2018年7月、2019年5月和2020年6月两地数据也同时出现的破年变变化,这些数据变化都可以推断为地表介质的运动造成本地区介质电性结构的变化 [9],只是由于能量的累积问题,有些此类变化并没有能够转化为地震。
在对地电阻率东西测向数据进行分析时也是采用以上相同的推论依据进行判定,但是在对典型事件的分析中,不论是北南向还是东西向数据如图4所示,固原和海原的数据变化趋势基本都是呈反向变化,这一现象需要我们在以后的工作中做进一步分析。
4. 宁南地区地电阻率差异性分析
由于宁夏南部地区属于阴湿和半阴湿地区,因此全年降雨相对较多,特别是在夏季降水量较大,降雨较为频繁,因此夏季降雨等环境干扰因素对地电阻率观测会产生一定程度的影响 [5] [10] [11] [12],造成夏天数据的离散度和波动性较大 [13];冬季宁南地区降雪量较少,雪水不易融化入地,因此对地电阻率的干扰较小,冬季数据也较为稳定 [14]。因此我们对两台站夏季和冬季的地电阻率数据变化情况做一些对比分析,研究它们之间的一些差异性变化,从而为数据的变化分析提供依据。
Figure 2. Comparative analysis of the annual variation of georesistivity in NS and EW directions. (Upper) NS direction; (Below) EW direction
图2. 北南向和东西向地电阻率年变幅度对比分析图。(Upper) 北南向;(Below) 东西向
Figure 3. Comparative analysis of NS georesistivity correlation in southern Ningxia
图3. 宁南地区北南向地电阻率相关性对比分析图
Figure 4. Comparative analysis of EW georesistivity correlation in southern Ningxia
图4. 宁南地区东西向地电阻率相关性对比分析图
4.1. 加速度分析
首先我们运用二阶差分函数建立一个加速度计算模型 [15],以此来分析研究固原和海原地电阻率数据之间及夏季和冬季之间数据变化的差异性,对比量化两者之间的差异指标 [16];在此我们定义X(t)是一个基于地电阻率
取样均值
的时间t相关离散函数:
离散函数
;
为随机变量;
为一阶差分函数;
为二阶差分函数。
通过计算出地电阻率测值
变化的速度为函数Y(t)值,地电阻率测值
变化的加速度为函数Z(t)值,产出的加速度即表明一定时间段内岩石介质的电性结构是否稳定,从而为介质的变化提供依据。因此本课题对固原和海原北南和东西两测向的时测值做加速度分析,得出表2数据内容。表2数据显示每年固原和海原各测道地电阻率数据加速度值固原要显著低于海原,固原NS测向加速度区间为0.026 Ω∙m~0.040 Ω∙m,海原NS测向加速度区间为0.026 Ω∙m~0.049 Ω∙m,除个别年份外固原比海原普遍都低于20%以上,最大为58.0%;EW测向固原加速度值全部低于海原,固原EW测向加速度区间为0.023 Ω∙m~0.042 Ω∙m,海原EW测向加速度区间为0.038 Ω∙m~0.059 Ω∙m,前者比后者变化加速度值平均降低了52.3%,最大为77.8%。
Table 2. Guyuan-Haiyuan georesistivity data acceleration value (unit: Ω∙m)
表2. 固原–海原地电阻率数据加速度值(单位:Ω∙m)
对固原和海原北南向地电阻率数据进行对比加速度分析(如图5所示),在相同计算参数及绘图参数条件下,可以由下图计算结果看到固原地电阻率北南向的整体变化离散度要小于海原地电阻率北南向数据。对固原和海原东西向地电阻率数据进行对比加速度分析,在相同计算参数及绘图参数条件下,可以由下图计算结果看到海原地电阻率东西向的整体变化离散度要大于固原地电阻率东西向数据。
对于地电阻率冬季和夏季的数据处理分析 [17],我们是提取了5月至9月为夏季数据区间,11月至次年3月为冬季数据区间来对地电阻率区间内的数据进行加速度计算求平均值(如表3和表4所示)。
对于固原和海原地电阻率冬夏季数据加速度值进行计算分析后,由表可以看出固原地电阻率不管是冬季还是夏季加速度值都是显著低于海原地电阻率,NS测向冬季固原比海原低11.2%,夏季低20.5%,EW测向冬季固原比海原低35.6%,夏季低46.7%。由以上计算分析得出(如图6和图7所示),固原地电阻率观测数据的离散度要显著低于海原地电阻率,即使是在固原地电阻率观测环境较海原复杂,干扰因素更多的情况下,也要高出20至40个百分点。因此,推断出固原地电阻率对于异常信息的捕获能力要高于海原地电阻率观测数据,从而具有更好的分析异常变化的研判能力。
Figure 5. Guyuan-Haiyuan NS and EW directions of georesistivity data acceleration analysis. (Upper) NS direction; (Below) EW direction
图5. 固原–海原北南向和东西向地电阻率加速度分析图。(Upper) 北南向;(Below) 东西向
Table 3. Guyuan winter and summer georesistivity data acceleration value (unit: Ω∙m)
表3. 固原冬夏季地电阻率数据加速度值(单位:Ω∙m)
Figure 6. Analysis of NS direction georesistivity acceleration in winter and summer in southern Ningxia. (Left) Data of Guyuan; (Right) Data of Haiyuan
图6. 宁南地区冬夏季地电阻率北南向加速度分析图。(左) 固原数据;(右) 海原数据
Figure 7. Analysis of EW direction georesistivity acceleration in winter and summer in southern Ningxia. (Left) Data of Guyuan; (Right) Data of Haiyuan
图7. 宁南地区冬夏季地电阻率东西向加速度分析图。(左) 固原数据;(右) 海原数据
Table 4. Haiyuan winter and summer georesistivity data acceleration value (unit: Ω∙m)
表4. 海原冬夏季地电阻率数据加速度值(单位:Ω∙m)
4.2. 地电阻率测值比分析
本课题首先对固原和海原地电阻率北南和东西向数据做测值比分析,计算公为
(
为测值比,
为北南向地电阻率,
为东西向地电阻率), 以此来分析研判该地区地层介质地电阻率的方位性变化,并通过相同地点地电阻率的正交方向分量变化,拟合成地电阻率虚拟方位角的变化。
图8左1图是2014年到现在固原地电阻率计算出的测值比结果,由图可以看出数据一致在上升,这主要是受东西向数据自2016年下半年以来的持续性下降所致。图8右1图是海原地电阻率计算出的测值比结果,由图可以看出数据虽然出现周期性变化,但整体变化趋势是呈现线性的下降变化,出现周期性变化是由于北南向数据和东西向数据呈反向变化所致;数据的持续下降推断出地层的表面介质在东西方向上承受的压力在逐渐增大,或者推断为在南北方向上承受的压力在减小。
图8右2图是固原和海原北南向地电阻率的测值比,变化整体呈上升趋势,出现正弦波式变化也是因为海原地电阻率北南向出现与其它测道呈反向年变形态所致。图8左2图是固原和海原东西向地电阻率的测值比,变化整体自2016年下半年开始呈线性下降态势,受固原地电阻率东西向数据下降变化所致。
Figure 8. Analysis of the measured value ratio of georesistivity in southern Ningxia
图8. 宁南地电阻率测值比分析图
对于宁南地区地电阻率数据的测值比计算结果显示,在宁南地区及邻区中强地震震前数据都会出现显著破年变变化,测值比震前值较正常值变化幅度超出范围为0.0037~0.0060,对于地方震可能对异常信息的提取更为显著,因此,测值比为震前异常信息的研判提供了一个依据。
5. 结论
1) 通过对宁夏南部地区地电阻率数据整体年变形态及相关性分析研究,认为在年变形态上固原和海原地电阻率较为一致,相关地电阻率测项之间的数据相关性较高,通过数据之间的相关性分析可以推断出相同测项数据在中强地震震前2至3个月内相关性会显著提高10个百分点左右,对于震前异常信息的提取提供了一个很好的依据。通过对相邻地区地电阻率相关性的分析,了解数据之间存在的异同,在一定程度上提升了对异常信息研判的能力。
2) 通过对固原和海原地电阻率冬夏季数据加速度值进行计算分析,得出固原地电阻率不管是冬季还是夏季加速度值都是显著低于海原地电阻率,NS测向冬季固原比海原低11.2%,夏季低20.5%,EW测向冬季固原比海原低35.6%,夏季低46.7%。从整体数据情况来看,固原地电阻率观测数据的离散度要显著低于海原地电阻率,即使是在固原地电阻率观测环境较海原复杂,干扰因素更多的情况下,也要高出20至40个百分点。因此,推断出固原地电阻率对于异常信息的捕获能力要高于海原地电阻率观测数据,从而具有更好的分析异常变化的研判能力。
3) 通过对宁南地区地电阻率数据的测值比计算,结果显示在宁南地区及邻区中强地震震前数据都会出现显著破年变变化,测值比震前值较正常值变化幅度超出范围为0.0037~0.0060,对于地方震可能对异常信息的提取更为显著,因此测值比为震前异常信息的研判提供了一个依据。
基金项目
宁夏地震局青年基金课题和宁夏地震局地下流体监测与应用研究创新团队资助。