0 引言

地震预测预报是当今世界的科学难题，近年来国内外学者利用电磁学方法研究地震的孕育与发生过程，记录到震前大量电磁异常，尤其是甚低频（VLF）到超低频（ELF）电磁异常信号，这些信号出现在震前几天或几小时，提供了明显的短临地震信息。因此，在大地震短期预报中，利用地震的电磁现象进行研究逐渐被认同，电磁方法成为一种有前途的地震研究方法（高曙德等，2010，2013；张建国等，2011，2013；孙维怀等，2012）。

“十五”期间，中国地震局在全国安装12套极低频电磁观测仪器。红山基准地震台（以下简称红山台）于2007年5月安装架设ELF电磁观测仪器，2009年8月按照统一规范进行连续观测，仪器总体运行良好。本文通过对该仪器观测资料的分析处理，探讨河北及国内大地震前后观测的电磁场功率谱密度及视电阻率中短期异常和变化规律。

1 地质构造背景

红山基准地震台（以下简称红山台）为河北省地磁基准台，地处隆尧县城西北9 km处，海拔37 m。红山台所在地区构造上属于邢台地堑区，东临沧县隆起，北接冀中坳陷，西为太行山隆起，南为隆尧断裂。邢台地堑总体构造以NNE向为主，内部发育一系列NNE—NE向断裂带，局部发育NW—NWW向断裂。受板块构造影响，华北地区的主压应力方向以NEE—SWW为主导，构造应力张量结构以走滑性质为主，兼有一定的正断性质。红山台周围地质构造复杂，基层断裂纵横交错。特别是东南部束鹿断陷为新生代断陷盆地，断陷内新构造运动发育强烈，第三系以来的沉积层厚达700—2 000 m。台站基岩为震旦纪石英砂岩，部分基岩出露，广布第四纪冲洪积相堆积（徐继山等，2011）。红山台地质构造见图1。

图1 红山台地质构造 Fig.1 Geological structure of Hongshan Seismic Station

2 数据观测及处理

2.1 数据观测

红山台使用俄罗斯ACF-4M型ELF电磁观测仪进行电磁场观测。电极采用正NS—EW向布设，观测4个分量的电磁场信号：2个分量水平电场（EW向Ex和NS向Ey）、2个分量水平磁场（NS向的Hx和EW向的Hy），目前主要观测天然电磁场，观测频率范围为0.1—800 Hz，电场传感器采用Pb-PbCl2不极化电极，与磁场传感器（感应线圈）共同埋入地下，连接线采用80 cm深度埋地敷设（范晔等，2013）。台站观测仪器布设见图2。

图2 台站观测仪器布设 Fig.2 Layout of observation instruments at the station

红山台每日分2个时段进行电磁观测，记为01、02，每个时段记录3个观测频段，记为 D1、D2、D3，其中 D3为高频段（1 — 800 Hz），D2为中频段（1 —400 Hz），D1为低频段（0.1 — 40 Hz），各频段采样率分别为 3 200 Hz、1 600 Hz、160 Hz。各频段每日观测时间统计见表 1。

4)鼓励游客在淡季进行错峰旅游体验，利用AR的虚实结合技术，弥补景点的季节性景色缺憾，对于低预算游客来说，既可以获得性价比较高的旅游资源，也可以欣赏到当季的“美景”。

表1 各频段每日观测时间 Table 1 Daily observation time list of each frequency bands

观测顺序 观测频段 观测时刻 观测时长/s 观测顺序 观测频段 观测时刻 观测时长/s 1 D3 00:01:00 1 200 4 D3 05:35:00 1 200 2 D2 00:25:00 3 600 5 D2 06:00:00 3 600 3 D1 01:30:00 14 400 6 D1 07:05:00 14 400

2.2 视电阻率、功率谱计算

ELF观测仪具有计算机网络设置功能， 使用专用软件设置仪器的各项参数、自动记录开机时间和记录长度。记录数据传输到计算机后，使用软件进行处理，对时间序列进行FFT（快速傅式变换）变换，实现时间域到频率域的转换，得到电场、磁场的自功率谱S(E)、S(H)及视电阻率ρxy和ρyx和阻抗相位Z（陈小斌等，2006；范晔等，2013）。其中Z和ρ由下式计算得到。

