2017年九寨沟7级地震前电磁扰动异常分析

第３９卷　第３期地 震Ｖｏｌ．３９，Ｎｏ．３ ２０１９年７月ＥＡＲＴＨＱＵＡＫＥ　Ｊｕｌ．，２０１９　２０１７年九寨沟７级地震前电磁扰动异常分析＊
黄继攀１，２，３，王新安１，２，３，雍珊珊１，２，王　墨１，３，刘　众１，３
（１．北京大学深圳研究生院集成微系统实验室，广东深圳　５１８０５５；２．北京大学深圳地震
监测预测研究中心，广东深圳　５１８０５５；３．北京大学信息科学技术学院，北京　１００８７１）
摘要：地震前的电磁扰动现象曾被多次观测到，但需进一步研究。

在中国地震局的支持下，可
以监测电磁扰动数据的多分量地震监测系统———ＡＥＴＡ在中国境内已经布设了２００余套。

２０１７年８月８日，四川省九寨沟县发生７级地震。

震中周边５００ｋｍ范围内当时已经布设２２套
ＡＥＴＡ设备。

基于滑动四分位距法（时间窗为２７ｄ）分析了这些台站震前３０ｄ的数据。

１１个台
站出现了异常，其中震中周边２００ｋｍ内的６个台站有４个发生异常。

同时，正异常明显多于
负异常。

另外，大部分异常发生在夜间（２０时～６时）。

关键词：九寨沟地震；电磁扰动；时间序列分析；异常
中图分类号：Ｐ３１５．６２ 文献标识码：Ａ 文章编号：１０００－３２７４（２０１９）０３－００８４－１１
引言
大地震是一种具有高度破坏性的自然灾害，经常造成巨大的人员伤亡和财产损失。

现代地震学中，已经有不少对地壳形变、应变应力、重力、电离层、氡气、地磁等的观测工作［１～７］。

考虑到地震的复杂性，应该观测更多的地震前兆信号。

电磁扰动虽然有前人研究，但基本上是研究了个别地震的个别台站的电磁扰动现象［８～１２］。

地声方面无论是古代文献［１３］还是地震经历者均有较多震例，但近年来相关的研究比较少，研究中震例和台站也较少［１４，１５］。

自２０１２年起，我们开始研制多分量地震监测系统ＡＥＴＡ（Ａｃｏｕｓｔｉｃ　ａｎｄ　ｅｌｅｃｔｒｏｍａｇｎｅｔ－ｉｃ　ｔｅｓｔｉｎｇ　ａｌｌ　ｉｎ　ｏｎｅ　ｓｙｓｔｅｍ），可以同时监测电磁扰动和地声信号［１６～１９］，可形成较密的前兆台网，让单震例下多台站的电磁扰动异常相互验证，同时实现了一定规模的地声前兆监测台网。

２０１７年８月８日，中国四川省九寨沟县发生ＭＳ７．０地震，震中为１０３．８２°Ｅ，３３．２０°Ｎ，震源深度２０ｋｍ。

地震发生时，我们在中国地震局的支持下已经在震中周边５００ｋｍ范围内布设了２２套ＡＥＴＡ设备。

本文将用滑动四分位距法分析这些台站的电磁扰动数据异常。

＊收稿日期：２０１８－１１－０５；修改回日期：２０１９－０３－１８
基金项目：深圳市科技创新计划项目（ＫＪＹＹ２０１７０７２１１５１９５５８４９；ＪＣＹＪ２０１７０４１２１５１１５９４６１）
作者简介：黄继攀（１９９１－），男，河南商丘人，博士研究生，主要从事地震监测预测研究。

通讯作者：王新安，教授。

Ｅ－ｍａｉｌ：ａｎｘｉｎｗａｎｇ＠ｐｋｕ．ｅｄｕ．ｃｎ
图１　ＡＥＴＡ系统框图
１　系统简介
首先，对多分量地震监测系统ＡＥ－
ＴＡ做个简介。

它有２个子系统，数据采集系统与数据分析系统。

（１
）数据采集系统：包括埋于地下２ｍ的电磁扰动探头和地声探头，以及地面机柜中的数据预处理终端；
（２
）数据分析系统：包括布设在阿里云上的云服务器、数据分析客户端和数据分析网页。

整个设备有２种———电磁扰动传感器和地声传感器。

电磁扰动传感器与法国Ｄｅｍｅｔｅｒ
卫星上ＩＭＳＣ传感器［２
０，２１］和中国电磁卫星上磁传感器［２２］
原理相同，基于法拉第电磁感应定律，感应出垂直方向变化的磁场，得到感应电动势，进而经过放大、滤波、模数转换等处理。

其性能指标如下：
（１）０．１Ｈｚ～１０ｋＨｚ，０．１ｎＴ～１　
０００ｎＴ电磁信号，灵敏度大于２０ｍＶ／ｎＴ＠０．１Ｈｚ～１０ｋＨｚ，１８ｂｉｔ分辨率，低频采样率５００Ｈｚ，高频采样率３０ｋＨｚ；（２）尺寸：直径９０ｍｍ，长度８５０ｍｍ；
（３
）垂直安装。

地声传感器基于压电薄膜传感器［
２３］
，地声信号经过压电薄膜感应出电荷，进而经过电荷放大器、滤波器、放大电路和ＡＤ转换等处理。

其性能指标如下：
（１）０．１Ｈｚ～５０ｋＨｚ，灵敏度３ＬＳＢ／ｐ
ａ＠０．１Ｈｚ～５０ｋＨｚ，１８ｂｉｔ分辨率，低频采样率５００Ｈｚ，高频采样率１５０ｋＨｚ
；（２）尺寸：直径２６５ｍｍ，高度６８ｍｍ；（３
）水平安装。

图２　探头实际安装情形（ａ）和数据处理终端（ｂ
） 在设备可靠性方面，Ａ
ＥＴＡ设备作为地震前兆监测的新型仪器，研发完成
后由卓翼科技［２４，２５］
进行生产。

ＡＥＴＡ设
备在卓翼科技出厂前进行了性能测试、老化测试和防水测试等，设备的可靠性和稳定性得到了保证。

