自然电场的区域性变化特征

谭大诚，赵家骝，刘小凤，范莹莹，刘君，陈军营

1中国地震局兰州地震研究所，兰州 730000

2兰州地球物理国家野外科学观测研究站，兰州 730000

3中国地震局地震预测研究所，北京 100036

1 引言

人类在19世纪初已经观测到地球表面存在电流，1940年S.Chapman和J.Bartels提出了大地电场、自然电场的概念（孙正江和王华俊，1984）；1984年，希腊学者提出从地电场观测资料中提取SES（Seismic Electric Signals）信号预报地震的VAN法，其原理主要应用了自然电场数据（马钦忠等，2004；马钦忠，2008），但迄今也存在争议（Geller，1996；黄清华，2005；黄清华和林玉峰，2010）.赵国泽等（2009）、汤吉等（2010）开展的 MT测深反演主要应用大地电场数据，而频谱分析（黄清华和刘涛，2006；范莹莹等，2010）、常规波形分析（钱复业和赵玉林，2005；仇勇海等，2008）、极化方位计算（毛桐恩等，1999）、长短极距比值计算（田山等，2009）、数据处理（安张辉等，2010）等应用的地电场数据归类相对模糊.

近年，黄清华和林玉峰（2010），Huang and Lin（2010），Huang等（2011），Ren等（2012）对SES信号数值模拟与分析以及对DC-ULF电磁信号微观机理与传输通道的研究、马钦忠和钱家栋（2003）提出的非均匀电性结构导致地表电场复杂、其他学者开展的数据分析（钱复业和赵玉林，2005；仇勇海等，2008）与观测方法（席继楼等，2013）等研究，使得自然电场局部场地的机理、特征与分析更深入，但应用自然电场动电效应、接触-扩散和氧化还原等微观机理（孙正江和王华俊，1984；钱复业和赵玉林，2005；Huang et al.，2011），研究自然电场的区域变化特征、探寻微观机理与宏观现象的关联等方面，目前成果依然很少.

赵国泽和陆建勋（2003）指出视电阻率日变化与固体潮存在关联；黄清华（2005），Huang等（2011）认为从复杂电磁环境中提取相对较弱的地震电磁信号，有必要探寻数学、信号处理与地震电磁物理过程结合的物理解析方法；2006—2007年，黄清华和刘涛（2006），张学民等（2007）分析了地电场固体潮汐现象的频谱特征及震前异常.基于徐文耀（1992），赵旭东等（2008）对空间Sq电流、黄清华等对地电场固体潮汐波的研究成果，谭大诚等（2010）将地电场潮汐波分类为TGF-A、TGF-B两类，探讨了自然电场、大地电场物理解析初步原理（谭大诚等，2012）.

本文依据自然电场、大地电场起源分类原则（孙正江和王华俊，1984），充实了地电场构成解析原理和方法（谭大诚等，2012），通过对南北地震带附近有代表性的13个台站、66个月的自然电场分析，总结了该地区自然电场的场地、区域中长期变化特征，解析了其中2个场地短期裂隙水主体渗流方位逐日变化（谭大诚等，2011，2013）与自然电场变化的关联，探讨了松山场地2年中自然电场变化与裂隙水渗流、主压应力、裂隙结构变化的微观机理过程.

2 地电场构成解析基础与方法

自然电场源于固体地球内部的物理、化学作用，场源主要包括地下矿体、岩体及界面、应力应变、地下流体等，这导致其具有局部性特征.依据自然电场形成的不同微观机理与过程，可分类为过滤电场、接触-扩散电场和氧化还原电场等（孙正江和王华俊，1984）.通常，认为自然电场相对稳定（孙正江和王华俊，1984），地震电磁前兆可能更易表现在自然电场变化中（马钦忠等，2004；马钦忠，2008）.

大地电场源于固体地球外部的空间电流系和潮汐力（孙正江和王华俊，1984；黄清华和刘涛，2006；谭大诚等，2010，2011，2012，2013），场源主要包括电离层、月球和太阳的潮汐力、太阳风等，这使其具有全球性或区域性特征.依据大地电场形成的不同微观机理与过程，可分类为TGF-A、TGF-B、地电暴（孙正江和王华俊，1984；谭大诚等，2010）等.依据谭大诚等（2010，2011，2012，2013）近年的研究，TGFA、TGF-B潮汐地电场和地电暴的强度受场地岩性结构、裂隙和裂隙水、应力应变等因素的明显影响，因此地震电磁前兆也能够表现在场地的大地电场变化中（黄清华和刘涛，2006；张学民等，2007；谭大诚等，2010，2011，2012，2013）.

