基于褶积滤波的山丹地电阻率归一化速率异常特征分析.pdf

利用褶积滤波法对山丹地电阻率进行去除降雨校正处理,再分别对原始数据和去降雨后数据在进行归一化月速率分析前作不同处理方式组合的结果进行映震效能统计,从而选出归一化速率分析前数据最优处理方式。结果显示,N S测道适用于去倾、去周期处理,EW测道和NW测道均适用于去周期处理,且去除降雨后各测道映震效能均优于去降雨前,说明降雨对山丹地电阻率异常提取有较大干扰。由于2 0 2 2年祁连山地震带表现空前活跃,至8月逐渐恢复平静,因此选择2 0 1 42 0 2 2-0 8祁连山地震带及其附近中强地震活动前山丹地电阻率表现出的速率超限异常总结映震指标,利用2 0 2 2-0 8 1 2山丹地电阻率的速率超限异常对逐渐平静的祁连山地震带地震危险性进行预测,分析山丹观测站附近地震活动水平是否真的降低,结果显示,山丹地电阻率观测值仍然存在速率超限异常,说明祁连山地震带的地震危险性仍然存在,地震活动水平并未减弱。关键词:山丹地电阻率;褶积滤波;归一化月速率分析;祁连山地震带中图分类号:P 3 1 5;P 3 1 9 文献标识码:A 地电阻率是目前国内用于地震预报的定点前兆观测重要手段之一。随着地电仪的改进,台站数字化仪器稳定性有所提升,但地表的气象变化、人为施工等环境因素对地电阻率数据造成干扰,尤其是降雨量、地下水等干扰严重,常常在降雨后表现为地电阻率测值下降,且测值下降幅度随降雨量的增加而增大,随着地面降雨逐渐蒸发,地电阻率数值也逐渐恢复原来水平1-3。因此为使地电阻率资料分析结果可靠且严谨,在跟踪分析各台站地电阻率数据时,有必要对地电阻率数据进行去降雨因素分析。去除降雨量对地电阻率影响的方法包括相关分析校正2及褶积滤波法1,4-8,其中褶积滤波法应用效果最好且应用最为广泛。月速率法是描述地电阻率异常的常用因子之一,适用于中短期异常判定,在地电阻率异常分析中应用广泛5,9。本文选取祁连山地震带地电阻率观测场地效果最佳的山丹台数据,利用褶积滤波去除降雨影响,依据地电阻率观测数据和台站附近中强地震震例确定归一化速率分析最优数据处理方式,并对比去除降雨前后归一化速率超限异常的映震效果。2 0 2 2-0 1祁连山地震带中强地震活动进入空前活跃期,8月才逐渐趋于平静。结合此前中强地震震例总结得到的超限异常映震指标,分析8月之后山丹地电阻率的超限异常速率,对祁连山地震带地震活动预测有重要意义,同时为祁连山地震带地震危险性评价提供科学依据。1 山丹观测站基本情况山丹地震台位于山丹县城西8k m处,海拔17 4 8m;台站观测场岩性为砂岩,测区地表为旷野型戈壁沙土堆,覆盖层较深,地表几乎没有植被。研究区发育河西走廊区域性主干大断裂,切割地壳较深,走向呈NWW向,活动性较强(图1)。山丹地震台地电阻率观测采用对称四极法(图2)布极,布设有N S、EW、N 4 5 W共3条测道。其中,N S测道和N 4 5 W测道供电极极距A B=1.0 0 0k m,测量极极距MN=0.2 0 0k m,装置系数K=(AMAN)MN=3.7 7 0k m;EW测道供电极极距A B=1.0 0 0k m,测量极极距MN=0.2 2 6k m,装置系数K=(AMAN)MN=3.2 9 0k m。2 褶积滤波去降雨分析假设B是一个系统,响应函数为H(t),t时 第4 3卷第9期李 娜等:基于褶积滤波的山丹地电阻率归一化速率异常特征分析图1 山丹地电阻率观测站附近构造特征F i g.1 T e n t o n i cc h a r a c t e r i s t i c so fS h a n d a nr e s i s t i v i t ys t a t i o n图2 山丹地电阻率布极F i g.2 M e a s u r e m e n t e l e c t r o d ed i s t r i b u t i o no fS h a n d a nr e s i s t i v i t ys t a t i o n刻输入X(t),则输出Y(t)=H(t)X(t),二者的褶积为t=0H()X(t-)。本文将土层看作一个系统,降雨量CR(t)作为输入信号,降雨量对地电阻率s的影响系统函数为R(t),则降雨量对s的影响可表示为:s(t)=R(t)CR(t)=t=0R()CR(t-)(1)考虑到降雨对地电阻率的影响特点,R(t)为褶积算法的系统函数且满足式(2)特征,其中由于降雨下渗过程体现降雨的滞后效应,根据降雨渗透过程对地电阻率影响的研究结果1,6-7,选择M=1d,表示降雨即时影响时间为1d,N=6m表示降雨的滞后影响6个月:R(t)=0,t0A0+A1t,0tMA2+A3t+A4t-1+A5t-2,MN (2)联合式(1)和(2)可得s(t)=K X,其中,X=A0,A1,A2,A3,A4,A5TK=M=0CR(t1-)M=0 CR(t1-)M=0CR(t2-)M=0 CR(t2-)N=MCR(t1-)N=M CR(t1-)N=MCR(t2-)N=M CR(t2-)N=M-1CR(t1-)N=M-2CR(t1-)N=M-1CR(t2-)N=M-2CR(t2-)M=0CR(tn-)M=0 CR(tn-)N=MCR(tn-)N=M CR(tn-)N=M-1CR(tn-)N=M-2CR(tn-)选取多个时段的降雨数据和地电阻率整点值或日均值观测数据,给出s(t)的变化,得到降雨对地电阻率的影响特征曲线,经残差计算即可获得降雨对地电阻率的短期变化影响。再选用月均值通过褶积滤波去除降雨影响,便可分析降雨对地电阻率年变趋势的影响。选取观测站每年降雨量较为集中时段及集中降雨开始前1d各观测站不同测道观测日值,利用褶积滤波去除短期集中降雨对地电阻率的影响,再分析降雨量对山丹电阻率年变和趋势性变化的响应。2.1 降雨对地电阻率的短期变化影响选取2 0 1 42 0 2 2-0 8山丹观测站同台观测的气象 三要素观测 资料中降雨 量日均值数 据及2 0 1 4-0 6-1 10 7-2 2、2 0 1 4-0 8-0 20 9-1 1、2 0 1 5-0 5-2 70 7-2 2、2 0 1 5-0 9-0 22 1、2 0 1 6-0 4-1 50 5-0 1、2 0 1 6-0 5-1 10 6-0 3、2 0 1 6-0 6-2 00 7-2 1、2 0 1 6-0 8-1 02 9、2 0 1 7-0 4-0 72 4、2 0 1 7-0 5-0 32 1、2 0 1 7-0 6-0 30 7-3 1、2 0 1 7-0 8-0 70 9-0 4、2 0 1 7-1 0-1 73 1、2 0 1 8-0 6-1 60 7-2 6、2 0 1 8-0 8-0 40 9-1 7、2 0 1 8-1 0-2 63 0、2 0 1 9-0 4-0 80 5-0 9、2 0 1 9-0 6-0 90 7-2 3、2 0 1 9-0 8-0 50 9-1 9、2 0 2 0-0 5-0 50 6-2 2、2 0 2 0-0 7-1 70 9-2 1、2 0 2 1-0 4-1 30 5-2 2、2 0 2 1-0 6-0 60 7-2 1、2 0 2 1-0 8-0 20 9-1 4、2 0 2 2-0 5-2 80 6-2 0、2 0 2 2-0 7-1 30 8-2 8共2 6个时段的降雨数据和降雨时段前1dN S测道、EW测道及NW测道地电阻率数据,分别给出系数矩阵K和降雨影响地电阻率变化量s(t),根据s(t)=K X得出参数矩阵X,再将X作用于山丹地电阻率历年降雨观测数据日均值,给出降雨对地电阻率观测影响量s(t)的短期变化,结果如图3所示。图中,蓝色代表原始曲线,红色代表通过褶积滤波法去除降雨对各观测站地电阻率影响后的残差变化,可以看出,降159大 地 测 量 与 地 球 动 力 学2 0 2 3年9月雨多集中在夏季。山丹地电阻率EW测道和NW测道褶积拟合基本与实际观测曲线重合,N S测道拟合后残差数值高于实际观测值,说明短期降雨可能造成地电阻率降低,降雨结束后测值恢复。图3 日降雨量对山丹地电阻率短期变化影响F i g.3 E f f e c to fd a i l yr a i n f a l l o ns h o r t-t e r mv a r i a t i o no f t h eS h a n d a nr e s i s t i v i t y2.2 降雨对地电阻率的季节性影响选取2 0 1 42 0 2 2-0 8山丹观测站气象三要素中降雨量月均值数据,利用褶积滤波计算降雨对山丹地电阻率年变和趋势变化的响应,结果如图4所示。由图可知,降雨对电阻率季节性的最大影响出现在每年最高降雨之后,山丹地电阻率在剔除月降雨量对低电阻率的影响值后,原始观测曲线“冬高夏低”的年变规律均被削弱,但不改变长趋势变化形态,因此地电阻率去除降雨季节性周期影响之后变化曲线较为平稳,年变形态相比原始观测曲线不再清晰,趋势变化依旧存在,说明降雨量对地电阻率的年变形态有一定的影响,但对趋势性变化影响不大。图4 月降雨量对山丹地电阻率年变形态的影响F i g.4 E f f e c to fm o n t h l yr a i n f a l l o na n n u a l v a r i a t i o no f t h eS h a n d a nr e s i s t i v i t y3 归一化月速率分析月速率分析的原理为:设定各个观测站地电阻率观测序列=1,2,N,以滑动步长计算归一化月速率的序列为Si=RiKin-1(i=n,n+1,N),其中Ri为序列曲线斜率,Ki为自相关系数,分母为分子乘积关于时间序列的均方差,N为选取资料长度,滑动步长n视情况选定。若最终处理结果Si绝对值超过2.4,则视为异常5,1 0-1 3。采用G S E P软件1 4进行归一化月速率分析时,数据选用月均值,推荐滑动步长为41 2,优势滑动步长为681 5。采用不同滑动步长分别对去干扰数据进行归一化分析发现,选用步长越大,归一化异常形态越不显著,滑动步长为68时,异常形态基本一致,但随着步长增大,异常超限部分变小,从异常识别和指标提取的角度,选取步长6更为合理。月相对变化速率进行地电阻率前兆异常提取的关键是对地电阻率资料进行去倾、去年变及排除周期变化,因此在分别对不同测道数据采用去年变、排除周期变化和去倾3种不同处理方式单独及组合处理数据后再进行月速率分析,并选择最优分析方法。4 去降雨后地电阻率映震效能分析选取2 0 1 42 0 2 1年山丹地电阻率月均值及祁连山地震带1 0次5级 以 上 地 震 震 例,包 括259 第4 3卷第9期李 娜等:基于褶积滤波的山丹地电阻率归一化速率异常特征分析2 0 1 4年青海乌兰5.1级、2 0 1 5年青海祁连5.2级、2 0 1 6年青海门源6.4级、2 0 1 7年内蒙古阿拉善左旗5级、2 0 1 9年青海茫崖5级、甘州5级、2 0 2 1年青海茫崖5.8级、阿克塞5.5级及2 0 1 7年九寨沟7级、2 0 2 1年青海玛多7.4级地震,对山丹地电阻率原始数据和去除降雨后数据分别进行归一化速率计算。月速率分析方法是异常分辨率较高的中短期预测方法1 4-1 5,超限异常的预测有效时间自超限异常开始第2个月起,到异常结束后6个月止,根据异常预测时段和发震时间分析地震与超限异常的对应情况。表1分别统计了异常结束后半年内有震次数(发震)、异常结束后无震次数(虚报)、震前半年内无异常次数(漏报),结果显示,去除降雨后数据异常映震次数显著增多,异常的发震率增加,虚报率降低,多种组合方式获得的异常对应结果基本一致。表1 山丹地电阻率原始数据与去降雨数据归一化处理方式映震情况对比T a b.1 T h ec o m p a r i s o no f s e i s m i c r e s p o n s eo b t a i n e db yd i f f e r e n t t r e a t m e n tm e t h o d so fn o r m a l i z e da n a l y s i s f o rS h a n d a no r i g i n a l r e s i s t i v i t yd a t aa n dt h ed a t aa f t e r r e m o v i n gr a i n f a l l数据类型归一化前处理方式N S测道/次EW测道/次NW测道/次发震虚报漏报发震虚报漏报发震虚报漏报原始数据不作处理01 00191191去年变282462462去周期550553553去倾01 0001 0201 01去年变+去周期460552552去年变+去倾371461460去倾+去周期551463370去年变+去周期+去倾551370281去降雨数据不作处理282280191去年变193461282去周期371731461去倾28219001 01去年变+去周期462733461去年变+去倾461460192去倾+去周期461373191去年变+去周期+去倾461373191对比原始数据和去除降雨后数据归一化结果,N S测道和NW测道在去除降雨影响后,经归一化速率分析的超限异常对应地震效果有所下降,EW测道去除降雨后的归一化速率分析超限异常映震效果显著增强。对于原始数据,N S测道仅需去周期处理即可有较好的映震效能,而去除降雨后需同时去除年变和周期才能有较好的映震效能;EW测道去除降雨前后的数据在进行去年变+去周期的相同处理后,剔除降雨的映震效能更好,说明EW测道受降雨影响更为显著;NW测道剔除降雨前后的数据在进行归一化分析时的最优处理方式均为去年变+去周期,但明显在去除降雨后虚报率有所下降,说明对于山丹地电阻率而言,降雨对地电阻率的异常提取产生了一定影响。去降雨之后山丹地电阻率进行归一化月速率分析最优结果为:去除降雨后的山丹地电阻率N S测道归一化处理既可选择去倾+去年变,也可选择去倾+去周期;EW测道和NW测道归一化处理则选择去周期的方式。将山丹地电阻率数据研究时段和震例分析时段均延长至2 0 2 2-0 8,根据选取的归一化前最佳处理方式对山丹地电阻率进行计算,结果如图5和6所示,超限时段和有效预测时间内发震情况统计结果如表2所示。结果发现,N S测道去倾+去周期后共有4次超限异常,在有效预测期内对应3次地震,去倾+去年变后共有5次超限异常,有效预测期内对应地震1次;EW测道去降雨+去周期后共有8次超限异常,有效预测期内发震1 2次,其中2 0 2 0年年底出现的超限异常时段内共发生3次地震,2 0 2 1年底超限异常时段内共发生5次地震;NW测道去降雨+去周期后共出现超限异常3次,发生4次地震,其中2 0 2 1年年初超限 异 常 有 效 预 测 期 内 发 生3次 地 震。除2 0 1 42 0 2 1年1 0次地震,2 0 2 2年门源6.9级地震前,EW测道也出现超限异常。图5 山丹去降雨后N S测道归一化结果F i g.5 N o r m a l i z e dr e s u l t sa f t e r r e m o v i n gr a i n f a l lo fS h a n d a nN Sm e a s u r i n gw a y山丹地电阻率3个测道中映震效果最好的是EW测道,N S测道和NW测道映震效果相当,测359大 地 测 量 与 地 球 动 力 学2 0 2 3年9月图6 山丹去降雨、去周期后归一化结果F i g.6 N o r m a l i z e dr e s u l t sa f t e r r e m o v i n gr a i n f a l la n dr e m o v i n gt h ec y c l eo fS h a n d a ns t a t i o n道间异常变化特征并不相同,表现出各向异性,与地电阻率观测场地的地层岩性及层状岩石含水量影响有关1 6。因此归一化速率分析结果显示,并不是所有地震前3个测道均能同步出现异常特征,且由于各测道归一化前数据处理方式不同,映震效果差异也较大。地电阻率各向异性变化存在一个高值易发震变化区间,对台站周围5 0 0k m范围内5级以上 中强地震具 有明显的异 常特征1 7,而满足3 0 0k m范围内5.0MS5.9、4 0 0k m范围内6.0MS6.9及5 0 0k m内MS7.0的地震视为与异常对应1 8,由此看来,表2中山丹地电阻率N S测道在茫崖5.0级地震前、EW测道在茫崖5.0级、5.8级及九寨沟7.0级地震前出现的异常属于虚报,并非前兆异常。去除超出5 0 0k m范围的地震后,山丹地电阻表2 去降雨后山丹地电阻率超限异常映震情况T a b.2 A b n o r m a l s e i s m i c r e s p o n s eo f r e s i s t i v i t yo v e r l i m i t a t S h a n d a ns t a t i o na f t e rr e m o v i n gr a i n f a l l测道处理方式超限时段异常持续时间/月有效预测时间映震时间地震震中距/k m是否对应N S去倾+去周期去倾+去年变2 0 1 9-0 212 0 1 9-0 30 82 0 1 9-0 3-2 8茫崖5.0级8 8 3.8 0超范围,虚报2 0 1 9-0 712 0 1 9-0 82 0 2 0-0 12 0 1 9-0 9-1 6甘州5.0级6 7.0 6是2 0 1 9-1 01 122 0 1 9-1 12 0 2 0-0 5-虚报2 0 2 0-1 01 232 0 2 0-1 12 0 2 1-0 52 0 2 1-0 5-2 2玛多7.4级4 9 5.4 6是2 0 1 8-0 80 922 0 1 8-0 92 0 1 9-0 32 0 1 9-0 3-2 8茫崖5.0级8 8 3.8 0超范围,虚报2 0 1 8-1 22 0 1 9-0 232 0 1 9-0 10 8-虚报2 0 1 9-1 012 0 1 9-1 02 0 2 0-0 4-虚报2 0 2 0-0 612 0 2 0-0 61 2-虚报2 0 2 0-0 91 022 0 2 0-1 02 0 2 1-0 4-虚报EW去周期2 0 1 5-1 22 0 1 6-0 452 0 1 6-0 11 02 0 1 6-0 91 242 0 1 6-1 02 0 1 7-0 62 0 1 7-0 71 152 0 1 7-0 82 0 1 8-0 52 0 1 8-0 912 0 1 8-1 02 0 1 9-0 32 0 1 9-0 10 442 0 1 9-0 21 02 0 2 0-1 12 0 2 1-0 462 0 2 0-1 22 0 2 1-1 02 0 2 1-1 12 0 2 2-0 342 0 2 1-1 22 0 2 2-0 92 0 2 2-0 812 0 2 2-0 92 0 2 3-0 22 0 1 6-0 1-2 1门源6.4级1 3 2.8 0是-虚报2 0 1 8-0 8-0 8九寨沟7.0级6 6 8.5 8超范围,虚报2 0 1 9-0 3-2 8茫崖5.0级8 8 3.8 0超范围,虚报2 0 1 9-0 9-1 6甘州5.0级6 7.0 6是2 0 2 1-0 5-2 2玛多7.4级4 9 5.4 6是2 0 2 1-0 6-1 6茫崖5.8级6 3 2.9 5超范围,虚报2 0 2 1-0 8-2 6阿克塞5.5级4 7 9.6 0是2 0 2 2-0 1-0 8门源6.9级1 1 5.4 6是2 0 2 2-0 1-2 3德令哈5.8级3 2 0.6 9是2 0 2 2-0 3-1 5肃南5.1级2 9 3.3 8是2 0 2 2-0 3-2 7德令哈6.0级3 2 3.3 1是2 0 2 2-0 4-1 5德令哈5.4级3 2 3.0 6是待对应NW去周期2 0 1 6-0 212 0 1 6-0 30 82 0 1 8-0 70 932 0 1 8-0 82 0 1 9-0 32 0 2 1-0 212 0 2 1-0 30 9-虚报2 0 1 9-0 3-2 8茫崖5.0级8 8 3.8 0超范围,虚报2 0 2 1-0 5-2 2玛多7.4级4 9 5.4 6是2 0 2 1-0 6-1 6茫崖5.8级6 3 2.9 5超范围,虚报2 0 2 1-0 8-2 6阿克塞5.5级4 7 9.6 0是率N S测道经去倾+去周期处理后映震效果较去倾+去年变好。山丹地电阻率N S测道和EW测道均在甘州5.0级、玛多7.4级地震前出现超限异常且与其对应,NW测道也在玛多7.4级地震前出现异常且对应,EW测道和NW测道则在阿克塞5.5级地震前同时出现超限异常且与地震对应。EW测道超限异常持续时间为46个月,且对应多次中强地震。山丹地电阻率超限异常对祁连山地震带的预测指标如表3所示,其中N S测道速率超限异常持续12个月时,在异常开始至异常结束后半年内存在发生5级左右地震的可能;EW测道速率超限异常持续46个月时,异459 第4 3卷第9期李 娜等:基于褶积滤波的山丹地电阻率归一化速率异常特征分析常开始至结束后半年内祁连山地震带存在发生6级左右地震的可能;NW测道速率超限异常持续1个月时,异常开始至结束后半年内祁连山地震带存在发生56级地震的可能。表3 去降雨后山丹地电阻率超限异常预测指标T a b.3 P r e d i c t i o n i n d e xo f r e s i s t i v i t yo v e r l i m i t a n o m a l ya tS h a n d a ns t a t i o na f t e r r e m o v i n gr a i n f a l l测道异常形态异常超限时间/月预测震级震例N S速率超限125级左右甘州5.0级EW速率超限466级左右甘州5.0级、门源6.4级、门源6.9级、阿克塞5.5级、玛多7.4级、德令哈5.8级、德令哈6.0级、德令哈5.4级、肃南5.1级NW速率超限156级阿克塞5.5级37级左右玛多7.4级5 结 语通过对2 0 1 42 0 2 2-0 8山丹地电阻率观测数据进行褶积滤波去除降雨影响后,利用归一化月速率分析结果的映震效果显著好于原始数据映震效果,根据去除降雨后的电阻率数据映震指标可知,山丹地电阻率3个测道的结果均可作为祁连山地震带中强地震前的短期异常指标,根据震例情况,对台站周边3 0 0k m范围内56级地震的指示意义更强。由2 0 2 2-0 81 2山丹地电阻率各测道速率超限异常可知,EW测道2 0 2 2-0 8以来持续5个月出现超限异常,NW测道2 0 2 2-1 2出现超限异常,说明2 0 2 3年祁连山地震带中西段发生中强地震的危险性仍未得到缓解,需重点关注速率超限异常。后期应结合震例指标情况,根据速率超限异常结束时间来判定未来发震的优势时间,同时密切关注研究区地震活动性及应力演化特征,若出现地震活动显著增强及应力积累特征,则可进一步验证山丹地电阻率指标的可靠性及可行性。致谢:感谢河北省地震局张国苓高级工程师和陕西省地震局石富强副研究员为本文褶积滤波降雨分析给予指导和帮助。参考文献1 张学民,王志贤,臧明珍,等.降雨对地电阻率干扰的分析J.华北地震科学,1 9 9 6,1 4(4):7 1-7 5(Z h a n gX u e m i n,W a n gZ h i x i a n,Z a n gM i n g z h e n,e t a l.A n a l y s i so f I n t e r f e r-e n c eo f R a i n f a l lo n E a r t h R e s i s t i v i t yJ.N o r t h C h i n aE a r t h q u a k eS c i e n c e s,1 9 9 6,1 4(4):7 1-7 5)2 张伟峰,王远,潘洁,等.嘉山台地电阻率与降雨关系探究J.华南地震,2 0 1 5,3 5(1):8 4-8 9(Z h a n gW e i f e n g,W a n gY u a n,P a nJ i e,e ta l.E x p l o r a t i o no nt h eR e l a t i o n s h i po ft h eR a i n f a l l a n dE a r t hR e s i s t i v i t yi nJ i a s h a nS e i s m i cS t a t i o nJ.S o u t hC h i n a J o u r n a l o fS e i s m o l o g y,2 0 1 5,3 5(1):8 4-8 9)3 刘允秀,陈华静,程瑞年,等.地电阻率与地下水位、大气降水关系研究J.中国地震,1 9 9 9,1 5(2):1 8 4-1 8 9(L i uY u n x i u,C h e nH u a j i n g,C h e n gR u i n i a n,e ta l.T h eS t u d yo fR e l a t i o n s h i pb e t w e e nE a r t hR e s i s t i v i t ya n dG r o u n d w a t e rL e v e l,P r e c i pI t a t i o nJ.E a r t h q u a k eR e s e a r c hi nC h i n a,1 9 9 9,1 5(2):1 8 4-1 8 9)4 宋晓磊,黄清华.应用褶积滤波和多元回归分析降雨对地电阻率的干扰C.中国地球物理学会第2 1届年会,长春,2 0 0 5(S o n gX i a o l e i,H u a n gQ i n g h u a.A p p l i c a t i o no fC o n v o l u-t i o nF i l t e r i n ga n dM u l t i p l eR e g r e s s i o nA n a l y s i s o n t h e I n t e r f e r-e n c eo fR a i n f a l lo nG r o u n dR e s i s t i v i t yC.T h e2 1 s tA n n u a lM e e t i n go fC h i n e s eG e o p h y s i c s,C h a n g c h u n,2 0 0 5)5 李姜,乔子云,张国苓,等.华北地区地电阻率归一化速率异常分析J.地震地磁观测与研究,2 0 1 9,4 0(6):6 2-7 1(L i J i a n g,Q i a oZ i y u n,Z h a n gG u o l i n g,e t a l.A n o m a l i e so fN o r m a l i z e dR a t eo fG e o e l e c t r i cR e s i s t i v i t yi nN o r t hC h i n aJ.S e i s m o l o g i c a la n dG e o m a g n e t i cO b s e r v a t i o na n dR e-s e a r c h,2 0 1 9,4 0(6):6 2-7 1)6 张国苓,茅远哲,乔子云,等.降雨对昌黎台地电阻率的干扰分析J.地震,2 0 1 9,3 9(4):1 1 0-1 1 7(Z h a n gG u o l i n g,M a oY u a n z h e,Q i a oZ i y u n,e t a l.Q u a n t i t a t i v eA n a l y s i so fR a i n f a l lE f f e c t so nE a r t hR e s i s t i v i t y i nC h a n g l i S e i s m i cS t a-t i o nJ.E a r t h q u a k e,2 0 1 9,3 9(4):1 1 0-1 1 7)7 石富强,张国强,方炜,等.陕西周至地电台地电阻率年变特征分析J.地震学报,2 0 1 4,3 6(6):11 1 3-11 2 3(S h iF u q i a n g,Z h a n gG u o q i a n g,F a n gW e i,e t a l.A n n u a lV a r i a-t i o nC h a r a c t e r i s t i c so fG e o r e s i s t i v i t ya tZ h o u z h iG e o e l e c t r i cS t a t i o n,S h a a n x iP r o v i n c eJ.A c t aS e i s m o l o g i c aS i n i c a,2 0 1 4,3 6(6):11 1 3-11 2 3)8 番绍辉,缪素秋,黎朕灵,等.降雨对腾冲地震台地电阻率的干扰分析J.高原地震,2 0 2 1,3 3(3):3 3-4 0(F a nS h a o-h u i,M i a oS u q i u,L iZ h e n l i n g,e ta l.A n a l y s i so f I n f l u e n c eo fR a i n f a l l o nG e o e l e c t r i cR e s i s t i v i t yo fT e n g c h o n gS e i s m i c S t a t i o nJ.P l a t e a uE a r t h q u a k eR e s e a r c h,2 0 2 1,3 3(3):3 3-4 0)9 张彩艳,闫勋,薛汶汶,等.基于归一化月速率法的甘肃四四地震台地电阻率异常分析J.山西地震,2 0 2 2(1):3 6-3 9(Z h a n gC a i y a n,Y a nX u n,X u eW e n w e n,e t a l.A n a l-y s i so fE a r t h R e s i s t i v i t y A n o m a l ya tG a n s u4 0 4S e i s m i cS t a t i o nB a s e do n N o r m a l i z e d M o n t h l yR a t e M e t h o dJ.E a r t h q u a k eR e s e a r c h i nS h a n x i,2 0 2 2(1):3 6-3 9)1 0杜学彬,阮爱国,范世宏,等.强震近震中区地电阻率变化速率的各向异性J.地震学报,2 0 0 1,2 3(3):2 8 9-2 9 7(D uX u e b i n,R u a nA i g u o,F a nS h i h o n g,e ta l.A n i s o t r o p yo ft h eV a r i a t i o nR a t eo fA p p a r e n tR e s i s t i v i t yn e a r t h eE p i c e n-t r a lR e g i o no fS t r o n gE a r t h q u a k e sJ.A c t aS e i s m o l o g i c aS i n i c a,2 0 0 1,2 3(3):2 8 9-2 9 7)1 1杜学彬.在地震预报中的两类视电阻率变化J.中国科学:地球科学,2 0 1 0,4 0(1 0):13 2 1-13 3 0(D uX u e b i n.T w oT y p e so fC h a n g e si nA p p a r e n tR e s i s t i v i t yi nE a r t h-559大 地 测 量 与 地 球 动 力 学2 0 2 3年9月q u a k eP r e d i c t i o nJ.S c i e n c eC h i n a:E a r t hS c i e n c e,2 0 1 0,4 0(1 0):13 2 1-13 3 0)1 2杜学彬,孙君嵩,陈军营.地震预测中的地电阻率数据处理方法J.地震学报,2 0 1 7,3 9(4):5 3 1-5 4 8(D uX u e b i n,S u nJ u n s o n g,C h e nJ u n y i n g.P r o c e s s i n g M e t h o d sf o rt h eO b s e r v a t i o nD a t ao fG e o r e s i s t i v i t y i nE a r t h q u a k eP r e d i c t i o nJ.A c t aS e i s m o l o g i c aS i n i c a,2 0 1 7,3 9(4):5 3 1-5 4 8)1 3李相平,张建国,李晓鹏,等.强震前视电阻率月速率变化特征分析 J.华 北 地 震 科 学,2 0 1 6,3 4(3):6 8-7 2(L iX i a n g p i n g,Z h a n gJ i a n g u o,L iX i a o p e n g,e ta l.R e s e a r c ho nN o r m a l i z e dV a r i a t i o nC h a r a c t e r i s t i c so fA p p a r e n tR e s i s-t i v i t yb e f o r eS t r o n gE a r t h q u a k e sJ.N o r t hC h i n aE a r t h-q u a k eS c i e n c e s,2 0 1 6,3 4(3):6 8-7 2)1 4郝臻,杜学彬,王静波.地电学地震预报方法软件系统(G S E P)J.西北地震学报,2 0 0 0,2 2(2):1 5 4-1 5 9(H a oZ h e n,D uX u e b i n,W a n gJ i n g b o.G e o e l e c t r i c i t yS o f t w a r ef o rE a r t h q u a k eP r e d i c t i o nJ.N o r t h w e