电磁扰动传感器的频响测试、灵敏度测试、噪声水平测试和一致性测试在中科院河北怀来遥感综合试验站磁屏蔽室完成（２０１６年１２月）
，相关指标达到了要求。

这两种传感器埋于地下２ｍ。

ＡＥ－
ＴＡ设备绝大多数安装于县级和乡镇级
的中国地震局台站，周围的环境足够安静，尽可能地降低了诸如施工、高压电网等人类活动对设备的干扰。

图２展示了具体的安装情形和摆放ＡＥＴＡ终端的机柜。

５
８　３期 黄继攀等：
２０１７年九寨沟７级地震前电磁扰动异常分析　
截至２０１８年１２月，在中国地震局的支持下，已经在川滇地区、首都圈、广东和台湾地区布设了２００余套设备，图３是当前的仪器分布图。

图３　ＡＥＴＡ在中国分布
２　数据与分析方法
２．１　数据
本文的数据分析对象是九寨沟周边５００ｋｍ范围内的２２个台站的电磁扰动数据。

这２２个台站的位置详情见表１。

数据可以从我们的数据分析客户端或者网页（ｈｔｔｐ：∥ｗｗｗ．ａｅ－ｔａ．ｃｎ）获取。

使用的数据分量是低频（０．１～２００Ｈｚ
）电磁扰动数据，精度是１ｈ。

图４给出了地震前后６０天共１２０ｄ的原始数据，横坐标中０对应着地震发生当天（２０１７－０８－０８）。

可以看到平静状态的数据如ＳＦ和ＫＤＧＺ，每天较为规律的波动，如果波动幅值较小，长期来看近似一根直线。

平静状态下，不同台站的电磁扰动信号的幅值和形态虽不完全相同，但普遍存在日变特征信号（图４中幅值较大的ＪＺＧ、ＳＰ、ＭＸ、ＬＤ等台站可以明显地看到）。

２．２　分析方法
在地震前兆数据分析领域中，滑动四分位距法是一种流行的异常分析方法［
４，２６～２８］。

根据前人研究［
９，１０，１２，２９，３０］
，震前电磁扰动异常多在地震发生时到震前一个月之间，这里分析地震前后３０ｄ的异常值。

由于地下电磁扰动信号很可能受到太阳活动的影响［３１］
，这里采用太阳辐射周期２７ｄ［３
１，３２
］作为滑动时间窗口。

每次滑动步长为１ｄ。

在每个滑动窗口内，同一时刻的数据（比如每天７时的数据）组成一个数组，组内进行升序排序，处在１／４位置的数称之为Ｑ１，１／２位置的数为中位数Ｍｅｄｉａｎ，３／４位置的数为Ｑ３。

这样，滑动四分位距
ＩＱＲ便可以被计算，公式为
ＩＱＲ＝Ｑ３－Ｑ１（１）ＩＱＲ约等于１．３４倍的标准差。

在这里，使用２倍的ＩＱＲ作为异常检测阈值，相当于２．７倍的标准差，具有９９．３％的检测异常置信概率。

这样，采用式（２）和（３）便可以确定正常数据的上下限。

滑动窗口后那一天的数据（第２８ｄ的数据）ＥＭＤ如果在上下限之内，则认为
６
８ 地 震 ３
９卷　
表１　２２个ＡＥＴＡ台站的位置信息
台站名台站缩写经纬度东经／（°）北纬／（°）震中距
／ｋ
ｍ布设时间
台站名
台站缩写经纬度东经／（°）北纬／（°）震中距
／ｋｍ
布设时间
九寨沟ＪＺＧ（１０４．２５，３３．２６）４０．５５　２０１７－０６－
１０名山ＭＳ（１０３．１５，３０．１９）３４０．６４　２０１７－０６－
０４松潘ＳＰ（１０３．６０，３２．６５）６４．５１　２０１７－０６－１２宝兴ＢＸ（１０２．８２，３０．２５）３４１．３８　２０１７－０６－０５平武ＰＷ（１０４．５５，３２．４１）１１１．２３　２０１７－０６－０８康定姑咱ＫＤＧＺ（１０２．１７，３０．１２）３７６．３９　２０１６－１２－
２１青川ＱＣ（１０５．２３，３２．５９）１４８．０９　２０１７－０６－
０７泸定ＬＤ（１０２．２３，２９．９２）３９４．５９　２０１６－１２－
２２茂县ＭＸ
（１０３．８５，３１．６９）１６７．９３　２０１７－０６－
１３峨眉山ＥＭＳ（１０３．５０，２９．５９）４０２．５６　２０１７－０６－０５汶川ＷＣ（１０３．５９，３１．４８）１９２．４７　２０１７－０６－１３犍为ＱＷ（１０３．９４，２９．２１）４４３．８１　２０１７－０６－
０８马尔康ＭＥＫ（１０２．２２，３１．９０）２０８．２９　２０１７－０６－１３金口河ＪＫＨ（１０３．０８，２９．２２）４４８．１１　２０１７－０６－
０６什邡ＳＦ（１０４．１６，３１．１３）２３２．３９　２０１７－０６－０５石棉挖角ＳＭＷＪ（１０２．３４，２９．２３）４６３．３２　２０１６－１２－１９都江堰ＤＪＹ（１０３．６５，３０．９８）２４７．３７　２０１６－１２－２７石棉山洞ＳＭＳＤ（１０２．３４，２９．２３）４６３．３２　２０１６－１２－１８金川ＪＣ（１０２．０７，３１．４７）２５３．０５　２０１７－０６－
１２沐川ＭＣ（１０３．９０，２８．９６）４７１．５３　２０１７－０６－１１小金县
ＸＪＸ
（１０２．３６，３１．００）２８０．６３　２０１７－０６－
０５马边
ＭＢ
（１０３．５４，２８．８３）４８６．６５　２０１７－０４－
１２图４　各个台站地震前后６０ｄ的原始数据序列
７
８　３期 黄继攀等：
２０１７年九寨沟７级地震前电磁扰动异常分析　
是正常数据，超出上下限的数据则被认为是异常数据，超出的部分称为异常值ΔＥＭＤ。

