笔者选取宝昌台2012年4月1～30日的地电阻率整点值资料进行分析,发现EW测道具有明显的日变形态,但NS测道该形态并不十分明显。在此采用小波方法对地电阻率数据进行去噪处理。去噪后的地电阻率NS测道整点值曲线,与去噪前相比,对随机波动的抗干扰能力强,且具有明显的日变形态,将其与去噪后的地电EW向整点值曲线进行比较,注意到两者基本重合(图3),说明宝昌地电两测道整点值曲线日变形态基本一致,同时也注意到两测道出现日变畸变时段也较为一致。

图4显示的是宝昌台地电阻率(NS)整点值日变最大值和最小值对应时刻随时间变化关系,2008年1月至2009年6月最大值和最小值时刻虽总体分属于0～11 h、12～23 h两个区间,但是较为离散,2009年7月之后,最大值和最小值时刻趋于集中,其中最大值时刻集中于C区间,而最小值时刻集中于A和B区间,这说明2009年7月之后地电阻率日变形态趋于稳定。

图3 宝昌台地电阻率整点值曲线(2012年04月)
Fig.3 Hour value curve of earth resistivity recorded by Baochang Station(Apr.,2012)

图4 宝昌台地电阻率(NS)日变最大值和最小值对应时刻随时间变化关系
Fig.4 The corresponding moment of the maximum and the minimun daily of earth resistivity rate recorded by Baochang Station varing with time

宝昌台地电阻率整点值曲线,类似于大地电场等,也存在规则的日变形态(李飞等,2011),可能是由一个或者多个具有周期特征的影响因素对地电阻率独立或者叠加施加作用产生的。

首先选取宝昌台地电场、锡林浩特地磁及宝昌台固体潮汐理论值,与同时段的宝昌地电阻率两测道整点值进行对比,结果显示:(1)宝昌台地电阻率日变形态与地电场、地磁形态相差甚远;(2)将宝昌台地电阻率整点值(去噪)与该台固体潮汐理论值以及滤除半日潮汐波之后的固体潮汐理论值进行逐日相关性分析,二者之间具有一定相关性,但关系不显著。

气象因素中,具有短周期性变化特征,且对地电阻率变化具有明显影响的物理量,仅有温度。

(1)通过相关系数方法,定量分析宝昌地电阻率与测区内气温之间的相关性特征,以寻找出制约地电阻率变化的主要因素。选取宝昌台地电阻率和温度2008年1月1日至2012年4月30日整点值数据,首先利用小波方法对地电阻率进行去噪处理,在此基础之上逐日计算地电阻率与温度之间的简单相关系数,2012年4月30日温度缺数未参与计算。地电阻率NS测道与温度的相关性结果(图5)显示,相关系数r为负数,且〖JB<1|〗r〖JB>1|〗≥0.60的天数占参与计算的总天数的49.34%,显著性水平p≤0.05的天数占参与计算的总天数的70.84%; 地电阻率EW测道与温度的相关性结果显示,相关系数r为负数,且〖JB<1|〗r〖JB>1|〗≥0.60的天数占参与计算的总天数的37.63%,显著性水平p≤0.05的天数占参与计算的总天数的68.31%。统计结果表明,宝昌台地电阻率与测区温度之间的日相关主要是以负相关为主,且相关关系显著。

(2)仅从数学角度定量研究,结果显示温度是导致宝昌台地电日变的主要原因,但这是否是真实情况的反应,还需作进一步讨论:① 宝昌台地电阻率自2006年开始采用ZD8B地电仪,观测室内温度为18～22 ℃。地电外线路主要按照NS、EW两测道布线,其中两测道供电极距和测量极距均相同,分别为560、80 m; 2005年9月对外线路进行了改造,外线路全部更换为铜芯屏蔽电缆,电极全部更换为新的铅电极,电极埋深3.5 m; 改造初期由于用于固定电缆的钢绞线与电线杆之间未用瓷瓶隔开,导致下雨时地电阻率出现初期上升变化,之后才下降的现象,2007年5月采用瓷瓶隔开后,该现象不明显。②自瓷瓶安装之后,至今没有清洁过,这可能是导致地电阻率日变的主要原因:由于瓷瓶未定期清洁,外线系统在对地绝缘方面稳定性差,每天进入晚间之后随着温度的降低,空气中冷凝水附着在瓷瓶上造成钢绞线与大地连通,电阻率整体上升,进入白天之后温度升高瓷瓶绝缘性增加,地电阻率整体下降。为验证这一分析是否准确,2012年8月26、27日对NS测道的瓷瓶进行了清洁,之后地电阻率NS向日变虽仍存在,但幅度有所减小(图6)。

