基金项目:上海市科学技术委员会科研计划项目(16dz1205800,23DZ1200200)
作者简介:马钦忠(1959-),男,研究员,主要从事地震电磁学研究。E-mail: mqz1234@sina.com
1.上海市地震局,上海 200062;2.上海佘山地球物理国家野外科学观测站,上海 200062;3.福建省地震局漳州地震监测中心站,福建 漳州 363118
1.Shanghai Earthquake Agency, Shanghai 200062, China;2.Shanghai Sheshan National Geophysical Observatory, Shanghai 200062, China;3.Zhangzhou Earthquake Monitoring Center Station, Fujian Earthquake Agency,Zhangzhou 363118, China
Fujian-Taiwan region; Geoelectric field observation; Fingerprint information method; Earthquake prediction
DOI: 10.13512/j.hndz.2026.04.04
闽台地区、特别是台湾属于地震多发地带之一。台湾地处东面的琉球弧和南面的吕宋弧两个沟—弧体系的结合带,复杂的构造环境以及菲律宾亚板块与欧亚大陆的相互作用,使台湾成为全球研究弧—陆碰撞最受关注的地区[1]。因此该区域的地震频发,是全球地震活动最为活跃的地区之一。在研究地震预测特别是短临预测时所面临的主要难点问题之一就是地震是小概率事件,所以要检验地震地电场用于地震短临预测的有效性最好方式就是选择地震多发区域进行连续观测和探索预测方法。而我国闽台地区就是很好的实验观测区域。在台湾的东海岸,中强地震和大地震出现的频率非常高,在岛内区域地震同样活跃,据有记录以来,花莲海域周边300 km已经发生过49次7级以上地震。1999年发生在台湾中部的7.7级集集地震造成了非常严重的灾难。所以选择地震监测预报研究的区域和用于预报地震活动的监控方法是非常重要的。
利用地震前地电场异常信息进行地震三要素(时间、地点、强度)的预测,在国内外许多研究文献中已有报道,且积累了大量的研究震例[2-21]。1966年邢台M7.2地震后,我国开展了地电场观测研究,观测技术从早期的模拟观测(土地电),到后来的数字化地电场观测,记录到了典型的震例[22]。我国地电场观测研究始于1966年3月22日我国邢台发生M7.2地震后,自此无论是先期的模拟观测(土地电)还是后来的数字化地电场观测研究都已经积累了许多震例[22]。上世纪90年代我国研制了ZD9型观测频带范围为DC-0.005 Hz的数字化大地电场仪[23],固体不极化电极[24-25],这些技术应用极大的提高了观测数据的质量[24, 26]。截止目前已经建设120多个地电场观测台站,其中电极敷设方式:长极距在300~400 m范围、短极距在150~250 m的长度,分别在南北、东西及北东(北西)方向共布设了六个测道的测量装置,分别在三个方向上都有长短不一的两个测道,构成了地电场多极距观测系统,在相同方向上两个测道比测的方式实现排除电极极化等影响,为剔除与地震无关的噪声并提取地震电信号提供判定依据[15, 19, 27]。我们近年开展的针对地电场点阵观测系统的研究(仪器微型化、低功耗、低成本化、数据无线传输、电极布设的小极距化),使得操作方便易行的小极距地电场观测取得了较好的效果,为野外临时观测提供了极大的方便[22]。
在福建省地震局“十五”项目建设期间在泉州市地震台建立了泉州地电场观测,但只有一个地电场台站,还是太少,而且也没有ULF频段的电磁扰动观测。为了能够验证福建沿海地震地电场观测的有效性,我们分别在福建漳洲台、东山台、平潭台建立了临时性的地电场观测,与泉州地电场台站一起构成了福建沿海地电场观测网,以期检验该观测台网对闽台地区地震前的异常信息观测的有效性。本文是对这方面工作的初步总结,以期为今后在本地区开展地震电磁观测提供十分有意义的参考。
目前中国地震系统敷设地电场观测仪器是ZD9系列及后期的GEF系列大地电场仪[23,27-29]。建台中,电极布设方式按照“多方向、多极距”观测方法布设测量线路,即在NS、EW、NE或NW方向上布设2个以上极距的测量线路,长极距长度为300~400 m之间,短极距在150~200 m之间,构成地电场的多极距观测系统。图1(c)为在漳州台使用的ZD9仪器。
图1 福建沿海地电场观测台站示意图Fig.1 Schematic diagram of the geoelectric field observatries along Fujian coast
为了开展地震地电场的阵列观测或点阵观测以及流动观测活动,我们进行了地电场仪微型化改造,形成了微型地电场仪,探索了微型化、低功耗、高采样率、数据无线传输、低成本、高密度的流动性地电场观测技术问题,2016年在上海市科委项目《短临地震地电观测新技术研究》支持下本项目研制地电场点阵观测装置,将该仪器命名为“DCDZ-01地电场点阵观测系统”,实现了小极距地电场观测目标。该观测仪的测量分辨力优于1 μV,频带范围为DC-0.005 Hz,采样率为16 Hz,每天输出一次数据文件到服务器,数据可以用手机或无线联网的PC接收;观测数据采集后通过物联网2G网络直接发送至云端,设备通过网络每日自动校时,实现了数据采集到存储的基本观测功能。野外流动观测时用锂电池供电可以持续观测30天;电池电量报警、接收开路等报警信息可在手机端收到短信通知。
由于闽台地区属于地震多发地区,尤其是台海地区地震频发,为了观测并获得更多的震前地电场异常信号,我们分别在漳州地磁台、东山地震台、平潭地震台内建立了临时地电场观测,与泉州台一起构成了面对台海地区的福建省沿海较为均匀分布的地电场观测站。
福建省地震局在泉州已经建立了一套ZD9大地电场仪观测系统,该系统在NS、EW、NW方向都有两个测道,NS和EW向长极距为300 m、短极距为150 m,NW向长极距和短极距长度分别为424 m和212 m。泉州地电场台台址位于南安市洪濑镇厝斗产田村。观测区坐标:北纬 25.02°、东经118.51°,高程约41 m。台址坐落在葵山脚下一片较为开阔的平地边缘,其南北均为低矮的丘陵台地,西面为开阔的冲击平原,东面为狭长平地,直达葵山脚下。该台址位于闽东滨海断隆带中段西侧,是北北东—北东向长乐—诏安断裂带与北西向永安—晋江断裂带的交汇地带。台址处于燕山期花岗岩为底之上的较厚洪积、淤积地层,周边石英岩脉成片出露,场址内地层均匀分布,石英砂层厚达6~7 m,石英砂砾层厚达10 m。
我们在2017年4月27日—2022年12月14日期间使用ZD9大地电场仪在漳州地震台建立了地电场观测,其电极布设方式如图1中间第三图所示,数据传输通过该台内网传回。该台位于龙海市九湖镇衍后村石狮岩,地理坐标:东经117.37°、北纬24. 28°,海拔42.3 m。台站基岩为花岗闪长岩。漳州台周边地区地质构造复杂,台站位于几条断裂带交汇合围形成的区域内,距东南一侧北东走向的长乐—诏安断裂带仅几十公里,西北部和西南部分别有北东向政和—海丰断裂带和北西向上杭—云霄断裂带经过,东北部则是九龙江断裂。
在2018年8月27日—2021年10月1日3期间在福建省东山地震台使用微型地电场仪建立了地电场观测,其电极布设如图1(b)所示,布极区在该台站办公楼前的花草园地里,电极距都是小极距布设,极距长度都在20 m内,因此操作简单,数据传输通过2G网络无线传回。
2019年9月12日—2021年7月4日期间在福建省平潭地震台使用微型地电场仪建立了地电场观测,其电极布设如图1中间第2图所示,在不同方向极距长度分别为15 m、16 m、20 m,布极区在该台站台地上面的草地里。观测数据通过2G网络无线传输。
动态跟踪每个台站地震电信号及其变化特征。具体震前地电场“指纹”信息跟踪原理如下:
第一步:收集每个台站观测到的各种异常信号。
第二步:落实和识别源自各种源的异常信号:
(a)电极极化噪声(由于电极性能的改变或埋设条件不良而产生);
(b)附近的电磁干扰噪声,如工业游散电流流过测区在测量电极上引起的附加电信号;
(c)台站附近地下介质电性变化导致地电场的局部畸变引起的噪声(降雨和温度变化引起的);
(d)雷电噪声;
(e)电离层电磁扰动噪声(该类噪声特点是其引起的地电暴具有广域性,即在全国地电场台站基本上都能在同时段观测到,因此通过多台和远距离参考台对比观测可识别;也可参考地磁指数Kp≥5.0的条件加以识别);
(f)高压输电接地极入地大电流信号。
第三步:分类和归类源自各种信号源的异常信号集合。
第四步:研究源自不同信号源的每种异常信号的特点。
第五步:在上述工作基础上挑选出疑似与地震有关的异常信号。
第六步:用实际发生地震信息对所选出的疑似与地震有关的异常信号进行比对分析和确认。
第七步:经过比对分析结果后,针对新出现的、特征能够相吻合的确定性地震电信号进行更深入研究,以获得对未来地震三要素的判定。
上述工作的基础就是积累并形成各个台站地电场异常信息库,该库收集了上述(a)、(b)、(c)、(d)、(e)、(f)类的噪声信号以及除此之外的异常信号,并在此基础上进行模式识别处理,模式识别是指对表征事物或现象的各种形式的信息进行处理和分析,以对事物或现象进行描述、辨认、分类和解释的过程。
通过上述对地电场观测数据的分析处理,可以获得信度高的地震电信号,并依此判定未来地震震级、发震时间和地点三要素,信息越可靠、分析越细致则分析结果的误差越小。关于判定引起地电场变化是远场、还是近场的基本原理已有一些研究[19,30]。本文下节给出一些基于福建沿海地电场台站特征信息库对比分析所得到的震前地震电信号特征图,由于篇幅所限,许多特征图在此省略。
该台在观测运行期间出现几次数据传输中断现象(服务器故障)。故所选出的地震电信号并没有完全覆盖所有地震。
图2所示的是东山台在2018年11月1日所观测到的地震电信号,由该图可以看到,该异常信号的形态和后面图3~图5的信号形态有很大的差异性,持续时间长,持时达8个多小时。25天之后于2018年11月26日在台湾海峡发生了M6.2地震,该地震震中位置距离东山观测点只有约126 km,并在福建省引起了强烈震感。
图2 2018年11月1日东山台地震电信号特征图Fig.2 Characteristic of seismic electric signals at Dongshan Station on November 1, 2018.
图3所示的是2019年3月8日14:30—15:30期间观测到的地震电信号。在41天后台湾花莲县海边发生了M6.7地震。
图4所示的是东山台于2020年5月11日观测到的地震电信号特征图,该异常信号出现的时间是夜里临晨00:40—02:00期间,异常信号出现后的第34天在台湾宜兰县海域发生了M5.5地震。
漳州台地电场观测时段在2017年4月27日—2022年12月14日,期间出现多次故障,包括观测场地因台站改造而中断、仪器故障而维修和更换、网线出现故障等。因此,观测数据不完整。2018年10月1日至11月23日因故观测中断,11月24日才恢复观测,所以很遗憾在此期间不知道11月26日台湾海峡发生M6.2地震之前的这段时间异常信号的情况。图5所示的是2018年11月25日漳州台恢复观测后观测到的疑似地震电信号特征图,该异常信号的特点是在NS方向和EW方向的短极距上观测到同步信号,具体持续时间是17:00开始大幅度下降,到17:30下降幅度达到最大值,然后在17:30—18:00期间快速恢复,然后继续逐渐恢复直到24:00为止。在其它测道上没有出现,参考图1中的漳州台电极布设图可以看到,不可能是中心电极区域电极极化噪声,因为同测向长极距没有出现该信号。第二天,既2018年11月26日在距离漳州台160 km处的台湾海峡发生了M6.2地震,福建省及广东省东部震感强烈。
图3 2019年3月8日东山台地震电信号特征图Fig.3 Characteristic of seismic electric signals at Dongshan Station on March 8, 2019.
图4 2020年5月11日东山台地震电信号特征图Fig.4 Characteristic of seismic electric signals at Dongshan Station on May 11, 2020.
图5 2018年11月25日漳州台秒种值疑似地震电信号特征图Fig.5 Characteristic of suspected seismic electric signals of second-sampling values at Zhangzhou Station on November 25, 2018.
图6所示的是2021年4月2日漳州台地震电信号特征图,该信号出现在夜里临晨01:50—02:10期间,由该图可以看到依据多极距原理可排除电极极化噪声问题,而夜里又是噪声水平最低的时间。16天之后,即2021年4月18日在台湾花莲县相继发生了M5.6地震和M6.1地震,震中距约390 km。
图6 2021年4月2日漳州台地震电信号特征图Fig.6 Characteristic of seismic electric signals at Zhangzhou Station on April 2, 2021.
该台属于福建省地震局建设,正常观测运行情况良好,由于该台地电场观测在2018年前后进行过一次搬迁改造,因此搬迁前后所观测到的地震电信号特征变化也较大。图7所示的是2018年11月3日泉州台地电场特征图,其特征是当日20:30前后在6个测道上都观测到了异常信号,持续时间约10 min左右。23天之后,既2018年11月26日台湾海峡发生了M6.2地震,震中距约190 km。
图8所示的是2024年2月29日—3月17日期间泉州台地电场特征图,分别在2024年2月29日13:00—19:00期间、3月10日10:00—3月11日19:00期间和3月17日00:00—22:00期间在NS方向长极距上观测到了异常信号,这些异常信号持续时间分别为:5小时30分钟、22小时和31小时50分钟。从观测系统特点来分析,首先排除电极极化噪声的因素,因为NS方向长极距两端的任何一个电极极化的话都必须在六个测道曲线中出现两条与同一电极相关的含噪声曲线,但这一情况并未出现;其次,也不是仪器故障,因为EW、NE方向长、短极距测道以及NS方向短极距测道都观测正常;再则,也排除雷电噪声和高压输电噪声,因为另外五个测道上都没有观测到该异常;也排除地电暴噪声,因为该信号不具备广域性。3月17日异常电信号结束后的第17天,即4月3日7时58分在台湾花莲县海域(23.81°N,121.74°E)发生7.3级地震,震源深度12 km。因此,我们认为图9所示的异常是由该地震在孕震过程中产生的地震电信号。根据以往泉州地电场台搬迁之前的资料分析,台湾发生中强地震前泉州台所观测到的地震电信号在六个测道上都可观测到,但该台站在搬迁后在此次7.3级大地震前只在NS方向的长极距上观测到了很显著的地震电信号而在其余五个测道上没有记录到该信号,推断这是由于台站台址下方地下电性结构的复杂性所致。图8所示4月份出现的异常是干扰,这可由地电场多极距观测原理所证明。
图7 2018年11月3日20:00—21:00期间泉州台地电场特征图Fig.7 Characteristic of seismic electric signals at Quanzhou Station from 20:00 to 21:00 on November 3, 2018.
图8 2024年2月29日—3月17日期间泉州台地电场特征图Fig.8 Characteristic of seismic electric signals at Quanzhou Station from November 3 to March 17, 2018
由于经费所限,我们在该台开展的地电场观测时段最短,观测时间为2019年9月12日—2021年7月4日。图9上图是平潭台于2020年11月29日至12月3日期间所观测到的地震电信号,该异常信号持续时间较长,最大幅度93 mV,7天之后,既2020年12月10日在台湾宜兰县海域发生了M5.8地震,震中距约为350 km。图9中图所示的是平潭台在2021年4月8日—2021年4月9日期间观测到的地震电信号,最大幅度出现在9日临晨01:00前后,持续时间约12个小时。9天之后,既2020年4月18日在台湾花莲县先后发生M6.1地震和M5.6地震。图9下图是平潭台在2020年2月19—20日观测到的异常变化,最大幅度19 mV,5天后台湾宜兰县发生了4.9级地震。
图2~图9显示了一小部分福建沿海4个地电场台站观测到的与闽台地区地震相关性的异常信号,表1列出了部分福建沿海地区东山台、漳州台、泉州台和平潭台地电场观测到的一些异常信号与闽台地区地震关系,还有许多这样的关系因篇幅所限在此略去。
图9 平潭台于2020年11月29日—12月3日(a)、2021年4月8—9日(b)、2020年2月19—20日(c)观测到的地震电信号特征图Fig.9 Characteristic of seismic electric signals at Pingtan station from November 29 to December 3, 2020(a), from April 8 to 9, 2021(b), from February 19 to 20, 2020(c)
表1 福建沿海地电场观测一些异常信号与闽台地区地震关系表Table 1 Relationship between anomalous geoelectric signals observed at Fujian coastal stations and earthquakes in Fujian-Taiwan region