ＵＢ＝Ｍｅｄｉａｎ＋２ＩＱＲ
（２
）ＬＢ＝Ｍｅｄｉａｎ－２ＩＱＲ
（３）Δ
ＥＭＤ＝ＥＭＤ－ＵＢ　（ＥＭＤ＞ＵＢ）　０　（ＬＢ＜ＥＭＤ＜ＵＢ）ＥＭＤ－ＬＢ　（ＥＭＤ＜ＬＢ烅烄
烆
）（４
）３　结果与讨论
３．１　结果
使用滑动四分位距法检测到的异常值，如图５所示。

同样，横坐标中０对应着地震发生当天。

从图５中，可以看到异常发生在ＪＺＧ、ＱＣ、ＭＸ、ＪＣ、ＥＭＳ、ＪＫＨ、ＭＣＳＰ、ＤＪＹ和
ＱＷ台站。

其中，九寨沟台站在震前３０ｄ的２０～２４时发生轻度异常，震前８ｄ的０～４时发生４个小时的负异常、１３～１８时发生５个小时的正异常，震前７ｄ到震后６ｄ每天的５～６时发生持续异常，形成了条带；青川台站在震前２７ｄ的６～８时出现明显异常，震前２０ｄ的２０～２４时出现明显异常；茂县台站在震前９ｄ、６ｄ、５ｄ、３ｄ和２ｄ的７～２０时出现约１０～２０ｍＶ异常；金川台站在震前８ｄ到震后７ｄ出现持续的负异常，幅值在－６０～－１０
ｍＶ之间；峨眉山台站在震前２５ｄ、１６ｄ、１５ｄ、１１ｄ、５ｄ和０ｄ的０～６时间出现正异常，震前１７ｄ和震前１ｄ的２０～２４时间出现正异常；金河口台站在震前２３ｄ和震前１９ｄ的０～６时间出现正异常；沐川台站的震前２４ｄ和１８ｄ的１０～２４时出现正异常，震前１２ｄ和震
前１ｄ的６～８时也出现正异常。

更清晰的异常情况，如表２所示。

表２　电磁扰动数据异常总结表
台站异常时间／ｄ异常时间（ＵＴＣ＋８ｈ
）异常幅值（ｍＶ）
ＪＺＧ－３０～－２８，－８～＋４　
０：００～６：００，１３：００～１８：００，２１：００～２４：００＋１０～＋４０，－２０～－５
ＱＣ－２６，－２０　６：００～８：００，２０：００～２４：００＋１０～＋３０
ＭＸ－９，－６～－５，－３～－２　
７：００～２０：００
＋１０～＋２０ＪＣ－８～－４，－１～＋７　０：００～２４：００
－６０～－１０ＥＭＳ－２５，－１６～－１５，－１１，－５，－１～０
０：００～６：００，２０：００～２４：００＋１０～＋４０ＪＫＨ－２３～－１９　０：００～６：００＋１０～＋４０ＭＣ－３０～－１８，－１２～－１　１０：００～２４：００，６：００～８：００
＋１０～＋２０ＳＰ－３０～－２９　
０：００～５：００，０：００～２４：００＋１０～＋１５ＤＪＹ－２８，－２６，－２３，－２０，－１６　
０：００～６：００，１８：００～２４：００＋１０～＋１５
ＱＷ
－２５～－１６　
０：００～５：００，２０：００～４：００
＋５～１５，－１５～－５
３．２　讨论
从表２可以看到，大部分异常发生在夜间（２０～６时）。

同时，在异常的台站中，正异常远远多于负异常，这很可能是这一区域地震前的电磁扰动异常规律。

地震通常被认为是板块运动等构造活动导致的，为了探究上述异常是否与活动地块边
界有关系，绘制图７。

图７中地块边界数据、板块移动方向和速度来自张培震等［３３，３４］
的活动构造数据。

从图７中可以看到九寨沟地震发生在巴颜喀拉地块的东侧边界上，结合表１和表２，震中周边２００ｋｍ内的６个台站有４个表现出了正异常。

华南一级地块边缘的ＥＭＳ、
８
８ 地 震 ３
９卷　
图５　基于滑动四分位距法检测到的低频电磁扰动异常
白色区域代表数据缺失
ＭＣ、ＪＨＫ和ＱＷ一带也表现出了正异常。

从图６上来看，异常集中出现在九寨沟震中周
边区域和峨眉山区域，已经用红色虚线圆圈进行表示。

关于电磁扰动异常的机理，目前尚未定论，动电效应和微破裂机制被认为是最可能的
机制［
１１，３５～３８］。

九寨沟震中周边区域部分台站异常的出现，很可能是动电效应或者微破裂导致的局部地下电性结构发生改变的结果。

峨眉山区域的几个台站处于在震前也出现异常，可能是地下电性结构发生改变。

９
８　３期 黄继攀等：
２０１７年九寨沟７级地震前电磁扰动异常分析　
图６　异常台站分布
３．３　地声信号和长期的电磁扰动信号变化情况
ＡＥＴＡ台站同时有电磁扰动探头和地声探头，九寨沟地震前可能地声也发生了异常。

于是对地声数据进行分析，分析对象为原始数据的均值、振铃计数和峰值频率，发现震前不少台站出现异常，如图７所示，红色五角星表示地震发生当天。

从图７中可以看到，ＪＺＧ台站震前９ｄ开始出现地声异常一直持续到地震发生当天，附近的ＰＷ和ＱＣ在震前１１ｄ和震前１ｄ均出现地声异常等。

为了进一步地理解震前信号的变化，对几个变化较为平稳信度较大的台站（ＪＺＧ和ＭＸ）的电磁扰动信号进行了较长时间（２０１７年６月８日～２０１８年７月３日）
的研究和分析，得到图８。

图８中左侧红色虚线对应着九寨沟７级地震发生时刻。