按照地电场起源分类原则，自然电场、大地电场可部分具有同一种微观机理过程，例如TGF-A、TGF-B潮汐地电场形成过程中，裂隙水或裂隙水中电荷发生了周期性渗流或移动（谭大诚等，2010，2011）；在构造活动中，应力应变导致含水岩体裂隙扩展、破碎过程也存在裂隙水和裂隙水电荷的渗流或移动，其形成的地电场属自然电场.因此，应用大地电场观测数据，可能探讨场地自然电场宏观变化的微观机理.

在地电场观测数据中，以ESP代表自然电场（spontaneous field），ET代 表 大 地 电 场 （telluric field），干扰成分用Er表示，则地电场E的构成可写成：

2008—2012年，中国大陆100余个台站的地电场日变形态基本有TGF-A、TGF-B和无日变波形三种，前两种波形约占总数的三分之二（谭大诚等，2010，2011，2012），无日变波形地电场在岩性和电性结构复杂、构造活动剧烈的南北地震带附近较集中（谭大诚等，2012）.由此，可推论ET稳定的主要成分是潮汐地电场，TGF-A、TGF-B日变波形分别与固体潮、Sq电流的波形对应，其前10阶潮汐谐波的周期是23～24h、12h、7.9h、6h、4.8h、4h、3.4h、3h、2.7h、2.4h（黄清华和刘涛，2006；赵旭东等，2008；谭大诚等，2010，2011，2012，2013）.

图1 基于潮汐地电场计算裂隙水主体渗流α角示意图Fig.1 Sketch of calculatingαof the preferred orientation of crack water seepage based on tidal geoelectric field

在电磁环境和观测系统较理想时，Er在较小范围内随机波动.设地电场观测给出的分钟值为Ei，应用一天观测数据计算ESP的表达式可简化为（谭大诚等，2012）：

应用（2）式计算ESP时，因潮汐谐波的周期特点，潮汐地电场的主要影响基本消除（谭大诚，2012）.对无日变波形的地电场，FFT频谱分析表明潮汐谐波存在，多数情况下还是主要谐波成分（谭大诚等，2012），（2）式仍成立.

设第j天的自然电场为ESP（j），则ESP逐日跳变可用（3）式描述（谭大诚等，2012）.

场地岩体裂隙主要在构造作用下形成，小尺度岩石块体内的微裂隙分布、裂隙水渗流会十分复杂，但宏观大尺度岩体内的地下水渗流主要沿裂隙进行（周志芳等，2007；陈颙等，2009；谭大诚等，2011，2013）.数百米量级的地电场观测场地属宏观大尺度，观测到的裂隙水渗流电场方向则接近裂隙水主体渗流方位（谭大诚等，2011，2013）.

对图1a装置系统，在场地裂隙水周期性主体渗流方位如图1b所示时，则NS、NW方向地电场潮汐波相关性更高（谭大诚等，2011）.设潮汐地电场（大地电场主要成分）为ET，裂隙水主体渗流方位以北偏东角度α代表，应用日变波形前10阶潮汐谐波振幅和，逐日计算场地岩体α角的方法推导如下.

中国地电场台站的NS、EW方向的极距通常相等，由图1c得到：

故

潮汐地电场ET是周期性变化场，原理上一天中任意时刻的ET（NW）、ET（NS）均可以用于（6）式计算.在潮汐波峰、谷前后，潮汐电场强度最大、稳定性相对好，故（6）式宜取峰值或谷值计算.

在图1d中，选取潮汐波零值为X轴，设ET（NS）峰、谷值分别为 A1、A2，ET（NW）峰、谷值分别为 B1、B2，通常A1≠A2，B1≠B2，于是：

则

将ET（NS）、ET（NW）峰谷值分别以ANS、ANW表 示，由（7）、（8）式可得：

对日变波形清晰的地电场，（9）式已可计算.在日变波形畸变或难于识辨时，ANS、ANW取值存在较大误差，而前10阶潮汐谐波周期稳定，日变波形幅度主要取决于潮汐谐波的前5阶幅度（谭大诚等，2011，2012，2013）.

图2（a，e）是泸沽湖台TGF-A波形、古丰台TGF-B波形的分钟值曲线.图2（b，f）是其日变峰谷值A、前10阶潮汐谐波振幅和AT的逐日取值曲线，在 NS、NE（NW）两个方向，ANS、ANE（NW）变化趋势一致，AT（NS）、AT（NE／NW）变化趋势一致，且 A、AT变化趋势相近，这说明图1模型与原理基本符合这两个场地的实际情况.