红山台观测到的天然电磁场的电场和磁场分量的自功率谱曲线变化趋势整体具有较好的一致性，见图3。与磁场数据相比，电场数据受干扰影响普遍较大，主要是因为，电场信号受台站所在区域地质构造、观测场地环境、人为干扰等因素影响较大，具有较强的局部效应，而磁场信号则具有更好的信噪比和区域可对比性（范晔等，2013）。由图3中1 Hz、39 Hz、178 Hz的磁场NS和EW道自功率谱曲线可观察到明显的背景变化规律。红山台地处北半球中纬度地区，随着春夏秋冬的季节变化，磁场曲线具有夏季升高、冬季回落的规律，且以1年为周期，高频178 Hz的磁场强度明显弱于低频1 Hz、39 Hz的磁场强度，因此磁场强度随着频率增加缓慢减弱。电场数据同样呈现夏强冬弱的变化规律，但变化规律没有磁场数据显著。

Sphinx同样使用的CMU语言模型工具。其中ARPA文件也使用mkbingram转换为二进制格式，以提高解码的效率[15]。

其中，RE(t)为时间域信号。

那年部队来征兵，招财报名，所幸的是老支书庆高伯卸职去公社林场了，不然的话会是第一个就被刷下来。他在部队几年，是站岗、放哨？还是练武、杀鬼子？谁也不知道。他娘桃娭毑也不知道。他虽然小时候读过三个一年级，两个二年级，后来还混到四年级，可他不会写信，就是写了信寄回来，桃娭毑也不会看，她扁担倒地，一字不识。

频率域功率谱SE(f)可由下式计算得到。

其中，Ex为电场水平强度，Hy为垂直于Ex的磁场强度，T为周期。

3 观测记录变化规律

天然电磁场信号微弱，易受到各种干扰的影响。检查红山台2014年3月7日—2015年12月31日观测日志，查看每日观测和数据缺失记录，去除人为干扰、雷雨等自然环境干扰因素的影响，确保电磁数据在仪器正常运行状态下进行记录和处理。

使用Matlab编写软件，绘制每日各频率段电磁场分量的自功率谱叠加曲线，从中挑选舒曼谐振背景变化规律明显且质量较好的数据，同时在D1—D3频段中分别筛选代表频率 1 Hz、8 Hz、39 Hz、82 Hz、178 Hz、580 Hz；编写频率提取软件，抽取选定代表频段数据，绘制随时间变化的电磁场分量自功率谱和视电阻率观测曲线。

多层次分重点培训,主要内容:医院感染诊断标准及报告处置流程、手卫生、消毒隔离、标准预防等知识。通过培训,增强医务人员的感控意识。

其中，在企业内部精细化管理工作当中需要设置专门的审计人员，对企业相关的财务业务进行审计。但是大部分企业都没有设置专门的审计部门，财务审计工作都是由企业原有的财务部门工作人员来完成的，但是企业财务管理人员并不是专业的审计人员，审计工作实施的专业性、规范性较差，而且财务人员没有掌握高效的审计方法和理念，审计工作完成的效率也受到了影响，企业内部审计工作的实际效果没有真正的发挥出来。

表2 红山台代表频率电磁场强度分布范围 Table 2 Distribution range of electromagnetic field instensity of the selective frequency bands at Hongshan Seismic Station

时间 测项舒曼谐振频率1 Hz 8 Hz 39 Hz 82 Hz 178 Hz 580 Hz 2014 年 H/((A/m)2/Hz)10-13—10-11 10-14—10-11 10-14—10-10 10-14—10-11 10-15—10-12 10-16—10-13 E/((V/m)2/Hz)10-13—10-11 10-14—10-11 10-13—10-10 10-14—10-11 10-15—10-12 10-15—10-13 2015 年 H/((A/m)2/Hz)10-13—10-10 10-14—10-11 10-14—10-10 10-14—10-11 10-15—10-12 10-16—10-13 E/((V/m)2/Hz)10-13—10-11 10-14—10-11 10-13—10-10 10-14—10-11 10-15—10-11 10-16—10-12

近年来，随着城市化进程的不断加快，政府储备土地普遍存在管理难、利用难等老大难问题，丽水市委、市政府高度重视，审时度势，富有创意地提出并成立了全省首家丽水市储备土地经营管理有限公司。该公司于2017年3月注册成立，7月份完成人员招聘、队伍建设等工作，通过建队伍、搞调研、定制度、抓管理、拓业务，一年来，实现了良好起步，特别是在推进丽水市国家文明城市创建、闲置储备资产管理等方面取得了阶段性成效。

