宝昌台地电阻率变化特征.doc

1、1宝昌台地电阻率变化特征摘要：对宝昌台地电阻率日变、年变和长期变化 3 种典型变化进行了分析，结果显示：整点值曲线存在明显的日变形态，可能是由外线系统对地绝缘不良造成的；地电阻率年变主要是由气象因素年变引起的，其中温度和降雨量贡献基本相同；地电阻率 NS 测道存在长期下降变化，而 EW 测道没有，这一现象与测区内水位、温度和降雨量等长期变化无关，主要是由于受到台站所在区域近 EW 向应力长期、持续作用引起的，并不对应地震。 关键词：宝昌台；地电阻率；特征变化；机理 中图分类号：P319 文献标识码：A 文章编号：1000-0666（2013）03-0358-06 0 引言 由于气象、环境等因素

2、影响，地电阻率观测数据通常情况下并不是一条平直的曲线，而是具有趋势变化、年变和脉冲等多种变化形态的曲线。与地震孕育、发生有关系的地震前兆异常因上述变化的存在，通常被弱化甚至被掩盖，因此分析地电阻率曲线典型变化特征，并研究其产生机理，对于识别和提取前兆异常具有非常重要的意义。宝昌台（4190N，11527E）位于内蒙古太仆寺旗，地貌为低山丘陵地带，覆盖层较浅（内蒙古自治区地震局，2006） ，该台地电阻率自 1980 年投入观测，对曾发生于晋冀蒙交界地区的地震，在震前有不同程度异常反映2（高立新，黄根喜，1999；黄根喜，2000；戴勇等，2009） 。本文主要通过小波变换、相关分析等方法研究宝

3、昌台地电阻率存在的典型变化特征，并从气象因素、区域应力场等角度解释其产生机理。 1 地电阻率变化特征 宝昌台地电阻率变化主要有非周期性变化和周期性变化。典型的非周期性变化主要为地电阻率 NS 测道，自 1993 年至今呈现长期下降变化（图 1） 。而处于同一测区的地电阻率 EW 测道观测数据却没有出现类似NS 测道的显著的趋势下降，从这一角度来说，宝昌台地电阻率长期变化具有各向异性变化。 为准确获取宝昌台地电阻率存在的周期性成分，采用快速傅里叶变换方法（FFT） （万永革，2007）对 2008 年 1 月 1 日至 2012 年 4 月 30 日的地电阻率整点值数据和 1993 年 1 月

4、1 日至 2012 年 4 月 30 日的日均值数据分别进行振幅谱分析。结果显示，地电阻率两测道整点值数据对应振幅谱都存在一个变化幅度大的单峰值形态（图 2） ，其峰值频率一致，都为 115710-5 Hz（对应周期 24 h） ，这说明宝昌台地电阻率整点值曲线存在典图 1 宝昌台地电阻率日均值曲线 Fig1Daily mean value curve of earth resistivity rate recorded by Baochang Station 型的日变特征；地电阻率两测道日均值数据对应振幅谱也都存在一个变化幅度大的单峰值形态，其峰值频率一致，都为 311510-8 Hz（对应

5、周期是 372 d） ，这说明宝昌台地电阻率日均值曲线存在典型的年变特征。图 2 宝昌台地电阻率整点值振幅谱 32012 年 4 月 130 日的地电阻率整点值资料进行分析，发现 EW 测道具有明显的日变形态，但 NS 测道该形态并不十分明显。在此采用小波方法对地电阻率数据进行去噪处理。去噪后的地电阻率 NS 测道整点值曲线，与去噪前相比，对随机波动的抗干扰能力强，且具有明显的日变形态，将其与去噪后的地电 EW 向整点值曲线进行比较，注意到两者基本重合（图 3） ，说明宝昌地电两测道整点值曲线日变形态基本一致，同时也注意到两测道出现日变畸变时段也较为一致。 2 影响因素分析 宝昌台地电阻率整点