图2（c，g）是场地 ANE（NW）与 ANS、AT（NE／NW）与AT（NS）比值的逐日变化曲线，在这两个场地存在：

其中Ai是第i阶潮汐谐波振幅，由（9）、（10）式得到：

在观测数据FFT结果中，稳定的直流成分已排除（谭大诚等，2011），因此依据潮汐地电场频谱特征按照（11）式计算α时，ESP的影响基本消除.应用该式计算α的误差，理论上主要源于图1b模型与实际岩体的吻合度，这取决于场地岩体裂隙的发展阶段（周志芳等，2007；陈颙等，2009；谭大诚等，2011，2013）.一般来说，岩体裂隙的优势方位愈明显，计算α的误差应愈小、稳定性愈好.

在2008年汶川MS8.0地震前，泸沽湖、盐源两场地无明显干扰、观测系统正常.依据地电场前兆信息可信度的统计分析RTL（Region-Time-Length）方法（Huang，2006），应用（11）式计算的两场地α的变异相互对应，其变异在时间、空间上与此次地震相关，且异常显著（谭大诚等，2013）.

应用（9）、（11）式于泸沽湖、古丰场地，图2（d，h）中α值逐日差异均值分别为1.86°和1.38°，可见两种计算方法的结果通常较稳定且相当接近，在多个场地应用这两式计算对比也没有明显差异（谭大诚等，2013）.依据场地的主压应力σ1方位理论推算，在剪裂时这两式的计算结果通常会偏小5°左右（周志芳等，2007；陈颙等，2009；谭大诚等，2011，2013）.

应当指出，在装置系统变更或裂隙水主体渗流方位更接近EW 向时，（9）、（11）表达式需相应调整；无潮汐日变的地电场，潮汐谐波通常也是其基本谐波（谭大诚等，2012），这时（9）式中峰谷值已无法客观提取，但（11）式依然可逐日计算.

3 自然电场区域变化特征

2008—2013年，大震、强震频繁发生在南北地震带附近，区域构造活动剧烈.因此，本节在总结青藏高原东缘海原、渭河、龙门山和盐源断裂带附近ESP区域特征中，采用了对多个场地数据长期变化的“趋势性改变”进行对比，并将各场地1、2天左右或短时不确定的大幅突跳数据删除，以显示出ESP变化趋势.

3.1 主断裂带附近自然电场区域变化特征

图3a是上述断裂带附近有代表性的13个台站分布图；图3b绘制了海原断裂带东西两端的景泰和固原、渭河断裂带附近西安和宝鸡的自然电场变化曲线.可见：①同一场地ESP、ΔESP变化具有方向性和局部场地特征.②在同一断裂带附近的不同场地，ESP发生大幅跃变、突跳或渐变的趋势性转折点多具有时间同步或准同步性，同时段数据稳定性大致类似，但形态变化不具相关性.此外，ΔESP跳变曲线反映了ESP突跳变化，对其渐变过程反映弱.

图2 泸沽湖与古丰台地电场日变波形峰谷值与前10阶潮汐谐波振幅和（2009-03-15—04-15）Fig.2 Peak-to-trough values and amplitude sum of the first 10order tidal harmonic waves of diurnal variations of the geoelectric field at Lugu Lake and Gufeng stations（2009-03-15—04-15）

图3c中，2012年6月初，成都和汉王ESP、ΔESP大幅跃变、突跳等趋势性转折点时间对应较好；10月后，甘孜ESP（NS）出现与这两台时间准同步的大幅突跳、阶跃.图3d中，西昌台阵泸沽湖、盐源和小庙台在2008—2009、2012—2013年间的数据变动大.总体上，将龙门山、盐源断裂带附近各场地的ESP、ΔESP数据曲线分别对比，其特征与图3b中的①、②结论基本相同.

应指出，2008—2013年中期，上述台站观测系统基本没有出现长时间故障，成都、汉王、宝鸡、小庙等场地存在相对复杂的电磁环境，图3（b—d）表明其ESP、ΔESP趋势性转变仍然与相邻场地具有较好的时间同步性；图3（b，d）的小黑色箭头表明电极更换可能导致ESP出现阶跃；2010年5月初前后，小庙台布极区内埋设金属管道的影响在4个月内较明显.

3.2 局部区域内自然电场变化特征

