中国大陆西部强震前超低频磁场时空变化特征

冯丽丽, 冯志生, 樊文杰, 管贻亮, 贺曼秋, 李霞, 何畅, 廖晓峰, 艾萨・伊斯马伊力, 袁文秀, 李莎

冯丽丽,冯志生,樊文杰,管贻亮,贺曼秋,李霞,何畅,廖晓峰,艾萨・伊斯马伊力,袁文秀,李莎. 2021. 中国大陆西部强震前超低频磁场时空变化特征. 地震学报,43(3):359−375. DOI: 10.11939/jass.20200093
引用本文: 冯丽丽,冯志生,樊文杰,管贻亮,贺曼秋,李霞,何畅,廖晓峰,艾萨・伊斯马伊力,袁文秀,李莎. 2021. 中国大陆西部强震前超低频磁场时空变化特征. 地震学报,43(3):359−375. DOI: 10.11939/jass.20200093
Feng L L,Feng Z S,Fan W J,Guan Y L,He M Q,Li X,He C,Liao X F,Aisa Y,Yuan W X,Li S. 2021. Spatio-temporal variation characteristic of the ultra-low frequency magnetic field prior to strong earthquakes of western Chinese mainland. Acta Seismologica Sinica43(3):359−375. DOI: 10.11939/jass.20200093
Citation: Feng L L,Feng Z S,Fan W J,Guan Y L,He M Q,Li X,He C,Liao X F,Aisa Y,Yuan W X,Li S. 2021. Spatio-temporal variation characteristic of the ultra-low frequency magnetic field prior to strong earthquakes of western Chinese mainland. Acta Seismologica Sinica43(3):359−375. DOI: 10.11939/jass.20200093

中国大陆西部强震前超低频磁场时空变化特征

基金项目: 地震科技星火计划(XH20062)、中国地震局青海格尔木青藏高原内部地球动力学野外科学观测研究站专项课题“青藏高原岩石圈磁场特征分析”和中国地震局监测预报司震情跟踪定向工作任务(2020020401)联合资助
详细信息
    通讯作者:

    冯志生: e-mail:fengzs2001@sohu.com

  • 中图分类号: P318.5