6、值曲线，类似于大地电场等，也存在规则的日变形态（李飞等，2011） ，可能是由一个或者多个具有周期特征的影响因素对地电阻率独立或者叠加施加作用产生的。 首先选取宝昌台地电场、锡林浩特地磁及宝昌台固体潮汐理论值，与同时段的宝昌地电阻率两测道整点值进行对比，结果显示：（1）宝昌台地电阻率日变形态与地电场、地磁形态相差甚远；（2）将宝昌台地电阻率整点值（去噪）与该台固体潮汐理论值以及滤除半日潮汐波之后的固体潮汐理论值进行逐日相关性分析，二者之间具有一定相关性，但关系不显著。 气象因素中，具有短周期性变化特征，且对地电阻率变化具有明显影响的物理量，仅有温度。 （1）通过相关系数方法，定量分析宝昌地电阻

7、率与测区内气温之间的相关性特征，以寻找出制约地电阻率变化的主要因素。选取宝昌台地电阻率和温度 2008 年 1 月 1 日至 2012 年 4 月 30 日整点值数据，首先利4用小波方法对地电阻率进行去噪处理，在此基础之上逐日计算地电阻率与温度之间的简单相关系数，2012 年 4 月 30 日温度缺数未参与计算。地电阻率 NS 测道与温度的相关性结果（图 5）显示，相关系数 r 为负数，且 r060 的天数占参与计算的总天数的 4934%，显著性水平 p005 的天数占参与计算的总天数的 7084%；地电阻率 EW 测道与温度的相关性结果显示，相关系数 r 为负数，且 r060 的天数占参与计

8、算的总天数的 3763%，显著性水平 p005 的天数占参与计算的总天数的 6831%。统计结果表明，宝昌台地电阻率与测区温度之间的日相关主要是以负相关为主，且相关关系显著。 （2）仅从数学角度定量研究，结果显示温度是导致宝昌台地电日变的主要原因，但这是否是真实情况的反应，还需作进一步讨论： 宝昌台地电阻率自 2006 年开始采用 ZD8B 地电仪，观测室内温度为 1822 。地电外线路主要按照 NS、EW 两测道布线，其中两测道供电极距和测量极距均相同，分别为 560、80 m；2005 年 9 月对外线路进行了改造，外线路全部更换为铜芯屏蔽电缆，电极全部更换为新的铅电极，电极埋深 35 m

9、；改造初期由于用于固定电缆的钢绞线与电线杆之间未用瓷瓶隔开，导致下雨时地电阻率出现初期上升变化，之后才下降的现象，2007 年 5 月采用瓷瓶隔开后，该现象不明显。自瓷瓶安装之后，至今没有清洁过，这可能是导致地电阻率日变的主要原因：由于瓷瓶未定期清洁，外线系统在对地绝缘方面稳定性差，每天进入晚间之后随着温度的降低，空气中冷凝水附着在瓷瓶上造成钢绞线与大地连通，电阻率整体上升，进入白天之后温度升高瓷瓶绝缘性增加，地电阻率整体下降。为验证这一分5析是否准确，2012 年 8 月 26、27 日对 NS 测道的瓷瓶进行了清洁，之后地电阻率 NS 向日变虽仍存在，但幅度有所减小（图 6） 。图 5 宝

10、昌台地电阻率 NS 测道与温度整点值 3 年变特征 宝昌台地电阻率 NS、EW 两测道都具有明显的年变，一般情况下，3月最高、9 月最低，即具有典型的“夏低冬高”的年变形态。众多学者研究表明，地电阻率年变主要是由气象因素等周期性年变引起的（徐世浙，1985；刘允秀，陈华静，1999；李 飞，姚伟中，2004；张学民，2004） 。为定量研究宝昌台地电阻率年变是由气象因素中哪一个物理量主导的，在此选取宝昌台地电阻率 1993 年 1 月至 2008 年 12 月的月均值数据和宝昌台测区同期温度、降雨量数据进行了相关分析。统计结果显示，宝昌台地电阻率与测区温度、降雨量之间均存在较为显著的年相关，且

11、以负相关为主，温度、降雨量在对地电阻率产生年变的贡献上基本相同（表1） 。部分年相关系数较低，甚至出现正相关现象，这可能是由于地震前兆异常或者短期干扰等对地电阻率年变产生影响，导致年变畸变。 同处于内蒙古中部的乌加河台地电阻率曲线具有“夏高冬低”的反向年变特征，正好与宝昌台地电阻率曲线具有的年变特征相反。乌加河台、宝昌台干扰模型虽然同属 a 类干扰（即一层电阻率发生季节性变化的模型） ，但两台站地下电性结构不同，其中乌加河台地电断面为 KQQ 型（中间层电阻率高于表层及底层电阻率） （国家地震局预测预防司，1998） ，宝昌台地电断面为 H 型（中间层电阻率低于表层及底层电阻率）（国家地震局科