青藏高原东缘2个地电场台阵建于2007年，西昌台阵ESP变化特点在3.1小节已分析，天祝台阵由古丰、景泰、松山、红沙湾和黄羊川5个台组成，相邻间距不超过50km.

图4a是天祝台阵地电场台、邻区断层分布图.图3a表明景泰、古丰分别接近海原和祁连山主断裂带，图4a表明其他3个台附近还有较小的断裂带.

图4b绘出了古丰、松山和红沙湾场地ESP和ΔESP曲线，景泰场地曲线见图3b.2008—2012年，这四个场地ESP显示出年变形态，但不同方向、不同时段，ESP年变波形及稳定性不同.2009年2月、2010年10月、2012年2月，古丰ESP、ΔESP趋势转折前后，松山、红沙湾场地数据有相应剧变对应.可见，在这一局部区域，各场地ESP、ΔESP变化特征与图3b中的①、②结论同样基本相同.

图3 主断层附近典型台站ESP、ΔESP趋势转折（2008-01-01—2013-06-30）Fig.3 Trend turnings of ESPandΔESPat representative stations adjacent to major faults（2008-01-01—2013-06-30）

由图3、图4可知，南北地震带附近多数场地ESP不具有清晰年变波形，天祝台阵内多数场地ESP年变波形不稳定；在ESP的变化形态中，持续数月较平稳的升降、数日内的快速跃变、阶跃、突跳等现象常见；通常，场地ESP、ΔESP变化表现出方向性和局部场地特征.同一主断裂带附近、局部区域内，不同场地ESP大幅跃变、突跳或渐变等趋势转折点多具有时间同步或准同步性，同一时段数据稳定性可能类似，但形态变化存在场地因素而不具相关性.

4 自然电场变化机理初步解析

4.1 自然电场跃变、突跳的可能性

场地ESP不具有明显的短时周期性，震情复杂区域其跃变、突跳等变化的可信度较难判断.在图3、图4中，同一主断裂带附近或局部区域内，多台ESP跃变、突跳等趋势性转折点时间同步或准同步性表明了其现象具有相当的可靠性.对前几年的震例分析中，这种大幅阶跃、突跳也常见（谭大诚等，2012）；图5（a，b）分别是新疆温泉、四川甘孜台站的体应变与地电场观测曲线，在数百米范围内，场地应变变异或部分波形转折处，地电场出现了明显的跃变、突跳.20世纪80年代后，郭自强等（1999）、郝锦绮等（2003），钱书清等（1996）在岩石破裂实验中也检测到电磁信号跃变、突跳现象，只是这些实验检测的电磁信号频率更高.

图5中的地电场E曲线以分钟观测值绘制，在其发生跃变、突跳时，周围其他台地电场没有类似变异现象，这期间温泉NS、甘孜EW方向的地电场日变波形基本清楚.可见，这两台E曲线的跃变、突跳是ESP变化导致.

应指出，ESP源于地下介质的物理、化学过程，其机理不仅与应变关联，也与流体渗流、流体与介质的化学作用等相关，因此场地ESP的跃变、突跳并不能始终与应变变异对应，但图5所示情况表明场地ESP在构造活动中能够发生跃变、突跳等现象.

4.2 自然电场变化机理初步解析

图4 天祝台阵台站分布及ESP、ΔESP趋势转折（2008-01-01—2013-06-30）Fig.4 Distribution of stations and the trend turnings of ESPandΔESPof the Tianzhu array（2008-01-01—2013-06-30）

图5 地电场（上）与体应变（下）对应变化Fig.5 Corresponding variations of the spontaneous field（upper）and the volume strain（lower）

应用主断裂带附近或局部区域内多台ESP、ΔESP变化的对应关系，基本可识别其可靠性，但不能认识其变化的微观机理.地电场是自然电场ESP、大地电场ET的叠加.ET虽起源于空间电流、潮汐力，但在具体观测场地，各方向ET的日变波形特征与ESP都受场地水文地质、应力应变、流体渗流等因素的影响（陈颙等，2009；谭大诚等，2011，2012，2013）.因此，同场地的ESP、ET的变化应存在关联性.

以南北地震带南段的盐源、北段的松山台站为例，分别选取这2个台站ESP变化较剧烈的一段时间，依照（11）式解析其场地裂隙水主体渗流α角，结果如图6所示.