这段时间内的２０１７年１１月７日九寨沟发生４．５级强余震，２０１７年１１月１０日距离ＭＸ台站仅仅２０ｋｍ的绵竹发生４．２级地震，因时间较近图中用红色虚线仅标注了２０１７年１１月１０日绵竹地震时刻。

从图８中可以看到，距离九寨沟７级地震震中较近的ＪＺＧ和ＭＸ台站，在２０１７年８月地震发生时，处于较长期正常数据数值偏低的时期，但在短期内又达到了高的峰值；ＪＺＧ和ＭＸ两台站的信号在２０１７年１１月初的２次地震时同样出现了由低值到高值的变化，与九寨沟７级地震时候的波形具有一定的相似性。

０
９ 地 震 ３
９卷　
图７　
九寨沟地震震前地声异常统计图示
图８　几个台站一年以上电磁扰动原始数据
４　结语
本文进一步印证了大地震周边地区多处会出现电磁扰动异常，多分量地震监测系统
ＡＥＴＡ监测到了可能与九寨沟大震相关的低频电磁扰动异常，积累了同一大震下多个台站电磁扰动的数据和数据异常经验。

ＡＥＴＡ台站的密集布设有望进一步促进地震的监测预测研究，更多关于电磁扰动异常的研究值得去做。

