摘要:2019年6月17日在四川长宁县(104.9°E、28.34°N)发生6.0级地震。本文以该地震为例,基于中国地磁台网提供的地磁台站资料,利用帕金森矢量方法,计算出每个台站在不同频段的帕金森矢量方位角异常之后,将同一频段的多个台站的异常进行叠加,提取和分析该地震前震中周围地磁场变化特征。结果显示,震中周围多个台站短临地磁异常时空特征明显,通过异常综合叠加分析可以较好地估计震中方位。该研究结果可为今后短临地震预报提供一定参考价值。
根据电磁感应理论,理想导体在外加均匀垂直磁场的情况下,其感应产生的垂直磁场将与外加场相互抵消,其表面并无垂直磁场分量。因此,Parkinson等[1]学者曾提出,在假设均匀天然场源的情况下,地球表面观测到地磁日变垂直分量不为零的现象时需要考虑层状地球电导模型中存在的横向电性不均匀体。地磁加卸载响应比[2]、地磁转换函数[3]、地磁日变幅逐日比[4]、地磁谐波振幅比[5]和地磁低点位移[6]等方法都是基于此物理原理。这些方法在研究地震孕育和发生过程的震磁前兆异常中发挥了重要作用[7,8]。
前人研究指出,在地磁转换函数法中通常将转换函数表示成地磁感应矢量[9],并且其导出的帕金森矢量可指明地下感应电流相对集中或相对高导的区域[1],即中强地震发生前后地下的电磁参数(电导率)会发生变化,从而引起地磁场的异常变化。之后中国地震学者对其进行了广泛研究,并且发现帕金森矢量、转换函数及地下电性结构变化与中强地震的发生有较好的对应关系[10,11,12,13]。Chen等[14,15]利用帕金森矢量方法来评估中国台湾地区中强地震前地下电导率变化特征,据此通过地磁台站数据定出震中大致的方位。
本文以2019年四川长宁6.0级地震为例,基于中国地磁台网提供的震中附近地磁台站的连续波形数据,利用帕金森矢量方法和Chen等[14,15]的改进方法,获取了震中周围地磁感应矢量的时间变化序列,提取并分析了其中的长期背景变化特征,讨论了地下岩石出现电导率增大的现象及其地磁感应矢量空间分布特征,定出震中方位,从而分析研究此次地震前后地磁场变化异常特征。
1、发震构造与资料选取
1.1发震构造情况
据中国地震台网测定,北京时间2019年6月17日,在青藏高原东缘四川盆地南缘宜宾市长宁县(104.9°E、28.34°N)发生6.0级地震,震源深度约19km,其后5d内相继多次发生5级左右强余震。根据矩张量反演结果、序列空间分布、震区构造特征以及灾害现场调查,确定此次长宁6.0级地震震中位于长宁—双河复式大背斜附近,其发生可能是由长宁—双河复式大背斜附近的白象岩—狮子滩背斜和双河场褶皱及其伴生的断层活动控制,震源机制显示为主逆冲型地震,发震断层整体呈NE—SW向挤压为主、兼具一定NW—SE向拉张分量的构造变形特征(图1)。震中处于华南地块内部现今构造活动相对较弱的四川盆地南缘盆山结合地带,震中100km范围内历史上发生过两次6级以上地震,本次长宁6.0级地震是继1974年5月11日云南大关7.1级地震后川东南地区发生的再一次强地震活动,造成重大人员伤亡及大量房屋财产受损[16],因此,研究该地震的地震前兆异常,对该地区地震预测研究及防灾减灾具有重要意义。
图1震中区构造及地磁台站分布示意
1.2资料选取
本文使用的地磁台站数据全部由中国地磁台网中心提供,震中300km范围内台站的分布及主要的断裂背景见图1。值得一提的是,考虑地磁资料易受人为和环境干扰,我们只使用震中周围秒采样地磁台站夜间23:00~5:00的数据,仪器类型均为GM4类型磁通门磁力仪(表1)。
表1震中周围地磁观测项目概况
2、计算方法与原理
Parkinson等[1]提出,地磁短周期变化的水平分量变化量ΔH、垂直分量变化量ΔZ和磁偏角变化量ΔD之间存在线性关系,这种关系也适用于ΔZ、ΔX(北分量)和ΔY(东分量)。
Z(f)=A(f)X(f)+B(f)Y(f) (1)
式中,X、Y、Z是地磁场三分量,A、B是地磁转换系数。
对地面各测点记录到的天然磁场信号进行分析得到的转换函数A和B是地点和时间(λ,φ)、深度h、频率ω或周期T的函数。如果垂直分量变化量Z有外源场部分,则转换函数A和B一般为复数,称为复转换函数:
{A(f)=Ar(f)+iAu(f)B(f)=Br(f)+iBu(f) (2)
式中,Ar(f)和Au(f)、Bu(f)分别表示转换函数A(f)和B(f)实部和虚部。那么,帕金森矢量的方位角和幅度随之定义为
Pa=arctanBr(f)Ar(f)+180°Pm=[Ar(f)2+Br(f)2]1/2 (3)
式中,Pa是帕金森矢量的方向,指示地下电导率高的方向(0度代表正北方向);Pm是帕金森矢量的大小,反映地下电性横向不均匀的程度,分析和研究A、B和Pa、Pm随时间的变化,可了解和监测地震活动区地下介质电磁参数的变化特征[13,14]。
3、异常特征分析与讨论
本文分别选取震中周围300km以内的四川省马边台、重庆市涪陵台、贵州省贵阳台的2018—2019年UTC23:00~05:00的秒采样全波形数据,其中我们将2018全年的地磁数据计算所得到的帕金森矢量方位角作为正常背景值(以图2中马边台为例),以15d为移动窗长,1d为步长,计算窗口内帕金森矢量方位角分布比例,再扣除背景值消除长期变化并做归一化处理,将其作为窗口内最后1d的方位角异常,计算地震前后共45d,各台站在0.005~0.01Hz帕金森矢量方位角异常的日变化。值得注意的是,地球磁场由外源场和内源场两部分组成,在讨论震磁效应时着重研究内源场及其变化,需要消除外源场的影响,而外源场最明显的特征即为周期变化[17]。因此,在对数据进行分析的过程中,消除年变是非常必要的,本文提取震磁异常前,首先使用傅里叶拟合方法,扣除其长周期变化,从而达到消除外源场变化影响,同时提高异常识别的准确性。在识别震磁信息时,排除变化磁场本身的背景场及外空场的影响等干扰因素。
图22018年马边台三分量(D、H、Z)地磁场日变化
图3分别给出了本文选取的3个地磁台(马边台、涪陵台、贵阳台)0.005~0.01Hz频段的转换函数长期变化及其导出的帕金森矢量及背景值。反映了研究区3个台站的帕金森矢量在电性不均匀结构的响应,并且该分布具有明显的方向性。马边台、涪陵台和贵阳台的帕金森矢量的方向分别集中在103°、242°和275°左右,这给我们研究长宁6.0地震地磁异常指明了方向。
图3研究区地磁台站帕金森矢量背景值
五角星为地震;红线为帕金森矢量
图4给出了涪陵台在长宁6.0级地震前后异常变化曲线。震前21d开始在240°左右的方向显示比较明显的短临前兆异常,异常持续到地震当天,震后快速消失。追加往前更长时间或震后该方向没有显示类似变化,且异常方向跟该台站帕金森矢量方向基本一致,且其他方向未出现明显的高频异常变化,因此异常信度较高,推测为本次地震的短临异常信号;震中较近的马边台、贵阳台则未出现类似异常,这也反映了地震前兆异常的复杂性。但是从单台异常很难圈定震中位置,因此我们对研究区不同方向的3个地磁台站帕金森矢量方向进行异常叠加,以对地震前后异常演化进行分析(图5)。
图4涪陵台震前30天和震后15天异常变化
五角星为地震;红线为帕金森矢量
选取某一天的异常分布,将平面以10°间隔分成36份。将台站某一方向和其背向异常值较高的值同时作为这两个方向的异常值,得到单一台站的异常空间分布。将不同方向所有台站异常叠加,取以各台站为中心300km范围内至少有3个台站异常值相交的区域。图5为长宁6.0级地震前30d和震后15d所选取三个台站的异常叠加综合图,其中色标表示异常大小。从图5不难看到,震前14d图中几乎看不到明显的高异常区或高异常区特别地零星分布,但是震前13d开始这些零星分布的异常在震中附近汇聚,呈现很明显的高异常集中区,持续到震前1d,尤其是震前3~1d高异常区域特别集中,长宁6.0级地震就发生在这个高异常集中区内。长宁6.0级地震后,震中周围仍然没有异常区域。这些结果为长宁6.0级地震震源附近存在相对高电导率提供了积极的证据。
4、结论
本文使用马边台、涪陵台和贵阳地磁台站三分量磁力仪的秒资料,计算出帕金森矢量并借以找出2019年6月17日长宁6.0级地震前,地下电磁场(电导率)的变化。结合前人对地震前地磁场(电导率)异常的研究[14,15,17],笔者可以得到以下结论:
(1)从单台的帕金森矢量分析结果看,涪陵台此次地震前这种异常变化特别明显,异常特性为短临异常,即震前22d以前帕金森矢量所指向的频段几乎无异常,到震前22d开始可以看到明显的高异常出现到消失的现象,且震中对台站指向的方位角基本吻合,地震异常信度高。
图5长宁6.0级地震前30d和震后15d所选取三个台站的异常叠加综合示意图
(2)从不同方向不同角度的多个地磁台站帕金森矢量方向进行异常叠加图看,2019年长宁6.0级地震前后异常演化特征较明显,震前13d开始逐渐出现感应矢量高幅度异常,尤其是震前3d异常幅度达到最高值,震后立即消失,且震中就落到高异常叠加区。本文采用帕金森矢量方位角异常指示的位置与实际震中的位置最接近。
(3)由于地下构造活动的复杂性以及地表环境的干扰等因素,并不是所有的地震在震前都会发生这样的前兆异常,同时也并不是所有的异常之后都会有相应的地震发生[18]。但是,本文的研究结果显示,无论在空间上还是时间上该地震前后磁力场变化都很显著,即震中周围地磁台站表现出最初无异常,到出现高频异常再到异常消失的过程。因此,帕金森矢量在提取本次长宁6.0级地震前兆信号中的积极作用是值得肯定的。
笔者综合分析认为,基于地磁短周期变化,利用地磁秒数据计算分析,在0.005~0.01Hz频段范围内,地震前后20d左右帕金森矢量方位角分布出现明显异常,并指向震中附近。对于目前地震预报难题,前兆异常需要多方面的综合研究,本文研究结果可作为一种地磁方面的参考。
参考文献:
[2]尹祥础,尹灿.非线性系统失稳的前兆与地震预报——响应比理论及其应用[J].中国科学B辑,1991,21(5):512-518.
[3]龚绍京,陈化然,张翠芬,等.地磁水平场转换函数在唐山地震前的异常反应[J].地震学报,1997,19(1):51-58.
[4]冯志生,王建宇,蒋延林,等.地磁垂直分量日变幅逐日比及其与地震关系的探讨[J].华南地震,2001,21(2):20-27.
[5]冯志生,梅卫萍,张秀霞,等.中强震前地磁谐波振幅比的趋势性变化特征初步研究[J].西北地震学报,2004,26(1):50-56.
[6]丁鉴海,车时,余素荣,等.地磁日变地震预报方法及其震例研究[M].北京:地震出版社,2009:1-3.
[7]戴勇,高立新,姚丽,等.2013年岷县—漳县6.6级地震地磁垂直分量日变化低点时间空间分布特征[J].地震,2019,39(3):106-114.
[8]倪晓寅,赖见深,朱培育.汶川8.0级地震前地磁日变化异常综合分析[J].中国地震,2019,35(1):65-75.
[9]王桥,黄清华,等.华北地磁感应矢量时空特征分析[J].地球物理学报,2016,59(1):215-228.
[10]龚绍京.广东省地磁台的帕金森矢量及广州台的系数在河源地震前后的时间变化[J].地震研究,1987,10(5):575-582.
[11]于明鑫,曾小苹,董玉兰,等.大连地磁短周期变化特征的初步研究[J].地震地磁观测与研究,1990,(6):35-41+21.
[12]曾小苹,林云芳,朱忠杰,等.用地磁转换函数研究震源区介质的电性变化[J].地震学报,1995,(3):335-339.
[13]李琪,杨冬梅,张素琴,等.地磁转换函数方法之典型震例及数模对比[J].地震,2005,(3):85-90.
[16]易桂喜,龙锋,梁明剑,等.2019年6月17日四川长宁MS6.0地震序列震源机制解与发震构造分析[J].地球物理学报,2019,62(9):3432-3447.
[17]储贻晟.日本地区地震前电导率异常研究[D].国立中正大学,2017.
[18]张国民,傅征祥,桂燮泰.地震预报引论[M].北京:科学出版社,2001,357-362.
艾萨·伊斯马伊力,陈界宏,毛志强.利用帕金森矢量检测2019年长宁6.0级地震的地磁场变化[J].地震,2020,40(02):91-99.
基金:中国地震局震情跟踪青年项目(2019010409,2020010422);中国地震局监测、预报、科研三结合项目(3JH-201901097).
分享:
土壤磁学是将现代岩石磁学理论、技术方法应用到土壤学而产生的分支学科,通过研究土壤中磁性矿物的类型、含量及颗粒磁畴状态等特征,探索其在环境系统中的溶解、迁移和转化机制,利用物质在磁性特征上的联系及其反映的环境内涵,进而研究不同时空尺度的环境问题、环境过程和作用机制.随着中国黄土高原黄土/古土壤的磁学研究日趋成熟,南方土壤的磁性研究逐渐受到学者重视.
2020-09-05由于电离层层析反演存在着不适定问题,因此许多国内外学者提出了各种方法来解决。文献[1,2]针对GPS电离层层析反演中的不适定问题,先后提出了多种层析算法,取得了较好的反演结果。虽然算法的优化改进可在一定程度上解决不适定的问题,但并不能从根本上消除该问题所带来的负面影响。目前,越来越多的学者通过增加数据源的思路来解决该问题。
2020-07-09在太阳风-磁层-电离层耦合系统中,随着太阳风和行星际磁场的变化,近地空间的电流系统也会发生改变.在高纬地区,电离层电流的快速变化会引起地球表面的磁场发生扰动.地磁扰动是由于日地系统中太阳风能量、动量和物质输运至地球并与近地空间环境相互作用产生的一种地磁活动现象,地磁扰动主要包括磁暴、地磁亚暴和地磁脉动等变化磁场表征现象。
2020-07-09通过实际干扰数据分析可以看出,核自适应滤波算法作为一种非线性滤波方法,能通过实时调节滤波参数从而抑制干扰,为一种完全自适应的无需人工干预的滤波算法。相比于其他算法,该算法对于地磁数据的处理更具有适用性。本文首次将核自适应滤波算法应用于地磁干扰抑制研究,以期为地磁干扰抑制提供新方法。
2020-07-09在多系统精密单点定位中,各系统的观测值都需进行电离层高阶项改正。本文通过计算2014年与2017年不同纬度的3个MGEX测站中GPS、GLONASS、Galileo以及Beidou(BDS)卫星导航定位系统观测值的电离层高阶项延迟,分析电离层高阶项延迟对多系统精密单点定位结果的影响。
2020-07-08中国幅员辽阔,北方多地地磁台站位于高纬度、高海拔地区,冬季温度低,可能达不到地磁等观测仪器正常工作所需环境温度,影响观测数据质量,台站人员应实时监测并及时处理观测数据异常现象,应具备相应的职业素养,明确地磁台站记录数据干扰类型和表现形态,能快速准确地进行数据预处理,确保地磁观测仪器正常工作。
2020-07-06干扰使得观测数据精度降低,对观测资料造成很多大的影响,故在日常观测中对各种干扰的识别尤为重要。本文选取2015年1月~2018年12月的观测数据,分别统计观测环境因素、观测系统影响及人为干扰等各类干扰情况,分析各类干扰对地磁观测数据的影响特征,总结得出都兰地震台地磁FHD观测数据的干扰因素及干扰特征,并提出各类干扰的应对措施,以提高数据的观测质量。
2020-07-06固原4.6级地震是1998年7月29日海原4.9级地震后西海固地区发生的最大一次地震,也是海原4.9级地震以来西海固地区发生的唯一4级以上地震。本文主要对固原4.6级地震发生前的异常和预测指标应用等进行回顾与总结,以期为宁夏地区中等有感及破坏性地震的预测和震情监视跟踪工作提供借鉴。
2020-07-06通过相关方法将地磁数据中的地磁异常提取出来,再对提取出的异常数据进行相关分析。目前,地震预测以观测事实归纳总结为主,辅以物理模式演绎法,将一些地磁场影响因素通过数学物理建模进行定量的扣除。本文中通过相关系数分析、回归分析,确定温度对温泉地磁场子夜均值影响的数学表达式,从而剔除温度对地磁场数据分析中的影响。
2020-07-06在勘探过程中,钻探所得岩心是研究储层裂缝最基础也是最佳的工具,其最能直观反映储层的裂缝形态[8]。如果获得了定向煤心,便可通过CT扫描等方法来对煤样内部空间裂缝进行观测和定量描述,借此来确定煤层裂缝的空间展布,从而进一步服务于煤层气田开发。本文结合了油气勘探中岩心定向方法,提出了将古地磁岩心定向技术用于煤层气勘探开发中的新思路。
2020-07-06人气:2007
人气:1311
人气:1293
人气:1215
人气:1085
我要评论
期刊名称:大学物理
期刊人气:1310
主管单位:中国科学技术协会
主办单位:中国物理学会
出版地方:北京
专业分类:科学
国际刊号: 1000-0712
国内刊号:11-1910/O4
邮发代号:82-320
创刊时间:1982年
发行周期:月刊
期刊开本:大16开
见刊时间:10-12个月
影响因子:0.735
影响因子:0.645
影响因子:1.369
影响因子:0.874
影响因子:0.385
400-069-1609
您的论文已提交,我们会尽快联系您,请耐心等待!
你的密码已发送到您的邮箱,请查看!