2012年9月7日云南彝良发生MS5.6地震（104.03°E，27.33°N），图6a是震前30余天的盐源场地ESP（EW）、α变化曲线，上、下两曲线变化形态具有时间对应性，说明场地ESP（EW）的临震变化可能与裂隙水α角不稳定相关.在甘肃天祝附近，2010年10月11日发生 MS3.4（102.2°E，37.5°N）、11月20日发生MS3.6（102.9°E，37.5°N）地震，图6b是松山场地ESP（EW）、α变化曲线，在ESP（EW）剧变期间，α角也发生了间断性跳变.

在图6b和图4b中，松山场地ESP和ΔESP自2010年8月底出现剧烈变化，图7a是其3个方位地电场E在剧变前后的分钟值原始曲线.8月28日，EW方位地电场EEW出现大幅度跃变，而ENS、ENW却稳定.图7b是该台持续2年的ESP（EW）曲线，明显跃变一次发生在2009年4月8、9日，另一次始于2010年8月28日.

剪切破裂与主压应力夹角小于45°，当围压十分高时，默雷尔（Murrell）的结论是破裂面与最大主应力的夹角趋于45°（陈颙等，2009）.松山场地岩石是含水度低的碎屑岩，NS、NW方向潮汐波相关性更好，场地α角在2009年3月前基本稳定在图7（c，d）所示的 N160°E（N20°W）左右（谭大诚等，2013），基本符合该区域主压应力σ1的方位（约N20°E）（刘百篪，1979；许忠淮等，1989；张培震等，2002；崔效锋等，2006；徐纪人等，2008）.

2009年3月中旬，场地α角剧烈变化，4月底稳定在N180°E（NS方位）附近；2010年10月10日及后期多次、大幅度变化均使α角接近N90°E（EW方位）.

图7d左图内的绿色虚线标示出了2009年3月前裂隙水主体渗流α角（N20°W），则理论上σ1方位接近N25°E；2009年4月后，α角如图7d中图的绿线L1方位，场地σ1可能已偏转至N45°E附近；如果推论2010年8月28日场地岩体发生了沿另一共轭剪裂面L2方位的裂隙扩展，如图7d右图的红色虚线方位，则当天地电场E在EW方向就应剧变.由于观测曲线在NS、NW方位没有明显变化，可进一步推测场地出现了2个共轭剪裂面L1、L2并存的状况，NS方位裂隙水主体渗流受制于L1剪裂面，EW方位裂隙水主体渗流受制于L2剪裂面.

图7e是松山、古丰场地2009—2011年α角的变化曲线（取每月5、15、25日计算），两台分布、E变化曲线见图4.松山场地在2009年3、4月间，α角由N20°W偏转至近NS方位，而古丰场地这一现象则发生在2010年10月，与松山场地同时期α角的大幅变化对应.这可能反映了陈颙等（2009）指出的岩石裂隙膨胀过程存在局部性.应说明，该区域各场地近年基本都有这一趋势，但时间不一致.

局部场地σ1变化可能复杂（崔效锋等，2006），而应用潮汐地电场频谱特征及α角，可能对场地岩体裂隙水渗流、裂隙及σ1方位变化进行逐日跟踪.这对ESP变化的机理理解、应用地电场观测资料研究区域震情的发展，以及多学科的交叉分析可能重要.

图6 （a）盐源和（b）松山ESP（EW）（上）与裂隙水主体渗流方位（下）的对应变化Fig.6 Corresponding variations of ESP（EW）（upper）and preferred orientation of crack water seepage（lower）at（a）Yanyuan and（b）Songshan stations

图7 松山ESP（EW）变化的物理解析Fig.7 Physical analysis of variations of ESP（EW）at Songshan station

5 结论与展望

（1）自然电场ESP、大地电场ET可能具有同一种微观机理过程，但其不是各自的全部微观机理过程.以年尺度分析，ESP的稳定性因区域、场地及方位而异，不同场地ESP年变波形差异可能大.南北地震带附近多数场地ESP不具有清晰的年变波形，天祝台阵局部范围内近几年的年变波形不稳定.

（2）在ESP变化形态中，持续数月较平稳的升降、数日内的快速跃变、阶跃、突跳等现象常见.同一主断裂带附近、局部区域内，ESP大幅跃变、突跳或渐变等趋势转折点多具有时间同步或准同步性，同一时段数据稳定性可能类似，但形态变化存在场地因素而不具相关性.

（3）ET相对稳定的主体成分是潮汐地电场，依据其频谱特征逐日计算裂隙水主体渗流α角，这既是ET在局部场地的定量应用，也是探寻ESP变化机理的方法之一，其应用可能揭示场地岩体裂隙水渗流、裂隙、主压应力与地电场内在关联性.

