兴济地震台大地电场方位变化特征.pdf
《兴济地震台大地电场方位变化特征.pdf》由会员分享,可在线阅读,更多相关《兴济地震台大地电场方位变化特征.pdf(6页珍藏版)》请在咨信网上搜索。
1、李瑞卿,郭学增,李凤,等.兴济地震台大地电场方位变化特征 J.华北地震科学,2023,41(3):91-96.doi:10.3969/j.issn.10031375.2023.03.013.LI Ruiqing,GUO Xuezeng,LI Feng,et al.Azimuth Variation Characteristics of Geoelectric Field at Xingji Seismic StationJ.North ChinaEarthquake Sciences,2023,41(3):91-96.doi:10.3969/j.issn.10031375.2023.03.01
2、3.兴济地震台大地电场方位变化特征李瑞卿,郭学增,李 凤,张 蕾,信世民,尹宏伟(石家庄地震监测中心站,石家庄 050021)摘要:根据地电场成分分类取空间电磁环境较为平静且地电场观测数据日变形态较为稳定的时段对兴济台大地电场 ET进行提取,然后利用正交向两测道大地电场合成总的大地电场 ET 总。对总大地电场方位变化进行对比发现,兴济台午前午后大地电场方位角图像较平凉、大柏舍数据离散度高,无午前午后 3 个阶段的稳定形态变化特征;通过谐波分析及 F 检验发现,兴济台大地电场前 10 阶谐波拟合程度相对较低、高阶谐波成分显著,说明兴济台地电场成分复杂,可能是造成一天中大地电场方位角数据离散的原因
3、。关键词:大地电场方位角;谐波分析;拟合程度;F 检验中图分类号:P315.722 文献标志码:A 文章编号:10031375(2023)03009106doi:10.3969/j.issn.10031375.2023.03.013 0 引言地电场是重要的地球物理场之一,其不仅能够反映空间电磁场的变化,还能反映与之有关的地壳和地表的电性结构及其变化。20 世纪中叶 Chapman提出,地电场由大地电场和自然电场组成1。地电场观测数据是观测方法、技术、场地环境、电磁环境等的综合体现2。随着对地电场的深入研究,有学者将自然电场 Esp与大地电场 ET进行提取及分析研究3-4,研究表明地电场波形日变
4、化可能源于大地电场,其起源主要有空间电流系以及潮汐作用5-8;有学者对地电场稳定的峰-谷日变波形的研究探索出大地电场日变化的潮汐机理说5-6,认为地电场日变波形的物理过程是固体潮或空间电流导致裂隙水周期性渗流产生的过滤电场进而影响了地电场观测波形变化。裂隙水的渗透流动与地下裂隙结构密切相关,通过大地电场计算的优势方位角在一定程度上反映了地下裂隙的优势方位5-7。目前,许多震后分析研究表明优势方位角在震前存在不同形态的异常变化9-12。另外,基于谐波拟合计算的优势方位角对地电暴、雷电、高压直流输电等地电场干扰信号具有一定的抗干扰能力13-14。因此,大地电场优势方位角的研究对地震预测有一定的意义
5、。地电场是矢量场,对于地电场的矢量变化研究发现,地电场在很多情况下是线性极化的,在地震发生前后地电场的极化方位会发生变化15-17。根据地电场极化的线性特征探究出极化斜率法18、垂直极化投影法19,并对地电场数据极化现象进行分析。大地电场是地电场主要构成之一,谭大诚等7研究表明,源于潮汐作用和空间电流系的大地电场其强度和方向受岩石裂隙结构影响,对大地电场矢量特性进行分析发现,平凉台地电场属于午前午后有稳定正弦波形的 TGF-B 型,其大地电场方位角午前午后存在 3 个阶段的稳定变化7。在对特殊布极方位的兴济台地电场研究过程中发现,同属于 TGF-B 型的兴济台大地电场方位角在午前午后数据离散,
6、方位角不稳定。因此,需对这种不同的大地电场方位角变化形态研究分析,进一步了解兴济台地电场观测数据变化特征,为其他台站地电场观测数据分析及应用提供参考。1 大地电场 ET的提取及大地电场方位角计算地电场是较为稳定的自然电场 Esp和变化电场的叠加,其中变化的电场包括大地电场 ET和干扰收稿日期:2022-12-26基金项目:中国地震局监测、预报、科研三结合课题(3JH-202201056)第一作者简介:李瑞卿(1990),男,工程师,主要从事地电观测与研究工作E-mail:第 41 卷第 3 期华 北 地 震 科 学Vol.41,No.32023 年 7 月NORTH CHINA EARTHQU
7、AKE SCIENCESJul.,2023成分 Er(式 1)。当空间电磁环境较为平静、且地电场观测数据日变形态较为稳定时,视为干扰成分Er较少,那么地电场近似于自然电场 Esp和大地电场 ET的叠加。对地电场 E 进行 FFT 变换,前 10 阶谐波周期顺序是 2324 h、12 h、7.9 h、6 h、4.8 h、4 h、3.4 h、3 h、2.7 h、2.4 h3,将一整天 1 440 个分钟值数据求平均值(即日均值)可以消除大部分周期变化的大地电场 ET成分,进而得到稳定的自然电场 Esp(式 2)。使用地电场分钟值数据减去自然电场 Esp即可得到大地电场 ET分钟值数据(式 3)2-
8、3。E=Esp+ET+Er(1)Esp=1439i=0Ei1440(2)ETi=EiEsp(3)提取各测向大地电场成分后可利用三角函数求解分钟值的总大地电场方位角。根据地电场强度分量值的方向约定2(以电流向北、向东为正,即地电场矢量“南北分量”的北指向南为正,“东西分量”的东指向西为正,“北东/北西分量”的偏北指向偏南为正),可以构建坐标系:以电场中心点为坐标系原点,指向正东方向为 x 轴正向,指向正北方向为 y 轴正向。根据大地电场值的正负判断大地电场方位,进而使用正交向大地电场及三角函数求解总大地电场方位角,此方位角 为 y 轴顺时针偏转的角度(图 1)。则一天中总大地电场 ET 方位角
9、计算公式为:i=180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第一象限)360180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第二象限)180+180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第三象限)180180arcatan(|ETEWi|ETNSi|)(ETEWi0,ETNSi0,坐标系第四象限)(4)由于兴济台地电场布极方式较为特殊,布极方式为沿平行于沧东断裂带走向的 N30E 布极和垂直沧东断裂带走向的 N60W 布极(图 2)。对于布极方式比较特殊的兴济台,可以将N6
10、0W 及 N30E 两测向大地电场在东西、南北方向进行分解,求得东西、南北方向的大地电场,即:ET南北=ETN30Ecos30+ETN60Wcos60ET东西=ETN30Esin30+ETN60Wsin60(5)然后,再根据式 4 进行求解兴济台总大地电场方位角(式 5 与直接使用正交向 N30E 和 N60W 三角函数合成得到的总大地电场位置一致,此处为统一 y 轴指向正北的坐标系进行了转换)。2 大地电场 ET方位角对比谭大诚等7研究表明,平凉台大地电场方位角存在午前午后 3 个稳定的变化阶段,其相对稳定时刻的 3 个阶段为:()6:0010:00 大地电场指向第三象限;()其后大地电场方
11、向发生反转,在 10:0016:00 大地电场指向第一象限;()在 16:3019:30大地电场再次指向第三象限(图 3)。变化的大地电场潮汐机理认为:在固体潮汐作用或 Sq电流电磁感应下,裂隙水中电荷以周期性沿裂隙往返移动,进而造成大地电场 ET在不同阶段存在不同方位,此方位沿裂隙优势方向呈近 180反转。对兴济、平凉、大柏舍、昌黎 2021 年 10 月 28 日数据进行总大地电场方位角求解。兴济台总大地电场方位角散点图离散度高,而且没有午前午后3 个稳定阶段(图 4a)。平凉台总大地电场午前午后方位角数据稳定,其相对稳定时刻的 3 个阶段分别为:()7:3011:00 大地电场指向第三象
12、限,方位角为 247左右;()11:3015:00 大地电场指向第 A1B1A2B2ONyET总x 图 1 大地电场方位角示意图Fig.1 Azimuth schematic diagram of the geoelectric field92华 北 地 震 科 学41 卷 一象限,方位角为 67左右;()16:0018:30大地电场再次指向第三象限,方位角为 247左右(图 4b)。大柏舍台大地电场方位角与平凉台类似经历的 3 个阶段角度分别集中在 28、197、14左右(图 4c),昌黎台与兴济台类似一天中的大地电场方位角离散度高,无法辨别其稳定的形态(图 4d)。3 谐波分析拟合程度及
13、F 检验为了分析提取的大地电场成分,引入谐波分析。yt=y+n/2k=1Ckcos(2ktnk)=y+n/2k=1Akcos(2ktn)+Bksin(2ktn)(6)k=arcatan(BkAk),Ak0arcatan(BkAk),Ak02;Ak=0(7)AkBk其中,、由(8)式计算(k 表示谐波的阶数,同下)即Ak=2nnt=1ytcos(2ktn),Bk=2nnt=1ytsin(2ktn)(8)则各阶谐波周期及振幅为Tk=nk,Ck=A2k+B2k(9)yn/2k=1Ckcos(2ktnk)式 6 中:项为稳定的自然电场 Esp,即地电场日均值;项表示变化的电场:大地电场 ET及干扰电场
14、 Er。116.0E116.5E117.0E117.5E118.0E39.0N39.0N39.0N39.0N图例兴济台直辖市断裂NB1B2A2A1O1、O2A1O1=300 mB1O1=300 mA1B1=424 mA2O2=150 mB2O2=150 mA2B2=212 mNa 构造地质 b 测区布极 图 2 兴济台构造地质和测区布极Fig.2 Geological Structure of Xingji Seismic Station and distribution of surveying area A1B1ONxyET(+)ET()第一象限 图 3 平凉台大地电场方位变化示意图Fig
15、.3 Schematic diagram of azimuth variation of the geoelectric fieldon Pingliang Station3 期李瑞卿,等:兴济地震台大地电场方位变化特征93 理论上,使用谐波对地电场数据进行拟合,拟合使用的谐波阶数越多,拟合结果约接近真实的地电场观测变化。根据潮汐电场理论,地电场日变波形主体基本由前 10 阶谐波成分构成,使用谐波分析拟合程度及 F 检验分析兴济、平凉、大柏舍、昌黎地电场,可进一步分析大地电场成分。3.1谐波拟合程度Sk=12C2kni=1Sk=ni=112C2k选取电场空间电磁环境较为平静、且地电场观测数据日
16、变形态较为稳定的 2021 年 10 月 16 日数据为例进行分析。各测向谐波拟合的解释率可以用各谐波的方差表示,即振幅平方的一半。用前 n 阶谐波来拟合观测数据则总的方差解释率为,则拟合程度 为解释率与观测k=ni=1Sk/序列(1 440 个分钟值)方差的比值:。兴济台 3 个测向前 10 阶谐波拟合程度均较低,大柏舍、平凉前 5 阶拟合程度就已经达到 90%以上。昌黎台绘制的一天中总大地电场方位角变化与兴济台相似,数据离散度高,其拟合程度同样较低(表 1)。3.2F 检验谐波函数可以检测数据样本的周期,第 k 阶对应周期 T=1 440/k,前 10 阶谐波拟合分别对应周期为 2324
17、h、12 h、7.9 h、6 h、4.8 h、4 h、3.4 h、3 h、2.7 h、和 2.4 h 的潮汐谐波。可构造 F 检验分析谐波拟合周期的显著程度。构造统计函数 F 检验:Fk=12C2k/2(s2C2k/2)/(n21)(10)表 1 各测向前 5、10、15、20、100 阶谐波拟合程度Table 1 The fitting degree of forward 5,10,15,20 and 100 order harmonics of each measurement 台站测向5/%10/%15/%20/%100/%方差兴济台N75E58.563.866.469.485.60.4
- 配套讲稿:
如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。