由上节可以看到近些年在台湾海峡区域发生的最大地震是2018年11月26日发生的M6.2地震,而在此地震发生前在福建沿海东山台、漳州台、泉州台都提前观测到了地电场异常信号,东山台观测到的时间是2018年11月1日,泉州台观测到的时间是2018年11月3日,漳州台观测到的时间是2018年11月25日(由于故障从10月1日至11月23日期间最关键的时段观测数据丢失)。从国内外大量研究震例来看,震前观测到的电磁异常信号很少是各个台站同时观测到,都是在不同的时间观测到的。这是因为孕震体是一个区域,当应力积累过程是发生在不同的应力积累区或凹凸体时,应力积累加速区就会产生并传播电磁信号[6,31],这种信号因为地下介质电性非均匀性和各向异性的影响从而使得在地表有些观测点能够被观测到,而另一些观测点观测不到这样的信号,这就是所谓的“敏感点效应”[32-33]。因此,对于2018年11月26日台湾海峡M6.2地震,东山台、漳州台、泉州台观测到的异常信号时间不一致,可以认为是因为信号源或凹凸体不一致导致的。由图2~图4可以看到,在东山台观测到的地震电信号波形形态上来看,海峡M6.2地震前地震电信号形态与台湾地区地震电信号的波形形态上存在着较大差异,前者呈现连续下降然后上升恢复到正常的水平,而后者呈现出尖峰脉冲组的形态,这种特征差异性应该可以作为短临地点预测的参考因素。
在本文所总结的地震电信号的可信度上,如果能够通过实践的检验则可以大大提高可信度。在这方面的检验,我们大都通过检测出地震电信号,做出判断等待地震的发生,然后作为下次判断的参考依据保存在信息库里。例如,图 10所示的是东山台于2020年10月11日12:30—14:30期间观测到的地震电信号特征图,异常信号出现后,笔者于2020年10月12日在图中用红色字体给出了基本判断:28天后,当在11月8日左右台湾宜兰至花莲区域将会发生M5.0左右地震的标注,以等待后续验证。结果在2020年10月17日在台湾台南市发生M4.9;2020年11月6日在台湾台东县海域发生了M4.9地震。地震发生后,再将地震三要素(时间、地点、震级)发生的具体信息标注在图上,这样就形成了地震电信号特征信息库里的一个元素,以方便未来比对信息之用。
图 10 2020年10月11日东山台地震电信号特征图Fig.10 Characsteristic of seismic electric signals observed at Dongshan Station on October 11, 2020
利用表1给出的地震电信号的异常信息,是否可以做出有效的短临地震预测是检验闽台地区基于地震地电场方法短临预测地震的最佳方式之一。我们认为这种方法在本文所监测的区域里是可行的。例如2019年3月14日笔者与一位专家在微信里交流时曾经给出了台湾地区地震短临预测意见,该意见认为:“关于台湾近期震情,我觉得发生6.0级左右的地震才对,4.6级地震太小了。下图是东山岛的观测,请参考。”结果在2019年4月18日台湾花莲县海域发生了6.7级地震,具体如图 11所示。在该意见中所指的东山岛观测的图即图3。
图 11 2019年3月14日微信预测意见及实际发生地震图Fig.11 WeChat prediction comments on March 14, 2019 and corresponding actual earthquake occurrence
在2020年6月3日笔者在一个三人地震群里给出了预测意见:“震情交流:未来三周内在台湾花莲至宜兰及周边地区发生M5.6左右地震的可能性很大,临震时间很可能是6月6日至6月16日期间。5月25日的4.8级地震难以解释异常信号。”实际上,2020年6月14日在台湾宜兰海域发生5.5级地震,如图 12所示。
图 12 2020年6月3日微信群预测意见及实际发生地震图Fig.12 WeChat prediction comments on June 3, 2020 and corresponding actual earthquake occurrence
在分析表1中的一些地震电信号特征后笔者认为很快在台湾地区将要发生一次中强地震,2021年12月27日笔者在一个三人地震群里给出了预测意见:“应在两周内台湾会发生6.3级左右的地震”。实际上,在预测意见给出的6天后,即2022年1月3日,台湾花莲县海域发生了6.4级地震,具体参见图 13。
图 13 2021年12月27日微信群预测意见及实际发生地震图Fig.13 WeChat prediction comments on December 27, 2021 and corresponding actual earthquake occurrence
由上可以看到,随着科技发展,手机微信通讯和交流平台的广泛应用,使得资讯通讯的时间都可以在手机上标示并截留,为我们进行地震短临预测意见的交流和保留证据提供了很好的平台。这能够证明地震预测结果除了填写地震预测卡并向当地地震局或中国地震台网中心邮寄外,还多了一个为我们对地震短临预测实践探索活动留取证据的方式。上述实践中地震短临预测都比较成功,说明在福建沿海进行地电场观测并用于闽台地区地震短临预测可以做出有效的预测。
另外,在表1给出的震例中,其中最大的地震是发生在2024年4月3日台湾花莲县海域的M7.3地震,该地震给台湾地区造成了很大损失,同时也引起了中国大陆东南沿海地区如厦门、福州、泉州、宁德等地的强烈震感。上海、浙江、江苏、江西等地有感。很可惜的是由于经费短缺,我们在东山台、漳州台、平潭台建立的地电场观测运行后续难以为继,只好停测。所以这次大地震只在泉州台检测到了地震电信号,因为信息太少,难以对其三要素进行准确预测,也错失了这次短临预测地震的实践机会。相信未来本区地电场观测台站数增多后基于本文方法的地震短临预测会取得更好结果。
由本文可以看到,福建沿海地电场观测对闽台地区地震的发生是可以在短临预测阶段观测到地震电信号的,在我们临时建立的东山台、漳洲台、平潭台观测期间结合泉州台地电场观测,得到了一些地震特别是台湾地区强地震前的异常电信号,几个较好的震前短临地震预测意见很好地支持了所得到的地震电信号的可信度。这几个震前短临地震预测意见也再次证明了地震地电场指纹信息法的有效性。在上述3个临时性的地电场观测台运行停止后只由一个泉州地电场台站的观测资料分析结果进行较准确的判断2024年4月3日台湾7.3级地震三要素还是较为困难的,由此可见如果能够在福建沿海地区多建几个观测质量高的地电场观测台站则对闽台地区地震震情(尤其是短临阶段的震情)研判可以起到非常重要的作用。
另外,在我们临时性的地电场观测期间,通过东山台和平潭台使用微型地电场仪的小极距地电场观测可以看到,小极距地电场观测系统对观测场地要求比ZD9大地电场仪要低许多,电极布设方式操作简单易行,因此建立观测站点的成本也非常的低;在地电场的阵列观测或点阵观测中,我们在前期的工作中得到了体现,本文工作再次证明了小极距地电场点阵观测的快速便捷之处。
致谢:福建省地震局林树研究员、福建省地震局监测预报处、平潭地震台诸同志在本项工作中给予了大力支持,在此深表感谢!