参考文献：
［１］　马瑾，Ｓｈｅｒｍａｎ　Ｓ　
Ｉ，郭彦双．地震前亚失稳应力状态的识别———以５°拐折断层变形温度场演化的实验为例［Ｊ］．中国科学：地球科学，２０１２，４２（５）：６３３－
６４５．［２］　Ｍａｔｓｕｏ　Ｋ，Ｈｅｋｉ　Ｋ．Ｃｏｓｅｉｓｍｉｃ　ｇｒａｖｉｔｙ　
ｃｈａｎｇｅｓ　ｏｆ　ｔｈｅ　２０１１Ｔｏｈｏｋｕ－Ｏｋｉ　ｅａｒｔｈｑｕａｋｅ　ｆｒｏｍ　ｓａｔｅｌｌｉｔｅｇｒａｖｉｍｅｔｒｙ［Ｊ］．Ｇｅｏｐｈｙ
ｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｌｅｔｔｅｒｓ，２０１１，３８（７）：１１３－１２０．［３］　Ｋｉｍ　Ｊ　Ｗ，Ｈａｎ　Ｙ　Ｊ，Ｋｉｍ　Ｒ，ｅｔ　ａｌ．Ｉｎｖｅｓｔｉｇａｔｉｏｎ　ｏｆ　ｔｈｅ　ｒｅｌａｔｉｏｎｓｈｉｐ　
ｂｅｔｗｅｅｎ　ｅａｒｔｈｑｕａｋｅｓ　ａｎｄ　ｉｎｄｏｏｒｒａｄｏｎ　ｃｏｎｃｅｎｔｒａｔｉｏｎｓ　ａｔ　ａ　ｂｕｉｌｄｉｎｇ　ｉｎ　Ｇｙｅｏｎｇｊｕ，Ｋｏｒｅａ［Ｊ］．Ｎｕｃｌｅａｒ　Ｅｎｇｉｎｅｅｒｉｎｇ　＆Ｔｅｃｈｎｏｌｏｇｙ，２０１８，５０：５１２－
５１８．［４］　Ｌｉｕ　Ｊ　Ｙ，Ｃｈｅｎ　Ｙ　Ｉ，Ｈｕａｎｇ　
Ｈ　Ｋ，ｅｔ　ａｌ．Ｉｏｎｏｓｐｈｅｒｉｃ　ｆｏＦ２ａｎｄ　ＴＥＣ　ａｎｏｍａｌｏｕｓ　ｄａｙｓ　ａｓｓｏｃｉａｔｅｄ　ｗｉｔｈ　Ｍ≥５．０ｅａｒｔｈｑｕａｋｅｓ　ｉｎ　Ｔａｉｗａｎ　ｄｕｒｉｎｇ　１９９７—１９９９［Ｊ］．Ｔｅｒｒｅｓｔｒｉａｌ　Ａｔｍｏｓｐｈｅｒｉｃ　＆Ｏｃｅａｎｉｃ　Ｓｃｉｅｎｃｅｓ，２００４，１５（３）：３７１－
３８３．［５］　晏锐，蒋长胜，张浪平．汶川８．
０级地震前水氡浓度的临界慢化现象研究［Ｊ］．地球物理学报，２０１１，５４（７）：１　８１７－１　
８２６．［６］　江在森，方颖，武艳强，等．汶川８．
０级地震前区域地壳运动与变形动态过程［Ｊ］．地球物理学报，２００９，５２（２）：５０５－
５１８．１
９　３期 黄继攀等：
２０１７年九寨沟７级地震前电磁扰动异常分析　
［７］　丁鉴海，刘杰，余素荣，等．地磁日变化异常与强震的关系［Ｊ］．地震学报，２００４，２６（ｓ１）：８０－８８．［８］　Ｃ　Ｆ　Ａ，Ａ　Ｂ，Ｒ　Ｍ　Ｐ，ｅｔ　ａｌ．Ｌｏｗ－ｆｒｅｑｕｅｎｃｙ　
ｍａｇｎｅｔｉｃ　ｆｉｅｌｄ　ｍｅａｓｕｒｅｍｅｎｔｓ　ｎｅａｒ　ｔｈｅ　ｅｐｉｃｅｎｔｅｒ　ｏｆ　ｔｈｅ　ＭＳ７．１Ｌｏｍａ　Ｐｒｉｅｔａ　Ｅａｒｔｈｑｕａｋｅ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｌｅｔｔｅｒｓ，１９９０，９（１７）：１　４６５－１　４６８．［９］　Ｐｏｔｉｒａｋｉｓ　Ｓ　Ｍ，Ｓｃｈｅｋｏｔｏｖ　Ａ，Ａｓａｎｏ　Ｔ，ｅｔ　ａｌ．Ｎａｔｕｒａｌ　ｔｉｍｅ　ａｎａｌｙｓｉｓ　ｏｎ　ｔｈｅ　ｕｌｔｒａ－ｌｏｗ　ｆｒｅｑｕｅｎｃｙ　
ｍａｇ－ｎｅｔｉｃ　ｆｉｅｌｄ　ｖａｒｉａｔｉｏｎｓ　ｐｒｉｏｒ　ｔｏ　ｔｈｅ　２０１６Ｋｕｍａｍｏｔｏ（Ｊａｐａｎ）ｅａｒｔｈｑｕａｋｅｓ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ａｓｉａｎ　ＥａｒｔｈＳｃｉｅｎｃｅｓ，２０１８，１５４：４１９－
４２７．［１０］　Ｐｏｔｉｒａｋｉｓ　Ｓ　Ｍ，Ｈａｙａｋａｗａ　Ｍ，Ｓｃｈｅｋｏｔｏｖ　Ａ．Ｆｒａｃｔａｌ　ａｎａｌｙｓｉｓ　ｏｆ　ｔｈｅ　ｇ
ｒｏｕｎｄ－ｒｅｃｏｒｄｅｄ　ＵＬＦ　ｍａｇｎｅｔｉｃｆｉｅｌｄｓ　ｐｒｉｏｒ　ｔｏ　ｔｈｅ　１１Ｍａｒｃｈ　２０１１Ｔｏｈｏｋｕ　ｅａｒｔｈｑｕａｋｅ（ＭＷ＝９）：ｄｉｓｃｒｉｍｉｎａｔｉｎｇ　ｐｏｓｓｉｂｌｅ　ｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｓ　ｆｒｏｍ　ｓｐ
ａｃｅ－ｓｏｕｒｃｅｄ　ｄｉｓｔｕｒｂａｎｃｅｓ［Ｊ］．Ｎａｔｕｒａｌ　Ｈａｚａｒｄｓ，２０１７，８５（１）：５９－８６．［１１］　张建国，焦立果，刘晓灿，等．汶川ＭＳ
８．０级地震前后ＵＬＦ电磁辐射频谱特征研究［Ｊ］．地球物理学报，２０１３，５６（４）：１　２５３－１　
２６１．［１２］　李美，谭捍东，王中平，等．用电子－空穴理论推算汶川ＭＳ
８．０地震电磁异常的“能量源”大小［Ｊ］．地震学报，２０１５（５）：８４２－
８５２．［１３］　高建国．中国最大地震的地声研究［Ｊ］．应用声学，１９８６，５（４）：１８－
２３．［１４］　Ｈｉｌｌ　Ｄ　Ｐ，Ｆｉｓｃｈｅｒ　Ｆ　Ｇ，Ｌａｈｒ　Ｋ　Ｍ，ｅｔ　ａｌ．Ｅａｒｔｈｑｕａｋｅ　ｓｏｕｎｄｓ　ｇｅｎｅｒａｔｅｄ　ｂｙ　
ｂｏｄｙ－ｗａｖｅ　ｇｒｏｕｎｄ　ｍｏｔｉｏｎ［Ｊ］．Ｂｕｌｌ．Ｓｅｉｓｍｏｌ．Ｓｏｃ．Ａｍ．；（Ｕｎｉｔｅｄ　Ｓｔａｔｅｓ），１９７６，６６（４）：１　１５９－１　
１７２．［１５］　Ｘｉａ　Ｙ，Ｌｉｕ　Ｊ　Ｙ　Ｔ，Ｃｕｉ　Ｘ，ｅｔ　ａｌ．Ａｂｎｏｒｍａｌ　ｉｎｆｒａｓｏｕｎｄ　ｓｉｇｎａｌｓ　ｂｅｆｏｒｅ　９２　Ｍ≥７．０ｗｏｒｌｄｗｉｄｅ　ｅａｒｔｈ－ｑｕａｋｅｓ　ｄｕｒｉｎｇ　２００２—２００８［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ａｓｉａｎ　Ｅａｒｔｈ　Ｓｃｉｅｎｃｅｓ，２０１１，４１（４）：４３４－４４１．［１６］　庞瑞涛，雍珊珊，王新安，等．地震监测系统的电磁信号的采集设计与实现［
Ｊ］．计算机技术与发展，２０１８，２８（２）：２７－
３０．［１７］　雍珊珊，王新安，庞瑞涛，等．多分量地震监测系统Ａ
ＥＴＡ的感应式磁传感器磁棒研制［Ｊ］．北京大学学报（自然科学版），２０１８，５４（３）：４９５－
５０１．［１８］　王新安，雍珊珊，徐伯星，等．多分量地震监测系统Ａ
ＥＴＡ的研究与实现［Ｊ］．北京大学学报：自然科学版，２０１８，５４（３）：４８７－
４９４．［１９］　林科，王新安，张兴，等．一种适用于大地震临震预测的地声监测系统［
Ｊ］．华南地震，２０１３，３３（４）：５４－
６２．［２０］　Ｐａｒｒｏｔ　Ｍ，Ｂｅｎｏｉｓｔ　Ｄ，Ｂｅｒｔｈｅｌｉｅｒ　Ｊ　Ｊ，ｅｔ　ａｌ．Ｔｈｅ　ｍａｇｎｅｔｉｃ　ｆｉｅｌｄ　ｅｘｐ
ｅｒｉｍｅｎｔ　ＩＭＳＣ　ａｎｄ　ｉｔｓ　ｄａｔａ　ｐｒｏｃｅｓｓ－ｉｎｇ　ｏｎｂｏａｒｄ　ＤＥＭＥＴＥＲ：Ｓｃｉｅｎｔｉｆｉｃ　ｏｂｊｅｃｔｉｖｅｓ，ｄｅｓｃｒｉｐｔｉｏｎ　ａｎｄ　ｆｉｒｓｔ　ｒｅｓｕｌｔｓ［Ｊ］．Ｐｌａｎｅｔａｒｙ　
＆ＳｐａｃｅＳｃｉｅｎｃｅ，２００６，５４（５）：４４１－４５５．［２１］　汤吉，赵国泽，陈小斌，等．地震电磁卫星载荷及现状［Ｊ］．地球物理学进展，２００７，２２（３）：６７９－
６８６．
［２２］　Ｃａｏ　Ｊ　Ｂ，Ｚｅｎｇ　
Ｌ，Ｚｈａｎ　Ｆ，ｅｔ　ａｌ．Ｔｈｅ　ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ　ｗａｖｅ　ｅｘｐｅｒｉｍｅｎｔ　ｆｏｒ　ＣＳＥＳ　ｍｉｓｓｉｏｎ：Ｓｅａｒｃｈ　ｃｏｉｌｍａｇ
ｎｅｔｏｍｅｔｅｒ［Ｊ］．Ｓｃｉｅｎｃｅ　Ｃｈｉｎａ　Ｔｅｃｈｎｏｌｏｇｉｃａｌ　Ｓｃｉｅｎｃｅｓ，２０１８，６１（５）：１－６．［２３］　赵洪利．
ＰＶＤＦ压电薄膜传感器的研究及其在ＰＣＢＡ模态分析中的应用［Ｄ］．镇江：江苏大学，２００９．
［２４］　卓翼科技．关于卓翼［ＥＢ／ＯＬ］．［２０１８－１１－２９］．ｈｔｔｐ：∥ｗｗｗ．ｚｏｗｅｅ．ｃｏｍ．ｃｎ／ａｂｏｕｔ．ａｓｐ
ｘ．［２５］　百度百科．深圳卓翼科技股份有限公司［ＥＢ／ＯＬ］．［２０１８－１１－２９］．ｈｔｔｐ
ｓ：∥ｂａｉｋｅ．ｂａｉｄｕ．ｃｏｍ／ｉｔｅｍ／深圳卓翼科技股份有限公司／８８１４４９９？ｆｒ＝ａｌａｄｄｉｎ．
［２６］　Ｇｕｏ　Ｊ，Ｌｉ　Ｗ，Ｙｕ　Ｈ，ｅｔ　ａｌ．Ｉｍｐｅｎｄｉｎｇ　ｉｏｎｏｓｐｈｅｒｉｃ　ａｎｏｍａｌｙ　ｐｒｅｃｅｄｉｎｇ　
ｔｈｅ　Ｉｑｕｉｑｕｅ　ＭＷ８．２ｅａｒｔｈｑｕａｋｅｉｎ　Ｃｈｉｌｅ　ｏｎ　２０１４Ａｐｒｉｌ　１［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ　Ｊｏｕｒｎａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ，２０１５，２０３（３）：１４６１－１４７０．［２７］　Ｌｉｕ　Ｊ　Ｙ，Ｃｈｅｎ　Ｙ　Ｉ，Ｃｈｅｎ　Ｃ　Ｈ，ｅｔ　ａｌ．Ｓｅｉｓｍｏｉｏｎｏｓｐ
ｈｅｒｉｃ　ＧＰＳ　ｔｏｔａｌ　ｅｌｅｃｔｒｏｎ　ｃｏｎｔｅｎｔ　ａｎｏｍａｌｉｅｓ　ｏｂ－ｓｅｒｖｅｄ　ｂｅｆｏｒｅ　ｔｈｅ　１２Ｍａｙ　２００８　ＭＷ７．９Ｗｅｎｃｈｕａｎ　ｅａｒｔｈｑｕａｋｅ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｇｅｏｐｈｙｓｉｃａｌ　ＲｅｓｅａｒｃｈＳｐａｃｅ　Ｐｈｙ
ｓｉｃｓ，２００９，１１４（Ａ４）：２３１－２６１．２
９ 地 震 ３
９卷　
［２８］　Ｎｅｎｏｖｓｋｉ　Ｐ　Ｉ，Ｐｅｚｚｏｐａｎｅ　Ｍ，Ｃｉｒａｏｌｏ　Ｌ，ｅｔ　ａｌ．Ｌｏｃａｌ　ｃｈａｎｇｅｓ　ｉｎ　ｔｈｅ　ｔｏｔａｌ　ｅｌｅｃｔｒｏｎ　ｃｏｎｔｅｎｔ　ｉｍｍｅｄｉａｔｅｌｙ
ｂｅｆｏｒｅ　ｔｈｅ　２００９Ａｂｒｕｚｚｏ　ｅａｒｔｈｑｕａｋｅ［Ｊ］．Ａｄｖａｎｃｅｓ　ｉｎ　Ｓｐ
ａｃｅ　Ｒｅｓｅａｒｃｈ，２０１５，５５（１）：２４３－２５８．［２９］　Ｈａｔｔｏｒｉ　Ｋ．ＵＬＦ　ｇｅｏｍａｇ
ｎｅｔｉｃ　ｃｈａｎｇｅｓ　ａｓｓｏｃｉａｔｅｄ　ｗｉｔｈ　ｌａｒｇｅ　ｅａｒｔｈｑｕａｋｅｓ［Ｊ］．Ｔｅｒｒｅｓｔｒｉａｌ　Ａｔｍｏｓｐｈｅｒｉｃ＆Ｏｃｅａｎｉｃ　Ｓｃｉｅｎｃｅｓ，２００４，１５（３）：３２９－
３６０．［３０］　马钦忠，尹京苑，顾学章．崇明异常电磁扰动与台湾７．５级强震［Ｊ］．地震，２００３，２３（４）：４９－５６．［３１］　徐文耀．地球电磁现象物理学［
Ｍ］．合肥：中国科学技术大学出版社，２００９．［３２］　Ｒｉｃｈ　Ｆ　Ｊ，Ｓｕｌｔａｎ　Ｐ　Ｊ，Ｂｕｒｋｅ　Ｗ　Ｊ．Ｔｈｅ　２７－ｄａｙ　
ｖａｒｉａｔｉｏｎｓ　ｏｆ　ｐｌａｓｍａ　ｄｅｎｓｉｔｉｅｓ　ａｎｄ　ｔｅｍｐｅｒａｔｕｒｅｓ　ｉｎ　ｔｈｅｔｏｐｓｉｄｅ　ｉｏｎｏｓｐｈｅｒｅ［Ｊ］．Ｊｏｕｒｎａｌ　ｏｆ　Ｇｅｏｐｈｙｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｓｐａｃｅ　Ｐｈｙｓｉｃｓ，２００３，１０８（Ａ７）：１　２７９－１　
３０４．［３３］　张培震，邓起东，张国民，等．中国大陆的强震活动与活动地块［
Ｊ］．中国科学：地球科学，２００３，３３（ｓ１）：１２－
２０．［３４］　邓起东，张培震，冉勇康，等．中国活动构造与地震活动［Ｊ］．地学前缘，２００３，１０（ｓ１）：６６－７３．［３５］　Ｒｅｎ　Ｈ，Ｗｅｎ　Ｊ，Ｈｕａｎｇ　
Ｑ，ｅｔ　ａｌ．Ｅｌｅｃｔｒｏｋｉｎｅｔｉｃ　ｅｆｆｅｃｔ　ｃｏｍｂｉｎｅｄ　ｗｉｔｈ　ｓｕｒｆａｃｅ－ｃｈａｒｇｅ　ａｓｓｕｍｐｔｉｏｎ：ａｐ
ｏｓｓｉｂｌｅ　ｇｅｎｅｒａｔｉｏｎ　ｍｅｃｈａｎｉｓｍ　ｏｆ　ｃｏｓｅｉｓｍｉｃ　ＥＭ　ｓｉｇｎａｌｓ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ　Ｊｏｕｒｎａｌ　Ｉｎｔｅｒｎａｔｉｏｎａｌ，２０１５，２００（２）：８３５－
８４８．［３６］　Ｍｏｌｃｈａｎｏｖ　Ｏ　Ａ，Ｈａｙａｋａｗａ　Ｍ．Ｇｅｎｅｒａｔｉｏｎ　ｏｆ　ＵＬＦ　ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ　ｅｍｉｓｓｉｏｎｓ　ｂｙ　
ｍｉｃｒｏｆｒａｃｔｕｒｉｎｇ［Ｊ］．Ｇｅｏｐｈｙ
ｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｌｅｔｔｅｒｓ，２０１３，２２（２２）：３　０９１－３　０９４．［３７］　Ｓｉｍｐｓｏｎ　Ｊ　Ｊ，Ｔａｆｌｏｖｅ　Ａ．Ｅｌｅｃｔｒｏｋｉｎｅｔｉｃ　ｅｆｆｅｃｔ　ｏｆ　ｔｈｅ　Ｌｏｍａ　Ｐｒｉｅｔａ　ｅａｒｔｈｑｕａｋｅ　ｃａｌｃｕｌａｔｅｄ　ｂｙ　
ａｎ　ｅｎｔｉｒｅ－Ｅａｒｔｈ　ＦＤＴＤ　ｓｏｌｕｔｉｏｎ　ｏｆ　Ｍａｘｗｅｌｌ　ｅｑｕａｔｉｏｎｓ［Ｊ］．Ｇｅｏｐｈｙ
ｓｉｃａｌ　Ｒｅｓｅａｒｃｈ　Ｌｅｔｔｅｒｓ，２００５，３２（９）：Ｌ９３０２．［３８］　Ｍｅｒｚｅｒ　Ｍ，Ｋｌｅｍｐｅｒｅｒ　Ｓ　Ｌ．Ｍｏｄｅｌｉｎｇ　Ｌｏｗ－ｆｒｅｑｕｅｎｃｙ　
Ｍａｇｎｅｔｉｃ－ｆｉｅｌｄ　Ｐｒｅｃｕｒｓｏｒｓ　ｔｏ　ｔｈｅ　Ｌｏｍａ　ＰｒｉｅｔａＥａｒｔｈｑｕａｋｅ　ｗｉｔｈ　ａ　Ｐｒｅｃｕｒｓｏｒｙ　Ｉｎｃｒｅａｓｅ　ｉｎ　Ｆａｕｌｔ－ｚｏｎｅ　Ｃｏｎｄｕｃｔｉｖｉｔｙ［Ｊ］．Ｐｕｒｅ　＆Ａｐｐｌｉｅｄ　Ｇｅｏｐｈｙｓｉｃｓ，１９９７，１５０（２）：２１７－
２４８．Ａｎａｌｙｓｉｓ　ｏｆ　Ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ　Ｄｉｓｔｕｒｂａｎｃｅ　Ａｎｏｍａｌｙ　
ｐｒｉｏｒｔｏ　ｔｈｅ　２０１７　ＭＳ７．０Ｊｉｕｚｈａｉｇｏｕ　Ｅａｒｔｈｑ
ｕａｋｅＨＵＡＮＧ　Ｊｉ－ｐａｎ１，２，
３，ＷＡＮＧ　Ｘｉｎ－ａｎ１，２，
３，ＹＯＮＧ　Ｓｈａｎ－ｓｈａｎ１，
２，ＷＡＮＧ　Ｍｏ１，
３，ＬＩＵ　Ｚｈｏｎｇ
１，
３
（１．Ｔｈｅ　Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｆ　Ｉｎｔｅｇｒａｔｅｄ　Ｍｉｃｒｏ－ｓｙｓｔｅｍｓ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｅｎｇｉｎｅｅｒｉｎｇ　Ａｐｐｌｉｃａｔｉｏｎｓ，Ｐｅｋｉｎｇ　ＵｎｉｖｅｒｓｉｔｙＳｈｅｎｚｈｅｎ　Ｇｒａｄｕａｔｅ　Ｓｃｈｏｏｌ，Ｇｕａｎｇｄｏｎｇ　Ｓｈｅｎｚｈｅｎ　５１８０５５，Ｃｈｉｎａ；２．Ｔｈｅ　Ｒｅｓｅａｒｃｈ　Ｃｅｎｔｅｒ　ｏｆ　ＥａｒｔｈｑｕａｋｅＭｏｎｉｔｏｒｉｎｇ　ａｎｄ　Ｆｏｒｅｃａｓｔｉｎｇ　Ｔｅｃｈｎｏｌｏｇｙ，Ｐｅｋｉｎｇ　Ｕｎｉｖｅｒｓｉｔｙ　
Ｓｈｅｎｚｈｅｎ　Ｇｒａｄｕａｔｅ　Ｓｃｈｏｏｌ，ＧｕａｎｇｄｏｎｇＳｈｅｎｚｈｅｎ　５１８０５５，Ｃｈｉｎａ；３．Ｓｃｈｏｏｌ　ｏｆ　Ｅｌｅｃｔｒｏｎｉｃｓ　Ｅｎｇｉｎｅｅｒｉｎｇ　
ａｎｄ　Ｃｏｍｐｕｔｅｒ　Ｓｃｉｅｎｃｅ，Ｐｅｋｉｎｇ　Ｕｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ　
１００８７１，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ　ｄｉｓｔｕｒｂａｎｃｅ　ｐｈｅｎｏｍｅｎａ　ｂｅｆｏｒｅ　ｅａｒｔｈｑｕａｋｅｓ　ｈａｖｅ　ｂｅｅｎ　ｏｂｓｅｒｖｅｄｆｏｒ　ｍａｎｙ　ｔｉｍｅｓ，ｈｏｗｅｖｅｒ　ｔｈｅ　ｃｏｒｒｅｌａｔｉｏｎ　ｓｈｏｕｌｄ　ｂｅ　ｆｕｒｔｈｅｒ　ｓｔｕｄｉｅｄ．Ｗｅ　ｄｅｓｉｇｎｅｄ　ａ　ｍｕｌｔｉ－ｃｏｍｐｏｎｅｎｔ　ｓｅｉｓｍｉｃ　ｍｏｎｉｔｏｒｉｎｇ　ｓｙｓｔｅｍ－ＡＥＴＡ，ｔｏ　ｍｏｎｉｔｏｒ　ｔｈｅ　ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ　ｄｉｓｔｕｒｂａｎｃｅ．Ａｔ　ｐｒｅｓｅｎｔ　ｍｏｒｅ　ｔｈａｎ　２００ＡＥＴＡ　ｄｅｖｉｃｅｓ　ｈａｖｅ　ｂｅｅｎ　ｄｅｐｌｏｙｅｄ　ｉｎ　Ｃｈｉｎａ　ｕｎｄｅｒ　ｔｈｅ　ｓｕｐｐ
ｏｒｔ　ｏｆｔｈｅ　Ｃｈｉｎｅｓｅ　Ｅａｒｔｈｑｕａｋｅ　Ａｄｍｉｎｉｓｔｒａｔｉｏｎ．Ｏｎ　Ａｕｇｕｓｔ　８，２０１７，ａ　ｌａｒｇｅ　ＭＳ７．０ｅａｒｔｈｑｕａｋｅｏｃｃｕｒｒｅｄ　ｉｎ　Ｊｉｕｚｈａｉｇｏｕ　ｃｏｕｎｔｙ　ｏｆ　Ｓｉｃｈｕａｎ　Ｐｒｏｖｉｎｃｅ，Ｃｈｉｎａ．Ｔｈｅｒｅ　ａｒｅ　２２ｓｅｔｓ　ｏｆ　ｄｅｖｉｃｅｓｉｎｓｔａｌｌｅｄ　ａｔ　ｓｅｉｓｍｉｃ　ｓｔａｔｉｏｎｓ　ｗｉｔｈｉｎ　ｔｈｅ　５００ｋｍ　ｒａｎｇｅ　ｏｆ　ｔｈｅ　ｅｐｉｃｅｎｔｅｒ．Ｗｅ　ａｎａｌｙ
ｚｅｄ　ｔｈｅ３
９　３期 黄继攀等：
２０１７年九寨沟７级地震前电磁扰动异常分析　
ｏｂｓｅｒｖａｔｉｏｎ　ｄａｔａ　ｏｆ　３０ｄａｙｓ　ｂｅｆｏｒｅ　ａｎｄ　ａｆｔｅｒ　ｔｈｅ　ｅａｒｔｈｑｕａｋｅ　ｕｓｉｎｇ　ｔｈｅ　ｓｌｉｄｉｎｇ　
ｉｎｔｅｒｑｕａｒｔｉｌｅｍｅｔｈｏｄ　ａｎｄ　ａ　２７－ｄａｙ　
ｗｉｎｄｏｗ．Ｔｈｅ　ｒｅｓｕｌｔ　ｓｈｏｗｓ　ｔｈａｔ　ａｂｎｏｒｍａｌ　ｓｉｇｎａｌ　ａｐｐｅａｒｅｄ　ａｔ　１１ｓｔａ－ｔｉｏｎｓ．Ａｍｏｎｇ　
ｔｈｅｍ，４ｏｕｔ　ｏｆ　６ｓｔａｔｉｏｎｓ　ｗｉｔｈｉｎ　２００ｋｍ　ａｒｏｕｎｄ　ｔｈｅ　ｅｐｉｃｅｎｔｅｒ　ｈａｖｅ　ａｎｏｍａ－ｌｉｅｓ．Ａｎｄ　ｔｈｅ　ｐｏｓｉｔｉｖｅ　ａｎｏｍａｌｙ　
ｉｓ　ｍｕｃｈ　ｍｏｒｅ　ｔｈａｎ　ｔｈｅ　ｎｅｇａｔｉｖｅ　ａｎｏｍａｌｙ．Ｂｅｓｉｄｅｓ，ｍｏｓｔａｎｏｍａｌｉｅｓ　ｏｃｃｕｒｒｅｄ　ｄｕｒｉｎｇ　
ｎｉｇｈｔ　ｔｉｍｅ（２０：００～６：００）．Ｋｅｙ　
ｗｏｒｄｓ：　Ｔｈｅ　２０１７Ｊｉｕｚｈａｉｇｏｕ　Ｅａｒｔｈｑｕａｋｅ；Ｅｌｅｃｔｒｏｍａｇｎｅｔｉｃ　Ｄｉｓｔｕｒｂａｎｃｅ；Ｔｉｍｅ－ｓｅｒｉｅｓａｎａｌｙｓｉｓ；Ａｎｏｍａｌｙ４９ 地 震 ３
９卷　。