图5 宝昌台地电阻率NS测道与温度整点值逐日相关统计结果
Fig.5 Daily correlation stastics of the earth resistivity rate in NS-direction and tempreture hour value recorded by Baochang Station

图6 宝昌台地电阻率NS测道整点值曲线
Fig.6 Hour value curve of earth resistivity rate in NS-direction recorded by Baochang Station

位于宝昌台测区南端有一口水井,记录了2004年4月至2008年7月地下潜水位埋深数据,宝昌台测区潜水位埋深平均为6.21 m,每年因灌溉等有波动,变化幅度为14%。依据宝昌台钻孔剖面图及电测深曲线,该台地下电性结构可分3层,第一层主要是高电阻率的砂碎石层(0～8.5 m),第二层是由精砂、细砂和砂土组成的低电阻率层(8.5～62.5 m),第三层是由砂砾岩及石英斑岩组成的高电阻率层(62.5 m以下)。宝昌台水位虽逐年趋于下降,但这并不是导致地电阻率NS测道下降的主要因素,因为宝昌台水位主要在第一层内变化,对决定该台地电阻率长期变化的第二、三层电阻率影响较小。

为确定宝昌台地电长期变化是否由降雨量、温度等气象因素引起的,在此将宝昌台地电阻率NS测道1993年1月至2008年12月月均值数据ρ与同期的降雨量Y、温度W进行了线性回归,得到的回归关系为ρ=a₀+a₁Y+a₂W,其中a₀=142.626 7,a₁=-0.001 0,a₂=-0.052 1,并依据拟合关系给出了仅由降雨量、温度引起的地电阻率NS测道变化曲线(图7)。结果显示,由气象因素引起的地电阻率变化没有明显的趋势下降特征。这说明宝昌地电NS测道长期变化有可能是由于区域应力作用产生的。

宝昌台所在的华北地区构造复杂,断裂纵横交错,形成了一系列大的构造带。徐菊生等(1999)利用GPS观测结果,通过反复调整应力边界条件,计算得到了华北地区的应力、应变、位移数据,结果表明,华北地区现今构造应力场的基本特征是优势压应力方向为北东东。本文采用《中国震例》提供的资料(张肇诚,1990; 张肇诚等,2000; 陈棋福等,2003),统计了宝昌台周边区域典型地震震源机制解特征(表2),结果表明地震震源机制反映出的主压应力方位均为NEE-SWW。因此可确定宝昌台所在区域受到的主压应力是近EW向的。

在地震孕育发展阶段,直至地震发生过程中,通常可以观察到位于孕震区的地电存在各向异性变化现象,很多学者通过实验或者模型,对该现象进行了验证和解释(毛桐恩等,1999; 冯志生等,2004),笔者借用杜学彬等(2001)解释强震近震中区地电变化速率的各向异性理论,来分析宝昌地电NS测道长期变化特征。宝昌台所在区域由于长期受到近EW向主压应力持续作用,造成介质内部裂隙走向逐渐沿主压应力方向优势排列,导电通道连通、导电流体(如水)进入或重新分布,造成以主压应力方向真电阻率持续下降为主的真电阻率变化各向异性,导致以垂直主压应力方向的宝昌台地电NS测道持续下降为主的视各向异性。

图7 宝昌台地电阻率NS测道理论值曲线
Fig.7 Theoretical value of earth resistivity rate in NS direction at Baochang Station

表2 宝昌台周边区域典型地震应力特征
Tab.2 Typical seismic stress characteristic in Baochang Station and its adjacent area