总体上，本文分析的场地、数据有限，全面、深入解析地电场台网主布局区和局部场地的数据应是学科发展的需要.

致谢 感谢中国地震局地震预测研究所钱家栋研究员、中国地震局地质研究所赵国泽研究员、北京大学黄清华教授的支持；感谢甘肃、四川、新疆、宁夏、陕西等省地震局提供的地电场数据；感谢审稿专家和编辑部的意见和建议.

An Z H，Yuan L H，Li N，et al.2010.Discussion on the application of HHT method to Geo-electric field data processing.Progress in Geophysics（in Chinese），25（2）：525-532，doi：10.3969／j.issn.1004-2903.

Chen Y，Huang T F，Liu E R.2009.Rock Physics（in Chinese）.Hefei：University of Science and Technology of China Press.

Chou Y H，Liu C S，Dai Q W.2008.The Application of the Spontaneous Field Method to Earthquake Prediction（in Chinese）.Changsha：Central South University Press，16-60.

Cui X F，Xie F R，Zhang H Y.2006.Recent tectonic stress field zoning in Sichuan-Yunnan region and its dynamic interest.Acta Seismologica Sinica （in Chinese），28（5）：451-461.

Fan Y Y，Du X B，Zlotnicki J，et al.2010.The electromagnetic phenomena before the Ms8.0Wenchuan earthquake.Chinese J.Geophys.（in Chinese），53（12）：2887-2898，doi：10.3969／j.issn.0001-5733.

Geller R J.1996.Debate on evaluation of the VAN Method：Editor′s introduction.Geophys.Res.Lett.，23（11）：1291-1293.

Guo Z Q，Guo Z Q，Qian S Q，et al.1999.Electro-acoustic effect in rock fracturing.Chinese J.Geophys.（in Chinese），42（1）：74-83.

Hao J Q，Qian S Q，Gao J T，et al.2003.ULF electric and magnetic anomalies accompanying the cracking of rock sample.Acta Seismologica Sinica （in Chinese），25（1）：102-111.

Huang Q，Maggipinto T，Eftaxias K.2011.Rethinking earthquakerelated DC-ULF electromagnetic phenomena：towards a physics-based approach.Nat.Hazards Earth Syst.Sci.，11（11）：2941-2949，doi：10.5194／nhess-11-2941-2011.

Huang Q H.2005.The state-of-the-art in seismic electromagnetic observation.Recent Developments in World Seismology（in Chinese），（11）：2-5.

Huang Q H，Liu T.2006.Earthquake and tide response of geoelectric potential field at the Niijima station.Chinese J.Geophys.（in Chinese），49（6）：1745-1754.

Huang Q H.2006.Search for reliable precursors：A case study of the seismic quiescence of the 2000western Tottori prefecture earthquake.J.Geophys.Res.，111（B4）：B04301，doi：10.1029／2005JB003982.

Huang Q H，Lin Y F.2010.Numerical simulation of selectivity of seismic electric signal and its possible influences.Chinese J.Geophys.（in Chinese），53（3）：535-543，doi：10.3969／j.issn.0001-5733.

Huang Q H，Lin Y F.2010.Selectivity of seismic electric signal（SES）of the 2000Izu earthquake swarm：a 3DFEM numerical simulation model.Proc.Japan Acad.，86（3）：257-264.

Liu B C.1979.Dynamical model of the stress adjustment of the continental earthquakes in China.Seismology and Geology（in Chinese），1（3）：24-40.

Ma Q Z，Qian J D.2003.The influence of inhomogeneous geoelectric structure on the signals of geoelectric field.Earthquake（in Chinese），23（1）：1-7.

Ma Q Z，Feng Z S，Song Z P，et al.2004.Study on the variation characteristics of the geoelectric field preceding earthquakes.Acta Seismologica Sinica （in Chinese），26（3）：304-312.

Ma Q Z.2008.Multi-dipole observation system and study on the abnormal variation of the geoelectric field observed at Capital Network before the 2006Wen′an，Hebei of China，Ms5.1 earthquake.Acta Seismologica Sinica （in Chinese），30（6）：615-625.

Mao T E，Xi J L，Wang Y Q，et al.1999.The variation characteristics of the telluric field in the process of earthquake.Chinese J.Geophys.（in Chinese），42（4）：520-528.

Qian F Y，Zhao Y L.2005.Study on geoelectric field method for short-term and impending earthquake prediction.Earthquake（in Chinese），25（2）：33-40.