由电磁场功率谱曲线得到红山台选定频段电磁场强度日常分布范围，见表2，去除电磁场强度日常分布范围外的明显异常值，尽量保证观测数据的可靠性。

图 3 1 Hz、39 Hz、178 Hz 电磁场自功率谱 Fig.3 1 Hz、39 Hz、178 Hz electromagnetic self-power spectrum

4 异常分析

极低频电磁异常信息一般出现在地震发生前数日至数十日，在震中距几百千米范围内均可记录到异常信息（郝建国等，1995）。河北及邻区小震频发，大震缺少，选取2014—2015年河北及邻区3级以上且红山台与震中距离200 km左右的4次地震，以及距震中均1 000 km以上的云南鲁甸6.5级地震、尼泊尔8.1级地震，对记录6次地震的ELF资料，计算自功率谱，分析视电阻率变化。参与分析的地震参数见表3。

表3 2014年3月7日至2015年12月31日红山台ELF资料计算功率谱与视电阻率变化与对应地震统计 Table 3 ELF data from March 7,2014 to December 31,2015 and corresponding eathquakes

发震时间 φN/(°) λE/(°)震中位置 M 震源深度/km 台站与震中距离/km 2014-08-03T16:30 27.1 103.3 云南鲁甸 6.5 12 ＞1 000 2015-01-18T11:01 35.64 115.38 河南范县 3.7 7 216 2015-04-19T18:21 38.9 116.3 河北文安 3 29 210 2015-04-25T14:11 28.2 84.7 尼泊尔 8.1 20 ＞1 000 2015-06-02T18:51 37.80 112.54 山西太原 3.1 21 192 2015-12-10T13:27 38.72 113.01 山西原平 3.1 9 193

4.1 自功率谱变化

绘制红山台2014—2015年电磁场观测数据1 Hz、39 Hz、178 Hz的自功率谱变化曲线，见图4（图中蓝色曲线表示H1、E1、红色曲线表示H2、E2）。从图4可见：①在2014年8月3日鲁甸6.5级地震发生前10天内，电场和磁场功率谱在不同频率段均有显著变化，变化幅度分别高于正常值约1个数量级，电场功率谱在各频率段均存在1个显著的台阶变化，地震发生后，数据恢复正常，台阶消失；②2015年4月19日文安3.0级地震、4月25日尼泊尔8.1级地震发生前十几天内，各频率段电场和磁场功率谱变化显著，变化幅度分别高于正常值1—6个数量级不等，尤其电场功率谱变化显著，明显的向上突跳和阶跃，为脉冲集丛现象，而磁场功率谱变化为向下的阶跃和突跳，在地震发生后，数据变化恢复正常，查看国家地磁台网中心磁暴记录，以上地震发生时段无太阳活动和磁暴；③表3中给出的河南范县、山西太原及原平3次地震在1 Hz、39 Hz、178 Hz电场、磁场自功率谱曲线中并未有明显震前异常，可能是因为震级较小，能量衰减较快，在受干扰较大的自功率谱曲线中无异常出现。对于鲁甸6.5级、尼泊尔8.1级强震，ELF观测资料的震前电磁效应明显。

图4 2014年3月至2015年12月红山台电磁观测数据不同频率自功率谱曲线 (a)1 Hz；(b)39 Hz；(c)178 Hz Fig.4 Different self-power spectrum densities of electromagnetic observations data at Hongshan Seismic Station from March 2014 to December 2015

4.2 视电阻率变化

本文利用红山台2014年3月7日—2015年12月31日ELF观测的天然电磁场数据，采用电磁场功率谱密度对背景变化规律及观测资料异常特征进行分析，并通过视电阻率变化辅佐分析，得到以下结论。

图 5 1 Hz、39 Hz、178 Hz 频率 Ryx方向的视电阻率映震曲线 Fig.5 The relationship between the Ryx apparent risistivities of 1 Hz、39 Hz、178 Hz and earthquake

5 结论

红山台为电磁场定点观测台站，观测点附近各地层电阻率基本稳定，计算Rxy、Ryx方向的视电阻率，并绘制观测曲线（图略），可知Ryx方向曲线较为平滑，Rxy方向显示突跳较多。因此，分别绘制1 Hz、39 Hz、178 Hz 视电阻率Ryx方向曲线，分析视电阻率变化特征。从趋肤公式可知，1 Hz、39 Hz、178 Hz分别对应红山台地下由深部至浅部的视电阻率，即由高阻至低阻，越向上至地表，松散第四纪沉积物电阻率越低。绘制3个频率段Ryx方向的ELF视电阻率变化曲线，见图5，可以看出，鲁甸6.5级地震、尼泊尔8.1级地震和文安3.0级地震与不同频率段视电阻率变化均有较好的对应关系，视电阻率向上或向下明显跳动，变化的异常幅度达1—2个数量级，应与震级大小、震中距、红山地质构造环境等相关；39 Hz、178 Hz（中高频段）的视电阻率变化较1 Hz（低频段）更显著，充分说明视电阻率具有方向性和频率所选性。