12、技监测司，1988） 。依据钱复业等（1987）的研究成果，当乌加河台供电极距在 1802 000 m 之间时，a 类干扰 K 型地电断面出6现明显反常的季节变化（夏高冬低） ，而当宝昌台供电极距为 560 m 时，a 类干扰 H 型地电断面出现明显正常的季节变化。表 1 宝昌台地电阻率与温度、降雨量整点值 4 长期变化特征 宝昌台地电阻率 NS 测道自 1993 年至今一直存在长期下降变化，而EW 测道没有出现长期单调变化。全国地电台网中有很多地电台站存在着类似于宝昌台地电阻率 NS 测道的长期单调变化的现象，如天津宝坻台存在单调下降变化，山西大同台存在单调上升变化等。虽然这一变化和形变数据

13、由于零漂造成的单调变化，在形态上类似，但机理上不一样，该现象的出现并不是仪器本身造成的。对于地电阻率长期变化，赵和云（1994）总结了 26 年来全国近百个台站的地电阻率观测资料，认为“斜”型趋势变化一般不对应地震。在对地电阻率存在的长期单调变化现象异常信度判定，以及是否将其纳入会商判定依据，存在着不同的意见。因此，对地电长期变化产生的机理研究，具有十分重要的意义。本文采用拟合等方法，从气象、应力等角度，结合前人提出的各向异性理论，来解释宝昌台地电阻率 NS 测道长期下降变化机理。 41 地下水位、气象因素对地电阻率影响分析 位于宝昌台测区南端有一口水井，记录了 2004 年 4 月至 200

14、8 年 7月地下潜水位埋深数据，宝昌台测区潜水位埋深平均为 621 m，每年因灌溉等有波动，变化幅度为 14%。依据宝昌台钻孔剖面图及电测深曲线，该台地下电性结构可分 3 层，第一层主要是高电阻率的砂碎石层（085 m） ，第二层是由精砂、细砂和砂土组成的低电阻率层（85625 m） ，第7三层是由砂砾岩及石英斑岩组成的高电阻率层（625 m 以下） 。宝昌台水位虽逐年趋于下降，但这并不是导致地电阻率 NS 测道下降的主要因素，因为宝昌台水位主要在第一层内变化，对决定该台地电阻率长期变化的第二、三层电阻率影响较小。 为确定宝昌台地电长期变化是否由降雨量、温度等气象因素引起的，在此将宝昌台地电阻

15、率 NS 测道 1993 年 1 月至 2008 年 12 月月均值数据 与同期的降雨量 Y、温度 W 进行了线性回归，得到的回归关系为=a0+a1Y+a2W，其中 a0=142626 7，a1=-0001 0，a2=-0052 1，并依据拟合关系给出了仅由降雨量、温度引起的地电阻率 NS 测道变化曲线（图7） 。结果显示，由气象因素引起的地电阻率变化没有明显的趋势下降特征。这说明宝昌地电 NS 测道长期变化有可能是由于区域应力作用产生的。42 应力场对宝昌地电阻率影响分析 宝昌台所在的华北地区构造复杂，断裂纵横交错，形成了一系列大的构造带。徐菊生等（1999）利用 GPS 观测结果，通过反复

16、调整应力边界条件，计算得到了华北地区的应力、应变、位移数据，结果表明，华北地区现今构造应力场的基本特征是优势压应力方向为北东东。本文采用中国震例提供的资料（张肇诚，1990；张肇诚等，2000；陈棋福等，2003） ，统计了宝昌台周边区域典型地震震源机制解特征（表 2） ，结果表明地震震源机制反映出的主压应力方位均为 NEE-SWW。因此可确定宝昌台所在区域受到的主压应力是近 EW 向的。 在地震孕育发展阶段，直至地震发生过程中，通常可以观察到位于8孕震区的地电存在各向异性变化现象，很多学者通过实验或者模型，对该现象进行了验证和解释（毛桐恩等，1999；冯志生等，2004） ，笔者借用杜学彬等