Qian S Q，Ren K X，LüZ.1996.Experimental study on the characteristic of VLF MF HF and VHF seismo-electromagnetic radiation with rock burst.Acta Seismologica Sinica（in Chinese），18（3）：346-351.

Ren H X，Chen X F，Huang Q H.2012.Numerical simulation of coseismic electromagnetic fields associated with seismic waves due to finite faulting in porous media.Geophys.J.Int.，188（3）：925-944.

Sun Z J，Wang H J.1984.Introduction on the Geoelectric Subject（in Chinese）.Beijing：Seismological Press.

Tan D C，Zhao J L，Xi J L，et al.2010.A study on feature and mechanism of the tidal geoelectrical field.Chinese J.Geophys.（in Chinese），53（3）：544-555，doi：10.3969／j.issn.0001-5733.

Tan D C，Wang L W，Zhao J L，et al.2011.Influence factors of harmonic waves and directional waves of tidal geoelectrical field.Chinese J.Geophys.（in Chinese），54（7）：1842-1853，doi：10.3969／j.issn.0001-5733.

Tan D C，Zhao J L，Xi J L，et al.2012.The variation of waveform and analysis of composition for the geoelectrical field before moderate or strong earthquakes in Qinghai-Tibetan plateau regions.Chinese J.Geophys.（in Chinese），55（3）：875-885，doi：10.6038／j.issn.0001-5733.

Tan D C，Xi J L，Zhang H，et al.2013.Hydrogeologic factors of geoelectric field and diurnal computation of preferred orientation of crack water seepage.Acta Seismologica Sinica（in Chinese），35（1）：36-49.

Tang J，Zhan Y，Wang L F，et al.2010.Electromagnetic coseismic effect associated with aftershock of Wenchuan Ms8.0 earthquake.Chinese J.Geophys.（in Chinese），53（3）：526-534，doi：10.3969／j.issn.0001-5733.

Tian S，Wang J G，Xu X G，et al.2009.Research on the technique of extracting seismic precursory anomaly from telluric electric field observation.Acta Seismologica Sinica （in Chinese），31（4）：424-431.

Xi J L，Song Y R，Hu M Z，et al.2013.Research on the observation methods and techniques of omni-directional spontaneous electric field.Acta Seismologica Sinica （in Chinese），35（1）：94-107.

Xu J R，Zhao Z X，Ishikawa Y Z.2008.Regional characteristics of crustal stress field and tectonic motions in and around Chinese mainland.Chinese J.Geophys.（in Chinese），51（3）：770-781.

Xu W Y.1992.Day-to-day variability of the Sq dynamo currents and Sqindex.Chinese J.Geophys.（in Chinese），35（6）：676-683.

Xu Z H，Wang S Y，Huang Y R，et al.1989.The tectonic stress field of Chinese continent deduced from a great number of earthquakes.Chinese J.Geophys.（in Chinese），32（6）：636-647.

Zhang P Z，Wang Q，Ma Z J.2002.GPS velocity field and active crustal blocks of contemporary tectonic deformation in China.EarthScienceFrontiers（in Chinese），9（2）：430-441.

Zhang X M，Zhai Y Z，Guo X Z，et al.2007.Tidal wave anomalies of geoelectrical field before remote earthquakes.Acta Seismologica Sinica （in Chinese），29（1）：48-58.

Zhao G Z，Lu J X.2003.Monitoring and analysis of earthquake phenomena by artificial SLF waves.Engineering Science（in Chinese），5（10）：27-33.

Zhao G Z，Chen X B，Xiao Q B，et al.2009.Generation mechanism of Wenchuan strong earthquake of Ms8.0inferred from EM measurements in three levels.Chinese J.Geophys.（in Chinese），52（2）：553-563.

Zhao X D，Du A M，Xu W Y，et al.2008.The origin of the prenoon-postnoon asymmetry for Sq current system.Chinese J.Geophys.（in Chinese），51（3）：643-649.

Zhou Z F，Wang J G，Huang Y.2007.Theory on Dynamics of Fluids in Fractured Medium （in Chinese）.Beijing：Higher Education Press.

附中文参考文献

安张辉，元丽华，李宁等.2010.HHT方法在地电场数据处理中的应用.地球物理学进展，25（2）：525-532，doi：10.3969／j.issn.1004-2903.

陈颙，黄庭芳，刘恩儒.2009.岩石物理学.合肥：中国科技大学出版社.

