A method for locating interference source by long-and-short polar distance ratio and its practical application in geoelectric field observation
-
摘要: 提出长短极距比值干扰源定位方法用来定位地电观测中出现的地表干扰源。该方法基于均匀半空间模型且地表干扰为点电源，利用地表不同位置的干扰源所引起的同测向长短极距的电场变化差异，联合多测向比值来定位干扰源位置，并应用到徐庄子台漏电干扰一例和高邮台偶极源试验中。结果表明：利用长短极距比值法，能有效定位干扰点源，且能定位偶极干扰源的其中一个电极。Abstract: The long and short pole distance ratio method in geoelectric observation is employed here to locate the near-surface disturbance sources, on the assumption of a uniform half-space model and that the near-surface disturbance is a point source. The disturbance sources of different locations can cause different increments in the long-pole-distance and short-pole-distance geoelectric field at the same direction, which is the key mind of this method. By combining multi-directional long-and-short pole distance ratios in geoelectric field, we successfully applied this method into locating the disturbance sources of the leakage at Xuzhuangzi station and the dipole source test at Gaoyou station. The results show that the point source producing the interference in geoelectric field at Xuzhuangzi station can be effectively located, and one electrode of the dipole source can be discovered in the Gaoyou station test as well.
-
引言
20世纪六、七十年代，美国、苏联、日本和我国开始研究地震地电学，研究表明电磁现象与地震事件存在一定相关性（Miyakoshi，1986），近年来理论方法方面得到快速发展，而技术方面并没有明显突破（Helman，2020）。我国从20世纪末期开始逐步建设地电场观测台站，目前已经发展为国际上少有的拥有大规模、规范化的地电场观测网的国家，可以实现全面覆盖、时间上连续以及地震危险区重点监测的地电场观测（中国地震学会，2007）。迄今为止，已经记录了丰富的地电场观测资料，促进了利用地电场优势方位角等方法进行地震预报的研究（谭大诚等，2014）。但随着社会的发展，地电场观测台站附近的工农业用电环境越来越复杂，用电过程中若发生对地漏电，势必对地电场观测造成干扰，在原始数据上可能表现为较大的扰动、突变和失真，进而可能影响地电场观测数据的分析和应用。
在地电场观测过程中，常见的干扰有电极极化、雷击、高压直流输电以及测区内外漏电等干扰（李伟等，2014）。当漏电干扰发生时，必须找到干扰源才能消除该干扰，因此如何在短时间内对电磁干扰源进行定位，降低漏电发生对地电场观测的影响，在地电场干扰调查中至关重要。目前，国际上还没有研究者针对干扰源对地电场观测的影响进行专门分析，而国内地电场观测中漏电干扰的定位研究较少，查阅资料仅发现两例实际漏电干扰定位的研究：陈志刚等（2012）利用干扰幅度的矢量合成初步定位干扰源的方位和距离，再结合现场排查来定位干扰源。其方法的优势为应用简便，但该方法的计算结果与现场排查找到的干扰源距离存在较大误差；王福才等（2013）和蒋延林等（2015）结合实地调查、重点走访和分片区停电来定位干扰源。该方法能确保找到干扰源，但需要多部门协调配合，费时费力，难以推广应用。建立偶极源模型进行定位试验的有两例：卫清等（2017）利用多元牛顿迭代算法对漏电位置进行反演，张宇等（2017）提出动态枚举定位求解算法。这两种方法能定位偶极干扰源的两个电极，且精度高，误差仅为数米。但这两个方法的使用有较为严格的限制，要求一个共用电极和另外五个电极，多数地电场观测台站并未与地电阻率测量同场地观测，这种情形下这两种方法均不可用，且这两种方法目前都仅限于试验阶段。
VAN法认为可以用同方向不同极距的地电场比值是否为1来区别近场信号和地震电磁场信号（Varotsos，Alexopoulos，1984）。田山等（2009）认为可以利用长短极距比值的稳定性来分析地震地电场异常信息。这些研究表明同测向长短极距比值具有可应用性，实际地电场观测中，通常会利用同测向长短极距比值是否稳定来分析是否有近场干扰（马钦忠等，2016）。由于不同位置的点源所引起的干扰电场的长短极距比值大小不同，本文提出长短极距比值干扰源定位方法，将漏电干扰源视为地表点源，利用不同测向的比值就可以定位出点源的具体位置。利用该方法定位了徐庄子台漏电干扰源，验证了该方法在地电场观测中能准确定位可视为点源的漏电干扰源，并利用高邮地震台人工干扰源试验观测数据定位干扰源位置，验证了当存在距离近且极距大的偶极源干扰时，能准确定位其中一个电极。
1. 基本原理
我国在进行地电场观测时，通常采用多方向、多极距测量方式，一般为三个测向：东西、南北、北东或北西，在每个观测方向布设长、短两个极距的测道，一共设置6个测道，台站布极如图1所示。
![]()
图 1 我国地电场观测布极方式示意图（a） 十字形布设；（b） L形布设Figure 1. Schematic diagram of the electrode arrangement for our geoelectric field observation（a） Cruciform layout；（b） L-shaped layout建立电阻率均匀各向同性的半空间模型，图2左侧为漏电干扰源与测量电极坐标示意图，xoy平面直角坐标系为地表，其中：红点表示电量为q的干扰源，蓝点表示测量电极，以o点为坐标原点，向东为x轴正向，向北为y轴正向建立坐标系。
![]()
图 2 不同的电极布设坐标示意图（左）、地电场长短极距比值分布图（中）和比值确定时的等值线图（右）（a） 十字形东西向；（b） L形东西向；（c） L形北东向Figure 2. Schematic diagram of different electrode layout coordinates （left），distribution of the long and short pole distance ratio for geoelectric field （middle），and contour map when the ratio is determined （right）（a） Cross-shaped east-west direction；（b） L-shaped east-west direction； （c） L-shaped north-east direction干扰源在测量电极处产生的电势为：
$$ U ＝ \frac{{kq}}{r}，$$ (1) 式中，$k$为在均匀半无限空间中的比例系数，为常数，其值与均匀介质的介电常数相关。
电场强度的计算公式为：
$$ E ＝ \frac{{{U_1} － {U_2}}}{L}\text{，} $$ (2) 式中，${U_1}$为干扰源在测量电极1处产生的电势，${U_2}$为干扰源在测量电极2处产生的电势，$L$为两测量电极间的距离。将式（1）代入式（2）有：
$$ E ＝ \left(\frac{{kq}}{{{r_1}}} － \frac{{kq}}{{{r_2}}}\right) {\text{×}} \frac{1}{L} \text{，} $$ (3) 式中${r_1}$和${r_2}$分别为干扰源到两个测量电极的距离。
将长极距的测量电极坐标$ （ {x_{{\rm{L}}1}}， {y_{{\rm{L}}1}} ） $，$ （ {x_{{\rm{L}}2}}， {y_{{\rm{L}}2}}$）和干扰源坐标$ （ {x}_{0}， {y}_{0}$）带入式（3）：
$$ {E_{\rm{L}}} ＝ \left( {\frac{{kq}}{{\sqrt {{{ （ {x_0} － {x_{{\rm{L}}1}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{L}}1}} ） }^2}} }} － \frac{{kq}}{{\sqrt {{{ （ {x_0} － {x_{{\rm{L}}2}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{L}}2}} ） }^2}} }}} \right) {\text{×}} \frac{1}{{\sqrt {{{ （ {x_{{\rm{L}}1}} － {x_{{\rm{L}}2}} ） }^2} ＋ {{ （ {y_{{\rm{L}}1}} － {y_{{\rm{L}}2}} ） }^2}} }} ； $$ (4) 将短极距的测量电极坐标$ （ {x_{{\rm{S}}1}}， {y_{{\rm{S}}1}} ） $、$ （ {x_{{\rm{S}}2}}， {y_{{\rm{S}}2}} ） $和干扰源坐标$ （ {x}_{0}， {y}_{0} $）带入式（3）：
$$ {E_{\rm{S}}} ＝ \left( {\frac{{kq}}{{\sqrt {{{ （ {x_0} － {x_{{\rm{S}}1}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{S}}1}} ） }^2}} }} － \frac{{kq}}{{\sqrt {{{ （ {x_0} － {x_{{\rm{S}}2}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{S}}2}} ） }^2}} }}} \right) {\text{×}} \frac{1}{{\sqrt {{{ （ {x_{{\rm{S}}1}} － {x_{{\rm{S}}2}} ） }^2} ＋ {{ （ {y_{{\rm{S}}1}} － {y_{{\rm{S}}2}} ） }^2}} }} \text{.} $$ (5) 计算上述长极距测得的电场值和短极距测得的电场值的比值：
$$ \dfrac{{{E_{\rm{L}}}}}{{{E_{\rm{S}}}}} ＝ \dfrac{{\left( {\dfrac{1}{{\sqrt {{{ （ {x_0} － {x_{{\rm{L}}1}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{L}}1}} ） }^2}} }} － \dfrac{1}{{\sqrt {{{ （ {x_0} － {x_{{\rm{L}}2}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{L}}2}} ） }^2}} }}} \right) {\text{×}} \dfrac{1}{{\sqrt {{{ （ {x_{{\rm{L}}1}} － {x_{{\rm{L}}2}} ） }^2} ＋ {{ （ {y_{{\rm{L}}1}} － {y_{{\rm{L}}2}} ） }^2}} }}}}{{\left( {\dfrac{1}{{\sqrt {{{ （ {x_0} － {x_{{\rm{S}}1}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{S}}1}} ） }^2}} }} － \dfrac{1}{{\sqrt {{{ （ {x_0} － {x_{{\rm{S}}2}} ） }^2} ＋ {{ （ {y_0} － {y_{{\rm{S}}2}} ） }^2}} }}} \right) {\text{×}} \dfrac{1}{{\sqrt {{{ （ {x_{{\rm{S}}1}} － {x_{{\rm{S}}2}} ） }^2} ＋ {{ （ {y_{{\rm{S}}1}} － {y_{{\rm{S}}2}} ） }^2}} }}}} \text{.} $$ (6) 从式（6）可以看出，地电场长短极距比值EL/ES与干扰源的电量q无关，仅与各测量电极的坐标和干扰源的坐标$ （ {x_0}， {y_0} ） $有关。地电场台站中各测向的长短极距测量电极的坐标固定，所以只需知道干扰源的坐标，将其代入式（6）就可以计算出该干扰源引起的各测向地电场长短极距增量的比值。
我们将o点作为坐标原点，正北为y轴正方向，正东为x轴正方向，建立坐标系，固定长短极距测量电极的坐标。以图2a为例，此时东西向长极距测量电极A1，A2和短极距测量电极坐标A1 ′，A2 ′已知，那么某一处存在干扰源时，将该干扰源坐标带入式（6），便可计算干扰引起的地电场长短极距增量的比值。将原点附近2 km范围内按4 m间隔进行网格划分，计算仅每个网格点处存在干扰源时引起的东西向长短极距测道所产生的增量的比值，并将比值以不同颜色表示，得到图2a所示的地电场长短极距比值分布图。我们采用Python语言编写该程序，计算并绘制此图用时1.27 s，说明该方法计算速度快。由于该方法在干扰源的计算公式推导过程中时并未假设长短极距测道同方向，所以该方法并不仅限于同方向的两测道的比值计算，不同方向不同测道的数据也可用该方法计算干扰源的位置。当漏电干扰出现时，我们选取干扰前后的数据作差，分别计算干扰导致长极距和短极距的地电场增量，然后计算地电场长短极距增量的比值，在分布图中只画出该比值的等值线，那么该等值线上每一个坐标点都是干扰源可能的坐标点。
图2a左图为十字形东西向长短极距电极布设方式坐标示意图，图2a中间图是将图2a左图中的坐标点带入式（6）计算得到的十字形东西向地电场长短极距比值分布图，可以看出比值关于x轴和y轴对称，在两条45°对角线上比值接近于1，2 km以外存在干扰源时，比值接近于1。当确定比值为1.2时，即从图2a中间图中仅选出值为1.2的等值线（图2a右图），图中红线上每一个坐标点都可能是干扰源所在的坐标点。
图2b显示：地电场长短极距比值关于x轴轴对称，在两条30—40°的对角线上比值接近于1。当确定比值为0.8时，即从图2b中间图仅选出值为0.8的等值线（图2b右图），图中红线上每一个坐标点都可能是干扰源所在的坐标点。
图2c显示：地电场长短极距比值关于北西对角线对称，在北东对角线上比值接近于1。当确定比值为0.9时，即从图2c中间图仅选出值为0.9的等值线（图2c右图），图中蓝线上每一个坐标点都可能是干扰源所在的坐标点。
实际地电场观测中发生漏电干扰时，我们计算不同测向的地电场长短极距比值，将不同测线计算得到的可能的坐标点画在同一个图中，交点便是干扰源可能存在的位置。假设地电场观测布极方式为如图2b和2c的L形布设，且东西向比值为0.8，北东向比值为0.9，将两测向的等值线画在一个坐标系中（图3），可以看到有两个交点，那么这两个交点即为干扰源可能的位置。
![]()
图 3 长短极距比定位漏电干扰源示意图Figure 3. Schematic diagram of locating leakage interference source by long-and-short polar distance ratio从图3可以看出交点可能不止一个，那么需要利用以下两个方法排除不可能的坐标点：① 通过观察其它测向是否有干扰来排除近处的交点。② 结合各测向的地电场干扰的增减同步性来判断，即南北向和东西向地电场值同增同减时，干扰源在两测线所夹的象限或者斜对面的象限；南北向和东西向地电场值增减相反时，干扰源在另外两象限。
漏电电流分为直流电和交流电，本文仅讨论直流漏电干扰。在理想状况下，漏电发生时，电场表现为方波状，即：漏电前，电场保持不变；漏电时刻，电场突增；漏电持续时，电场维持不变；漏电消失时刻，电场突降。方波从低值到高值的增量是漏电干扰产生的电场值。但由于数据采集只能采集一定频率范围内的信号，漏电前后会出现吉布斯现象，即在方波脉冲的转折处会有明显的震荡。除了震荡的转折处，在其它时刻我们可以得到稳定的低值和高值，所以为了提高结果的准确性，应该选数据相对稳定期间的低值和高值来计算增量并以此作为漏电干扰产生的电场值。日常地电场数据采集过程中，分钟值是一分钟之内进行多次采样后平均得来，虽然在一定程度上消除了漏电前后吉布斯现象的影响，但由于一分钟之内地电场变化较大，数据表现出不稳定性，在干扰出现前后，无法确定稳定的高值和低值，所以选择不同时刻的数据会导致干扰源的定位结果不同。我们认为应该通过提高采样率来避免地电场变化对干扰源定位的影响，并选择震荡转折以外稳定的高值和低值来计算增量，以避免吉布斯现象的影响，从而提高该方法的准确性。
2. 点干扰源的定位应用
徐庄子台地电场电极布设方式为L型，各测量电极位置如图4所示，建立以o点为坐标原点，向东为x轴正向，向北为y轴正向的坐标轴，东西向测线和南北向测线所夹的象限为第三象限。2009年11月，徐庄子台的地电场观测数据出现方波干扰。图5为2010年11月16日徐庄子台地电场分钟数据对比曲线，可以看出在18时和21时左右，南北向和东西向同时发生了方波状漏电干扰，且长短极距测得的地电场干扰幅度成一定比值，而北西向没有明显干扰。
![]()
图 4 徐庄子地震台地电场观测系统电极布置方式示意图 （单位：m）Figure 4. The sketch map of layout patterns for electrodes in geoelectric field observatory systems at Xuzhuangzi seismic station （unit：m）![]()
图 5 徐庄子台2010年11月16日地电场分钟数据对比曲线Figure 5. Comparison curves of minute data of geoelectric field at Xuzhuangzi station on November 16，201011月出现了多次干扰事件，每次干扰事件各测向的长短极距增量比值差别不大，表现在干扰源坐标上的差别仅为几米，可以忽略不计，因此本文直接采用了陈志刚等（2012）在文章中随机挑选的11月中的5次干扰事件的东西向和南北向的长短极距的地电场值（表1），分别计算了这5次干扰事件的长短极距比，并求得它们的平均值。由于北西向地电场数据未表现出明显干扰，文中只使用了南北和东西测向的值。
表 1 徐庄子台干扰电场长短极距比Table 1. Long-and-short polar distance ratio of interference electric field at Xuzhuangzi station利用东西向和南北向长短极距比的平均值画图，如图6，黑线和红线有五个交叉点，因此需要结合其它条件进行排除。东西向测线和南北向测线所夹的象限为第三象限，从图5可以看出，东西向和南北向的电场干扰同增同减，干扰源应在两测线之间象限或者斜对面的象限，即第三象限或第一象限。因北西向地电场不受干扰，最终认定干扰源位于红圈所示的这一交叉点，其坐标为（250，400），方位为北偏东32°，距离o点471.7 m。最后经实地确认，干扰源为台站东北方、距离台站500 m的木材加工厂的电器供暖设备漏电，与文中计算得到的471.7 m相差约30 m，且方位与最终排查得到的干扰源方位相近。
![]()
图 6 徐庄子台干扰源定位结果示意图Figure 6. Schematic diagram of the location result of the interference source at Xuzhuangzi station陈志刚等（2012）的文章中利用矢量合成方法计算的方位角为北偏东43°，本文计算结果为北偏东32°。我们认为利用矢量合成的方法初步确定干扰源的方位和距离不可靠：一方面，出现干扰源时地电场值是有增有减的，不能因此推论干扰源在东北方向；另一方面，漏电干扰为近源，徐庄子台的电极布设方式为L型，不能简单的将东西向和南北向的干扰电场近似为o点的电场。
3. 偶极干扰源的定位应用
2016年6月卫清等（2017）和张宇等（2017）在高邮台地电观测场地进行人为对地供电试验，模拟偶极源漏电对地电场观测的干扰，用来验证牛顿迭代反演法和动态枚举法是否能准确定位干扰源。由于干扰源离测量电极很近、且极距大，因此不可看作点源，需要找出两个电极所在的位置。我们利用该文章中的数据研究长短极距比值法在偶极干扰源试验中的计算效果。图7为高邮地电观测布极图，建立以o点为坐标原点，向东为x轴正向，向北为y轴正向的坐标轴。将地电阻率供电电极${A_1}$和${B_1}$设为偶极干扰源，东西测向的长极距测道为a1o、短极距测道为a2o；南北测向的长极距测道为b1o、短极距测道为b2o；斜测向的长极距测道为M3o、短极距测道为N1o。我们需要利用这三个测向的长短极距比值来定位干扰源${A_1}$和${B_1}$的坐标。
![]()
图 7 高邮地震台地电场及地电阻率同场观测电极布极图（单位：m）Figure 7. The sketch map of layout patterns for electrodes in geoelectric field and earth resistivity observatory systems at Gaoyou seismic station （unit：m）利用张宇等（2017）文章中提供的人工干扰源试验电位差值ΔV和对应测道的极距长度L带入电场强度公式E＝ΔV/L，计算得到长短极距的地电场及其比值（表2）。由于是试验数据，采样率高，数据质量高，电场值随时间的变化呈较为标准的方波形态，选取方波从低值到高值突变时刻的数据。从三次试验的各测向长短极距比值大小基本一致可以看出，在采样率高、数据质量高的情况下，长短极距比值稳定也证明了长短极距比值能反应干扰源的位置。
表 2 偶极源试验中干扰电场长短极距比Table 2. Long-and-short polar distance ratio of interference electric field in dipole source test图8为将表2中的三个测向的长短极距比平均值的等值线画在一幅图中，其中三个交点为干扰源可能的坐标点。实际的${A_1}$电极的坐标为（303，−420），计算得到的坐标为（291，−410），相差约20 m。实际的${B_1}$电极的坐标为（303，579），该电极在另外两个交点（288，654）和（306，182）之间，说明该方法用于计算偶极干扰源时也能有效定位其中一个电极。也说明在地电阻率与地电场同地观测的情况下，出现漏电干扰时，可以记录地电阻率法的测量电极的电位值，计算长短极距比作为另一测向数据，从而更准确地定位干扰源。
![]()
图 8 偶极干扰源试验定位结果示意图Figure 8. Schematic diagram of the location results of the dipole interference source test4. 结论
尽管对地漏电在地电场观测过程中并不常见，但漏电一旦发生，将严重影响地电场观测数据，所以快速定位干扰源，排除干扰，对地电场观测意义重大。本工作提出长短极距比值法定位电磁干扰源，并应用到两个实例中，得到了较好的应用效果。
本方法的实际应用结果表明：① 当干扰源极距较小可以看作点源时，计算得到的干扰源位置较为准确，可以为排查干扰节省成本；② 当干扰源距离测量电极很近，且干扰源极距较大时，干扰源不可看作点源，该方法可以有效定位其中一个电极，再通过现场排查来确定干扰源的两极。长短极距比值法利用数据库提供的各测向的地电场分钟值就能定位漏电干扰源，应用简便，投入实际应用中可以大大缩小排查范围。
-
![]()
图 1 我国地电场观测布极方式示意图
（a） 十字形布设；（b） L形布设
Figure 1. Schematic diagram of the electrode arrangement for our geoelectric field observation
（a） Cruciform layout；（b） L-shaped layout
![]()
图 2 不同的电极布设坐标示意图（左）、地电场长短极距比值分布图（中）和比值确定时的等值线图（右）
（a） 十字形东西向；（b） L形东西向；（c） L形北东向
Figure 2. Schematic diagram of different electrode layout coordinates （left），distribution of the long and short pole distance ratio for geoelectric field （middle），and contour map when the ratio is determined （right）
（a） Cross-shaped east-west direction；（b） L-shaped east-west direction； （c） L-shaped north-east direction
![]()
图 3 长短极距比定位漏电干扰源示意图
Figure 3. Schematic diagram of locating leakage interference source by long-and-short polar distance ratio
![]()
图 4 徐庄子地震台地电场观测系统电极布置方式示意图 （单位：m）
Figure 4. The sketch map of layout patterns for electrodes in geoelectric field observatory systems at Xuzhuangzi seismic station （unit：m）
![]()
图 5 徐庄子台2010年11月16日地电场分钟数据对比曲线
Figure 5. Comparison curves of minute data of geoelectric field at Xuzhuangzi station on November 16，2010
![]()
图 6 徐庄子台干扰源定位结果示意图
Figure 6. Schematic diagram of the location result of the interference source at Xuzhuangzi station
![]()
图 7 高邮地震台地电场及地电阻率同场观测电极布极图（单位：m）
Figure 7. The sketch map of layout patterns for electrodes in geoelectric field and earth resistivity observatory systems at Gaoyou seismic station （unit：m）
![]()
图 8 偶极干扰源试验定位结果示意图
Figure 8. Schematic diagram of the location results of the dipole interference source test
表 1 徐庄子台干扰电场长短极距比
Table 1 Long-and-short polar distance ratio of interference electric field at Xuzhuangzi station
下载: 导出CSV
表 2 偶极源试验中干扰电场长短极距比
Table 2 Long-and-short polar distance ratio of interference electric field in dipole source test
下载: 导出CSV
-
陈志刚,田山,徐学恭,马朝晖,陈嵩. 2012. 地电场观测漏电干扰的排查一例[J]. 华北地震科学,30(4):44–48. Chen Z G,Tian S,Xu X G,Ma Z H,Chen S. 2012. Study on electricity-leaking interference in geoelectric field observation at Xuzhuangzi seismic station[J]. North China Earthquake Sciences,30(4):44–48 (in Chinese).
蒋延林,赵卫红,王福才,张骞,朱晔. 2015. 一种地电场观测环境干扰的分析和试验研究[J]. 中国地震,31(1):158–164. Jiang Y L,Zhao W H,Wang F C,Zhang Q,Zhu Y. 2015. Analysis and experimental study of the environment interference in the geo-electric field observation[J]. Earthquake Research in China,31(1):158–164 (in Chinese).
李伟,马钦忠,王冠玥. 2014. 地电场多极距观测的异常识别和分析[J]. 中国地震,30(1):91–101. Li W,Ma Q Z,Wang G Y. 2014. Abnormal signal recognition by multi-dipole observation of geoelectric field[J]. Earthquake Research in China,30(1):91–101 (in Chinese).
马钦忠,钱家栋,李伟,赵文舟,方国庆. 2016. 源自多个大电流源的华东地区地电场空间变化特征[J]. 地球物理学报,59(7):2598–2614. Ma Q Z,Qian J D,Li W,Zhao W Z,Fang G Q. 2016. Characteristics of the spatial variation of geoelectric field signals recorded at the stations in Huadong area in China when 4 heavy currents are injected[J]. Chinese Journal of Geophysics,59(7):2598–2614 (in Chinese).
谭大诚,赵家骝,刘小凤,范莹莹,刘君,陈军营. 2014. 自然电场的区域性变化特征[J]. 地球物理学报,57(5):1588–1598. Tan D C,Zhao J L,Liu X F,Fan Y Y,Liu J,Chen J Y. 2014. Features of regional variations of the spontaneous field[J]. Chinese Journal of Geophysics,57(5):1588–1598 (in Chinese).
田山,王建国,徐学恭,董洪军,崔晓峰,马骥. 2009. 大地电场观测地震前兆异常提取技术研究[J]. 地震学报,31(4):424–431. Tian S,Wang J G,Xu X G,Dong H J,Cui X F,Ma J. 2009. Research on the technique of extracting seismic precursory anomaly from telluric electric field observation[J]. Acta Seismologica Sinica,31(4):424–431 (in Chinese).
王福才,蒋延林,张骞,赵卫红,朱晔,薛家富. 2013. 地电观测场地干扰的一种测试查找方法[J]. 西北地震学报,35(增刊):126–134. Wang F C,Jiang Y L,Zhang Q,Zhao W H,Zhu Y,Xue J F. 2013. A test and search method for the interferences in geoelectric observation field[J]. China Earthquake Engineering Journal,35(S1):126–134 (in Chinese).
卫清,颜蕊,韩秀红,温晋. 2017. 基于牛顿迭代法的地电场观测人为漏电干扰源反演及其应用[J]. 地球物理学进展,32(4):1490–1495. Wei Q,Yan R,Han X H,Wen J. 2017. Inversion for artificial electricity-leaking in geo-electric field observations based on Newton iteration and its application[J]. Progress in Geophysics,32(4):1490–1495 (in Chinese).
张宇,王兰炜,赵家骝,蒋延林,胡哲,张兴国,王福才. 2017. 地电观测中电磁干扰源定位方法研究[J]. 地震学报,39(3):367–373. Zhang Y,Wang L W,Zhao J L,Jiang Y L,Hu Z,Zhang X G,Wang F C. 2017. Location method of electromagnetic interference sources in geo-electric field observation[J]. Acta Seismologica Sinica,39(3):367–373 (in Chinese).
中国地震学会. 2007. 中国地震科技事业的发展[J]. 国际地震动态,(8):23–37. doi: 10.3969/j.issn.0253-4975.2007.08.010 Seismological Society of China. 2007. The development of the Chinese seismological science and technology[J]. Recent Developments in World Seismology,(8):23–37 (in Chinese).
Helman D S. 2020. Seismic electric signals (SES) and earthquakes:A review of an updated VAN method and competing hypotheses for SES generation and earthquake triggering[J]. Phys Earth Planet Inter,302:106484. doi: 10.1016/j.pepi.2020.106484
Miyakoshi J. 1986. Anomalous time variation of the self-potential in the fractured zone of an active fault preceding the earthquake occurrence[J]. J Geomagn Geoelectr,38(10):1015–1030. doi: 10.5636/jgg.38.1015
Varotsos P,Alexopoulos K. 1984. Physical properties of the variations of the electric field of the earth preceding earthquakes,I[J]. Tectonophysics,110(1/2):73–98.
-
期刊类型引用(1)
1. 郭明瑞，张宇，张兴国，王锡娇，张慧，曾涛，杨世杰，卢启明，庞雪华. 琼中地电场观测干扰排查实验分析. 地震地磁观测与研究. 2024(05): 98-106 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 317
- HTML全文浏览量: 93
- PDF下载量: 64
- 被引次数: 1