17、（2001）解释强震近震中区地电变化速率的各向异性理论，来分析宝昌地电 NS 测道长期变化特征。宝昌台所在区域由于长期受到近EW 向主压应力持续作用，造成介质内部裂隙走向逐渐沿主压应力方向优势排列，导电通道连通、导电流体（如水）进入或重新分布，造成以主压应力方向真电阻率持续下降为主的真电阻率变化各向异性，导致以垂直主压应力方向的宝昌台地电 NS 测道持续下降为主的视各向异性。图 7宝昌台地电阻率 NS 测道理论值曲线 Fig7Theoretical value of earth resistivity rate in NS direction at Baochang Station 表 2 宝

18、昌台周边区域典型地震应力特征 Tab2Typical seismic stress characteristic in Baochang Station and its adjacent area 华北地区主要地震事件 1 主压应力或主张应力轴方位 1 备注 1976 年4 月 6 日 内蒙古和林格尔 MS63 地震 1NEE 向主压应力，NNW 向主张应力场 1 右旋走滑为主的逆断层 1989 年 10 月 19 日 山西大同阳高 MS61 地震 1 主压应力优势方向为 NEE，主张应力优势方向为 NWW，两个主应力轴方向都是近水平的，与华北 现代应力场的基本特征一致 1998 年 1 月

19、10 日张北 MS62 地震；1999 年 3 月 11 日 9张北 MS55 地震 1 主压应力轴方位 NEE-SWW 向、近于水平 1 最大主压应力轴与张北地震一致；张北地震的 序列的破裂面为 NWW、NNE 方向组合，恰好 是华北构造应力场最大剪切应力方向 1999 年 1 月 29 日 锡林浩特 MS52 地震 1NEE 向主压应力、NNW 主张应力 1 右旋走滑兼正断错动 1999 年 11 月 1 日 大同 MS52 地震 1 主压应力方向为 NE66，主张应力方向为 NW18；主压应力轴仰角近水平，主张应力轴仰角为 321 震源错动以走滑为主 5 结论与讨论 宝昌台位于晋冀蒙交界

20、地区，自投入观测以来地电阻率连续、稳定，典型变化特征明显，且对台站所在区域的地震具有较好的前兆异常反映。宝昌台地电阻率既有时间段较长的模拟数据，也有质量较为理想的整点值数据，气象要素等辅助观测较为完善。宝昌台地电阻率在全国地电台网的地电资料中具有代表性，研究该台地电阻率的典型变化特征，对于今后研究全国其它地电阻率变化特征具有重要意义。 （1）通过快速傅里叶变换方法，分别得到宝昌台地电阻率整点值、日均值对应的振幅谱。宝昌台地电阻率不仅具有早已熟知的年变特征，还具有日变特征。 （2）宝昌台地电阻率 EW 测道具有较为清晰的日变特征，但是 NS 测道由于噪声较大，日变特征不明显。通过小波方法对宝昌台

21、地电阻率进10行了去噪处理，发现两测道都具有清晰的日变形态，且基本重合。将地电阻率去噪后的结果，与同期的温度进行相关性分析，发现二者之间具有显著相关性。进一步分析认为，该现象可能是由于外线路系统对地绝缘不良造成的。 （3）宝昌台地电阻率具有明显的“夏低冬高”年变特征，这主要是由温度和降雨量等气象因素年变引起的。宝昌台和乌加河台地电阻率年变形态相反，这主要是由两台站地下电性结构差异引起的。 （4）宝昌台地电阻率长期变化存在各向异性特征，其中 NS 测道具有下降变化，而 EW 测道则没有。将宝昌台地电阻率 NS 测道 19932008年月均值数据与该台所在区域温度、降雨量同期月均值数据进行定量分析

22、，发现 NS 测道具有的长期变化与测区的水位、温度和降雨量的长期变化无关，主要是由于受到台站所在区域近 EW 向应力场长期、持续作用引起的。 参考文献： 陈棋福，郑大林，高荣胜 2003 中国震例（19971999）M.北京：地震出版社 戴勇，闫海宾，丁风和，等 2009 张北 62 级地震宝昌地震台地电异常特征及机制分析J.地震地磁观测与研究，30（1）：51-55 杜学彬，阮爱国，范世宏，等 2001 强震近震中区地电阻率变化速率的各向异性J.地震学报，23（3）：289-297 冯志生，杨建军，梅卫萍，等 2004 华东地区地电阻率各向异性度的地震前兆异常特征初步研究J.地震学报，26（2）：223-227