仇勇海，刘春生，戴前伟.2008.自然电场法预测地震.长沙：中南大学出版社，16-60.

崔效锋，谢富仁，张红艳.2006.川滇地区现代构造应力场分区及动力学意义.地震学报，28（5）：451-461.

范莹莹，杜学彬，Zlotnicki J等.2010.汶川Ms8.0大震前的电磁现象.地球物理学报，53（12）：2887-2898，doi：10.3969／j.issn.0001-5733.

郭自强，郭子祺，钱书清等.1999.岩石破裂中的电声效应.地球物理学报，42（1）：74-83.

郝锦绮，钱书清，高金田等.2003.岩石破裂过程中的超低频电磁异常.地震学报，25（1）：102-111.

黄清华.2005.地震电磁观测研究简述.国际地震动态，（11）：2-5.

黄清华，刘涛.2006.新岛台地电场的潮汐响应与地震.地球物理学报，49（6）：1745-1754.

黄清华，林玉峰.2010.地震电信号选择性数值模拟及可能影响因素.地球物理学报，53（3）：535-543，doi：10.3969／j.issn.0001-5733.

刘百篪.1979.中国大陆地震的应力调整场动态模型.地震地质，1（3）：24-40.

马钦忠，钱家栋.2003.地下介质非均匀结构对地电场信号的影响.地震，23（1）：1-7.

马钦忠，冯志生，宋治平等.2004.崇明与南京台震前地电场变化异常分析.地震学报，26（3）：304-312.

马钦忠.2008.地电场多极距观测装置系统与文安Ms5.1地震前首都圈地电场异常研究.地震学报，30（6）：615-625.

毛桐恩，席继楼，王燕琼等.1999.地震过程中的大地电场变化特征.地球物理学报，42（4）：520-528.

钱复业，赵玉林.2005.地电场短临预报方法研究.地震，25（2）：33-40.

钱书清，任克新，吕智.1996.伴随岩石破裂VLF，MF，HF和VHF电磁辐射特性的实验研究.地震学报，18（3）：346-351.

孙正江，王华俊.1984.地电概论.北京：地震出版社.

谭大诚，赵家骝，席继楼等.2010.潮汐地电场特征及机理研究.地球物理学报，53（3）：544-555，doi：10.3969／j.issn.0001-5733.

谭大诚，王兰炜，赵家骝等.2011.潮汐地电场谐波和各向波形的影响要素.地球物理学报，54（7）：1842-1853，doi：10.3969／j.issn.0001-5733.

谭大诚，赵家骝，席继楼等.2012.青藏高原中强地震前的地电场变异及构成解析.地球物理学报，55（3）：875-885，doi：10.6038／j.issn.0001-5733.

谭大诚，席继楼，张慧等.2013.地电场水文地质因素及裂隙水主体渗流方向逐日计算.地震学报，35（1）：36-49.

汤吉，詹艳，王立风等.2010.汶川地震强余震的电磁同震效应.地球物理学报，53（3）：526-534，doi：10.3969／j.issn.0001-5733.

田山，王建国，徐学恭等.2009.大地电场观测地震前兆异常提取技术研究.地震学报，31（4）：424-431.

席继楼，宋艳茹，胡明朝等.2013.全方位自然电场观测方法和观测技术研究.地震学报，35（1）：94-107.

徐纪人，赵志新，石川有三.2008.中国大陆地壳应力场与构造区域特征研究.地球物理学报，51（3）：770-781.

徐文耀.1992.Sq发电机电流的逐日变化和Sq指数.地球物理学报，35（6）：676-683.

许忠淮，汪素云，黄雨蕊等.1989.由大量的地震资料推断的我国大陆构造应力场.地球物理学报，32（6）：636-647.

张培震，王琪，马宗晋.2002.中国大陆现今构造运动的GPS速度场与活动地块.地学前缘，9（2）：430-441.

张学民，翟彦忠，郭学增等.2007.远震前的地电场潮汐波异常.地震学报，29（1）：48-58.

赵国泽，陆建勋.2003.利用人工源超低频电磁波监测地震的试验与分析.中国工程科学，5（10）：27-33.

赵国泽，陈小斌，肖骑彬等.2009.汶川Ms8.0级地震成因三“层次”分析——基于深部电性结构.地球物理学报，52（2）：553-563.

赵旭东，杜爱民，徐文耀等.2008.Sq电流系午前午后不对称性现象的来源.地球物理学报，51（3）：643-649.

周志芳，王锦国，黄勇.2007.裂隙介质水动力学原理.北京：高等教育出版社.