（1）电磁场强度呈现夏季强冬季弱的变化趋势，电场信号受干扰影响较大，磁场信号则具有更好的信噪比和区域可对比性。

（2）红山台各地层电阻率基本稳定，1 Hz视电阻率高于39 Hz，39 Hz视电阻率高于178 Hz，从趋肤公式可知，1 Hz、39 Hz、178 Hz分别对应红山台地下由深部至浅部视电阻率，即由高阻至低阻，越向上至地表，松散第四纪沉积物电阻率较低。

（3）电磁场自功率谱及视电阻率变化曲线中出现的一些突跳并未伴随地震发生，有些地震未发现明显异常变化，可能原因为地震震级较小，电磁信号较弱，有效信号淹没在环境噪声中。台站ELF记录天然电磁场信号，场源主要来自雷电产生的电磁场和高空电离层扰动产生的电磁场以及各类人为干扰、环境干扰，也有来自地震前的电磁辐射，但单个地震台站无法确定观测的异常是否来自人为干扰、雷电信号或震前异常，因此，为了从复杂的电磁环境中提取相对较微弱的地震电磁信号，需要加强对信号提取方法的研究，需要采用同台站多手段或其他台站相同观测手段进行对比分析及异常跟踪，并考虑传统的数学分析方法与一些信号处理分析方法的结合使用。

（4）选取的1—178 Hz电磁场自功率谱在2014年8月3日鲁甸6.5级地震，2015年4月19日文安3.0级地震和4月25日尼泊尔8.1级地震发生前十几天内均发生显著变化，包括台阶、向下的阶跃变化以及脉冲集丛现象，异常幅度为正常月份变化的1—6个数量级。且地震发生后数据恢复正常变化。查看国家地磁台网中心磁暴记录，以上地震发生时段均无太阳活动和磁暴。

（5）鲁甸6.5级地震、文安3.0级地震和尼泊尔8.1级地震前，视电阻率均有显著变化，视电阻率向上或向下明显跳动，变化的异常幅度达到1—2个数量级，与地震震级大小、震中距、红山地质构造环境等相关。中高频段变化显著，充分说明视电阻率具有方向性和频率所选性，与电磁场自功率谱变化对应性较好。

在开展电磁研究时，中国地震局地质研究所赵国泽研究员、汤吉研究员、王立凤老师给予指导和帮助，红山基准地震台提供ELF数据，河北省局测震台网提供地震目录，在此一并表示感谢。

参考文献

陈小斌，赵国泽，马霄.大地电磁三维模型的一维二维反演近似问题研究[J].工程地球物理学报，2006，3（1）：9-15.

范晔，汤吉，赵国泽，等.电磁监测台站观测的舒曼谐振背景变化[J].地球物理学报，2013，56（7）：2 369-2 377.

高曙德，汤吉，杜学彬，等.汶川8.0 级地震前后电磁场的变化特征[J].地球物理学报，2010，53（3）：512- 525.

高曙德，汤吉，孙维怀.盈江5.8级和缅甸7.2级地震前电磁异常[J].地球物理学报，2013，56（5）：1 538-1 548.

郝建国，张云福，潘怀文，等.震前极低频电磁异常及其频谱特征[J].地震学报，1995，17（1）：81-88.

孙维怀，汤吉，张平，杨学慧.云南ELF电磁观测资料的异常特征分析[J].地震地质，2012，34（3）：488-499.

许继山，彭建兵，马学军，等.邢台市隆尧地裂缝发育特征及成因分析[J].工程地质学报，2011，20（2）：160-169.

张建国，焦立果，刘晓灿，马新欣.汶川MS 8.0级地震前后ULF电磁辐射频谱特征研究[J].地球物理学报，2013，56（4）：1 253-1 261.

张建国，刘晓灿，马新欣，焦立果.玉树7.1级地震前后电磁辐射小波变换分析[J].大地测量与地球动力学，2011，31（6）：32-38.