基于褶积滤波的山丹地电阻率归一化速率异常特征分析.pdf
《基于褶积滤波的山丹地电阻率归一化速率异常特征分析.pdf》由会员分享,可在线阅读,更多相关《基于褶积滤波的山丹地电阻率归一化速率异常特征分析.pdf(7页珍藏版)》请在咨信网上搜索。
1、利用褶积滤波法对山丹地电阻率进行去除降雨校正处理,再分别对原始数据和去降雨后数据在进行归一化月速率分析前作不同处理方式组合的结果进行映震效能统计,从而选出归一化速率分析前数据最优处理方式。结果显示,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山丹地电阻率的速率超限
2、异常对逐渐平静的祁连山地震带地震危险性进行预测,分析山丹观测站附近地震活动水平是否真的降低,结果显示,山丹地电阻率观测值仍然存在速率超限异常,说明祁连山地震带的地震危险性仍然存在,地震活动水平并未减弱。关键词:山丹地电阻率;褶积滤波;归一化月速率分析;祁连山地震带中图分类号:P 3 1 5;P 3 1 9 文献标识码:A 地电阻率是目前国内用于地震预报的定点前兆观测重要手段之一。随着地电仪的改进,台站数字化仪器稳定性有所提升,但地表的气象变化、人为施工等环境因素对地电阻率数据造成干扰,尤其是降雨量、地下水等干扰严重,常常在降雨后表现为地电阻率测值下降,且测值下降幅度随降雨量的增加而增大,随着地
3、面降雨逐渐蒸发,地电阻率数值也逐渐恢复原来水平1-3。因此为使地电阻率资料分析结果可靠且严谨,在跟踪分析各台站地电阻率数据时,有必要对地电阻率数据进行去降雨因素分析。去除降雨量对地电阻率影响的方法包括相关分析校正2及褶积滤波法1,4-8,其中褶积滤波法应用效果最好且应用最为广泛。月速率法是描述地电阻率异常的常用因子之一,适用于中短期异常判定,在地电阻率异常分析中应用广泛5,9。本文选取祁连山地震带地电阻率观测场地效果最佳的山丹台数据,利用褶积滤波去除降雨影响,依据地电阻率观测数据和台站附近中强地震震例确定归一化速率分析最优数据处理方式,并对比去除降雨前后归一化速率超限异常的映震效果。2 0 2
4、 2-0 1祁连山地震带中强地震活动进入空前活跃期,8月才逐渐趋于平静。结合此前中强地震震例总结得到的超限异常映震指标,分析8月之后山丹地电阻率的超限异常速率,对祁连山地震带地震活动预测有重要意义,同时为祁连山地震带地震危险性评价提供科学依据。1 山丹观测站基本情况山丹地震台位于山丹县城西8k m处,海拔17 4 8m;台站观测场岩性为砂岩,测区地表为旷野型戈壁沙土堆,覆盖层较深,地表几乎没有植被。研究区发育河西走廊区域性主干大断裂,切割地壳较深,走向呈NWW向,活动性较强(图1)。山丹地震台地电阻率观测采用对称四极法(图2)布极,布设有N S、EW、N 4 5 W共3条测道。其中,N S测道
5、和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 n
6、r 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
7、)考虑到降雨对地电阻率的影响特点,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
8、=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各观测站不同测道观测日值,利用褶积滤波去除短期集中降雨对地电阻率的影响,再分析降雨量对山
9、丹电阻率年变和趋势性变化的响应。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
10、、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-
11、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月雨多集中在夏季。山丹地电阻率E
12、W测道和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山丹观测站气象三要素中降雨量月均值数据,利用褶积滤波计算降雨对山丹地电阻率年变和趋势变化的
13、响应,结果如图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
14、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。采用不同滑动步长分别对去干扰数据进行归一化分析发现,选用步长越大,
15、归一化异常形态越不显著,滑动步长为68时,异常形态基本一致,但随着步长增大,异常超限部分变小,从异常识别和指标提取的角度,选取步长6更为合理。月相对变化速率进行地电阻率前兆异常提取的关键是对地电阻率资料进行去倾、去年变及排除周期变化,因此在分别对不同测道数据采用去年变、排除周期变化和去倾3种不同处理方式单独及组合处理数据后再进行月速率分析,并选择最优分析方法。4 去降雨后地电阻率映震效能分析选取2 0 1 42 0 2 1年山丹地电阻率月均值及祁连山地震带1 0次5级 以 上 地 震 震 例,包 括259 第4 3卷第9期李 娜等:基于褶积滤波的山丹地电阻率归一化速率异常特征分析2 0 1 4
16、年青海乌兰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分别统计了异常结束后半年内有震次数(发震)、异常结束后无震次数(
17、虚报)、震前半年内无异常次数(漏报),结果显示,去除降雨后数据异常映震次数显著增多,异常的发震率增加,虚报率降低,多种组合方式获得的异常对应结果基本一致。表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
18、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去年变19346128
19、2去周期371731461去倾28219001 01去年变+去周期462733461去年变+去倾461460192去倾+去周期461373191去年变+去周期+去倾461373191对比原始数据和去除降雨后数据归一化结果,N S测道和NW测道在去除降雨影响后,经归一化速率分析的超限异常对应地震效果有所下降,EW测道去除降雨后的归一化速率分析超限异常映震效果显著增强。对于原始数据,N S测道仅需去周期处理即可有较好的映震效能,而去除降雨后需同时去除年变和周期才能有较好的映震效能;EW测道去除降雨前后的数据在进行去年变+去周期的相同处理后,剔除降雨的映震效能更好,说明EW测道受降雨影响更为显著;N
20、W测道剔除降雨前后的数据在进行归一化分析时的最优处理方式均为去年变+去周期,但明显在去除降雨后虚报率有所下降,说明对于山丹地电阻率而言,降雨对地电阻率的异常提取产生了一定影响。去降雨之后山丹地电阻率进行归一化月速率分析最优结果为:去除降雨后的山丹地电阻率N S测道归一化处理既可选择去倾+去年变,也可选择去倾+去周期;EW测道和NW测道归一化处理则选择去周期的方式。将山丹地电阻率数据研究时段和震例分析时段均延长至2 0 2 2-0 8,根据选取的归一化前最佳处理方式对山丹地电阻率进行计算,结果如图5和6所示,超限时段和有效预测时间内发震情况统计结果如表2所示。结果发现,N S测道去倾+去周期后共
21、有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
22、 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
23、 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中山丹地
24、电阻率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测道处理方式超限时段异常持续时间/月有效预测时间映震时
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 滤波 山丹 电阻率 归一化 速率 异常 特征 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。