Spatio-temporal variation characteristic of the ultra-low frequency magnetic field prior to strong earthquakes of western Chinese mainland

  • 摘要: 对2015年至2019年期间中国大陆西部的磁通门磁力仪秒采样观测资料开展了5—100 s频段的地磁垂直强度极化分析,并运用一些数学方法对分析结果进行了处理。结果显示,极化高值在经向和纬向均无明显的形态和幅值变化,且极化高值与地磁外源场扰动无关。在此基础上筛选出18次极化高值异常事件,利用插值方法得到了极化高值异常的空间分布图,并分析了18次高值异常与中国大陆西部及周边15次强震的时空关系。分析结果表明:地磁场出现多台同步极化高值现象后的半年内,高值区可能发生M6.0以上强震;高值现象出现后,多个高值区均有可能发生强震;后续强震的震级与高值区面积呈正相关。
    Abstract: Based on the one second sampling data from fluxgate magnetometers in western Chinese mainland from 2015 to 2019, this paper carried out the vertical component polarization analysis of the frequency band between 5 s to 100 s, and then processed the analysis results by some mathematical methods. The results show that the high polarization value has no obvious shape and amplitude change in both meridional and zonal directions, and the high polarization value has nothing to do with geomagnetic field disturbance. On this basis, 18 high-value anomalies were screened out, and their spatial distribution map was obtained by interpolation method. The results show that within half a year after the synchronous appearance of multiple high-value anomalies of the vertical component polarization of geomagnetic field, the high-value region may have a strong earthquake with magnitude over M6.0. After the appearance of high value phenomenon, several high value regions are likely to have strong earthquakes; furthermore the magnitude of subsequent strong earthquake seems to be positively correlated with the area of the high-value zone.
  • 近年来利用卫星技术观测地震前后的电离层异常变化受到了地震学界的广泛关注, 并被认为是研究地震电磁前兆的有效途径之一(Larkina et al, 1983Chmyrev et al, 1989Molchanov et al, 1993Parrot, 1994). 尤其自2004年首颗专门用于地震观测的DEMETER卫星发射和应用以来, 地震电离层前兆观测研究积累了更多的经验和数据(Parrot et al, 2006aParrot, 2011Malchanov et al, 2006Bhattacharya et al, 2007张学民等, 2009, 张学民, 2011泽仁志玛等, 2010安张辉等, 2010欧阳新艳等, 2011).

    2008年5月12日汶川MS8.0地震、 2010年1月12日海地MS7.3地震和2010年2月27日智利MS8.8地震是近几年最受关注的三大强震, 利用DEMETER卫星数据对这3次强震的研究已有不少文献报道, 其中汶川地震的研究最多, 发现的电离层异常现象涉及到VLF电磁场(张学民等, 2009b曾中超等, 2009), ELF电场(张蓓等, 2010), LF电场(朱涛, 王兰炜, 2011), 以及相关等离子体密度和温度等参量(安张辉等, 2010武安绪等, 2011张学民等, 2009c董建等, 2009). 智利地震前后的电离层异常现象, 主要有电磁场及部分等离子体参量的同步扰动(Zhang et al, 2011刘静等, 2011Pisa et al, 2011)、 高能粒子通量的异常变化(黄建平等, 2010). 海地地震的研究较少, 仅焦其松等(2011)发现了地震发生当天多个离子参量的剧烈扰动.

    以上研究为3次强震前后电离层异常现象的探索提供了一定的研究经验和基础, 但均是针对单次地震的研究, 未关注同一参量在3次地震前后是否有相似或不同的表现. 本文利用DEMETER卫星数据, 研究对异常反应较敏感的离子温度(Ti)和VLF电磁场频谱在汶川、 海地和智利等3次强震前后所表现出的异常现象、 特征及其之间的异同.

    法国地震电磁探测卫星(detection of electro-magnetic emission from earthquake regions, 简称为DEMETER卫星), 自2004年6月29日发射后, 在轨期间共积累了近6年半的电离层观测资料. 卫星采用准太阳同步圆形轨道, 轨道倾角98.3°, 高度710 km, 其主要载荷包括电场探测仪(ICE)、 感应式磁力仪(IMSC)、 等离子体分析仪(IAP)、 朗缪尔探针(ISL)和高能粒子探测仪(IDP)(Thibery et al, 2006), 用来探测不同频段的电磁场和电离层等离子体参数.

    总结卫星数据震例研究结果(Gokhberg et al, 1983Serebryakova et al, 1992Parrot et al, 2006bBhattacharya et al, 2009), 电磁场在地震前后表现出的异常现象较多, 在DEMETER卫星记录的电场和磁场数据中, 仅VLF频段范围同时记录的巡查和详查两种模式(Lagoutte, 2006)数据最全, 因此本文首先选择了VLF频段的电磁场频谱数据进行分析. 此外, 笔者曾做过上述3次地震多参量异常现象的研究(颜蕊, 2010), 发现多参量典型异常现象出现较多的日期与Ti时间序列异常现象出现的日期十分吻合, 因此同时选择 了Ti作为等离子体参量进行研究. 电场探测仪(ICE)观测的VLF电场频谱数据提供了分辨率为19.531 25 Hz、 频段范围为19.531 25—20 000 Hz的功率谱密度, 单位为(μV2·m-2)/Hz, 共1 024个频点(Berthelier et al, 2006b). 感应式磁力仪(IMSC)观测的VLF磁场频谱频点分配与电场类似: 以19.531 25 Hz为间隔, 划分19.531 25—20 000 Hz的频段范围, 单位为nT2/Hz (Parrot et al, 2006b). 等离子体分析仪(IAP)直接提供了离子温度(Ti)数据(Berthelier et al, 2006a).

    汶川、 海地和智利3次地震参数见表1. 依据前人研究结果(朱荣, 2007欧阳新艳, 2008), 本文选择震前9天、 地震当天及震后2天, 共12天的临震时段进行研究(表2); 空间选择以震中为中心, 40°(纬度)×40°(经度)范围内的数据, 轨道与震中分布如图1所示. 观测时段均选择了夜侧数据(升轨)以尽量减少太阳活动对数据的影响, 因此所选数据地方时均为20:00—24:00, 即3次地震发生当天数据均为地震发生后记录的.

    表  1  汶川、海地、智利地震相关信息表
    Table  1.  Information on the Wenchuan,Haiti and Chile earthquakes
    下载: 导出CSV 
    | 显示表格
    表  2  汶川、海地和智利地震选取日期、轨道
    Table  2.  Date and orbits selected for the Wenchuan,Haitiand Chile earthquakes
    下载: 导出CSV 
    | 显示表格
    图  1  汶川、 海地、 智利地震震中与选择轨道关系
    Figure  1.  Relative epicenter location and satellite orbits for the Wenchuan, Haiti and Chile earthquakes

    同时, 通过ΣKpDst两个指数排除了地磁活动的影响.Kp指数是单个地磁台用来描述每日每3小时内的地磁扰动强度的指数,ΣKp指数则表示每天8个Kp指数之和,如果ΣKp>30, 则认为当天地磁活动强烈; Dst是描述磁暴时变化指数, 一般认为Dst指数高于20 nT时为平静期(徐文耀, 2003). 分析发现, 三次地震的Dst指数均高于20 nT、 ΣKp均小于20(图2)(DstKp数据均来自京都世界地磁数据中心).因此震例研究时段均处于地磁平静期. 此外, 还需根据卫星辅助文件DATA_RELATED_EVENT记录, 对受影响的轨道进行剔除, 本文中3次地震选取的轨道数据均不存在此种干扰, 因此不需剔除.

    图  2  汶川、海地、智利震例研究选取时段与地磁指数Dst和ΣKp
    Figure  2.  Geomagnetic indices Dst stand ΣKpduring the time period selected for Wenchuan,Haiti and Chile earthquakes

    强震前后Ti的研究, 首先提取时间序列中出现的异常现象, 判断异常出现时间, 然后做时空演化分析, 找出异常出现的位置.

    图34, 5分别给出了汶川、 海地和智利地震前后12天Ti 的时间序列, 横轴为日期,纵轴为Ti值. 选取数据为震中周围20°而非整轨数据, 因此区域内若存在温度异常变化, 则与区域内地震关系较大. 图67, 8对应3次强震前后Ti的时空演化图, 将研究时段内的数据按每3天分割, 绘成等值线. 其中图6a对应汶川地震前9—7天(2008-05-03—05)的数据, 图6b为震前6—4天(2008-05-06—08)的数据, 图6c为震前3—1天(2008-05-09—11)的数据, 图6d为是汶川地震当天(2008-05-12) 至震后2天(2008-05-13, 14)的数据. 同样图7a,b,c,d分别对应海地地震2010-01-03—05, 2010-01-06—08, 2010-01-09—11, 2010-01-12—14的数据; 图8a,b,c,d分别对应海地地震2010-02-18—20, 2010-02-21—23, 2010-02-24—26, 2010-02-27—03-01的数据.

    图  3  2008年汶川地震前后离子温度(Ti)时间序列显示图. 椭圆内为异常区
    Figure  3.  The time series of ion temperature (Ti) before and after the 2008 Wenchuan earthquake. Ellipse represents anomalous area
    图  4  2010年海地地震前后离子温度(Ti)时间序列图. 椭圆内为温度异常升高区
    Figure  4.  The time series of ion temperature (Ti) before and after the 2010 Haiti earthquake. Ellipse represents area of abnormal Ti rise
    图  5  2010年智利地震前后离子温度(Ti)时间序列图. 小椭圆内为震前异常区, 大椭圆内为震后温度升高区
    Figure  5.  The time series of ion temperature (Ti) before and after the 2010 Chile earthquake. Small ellipse represents area of abnormal disturbance, big ellipse stands for area of abnormal Ti rise
    图  6  汶川地震Ti时空演化图(图的周期均为3天)红色十字表示震中, 白色矩形框内为异常区
    Figure  6.  Ti spatio-temporal evolution of Wenchuan earthquake. The period of all maps is 3 days. Red cross stands for the epicenter, white rectangle represents anomalous area
    图  7  海地地震Ti时空演化图(图的周期均为3天) 红色十字符表示震中, 白色矩形框为异常区
    Figure  7.  Ti spatio-temporal evolution of Haiti earthquake. The period of all maps is 3 days. Red cross stands for the epicenter, white rectangle represents anomalous area
    图  8  智利Ti时空演化图(图的周期均为3天)红色十字表示震中, 白色矩形框为异常区
    Figure  8.  Ti spatio-temporal evolution of the Chili earthquake. The period of all maps is 3 days.Red cross stands for the epicenter, white rectangle represents anomalous area

    图3可以看出, 汶川地震前后Ti不太稳定, 温度变化波动较大, 以2008-05-09的突然升高最为明显(椭圆框内), 高于1 500 K, 平时Ti一般低于1200 K; 对应图6c(2008-05-09—11)的时空演化可以看出, 震中北偏西方向有明显的异常高值区(白色矩形框内). 可初步判断, 汶川震前3天(2010-05-09), Ti在震中北偏西方向出现异常高值, 疑似为地震前兆. 海地地震前后(图4), Ti比较稳定, 由于震中所处纬度相对较低, 因此整体值都比较高, 多在2 250 K左右, 地震发生当天(2010-01-12), 即震后4—5小时, Ti突然飙升, 幅度达50%(椭圆框内); 对应图7d震后1—2天演化图可见, 震中东南方向(白色框内)出现明显高值, 可见地震发生对Ti的影响还是非常大的. 智利地震前后, Ti变化相对汶川地震较稳定, 但相对海地地震较动荡, 多集中在800—1 200 K之间(图5); 震前9天(2010-02-18), Ti表现出剧烈扰动, 值比平时略微升高, 达到1 500 K左右(小椭圆框内); 震后发生了明显的飙升, 达到3 000 K以上(大椭圆框内); 对应图8a, d可以看出, 震前(2010-02-18)的扰乱在震中北东方向表现出了扰乱团, 震后震中正北方向有明显的高值区(白色矩形框内).

    地震电磁扰动的频率范围从直流一直到高频均有覆盖, 由于低频电磁波(ULF, ELF和VLF)在传播过程中的低衰减性, 地震前后更容易被低轨卫星观测到. 例如, OGO-6卫星观测到的100, 216和467 Hz磁场强度增强现象(Gokhberg et al, 1983); Cosmos-1809 卫星资料在140和450 Hz频段观测到电磁信号增强现象(Serebryakova et al, 1992); Parrot(1994)(2011)及Parrot等(2006a)利用Aureol-3卫星的电磁场观测数据发现了电场(72 Hz) 和磁场(140 Hz)分别出现信号增强现象, 并且通过DEMETER卫星统计多次地震, 观测到130 Hz频段的电磁辐射, 极低频段(<100 Hz)的电场静电紊流现象, 以及多次截止频率下降等现象; DEMETER卫星电磁场数据350和400 Hz附近在震前1—2天也出现辐射增强现象(Bhattacharya et al, 2009张学民等, 2009d). 可见, 前人研究的震前地震电磁场异常现象大都集中在500 Hz以下的低频段, 并且多表现为信号增强. 仅有Nemec等(2009)通过改进的统计方法得到1.7 kHz电场频谱震前0—4小时有小幅度的减弱现象.

    本文对VLF电磁场频谱的研究首先提取出单频点曲线, 从单频点曲线中筛出震中上空或震中磁共轭点上空有明显异常扰动的频点, 然后以此频点为中心确定异常出现频段范围. 文中重点选取了500 Hz以下幅度有增强现象的电磁场异常现象进行分析. 从异常出现的时间、 空间以及前人研究经验来看, 与地震的发生更相关, 具体见图9-图11. 图中横坐标为纬度, 纵坐标对应不同单频点谱值, 垂直虚线标注震中位置.

    图  9  汶川地震前后VLF电场与磁场频谱曲线(垂直虚线标注地震位置, 箭头标注异常突升位置)(a) 2008年5月13日电场频谱曲线; (b) 2008年5月8日磁场频谱曲线
    Figure  9.  Single frequency spectrum curves of VLF electromagnetic field before and after the Wenchuan earthquake. Vertical dotted lines mark earthquake location; arrows mark abnormal sudden rise(a) The electric field frequency spectrum curve on May 13, 2008; (b) The magnetic field frequency spectrum curve on May 8, 2008
    图  10  2010年1月12日海地地震当天VLF电场与磁场频谱曲线(垂直虚线标注地震位置, 箭头标注异常突起)(a) 电场频谱曲线; (b) 磁场频谱曲线
    Figure  10.  Single frequency spectrum curve of VLF electromagnetic field on the day of Haiti earthquake. Vertical dotted lines mark earthquake location, arrows mark abnormal sudden rise(a) The electric field frequency spectrum curve on January 12, 2010; (b) The magnetic field frequency spectrum curve on January 12, 2010
    图  11  智利地震前VLF磁场频谱单频点曲线(a) 2010-02-23电场频谱曲线; (b) 同天磁共轭点位置处电场频谱曲线(垂直虚线标注地震位置, 箭头标注异常突起)
    Figure  11.  Single frequency spectrum curve of VLF magnetic field before Chile earthquake (a) is the magnetic field spectrum curve on February 23 around the epicenter, and (b) is the curve around the conjugate point of epicenter in the same day. Vertical dotted lines mark earthquake location, arrows mark abnormal sudden rise

    图9为汶川地震前后VLF电磁场单频点曲线. 其中图9a为汶川震后(2008-05-13)电场单频点的频谱曲线, 自上至下频点分别对应19, 39, 60, 80和100 Hz. 可以明显看出, 在震中以北2°范围内, 集中在39—80 Hz频段范围内有异常突升(箭头处), 比平时背景值升高达3个lg(μV2 · m-2 · Hz-1)量级, 即高达1 000 (μV2 · m-2)/Hz; 80 Hz以上频段, 异常突升逐渐不明显. 图9b为汶川地震前4天(5月8日)磁场单频点频谱曲线, 自上至下对应100, 120, 137, 195和215 Hz. 可以看出, 在震中以北5°范围内磁场功率谱值有明显的突升(箭头处), 约1—2个lg(nT2 · Hz-1) 量级, 此突升在100—120 Hz频段范围内非常明显; 但高于120 Hz频段, 背景值逐渐增大, 则此异常不再突出.

    海地地震发生当天(2010>-01-12, 震后约4—5小时)的观测数据, VLF电场和磁场均发现了异常反应(图10). 图10a为电场频谱数据, 在距离震中以南5°范围内, 19—160 Hz频段范围出现了电场频谱突升(箭头处), 比正常谱值增大约3—4个lg(μV2 · m-2 · Hz-1) 量级; 图10b为磁场频谱数据, 在120—480 Hz频段范围之间, 增高约1—2个lg(nT2 · Hz-1) 量级. 由图10可见, 海地地震1月12日的VLF电场与磁场频谱异常突升的位置相同, 可能均为地震发生后能量释放所引起的共同反应. 可见, 地震的发生会引起对应空间电磁场的同步扰动, 但影响的频段范围不同, 仅有一定的频段交集.

    图11为智利地震震前VLF磁场单频点频谱曲线: 震前4天(2010-02-23)VLF磁场频谱在震中偏南5°内出现突升(图11a箭头处), 升幅比正常谱值高约1个lg(nT2 · Hz-1) 量级, 集中在100—160 Hz频段范围内. 同时, 震中磁共轭点位置也出现同步异常, 表现为磁共轭点以南1°以内发生突升(图11b箭头处), 升高约为2—3个lg(nT2 · Hz-1) 量级. 但共轭点位置处的异常频率与震中点异常频率不一致, 不好确定是同一信号源造成的, 因此此处仅是将出现的客观现象描述, 并不能确定此种效应为共轭效应.

    1) 笔者曾做过此3次强震多参量异常现象的分析研究(颜蕊, 2010), 发现多参量出现同步扰动较多的日期与Ti长时间序列异常现象出现的日期多比较吻合, 如汶川地震震前3天(2008-05-09)及智利地震震前9天(2010-02-18), 因此初步判断Ti的变化与地震发生的相关性较大, 尤其是2008年5月9日汶川地震前的异常现象, 在很多研究中都发现了其它参量的同步扰动(张学民等, 2009c曾中超等, 2009董建等, 2009杨剑等, 2011). 利用Ti可在后续工作中更多积累震例, 以及延长时间序列等进行更深入地研究. 此外, 根据Ti的时间序列可判别异常出现的时间, 通过时空演化分析, 则可以确定异常出现的空间位置, 是一种比较好的判别前兆时空位置的方法.

    2) 就此3次强震前后VLF电磁场频谱异常研究来看, VLF电磁场单频点频谱异常大部分出现在震中上空, 相对Ti来说, 电磁场异常出现的位置与震中位置相关性更大; 震前的VLF电场和磁场频谱异常无论时空位置还是频段几乎都不同步, 不同地震有不同的电磁异常现象反映, 可能会有电场异常, 也可能会有磁场异常, 也可能两种异常都有反应, 不容易确定; 但海地地震当天(2010-01-12, 震后约4—5小时)的VLF电场和磁场频谱均出现了同步异常, 异常出现的频段也有交叉. 本文认为: 震前的电场和磁场扰动并不是同步的, 两者之间需要经过一个相互耦合的过程, 因此造成了电磁场异常现象出现的时空及频段都存在一定的差异, 但两者之间的相互作用是不可分割的. 这种耦合反应在震前震后会有不同, 与地震孕育时和地震发生后参量的改变不同有关.

    3) 由文中数据结果可以看出, 汶川、 海地和智利3次地震震后都出现了明显的异常现象. 尤其是海地地震当天(2010-01-12, 震后约4—5小时)的VLF电场和磁场频谱均出现了同步异常, 其异常出现的位置和时段是一致的, 异常出现的频段也有交叉. 本文认为: 海地地震震后的异常强度很大, 表现激烈, 更多可能是地震发生后能量释放造成的电离层电磁状态的紊乱, 可见海地地震发生后能量释放对电离层有较大的影响, 能量释放也更集中. 虽然震后的异常现象不能说明地震前兆问题, 但是可帮助积累地震引起的各参量的异常现象的经验, 有助于对震前与地震有关的异常现象的判断.

    4) 通过本文的3个参量研究发现: 汶川、 海地和智利3次强震前后并未出现相似的异常现象; 由于发震机制等各种情况的不同, 其地震前后的表现各不相同, 所以共有的规律较难探索. 相比之下, 3次强震中海地地震震级小一些, 所以震前未表现出明显的异常现象, 也许正是因为震前未有更多能量释放, 才将能量一直累积在震后, 所以震后的表现更明显; 汶川和智利地震震级较大, 震前异常表现也更明显, 震后也有一定的反映. 当然, 震前震后异常与地震的发震机理、 震级和震源深度等综合因素紧密相关, 不易从单方面解释; 对岩石圈-大气层-电离层耦合机理(lithosphere-atmosphere-ionosphere, 简写为LAI)的认识不足也是目前很多现象不能得到更好解释的重要原因.

    Pulinets和Boyarchuk(2004)综合了多位科学家的资料和观点认为, 地震孕育的不同阶段, 电离层表现出不同的异常特征, 离子密度的异常出现比较早, 而VLF电磁信号在临震阶段则更为突出, 反映临震阶段电磁信号辐射增强, 这应该也与电离层中离子密度和温度的增加是联系在一起的. 本文临震前的离子温度和电磁场信号的增强与Pulinets和Boyarchuk的这一观点吻合. 但地震孕育与电离层异常信号之间的耦合不是单一模型能够完全解释的, 而是一个复杂的过程. 地震的孕育激发电离层异常的机制目前还不明确, 这就给地震异常的识别带来极大困难. 因此需要对不同地震事件期间的电离层现象进行更为广泛和深入的分析, 才能给出更确切的结果. 本文仅是将这些观测到的现象客观记录了下来, 对客观出现的震前异常现象进行总结和分析, 并在分析过程中尽可能地排除了太阳、 地磁等已知因素引起的电离层异常的可能. 当然目前还只是很初步的研究分析结果, 对这些信息的起源及空间电磁信息源场的判定等都还需要更多事件和经验的积累.

    1) 本文基于法国DEMETER卫星观测的离子温度(Ti)以及VLF电磁场单频点频谱数据, 探索了2008年5月12日汶川MS8.0, 2010年1月12日海地MS7.3和2010年2月27日智利MS8.8等3次强震前后与地震有关的电离层异常现象. 结果发现, 汶川地震前3天(2008-05-09)在震中北偏西方向离子温度明显升高, 震前4天(2008-05-08)VLF磁场低于200 Hz频段范围的频谱有明显突升. 智利地震前9天(2010-02-18)震中北东方向离子温度(Ti)有剧烈扰动; 震前4天(-02-23)VLF磁场100—160 Hz频段范围内突升. 海地地震震前没有观测到明显的异常现象, 但地震发生当天的数据(震后约4—5小时)中Ti, 40—160Hz频段电场频谱, 以及120—480 Hz磁场频谱均有明显突升, 应为地震发生后能量释放所引起的. 汶川和智利地震震前明显的异常现象, 虽然目前没有直接的证据表明这些变化是由地震孕育引起的, 但在数据处理中已尽可能排除了太阳、 地磁等因素的影响, 并且与前人(Gokhberg et al, 1983Serebryakova et al, 1992Parrot et al, 19942006aBhattacharya et al,2009)研究结果吻合. 因此, 本文研究的电磁辐射信号增强及温度的增高与地震孕育发生的关系较大.

    2) 本文研究的可能与地震有关的3个参量异常都具有很强的局部性, 多出现在震中区上空周围15°范围内, VLF电磁场单频点频谱的异常更多地出现在震中上空5°范围以内. 这种区域性特征与震源区孕震过程相关, 也就是说是孕震区的活动影响了电离层的各参数变化, 孕震区是作为一个扰动信号的主动源存在的, 这对追溯信号来源和分析异常产生机理是非常重要的.

    3) 由于电离层异常在地震前后出现的时间不确定, 卫星获取的数据时空覆盖也不够全面, 所以距离震中最近的轨道并不一定表现出最明显的异常, 距离发震时间最近的轨道也未必表现出最明显的异常, 所以卫星获取的数据, 其异常的出现在大多数情况下都比较随机, 较难寻找规律. 并且震例积累越多, 地震所表现出的异常现象也越多, 规律越难寻找. 若将来发射多颗同类卫星的星座, 弥补卫星获取时空密度上的不足, 会有助于解决此问题.

    作者衷心感谢法国DEMETER数据中心为本文提供的1级观测数据; 感谢中国地震局地震预测研究所同行给予的帮助; 感谢评审专家提出的修改意见和建议.

  • 图  1   2011年至2019年国家地磁台网中心磁通门磁力仪数量

    Figure  1.   The number of fluxgate magnetometers in National Geomagnetic Networks Center of China from 2011 to 2019

    图  2   2019年中国西部磁通门磁力仪(秒采样)分布

    Figure  2.   Distribution of fluxgate magnetometers (second sampling) in western Chinese mainland in 2019

    图  3   喀什台2017年YZH1YZH2的对比

    Figure  3.   Comparison of YZH1 with YZH2 results of Kashi station in 2017

    图  4   沿经线和纬线分布的部分台站及其YZHYZH1 2017年结果

    Figure  4.   Location of some stations and their YZH and YZH1 results in 2017

    图  5   喀什台2017年YZH1序列与Kp (a)和Dst (b)指数曲线的对比

    Figure  5.   Comparison of YZH1 series with Kp (a) and Dst (b) curves in Kashgar station in 2017

    图  6   YZH1高值异常在12个月中的频次分布

    Figure  6.   The frequency distribution of YZH1 high value anomalies in 12 months

    图  7   YZH1高值异常频次与地震频次图

    Figure  7.   The frequency of YZH1 high value anomalies and earthquake frequency

    图  8   2015年至2019年部分台站YZH1时序曲线(图中数字序号与表2中的异常序号相对应)

    Figure  8.   YZH1 time series curves of some stations from 2015 to 2017 (The numbers correspond with the series numbers of the anomalies in Table 2

    图  9   2015年至2016年磁场高值日YZH2的空间分布图

    Figure  9.   Spatial distribution of YZH2 on the magnetic field high value day from 2015 to 2016

    图  10   2017年至2019年磁场高值日YZH2的空间分布图

    Figure  10.   Spatial distribution of YZH2 on the magnetic field high value day from 2017 to 2019

    图  11   喀什台、静海台和高邮台的空间分布(a)以及2017年极化YZH1时序曲线对比图(b)

    图中红色竖线标示中国大陆西部YZH1出现高值异常的时间

    Figure  11.   The locations of the stations Kashi,Jinghai and Gaoyou (a) and YZH1 time series curves in 2017 (b)

    The red vertical line indicates the time when the YZH1 high value anomaly appeared in western Chinese mainland

    图  12   2015—2019年6组高值异常叠加前、后的YZH2与后续6个月内地震的空间关系

    Figure  12.   Spatial relationship of six groups of high YZH2 anomalies before and after superposition as well as the distribution of earthquakes in the following six months from 2015 to 2019

    图  12   2015—2019年6组高值异常叠加前、后的YZH2与后续6个月内地震的空间关系

    Figure  12.   Spatial relationship of six groups of high YZH2 anomalies before and after superposition as well as the distribution of earthquakes in the following six months from 2015 to 2019

    图  13   叠加后的YZH2高值区面积与后续6个月内地震震级的关系

    Figure  13.   The relationship between YZH2 high-value areas of the superposition of several high-value anomalies and the magnitude of strong earthquakes in six months

    表  1   2015年至2020年1月中国大陆西部及周边M≥6.0强震信息

    Table  1   Information of strong earthquakes (M≥6.0) during the period from 2015 to January of 2020 in western Chinese mainland and its surrounding areas

    发震时刻
     年-月-日 时:分:秒
    M北纬/°东经/°震源深度/km参考位置
    2020-01-19 21:27:55 6.4 39.83 77.21 16 新疆伽师县
    2019-11-21 07:50:45 6.0 19.50 101.35 10 老挝
    2019-06-17 22:55:43 6.0 28.34 104.9 16 四川长宁县
    2019-04-24 04:15:48 6.3 28.40 94.61 10 西藏墨脱县
    2017-11-18 06:34:19 6.9 29.75 95.02 10 西藏米林县
    2017-08-09 07:27:52 6.6 44.27 82.89 11 新疆精河县
    2017-08-08 21:19:46 7.0 33.2 103.82 20 四川九寨沟县
    2016-12-08 13:15:03 6.2 43.83 86.35 6 新疆呼图壁县
    2016-11-25 22:24:30 6.7 39.27 74.04 10 新疆阿克陶县
    2016-10-17 15:14:49 6.2 32.81 94.93 9 青海杂多县
    2016-06-26 19:17:11 6.7 39.43 73.4 10 吉尔吉斯斯坦
    2016-01-21 01:13:13 6.4 37.68 101.62 10 青海门源县
    2015-12-07 15:50:03 7.4 38.2 72.9 30 塔吉克斯坦
    2015-07-03 09:07:46 6.5 37.6 78.2 10 新疆皮山县
    2015-04-25 14:11:26 8.1 28.2 84.7 20 尼泊尔
    下载: 导出CSV

    表  2   2015年至2019年中国大陆西部$Y_{ZH2}$高值异常信息

    Table  2   Information of high YZH2 value anomalies in western Chinese mainland from 2015 to 2019

    序号异常出现时间
     年-月-日
    异常持续时间/d异常台站最多日(异常日)
     年-月-日
    异常台站数总台数YZH2
    12015-01-1662015-01-1828360.73
    22015-02-0992015-02-1314360.39
    32015-07-1752015-07-188360.29
    42015-09-2652015-09-2913270.47
    52015-11-18112015-11-2321270.92
    62016-01-2782016-01-3017410.75
    72016-08-1492016-08-1715330.92
    82016-09-1192016-09-1514330.88
    92016-12-0162016-12-0214370.75
    102017-01-1352017-01-158380.41
    112017-02-1352017-02-1410370.49
    122017-03-1192017-03-1814340.66
    132017-08-2462017-08-2822400.65
    142017-10-2962017-10-3014400.70
    152018-04-1632018-04-1711430.62
    162018-06-1062018-06-1217430.44
    172019-09-2232019-09-2413541.15
    182019-11-1252019-11-1226560.85
    下载: 导出CSV

    表  3   单次高值异常及短期内多次高值异常叠加后发生强震的高值区面积

    Table  3   The high value areas of single high-value anomalies and the superposition of several high-value anomalies in a short period which followed with strong earthquakes

    序号出现高值日期
    年-月-日
    震后出现高值区的
    地震(提前的天数)
    半年内高值区强震高值区面积
    /(104 km2
    高值区总面积
    /(104 km2
    1 2015-01-18     无 2015-07-03皮山M6.5 251.6 313
    2 2015-02-13 2015-04-25尼泊尔M8.1
    3 2015-07-18 2015-07-03皮山M6.5
    (15)
    2015-12-07塔吉克斯坦M7.4 13 20.6
    4 2015-09-29     无 2016-01-21门源M6.4 20 144
    5 2015-11-23 2015-12-07塔吉克斯坦M7.4 123.7
    6 2016-01-30 2016-01-21门源M6.4 (9) 2016-06-26吉尔吉斯斯坦M6.7 17.5 84.7
    7 2016-08-17     无 2016-12-08呼图壁M6.2 24.1 82.2
    8 2016-09-15 2016-10-17杂多M6.2 58
    9 2016-12-02     无 2016-12-08呼图壁M6.2 18.2 80.8
    10 2017-01-13 2016-12-08呼图壁M6.2 (37)     无 47.2 47.2
    11 2017-02-14     无 2017-08-08九寨沟M7.0 41.6 84.4
    12 2017-03-18 2017-08-09精河M6.6 30.5
    13 2017-08-28 2017-08-08九寨沟M7.0 (19)
    2017-08-09精河M6.6 (20)
        无 39.5
    14 2017-10-31     无 16.1
    15 2018-04-17     无     无 32.4
    16 2018-06-12
    17 2019-09-24     无 2019-11-21老挝M6.0 19.2 72.3
    18 2019-11-12
    注:① 指地震震中所在高值区的面积;② 指一次高值现象出现的多个高值区面积之和。
    下载: 导出CSV
  • 杜爱民,周志坚,徐文耀,杨少峰. 2004. 新疆和田ML7.1地震前ULF电磁辐射的激发机理[J]. 地球物理学报,47(5):832–837. doi: 10.3321/j.issn:0001-5733.2004.05.014

    Du A M,Zhou Z J,Xu W Y,Yang S F. 2004. Generation mechanisms of ULF electromagnetic emissions before the ML=7.1 earthquake at Hotan of Xinjiang[J]. Chinese Journal of Geophysics,47(5):832–837 (in Chinese). doi: 10.3321/j.issn:0001-5733.2004.05.014(inChinese)

    冯志生,李琪,卢军,李鸿宇,居海华,孙海军,杨福喜,张翼. 2010. 基于磁通门秒值数据的地震ULF磁场可靠信息提取研究[J]. 华南地震,30(2):1–7. doi: 10.3969/j.issn.1001-8662.2010.02.001

    Feng Z S,Li Q,Lu J,Li H Y,Ju H H,Sun H J,Yang F X,Zhang Y. 2010. The seismic ULF geomagnetic reliable information exaction based on fluxgate magnetometer data of second value[J]. South China Journal of Seismology,30(2):1–7 (in Chinese).

    何畅,冯志生. 2017. 极化方法在成都地磁台的应用研究[J]. 地震学报,39(4):558–564. doi: 10.11939/jass.2017.04.010

    He C,Feng Z S. 2017. Application of polarization method to geomagnetic data from the station Chengdu[J]. Acta Seismologica Sinica,39(4):558–564 (in Chinese).

    李琪,杨星,蔡绍平. 2015. 极化方法应用于地磁台阵的震例分析[J]. 震灾防御技术,10(2):412–417. doi: 10.11899/zzfy20150222

    Li Q,Yang X,Cai S P. 2015. Case study of applying polarization method to geomagnetic array data[J]. Technology for Earthquake Disaster Prevention,10(2):412–417 (in Chinese).

    马宗晋. 1980. 华北地壳的多(应力集中)点场与地震[J]. 地震地质,2(1):39–47.

    Ma Z J. 1980. Multipoints (concentrated) stress field and earthquakes in North China crust[J]. Seismology and Geology,2(1):39–47 (in Chinese).

    梅世蓉. 1996. 地震前兆场物理模式与前兆时空分布机制研究(二):强震孕育时应力、应变场的演化与地震活动、地震前兆的关系[J]. 地震学报,18(1):1–10.

    Mei S R. 1996. Study on physical model of seismic precursors and time-space distribution mechanism of precursors(2)[J]. Acta Seismologica Sinica,18(1):1–10 (in Chinese).

    潘晖,张建国,杨冬梅,李美. 2014. 地震地磁低点位移成因的一种可能机制[J]. 大地测量与地球动力学,34(4):83–87. doi: 10.14075/j.jgg.2014.04.029

    Pan H,Zhang J G,Yang D M,Li M. 2014. A possible mechanism of the contributing factor of geomagnetic lowest shift[J]. Journal of Geodesy and Geodynamics,34(4):83–87 (in Chinese).

    汤吉,詹艳,王立凤,徐建郎,赵国泽,陈小斌,董泽义,肖骑彬,王继军,蔡军涛,徐光晶. 2008. 5月12日汶川8.0级地震强余震观测的电磁同震效应[J]. 地震地质,30(3):739–745. doi: 10.3969/j.issn.0253-4967.2008.03.012

    Tang J,Zhan Y,Wang L F,Xu J L,Zhao G Z,Chen X B,Dong Z Y,Xiao Q B,Wang J J,Cai J T,Xu G J. 2008. Coseismic signal associated with aftershock of the MS8.0 Wenchuan earthquake[J]. Seismology and Geology,30(3):739–745 (in Chinese). doi: 10.3969/j.issn.0253-4967.2008.03.012(inChinese)

    汤吉,詹艳,王立凤,董泽义,赵国泽,徐建郎. 2010. 汶川地震强余震的电磁同震效应[J]. 地球物理学报,53(3):526–534. doi: 10.3969/j.issn.0001-5733.2010.03.006

    Tang J,Zhan Y,Wang L F,Dong Z Y,Zhao G Z,Xu J L. 2010. Electromagnetic coseismic effect associated with aftershock of Wenchuan MS8.0 earthquake[J]. Chinese Journal of Geophysics,53(3):526–534 (in Chinese). doi: 10.3969/j.issn.0001-5733.2010.03.006(inChinese)

    姚休义,滕云田. 2017. 2010年元谋MS5.2地震前地磁极化异常研究[J]. 地震研究,40(3):444–448. doi: 10.3969/j.issn.1000-0666.2017.03.021

    Yao X Y,Teng Y T. 2017. Variation of geomagnetic polarization value associated with the 2010 Yuanmou MS5.2 earthquake[J]. Journal of Seismological Research,40(3):444–448 (in Chinese).

    姚休义,冯志生. 2018. 地震磁扰动分析方法研究进展[J]. 地球物理学进展,33(2):511–520. doi: 10.6038/pg2018BB0063

    Yao X Y,Feng Z S. 2018. Review on the recent development of analysis methods on magnetic disturbance associated with earthquakes[J]. Progress in Geophysics,33(2):511–520 (in Chinese).

    袁桂平,张学民,吴迎燕,赵旭东. 2015. 汶川8.0级地震前地磁低点位移与内外源Sq等效电流体系关系的研究[J]. 地震,35(3):102–112. doi: 10.3969/j.issn.1000-3274.2015.03.011

    Yuan G P,Zhang X M,Wu Y Y,Zhao X D. 2015. Minimum point shift of the geomagnetic vertical component in diurnal variation and the internal-external equivalent current system Sq before the 2008 Wenchuan MS8.0 earthquake[J]. Earthquake,35(3):102–112 (in Chinese).

    Chang X T,Zou B,Guo J Y,Zhu G B,Li W,Li W D. 2017. One sliding PCA method to detect ionospheric anomalies before strong earthquakes:Cases study of Qinghai,Honshu,Hotan and Nepal earthquakes[J]. Adv Space Res,59(8):2058–2070. doi: 10.1016/j.asr.2017.02.007

    Cianchini G,De Santis A,Barraclough D R,Wu L X,Qin K. 2012. Magnetic transfer function entropy and the 2009 MW=6.3 L’Aquila earthquake (Central Italy)[J]. Nonlinear Proc Geophys,19(4):401–409. doi: 10.5194/npg-19-401-2012

    Currie J L,Waters C L. 2014. On the use of geomagnetic indices and ULF waves for earthquake precursor signatures[J]. J Geophys Res,119(2):992–1003. doi: 10.1002/2013JA019530

    Gladychev V,Baransky L,Schekotov A,Fedorov E,Pokhotelov O,Andreevsky S,Rozhnoi A,Khabazin Y,Belyaev G,Gorbatikov A,Gordeev E,Chebrov V,Sinitsin V,Lutikov A,Yunga S,Kosarev G,Surkov V,Molchanov O,Hayakawa M,Uyeda S,Nagao T,Hattori K,Noda Y. 2001. Study of electromagnetic emissions associated with seismic activity in Kamchatka region[J]. Nat Hazards Earth Syst Sci,1(3):127–136. doi: 10.5194/nhess-1-127-2001

    Han P,Hattori K,Xu G J,Ashida R,Chen C H,Febriani F,Yamaguchi H. 2015. Further investigations of geomagnetic diurnal variations associated with the 2011 off the Pacific coast of Tohoku earthquake (MW9.0)[J]. J Asian Earth Sci,114:321–326. doi: 10.1016/j.jseaes.2015.02.022

    Hattori K. 2004. ULF geomagnetic changes associated with large earthquakes[J]. Terr Atmos Ocean Sci,15(3):329–360. doi: 10.3319/TAO.2004.15.3.329(EP)

    Hattori K,Han P,Yoshino C,Febriani F,Yamaguchi H,Chen C H. 2013. Investigation of ULF seismo-magnetic phenomena in Kanto,Japan during 2000−2010:Case studies and statistical studies[J]. Surv Geophys,34(3):293–316. doi: 10.1007/s10712-012-9215-x

    Hayakawa M,Kawate R,Molchanov O A,Yumoto K. 1996. Results of ultra-low-frequency magnetic field measurements during the Guam earthquake of 8 August 1993[J]. Geophys Res Lett,23(3):241–244. doi: 10.1029/95GL02863

    Hayakawa M,Itoh T,Hattori K,Yumoto K. 2000. ULF electromagnetic precursors for an earthquake at Biak,Indonesia on February 17,1996[J]. Geophys Res Lett,27(10):1531–1534. doi: 10.1029/1999GL005432

    Huang Q H,Ikeya M. 1998. Seismic electromagnetic signals (SEMS) explained by a simulation experiment using electromagnetic waves[J]. Phys Earth Planet Inter,109(3/4):107–114. doi: 10.1016/S0031-9201(98)00135-6

    Ismaguilov V S,Kopytenko Y A,Hattori K,Voronov P M,Molchanov O A,Hayakawa M. 2001. ULF magnetic emissions connected with under sea bottom earthquakes[J]. Nat Hazards Earth Syst Sci,1(1/2):23–31. doi: 10.5194/nhess-1-23-2001

    Karakelian D,Klemperer S L,Fraser-Smith A C,Thompson G A. 2002. Ultra-low frequency electromagnetic measurements associated with the 1998 MW5.1 San Juan Bautista,California earthquake and implications for mechanisms of electromagnetic earthquake precursors[J]. Tectonophysics,359(1/2):65–79.

    Li Q,Zhu P Y,Mamatemin A,Xu X G. 2011. Detection of ULF electromagnetic emissions as a precursor to two earthquakes in China[J]. Earthquake Science,24(6):601–607. doi: 10.1007/s11589-011-0822-2

    Masci F. 2011. On the seismogenic increase of the ratio of the ULF geomagnetic field components[J]. Phys Earth Planet Inter,187(1/2):19–32. doi: 10.1016/j.pepi.2011.05.001

    Matsushima M,Honkura Y,Oshiman N,Bariş Ş,Tunçer M K,Tank S B,Çelik C,Takahashi F,Nakanishi M,Yoshimura R,Pektaş R,Komut T,Tolak E,Ito A,Iio Y,Işikara A M. 2002. Seismoelectromagnetic effect associated with the Izmit earthquake and its aftershocks[J]. Bull Seismol Soc Am,92(1):350–360. doi: 10.1785/0120000807

    Molchanov O A,Kopytenko Y A,Voronov P M,Kopytenko E A,Matiashvili T G,Fraser-Smith A C,Bernardi A. 1992. Results of ULF magnetic field measurements near the epicenters of the Spitak (MS=6.9) and Loma Prieta (MS=7.1) earthquakes:Comparative analysis[J]. Geophys Res Lett,19(14):1495–1498. doi: 10.1029/92GL01152

    Molchanov O A,Schekotov A,Fedorov E,Belyaev G,Gordeev E. 2003. Presismic ULF electromagnetic effect from observation at Kamchatka[J]. Nat Hazards Earth Syst Sci,3(3/4):203–209. doi: 10.5194/nhess-3-203-2003

    Nagao T,Enomoto Y,Fujinawa Y,Hata M,Hayakawa M,Huang Q,Izutsu J,Kushida Y,Maeda K,Oike K,Uyeda S,Yoshino T. 2002. Electromagnetic anomalies associated with 1995 Kobe earthquake[J]. J Geodynam,33(4/5):401–411.

    Prattes G,Schwingenschuh K,Eichelberger H U,Magnes W,Boudjada M,Stachel M,Vellante M,Wesztergom V,Nenovski P. 2008. Multi-point ground-based ULF magnetic field observations in Europe during seismic active periods in 2004 and 2005[J]. Nat Hazard Earth Syst Sci,8(3):501–507. doi: 10.5194/nhess-8-501-2008

    Thomas J N,Love J J,Johnston M J S,Yumoto K. 2009. On the reported magnetic precursor of the 1993 Guam earthquake[J]. Geophys Res Lett,36(16):L16301. doi: 10.1029/2009GL039020

    Uyeda S,Hayakawa M,Nagao T,Molchanov O,Hattori K,Orihara Y,Gotoh K,Akinaga Y,Tanaka H. 2002. Electric and magnetic phenomena observed before the volcano seismic activity in 2000 in the Izu island region,Japan[J]. Proc Natl Acad Sci USA,99(11):7352–7355. doi: 10.1073/pnas.072208499

    Wang Z D,Chen B,Yuan J H,Yang F X,Jia L,Wang C. 2018. Localized geomagnetic field anomalies in an underground gas storage[J]. Phys Earth Planet Inter,283:92–97. doi: 10.1016/j.pepi.2018.08.009

    Wessel P,Smith W H F. 1995. New version of the generic mapping tools[J]. Eos Trans AGU,76(33):329.

  • 期刊类型引用(13)

    1. 赵彬彬,余怀忠,贾东辉,尼鲁帕尔·买买吐孙,钱才. 新疆及周边地震前电离层扰动异常特征统计分析. 华南地震. 2024(S1): 46-48 . 百度学术
    2. 李忠,杨百一,黄建平,陈朝阳,杨旭明,韩莹,张富志,鲁恒新. 空基电场数据处理方法综述. 地球物理学进展. 2023(02): 562-577 . 百度学术
    3. 赵彬彬,钱才,尼鲁帕尔·买买吐孙,张小涛,余怀忠. 2022年1月8日青海门源6.9级地震前电离层异常特征分析. 地震工程学报. 2022(03): 691-699 . 百度学术
    4. 杨静,胡耀垓,赵正予. 汶川地区夏季电离层VLF背景分析. 地球物理学进展. 2021(02): 531-539 . 百度学术
    5. 何苗,吴立新,崔静,王威,齐源,毛文飞,苗则朗,陈必焰,申旭辉. 汶川地震前多圈层短—临遥感异常回顾及其时空关联性. 遥感学报. 2020(06): 681-700 . 百度学术
    6. 宋程,张永仙,夏彩韵,张小涛. 2008年汶川M_S8.0地震电离层图像信息长时程演化初探. 地震. 2019(02): 135-146 . 百度学术
    7. 夏彩韵,张小涛. 利用改进型图像信息方法研究汶川8.0级地震震前电离层扰动异常特征. 防灾减灾学报. 2019(04): 25-30 . 百度学术
    8. 蔡润,武震,谭大诚,云欢,郭鹏. 地震前的电磁异常综述. 华南地震. 2018(01): 1-16 . 百度学术
    9. 郑香脂,张爱兵,关燚炳,刘超,王文静,田峥,孔令高,孙越强. 电磁监测试验卫星阻滞势分析器探测技术. 物理学报. 2017(07): 453-462 . 百度学术
    10. 郑香脂,张爱兵,关燚炳,刘超,孙越强,王文静,田峥,孔令高,丁建京. 电磁监测试验卫星离子漂移计探测技术. 物理学报. 2017(20): 295-303 . 百度学术
    11. 宋奕瑶,李忠,安建琴,张丽娟,史鹏飞. 基于SOM网的DEMETER卫星电场数据聚类异常分析. 软件导刊. 2016(08): 140-142 . 百度学术
    12. 欧阳新艳,申旭辉. DEMETER卫星观测的ULF电场波形干扰分析. 地震地磁观测与研究. 2015(04): 19-25 . 百度学术
    13. 欧阳新艳,申旭辉. DEMETER卫星观测的超低频电场扰动预处理方法及震例应用分析. 地震学报. 2015(05): 820-829 . 本站查看

    其他类型引用(8)

图(14)  /  表(3)
计量
  • 文章访问数:  831
  • HTML全文浏览量:  474
  • PDF下载量:  63
  • 被引次数: 21
出版历程
  • 收稿日期:  2020-06-04
  • 修回日期:  2021-03-01
  • 网络出版日期:  2021-07-06
  • 发布日期:  2021-05-14

目录

/

返回文章
返回